FUBAR: a fast, unconstrained bayesian approximation for inferring selection

Complete sequences of 4 viral hemorrhagic septicemia virus IVb isolates and their virulence in northern pike fry

Four viral hemorrhagic septicemia virus (VHSV) genotype IVb isolates were sequenced, their genetic variation explored, and comparative virulence assayed with experimental infections of northern pike Esox lucius fry. In addition to the type strain MI03, the complete 11183 bp genome of the first round goby Neogobius melanostomus isolate from the St. Lawrence River, and the 2013 and 2014 isolates from gizzard shad Dorosoma cepedianum die-offs in Irondequoit Bay, Lake Ontario and Dunkirk Harbor, Lake Erie were all deep sequenced on an Illumina platform. Mutations documented in the 11 yr since the MI03 index case from Lake St. Clair muskellunge Esox masquinongy showed 87 polymorphisms among the 4 isolates. Twenty-six mutations were non-synonymous and located at 18 different positions within the matrix protein, glycoprotein, non-virion protein, and RNA polymerase genes. The same 4 isolates were used to infect northern pike fry by a single 1 h bath exposure. Cumulative percent mortality varied from 42.5 to 62.5%. VHSV was detected in 57% (41/72) of the survivors at the end of the 21-d trial, suggesting that the virus was not rapidly cleared. Lesions were observed in many of the moribund and dead northern pike, such as hemorrhaging in the skin and fins, as well as hydrocephalus. Mean viral load measured from the trunk and visceral tissues of MI03-infected pike was significantly higher than the quantities detected in fish infected with the most recent isolates of genotype IVb, but there were no differences in cumulative mortality observed.

Insights into the Impact of CD8+ Immune Modulation on Human Immunodeficiency Virus Evolutionary Dynamics in Distinct Anatomical Compartments by Using Simian Immunodeficiency Virus-Infected Macaque Models of AIDS Progression

A thorough understanding of the role of human immunodeficiency virus (HIV) intrahost evolution in AIDS pathogenesis has been limited by the need for longitudinally sampled viral sequences from the vast target space within the host, which are often difficult to obtain from human subjects. CD8⁺ lymphocyte-depleted macaques infected with simian immunodeficiency virus (SIV) provide an increasingly utilized model of pathogenesis due to clinical manifestations similar to those for HIV-1 infection and AIDS progression, as well as a characteristic rapid disease onset. Comparison of this model with SIV-infected non-CD8⁺ lymphocyte-depleted macaques also provides a unique opportunity to investigate the role of CD8⁺ cells in viral evolution and population dynamics throughout the duration of infection. Using several different phylogenetic methods, we analyzed viral gp120 sequences obtained from extensive longitudinal sampling of multiple tissues and enriched leukocyte populations from SIVmac251-infected macaques with or without CD8⁺ lymphocyte depletion. SIV evolutionary and selection patterns in non-CD8⁺ lymphocyte-depleted animals were characterized by sequential population turnover and continual viral adaptation, a scenario readily comparable to intrahost evolutionary patterns during human HIV infection in the absence of antiretroviral therapy. Alternatively, animals that were depleted of CD8⁺ lymphocytes exhibited greater variation in population dynamics among tissues and cell populations over the course of infection. Our findings highlight the major role for CD8⁺ lymphocytes in prolonging disease progression through continual control of SIV subpopulations from various anatomical compartments and the potential for greater independent viral evolutionary behavior among these compartments in response to immune modulation.IMPORTANCE Although developments in combined antiretroviral therapy (cART) strategies have successfully prolonged the time to AIDS onset in HIV-1-infected individuals, a functional cure has yet to be found. Improvement of drug interventions for a virus that is able to infect a wide range of tissues and cell types requires a thorough understanding of viral adaptation and infection dynamics within this target milieu. Although it is difficult to accomplish in the human host, longitudinal sampling of multiple anatomical locations is readily accessible in the SIV-infected macaque models of neuro-AIDS. The significance of our research is in identifying the impact of immune modulation, through differing immune selective pressures, on viral evolutionary behavior in a multitude of anatomical compartments. The results provide evidence encouraging the development of a more sophisticated model that considers a network of individual viral subpopulations within the host, with differing infection and transmission dynamics, which is necessary for more effective treatment strategies.

Biological and phylogenetic characteristics of West African lineages of West Nile virus

The West Nile virus (WNV), isolated in 1937, is an arbovirus (arthropod-borne virus) that infects thousands of people each year. Despite its burden on global health, little is known about the virus’ biological and evolutionary dynamics. As several lineages are endemic in West Africa, we obtained the complete polyprotein sequence from three isolates from the early 1990s, each representing a different lineage. We then investigated differences in growth behavior and pathogenicity for four distinct West African lineages in arthropod (Ap61) and primate (Vero) cell lines, and in mice. We found that genetic differences, as well as viral-host interactions, could play a role in the biological properties in different WNV isolates in vitro, such as: (i) genome replication, (ii) protein translation, (iii) particle release, and (iv) virulence. Our findings demonstrate the endemic diversity of West African WNV strains and support future investigations into (i) the nature of WNV emergence, (ii) neurological tropism, and (iii) host adaptation.

The West Nile virus (WNV) can cause severe neurological diseases including meningitis, encephalitis, and acute flaccid paralysis. Differences in WNV genetics could play a role in the frequency of neurological symptoms from an infection. For the first time, we observed how geographically similar but genetically distinct lineages grow in cellular environments that agree with the transmission chain of West Nile virus—vertebrate-arthropod-vertebrate. We were able to connect our in vitro and in vivo results with relevant epidemiological and molecular data. Our findings highlight the existence of West African lineages with higher virulence and replicative efficiency in vitro and in vivo compared to lineages similar to circulating strains in the United States and Europe. Our investigation of four West African lineages of West Nile virus will help us better understand the biology of the virus and assess future epidemiological threats.

The complete chloroplast genome sequence of the CAM epiphyte Spanish moss (Tillandsia usneoides, Bromeliaceae) and its comparative analysis

Spanish moss (Tillandsia usneoides) is an epiphytic bromeliad widely distributed throughout tropical and warm temperate America. This plant is highly adapted to extreme environmental conditions. Striking features of this species include specialized trichomes (scales) covering the surface of its shoots aiding the absorption of water and nutrients directly from the atmosphere and a specific photosynthesis using crassulacean acid metabolism (CAM). Here we report the plastid genome of Spanish moss and present the comparison of genome organization and sequence evolution within Poales. The plastome of Spanish moss has a quadripartite structure consisting of a large single copy (LSC, 87,439 bp), two inverted regions (IRa and IRb, 26,803 bp) and short single copy (SSC, 18,612 bp) region. The plastid genome had 37.2% GC content and 134 genes with 88 being unique protein-coding genes and 20 of these are duplicated in the IR, similar to other reported bromeliads. Our study shows that early diverging lineages of Poales do not have high substitution rates as compared to grasses, and plastid genomes of bromeliads show structural features considered to be ancestral in graminids. These include the loss of the introns in the clpP and rpoC1 genes and the complete loss or partial degradation of accD and ycf genes in the Graminid clade. Further structural rearrangements appeared in the graminids lacking in Spanish moss, which include a 28-kb inversion between the trnG-UCC-rps14 region and 6-kb in the trnG-UCC-psbD, followed by a third <1kb inversion in the trnT sequence.

Host Genetic Variation Does Not Determine Spatio-Temporal Patterns of European Bat 1 Lyssavirus

The majority of bat rabies cases in Europe are attributed to European bat 1 lyssavirus (EBLV-1), circulating mainly in serotine bats (Eptesicus serotinus). Two subtypes have been defined (EBLV-1a and EBLV-1b), each associated with a different geographical distribution. In this study, we undertake a comprehensive sequence analysis based on 80 newly obtained EBLV-1 nearly complete genome sequences from nine European countries over a 45-year period to infer selection pressures, rates of nucleotide substitution, and evolutionary time scale of these two subtypes in Europe. Our results suggest that the current lineage of EBLV-1 arose in Europe ∼600 years ago and the virus has evolved at an estimated average substitution rate of ∼4.19×10-5 subs/site/year, which is among the lowest recorded for RNA viruses. In parallel, we investigate the genetic structure of French serotine bats at both the nuclear and mitochondrial level and find that they constitute a single genetic cluster. Furthermore, Mantel tests based on interindividual distances reveal the absence of correlation between genetic distances estimated between viruses and between host individuals. Taken together, this indicates that the genetic diversity observed in our E. serotinus samples does not account for EBLV-1a and -1b segregation and dispersal in Europe.

