Nucleotide sequence analysis of rotavirus gene 11 from two tissue culture-adapted ATCC strains, RRV and Wa

[Variability of genes encoding nonstructural proteins of rotavirus А (Reoviridae: Rotavirus: Rotavirus A) genotype G9P[8] during the period of dominance in the territory of Nizhny Novgorod (central part of Russia) (2011-2020)]

INTRODUCTION

In Russia, rotavirus A is the main cause of severe viral gastroenteritis in young children. The molecular features that allow a rotavirus of a particular genotype to gain an evolutionary advantage remain unclear, therefore, the study of the genetic diversity of rotaviruses based on genes encoding nonstructural proteins (NSPs) responsible for the reproduction of the virus in the cell is an urgent task.

OBJECTIVE

To study the genetic diversity of rotaviruses of genotype G9P[8], which dominated Nizhny Novgorod in 20112020, based on genes encoding nonstructural proteins.

MATERIALS AND METHODS

Rotavirus-positive samples were subjected to PCR-genotyping and sequencing of NSP1 NSP5 genes. Phylogenetic analysis was carried out in the MEGA X program.

RESULTS

In the period 20112020, G9P[8] rotaviruses with four variants of the NSP2 gene were co-circulating in Nizhny Novgorod. New alleles were noted in 2012 (N1-a-III), 2016 (N1-a-IV) and in 2019 (N1-a-II). The appearance of new variants of other genes occurred in 2014 (E1-3, NSP4), 2018 (T1-a3-III, NSP3) and in 2019 (A1-b-II, NSP1). NSP2 gene had the most variable amino acid sequence (16 substitutions), 2 to 7 substitutions were observed in NSP1, NSP3 and NSP4, NSP5 was conservative.

DISCUSSION

The results obtained are consistent with the literature data and indicate the participation of NSP genes in maintaining the heterogeneity of the rotavirus population.

CONCLUSION

Until 2018, the genetic diversity of rotaviruses in Nizhny Novgorod was determined by the circulation of strains carrying several alleles of the NSP2 gene and conservative genes NSP1, NSP3NSP5. By the end of the study period, new variants of the genotype G9P[8] were formed in the population, carrying previously unknown combinations of alleles of nonstructural genes.

Characterization of a second open reading frame in genome segment 10 of bluetongue virus

Viruses have often evolved overlapping reading frames in order to maximize their coding capacity. Until recently, the segmented dsRNA genome of viruses of the Orbivirus genus was thought to be monocistronic, but the identification of the bluetongue virus (BTV) NS4 protein changed this assumption. A small ORF in segment 10, overlapping the NS3 ORF in the +1 position, is maintained in more than 300 strains of the 27 different BTV serotypes and in more than 200 strains of the phylogenetically related African horse sickness virus (AHSV). In BTV, this ORF (named S10-ORF2 in this study) encodes a putative protein 50–59 residues in length and appears to be under strong positive selection. HA- or GFP-tagged versions of S10-ORF2 expressed from transfected plasmids localized within the nucleoli of transfected cells, unless a putative nucleolar localization signal was mutated. S10-ORF2 inhibited gene expression, but not RNA translation, in transient transfection reporter assays. In both mammalian and insect cells, BTV S10-ORF2 deletion mutants (BTV8ΔS10-ORF2) displayed similar replication kinetics to wt virus. In vivo, S10-ORF2 deletion mutants were pathogenic in mouse models of disease. Although further evidence is required for S10-ORF2 expression during infection, the data presented provide an initial characterization of this ORF.

FRESCo: finding regions of excess synonymous constraint in diverse viruses

Background

The increasing availability of sequence data for many viruses provides power to detect regions under unusual evolutionary constraint at a high resolution. One approach leverages the synonymous substitution rate as a signature to pinpoint genic regions encoding overlapping or embedded functional elements. Protein-coding regions in viral genomes often contain overlapping RNA structural elements, reading frames, regulatory elements, microRNAs, and packaging signals. Synonymous substitutions in these regions would be selectively disfavored and thus these regions are characterized by excess synonymous constraint. Codon choice can also modulate transcriptional efficiency, translational accuracy, and protein folding.

