CodonTest: modeling amino acid substitution preferences in coding sequences

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.

Genome wide and evolutionary analysis of heat shock protein 70 proteins in tomato and their role in response to heat and drought stress

Heat shock protein 70 (HSP70) proteins play a crucial role in mitigating the detrimental effects of abiotic stresses in plants. In the present study, 21 full length non-redundant SlHSP70 genes were detected and characterized in tomato (Solanum lycopersicum L.). The SlHSP70 genes were classified into four groups based on phylogenetic analysis. Similarities were observed in gene features and motif structures of SlHSP70s belonging to the same group. SlHSP70 genes were unevenly and unequally mapped on 11 chromosomes. Segmental and tandem duplication are the main events that have contributed to the expansion of the SlHSP70 genes. A large number of groups and sub-groups were generated during comparative analysis of HSP70 genes in multiple plant species including tomato. These findings indicated a common ancestor which created diverse sub-groups prior to a mono-dicot split. The selection pressure on specific codons was identified through a maximum-likelihood approach and we found some important coding sites in the coding region of all groups. Diversifying positive selection was indirectly associated with evolutionary changes in SlHSP70 proteins and suggests that gene evolution modulated the tomato domestication event. In addition, expression analysis using RNA-seq revealed that 21 SlHSP70 genes were differentially expressed in response to drought and heat stress. SlHSP70-5 was down-regulated by heat treatment and up-regulated by drought stress. Furthermore, the expression of some of the duplicate genes was partially redundant, while others showed functional diversity. Our results indicate the diverse role of HSP70 gene family in S. lycopersicum under drought and heat stress conditions and open the gate for further investigation of HSP70 gene family functions, especially under drought and heat stress.

Identification, evolutionary analysis and functional diversification of RAV gene family in cotton (G. hirsutum L.)

Genome wide analysis, expression pattern analysis, and functional characterization of RAV genes highlight their roles in roots, stem development and hormonal response. RAV (Related to ABI3 and VP1) gene family members have been involved in tissues/organs growth and hormone signaling in various plant species. Here, we identified 247 RAVs from 12 different species with 33 RAV genes from G. hirsutum. Phylogenetic analysis classified RAV genes into four distinct groups. Analysis of gene structure showed that most GhRAVs lack introns. Motif distribution pattern and protein sequence logos indicated that GhRAV genes were highly conserved during the process of evolution. Promotor cis-acting elements revealed that promotor regions of GhRAV genes encode numerous elements related to plant growth, abiotic stresses and phytohormones. Chromosomal location information showed uneven distribution of 33 GhRAV genes on different chromosomes. Collinearity analysis identified 628 and 52 orthologous/ paralogous gene pairs in G. hirsutum and G. barbadense, respectively. Ka/Ks values indicated that GhRAV and GbRAV genes underwent strong purifying selection pressure. Selecton model and codon model selection revealed that GhRAV amino acids were under purifying selection and adaptive evolution exists among GhRAV proteins. Three dimensional structure of GhRAVs indicated the presence of numerous alpha helix and beta-barrels. Expression level revealed that some GhRAV genes exhibited high expression in roots (GhRAV3, GhRAV4, GhRAV11, GhRAV18, GhRAV20 and GhRAV30) and stem (GhRAV3 and GhRAV18), indicating their potential role in roots and stem development. GhRAV genes can be regulated by phytohormonal stresses (BL, JA and IAA). Our study provides a reference for future studies related to the functional analysis of GhRAVs in cotton.

Weak selection on synonymous codons substantially inflates dN/dS estimates in bacteria

Synonymous codon substitutions are not always selectively neutral as revealed by several types of analyses, including studies of codon usage patterns among genes. We analyzed codon usage in 13 bacterial genomes sampled from across a large order of bacteria, Enterobacterales, and identified presumptively neutral and selected classes of synonymous substitutions. To estimate substitution rates, given a neutral/selected classification of synonymous substitutions, we developed a flexible [Formula: see text] substitution model that allows multiple classes of synonymous substitutions. Under this multiclass synonymous substitution (MSS) model, the denominator of [Formula: see text] includes only the strictly neutral class of synonymous substitutions. On average, the value of [Formula: see text] under the MSS model was 80% of that under the standard codon model in which all synonymous substitutions are assumed to be neutral. The indication is that conventional [Formula: see text] analyses overestimate these values and thus overestimate the frequency of positive diversifying selection and underestimate the strength of purifying selection. To quantify the strength of selection necessary to explain this reduction, we developed a model of selected compensatory codon substitutions. The reduction in synonymous substitution rate, and thus the contribution that selection makes to codon bias variation among genes, can be adequately explained by very weak selection, with a mean product of population size and selection coefficient, [Formula: see text].

Contrast-FEL-A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches

A number of evolutionary hypotheses can be tested by comparing selective pressures among sets of branches in a phylogenetic tree. When the question of interest is to identify specific sites within genes that may be evolving differently, a common approach is to perform separate analyses on subsets of sequences and compare parameter estimates in a post hoc fashion. This approach is statistically suboptimal and not always applicable. Here, we develop a simple extension of a popular fixed effects likelihood method in the context of codon-based evolutionary phylogenetic maximum likelihood testing, Contrast-FEL. It is suitable for identifying individual alignment sites where any among the K≥2 sets of branches in a phylogenetic tree have detectably different ω ratios, indicative of different selective regimes. Using extensive simulations, we show that Contrast-FEL delivers good power, exceeding 90% for sufficiently large differences, while maintaining tight control over false positive rates, when the model is correctly specified. We conclude by applying Contrast-FEL to data from five previously published studies spanning a diverse range of organisms and focusing on different evolutionary questions.

Adaptive evolution of peptidoglycan recognition protein family regulates the innate signaling against microbial pathogens in vertebrates

The innate immune system is the first line of defense in vertebrates against microbial pathogens. This defense system depends on the peptidoglycan pathogen recognition of receptors (PGRPs) existing in both invertebrates and vertebrates. Although some studies revealed the structural and functional differences between them, however, the evolutionary history and the selection pressures on these genes during adaptive evolution are poorly understood. In this study, we examined four (PGLYRP1, PGLYRP2, PGLYRP3, and PGLYRP4) genes of 127 vertebrates' species, conserved across vertebrates to evaluate positive selection pressure drives by adaptive evolution. The codons under positive selection were recognized through likelihood tests by comparing different models based on ω ratios in these genes across the vertebrate species. The positive selection test used two sets of models M1a vs. M2a and M7 vs. M8. The results showed that the test of these genes in M1a vs. M2a was not significant with the likelihood value 2ΔlnL = 0, while the likelihood ratios (2ΔlnL) were 2ΔlnL = 12.386, 2ΔlnL = 4.9283, 2ΔlnL = 24.031, and 2ΔlnL = 103.39 for PGLYRP1, PGLYRP2, PGLYRP3, and PGLYRP4 in M7 vs. M8, respectively. Our study identified the evidence of robust positive selection for these four genes across the vertebrates. These protuberant changes in PGRPs evolution of vertebrates reveal their role in innate immunity. Our study provides an insight based on PGRP genes to understand the evolution of host and pathogens interaction that leads to the progress of the novel conducts for immune diseases that include proteins linked to the recognition of pathogens.

Adaptive Molecular Evolution of AKT3 Gene for Positive Diversifying Selection in Mammals

The V-Akt Murine Thymoma Viral Oncogene Homolog 3 (AKT3) gene is of the serine/threonine-protein kinase family and influences the production of milk fats and cholesterol by acting on the sterol administrative area restricting protein (SREBP). The AKT3 gene is highly preserved in animals, and during lactation in cattle, its expression increases. The AKT3 gene is expressed in the digestive system, mammary gland, and immune cells. A phylogenetic investigation was performed to clarify the evolutionary role of AKT3, by maximum probability. The AKT3 gene sequence data of various mammalian species was evident even with animals undergoing breeding selection. From 39 mammalian species studied, there was a signal of positive diversifying selection with Hominidae at 13Q, 16G, 23R, 24P, 121P, 294K, 327V, 376L, 397K, 445T, and 471F among other codon sites of the AKT3 gene. These sites were codes for amino acids such as arginine, proline, lysine, and leucine indicating major roles for the function of immunological proteins, and in particular, the study highlighted the importance of changes in gene expression of AKT3 on immunity.