The Influence of HIV on the Evolution of Mycobacterium tuberculosis

HIV significantly affects the immunological environment during tuberculosis coinfection, and therefore may influence the selective landscape upon which M. tuberculosis evolves. To test this hypothesis whole genome sequences were determined for 169 South African M. tuberculosis strains from HIV-1 coinfected and uninfected individuals and analyzed using two Bayesian codon-model based selection analysis approaches: FUBAR which was used to detect persistent positive and negative selection (selection respectively favoring and disfavoring nonsynonymous substitutions); and MEDS which was used to detect episodic directional selection specifically favoring nonsynonymous substitutions within HIV-1 infected individuals. Among the 25,251 polymorphic codon sites analyzed, FUBAR revealed that 189-fold more were detectably evolving under persistent negative selection than were evolving under persistent positive selection. Three specific codon sites within the genes celA2b, katG, and cyp138 were identified by MEDS as displaying significant evidence of evolving under directional selection influenced by HIV-1 coinfection. All three genes encode proteins that may indirectly interact with human proteins that, in turn, interact functionally with HIV proteins. Unexpectedly, epitope encoding regions were enriched for sites displaying weak evidence of directional selection influenced by HIV-1. Although the low degree of genetic diversity observed in our M. tuberculosis data set means that these results should be interpreted carefully, the effects of HIV-1 on epitope evolution in M. tuberculosis may have implications for the design of M. tuberculosis vaccines that are intended for use in populations with high HIV-1 infection rates.

No Substantial Evidence for Sexual Transmission of Minority HIV Drug Resistance Mutations in Men Who Have Sex with Men

During primary HIV infection, the presence of minority drug resistance mutations (DRM) may be a consequence of sexual transmission, de novo mutations, or technical errors in identification. Baseline blood samples were collected from 24 HIV-infected antiretroviral-naive, genetically and epidemiologically linked source and recipient partners shortly after the recipient's estimated date of infection. An additional 32 longitudinal samples were available from 11 recipients. Deep sequencing of HIV reverse transcriptase (RT) was performed (Roche/454), and the sequences were screened for nucleoside and nonnucleoside RT inhibitor DRM. The likelihood of sexual transmission and persistence of DRM was assessed using Bayesian-based statistical modeling. While the majority of DRM (>20%) were consistently transmitted from source to recipient, the probability of detecting a minority DRM in the recipient was not increased when the same minority DRM was detected in the source (Bayes factor [BF] = 6.37). Longitudinal analyses revealed an exponential decay of DRM (BF = 0.05) while genetic diversity increased. Our analysis revealed no substantial evidence for sexual transmission of minority DRM (BF = 0.02). The presence of minority DRM during early infection, followed by a rapid decay, is consistent with the "mutation-selection balance" hypothesis, in which deleterious mutations are more efficiently purged later during HIV infection when the larger effective population size allows more efficient selection. Future studies using more recent sequencing technologies that are less prone to single-base errors should confirm these results by applying a similar Bayesian framework in other clinical settings.IMPORTANCE The advent of sensitive sequencing platforms has led to an increased identification of minority drug resistance mutations (DRM), including among antiretroviral therapy-naive HIV-infected individuals. While transmission of DRM may impact future therapy options for newly infected individuals, the clinical significance of the detection of minority DRM remains controversial. In the present study, we applied deep-sequencing techniques within a Bayesian hierarchical framework to a cohort of 24 transmission pairs to investigate whether minority DRM detected shortly after transmission were the consequence of (i) sexual transmission from the source, (ii) de novo emergence shortly after infection followed by viral selection and evolution, or (iii) technical errors/limitations of deep-sequencing methods. We found no clear evidence to support the sexual transmission of minority resistant variants, and our results suggested that minor resistant variants may emerge de novo shortly after transmission, when the small effective population size limits efficient purge by natural selection.

Ancient Evolution and Dispersion of Human Papillomavirus 58 Variants

Human papillomavirus 58 (HPV58) is found in 10 to 18% of cervical cancers in East Asia but is rather uncommon elsewhere. The distribution and oncogenic potential of HPV58 variants appear to be heterogeneous, since the E7 T20I/G63S variant is more prevalent in East Asia and confers a 7- to 9-fold-higher risk of cervical precancer and cancer. However, the underlying genomic mechanisms that explain the geographic and carcinogenic diversity of HPV58 variants are still poorly understood. In this study, we used a combination of phylogenetic analyses and bioinformatics to investigate the deep evolutionary history of HPV58 complete genome variants. The initial splitting of HPV58 variants was estimated to occur 478,600 years ago (95% highest posterior density [HPD], 391,000 to 569,600 years ago). This divergence time is well within the era of speciation between Homo sapiens and Neanderthals/Denisovans and around three times longer than the modern Homo sapiens divergence times. The expansion of present-day variants in Eurasia could be the consequence of viral transmission from Neanderthals/Denisovans to non-African modern human populations through gene flow. A whole-genome sequence signature analysis identified 3 amino acid changes, 16 synonymous nucleotide changes, and a 12-bp insertion strongly associated with the E7 T20I/G63S variant that represents the A3 sublineage and carries higher carcinogenetic potential. Compared with the capsid proteins, the oncogenes E7 and E6 had increased substitution rates indicative of higher selection pressure. These data provide a comprehensive evolutionary history and genomic basis of HPV58 variants to assist further investigation of carcinogenic association and the development of diagnostic and therapeutic strategies.IMPORTANCE Papillomaviruses (PVs) are an ancient and heterogeneous group of double-stranded DNA viruses that preferentially infect the cutaneous and mucocutaneous epithelia of vertebrates. Persistent infection by specific oncogenic human papillomaviruses (HPVs), including HPV58, has been established as the primary cause of cervical cancer. In this work, we reveal the complex evolutionary history of HPV58 variants that explains the heterogeneity of oncogenic potential and geographic distribution. Our data suggest that HPV58 variants may have coevolved with archaic hominins and dispersed across the planet through host interbreeding and gene flow. Certain genes and codons of HPV58 variants representing higher carcinogenic potential and/or that are under positive selection may have important implications for viral host specificity, pathogenesis, and disease prevention.

Characterization of HIV diversity, phylodynamics and drug resistance in Washington, DC

BACKGROUND

Washington DC has a high burden of HIV with a 2.0% HIV prevalence. The city is a national and international hub potentially containing a broad diversity of HIV variants; yet few sequences from DC are available on GenBank to assess the evolutionary history of HIV in the US capital. Towards this general goal, here we analyze extensive sequence data and investigate HIV diversity, phylodynamics, and drug resistant mutations (DRM) in DC.

METHODS

Molecular HIV-1 sequences were collected from participants infected through 2015 as part of the DC Cohort, a longitudinal observational study of HIV+ patients receiving care at 13 DC clinics. Sequences were paired with Cohort demographic, risk, and clinical data and analyzed using maximum likelihood, Bayesian and coalescent approaches of phylogenetic, network and population genetic inference. We analyzed 601 sequences from 223 participants for int (~864 bp) and 2,810 sequences from 1,659 participants for PR/RT (~1497 bp).

RESULTS

Ninety-nine and 94% of the int and PR/RT sequences, respectively, were identified as subtype B, with 14 non-B subtypes also detected. Phylodynamic analyses of US born infected individuals showed that HIV population size varied little over time with no significant decline in diversity. Phylogenetic analyses grouped 13.5% of the int sequences into 14 clusters of 2 or 3 sequences, and 39.0% of the PR/RT sequences into 203 clusters of 2-32 sequences. Network analyses grouped 3.6% of the int sequences into 4 clusters of 2 sequences, and 10.6% of the PR/RT sequences into 76 clusters of 2-7 sequences. All network clusters were detected in our phylogenetic analyses. Higher proportions of clustered sequences were found in zip codes where HIV prevalence is highest (r = 0.607; P<0.00001). We detected a high prevalence of DRM for both int (17.1%) and PR/RT (39.1%), but only 8 int and 12 PR/RT amino acids were identified as under adaptive selection. We observed a significant (P<0.0001) association between main risk factors (men who have sex with men and heterosexuals) and genotypes in the five well-supported clusters with sufficient sample size for testing.

