Evidence of recombination among early-vaccination era measles virus strains

The first report and biological characterization of Avian Orthoavulavirus 16 in wild migratory waterfowl and domestic poultry in China reveal a potential threat to birds

RESEARCH HIGHLIGHTS

First confirmation of AOAV-16 in domestic and wild birds in China.AOAV-16 are low virulent viruses for chickens.Co-circulation/co-infection of AOAV-16 and H9N2 subtype AIV enhanced pathogenicity.Different intergenic sequences and recombination events exist within AOAV-16.

Evolutionary Analysis of a Parrot Bornavirus 2 Detected in a Sulphur-Crested Cockatoo (Cacatua galerita) Suggests a South American Ancestor

Parrot bornavirus (PaBV) is an RNA virus that causes Proventricular Dilatation Disease (PDD), neurological disorders, and death in Psittaciformes. Its diversity in South America is poorly known. We examined a Cacatua galerita presenting neuropathies, PDD, and oculopathies as the main signs. We detected PaBV through reverse transcription polymerase chain reaction (RT-PCR) and partial sequencing of the nucleoprotein (N) and matrix (M) genes. Maximum likelihood and Bayesian phylogenetic inferences classified it as PaBV-2. The nucleotide identity of the sequenced strain ranged from 88.3% to 90.3% against genotype PaBV-2 and from 80.2% to 84.4% against other genotypes. Selective pressure analysis detected signs of episodic diversifying selection in both the N and M genes. No recombination events were detected. Phylodynamic analysis estimated the time to the most recent common ancestor (TMRCA) as the year 1758 for genotype PaBV-2 and the year 1049 for the Orthobornavirus alphapsittaciforme species. Substitution rates were estimated at 2.73 × 10-4 and 4.08 × 10-4 substitutions per year per site for N and M, respectively. The analysis of population dynamics showed a progressive decline in the effective population size during the last century. Timescale phylogeographic analysis revealed a potential South American ancestor as the origin of genotypes 1, 2, and 8. These results contribute to our knowledge of the evolutionary origin, diversity, and dynamics of PaBVs in South America and the world. Additionally, it highlights the importance of further studies in captive Psittaciformes and the potential impact on endangered wild birds.

A curated dataset of peste des petits ruminants virus sequences for molecular epidemiological analyses

Peste des petits ruminants (PPR) is a highly contagious and devastating viral disease infecting predominantly sheep and goats. Tracking outbreaks of disease and analysing the movement of the virus often involves sequencing part or all of the genome and comparing the sequence obtained with sequences from other outbreaks, obtained from the public databases. However, there are a very large number (>1800) of PPRV sequences in the databases, a large majority of them relatively short, and not always well-documented. There is also a strong bias in the composition of the dataset, with countries with good sequencing capabilities (e.g. China, India, Turkey) being overrepresented, and most sequences coming from isolates in the last 20 years. In order to facilitate future analyses, we have prepared sets of PPRV sequences, sets which have been filtered for sequencing errors and unnecessary duplicates, and for which date and location information has been obtained, either from the database entry or from other published sources. These sequence datasets are freely available for download, and include smaller datasets which maximise phylogenetic information from the minimum number of sequences, and which will be useful for simple lineage identification. Their utility is illustrated by uploading the data to the MicroReact platform to allow simultaneous viewing of lineage date and geographic information on all the viruses for which we have information. While preparing these datasets, we identified a significant number of public database entries which contain clear errors, and propose guidelines on checking new sequences and completing metadata before submission.

Homologous recombination as a mechanism of genetic changes in bovine parainfluenza-3 virus

Bovine parainfluenza-3 virus (BPIV-3) is one of the main viruses associated with bovine respiratory disease complex (BRDC) worldwide. BPIV-3 infect the bovine respiratory tract causing from subclinical infections to severe pneumonia with significant economic losses in the cattle industry. BPIV-3 is a RNA virus with high genetic variability, nevertheless, the contribution of recombination events to its variability has not been assessed so far. In this study the 25 complete genome sequences (CGS) reported so far and 215 partial sequences of different viral genes of BPIV-3 were analyzed to determine their genotypes and subgenotypes, distribution, and the existence of potential recombination events. Based on the analysis of the HN, M, N, and P genes one hypothetical subgenotype was found (subgenotype A4). Four recombination events between sequences of swine and cattle were detected by RDP4 analysis in conjunction with phylogenetic incongruences in the L gene. In addition, 9 sequences reported from Argentina were found to be miss-classified. These results reveal that homologous recombination events have a relevant role in the evolution of BPIV-3 and highlight the importance of implement advanced molecular characterization to better understand the variability and evolution of BPIV-3 as a component of BRDC.

