Signatures of Adaptation in Mitochondrial Genomes of Palearctic Subterranean Voles (Arvicolinae, Rodentia)

This study evaluates signatures of selection in the evolution of the mitochondrial DNA of voles, subfamily Arvicolinae, during the colonization of subterranean environments. The comparative sequence analysis of mitochondrial protein-coding genes of eight subterranean vole species (Prometheomys schaposchnikowi, three species of the genus Ellobius: Ellobius talpinus, Ellobius fuscocapillus and Ellobius lutescens, two species of the genus Terricola: Terricola subterraneus and Terricola daghestanicus, Lasiopodomys mandarinus, and Hyperacrius fertilis) and their closest aboveground relatives was applied using codon-substitution models. The highest number of selection signatures was detected in genes ATP8 and CYTB. The relaxation of selection was observed in most mitochondrial DNA protein-coding genes for subterranean species. The largest amount of relaxed genes is discovered in mole voles (genus Ellobius). The number of selection signatures was found to be independent of the evolutionary age of the lineage but fits the degree of specialization to the subterranean niche. The common trends of selective pressures were observed among the evolutionary ancient and highly specialized subterranean rodent families and phylogenetically young lineages of voles. It suggests that the signatures of adaptation in individual mitochondrial protein-coding genes associated with the colonization of the subterranean niche may appear within a rather short evolutionary timespan.

Signatures of Adaptation in Mitochondrial Genomes of the Palearctic Subterranean Voles (Arvicolinae, Rodentia)

The current study evaluates the selection signals in the evolution of mitochondrial DNA of voles, subfamily Arvicolinae, during the colonization of subterranean environments. The comparative sequence analysis of mitochondrial protein-coding genes of eight subterranean vole species (Prometheomys schaposchnikowi, three species of the genus Ellobius: E. talpinus, E. fuscocapillus and E. lutescens, two species of the genus Terricola: T. subterraneus and T. daghestanicus, Lasiopodomys mandarinus and Hyperacrius fertilis) and their closest aboveground relatives using codon-substitution models was applied. The highest number of selection signatures was detected in genes ATP8 and CYTB. The relaxation of selection was observed in most mtDNA protein-coding genes. In mole voles (genus Ellobius) the signatures of adaptive evolution of mitochondrial genes related to subterranean niche were most pronounced. The number of selection signatures was found to be independent of the evolutionary age of the lineage but fits the degree of specialization to the subterranean niche. The common trends of selective pressures were observed among the evolutionary ancient and highly specialized subterranean rodent families and phylogenetically young lineages of voles. It suggests that the signatures of adaptations in individual mitochondrial protein-coding genes associated with the colonization of the subterranean niche may appear within a rather short evolutionary timespan.

The Adaptive Evolution and Gigantism Mechanisms of the Hadal “Supergiant” Amphipod Alicella gigantea

Hadal trenches are commonly referred to as the deepest areas in the ocean and are characterized by extreme environmental conditions such as high hydrostatic pressures and very limited food supplies. Amphipods are considered the dominant scavengers in the hadal food web. Alicella gigantea is the largest hadal amphipod and, as such, has attracted a lot of attention. However, the adaptive evolution and gigantism mechanisms of the hadal “supergiant” remain unknown. In this study, the whole-body transcriptome analysis was conducted regarding the two hadal amphipods, one being the largest sized species A. gigantea from the New Britain Trench and another the small-sized species Bathycallisoma schellenbergi from the Marceau Trench. The size and weight measurement of the two hadal amphipods revealed that the growth of A. gigantea was comparatively much faster than that of B. schellenbergi. Phylogenetic analyses showed that A. gigantea and B. schellenbergi were clustered into a Lysianassoidea clade, and were distinct from the Gammaroidea consisting of shallow-water Gammarus species. Codon substitution analyses revealed that “response to starvation,” “glycerolipid metabolism,” and “meiosis” pathways were enriched among the positively selected genes (PSGs) of the two hadal amphipods, suggesting that hadal amphipods are subjected to intense food shortage and the pathways are the main adaptation strategies to survive in the hadal environment. To elucidate the mechanisms underlying the gigantism of A. gigantea, small-sized amphipods were used as the background for evolutionary analysis, we found the seven PSGs that were ultimately related to growth and proliferation. In addition, the evolutionary rate of the gene ontology (GO) term “growth regulation” was significantly higher in A. gigantea than in small-sized amphipods. By combining, those points might be the possible gigantism mechanisms of the hadal “supergiant” A. gigantea.

The A756T Mutation of the ERG11 Gene Associated With Resistance to Itraconazole in Candida Krusei Isolated From Mycotic Mastitis of Cows

