Neutral Drift Research Articles

Homeologue-specific expression divergence in the recently formed tetraploid Capsella bursa-pastoris (Brassicaceae).

Following allopolyploid formation, extensive genome evolution occurs, with the eventual loss of many homeologous gene copies. Although this process of diploidization has occurred many times independently, the evolutionary forces determining the probability and rate of gene loss remain poorly understood. Here, we conduct genome and transcriptome sequencing in a broad sample of Chinese accessions of Capsella bursa-pastoris, a recently formed allotetraploid. Our whole genome data reveal three groups of these accessions: an Eastern group from low-altitude regions, a Western group from high-altitude regions, and a much more differentiated Northwestern group. Population differentiation in total expression was limited among closely related populations; by contrast, the relative expression of the two homeologous copies closely mirrors the genome-wide SNP divergence. Consistent with this, we observe a negative correlation between expression changes in the two homeologues. However, genes showing population genomic evidence for adaptive evolution do not show an enrichment for expression divergence between homeologues, providing no clear evidence for adaptive shifts in relative gene expression. Overall, these patterns suggest that neutral drift may contribute to the population differentiation in the expression of the homeologues, and drive eventual gene loss over longer periods of time.

Shuffling the Neutral Drift of Unspecific Peroxygenase in Saccharomyces cerevisiae.

Unspecific peroxygenase (UPO) is a highly promiscuous biocatalyst, and its selective mono(per)oxygenase activity makes it useful for many synthetic chemistry applications. Among the broad repertory of library creation methods for directed enzyme evolution, genetic drift allows neutral mutations to be accumulated gradually within a polymorphic network of variants. In this study, we conducted a campaign of genetic drift with UPO in Saccharomyces cerevisiae, so that neutral mutations were simply added and recombined in vivo With low mutational loading and an activity threshold of 45% of the parent's native function, mutant libraries enriched in folded active UPO variants were generated. After only eight rounds of genetic drift and DNA shuffling, we identified an ensemble of 25 neutrally evolved variants with changes in peroxidative and peroxygenative activities, kinetic thermostability, and enhanced tolerance to organic solvents. With an average of 4.6 substitutions introduced per clone, neutral mutations covered approximately 10% of the protein sequence. Accordingly, this study opens new avenues for UPO design by bringing together neutral genetic drift and DNA recombination in vivoIMPORTANCE Fungal peroxygenases resemble the peroxide shunt pathway of cytochrome P450 monoxygenases, performing selective oxyfunctionalizations of unactivated C-H bonds in a broad range of organic compounds. In this study, we combined neutral genetic drift and in vivo DNA shuffling to generate highly functional peroxygenase mutant libraries. The panel of neutrally evolved peroxygenases showed different activity profiles for peroxygenative substrates and improved stability with respect to temperature and the presence of organic cosolvents, making the enzymes valuable blueprints for emerging evolution campaigns. This association of DNA recombination and neutral drift is paving the way for future work in peroxygenase engineering and, from a more general perspective, to any other enzyme system heterologously expressed in S. cerevisiae.

Constraint on ion–neutral drift velocity in the Class 0 protostar B335 from ALMA observations

Aims. Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal magnetohydrodynamics (MHD) effects proposed to enable the formation of large-scale Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare here gas kinematics traced by ionized and neutral molecular lines and discuss the implication on ambipolar diffusion. Methods. We analyzed the data of the H13CO+ (3–2) and C18O (2–1) emission in the Class 0 protostar B335 obtained with our ALMA observations. We constructed kinematical models to fit the velocity structures observed in the H13CO+ and C18O emission and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. Results. A central compact (~1′′–2′′) component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. With our kinematical models, the infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85 ± 0.2 km s−1 at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km s−1 at a radius of 100 au in B335. Conclusions. The Hall parameter for H13CO+ is estimated to be ≫1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between the ionized and neutral gas could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. This result is consistent with the expectation from the MHD simulations with a typical ambipolar diffusivity and those without ambipolar diffusion. On the other hand, the high ambipolar drift velocity of 0.5–1.0 km s−1 on a 100 au scale predicted in the MHD simulations with an enhanced ambipolar diffusivity by removing small dust grains, where the minimum grain size is 0.1 μm, is not detected in our observations. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0. ′′1) are needed to examine this possibility and ambipolar diffusion on a smaller scale.