Adaptation to host-specific bacterial pathogens drive rapid evolution of novel PhoP/PhoQ regulation pathway modulating the virulence

The presence of the PhoP-PhoQ system is usually different in various bacterial groups, suggesting that PhoP can control the expression of different genes in species. However, little is known about the evolution of the PhoP-PhoQ system among bacterial pathogens. Here, we study the evolution of PhoP and PhoQ regulation in 15 species of Enterobacteriaceae family. We have determined that the regulatory objectives adopted by PhoP and PhoQ are mainly different, due to the result of horizontal gene transfer events and even the change in the genetic content between closely related species. We have compared many possibilities tests (M1 vs. M2 and M7 with M8) to determine the positive selection. Estimating parameters at M1 and M2, with positive selection in M2 of the two proteins. The proportions of positive selection sites significant with ω = 4.53076 for PhoP and ω = 4.21041 PhQ. M8 was significant for PhoP and PhQ proteins. To further confirm the positive selection results, we used the Selecton server to confer positive selection on individual sites using the Mechanistic-Empirical Combination model, and we noticed that several sites had been identified under selection pressure during the evolution. There was a strong indication for the positive selection in bacterial genes of PhoP and PhoQ showed the results. By the use of REL and IFEL, the positive selection for PhoP was detected 14 and 11 sites respectively at different codon positions. The positively selected sites of amino acids such as Arginine, Alanine, Lysine, and Leucine are more important for the production of signals. Our results suggest that the positive selection of PhoP-PhoQ genes in host adaptation during evolution raises an intriguing possibility causes subtle variations in actions of PhoP-PhoQ and also increases the opportunities that cause modification in protein structure for the evolution of increasing pathogenicity in bacterial pathogens.

Non-neutral evolution of H3.3-encoding genes occurs without alterations in protein sequence

Histone H3.3 is a developmentally essential variant encoded by two independent genes in human (H3F3A and H3F3B). While this two-gene arrangement is evolutionarily conserved, its origins and function remain unknown. Phylogenetics, synteny and gene structure analyses of H3.3 genes from 32 metazoan genomes indicate independent evolutionary paths for H3F3A and H3F3B. While H3F3B bears similarities with H3.3 genes in distant organisms and with canonical H3 genes, H3F3A is sarcopterygian-specific and evolves under strong purifying selection. Additionally, H3F3B codon-usage preferences resemble those of broadly expressed genes and 'cell differentiation-induced' genes, while codon-usage of H3F3A resembles that of 'cell proliferation-induced' genes. We infer that H3F3B is more similar to the ancestral H3.3 gene and likely evolutionarily adapted for a broad expression pattern in diverse cellular programs, while H3F3A adapted for a subset of gene expression programs. Thus, the arrangement of two independent H3.3 genes facilitates fine-tuning of H3.3 expression across cellular programs.

Identification of a Novel Gammaherpesvirus in Canada lynx (Lynx canadensis)

Gammaherpesviruses (GHVs) infect many animal species and are associated with lymphoproliferative disorders in some. Previously, we identified several novel GHVs in North American felids; however, a GHV had never been identified in Canada lynx (Lynx canadensis). We, therefore, hypothesized the existence of an unidentified GHV in lynx. Using degenerate nested and subsequently virus-specific PCR, we amplified and sequenced 3.4 kb of DNA from a novel GHV in lynx, which we named Lynx canadensis gammaherpesvirus 1 (LcaGHV1). Phylogenetic analysis determined that LcaGHV1 is a distinct GHV species belonging to the genus Percavirus. We then estimated the prevalence of LcaGHV1 in lynx by developing a PCR-based assay and detected LcaGHV1 DNA in 36% (95% CI: 22-53%) of lynx spleen DNA samples from Maine, USA and 17% (95% CI: 8-31%) from Newfoundland, Canada. The LcaGHV1 DNA sequences from Maine and Newfoundland lynx were nearly identical to each other (two nucleotide substitutions in 3.4 kb), suggesting that the unique lynx subspecies present on the island of Newfoundland (Lynx canadensis subsolanus) is infected with virus that very closely resembles virus found in mainland lynx. The potential ecologic and pathologic consequences of this novel virus for Canada lynx populations warrant further study.

Improved inference of site-specific positive selection under a generalized parametric codon model when there are multinucleotide mutations and multiple nonsynonymous rates

BACKGROUND

An excess of nonsynonymous substitutions, over neutrality, is considered evidence of positive Darwinian selection. Inference for proteins often relies on estimation of the nonsynonymous to synonymous ratio (ω = dN/dS) within a codon model. However, to ease computational difficulties, ω is typically estimated assuming an idealized substitution process where (i) all nonsynonymous substitutions have the same rate (regardless of impact on organism fitness) and (ii) instantaneous double and triple (DT) nucleotide mutations have zero probability (despite evidence that they can occur). It follows that estimates of ω represent an imperfect summary of the intensity of selection, and that tests based on the ω > 1 threshold could be negatively impacted.

RESULTS

We developed a general-purpose parametric (GPP) modelling framework for codons. This novel approach allows specification of all possible instantaneous codon substitutions, including multiple nonsynonymous rates (MNRs) and instantaneous DT nucleotide changes. Existing codon models are specified as special cases of the GPP model. We use GPP models to implement likelihood ratio tests for ω > 1 that accommodate MNRs and DT mutations. Through both simulation and real data analysis, we find that failure to model MNRs and DT mutations reduces power in some cases and inflates false positives in others. False positives under traditional M2a and M8 models were very sensitive to DT changes. This was exacerbated by the choice of frequency parameterization (GY vs. MG), with rates sometimes > 90% under MG. By including MNRs and DT mutations, accuracy and power was greatly improved under the GPP framework. However, we also find that over-parameterized models can perform less well, and this can contribute to degraded performance of LRTs.

CONCLUSIONS

We suggest GPP models should be used alongside traditional codon models. Further, all codon models should be deployed within an experimental design that includes (i) assessing robustness to model assumptions, and (ii) investigation of non-standard behaviour of MLEs. As the goal of every analysis is to avoid false conclusions, more work is needed on model selection methods that consider both the increase in fit engendered by a model parameter and the degree to which that parameter is affected by un-modelled evolutionary processes.

Evolutionary Analysis Provides Insight Into the Origin and Adaptation of HCV

Hepatitis C virus (HCV) belongs to the Hepacivirus genus and is genetically heterogeneous, with seven major genotypes further divided into several recognized subtypes. HCV origin was previously dated in a range between ∼200 and 1000 years ago. Hepaciviruses have been identified in several domestic and wild mammals, the largest viral diversity being observed in bats and rodents. The closest relatives of HCV were found in horses/donkeys (equine hepaciviruses, EHV). However, the origin of HCV as a human pathogen is still an unsolved puzzle. Using a selection-informed evolutionary model, we show that the common ancestor of extant HCV genotypes existed at least 3000 years ago (CI: 3192–5221 years ago), with the oldest genotypes being endemic to Asia. EHV originated around 1100 CE (CI: 291–1640 CE). These time estimates exclude that EHV transmission was mainly sustained by widespread veterinary practices and suggest that HCV originated from a single zoonotic event with subsequent diversification in human populations. We also describe a number of biologically important sites in the major HCV genotypes that have been positively selected and indicate that drug resistance-associated variants are significantly enriched at positively selected sites. HCV exploits several cell-surface molecules for cell entry, but only two of these (CD81 and OCLN) determine the species-specificity of infection. Herein evolutionary analyses do not support a long-standing association between primates and hepaciviruses, and signals of positive selection at CD81 were only observed in Chiroptera. No evidence of selection was detected for OCLN in any mammalian order. These results shed light on the origin of HCV and provide a catalog of candidate genetic modulators of HCV phenotypic diversity.

Hyperendemic dengue transmission and identification of a locally evolved DENV-3 lineage, Papua New Guinea 2007-2010

Background