DISCUSSION

Pairing molecular data with clinical and demographic data provided novel insights into HIV population dynamics in Washington, DC. Identification of populations and geographic locations where clustering occurs can inform and complement active surveillance efforts to interrupt HIV transmission.

Mechanistic insights into molecular evolution of species-specific differential glycosaminoglycan binding surfaces in growth-related oncogene chemokines

Chemokines are chemotactic cytokines involved in leucocyte trafficking to infected tissue. Growth-related oncogene (GRO) chemokines namely CXCL1, CXCL2 and CXCL3 are neutrophil activating chemokines sharing a conserved three-dimensional structure, but encompassing functional diversity due to gene duplication and evolutionary events. However, the evolutionary mechanisms including selection pressures involved in diversification of GRO genes have not yet been characterized. Here, we performed comprehensive evolutionary analysis of GRO genes among different mammalian species. Phylogenetic analysis illustrated a species-specific evolution pattern. Selection analysis evidenced that these genes have undergone concerted evolution. Seventeen positively selected sites were obtained, although the majority of the protein is under purifying selection. Interestingly, these positively selected sites are more concentrated on the C-terminal/glycosaminoglycan (GAG) binding and dimerization segment compared to receptor binding domain. Substitution rate analysis confirmed the C-terminal domain of GRO genes as the highest substituted segment. Further, structural analysis established that the nucleotide alterations in the GAG binding domain are the source of surface charge modulation, thus generating the differential GAG binding surfaces and multiple binding sites as per evolutionary pressure, although the helical surface is primordial for GAG binding. Indeed, such variable electrostatic surfaces are crucial to regulate chemokine gradient formation during a host's defence against pathogens and also explain the significance of chemokine promiscuity.

Balancing selection maintains diversity in a cold tolerance gene in broadly distributed live oaks

Cold poses major physiological challenges to plants, especially long-lived trees. In trees occurring along variable temperature clines, the expected direction and consequences of selection on cold acclimation ability and freezing tolerance are not straightforward. Here we estimated selection in cold acclimation genes at two evolutionary timescales in all seven species of the American live oaks (Quercus subsection Virentes). Two cold response candidate genes were chosen: ICE1, a key gene in the cold acclimation pathway, and HOS1, which modulates cold response by negatively regulating ICE1. Two housekeeping genes, GAPDB and CHR11, were also analyzed. At the shallow evolutionary timescale, we demonstrate that HOS1 experienced recent balancing selection in the two most broadly distributed species, Q. virginiana and Q. oleoides. At a deeper evolutionary scale, a codon-based model of evolution revealed the signature of negative selection in ICE1. In contrast, three positively selected codons have been identified in HOS1, possibly a signature of the diversification of Virentes into warmer climates from a freezing adapted lineage of oaks. Our findings indicate that evolution has favored diversity in cold tolerance modulation through balancing selection in HOS1 while maintaining core cold acclimation ability, as evidenced by purifying selection in ICE1.

REST, a master regulator of neurogenesis, evolved under strong positive selection in humans and in non human primates

The transcriptional repressor REST regulates many neuronal genes by binding RE1 motifs. About one third of human RE1s are recently evolved and specific to primates. As changes in the activity of a transcription factor reverberate on its downstream targets, we assessed whether REST displays fast evolutionary rates in primates. We show that REST was targeted by very strong positive selection during primate evolution. Positive selection was also evident in the human lineage, with six selected sites located in a region that surrounds a VNTR in exon 4. Analysis of expression data indicated that REST brain expression peaks during aging in humans but not in other primates. Because a REST coding variant (rs3796529) was previously associated with protection from hippocampal atrophy in elderly subjects with mild cognitive impairment (MCI), we analyzed a cohort of Alzheimer disease (AD) continuum patients. Genotyping of two coding variants (rs3796529 and rs2227902) located in the region surrounding the VNTR indicated a role for rs2227902 in modulation of hippocampal volume loss, indirectly confirming a role for REST in neuroprotection. Experimental studies will be instrumental to determine the functional effect of positively selected sites in REST and the role of REST variants in neuropreservation/neurodegeneration.

Genetic and codon usage bias analyses of polymerase genes of equine influenza virus and its relation to evolution

BACKGROUND

Equine influenza is a major health problem of equines worldwide. The polymerase genes of influenza virus have key roles in virus replication, transcription, transmission between hosts and pathogenesis. Hence, the comprehensive genetic and codon usage bias of polymerase genes of equine influenza virus (EIV) were analyzed to elucidate the genetic and evolutionary relationships in a novel perspective.

RESULTS

The group - specific consensus amino acid substitutions were identified in all polymerase genes of EIVs that led to divergence of EIVs into various clades. The consistent amino acid changes were also detected in the Florida clade 2 EIVs circulating in Europe and Asia since 2007. To study the codon usage patterns, a total of 281,324 codons of polymerase genes of EIV H3N8 isolates from 1963 to 2015 were systemically analyzed. The polymerase genes of EIVs exhibit a weak codon usage bias. The ENc-GC3s and Neutrality plots indicated that natural selection is the major influencing factor of codon usage bias, and that the impact of mutation pressure is comparatively minor. The methods for estimating host imposed translation pressure suggested that the polymerase acidic (PA) gene seems to be under less translational pressure compared to polymerase basic 1 (PB1) and polymerase basic 2 (PB2) genes. The multivariate statistical analysis of polymerase genes divided EIVs into four evolutionary diverged clusters - Pre-divergent, Eurasian, Florida sub-lineage 1 and 2.

CONCLUSIONS

Various lineage specific amino acid substitutions observed in all polymerase genes of EIVs and especially, clade 2 EIVs underwent major variations which led to the emergence of a phylogenetically distinct group of EIVs originating from Richmond/1/07. The codon usage bias was low in all the polymerase genes of EIVs that was influenced by the multiple factors such as the nucleotide compositions, mutation pressure, aromaticity and hydropathicity. However, natural selection was the major influencing factor in defining the codon usage patterns and evolution of polymerase genes of EIVs.

Enter the Dragon: The Dynamic and Multifunctional Evolution of Anguimorpha Lizard Venoms

While snake venoms have been the subject of intense study, comparatively little work has been done on lizard venoms. In this study, we have examined the structural and functional diversification of anguimorph lizard venoms and associated toxins, and related these results to dentition and predatory ecology. Venom composition was shown to be highly variable across the 20 species of Heloderma, Lanthanotus, and Varanus included in our study. While kallikrein enzymes were ubiquitous, they were also a particularly multifunctional toxin type, with differential activities on enzyme substrates and also ability to degrade alpha or beta chains of fibrinogen that reflects structural variability. Examination of other toxin types also revealed similar variability in their presence and activity levels. The high level of venom chemistry variation in varanid lizards compared to that of helodermatid lizards suggests that venom may be subject to different selection pressures in these two families. These results not only contribute to our understanding of venom evolution but also reveal anguimorph lizard venoms to be rich sources of novel bioactive molecules with potential as drug design and development lead compounds.

Superoxide dismutase 1 is positively selected to minimize protein aggregation in great apes

Positive (adaptive) selection has recently been implied in human superoxide dismutase 1 (SOD1), a highly abundant antioxidant protein with energy signaling and antiaging functions, one of very few examples of direct selection on a human protein product (exon); the molecular drivers of this selection are unknown. We mapped 30 extant SOD1 sequences to the recently established mammalian species tree and inferred ancestors, key substitutions, and signatures of selection during the protein's evolution. We detected elevated substitution rates leading to great apes (Hominidae) at ~1 per 2 million years, significantly higher than in other primates and rodents, although these paradoxically generally evolve much faster. The high evolutionary rate was partly due to relaxation of some selection pressures and partly to distinct positive selection of SOD1 in great apes. We then show that higher stability and net charge and changes at the dimer interface were selectively introduced upon separation from old world monkeys and lesser apes (gibbons). Consequently, human, chimpanzee and gorilla SOD1s have a net charge of -6 at physiological pH, whereas the closely related gibbons and macaques have -3. These features consistently point towards selection against the malicious aggregation effects of elevated SOD1 levels in long-living great apes. The findings mirror the impact of human SOD1 mutations that reduce net charge and/or stability and cause ALS, a motor neuron disease characterized by oxidative stress and SOD1 aggregates and triggered by aging. Our study thus marks an example of direct selection for a particular chemical phenotype (high net charge and stability) in a single human protein with possible implications for the evolution of aging.

