The effect of species representation on the detection of positive selection in primate gene data sets

The Diverse Evolutionary Histories of Domesticated Metaviral Capsid Genes in Mammals

Selfish genetic elements comprise significant fractions of mammalian genomes. In rare instances, host genomes domesticate segments of these elements for function. Using a complete human genome assembly and 25 additional vertebrate genomes, we re-analyzed the evolutionary trajectories and functional potential of capsid (CA) genes domesticated from Metaviridae, a lineage of retrovirus-like retrotransposons. Our study expands on previous analyses to unearth several new insights about the evolutionary histories of these ancient genes. We find that at least five independent domestication events occurred from diverse Metaviridae, giving rise to three universally retained single-copy genes evolving under purifying selection and two gene families unique to placental mammals, with multiple members showing evidence of rapid evolution. In the SIRH/RTL family, we find diverse amino-terminal domains, widespread loss of protein-coding capacity in RTL10 despite its retention in several mammalian lineages, and differential utilization of an ancient programmed ribosomal frameshift in RTL3 between the domesticated CA and protease domains. Our analyses also reveal that most members of the PNMA family in mammalian genomes encode a conserved putative amino-terminal RNA-binding domain (RBD) both adjoining and independent from domesticated CA domains. Our analyses lead to a significant correction of previous annotations of the essential CCDC8 gene. We show that this putative RBD is also present in several extant Metaviridae, revealing a novel protein domain configuration in retrotransposons. Collectively, our study reveals the divergent outcomes of multiple domestication events from diverse Metaviridae in the common ancestor of placental mammals.

The diverse evolutionary histories of domesticated metaviral capsid genes in mammals

Selfish genetic elements and their remnants comprise at least half of the human genome. Active transposons duplicate by inserting copies at new sites in a host genome. Following insertion, transposons can acquire mutations that render them inactive; the accrual of additional mutations can render them unrecognizable over time. However, in rare instances, segments of transposons become useful for the host, in a process called gene domestication. Using the first complete human genome assembly and 25 additional vertebrate genomes, we analyzed the evolutionary trajectories and functional potential of genes domesticated from the capsid genes of Metaviridae, a retroviral-like retrotransposon family. Our analysis reveals four families of domesticated capsid genes in placental mammals with varied evolutionary outcomes, ranging from universal retention to lineage-specific duplications or losses and from purifying selection to lineage-specific rapid evolution. The four families of domesticated capsid genes have divergent amino-terminal domains, inherited from four distinct ancestral metaviruses. Structural predictions reveal that many domesticated genes encode a previously unrecognized RNA-binding domain retained in multiple paralogs in mammalian genomes both adjacent to and independent from the capsid domain. Collectively, our study reveals diverse outcomes of domestication of diverse metaviruses, which led to structurally and evolutionarily diverse genes that encode important, but still largely-unknown functions in placental mammals. (207).

Unconventional conservation reveals structure-function relationships in the synaptonemal complex

Functional requirements constrain protein evolution, commonly manifesting in a conserved amino acid sequence. Here, we extend this idea to secondary structural features by tracking their conservation in essential meiotic proteins with highly diverged sequences. The synaptonemal complex (SC) is a ~100-nm-wide ladder-like meiotic structure present in all eukaryotic clades, where it aligns parental chromosomes and regulates exchanges between them. Despite the conserved ultrastructure and functions of the SC, SC proteins are highly divergent within Caenorhabditis. However, SC proteins have highly conserved length and coiled-coil domain structure. We found the same unconventional conservation signature in Drosophila and mammals, and used it to identify a novel SC protein in Pristionchus pacificus, Ppa-SYP-1. Our work suggests that coiled-coils play wide-ranging roles in the structure and function of the SC, and more broadly, that expanding sequence analysis beyond measures of per-site similarity can enhance our understanding of protein evolution and function.

Pervasive positive selection on virus receptors driven by host-virus conflicts in mammals

