Search

Bacteria, Archaea, and Eukarya all share a common set of metabolic reactions. This implies that the function and topology of central metabolism has been evolving under purifying selection over deep time. Central metabolism may similarly evolve under purifying selection during long-term evolution experiments, although it is unclear how long such experiments would have to run (decades, centuries, millennia) before signs of purifying selection on metabolism appear. I hypothesized that central and superessential metabolic enzymes would show evidence of purifying selection in the long-term evolution experiment with Escherichia coli (LTEE). I also hypothesized that enzymes that specialize on single substrates would show stronger evidence of purifying selection in the LTEE than generalist enzymes that catalyze multiple reactions. I tested these hypotheses by analyzing metagenomic time series covering 62,750 generations of the LTEE. I find mixed support for these hypotheses, because the observed patterns of purifying selection are idiosyncratic and population-specific. To explain this finding, I propose the Jenga hypothesis, named after a children's game in which blocks are removed from a tower until it falls. The Jenga hypothesis postulates that loss-of-function mutations degrade costly, redundant, and non-essential metabolic functions. Replicate populations can therefore follow idiosyncratic trajectories of lost redundancies, despite purifying selection on overall function. I tested the Jenga hypothesis by simulating the evolution of 1,000 minimal genomes under strong purifying selection. As predicted, the minimal genomes converge to different metabolic networks. Strikingly, the core genes common to all 1,000 minimal genomes show consistent signatures of purifying selection in the LTEE. [ABSTRACT FROM AUTHOR]

Background: Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites. Results: Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in "cryptic" splicing regions beyond canonical donor or acceptor splice sites. Conclusion: Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases. [ABSTRACT FROM AUTHOR]

MADS domain transcription factors play roles throughout the whole lifecycle of plants from seeding to flowering and fruit-bearing. However, systematic research into MADS-box genes of the economically important vegetable crop pepper (Capsicum spp.) is still lacking. We identified 174, 207, and 72 MADS-box genes from the genomes of C. annuum, C. baccatum, and C. chinense, respectively. These 453 MADS-box genes were divided into type I (Mα, Mβ, Mγ) and type II (MIKC* and MIKCC) based on their phylogenetic relationships. Collinearity analysis identified 144 paralogous genes and 195 orthologous genes in the three Capsicum species, and 70, 114, and 10 MADS-box genes specific to C. annuum, C. baccatum, and C. chinense, respectively. Comparative genomic analysis highlighted functional differentiation among homologous MADS-box genes during pepper evolution. Tissue expression analysis revealed three main expression patterns: highly expressed in roots, stems, leaves, and flowers (CaMADS93/CbMADS35/CcMADS58); only expressed in roots; and specifically expressed in flowers (CaMADS26/CbMADS31/CcMADS11). Protein interaction network analysis showed that type II CaMADS mainly interacted with proteins related to flowering pathway and flower organ development. This study provides the basis for an in-depth study of the evolutionary features and biological functions of pepper MADS-box genes. [ABSTRACT FROM AUTHOR]

Background: Mycobacteroides abscessus complex (MABC), an emerging pathogen, causes human infections resistant to multiple antibiotics. In this study, the genome data of 1,581 MABC strains were downloaded from NCBI database for phylogenetic relatedness inference, resistance profile identification and the estimation of evolutionary pressure on resistance genes in silico. Results: From genes associated with resistance to 28 antibiotic classes, 395 putative proteins (ARPs) were identified, based on the information in two antibiotic resistance databases (CARD and ARG-ANNOT). The ARPs most frequently identified in MABC were those associated with resistance to multiple antibiotic classes, beta-lactams and aminoglycosides. After excluding ARPs that had undergone recombination, two ARPs were predicted to be under diversifying selection and 202 under purifying selection. This wide occurrence of purifying selection suggested that the diversity of commonly shared ARPs in MABC have been reduced to achieve stability. The unequal distribution of ARPs in members of the MABC could be due to horizontal gene transfer or ARPs pseudogenization events. Most (81.5%) of the ARPs were observed in the accessory genome and 72.2% ARPs were highly homologous to proteins associated with mobile genetic elements such as plasmids, prophages and viruses. On the other hand, with TBLASTN search, only 18 of the ARPs were identified as pseudogenes. Conclusion: Altogether, our results suggested an important role of horizontal gene transfer in shaping the resistome of MABC. [ABSTRACT FROM AUTHOR]

The evolutionary dynamics and phylogenetic utility of mitochondrial genomes (mitogenomes) have been of particular interest to systematists and evolutionary biologists. However, certain mitochondrial features, such as the molecular evolution of the control region in insects, remain poorly explored due to technological constraints. Using a combination of long- and short-read sequencing data, we assembled ten complete mitogenomes from ten Hercules beetles. We found large-sized mitogenomes (from 24 to 28 kb), which are among the largest in insects. The variation in genome size can be attributed to copy-number evolution of tandem repeats in the control region. Furthermore, one type of tandem repeat was found flanking the conserved sequence block in the control region. Importantly, such variation, which made up around 30% of the size of the mitogenome, may only become detectable should long-read sequencing technology be applied. We also found that, although different mitochondrial loci often inferred different phylogenetic histories, none of the mitochondrial loci statistically reject a concatenated mitochondrial phylogeny, supporting the hypothesis that all mitochondrial loci share a single genealogical history. We on the other hand reported statistical support for mito-nuclear phylogenetic discordance in 50% of mitochondrial loci. We argue that long-read DNA sequencing should become a standard application in the rapidly growing field of mitogenome sequencing. Furthermore, mitochondrial gene trees may differ even though they share a common genealogical history, and ND loci could be better candidates for phylogenetics than the commonly used COX1. [ABSTRACT FROM AUTHOR]

Sugar plays an important role in regulating the flowering of plants. However, studies of genes related to flowering regulation by the sugar pathway of Brassicaceae plants are scarce. In this study, we performed a comprehensive comparative genomics analysis of the flowering genes in the sugar pathway from seven members of the Brassicaceae, including: Arabidopsis thaliana, Arabidopsis lyrata, Astelia pumila, Camelina sativa, Brassica napus, Brassica oleracea, and Brassica rapa. We identified 105 flowering genes in the sugar pathway of these plants, and they were categorized into nine groups. Protein domain analysis demonstrated that the IDD8 showed striking structural variations in different Brassicaceae species. Selection pressure analysis revealed that sugar pathway genes related to flowering were subjected to strong purifying selection. Collinearity analysis showed that the identified flowering genes expanded to varying degrees, but SUS4 was absent from the genomes of Astelia pumila, Camelina sativa, Brassica napus, Brassica oleracea, and Brassica rapa. Tissue-specific expression of ApADG indicated functional differentiation. To sum up, genome-wide identification revealed the expansion, contraction, and diversity of flowering genes in the sugar pathway during Brassicaceae evolution. This study lays a foundation for further study on the evolutionary characteristics and potential biological functions of flowering genes in the sugar pathway of Brassicaceae. [ABSTRACT FROM AUTHOR]