Distinct selective forces and Neanderthal introgression shaped genetic diversity at genes involved in neurodevelopmental disorders

In addition to high intelligence, humans evolved specialized social-cognitive skills, which are specifically affected in children with autism spectrum disorder (ASD). Genes affected in ASD represent suitable candidates to study the evolution of human social cognition. We performed an evolutionary analysis on 68 genes associated to neurodevelopmental disorders; our data indicate that genetic diversity was shaped by distinct selective forces, including natural selection and introgression from archaic hominins. We discuss the possibility that segregation distortion during spermatogenesis accounts for a subset of ASD mutations. Finally, we detected modern-human-specific alleles in DYRK1A and TCF4. These variants are located within regions that display chromatin features typical of transcriptional enhancers in several brain areas, strongly suggesting a regulatory role. These SNPs thus represent candidates for association with neurodevelopmental disorders, and await experimental validation in future studies.

Retention of functional variation despite extreme genomic erosion: MHC allelic repertoires in the Lynx genus

BACKGROUND

Demographic bottlenecks erode genetic diversity and may increase endangered species' extinction risk via decreased fitness and adaptive potential. The genetic status of species is generally assessed using neutral markers, whose dynamic can differ from that of functional variation due to selection. The MHC is a multigene family described as the most important genetic component of the mammalian immune system, with broad implications in ecology and evolution. The genus Lynx includes four species differing immensely in demographic history and population size, which provides a suitable model to study the genetic consequences of demographic declines: the Iberian lynx being an extremely bottlenecked species and the three remaining ones representing common and widely distributed species. We compared variation in the most variable exon of the MHCI and MHCII-DRB loci among the four species of the Lynx genus.

RESULTS

The Iberian lynx was characterised by lower number of MHC alleles than its sister species (the Eurasian lynx). However, it maintained most of the functional genetic variation at MHC loci present in the remaining and genetically healthier lynx species at all nucleotide, amino acid, and supertype levels.

CONCLUSIONS

Species-wide functional genetic diversity can be maintained even in the face of severe population bottlenecks, which caused devastating whole genome genetic erosion. This could be the consequence of divergent alleles being retained across paralogous loci, an outcome that, in the face of frequent gene conversion, may have been favoured by balancing selection.

Evolution of nonspectral rhodopsin function at high altitudes

High-altitude environments present a range of biochemical and physiological challenges for organisms through decreases in oxygen, pressure, and temperature relative to lowland habitats. Protein-level adaptations to hypoxic high-altitude conditions have been identified in multiple terrestrial endotherms; however, comparable adaptations in aquatic ectotherms, such as fishes, have not been as extensively characterized. In enzyme proteins, cold adaptation is attained through functional trade-offs between stability and activity, often mediated by substitutions outside the active site. Little is known whether signaling proteins [e.g., G protein-coupled receptors (GPCRs)] exhibit natural variation in response to cold temperatures. Rhodopsin (RH1), the temperature-sensitive visual pigment mediating dim-light vision, offers an opportunity to enhance our understanding of thermal adaptation in a model GPCR. Here, we investigate the evolution of rhodopsin function in an Andean mountain catfish system spanning a range of elevations. Using molecular evolutionary analyses and site-directed mutagenesis experiments, we provide evidence for cold adaptation in RH1. We find that unique amino acid substitutions occur at sites under positive selection in high-altitude catfishes, located at opposite ends of the RH1 intramolecular hydrogen-bonding network. Natural high-altitude variants introduced into these sites via mutagenesis have limited effects on spectral tuning, yet decrease the stability of dark-state and light-activated rhodopsin, accelerating the decay of ligand-bound forms. As found in cold-adapted enzymes, this phenotype likely compensates for a cold-induced decrease in kinetic rates-properties of rhodopsin that mediate rod sensitivity and visual performance. Our results support a role for natural variation in enhancing the performance of GPCRs in response to cold temperatures.

Hypoxia Inducible Factor (HIF) transcription factor family expansion, diversification, divergence and selection in eukaryotes

Hypoxia inducible factor (HIF) transcription factors are crucial for regulating a variety of cellular activities in response to oxygen stress (hypoxia). In this study, we determine the evolutionary history of HIF genes and their associated transactivation domains, as well as perform selection and functional divergence analyses across their four characteristic domains. Here we show that the HIF genes are restricted to metazoans: At least one HIF-α homolog is found within the genomes of non-bilaterians and bilaterian invertebrates, while most vertebrate genomes contain between two and six HIF-α genes. We also find widespread purifying selection across all four characteristic domain types, bHLH, PAS, NTAD, CTAD, in HIF-α genes, and evidence for Type I functional divergence between HIF-1α, HIF-2α /EPAS, and invertebrate HIF genes. Overall, we describe the evolutionary histories of the HIF transcription factor gene family and its associated transactivation domains in eukaryotes. We show that the NTAD and CTAD domains appear de novo, without any appearance outside of the HIF-α subunits. Although they both appear in invertebrates as well as vertebrate HIF- α sequences, there seems to have been a substantial loss across invertebrates or were convergently acquired in these few lineages. We reaffirm that HIF-1α is phylogenetically conserved among most metazoans, whereas HIF-2α appeared later. Overall, our findings can be attributed to the substantial integration of this transcription factor family into the critical tasks associated with maintenance of oxygen homeostasis and vascularization, particularly in the vertebrate lineage.

HIV Trafficking Between Blood and Semen During Early Untreated HIV Infection

BACKGROUND

Understanding the dynamics of HIV across anatomic compartments is important to design effective eradication strategies. In this study, we evaluated viral trafficking between blood and semen during primary HIV infection in 6 antiretroviral-naive men who have sex with men.

METHODS

Deep sequencing data of HIV env were generated from longitudinal blood plasma, peripheral blood mononuclear cells, and seminal plasma samples. The presence or absence of viral compartmentalization was assessed using tree-based Slatkin-Maddison and distance-based Fst methods. Phylogeographic analyses were performed using a discrete Bayesian asymmetric approach of diffusion with Markov jump count estimation to evaluate the gene flow between blood and semen during primary HIV infection. Levels of DNA from human herpesviruses and selected inflammatory cytokines were also measured on genital secretions collected at baseline to evaluate potential correlates of increased viral migration between anatomic compartments.

RESULTS

We detected varying degrees of compartmentalization in all 6 individuals evaluated. None of them maintained viral compartmentalization between blood and seminal plasma throughout the analyzed time points. Phylogeographic analyses revealed that the HIV population circulating in blood plasma populated the seminal compartment during the earliest stages of infection. In our limited data set, we found no association between local inflammation or herpesvirus shedding at baseline and viral trafficking between semen and blood.

CONCLUSIONS

The early spread of virus from blood plasma to genital tract and the complex viral interplay between these compartments suggest that viral eradication efforts will require monitoring viral subpopulations in anatomic sites and viral trafficking during the course of infection.

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.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

PvGAMA reticulocyte binding activity: predicting conserved functional regions by natural selection analysis

Background

Adhesin proteins are used by Plasmodium parasites to bind and invade target cells. Hence, characterising molecules that participate in reticulocyte interaction is key to understanding the molecular basis of Plasmodium vivax invasion. This study focused on predicting functionally restricted regions of the P. vivax GPI-anchored micronemal antigen (PvGAMA) and characterising their reticulocyte binding activity.

Results

The pvgama gene was initially found in P. vivax VCG-I strain schizonts. According to the genetic diversity analysis, PvGAMA displayed a size polymorphism very common for antigenic P. vivax proteins. Two regions along the antigen sequence were highly conserved among species, having a negative natural selection signal. Interestingly, these regions revealed a functional role regarding preferential target cell adhesion.

Conclusions

To our knowledge, this study describes PvGAMA reticulocyte binding properties for the first time. Conserved functional regions were predicted according to natural selection analysis and their binding ability was confirmed. These findings support the notion that PvGAMA may have an important role in P. vivax merozoite adhesion to its target cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2183-8) contains supplementary material, which is available to authorized users.

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.

