FUBAR: a fast, unconstrained bayesian approximation for inferring selection

Diversifying Selection Analysis Predicts Antigenic Evolution of 2009 Pandemic H1N1 Influenza A Virus in Humans

Although a large number of immune epitopes have been identified in the influenza A virus (IAV) hemagglutinin (HA) protein using various experimental systems, it is unclear which are involved in protective immunity to natural infection in humans. We developed a data mining approach analyzing natural H1N1 human isolates to identify HA protein regions that may be targeted by the human immune system and can predict the evolution of IAV. We identified 16 amino acid sites experiencing diversifying selection during the evolution of prepandemic seasonal H1N1 strains and found that 11 sites were located in experimentally determined B-cell/antibody (Ab) epitopes, including three distinct neutralizing Caton epitopes: Sa, Sb, and Ca2 [A. J. Caton, G. G. Brownlee, J. W. Yewdell, and W. Gerhard, Cell 31:417-427, 1982, http://dx.doi.org/10.1016/0092-8674(82)90135-0]. We predicted that these diversified epitope regions would be the targets of mutation as the 2009 H1N1 pandemic (pH1N1) lineage evolves in response to the development of population-level protective immunity in humans. Using a chi-squared goodness-of-fit test, we identified 10 amino acid sites that significantly differed between the pH1N1 isolates and isolates from the recent 2012-2013 and 2013-2014 influenza seasons. Three of these sites were located in the same diversified B-cell/Ab epitope regions as identified in the analysis of prepandemic sequences, including Sa and Sb. As predicted, hemagglutination inhibition (HI) assays using human sera from subjects vaccinated with the initial pH1N1 isolate demonstrated reduced reactivity against 2013-2014 isolates. Taken together, these results suggest that diversifying selection analysis can identify key immune epitopes responsible for protective immunity to influenza virus in humans and thereby predict virus evolution.

IMPORTANCE

The WHO estimates that approximately 5 to 10% of adults and 20 to 30% of children in the world are infected by influenza virus each year. While an adaptive immune response helps eliminate the virus following acute infection, the virus rapidly evolves to evade the established protective memory immune response, thus allowing for the regular seasonal cycles of influenza virus infection. The analytical approach described here, which combines an analysis of diversifying selection with an integration of immune epitope data, has allowed us to identify antigenic regions that contribute to protective immunity and are therefore the key targets of immune evasion by the virus. This information can be used to determine when sequence variations in seasonal influenza virus strains have affected regions responsible for protective immunity in order to decide when new vaccine formulations are warranted.

Gene-wide identification of episodic selection

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

Molecular evolution of candidate genes involved in post-mating-prezygotic reproductive isolation

Traits involved in post-copulatory interactions between the sexes may evolve rapidly as a result of sexual selection and/or sexual conflict, leading to post-mating-prezygotic (PMPZ) reproductive isolating barriers between diverging lineages. Although the importance of PMPZ isolation is recognized, the molecular basis of such incompatibilities is not well understood. Here, we investigate molecular evolution of a subset of Drosophila mojavensis and Drosophila arizonae reproductive tract genes. These include genes that are transcriptionally regulated by conspecific mating in females, many of which are misregulated in heterospecific crosses, and a set of male genes whose transcripts are transferred to females during mating. As a group, misregulated female genes are not more divergent and do not appear to evolve under different selection pressures than other female reproductive genes. Male transferred genes evolve at a higher rate than testis-expressed genes, and at a similar rate compared to accessory gland protein genes, which are known to evolve rapidly. Four of the individual male transferred genes show patterns of divergent positive selection between D. mojavensis and D. arizonae. Three of the four genes belong to the sperm-coating protein-like family, including an ortholog of antares, which influences female fertility and receptivity in Drosophila melanogaster. Synthesis of these molecular evolutionary analyses with transcriptomics and predicted functional information makes these genes candidates for involvement in PMPZ reproductive incompatibilities between D. mojavensis and D. arizonae.

An evolutionarily-unique heterodimeric voltage-gated cation channel found in aphids

•

Aphids have a unique heterodimeric voltage-gated sodium channel.

•

The aphid channel has an atypical ion-selectivity filter (DENS rather than DEKA).

•

The channel’s novel selectivity filter may result in a loss of sodium selectivity.

•

This is the only identifiable voltage-gated sodium channel in aphid genome(s).

•

This channel has most likely arisen by gene fission or gene duplication.

We describe the identification in aphids of a unique heterodimeric voltage-gated sodium channel which has an atypical ion selectivity filter and, unusually for insect channels, is highly insensitive to tetrodotoxin. We demonstrate that this channel has most likely arisen by adaptation (gene fission or duplication) of an invertebrate ancestral mono(hetero)meric channel. This is the only identifiable voltage-gated sodium channel homologue in the aphid genome(s), and the channel’s novel selectivity filter motif (DENS instead of the usual DEKA found in other eukaryotes) may result in a loss of sodium selectivity, as indicated experimentally in mutagenised Drosophila channels.

The contrasting phylodynamics of human influenza B viruses

A complex interplay of viral, host, and ecological factors shapes the spatio-temporal incidence and evolution of human influenza viruses. Although considerable attention has been paid to influenza A viruses, a lack of equivalent data means that an integrated evolutionary and epidemiological framework has until now not been available for influenza B viruses, despite their significant disease burden. Through the analysis of over 900 full genomes from an epidemiological collection of more than 26,000 strains from Australia and New Zealand, we reveal fundamental differences in the phylodynamics of the two co-circulating lineages of influenza B virus (Victoria and Yamagata), showing that their individual dynamics are determined by a complex relationship between virus transmission, age of infection, and receptor binding preference. In sum, this work identifies new factors that are important determinants of influenza B evolution and epidemiology.

The 'emergence' of turnip mosaic virus was probably a 'gene-for-quasi-gene' event

Turnip mosaic potyvirus is a virus of brassicas that emerged from a lineage of monocotyledon-infecting potyviruses about 1000 years ago. In vivo and in silico studies all indicate that sites, primarily in its protein 3 (P3) and cylindrical inclusion protein (CI) genes, but also its small 6 kDa 2 protein (6K2) and genome-linked viral protein (VPg) genes, control host specificity in a dynamic way. It is most likely that non-unique combinations of transient viral genomic single nucleotide polymorphisms (SNPs), not all of them non-synonymous, allowed the host switch to occur. These SNPs were probably ephemeral and replaced over time by other combinations as the population subsequently diverged within, and adapted to, the brassica host population.

Phylodynamic analysis of clinical and environmental Vibrio cholerae isolates from Haiti reveals diversification driven by positive selection

Phylodynamic analysis of genome-wide single-nucleotide polymorphism (SNP) data is a powerful tool to investigate underlying evolutionary processes of bacterial epidemics. The method was applied to investigate a collection of 65 clinical and environmental isolates of Vibrio cholerae from Haiti collected between 2010 and 2012. Characterization of isolates recovered from environmental samples identified a total of four toxigenic V. cholerae O1 isolates, four non-O1/O139 isolates, and a novel nontoxigenic V. cholerae O1 isolate with the classical tcpA gene. Phylogenies of strains were inferred from genome-wide SNPs using coalescent-based demographic models within a Bayesian framework. A close phylogenetic relationship between clinical and environmental toxigenic V. cholerae O1 strains was observed. As cholera spread throughout Haiti between October 2010 and August 2012, the population size initially increased and then fluctuated over time. Selection analysis along internal branches of the phylogeny showed a steady accumulation of synonymous substitutions and a progressive increase of nonsynonymous substitutions over time, suggesting diversification likely was driven by positive selection. Short-term accumulation of nonsynonymous substitutions driven by selection may have significant implications for virulence, transmission dynamics, and even vaccine efficacy.

IMPORTANCE

Cholera, a dehydrating diarrheal disease caused by toxigenic strains of the bacterium Vibrio cholerae, emerged in 2010 in Haiti, a country where there were no available records on cholera over the past 100 years. While devastating in terms of morbidity and mortality, the outbreak provided a unique opportunity to study the evolutionary dynamics of V. cholerae and its environmental presence. The present study expands on previous work and provides an in-depth phylodynamic analysis inferred from genome-wide single nucleotide polymorphisms of clinical and environmental strains from dispersed geographic settings in Haiti over a 2-year period. Our results indicate that even during such a short time scale, V. cholerae in Haiti has undergone evolution and diversification driven by positive selection, which may have implications for understanding the global clinical and epidemiological patterns of the disease. Furthermore, the continued presence of the epidemic strain in Haitian aquatic environments has implications for transmission.

454 screening of individual MHC variation in an endemic island passerine

Genes of the major histocompatibility complex (MHC) code for receptors that are central to the adaptive immune response of vertebrates. These genes are therefore important genetic markers with which to study adaptive genetic variation in the wild. Next-generation sequencing (NGS) has increasingly been used in the last decade to genotype the MHC. However, NGS methods are highly prone to sequencing errors, and although several methodologies have been proposed to deal with this, until recently there have been no standard guidelines for the validation of putative MHC alleles. In this study, we used the 454 NGS platform to screen MHC class I exon 3 variation in a population of the island endemic Berthelot's pipit (Anthus berthelotii). We were able to characterise MHC genotypes across 309 individuals with high levels of repeatability. We were also able to determine alleles that had low amplification efficiencies, whose identification within individuals may thus be less reliable. At the population level we found lower levels of MHC diversity in Berthelot's pipit than in its widespread continental sister species the tawny pipit (Anthus campestris), and observed trans-species polymorphism. Using the sequence data, we identified signatures of gene conversion and evidence of maintenance of functionally divergent alleles in Berthelot's pipit. We also detected positive selection at 10 codons. The present study therefore shows that we have an efficient method for screening individual MHC variation across large datasets in Berthelot's pipit, and provides data that can be used in future studies investigating spatio-temporal patterns and scales of selection on the MHC.

Heterogeneous genetic diversity pattern in Plasmodium vivax genes encoding merozoite surface proteins (MSP) -7E, -7F and -7L

Background

The msp-7 gene has become differentially expanded in the Plasmodium genus; Plasmodium vivax has the highest copy number of this gene, several of which encode antigenic proteins in merozoites.

Methods

DNA sequences from thirty-six Colombian clinical isolates from P. vivax (pv) msp-7E, −7F and -7L genes were analysed for characterizing and studying the genetic diversity of these pvmsp-7 members which are expressed during the intra-erythrocyte stage; natural selection signals producing the variation pattern so observed were evaluated.

Results

The pvmsp-7E gene was highly polymorphic compared to pvmsp-7F and pvmsp-7L which were seen to have limited genetic diversity; pvmsp-7E polymorphism was seen to have been maintained by different types of positive selection. Even though these copies seemed to be species-specific duplications, a search in the Plasmodium cynomolgi genome (P. vivax sister taxon) showed that both species shared the whole msp-7 repertoire. This led to exploring the long-term effect of natural selection by comparing the orthologous sequences which led to finding signatures for lineage-specific positive selection.

Conclusions

The results confirmed that the P. vivax msp-7 family has a heterogeneous genetic diversity pattern; some members are highly conserved whilst others are highly diverse. The results suggested that the 3′-end of these genes encode MSP-7 proteins’ functional region whilst the central region of pvmsp-7E has evolved rapidly. The lineage-specific positive selection signals found suggested that mutations occurring in msp-7s genes during host switch may have succeeded in adapting the ancestral P. vivax parasite population to humans.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-495) contains supplementary material, which is available to authorized users.