The conventionally clear distinction between exons and introns in eukaryotic genes is actually blurred. To illustrate this point, consider sequences that are retained in mature mRNAs about 50% of the time: how should they be classified? Moreover, although it is clear that RNA splicing influences gene expression levels and is an integral part of interdependent cellular networks, introns continue to be regarded as accidental insertions; exogenous sequences whose evolutionary origin is independent of mRNA-associated processes and somewhat still elusive. Here, we present evidence that aids to resolve this disconnect between conventional views about introns and current knowledge about the role of RNA splicing in the eukaryotic cell. We first show that coding sequences flanked by cryptic splice sites are negatively selected on a genome-wide scale in Paramecium. Then, we exploit selection intensity to infer splicing-related evolutionary dynamics. Our analyses suggest that intron gain begins as a splicing error, involves a transient phase of alternative splicing, and is preferentially completed at the 5' end of genes, which through intron gain can become highly expressed. We conclude that relaxed selective constraints may promote biological complexity in Paramecium and that the relationship between exons and introns is fluid on an evolutionary scale. [ABSTRACT FROM AUTHOR]

Abstract Detecting selective sweeps driven by strong positive selection and localizing the targets of selection in the genome play a major role in modern population genetics and genomics. Most of these analyses are based on the classical model of genetic hitchhiking proposed by Maynard Smith and Haigh (1974, Genetical Research, 23, 23). Here, we consider extensions of the classical two‐locus model. Introducing mutation at the strongly selected site, we analyze the conditions under which soft sweeps may arise. We identify a new parameter (the ratio of the beneficial mutation rate to the selection coefficient) that characterizes the occurrence of multiple‐origin soft sweeps. Furthermore, we quantify the hitchhiking effect when the polymorphism at the linked locus is not neutral but maintained in a mutation‐selection balance. In this case, we find a smaller relative reduction of heterozygosity at the linked site than for a neutral polymorphism. In our analysis, we use a semi‐deterministic approach; i.e., we analyze the frequency process of the beneficial allele in an infinitely large population when its frequency is above a certain threshold; however, for very small frequencies in the initial phase after the onset of selection we rely on diffusion theory.

ABSTRACT Long-standing cyanobacterial harmful algal blooms (CyanoHABs) are known to result from synergistic interaction between elevated nutrients and superior ecophysiology of cyanobacteria. However, it remains to be determined whether CyanoHABs are a result of positive selection by eutrophic waters. To address this, we conducted molecular evolutionary analyses on the genomes of 9 bloom-forming cyanobacteria, combined with pangenomics and metatranscriptomics. The results showed no positive selection by water eutrophication. Instead, all homologous genes in the species are under strong purifying selection based on the ratio of divergence at nonsynonymous and synonymous sites (dN/dS) and phylogeny. The dN/dS

Analyses of genetic variation in many taxa have established that neutral genetic diversity is shaped by natural selection at linked sites. Whether the mode of selection is primarily the fixation of strongly beneficial alleles (selective sweeps) or purifying selection on deleterious mutations (background selection) remains unknown, however. We address this question in humans by fitting a model of the joint effects of selective sweeps and background selection to autosomal polymorphism data from the 1000 Genomes Project. After controlling for variation in mutation rates along the genome, a model of background selection alone explains ~60% of the variance in diversity levels at the megabase scale. Adding the effects of selective sweeps driven by adaptive substitutions to the model does not improve the fit, and when both modes of selection are considered jointly, selective sweeps are estimated to have had little or no effect on linked neutral diversity. The regions under purifying selection are best predicted by phylogenetic conservation, with ~80% of the deleterious mutations affecting neutral diversity occurring in non-exonic regions. Thus, background selection is the dominant mode of linked selection in humans, with marked effects on diversity levels throughout autosomes.

Mitochondria harbor an independent genome, called mitochondrial DNA (mtDNA), which contains essential metabolic genes. Although mtDNA mutations occur at high frequency, they are inherited infrequently, indicating that germline mechanisms limit their accumulation. To determine how germline mtDNA is regulated, we examined the control of mtDNA quantity and quality in C. elegans primordial germ cells (PGCs). We show that PGCs combine strategies to generate a low point in mtDNA number by segregating mitochondria into lobe-like protrusions that are cannibalized by adjacent cells, and by concurrently eliminating mitochondria through autophagy, reducing overall mtDNA content twofold. As PGCs exit quiescence and divide, mtDNAs replicate to maintain a set point of ~200 mtDNAs per germline stem cell. Whereas cannibalism and autophagy eliminate mtDNAs stochastically, we show that the kinase PTEN-induced kinase 1 (PINK1), operating independently of Parkin and autophagy, preferentially reduces the fraction of mutant mtDNAs. Thus, PGCs employ parallel mechanisms to control both the quantity and quality of the founding population of germline mtDNAs.

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied proand eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection. [ABSTRACT FROM AUTHOR]

Recent plant genomic studies provide fine-grained details on the evolutionary consequences of adaptive introgression during crop domestication. Modern genomic approaches and analytical methods now make it possible to better separate the introgression signal from the demographic signal thus providing a more comprehensive and complex picture of the role of introgression in local adaptation. Adaptive introgression has been fundamental for crop expansion and has involved complex patterns of gene flow. In addition to providing new and more favorable alleles of large effect, introgression during the early stages of domestication also increased allelic diversity at adaptive loci. Previous studies have largely underestimated the effect of such increased diversity following introgression. Recent genomic studies in wheat, potato, maize, grapevine, and ryegrass show that introgression of multiple genes, of as yet unknown effect, increased the effectiveness of purifying selection, and promoted disruptive or fluctuating selection in early cultivars and landraces. Historical selection processes associated with introgression from crop wild relatives provide an instructive analog for adaptation to current climate change and offer new avenues for crop breeding research that are expected to be instrumental for strengthening food security in the coming years. [ABSTRACT FROM AUTHOR]