Evolutionary dynamics of paroxysmal nocturnal hemoglobinuria.

Paroxysmal nocturnal hemoglobinuria (PNH) is an acquired clonal blood disorder characterized by hemolysis and a high risk of thrombosis, that is due to a deficiency in several cell surface proteins that prevent complement activation. Its origin has been traced to a somatic mutation in the PIG-A gene within hematopoietic stem cells (HSC). However, to date the question of how this mutant clone expands in size to contribute significantly to hematopoiesis remains under debate. One hypothesis posits the existence of a selective advantage of PIG-A mutated cells due to an immune mediated attack on normal HSC, but the evidence supporting this hypothesis is inconclusive. An alternative (and simpler) explanation attributes clonal expansion to neutral drift, in which case selection neither favours nor inhibits expansion of PIG-A mutated HSC. Here we examine the implications of the neutral drift model by numerically evolving a Markov chain for the probabilities of all possible outcomes, and investigate the possible occurrence and evolution, within this framework, of multiple independently arising clones within the HSC pool. Predictions of the model agree well with the known incidence of the disease and average age at diagnosis. Notwithstanding the slight difference in clonal expansion rates between our results and those reported in the literature, our model results lead to a relative stability of clone size when averaging multiple cases, in accord with what has been observed in human trials. The probability of a patient harbouring a second clone in the HSC pool was found to be extremely low (). Thus our results suggest that in clinical cases of PNH where two independent clones of mutant cells are observed, only one of those is likely to have originated in the HSC pool.

INTEGRATING MACROECOLOGY AND QUANTITATIVE GENETICS: EVOLUTION OF BODY SIZE AND BRAIN SIZE UNDER ISLAND RULE

Island rule proposes that a negative correlation between ancestral body size in continents and the descendent body size in islands exists, and this pattern have been widely studied in a macroecological and comparative perspectives. However, there are doubts about what mechanisms underlie body size evolution in islands. Here we review methodological and theoretical framework on evolutionary quantitative genetics, showing their application on body and brain size evolution in islands, using Hippopotamus dwarfism as example. In our analyses we started by generating 10,000 combinations of model parameters (generation time, effective size and heritability) and tested by Mutation-Drift Equilibrium model if body size dwarfism is a consequence of neutral drift or directional selection. We found that 99.9 % of simulations rejected neutral model. Then, we estimated the strength of directional selection necessary to differentiate the island species and found that a relatively low proportion of population (0.01 - 0.2%) should be selectively killed to decrease body size. Our results also showed that one unit decrease in body size would increase, on average, fitness by 4% in each generation, so directional selection is a plausible explanation to island rule. Finally, we also simulated the evolution of brain size of dwarfed Hippopotamus as a consequence of body size evolution alone. Our estimates of Expected Brain Size (EBS = 484 A 64 cc) were larger than the observed brain size (equal to 380 cc), which suggests the need to estimate directional selection acting on brain size independently of body size evolution. This supports the overall idea that brain size reduction is advantageous in island environments under a scenario reduction in resources, due to the high energetic budget of brain. Our analyses using evolutionary quantitative genetic support that Island Rule as a parsimonious adaptive explanation for the reduction in brain and body sizes and illustrates how to couple evolutionary analyses at population level to better understand macroecological patterns.

Mature enteroendocrine cells contribute to basal and pathological stem cell dynamics in the small intestine.

