Molecular biology and evolution最新文献

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is an ancient protein critical for CO2-fixation and global biogeochemistry. Form-I RuBisCO complexes uniquely harbor small subunits that form a hexadecameric complex together with their large subunits. The small subunit protein is thought to have significantly contributed to RuBisCO's response to the atmospheric rise of O2 ∼2.5 billion years ago, marking a pivotal point in the enzyme's evolutionary history. Here, we performed a comprehensive evolutionary analysis of extant and ancestral RuBisCO sequences and structures to explore the impact of the small subunit's earliest integration on the molecular dynamics of the overall complex. Our simulations suggest that the small subunit restricted the conformational flexibility of the large subunit early in its history, impacting the evolutionary trajectory of the Form-I RuBisCO complex. Molecular dynamics investigations of CO2 and O2 gas distribution around predicted ancient RuBisCO complexes suggest that a proposed "CO2-reservoir" role for the small subunit is not conserved throughout the enzyme's evolutionary history. The evolutionary and biophysical response of RuBisCO to changing atmospheric conditions on ancient Earth showcase multi-level and trackable responses of enzymes to environmental shifts over long timescales.

Plant cells have two major organelles with their own genomes: chloroplasts and mitochondria. While chloroplast genomes tend to be structurally conserved, the mitochondrial genomes of plants, which are much larger than those of animals, are characterized by complex structural variation. We introduce TIPPo, a user-friendly, reference-free assembly tool that uses PacBio high-fidelity long-read data and that does not rely on genomes from related species or nuclear genome information for the assembly of organellar genomes. TIPPo employs a deep learning model for initial read classification and leverages k-mer counting for further refinement, significantly reducing the impact of nuclear insertions of organellar DNA on the assembly process. We used TIPPo to completely assemble a set of 54 complete chloroplast genomes. No other tool was able to completely assemble this set. TIPPo is comparable with PMAT in assembling mitochondrial genomes from most species but does achieve even higher completeness for several species. We also used the assembled organelle genomes to identify instances of nuclear plastid DNA (NUPTs) and nuclear mitochondrial DNA (NUMTs) insertions. The cumulative length of NUPTs/NUMTs positively correlates with the size of the nuclear genome, suggesting that insertions occur stochastically. NUPTs/NUMTs show predominantly C:G to T:A changes, with the mutated cytosines typically found in CG and CHG contexts, suggesting that degradation of NUPT and NUMT sequences is driven by the known elevated mutation rate of methylated cytosines. Small interfering RNA loci are enriched in NUPTs and NUMTs, consistent with the RdDM pathway mediating DNA methylation in these sequences.

Regeneration, the ability to restore body parts after injury, is widespread in metazoans; however, the underlying molecular and cellular mechanisms involved in this process remain largely unknown, and its evolutionary history is consequently unresolved. Recently, reactive oxygen species (ROS) have been shown in several metazoan models to be triggers of apoptosis and cell proliferation that drive regenerative success. However, it is not known whether the contribution of ROS to regeneration relies on conserved mechanisms. Here we performed a comparative genomic analysis of ROS metabolism actors across metazoans, and carried out a comparative study of the deployment and roles of ROS during regeneration in two different metazoan models: the annelid Platynereis dumerilii and the cnidarian Nematostella vectensis. We established that the vast majority of metazoans encode a core redox kit allowing for the production and detoxification of ROS, and overall regulation of ROS levels. However, the precise composition of the redox arsenal can vary significantly from species to species, suggesting that evolutionary constraints apply to ROS metabolism functions rather than precise actors. We found that while ROS are necessary for regeneration in both Platynereis and Nematostella, the two species deploy different enzymatic activities controlling ROS dynamics, and display distinct effects of ROS signaling on injury-induced apoptosis and cell proliferation. We conclude that, while ROS are a common feature of metazoan regeneration, their production and contribution to this phenomenon may depend on different molecular mechanisms highlighting the overall plasticity of the machinery.