Viruses hijack cellular proteins known as viral receptors to initiate their infection. Viral receptors are subject to two conflicting directional forces, namely negative selection to maintain their cellular function and positive selection resulted from everchanging host-virus arms race. Much remains unclear how viral receptors evolved in mammals, and whether viral receptors from different mammal groups experienced different strength of natural selection. Here, we perform evolutionary analyses of 92 viral receptors in five major orders of mammals, including Carnivora, Cetartiodactyla, Chiroptera, Primates, and Rodentia. In all the five mammal orders, signals of positive selection are detected for a high proportion of viral receptors (from 41% in Carnivora to 65% in Rodentia). Many positively selected residues overlap host-virus interaction interface. Compared with control genes, we find viral receptors underwent elevated rate of adaptive evolution in all the five mammal orders, suggesting that host-virus conflicts are the main driver of the adaptive evolution of viral receptors in mammals. Interestingly, the overall strength of natural selection acting on viral receptors driven by host-virus arms race is largely homogenous and correlated among different mammal orders with bats and rodents, zoonosis reservoirs of importance, unexceptional. Taken together, our findings indicate host-virus conflicts have driven the elevated rate of adaptive evolution in viral receptors across mammals, and might have important implications in zoonosis surveillance and prediction. Importance Viral receptors are cellular proteins hijacked by viruses to help their infections. A complete picture on the evolution of viral receptors in mammals is still lacking. Here, we perform a comprehensive evolutionary analysis of the evolution of 92 viral receptors in five mammal orders, including Carnivora, Cetartiodactyla, Chiroptera, Primates, and Rodentia. We find that positive selection pervasively occurred during the evolution of viral receptors, and viral receptors exhibit at an elevated rate of adaptive evolution than control genes in all the five mammal orders, suggesting host-virus conflicts are a major driver of the adaptive evolution of viral receptors. Interestingly, the strength of positive selection acting on viral receptors is similar among the five mammal orders. Our study might have important implications in understanding the evolution of host-virus interaction.

Hearing Sensitivity of Primates: Recurrent and Episodic Positive Selection in Hair Cells and Stereocilia Protein-Coding Genes

The large spectrum of hearing sensitivity observed in primates results from the impact of environmental and behavioral pressures to optimize sound perception and localization. Although evidence of positive selection in auditory genes has been detected in mammals including in Hominoids, selection has never been investigated in other primates. We analyzed 123 genes highly expressed in the inner ear of 27 primate species and tested to what extent positive selection may have shaped these genes in the order Primates tree. We combined both site and branch-site tests to obtain a comprehensive picture of the positively selected genes (PSGs) involved in hearing sensitivity, and drew a detailed description of the most affected branches in the tree. We chose a conservative approach, and thus focused on confounding factors potentially affecting PSG signals (alignment, GC-biased gene conversion, duplications, heterogeneous sequencing qualities). Using site tests, we showed that around 12% of these genes are PSGs, an α selection value consistent with average human genome estimates (10-15%). Using branch-site tests, we showed that the primate tree is heterogeneously affected by positive selection, with the black snub-nosed monkey, the bushbaby, and the orangutan, being the most impacted branches. A large proportion of these genes is inclined to shape hair cells and stereocilia, which are involved in the mechanotransduction process, known to influence frequency perception. Adaptive selection, and more specifically recurrent adaptive evolution, could have acted in parallel on a set of genes (ADGRV1, USH2A, PCDH15, PTPRQ, and ATP8A2) involved in stereocilia growth and the whole complex of bundle links connecting them, in species across different habitats, including high altitude and nocturnal environments.

PoSeiDon: a Nextflow pipeline for the detection of evolutionary recombination events and positive selection

SUMMARY

PoSeiDon is an easy-to-use pipeline that helps researchers to find recombination events and sites under positive selection in protein-coding sequences. By entering homologous sequences, PoSeiDon builds an alignment, estimates a best-fitting substitution model and performs a recombination analysis followed by the construction of all corresponding phylogenies. Finally, significantly positive selected sites are detected according to different models for the full alignment and possible recombination fragments. The results of PoSeiDon are summarized in a user-friendly HTML page providing all intermediate results and the graphical representation of recombination events and positively selected sites.

AVAILABILITY AND IMPLEMENTATION

PoSeiDon is freely available at https://github.com/hoelzer/poseidon. The pipeline is implemented in Nextflow with Docker support and processes the output of various tools.

Positive natural selection in primate genes of the type I interferon response

Background