A comprehensive global perspective on phylogenomics and evolutionary dynamics of Small ruminant morbillivirus

A string of complete genome sequences of Small ruminant morbillivirus (SRMV) have been reported from different parts of the globe including Asia, Africa and the Middle East. Despite individual genome sequence-based analysis, there is a paucity of comparative genomic and evolutionary analysis to provide overarching and comprehensive evolutionary insights. Therefore, we first enriched the existing database of complete genome sequences of SRMVs with Pakistan-originated strains and then explored overall nucleotide diversity, genomic and residue characteristics, and deduced an evolutionary relationship among strains representing a diverse geographical region worldwide. The average number of pairwise nucleotide differences among the whole genomes was found to be 788.690 with a diversity in nucleotide sequences (0.04889 ± S.D. 0.00468) and haplotype variance (0.00001). The RNA-dependent-RNA polymerase (L) gene revealed phylogenetic relationship among SRMVs in a pattern similar to those of complete genome and the nucleoprotein (N) gene. Therefore, we propose another useful molecular marker that may be employed for future epidemiological investigations. Based on evolutionary analysis, the mean evolution rate for the complete genome, N, P, M, F, H and L genes of SRMV was estimated to be 9.953 × 10^–4, 1.1 × 10^–3, 1.23 × 10^–3, 2.56 × 10^–3, 2.01 × 10^–3, 1.47 × 10^–3 and 9.75 × 10^–4 substitutions per site per year, respectively. A recombinant event was observed in a Pakistan-originated strain (KY967608) revealing Indian strains as major (98.1%, KR140086) and minor parents (99.8%, KT860064). Taken together, outcomes of the study augment our knowledge and current understanding towards ongoing phylogenomic and evolutionary dynamics for better comprehensions of SRMVs and effective disease control interventions.

Genetic and evolutionary analysis of a new Asia-4 lineage and naturally recombinant canine distemper virus strains from Thailand

Canine distemper virus (CDV), a pathogen causing fatal disease in a wide range of carnivores, can be classified into several geographically-related lineages. It is unclear how genetic recombination contributed to the evolution and emergence of the novel CDV strains and the evolutions of these strains are not fully yet investigated. In this study, the complete genome sequences of eight CDV viruses, isolated from domestic dogs in Thailand, were investigated. Interestingly, most of the identified CDV strains (CDV1-3, -5, -8 TH/2014) clustered as a novel Asia-4 lineage, while the CDV4, -6, -7 TH/2014 belonged to the Asia-1 lineage. Recombination analysis revealed that the CDV4 TH/2014 is a putative recombinant virus from the Asia-1 and America-2 parent viruses. In contrast, no recombination events were detected in the Asia-4 lineage, indicating that it is a distinctive lineage. Evolutionary analysis suggested that the CDV Asia-4 lineage had emerged since 1924 and shared common ancestor with the America-2 lineage. Pressure analysis revealed that CDV nucleotides were under negative selection pressure for its rapid adaptation. These findings demonstrate the evolution of CDV Asia-4 lineage and identified the Asia-1 recombination event. The information regarding genetic diversity of CDVs is essential for further CDV’s research and monitoring.

Constraints on the Genetic and Antigenic Variability of Measles Virus

Antigenic drift and genetic variation are significantly constrained in measles virus (MeV). Genetic stability of MeV is exceptionally high, both in the lab and in the field, and few regions of the genome allow for rapid genetic change. The regions of the genome that are more tolerant of mutations (i.e., the untranslated regions and certain domains within the N, C, V, P, and M proteins) indicate genetic plasticity or structural flexibility in the encoded proteins. Our analysis reveals that strong constraints in the envelope proteins (F and H) allow for a single serotype despite known antigenic differences among its 24 genotypes. This review describes some of the many variables that limit the evolutionary rate of MeV. The high genomic stability of MeV appears to be a shared property of the Paramyxovirinae, suggesting a common mechanism that biologically restricts the rate of mutation.