Dengue is endemic in the Western Pacific and Oceania and the region reports more than 200,000 cases annually. Outbreaks of dengue and severe dengue occur regularly and movement of virus throughout the region has been reported. Disease surveillance systems, however, in many areas are not fully established and dengue incidence is underreported. Dengue epidemiology is likely least understood in Papua New Guinea (PNG), where the prototype DENV-2 strain New Guinea C was first isolated by Sabin in 1944 but where routine surveillance is not undertaken and little incidence and prevalence data is available.

Methodology/Principal findings

Serum samples from individuals with recent acute febrile illness or with non-febrile conditions collected between 2007–2010 were tested for anti-DENV neutralizing antibody. Responses were predominantly multitypic and seroprevalence increased with age, a pattern indicative of endemic dengue. DENV-1, DENV-2 and DENV-3 genomes were detected by RT-PCR within a nine-month period and in several instances, two serotypes were identified in individuals sampled within a period of 10 days. Phylogenetic analysis of whole genome sequences identified a DENV-3 Genotype 1 lineage which had evolved on the northern coast of PNG which was likely exported to the western Pacific five years later, in addition to a DENV-2 Cosmopolitan Genotype lineage which had previously circulated in the region.

Conclusions/Significance

We show that dengue is hyperendemic in PNG and identify an endemic, locally evolved lineage of DENV-3 that was associated with an outbreak of severe dengue in Pacific countries in subsequent years, although severe disease was not identified in PNG. Additional studies need to be undertaken to understand dengue epidemiology and burden of disease in PNG.

Dengue virus (DENV) was first identified in Papua New Guinea (PNG) in 1944. Dengue is currently assumed to be an endemic disease in PNG although there is little incidence or prevalence data, and the evidence consensus for dengue presence is low. Routine surveillance is not undertaken and dengue is not a notifiable disease. Severe dengue is rarely identified by local clinicians and the reasons for this are unclear but may be related to poor recognition of dengue and a low index of suspicion, despite high incidence and prevalence rates in neighbouring countries. For example, Indonesia shares borders with PNG and regularly reports outbreaks of severe dengue and transmission of multiple DENV serotypes. DENV infection is identified in travellers from PNG however there are no data on locally circulating strains and how they may compare to viruses associated with severe dengue epidemics in other countries in the Asia Pacific region. We identified evidence for previous infection with all four DENV serotypes among people living on the northern coast of PNG, in Madang, and on Lihir Island in the Bismarck Archipelago off the northeastern coast. We also detected DENV-1, DENV-2, and DENV-3 virus in febrile patients, and we describe the first whole genome sequences of endemically circulating DENV since the prototype 1944 DENV-2_{New Guinea C} strain was characterized. Of note, severe dengue was not diagnosed in any patient infected with these viruses in PNG although introduction of the PNG DENV-3 strain into the Solomon Islands five years later resulted in a large outbreak of severe dengue with hospitalizations and deaths in that country. Dengue epidemiology and burden of disease should be investigated in PNG.

Origin and dispersal of Hepatitis E virus

Hepatitis E virus (HEV, genus Orthohepevirus) is a common cause of hepatitis worldwide. Human-infecting HEV strains (Orthohepevirus A) include human-restricted and enzootic genotypes. Viruses in the Orthohepevirus A species also infect rabbits (HEV-3ra), camels, and swine. Using a selection-informed method, we dated the origin of the Orthohepevirus genus at least 21 million years ago, whereas the Orthohepevirus A species originated in Asia, most likely from a human-infecting ancestor that existed ~4500 to 6800 years ago. In this period, the appearance of large human settlements probably facilitated HEV emergence and spread. The earliest events in Orthohepevirus A evolutionary history involved the separation of the enzootic and human-restricted genotypes, as well as the split of the camel-infecting genotypes, which occurred during the time-frame of camel domestication. The place and timing of HEV-3ra divergence also correspond to the circumstances of rabbit domestication. This study clarifies the origin and historical events underlying HEV dispersal.

Molecular evolutionary analysis of type-1 human astroviruses identifies putative sites under selection pressure on the capsid protein

Human astroviruses (HAstV) are important enteric pathogens that can be classified into eight sero/genotypes (HAstV-1 to -8). Although the various HAstV types show global spread, type-1 strains tend to be predominant. Molecular analysis of the genomic region encoding the capsid protein (ORF2) has revealed discrete sequence variation, with different lineages within each HAstV type and at least three major lineages have been identified within HAstV-1. Longitudinal epidemiological surveillance has revealed temporal shift of the various HAstV-1 lineages. Metadata analysis of HAstV-1 sequences available in the databases also revealed temporal shifts of the circulation of HAstV-1 lineages, suggesting possible antigenic-related mechanisms of selection at the sub-genotype level. By comparison of HAstV-1 capsid sequences, lineage-defining residues under positive selection were identified. Structural analysis of HAstV-1 capsid allowed identifying at least six residues exposed on the virion surface. Two residues were located in the VP34 (shell region) whilst four residues were mapped in the VP25/27 (protruding region) of HAstV capsid protein, in proximity of the putative receptor binding S site. These findings suggest that mechanisms similar to those observed and/or hypothesized for other enteric viruses are also shaping the evolution of HAstVs, with intra-typic diversification being a possible mechanism to decrease the antigenic pressure to which these viruses are exposed.

Evolutionary Analysis Predicts Sensitive Positions of MMP20 and Validates Newly- and Previously-Identified MMP20 Mutations Causing Amelogenesis Imperfecta

Amelogenesis imperfecta (AI) designates a group of genetic diseases characterized by a large range of enamel disorders causing important social and health problems. These defects can result from mutations in enamel matrix proteins or protease encoding genes. A range of mutations in the enamel cleavage enzyme matrix metalloproteinase-20 gene (MMP20) produce enamel defects of varying severity. To address how various alterations produce a range of AI phenotypes, we performed a targeted analysis to find MMP20 mutations in French patients diagnosed with non-syndromic AI. Genomic DNA was isolated from saliva and MMP20 exons and exon-intron boundaries sequenced. We identified several homozygous or heterozygous mutations, putatively involved in the AI phenotypes. To validate missense mutations and predict sensitive positions in the MMP20 sequence, we evolutionarily compared 75 sequences extracted from the public databases using the Datamonkey webserver. These sequences were representative of mammalian lineages, covering more than 150 million years of evolution. This analysis allowed us to find 324 sensitive positions (out of the 483 MMP20 residues), pinpoint functionally important domains, and build an evolutionary chart of important conserved MMP20 regions. This is an efficient tool to identify new- and previously-identified mutations. We thus identified six functional MMP20 mutations in unrelated families, finding two novel mutated sites. The genotypes and phenotypes of these six mutations are described and compared. To date, 13 MMP20 mutations causing AI have been reported, making these genotypes and associated hypomature enamel phenotypes the most frequent in AI.

Adaptive molecular evolution of MC1R gene reveals the evidence for positive diversifying selection in indigenous goat populations

Detecting signatures of selection can provide a new insight into the mechanism of contemporary breeding and artificial selection and further reveal the causal genes associated to the phenotypic variation. However, the signatures of selection on genes entailing for profitable traits between Chinese commercial and indigenous goats have been poorly interpreted. We noticed footprints of positive selection at MC1R gene containing SNPs genotyped in five Chinese native goat breeds. An experimental distribution of F_ST was built based on approximations of F_ST for each SNP across five breeds. We identified selection using the high F_ST outlier method and found that MC1R candidate gene show evidence of positive selection. Furthermore, adaptive selection pressure on specific codons was determined using different codon based on maximum‐likelihood methods; signature of positive selection in mammalian MC1R was explored in individual codons. Evolutionary analyses were inferred under maximum likelihood models, the HyPhy package implemented in the DATAMONKEY Web Server. The results of codon selection displayed positive diversifying selection at the sites were mainly involved in development of genetic variations in coat color in various mammalian species. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat MC1R provides new insights that the gene evolution may have been modulated by domestication events in goats.

Purifying selection and concerted evolution of RNA-sensing toll-like receptors in migratory waders

