Extreme variations in the ratios of non-synonymous to synonymous nucleotide substitution rates in signal peptide evolution.

Nucleotide sequences encoding signal peptides from the precursors of alpha-amylase/trypsin inhibitors from cereals are homologous to those corresponding to the precursors of thaumatin II and of plastocyanins. Non-synonymous (KA) and synonymous (KS) rates of nucleotide substitutions have been calculated for all possible binary combinations. Extreme variation in KA/KS ratios has been observed; from the 0.167 average found within the plastocyanin family to an average of 1.90 calculated for the inhibitors/thaumatin II transition. A similar calculation has been carried out for the signal peptide sequences of thionins, which are unrelated to those of the alpha-amylase/trypsin inhibitor family, and an average KA/KS of 0.12 has been obtained. This variation can be largely explained in terms of an empirical index of stability related to amino acid composition and seems to be independent of functional constraints.     

Weak preservation of local neutral substitution rates across mammalian genomes

Background  

The rate at which neutral (non-functional) bases undergo substitution is highly dependent on their location within a genome. However, it is not clear how fast these location-dependent rates change, or to what extent the substitution rate patterns are conserved between lineages. To address this question, which is critical not only for understanding the substitution process but also for evaluating phylogenetic footprinting algorithms, we examine ancestral repeats: a predominantly neutral dataset with a significantly higher genomic density than other datasets commonly used to study substitution rate variation. Using this repeat data, we measure the extent to which orthologous ancestral repeat sequences exhibit similar substitution patterns in separate mammalian lineages, allowing us to ascertain how well local substitution rates have been preserved across species.     

Efficient biased estimation of evolutionary distances when substitution rates vary across sites

This paper deals with phylogenetic inference when the variability of substitution rates across sites (VRAS) is modeled by a gamma distribution. We show that underestimating VRAS, which results in underestimates for the evolutionary distances between sequences, usually improves the topological accuracy of phylogenetic tree inference by distance-based methods, especially when the molecular clock holds. We propose a method to estimate the gamma shape parameter value which is most suited for tree topology inference, given the sequences at hand. This method is based on the pairwise evolutionary distances between sequences and allows one to reconstruct the phylogeny of a high number of taxa (>1,000). Simulation results show that the topological accuracy is highly improved when using the gamma shape parameter value given by our method, compared with the true (unknown) value which was used to generate the data. Furthermore, when VRAS is high, the topological accuracy of our distance-based method is better than that of a maximum likelihood approach. Finally, a data set of Maoricicada species sequences is analyzed, which confirms the advantage of our method.     

A targeted amplicon sequencing panel to simultaneously identify mosquito species and Plasmodium presence across the entire Anopheles genus

Anopheles is a diverse genus of mosquitoes comprising over 500 described species, including all known human malaria vectors. While a limited number of key vector species have been studied in detail, the goal of malaria elimination calls for surveillance of all potential vector species. Here, we develop a multilocus amplicon sequencing approach that targets 62 highly variable loci in the Anopheles genome and two conserved loci in the Plasmodium mitochondrion, simultaneously revealing both the mosquito species and whether that mosquito carries malaria parasites. We also develop a cheap, nondestructive, and high-throughput DNA extraction workflow that provides template DNA from single mosquitoes for the multiplex PCR, which means specimens producing unexpected results can be returned to for morphological examination. Over 1000 individual mosquitoes can be sequenced in a single MiSeq run, and we demonstrate the panel’s power to assign species identity using sequencing data for 40 species from Africa, Southeast Asia, and South America. We also show that the approach can be used to resolve geographic population structure within An. gambiae and An. coluzzii populations, as the population structure determined based on these 62 loci from over 1000 mosquitoes closely mirrors that revealed through whole genome sequencing. The end-to-end approach is quick, inexpensive, robust, and accurate, which makes it a promising technique for very large-scale mosquito genetic surveillance and vector control.     

Partitioning the variation in mammalian substitution rates

We have used analysis of variance to partition the variation in synonymous and amino acid substitution rates between three effects (gene, lineage, and a gene-by-lineage interaction) in mammalian nuclear and mitochondrial genes. We find that gene effects are stronger for amino acid substitution rates than for synonymous substitution rates and that lineage effects are stronger for synonymous substitution rates than for amino acid substitution rates. Gene-by-lineage interactions, equivalent to overdispersion corrected for lineage effects, are found in amino acid substitutions but not in synonymous substitutions. The variance in the ratio of amino acid and synonymous substitution rates is dominated by gene effects, but there is also a significant gene-by-lineage interaction.     

Age-specific metabolic rates and mortality rates in the genus Drosophila

