Evolutionary innovations generate phenotypic and species diversity. Elucidating the genomic processes underlying such innovations is central to understanding biodiversity. In this study, we addressed the genomic basis of evolutionary novelties in the glassy-winged sharpshooter (Homalodisca vitripennis, GWSS), an agricultural pest. Prominent evolutionary innovations in leafhoppers include brochosomes, proteinaceous structures that are excreted and used to coat the body, and obligate symbiotic associations with two bacterial types that reside within cytoplasm of distinctive cell types. Using PacBio long-read sequencing and Dovetail Omni-C technology, we generated a chromosome-level genome assembly for the GWSS and then validated the assembly using flow cytometry and karyotyping. Additional transcriptomic and proteomic data were used to identify novel genes that underlie brochosome production. We found that brochosome-associated genes include novel gene families that have diversified through tandem duplications. We also identified the locations of genes involved in interactions with bacterial symbionts. Ancestors of the GWSS acquired bacterial genes through horizontal gene transfer (HGT), and these genes appear to contribute to symbiont support. Using a phylogenomics approach, we inferred HGT sources and timing. We found that some HGT events date to the common ancestor of the hemipteran suborder Auchenorrhyncha, representing some of the oldest known examples of HGT in animals. Overall, we show that evolutionary novelties in leafhoppers are generated by the combination of acquiring novel genes, produced both de novo and through tandem duplication, acquiring new symbiotic associations that enable use of novel diets and niches, and recruiting foreign genes to support symbionts and enhance herbivory.

The origin of nitrogen fixation is an important issue in evolutionary biology. While nitrogen is required by all living organisms, only a small fraction of bacteria and archaea can fix nitrogen. The prevailing view is that nitrogen fixation first evolved in archaea and was later transferred to bacteria. However, nitrogen-fixing (Nif) bacteria are far larger in number and far more diverse in ecological niches than Nif archaea. We, therefore, propose the bacteria-first hypothesis, which postulates that nitrogen fixation first evolved in bacteria and was later transferred to archaea. As >30,000 prokaryotic genomes have been sequenced, we conduct an in-depth comparison of the two hypotheses. We first identify the six genes involved in nitrogen fixation in all sequenced prokaryotic genomes and then reconstruct phylogenetic trees using the six Nif proteins individually or in combination. In each of these trees, the earliest lineages are bacterial Nif protein sequences and in the oldest clade (group) the archaeal sequences are all nested inside bacterial sequences, suggesting that the Nif proteins first evolved in bacteria. The bacteria-first hypothesis is further supported by the observation that the majority of Nif archaea carry the major bacterial Mo (molybdenum) transporter (ModABC) rather than the archaeal Mo transporter (WtpABC). Moreover, in our phylogeny of all available ModA and WtpA protein sequences, the earliest lineages are bacterial sequences while archaeal sequences are nested inside bacterial sequences. Furthermore, the bacteria-first hypothesis is supported by available isotopic data. In conclusion, our study strongly supports the bacteria-first hypothesis.

Sirtuins are a family of proteins that protect against cellular injury and aging; understanding their evolution should reveal fundamental mechanisms governing longevity. "Early-branching" animals such as sea sponges and jellyfish have been understudied in previous analyses of sirtuin diversity. These organisms not only hold important positions at the base of the evolutionary tree, but also have unique aging dynamics that defy convention, such as quasi-immortality and high regenerative capacity. In this study, we survey the evolution of sirtuin proteins in animals, with a focus on the oldest living lineages. We describe previously unrecognized expansions of "Class IV" and "Class I" sirtuins around the origin of animals, raising the number of sirtuin families in the last common ancestor to at least nine. Most of these undescribed sirtuins have been lost in vertebrates and other bilaterian animals. Our work also clarifies the evolution of PNC1 and NAMPT enzymes that carry out the rate-limiting step in sirtuin-related NAD+ biosynthesis. The genes for PNC1 and NAMPT enzymes were both present in the first animals, with the genes being lost a minimum of 11 and 13 times, respectively, over the course of animal evolution. We propose that species with these ancestral gene repertoires are ideal model organisms for studying the genetic regulation of animal longevity and will provide clues to increasing longevity in humans.