Migratory birds encounter a broad range of pathogens during their journeys, making them ideal models for studying immune gene evolution. Despite the potential value of these species to immunoecology and disease epidemiology, previous studies have typically focused on their adaptive immune gene repertoires. In this study, we examined the evolution of innate immune genes in three long-distance migratory waders (order Charadriiformes). We analysed two parts of the extracellular domains of two Toll-like receptors (TLR3 and TLR7) involved in virus recognition in the Sanderling (Calidris alba), Red-necked Stint (Calidris ruficollis), and Ruddy Turnstone (Arenaria interpres). Our analysis was extended to 50 avian species for which whole-genome sequences were available, including two additional waders. We found that the inferred relationships among avian TLR3 and TLR7 do not match the whole-genome phylogeny of birds. Further analyses showed that although both loci are predominantly under purifying selection, the evolution of the extracellular domain of avian TLR3 has also been driven by episodic diversifying selection. TLR7 was found to be duplicated in all five wader species and in two other orders of birds, Cuculiformes and Passeriformes. The duplication is likely to have occurred in the ancestor of each order, and the duplicated copies appear to be undergoing concerted evolution. The phylogenetic relationships of wader TLR7 matched those of the five wader species, but that of TLR3 did not. Instead, the tree inferred from TLR3 showed potential associations with the species' ecology, including migratory behaviour and exposure to pathogens. Our study demonstrates the importance of combining immunological and ecological knowledge to understand the impact of immune gene polymorphism on the evolutionary ecology of infectious diseases.

Functional evolution of the OAS1 viral sensor: Insights from old world primates

Infections with viral pathogens impose considerable selective pressure on host defensive genes. Those genes at the forefront, responsible for identifying and binding exogenous molecular viral components, will carry the hallmarks of this struggle. Oligoadenylate synthetase (OAS) enzymes play a major role in the innate defense against a large number of viruses by acting as sensors of viral infections. Following their up-regulation by the interferon pathway, OASs bind viral dsRNA and then signal ribonuclease L (RNase L) to degrade RNA, shutting down viral and host protein synthesis. We have investigated the evolution of OAS1 in twenty-two Old World monkey species. We identified a total of 35 codons with the earmarks of positive selection and we performed a comprehensive analysis of their functional significance using in silico modeling of the OAS1 protein. Subdividing OAS1 into functional domains revealed intense purifying selection in the active domain but significant positive directional selection in the RNA-binding domain (RBD), the region where OAS1 binds viral dsRNA. The modeling analysis revealed a concentration of rapidly evolving residues in one region of the RBD suggestive of the sub-functionalization of different regions of the RBD. This analysis also identified several positively selected residues circumscribing the entry to the active site suggesting adaptive evasion of viral antagonism and/or selection for production of oligoadenylate of different length.

An Analysis of Single Nucleotide Substitution in Genetic Codons - Probabilities and Outcomes

Background: Single nucleotide substitutions (SNS) in genetic codon are of prime importance due to their ability to alter an amino acid sequence as a result. Given the nature of genetic code, any SNS is expected to change the protein sequence randomly into any of the 64 possible codons. In this paper, we present a theoretical analysis of how single nucleotide substitutions in genetic codon may affect resulting amino acid residue and what is the most likely amino acid that will get selected as a result. Methods: A probability matrix was developed showing possible changes and routes likely being followed as a result of base substitution mutation causing changes at the translational level for the amino acid being encoded. Results: We observe that in event of single base pair substitution in a given amino acid; a chosen set of amino acids is theoretically more probable to be resulted suggesting a directional rather than a random change. This study also indicates that for a given amino acid coded by a number of synonymous codons, all synonymous codons will result into same list of amino acids in case of all possible SNS at three positions. Conclusion: The present work has resulted into development of a theoretical probability matrix which can be used to predict changes in amino acid residues in a protein sequence caused by single base substitutions.

Trends in substitution models of molecular evolution

Substitution models of evolution describe the process of genetic variation through fixed mutations and constitute the basis of the evolutionary analysis at the molecular level. Almost 40 years after the development of first substitution models, highly sophisticated, and data-specific substitution models continue emerging with the aim of better mimicking real evolutionary processes. Here I describe current trends in substitution models of DNA, codon and amino acid sequence evolution, including advantages and pitfalls of the most popular models. The perspective concludes that despite the large number of currently available substitution models, further research is required for more realistic modeling, especially for DNA coding and amino acid data. Additionally, the development of more accurate complex models should be coupled with new implementations and improvements of methods and frameworks for substitution model selection and downstream evolutionary analysis.

The heptad repeat region is a major selection target in MERS-CoV and related coronaviruses

Middle East respiratory syndrome coronavirus (MERS-CoV) originated in bats and spread to humans via zoonotic transmission from camels. We analyzed the evolution of the spike (S) gene in betacoronaviruses (betaCoVs) isolated from different mammals, in bat coronavirus populations, as well as in MERS-CoV strains from the current outbreak. Results indicated several positively selected sites located in the region comprising the two heptad repeats (HR1 and HR2) and their linker. Two sites (R652 and V1060) were positively selected in the betaCoVs phylogeny and correspond to mutations associated with expanded host range in other coronaviruses. During the most recent evolution of MERS-CoV, adaptive mutations in the HR1 (Q/R/H1020) arose in camels or in a previous host and spread to humans. We determined that different residues at position 1020 establish distinct inter- and intra-helical interactions and affect the stability of the six-helix bundle formed by the HRs. A similar effect on stability was observed for a nearby mutation (T1015N) that increases MERS-CoV infection efficiency in vitro. Data herein indicate that the heptad repeat region was a major target of adaptive evolution in MERS-CoV-related viruses; these results are relevant for the design of fusion inhibitor peptides with antiviral function.

Patterns of natural selection acting on the mitochondrial genome of a locally adapted fish species

Background

Mitochondrial DNA (mtDNA) is frequently used in population genetic studies and is usually considered as a neutral marker. However, given the functional importance of the proteins encoded by the mitochondrial genome, and the prominent role of mitochondria in cellular energy production, the assumption of neutrality is increasingly being questioned.

Results

We tested for evidence of selection on the mitochondrial genome of the Atlantic salmon, which is a locally adapted and widely farmed species and is distributed across a large latitudinal cline. We analysed 20 independent regions of the salmon mtDNA that represented nine genes (ND1, ND2, ND3, COX1, COX2, ATP6, ND4, ND5, and CYTB). These 20 mtDNA regions were sequenced using a 454 approach from samples collected across the entire European range of this species. We found evidence of positive selection at the ND1, ND3 and ND4 genes, which is supported by at least two different codon-based methods and also by differences in the chemical properties of the amino acids involved. The geographical distribution of some of the mutations indicated to be under selection was not random, and some mutations were private to artic populations. We discuss the possibility that selection acting on the Atlantic salmon mtDNA genome might be related to the need for increased metabolic efficiency at low temperatures.

Conclusions

The analysis of sequences representing nine mitochondrial genes that are involved in the OXPHOS pathway revealed signatures of positive selection in the mitochondrial genome of the Atlantic salmon. The properties of the amino acids involved suggest that some of the mutations that were identified to be under positive selection might have functional implications, possibly in relation to metabolic efficiency. Experimental evidence, and better understanding of regional phylogeographic structuring, are needed to clarify the potential role of selection acting on the mitochondrial genome of Atlantic salmon and other locally adapted fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0138-0) contains supplementary material, which is available to authorized users.

Selection on MHC class II supertypes in the New Zealand endemic Hochstetter's frog

BACKGROUND