The Type I interferon response is an important first-line defense against viruses. In turn, viruses antagonize (i.e., degrade, mis-localize, etc.) many proteins in interferon pathways. Thus, hosts and viruses are locked in an evolutionary arms race for dominance of the Type I interferon pathway. As a result, many genes in interferon pathways have experienced positive natural selection in favor of new allelic forms that can better recognize viruses or escape viral antagonists. Here, we performed a holistic analysis of selective pressures acting on genes in the Type I interferon family. We initially hypothesized that the genes responsible for inducing the production of interferon would be antagonized more heavily by viruses than genes that are turned on as a result of interferon. Our logic was that viruses would have greater effect if they worked upstream of the production of interferon molecules because, once interferon is produced, hundreds of interferon-stimulated proteins would activate and the virus would need to counteract them one-by-one.

Results

We curated multiple sequence alignments of primate orthologs for 131 genes active in interferon production and signaling (herein, “induction” genes), 100 interferon-stimulated genes, and 100 randomly chosen genes. We analyzed each multiple sequence alignment for the signatures of recurrent positive selection. Counter to our hypothesis, we found the interferon-stimulated genes, and not interferon induction genes, are evolving significantly more rapidly than a random set of genes. Interferon induction genes evolve in a way that is indistinguishable from a matched set of random genes (22% and 18% of genes bear signatures of positive selection, respectively). In contrast, interferon-stimulated genes evolve differently, with 33% of genes evolving under positive selection and containing a significantly higher fraction of codons that have experienced selection for recurrent replacement of the encoded amino acid.

Conclusion

Viruses may antagonize individual products of the interferon response more often than trying to neutralize the system altogether.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01783-z.

Rapid Evolution of HERC6 and Duplication of a Chimeric HERC5/6 Gene in Rodents and Bats Suggest an Overlooked Role of HERCs in Mammalian Immunity

Studying the evolutionary diversification of mammalian antiviral defenses is of main importance to better understand our innate immune repertoire. The small HERC proteins are part of a multigene family, including HERC5 and HERC6, which have probably diversified through complex evolutionary history in mammals. Here, we performed mammalian-wide phylogenetic and genomic analyses of HERC5 and HERC6, using 83 orthologous sequences from bats, rodents, primates, artiodactyls, and carnivores—the top five representative groups of mammalian evolution. We found that HERC5 has been under weak and differential positive selection in mammals, with only primate HERC5 showing evidences of pathogen-driven selection. In contrast, HERC6 has been under strong and recurrent adaptive evolution in mammals, suggesting past and widespread genetic arms-races with viral pathogens. Importantly, the rapid evolution of mammalian HERC6 spacer domain suggests that it might be a host-pathogen interface, targeting viral proteins and/or being the target of virus antagonists. Finally, we identified a HERC5/6 chimeric gene that arose from independent duplication in rodent and bat lineages and encodes for a conserved HERC5 N-terminal domain and divergent HERC6 spacer and HECT domains. This duplicated chimeric gene highlights adaptations that potentially contribute to rodent and bat immunity. Our findings open new research avenues on the functions of HERC6 and HERC5/6 in mammals, and on their implication in antiviral innate immunity.

A Burst of Genetic Innovation in Drosophila Actin-Related Proteins for Testis-Specific Function

Many cytoskeletal proteins perform fundamental biological processes and are evolutionarily ancient. For example, the superfamily of actin-related proteins (Arps) specialized early in eukaryotic evolution for diverse cellular roles in the cytoplasm and the nucleus. Despite its strict conservation across eukaryotes, we find that the Arp superfamily has undergone dramatic lineage-specific diversification in Drosophila. Our phylogenomic analyses reveal four independent Arp gene duplications that occurred in the common ancestor of the obscura group of Drosophila and have been mostly preserved in this lineage. All four obscura-specific Arp paralogs are predominantly expressed in the male germline and have evolved under positive selection. We focus our analyses on the divergent Arp2D paralog, which arose via a retroduplication event from Arp2, a component of the Arp2/3 complex that polymerizes branched actin networks. Computational modeling analyses suggest that Arp2D can replace Arp2 in the Arp2/3 complex and bind actin monomers. Together with the signature of positive selection, our findings suggest that Arp2D may augment Arp2's functions in the male germline. Indeed, we find that Arp2D is expressed during and following male meiosis, where it localizes to distinct locations such as actin cones-specialized cytoskeletal structures that separate bundled spermatids into individual mature sperm. We hypothesize that this unprecedented burst of genetic innovation in cytoskeletal proteins may have been driven by the evolution of sperm heteromorphism in the obscura group of Drosophila.

