Molecular evolution of the Newcastle disease virus matrix protein gene and phylogenetic relationships among the paramyxoviridae

Matrix (M) gene sequences for recent field isolates and older reference Newcastle disease viruses (NDV) were examined to determine phylogenetic relationships and population trends among these viruses. Overall, the M gene has a majority of synonymous nucleotide sequence substitutions occurring among NDV isolates. However, several predicted amino acid changes in the M protein of specific NDV isolates have occurred that correlate to phylogenetic relationships. Nucleotide substitutions in these codons have a greater number of nonsynonymous base changes. The NDV isolates arising since the 1970s belong to a population of viruses that expanded worldwide at an exponential rate. These viruses may have their origins in free-living birds, are present worldwide, and continue to circulate causing disease in poultry. A specific NDV lineage composed of virulent isolates obtained in the US prior to 1970 appears to no longer exists among free-living birds or commercial poultry. However, "vaccine-like" viruses are common in the US and continue to circulate among commercial poultry. Based on M protein amino acid sequences, NDV separates as a clade most closely related to morbilliviruses and not with their current designated category, the rubulaviruses among the Paramyxoviridae. Consequently, avian paramyxoviruses should have their own taxonomic subfamily among the Paramyxovirinae.     

Identification of conserved and variable regions in the envelope glycoprotein sequences of two feline immunodeficiency viruses isolated in Zurich, Switzerland.

The nucleotide sequences of the envelope (env) coding regions of two strains of the feline immunodeficiency virus isolated in Zurich, Switzerland (FIVZ1, FIVZ2) have been analysed. In addition, the complete sequence of the FIVZ1 isolate has been determined. Comparisons have been made with the previously published sequences of three North American isolates (PPR and the Petaluma strains FIV34TF10 and FIV14). The isolate FIVZ1 was very similar to the Petaluma strains of FIV and may represent a clonal derivative acquired by 'contamination'. Overall there are between 2.6% and 15.1% amino acid changes in the env gene products of the five isolates. Of the Zurich isolates, FIVZ2 exhibited the greatest divergence to the other viruses and based on its genotype, phenotype and origins probably represents a new isolate of FIV. Possibly the viruses diverged only recently from a common ancestor. Some 31 of the 33 cysteine residues and 17 of the 21 potential N-linked glycosylation sites of the FIV34TF10 env gene product were conserved among all five isolates. The open reading frame 3 (ORF3, or D) which overlaps the env gene (but is encoded in a different frame) has an ATG codon downstream of a potential splice acceptor site in all five isolates, supporting the view that it encodes a viral gene product. In ORF3 of FIVZ1 a stop codon was located 16 amino acids upstream of the stop codon of ORF3 of the other isolates. The ORF4 (or G) of isolate FIVZ2, thought to be the second coding exon of an FIV rev-like gene, contained a nucleotide deletion in amino acid 45 of ORF4, resulting in a--1 frameshift at this position. Comparison of the LTR sequences of the five isolates identified conserved promoter/enhancer elements. A potential stem-loop structure was identified in the R region of the LTRs of all the isolates, despite the heterogeneity of nucleotide sequences in that region. Such structures (TAR) are present in analogous regions of other lentiviruses and are responsible for tat-mediated trans-activation.     

Analysis of matrix protein gene nucleotide sequence diversity among Newcastle disease virus isolates demonstrates that recent disease outbreaks are caused by viruses of psittacine origin

Nucleotide sequence analysis was completed for isolates of Newcastle disease virus (NDV; avian paramyxovirus 1) from 1992 outbreaks in cormorants and turkeys. These isolates were of the neurotropic velogenic type. The cormorant and turkey NDV isolates had the fusion protein cleavage sequence¹⁰⁹SRGRRQKR/FVG¹¹⁹, as opposed to the consensus sequence¹⁰⁹SGGRRQKR/FIG¹¹⁹ of most known velogenic NDV isolates. The R for G substitution at position 110 may be unique for the cormorant and turkey isolates. For comparative purposes, nucleotide sequencing and analysis of the conserved matrix protein gene coding region were completed for isolates representing all pathotypes. Phylogenetic relationships demonstrated that there are two major groups of NDV isolates. One group includes viruses found in North America and worldwide, such as B1, LaSota, Texas/GB, and Beaudette/C. The second group contains isolates, such as Ulster/2C, Australia/Victoria, and Herts/33, considered exotic to North America. Within this second group are viruses of psittacine origin. The viruses from 1992 outbreaks of Newcastle disease in North America, and an isolate thought to have caused the major outbreak in southern California during the 1970s, are most closely related to an NDV isolate of psittacine origin.     