Mitogenomics of recombinant mitochondrial genomes of Baltic Sea Mytilus mussels

Recombination in the control region (CR) of Mytilus mitochondrial DNA (mtDNA) was originally reported based on the relatively short, sequenced fragments of mitochondrial genomes. Recombination outside the CR has been reported recently with the suggestion that such processes are common in Mytilus. We have fully sequenced a set of 11 different mitochondrial haplotypes representing the high diversity of paternally inherited mitochondrial genomes of Baltic Sea Mytilus mussels, including the haplotype close to the native Mytilus trossulus mitochondrial genome, which was thought to have been entirely eliminated from this population. Phylogenetic and comparative analysis showed that the recombination is limited to the vicinity of the CR in all sequenced genomes. Coding sequence comparison indicated that all paternally inherited genomes showed increased accumulation of nonsynonymous substitutions, including the genomes which switched their transmission route very recently. The acquisition of certain CR sequences through recombination with highly divergent paternally inherited genomes seems to precede and favor the switch, but it is not a prerequisite for this process. Interspecies hybridization in the Baltic Sea during the recent 10,000 years created conditions for both structural and evolutionary mitochondrial instability which resulted in the observed variation and dynamics of mtDNA in Baltic Sea Mytilus mussels. In conclusion, the data shows that the effects of mitochondrial recombination are limited to the CR of few phylogenetic lineages.

Molecular diversity of Chickpea chlorotic dwarf virus in Sudan: high rates of intra-species recombination - a driving force in the emergence of new strains

In Sudan Chickpea chlorotic dwarf virus (CpCDV, genus Mastrevirus, family Geminiviridae) is an important pathogen of pulses that are grown both for local consumption, and for export. Although a few studies have characterised CpCDV genomes from countries in the Middle East, Africa and the Indian subcontinent, little is known about CpCDV diversity in any of the major chickpea production areas in these regions. Here we analyse the diversity of 146 CpCDV isolates characterised from pulses collected across the chickpea growing regions of Sudan. Although we find that seven of the twelve known CpCDV strains are present within the country, strain CpCDV-H alone accounted for ∼73% of the infections analysed. Additionally we identified four new strains (CpCDV-M, -N, -O and -P) and show that recombination has played a significant role in the diversification of CpCDV, at least in this region. Accounting for observed recombination events, we use the large amounts of data generated here to compare patterns of natural selection within protein coding regions of CpCDV and other dicot-infecting mastrevirus species.

Whole genome phylogenetic investigation of a West Nile virus strain isolated from a tick sampled from livestock in north eastern Kenya

Background

West Nile virus (WNV) has a wide geographical distribution and has been associated to cause neurological disease in humans and horses. Mosquitoes are the traditional vectors for WNV; however, the virus has also been isolated from tick species in North Africa and Europe which could be a means of introduction and spread of the virus over long distances through migratory birds. Although WNV has been isolated in mosquitoes in Kenya, paucity of genetic and pathogenicity data exists. We previously reported the isolation of WNV from ticks collected from livestock and wildlife in Ijara District of Kenya, a hotspot for arbovirus activity. Here we report the full genome sequence and phylogenetic investigation of their origin and relation to strains from other regions.

Methods

A total of 10,488 ticks were sampled from animal hosts, classified to species and processed in pools of up to eight ticks per pool. Virus screening was performed by cell culture, RT-PCR and sequencing. Phylogenetic analysis was carried out to determine the evolutionary relationships of our isolate.

Results

Among other viruses, WNV was isolated from a pool of Rhipicephalus pulchellus sampled from cattle, sequenced and submitted to GenBank (Accession number: KC243146). Comparative analysis with 27 different strains revealed that our isolate belongs to lineage 1 and clustered relatively closely to isolates from North Africa and Europe, Russia and the United States. Overall, Bayesian analysis based on nucleotide sequences showed that lineage 1 strains including the Kenyan strain had diverged 200 years ago from lineage 2 strains of southern Africa. Ijara strain collected from a tick sampled on livestock was closest to another Kenyan strain and had diverged 20 years ago from strains detected in Morocco and Europe and 30 years ago from strains identified in the USA.

Conclusion

To our knowledge, this is the first characterized WNV strain isolated from R. pulchellus. The epidemiological role of this tick in WNV transmission and dissemination remains equivocal but presents tick verses mosquito virus transmission has been neglected. Genetic data of this strain suggest that lineage 1 strains from Africa could be dispersed through tick vectors by wild migratory birds to Europe and beyond.

