Search

Pathogens were a major driver of genetic adaptation during human evolution. Viruses in particular were a dominant driver of adaptation in the thousands of proteins that physically interact with viruses (VIPs for Virus-Interacting Proteins). This however poses a conundrum. The best understood cases of virus-driven adaptation in specialized immune antiviral factors or in host viral receptors are numerically vastly insufficient to explain abundant adaptations in VIPs. What adaptive mechanisms can then at least partly close this gap? VIPs tend to be broadly conserved proteins with conserved host native molecular functions. Because many amino acid changes in a protein can alter its stability –the balance between the folded and unfolded forms of a protein– without destroying conserved native activities, here we ask if stability evolution was an important mechanism of virus-driven human protein adaptation. Using predictions of protein stability changes based on Alphafold 2 structures and validated by multiple lines of evidence, we find that amino acid changes that altered stability experienced highly elevated adaptative evolution in VIPs, compared to changes with a weaker impact on stability. We further find that RNA viruses, rather DNA viruses, predominantly drove strong adaptation through stability changes in VIPs. We also observe that stability in immune antiviral VIPs preferentially evolved under directional selection. Conversely, stability in proviral VIPs needed by viruses evolved under compensatory evolution following viral epidemics. Together, these results suggest that stability evolution, and thus functional host protein abundance evolution, was a prominent mechanism of host protein adaptation during viral epidemics.

Advances in genome sequencing have improved our understanding of the genetic basis of human diseases, and thousands of human genes have been associated with different diseases. Recent genomic adaptation at disease genes has not been well characterized. Here, we compare the rate of strong recent adaptation in the form of selective sweeps between mendelian, non-infectious disease genes and non-disease genes across distinct human populations from the 1000 Genomes Project. We find that mendelian disease genes have experienced far less selective sweeps compared to non-disease genes especially in Africa. Investigating further the possible causes of the sweep deficit at disease genes, we find that this deficit is very strong at disease genes with both low recombination rates and with high numbers of associated disease variants, but is almost non-existent at disease genes with higher recombination rates or lower numbers of associated disease variants. Because segregating recessive deleterious variants have the ability to interfere with adaptive ones, these observations strongly suggest that adaptation has been slowed down by the presence of interfering recessive deleterious variants at disease genes. These results suggest that disease genes suffer from a transient inability to adapt as fast as the rest of the genome.

Advances in genome sequencing have dramatically improved our understanding of the genetic basis of human diseases, and thousands of human genes have been associated with different diseases. Despite our expanding knowledge of gene-disease associations, and despite the medical importance of disease genes, their evolution has not been thoroughly studied across diverse human populations. In particular, recent genomic adaptation at disease genes has not been well characterized, even though multiple evolutionary processes are expected to connect disease and adaptation at the gene level. Understanding the relationship between disease and adaptation at the gene level in the human genome is severely hampered by the fact that we don’t even know whether disease genes have experienced more, less, or as much adaptation as non-disease genes during recent human evolution. Here, we compare the rate of strong recent adaptation in the form of selective sweeps between disease genes and non-disease genes across 26 distinct human populations from the 1,000 Genomes Project. We find that disease genes have experienced far less selective sweeps compared to non-disease genes during recent human evolution. This sweep deficit at disease genes is particularly visible in Africa, and less visible in East Asia or Europe, likely due to more intense genetic drift in the latter populations creating more spurious selective sweeps signals. Investigating further the possible causes of the sweep deficit at disease genes, we find that this deficit is very strong at disease genes with both low recombination rates and with high numbers of associated disease variants, but is inexistant at disease genes with higher recombination rates or lower numbers of associated disease variants. Because recessive deleterious variants have the ability to interfere with adaptive ones, these observations strongly suggest that adaptation has been slowed down by the presence of interfering recessive deleterious variants at disease genes. These results clarify the evolutionary relationship between disease genes and recent genomic adaptation, and suggest that disease genes suffer not only from a higher load of segregating deleterious mutations, but also an inability to adapt as much, and/or as fast as the rest of the genome.

SummaryMany crops are polyploid or have a polyploid ancestry. Recent phylogenetic analyses have found that polyploidy often preceded the domestication of crop plants. One explanation for this observation is that increased genetic diversity following polyploidy may have been important during the strong artificial selection that occurs during domestication.To test the connection between domestication and polyploidy, we identified and examined candidate genes associated with the domestication of the diverse crop varieties ofBrassica rapa. Like all “diploid” flowering plants,B. rapahas a diploidized paleopolyploid genome and experienced many rounds of whole genome duplication (WGD). We analyzed transcriptome data of more than 100 cultivatedB. rapaaccessions.Using a combination of approaches, we identified more than 3,000 candidate genes associated with the domestication of four majorB. rapacrop varieties. Consistent with our expectation, we found that the candidate genes were significantly enriched with genes derived from the Brassiceae mesohexaploidy. We also observed that paleologs were significantly more diverse than non-paleologs.Our analyses find evidence for that genetic diversity derived from ancient polyploidy played a key role in the domestication ofB. rapaand provide support for its importance in the success of modern agriculture.

The diversity of life is shaped by rates of speciation and extinction, and so estimating these rates correctly is crucial for understanding diversity patterns among clades, regions, and habitats. In 2011, Morlon and collaborators developed a promising likelihood-based approach to estimate speciation and extinction and to infer the model describing how these rates change over time based on AICc. This approach is now implemented in an R package (RPANDA). Here, we test the accuracy of this approach under simulated conditions, to evaluate its ability to correctly estimate rates of speciation, extinction, and diversification (speciation—extinction) and to choose the correct underlying model of diversification (e.g. constant or changing rates of speciation and extinction over time). We found that this likelihood-based approach frequently picked the incorrect model. For example, with changing speciation rates over time, the correct model was chosen in only ∼10 per cent of replicates. There were significant relationships between true and estimated speciation rates using this approach, but relationships were weak when speciation rates were constant within clades. Relationships were consistently weak between true and estimated rates of extinction and of diversification. Overall, we suggest that results from this approach should be interpreted with considerable caution.

Chinese cymbidiums (Cymbidium sp.) are important ornamental plants because of their foliage, flower shape, and fragrance. Well-known Chinese cymbidiums mainly include Cymbidium goeringii, Cymbidium faberi, Cymbidium ensifolium, Cymbidium kanran, and Cymbidium sinense. The population genetics of Chinese cymbidiums can be efficiently analyzed using small-scale marker panels with high discriminatory power. In this study, we tested several genic simple sequence repeats (SSRs) and built six genic SSR panels. The panels included several robust markers, which can rapidly assign Chinese cymbidium accessions to their source species. Fifty-three accessions of Chinese cymbidiums were analyzed using 25 markers, which exhibited polymorphism among five species. These markers were ranked according to their discriminatory scores (D scores). The program selected six markers to build an “overall” panel for all Cymbidium classifications and yielded 95.16% population assignment accuracy. Considering one species as the “critical” population and the four other species as one population, we built five genic SSR panels: C. ensifolium panel (four markers, 98.05% accuracy), C. faberi panel (six markers, 95.90% accuracy), C. goeringii panel (six markers, 95.15% accuracy), C. sinense panel (six markers, 96.35% accuracy), and C. kanran panel (five markers, 96.10% accuracy). Genetic distance matrices calculated using the “overall” panels and those derived with the 25 markers were compared. Results showed a high correlation (R = 0.807) with statistical significance (P = 0.042). Moreover, “all panels” revealed higher genetic variations among populations than “all markers.” Hence, the developed panels are suitable for efficient population classification of Chinese cymbidiums.