Determining the impact of mutations on the thermodynamic stability of proteins is essential for a wide range of applications such as rational protein design and genetic variant interpretation. Since protein stability is a major driver of evolution, evolutionary data are often used to guide stability predictions. Many state-of-the-art stability predictors extract evolutionary information from multiple sequence alignments of proteins homologous to a query protein, and leverage it to predict the effects of mutations on protein stability. To evaluate the power and the limitations of such methods, we used the massive amount of stability data recently obtained by deep mutational scanning to study how best to construct multiple sequence alignments and optimally extract evolutionary information from them. We tested different evolutionary models and found that, unexpectedly, independent-site models achieve similar accuracy to more complex epistatic models. A detailed analysis of the latter models suggests that their inference often results in noisy couplings, which do not appear to add predictive power over the independent-site contribution, at least in the context of stability prediction. Interestingly, by combining any of the evolutionary features with a simple structural feature, the relative solvent accessibility of the mutated residue, we achieved similar prediction accuracy to supervised, machine learning-based, protein stability change predictors. Our results provide new insights into the relationship between protein evolution and stability, and show how evolutionary information can be exploited to improve the performance of mutational stability prediction.

For a long time, snakes were presented as a textbook example of a group with gradual differentiation of homologous ZZ/ZW sex chromosomes. However, recent advances revealed that the ZZ/ZW sex chromosomes characterize only caenophidian snakes and certain species of boas and pythons have nonhomologous XX/XY sex chromosomes. We used genome coverage analysis in four non-caenophidian species to identify their sex chromosomes, and we examined the homology of sex chromosomes across phylogenetically informative snake lineages. We identified sex chromosomes for the first time in 13 species of non-caenophidian snakes, providing much deeper insights into the evolutionary history of snake sex chromosomes. The evolution of sex chromosomes in snakes is more complex than previously thought. Snakes may have had ancestral XX/XY sex chromosomes, which are still present in a blind snake and some boas, and there were several transitions to derived XX/XY sex chromosomes with different gene content and two or even three transitions to ZZ/ZW sex chromosomes. However, we discuss more alternative scenarios. In any case, we document that (1) some genomic regions were likely repeatedly co-opted as sex chromosomes in phylogenetically distant lineages, even with opposite types of heterogamety; (2) snake lineages differ greatly in the rate of differentiation of sex chromosomes; (3) snakes likely originally possessed sex chromosomes prone to turnovers. The sex chromosomes became evolutionarily highly stable once their differentiation progressed in the megadiverse caenophidian snakes. Snakes thus provide an ideal system for studying the evolutionary factors that drive unequal rates of differentiation, turnovers and stability of sex chromosomes.

The impact of gene loss on the diversification of taxa and the emergence of evolutionary innovations remains poorly understood. Here, our investigation on the evolution of the Fibroblast Growth Factors (FGFs) in appendicularian tunicates as a case study reveals a scenario of "less, but more" characterized by massive losses of all Fgf gene subfamilies, except for the Fgf9/16/20 and Fgf11/12/13/14, which in turn underwent two bursts of duplications. Through phylogenetic analysis, synteny conservation, and gene and protein structure, we reconstruct the history of appendicularian Fgf genes, highlighting their paracrine and intracellular functions. An exhaustive analysis of developmental Fgf expression in Oikopleura dioica allows us to identify four associated evolutionary patterns characterizing the "less, but more" conceptual framework: conservation of ancestral functions; function shuffling between paralogs linked to gene losses; innovation of new functions after the duplication bursts; and function extinctions linked to gene losses. Our findings allow us to formulate novel hypotheses about the impact of Fgf losses and duplications on the transition from an ancestral ascidian-like biphasic lifestyle to the fully free-living appendicularians. These hypotheses include massive co-options of Fgfs for the development of the oikoblast and the tail fin; recruitment of Fgf11/12/13/14s into the evolution of a new mouth, and their role modulating neuronal excitability; the evolutionary innovation of an anterior tail FGF signaling source upon the loss of retinoic acid signaling; and the potential link between the loss of Fgf7/10/22 and Fgf8/17/18 and the loss of drastic metamorphosis and tail absorption in appendicularians, in contrast to ascidians.