The human Respiratory Syncytial Virus (hRSV) is the main causative agent of acute respiratory infections (ARI), such as pneumonia and bronchiolitis. One of the factors that lead to success in viral replication is the interaction of the M2-2 protein with the ribosomal complex. This interaction is responsible for the phase change of viral activity, acting as an inhibitor or inducer of viral replication, according to the concentration of mRNA. Based on the importance of M2-2 gene and protein have to viral physiology, we performed here evaluations of genetic diversity, phylogenetic reconstructions, phylodynamics, and selection test. Our results suggested an alternative way of classifying this virus in clades A and B, based on a new phylogenetic marker, the M2-2 gene. Therefore, our study is the first one to investigate the dynamics of the evolutionary diversification process of hRSV from the perspective of the M2-2 viral gene. In our study was also identified that the M2-2 gene is under the effect of purifying selection originated by population genetic bottlenecks. Therefore, the M2-2 gene demonstrated an interesting potential to be applied in evolutionary studies involving hRSV, recovering phylogenetic signals and traits of natural selection under the evolution of this virus.

The evolution and diversity of the supergroup Amoebozoa is complex and poorly understood. The supergroup encompasses predominantly amoeboid lineages characterized by extreme diversity in phenotype, behavior and genetics. The study of natural selection, a driving force of diversification, within and among species of Amoebozoa will play a crucial role in understanding the evolution of the supergroup. In this study, we searched for traces of natural selection based on a set of highly conserved protein-coding genes in a phylogenetic framework from a broad sampling of amoebozoans. Using these genes, we estimated substitution rates and inferred patterns of selective pressure in lineages and sites with various models. We also examined the effect of selective pressure on codon usage bias and potential correlations with observed biological traits and habitat. Results showed large heterogeneity of selection across lineages of Amoebozoa, indicating potential species-specific optimization of adaptation to their diverse ecological environment. Overall, lineages in Tubulinea had undergone stronger purifying selection with higher average substitution rates compared to Discosea and Evosea. Evidence of adaptive evolution was observed in some representative lineages and in a gene (Rpl7a) within Evosea, suggesting potential innovation and beneficial mutations in these lineages. Our results revealed that members of the fast-evolving lineages, Entamoeba and Cutosea, all underwent strong purifying selection but had distinct patterns of codon usage bias. For the first time, this study revealed an overall pattern of natural selection across the phylogeny of Amoebozoa and provided significant implications on their distinctive evolutionary processes.

ABSTRACT Candidatus (Ca.) Liberibacter taxa are economically important bacterial plant pathogens that are not culturable; however, genome-enabled insights can help us develop a deeper understanding of their host-microbe interactions and evolution. The draft genome of a recently identified Liberibacter taxa, Ca. Liberibacter capsica, was curated and annotated here with a total draft genome size of 1.1 MB with 1,036 proteins, which is comparable to other Liberibacter species with complete genomes. A total of 459 orthologous clusters were identified among Ca. L. capsica, Ca. L. asiaticus, Ca. L. psyllaurous, Ca. L. americanus, Ca. L. africanus, and L. crescens, and these genes within these clusters consisted of housekeeping and environmental response functions. We estimated the rates of molecular evolution for each of the 443 one-to-one ortholog clusters and found that all Ca. L. capsica orthologous pairs were under purifying selection when the synonymous substitutions per synonymous site (dS) were not saturated. These results suggest that these genes are largely maintaining their conserved functions. We also identified the most divergent single-copy orthologous proteins in Ca. L. capsica by analyzing the ortholog pairs that represented the highest nonsynonymous substitutions per nonsynonymous site (dN) values for each pairwise comparison. From these analyses, we found that 21 proteins which are known to be involved in pathogenesis and host-microbe interactions, including the Tad pilus complex, were consistently divergent between Ca. L. capsica and the majority of other Liberibacter species. These results further our understanding of the evolutionary genetics of Ca. L. capsica and, more broadly, the evolution of Liberibacter. IMPORTANCE “Candidatus” (Ca.) Liberibacter taxa are economically important plant pathogens vectored by insects; however, these host-dependent bacterial taxa are extremely difficult to study because they are unculturable. Recently, we identified a new Ca. Liberibacter lineage (Ca. Liberibacter capsica) from a rare insect metagenomic sample. In this current study, we report that the draft genome of Ca. Liberibacter capsica is similar in genome size and protein content compared to the other Ca. Liberibacter taxa. We provide evidence that many of their shared genes, which encode housekeeping and environmental response functions, are evolving under purifying selection, suggesting that these genes are maintaining similar functions. Our study also identifies 21 proteins that are rapidly evolving amino acid changes in Ca. Liberibacter capsica compared to the majority of other Liberibacter taxa. Many of these proteins represent key genes involved in Liberibacter-host interactions and pathogenesis and are valuable candidate genes for future studies.

ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to represent a global health emergency as a highly transmissible, airborne virus. An important coronaviral drug target for treatment of COVID-19 is the conserved main protease (Mpro). Nirmatrelvir is a potent Mpro inhibitor and the antiviral component of Paxlovid. The significant viral sequencing effort during the ongoing COVID-19 pandemic represented a unique opportunity to assess potential nirmatrelvir escape mutations from emerging variants of SARS-CoV-2. To establish the baseline mutational landscape of Mpro prior to the introduction of Mpro inhibitors, Mpro sequences and its cleavage junction regions were retrieved from ~4,892,000 high-quality SARS-CoV-2 genomes in the open-access Global Initiative on Sharing Avian Influenza Data (GISAID) database. Any mutations identified from comparison to the reference sequence (Wuhan-Hu-1) were catalogued and analyzed. Mutations at sites key to nirmatrelvir binding and protease functionality (e.g., dimerization sites) were still rare. Structural comparison of Mpro also showed conservation of key nirmatrelvir contact residues across the extended Coronaviridae family (α-, β-, and γ-coronaviruses). Additionally, we showed that over time, the SARS-CoV-2 Mpro enzyme remained under purifying selection and was highly conserved relative to the spike protein. Now, with the emergency use authorization (EUA) of Paxlovid and its expected widespread use across the globe, it is essential to continue large-scale genomic surveillance of SARS-CoV-2 Mpro evolution. This study establishes a robust analysis framework for monitoring emergent mutations in millions of virus isolates, with the goal of identifying potential resistance to present and future SARS-CoV-2 antivirals. IMPORTANCE The recent authorization of oral severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antivirals, such as Paxlovid, has ushered in a new era of the COVID-19 pandemic. The emergence of new variants, as well as the selective pressure imposed by antiviral drugs themselves, raises concern for potential escape mutations in key drug binding motifs. To determine the potential emergence of antiviral resistance in globally circulating isolates and its implications for the clinical response to the COVID-19 pandemic, sequencing of SARS-CoV-2 viral isolates before, during, and after the introduction of new antiviral treatments is critical. The infrastructure built herein for active genetic surveillance of Mpro evolution and emergent mutations will play an important role in assessing potential antiviral resistance as the pandemic progresses and Mpro inhibitors are introduced. We anticipate our framework to be the starting point in a larger effort for global monitoring of the SARS-CoV-2 Mpro mutational landscape.