The New Zealand native frogs, family Leiopelmatidae, are among the most archaic in the world. Leiopelma hochstetteri (Hochstetter's frog) is a small, semi-aquatic frog with numerous, fragmented populations scattered across New Zealand's North Island. We characterized a major histocompatibility complex (MHC) class II B gene (DAB) in L. hochstetteri from a spleen transcriptome, and then compared its diversity to neutral microsatellite markers to assess the adaptive genetic diversity of five populations ("evolutionarily significant units", ESUs).

RESULTS

L. hochstetteri possessed very high MHC diversity, with 74 DAB alleles characterized. Extremely high differentiation was observed at the DAB locus, with only two alleles shared between populations, a pattern that was not reflected in the microsatellites. Clustering analysis on putative peptide binding residues of the DAB alleles indicated four functional supertypes, all of which were represented in 4 of 5 populations, albeit at different frequencies. Otawa was an exception to these observations, with only two DAB alleles present.

CONCLUSIONS

This study of MHC diversity highlights extreme population differentiation at this functional locus. Supertype differentiation was high among populations, suggesting spatial and/or temporal variation in selection pressures. Low DAB diversity in Otawa may limit this population's adaptive potential to future pathogenic challenges.

Transmitted/Founder Viruses Rapidly Escape from CD8+ T Cell Responses in Acute Hepatitis C Virus Infection

The interaction between hepatitis C virus (HCV) and cellular immune responses during very early infection is critical for disease outcome. To date, the impact of antigen-specific cellular immune responses on the evolution of the viral population establishing infection and on potential escape has not been studied. Understanding these early host-virus dynamics is important for the development of a preventative vaccine. Three subjects who were followed longitudinally from the detection of viremia preseroconversion until disease outcome were analyzed. The evolution of transmitted/founder (T/F) viruses was undertaken using deep sequencing. CD8(+) T cell responses were measured via enzyme-linked immunosorbent spot (ELISpot) assay using HLA class I-restricted T/F epitopes. T/F viruses were rapidly extinguished in all subjects associated with either viral clearance (n = 1) or replacement with viral variants leading to establishment of chronic infection (n = 2). CD8(+) T cell responses against 11 T/F epitopes were detectable by 33 to 44 days postinfection, and 5 of these epitopes had not previously been reported. These responses declined rapidly in those who became chronically infected and were maintained in the subject who cleared infection. Higher-magnitude CD8(+) T cell responses were associated with rapid development of immune escape variants at a rate of up to 0.1 per day. Rapid escape from CD8(+) T cell responses has been quantified for the first time in the early phase of primary HCV infection. These rapid escape dynamics were associated with higher-magnitude CD8(+) T cell responses. These findings raise questions regarding optimal selection of immunogens for HCV vaccine development and suggest that detailed analysis of individual epitopes may be required.

IMPORTANCE

A major limitation in our detailed understanding of the role of immune response in HCV clearance has been the lack of data on very early primary infection when the transmitted viral variants successfully establish the acute infection. This study was made possible through the availability of specimens from a unique cohort of asymptomatic primary infection cases in whom the first available viremic samples were collected approximately 3 weeks postinfection and at regular intervals thereafter. The study included detailed examination of both the evolution of the viral population and the host cellular immune responses against the T/F viruses. The findings here provide the first evidence of host cellular responses targeting T/F variants and imposing a strong selective force toward viral escape. The results of this study provide useful insight on how virus escapes the host response and consequently on future analysis of vaccine-induced immunity.

Less is more: an adaptive branch-site random effects model for efficient detection of episodic diversifying selection

Over the past two decades, comparative sequence analysis using codon-substitution models has been honed into a powerful and popular approach for detecting signatures of natural selection from molecular data. A substantial body of work has focused on developing a class of "branch-site" models which permit selective pressures on sequences, quantified by the ω ratio, to vary among both codon sites and individual branches in the phylogeny. We develop and present a method in this class, adaptive branch-site random effects likelihood (aBSREL), whose key innovation is variable parametric complexity chosen with an information theoretic criterion. By applying models of different complexity to different branches in the phylogeny, aBSREL delivers statistical performance matching or exceeding best-in-class existing approaches, while running an order of magnitude faster. Based on simulated data analysis, we offer guidelines for what extent and strength of diversifying positive selection can be detected reliably and suggest that there is a natural limit on the optimal parametric complexity for "branch-site" models. An aBSREL analysis of 8,893 Euteleostomes gene alignments demonstrates that over 80% of branches in typical gene phylogenies can be adequately modeled with a single ω ratio model, that is, current models are unnecessarily complicated. However, there are a relatively small number of key branches, whose identities are derived from the data using a model selection procedure, for which it is essential to accurately model evolutionary complexity.

Gene-wide identification of episodic selection

We present BUSTED, a new approach to identifying gene-wide evidence of episodic positive selection, where the non-synonymous substitution rate is transiently greater than the synonymous rate. BUSTED can be used either on an entire phylogeny (without requiring an a priori hypothesis regarding which branches are under positive selection) or on a pre-specified subset of foreground lineages (if a suitable a priori hypothesis is available). Selection is modeled as varying stochastically over branches and sites, and we propose a computationally inexpensive evidence metric for identifying sites subject to episodic positive selection on any foreground branches. We compare BUSTED with existing models on simulated and empirical data. An implementation is available on www.datamonkey.org/busted, with a widget allowing the interactive specification of foreground branches.

Reduced evolutionary rates in HIV-1 reveal extensive latency periods among replicating lineages

Background

HIV-1 can persist for the duration of a patient’s life due in part to its ability to hide from the immune system, and from antiretroviral drugs, in long-lived latent reservoirs. Latent forms of HIV-1 may also be disproportionally involved in transmission. Thus, it is important to detect and quantify latency in the HIV-1 life cycle.

Results

We developed a novel molecular clock–based phylogenetic tool to investigate the prevalence of HIV-1 lineages that have experienced latency. The method removes alternative sources that may affect evolutionary rates, such as hypermutation, recombination, and selection, to reveal the contribution of generation-time effects caused by latency. Our method was able to recover latent lineages with high specificity and sensitivity, and low false discovery rates, even on relatively short branches on simulated phylogenies. Applying the tool to HIV-1 sequences from 26 patients, we show that the majority of phylogenetic lineages have been affected by generation-time effects in every patient type, whether untreated, elite controller, or under effective or failing treatment. Furthermore, we discovered extensive effects of latency in sequence data (gag, pol, and env) from reservoirs as well as in the replicating plasma population. To better understand our phylogenetic findings, we developed a dynamic model of virus-host interactions to investigate the proportion of lineages in the actively replicating population that have ever been latent. Assuming neutral evolution, our dynamic modeling showed that under most parameter conditions, it is possible for a few activated latent viruses to propagate so that in time, most HIV-1 lineages will have been latent at some time in their past.

Conclusions

These results suggest that cycling in and out of latency plays a major role in the evolution of HIV-1. Thus, no aspect of HIV-1 evolution can be fully understood without considering latency - including treatment, drug resistance, immune evasion, transmission, and pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-014-0081-0) contains supplementary material, which is available to authorized users.

A generalized mechanistic codon model

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

Mutations associated with severity of the pandemic influenza A(H1N1)pdm09 in humans: a systematic review and meta-analysis of epidemiological evidence

Mutations in the haemagglutinin (HA), non-structural protein 1 (NS1) and polymerase basic protein 2 (PB2) of influenza viruses have been associated with virulence. This study investigated the association between mutations in these genes in influenza A(H1N1)pdm09 virus and the risk of severe or fatal disease. Searches were conducted on the MEDLINE, EMBASE and Web of Science electronic databases and the reference lists of published studies. The PRISMA and STROBE guidelines were followed in assessing the quality of studies and writing-up. Eighteen (18) studies, from all continents, were included in the systematic review (recruiting patients 0 - 77 years old). The mutation D222G was associated with a significant increase in severe disease (pooled RD: 11 %, 95 % CI: 3.0 % - 18.0 %, p = 0.004) and the risk of fatality (RD: 23 %, 95 % CI: 14.0 %-31.0 %, p = < 0.0001). No association was observed between the mutations HA-D222N, D222E, PB2-E627K and NS1-T123V and severe/fatal disease. The results suggest that no virus quasispecies bearing virulence-conferring mutations in the HA, PB2 and NS1 predominated. However issues of sampling bias, and bias due to uncontrolled confounders such as comorbidities, and viral and bacterial coinfection, should be born in mind. Influenza A viruses should continue to be monitored for the occurrence of virulence-conferring mutations in HA, PB2 and NS1. There are suggestions that respiratory virus coinfections also affect virus virulence. Studies investigating the role of genetic mutations on disease outcome should make efforts to also investigate the role of respiratory virus coinfections.