Biological nitrogen fixation (BNF) by cyanobacteria is of significant importance for the Earth's biogeochemical nitrogen cycle but is restricted to a few genera that do not form monophyletic group. To explore the evolutionary trajectory of BNF and investigate the driving forces of its evolution, we analyze 650 cyanobacterial genomes and compile the database of diazotrophic cyanobacteria based on the presence of nitrogen fixation gene clusters (NFGCs). We report that 266 of 650 examined genomes are NFGC-carrying members, and these potentially diazotrophic cyanobacteria are unevenly distributed across the phylogeny of Cyanobacteria, that multiple independent losses shaped the scattered distribution. Among the diazotrophic cyanobacteria, two types of NFGC exist, with one being ancestral and abundant, which have descended from diazotrophic ancestors, and the other being anaerobe-like and sparse, possibly being acquired from anaerobic microbes through horizontal gene transfer. Interestingly, we illustrate that the origin of BNF in Cyanobacteria coincide with two major evolutionary events. One is the origin of multicellularity of cyanobacteria, and the other is concurrent genetic innovations with massive gene gains and expansions, implicating their key roles in triggering the evolutionary transition from nondiazotrophic to diazotrophic cyanobacteria. Additionally, we reveal that genes involved in accelerating respiratory electron transport (coxABC), anoxygenic photosynthetic electron transport (sqr), as well as anaerobic metabolisms (pfor, hemN, nrdG, adhE) are enriched in diazotrophic cyanobacteria, representing adaptive genetic signatures that underpin the diazotrophic lifestyle. Collectively, our study suggests that multicellularity, together with concurrent genetic adaptations contribute to the evolution of diazotrophic cyanobacteria.

Deleterious genetic variation is maintained in populations at low frequencies. Under a model of stabilizing selection, rare (and presumably deleterious) genetic variants are associated with increase or decrease in gene expression from some intermediate optimum. We investigate this phenomenon in a population of largely Oryza sativa ssp. indica rice landraces under normal unstressed wet and stressful drought field conditions. We include single nucleotide polymorphisms, insertion/deletion mutations, and structural variants in our analysis and find a stronger association between rare variants and gene expression outliers under the stress condition. We also show an association of the strength of this rare variant effect with linkage, gene expression levels, network connectivity, local recombination rate, and fitness consequence scores, consistent with the stabilizing selection model of gene expression.

The emergence of the placenta is a revolutionary event in the evolution of therian mammals, to which some LTR retroelement-derived genes, such as PEG10, RTL1, and syncytin, are known to contribute. However, therian genomes contain many more LTR retroelement-derived genes that may also have contributed to placental evolution. We conducted large-scale evolutionary genomic and transcriptomic analyses to comprehensively search for LTR retroelement-derived genes whose origination coincided with therian placental emergence and that became consistently expressed in therian placentae. We identified NYNRIN as another Ty3/Gypsy LTR retroelement-derived gene likely to contribute to placental emergence in the therian stem lineage. NYNRIN knockdown inhibited the invasion of HTR8/SVneo invasive-type trophoblasts, whereas the knockdown of its nonretroelement-derived homolog KHNYN did not. Functional enrichment analyses suggested that NYNRIN modulates trophoblast invasion by regulating epithelial-mesenchymal transition and extracellular matrix remodeling and that the ubiquitin-proteasome system is responsible for the functional differences between NYNRIN and KHNYN. These findings extend our knowledge of the roles of LTR retroelement-derived genes in the evolution of therian mammals.