Purifying selection and genetic drift shaped Pleistocene evolution of the mitochondrial genome in an endangered Australian freshwater fish

Genetic variation in mitochondrial genes could underlie metabolic adaptations because mitochondrially encoded proteins are directly involved in a pathway supplying energy to metabolism. Macquarie perch from river basins exposed to different climates differ in size and growth rate, suggesting potential presence of adaptive metabolic differences. We used complete mitochondrial genome sequences to build a phylogeny, estimate lineage divergence times and identify signatures of purifying and positive selection acting on mitochondrial genes for 25 Macquarie perch from three basins: Murray-Darling Basin (MDB), Hawkesbury-Nepean Basin (HNB) and Shoalhaven Basin (SB). Phylogenetic analysis resolved basin-level clades, supporting incipient speciation previously inferred from differentiation in allozymes, microsatellites and mitochondrial control region. The estimated time of lineage divergence suggested an early- to mid-Pleistocene split between SB and the common ancestor of HNB+MDB, followed by mid-to-late Pleistocene splitting between HNB and MDB. These divergence estimates are more recent than previous ones. Our analyses suggested that evolutionary drivers differed between inland MDB and coastal HNB. In the cooler and more climatically variable MDB, mitogenomes evolved under strong purifying selection, whereas in the warmer and more climatically stable HNB, purifying selection was relaxed. Evidence for relaxed selection in the HNB includes elevated transfer RNA and 16S ribosomal RNA polymorphism, presence of potentially mildly deleterious mutations and a codon (ATP6113) displaying signatures of positive selection (ratio of nonsynonymous to synonymous substitution rates (dN/dS) >1, radical change of an amino-acid property and phylogenetic conservation across the Percichthyidae). In addition, the difference could be because of stronger genetic drift in the smaller and historically more subdivided HNB with low per-population effective population sizes.

The herpevac trial for women: Sequence analysis of glycoproteins from viruses obtained from infected subjects

The Herpevac Trial for Women revealed that three dose HSV-2 gD vaccine was 58% protective against culture-positive HSV-1 genital disease, but it was not protective against HSV-2 infection or disease. To determine whether vaccine-induced immune responses had selected for a particular gD sequence in strains infecting vaccine recipients compared with viruses infecting control subjects, genetic sequencing studies were carried out on viruses isolated from subjects infected with HSV-1 or HSV-2. We identified naturally occurring variants among the gD sequences obtained from 83 infected subjects. Unique or low frequency amino acid substitutions in the ectodomain of gD were found in 6 of 39 HSV-1-infected subjects and in 7 of 44 HSV-2-infected subjects. However, no consistent amino acid change was identified in isolates from gD-2 vaccine recipients compared with infected placebo recipients. gC and gE surround and partially shield gD from neutralizing antibody, and gB also participates closely in the viral entry process. Therefore, these genes were sequenced from a number of isolates to assess whether sequence variation may alter protein conformation and influence the virus strain’s capacity to be neutralized by vaccine-induced antibody. gC and gE genes sequenced from HSV-1-infected subjects showed more variability than their HSV-2 counterparts. The gB sequences of HSV-1 oral isolates resembled each other more than they did gB sequences rom genital isolates. Overall, however, comparison of glycoprotein sequences of viral isolates obtained from infected subjects did not reveal any singular selective pressure on the viral cell attachment protein or surrounding glycoproteins due to administration of gD-2 vaccine.

Insights into intercontinental spread of Zika virus

The epidemic of Zika virus (ZIKV) infection in South America has led to World Health Organization's declaration of a Public Health Emergency of International Concern. To further inform effective public health policy, an understanding of ZIKV's transmission mechanisms is crucial. To characterize the intercontinental transmission of ZIKV, we compiled and analyzed more than 250 gene sequences together with their sequence-related geographic and temporal information, sampled across 27 countries spanning from 1947 to 2016. After filtering and selecting appropriate sequences, extensive phylogenetic analyses were performed. Although phylogeographic reconstruction supported the transmission route of the virus in Africa, South-eastern Asia, Oceania and Latin America, we discovered that the Eastern Africa origin of ZIKV was disputable. On a molecular level, purifying selection was found to be largely responsible for the evolution of non-structural protein 5 and envelope protein E. Our dataset and ancestral sequences reconstruction analysis captured previously unidentified amino acid changes during evolution. Finally, based on the estimation of the time to the most recent common ancestors for the non-structural protein 5 gene, we hypothesized potential specific historic events that occurred in the 1940s and might have facilitated the spread of Zika virus from Africa to South-eastern Asia. Our findings provide new insights into the transmission characteristics of ZIKV, while further genetic and serologic studies are warranted to support the design of tailored prevention strategies.

Long-wavelength sensitive opsin (LWS) gene variability in Neotropical cichlids (Teleostei: Cichlidae)

Cichlid fishes are an important group in evolutionary biology due to their fast speciation. This group depends widely of vision for feeding and reproduction. During the evolutionary process it plays a significant role in interspecific and intraspecific recognition and in its ecology. The molecular basis of vision is formed by the interaction of the protein opsin and retinal chromophore. Long-wavelength sensitive opsin (LWS) gene is the most variable among the opsin genes and it has an ecological significance. Current assay identifies interspecific variation of Neotropical cichlids that would modify the spectral properties of the LWS opsin protein and codons selected. Neotropical species present more variable sites for LWS gene than those of the African lakes species. The LWS opsin gene in Crenicichla britskii has a higher amino acid similarity when compared to that in the African species, but the variable regions do not overlap. Neotropical cichlids accumulate larger amounts of variable sites for LWS opsin gene, probably because they are spread over a wider area and submitted to a wider range of selective pressures by inhabiting mainly lotic environments. Furthermore, the codons under selection are different when compared to those of the African cichlids.

Genotype-Specific Evolution of Hepatitis E Virus

Hepatitis E virus (HEV) is the most common cause of acute viral hepatitis globally. HEV comprises four genotypes with different geographic distributions and host ranges. We utilize this natural case-control study for investigating the evolution of zoonotic viruses compared to single-host viruses, using 244 near-full-length HEV genomes. Genome-wide estimates of the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS ratio) located a region of overlapping reading frames, which is subject to positive selection in genotypes 3 and 4. The open reading frames (ORFs) involved have functions related to host-pathogen interaction, so genotype-specific evolution of these regions may reflect their fitness. Bayesian inference of evolutionary rates shows that genotypes 3 and 4 have significantly higher rates than genotype 1 across all ORFs. Reconstruction of the phylogenies of zoonotic genotypes demonstrates significant intermingling of isolates between hosts. We speculate that the genotype-specific differences may result from cyclical adaptation to different hosts in genotypes 3 and 4.IMPORTANCE Hepatitis E virus (HEV) is increasingly recognized as a pathogen that affects both the developing and the developed world. While most often clinically mild, HEV can be severe or fatal in certain demographics, such as expectant mothers. Like many other viral pathogens, HEV has been classified into several distinct genotypes. We show that most of the HEV genome is evolutionarily constrained. One locus of positive selection is unusual in that it encodes two distinct protein products. We are the first to detect positive selection in this overlap region. Genotype 1, which infects humans only, appears to be evolving differently from genotypes 3 and 4, which infect multiple species, possibly because genotypes 3 and 4 are unable to achieve the same fitness due to repeated host jumps.

Modular Organization of Residue-Level Contacts Shapes the Selection Pressure on Individual Amino Acid Sites of Ribosomal Proteins