Lgr5-expressing intestinal stem cells (ISCs) maintain continuous and rapid generation of the intestinal epithelium. Here, we present evidence that dedifferentiation of committed enteroendocrine cells (EECs) contributes to maintenance of the epithelium under both basal conditions and in response to injury. Lineage-tracing studies identified a subset of EECs that reside at +4 position for more than 2 wk, most of which were BrdU-label-retaining cells. Under basal conditions, cells derived from these EECs grow from the bottom of the crypt to generate intestinal epithelium according to neutral drift kinetics that is consistent with dedifferentiation of mature EECs to ISCs. The lineage tracing of EECs demonstrated reserve stem cell properties in response to radiation-induced injury with the generation of reparative EEC-derived epithelial patches. Finally, the enterochromaffin (EC) cell was the predominant EEC type participating in these stem cell dynamics. These results provide novel insights into the +4 reserve ISC hypothesis, stem cell dynamics of the intestinal epithelium, and in the development of EC-derived small intestinal tumors. NEW & NOTEWORTHY The current manuscript demonstrating that a subset of mature enteroendocrine cells (EECs), predominantly enterochromaffin cells, dedifferentiates to fully functional intestinal stem cells (ISCs) is novel, timely, and important. These cells dedifferentiate to ISCs not only in response to injury but also under basal homeostatic conditions. These novel findings provide a mechanism in which a specified cell can dedifferentiate and contribute to normal tissue plasticity as well as the development of EEC-derived intestinal tumors under pathologic conditions.

Nucleoside analogs to manage sequence divergence in nucleic acid amplification and SNP detection.

Described here are the synthesis, enzymology and some applications of a purine nucleoside analog (H) designed to have two tautomeric forms, one complementary to thymidine (T), the other complementary to cytidine (C). The performance of H is compared by various metrics to performances of other ‘biversal’ analogs that similarly rely on tautomerism to complement both pyrimidines. These include (i) the thermodynamic stability of duplexes that pair these biversals with various standard nucleotides, (ii) the ability of the biversals to support polymerase chain reaction (PCR), (iii) the ability of primers containing biversals to equally amplify targets having polymorphisms in the primer binding site, and (iv) the ability of ligation-based assays to exploit the biversals to detect medically relevant single nucleotide polymorphisms (SNPs) in sequences flanked by medically irrelevant polymorphisms. One advantage of H over the widely used inosine ‘universal base’ and ‘mixed sequence’ probes is seen in ligation-based assays to detect SNPs. The need to detect medically relevant SNPs within ambiguous sequences is especially important when probing RNA viruses, which rapidly mutate to create drug resistance, but also suffer neutral drift, the second obstructing simple methods to detect the first. Thus, H is being developed to detect variants of viruses that are rapidly mutating.

Computer simulation of neutral drift among limbal epithelial stem cells of mosaic mice

The use of mice that are mosaic for reporter gene expression underlies many lineage-tracing studies in stem cell biology. For example, using mosaic LacZ reporter mice, it was shown that limbal epithelial stem cells (LESCs) around the periphery of the cornea maintain radial sectors of the corneal epithelium and that radial stripe numbers declined with age. Originally, the corneal results were interpreted as progressive, age-related loss or irreversible inactivation of some LESC clones. In this study we used computer simulations to show that these results could also be explained by stochastic replacement of LESCs by neighbouring LESCs, leading to neutral drift of LESC populations. This was shown to reduce the number of coherent clones of LESCs and hence would coarsen the mosaic pattern in the corneal epithelium without reducing the absolute number of LESCs. Simulations also showed that corrected stripe numbers declined more slowly when LESCs were grouped non-randomly and that mosaicism was rarely lost unless simulated LESC numbers were unrealistically low. Possible reasons why age-related changes differ between mosaic corneal epithelia and other systems, such as adrenal cortices and intestinal crypts, are discussed.

Mo1063 - Neutral Drift Dynamics of Intestinal Stem Cells: Simulation Based on the Branching Process Model with a Symmetric/Asymmetric Cell Division Mode

Mo1063 - Neutral Drift Dynamics of Intestinal Stem Cells: Simulation Based on the Branching Process Model with a Symmetric/Asymmetric Cell Division Mode

Neutral Theory, Transposable Elements, and Eukaryotic Genome Evolution