RNA-dependent RNA polymerase gene analysis of worldwide Newcastle disease virus isolates representing different virulence types and their phylogenetic relationship with other members of the paramyxoviridae

Nucleotide sequence was determined for the RNA-dependent RNA polymerase (L) gene of 16 Newcastle disease virus (NDV) isolates from diverse geographic and chronological origins. The observed consensus amino acid sequence conformed to the six domains previously identified among paramyxovirus L proteins, and the putative 749QGDNQ753 active site was strictly conserved among all isolates. Analysis of predicted amino acid sequences allowed us to identify a sequencing error in the previously reported L genes for NDV. The correct sequences reported herein provided a more accurate alignment with predicted l-amino acid sequences of other paramyxoviruses. Comparison of L gene coding sequences among isolates revealed that synonymous substitutions dominated non-synonymous substitutions, as observed previously with other NDV genes. However, the overall substitution rate was lower than other genes examined making the L gene the most conserved of the NDV genome. Phylogenetic analysis of L genes among NDV isolates was consistent with previous results that suggested the existence of two major lineages. One group contained strains isolated in North America prior to 1970 and included virulent and vaccine strains, while the second group included virulent viruses isolated worldwide. A comparison of the NDV L coding sequences to other Paramyxoviridae illustrated the unique clustering of the avian-specific paramyxoviruses, further justifying the newly created Avulavirus genus.     

Mutational analysis of the Sendai virus V protein: importance of the conserved residues for Zn binding,virus pathogenesis,and efficient RNA editing

The V protein of Sendai virus (SeV) is nonessential for virus replication in cell culture but indispensable for viral pathogenicity in mice. At the C terminus of the V protein, there are amino acid residues conserved among the members of the Paramyxovinae subfamily that are clustered in three regions: region I, just downstream of the RNA editing site; and regions II and III, cysteine-rich zinc-finger-like regions. In the present study, we introduced mutations into the conserved amino acids and generated nine mutant viruses. All of the viruses had impaired virus replication in mouse lungs and attenuated virulence in mice. Furthermore, the C-terminal polypeptides fused with glutathione-S-transferase with a mutation in region I, II, or III all had impaired Zn binding in a (65)Zn-binding assay in solution. These results demonstrate that the conserved amino acids are important for V protein function, probably via protein conformation dependent on Zn binding. One mutant, SeV V-H(318)N, had inefficient RNA editing, indicating that the nucleotide that is a part of the codon encoding histidine at position 318 is conserved for the RNA editing machinery. In addition, to determine the function of the C-terminal extension of the V protein, which is not translated in recent virulent field isolates, a translational stop codon was introduced to generate the corresponding short V protein. The mutant virus showed similar virus propagation and pathogenicity, indicating that C-terminal extension of the V protein is not relevant to virus pathogenesis.     

Use of a heteroduplex mobility assay to detect differences in the fusion protein cleavage site coding sequence among Newcastle disease virus isolates

Newcastle disease virus (NDV) is an economically important pathogen of poultry that may cause clinical disease that ranges from a mild respiratory syndrome to a virulent form with high mortality, depending on an isolate's pathotype. Infections with virulent NDV strains are required to be reported by member nations to the Office of International Epizootes (OIE). The primary determinant for virulence among NDV isolates is the presence or absence of dibasic amino acids in the fusion (F) protein cleavage activation site. Along with biological virulence determinations as the definitive tests, OIE accepts reporting of the F protein cleavage site sequence of NDV isolates as a virulence criterion. Nucleotide sequence data for many NDV isolates recently isolated from infected chickens and other avian species worldwide have been deposited in GenBank. Consequently, viral genomic information surrounding the F protein cleavage site coding sequence was used to develop a heteroduplex mobility assay (HMA) to aid in further identification of molecular markers as predictors of NDV virulence. Using common vaccine strains as a reference, we were able to distinguish virulent viruses among NDV isolates that correlated with phylogenetic analysis of the nucleotide sequence. This technique was also used to examine NDV isolates not previously characterized. We were able to distinguish vaccine-like viruses from other isolates potentially virulent for chickens. This technique will help improve international harmonization of veterinary biologics as set forth by the OIE and the Veterinary International Cooperation on Harmonization of Technical Requirements of Veterinary Medicinal Products. Ultimately, the HMA could be used for initial screening among a large number of isolates and rapid identification of potentially virulent NDV that continue to threaten commercial poultry worldwide.     