Understanding the molecular evolution of macromolecular complexes in the light of their structure, assembly, and stability is of central importance. Here, we address how the modular organization of native molecular contacts shapes the selection pressure on individual residue sites of ribosomal complexes. The bacterial ribosomal complex is represented as a residue contact network where nodes represent amino acid/nucleotide residues and edges represent their van der Waals interactions. We find statistically overrepresented native amino acid-nucleotide contacts (OaantC, one amino acid contacts one or multiple nucleotides, internucleotide contacts are disregarded). Contact number is defined as the number of nucleotides contacted. Involvement of individual amino acids in OaantCs with smaller contact numbers is more random, whereas only a few amino acids significantly contribute to OaantCs with higher contact numbers. An investigation of structure, stability, and assembly of bacterial ribosome depicts the involvement of these OaantCs in diverse biophysical interactions stabilizing the complex, including high-affinity protein-RNA contacts, interprotein cooperativity, intersubunit bridge, packing of multiple ribosomal RNA domains, etc. Amino acid-nucleotide constituents of OaantCs with higher contact numbers are generally associated with significantly slower substitution rates compared with that of OaantCs with smaller contact numbers. This evolutionary rate heterogeneity emerges from the strong purifying selection pressure that conserves the respective amino acid physicochemical properties relevant to the stabilizing interaction with OaantC nucleotides. An analysis of relative molecular orientations of OaantC residues and their interaction energetics provides the biophysical ground of purifying selection conserving OaantC amino acid physicochemical properties.

Comparative analysis of protein evolution in the genome of pre-epidemic and epidemic Zika virus

Zika virus (ZIKV) causes microcephaly in congenital infection, neurological disorders, and poor pregnancy outcome and no vaccine is available for use in humans or approved. Although ZIKV was first discovered in 1947, the exact mechanism of virus replication and pathogenesis remains unknown. Recent outbreaks of Zika virus in the Americas clearly suggest a human-mosquito cycle or urban cycle of transmission. Understanding the conserved and adaptive features in the evolution of ZIKV genome will provide a hint on the mechanism of ZIKV adaptation to a new cycle of transmission. Here, we show comprehensive analysis of protein evolution of ZIKV strains including the current 2015-16 outbreak. To identify the constraints on ZIKV evolution, selection pressure at individual codons, immune epitopes and co-evolving sites were analyzed. Phylogenetic trees show that the ZIKV strains of the Asian genotype form distinct cluster and share a common ancestor with African genotype. The TMRCA (Time to the Most Recent Common Ancestor) for the Asian lineage and the subsequently evolved Asian human strains was calculated at 88 and 34years ago, respectively. The proteome of current 2015/16 epidemic ZIKV strains of Asian genotype was found to be genetically conserved due to genome-wide negative selection, with limited positive selection. We identified a total of 16 amino acid substitutions in the epidemic and pre-epidemic strains from human, mosquito, and monkey hosts. Negatively selected amino acid sites of Envelope protein (E-protein) (positions 69, 166, and 174) and NS5 (292, 345, and 587) were located in central dimerization domains and C-terminal RNA-directed RNA polymerase regions, respectively. The predicted 137 (92 CD4 TCEs; 45 CD8 TCEs) immunogenic peptide chains comprising negatively selected amino acid sites can be considered as suitable target for sub-unit vaccine development, as these sites are less likely to generate immune-escape variants due to strong functional constrains operating on them. The targeted changes at the amino acid level may contribute to better adaptation of ZIKV strains to human-mosquito cycle or urban cycle of transmission.

The diversification of begomovirus populations is predominantly driven by mutational dynamics

Begomoviruses (single-stranded DNA plant viruses) are responsible for serious agricultural threats. Begomovirus populations exhibit a high degree of within-host genetic variation and evolve as quickly as RNA viruses. Although the recombination-prone nature of begomoviruses has been extensively demonstrated, the relative contribution of recombination and mutation to the genetic variation of begomovirus populations has not been assessed. We estimated the genetic variability of begomovirus datasets from around the world. An uneven distribution of genetic variation across the length of the cp and rep genes due to recombination was evident from our analyses. To estimate the relative contributions of recombination and mutation to the genetic variability of begomoviruses, we mapped all substitutions over maximum likelihood trees and counted the number of substitutions on branches which were associated with recombination (η_r) and mutation (η_μ). In addition, we also estimated the per generation relative rates of both evolutionary mechanisms (r/μ) to express how frequently begomovirus genomes are affected by recombination relative to mutation. We observed that the composition of genetic variation in all begomovirus datasets was dominated by mutation. Additionally, the low correlation between the estimates indicated that the relative contributions of recombination and mutation are not necessarily a function of their relative rates. Our results show that, although a considerable fraction of the genetic variation levels could be assigned to recombination, it was always lower than that due to mutation, indicating that the diversification of begomovirus populations is predominantly driven by mutational dynamics.

Porphyromonas gingivalis Uses Specific Domain Rearrangements and Allelic Exchange to Generate Diversity in Surface Virulence Factors

Porphyromonas gingivalis is a keystone pathogen of chronic periodontitis. The virulence of P. gingivalis is reported to be strain related and there are currently a number of strain typing schemes based on variation in capsular polysaccharide, the major and minor fimbriae and adhesin domains of Lys-gingipain (Kgp), amongst other surface proteins. P. gingivalis can exchange chromosomal DNA between strains by natural competence and conjugation. The aim of this study was to determine the genetic variability of P. gingivalis strains sourced from international locations over a 25-year period and to determine if variability in surface virulence factors has a phylogenetic basis. Whole genome sequencing was performed on 13 strains and comparison made to 10 previously sequenced strains. A single nucleotide polymorphism-based phylogenetic analysis demonstrated a shallow tri-lobed phylogeny. There was a high level of reticulation in the phylogenetic network, demonstrating extensive horizontal gene transfer between the strains. Two highly conserved variants of the catalytic domain of the major virulence factor the Kgp proteinase (Kgp_catI and Kgp_catII) were found. There were three variants of the fourth Kgp C-terminal cleaved adhesin domain. Specific variants of the cell surface proteins FimA, FimCDE, MfaI, RagAB, Tpr, and PrtT were also identified. The occurrence of all these variants in the P. gingivalis strains formed a mosaic that was not related to the SNP-based phylogeny. In conclusion P. gingivalis uses domain rearrangements and genetic exchange to generate diversity in specific surface virulence factors.

Identification of positive selection in genes is greatly improved by using experimentally informed site-specific models

Background

Sites of positive selection are identified by comparing observed evolutionary patterns to those expected under a null model for evolution in the absence of such selection. For protein-coding genes, the most common null model is that nonsynonymous and synonymous mutations fix at equal rates; this unrealistic model has limited power to detect many interesting forms of selection.

Results

I describe a new approach that uses a null model based on experimental measurements of a gene’s site-specific amino-acid preferences generated by deep mutational scanning in the lab. This null model makes it possible to identify both diversifying selection for repeated amino-acid change and differential selection for mutations to amino acids that are unexpected given the measurements made in the lab. I show that this approach identifies sites of adaptive substitutions in four genes (lactamase, Gal4, influenza nucleoprotein, and influenza hemagglutinin) far better than a comparable method that simply compares the rates of nonsynonymous and synonymous substitutions.

Conclusions

As rapid increases in biological data enable increasingly nuanced descriptions of the constraints on individual protein sites, approaches like the one here can improve our ability to identify many interesting forms of selection in natural sequences.

Reviewers

This article was reviewed by Sebastian Maurer-Stroh, Olivier Tenaillon, and Tal Pupko. All three reviewers are members of the Biology Direct editorial board.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0172-z) contains supplementary material, which is available to authorized users.

Unexpected differences in the population genetics of phasmavirids (Bunyavirales) from subarctic ponds

Little is known of the evolution of RNA viruses in aquatic systems. Here, we assess the genetic connectivity of two bunyaviruses (Kigluaik phantom orthophasmavirus or KIGV and Nome phantom orthophasmavirus or NOMV) with zooplanktonic hosts from subarctic ponds. We expected weak genetic structure among populations as the hosts (phantom midges) have a terrestrial winged dispersal stage. To test whether their respective viruses mirror this structure, we collected and analyzed population datasets from 21 subarctic freshwater ponds and obtained sequences from all four genes in the viral genomes. Prevalence averaged 66 per cent for 514 host specimens and was not significantly different between recently formed thaw ponds and glacial ponds. Unexpectedly, KIGV from older ponds showed pronounced haplotype divergence with little evidence of genetic connectivity. However, KIGV populations from recent thaw ponds appeared to be represented by a closely related haplotype group, perhaps indicating a genotypic dispersal bias. Unlike KIGV, NOMV had modest structure and diversity in recently formed thaw ponds. For each virus, we found elevated genetic diversity relative to the host, but similar population structures to the host. Our results suggest that non-random processes such as virus-host interactions, genotypic bias, and habitat effects differ among polar aquatic RNA viruses.