Influence of vaccine strains on the evolution of canine distemper virus

Canine distemper virus (CDV) is a major dog pathogen belonging to the genus Morbillivirus of the family Paramyxoviridae. CDV causes disease and high mortality in dogs and wild carnivores. Although homologous recombination has been demonstrated in many members of Paramyxoviridae, these events have rarely been reported for CDV. To detect potential recombination events, the complete CDV genomes available in GenBank up to June 2015 were screened using distinct algorithms to detect genetic conversions and incongruent phylogenies. Eight putative recombinant viruses derived from different CDV genotypes and different hosts were detected. The breakpoints of the recombinant strains were primarily located on fusion and hemagglutinin glycoproteins. These results suggest that homologous recombination is a frequent phenomenon in morbillivirus populations under natural replication, and CDV vaccine strains might play an important role in shaping the evolution of this virus.

Development of a helper cell-dependent form of peste des petits ruminants virus: a system for making biosafe antigen

Peste des petits ruminants (PPR) is a viral disease of sheep and goats that is spreading through many countries in the developing world. Work on the virus is often restricted to studies of attenuated vaccine strains or to work in laboratories that have high containment facilities. We have created a helper cell dependent form of PPR virus by removing the entire RNA polymerase gene and complementing it with polymerase made constitutively in a cell line. The resultant L-deleted virus grows efficiently in the L-expressing cell line but not in other cells. Virus made with this system is indistinguishable from normal virus when used in diagnostic assays, and can be grown in normal facilities without the need for high level biocontainment. The L-deleted virus will thus make a positive contribution to the control and study of this important disease.

Epitope dampening monotypic measles virus hemagglutinin glycoprotein results in resistance to cocktail of monoclonal antibodies

The measles virus (MV) is serologically monotypic. Life-long immunity is conferred by a single attack of measles or following vaccination with the MV vaccine. This is contrary to viruses such as influenza, which readily develop resistance to the immune system and recur. A better understanding of factors that restrain MV to one serotype may allow us to predict if MV will remain monotypic in the future and influence the design of novel MV vaccines and therapeutics. MV hemagglutinin (H) glycoprotein, binds to cellular receptors and subsequently triggers the fusion (F) glycoprotein to fuse the virus into the cell. H is also the major target for neutralizing antibodies. To explore if MV remains monotypic due to a lack of plasticity of the H glycoprotein, we used the technology of Immune Dampening to generate viruses with rationally designed N-linked glycosylation sites and mutations in different epitopes and screened for viruses that escaped monoclonal antibodies (mAbs). We then combined rationally designed mutations with naturally selected mutations to generate a virus resistant to a cocktail of neutralizing mAbs targeting four different epitopes simultaneously. Two epitopes were protected by engineered N-linked glycosylations and two epitopes acquired escape mutations via two consecutive rounds of artificial selection in the presence of mAbs. Three of these epitopes were targeted by mAbs known to interfere with receptor binding. Results demonstrate that, within the epitopes analyzed, H can tolerate mutations in different residues and additional N-linked glycosylations to escape mAbs. Understanding the degree of change that H can tolerate is important as we follow its evolution in a host whose immunity is vaccine induced by genotype A strains instead of multiple genetically distinct wild-type MVs.

Interclade recombination in porcine parvovirus strains

A detailed analysis of the Ns1/Vp1Vp2 genome region of the porcine parvovirus (PPV) strains isolated from vaccinated animals was performed. We found many inconsistencies in the phylogenetic trees of these viral isolates, such as low statistical support and strains with long branches in the phylogenetic trees. Thus, we used distance-based and phylogenetic methods to distinguish de facto recombinants from spurious recombination signals. We found a mosaic virus in which the Ns1 gene was acquired from one PPV clade and the Vp1Vp2 gene was acquired from a distinct phylogenetic clade. We also described the interclade mosaic structure of the Vp1Vp2 gene of a reference strain. If recombination is an adaptive mechanism over the course of PPV evolution, we would likely observe increasing numbers of chimeric strains over time. However, when the PPV sequences isolated from 1964 to 2011 were analysed, only two chimeric strains were detected. Thus, PPV recombination is an independent event, resulting from close contact between animals housed in high-density conditions.