Results

We developed a phylogenetic codon model-based framework, FRESCo, designed to find regions of excess synonymous constraint in short, deep alignments, such as individual viral genes across many sequenced isolates. We demonstrated the high specificity of our approach on simulated data and applied our framework to the protein-coding regions of approximately 30 distinct species of viruses with diverse genome architectures.

Conclusions

FRESCo recovers known multifunctional regions in well-characterized viruses such as hepatitis B virus, poliovirus, and West Nile virus, often at a single-codon resolution, and predicts many novel functional elements overlapping viral genes, including in Lassa and Ebola viruses. In a number of viruses, the synonymously constrained regions that we identified also display conserved, stable predicted RNA structures, including putative novel elements in multiple viral species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0603-7) contains supplementary material, which is available to authorized users.

Phylogenetic analysis of a G6P[5] bovine rotavirus strain isolated in a neonatal diarrhea outbreak in a beef cattle herd vaccinated with G6P[1] and G10P[11] genotypes

The aim of this study was to perform the molecular characterization of the eleven genes of a G6P[5] bovine group A rotavirus (RVA) strain detected in a diarrhea outbreak from a vaccinated beef cattle herd. The outbreak affected 80 % of calves between 15–30 days old. RVA was identified by RT-PCR in 12 (70.6 %) out of 17 diarrheic fecal samples evaluated. The rotavirus wild-type strain had the genotype constellation G6(IV)-P[5](IX)-I2c-R2-C2-M2-A3-N2-T6-E2e-H3a. This study confirms the importance of homotypic immunity against the bovine RVA P[5] genotype in neonatal diarrhea in cattle herds that are regularly vaccinated against rotaviruses.

Whole-genomic analysis of rotavirus strains: current status and future prospects

Studies on genetic diversity of rotaviruses have been primarily based on the genes encoding the antigenically significant VP7 and VP4 proteins. Since the rotavirus genome has 11 segments of RNA that are vulnerable to reassortment events, analyses of the VP7 and VP4 genes may not be sufficient to obtain conclusive data on the overall genetic diversity, or true origin of strains. In the last few years following the advent of the whole-genome-based genotype classification system, the whole genomes of at least 167 human group A rotavirus strains have been analyzed, providing a plethora of new and important information on the complex origin of strains, inter- and intra-genogroup reassortment events, animal–human reassortment events, zoonosis, and genetic linkages involving different group A rotavirus gene segments. In addition, the whole genomes of a limited number of human group B, C and novel group rotavirus strains have been analyzed. This article briefly reviews the available data on whole-genomic analysis of human rotavirus strains. The significance and future prospects of whole-genome-based studies are also discussed.

Molecular characterization of rare bovine group A rotavirus G15P[11] and G15P[21] strains from eastern India: identification of simian SA11-like VP6 genes in G15P[21] strains

During a surveillance study (November 2001-March 2005), one rare G15P[11] and two rare G15P[21] bovine group A rotavirus strains were detected in diarrhoeic calves in Eastern India. Sequence analysis of the VP8*, VP6, NSP4 and NSP5 genes of the G15P[11] strain confirmed its bovine origin. Although the NSP4 and NSP5 genes of the two G15P[21] strains were of bovine origin, their VP6 genes shared higher nucleotide and amino acid identities with simian strain SA11 (92.5-93.1% and 98.5-98.7%) than bovine strains (88.5-88.9% and 97-97.2%), and by phylogenetic analysis, exhibited clustering with SA11, distantly related to bovine strains. All these pointed towards a possible reassortment event of VP6 gene between bovine and simian (SA11-like) strains. Therefore, the present study provided molecular evidence for bovine origin of G15 strains and revealed a rare instance of genetic diversity in the bovine VP6 gene, otherwise conserved in group A rotavirus strains from cattle.

Molecular characterization of a porcine Group A rotavirus strain with G12 genotype specificity

