Adaptive Amino Acid Changes Research Articles

Contrasting Patterns in the Early Stage of SARS-CoV-2 Evolution between Humans and Minks.

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor-binding domain and receptor-binding motif. An excess of high-frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low-frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies that the virus may not be preadapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.

Large Scale Analyses and Visualization of Adaptive Amino Acid Changes Projects

When changes at few amino acid sites are the target of selection, adaptive amino acid changes in protein sequences can be identified using maximum-likelihood methods based on models of codon substitution (such as codeml). Although such methods have been employed numerous times using a variety of different organisms, the time needed to collect the data and prepare the input files means that tens or hundreds of coding regions are usually analyzed. Nevertheless, the recent availability of flexible and easy to use computer applications that collect relevant data (such as BDBM) and infer positively selected amino acid sites (such as ADOPS), means that the entire process is easier and quicker than before. However, the lack of a batch option in ADOPS, here reported, still precludes the analysis of hundreds or thousands of sequence files. Given the interest and possibility of running such large-scale projects, we have also developed a database where ADOPS projects can be stored. Therefore, this study also presents the B+ database, which is both a data repository and a convenient interface that looks at the information contained in ADOPS projects without the need to download and unzip the corresponding ADOPS project file. The ADOPS projects available at B+ can also be downloaded, unzipped, and opened using the ADOPS graphical interface. The availability of such a database ensures results repeatability, promotes data reuse with significant savings on the time needed for preparing datasets, and effortlessly allows further exploration of the data contained in ADOPS projects.

Differential Contribution of Subunit Interfaces to α9α10 Nicotinic Acetylcholine Receptor Function.

Nicotinic acetylcholine receptors can be assembled from either homomeric or heteromeric pentameric subunit combinations. At the interface of the extracellular domains of adjacent subunits lies the acetylcholine binding site, composed of a principal component provided by one subunit and a complementary component of the adjacent subunit. Compared with neuronal nicotinic acetylcholine cholinergic receptors (nAChRs) assembled from α and β subunits, the α9α10 receptor is an atypical member of the family. It is a heteromeric receptor composed only of α subunits. Whereas mammalian α9 subunits can form functional homomeric α9 receptors, α10 subunits do not generate functional channels when expressed heterologously. Hence, it has been proposed that α10 might serve as a structural subunit, much like a β subunit of heteromeric nAChRs, providing only complementary components to the agonist binding site. Here, we have made use of site-directed mutagenesis to examine the contribution of subunit interface domains to α9α10 receptors by a combination of electrophysiological and radioligand binding studies. Characterization of receptors containing Y190T mutations revealed unexpectedly that both α9 and α10 subunits equally contribute to the principal components of the α9α10 nAChR. In addition, we have shown that the introduction of a W55T mutation impairs receptor binding and function in the rat α9 subunit but not in the α10 subunit, indicating that the contribution of α9 and α10 subunits to complementary components of the ligand-binding site is nonequivalent. We conclude that this asymmetry, which is supported by molecular docking studies, results from adaptive amino acid changes acquired only during the evolution of mammalian α10 subunits.

Viruses are a dominant driver of protein adaptation in mammals.

Viruses interact with hundreds to thousands of proteins in mammals, yet adaptation against viruses has only been studied in a few proteins specialized in antiviral defense. Whether adaptation to viruses typically involves only specialized antiviral proteins or affects a broad array of virus-interacting proteins is unknown. Here, we analyze adaptation in ~1300 virus-interacting proteins manually curated from a set of 9900 proteins conserved in all sequenced mammalian genomes. We show that viruses (i) use the more evolutionarily constrained proteins within the cellular functions they interact with and that (ii) despite this high constraint, virus-interacting proteins account for a high proportion of all protein adaptation in humans and other mammals. Adaptation is elevated in virus-interacting proteins across all functional categories, including both immune and non-immune functions. We conservatively estimate that viruses have driven close to 30% of all adaptive amino acid changes in the part of the human proteome conserved within mammals. Our results suggest that viruses are one of the most dominant drivers of evolutionary change across mammalian and human proteomes.

Conservation and Innovation of APOBEC3A Restriction Functions during Primate Evolution.