The diversity of Central Asians has been shaped by multiple migrations and cultural diffusion. Although ancient DNA studies have revealed the demographic changes of the Central Asian since the Bronze Age, the contribution of the ancient populations to the modern Central Asian remains opaque. Herein, we performed high-coverage sequencing of 131 whole genomes of Indo-European-speaking Tajik and Turkic-speaking Kyrgyz populations to explore their genomic diversity and admixture history. By integrating the ancient DNA data, we revealed more details of the origins and admixture history of Central Asians. We found that the major ancestry of present-day Tajik populations can be traced back to the admixture of the Bronze Age Bactria-Margiana Archaeological Complex and Andronovo-related populations. Highland Tajik populations further received additional gene flow from the Tarim mummies, an isolated ancient North Eurasian-related population. The West Eurasian ancestry of Kyrgyz is mainly derived from Historical Era populations in Xinjiang of China. Furthermore, the recent admixture signals detected in both Tajik and Kyrgyz are ascribed to the expansions of Eastern Steppe nomadic pastoralists during the Historical Era.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Bacteria and lytic viruses (phages) engage in highly dynamic coevolutionary interactions over time, yet we have little idea of how transient selection by phages might shape the future evolutionary trajectories of their host populations. To explore this question, we generated genetically diverse phage-resistant mutants of the bacterium Pseudomonas syringae. We subjected the panel of mutants to prolonged experimental evolution in the absence of phages. Some populations re-evolved phage sensitivity, whereas others acquired compensatory mutations that reduced the costs of resistance without altering resistance levels. To ask whether these outcomes were driven by the initial genetic mechanisms of resistance, we next evolved independent replicates of each individual mutant in the absence of phages. We found a strong signature of historical contingency: some mutations were highly reversible across replicate populations, whereas others were highly entrenched. Through whole-genome sequencing of bacteria over time, we also found that populations with the same resistance gene acquired more parallel sets of mutations than populations with different resistance genes, suggesting that compensatory adaptation is also contingent on how resistance initially evolved. Our study identifies an evolutionary ratchet in bacteria-phage coevolution and may explain previous observations that resistance persists over time in some bacterial populations but is lost in others. We add to a growing body of work describing the key role of phages in the ecological and evolutionary dynamics of their host communities. Beyond this specific trait, our study provides a new insight into the genetic architecture of historical contingency, a crucial component of interpreting and predicting evolution.

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor-binding domain and receptor-binding motif. An excess of high-frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low-frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies that the virus may not be preadapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.

Bacterial adaptation to stressful environments often produces evolutionary constraints whereby increases in resistance are associated with reduced fitness in a different environment. The exploitation of this resistance-cost trade-off has been proposed as the basis of rational antimicrobial treatment strategies designed to limit the evolution of drug resistance in bacterial pathogens. Recent theoretical, laboratory, and clinical studies have shown that fluctuating selection can maintain drug efficacy and even restore drug susceptibility, but can also increase the rate of adaptation and promote cross-resistance to other antibiotics. In this paper, we combine mathematical modeling, experimental evolution, and whole-genome sequencing to follow evolutionary trajectories towards β-lactam resistance under fluctuating selective conditions. Our experimental model system consists of eight populations of Escherichia coli K12 evolving in parallel to a serial dilution protocol designed to dynamically control the strength of selection for resistance. We implemented adaptive ramps with mild and strong selection, resulting in evolved populations with similar levels of resistance, but with different evolutionary dynamics and diverging genotypic profiles. We found that mutations that emerged under strong selection are unstable in the absence of selection, in contrast to resistance mutations previously selected in the mild selection regime that were stably maintained in drug-free environments and positively selected for when antibiotics were reintroduced. Altogether, our population dynamics model and the phenotypic and genomic analysis of the evolved populations show that the rate of resistance adaptation is contingent upon the strength of selection, but also on evolutionary constraints imposed by prior drug exposures.

Unisexual animals are commonly found in some polyploid species complexes, and most of these species have had a long evolutionary history. However, their method for avoiding genomic decay remains unclear. The polyploid Carassius complex naturally comprises the sexual amphidiploid C. auratus (crucian carp or goldfish) (AABB) and the gynogenetic amphitriploid C. gibelio (gibel carp) (AAABBB). Recently, we developed a fertile synthetic amphitetraploid (AAAABBBB) male from C. gibelio by incorporating a C. auratus genome. In this study, we generated novel amphitriploids (AAABBB) by backcrossing the amphitetraploid male with the amphidiploid C. auratus. Whole-genome resequencing revealed the genomic changes, including recombination and independent assortment between homologs of C. gibelio and C. auratus. The fertility, sex determination system, oocyte development, and fertilization behaviors of the novel amphitriploids were investigated. Approximately 80% of the novel amphitriploid females recovered the unisexual gynogenesis ability. Intriguingly, two types of primary oocyte (with and without homolog synapsis) were discovered, and their distinct development fates were observed. Type I oocytes entered apoptosis due to improper synaptonemal complex assembly and incomplete double-strand break repair, whereas subsequent type II oocytes bypassed meiosis through an alternative ameiotic pathway to develop into mature eggs. Moreover, gynogenesis was stabilized in their offspring, and a new array of diverse gynogenetic amphitriploid clones was produced. These revealed genomic changes and detailed cytological data provide comprehensive evidence that changes in ploidy drive unisexual and sexual reproduction transition, thereby resulting in genomic diversity and allowing C. gibelio avoid genomic decay.