Analysis of codon usage in Newcastle disease virus

In this study, the relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents, and dinucleotide were used to investigate codon usage pattern of each protein-coding gene and genome among 31 Newcastle disease virus (NDV) isolates. The result shows that the overall extent of codon usage bias in NDV is low (mean ENC = 56.15 > 40). The good correlation between the (C + G)₁₂% and (G + C)₃% suggests that the mutational pressure, rather than natural selection, is the main factor that determines the codon usage bias and base component in NDV. It is observed that synonymous codon usage pattern in NDV genes is gene function and geography specific, but not host specific. By contrasting synonymous codon usage patterns of different NDV isolates, we suggest that more than one genotype of NDV circulates in waterfowl in USA; and gene length has no significant effect on the variations of synonymous codon usage in these virus genes. CpG under-represented is a characteristic for NDV to fit in its host. These results not only provide an insight into the variation of codon usage pattern among the genomes of NDV, but also may help in understanding the processes governing the evolution of NDV.     

Nucleotide sequence analysis of the Newcastle disease virus nucleocapsid protein gene and phylogenetic relationships among the Paramyxoviridae

The nucleocapsid (N) protein genes from 24 Newcastle disease virus (NDV) isolates representing various pathotypes with different geographical and chronological origins were cloned and sequenced. The N-terminal region of the N protein to residue 401 was highly conserved among isolates with several conservative substitutions occurring that correlated with phylogenetic relationships. Variability of the N protein was detected in the C-terminal portion similar to what has been reported for other members of the Paramyxovirinae. Amino acids previously identified as invariant or highly conserved in N proteins of other paramyxoviruses were also present in the NDV protein. Phylogenetic analysis of N gene coding sequences among NDV isolates again demonstrated the existence of two major groups. One clade contained viruses that included vaccine and virulent strains isolated in the USA prior to 1970 while a second clade included vaccine and virulent viruses isolated worldwide. Comparison of N protein amino acid sequences among members of the Paramyxoviridae resulted in NDV and avian paramyxovirus 6 separating as a cluster distinct from the Rubulavirus genus. This provides further support for avian paramyxoviruses being considered for their own genus among the Paramyxovirinae.     

Glycoprotein evolution of vesicular stomatitis virus New Jersey

A T1 ribonuclease fingerprinting study of a large number of virus isolates had previously demonstrated that considerable genetic variability existed among natural isolates of the vesicular stomatitis virus (VSV) New Jersey (NJ) serotype [S.T. Nichol (1988) J. Virol. 62, 572-579]. Based on these results, 34 virus isolates were chosen as representing the extent of genetic diversity within the VSV NJ serotype. We report the entire glycoprotein (G) gene nucleotide sequence and the deduced amino acid sequence for each of these viruses. Up to 19.8% G gene sequence differences could be seen among NJ serotype isolates. Analysis of the distribution of nucleotide substitutions relative to nucleotide codon position revealed that third position changes were distributed randomly throughout the gene. Third base changes constituted 84% of the observed nucleotide substitutions and affected 89% of the third base positions located in the G gene. Only three short oligonucleotide stretches of complete sequence conservation were observed. The remaining nucleotide changes located in the first and second positions were not distributed randomly, indicating that most of the amino acids coded by the G gene cannot be altered without reducing the fitness of the VSV NJ serotype viruses. Despite these constraints, up to 8.5% amino acid differences were observed between virus isolates. These differences were located throughout the G protein including regions adjacent to defined major antibody neutralization epitopes. Apparent clusters of amino acid substitutions were present in the hydrophobic signal sequence, transmembrane domain, and within the cytoplasmic domain of the G protein. A maximum parsimony analysis of the G gene nucleotide sequences allowed construction of a phylogram indicating the evolutionary relationship of these viruses. The VSV NJ serotype appears to contain at least three distinct lineages or subtypes. All recent virus isolates from the United States and Mexico are within subtype I and appear to have evolved from an ancestor more closely related to the Hazelhurst historic strain than other older strains. The implications of these findings for the evolution, epizootiology, and classification of these viruses are discussed.     

Genotypic characteristics of HTLV-II isolates from Amerindian and non-Indian populations