LINE-1 (long interspersed element-1) retroelements are the only active autonomous endogenous retroelements in human genomes. Their retrotransposition activity has created close to 50% of the current human genome. Due to the apparent costs of this proliferation, host genomes have evolved multiple mechanisms to curb LINE-1 retrotransposition. Here, we investigate the evolution and function of the LINE-1 restriction factor APOBEC3A, a member of the APOBEC3 cytidine deaminase gene family. We find that APOBEC3A genes have evolved rapidly under diversifying selection in primates, suggesting changes in APOBEC3A have been recurrently selected in a host-pathogen "arms race." Nonetheless, in contrast to previous reports, we find that the LINE-1 restriction activity of APOBEC3A proteins has been strictly conserved throughout simian primate evolution in spite of its pervasive diversifying selection. Based on these results, we conclude that LINE-1s have not driven the rapid evolution of APOBEC3A in primates. In contrast to this conserved LINE-1 restriction, we find that a subset of primate APOBEC3A genes have enhanced antiviral restriction. We trace this gain of antiviral restriction in APOBEC3A to the common ancestor of a subset of Old World monkeys. Thus, APOBEC3A has not only maintained its LINE-1 restriction ability, but also evolved a gain of antiviral specificity against other pathogens. Our findings suggest that while APOBEC3A has evolved to restrict additional pathogens, only those adaptive amino acid changes that leave LINE-1 restriction unperturbed have been tolerated.

Adaptive amino acid substitutions enhance the virulence of a novel human H7N9 influenza virus in mice

To identify molecular features that confer enhanced H7N9 virulence in mammals, we independently generated three mouse-adapted variants of A/Shanghai/2/2013 (H7N9) by serial passage in mice. The mouse lethal doses (MLD50) of the mouse-adapted variants were reduced >1000–100000-fold when compared to the parental virus. Adapted variants displayed enhanced replication kinetics in vivo, and were capable of replicating in multiple organs. Analysis of adapted viral genomes revealed a total of 14 amino acid changes among the three variant viruses in the PA (T97I, K328R, P332T, and Q556R), HA (H3 numbering; A107T, R220I, L226Q, and R354K), NP (A284T and M352I), NA (M26I, N142S, and G389D), and M1 (M128R) proteins. Notably, many of these adaptive amino acid changes have been identified in naturally occurring H7 isolates. Our results identify amino acid substitutions that collectively enhance the ability of a human H7N9 virus to replicate and cause severe disease in mice.

Adaptive Mutations Enhance Assembly and Cell-to-Cell Transmission of a High-Titer Hepatitis C Virus Genotype 5a Core-NS2 JFH1-Based Recombinant.

Recombinant hepatitis C virus (HCV) clones propagated in human hepatoma cell cultures yield relatively low infectivity titers. Here, we adapted the JFH1-based Core-NS2 recombinant SA13/JFH1C3405G,A3696G (termed SA13/JFH1orig), of the poorly characterized genotype 5a, to Huh7.5 cells, yielding a virus with greatly improved spread kinetics and an infectivity titer of 6.7 log10 focus-forming units (FFU)/ml. We identified several putative adaptive amino acid changes. In head-to-head infections at fixed multiplicities of infection, one SA13/JFH1orig mutant termed SA13/JFH1Core-NS5B, containing 13 amino acid changes (R114W and V187A [Core]; V235L [E1]; T385P [E2]; L782V [p7]; Y900C [NS2]; N2034D, E2238G, V2252A, L2266P, and I2340T [NS5A]; A2500S and V2841A [NS5B]), displayed fitness comparable to that of the polyclonal high-titer adapted virus. Single-cycle virus production assays in CD81-deficient Huh7-derived cells demonstrated that these changes did not affect replication but increased HCV assembly and specific infectivity as early as 24 h posttransfection. Infectious coculture assays in Huh7.5 cells showed a significant increase in cell-to-cell transmission for SA13/JFH1Core-NS5B viruses as well as viruses with only p7 and nonstructural protein mutations. Interestingly, the E2 hypervariable region 1 (HVR1) mutation T385P caused (i) increased sensitivity to neutralizing patient IgG and human monoclonal antibodies AR3A and AR4A and (ii) increased accessibility of the CD81 binding site without affecting the usage of CD81 and SR-BI. We finally demonstrated that SA13/JFH1orig and SA13/JFH1Core-NS5B, with and without the E2 mutation T385P, displayed similar biophysical properties following iodixanol gradient ultracentrifugation. This study has implications for investigations requiring high virus concentrations, such as studies of HCV particle composition and development of whole-virus vaccine antigens. Hepatitis C virus (HCV) is a major global health care burden, affecting more than 150 million people worldwide. These individuals are at high risk of developing severe end-stage liver diseases. No vaccine exists. While it is possible to produce HCV particles resembling isolates of all HCV genotypes in human hepatoma cells (HCVcc), production efficacy varies. Thus, for several important studies, including vaccine development, in vitro systems enabling high-titer production of diverse HCV strains would be advantageous. Our study offers important functional data on how cell culture-adaptive mutations identified in genotype 5a JFH1-based HCVcc permit high-titer culture by affecting HCV genesis through increasing virus assembly and HCV fitness by enhancing the virus specific infectivity and cell-to-cell transmission ability, without influencing the biophysical particle properties. High-titer HCVcc like the one described in this study may be pivotal in future vaccine-related studies where large quantities of infectious HCV particles are necessary.

