Nucleotide Substitution Rates Research Articles

Exploring the Limits and Causes of Plastid Genome Expansion in Volvocine Green Algae.

Plastid genomes are not normally celebrated for being large. But researchers are steadily uncovering algal lineages with big and, in rare cases, enormous plastid DNAs (ptDNAs), such as volvocine green algae. Plastome sequencing of five different volvocine species has revealed some of the largest, most repeat-dense plastomes on record, including that of Volvox carteri (∼525 kb). Volvocine algae have also been used as models for testing leading hypotheses on organelle genome evolution (e.g., the mutational hazard hypothesis), and it has been suggested that ptDNA inflation within this group might be a consequence of low mutation rates and/or the transition from a unicellular to multicellular existence. Here, we further our understanding of plastome size variation in the volvocine line by examining the ptDNA sequences of the colonial species Yamagishiella unicocca and Eudorina sp. NIES-3984 and the multicellular Volvox africanus, which are phylogenetically situated between species with known ptDNA sizes. Although V. africanus is closely related and similar in multicellular organization to V. carteri, its ptDNA was much less inflated than that of V. carteri. Synonymous- and noncoding-site nucleotide substitution rate analyses of these two Volvox ptDNAs suggest that there are drastically different plastid mutation rates operating in the coding versus intergenic regions, supporting the idea that error-prone DNA repair in repeat-rich intergenic spacers is contributing to genome expansion. Our results reinforce the idea that the volvocine line harbors extremes in plastome size but ultimately shed doubt on some of the previously proposed hypotheses for ptDNA inflation within the lineage.

Purifying Selective Pressure Suggests the Functionality of a Vitamin B12 Biosynthesis Pathway in a Global Population of Mycobacterium tuberculosis.

Mycobacterium tuberculosis is one of the deadliest and most challenging pathogens to study in current microbiological research. One of the issues that remains to be resolved is the importance of cobalamin in the metabolism of M. tuberculosis. The functionality of a vitamin B12 biosynthesis pathway in M. tuberculosis is under dispute, and the ability of this pathogen to scavenge vitamin B12 from the host is unknown. Here, we quantified the ratios of nonsynonymous and synonymous nucleotide substitution rates (dN/dS) in the genes involved in vitamin B12 biosynthesis and transport and in genes encoding cobalamin-dependent enzymes in nearly four thousand strains of M. tuberculosis. We showed that purifying selection is the dominant force acting on cobalamin-related genes at the levels of individual codons, genes and groups of genes. We conclude that cobalamin-related genes may not be essential but are adaptive for M. tuberculosis in clinical settings. Furthermore, the cobalamin biosynthesis pathway is likely to be functional in this species.

Possible risks posed by single-stranded DNA viruses of pigs associated with xenotransplantation.

Routine large‐scale xenotransplantation from pigs to humans is getting closer to clinical reality owing to several state‐of‐the‐art technologies, especially the ability to rapidly engineer genetically defined pigs. However, using pig organs in humans poses risks including unwanted cross‐species transfer of viruses and adaption of these pig viruses to the human organ recipient. Recent developments in the field of virology, including the advent of metagenomic techniques to characterize entire viromes, have led to the identification of a plethora of viruses in many niches. Single‐stranded DNA (ssDNA) viruses are the largest group prevalent in virome studies in mammals. Specifically, the ssDNA viral genomes are characterized by a high rate of nucleotide substitution, which confers a proclivity to adapt to new hosts and cross‐species barriers. Pig‐associated ssDNA viruses include torque teno sus viruses (TTSuV) in the Anelloviridae family, porcine parvoviruses (PPV), and porcine bocaviruses (PBoV) both in the family of Parvoviridae, and porcine circoviruses (PCV) in the Circoviridae family, some of which have been confirmed to be pathogenic to pigs. The risks of these viruses for the human recipient during xenotransplantation procedures are relatively unknown. Based on the scant knowledge available on the prevalence, predilection, and pathogenicity of pig‐associated ssDNA viruses, careful screening and monitoring are required. In the case of positive identification, risk assessments and strategies to eliminate these viruses in xenotransplantation pig stock may be needed.

