Two decades of suspect evidence for adaptive molecular evolution – Negative selection confounding positive selection signals

Characterization of cancer-driving nucleotides (CDNs) across genes, cancer types, and patients

A central goal of cancer genomics is to identify, in each patient, all the cancer-driving mutations. Among them, point mutations are referred to as cancer-driving nucleotides (CDNs), which recur in cancers. The companion study shows that the probability of i recurrent hits in n patients would decrease exponentially with i; hence, any mutation with i ≥ 3 hits in The Cancer Genome Atlas (TCGA) database is a high-probability CDN. This study characterizes the 50-150 CDNs identifiable for each cancer type of TCGA (while anticipating 10 times more undiscovered ones) as follows: (i) CDNs tend to code for amino acids of divergent chemical properties. (ii) At the genic level, far more CDNs (more than fivefold) fall on noncanonical than canonical cancer-driving genes (CDGs). Most undiscovered CDNs are expected to be on unknown CDGs. (iii) CDNs tend to be more widely shared among cancer types than canonical CDGs, mainly because of the higher resolution at the nucleotide than the whole-gene level. (iv) Most important, among the 50-100 coding region mutations carried by a cancer patient, 5-8 CDNs are expected but only 0-2 CDNs have been identified at present. This low level of identification has hampered functional test and gene-targeted therapy. We show that, by expanding the sample size to 105, most CDNs can be identified. Full CDN identification will then facilitate the design of patient-specific targeting against multiple CDN-harboring genes.

The theory of massively repeated evolution and full identifications of cancer-driving nucleotides (CDNs)

Tumorigenesis, like most complex genetic traits, is driven by the joint actions of many mutations. At the nucleotide level, such mutations are cancer-driving nucleotides (CDNs). The full sets of CDNs are necessary, and perhaps even sufficient, for the understanding and treatment of each cancer patient. Currently, only a small fraction of CDNs is known as most mutations accrued in tumors are not drivers. We now develop the theory of CDNs on the basis that cancer evolution is massively repeated in millions of individuals. Hence, any advantageous mutation should recur frequently and, conversely, any mutation that does not is either a passenger or deleterious mutation. In the TCGA cancer database (sample size n=300-1000), point mutations may recur in i out of n patients. This study explores a wide range of mutation characteristics to determine the limit of recurrences (i*) driven solely by neutral evolution. Since no neutral mutation can reach i*=3, all mutations recurring at i≥3 are CDNs. The theory shows the feasibility of identifying almost all CDNs if n increases to 100,000 for each cancer type. At present, only <10% of CDNs have been identified. When the full sets of CDNs are identified, the evolutionary mechanism of tumorigenesis in each case can be known and, importantly, gene targeted therapy will be far more effective in treatment and robust against drug resistance.

Bona Fide Plant Steroid Receptors are Innovated in Seed Plants and Angiosperms through Successive Whole-Genome Duplication Events

Whole-genome duplication (WGD) events are widespread in plants and animals, thus their long-term evolutionary contribution has long been speculated, yet a specific contribution is difficult to verify. Here, we show that ɛ-WGD and ζ-WGD contribute to the origin and evolution of bona fide brassinosteroid (BR) signaling through the innovation of active BR biosynthetic enzymes and active BR receptors from their respective ancestors. We found that BR receptors BRI1 (BR INSENSITIVE 1) and BRL1/3 (BRI1-LIKES 1/3) derived by ɛ-WGD and ζ-WGD, which occurred in the common ancestor of angiosperms and seed plants, respectively, while orphan BR receptor BRL2 first appeared in stomatophytes. Additionally, CYP85A enzymes synthesizing the bioactive BRs derived from a common ancestor of seed plants, while its sister enzymes CYP90 synthesizing BR precursors presented in all land plants, implying possible ligand-receptor coevolution. Consistently, the island domains (IDs) responsible for BR perception in BR receptors were most divergent among different receptor branches, supporting ligand-driven evolution. As a result, BRI1 was the most diversified BR receptor in angiosperms. Importantly, relative to the BR biosynthetic DET2 gene presented in all land plants, BRL2, BRL1/3 and BRI1 had high expression in vascular plants ferns, gymnosperms and angiosperms, respectively. Notably, BRI1 is the most diversified BR receptor with the most abundant expression in angiosperms, suggesting potential positive selection. Therefore, WGDs initiate a neofunctionalization process diverged by ligand-perception and transcriptional expression, which might optimize both BR biosynthetic enzymes and BR receptors, likely contributing to the evolution of land plants, especially seed plants and angiosperms.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Intra- vs. Interhost Evolution of SARS-CoV-2 Driven by Uncorrelated Selection-The Evolution Thwarted