Isolation and molecular identification of Sunshine virus, a novel paramyxovirus found in Australian snakes

This paper describes the isolation and molecular identification of a novel paramyxovirus found during an investigation of an outbreak of neurorespiratory disease in a collection of Australian pythons. Using Illumina® high-throughput sequencing, a 17,187 nucleotide sequence was assembled from RNA extracts from infected viper heart cells (VH2) displaying widespread cytopathic effects in the form of multinucleate giant cells. The sequence appears to contain all the coding regions of the genome, including the following predicted paramyxoviral open reading frames (ORFs): 3'--Nucleocapsid (N)--putative Phosphoprotein (P)--Matrix (M)--Fusion (F)--putative attachment protein--Polymerase (L)--5'. There is also a 540 nucleotide ORF between the N and putative P genes that may be an additional coding region. Phylogenetic analyses of the complete N, M, F and L genes support the clustering of this virus within the family Paramyxoviridae but outside both of the current subfamilies: Paramyxovirinae and Pneumovirinae. We propose to name this new virus, Sunshine virus, after the geographic origin of the first isolate--the Sunshine Coast of Queensland, Australia.

Evidence for inter- and intra-clade recombinations in rabies virus

Homologous recombination is considered rare in negative-strand RNA viruses and has not been reported for rabies virus. In this study, full-length genomes of 44 rabies virus strains were analyzed for potential homologous recombination events. Phylogenetic analysis classified these strains into three clades. By applying six different recombination detection methods, one inter-clade and one intra-clade potential recombination events were identified with high confidence values. Software-predicted recombination break points of the two events were all located within the polymerase gene. This report presents the first evidence suggesting the possibility of homologous recombination in rabies virus, which could provide valuable insights for understanding the diversity and evolution of rabies virus as well as other negative-strand RNA viruses.

Homologous recombination in negative sense RNA viruses

Recombination is an important process that influences biological evolution at many different levels. More and more homologous recombination events have been reported among negative sense RNA viruses recently. While sporadic authentic examples indicate that homologous recombination does occur, recombination seems to be generally rare or even absent in most negative sense RNA viruses, and most of the homologous recombination events reported in the literature were likely generated artificially due to lab contamination or inappropriate bioinformatics methods. Homologous recombination in negative sense RNA viruses should be reported with caution in the future, and only after stringent quality control efforts. Moreover, co-infection experiments should be performed to confirm whether recombination can occur.

Respiratory tract immunization of non-human primates with a Newcastle disease virus-vectored vaccine candidate against Ebola virus elicits a neutralizing antibody response

We previously developed a respiratory tract vaccine candidate against Ebola virus (EBOV) based on human parainfluenza virus type 3 (HPIV3), a respiratory paramyxovirus, expressing the EBOV GP envelope protein (HPIV3/GP) from an added gene. Two doses of this vaccine candidate delivered by the intranasal and intratracheal route protected monkeys against intraperitoneal challenge with EBOV; however, concerns exist that the vaccine candidate may have reduced immunogenicity in the adult human population due to pre-existing immunity against HPIV3. Here we developed a new vaccine candidate (NDV/GP) based on Newcastle disease virus (NDV), an avian paramyxovirus that is antigenically distinct from human viral pathogens and is highly attenuated in monkeys. Following one intranasal and intratracheal inoculation of Rhesus monkeys with NDV/GP, titers of EBOV-specific antibodies in respiratory tract secretions and serum samples determined by ELISA, as well as serum EBOV-neutralizing antibodies, were undetectable or low compared to those induced by HPIV3/GP. A second immunization resulted in a substantial boost in serum IgG ELISA titers, yet the titers remained lower than those induced by a second dose of HPIV3/GP. In contrast, the ELISA IgA titers in respiratory tract secretions and, more importantly, the serum EBOV-neutralizing antibody titers were equal to those induced after the second dose of HPIV3/GP. These data suggest that NDV/GP can be effective for immunization against EBOV alone, or in combination with either HPIV3/GP or another vaccine platform in a heterologous prime-boost regimen.