Molecular evolution of psbA gene in ferns: unraveling selective pressure and co-evolutionary pattern

BackgroundThe photosynthetic oxygen-evolving photo system II (PS II) produces almost the entire oxygen in the atmosphere. This unique biochemical system comprises a functional core complex that is encoded by psbA and other genes. Unraveling the evolutionary dynamics of this gene is of particular interest owing to its direct role in oxygen production. psbA underwent gene duplication in leptosporangiates, in which both copies have been preserved since. Because gene duplication is often followed by the non-fictionalization of one of the copies and its subsequent erosion, preservation of both psbA copies pinpoint functional or regulatory specialization events. The aim of this study was to investigate the molecular evolution of psbA among fern lineages.ResultsWe sequenced psbA , which encodes D1 protein in the core complex of PSII, in 20 species representing 8 orders of extant ferns; then we searched for selection and convolution signatures in psbA across the 11 fern orders. Collectively, our results indicate that: (1) selective constraints among D1 protein relaxed after the duplication in 4 leptosporangiate orders; (2) a handful positively selected codons were detected within species of single copy psbA, but none in duplicated ones; (3) a few sites among D1 protein were involved in co-evolution process which may intimate significant functional/structural communications between them.ConclusionsThe strong competition between ferns and angiosperms for light may have been the main cause for a continuous fixation of adaptive amino acid changes in psbA , in particular after its duplication. Alternatively, a single psbA copy may have undergone bursts of adaptive changes at the molecular level to overcome angiosperms competition. The strong signature of positive Darwinian selection in a major part of D1 protein is testament to this. At the same time, species own two psbA copies hardly have positive selection signals among the D1 protein coding sequences. In this study, eleven co-evolving sites have been detected via different molecules, which may be more important than others.

Evolution of cagA Oncogene of Helicobacter pylori through Recombination

Helicobacter pylori is a gastric pathogen that infects half the human population and causes gastritis, ulcers, and cancer. The cagA gene product is a major virulence factor associated with gastric cancer. It is injected into epithelial cells, undergoes phosphorylation by host cell kinases, and perturbs host signaling pathways. CagA is known for its geographical, structural, and functional diversity in the C-terminal half, where an EPIYA host-interacting motif is repeated. The Western version of CagA carries the EPIYA segment types A, B, and C, while the East Asian CagA carries types A, B, and D and shows higher virulence. Many structural variants such as duplications and deletions are reported. In this study, we gained insight into the relationships of CagA variants through various modes of recombination, by analyzing all known cagA variants at the DNA sequence level with the single nucleotide resolution. Processes that occurred were: (i) homologous recombination between DNA sequences for CagA multimerization (CM) sequence; (ii) recombination between DNA sequences for the EPIYA motif; and (iii) recombination between short similar DNA sequences. The left half of the EPIYA-D segment characteristic of East Asian CagA was derived from Western type EPIYA, with Amerind type EPIYA as the intermediate, through rearrangements of specific sequences within the gene. Adaptive amino acid changes were detected in the variable region as well as in the conserved region at sites to which no specific function has yet been assigned. Each showed a unique evolutionary distribution. These results clarify recombination-mediated routes of cagA evolution and provide a solid basis for a deeper understanding of its function in pathogenesis.

Adaptive Mutations in a Human Immunodeficiency Virus Type 1 Envelope Protein with a Truncated V3 Loop Restore Function by Improving Interactions with CD4