Abstract Background Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection. Results Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species. Conclusions Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

[Display omitted] • There is an increase in the transition transversion bias, nucleotide diversity and purifying selection in the B lineage on a temporal scale. • The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain leading to a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. • The binding site 501 on the spike protein in B lineage coevolved with other amino acid sites. The descendants of the B lineage are the most predominant variants among the SARS-CoV-2 virus due to the incorporation of new mutations augmenting the infectivity of the virus. There is a substantial increase in the transition transversion bias, nucleotide diversity and purifying selection on the spike protein in the descendants of the B lineage of the SARS-CoV-2 virus on a temporal scale. A strong bias for C-to-U substitutions is found in the genes encoding spike protein in this lineage. The positive selection has operated on the spike gene of B lineages and its sub-lineages. The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain ultimately reflected in a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. The intensity of purifying selection on the multiple sites of the spike gene has increased substantially in the sub-lineages of B.1 in a timescale. The binding site 501 on the spike protein in B lineage is found to coevolve with other amino acid sites. This study sheds light on the evolutionary trajectory of the B lineage into highly infectious descendants in the recent past under the influence of positive and purifying selection exerted by natural immunity and vaccination of the host. [ABSTRACT FROM AUTHOR]

Cyanobacterial harmful algal blooms (CyanoHABs) are longstanding aquatic hazards worldwide, of which the mechanism is not yet fully understood, i.e., the process in which cyanobacteria establish dominance over coexisting algae in the same eutrophic waters. The dominance of CyanoHABs represents a deviation from their low abundance under conventional evolution in the oligotrophic state, which has been the case since the origin of cyanobacteria on early Earth. To piece together a comprehensive mechanism of CyanoHABs, we revisit the origin and adaptive radiation of cyanobacteria in oligotrophic Earth, demonstrating ubiquitous adaptive radiation enabled by corresponding biological functions under various oligotrophic conditions. Next, we summarize the biological functions (ecophysiology) which drive CyanoHABs and ecological evidence to synthesize a working mechanism at the population level (the special mechanism) for CyanoHABs: CyanoHABs are the consequence of the synergistic interaction between superior cyanobacterial ecophysiology and elevated nutrients. Interestingly, these biological functions are not a result of positive selection by water eutrophication, but an adaptation to a longstanding oligotrophic state as all the genes in cyanobacteria are under strong negative selection. Last, to address the relative dominance of cyanobacteria over coexisting algae, we postulate a “general” mechanism of CyanoHABs at the community level from an energy and matter perspective: cyanobacteria are simpler life forms and thus have lower per capita nutrient demand for growth than coexisting eukaryotic algae. We prove this by comparing cyanobacteria and eukaryotic algae in cell size and structure, genome size, size of genome-scale metabolic networks, cell content, and finally the golden standard—field studies with nutrient supplementation in the same waters. To sum up, the comprehensive mechanism of CyanoHABs comprises a necessary condition, which is the general mechanism, and a sufficient condition, which is the special mechanism. One prominent prediction based on this tentative comprehensive mechanism is that eukaryotic algal blooms will coexist with or replace CyanoHABs if eutrophication continues and goes over the threshold nutrient levels for eukaryotic algae. This two-fold comprehensive mechanism awaits further theoretic and experimental testing and provides an important guide to control blooms of all algal species.

Evolution is change over time. Although neutral changes promoted by drift effects are most reliable for phylogenetic reconstructions, selection-relevant changes are of only limited use to reconstruct phylogenies. On the other hand, comparative analyses of neutral and selected changes of protein-coding DNA sequences (CDS) retrospectively tell us about episodic constrained, relaxed, and adaptive incidences. The ratio of sites with nonsynonymous (amino acid altering) versus synonymous (not altering) mutations directly measures selection pressure and can be analysed by using the Phylogenetic Analysis by Maximum Likelihood (PAML) software package. We developed a CDS extractor for compiling protein-coding sequences (CDS-extractor) and parallel PAML (paPAML) to simplify, amplify, and accelerate selection analyses via parallel processing, including detection of negatively selected sites. paPAML compiles results of site, branch-site, and branch models and detects site-specific negative selection with the output of a codon list labelling significance values. The tool simplifies selection analyses for casual and inexperienced users and accelerates computing speeds up to the number of allocated computer threads. We then applied paPAML to examine the evolutionary impact on a new GINS Complex Subunit 3 exon, and neutrophil-associated as well as lysin and apolipoprotein genes. Compared with codeml (PAML version 4.9j) and HyPhy (HyPhy FEL version 2.5.26), all paPAML test runs performed with 10 computing threads led to identical selection pressure results, whereas the total selection analysis via paPAML, including all model comparisons, was about 3 to 5 times faster than the longest running codeml model and about 7 to 15 times faster than the entire processing time of these codeml runs. [ABSTRACT FROM AUTHOR]