No observed effect of homologous recombination on influenza C virus evolution

The occurrence of homologous recombination in influenza viruses has been under some debate recently. To determine the extent of homologous recombination in influenza C virus, recombination analyses of all available gene sequences of influenza C virus were carried out. No recombination signal was found. With the previous evidence in influenza A and B viruses, it seems that homologous recombination has minimal or no effect on influenza virus evolution.

Isolation and analysis of two naturally-occurring multi-recombination Newcastle disease viruses in China

Two Newcastle disease viruses (NDV), designated QG/Hebei/07 and XD/Shandong/08, were isolated from infected chicken flocks in China in 2007 and 2008, respectively. The results of phylogenetic and recombination analyses on complete NDV genome sequences (excluding terminal segments) show that the QG/Hebei/07 isolate had evidence of recombination in the M and F genes, and recombination in the XD/Shandong/08 isolate in the F, L genes and the non-coding region between the HN and L genes. These two naturally-occurring recombinants we found to be descended from at least three putative parents from vaccine and circulating virus lineages. Moreover, we found that evidence that homologous recombination also occurred between NDV viruses of chicken and swine lineages, while the major putative parent is likely to have been derived from the chicken avirulent vaccine lineage. This study suggests that homologous recombination can occur in all coding and non-coding regions of the NDV genome and a live vaccine strain is capable of recombination with circulating viruses resulting in significant genetic change. The potential role of swine-origin viruses in the evolution of virulent NDV warrants further investigation.

Reassessing conflicting evolutionary histories of the Paramyxoviridae and the origins of respiroviruses with Bayesian multigene phylogenies

The evolution of paramyxoviruses is still poorly understood since past phylogenetic studies have revealed conflicting evolutionary signals among genes, and used varying methods and datasets. Using Bayesian phylogenetic analysis of full length single and concatenated sequences for the 6 genes shared among paramyxovirus genera, we reassess the ambiguous evolutionary relationships within the family, and examine causes of varying phylogenetic signals among different genes. Relative to a pneumovirus outgroup, the concatenated gene phylogeny, splits the Paramyxovirinae into two lineages, one comprising the avulaviruses and rubulaviruses, and a second containing the respiroviruses basal to the henipaviruses, and morbilliviruses. Phylogenies for the matrix (M), RNA dependent RNA polymerase (L) and the fusion (F) glycoprotein genes, are concordant with the topology from the concatenated dataset. In phylogenies derived from the nucleocapsid (N) and phosphoprotein (P) genes, the respiroviruses form the most basal genus of the Paramyxovirinae subfamily, with the avulaviruses and rubulaviruses in one lineage, and the henipaviruses, and morbilliviruses in a second. The phylogeny of the hemagglutinin (H) gene places the respiroviruses basal to the avula-rubulavirus group, but the relationship of this lineage with henipa and morbillviruses is not resolved. Different genes may be under varying evolutionary pressures giving rise to these conflicting signals. Given the level of conservation in the M and L genes, we suggest that together with F gene, these or concatenated datasets for all six genes are likely to reveal the most reliable phylogenies at a family level, and should be used for future phylogenetic studies in this group. Split decomposition analysis suggests that recombination within genera, may have a contributed to the emergence of dolphin morbillivirus, and several species within respiroviruses. A partial L gene alignment, resolves the relationship of 25 unclassified paramxyoviruses into 4 clades (Chiopteran-, Salmon-, Rodentian- and Ophidian paramyxoviruses) which group with rubula-, respiro-, morbilliviruses, and within the paramxyovirinae respectively.

Coinfection rates in Φ6 bacteriophage are enhanced by virus-induced changes in host cells