Porcine circovirus type 2 (PCV2) evolution before and after the vaccination introduction: A large scale epidemiological study

Since their commercialization, vaccines against Porcine circovirus type 2 (PCV2) have been the cornerstone control strategy. Nevertheless, the periodic emergence of new genotype waves and the recent reports of vaccine failure outbreaks have raised the question if widespread vaccination strategies could have driven viral evolution and affected different genotype fitness. To investigate this issue an in-deep analysis, based on a bioinformatics and biostatistics approach, has been implemented. ORF2 sequences from vaccinated and non-vaccinated populations (i.e. domestic pigs before and after vaccine introduction and wild boars) were considered. The action of selective forces on PCV2 strains has been analyzed and compared among groups. Remarkable differences were found in the selective forces acting on viral populations circulating in different “immune environments”. Particularly for PCV2a, a directional selection promoting a change in the viral capsid away from the vaccine specific antigenic determinants has been detected after vaccine introduction. Involved amino acids were previously reported to be part of viral epitopes whose variability is responsible of immune escape. Our findings support a change in PCV2 evolutionary pattern after widespread vaccination introduction and stress once more the compulsoriness of a continuous monitoring of PCV2 epidemiology to promptly act in response to the emergence of possible vaccine-escaping mutants.

Differential Infection Patterns and Recent Evolutionary Origins of Equine Hepaciviruses in Donkeys

The hepatitis C virus (HCV) is a major human pathogen. Genetically related viruses in animals suggest a zoonotic origin of HCV. The closest relative of HCV is found in horses (termed equine hepacivirus [EqHV]). However, low EqHV genetic diversity implies relatively recent acquisition of EqHV by horses, making a derivation of HCV from EqHV unlikely. To unravel the EqHV evolutionary history within equid sister species, we analyzed 829 donkeys and 53 mules sampled in nine European, Asian, African, and American countries by molecular and serologic tools for EqHV infection. Antibodies were found in 278 animals (31.5%), and viral RNA was found in 3 animals (0.3%), all of which were simultaneously seropositive. A low RNA prevalence in spite of high seroprevalence suggests a predominance of acute infection, a possible difference from the mostly chronic hepacivirus infection pattern seen in horses and humans. Limitation of transmission due to short courses of infection may explain the existence of entirely seronegative groups of animals. Donkey and horse EqHV strains were paraphyletic and 97.5 to 98.2% identical in their translated polyprotein sequences, making virus/host cospeciation unlikely. Evolutionary reconstructions supported host switches of EqHV between horses and donkeys without the involvement of adaptive evolution. Global admixture of donkey and horse hepaciviruses was compatible with anthropogenic alterations of EqHV ecology. In summary, our findings do not support EqHV as the origin of the significantly more diversified HCV. Identification of a host system with predominantly acute hepacivirus infection may enable new insights into the chronic infection pattern associated with HCV.

IMPORTANCE

The evolutionary origins of the human hepatitis C virus (HCV) are unclear. The closest animal-associated relative of HCV occurs in horses (equine hepacivirus [EqHV]). The low EqHV genetic diversity implies a relatively recent acquisition of EqHV by horses, limiting the time span for potential horse-to-human infections in the past. Horses are genetically related to donkeys, and EqHV may have cospeciated with these host species. Here, we investigated a large panel of donkeys from various countries using serologic and molecular tools. We found EqHV to be globally widespread in donkeys and identify potential differences in EqHV infection patterns, with donkeys potentially showing enhanced EqHV clearance compared to horses. We provide strong evidence against EqHV cospeciation and for its capability to switch hosts among equines. Differential hepacivirus infection patterns in horses and donkeys may enable new insights into the chronic infection pattern associated with HCV.

Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15^th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

Zoonoses account for most recently emerged infectious diseases of humans, although little is known about the evolutionary mechanisms involved in cross-species virus transmission. Understanding the evolutionary patterns and processes that underpin such cross-species transmission is of importance for predicting the spread of zoonotic infections, and hence to their ultimate control. We present a large-scale and detailed reconstruction of the evolutionary history of rabies virus (RABV) in domestic and wildlife animal species. RABV is of particular interest as it is capable of infecting many mammals but, paradoxically, is only maintained in distinct epidemiological cycles associated with animal species from the orders Carnivora and Chiroptera. We show that bat-related RABV and dog-related RABV have experienced very different evolutionary dynamics, and that host jumps are sometimes characterized by significant increases in evolutionary rate. Among Carnivora, the association between RABV and particular host species most likely arose from a combination of the historical human-mediated spread of the virus and jumps into new primary host species. In addition, we show that changes in host species are associated with multiple evolutionary pathways including the occurrence of host-specific parallel evolution. Overall, our data indicate that the establishment of dog-related RABV in new carnivore hosts may only require subtle adaptive evolution.

Functional and Evolutionary Analyses Identify Proteolysis as a General Mechanism for NLRP1 Inflammasome Activation

Inflammasomes are cytosolic multi-protein complexes that initiate immune responses to infection by recruiting and activating the Caspase-1 protease. Human NLRP1 was the first protein shown to form an inflammasome, but its physiological mechanism of activation remains unknown. Recently, specific variants of mouse and rat NLRP1 were found to be activated upon N-terminal cleavage by the anthrax lethal factor protease. However, agonists for other NLRP1 variants, including human NLRP1, are not known, and it remains unclear if they are also activated by proteolysis. Here we demonstrate that two mouse NLRP1 paralogs (NLRP1AB6 and NLRP1BB6) are also activated by N-terminal proteolytic cleavage. We also demonstrate that proteolysis within a specific N-terminal linker region is sufficient to activate human NLRP1. Evolutionary analysis of primate NLRP1 shows the linker/cleavage region has evolved under positive selection, indicative of pathogen-induced selective pressure. Collectively, these results identify proteolysis as a general mechanism of NLRP1 inflammasome activation that appears to be contributing to the rapid evolution of NLRP1 in rodents and primates.

Adaptive molecular evolution of the two-pore channel 1 gene TPC1 in the karst-adapted genus Primulina (Gesneriaceae)

BACKGROUND AND AIMS

Limestone karst areas possess high floral diversity and endemism. The genus Primulina, which contributes to the unique calcicole flora, has high species richness and exhibit specific soil-based habitat associations that are mainly distributed on calcareous karst soils. The adaptive molecular evolutionary mechanism of the genus to karst calcium-rich environments is still not well understood. The Ca²⁺-permeable channel TPC1 was used in this study to test whether its gene is involved in the local adaptation of Primulina to karst high-calcium soil environments.

METHODS

Specific amplification and sequencing primers were designed and used to amplify the full-length coding sequences of TPC1 from cDNA of 76 Primulina species. The sequence alignment without recombination and the corresponding reconstructed phylogeny tree were used in molecular evolutionary analyses at the nucleic acid level and amino acid level, respectively. Finally, the identified sites under positive selection were labelled on the predicted secondary structure of TPC1.

KEY RESULTS

Seventy-six full-length coding sequences of Primulina TPC1 were obtained. The length of the sequences varied between 2220 and 2286 bp and the insertion/deletion was located at the 5' end of the sequences. No signal of substitution saturation was detected in the sequences, while significant recombination breakpoints were detected. The molecular evolutionary analyses showed that TPC1 was dominated by purifying selection and the selective pressures were not significantly different among species lineages. However, significant signals of positive selection were detected at both TPC1 codon level and amino acid level, and five sites under positive selective pressure were identified by at least three different methods.

CONCLUSIONS

The Ca²⁺-permeable channel TPC1 may be involved in the local adaptation of Primulina to karst Ca²⁺-rich environments. Different species lineages suffered similar selective pressure associated with calcium in karst environments, and episodic diversifying selection at a few sites may play a major role in the molecular evolution of Primulina TPC1.

HIV-1 drug resistance mutations emerging on darunavir therapy in PI-naive and -experienced patients in the UK

BACKGROUND

Darunavir is considered to have a high genetic barrier to resistance. Most darunavir-associated drug resistance mutations (DRMs) have been identified through correlation of baseline genotype with virological response in clinical trials. However, there is little information on DRMs that are directly selected by darunavir in clinical settings.

OBJECTIVES

We examined darunavir DRMs emerging in clinical practice in the UK.