The Complete Mitochondrial Genome of Glyptothorax macromaculatus Provides a Well-Resolved Molecular Phylogeny of the Chinese Sisorid Catfishes.

Previous phylogenetic analyses of the Chinese sisorid catfishes have either been poorly resolved or have not included all the 12 sisorid genera. Here, we successfully assembled the first complete mitochondrial genome of the sisorid fish Glyptothorax macromaculatus. Based on this novel mitochondrial genome and previously published mitochondrial genomes in the Sisoridae, we generated maximum likelihood and Bayesian phylogenies. We dated our preferred topology using fossil calibration points. We also tested the protein-coding genes in the mitochondrial genomes of the glyptosternoid fishes for signals of natural selection by comparing the nucleotide substitution rate along the branch ancestral to the glyptosternoid fishes to other branches in our topology. The mitochondrial sequence structure of G. macromaculatus was similar to those known from other vertebrates, with some slight differences. Our sisorid phylogenies were well-resolved and well-supported, with exact congruence between the different phylogenetic methods. This robust phylogeny clarified the relationships among the Chinese sisorid genera and strongly supported the division of the family into three main clades. Interestingly, the glyptosternoid divergence time predicted by our molecular dating analysis coincided with the uplift of the Tibetan Plateau, suggesting that geology may have influenced speciation in the Sisoridae. Among the mitochondrial protein-coding genes, atp8 may have most rapidly evolved, and atp6 may have been subjected to positive selection pressure to adapt to high elevations. In summary, this study provided novel insights into the phylogeny, evolution and high-altitude adaptions of the Chinese sisorid fishes.

Divergence time, historical biogeography and evolutionary rate estimation of the order Bangiales (Rhodophyta) inferred from multilocus data

Bangiales is the only order of the Bangiophyceae and has been suggested to be monophyletic. This order contains approximately 190 species and is distributed worldwide. Previous molecular studies have produced robust phylogenies among the red algae, but the divergence times, historical biogeography and evolutionary rates of Bangiales have rarely been studied. Phylogenetic relationships within the Bangiales were examined using the concatenated gene sets from all available organellar genomes. This analysis has revealed the topology (((Bangia, Porphyra) Pyropia) Wildemania). Molecular dating indicates that Bangiales diversified approximately 246.40 million years ago (95% highest posterior density (HPD)=194.78–318.24 Ma, posterior probability (PP)=0.99) in the Late Permian and Early Triassic, and that the ancestral species most likely originated from eastern Gondwanaland (currently New Zealand and Australia) and subsequently began to spread and evolve worldwide. Based on pairwise comparisons, we found a slower rate of nucleotide substitutions and lower rates of diversification in Bangiales relative to Florideophyceae. Compared with Viridiplantae (green algae and land plants), the evolutionary rates of Bangiales and other Rhodophyte groups were found to be dramatically faster, by more than 3-fold for plastid genome (ptDNA) and 15-fold for mitochondrial genome (mtDNA). In addition, an average 2.5-fold lower dN/dS was found for the algae than for the land plants, which indicates purifying selection of the algae.

BA9 lineage of respiratory syncytial virus from across the globe and its evolutionary dynamics.