Baird's tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as 'endangered' on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG's). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 proteincoding genes (PCG's), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (OL), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (OH). A majority of the PCG's and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG's show values<1, suggesting that all these PCGs experience strong purifying selection. A maximum-likelihood phylogenetic analysis supports the monophyly of the genus Tapirus and the order Perissodactyla. The complete annotation and analysis of the mitochondrial genome of T. bairdii will contribute to a better understanding of the population genomic diversity and structure of this species, and it will assist in the conservation and protection of its dwindling populations. [ABSTRACT FROM AUTHOR]

Feather colours are used by avian species for defense, adaptation and signaling. Melanocortin-1 receptor (MC1R) gene is one of the genes responsible for feather colour. This study identified selection signatures in MC1R gene of Nigerian indigenous turkeys (NIT) using British United turkeys (BUT) as control breed to investigate the evolutionary processes that have shaped NIT with various feather colours. Complete MC1R gene of 146 NIT (76 males and 70 females) and 32 BUT (18 males and 14 females) were sequenced. Transition/transversion and codon usage biases were predicted using MEGA v6 software. The selective force acting on the gene was predicted using HyPhy software. The FST values were estimated using Arlequin v3.5. The highest transition/transversion bias was predicted for white BUT (1.00) while the lowest was predicted for black NIT (0.50). Negative dN-dS values, indicative of purifying selection, were observed in MC1R gene of all the turkeys. The highest pairwise FST was observed between the MC1R gene of white BUT and black NIT while the least was observed between lavender NIT and white NIT. No recombination event was observed in black NIT and white BUT. The relative synonymous codon usage was the same among different colours for some codons. Presence of purifying selection in MC1R gene of all the turkeys with different feather colours confirms that the gene plays role in many biological processes such as feather colouration, behaviour, pain perception, immunity, growth and adaptation. The results also suggested that the genetic mechanisms generating different feather colours in turkeys are conserved. [ABSTRACT FROM AUTHOR]

"Junk DNA" is a popular yet controversial concept that states that organisms carry in their genomes DNA that has no positive impact on their fitness. Nonetheless, biochemical functions have been identified for an increasing fraction of DNA elements traditionally seen as "Junk DNA". These findings have been interpreted as fundamentally undermining the "Junk DNA" concept. Here, we reinforce previous arguments that this interpretation relies on an inadequate concept of biological function that does not consider the selected effect of a given genomic structure, which is central to the "Junk DNA" concept. Next, we suggest that another (though ignored) confounding factor is that the discussion about biological functions includes two different dimensions: a horizontal, ecological dimension that reflects how a given genomic element affects fitness in a specific time, and a vertical, temporal dimension that reflects how a given genomic element persisted along time. We suggest that "Junk DNA" should be used exclusively relative to the horizontal dimension, while for the vertical dimension, we propose a new term, "Spam DNA", that reflects the fact that a given genomic element may persist in the genome even if not selected for on their origin. Importantly, these concepts are complementary. An element can be both "Spam DNA" and "Junk DNA", and "Spam DNA" can also be recruited to perform evolved biological functions, as illustrated in processes of exaptation or constructive neutral evolution. [ABSTRACT FROM AUTHOR]

Transposable elements are a major component of most eukaryotic genomes. Here, we present a new approach which allows us to study patterns of natural selection in the evolution of transposable elements over short time scales. The method uses the alignment of all elements with intact gag/pol genes of a transposable element family from a single genome. We predict that the ratio of nonsynonymous to synonymous variants in the alignment should decrease as a function of the frequency of the variants, because elements with nonsynonymous variants that reduce transposition will have fewer progeny. We apply our method to Sirevirus long-terminal repeat retrotransposons that are abundant in maize and other plant species and show that nonsynonymous to synonymous variants declines as variant frequency increases, indicating that negative selection is acting strongly on the Sirevirus genome. The asymptotic value of nonsynonymous to synonymous variants suggests that at least 85% of all nonsynonymous mutations in the transposable element reduce transposition. Crucially, these patterns in nonsynonymous to synonymous variants are only predicted to occur if the gene products from a particular transposable element insertion preferentially promote the transposition of the same insertion. Overall, by using large numbers of intact elements, this study sheds new light on the selective processes that act on transposable elements. [ABSTRACT FROM AUTHOR]

Invariant sites are a common feature of amino acid sequence evolution. The presence of invariant sites is frequently attributed to the need to preserve function through site-specific conservation of amino acid residues. Amino acid substitution models without a provision for invariant sites often fit the data significantly worse than those that allow for an excess of invariant sites beyond those predicted by models that only incorporate rate variation among sites (e.g. a Gamma distribution). An alternative is epistasis between sites to preserve residue interactions that can create invariant sites. Through computer-simulated sequence evolution, we evaluated the relative effects of site-specific preferences and site-site couplings in the generation of invariant sites and the modulation of the rate of molecular evolution. In an analysis of ten major families of protein domains with diverse sequence and functional properties, we find that the negative selection imposed by epistasis creates many more invariant sites than site-specific residue preferences alone. Further, epistasis plays an increasingly larger role in creating invariant sites over longer evolutionary periods. Epistasis also dictates rates of domain evolution over time by exerting significant additional purifying selection to preserve site couplings. These patterns illuminate the mechanistic role of epistasis in the processes underlying observed site invariance and evolutionary rates. [ABSTRACT FROM AUTHOR]

The statistical associations between mutations, collectively known as linkage disequilibrium, encode important information about the evolutionary forces acting within a population. Yet in contrast to single-site analogues like the site frequency spectrum, our theoretical understanding of linkage disequilibrium remains limited. In particular, little is currently known about how mutations with different ages and fitness costs contribute to expected patterns of linkage disequilibrium, even in simple settings where recombination and genetic drift are the major evolutionary forces. Here, I introduce a forward-time framework for predicting linkage disequilibrium between pairs of neutral and deleterious mutations as a function of their present-day frequencies. I show that the dynamics of linkage disequilibrium become much simpler in the limit that mutations are rare, where they admit a simple heuristic picture based on the trajectories of the underlying lineages. I use this approach to derive analytical expressions for a family of frequency-weighted linkage disequilibrium statistics as a function of the recombination rate, the frequency scale, and the additive and epistatic fitness costs of the mutations. I find that the frequency scale can have a dramatic impact on the shapes of the resulting linkage disequilibrium curves, reflecting the broad range of time scales over which these correlations arise. I also show that the differences between neutral and deleterious linkage disequilibrium are not purely driven by differences in their mutation frequencies and can instead display qualitative features that are reminiscent of epistasis. I conclude by discussing the implications of these results for recent linkage disequilibrium measurements in bacteria. This forward-time approach may provide a useful framework for predicting linkage disequilibrium across a range of evolutionary scenarios. [ABSTRACT FROM AUTHOR]

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex , Daphnia pulicaria , and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa , where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion. [ABSTRACT FROM AUTHOR]

Abstract Human‐driven habitat fragmentation and loss have led to a proliferation of small and isolated plant and animal populations with high risk of extinction. One of the main threats to extinction in these populations is inbreeding depression, which is primarily caused by recessive deleterious mutations becoming homozygous due to inbreeding. The typical approach for managing these populations is to maintain high genetic diversity, increasingly by translocating individuals from large populations to initiate a “genetic rescue.” However, the limitations of this approach have recently been highlighted by the demise of the gray wolf population on Isle Royale, which declined to the brink of extinction soon after the arrival of a migrant from the large mainland wolf population. Here, we use a novel population genetic simulation framework to investigate the role of genetic diversity, deleterious variation, and demographic history in mediating extinction risk due to inbreeding depression in small populations. We show that, under realistic models of dominance, large populations harbor high levels of recessive strongly deleterious variation due to these mutations being hidden from selection in the heterozygous state. As a result, when large populations contract, they experience a substantially elevated risk of extinction after these strongly deleterious mutations are exposed by inbreeding. Moreover, we demonstrate that, although genetic rescue is broadly effective as a means to reduce extinction risk, its effectiveness can be greatly increased by drawing migrants from small or moderate‐sized source populations rather than large source populations due to smaller populations harboring lower levels of recessive strongly deleterious variation. Our findings challenge the traditional conservation paradigm that focuses on maximizing genetic diversity in small populations in favor of a view that emphasizes minimizing strongly deleterious variation. These insights have important implications for managing small and isolated populations in the increasingly fragmented landscape of the Anthropocene.

Baird’s tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as ‘endangered’ on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG’s). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 protein-coding genes (PCG’s), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (OL), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (OH). A majority of the PCG’s and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG’s show values

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1–Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.