DGINN, an automated and highly-flexible pipeline for the detection of genetic innovations on protein-coding genes

Adaptive evolution has shaped major biological processes. Finding the protein-coding genes and the sites that have been subjected to adaptation during evolutionary time is a major endeavor. However, very few methods fully automate the identification of positively selected genes, and widespread sources of genetic innovations such as gene duplication and recombination are absent from most pipelines. Here, we developed DGINN, a highly-flexible and public pipeline to Detect Genetic INNovations and adaptive evolution in protein-coding genes. DGINN automates, from a gene's sequence, all steps of the evolutionary analyses necessary to detect the aforementioned innovations, including the search for homologs in databases, assignation of orthology groups, identification of duplication and recombination events, as well as detection of positive selection using five methods to increase precision and ranking of genes when a large panel is analyzed. DGINN was validated on nineteen genes with previously-characterized evolutionary histories in primates, including some engaged in host-pathogen arms-races. Our results confirm and also expand results from the literature, including novel findings on the Guanylate-binding protein family, GBPs. This establishes DGINN as an efficient tool to automatically detect genetic innovations and adaptive evolution in diverse datasets, from the user's gene of interest to a large gene list in any species range.

SARS-CoV-2, an evolutionary perspective of interaction with human ACE2 reveals undiscovered amino acids necessary for complex stability

The emergence of SARS-CoV-2 has resulted in nearly 1,280,000 infections and 73,000 deaths globally so far. This novel virus acquired the ability to infect human cells using the SARS-CoV cell receptor hACE2. Because of this, it is essential to improve our understanding of the evolutionary dynamics surrounding the SARS-CoV-2 hACE2 interaction. One way theory predicts selection pressures should shape viral evolution is to enhance binding with host cells. We first assessed evolutionary dynamics in select betacoronavirus spike protein genes to predict whether these genomic regions are under directional or purifying selection between divergent viral lineages, at various scales of relatedness. With this analysis, we determine a region inside the receptor-binding domain with putative sites under positive selection interspersed among highly conserved sites, which are implicated in structural stability of the viral spike protein and its union with human receptor ACE2. Next, to gain further insights into factors associated with recognition of the human host receptor, we performed modeling studies of five different betacoronaviruses and their potential binding to hACE2. Modeling results indicate that interfering with the salt bridges at hot spot 353 could be an effective strategy for inhibiting binding, and hence for the prevention of SARS-CoV-2 infections. We also propose that a glycine residue at the receptor-binding domain of the spike glycoprotein can have a critical role in permitting bat SARS-related coronaviruses to infect human cells.

Positive selection in dNTPase SAMHD1 throughout mammalian evolution

The vertebrate protein SAMHD1 is highly unusual in having roles in cellular metabolic regulation, antiviral restriction, and regulation of innate immunity. Its deoxynucleoside triphosphohydrolase activity regulates cellular dNTP concentration, reducing levels below those required by lentiviruses and other viruses to replicate. To counter this threat, some primate lentiviruses encode accessory proteins that bind SAMHD1 and induce its degradation; in turn, positive diversifying selection has been observed in regions bound by these lentiviral proteins, suggesting that primate SAMHD1 has coevolved to evade these countermeasures. Moreover, deleterious polymorphisms in human SAMHD1 are associated with autoimmune disease linked to uncontrolled DNA synthesis of endogenous retroelements. Little is known about how evolutionary pressures affect these different SAMHD1 functions. Here, we examine the deeper history of these interactions by testing whether evolutionary signatures in SAMHD1 extend to other mammalian groups and exploring the molecular basis of this coevolution. Using codon-based likelihood models, we find positive selection in SAMHD1 within each mammal lineage for which sequence data are available. We observe positive selection at sites clustered around T592, a residue that is phosphorylated to regulate SAMHD1 activity. We verify experimentally that mutations within this cluster affect catalytic rate and lentiviral restriction, suggesting that virus-host coevolution has required adaptations of enzymatic function. Thus, persistent positive selection may have involved the adaptation of SAMHD1 regulation to balance antiviral, metabolic, and innate immunity functions.

Biased Inference of Selection Due to GC-Biased Gene Conversion and the Rate of Protein Evolution in Flycatchers When Accounting for It