Candida krusei (C. krusei) has been recently recognized as an important pathogen involved in mycotic mastitis of cows. The phenotypic and molecular characteristics of 15 C. krusei clinical isolates collected from cows with clinical mastitis in three herds of Yinchuan, Ningxia, were identified by matrix-assisted laser desorption ionization–time of flight mass spectrometry. In addition to sequencing analysis, the ERG11 gene that encodes 14α-demethylases, the expression of the ERG11 gene, and efflux transporters ABC1 and ABC2 in itraconazole-susceptible (S), itraconazole-susceptible dose dependent (SDD), and itraconazole-resistant (R) C. krusei isolates was also quantified by a quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) assay. Sequencing analysis revealed three synonymous codon substitutions of the ERG11 gene including T939C, A756T, and T642C in these C. krusei clinical isolates. Among them, T642C and T939C mutations were detected in itraconazole-resistant and -susceptible C. krusei isolates, but the A756T substitution was found only in itraconazole-resistant isolates. Importantly, the expression of the ERG11 gene in itraconazole-resistant isolates was significantly higher compared with itraconazole-SDD and itraconazole-susceptible isolates (p = 0.052 and p = 0.012, respectively), as determined by the qRT-PCR assay. Interestingly, the expression of the ABC2 gene was also significantly higher in itraconazole-resistant isolates relative to the itraconazole-SDD and itraconazole-susceptible strains. Notably, the expression of ERG11 was positively associated with resistance to itraconazole (p = 0.4177 in SDD compared with S, p = 0.0107 in SDD with R, and p = 0.0035 in S with R, respectively). These data demonstrated that mutations of the ERG11 gene were involved in drug resistance in C. krusei. The A756T synonymous codon substitution of the ERG11 gene was correlated with an increased expression of drug-resistant genes including ERG11 and ABC2 in itraconazole-resistant C. krusei isolates examined in this study.

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

We use first principles of population genetics to model the evolution of proteins under persistent positive selection (PPS). PPS may occur when organisms are subjected to persistent environmental change, during adaptive radiations, or in host-pathogen interactions. Our mutation-selection model indicates protein evolution under PPS is an irreversible Markov process, and thus proteins under PPS show a strongly asymmetrical distribution of selection coefficients among amino acid substitutions. Our model shows the criteria ω > 1 (where ω is the ratio of non-synonymous over synonymous codon substitution rates) to detect positive selection is conservative and indeed arbitrary, because in real proteins many mutations are highly deleterious and are removed by selection even at positively-selected sites. We use a penalized-likelihood implementation of our model to successfully detect PPS in plant RuBisCO and influenza HA proteins. By directly estimating selection coefficients at protein sites, our inference procedure bypasses the need for using ω as a surrogate measure of selection and improves our ability to detect molecular adaptation in proteins.

Extra base hits: Widespread empirical support for instantaneous multiple-nucleotide changes

Despite many attempts to introduce evolutionary models that permit substitutions to instantly alter more than one nucleotide in a codon, the prevailing wisdom remains that such changes are rare and generally negligible or are reflective of non-biological artifacts, such as alignment errors. Codon models continue to posit that only single nucleotide change have non-zero rates. Here, we develop and test a simple hierarchy of codon-substitution models with non-zero evolutionary rates for only one-nucleotide (1H), one- and two-nucleotide (2H), or any (3H) codon substitutions. Using over 42, 000 empirical alignments, we find widespread statistical support for multiple hits: 61% of alignments prefer models with 2H allowed, and 23%—with 3H allowed. Analyses of simulated data suggest that these results are not likely to be due to simple artifacts such as model misspecification or alignment errors. Further modeling reveals that synonymous codon island jumping among codons encoding serine, especially along short branches, contributes significantly to this 3H signal. While serine codons were prominently involved in multiple-hit substitutions, there were other common exchanges contributing to better model fit. It appears that a small subset of sites in most alignments have unusual evolutionary dynamics not well explained by existing model formalisms, and that commonly estimated quantities, such as dN/dS ratios may be biased by model misspecification. Our findings highlight the need for continued evaluation of assumptions underlying workhorse evolutionary models and subsequent evolutionary inference techniques. We provide a software implementation for evolutionary biologists to assess the potential impact of extra base hits in their data in the HyPhy package and in the Datamonkey.org server.

Holocentric Chromosomes Probably Do Not Prevent Centromere Drive in Cyperaceae