Much research has focused on visual system evolution in bony fishes. The capacity of visual systems to perceive and respond to external signals is integral to evolutionary success. However, integrated research on the mechanisms of adaptive evolution based on corneal structure and related genes remains limited. In this study, scanning electron microscopy (SEM) was used to assess the microstructure and adaptation of corneal epithelial cells. Then, the evolution of the cornea-related COL8A2 gene was investigated. We found various projections (microridges, microplicae, microholes, and microvilli) on the corneal epithelial cells of amphibious mudskippers. Compared with those of fully aquatic fishes, these microstructures were considered adaptations to the variable environments experienced by amphibious mudskippers, as they can resist dryness in terrestrial environments and infection in aquatic environments. Moreover, strong purifying selection was detected for COL8A2. In addition, some specific amino acid substitution sites were also identified in the COL8A2 sequence in mudskippers. Interestingly, the evolutionary rate of the COL8A2 gene was significantly and positively correlated with maximum diving depth in our dataset. Specifically, with increasing diving depth, the evolutionary rate of the COL8A2 gene seemed to gradually accelerate. The results indicated that the cornea of bony fishes has evolved through adaptation to cope with the different diving depths encountered during the evolutionary process, with the corneal evolution of the amphibious mudskipper group showing a unique pattern.

The Out-of-Africa (OOA) dispersal ∼ 50,000 y ago is characterized by a series of founder events as modern humans expanded into multiple continents. Population genetics theory predicts an increase of mutational load in populations undergoing serial founder effects during range expansions. To test this hypothesis, we have sequenced full genomes and high-coverage exomes from seven geographically divergent human populations from Namibia, Congo, Algeria, Pakistan, Cambodia, Siberia, and Mexico. We find that individual genomes vary modestly in the overall number of predicted deleterious alleles. We show via spatially explicit simulations that the observed distribution of deleterious allele frequencies is consistent with the OOA dispersal, particularly under a model where deleterious mutations are recessive. We conclude that there is a strong signal of purifying selection at conserved genomic positions within Africa, but that many predicted deleterious mutations have evolved as if they were neutral during the expansion out of Africa. Under a model where selection is inversely related to dominance, we show that OOA populations are likely to have a higher mutation load due to increased allele frequencies of nearly neutral variants that are recessive or partially recessive.

Simple Summary: The evolution of distinct sex chromosomes, such as the X and Y chromosomes of humans, has long been of interest. Unlike most chromosomes, sex chromosomes do not fully recombine during meiosis. The non-recombining region is typically small early in their evolution, but additional regions stop recombining over time. Avian sex chromosomes, where males are ZZ and females ZW, show several evolutionary strata. Previous studies have demonstrated variation in barriers to recombination among bird species in the youngest strata. However, previous studies have sampled relatively few avian orders. Here, we examine one locus (ATP5F1A) from species in 22 different avian orders. Our results indicate that this locus has stopped recombining in almost all orders sampled. However, our results also show that the recombination stopped independently in each order of birds, highlighting the dynamic nature of avian sex chromosome evolution. The cessation of recombination is thought to have functional consequences because it is expected to weaken natural selection. We compared the proteins encoded by the recombining Z chromosome with those on the W chromosome. Most W chromosome proteins appeared to be functional, but they also appeared to be subject to weaker selection than the Z chromosome proteins. Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution. [ABSTRACT FROM AUTHOR]

Bacterial genomes are massively sequenced, and they provide valuable data to better know the complete set of genes of a species. The analysis of thousands of bacterial strains can identify both shared genes and those appearing only in the pathogenic ones. Current computational gene finders facilitate this task but often miss some existing genes. However, the present availability of different genomes from the same species is useful to estimate the selective pressure applied on genes of complete pangenomes. It may assist in evaluating gene predictions either by checking the certainty of a new gene or annotating it as a gene under positive selection. Here, we estimated the selective pressure of 19 271 genes that are part of the pangenome of the human opportunistic pathogen Acinetobacter baumannii and found that most genes in this bacterium are subject to negative selection. However, 23% of them showed values compatible with positive selection. These latter were mainly uncharacterized proteins or genes required to evade the host defence system including genes related to resistance and virulence whose changes may be favoured to acquire new functions. Finally, we evaluated the utility of measuring selection pressure in the detection of sequencing errors and the validation of gene prediction. [ABSTRACT FROM AUTHOR]

The rate of mutations varies >100-fold across the genome, altering the rate of evolution, and susceptibility to genetic diseases. The strongest predictor of mutation rate is the sequence itself, varying 75-fold between trinucleotides. The fact that DNA sequence drives its own mutation rate raises a simple but important prediction; highly mutable sequences will mutate more frequently and eliminate themselves in favor of sequences with lower mutability, leading to a lower equilibrium mutation rate. However, purifying selection constrains changes in mutable sequences, causing higher rates of mutation. We conduct a simulation using real human mutation data to test if 1) DNA evolves to a low equilibrium mutation rate and 2) purifying selection causes a higher equilibrium mutation rate in the genome's most important regions. We explore how this simple process affects sequence evolution in the genome, and discuss the implications for modeling evolution and susceptibility to DNA damage. [ABSTRACT FROM AUTHOR]

Mitochondrial genomes (mitogenomes) are involved in cellular energy metabolism and have been shown to undergo adaptive evolution in organisms with increased energy-consuming activities. The genetically selected high royal jelly-producing bees (RJBs, Apis mellifera ligustica) in China can produce 10 times more royal jelly, a highly nutritional and functional food, relative to unselected Italian bees (ITBs). To test for potential adaptive evolution of RJB mitochondrial genes, we sequenced mitogenomes from 100 RJBs and 30 ITBs. Haplotype network and phylogenetic analysis indicate that RJBs and ITBs are not reciprocally monophyletic but mainly divided into the RJB- and ITB-dominant sublineages. The RJB-dominant sublineage proportion is 6-fold higher in RJBs (84/100) than in ITBs (4/30), which is mainly attributable to genetic drift rather than positive selection. The RJB-dominant sublineage exhibits a low genetic diversity due to purifying selection. Moreover, mitogenome abundance is not significantly different between RJBs and ITBs, thereby rejecting the association between mitogenome copy number and royal jelly-producing performance. Our findings demonstrate low genetic diversity levels of RJB mitogenomes and reveal genetic drift and purifying selection as potential forces driving RJB mitogenome evolution. [ABSTRACT FROM AUTHOR]