We previously reported that a human immunodeficiency virus type 1 (HIV-1) clade B envelope protein with a severely truncated V3 loop regained function after passage in tissue culture. The adapted virus, termed TA1, retained the V3 truncation, was exquisitely sensitive to neutralization by the CD4 binding site monoclonal antibody b12 and by HIV-positive human sera, used CCR5 to enter cells, and was completely resistant to small molecule CCR5 antagonists. To examine the mechanistic basis for these properties, we singly and in combination introduced each of the 5 mutations from the adapted clone TA1 into the unadapted envelope. We found that single amino acid changes in the C3 region, the V3 loop, and in the fusion peptide were responsible for imparting near-normal levels of envelope function to TA1. T342A, which resulted in the loss of a highly conserved glycosylation site in C3, played the primary role. The adaptive amino acid changes had no impact on CCR5 antagonist resistance but made virus more sensitive to neutralization by antibodies to the CD4 binding site, modestly enhanced affinity for CD4, and made TA1 more responsive to CD4 binding. Specifically, TA1 was triggered by soluble CD4 more readily than the parental Env and, unlike the parental Env, could mediate entry on cells that express low levels of CD4. In contrast, TA1 interacted with CCR5 less efficiently and was highly sensitive to antibodies that bind to the CCR5 N terminus and ECL2. Therefore, enhanced utilization of CD4 is one mechanism by which HIV-1 can overcome mutations in the V3 region that negatively affect CCR5 interactions.

Mitochondrial Targeting Adaptation of the Hominoid-Specific Glutamate Dehydrogenase Driven by Positive Darwinian Selection

Many new gene copies emerged by gene duplication in hominoids, but little is known with respect to their functional evolution. Glutamate dehydrogenase (GLUD) is an enzyme central to the glutamate and energy metabolism of the cell. In addition to the single, GLUD-encoding gene present in all mammals (GLUD1), humans and apes acquired a second GLUD gene (GLUD2) through retroduplication of GLUD1, which codes for an enzyme with unique, potentially brain-adapted properties. Here we show that whereas the GLUD1 parental protein localizes to mitochondria and the cytoplasm, GLUD2 is specifically targeted to mitochondria. Using evolutionary analysis and resurrected ancestral protein variants, we demonstrate that the enhanced mitochondrial targeting specificity of GLUD2 is due to a single positively selected glutamic acid-to-lysine substitution, which was fixed in the N-terminal mitochondrial targeting sequence (MTS) of GLUD2 soon after the duplication event in the hominoid ancestor ∼18–25 million years ago. This MTS substitution arose in parallel with two crucial adaptive amino acid changes in the enzyme and likely contributed to the functional adaptation of GLUD2 to the glutamate metabolism of the hominoid brain and other tissues. We suggest that rapid, selectively driven subcellular adaptation, as exemplified by GLUD2, represents a common route underlying the emergence of new gene functions.

ZPS: visualization of recent adaptive evolution of proteins.

BackgroundDetection of adaptive amino acid changes in proteins under recent short-term selection is of great interest for researchers studying microevolutionary processes in microbial pathogens or any other biological species. However, independent occurrence of such point mutations within genetically diverse haplotypes makes it difficult to detect the selection footprint by using traditional molecular evolutionary analyses. The recently developed Zonal Phylogeny (ZP) has been shown to be a useful analytic tool for identifying the footprints of short-term positive selection. ZP separates protein-encoding genes into evolutionarily long-term (with silent diversity) and short-term (without silent diversity) categories, or zones, followed by statistical analysis to detect signs of positive selection in the short-term zone. However, successful broad application of ZP for analysis of large haplotype datasets requires automation of the relatively labor-intensive computational process.ResultsHere we present Zonal Phylogeny Software (ZPS), an application that describes the distribution of single nucleotide polymorphisms (SNPs) of synonymous (silent) and non-synonymous (replacement) nature along branches of the DNA tree for any given protein-coding gene locus. Based on this information, ZPS separates the protein variant haplotypes with silent variability (Primary zone) from those that have recently evolved from the Primary zone variants by amino acid changes (External zone). Further comparative analysis of mutational hot-spot frequencies and haplotype diversity between the two zones allows determination of whether the External zone haplotypes emerged under positive selection.ConclusionsAs a visualization tool, ZPS depicts the protein tree in a DNA tree, indicating the most parsimonious numbers of synonymous and non-synonymous changes along the branches of a maximum-likelihood based DNA tree, along with information on homoplasy, reversion and structural mutation hot-spots. Through zonal differentiation, ZPS allows detection of recent adaptive evolution via selection of advantageous structural mutations, even when the advantage conferred by such mutations is relatively short-term (as in the case of "source-sink" evolutionary dynamics, which may represent a major mode of virulence evolution in microbes).