Fitness landscape mapping and the prediction of evolutionary trajectories on these landscapes are major tasks in evolutionary biology research. Evolutionary dynamics is tightly linked to the landscape topography, but this relation is not straightforward. Here, we analyze a fitness landscape of a yeast tRNA gene, previously measured under four different conditions. We find that the wild type allele is sub-optimal, and 8-10% of its variants are fitter. We rule out the possibilities that the wild type is fittest on average on these four conditions or located on a local fitness maximum. Notwithstanding, we cannot exclude the possibility that the wild type might be fittest in some of the many conditions in the complex ecology that yeast lives at. Instead, we find that the wild type is mutationally robust ("flat"), while more fit variants are typically mutationally fragile. Similar observations of mutational robustness or flatness have been so far made in very few cases, predominantly in viral genomes.

Quantifying the magnitude of the global extinction crisis is important but remains challenging, as many extinction events pass unnoticed owing to our limited taxonomic knowledge of the world's organisms. The increasing rarity of many taxa renders comprehensive sampling difficult, further compounding the problem. Vertebrate lineages such as birds, which are thought to be taxonomically well understood, are therefore used as indicator groups for mapping and quantifying global extinction. To test whether extinction patterns are adequately gauged in well-studied groups, we implemented ancient-DNA protocols and retrieved whole genomes from the historic DNA of museum specimens in a widely known songbird radiation of shamas (genus Copsychus) that is assumed to be of least conservation concern. We uncovered cryptic diversity and an unexpected degree of hidden extinction and terminal endangerment. Our analyses reveal that >40% of the phylogenetic diversity of this radiation is already either extinct in the wild or nearly so, including the two genomically most distinct members of this group (omissus and nigricauda), which have so far flown under the conservation radar as they have previously been considered subspecies. Comparing the genomes of modern samples with those from roughly a century ago, we also found a significant decrease in genetic diversity and a concomitant increase in homozygosity affecting various taxa, including small-island endemics that are extinct in the wild as well as subspecies that remain widespread across the continental scale. Our application of modern genomic approaches demonstrates elevated levels of allelic and taxonomic diversity loss in a songbird clade that has not been listed as globally threatened, highlighting the importance of ongoing reassessments of extinction incidence even across well-studied animal groups. Key words: extinction, introgression, white-rumped shama, conservation.

Host-dependency factors have increasingly been targeted to minimize antiviral drug resistance. In this study, we have demonstrated that inhibition of p38 mitogen-activated protein kinase (a cellular protein) suppresses buffalopox virus (BPXV) protein synthesis by targeting p38-MNK1-eIF4E signaling pathway. In order to provide insights into the evolution of drug resistance, we selected resistant mutants by long-term sequential passages (P; n = 60) in the presence of p38 inhibitor (SB239063). The P60-SB239063 virus exhibited significant resistance to SB239063 as compared to the P60-Control virus. To provide mechanistic insights on the acquisition of resistance by BPXV-P60-SB239063, we generated p38-α and p38-ϒ (isoforms of p38) knockout Vero cells by CRISPR/Cas9-mediated genome editing. It was demonstrated that unlike the wild type (WT) virus which is dependent on p38-α isoform, the resistant virus (BPXV-P60-SB239063) switches over to use p38-ϒ so as to efficiently replicate in the target cells. This is a rare evidence wherein a virus was shown to bypass the dependency on a critical cellular factor under selective pressure of a drug.

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.

Erratum in Mol Biol Evol. 2022 Nov 3;39(11):

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.