Adaptation to derived habitats often occurs from standing genetic variation. The maintenance within ancestral populations of genetic variants favourable in derived habitats is commonly ascribed to long‐term antagonism between purifying selection and gene flow resulting from hybridization across habitats. A largely unexplored alternative idea based on quantitative genetic models of polygenic adaptation is that variants favoured in derived habitats are neutral in ancestral populations when their frequency is relatively low. To explore the latter, we first identify genetic variants important to the adaptation of threespine stickleback fish (Gasterosteus aculeatus) to a rare derived habitat—nutrient‐depleted acidic lakes—based on whole‐genome sequence data. Sequencing marine stickleback from six locations across the Atlantic Ocean then allows us to infer that the frequency of these derived variants in the ancestral habitat is unrelated to the likely opportunity for gene flow of these variants from acidic‐adapted populations. This result is consistent with the selective neutrality of derived variants within the ancestor. Our study thus supports an underappreciated explanation for the maintenance of standing genetic variation, and calls for a better understanding of the fitness consequences of adaptive variation across habitats and genomic backgrounds. [ABSTRACT FROM AUTHOR]

Whole-genome duplication (polyploidization) is among the most dramatic mutational processes in nature, so understanding how natural selection differs in polyploids relative to diploids is an important goal. Population genetics theory predicts that recessive deleterious mutations accumulate faster in allopolyploids than diploids due to the masking effect of redundant gene copies, but this prediction is hitherto unconfirmed. Here, we use the cotton genus (Gossypium), which contains seven allopolyploids derived from a single polyploidization event 1–2 Million years ago, to investigate deleterious mutation accumulation. We use two methods of identifying deleterious mutations at the nucleotide and amino acid level, along with whole-genome resequencing of 43 individuals spanning six allopolyploid species and their two diploid progenitors, to demonstrate that deleterious mutations accumulate faster in allopolyploids than in their diploid progenitors. We find that, unlike what would be expected under models of demographic changes alone, strongly deleterious mutations show the biggest difference between ploidy levels, and this effect diminishes for moderately and mildly deleterious mutations. We further show that the proportion of nonsynonymous mutations that are deleterious differs between the two coresident subgenomes in the allopolyploids, suggesting that homoeologous masking acts unequally between subgenomes. Our results provide a genome-wide perspective on classic notions of the significance of gene duplication that likely are broadly applicable to allopolyploids, with implications for our understanding of the evolutionary fate of deleterious mutations. Finally, we note that some measures of selection (e.g. d N /d S , πN/ π S) may be biased when species of different ploidy levels are compared. [ABSTRACT FROM AUTHOR]

High-throughput sequencing enables rapid genome sequencing during infectious disease outbreaks and provides an opportunity to quantify the evolutionary dynamics of pathogens in near real-time. One difficulty of undertaking evolutionary analyses over short timescales is the dependency of the inferred evolutionary parameters on the timespan of observation. Crucially, there are an increasing number of molecular clock analyses using external evolutionary rate priors to infer evolutionary parameters. However, it is not clear which rate prior is appropriate for a given time window of observation due to the time-dependent nature of evolutionary rate estimates. Here, we characterize the molecular evolutionary dynamics of SARS-CoV-2 and 2009 pandemic H1N1 (pH1N1) influenza during the first 12 months of their respective pandemics. We use Bayesian phylogenetic methods to estimate the dates of emergence, evolutionary rates, and growth rates of SARS-CoV-2 and pH1N1 over time and investigate how varying sampling window and data set sizes affect the accuracy of parameter estimation. We further use a generalized McDonald–Kreitman test to estimate the number of segregating nonneutral sites over time. We find that the inferred evolutionary parameters for both pandemics are time dependent, and that the inferred rates of SARS-CoV-2 and pH1N1 decline by ∼50% and ∼100%, respectively, over the course of 1 year. After at least 4 months since the start of sequence sampling, inferred growth rates and emergence dates remain relatively stable and can be inferred reliably using a logistic growth coalescent model. We show that the time dependency of the mean substitution rate is due to elevated substitution rates at terminal branches which are 2–4 times higher than those of internal branches for both viruses. The elevated rate at terminal branches is strongly correlated with an increasing number of segregating nonneutral sites, demonstrating the role of purifying selection in generating the time dependency of evolutionary parameters during pandemics. [ABSTRACT FROM AUTHOR]

Background: The 'Zacatuche', 'Teporingo', or Volcano rabbit (Romerolagus diazi) belongs to the family Leporidae, is an endemic species restricted to the Central part of the Trans-Mexican Volcanic Belt, and is considered 'endangered' by the IUCN Red List of Threatened Species. Methods and Results: This study reports, for the first time, the complete mitochondrial genome of R. diazi and examined the phylogenetic position of R. diazi among other closely related co-familiar species using mitochondrial protein-coding genes (PCGs). The mitogenome of R. diazi was assembled from short Illumina 150 bp pair-end reads with a coverage of 189x. The AT-rich mitochondrial genome of R. diazi is 17,400 bp in length and is comprised of 13 PCGs, two ribosomal RNA genes, and 22 transfer RNA genes. The gene order observed in the mitochondrial genome of R. diazi is identical to that reported for other leporids. Phylogenetic analyses based on PCGs support the basal position of Romerolagus within the Leporidae, at least when compared to the genera Oryctolagus and Lepus. Nonetheless, additional mitochondrial genomes from species belonging to the genera Bunolagus, Sylvilagus, and Pronolagus, among others, are needed before a more robust conclusion about the derived vs basal placement of Romerolagus within the family Leporidae can be reached based on mitochondrial PCGs. Conclusions: This is the first genomic resource developed for R. diazi and it represents a tool to improve our understanding about the ecology and evolutionary biology of this iconic and endangered species. [ABSTRACT FROM AUTHOR]