The rate of recombination impacts on rates of protein evolution for at least two reasons: it affects the efficacy of selection due to linkage and influences sequence evolution through the process of GC-biased gene conversion (gBGC). We studied how recombination, via gBGC, affects inferences of selection in gene sequences using comparative genomic and population genomic data from the collared flycatcher (Ficedula albicollis). We separately analyzed different mutation categories (“strong”-to-“weak,” “weak-to-strong,” and GC-conservative changes) and found that gBGC impacts on the distribution of fitness effects of new mutations, and leads to that the rate of adaptive evolution and the proportion of adaptive mutations among nonsynonymous substitutions are underestimated by 22–33%. It also biases inferences of demographic history based on the site frequency spectrum. In light of this impact, we suggest that inferences of selection (and demography) in lineages with pronounced gBGC should be based on GC-conservative changes only. Doing so, we estimate that 10% of nonsynonymous mutations are effectively neutral and that 27% of nonsynonymous substitutions have been fixed by positive selection in the flycatcher lineage. We also find that gene expression level, sex-bias in expression, and the number of protein–protein interactions, but not Hill–Robertson interference (HRI), are strong determinants of selective constraint and rate of adaptation of collared flycatcher genes. This study therefore illustrates the importance of disentangling the effects of different evolutionary forces and genetic factors in interpretation of sequence data, and from that infer the role of natural selection in DNA sequence evolution.

A metaanalysis of bat phylogenetics and positive selection based on genomes and transcriptomes from 18 species

This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation.

Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

Derivation of simian tropic HIV-1 infectious clone reveals virus adaptation to a new host

To replicate in a new host, lentiviruses must adapt to exploit required host factors and evade species-specific antiviral proteins. Understanding how host protein variation drives lentivirus adaptation allowed us to expand the host range of HIV-1 to pigtail macaques. We have previously derived a viral swarm (in the blood of infected animals) that can cause AIDS in this new host. To further exploit this reagent, we generated infectious molecular clones (IMCs) from the viral swarm. We identified clones with high replicative capacity in pigtail peripheral blood mononuclear cells (PBMC) in vitro and used in vivo replication to select an individual IMC, named stHIV-A19 (for simian tropic HIV-1 clone A19), which recapitulated the phenotype obtained with the viral swarm. Adaptation of HIV-1 in macaques led to the acquisition of amino acid changes in viral proteins, such as capsid (CA), that are rarely seen in HIV-1-infected humans. Using stHIV-A19, we show that these CA changes confer a partial resistance to the host cell inhibitor Mx2 from pigtail macaques, but that complete resistance is associated with a fitness defect. Adaptation of HIV-1 to a new host will lead to a more accurate animal model and a better understanding of virus-host interactions.

Evolution of Hepatitis B Virus Receptor NTCP Reveals Differential Pathogenicities and Species Specificities of Hepadnaviruses in Primates, Rodents, and Bats

Human hepatitis B virus (HBV) is a global health problem, affecting more than 250 million people worldwide. HBV-like viruses, named orthohepadnaviruses, also naturally infect nonhuman primates, rodents, and bats, but their pathogenicity and evolutionary history are unclear. Here, we determined the evolutionary history of the HBV receptors NTCP and GPC5 over millions of years of primate, rodent, and bat evolution. We use this as a proxy to understand the pathogenicity of orthohepadnaviruses in mammalian hosts and to determine the implications for species specificity. We found that NTCP, but not GPC5, has evolved under positive selection in primates (27 species), rodents (18 species), and bats (21 species) although at distinct residues. Notably, the positively selected codons map to the HBV-binding sites in primate NTCP, suggesting past genetic "arms races" with pathogenic orthohepadnaviruses. In rodents, the positively selected codons fall outside and within the presumed HBV-binding sites, which may contribute to the restricted circulation of rodent orthohepadnaviruses. In contrast, the presumed HBV-binding motifs in bat NTCP are conserved, and none of the positively selected codons map to this region. This suggests that orthohepadnaviruses may bind to different surfaces in bat NTCP. Alternatively, the patterns may reflect adaptive changes associated with metabolism rather than pathogens. Overall, our findings further point to NTCP as a naturally occurring genetic barrier for cross-species transmissions in primates, which may contribute to the narrow host range of HBV. In contrast, this constraint seems less important in bats, which may correspond to greater orthohepadnavirus circulation and diversity.IMPORTANCE Chronic infection with hepatitis B virus (HBV) is a major cause of liver disease and cancer in humans. Mammalian HBV-like viruses are also found in nonhuman primates, rodents, and bats. As for most viruses, HBV requires a successful interaction with a host receptor for replication. Cellular receptors are thus key determinants of host susceptibility as well as specificity. One hallmark of pathogenic virus-host relationships is the reciprocal evolution of host receptor and viral envelope proteins, as a result of their antagonistic interaction over time. The dynamics of these so-called "evolutionary arms races" can leave signatures of adaptive selection, which in turn reveal the evolutionary history of the virus-host interaction as well as viral pathogenicity and the genetic determinants of species specificity. Here, we show how HBV-like viruses have shaped the evolutionary history of their mammalian host receptor, as a result of their ancient pathogenicity, and decipher the genetic determinants of cross-species transmissions.