Among the multitude of papers published yearly in scientific journals, precious few publications may be worth looking back in half a century to appreciate the significance of the discoveries that would later become common knowledge and get a chance to shape a field or several adjacent fields. Here, Kimura's fundamental concept of neutral mutation-random drift, which was published 50 years ago, is re-examined in light of its pervasive influence on comparative genomics and, more specifically, on the contribution of transposable elements to eukaryotic genome evolution.

Hidden long evolutionary memory in a model biochemical network.

We introduce a minimal model for the evolution of functional protein-interaction networks using a sequence-based mutational algorithm, and apply the model to study neutral drift in networks that yield oscillatory dynamics. Starting with a functional core module, random evolutionary drift increases network complexity even in the absence of specific selective pressures. Surprisingly, we uncover a hidden order in sequence space that gives rise to long-term evolutionary memory, implying strong constraints on network evolution due to the topology of accessible sequence space.

Wnt ligands influence tumour initiation by controlling the number of intestinal stem cells

Many epithelial stem cell populations follow a pattern of stochastic stem cell divisions called 'neutral drift'. It is hypothesised that neutral competition between stem cells protects against the acquisition of deleterious mutations. Here we use a Porcupine inhibitor to reduce Wnt secretion at a dose where intestinal homoeostasis is maintained despite a reduction of Lgr5+ stem cells. Functionally, there is a marked acceleration in monoclonal conversion, so that crypts become rapidly derived from a single stem cell. Stem cells located further from the base are lost and the pool of competing stem cells is reduced. We tested whether this loss of stem cell competition would modify tumorigenesis. Reduction of Wnt ligand secretion accelerates fixation of Apc-deficient cells within the crypt leading to accelerated tumorigenesis. Therefore, ligand-based Wnt signalling influences the number of stem cells, fixation speed of Apc mutations and the speed and likelihood of adenoma formation.

Neutral dynamics and cell renewal of colonic crypts in homeostatic regime

The self renewal process in colonic crypts is the object of several studies. We present here a new compartment model with the following characteristics: (a) we distinguish different classes of cells: stem cells, six generations of transit amplifying cells and the differentiated cells; (b) in order to take into account the monoclonal character of crypts in homeostatic regimes we include symmetric divisions of the stem cells. We first consider the dynamic differential equations that describe the evolution of the mean values of the populations, but the small observed value of the total number of cells involved plus the huge dispersion of experimental data found in the literature leads us to study the stochastic discrete process. This analysis allows us to study fluctuations, the neutral drift that leads to monoclonality, and the effects of the fixation of mutant clones.

To the Final Goal: Can We Predict and Suggest Mutations for Protein to Develop Desired Phenotype?

Directed evolution of proteins is a good approach to develop desired phenotypes from existing proteins. Fully experimental protein evolution usually utilizes randomization of a given protein sequence by error-prone PCR or gene shuffling followed by high-throughput selection or timeconsuming screening method. However, these random methods create mutant library full of deleterious mutations. In addition, they need high-throughput screening or selection method to search for positive clones from an enormous size of mutant library. Construction of a mutant library while retaining the original function is important for efficient protein evolution because it greatly reduces time and effort for the identification of positive mutants. Therefore, researchers have tried to reduce the size of mutant library by minimizing the occurrence of deleterious mutants. Such efforts have led to the creation of a concept of ‘small but smart library’. For this goal, neutral drift theory has been applied. Although smart library greatly reduces the library size, it is still the beyond the capacity of low-throughput assay. In parallel, computational analysis of protein structure and efforts to discriminate mutatable residues from all residues of a given protein have been consistently pursued. Accumulated knowledge of protein evolution through random mutation and selection has improved our understanding of functions of amino acids in protein structure. Protein evolution by rational design is being developed based on such understanding. In this review, we describe how the use of semi-rationally designed library rather than completely random one has impacted the overall procedure of directed evolution. We also describe efforts made to evaluate the effect of single mutation. Such efforts will bring lazy boys to the final goal - computational mutation suggestion system.

Spatial and environmental correlates of intraspecific morphological variation in three species of passerine birds from the Purus–Madeira interfluvium, Central Amazonia

