Model-based detection and analysis of introgressed Neanderthal ancestry in modern humans, bioRxiv, 2017-12-02

AbstractGenetic evidence has revealed that the ancestors of modern human populations outside of Africa and their hominin sister groups, notably the Neanderthals, exchanged genetic material in the past. The distribution of these introgressed sequence-tracts along modern-day human genomes provides insight into the ancient structure and migration patterns of these archaic populations. Furthermore, it facilitates studying the selective processes that lead to the accumulation or depletion of introgressed genetic variation. Recent studies have developed methods to localize these introgressed regions, reporting long regions that are depleted of Neanderthal introgression and enriched in genes, suggesting negative selection against the Neanderthal variants. On the other hand, enriched Neanderthal ancestry in hair- and skin-related genes suggests that some introgressed variants facilitated adaptation to new environments. Here, we present a model-based method called diCal-admix and apply it to detect tracts of Neanderthal introgression in modern humans. We demonstrate its efficiency and accuracy through extensive simulations. We use our method to detect introgressed regions in modern human individuals from the 1000 Genomes Project, using a high coverage genome from a Neanderthal individual from the Altai mountains as reference. Our introgression detection results and findings concerning their functional implications are largely concordant with previous studies, and are consistent with weak selection against Neanderthal ancestry. We find some evidence that selection against Neanderthal ancestry was due to higher genetic load in Neanderthals, resulting from small effective population size, rather than Dobzhansky-Müller incompatibilities. Finally, we investigate the role of the X-chromosome in the divergence between Neanderthals and modern humans.

biorxiv evolutionary-biology 0-100-users 2017

Partially redundant actin genes in Chlamydomonas control flagellum-directed traffic and transition zone organization, bioRxiv, 2017-12-01

ABSTRACTFlagella of the unicellular green alga Chlamydomonas reinhardtii are nearly identical to cilia of mammalian cells and provide an excellent model to study ciliogenesis. These biflagellated cells have two actin genes one encoding a conventional actin (IDA5) and the other encoding a divergent novel actin-like protein (NAP1). Previously, we described a role for actin in the regulation of flagella-building intraflagellar transport machinery. Here, we probe how actin redundancy contributes to this process using a nap1 mutant Chlamydomonas strain. Disruption of a single actin allows normal or slower incorporation but complete flagellar assembly. However, when we disrupt both actins using Latrunculin B (LatB) treatment on the nap1 mutant background, we find flagellar growth from newly synthesized limiting flagellar proteins is actin-dependent. Upon total actin disruption during flagellar assembly, transmission electron microscopy identified an accumulation of Golgi-adjacent vesicles, suggesting impaired vesicular trafficking may be the mechanism by which actin supports flagellar growth from new flagellar proteins. We also find there is a mislocalization of a key transition zone gating and ciliopathy protein, NPHP-4. Extended (2 hour) treatment with LatB, a condition under which NAP1 is upregulated, restores NPHP-4 localization. This suggests NAP1 can perform the functions of conventional actin at the transition zone. Our experiments demonstrate that each stage of flagellar biogenesis requires redundant actin function to varying degrees, with an absolute requirement for these actins in transport of Golgi-adjacent vesicles and flagellar incorporation of newly synthesized proteins.

biorxiv cell-biology 0-100-users 2017

Biomolecular simulations under realistic macroscopic salt conditions, bioRxiv, 2017-11-30

Biomolecular simulations are typically performed in an aqueous environment where the number of ions remains fixed for the duration of the simulation, generally with either a minimally neutralizing ion environment or a number of salt pairs intended to match the macroscopic salt concentration. In contrast, real biomolecules experience local ion environments where the salt concentration is dynamic and may differ from bulk. The degree of salt concentration variability and average deviation from the macroscopic concentration remains, as yet, unknown. Here, we describe the theory and implementation of a Monte Carlo osmostat that can be added to explicit solvent molecular dynamics or Monte Carlo simulations to sample from a semigrand canonical ensemble in which the number of salt pairs fluctuates dynamically during the simulation. The osmostat reproduce the correct equilibrium statistics for a simulation volume that can exchange ions with a large reservoir at a defined macroscopic salt concentration. To achieve useful Monte Carlo acceptance rates, the method makes use of nonequilibrium candidate Monte Carlo (NCMC) moves in which monovalent ions and water molecules are alchemically transmuted using short nonequilibrium trajectories, with a modified Metropolis-Hastings criterion ensuring correct equilibrium statistics for an (Δ𝜇, 𝑁, 𝑝, 𝑇) ensemble. We demonstrate how typical protein (DHFR and the tyrosine kinase Src) and nucleic acid (Drew-Dickerson B-DNA dodecamer) systems exhibit salt concentration distributions that significantly differ from fixed-salt bulk simulations and display fluctuations that are on the same order of magnitude as the average.

biorxiv biophysics 0-100-users 2017

 

Created with the audiences framework by Jedidiah Carlson

Powered by Hugo