Genetic conflicts with Plasmodium parasites and functional constraints shape the evolution of erythrocyte cytoskeletal proteins

Plasmodium parasites exerted a strong selective pressure on primate genomes and mutations in genes encoding erythrocyte cytoskeleton proteins (ECP) determine protective effects against Plasmodium infection/pathogenesis. We thus hypothesized that ECP-encoding genes have evolved in response to Plasmodium-driven selection. We analyzed the evolutionary history of 15 ECP-encoding genes in primates, as well as of their Plasmodium-encoded ligands (KAHRP, MESA and EMP3). Results indicated that EPB42, SLC4A1, and SPTA1 evolved under pervasive positive selection and that episodes of positive selection tended to occur more frequently in primate species that host a larger number of Plasmodium parasites. Conversely, several genes, including ANK1 and SPTB, displayed extensive signatures of purifying selection in primate phylogenies, Homininae lineages, and human populations, suggesting strong functional constraints. Analysis of Plasmodium genes indicated adaptive evolution in MESA and KAHRP; in the latter, different positively selected sites were located in the spectrin-binding domains. Because most of the positively selected sites in alpha-spectrin localized to the domains involved in the interaction with KAHRP, we suggest that the two proteins are engaged in an arms-race scenario. This observation is relevant because KAHRP is essential for the formation of “knobs”, which represent a major virulence determinant for P. falciparum.

Smc5/6 Antagonism by HBx Is an Evolutionarily Conserved Function of Hepatitis B Virus Infection in Mammals

Chronic infection with hepatitis B virus (HBV) is a major cause of liver disease and cancer in humans. HBVs (family Hepadnaviridae) have been associated with mammals for millions of years. Recently, the Smc5/6 complex, known for its essential housekeeping functions in genome maintenance, was identified as an antiviral restriction factor of human HBV. The virus has, however, evolved to counteract this defense mechanism by degrading the complex via its regulatory HBx protein. Whether the antiviral activity of the Smc5/6 complex against hepadnaviruses is an important and evolutionarily conserved function is unknown. In this study, we used an evolutionary and functional approach to address this question. We first performed phylogenetic and positive selection analyses of the Smc5/6 complex subunits and found that they have been conserved in primates and mammals. Yet, Smc6 showed marks of adaptive evolution, potentially reminiscent of a virus-host "arms race." We then functionally tested the HBx proteins from six divergent hepadnaviruses naturally infecting primates, rodents, and bats. We demonstrate that despite little sequence homology, these HBx proteins efficiently degraded mammalian Smc5/6 complexes, independently of the host species and of the sites under positive selection. Importantly, all HBx proteins also rescued the replication of an HBx-deficient HBV in primary human hepatocytes. These findings point to an evolutionarily conserved requirement for Smc5/6 inactivation by HBx, showing that Smc5/6 antiviral activity has been an important defense mechanism against hepadnaviruses in mammals. It will be interesting to investigate whether Smc5/6 may further be a restriction factor of other, yet-unidentified viruses that may have driven some of its adaptation.IMPORTANCE Infection with hepatitis B virus (HBV) led to 887,000 human deaths in 2015. HBV has been coevolving with mammals for millions of years. Recently, the Smc5/6 complex, which has essential housekeeping functions, was identified as a restriction factor of human HBV antagonized by the regulatory HBx protein. Here we address whether the antiviral activity of Smc5/6 is an important evolutionarily conserved function. We found that all six subunits of Smc5/6 have been conserved in primates, with only Smc6 showing signatures of an "evolutionary arms race." Using evolution-guided functional analyses that included infections of primary human hepatocytes, we demonstrated that HBx proteins from very divergent mammalian HBVs could all efficiently antagonize Smc5/6, independently of the host species and sites under positive selection. These findings show that Smc5/6 antiviral activity against HBV is an important function in mammals. They also raise the intriguing possibility that Smc5/6 may restrict other, yet-unidentified viruses.