Electronic supplementary material

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

Next generation sequencing improves detection of drug resistance mutations in infants after PMTCT failure

BACKGROUND

Next generation sequencing (NGS) allows the detection of minor variant HIV drug resistance mutations (DRMs). However data from new NGS platforms after Prevention-of-Mother-to-Child-Transmission (PMTCT) regimen failure are limited.

OBJECTIVE

To compare major and minor variant HIV DRMs with Illumina MiSeq and Life Technologies Ion Personal Genome Machine (PGM) in infants infected despite a PMTCT regimen.

STUDY DESIGN

We conducted a cross-sectional study of NGS for detecting DRMs in infants infected despite a zidovudine (AZT) and Nevirapine (NVP) regimen, before initiation of combination antiretroviral therapy. Sequencing was performed on PCR products from plasma samples on PGM and MiSeq platforms. Bioinformatic analyses were undertaken using a codon-aware version of the Smith-Waterman mapping algorithm and a mixture multinomial error filtering statistical model.

RESULTS

Of 15 infants, tested at a median age of 3.4 months after birth, 2 (13%) had non-nucleoside reverse transcriptase inhibitor (NNRTI) DRMs (K103N and Y181C) by bulk sequencing, whereas PGM detected 4 (26%) and MiSeq 5 (30%). NGS enabled the detection of additional minor variant DRMs in the infant with K103N. Coverage and instrument quality scores were higher with MiSeq, increasing the confidence of minor variant calls.

CONCLUSIONS

NGS followed by bioinformatic analyses detected multiple minor variant DRMs in HIV-1 RT among infants where PMTCT failed. The high coverage of MiSeq and high read quality improved the confidence of identified DRMs and may make this platform ideal for minor variant detection.

Genetic and epidemiological insights into the emergence of peste des petits ruminants virus (PPRV) across Asia and Africa

Small ruminants are important components in the livelihood of millions of households in many parts of the world. The spread of the highly contagious peste des petits ruminants (PPR) disease, which is caused by an RNA virus, PPRV, across Asia and Africa remains a major concern. The present study explored the evolutionary and epidemiological dynamics of PPRV through the analyses of partial N-gene and F-gene sequences of the virus. All the four previously described PPRV lineages (I-IV) diverged from their common ancestor during the late-19(th) to early-20(th) century. Among the four lineages, PPRV-IV showed pronounced genetic structuring across the region; however, haplotype sharing among the geographic regions, together with the presence of multiple genetic clusters within a country, indicates the possibility of frequent mobility of the diseased individuals across the region. The gradual decline in the effective number of infections suggests a limited genetic variation, which could be attributed to the effective vaccination that has been practiced since 1990s. However, the movement of infected animals across the region likely contributes to the spread of PPRV-IV. No evidence of positive selection was identified from this study.

Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds

BACKGROUND

Toll-like receptors (TLR) are essential activators of the innate part of the vertebrate immune system. In this study, we analysed the interspecific variability of three TLR (bacterial-sensing TLR4 and TLR5 and viral-sensing TLR7) within the Galloanserae bird clade, investigated their phylogeny, assessed their structural conservation and estimated site-specific selection pressures.

RESULTS

Physiochemical properties varied according to the TLR analysed, mainly with regards to the surface electrostatic potential distribution. The predicted ligand-binding features (mainly in TLR4 and TLR5) differed between the avian proteins and their fish and mammalian counterparts, but also varied within the Galloanserae birds. We identified 20 positively selected sites in the three TLR, among which several are topologically close to ligand-binding sites reported for mammalian and fish TLR. We described 26, 28 and 25 evolutionarily non-conservative sites in TLR4, TLR5 and TLR7, respectively. Thirteen of these sites in TLR4, and ten in TLR5 were located in functionally relevant regions. The variability appears to be functionally more conserved for viral-sensing TLR7 than for the bacterial-sensing TLR. Amino-acid positions 268, 270, 343, 383, 444 and 471 in TLR4 and 180, 183, 209, 216, 264, 342 and 379 in TLR5 are key candidates for further functional research.

CONCLUSIONS

Host-pathogen co-evolution has a major effect on the features of host immune receptors. Our results suggest that avian and mammalian TLR may be differentially adapted to pathogen-derived ligand recognition. We have detected signatures of positive selection even within the Galloanserae lineage. To our knowledge, this is the first study to depict evolutionary pressures on Galloanserae TLR and to estimate the validity of current knowledge on TLR function (based on mammalian and chicken models) for non-model species of this clade.

Evolution and Structural Analyses of Glossina morsitans (Diptera; Glossinidae) Tetraspanins

Tetraspanins are important conserved integral membrane proteins expressed in many organisms. Although there is limited knowledge about the full repertoire, evolution and structural characteristics of individual members in various organisms, data obtained so far show that tetraspanins play major roles in membrane biology, visual processing, memory, olfactory signal processing, and mechanosensory antennal inputs. Thus, these proteins are potential targets for control of insect pests. Here, we report that the genome of the tsetse fly, Glossina morsitans (Diptera: Glossinidae) encodes at least seventeen tetraspanins (GmTsps), all containing the signature features found in the tetraspanin superfamily members. Whereas six of the GmTsps have been previously reported, eleven could be classified as novel because their amino acid sequences do not map to characterized tetraspanins in the available protein data bases. We present a model of the GmTsps by using GmTsp42Ed, whose presence and expression has been recently detected by transcriptomics and proteomics analyses of G. morsitans. Phylogenetically, the identified GmTsps segregate into three major clusters. Structurally, the GmTsps are largely similar to vertebrate tetraspanins. In view of the exploitation of tetraspanins by organisms for survival, these proteins could be targeted using specific antibodies, recombinant large extracellular loop (LEL) domains, small-molecule mimetics and siRNAs as potential novel and efficacious putative targets to combat African trypanosomiasis by killing the tsetse fly vector.

Adaptive evolution of a novel avian-origin influenza A/H7N9 virus

In China, the recent outbreak of novel influenza A/H7N9 virus has been assumed to be severe, and it may possibly turn brutal in the near future. In order to develop highly protective vaccines and drugs for the A/H7N9 virus, it is critical to find out the selection pressure of each amino acid site. In the present study, six different statistical methods consisting of four independent codon-based maximum likelihood (CML) methods, one hierarchical Bayesian (HB) method and one branch-site (BS) method, were employed to determine if each amino acid site of A/H7N9 virus is under natural selection pressure. Functions for both positively and negatively selected sites were inferred by annotating these sites with experimentally verified amino acid sites. Comprehensively, the single amino acid site 627 of PB2 protein was inferred as positively selected and it function was identified as a T-cell epitope (TCE). Among the 26 negatively selected amino acid sites of PB2, PB1, PA, HA, NP, NA, M1 and NS2 proteins, only 16 amino acid sites were identified to be involved in TCEs. In addition, 7 amino acid sites including, 608 and 609 of PA, 480 of NP, and 24, 25, 109 and 205 of M1, were identified to be involved in both B-cell epitopes (BCEs) and TCEs. Conversely, the function of positions 62 of PA, and, 43 and 113 of HA was unknown. In conclusion, the seven amino acid sites engaged in both BCEs and TCEs were identified as highly suitable targets, as these sites will be predicted to play a principal role in inducing strong humoral and cellular immune responses against A/H7N9 virus.

IMPACT_S: integrated multiprogram platform to analyze and combine tests of selection