A porcine Group A rotavirus strain (RU172) was detected and molecularly characterized during a surveillance study conducted for rotavirus infection in a pig farm located in a suburban area of Kolkata City, India. The G12 genotype specificity of RU172 was revealed by PCR-based genotyping assays following addition of a G12 type-specific primer (designed in our laboratory to pick up G12 isolates from field samples) and was confirmed by sequence analysis of the VP7-encoding gene. The RU172 strain exhibited maximum VP7 identities of 93.6% to 94.5% with human G12 strains at the deduced amino acid level. In spite of its G12 genotype nature, RU172 appeared to be distinct from human G12 rotaviruses and, on phylogenetic analysis, formed a separate lineage with human G12 strains. Among the other gene segments analyzed, RU172 belonged to NSP4 genotype B, had a NSP5 and VP6 of porcine origin, and shared maximum VP4 identities with porcine P[7] rotaviruses (94.3%-95.4% at the deduced amino acid level). Therefore, to the best of our knowledge, this is the first report of detection of an animal rotavirus strain with G12 genotype specificity. Detection of strains like RU172 provides vital insights into the genomic diversity of Group A rotaviruses of man and animals.

Hyperphosphorylation of the rotavirus NSP5 protein is independent of serine 67, [corrected] NSP2, or [corrected] the intrinsic insolubility of NSP5 is regulated by cellular phosphatases

The NSP5 protein is required for viroplasm formation during rotavirus infection and is hyperphosphorylated into 32- to 35-kDa isoforms. Earlier studies reported that NSP5 is not hyperphosphorylated without NSP2 coexpression or deleting the NSP5 N terminus and that serine 67 is essential for NSP5 hyperphosphorylation. In this report, we show that full-length NSP5 is hyperphosphorylated in the absence of NSP2 or serine 67 and demonstrate that hyperphosphorylated NSP5 is predominantly present in previously unrecognized cellular fractions that are insoluble in 0.2% sodium dodecyl sulfate. The last 68 residues of NSP5 are sufficient to direct green fluorescent protein into insoluble fractions and cause green fluorescent protein localization into viroplasm-like structures; however, NSP5 insolubility was intrinsic and did not require NSP5 hyperphosphorylation. When we mutated serine 67 to alanine we found that the NSP5 mutant was both hyperphosphorylated and insoluble, identical to unmodified NSP5, and as a result serine 67 is not required for NSP5 phosphorylation. Interestingly, treating cells with the phosphatase inhibitor calyculin A permitted the accumulation of soluble hyperphosphorylated NSP5 isoforms. This suggests that soluble NSP5 is constitutively dephosphorylated by cellular phosphatases and demonstrates that hyperphosphorylation does not direct NSP5 insolubility. Collectively these findings indicate that NSP5 hyperphosphorylation and insolubility are completely independent parameters and that analyzing insoluble NSP5 is essential for studies assessing NSP5 phosphorylation. Our results also demonstrate the involvement of cellular phosphatases in regulating NSP5 phosphorylation and indicate that in the absence of other rotavirus proteins, domains on soluble and insoluble NSP5 recruit cellular kinases and phosphatases that coordinate NSP5 hyperphosphorylation.

Comparative molecular characterization of gene segment 11-derived NSP6 from lamb rotavirus LLR strain used as a human vaccine in China

Sequence-length polymorphism is known for rotavirus genetic segment 11 (encodes non-structural protein, NSP6). With the exception of 11 strains that have the coding potential for a 98-residue NSP6, majority of the strains have the potential for a 92-residue NSP6. In nine strains, the coding potential for this protein is even shorter. This report focuses on the NSP6 gene nucleotide sequence of Lanzhou Lamb Rotavirus (LLR) strain and its comparative molecular characterization. The LLR strain is a G10 P12 type, which is in use as a licensed human vaccine in China. The LLR NSP6 was compared with 56 other rotaviral NSP6 sequences including a rhesus strain (RRV) available in the database. Analyses indicate that while RRV-NSP6 belongs to the majority (92-residue) group, the LLR NSP6 belongs to the 98-residue group. When the rotavirus NSP6 protein was expressed in cells as GFP fusion protein from human, simian and the LLR strains, they all demonstrated punctate cytoplasmic distribution and, contrary to the computer-aided prediction, the NSP6 did not undergo phosphorylation, which in itself is a novel observation for the rotavirus NSP6.