Molecular Adaptation during Adaptive Radiation in the Hawaiian Endemic Genus Schiedea

Background“Explosive” adaptive radiations on islands remain one of the most puzzling evolutionary phenomena. The rate of phenotypic and ecological adaptations is extremely fast during such events, suggesting that many genes may be under fairly strong selection. However, no evidence for adaptation at the level of protein coding genes was found, so it has been suggested that selection may work mainly on regulatory elements. Here we report the first evidence that positive selection does operate at the level of protein coding genes during rapid adaptive radiations. We studied molecular adaptation in Hawaiian endemic plant genus Schiedea (Caryophyllaceae), which includes closely related species with a striking range of morphological and ecological forms, varying from rainforest vines to woody shrubs growing in desert-like conditions on cliffs. Given the remarkable difference in photosynthetic performance between Schiedea species from different habitats, we focused on the “photosynthetic” Rubisco enzyme, the efficiency of which is known to be a limiting step in plant photosynthesis.ResultsWe demonstrate that the chloroplast rbcL gene, encoding the large subunit of Rubisco enzyme, evolved under strong positive selection in Schiedea. Adaptive amino acid changes occurred in functionally important regions of Rubisco that interact with Rubisco activase, a chaperone which promotes and maintains the catalytic activity of Rubisco. Interestingly, positive selection acting on the rbcL might have caused favorable cytotypes to spread across several Schiedea species.SignificanceWe report the first evidence for adaptive changes at the DNA and protein sequence level that may have been associated with the evolution of photosynthetic performance and colonization of new habitats during a recent adaptive radiation in an island plant genus. This illustrates how small changes at the molecular level may change ecological species performance and helps us to understand the molecular bases of extremely fast rate of adaptation during island adaptive radiations.

Did brain-specific genes evolve faster in humans than in chimpanzees?

One of the most distinctive characteristics of humans among primates is the size, organization and function of the brain. A recent study has proposed that there was widespread accelerated sequence evolution of genes functioning in the nervous system during human origins. Here we test this hypothesis by a genome-wide analysis of genes that are expressed predominantly or specifically in brain tissues and genes that have important roles in the brain, identified on the basis of five different definitions of brain specificity. Although there is little overlap among the five sets of brain-specific genes, none of them supports human acceleration. On the contrary, some datasets show significantly fewer nonsynonymous substitutions in humans than in chimpanzees for brain-specific genes relative to other genes in the genome. Our results suggest that the unique features of the human brain did not arise by a large number of adaptive amino acid changes in many proteins.

Evolution of Human Longevity and the Genetic Complexity Governing Aging Rate

The genetic complexity of processes governing the aging rate of humans was estimated by determining the maximum rate at which life-span has evolved along the hominid ancestral-descendant sequence. Maximum life-span potential was found to have increased approximately twofold over the past 3 million years, reaching a maximum rate of increase of 14 years per 100,000 years about 100,000 years ago. It is estimated that about 0.6% of the total functional genes have received substitutions leading to one or more adaptive amino acid changes during this 100,000-year time-period. This suggests that aging is not the result of the expression of a large number of independently acting processes. Instead, primary aging processes appear to exist in which only a few genetic changes are necessary to uniformly decrease the aging rate of many different physiological functions. Reproduced by permission. Richard G. Cutler, Evolution of Human Longevity and the Genetic Complexity Governing Aging Rate. Proc. Natl. Acad. Sci. U.S.A. 72 , 4664-4668 (1975).

Evolution of human longevity and the genetic complexity governing aging rate.

Genetic complexity of processes governing the aging rate of man was estimated by determining the maximum rate lifespan has evolved along the hominid ancestral-descendant sequence. Maximum lifespan potential was found to have increased approximately 2-fold over the past 3 million years, reaching a maximum rate of increase of 14 years per 100,000 years about 100,000 years ago. It is estimated that about 0.6% of the total functional genes have received substitutions leading to one or more adaptive amino acid changes during this 100,000-year time-period. This suggests that aging is not the result of an expression of a large number of independently acting processes. Instead, primary aging processes appear to exist where only a few genetic changes are necessary to decrease uniformly the aging rate of many different physiological functions.