Molecular evolution of the tissue-nonspecific alkaline phosphatase allows prediction and validation of missense mutations responsible for hypophosphatasia

ALPL encodes the tissue nonspecific alkaline phosphatase (TNSALP), which removes phosphate groups from various substrates. Its function is essential for bone and tooth mineralization. In humans, ALPL mutations lead to hypophosphatasia, a genetic disorder characterized by defective bone and/or tooth mineralization. To date, 275 ALPL mutations have been reported to cause hypophosphatasia, of which 204 were simple missense mutations. Molecular evolutionary analysis has proved to be an efficient method to highlight residues important for the protein function and to predict or validate sensitive positions for genetic disease. Here we analyzed 58 mammalian TNSALP to identify amino acids unchanged, or only substituted by residues sharing similar properties, through 220 millions years of mammalian evolution. We found 469 sensitive positions of the 524 residues of human TNSALP, which indicates a highly constrained protein. Any substitution occurring at one of these positions is predicted to lead to hypophosphatasia. We tested the 204 missense mutations resulting in hypophosphatasia against our predictive chart, and validated 99% of them. Most sensitive positions were located in functionally important regions of TNSALP (active site, homodimeric interface, crown domain, calcium site, …). However, some important positions are located in regions, the structure and/or biological function of which are still unknown. Our chart of sensitive positions in human TNSALP (i) enables to validate or invalidate at low cost any ALPL mutation, which would be suspected to be responsible for hypophosphatasia, by contrast with time consuming and expensive functional tests, and (ii) displays higher predictive power than in silico models of prediction.

Episodic positive selection in the evolution of avian toll-like receptor innate immunity genes

Toll-like receptors (TLRs) are a family of conserved pattern-recognition molecules responsible for initiating innate and acquired immune responses. Because they play a key role in host defence, these genes have received increasing interest in the evolutionary and population genetics literature, as their variation represents a potential target of adaptive evolution. However, the role of pathogen-mediated selection (i.e. episodic positive selection) in the evolution of these genes remains poorly known and has not been examined outside of mammals. A recent increase in the number of bird species for which TLR sequences are available has enabled us to examine the selective processes that have influenced evolution of the 10 known avian TLR genes. Specifically, we tested for episodic positive selection to identify codons that experience purifying selection for the majority of their evolution, interspersed with bursts of positive selection that may occur only in restricted lineages. We included up to 23 species per gene (mean = 16.0) and observed that, although purifying selection was evident, an average of 4.5% of codons experienced episodic positive selection across all loci. For four genes in which sequence coverage traversed both the extracellular leucine-rich repeat region (LRR) and transmembrane/intracellular domains of the proteins, increased positive selection was observed at the extracellular domain, consistent with theoretical predictions. Our results provide evidence that episodic positive selection has played an important role in the evolution of most avian TLRs, consistent with the role of these loci in pathogen recognition and a mechanism of host-pathogen coevolution.

Intragenic homogenization and multiple copies of prey-wrapping silk genes in Argiope garden spiders

Background

Spider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals.

Results

We present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region.

Conclusions

Using the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

Novel gammaherpesviruses in North American domestic cats, bobcats, and pumas: identification, prevalence, and risk factors

Gammaherpesviruses (GHVs) are a diverse and rapidly expanding group of viruses associated with a variety of disease conditions in humans and animals. To identify felid GHVs, we screened domestic cat (Felis catus), bobcat (Lynx rufus), and puma (Puma concolor) blood cell DNA samples from California, Colorado, and Florida using a degenerate pan-GHV PCR. Additional pan-GHV and long-distance PCRs were used to sequence a contiguous 3.4-kb region of each putative virus species, including partial glycoprotein B and DNA polymerase genes. We identified three novel GHVs, each present predominantly in one felid species: Felis catus GHV 1 (FcaGHV1) in domestic cats, Lynx rufus GHV 1 (LruGHV1) in bobcats, and Puma concolor GHV 1 (PcoGHV1) in pumas. To estimate infection prevalence, we developed real-time quantitative PCR assays for each virus and screened additional DNA samples from all three species (n = 282). FcaGHV1 was detected in 16% of domestic cats across all study sites. LruGHV1 was detected in 47% of bobcats and 13% of pumas across all study sites, suggesting relatively common interspecific transmission. PcoGHV1 was detected in 6% of pumas, all from a specific region of Southern California. The risk of infection for each host varied with geographic location. Age was a positive risk factor for bobcat LruGHV1 infection, and age and being male were risk factors for domestic cat FcaGHV1 infection. Further characterization of these viruses may have significant health implications for domestic cats and may aid studies of free-ranging felid ecology.

IMPORTANCE

Gammaherpesviruses (GHVs) establish lifelong infection in many animal species and can cause cancer and other diseases in humans and animals. In this study, we identified the DNA sequences of three GHVs present in the blood of domestic cats (Felis catus), bobcats (Lynx rufus), and pumas (Puma concolor; also known as mountain lions, cougars, and panthers). We found that these viruses were closely related to, but distinct from, other known GHVs of animals and represent the first GHVs identified to be native to these feline species. We developed techniques to rapidly and specifically detect the DNA of these viruses in feline blood and found that the domestic cat and bobcat viruses were widespread across the United States. In contrast, puma virus was found only in a specific region of Southern California. Surprisingly, the bobcat virus was also detected in some pumas, suggesting relatively common virus transmission between these species. Adult domestic cats and bobcats were at greater risk for infection than juveniles. Male domestic cats were at greater risk for infection than females. This study identifies three new viruses that are widespread in three feline species, indicates risk factors for infection that may relate to the route of infection, and demonstrates cross-species transmission between bobcats and pumas. These newly identified viruses may have important effects on feline health and ecology.

Dietary and flight energetic adaptations in a salivary gland transcriptome of an insectivorous bat

We hypothesized that evolution of salivary gland secretory proteome has been important in adaptation to insectivory, the most common dietary strategy among Chiroptera. A submandibular salivary gland (SMG) transcriptome was sequenced for the little brown bat, Myotis lucifugus. The likely secretory proteome of 23 genes included seven (RETNLB, PSAP, CLU, APOE, LCN2, C3, CEL) related to M. lucifugus insectivorous diet and metabolism. Six of the secretory proteins probably are endocrine, whereas one (CEL) most likely is exocrine. The encoded proteins are associated with lipid hydrolysis, regulation of lipid metabolism, lipid transport, and insulin resistance. They are capable of processing exogenous lipids for flight metabolism while foraging. Salivary carboxyl ester lipase (CEL) is thought to hydrolyze insect lipophorins, which probably are absorbed across the gastric mucosa during feeding. The other six proteins are predicted either to maintain these lipids at high blood concentrations or to facilitate transport and uptake by flight muscles. Expression of these seven genes and coordinated secretion from a single organ is novel to this insectivorous bat, and apparently has evolved through instances of gene duplication, gene recruitment, and nucleotide selection. Four of the recruited genes are single-copy in the Myotis genome, whereas three have undergone duplication(s) with two of these genes exhibiting evolutionary 'bursts' of duplication resulting in multiple paralogs. Evidence for episodic directional selection was found for six of seven genes, reinforcing the conclusion that the recruited genes have important roles in adaptation to insectivory and the metabolic demands of flight. Intragenic frequencies of mobile- element-like sequences differed from frequencies in the whole M. lucifugus genome. Differences among recruited genes imply separate evolutionary trajectories and that adaptation was not a single, coordinated event.

Urotensin-II system in genetic control of blood pressure and renal function