Biogeographic studies in Amazonia typically describe biodiversity across interfluvia, rarely within them, where geographic variability in morphological traits might be observed. We tested for intraspecific phenotypic variation in three bird species within the Purus–Madeira interfluvium (Central Amazon) and whether phenotypes were correlated with environmental heterogeneity or geographic distance among sites. We compared coloration indexes derived from reflectance spectra and morphometrics of up to five adult individuals of each sex among 11 sites within the interfluvium and contrasted them with proxies for geographic distance and environmental variation (tree basal area and bird community). Environmental heterogeneity was minimally spatially autocorrelated, and there were no obvious geographical barriers to dispersal in the study region. The null hypothesis was that we would see either no phenotypic variation or random variation that was not explained by the tested variables. Half of the cases analyzed showed intraspecific morphological variation. Coloration varied more frequently than morphometrics, and color was better explained by environmental heterogeneity, particularly in males, whereas brightness also varied with geographic distance. Geographic distance explained the only case of variation in morphometrics. Our results indicate that coloration, particularly plumage brightness, is more labile than morphometric traits and that plumage color might be under stronger effects of local adaptation than brightness, which also seems to be under effects of neutral drift and gene flow among populations. Higher frequencies of association between male coloration and the environment suggest a role of non-arbitrary mechanisms of sexual selection on the expression of male phenotypes, whereas arbitrary intersexual selection might explain the randomly distributed variation that is not explained by environmental heterogeneity or geographic distance. We revealed intraspecific phenotypic variation in a spatial extent usually not considered in biogeographic studies in the Amazon and demonstrate that both local adaptation and neutral drift are important to explain intraspecific trait diversification at this geographical scale.

Neutral Theory and Scale-Free Neural Dynamics

Avalanches of electrochemical activity in brain networks have been empirically reported to obey scale-invariant behavior --characterized by power-law distributions up to some upper cut-off-- both in vitro and in vivo. Elucidating whether such scaling laws stem from the underlying neural dynamics operating at the edge of a phase transition is a fascinating possibility, as systems poised at criticality have been argued to exhibit a number of important functional advantages. Here we employ a well-known model for neural dynamics with synaptic plasticity, to elucidate an alternative scenario in which neuronal avalanches can coexist, overlapping in time, but still remaining scale-free. Remarkably their scale-invariance does not stem from underlying criticality nor self-organization at the edge of a continuous phase transition. Instead, it emerges from the fact that perturbations to the system exhibit a neutral drift --guided by demographic fluctuations-- with respect to endogenous spontaneous activity. Such a neutral dynamics --similar to the one in neutral theories of population genetics-- implies marginal propagation of activity, characterized by power-law distributed causal avalanches. Importantly, our results underline the importance of considering causal information --on which neuron triggers the firing of which-- to properly estimate the statistics of avalanches of neural activity. We discuss the implications of these findings both in modeling and to elucidate experimental observations, as well as its possible consequences for actual neural dynamics and information processing in actual neural networks.

Muscle Stem Cells Exhibit Distinct Clonal Dynamics in Response to Tissue Repair and Homeostatic Aging

The clonal complexity of adult stem cell pools is progressively lost during homeostatic turnover in several tissues, suggesting a decrease in the number of stem cells with distinct clonal origins. The functional impact of reduced complexity on stem cell pools, and how different tissue microenvironments may contribute to such a reduction, are poorly understood. Here, we performed clonal multicolor lineage tracing of skeletal muscle stem cells (MuSCs) to address these questions. We found that MuSC clonal complexity is maintained during aging despite heterogenous reductions in proliferative capacity, allowing aged muscle to mount a clonally diverse, albeit diminished, response to injury. In contrast, repeated bouts of tissue repair cause a progressive reduction in MuSC clonal complexity indicative of neutral drift. Consistently, biostatistical modeling suggests that MuSCs undergo symmetric expansions with stochastic fate acquisition during tissue repair. These findings establish distinct principles that underlie stem cell dynamics during homeostatic aging and muscle regeneration.