Two or more viruses infecting the same host cell can interact in ways that profoundly affect disease dynamics and control, yet the factors determining coinfection rates are incompletely understood. Previous studies have focused on the mechanisms that viruses use to suppress coinfection, but recently the phenomenon of enhanced coinfection has also been documented. In the experiments described here, we explore the hypothesis that enhanced coinfection rates in the bacteriophage Φ6 are achieved by virus-induced upregulation of the Φ6 receptor, which is the bacterial pilus. First, we confirmed that coinfection enhancement in Φ6 is virus-mediated by showing that Φ6 attaches significantly faster to infected cells than to uninfected cells. Second, we explored the hypothesis that coinfection enhancement in Φ6 depends upon changes in the expression of an inducible receptor. Consistent with this hypothesis, the closely related phage, Φ12, that uses constitutively expressed lipopolysaccharide as its receptor, attaches to infected and uninfected cells at the same rate. Our results, along with the previous finding that coinfection in Φ6 is limited to two virions, suggest that viruses may closely regulate rates of coinfection through mechanisms for both coinfection enhancement and exclusion.

Cross-species recombination in the haemagglutinin gene of canine distemper virus

Canine distemper virus (CDV) has high prevalence in the world dog population and poses an important conservation threat to many carnivore species. In this study, extensive phylogenetic and recombination analyses were performed on all available complete haemagglutinin gene sequences and a strain (AF178038) isolated from giant panda was identified as putative recombinant. Interestingly, the mosaic was produced by recombination between genotypes European wildlife and Asia-1 and the recombination event involves viruses infecting different host species. This finding may have important implications for the evolution of CDV.

Pathogen evolution and disease emergence in carnivores

Emerging infectious diseases constitute some of the most pressing problems for both human and domestic animal health, and biodiversity conservation. Currently it is not clear whether the removal of past constraints on geographical distribution and transmission possibilities for pathogens alone are sufficient to give rise to novel host-pathogen combinations, or whether pathogen evolution is also generally required for establishment in novel hosts. Canine distemper virus (CDV) is a morbillivirus that is prevalent in the world dog population and poses an important conservation threat to a diverse range of carnivores. We performed an extensive phylogenetic and molecular evolution analysis on complete sequences of all CDV genes to assess the role of selection and recombination in shaping viral genetic diversity and driving the emergence of CDV in non-dog hosts. We tested the specific hypothesis that molecular adaptation at known receptor-binding sites of the haemagglutinin gene is associated with independent instances of the spread of CDV to novel non-dog hosts in the wild. This hypothesis was upheld, providing compelling evidence that repeated evolution at known functional sites (in this case residues 530 and 549 of the haemagglutinin molecule) is associated with multiple independent occurrences of disease emergence in a range of novel host species.

F gene recombination between genotype II and VII Newcastle disease virus

A velogenic Newcastle disease virus (NDV) strain, designated as SRZ03, was isolated from an egg layer flock with NDV vaccine immunization failure in China in 2003. Recombination was found in the F gene of SRZ03. Complete genome sequences analysis indicated that the N-terminal of SRZ03 F gene originated from a genotype II NDV strain, whereas the C-terminal of F gene and the rest of the genes originated from a prevalent velogenic genotype VII NDV strain. It provides us valuable information for understanding the recombination of nonsegmented negative-sense RNA viruses.

An exact nonparametric method for inferring mosaic structure in sequence triplets

Statistical tests for detecting mosaic structure or recombination among nucleotide sequences usually rely on identifying a pattern or a signal that would be unlikely to appear under clonal reproduction. Dozens of such tests have been described, but many are hampered by long running times, confounding of selection and recombination, and/or inability to isolate the mosaic-producing event. We introduce a test that is exact, nonparametric, rapidly computable, free of the infinite-sites assumption, able to distinguish between recombination and variation in mutation/fixation rates, and able to identify the breakpoints and sequences involved in the mosaic-producing event. Our test considers three sequences at a time: two parent sequences that may have recombined, with one or two breakpoints, to form the third sequence (the child sequence). Excess similarity of the child sequence to a candidate recombinant of the parents is a sign of recombination; we take the maximum value of this excess similarity as our test statistic Delta(m,n,b). We present a method for rapidly calculating the distribution of Delta(m,n,b) and demonstrate that it has comparable power to and a much improved running time over previous methods, especially in detecting recombination in large data sets.