Among the major goals of research in evolutionary biology are the identification of genes targeted by natural selection and understanding how various regimes of evolution affect the fitness of an organism. In particular, adaptive evolution enables organisms to adapt to changing ecological factors such as diet, temperature, habitat, predatory pressures and prey abundance. An integrative approach is crucial for the identification of non-synonymous mutations that introduce radical changes in protein biochemistry and thus in turn influence the structure and function of proteins. Performing such analyses manually is often a time-consuming process, due to the large number of statistical files generated from multiple approaches, especially when assessing numerous taxa and/or large datasets. We present IMPACT_S, an easy-to-use Graphical User Interface (GUI) software, which rapidly and effectively integrates, filters and combines results from three widely used programs for assessing the influence of selection: Codeml (PAML package), Datamonkey and TreeSAAP. It enables the identification and tabulation of sites detected by these programs as evolving under the influence of positive, neutral and/or negative selection in protein-coding genes. IMPACT_S further facilitates the automatic mapping of these sites onto the three-dimensional structures of proteins. Other useful tools incorporated in IMPACT_S include Jmol, Archaeopteryx, Gnuplot, PhyML, a built-in Swiss-Model interface and a PDB downloader. The relevance and functionality of IMPACT_S is shown through a case study on the toxicoferan-reptilian Cysteine-rich Secretory Proteins (CRiSPs). IMPACT_S is a platform-independent software released under GPLv3 license, freely available online from http://impact-s.sourceforge.net.

Molecular evolutionary and epidemiological dynamics of genotypes 1G and 2B of rubella virus

Rubella Virus (RV), which causes measles-like rashes in children, puts millions of infants at risk of congenital defects across the globe. Employing phylogenetic approaches to the whole genome sequence data and E1 glycoprotein sequence data, the present study reports the substitution rates and dates of emergence of all thirteen previously described rubella genotypes, and gains important insights into the epidemiological dynamics of two geographically widely distributed genotypes 1G and 2B. The overall nucleotide substitution rate of this non-vector-borne RV is in the order of 10-3 substitutions/site/year, which is considerably higher than the substitution rates previously reported for the vector-borne alphaviruses within the same family. Currently circulating strains of RV share a common ancestor that existed within the last 150 years, with 95% Highest Posterior Density values ranging from 1868 to 1926 AD. Viral strains within the respective genotypes began diverging between the year 1930 s and 1980 s. Both genotype 1G and 2B have shown a decline in effective number of infections since 1990 s, a period during which mass immunization programs against RV were adapted across the globe. Although both genotypes showed some extent of spatial genetic structuring, the analyses also depicted an inter-continental viral dispersal. Such a viral dispersal pattern could be related to the migration of infected individuals across the regions coupled with a low coverage of MMR vaccination.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

The loss of taste genes in cetaceans

BACKGROUND

Five basic taste modalities, sour, sweet, bitter, salt and umami, can be distinguished by humans and are fundamental for physical and ecological adaptations in mammals. Molecular genetic studies of the receptor genes for these tastes have been conducted in terrestrial mammals; however, little is known about the evolution and adaptation of these genes in marine mammals.

RESULTS

Here, all five basic taste modalities, sour, sweet, bitter, salt and umami, were investigated in cetaceans. The sequence characteristics and evolutionary analyses of taste receptor genes suggested that nearly all cetaceans may have lost all taste modalities except for that of salt.

CONCLUSIONS

This is the first study to comprehensively examine the five basic taste modalities in cetaceans with extensive taxa sampling. Our results suggest that cetaceans have lost four of the basic taste modalities including sour, sweet, umami, and most of the ability to sense bitter tastes. The integrity of the candidate salt taste receptor genes in all the cetaceans examined may be because of their function in Na(+) reabsorption, which is key to osmoregulation and aquatic adaptation.

The genetics of monarch butterfly migration and warning colouration

The monarch butterfly, Danaus plexippus, is famous for its spectacular annual migration across North America, recent worldwide dispersal, and orange warning coloration. Despite decades of study and broad public interest, we know little about the genetic basis of these hallmark traits. By sequencing 101 monarch genomes from around the globe, we uncover the history of the monarch's evolutionary origin and global dispersal, characterize the genes and pathways associated with migratory behavior, and identify the discrete genetic basis of warning coloration. The results rewrite our understanding of this classic system, showing that D. plexippus was ancestrally migratory and dispersed out of North America to occupy its broad distribution. We find the strongest signatures of selection associated with migration center on flight muscle function, resulting in greater flight efficiency among migratory monarchs, and that variation in monarch warning coloration is controlled by a single myosin gene not previously implicated in insect pigmentation.

Deep sequencing of norovirus genomes defines evolutionary patterns in an urban tropical setting

Norovirus is a highly transmissible infectious agent that causes epidemic gastroenteritis in susceptible children and adults. Norovirus infections can be severe and can be initiated from an exceptionally small number of viral particles. Detailed genome sequence data are useful for tracking norovirus transmission and evolution. To address this need, we have developed a whole-genome deep-sequencing method that generates entire genome sequences from small amounts of clinical specimens. This novel approach employs an algorithm for reverse transcription and PCR amplification primer design using all of the publically available norovirus sequence data. Deep sequencing and de novo assembly were used to generate norovirus genomes from a large set of diarrheal patients attending three hospitals in Ho Chi Minh City, Vietnam, over a 2.5-year period. Positive-selection analysis and direct examination of protein changes in the virus over time identified codons in the regions encoding proteins VP1, p48 (NS1-2), and p22 (NS4) under positive selection and expands the known targets of norovirus evolutionary pressure.

IMPORTANCE

The high transmissibility and rapid evolutionary rate of norovirus, combined with a short-lived host immune responses, are thought to be the reasons why the virus causes the majority of pediatric viral diarrhea cases. The evolutionary patterns of this RNA virus have been described in detail for only a portion of the virus genome and never for a virus from a detailed urban tropical setting. We provide a detailed sequence description of the noroviruses circulating in three Ho Chi Minh City hospitals over a 2.5-year period. This study identified patterns of virus change in known sites of host immune response and identified three additional regions of the virus genome under selection that were not previously recognized. In addition, the method described here provides a robust full-genome sequencing platform for community-based virus surveillance.

Low genetic diversity in the locus encoding the Plasmodium vivax P41 protein in Colombia's parasite population

Background

The development of malaria vaccine has been hindered by the allele-specific responses produced by some parasite antigens’ high genetic diversity. Such antigen genetic diversity must thus be evaluated when designing a completely effective vaccine. Plasmodium falciparum P12, P38 and P41 proteins have red blood cell binding regions in the s48/45 domains and are located on merozoite surface, P41 forming a heteroduplex with P12. These three genes have been identified in Plasmodium vivax and share similar characteristics with their orthologues in Plasmodium falciparum. Plasmodium vivax pv12 and pv38 have low genetic diversity but pv41 polymorphism has not been described.

Methods

The present study was aimed at evaluating the P. vivax p41 (pv41) gene’s polymorphism. DNA sequences from Colombian clinical isolates from pv41 gene were analysed for characterising and studying the genetic diversity and the evolutionary forces that produced the variation pattern so observed.

Results

Similarly to other members of the 6-Cys family, pv41 had low genetic polymorphism. pv41 3′-end displayed the highest nucleotide diversity value; several substitutions found there were under positive selection. Negatively selected codons at inter-species level were identified in the s48/45 domains; p41 would thus seem to have functional/structural constraints due to the presence of these domains.

Conclusions

In spite of the functional constraints of Pv41 s48/45 domains, immune system pressure seems to have allowed non-synonymous substitutions to become fixed within them as an adaptation mechanism; including Pv41 s48/45 domains in a vaccine should thus be carefully evaluated due to these domains containing some allele variants.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-388) contains supplementary material, which is available to authorized users.

The influence of secondary structure, selection and recombination on rubella virus nucleotide substitution rate estimates

BACKGROUND

Annually, rubella virus (RV) still causes severe congenital defects in around 100 000 children globally. An attempt to eradicate RV is currently underway and analytical tools to monitor the global decline of the last remaining RV lineages will be useful for assessing the effectiveness of this endeavour. RV evolves rapidly enough that much of this information might be inferable from RV genomic sequence data.

METHODS

Using BEASTv1.8.0, we analysed publically available RV sequence data to estimate genome-wide and gene-specific nucleotide substitution rates to test whether current estimates of RV substitution rates are representative of the entire RV genome. We specifically accounted for possible confounders of nucleotide substitution rate estimates, such as temporally biased sampling, sporadic recombination, and natural selection favouring either increased or decreased genetic diversity (estimated by the PARRIS and FUBAR methods), at nucleotide sites within the genomic secondary structures (predicted by the NASP method).