In viral evolution, a new mutation has to proliferate within the host (Stage I) in order to be transmitted and then compete in the host population (Stage II). We now analyze the intrahost single nucleotide variants (iSNVs) in a set of 79 SARS-CoV-2 infected patients with most transmissions tracked. Here, every mutation has two measures: 1) iSNV frequency within each individual host in Stage I; 2) occurrence among individuals ranging from 1 (private), 2-78 (public), to 79 (global) occurrences in Stage II. In Stage I, a small fraction of nonsynonymous iSNVs are sufficiently advantageous to rise to a high frequency, often 100%. However, such iSNVs usually fail to become public mutations. Thus, the selective forces in the two stages of evolution are uncorrelated and, possibly, antagonistic. For that reason, successful mutants, including many variants of concern, have to avoid being eliminated in Stage I when they first emerge. As a result, they may not have the transmission advantage to outcompete the dominant strains and, hence, are rare in the host population. Few of them could manage to slowly accumulate advantageous mutations to compete in Stage II. When they do, they would appear suddenly as in each of the six successive waves of SARS-CoV-2 strains. In conclusion, Stage I evolution, the gate-keeper, may contravene the long-term viral evolution and should be heeded in viral studies.

Genes and sites under adaptation at the phylogenetic scale also exhibit adaptation at the population-genetic scale

Adaptation in protein-coding sequences can be detected from multiple sequence alignments across species or alternatively by leveraging polymorphism data within a population. Across species, quantification of the adaptive rate relies on phylogenetic codon models, classically formulated in terms of the ratio of nonsynonymous over synonymous substitution rates. Evidence of an accelerated nonsynonymous substitution rate is considered a signature of pervasive adaptation. However, because of the background of purifying selection, these models are potentially limited in their sensitivity. Recent developments have led to more sophisticated mutation-selection codon models aimed at making a more detailed quantitative assessment of the interplay between mutation, purifying, and positive selection. In this study, we conducted a large-scale exome-wide analysis of placental mammals with mutation-selection models, assessing their performance at detecting proteins and sites under adaptation. Importantly, mutation-selection codon models are based on a population-genetic formalism and thus are directly comparable to the McDonald and Kreitman test at the population level to quantify adaptation. Taking advantage of this relationship between phylogenetic and population genetics analyses, we integrated divergence and polymorphism data across the entire exome for 29 populations across 7 genera and showed that proteins and sites detected to be under adaptation at the phylogenetic scale are also under adaptation at the population-genetic scale. Altogether, our exome-wide analysis shows that phylogenetic mutation-selection codon models and the population-genetic test of adaptation can be reconciled and are congruent, paving the way for integrative models and analyses across individuals and populations.

Canalization of Phenotypes-When the Transcriptome is Constantly but Weakly Perturbed

Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

Positive selection-driven fixation of a hominin-specific amino acid mutation related to dephosphorylation in IRF9

The arms race between humans and pathogens drives the evolution of the human genome. It is thus expected that genes from the interferon-regulatory factors family (IRFs), a critical family for anti-viral immune response, should be undergoing episodes of positive selection. Herein, we tested this hypothesis and found multiple lines of evidence for positive selection on the amino acid site Val129 (NP_006075.3:p.Ser129Val) of human IRF9. Interestingly, the ancestral reconstruction and population distribution analyses revealed that the ancestral state (Ser129) is conserved among mammals, while the derived positively selected state (Val129) was fixed before the “out-of-Africa” event ~ 500,000 years ago. The motif analysis revealed that this young amino acid (Val129) may serve as a dephosphorylation site of IRF9. Structural parallelism between homologous genes further suggested the functional effects underlying the dephosphorylation that may affect the immune activity of IRF9. This study provides a model in which a strong positive Darwinian selection drives a recent fixation of a hominin-specific amino acid leading to molecular adaptation involving dephosphorylation in an immune-responsive gene.