Urotensin-II controls ion/water homeostasis in fish and vascular tone in rodents. We hypothesised that common genetic variants in urotensin-II pathway genes are associated with human blood pressure or renal function. We performed family-based analysis of association between blood pressure, glomerular filtration and genes of the urotensin-II pathway (urotensin-II, urotensin-II related peptide, urotensin-II receptor) saturated with 28 tagging single nucleotide polymorphisms in 2024 individuals from 520 families; followed by an independent replication in 420 families and 7545 unrelated subjects. The expression studies of the urotensin-II pathway were carried out in 97 human kidneys. Phylogenetic evolutionary analysis was conducted in 17 vertebrate species. One single nucleotide polymorphism (rs531485 in urotensin-II gene) was associated with adjusted estimated glomerular filtration rate in the discovery cohort (p = 0.0005). It showed no association with estimated glomerular filtration rate in the combined replication resource of 8724 subjects from 6 populations. Expression of urotensin-II and its receptor showed strong linear correlation (r = 0.86, p<0.0001). There was no difference in renal expression of urotensin-II system between hypertensive and normotensive subjects. Evolutionary analysis revealed accumulation of mutations in urotensin-II since the divergence of primates and weaker conservation of urotensin-II receptor in primates than in lower vertebrates. Our data suggest that urotensin-II system genes are unlikely to play a major role in genetic control of human blood pressure or renal function. The signatures of evolutionary forces acting on urotensin-II system indicate that it may have evolved towards loss of function since the divergence of primates.

Superiority of a mechanistic codon substitution model even for protein sequences in phylogenetic analysis

Background

Nucleotide and amino acid substitution tendencies are characteristic of each species, organelle, and protein family. Hence, various empirical amino acid substitution rate matrices have needed to be estimated for phylogenetic analysis: JTT, WAG, and LG for nuclear proteins, mtREV for mitochondrial proteins, cpREV10 and cpREV64 for chloroplast-encoded proteins, and FLU for influenza proteins. On the other hand, in a mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the ratio of fixation depending on the type of amino acid replacement, mutation rates and the strength of selective constraint on amino acids can be tailored to each protein family with additional 11 parameters. As a result, in the evolutionary analysis of codon sequences it outperforms codon substitution models equivalent to empirical amino acid substitution matrices. Is it superior even for amino acid sequences, among which synonymous substitutions cannot be identified?

Results

Nucleotide mutations are assumed to occur independently of codon positions but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene with a linear function of a given estimate of selective constraints, which were estimated by maximizing the likelihood of an empirical amino acid or codon substitution frequency matrix, each of JTT, WAG, LG, and KHG. It is shown that the mechanistic codon substitution model with the assumption of equal codon usage yields better values of Akaike and Bayesian information criteria for all three phylogenetic trees of mitochondrial, chloroplast, and influenza-A hemagglutinin proteins than the empirical amino acid substitution models with mtREV, cpREV64, and FLU, which were designed specifically for those protein families, respectively. The variation of selective constraint across sites fits the datasets significantly better than variable codon mutation rates, confirming that substitution rate variations across sites detected by amino acid substitution models are caused primarily by the variation of selective constraint against amino acid substitutions rather than the variation of codon mutation rate.

Conclusions

The mechanistic codon substitution model is superior to amino acid substitution models even in the evolutionary analysis of protein sequences.

Neocortical development as an evolutionary platform for intragenomic conflict

Embryonic development in mammals has evolved a platform for genomic conflict between mothers and embryos and, by extension, between maternal and paternal genomes. The evolutionary interests of the mother and embryo may be maximized through the promotion of sex-chromosome genes and imprinted alleles, resulting in the rapid evolution of postzygotic phenotypes preferential to either the maternal or paternal genome. In eutherian mammals, extraordinary in utero maternal investment in the brain, and neocortex especially, suggests that convergent evolution of an expanded mammalian neocortex along divergent lineages may be explained, in part, by parent-of-origin-linked gene expression arising from parent-offspring conflict. The influence of this conflict on neocortical development and evolution, however, has not been investigated at the genomic level. In this hypothesis and theory article, we provide preliminary evidence for positive selection in humans in the regions of two platforms of intragenomic conflict—chromosomes 15q11-q13 and X—and explore the potential relevance of cis-regulated imprinted domains to neocortical expansion in mammalian evolution. We present the hypothesis that maternal- and paternal-specific pressures on the developing neocortex compete intragenomically to influence neocortical expansion in mammalian evolution.

Hippo pathway phylogenetics predicts monoubiquitylation of Salvador and Merlin/Nf2

Recently we employed phylogenetics to predict that the cellular interpretation of TGF-β signals is modulated by monoubiquitylation cycles affecting the Smad4 signal transducer/tumor suppressor. This prediction was subsequently validated by experiments in flies, frogs and mammalian cells. Here we apply a phylogenetic approach to the Hippo pathway and predict that two of its signal transducers, Salvador and Merlin/Nf2 (also a tumor suppressor) are regulated by monoubiquitylation. This regulatory mechanism does not lead to protein degradation but instead serves as a highly efficient “off/on” switch when the protein is subsequently deubiquitylated. Overall, our study shows that the creative application of phylogenetics can predict new roles for pathway components and new mechanisms for regulating intercellular signaling pathways.

Estimating empirical codon hidden Markov models

Empirical codon models (ECMs) estimated from a large number of globular protein families outperformed mechanistic codon models in their description of the general process of protein evolution. Among other factors, ECMs implicitly model the influence of amino acid properties and multiple nucleotide substitutions (MNS). However, the estimation of ECMs requires large quantities of data, and until recently, only few suitable data sets were available. Here, we take advantage of several new Drosophila species genomes to estimate codon models from genome-wide data. The availability of large numbers of genomes over varying phylogenetic depths in the Drosophila genus allows us to explore various divergence levels. In consequence, we can use these data to determine the appropriate level of divergence for the estimation of ECMs, avoiding overestimation of MNS rates caused by saturation. To account for variation in evolutionary rates along the genome, we develop new empirical codon hidden Markov models (ecHMMs). These models significantly outperform previous ones with respect to maximum likelihood values, suggesting that they provide a better fit to the evolutionary process. Using ECMs and ecHMMs derived from genome-wide data sets, we devise new likelihood ratio tests (LRTs) of positive selection. We found classical LRTs very sensitive to the presence of MNSs, showing high false-positive rates, especially with small phylogenies. The new LRTs are more conservative than the classical ones, having acceptable false-positive rates and reduced power.

Codon substitution models based on residue similarity and their applications

Codon models are now widely used to draw evolutionary inferences from alignments of homologous sequence data. Incorporating physicochemical properties of amino acids into codon models, two novel codon substitution models describing the evolution of protein-coding DNA sequences are presented based on the similarity scores of amino acids. To describe substitutions between codons a continue-time Markov process is used. Transition/transversion rate bias and nonsynonymous codon usage bias are allowed in the models. In our implementation, the parameters are estimated by maximum-likelihood (ML) method as in previous studies. Furthermore, instantaneous mutations involving more than one nucleotide position of a codon are considered in the second model. Then the two suggested models are applied to five real data sets. The analytic results indicate that the new codon models considering physicochemical properties of amino acids can provide a better fit to the data comparing with existing codon models, and then produce more reliable estimates of certain biologically important measures than existing methods.

Modeling coding-sequence evolution within the context of residue solvent accessibility

Background

Protein structure mediates site-specific patterns of sequence divergence. In particular, residues in the core of a protein (solvent-inaccessible residues) tend to be more evolutionarily conserved than residues on the surface (solvent-accessible residues).

Results

Here, we present a model of sequence evolution that explicitly accounts for the relative solvent accessibility of each residue in a protein. Our model is a variant of the Goldman-Yang 1994 (GY94) model in which all model parameters can be functions of the relative solvent accessibility (RSA) of a residue. We apply this model to a data set comprised of nearly 600 yeast genes, and find that an evolutionary-rate ratio ω that varies linearly with RSA provides a better model fit than an RSA-independent ω or an ω that is estimated separately in individual RSA bins. We further show that the branch length t and the transition-transverion ratio κ also vary with RSA. The RSA-dependent GY94 model performs better than an RSA-dependent Muse-Gaut 1994 (MG94) model in which the synonymous and non-synonymous rates individually are linear functions of RSA. Finally, protein core size affects the slope of the linear relationship between ω and RSA, and gene expression level affects both the intercept and the slope.

Conclusions

Structure-aware models of sequence evolution provide a significantly better fit than traditional models that neglect structure. The linear relationship between ω and RSA implies that genes are better characterized by their ω slope and intercept than by just their mean ω.

Evolutionary analysis suggests that AMTN is enamel-specific and a candidate for AI