When introduced to multiple distinct ranges, invasive species provide a compelling natural experiment for understanding the repeatability of adaptation. Ambrosia artemisiifolia is an invasive, noxious weed, and chief cause of hay fever. Leveraging over 400 whole-genome sequences spanning the native-range in North America and 2 invasions in Europe and Australia, we inferred demographically distinct invasion histories on each continent. Despite substantial differences in genetic source and effective population size changes during introduction, scans of both local climate adaptation and divergence from the native-range revealed genomic signatures of parallel adaptation between invasions. Disproportionately represented among these parallel signatures are 37 large haploblocks-indicators of structural variation-that cover almost 20% of the genome and exist as standing genetic variation in the native-range. Many of these haploblocks are associated with traits important for adaptation to local climate, like size and the timing of flowering, and have rapidly reformed native-range clines in invaded ranges. Others show extreme frequency divergence between ranges, consistent with a response to divergent selection on different continents. Our results demonstrate the key role of large-effect standing variants in rapid adaptation during range expansion, a pattern that is robust to diverse invasion histories.

Bacterial genomes primarily diversify via gain, loss, and rearrangement of genetic material in their flexible accessory genome. Yet the dynamics of accessory genome evolution are very poorly understood, in contrast to the core genome where diversification is readily described by mutations and homologous recombination. Here, we tackle this problem for the case of very closely related genomes. We comprehensively describe genome evolution within n=222 genomes of E. coli ST131, which likely shared a common ancestor around one hundred years ago. After removing putative recombinant diversity, the total length of the phylogeny is 6000 core genome substitutions. Within this diversity, we find 22 modifications to core genome synteny and estimate around 2000 structural changes within the accessory genome, i.e. one structural change for every 3 core genome substitutions. 63% of loci with structural diversity could be resolved into individual gain and loss events with ten-fold more gains than losses, demonstrating a dominance of gains due to insertion sequences and prophage integration. Our results suggest the majority of synteny changes and insertions in our dataset are likely deleterious and only persist for a short time before being removed by purifying selection.

We introduce the 12th version of the Molecular Evolutionary Genetics Analysis (MEGA) software. This latest version brings many significant improvements by reducing the computational time needed for selecting optimal substitution models and conducting bootstrap tests on phylogenies using maximum likelihood (ML) methods. These improvements are achieved by implementing heuristics that minimize likely unnecessary computations. Analyses of empirical and simulated datasets show substantial time savings by using these heuristics without compromising the accuracy of results. MEGA12 also implements an evolutionary sparse learning approach to identify fragile clades and associated sequences in evolutionary trees inferred through phylogenomic analyses. In addition, this version includes fine-grained parallelization for ML analyses, support for high-resolution monitors, and an enhanced Tree Explorer. MEGA12 can be downloaded from https://www.megasoftware.net.

Marine mammals have evolved unihemispheric slow wave sleep (USWS), a unique state during which one cerebral hemisphere sleeps while the other remains awake, to mitigate the fundamental conflict between sleep and wakefulness. However, the underlying mechanisms remain largely unclear. Here, we use a comparative phylogenetic approach to analyze genes associated with light-dependent circadian mechanisms, aiming to reconstruct the evolution of the circadian rhythm pathway in mammals and to identify adaptively changed components likely to have contributed to the development of USWS. Specifically, among eight genes with shared signals of positive selection in two USWS-specific lineages, seven genes showed direct evidence of affecting sleep and spontaneous movements. Both in vitro and in vivo experiments indicated that functional innovation in cetacean and non-phocid pinniped FBXL21, which was found to undergo positive selection, may be beneficial for decoupling sleep-wake patterns from daily rhythms to sustain continuous swimming. For cetaceans exhibiting only USWS, we identified 73 genes as rapidly evolving and 92 genes containing unique amino acid substitutions. Functional assays showed that a cetacean-specific mutation (F411Y) in NFIL3 led to a decrease in repressor activity and protein stability. Furthermore, convergent amino acid replacements detected in genes related to Ca2+ signaling and CREB phosphorylation suggest their crucial role in USWS adaptation. Overall, this study enhances our understanding of the evolutionary mechanisms underlying USWS and provides a foundation for investigating how circadian rhythm changes contribute to variations in sleep and circadian behavior.