RESULTS

We determine that RV nucleotide substitution rates range from 1.19 × 10(-3) substitutions/site/year in the E1 region to 7.52 × 10(-4) substitutions/site/year in the P150 region. We find that differences between substitution rate estimates in different RV genome regions are largely attributable to temporal sampling biases such that datasets containing higher proportions of recently sampled sequences, will tend to have inflated estimates of mean substitution rates. Although there exists little evidence of positive selection or natural genetic recombination in RV, we show that RV genomes possess pervasive biologically functional nucleic acid secondary structure and that purifying selection acting to maintain this structure contributes substantially to variations in estimated nucleotide substitution rates across RV genomes.

CONCLUSION

Both temporal sampling biases and purifying selection favouring the conservation of RV nucleic acid secondary structures have an appreciable impact on substitution rate estimates but do not preclude the use of RV sequence data to date ancestral sequences. The combination of uniformly high substitution rates across the RV genome and strong temporal structure within the available sequence data, suggests that such data should be suitable for tracking the demographic, epidemiological and movement dynamics of this virus during eradication attempts.

Positive selection and multiple losses of the LINE-1-derived L1TD1 gene in mammals suggest a dual role in genome defense and pluripotency

Mammalian genomes comprise many active and fossilized retroelements. The obligate requirement for retroelement integration affords host genomes an opportunity to 'domesticate' retroelement genes for their own purpose, leading to important innovations in genome defense and placentation. While many such exaptations involve retroviruses, the L1TD1 gene is the only known domesticated gene whose protein-coding sequence is almost entirely derived from a LINE-1 (L1) retroelement. Human L1TD1 has been shown to play an important role in pluripotency maintenance. To investigate how this role was acquired, we traced the origin and evolution of L1TD1. We find that L1TD1 originated in the common ancestor of eutherian mammals, but was lost or pseudogenized multiple times during mammalian evolution. We also find that L1TD1 has evolved under positive selection during primate and mouse evolution, and that one prosimian L1TD1 has 'replenished' itself with a more recent L1 ORF1 from the prosimian genome. These data suggest that L1TD1 has been recurrently selected for functional novelty, perhaps for a role in genome defense. L1TD1 loss is associated with L1 extinction in several megabat lineages, but not in sigmodontine rodents. We hypothesize that L1TD1 could have originally evolved for genome defense against L1 elements. Later, L1TD1 may have become incorporated into pluripotency maintenance in some lineages. Our study highlights the role of retroelement gene domestication in fundamental aspects of mammalian biology, and that such domesticated genes can adopt different functions in different lineages.

An evolutionary analysis of the Secoviridae family of viruses

The plant-infecting Secoviridae family of viruses forms part of the Picornavirales order, an important group of non-enveloped viruses that infect vertebrates, arthropods, plants and algae. The impact of the secovirids on cultivated crops is significant, infecting a wide range of plants from grapevine to rice. The overwhelming majority are transmitted by ecdysozoan vectors such as nematodes, beetles and aphids. In this study, we have applied a variety of computational methods to examine the evolutionary traits of these viruses. Strong purifying selection pressures were calculated for the coat protein (CP) sequences of nine species, although for two species evidence of both codon specific and episodic diversifying selection were found. By using Bayesian phylogenetic reconstruction methods CP nucleotide substitution rates for four species were estimated to range from between 9.29×10⁻³ to 2.74×10⁻³ (subs/site/year), values which are comparable with the short-term estimates of other related plant- and animal-infecting virus species. From these data, we were able to construct a time-measured phylogeny of the subfamily Comovirinae that estimated divergence of ninety-four extant sequences occurred less than 1,000 years ago with present virus species diversifying between 50 and 250 years ago; a period coinciding with the intensification of agricultural practices in industrial societies. Although recombination (modularity) was limited to closely related taxa, significant and often unique similarities in the protein domains between secovirid and animal infecting picorna-like viruses, especially for the protease and coat protein, suggested a shared ancestry. We discuss our results in a wider context and find tentative evidence to indicate that some members of the Secoviridae might have their origins in insects, possibly colonizing plants in a number of founding events that have led to speciation. Such a scenario; virus infection between species of different taxonomic kingdoms, has significant implications for virus emergence.

Key mutations in the cylindrical inclusion involved in lettuce mosaic virus adaptation to eIF4E-mediated resistance in lettuce

We previously showed that allelic genes mol¹ and mo1² used to protect lettuce crops against Lettuce mosaic virus (LMV) correspond to mutant alleles of the gene encoding the eukaryotic translation initiation factor 4E. LMV resistance-breaking determinants map not only to the main potyvirus virulence determinant, a genome-linked viral protein, but also to the C-terminal region of the cylindrical inclusion (CI), with a key role of amino acid at position 621. Here, we show that the propagation of several non-lettuce isolates of LMV in mo1¹ plants is accompanied by a gain of virulence correlated with the presence in the CI C terminus of a serine at position 617 and the accumulation of mutations at positions 602 or 627. Whole-genome sequencing of native and evolved isolates showed that no other mutation could be associated with adaptation to mo1 resistance. Site-directed mutagenesis pinpointed the key role in the virulence of the combination of mutations at positions 602 and 617, in addition to position 621. The impact of these mutations on the fitness of the virus was evaluated, suggesting that the durability of mo1 resistance in the field relies on the fitness cost associated with the resistance-breaking mutations, the nature of the mutations, and their potential antagonistic effects.

Increased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (Oryctolagus cuniculus)

The release of myxoma virus (MYXV) and Rabbit Haemorrhagic Disease Virus (RHDV) in Australia with the aim of controlling overabundant rabbits has provided a unique opportunity to study the initial spread and establishment of emerging pathogens, as well as their co-evolution with their mammalian hosts. In contrast to MYXV, which attenuated shortly after its introduction, rapid attenuation of RHDV has not been observed. By studying the change in virulence of recent field isolates at a single field site we show, for the first time, that RHDV virulence has increased through time, likely because of selection to overcome developing genetic resistance in Australian wild rabbits. High virulence also appears to be favoured as rabbit carcasses, rather than diseased animals, are the likely source of mechanical insect transmission. These findings not only help elucidate the co-evolutionary interaction between rabbits and RHDV, but reveal some of the key factors shaping virulence evolution.

The genetics of host-virus coevolution in invertebrates

Although viral infection and antiviral defence are ubiquitous, genetic data are currently unavailable from the vast majority of animal phyla-potentially biasing our overall perspective of the coevolutionary process. Rapid adaptive evolution is seen in some insect antiviral genes, consistent with invertebrate-virus 'arms-race' coevolution, but equivalent signatures of selection are hard to detect in viruses. We find that, despite the large differences in vertebrate, invertebrate, and plant immune responses, comparison of viral evolution fails to identify any difference among these hosts in the impact of positive selection. The best evidence for invertebrate-virus coevolution is currently provided by large-effect polymorphisms for host resistance and/or viral evasion, as these often appear to have arisen and spread recently, and can be favoured by virus-mediated selection.

Genetic differentiation among isolates of Teredinibacter turnerae, a widely occurring intracellular endosymbiont of shipworms