The Amerindian human T-cell lymphotropic virus type II isolate HTLV-II_G12 has been demonstrated to be an HTLV-IIb with several unique features, including several restriction enzyme site changes, a distinctive pre-gag region, a stop codon within thepol gene, and an extended Tax protein. In this study, HTLV-II isolates from Amerindian and non-Indian populations were characterized by restriction enzyme site analysis to determine the prevalent HTLV-II subtype. In addition, DNA amplification by the polymerase chain reaction and Southern blot analyses were used to probe for the HTLV-II_G12 pre-gag region. Our results showed that of 13 Guaymi Indian isolates subtyped, all were HTLV-IIb, and that approximately one third of 17 isolates had the unique pre-gag region. While other HTLV-II-infected groups contained both HTLV-IIa and HTLV-IIb isolates, none of these isolates showed evidence of the distinctive HTLV-II_G12 pre-gag region. Lastly, DNA sequence analysis was used to determine the prevalence of the stop codon within thepol gene open reading frame. These analyses revealed that the occurrence of a stop codon within this sequence appeared to be characteristic of most HTLV-IIb subtypes. These results further our understanding of the genetic variations and evolution of the HTLV-II viruses within the endemically infected Amerindian populations, as well as U.S. intravenous drug users and other non-Indian populations.     

Molecular Characterization of Virulent Newcastle Disease Virus Isolates from Chickens during the 1998 NDV Outbreak in Kazakhstan

Newcastle disease virus (NDV) infects domesticated and wild birds throughout the world and has the possibility to cause outbreaks in chicken flocks in future. To assess the evolutionary characteristics of 10 NDV strains isolated from chickens in Kazakhstan during 1998 we investigated the phylogenetic relationships among these viruses and viruses described previously. For genotyping, fusion (F) gene phylogenetic analysis (nucleotide number 47-421) was performed using sequences of Kazakhstanian isolates as compared to sequences of selected NDV strains from GenBank. Phylogenetic analysis showed that all newly characterized strains belonged to the genetic group designated as VIIb. All strains possessed a virulent fusion cleavage site (RRQRR/F) belonging to velogenic or mesogenic pathotypes with intracerebral pathogenicity indexes (ICPI) varying from 1.05 to 1.87.     

Molecular Cloning and Sequence Analysis of the Phosphoprotein (P) Gene of the Lapinized Rinderpest Virus

We determined the nucleotide sequence of the coding region for the phosphoprotein (P) gene of the L strain of rinderpest virus (RPV). The gene encodes two overlapping open reading frames of 1521 and 531 nucleotides. Use of the first ATG would produce a P polypeptide of 507 amino acids, while use of the second ATG would produce a C polypeptide of 177 amino acids. In addition, the insertion of an extra G residue at the editing site generates an alternative mRNA potentially encoding the V protein of RPV. Homology comparisons of the P, C and V proteins among various viruses suggest that RPV is closer to measles virus (MV) than to canine distemper virus (CDV). Alignment of the sequences unique to the V protein revealed that the cysteine residues are well conserved among RPV, MV and CDV, and form a zinc finger-like motif.     

Presence of Newcastle disease viruses of sub-genotypes Vc and VIn in backyard chickens and in apparently healthy wild birds from Mexico in 2017

Virulent Newcastle disease viruses (NDV) have been present in Mexico since 1946, and recently, multiple outbreaks have been reported in the country. Here, we characterized eleven NDV isolated from apparently healthy wild birds and backyard chickens in three different locations of Jalisco, Mexico in 2017. Total RNA from NDV was reverse-transcribed, and 1285 nucleotides, which includes 3/4 of the fusion gene, was amplified and sequenced using a long-read MinION sequencing method. The sequences were 99.99–100% identical to the corresponding region obtained using the Illumina MiSeq. Phylogenetic analysis using MinION sequences demonstrated that nine virulent NDV from wild birds belonged to sub-genotypes Vc and VIn, and two backyard chicken isolates were of sub-genotype Vc. The sub-genotype Vc viruses had nucleotide sequence identity that ranged from 97.7 to 98% to a virus of the same sub-genotype isolated from a chicken in Mexico in 2010. Three viruses from pigeons had 96.3–98.7% nucleotide identity to sub-genotype VIn pigeon viruses, commonly referred to as pigeon paramyxovirus, isolated in the USA during 2000–2016. This study demonstrates that viruses of sub-genotype Vc are still present in Mexico, and the detection of this sub-genotype in both chickens and wild birds suggests that transmission among these species may represent a biosecurity risk. This is the first detection and complete genome sequencing of genotype VI NDV from Mexico. In addition, the utilization of an optimized long-read sequencing method for rapid virulence and genotype identification using the Oxford nanopore MinION system is demonstrated.

     

Distinct patterns of evolution between respiratory syncytial virus subgroups A and B from New Zealand isolates collected over thirty-seven years