Evolution and Antiviral Specificities of Interferon-Induced Mx Proteins of Bats against Ebola, Influenza, and Other RNA Viruses

Bats serve as a reservoir for various, often zoonotic viruses, including significant human pathogens such as Ebola and influenza viruses. However, for unknown reasons, viral infections rarely cause clinical symptoms in bats. A tight control of viral replication by the host innate immune defense might contribute to this phenomenon. Transcriptomic studies revealed the presence of the interferon-induced antiviral myxovirus resistance (Mx) proteins in bats, but detailed functional aspects have not been assessed. To provide evidence that bat Mx proteins might act as key factors to control viral replication we cloned Mx1 cDNAs from three bat families, Pteropodidae, Phyllostomidae, and Vespertilionidae. Phylogenetically these bat Mx1 genes cluster closely with their human ortholog MxA. Using transfected cell cultures, minireplicon systems, virus-like particles, and virus infections, we determined the antiviral potential of the bat Mx1 proteins. Bat Mx1 significantly reduced the polymerase activity of viruses circulating in bats, including Ebola and influenza A-like viruses. The related Thogoto virus, however, which is not known to infect bats, was not inhibited by bat Mx1. Further, we provide evidence for positive selection in bat Mx1 genes that might explain species-specific antiviral activities of these proteins. Together, our data suggest a role for Mx1 in controlling these viruses in their bat hosts.IMPORTANCE Bats are a natural reservoir for various viruses that rarely cause clinical symptoms in bats but are dangerous zoonotic pathogens, like Ebola or rabies virus. It has been hypothesized that the interferon system might play a key role in controlling viral replication in bats. We speculate that the interferon-induced Mx proteins might be key antiviral factors of bats and have coevolved with bat-borne viruses. This study evaluated for the first time a large set of bat Mx1 proteins spanning three major bat families for their antiviral potential, including activity against Ebola virus and bat influenza A-like virus, and we describe here their phylogenetic relationship, revealing patterns of positive selection that suggest a coevolution with viral pathogens. By understanding the molecular mechanisms of the innate resistance of bats against viral diseases, we might gain important insights into how to prevent and fight human zoonotic infections caused by bat-borne viruses.

Functional Divergence in Teleost Cardiac Troponin Paralogs Guides Variation in the Interaction of TnI Switch Region with TnC

Gene duplication results in extra copies of genes that must coevolve with their interacting partners in multimeric protein complexes. The cardiac troponin (Tn) complex, containing TnC, TnI, and TnT, forms a distinct functional unit critical for the regulation of cardiac muscle contraction. In teleost fish, the function of the Tn complex is modified by the consequences of differential expression of paralogs in response to environmental thermal challenges. In this article, we focus on the interaction between TnI and TnC, coded for by genes that have independent evolutionary origins, but the co-operation of their protein products has necessitated coevolution. In this study, we characterize functional divergence of TnC and TnI paralogs, specifically the interrelated roles of regulatory subfunctionalization and structural subfunctionalization. We determined that differential paralog transcript expression in response to temperature acclimation results in three combinations of TnC and TnI in the zebrafish heart: TnC1a/TnI1.1, TnC1b/TnI1.1, and TnC1a/TnI1.5. Phylogenetic analysis of these highly conserved proteins identified functionally divergent residues in TnI and TnC. The structural and functional effect of these Tn combinations was modeled with molecular dynamics simulation to link divergent sites to changes in interaction strength. Functional divergence in TnI and TnC were not limited to the residues involved with TnC/TnI switch interaction, which emphasizes the complex nature of Tn function. Patterns in domain-specific divergent selection and interaction energies suggest that substitutions in the TnI switch region are crucial to modifying TnI/TnC function to maintain cardiac contraction with temperature changes. This integrative approach introduces Tn as a model of functional divergence that guides the coevolution of interacting proteins.