Molecular evolutionary analysis is an efficient method to predict and/or validate amino acid substitutions that could lead to a genetic disease and to highlight residues and motifs that could play an important role in the protein structure and/or function. We have applied such analysis to amelotin (AMTN), a recently identified enamel protein in the rat, mouse, and humans. An in silico search for AMTN provided 42 new mammalian sequences that were added to the 3 published sequences with which we performed the analysis using a dataset representative of all lineages (circa 220 million years of evolution), including 2 enamel-less species, sloth and armadillo. During evolution, of the 209 residues of human AMTN, 17 were unchanged and 34 had conserved their chemical properties. Substituting these important residues could lead to amelogenesis imperfecta (AI). Also, AMTN possesses a well-conserved signal peptide, 2 conserved motifs whose function is certainly important but unknown, and a putative phosphorylation site (SXE). In addition, the sequences of the 2 enamel-less species display mutations revealing that AMTN underwent pseudogenization, which suggests that AMTN is an enamel-specific protein.

Modeling HIV-1 drug resistance as episodic directional selection

The evolution of substitutions conferring drug resistance to HIV-1 is both episodic, occurring when patients are on antiretroviral therapy, and strongly directional, with site-specific resistant residues increasing in frequency over time. While methods exist to detect episodic diversifying selection and continuous directional selection, no evolutionary model combining these two properties has been proposed. We present two models of episodic directional selection (MEDS and EDEPS) which allow the a priori specification of lineages expected to have undergone directional selection. The models infer the sites and target residues that were likely subject to directional selection, using either codon or protein sequences. Compared to its null model of episodic diversifying selection, MEDS provides a superior fit to most sites known to be involved in drug resistance, and neither one test for episodic diversifying selection nor another for constant directional selection are able to detect as many true positives as MEDS and EDEPS while maintaining acceptable levels of false positives. This suggests that episodic directional selection is a better description of the process driving the evolution of drug resistance.

When exposed to treatment, HIV-1 and other rapidly evolving viruses have the capacity to acquire drug resistance mutations (DRAMs), which limit the efficacy of antivirals. There are a number of experimentally well characterized HIV-1 DRAMs, but many mutations whose roles are not fully understood have also been reported. In this manuscript we construct evolutionary models that identify the locations and targets of mutations conferring resistance to antiretrovirals from viral sequences sampled from treated and untreated individuals. While the evolution of drug resistance is a classic example of natural selection, existing analyses fail to detect the majority of DRAMs. We show that, in order to identify resistance mutations from sequence data, it is necessary to recognize that in this case natural selection is both episodic (it only operates when the virus is exposed to the drugs) and directional (only mutations to a particular amino-acid confer resistance while allowing the virus to continue replicating). The new class of models that allow for the episodic and directional nature of adaptive evolution performs very well at recovering known DRAMs, can be useful at identifying unknown resistance-associated mutations, and is generally applicable to a variety of biological scenarios where similar selective forces are at play.

Evolution of viral genomes: interplay between selection, recombination, and other forces

RNA viruses evolve very rapidly, often recombine, and are subject to strong host (immune response) and anthropogenic (antiretroviral drugs) selective forces. Given their compact and extensively sequenced genomes, comparative analysis of RNA viral data can provide important insights into the molecular mechanisms of adaptation, pathogenicity, immune evasion, and drug resistance. In this chapter, we present an example-based overview of recent advances in evolutionary models and statistical approaches that enable screening viral alignments for evidence of adaptive change in the presence of recombination, detecting bursts of directional adaptive evolution associated with the phenotypic changes, and detecting of coevolving sites in viral genes.

Advantages of a mechanistic codon substitution model for evolutionary analysis of protein-coding sequences

BACKGROUND

A mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the average fixation probability depending on the type of amino acid replacement, has advantages over nucleotide, amino acid, and empirical codon substitution models in evolutionary analysis of protein-coding sequences. It can approximate a wide range of codon substitution processes. If no selection pressure on amino acids is taken into account, it will become equivalent to a nucleotide substitution model. If mutation rates are assumed not to depend on the codon type, then it will become essentially equivalent to an amino acid substitution model. Mutation at the nucleotide level and selection at the amino acid level can be separately evaluated.

RESULTS

The present scheme for single nucleotide mutations is equivalent to the general time-reversible model, but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene in a linear function of a given estimate of selective constraints. Their good estimates are those calculated by maximizing the respective likelihoods of empirical amino acid or codon substitution frequency matrices. Akaike and Bayesian information criteria indicate that the present model performs far better than the other substitution models for all five phylogenetic trees of highly-divergent to highly-homologous sequences of chloroplast, mitochondrial, and nuclear genes. It is also shown that multiple nucleotide changes in infinitesimal time are significant in long branches, although they may be caused by compensatory substitutions or other mechanisms. The variation of selective constraint over sites fits the datasets significantly better than variable mutation rates, except for 10 slow-evolving nuclear genes of 10 mammals. An critical finding for phylogenetic analysis is that assuming variable mutation rates over sites lead to the overestimation of branch lengths.

Non-negative matrix factorization for learning alignment-specific models of protein evolution

Models of protein evolution currently come in two flavors: generalist and specialist. Generalist models (e.g. PAM, JTT, WAG) adopt a one-size-fits-all approach, where a single model is estimated from a number of different protein alignments. Specialist models (e.g. mtREV, rtREV, HIVbetween) can be estimated when a large quantity of data are available for a single organism or gene, and are intended for use on that organism or gene only. Unsurprisingly, specialist models outperform generalist models, but in most instances there simply are not enough data available to estimate them. We propose a method for estimating alignment-specific models of protein evolution in which the complexity of the model is adapted to suit the richness of the data. Our method uses non-negative matrix factorization (NNMF) to learn a set of basis matrices from a general dataset containing a large number of alignments of different proteins, thus capturing the dimensions of important variation. It then learns a set of weights that are specific to the organism or gene of interest and for which only a smaller dataset is available. Thus the alignment-specific model is obtained as a weighted sum of the basis matrices. Having been constrained to vary along only as many dimensions as the data justify, the model has far fewer parameters than would be required to estimate a specialist model. We show that our NNMF procedure produces models that outperform existing methods on all but one of 50 test alignments. The basis matrices we obtain confirm the expectation that amino acid properties tend to be conserved, and allow us to quantify, on specific alignments, how the strength of conservation varies across different properties. We also apply our new models to phylogeny inference and show that the resulting phylogenies are different from, and have improved likelihood over, those inferred under standard models.

Potent autologous and heterologous neutralizing antibody responses occur in HIV-2 infection across a broad range of infection outcomes

Few studies have explored the role of neutralizing antibody (NAb) responses in controlling HIV-2 viremia and disease progression. Using a TZM-bl neutralization assay, we assessed heterologous and autologous NAb responses from a community cohort of HIV-2-infected individuals with a broad range of disease outcomes in rural Guinea-Bissau. All subjects (n = 40) displayed exceptionally high heterologous NAb titers (50% inhibitory plasma dilution or 50% inhibitory concentration [IC(50)], 1:7,000 to 1:1,000,000) against 5 novel primary HIV-2 envelopes and HIV-2 7312A, whereas ROD A and 3 primary envelopes were relatively resistant to neutralization. Most individuals also showed high autologous NAb against contemporaneous envelopes (78% of plasma-envelope combinations in 69 envelopes from 21 subjects), with IC(50)s above 1:10,000. No association between heterologous or autologous NAb titer and greater control of HIV-2 was found. A subset of envelopes was found to be more resistant to neutralization (by plasma and HIV-2 monoclonal antibodies). These envelopes were isolated from individuals with greater intrapatient sequence diversity and were associated with changes in potential N-linked glycosylation sites but not CD4 independence or CXCR4 use. Plasma collected from up to 15 years previously was able to potently neutralize recent autologous envelopes, suggesting a lack of escape from NAb and the persistence of neutralization-sensitive variants over time, despite significant NAb pressure. We conclude that despite the presence of broad and potent NAb responses in HIV-2-infected individuals, these are not the primary forces behind the dichotomous outcomes observed but reveal a limited capacity for adaptive selection and escape from host immunity in HIV-2 infection.