Teredinibacter turnerae is a cultivable intracellular endosymbiont of xylotrophic (woodfeeding)bivalves of the Family Teredinidae (shipworms). Although T. turnerae has been isolated from many shipworm taxa collected in many locations, no systematic effort has been made to explore genetic diversity within this symbiont species across the taxonomic and geographical range of its hosts. The mode of symbiont transmission is unknown. Here, we examine sequence diversity in fragments of six genes (16S rRNA, gyrB, sseA, recA, rpoB and celAB) among 25 isolates of T. turnerae cultured from 13 shipworm species collected in 15 locations in the Atlantic, Pacific and Indian Oceans. While 16S rRNA sequences are nearly invariant between all examined isolates (maximum pairwise difference <0.26%), variation between examined protein-coding loci is greater (mean pairwise difference 2.2–5.9%). Phylogenetic analyses based on each protein-coding locus differentiate the 25 isolates into two distinct and well-supported clades. With five exceptions, clade assignments for each isolate were supported by analysis of alleles of each of the five protein-coding loci. These exceptions include (i) putative recombinant alleles of the celAB and gyrB loci in two isolates (PMS-535T.S.1b.3 and T8510), suggesting homologous recombination between members of the two clades; and (ii) evidence for a putative lateral gene transfer event affecting a second locus (recA) in three isolates (T8412, T8503 and T8513). These results demonstrate that T. turnerae isolates do not represent a homogeneous global population. Instead, they indicate the emergence of two lineages that, although distinct, likely experience some level of genetic exchange with each other and with other bacterial species.

Molecular evolutionary history of Sugarcane yellow leaf virus based on sequence analysis of RNA-dependent RNA polymerase and putative aphid transmission factor-coding genes

RNA-dependent RNA polymerase (RdRp) encoded by ORF2 and putative aphid transmission factor (PATF) encoded by ORF5 of Sugarcane yellow leaf virus (SCYLV) were detected in six sugarcane cultivars affected by yellow leaf using RT-PCR and real-time RT-PCR assays. Expression of both genes varied among infected plants, but overall expression of RdRp was higher than expression of PATF. Cultivar H87-4094 from Hawaii yielded the highest transcript levels of RdRp, whereas cultivar C1051-73 from Cuba exhibited the lowest levels. Sequence comparisons among 25 SCYLV isolates from various geographical locations revealed an amino acid similarity of 72.1-99.4 and 84.7-99.8 % for the RdRp and PATF genes, respectively. The 25 SCYLV isolates were separated into three (RdRp) and two (PATF) phylogenetic groups using the MEGA6 program that does not account for genetic recombination. However, the SCYLV genome contained potential recombination signals in the RdRp and PATF coding genes based on the GARD genetic algorithm. Use of this later program resulted in the reconstruction of phylogenies on the left as well as on the right sides of the putative recombination breaking points, and the 25 SCYLV isolates were distributed into three distinct phylogenetic groups based on either RdRp or PATF sequences. As a result, recombination reshuffled the affiliation of the accessions to the different clusters. Analysis of selection pressures exerted on RdRp and PATF encoded proteins revealed that ORF 2 and ORF 5 underwent predominantly purifying selection. However, a few sites were also under positive selection as assessed by various models such as FEL, IFEL, REL, FUBAR, MEME, GA-Branch, and PRIME.

DNA sequence polymorphisms of the pfmdr1 gene and association of mutations with the pfcrt gene in Indian Plasmodium falciparum isolates

Mutations in the Plasmodium falciparum multidrug resistance (pfmdr1) gene are known to provide compensatory fitness benefits to the chloroquine (CQ)-resistant malaria parasites and are often associated with specific mutations in the P. falciparum CQ resistant transporter (pfcrt) gene. Prevalence of the specific mutations in these two genes across different malaria endemic regions was mostly studies. However, reports on mutations in the pfmdr1 gene and their genetic associations with mutations in the pfcrt gene in Indian P. falciparum field isolates are scarce. We have sequenced a 560 bp region of pfmdr1 coding sequence in 64 P. falciparum isolates collected from different malaria endemic populations in India. Twenty out of these 64 isolates were laboratory cultured with known in vitro CQ sensitiveness (10 sensitive and 10 resistant). Three low frequency mutations (two non-synonymous and one synonymous) in the pfmdr1 gene were segregating in Indian isolates in addition to the predominant Y₈₆ and Y₁₈₄ ones, with high haplotype and nucleotide diversity in the field isolates in comparison to the cultured ones. No statistically significant genetic association between the mutations in the pfmdr1 and pfcrt gene could be detected; almost all observed associations were intragenic in nature. The results on the genetic diversity of the pfmdr1 gene were discussed in term of evolutionary perspectives in Indian P. falciparum, with possible future potential of gaining further insights on this gene in view of evolving malaria parasites resistant to artemisinin partner drugs.

Western equine encephalitis virus: evolutionary analysis of a declining alphavirus based on complete genome sequences

Western equine encephalitis virus (WEEV) is an arbovirus from the genus Alphavirus, family Togaviridae, which circulates in North America between birds and mosquitoes, occasionally causing disease in humans and equids. In recent decades, human infection has decreased dramatically; the last documented human case in North America occurred in 1994, and the virus has not been detected in mosquito pools since 2008. Because limited information exists regarding the evolution of WEEV, we analyzed the genomic sequences of 33 low-passage-number strains with diverse geographic and temporal distributions and performed comprehensive phylogenetic analyses. Our results indicated that WEEV is a highly conserved alphavirus with only approximately 5% divergence in its most variable genes. We confirmed the presence of the previously determined group A and B lineages and further resolved group B into three sublineages. We also observed an increase in relative genetic diversity during the mid-20th century, which correlates with the emergence and cocirculation of several group B sublineages. The estimated WEEV population size dropped in the 1990s, with only the group B3 lineage being sampled in the past 20 years. Structural mapping showed that the majority of substitutions in the envelope glycoproteins occurred at the E2-E2 interface. We hypothesize that an event occurred in the mid-20th century that resulted in the increased genetic diversity of WEEV in North America, followed by genetic constriction due to either competitive displacement by the B3 sublineage or stochastic events resulting from a population decline.

IMPORTANCE

Western equine encephalitis virus (WEEV) has caused several epidemics that resulted in the deaths of thousands of humans and hundreds of thousands of equids during the past century. During recent decades, human infection decreased drastically and the virus has not been found in mosquito pools since 2008. Because limited information exists regarding the evolution of WEEV, we analyzed 33 complete genome sequences and conducted comprehensive phylogenetic analyses. We confirmed the presence of two major lineages, one of which diverged into three sublineages. Currently, only one of those sublineages is found circulating in nature. Understanding the evolution of WEEV over the past century provides a unique opportunity to observe an arbovirus that is in decline and to better understand what factors can cause said decline.

Limited utility of residue masking for positive-selection inference

Errors in multiple sequence alignments (MSAs) can reduce accuracy in positive-selection inference. Therefore, it has been suggested to filter MSAs before conducting further analyses. One widely used filter, Guidance, allows users to remove MSA positions aligned with low confidence. However, Guidance's utility in positive-selection inference has been disputed in the literature. We have conducted an extensive simulation-based study to characterize fully how Guidance impacts positive-selection inference, specifically for protein-coding sequences of realistic divergence levels. We also investigated whether novel scoring algorithms, which phylogenetically corrected confidence scores, and a new gap-penalization score-normalization scheme improved Guidance's performance. We found that no filter, including original Guidance, consistently benefitted positive-selection inferences. Moreover, all improvements detected were exceedingly minimal, and in certain circumstances, Guidance-based filters worsened inferences.

Speciation and demographic history of Atlantic eels (Anguilla anguilla and A. rostrata) revealed by mitogenome sequencing

Processes leading to speciation in oceanic environments without obvious physical barriers remain poorly known. European and American eel (Anguilla anguilla and A. rostrata) spawn in partial sympatry in the Sargasso Sea. Larvae are advected by the Gulf Stream and other currents towards the European/North African and North American coasts, respectively. We analyzed 104 mitogenomes from the two species along with mitogenomes of other Anguilla and outgroup species. We estimated divergence time between the two species to identify major events involved in speciation. We also considered two previously stated hypotheses: one where the ancestral species was present in only one continent but was advected across the Atlantic by ocean current changes and another where population declines during Pleistocene glaciations led to increasing vicariance, facilitating speciation. Divergence time was estimated to ∼3.38 Mya, coinciding with the closure of the Panama Gateway that led to reinforcement of the Gulf Stream. This could have advected larvae towards European/North African coasts, in which case American eel would be expected to be the ancestral species. This scenario could, however, not be unequivocally confirmed by analyses of dN/dS, nucleotide diversity and effective population size estimates. Extended bayesian skyline plots showed fluctuations of effective population sizes and declines during glaciations, and thus also lending support to the importance of vicariance during speciation. There was evidence for positive selection at the ATP6 and possibly ND5 genes, indicating a role in speciation. The findings suggest an important role of ocean current changes in speciation of marine organisms.