PATIENTS AND METHODS

Baseline and post-exposure protease genotypes were compared for individuals in the UK Collaborative HIV Cohort Study who had received darunavir; analyses were stratified for PI history. A selection analysis was used to compare the evolution of subtype B proteases in darunavir recipients and matched PI-naive controls.

RESULTS

Of 6918 people who had received darunavir, 386 had resistance tests pre- and post-exposure. Overall, 2.8% (11/386) of these participants developed emergent darunavir DRMs. The prevalence of baseline DRMs was 1.0% (2/198) among PI-naive participants and 13.8% (26/188) among PI-experienced participants. Emergent DRMs developed in 2.0% of the PI-naive group (4 mutations) and 3.7% of the PI-experienced group (12 mutations). Codon 77 was positively selected in the PI-naive darunavir cases, but not in the control group.

CONCLUSIONS

Our findings suggest that although emergent darunavir resistance is rare, it may be more common among PI-experienced patients than those who are PI-naive. Further investigation is required to explore whether codon 77 is a novel site involved in darunavir susceptibility.

Strong and widespread action of site-specific positive selection in the snake venom Kunitz/BPTI protein family

S1 family of serine peptidases is the largest family of peptidases. They are specifically inhibited by the Kunitz/BPTI inhibitors. Kunitz domain is characterized by the compact 3D structure with the most important inhibitory loops for the inhibition of S1 peptidases. In the present study we analysed the action of site-specific positive selection and its impact on the structurally and functionally important parts of the snake venom Kunitz/BPTI family of proteins. By using numerous models we demonstrated the presence of large numbers of site-specific positively selected sites that can reach between 30-50% of the Kunitz domain. The mapping of the positively selected sites on the 3D model of Kunitz/BPTI inhibitors has shown that these sites are located in the inhibitory loops 1 and 2, but also in the Kunitz scaffold. Amino acid replacements have been found exclusively on the surface, and the vast majority of replacements are causing the change of the charge. The consequence of these replacements is the change in the electrostatic potential on the surface of the Kunitz/BPTI proteins that may play an important role in the precise targeting of these inhibitors into the active site of S1 family of serine peptidases.

Characterization of the antimicrobial peptide family defensins in the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii)

Defensins comprise a family of cysteine-rich antimicrobial peptides with important roles in innate and adaptive immune defense in vertebrates. We characterized alpha and beta defensin genes in three Australian marsupials: the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii) and identified 48, 34, and 39 defensins, respectively. One hundred and twelve have the classical antimicrobial peptides characteristics required for pathogen membrane targeting, including cationic charge (between 1+ and 15+) and a high proportion of hydrophobic residues (>30%). Phylogenetic analysis shows that gene duplication has driven unique and species-specific expansions of devil, koala, and tammar wallaby beta defensins and devil alpha defensins. Defensin genes are arranged in three genomic clusters in marsupials, whereas further duplications and translocations have occurred in eutherians resulting in four and five gene clusters in mice and humans, respectively. Marsupial defensins are generally under purifying selection, particularly residues essential for defensin structural stability. Certain hydrophobic or positively charged sites, predominantly found in the defensin loop, are positively selected, which may have functional significance in defensin-target interaction and membrane insertion.

Spatio-temporal analysis of Nova virus, a divergent hantavirus circulating in the European mole in Belgium

Over the last decade, the recognized host range of hantaviruses has expanded considerably with the discovery of distinct hantaviruses in shrews, moles and bats. Unfortunately, in-depth studies of these viruses have been limited. Here we describe a comprehensive analysis of the spatial distribution, genetic diversity and evolution of Nova virus, a hantavirus that has the European mole as its natural host. Our analysis demonstrated that Nova virus has a high prevalence and widespread distribution in Belgium. While Nova virus displayed relatively high nucleotide diversity in Belgium, amino acid changes were limited. The nucleocapsid protein was subjected to strong purifying selection, reflecting the strict evolutionary constraints placed upon Nova virus by its host. Spatio-temporal analysis using Bayesian evolutionary inference techniques demonstrated that Nova virus had efficiently spread in the European mole population in Belgium, forming two distinct clades, representing east and west of Belgium. The influence of landscape barriers, in the form of the main waterways, on the dispersal velocity of Nova virus was assessed using an analytical framework for comparing Bayesian viral phylogenies with environmental landscape data. We demonstrated that waterways did not act as an environmental resistance factor slowing down Nova virus diffusion in the mole population. With this study, we provide information about the spatial diffusion of Nova virus and contribute sequence information that can be applied in further functional studies.

Annotating Protein Functional Residues by Coupling High-Throughput Fitness Profile and Homologous-Structure Analysis

Identification and annotation of functional residues are fundamental questions in protein sequence analysis. Sequence and structure conservation provides valuable information to tackle these questions. It is, however, limited by the incomplete sampling of sequence space in natural evolution. Moreover, proteins often have multiple functions, with overlapping sequences that present challenges to accurate annotation of the exact functions of individual residues by conservation-based methods. Using the influenza A virus PB1 protein as an example, we developed a method to systematically identify and annotate functional residues. We used saturation mutagenesis and high-throughput sequencing to measure the replication capacity of single nucleotide mutations across the entire PB1 protein. After predicting protein stability upon mutations, we identified functional PB1 residues that are essential for viral replication. To further annotate the functional residues important to the canonical or noncanonical functions of viral RNA-dependent RNA polymerase (vRdRp), we performed a homologous-structure analysis with 16 different vRdRp structures. We achieved high sensitivity in annotating the known canonical polymerase functional residues. Moreover, we identified a cluster of noncanonical functional residues located in the loop region of the PB1 β-ribbon. We further demonstrated that these residues were important for PB1 protein nuclear import through the interaction with Ran-binding protein 5. In summary, we developed a systematic and sensitive method to identify and annotate functional residues that are not restrained by sequence conservation. Importantly, this method is generally applicable to other proteins about which homologous-structure information is available.

IMPORTANCE

To fully comprehend the diverse functions of a protein, it is essential to understand the functionality of individual residues. Current methods are highly dependent on evolutionary sequence conservation, which is usually limited by sampling size. Sequence conservation-based methods are further confounded by structural constraints and multifunctionality of proteins. Here we present a method that can systematically identify and annotate functional residues of a given protein. We used a high-throughput functional profiling platform to identify essential residues. Coupling it with homologous-structure comparison, we were able to annotate multiple functions of proteins. We demonstrated the method with the PB1 protein of influenza A virus and identified novel functional residues in addition to its canonical function as an RNA-dependent RNA polymerase. Not limited to virology, this method is generally applicable to other proteins that can be functionally selected and about which homologous-structure information is available.

Size polymorphism and low sequence diversity in the locus encoding the Plasmodium vivax rhoptry neck protein 4 (PvRON4) in Colombian isolates

Background

Designing a vaccine against Plasmodium vivax has focused on selecting antigens involved in invasion mechanisms that must have domains with low polymorphism for avoiding allele-specific immune responses. The rhoptry neck protein 4 (RON4) forms part of the tight junction, which is essential in the invasion of hepatocytes and/or erythrocytes; however, little is known about this locus’ genetic diversity.

Methods

DNA sequences from 73 Colombian clinical isolates from pvron4 gene were analysed for characterizing their genetic diversity; pvron4 haplotype number and distribution, as well as the evolutionary forces determining diversity pattern, were assessed by population genetics and molecular evolutionary approaches.

Results

ron4 has low genetic diversity in P. vivax at sequence level; however, a variable amount of tandem repeats at the N-terminal region leads to extensive size polymorphism. This region seems to be exposed to the immune system. The central region has a putative esterase/lipase domain which, like the protein’s C-terminal fragment, is highly conserved at intra- and inter-species level. Both regions are under purifying selection.

Conclusions

pvron4 is the locus having the lowest genetic diversity described to date for P. vivax. The repeat regions in the N-terminal region could be associated with immune evasion mechanisms while the central region and the C-terminal region seem to be under functional or structural constraint. Bearing such results in mind, the PvRON4 central and/or C-terminal portions represent promising candidates when designing a subunit-based vaccine as they are aimed at avoiding an allele-specific immune response, which might limit vaccine efficacy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1563-4) contains supplementary material, which is available to authorized users.