Early theories of aging suggested that organisms with relatively high metabolic rates would live shorter lives. Despite widespread tests of this 'rate of living' theory of aging, there is little empirical evidence to support the idea. A more fine-grained approach that examined age-related changes in metabolic rate over the life span could provide valuable insight into the relationship between metabolic rate and aging. Here we compare age-related metabolic rate (measured as CO2 production per hour) and age-related mortality rate among five species in the genus Drosophila. We find no evidence that longer-lived species have lower metabolic rates. In all five species, there is no clear evidence of an age-related metabolic decline. Metabolic rates are strikingly constant throughout the life course, with the exception of females of D. hydei, in which metabolic rates show an increase over the first third of the life span and then decline. We argue that some physiological traits may have been shaped by such strong selection over evolutionary time that they are relatively resistant to the decline in the force of selection that occurs within the life time of a single individual. We suggest that comparisons of specific traits that do not show signs of aging with those traits that do decline with age could provide insight into the aging process.     

Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene

Background

Synonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.

Methodology/Principle Findings

We found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.

Conclusions/Significance

We found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.     

Site-to-site variation of synonymous substitution rates

We develop a new model for studying the molecular evolution of protein-coding DNA sequences. In contrast to existing models, we incorporate the potential for site-to-site heterogeneity of both synonymous and nonsynonymous substitution rates. We demonstrate that within-gene heterogeneity of synonymous substitution rates appears to be common. Using the new family of models, we investigate the utility of a variety of new statistical inference procedures, and we pay particular attention to issues surrounding the detection of sites undergoing positive selection. We discuss how failure to model synonymous rate variation in the model can lead to misidentification of sites as positively selected.     

Conservation of physiological dysregulation signatures of aging across primates

Two major goals in the current biology of aging are to identify general mechanisms underlying the aging process and to explain species differences in aging. Recent research in humans suggests that one important driver of aging is dysregulation, the progressive loss of homeostasis in complex biological networks. Yet, there is a lack of comparative data for this hypothesis, and we do not know whether dysregulation is widely associated with aging or how well signals of homeostasis are conserved. To address this knowledge gap, we use unusually detailed longitudinal biomarker data from 10 species of nonhuman primates housed in research centers and data from two human populations to test the hypotheses that (a) greater dysregulation is associated with aging across primates and (b) physiological states characterizing homeostasis are conserved across primates to degrees associated with phylogenetic proximity. To evaluate dysregulation, we employed a multivariate distance measure, calculated from sets of biomarkers, that is associated with aging and mortality in human populations. Dysregulation scores positively correlated with age and risk of mortality in most nonhuman primates studied, and signals of homeostatic state were significantly conserved across species, declining with phylogenetic distance. Our study provides the first broad demonstration of physiological dysregulation associated with aging and mortality risk in multiple nonhuman primates. Our results also imply that emergent signals of homeostasis are evolutionarily conserved, although with notable variation among species, and suggest promising directions for future comparative studies on dysregulation and the aging process.     

Immune signatures in cutaneous melanoma correlate with survival independently of immunotherapy treatment

Immune checkpoint inhibitors (ICIs) have fundamentally improved survival from advanced cutaneous melanoma. Significant efforts have been made to understand the ICI response to identify ways to further improve outcomes. One such approach has been to investigate gene expression associated with response to ICI, which has identified various immune-related mRNA signatures, including a six-gene IFN-γ signature (IFN-γ⁶), an expanded immune signature (IFN-γ¹⁸), an effector T-cell gene signature (T_eff), and a T_eff-associated and IFN-γ-associated gene signature (T_eff + IFN-γ). Given that these signatures appear to reflect expression from T cells and the level of tumour-infiltrating immune cells has been associated with survival, we hypothesised that the prognostic value of the signatures is not limited to ICI treatment and investigated if they were associated with survival also in patients who never received ICI. The signatures were not present in melanoma cell lines when compared with tumour samples, confirming that the signatures were likely derived from the samples' non-tumour (immune) components. We acquired expression and survival data from five melanoma cohorts with a wide range of disease stages, treatments and metrics for survival, and correlated the expression signatures with survival. All four signatures were significantly associated (p < .05) with survival in four of five cohorts, with hazard ratios ranging from 0.69 to 0.92. We conclude that these immune signatures' association with survival is not specific to ICI-treated patients, but present in a number of settings.     

Meta-analysis of diazotrophic signatures across terrestrial ecosystems at the continental scale

Biological nitrogen fixation performed by diazotrophs forms a cornerstone of Earth's terrestrial ecosystem productivity. However, the composition, diversity and distribution of soil diazotrophs are poorly understood across different soil ecosystems. Furthermore, the biological potential of the key diazotroph species in relation to key environmental parameters is unknown. To address this, we used meta-analysis approach to merge together 39 independent diazotroph amplicon sequencing (nifH gene) datasets consisting of 1988 independent soil samples. We then employed multiple statistical analyses and machine-learning approaches to compare diazotroph community differences and indicator species between terrestrial ecosystems on a global scale. The distribution, composition and structure of diazotroph communities varied across seven different terrestrial ecosystems, with community composition exhibiting an especially clear effect. The Cyanobacteria were the most abundant taxa in crust ecosystems (accounting for ~45% of diazotrophs), while other terrestrial ecosystems were dominated by Proteobacteria, including Alpha-, Beta- and Gamma-Proteobacteria (accounting for ~70% of diazotrophs). Farmland ecosystems harboured the highest and crust ecosystems the lowest alpha and phylogenetic diversities. Azospirillum zeae, Skermanella aerolata and four Bradyrhizobium species were identified as key indicator species of potential diazotroph activity. Overall, diazotroph abundances and distribution were affected by multiple environmental parameters, including soil pH, nitrogen, organic carbon, C:N ratio and annual mean precipitation and temperature. Together, our findings suggest that based on the relative abundance and diversity of nifH marker gene, diazotrophs have adapted to a range of environmental niches globally.     