Abstract The emergence of SARS‐CoV‐2 has resulted in nearly 1,280,000 infections and 73,000 deaths globally so far. This novel virus acquired the ability to infect human cells using the SARS‐CoV cell receptor hACE2. Because of this, it is essential to improve our understanding of the evolutionary dynamics surrounding the SARS‐CoV‐2 hACE2 interaction. One way theory predicts selection pressures should shape viral evolution is to enhance binding with host cells. We first assessed evolutionary dynamics in select betacoronavirus spike protein genes to predict whether these genomic regions are under directional or purifying selection between divergent viral lineages, at various scales of relatedness. With this analysis, we determine a region inside the receptor‐binding domain with putative sites under positive selection interspersed among highly conserved sites, which are implicated in structural stability of the viral spike protein and its union with human receptor ACE2. Next, to gain further insights into factors associated with recognition of the human host receptor, we performed modeling studies of five different betacoronaviruses and their potential binding to hACE2. Modeling results indicate that interfering with the salt bridges at hot spot 353 could be an effective strategy for inhibiting binding, and hence for the prevention of SARS‐CoV‐2 infections. We also propose that a glycine residue at the receptor‐binding domain of the spike glycoprotein can have a critical role in permitting bat SARS‐related coronaviruses to infect human cells.

Abstract DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.

Abstract The endosperm cell walls of mature coffee seeds accumulate large amounts of mannan storage polysaccharides, which serve as nutrient reserve for embryo and contribute to beverage quality. Our study investigated the evolutionary patterns of key galactomannan (GM) biosynthesis genes using dN/dS ratio, synteny, and phylogenetic analysis and detected heterogeneity in rate of evolution among gene copies. Selection ratio index revealed evidence of positive selection in the branch editing gene Coffea canephora alpha (α) galactosidase (Cc‐alpha Gal) at Cc11_g15950 copy (ω = 1.12), whereas strong purifying selection on deleterious mutations was observed in the Coffea canephora uridine diphosphate (UDP)‐glucose 4′‐epimerase (Cc‐UG4E) and Coffea canephora mannose‐1P guanylytransferase (Cc‐MGT) genes controlling the crucial nucleotide carbon sugar building blocks flux in the pathway. Relatively low sequence diversity and strong syntenic linkages were detected in all GM pathway genes except in Cc‐alpha Gal, which suggests a correlation between selection pressure and nucleotide diversity or synteny analysis. In addition, phylogenetic analysis revealed independent evolution or expansion of GM pathway genes in different plant species, with no obvious inferable clustering patterns according to either gene family or congruent with evolutionary plants lineages tested due to high dynamic nature and specific biochemical cell wall modification requirements. Altogether, our study shows a significant high rate of evolutionary variation among GM pathway genes in the diploid C. canephora and demonstrates the inherent variation in evolution of gene copies and their potential role in understanding selection rates in a homogenously connected metabolic pathway.

Abstract Recombination's omnipresence in nature is one of the most intriguing problems in evolutionary biology. The question of why recombination exhibits certain general features is no less interesting than that of why it exists at all. One such feature is recombination's fitness dependence (FD). The so far developed population genetics models have focused on the evolution of FD recombination mainly in haploids, although the empirical evidence for this phenomenon comes mostly from diploids. Using numerical analysis of modifier models for infinite panmictic populations, we show here that FD recombination can be evolutionarily advantageous in diploids subjected to purifying selection. We ascribe this advantage to the differential rate of disruption of lower‐ versus higher‐fitness genotypes, which can be manifested in selected systems with at least three loci. We also show that if the modifier is linked to such selected system, it can additionally benefit from modifying this linkage in a fitness‐dependent manner. The revealed evolutionary advantage of FD recombination appeared robust to crossover interference within the selected system, either positive or negative. Remarkably, FD recombination was often favored in situations where any constant nonzero recombination was evolutionarily disfavored, implying a relaxation of the rather strict constraints on major parameters (e.g., selection intensity and epistasis) required for the evolutionary advantage of nonzero recombination formulated by classical models.

UDP- glucose: glycoprotein glucosyltransferase (UGGT) is a protein that operates as the gatekeeper for the endoplasmic reticulum (ER) quality control mechanism of glycoprotein folding. It is known that vertebrates and Caenorhabditis genomes harbor two uggt gene copies that exhibit differences in their properties.

Mitochondrial genomes (mitogenomes) are involved in cellular energy metabolism and have been shown to undergo adaptive evolution in organisms with increased energy-consuming activities. The genetically selected high royal jelly-producing bees (RJBs, Apis mellifera ligustica) in China can produce 10 times more royal jelly, a highly nutritional and functional food, relative to unselected Italian bees (ITBs). To test for potential adaptive evolution of RJB mitochondrial genes, we sequenced mitogenomes from 100 RJBs and 30 ITBs. Haplotype network and phylogenetic analysis indicate that RJBs and ITBs are not reciprocally monophyletic but mainly divided into the RJB- and ITB-dominant sublineages. The RJB-dominant sublineage proportion is 6-fold higher in RJBs (84/100) than in ITBs (4/30), which is mainly attributable to genetic drift rather than positive selection. The RJB-dominant sublineage exhibits a low genetic diversity due to purifying selection. Moreover, mitogenome abundance is not significantly different between RJBs and ITBs, thereby rejecting the association between mitogenome copy number and royal jelly-producing performance. Our findings demonstrate low genetic diversity levels of RJB mitogenomes and reveal genetic drift and purifying selection as potential forces driving RJB mitogenome evolution.

Macrocystis pyrifera (giant kelp), a haplodiplontic brown macroalga that alternates between a macroscopic diploid (sporophyte) and a microscopic haploid (gametophyte) phase, provides an ideal system to investigate how ploidy background affects the evolutionary history of a gene. In M. pyrifera, the same genome is subjected to different selective pressures and environments as it alternates between haploid and diploid life stages. We assembled M. pyrifera gene models using available expression data and validated 8,292 genes models using the model alga Ectocarpus siliculosus. Differential expression analysis identified gene models expressed in either or both the haploid and diploid life stages while functional annotation identified processes enriched in each stage. Genes expressed preferentially or exclusively in the gametophyte stage were found to have higher nucleotide diversity (π = 2.3 × 10–3 and 2.8 × 10–3, respectively) than those for sporophytes (π = 1.1 × 10–3 and 1 × 10–3, respectively). While gametophyte-biased genes show faster sequence evolution, the sequence evolution exhibits less signatures of adaptations when compared to sporophyte-biased genes. Our findings contrast the standing masking hypothesis, which predicts higher standing genetic variation at the sporophyte stage, and support the strength of expression theory, which posits that genes expressed more strongly are expected to evolve slower. We argue that the sporophyte stage undergoes more stringent selection compared with the gametophyte stage, which carries a heavy genetic load associated with broadcast spawning. Furthermore, using whole-genome sequencing, we confirm the strong population structure in wild M. pyrifera populations previously established using microsatellite markers, and estimate population genetic parameters, such as pairwise genetic diversity and Tajima’s D, important for conservation and domestication of M. pyrifera.