Relationships among porcine and human P[6] rotaviruses: Evidence that the different human P[6] lineages have originated from multiple interspecies transmission events

Porcine rotavirus strains (PoRVs) bearing human-like VP4 P[6] gene alleles were identified. Genetic characterization with either PCR genotyping or sequence analysis allowed to determine the VP7 specificity of the PoRVs as G3, G4, G5 and G9, and the VP6 as genogroup I, that is predictive of a subgroup I specificity. Sequence analysis of the VP8* trypsin-cleavage product of VP4 allowed PoRVs to be characterized further into genetic lineages within the P[6] genotype. Unexpectedly, the strains displayed significantly higher similarity (up to 94.6% and 92.5% at aa and nt level, respectively) to human M37-like P[6] strains (lineage I), serologically classifiable as P2A, or to the atypical Hungarian P[6] human strains (HRVs), designated as lineage V (up to 97.0% aa and 96.1% nt), than to the porcine P[6] strain Gottfried, lineage II (<85.1% aa and 82.2 nt), which is serologically classified as P2B. Interestingly, no P[6] PoRV resembling the original prototype porcine strain, Gottfried, was detected, while Japanase P[6] PoRV clustered with the atypical Japanase G1 human strain AU19. By analysis of the 10th and 11th genome segments, all the strains revealed a NSP4B genogroup (Wa-like) and a NSP5/6 gene of porcine origin. These findings strongly suggest interspecies transmission of rotavirus strains and/or genes, and may indicate the occurrence of at least 3 separate rotavirus transmission events between pigs and humans, providing convincing evidence that evolution of human rotaviruses is tightly intermingled with the evolution of animal rotaviruses.

Reduced expression of the rotavirus NSP5 gene has a pleiotropic effect on virus replication

Rotavirus RRV gene 11 encodes two non-structural proteins, NSP5 and NSP6. NSP5 is a phosphorylated non-structural protein that binds single- and double-stranded RNA in a non-specific manner. Transient expression of this protein in uninfected cells has provided evidence for its participation in the formation of electron-dense cytoplasmic structures, known as viroplasms, which are thought to be key structures for the replication of the virus. NSP6 is a protein of unknown function that seems not to be essential for virus replication in cell culture. To study the function of NSP5 in the context of a viral infection, the expression of RRV gene 11 was silenced by RNA interference. Reduction in the synthesis of NSP5, as shown by immunoblot and immunofluorescence assays, correlated with a reduction in the number and size of viroplasms and with an altered intracellular distribution of other viroplasm-associated proteins. Silencing of gene 11 also resulted in a reduced synthesis of viral RNA(+) and double-stranded RNA and of all viral proteins, as well as in a decreased production of infectious virus. A similar phenotype was observed when the NSP5 coding gene of the lapine rotavirus strain Alabama was silenced. The fact that the NSP5 gene of rotavirus Alabama lacks the AUG initiator codon for a complete NSP6 protein, suggests that the described phenotype in gene 11-silenced cells is mostly due to the absence of NSP5. The data presented in this work suggest that NSP5 is a key protein during the replication cycle of rotaviruses.

Sequence analysis of the VP7 and VP4 genes identifies a novel VP7 gene allele of porcine rotaviruses, sharing a common evolutionary origin with human G2 rotaviruses