Genomic Loss and Silencing on the Y Chromosomes of Rumex.

Across many unrelated lineages of plants and animals, Y chromosomes show a recurrent pattern of gene degeneration and loss, but the relative importance of inefficient selection, adaptive gene silencing, and neutral genetic drift in causing degeneration remain poorly understood. Here, we use next-generation genome and transcriptome sequencing to investigate patterns of ongoing Y chromosome degeneration in two annual plant species of Rumex (Polygonaceae) differing in their degree of degeneration and sex chromosome heteromorphism. We find evidence for both gene loss as well as silencing in these young plant sex chromosomes. Our analyses revealed significantly more gene deletion relative to silencing in R. rothschildianus, which has had a larger nonrecombining region for a longer period than R. hastatulus, consistent with this system being at a more advanced stage of degeneration. Intra- and interspecific comparisons of genomic coverage and heterozygosity indicated that loss of expression precedes gene deletion, implying that the final stages of mutation accumulation and gene loss may often occur neutrally. We found no evidence for adaptive silencing of genes that have lost expression. Our results suggest that the initial spread of deleterious regulatory variants and/or epigenetic silencing may be an important driver of early degeneration of Y chromosomes.

Genomic Signature of Adaptive Divergence despite Strong Nonadaptive Forces on Edaphic Islands: A Case Study of Primulina juliae.

Both genetic drift and divergent selection are expected to be strong evolutionary forces driving population differentiation on edaphic habitat islands. However, the relative contribution of genetic drift and divergent selection to population divergence has rarely been tested simultaneously. In this study, restriction-site associated DNA-based population genomic analyses were applied to assess the relative importance of drift and divergent selection on population divergence of Primulina juliae, an edaphic specialist from southern China. All populations were found with low standing genetic variation, small effective population size (NE), and signatures of bottlenecks. Populations with the lowest genetic variation were most genetically differentiated from other populations and the extent of genetic drift increased with geographic distance from other populations. Together with evidence of isolation by distance, these results support neutral drift as a critical evolutionary driver. Nonetheless, redundancy analysis revealed that genomic variation is significantly associated with both edaphic habitats and climatic factors independently of spatial effects. Moreover, more genomic variation was explained by environmental factors than by geographic variables, suggesting that local adaptation might have played an important role in driving population divergence. Finally, outlier tests and environment association analyses identified 31 single-nucleotide polymorphisms as candidates for adaptive divergence. Among these candidates, 26 single-nucleotide polymorphisms occur in/near genes that potentially play a role in adaptation to edaphic specialization. This study has important implications that improve our understanding of the joint roles of genetic drift and adaptation in generating population divergence and diversity of edaphic specialists.

Ecolottery: Simulating and assessing community assembly with environmental filtering and neutral dynamics inR

AbstractWe introduce theRpackageecolotterydedicated to quick and efficient simulation of communities undergoing local neutral dynamics with environmentally filtered immigration from a reference species pool (spatially implicit model). The package includes an Approximate Bayesian Computation (ABC) tool to estimate the parameters of these processes. We present the rationale of the approach and show examples of simulations and ABC analysis.The species in the reference pool differ in their abundances and trait values. Environmental filtering weights the probability of immigration success depending on trait values, while the descendants of established immigrants undergo neutral stochastic drift. The reference pool can be defined in a flexible way as representing, e.g. the composition of a broad biogeographical region, or available dispersers around local communities. The package provides a process‐based alternative to the use of randomization‐based null models.The package proposes a coalescent‐based simulation algorithm that presents significant advantages over alternative algorithms. It does not require simulating community dynamics from an initial state forward in time but does still allow measurement of the influence of environmental filtering. Because of its high calculation speed, this approach allows simulating many communities within a reasonable amount of time.Diverse patterns of taxonomic, functional and phylogenetic compositions can be generated. The package can be used to explore the outcome of ecological and evolutionary processes playing at local and regional scales, and to estimate the parameters of these processes based on observed patterns.