Respiratory syncytial virus (RSV) is an important pathogen of global significance. The BA9 is one of the most predominant lineages of the BA genotype of group B RSV that has acquired a 60bp duplication in its G protein gene. We describe the local and global evolutionary dynamics of the second hyper variable region in the C- terminal of the G protein gene of the BA9 lineage. A total of 418 sequences (including 31 study and 387 GenBank strains) from 29 different countries were used for phylogenetic analysis. This analysis showed that the study strains clustered with BA (BA9 and BA8) and SAB4 genotype of group B RSV. We performed time-scaled evolutionary clock analyses using Bayesian Markov chain Monte Carlo methods. We also carried out glycosylation, selection pressure, mutational, entropy and Network analyses of the BA9 lineage. The time to the most recent common ancestor (tMRCA) of the BA genotype and BA9 lineage were estimated to be the years 1995 (95% HPD; 1987–1997) and 2000 (95% HPD; 1998–2001), respectively. The nucleotide substitution rate of the BA genotype [(4.58×10−3 (95% HPD; 3.89–5.29×10−3) substitution/site/year] was slightly faster than the BA9 lineage [4.03×10−3 (95% HPD; 4.65–5.2492×10−3)]. The BA9 lineage was categorized into 3 sub lineages (I, II and III) based on the Bayesian and Network analyses. The local transmission pattern suggested that BA9 is the predominant lineage of BA viruses that has been circulating in India since 2002 though showing fluctuations in its effective population size. The BA9 lineage established its global distribution with report from 23 different countries over the past 16 years. The present study augments our understanding of RSV infection, its epidemiological dynamics warranting steps towards its overall global surveillance.

Plastome characteristics of Cannabaceae

Cannabaceae is an economically important family that includes ten genera and ca. 117 accepted species. To explore the structure and size variation of their plastomes, we sequenced ten plastomes representing all ten genera of Cannabaceae. Each plastome possessed the typical angiosperm quadripartite structure and contained a total of 128 genes. The Inverted Repeat (IR) regions in five plastomes had experienced small expansions (330–983 bp) into the Large Single-Copy (LSC) region. The plastome of Chaetachme aristata has experienced a 942-bp IR contraction and lost rpl22 and rps19 in its IRs. The substitution rates of rps19 and rpl22 decreased after they shifted from the LSC to IR. A 270-bp inversion was detected in the Parasponia rugosa plastome, which might have been mediated by 18-bp inverted repeats. Repeat sequences, simple sequence repeats, and nucleotide substitution rates varied among these plastomes. Molecular markers with more than 13% variable sites and 5% parsimony-informative sites were identified, which may be useful for further phylogenetic analysis and species identification. Our results show strong support for a sister relationship between Gironniera and Lozanell (BS = 100). Celtis, Cannabis-Humulus, Chaetachme-Pteroceltis, and Trema-Parasponia formed a strongly supported clade, and their relationships were well resolved with strong support (BS = 100). The availability of these ten plastomes provides valuable genetic information for accurately identifying species, clarifying taxonomy and reconstructing the intergeneric phylogeny of Cannabaceae.

Madagascar ground gecko genome analysis characterizes asymmetric fates of duplicated genes

BackgroundConventionally, comparison among amniotes – birds, mammals, and reptiles – has often been approached through analyses of mammals and, for comparison, birds. However, birds are morphologically and physiologically derived and, moreover, some parts of their genomes are recognized as difficult to sequence and/or assemble and are thus missing in genome assemblies. Therefore, sequencing the genomes of reptiles would aid comparative studies on amniotes by providing more comprehensive coverage to help understand the molecular mechanisms underpinning evolutionary changes.ResultsHerein, we present the whole genome sequences of the Madagascar ground gecko (Paroedura picta), a promising study system especially in developmental biology, and used it to identify changes in gene repertoire across amniotes. The genome-wide analysis of the Madagascar ground gecko allowed us to reconstruct a comprehensive set of gene phylogenies comprising 13,043 ortholog groups from diverse amniotes. Our study revealed 469 genes retained by some reptiles but absent from available genome-wide sequence data of both mammals and birds. Importantly, these genes, herein collectively designated as ‘elusive’ genes, exhibited high nucleotide substitution rates and uneven intra-genomic distribution. Furthermore, the genomic regions flanking these elusive genes exhibited distinct characteristics that tended to be associated with increased gene density, repeat element density, and GC content.ConclusionThis highly continuous and nearly complete genome assembly of the Madagascar ground gecko will facilitate the use of this species as an experimental animal in diverse fields of biology. Gene repertoire comparisons across amniotes further demonstrated that the fate of a duplicated gene can be affected by the intrinsic properties of its genomic location, which can persist for hundreds of millions of years.

Positive Selection Driving Cytoplasmic Genome Evolution of the Medicinally Important Ginseng Plant Genus Panax.

Panax L. (the ginseng genus) is a shade-demanding group within the family Araliaceae and all of its species are of crucial significance in traditional Chinese medicine. Phylogenetic and biogeographic analyses demonstrated that two rounds of whole genome duplications accompanying with geographic and ecological isolations promoted the diversification of Panax species. However, contributions of the cytoplasmic genomes to the adaptive evolution of Panax species remained largely uninvestigated. In this study, we sequenced the chloroplast and mitochondrial genomes of 11 accessions belonging to seven Panax species. Our results show that heterogeneity in nucleotide substitution rate is abundant in both of the two cytoplasmic genomes, with the mitochondrial genome possessing more variants at the total level but the chloroplast showing higher sequence polymorphisms at the genic regions. Genome-wide scanning of positive selection identified five and 12 genes from the chloroplast and mitochondrial genomes, respectively. Functional analyses further revealed that these selected genes play important roles in plant development, cellular metabolism and adaptation. We therefore conclude that positive selection might be one of the potential evolutionary forces that shaped nucleotide variation pattern of these Panax species. In particular, the mitochondrial genes evolved under stronger selective pressure compared to the chloroplast genes.

Analysis of evolutionary rate of HIV-1 subtype B using blood donor samples in Japan.

There are few reports on HIV-1 intra-host evolutionary rate in asymptomatic treatment-naïve patients. Here, the HIV-1 intra-host evolutionary rate was estimated based on HIV-1 RNA sequences from plasma samples of blood donors in Japan. Blood donors were assumed to have received no treatment for and have no symptoms of HIV-1 infection because they were healthy, and declared no risky behaviors of HIV-1 infection on a self-reported questionnaire or interview followed by donation. HIV-1 RNA was obtained from 85 plasma samples from 36 blood donors who donated blood multiple times and were HIV-1-positive. The C2V3C3 region which encodes for a part of the envelope protein, and the V3 loop in the C2V3C3 region were analyzed by RT-PCR and direct sequencing, and the sequences were compared. The nucleotide substitution rate was calculated by linear regression. All HIV-1 samples analyzed were classified as subtype B. The mean nucleotide substitution rate in C2V3C3 was calculated to be 6.2 × 10-3-1.8 × 10-2/site/year (V3: 4.5 × 10-3-2.3 × 10-2/site/year). The mean non-synonymous substitution rate in C2V3C3 was calculated to be 5.2 × 10-3-1.7 × 10-2/site/year (V3: 4.5 × 10-3-2.1 × 10-2/site/year). The mean synonymous substitution rate in C2V3C3 was calculated to be 1.1 × 10-4-2.3 × 10-3/site/year (V3: 2.9 × 10-3/site/year). Among HIV-1 subtype B RNA-positive blood donors in Japan, the nucleotide substitution rate in C2V3C3 was estimated to be higher than that of reported cases using HIV-1 samples mainly obtained from AIDS patients. Compared to AIDS patients, immune responses against HIV-1 are probably more effective in HIV-1 RNA-positive blood donors. Consequently, immune pressure presumably promotes mutation of the virus genome.

Molecular Characterization of Mx1 Gene in Native Indian Breeds of Chicken

ABSTRACTThe genetic polymorphism of Mx1 gene was explored in Indian chicken breeds. PCR-RFLP analysis in 102 bp fragment of partial intron 13 and partial exon 14 of Mx1 gene revealed two genotypes viz. RS and SS with two alleles viz. R and S both in Naked Neck and Tellicherry breeds of chicken. The homozygous genotype RR was not identified. When deduced amino acid sequences were compared, the asparagine amino acid was found to be substituted in “R” allele for serine in “S” allele. PCR-SSCP analysis of 284 bp fragment in 5′-UTR and partial promoter region revealed three genotypes viz. CC, CG, and CH with three different alleles viz. C, G, and H in Naked Neck breed of chicken and five genotypes viz. DI, JK, KK, KL, and KM with six different alleles viz. D, I, J, K, L, and M in Tellicherry breed of chicken. The homozygous genotypes viz. GG and HH in Naked Neck and DD, II, JJ, LL, and MM in Tellicherry chicken was not identified. The nucleotide substitution rate estimated to be in the range of 0.004–0.011. The identified genetic variation can be helpful for better insight to disease resistance property of the Mx1 gene.

PhyloMAd: efficient assessment of phylogenomic model adequacy.

Statistical phylogenetic inference plays an important role in evolutionary biology. The accuracy of phylogenetic methods relies on having suitable models of the evolutionary process. Various tools allow comparisons of candidate phylogenetic models, but assessing the absolute performance of models remains a considerable challenge. We introduce PhyloMAd, a user-friendly application for assessing the adequacy of commonly used models of nucleotide substitution and among-lineage rate variation. Our software implements a fast, likelihood-based method of model assessment that is tractable for analyses of large multi-locus datasets. PhyloMAd provides a means of informing model improvement, or selecting data to enhance the evolutionary signal in phylogenomic analyses. PhyloMAd, together with a manual, a tutorial and the source code, are freely available from the GitHub repository github.com/duchene/phylomad.

Molecular evolution and phylodynamics of hepatitis B virus infection circulating in Iran.

Previous local and national Iranian publications indicate that all Iranian hepatitis B virus (HBV) strains belong to HBV genotype D. The aim of this study was to analyze the evolutionary history of HBV infection in Iran for the first time, based on an intensive phylodynamic study. The evolutionary parameters, time to most recent common ancestor (tMRCA), and the population dynamics of infections were investigated using the Bayesian Monte Carlo Markov chain (BMCMC). The effective sample size (ESS) and sampling convergence were then monitored. After sampling from the posterior distribution of the nucleotide substitution rate and other evolutionary parameters, the point estimations (median) of these parameters were obtained. All Iranian HBV isolates were of genotype D, sub-type ayw2. The origin of HBV is regarded as having evolved first on the eastern border, before moving westward, where Isfahan province then hosted the virus. Afterwards, the virus moved to the south and west of the country. The tMRCA of HBV in Iran was estimated to be around 1894, with a 95% credible interval between the years 1701 and 1957. The effective number of infections increased exponentially from around 1925 to 1960. Conversely, from around 1992 onwards, the effective number of HBV infections has decreased at a very high rate. Phylodynamic inference clearly demonstrates a unique homogenous pattern of HBV genotype D compatible with a steady configuration of the decreased effective number of infections in the population in recent years, possibly due to the implementation of blood donation screening and vaccination programs. Adequate molecular epidemiology databases for HBV are crucial for infection prevention and treatment programs.

Salt and drought stress and ABA responses related to bZIP genes from V. radiata and V. angularis

Mung bean and adzuki bean are warm-season legumes widely cultivated in China. However, bean production in major producing regions is limited by biotic and abiotic stress, such as drought and salt stress. Basic leucine zipper (bZIP) genes play key roles in responses to various biotic and abiotic stresses. However, only several bZIP genes involved in drought and salt stress in legumes, especially Vigna radiata and Vigna angularis, have been identified. In this study, we identified 54 and 50 bZIP proteins from whole-genome sequences of V. radiata and V. angularis, respectively. First, we comprehensively surveyed the characteristics of all bZIP genes, including their gene structure, chromosome distribution and motif composition. Phylogenetic trees showed that VrbZIP and VabZIP proteins were divided into ten clades comprising nine known and one unknown subgroup. The results of the nucleotide substitution rate of the orthologous gene pairs showed that bZIP proteins have undergone strong purifying selection: V. radiata and V. angularis diverged 1.25 million years ago (mya) to 9.20 mya (average of 4.95 mya). We also found that many cis-acting regulatory elements (CAREs) involved in abiotic stress and plant hormone responses were detected in the putative promoter regions of the bZIP genes. Finally, using the quantitative real-time PCR (qRT-PCR) method, we performed expression profiling of the bZIP genes in response to drought, salt and abscisic acid (ABA). We identified several bZIP genes that may be involved in drought and salt responses. Generally, our results provided useful and rich resources of VrbZIP and VabZIP genes for the functional characterization and understanding of bZIP transcription factors (TFs) in warm-season legumes. In addition, our results revealed important and interesting data – a subset of VrbZIP and VabZIP gene expression profiles in response to drought, salt and ABA stress. These results provide gene expression evidence for the selection of candidate genes under drought and salt stress for future study.

Evolutionary dynamics of GII.17 norovirus.

BackgroundDuring the winter of 2014–2015, a rarely reported norovirus (NoV) genotype GII.17 was found to have increased its frequency in norovirus outbreaks in East Asia, surpassing the GII.4 NoV infections. GII.17 genotype has been detected for over three decades in the world. The aim of this study is to examine the evolutionary dynamics of GII.17 over the last four decades.MethodsNoV GII.17 sequences with complete or nearly complete VP1 were downloaded from GenBank and the phylogenetic analyses were then conducted.ResultsThe maximum likelihood analysis showed that GII.17 genotype could be divided into four different clades (Clades A–D). The strains detected after 2012, which could be the cause of the outbreaks, were separated into Clades C–D with their mean amino acid distance being 4.5%. Bayesian Markov chain Monte Carlo analyses indicated that the rate of nucleotide substitution per sites was 1.68 × 10−3 nucleotide substitutions/site/year and the time of the most recent common ancestor was 1840. The P2 subdomain of GII.17 was highly variable with 44% (56/128) amino acids variations including two insertions at positions 295–296 and one deletion at position 385 (Clades C and D) and one insertion at position 375 (Clade D). Variations existed in Epitopes A, B and D corresponding to GII.4 and human histo-blood group antigens binding site I in P2 subdomain.ConclusionThe novel GII.17 strains that caused outbreaks in 2013–2015 may have two new variants. The evolvement of HBGAs binding site and epitopes in P2 subdomain might contribute to the novel GII.17 strains predominance in some regions.

Dating the Cryptococcus gattii Dispersal to the North American Pacific Northwest.

The emergence of Cryptococcus gattii, previously regarded as a predominantly tropical pathogen, in the temperate climate of the North American Pacific Northwest (PNW) in 1999 prompted several questions. The most prevalent among these was the timing of the introduction of this pathogen to this novel environment. Here, we infer tip-dated timing estimates for the three clonal C.gattii populations observed in the PNW, VGIIa, VGIIb, and VGIIc, based on whole-genome sequencing of 134 C.gattii isolates and using Bayesian evolutionary analysis by sampling trees (BEAST). We estimated the nucleotide substitution rate for each lineage (1.59 × 10-8, 1.59 × 10-8, and 2.70 × 10-8, respectively) to be an order of magnitude higher than common neutral fungal mutation rates (2.0 × 10-9), indicating a microevolutionary rate (e.g., successive clonal generations in a laboratory) in comparison to a species' slower, macroevolutionary rate (e.g., when using fossil records). The clonal nature of the PNW C.gattii emergence over a narrow number of years would therefore possibly explain our higher mutation rates. Our results suggest that the mean time to most recent common ancestor for all three sublineages occurred within the last 60 to 100years. While the cause of C.gattii dispersal to the PNW is still unclear, our research estimates that the arrival is neither ancient nor very recent (i.e., <25 years ago), making a strong case for an anthropogenic introduction. IMPORTANCE The recent emergence of the pathogenic fungus Cryptococcus gattii in the Pacific Northwest (PNW) resulted in numerous investigations into the epidemiological and enzootic impacts, as well as multiple genomic explorations of the three primary molecular subtypes of the fungus that were discovered. These studies lead to the general conclusion that the subtypes identified likely emerged out of Brazil. Here, we conducted genomic dating analyses to determine the ages of the various lineages seen in the PNW and propose hypothetical causes for the dispersal events. Bayesian evolutionary analysis strongly suggests that these independent fungal populations in the PNW are all 60 to 100 years old, providing a timing that is subsequent to the opening of the Panama Canal, which allowed for more direct shipping between Brazil and the western North American coastline, a possible driving event for these fungal translocation events.

Complexity and recombination analysis of novel begomovirus associated with Spinach yellow vein disease in India

Leaves with vein yellowing were observed on Spinach (Spinacia oleracea) in Rajasthan province of India. The plants were severely stunted in growth and exhibited symptoms typical of Begomovirus infection. Here we report a comprehensive study of complexity and recombination of novel begomovirus associated with spinach yellow vein disease. It is probably a perpetuation of Papaya leaf curl virus in absence of its main host and thus raises concern about its spread to other crops. Analysis of samples collected in the survey indicates that SYVD infected plants are associated with a begomoviruses (novel species showed <91% nucleotide sequence similarity) and its associated betasatellites. Interestingly, most of the begomoviruses were found to be inter-species recombinants and the betasatellites possess high nucleotide variability, due to recombination. This suggested the existence of a complex recombination pattern to form “gene pool” of the crop-infecting begomovirus–betasatellites complexes. The Nucleotide diversity (π) and nucleotide substitution rates were determined for the DNA-A of begomoviruses SYVSV is 0.10382 and 2.436×10−3substitutions site−1year−1 respectively. Thus mutation and recombination are driving forces for the emergence and evolution of new begomoviruses. This strengthens the hypothesis for increase of host range of begomovirus and indicates that the high evolutionary rate of the begomoviruses is primarily due to frequent recombination and mutation.

Longitudinal Surveillance of Porcine Rotavirus B Strains from the United States and Canada and In Silico Identification of Antigenically Important Sites.

Rotavirus B (RVB) is an important swine pathogen, but control and prevention strategies are limited without an available vaccine. To develop a subunit RVB vaccine with maximal effect, we characterized the amino acid sequence variability and predicted antigenicity of RVB viral protein 7 (VP7), a major neutralizing antibody target, from clinically infected pigs in the United States and Canada. We identified genotype-specific antigenic sites that may be antibody neutralization targets. While some antigenic sites had high amino acid functional group diversity, nine antigenic sites were completely conserved. Analysis of nucleotide substitution rates at amino acid sites (dN/dS) suggested that negative selection appeared to be playing a larger role in the evolution of the identified antigenic sites when compared to positive selection, and was identified in six of the nine conserved antigenic sites. These results identified important characteristics of RVB VP7 variability and evolution and suggest antigenic residues on RVB VP7 that are negatively selected and highly conserved may be good candidate regions to include in a subunit vaccine design due to their tendency to remain stable.

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.

Global genetic variation and transmission dynamics of H9N2 avian influenza virus.

The H9N2 influenza viruses are extensively circulating in the poultry population, and variable genotypes can be generated through mutation, recombination and reassortment, which may be better adapted to infect a new host, resist drug treatment or escape immune pressure. The LPAI H9N2 viruses have the potential to evolve towards high levels of virulence in human. Some studies about the regional dispersal were reported, but global dissemination and the drivers of the virus are poorly understood, particularly at the genome scale. Here, we have analysed all eight gene segments of 168 H9N2 genomes sampled randomly aiming to provide a panoramic framework for better understanding the genesis and genetic variation of the viruses, and utilized phylogeography and spatial epidemiology approaches to uncover the effects of the genetic variation, predictors and spread of H9N2 viruses. We found that more frequent reassortment events involve segments PA, NP and NS, and 21 isolates have possible mosaic structure resulting from recombination events. Estimates of gene-specific global dN/dS ratios showed that all genes were subject to purifying selection. However, a total of 13 sites were detected under positive selection by at least two of three methods, which located within segments HA, NA, M2, NS1 and PA. Additionally, we inferred that NA segment has the highest rate of nucleotide substitution, and its tMRCA estimate is the youngest than the remaining segments' inference. About the spatial history, air transportation of human was identified as the predominant driver of global viral migration using GLM analysis, and economic factors and geographical distance were the modest predictors. Higher migration rates were estimated between five pairs of regions (>0.01) indicating the frequent migration of the viruses between discrete geographical locations. Further, our Markov jumps analysis showed that viral migration is more frequent between Southern China and Northern China, and high rate of gene flow was observed between America and East Asia. Moreover, the America together with Southeast Asia acted as the primary hubs of global transmission, forming the trunk of evolutionary tree. These findings suggested a complex interaction between virus evolution, epidemiology and human behaviour.