An Evolutionary View of the Arms Race between Protein Kinase R and Large DNA Viruses

To establish productive infections, viruses must counteract numerous cellular defenses that are poised to recognize viruses as nonself and to activate antiviral pathways. The opposing goals of host and viral factors lead to evolutionary arms races that can be illuminated by evolutionary and computational methods and tested in experimental models. Here we illustrate how this perspective has been contributing to our understanding of the interactions of the protein kinase R pathway with large DNA viruses.

Evolutionary Analyses Suggest a Function of MxB Immunity Proteins Beyond Lentivirus Restriction

Viruses impose diverse and dynamic challenges on host defenses. Diversifying selection of codons and gene copy number variation are two hallmarks of genetic innovation in antiviral genes engaged in host-virus genetic conflicts. The myxovirus resistance (Mx) genes encode interferon-inducible GTPases that constitute a major arm of the cell-autonomous defense against viral infection. Unlike the broad antiviral activity of MxA, primate MxB was recently shown to specifically inhibit lentiviruses including HIV-1. We carried out detailed evolutionary analyses to investigate whether genetic conflict with lentiviruses has shaped MxB evolution in primates. We found strong evidence for diversifying selection in the MxB N-terminal tail, which contains molecular determinants of MxB anti-lentivirus specificity. However, we found no overlap between previously-mapped residues that dictate lentiviral restriction and those that have evolved under diversifying selection. Instead, our findings are consistent with MxB having a long-standing and important role in the interferon response to viral infection against a broader range of pathogens than is currently appreciated. Despite its critical role in host innate immunity, we also uncovered multiple functional losses of MxB during mammalian evolution, either by pseudogenization or by gene conversion from MxA genes. Thus, although the majority of mammalian genomes encode two Mx genes, this apparent stasis masks the dramatic effects that recombination and diversifying selection have played in shaping the evolutionary history of Mx genes. Discrepancies between our study and previous publications highlight the need to account for recombination in analyses of positive selection, as well as the importance of using sequence datasets with appropriate depth of divergence. Our study also illustrates that evolutionary analyses of antiviral gene families are critical towards understanding molecular principles that govern host-virus interactions and species-specific susceptibility to viral infection.

Evolutionary analyses have the potential to reveal not only biochemical details about host-virus arms-races but also the nature of the pathogens that drove them. Primate MxB was recently shown to restrict the replication of primate lentiviruses, including HIV-1. However, we find that positive selection in primate MxB is incongruent with known molecular determinants of lentiviral restriction. This suggests that MxB has antiviral activity against a broader range of viruses than is currently appreciated. We also identified multiple losses of MxB in mammals, as well as rampant recombination between Mx paralogs, which has distorted gene orthology. Our study illustrates the importance of evolution-guided functional analyses of antiviral gene families.

Positive selection underlies the species-specific binding of Plasmodium falciparum RH5 to human basigin

Plasmodium falciparum, the causative agent of the deadliest form of malaria, is a member of the Laverania subgenus, which includes ape-infecting parasites. P. falciparum is thought to have originated in gorillas, although infection is now restricted to humans. Laverania parasites display remarkable host-specificity, which is partially mediated by the interaction between parasite ligands and host receptors. We analyse the evolution of BSG (basigin) and GYPA (glycophorin A) in primates/hominins, as well as of their Plasmodium-encoded ligands, PfRH5 and PfEBA175. We show that, in primates, positive selection targeted two sites in BSG (F27 and H102), both involved in PfRH5 binding. A population genetics-phylogenetics approach detected the strongest selection for the gorilla lineage: one of the positively selected sites (K191) is a major determinant of PfRH5 binding affinity. Analysis of RH5 genes indicated episodic selection on the P. falciparum branch; the positively selected W447 site is known to stabilize the interaction with human basigin. Conversely, we detect no selection in the receptor-binding region of EBA175 in the P. falciparum lineage. Its host receptor, GYPA, shows evidence of positive selection in all hominid lineages; selected codons include glycosylation sites that modulate PfEBA175 binding affinity. Data herein provide an evolutionary explanation for species-specific binding of the PfRH5-BSG ligand-receptor pair and support the hypothesis that positive selection at these genes drove the host shift leading to the emergence of P. falciparum as a human pathogen.