Discovery and evolution of bunyavirids in arctic phantom midges and ancient bunyavirid-like sequences in insect genomes

Bunyaviridae is a large family of RNA viruses chiefly comprised of vertebrate and plant pathogens. We discovered novel bunyavirids that are approximately equally divergent from each of the five known genera. We characterized novel genome sequences for two bunyavirids, namely, Kigluaik phantom virus (KIGV), from tundra-native phantom midges (Chaoborus), and Nome phantom virus (NOMV), from tundra-invading phantom midges, and demonstrated that these bunyavirid-like sequences belong to an infectious virus by passaging KIGV in mosquito cell culture, although the infection does not seem to be well sustained beyond a few passages. Virus and host gene sequences from individuals collected on opposite ends of North America, a region spanning 4,000 km, support a long-term, vertically transmitted infection of KIGV in Chaoborus trivittatus. KIGV-like sequences ranging from single genes to full genomes are present in transcriptomes and genomes of insects belonging to six taxonomic orders, suggesting an ancient association of this clade with insect hosts. In Drosophila, endogenous virus genes have been coopted, forming an orthologous tandem gene family that has been maintained by selection during the radiation of the host genus. Our findings indicate that bunyavirid-host interactions in nonbloodsucking arthropods have been much more extensive than previously thought.

IMPORTANCE

Very little is known about the viral diversity in polar freshwater ponds, and perhaps less is known about the effects that climate-induced habitat changes in these regions will have on virus-host interactions in the coming years. Our results show that at the tundra-boreal boundary, a hidden viral landscape is being altered as infected boreal phantom midges colonize tundra ponds. Likewise, relatively little is known of the deeper evolutionary history of bunyavirids that has led to the stark lifestyle contrasts between some genera. The discovery of this novel bunyavirid group suggests that ancient and highly divergent bunyavirid lineages remain undetected in nature and may offer fresh insight into host reservoirs, potential sources of emerging disease, and major lifestyle shifts in the evolutionary history of viruses in the family Bunyaviridae.

Characterization and expression of the CXCR1 and CXCR4 in miiuy croaker and evolutionary analysis shows the strong positive selection pressures imposed in mammal CXCR1

The innate immune system can recognize non-self, danger signals, and pathogen associated molecular patterns and provides a first line of antimicrobial host defense. Therefore, it plays an instructive role and is pretty important in vertebrates. In innate immune responses, CXCRs act as the main receptors of CXC chemokines and play a vital role in host defense and inflammation. In present study, we cloned two cDNA molecules of CXCR1 and CXCR4 in Miichthys miiuy (miiuy croaker). In these two genes, we found the most highly conserved DRY motif in the second intracellular loop adjacent to the third transmembrane domain. The expressions of CXCR1 and CXCR4 showed that they were ubiquitously expressed in ten normal tissues. After infection with Vibrio anguillarum and Vibrio harveyi, the expressions of CXCRs in the immune tissues were significantly regulated in most of tissues except that of CXCR1 in the kidney after V. harveyi injection. Evolutionary analysis showed that only the ancestral lineages of CXCR4 in amphibians underwent positive selection, indicating that the ancestors of amphibians boarded the land and had to further evolve to adapt to terrestrial environments. Multiple ML methods were implemented to detect the robust positively selected candidates for sites. In total, we detected 12 and 3 positively selected sites in the subsets of current mammal and fish CXCR1 genes, and only one site under positive selection was found in mammalian CXCR4 subsets. These positively selected sites were mainly located in the extracellular domains of CXCRs. The sliding window analysis and evolution test tended to favor positive selection acting on the N-terminal domain of CXCR1, which was the critical region for ligand/receptor signaling for neutrophils and receptor-ligand interaction, indicating that the N-terminal of CXCR1 in mammals underwent more positive selection than that of fish.

Eco-immunology of fish invasions: the role of MHC variation

The relationship between invaders and the pathogens encountered in their new environment can have a large effect on invasion success. Invaders can become free from their natural pathogens and reallocate costly immune resources to growth and reproduction, thereby increasing invasion success. Release from enemies and relaxation of selective pressures could render newly founded populations more variable at immune-related genes, such as the major histocompatibility complex (MHC), particularly when they have different origins. Using rainbow and brown trout, two of the world's most successful fish invaders, we tested the general hypothesis that invaders should display high intrapopulation immunogenetic diversity and interpopulation divergence, due to the interplay between genetic drift and successive waves of genetically divergent introductions. We analysed genetic diversity and signatures of selection at the MHC class II β immune-related locus. In both species, MHC diversity (allelic richness and heterozygosity) for southern hemisphere populations was similar to values reported for populations at their native range. However, MHC functional diversity was limited, and population immunogenetic structuring weaker than that observed using neutral markers. Depleted MHC functional diversity could reflect a decrease in immune response, immune-related assortative mating or selection for resistance to newly encountered parasites. Given that the role of MHC diversity in the survival of these populations remains unclear, depleted functional diversity of invasive salmonids could compromise their long-term persistence. A better understanding of the eco-immunology of invaders may help in managing and preventing the impact of biological invasions, a major cause of loss of biodiversity worldwide.

On the validity of evolutionary models with site-specific parameters

Evolutionary models that make use of site-specific parameters have recently been criticized on the grounds that parameter estimates obtained under such models can be unreliable and lack theoretical guarantees of convergence. We present a simulation study providing empirical evidence that a simple version of the models in question does exhibit sensible convergence behavior and that additional taxa, despite not being independent of each other, lead to improved parameter estimates. Although it would be desirable to have theoretical guarantees of this, we argue that such guarantees would not be sufficient to justify the use of these models in practice. Instead, we emphasize the importance of taking the variance of parameter estimates into account rather than blindly trusting point estimates – this is standardly done by using the models to construct statistical hypothesis tests, which are then validated empirically via simulation studies.

The evolution of Vp1 gene in enterovirus C species sub-group that contains types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99

Genus Enterovirus (Family Picornaviridae,) consists of twelve species divided into genetically diverse types by their capsid protein VP1 coding sequences. Each enterovirus type can further be divided into intra-typic sub-clusters (genotypes). The aim of this study was to elucidate what leads to the emergence of novel enterovirus clades (types and genotypes). An evolutionary analysis was conducted for a sub-group of Enterovirus C species that contains types Coxsackievirus A21 (CVA-21), CVA-24, Enterovirus C95 (EV-C95), EV-C96 and EV-C99. VP1 gene datasets were collected and analysed to infer the phylogeny, rate of evolution, nucleotide and amino acid substitution patterns and signs of selection. In VP1 coding gene, high intra-typic sequence diversities and robust grouping into distinct genotypes within each type were detected. Within each type the majority of nucleotide substitutions were synonymous and the non-synonymous substitutions tended to cluster in distinct highly polymorphic sites. Signs of positive selection were detected in some of these highly polymorphic sites, while strong negative selection was indicated in most of the codons. Despite robust clustering to intra-typic genotypes, only few genotype-specific ‘signature’ amino acids were detected. In contrast, when different enterovirus types were compared, there was a clear tendency towards fixation of type-specific ‘signature’ amino acids. The results suggest that permanent fixation of type-specific amino acids is a hallmark associated with evolution of different enterovirus types, whereas neutral evolution and/or (frequency-dependent) positive selection in few highly polymorphic amino acid sites are the dominant forms of evolution when strains within an enterovirus type are compared.

Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic

HLA-restricted immune escape mutations that persist following HIV transmission could gradually spread through the viral population, thereby compromising host antiviral immunity as the epidemic progresses. To assess the extent and phenotypic impact of this phenomenon in an immunogenetically diverse population, we genotypically and functionally compared linked HLA and HIV (Gag/Nef) sequences from 358 historic (1979-1989) and 382 modern (2000-2011) specimens from four key cities in the North American epidemic (New York, Boston, San Francisco, Vancouver). Inferred HIV phylogenies were star-like, with approximately two-fold greater mean pairwise distances in modern versus historic sequences. The reconstructed epidemic ancestral (founder) HIV sequence was essentially identical to the North American subtype B consensus. Consistent with gradual diversification of a "consensus-like" founder virus, the median "background" frequencies of individual HLA-associated polymorphisms in HIV (in individuals lacking the restricting HLA[s]) were ∼ 2-fold higher in modern versus historic HIV sequences, though these remained notably low overall (e.g. in Gag, medians were 3.7% in the 2000s versus 2.0% in the 1980s). HIV polymorphisms exhibiting the greatest relative spread were those restricted by protective HLAs. Despite these increases, when HIV sequences were analyzed as a whole, their total average burden of polymorphisms that were "pre-adapted" to the average host HLA profile was only ∼ 2% greater in modern versus historic eras. Furthermore, HLA-associated polymorphisms identified in historic HIV sequences were consistent with those detectable today, with none identified that could explain the few HIV codons where the inferred epidemic ancestor differed from the modern consensus. Results are therefore consistent with slow HIV adaptation to HLA, but at a rate unlikely to yield imminent negative implications for cellular immunity, at least in North America. Intriguingly, temporal changes in protein activity of patient-derived Nef (though not Gag) sequences were observed, suggesting functional implications of population-level HIV evolution on certain viral proteins.

Evolution of separate predation- and defence-evoked venoms in carnivorous cone snails

Venomous animals are thought to inject the same combination of toxins for both predation and defence, presumably exploiting conserved target pharmacology across prey and predators. Remarkably, cone snails can rapidly switch between distinct venoms in response to predatory or defensive stimuli. Here, we show that the defence-evoked venom of Conus geographus contains high levels of paralytic toxins that potently block neuromuscular receptors, consistent with its lethal effects on humans. In contrast, C. geographus predation-evoked venom contains prey-specific toxins mostly inactive at human targets. Predation- and defence-evoked venoms originate from the distal and proximal regions of the venom duct, respectively, explaining how different stimuli can generate two distinct venoms. A specialized defensive envenomation strategy is widely evolved across worm, mollusk and fish-hunting cone snails. We propose that defensive toxins, originally evolved in ancestral worm-hunting cone snails to protect against cephalopod and fish predation, have been repurposed in predatory venoms to facilitate diversification to fish and mollusk diets.

Pigeon circoviruses display patterns of recombination, genomic secondary structure and selection similar to those of beak and feather disease viruses

Pigeon circovirus (PiCV) has a ~2 kb genome circular ssDNA genome. All but one of the known PiCV isolates have been found infecting pigeons in various parts of the world. In this study, we screened 324 swab and tissue samples from Polish pigeons and recovered 30 complete genomes, 16 of which came from birds displaying no obvious pathology. Together with 17 other publicly available PiCV complete genomes sampled throughout the Northern Hemisphere and Australia, we find that PiCV displays a similar degree of genetic diversity to that of the related psittacine-infecting circovirus species, beak and feather disease virus (BFDV). We show that, as is the case with its pathology and epidemiology, PiCV also displays patterns of recombination, genomic secondary structure and natural selection that are generally very similar to those of BFDV. It is likely that breeding facilities play a significant role in the emergence of new recombinant PiCV variants and given that ~50 % of the domestic pigeon population is infected subclinically, all pigeon breeding stocks should be screened routinely for this virus.

Changes in diversification patterns and signatures of selection during the evolution of murinae-associated hantaviruses

In the last 50 years, hantaviruses have significantly affected public health worldwide, but the exact extent of the distribution of hantavirus diseases, species and lineages and the risk of their emergence into new geographic areas are still poorly known. In particular, the determinants of molecular evolution of hantaviruses circulating in different geographical areas or different host species are poorly documented. Yet, this understanding is essential for the establishment of more accurate scenarios of hantavirus emergence under different climatic and environmental constraints. In this study, we focused on Murinae-associated hantaviruses (mainly Seoul Dobrava and Hantaan virus) using sequences available in GenBank and conducted several complementary phylogenetic inferences. We sought for signatures of selection and changes in patterns and rates of diversification in order to characterize hantaviruses’ molecular evolution at different geographical scales (global and local). We then investigated whether these events were localized in particular geographic areas. Our phylogenetic analyses supported the assumption that RNA virus molecular variations were under strong evolutionary constraints and revealed changes in patterns of diversification during the evolutionary history of hantaviruses. These analyses provide new knowledge on the molecular evolution of hantaviruses at different scales of time and space.

Diversification of a single ancestral gene into a successful toxin superfamily in highly venomous Australian funnel-web spiders

BACKGROUND

Spiders have evolved pharmacologically complex venoms that serve to rapidly subdue prey and deter predators. The major toxic factors in most spider venoms are small, disulfide-rich peptides. While there is abundant evidence that snake venoms evolved by recruitment of genes encoding normal body proteins followed by extensive gene duplication accompanied by explosive structural and functional diversification, the evolutionary trajectory of spider-venom peptides is less clear.

RESULTS

Here we present evidence of a spider-toxin superfamily encoding a high degree of sequence and functional diversity that has evolved via accelerated duplication and diversification of a single ancestral gene. The peptides within this toxin superfamily are translated as prepropeptides that are posttranslationally processed to yield the mature toxin. The N-terminal signal sequence, as well as the protease recognition site at the junction of the propeptide and mature toxin are conserved, whereas the remainder of the propeptide and mature toxin sequences are variable. All toxin transcripts within this superfamily exhibit a striking cysteine codon bias. We show that different pharmacological classes of toxins within this peptide superfamily evolved under different evolutionary selection pressures.

CONCLUSIONS

Overall, this study reinforces the hypothesis that spiders use a combinatorial peptide library strategy to evolve a complex cocktail of peptide toxins that target neuronal receptors and ion channels in prey and predators. We show that the ω-hexatoxins that target insect voltage-gated calcium channels evolved under the influence of positive Darwinian selection in an episodic fashion, whereas the κ-hexatoxins that target insect calcium-activated potassium channels appear to be under negative selection. A majority of the diversifying sites in the ω-hexatoxins are concentrated on the molecular surface of the toxins, thereby facilitating neofunctionalisation leading to new toxin pharmacology.

Extraordinary MHC class II B diversity in a non-passerine, wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae)

The major histocompatibility complex (MHC) hosts the most polymorphic genes ever described in vertebrates. The MHC triggers the adaptive branch of the immune response, and its extraordinary variability is considered an evolutionary consequence of pathogen pressure. The last few years have witnessed the characterization of the MHC multigene family in a large diversity of bird species, unraveling important differences in its polymorphism, complexity, and evolution. Here, we characterize the first MHC class II B sequences isolated from a Rallidae species, the Eurasian Coot Fulica atra. A next-generation sequencing approach revealed up to 265 alleles that translated into 251 different amino acid sequences (β chain, exon 2) in 902 individuals. Bayesian inference identified up to 19 codons within the presumptive peptide-binding region showing pervasive evidence of positive, diversifying selection. Our analyses also detected a significant excess of high-frequency segregating sites (average Tajima's D = 2.36, P < 0.05), indicative of balancing selection. We found one to six different alleles per individual, consistent with the occurrence of at least three MHC class II B gene duplicates. However, the genotypes comprised of three alleles were by far the most abundant in the population investigated (49.4%), followed by those with two (29.6%) and four (17.5%) alleles. We suggest that these proportions are in agreement with the segregation of MHC haplotypes differing in gene copy number. The most widespread segregating haplotypes, according to our findings, would contain one single gene or two genes. The MHC class II of the Eurasian Coot is a valuable system to investigate the evolutionary implications of gene copy variation and extensive variability, the greatest ever found, to the best of our knowledge, in a wild population of a non-passerine bird.