Respiratory syncytial virus (RSV) is the most important cause of viral lower respiratory tract infections in infants and children worldwide. In New Zealand, infants with RSV disease are hospitalized at a higher rate than other industrialized countries, without a proportionate increase in known risk factors. The molecular epidemiology of RSV in New Zealand has never been described. Therefore, we analyzed viral attachment glycoprotein (G) gene sequences from 106 RSV subgroup A isolates collected in New Zealand between 1967 and 2003, and 38 subgroup B viruses collected between 1984 and 2004. Subgroup A and B sequences were aligned separately, and compared to sequences of viruses isolated from other countries during a similar period. Genotyping and clustering analyses showed RSV in New Zealand is similar and temporally related to viruses found in other countries. By quantifying temporal clustering, we found subgroup B viruses clustered more strongly than subgroup A viruses. RSV B sequences displayed more variability in stop codon usage and predicted protein length, and had a higher degree of predicted O-glycosylation site changes than RSV A. The mutation rate calculated for the RSV B G gene was significantly higher than for RSV A. Together, these data reveal that RSV subgroups exhibit different patterns of evolution, with subgroup B viruses evolving faster than A.     

Efficient expression of the 15-kDa form of infectious pancreatic necrosis virus VP5 by suppression of a UGA codon

Infectious pancreatic necrosis virus (IPNV), a member of the Birnaviridae family, encodes a nonstructural VP5 protein from a small open reading frame (ORF), which overlaps with a major ORF encoding pVP2, VP4 and VP3 proteins. In majority of the Sp strains of IPNV sequenced to date, VP5 gene codes for a 15-kDa protein. However, we have shown that in highly virulent strains, there is a premature in-frame stop codon (UGA) at nucleotide (nt) position 427, (preceding the 15-kDa stop codon at nt position 511) which could encode a 12-kDa protein. Using reverse genetics, we recovered recombinant rNVI15, rNVI15-15K and rNVI15-DeltaVP5 viruses (which could encode 12 or 15-kDa VP5 or lack the expression of VP5, respectively) and demonstrated that VP5 is dispensable for viral replication in vivo but is not involved in virulence (Santi, N., Song, H., Vakharia, V. N., Evensen, ?., 2005a. Infectious pancreatic necrosis virus VP5 is dispensable for virulence and persistence. J. Virol. 79, 9206-9216). Here, we utilized these viruses to investigate the gene expression of VP5 in vitro. Our results indicate that a 15-kDa VP5 is produced in rNVI15-infected cells, albeit at lower levels than in rNVI15-15K-infected cells, suggesting that the opal stop codon at nt 427 is suppressed. Furthermore, to examine translational suppression of the opal stop codon in VP5 gene, we constructed plasmids containing VP5-specific sequence and employed a yeast-based bicistronic dual-luciferase reporter system (Harger, J.W., Dinman, J.D., 2003. An in vivo dual-luciferase assay system for studying translational recoding in the yeast Saccharomyces cerevisiae. RNA 9, 1019-1024). Our results demonstrate that the VP5 sequence (with or without a stop codon) yielded approximately 13% termination suppression and the efficiency is directly related to the base immediately 3' of the termination codon, C>A>U>G.     

Newcastle disease outbreaks in Kazakhstan and Kyrgyzstan during 1998, 2000, 2001, 2003, 2004, and 2005 were caused by viruses of the genotypes VIIb and VIId

Newcastle disease virus (NDV) infects domesticated and wild birds throughout the world, and infections with virulent NDV strains continue to cause disease outbreaks in poultry and wild birds. To assess the evolutionary characteristics of 28 NDV strains isolated from chickens in Kazakhstan and Kyrgyzstan during 1998, 2000, 2001, 2003, 2004, and 2005, we investigated the phylogenetic relationships among these viruses and viruses described previously. For genotyping, fusion (F) gene phylogenetic analysis (nucleotide number 47–421) was performed using sequences of Kazakhstanian and Kyrgyzstanian isolates as compared to sequences of selected NDV strains from GenBank. Phylogenetic analysis demonstrated that the 14 newly characterized strains from years 1998 to 2001 belong to the NDV genotype VIIb, whereas the 14 strains isolated during 2003–2005 were of genotype VIId. All strains possessed a virulent fusion protein cleavage site (R-R-Q-R/K-R-F) and had intracerebral pathogenicity indexes in day-old chickens that ranged from 1.05 to 1.87, both properties typical of NDV strains classified in the mesogenic or velogenic pathotype.