A non-synonymous nucleotide substitution can account for one evolutionary route to sesquiterpene synthase activity in the TPS-b subgroup

Plant sesquiterpene and hemiterpene synthases in the monoterpene synthase dominated TPS-b subgroup are thought to have evolved independently from a monoterpene synthase ancestor. A TPS-b sesquiterpene synthase from apple (MdAFS1), which predominantly produces α-farnesene, can also synthesize the monoterpene (E)-β-ocimene. The dual activity offered a functional link to an ancestral MdAFS1 enzyme and a rational basis for investigation of the evolution of TPS-b sesquiterpene enzymes. Protein modelling and mutagenesis analysis of the MdAFS1 active site identified a non-synonymous nucleotide substitution that could account for the requisite shift in substrate specificity necessary for the emergence of its sesquiterpene activity during the evolution of the TPS-b enzymes.     

In silico investigation of the ATP7B gene: insights from functional prediction of non-synonymous substitution to protein structure

ATP7B is a copper-transporting ATPase that plays a key role in the regulation of copper homeostasis. Mutations in the ATP7B gene are causative for Wilson’s disease, and recent reports have suggested that genetic variants are associated with susceptibility to Alzheimer’s disease. Unfortunately, it is difficult to profile experimentally novel genetic variants in the ATP7B gene, because the human protein X-ray structure is not yet entirely understood. In order to investigate ATP7B non-synonymous substitutions, we used an in silico amino acid sequence-based approach. Specifically, we analyzed 337 ATP7B non-synonymous substitutions, which included Wilson’s disease-causing mutations (DVs) and non Wilson’s disease-causing variants (NDVs), with an algorithm that estimated a combined probability (cP_del) of an amino acidic change to be deleterious for the protein function. This approach appeared to reliably indentify the probability of DVs and NDVs to be deleterious and to profile still unknown gene variants. Specifically, after analyzing ATP7B protein domains with the cP_del method, we found results in line with the predicted–modeled domains and some new suggestions. In conclusion, a functional survey of amino acid changes in the ATP7B protein is provided herein, and we suggest that this bioinformatic method can furnish information about novel ATP7B mutations. Furthermore, the same approach can be applied to other uncharacterized proteins.     

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

Summary The current convention in estimating the number of substitutions per synonymous site (K _S ) and per nonsynonymous site (K _A ) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theK _S value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aK _S value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aK _A value close to the number of substitutions per nondegenerate site.     

Approximate methods for estimating the pattern of nucleotide substitution and the variation of substitution rates among sites

We propose two approximate methods (one based on parsimony and one on pairwise sequence comparison) for estimating the pattern of nucleotide substitution and a parsimony-based method for estimating the gamma parameter for variable substitution rates among sites. The matrix of substitution rates that represents the substitution pattern can be recovered through its relationship with the observable matrix of site pattern frequences in pairwise sequence comparisons. In the parsimony approach, the ancestral sequences reconstructed by the parsimony algorithm were used, and the two sequences compared are those at the ends of a branch in the phylogenetic tree. The method for estimating the gamma parameter was based on a reinterpretation of the numbers of changes at sites inferred by parsimony. Three data sets were analyzed to examine the utility of the approximate methods compared with the more reliable likelihood methods. The new methods for estimating the substitution pattern were found to produce estimates quite similar to those obtained from the likelihood analyses. The new method for estimating the gamma parameter was effective in reducing the bias in conventional parsimony estimates, although it also overestimated the parameter. The approximate methods are computationally very fast and appear useful for analyzing large data sets, for which use of the likelihood method requires excessive computation.      

Calculating independent contrasts for the comparative study of substitution rates

Phylogenetic comparative methods have been used to study patterns of correlated evolution between biological traits of all kinds, and are increasingly used to identify predictors of the rate of DNA substitution. Substitution rate differs from most traits studied in that it cannot be observed directly, but must be inferred from substitutions accrued over a long period of time. Studying a mathematical model of trait and rate evolution, we show that the special nature of substitution rates can lead to a severe loss of power for standard comparative methods. We further show how strategies designed to maximise power, by increasing the number of data points, can have the opposite effect when substitution rate is involved. We then propose two modifications of the standard methods that can mitigate these problems. First, we show how pre-existing tests for homogeneity of variance can be used to identify and exclude those data from which changes in substitution rate cannot be reliably inferred. Second, we show that power can be increased by comparing substitution rate with the time average of the predictor trait along the history of the lineage, and introduce a maximum likelihood estimator of this quantity.