A mutation-selection model of protein evolution under persistent positive selection

Increased Positive Selection in Highly Recombining Genes Does not Necessarily Reflect an Evolutionary Advantage of Recombination

It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.

Quantifying microbial guilds

The ecological role of microorganisms is of utmost importance due to their multiple interactions with the environment. However, assessing the contribution of individual taxonomic groups has proven difficult despite the availability of high throughput data, hindering our understanding of such complex systems. Here, we propose a quantitative definition of guild that is readily applicable to metagenomic data. Our framework focuses on the functional character of protein sequences, as well as their diversifying nature. First, we discriminate functional sequences from the whole sequence space corresponding to a gene annotation to then quantify their contribution to the guild composition across environments. In addition, we identify and distinguish functional implementations, which are sequence spaces that have different ways of carrying out the function. In contrast, we found that orthology delineation did not consistently align with ecologically (or functionally) distinct implementations of the function. We demonstrate the value of our approach with two case studies: the ammonia oxidation and polyamine uptake guilds from the Malaspina circumnavigation cruise, revealing novel ecological dynamics of the latter in marine ecosystems. Thus, the quantification of guilds helps us to assess the functional role of different taxonomic groups with profound implications on the study of microbial communities.

Comparative genomics analysis provide insights into evolution and stress responses of Lhcb genes in Rosaceae fruit crops

BACKGROUND

Light-harvesting chlorophyll a/b b evelopment of higher plants and in response to abiotic stress. Previous works has demonstrated that that Lhcb genes were involved in the phytochrome regulation and responded to the different light and temperature conditions in Poaceae (such as maize). However, the evolution and functions of Lhcb genes remains poorly characterized in important Rosaceae species.

RESULTS

In this investigation, we conducted a genome-wide analysis and identified a total of 212 Lhcb genes across nine Rosaceae species. Specifically, we found 23 Lhcb genes in Fragaria vesca, 20 in Prunus armeniaca, 33 in Malus domestica 'Gala', 21 in Prunus persica, 33 in Rosa chinensis, 29 in Pyrus bretschneideri, 18 in Rubus occidentalis, 20 in Prunus mume, and 15 in Prunus salicina. Phylogenetic analysis revealed that the Lhcb gene family could be classified into seven major subfamilies, with members of each subfamily sharing similar conserved motifs. And, the functions of each subfamily was predicted based on the previous reports from other species. The Lhcb proteins were highly conserved within their respective subfamilies, suggesting similar functions. Interestingly, we observed similar peaks in Ks values (0.1-0.2) for Lhcb genes in apple and pear, indicating a recent whole genome duplication event (about 30 to 45 million years ago). Additionally, a few Lhcb genes underwent tandem duplication and were located across all chromosomes of nine species of Rosaceae. Furthermore, the analysis of the cis-acting elements in the 2000 bp promoter region upstream of the pear Lhcb gene revealed four main categories: light response correlation, stress response correlation, hormone response correlation, and plant growth. Quantitative expression analysis demonstrated that Lhcb genes exhibited tissue-specific expression patterns and responded differently to low-temperature stress in Rosaceae species.

CONCLUSIONS

These findings shed light on the evolution and phylogeny of Lhcb genes in Rosaceae and highlight the critical role of Lhcb in pear's response to low temperatures. The results obtained provide valuable insights for further investigations into the functions of Lhcb genes in Rosaceae, and these functional genes will be used for further fruit tree breeding and improvement to cope with the current climate changes.

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses

Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

Harmonized cross-species cell atlases of trigeminal and dorsal root ganglia

Peripheral sensory neurons in the dorsal root ganglion (DRG) and trigeminal ganglion (TG) are specialized to detect and transduce diverse environmental stimuli including touch, temperature, and pain to the central nervous system. Recent advances in single-cell RNA-sequencing (scRNA-seq) have provided new insights into the diversity of sensory ganglia cell types in rodents, non-human primates, and humans, but it remains difficult to compare transcriptomically defined cell types across studies and species. Here, we built cross-species harmonized atlases of DRG and TG cell types that describe 18 neuronal and 11 non-neuronal cell types across 6 species and 19 studies. We then demonstrate the utility of this harmonized reference atlas by using it to annotate newly profiled DRG nuclei/cells from both human and the highly regenerative axolotl. We observe that the transcriptomic profiles of sensory neuron subtypes are broadly similar across vertebrates, but the expression of functionally important neuropeptides and channels can vary notably. The new resources and data presented here can guide future studies in comparative transcriptomics, simplify cell type nomenclature differences across studies, and help prioritize targets for future pain therapy development.

Beginner's Guide on the Use of PAML to Detect Positive Selection

The CODEML program in the PAML package has been widely used to analyze protein-coding gene sequences to estimate the synonymous and nonsynonymous rates (dS and dN) and to detect positive Darwinian selection driving protein evolution. For users not familiar with molecular evolutionary analysis, the program is known to have a steep learning curve. Here, we provide a step-by-step protocol to illustrate the commonly used tests available in the program, including the branch models, the site models, and the branch-site models, which can be used to detect positive selection driving adaptive protein evolution affecting particular lineages of the species phylogeny, affecting a subset of amino acid residues in the protein, and affecting a subset of sites along prespecified lineages, respectively. A data set of the myxovirus (Mx) genes from ten mammal and two bird species is used as an example. We discuss a new feature in CODEML that allows users to perform positive selection tests for multiple genes for the same set of taxa, as is common in modern genome-sequencing projects. The PAML package is distributed at https://github.com/abacus-gene/paml under the GNU license, with support provided at its discussion site (https://groups.google.com/g/pamlsoftware). Data files used in this protocol are available at https://github.com/abacus-gene/paml-tutorial.

Evolution of amino acid propensities under stability-mediated epistasis

Site-specific amino acid preferences are influenced by the genetic background of the protein. The preferences for resident amino acids are expected to, on average, increase over time because of replacements at other sites - a nonadaptive phenomenon referred to as the 'evolutionary Stokes shift'. Alternatively, decreases in resident amino acid propensity have recently been viewed as evidence of adaptations to external environmental changes. Using population genetics theory and thermodynamic stability-constraints, we show that nonadaptive evolution can lead to both positive and negative shifts in propensities following the fixation of an amino acid, emphasizing that the detection of negative shifts is not conclusive evidence of adaptation. Considering shifts in propensities over windows between substitutions at a focal site, we find that following ≈ 50% of substitutions the propensity for the new resident amino acid decreases over time, and both positive and negative shifts were comparable in magnitude. Preferences were often conserved via a significant negative autocorrelation in propensity changes-increases in propensities often followed by decreases, and vice versa. Lastly, we explore the underlying mechanisms that lead propensities to fluctuate. We observe that stabilizing replacements increase the mutational tolerance at a site and in doing so decrease the propensity for the resident amino acid. In contrast, destabilizing substitutions result in more rugged fitness landscapes that tend to favor the resident amino acid. In summary, our results characterize propensity trajectories under nonadaptive stability-constrained evolution against which evidence of adaptations should be calibrated.