During an epidemiological survey encompassing several porcine herds in Saragoza, Spain, the VP7 and VP4 of a rotavirus-positive sample, 34461-4, could not be predicted by using multiple sets of G- and P-type-specific primers. Sequence analysis of the VP7 gene revealed a low amino acid (aa) identity with those of well-established G serotypes, ranging between 58.33% and 88.88%, with the highest identity being to human G2 rotaviruses. Analysis of the VP4 gene revealed a P[23] VP4 specificity, as its VP8* aa sequence was 95.9% identical to that of the P14[23],G5 porcine strain A34, while analysis of the VP6 indicated a genogroup I, that is predictive of subgroup I specificity. Analysis of the 10th and 11th RNA segments revealed close identity to strains of porcine and human origin, respectively. The relatively low overall aa sequence conservation (<89% aa) to G2 human rotaviruses, the lack of N-glycosylation sites that are usually highly conserved in G2 rotaviruses, and the presence of several amino acid substitutions in the major antigenic hypervariable regions hampered an unambiguous classification of the porcine strain 34461-4 as G2 serotype on the basis of sequence analysis alone. The identification of a borderline, G2-like, VP7 gene allele in pigs, while reinforcing the hypotheses of a tight relationship in the evolution of human and animal rotaviruses, provides additional evidence for the wide genetic/antigenic diversity of group A rotaviruses.

Molecular analysis of the VP7, VP4, VP6, NSP4, and NSP5/6 genes of a buffalo rotavirus strain: identification of the rare P[3] rhesus rotavirus-like VP4 gene allele

We report the detection and molecular characterization of a rotavirus strain, 10733, isolated from the feces of a buffalo calf affected with diarrhea in Italy. Strain 10733 was classified as a P[3] rotavirus, as the VP8* trypsin cleavage product of the VP4 protein revealed a high amino acid identity (96.2%) with that of rhesus rotavirus strain RRV (P5B[3]), used as the recipient virus in the human-simian reassortant vaccine. Analysis of the VP7 gene product revealed that strain 10733 possessed G6 serotype specificity, a type common in ruminants, with an amino acid identity to G6 rotavirus strains ranging from 88 to 98%, to Venezuelan bovine strain BRV033, and Hungarian human strain Hun4. Phylogenetic analysis based on the VP7 gene of G6 rotaviruses identified at least four lineages and an apparent linkage between each lineage and the VP4 specificity, suggesting the occurrence of repeated interspecies transmissions and genetic reassortment events between ruminant and human rotaviruses. Moreover, strain 10733 displayed a bovine-like NSP4 and NSP5/6 and a subgroup I VP6 specificity, as well as a long electropherotype pattern. The detection of the rare P[3] genotype in ruminants provides additional evidence for the wide genetic and antigenic diversity of group A rotaviruses.

The N- and C-terminal regions of rotavirus NSP5 are the critical determinants for the formation of viroplasm-like structures independent of NSP2

Molecular events and the interdependence of the two rotavirus nonstructural proteins, NSP5 and NSP2, in producing viroplasm-like structures (VLS) were previously evaluated by using transient cellular coexpression of the genes for the two proteins, and VLS domains as well as the NSP2-binding region of NSP5 were mapped in the context of NSP2. Review of the previous studies led us to postulate that NSP2 binding of NSP5 may block the N terminus of NSP5 or render it inaccessible and that any similar N-terminal blockage may render NSP5 alone capable of producing VLS independent of NSP2. This possibility was addressed in this report by using two forms of NSP5-green fluorescent protein (GFP) chimeras wherein GFP is fused at either the N or the C terminus of NSP5 (GFP-NSP5 and NSP5-GFP) and evaluating their VLS-forming capability (by light and electron microscopy) and phosphorylation and multimerization potential independent of NSP2. Our results demonstrate that NSP5 alone can form VLS when the N terminus is blocked by fusion with a nonrotavirus protein (GFP-NSP5) but the C terminus is unmodified. Only GFP-NSP5 was able to undergo hyperphosphorylation and multimerization with the native form of NSP5, emphasizing the importance of an unmodified C terminus for these events. Deletion analysis of NSP5 mapped the essential signals for VLS formation to the C terminus and clearly suggested that hyperphosphorylation of NSP5 is not required for VLS formation. The present study emphasizes in general that when fusion proteins are used for functional studies, constructs that represent fusions at both the N and the C termini of the protein should be evaluated.