Centromere drive model describes an evolutionary process initiated by centromeric repeats expansion, which leads to the recruitment of excess kinetochore proteins and consequent preferential segregation of an expanded centromere to the egg during female asymmetric meiosis. In response to these selfish centromeres, the histone protein CenH3, which recruits kinetochore components, adaptively evolves to restore chromosomal parity and counter the detrimental effects of centromere drive. Holocentric chromosomes, whose kinetochores are assembled along entire chromosomes, have been hypothesized to prevent expanded centromeres from acquiring a selective advantage and initiating centromere drive. In such a case, CenH3 would be subjected to less frequent or no adaptive evolution. Using codon substitution models, we analyzed 36 CenH3 sequences from 35 species of the holocentric family Cyperaceae. We found 10 positively selected codons in the CenH3 gene [six codons in the N-terminus and four in the histone fold domain (HFD)] and six branches of its phylogeny along which the positive selection occurred. One of the positively selected codons was found in the centromere targeting domain (CATD) that directly interacts with DNA and its mutations may be important in centromere drive suppression. The frequency of these positive selection events was comparable to the frequency of positive selection in monocentric clades with asymmetric female meiosis. Taken together, these results suggest that preventing centromere drive is not the primary adaptive role of holocentric chromosomes, and their ability to suppress it likely depends on their kinetochore structure in meiosis.

Soluble Expression of hFGF19 without Fusion Protein through Synonymous Codon Substitutions and DsbC Co-Expression in E. coli

Human fibroblast growth factor 19 (hFGF19) is a difficult-to-express protein that is frequently fused with another protein for soluble expression. However, residual amino acids after cleavage with protease represent one of the major problems in therapeutic protein development. Here, we introduced synonymous codon substitutions in the N-terminal region encoding sequence of hFGF19 and co-expressed disulfide bond isomerase (ΔssDsbC) to functionally express hFGF19 without any fusion protein. Synonymous codon substitution significantly increased hFGF19 expression. Subsequent co-expression of ΔssDsbC with a selected variant of hFGF19 (scvhFGF19) further increased the proportion of soluble hFGF19 expression in Escherichia coli XL1-Blue. Both total and soluble scvhFGF19 expression increased remarkably in the alternative host, E. coli Origami 2 with mutated thioredoxin reductase and glutathione reductase. scvhFGF19 purification by anion exchange and heparin affinity chromatography resulted in a yield of 6.5 mg/L under normal induction conditions in flask culture. As such, a high cell density culture is expected to achieve an even higher yield. The biological activities of purified scvhFGF19 were assessed based on its ability to activate ERK1/2 signaling pathway in HepG2 hepatocarcinoma cells. In conclusion, the strategy described here may represent an efficient alternative process for the production of hFGF19 and/or related proteins.

A Bayesian Mutation–Selection Framework for Detecting Site-Specific Adaptive Evolution in Protein-Coding Genes

In recent years, codon substitution models based on the mutation–selection principle have been extended for the purpose of detecting signatures of adaptive evolution in protein-coding genes. However, the approaches used to date have either focused on detecting global signals of adaptive regimes—across the entire gene—or on contexts where experimentally derived, site-specific amino acid fitness profiles are available. Here, we present a Bayesian site-heterogeneous mutation–selection framework for site-specific detection of adaptive substitution regimes given a protein-coding DNA alignment. We offer implementations, briefly present simulation results, and apply the approach on a few real data sets. Our analyses suggest that the new approach shows greater sensitivity than traditional methods. However, more study is required to assess the impact of potential model violations on the method, and gain a greater empirical sense its behavior on a broader range of real data sets. We propose an outline of such a research program.

Consequences of Stability-Induced Epistasis for Substitution Rates

Do interactions between residues in a protein (i.e., epistasis) significantly alter evolutionary dynamics? If so, what consequences might they have on inference from traditional codon substitution models which assume site-independence for the sake of computational tractability? To investigate the effects of epistasis on substitution rates, we employed a mechanistic mutation-selection model in conjunction with a fitness framework derived from protein stability. We refer to this as the stability-informed site-dependent (S-SD) model and developed a new stability-informed site-independent (S-SI) model that captures the average effect of stability constraints on individual sites of a protein. Comparison of S-SI and S-SD offers a novel and direct method for investigating the consequences of stability-induced epistasis on protein evolution. We developed S-SI and S-SD models for three natural proteins and showed that they generate sequences consistent with real alignments. Our analyses revealed that epistasis tends to increase substitution rates compared with the rates under site-independent evolution. We then assessed the epistatic sensitivity of individual site and discovered a counterintuitive effect: Highly connected sites were less influenced by epistasis relative to exposed sites. Lastly, we show that, despite the unrealistic assumptions, traditional models perform comparably well in the presence and absence of epistasis and provide reasonable summaries of average selection intensities. We conclude that epistatic models are critical to understanding protein evolutionary dynamics, but epistasis might not be required for reasonable inference of selection pressure when averaging over time and sites.