key: cord-0003117-s3zfz24h
authors: Bhat, Sudipta; Kattoor, Jobin Jose; Malik, Yashpal Singh; Sircar, Shubhankar; Deol, Pallavi; Rawat, Vinita; Rakholia, Ritu; Ghosh, Souvik; Vlasova, Anastasia N.; Nadia, Touil; Dhama, Kuldeep; Kobayashi, Nobumichi
title: Species C Rotaviruses in Children with Diarrhea in India, 2010–2013: A Potentially Neglected Cause of Acute Gastroenteritis
date: 2018-02-17
journal: Pathogens
DOI: 10.3390/pathogens7010023
sha: f7658f780f7b559a2ddb1649957efda8a94f534b
doc_id: 3117
cord_uid: s3zfz24h

All over the world, children and adults are severely affected by acute gastroenteritis, caused by one of the emerging enteric pathogens, rotavirus C (RVC). At present, no extensive surveillance program is running for RVC in India, and its prevalence is largely unknown except cases of local outbreaks. Here, we intended to detect the presence of RVC in diarrheic children visiting or admitted to hospitals in Haldwani (state of Uttarakhand, India), a city located in the foothills of the Himalayas. During 2010–2013, we screened 119 samples for RVC by an RVC VP6 gene-specific RT-PCR. Of these, 38 (31.93%) were found positive, which is higher than the incidence rates reported so far from India. The phylogenetic analysis of the derived nucleotide sequences from one of the human RVC (HuRVC) isolates, designated as HuRVC/H28/2013/India, showed that the study isolate belongs to genotype I2, P2 and E2 for RVC structural genes 6 and 4 (VP6, and VP4) and non-structural gene 4 (NSP4), respectively. Furthermore, the VP6 gene of HuRVC/H28/2013/India shows the highest similarity to a recently-reported human-like porcine RVC (PoRVC/ASM140/2013/India, KT932963) from India suggesting zoonotic transmission. We also report a full-length NSP4 gene sequence of human RVC from India. Under the One-health platforms there is a need to launch combined human and animal RVC surveillance programs for a better understanding of the epidemiology of RVC infections and for implementing control strategies.

picobirnavirus in a child suffering with gastroenteritis [41] and further the presence of porcine RVC in India [14] . The previously discussed report suggests the possibilities of gene re-assortment events through interspecies transmission [42] . Molecular epidemiological studies support the detection of HuRVCs from neighboring countries of India like Bangladesh [43] , Russia [44] , Malaysia [45] , South Korea [46] , Thailand, Nepal [47] and Kenya [48] , supporting the importance of comprehensive surveillance in the Indian population. As mixed infections of different strain of RVs are a prerequisite for re-assortment events, co-surveillance of animal and human virus strains is vital to gain a better understanding of the relationships between co-circulating viruses.

During 2010-2013, 119 fecal samples were screened for RVC from Haldwani (state of Uttarakhand, India), a city known as the gate way of the Himalayan foothills, located at 29 • 13'0" N 79 • 31'0" E with an average elevation of 424 m (1391 feet), having a mixed urban and rural culture of various religions with diverse living systems. Hospital-based (Out Patient Department patients/hospitalized, less than one year of age) diarrheic samples were collected with the respective parent's consent. A 10% (w/v) fecal suspension was prepared with phosphate-buffered saline (pH 7.2) followed by clarification of debris by centrifugation at 6797× g for 5 min, and the supernatants were used for the viral RNA extraction.

Total RNA from suspended diarrheic stool samples was isolated using Qiazol reagent (Qiagen GmbH, Hilden, Germany) following the manufacturer's protocol and quantified on a Nanodrop spectrometer [49] . Extracted RNA (500 ng) was used to prepare a pool of first-strand cDNA by random-priming reverse transcription (RT) using recombinant MMLV-RT (Promega Corporation, Madison, WI, USA) and random hexamer (Qiagen GmbH, Hilden, Germany) at 37 • C. Enzymatic activity of the MMLV-RT was stopped by keeping the cDNA reaction mixture at 80 • C for 3 min.

Presence of RVC in diarrheic stool samples was detected by using the RT-PCR assay targeting the VP6 gene with primers RVC-VP6-DF; 5 -ARTCHGTTCTATGYGATTC-3 [14] and BMJ44; 5 -AGCCACATAG TTCACATTTC-3 [20] (Supplementary Data 1). The diagnostic primers were expected to amplify a 340-bp amplicon. PCR conditions were optimized as described previously [14] with SapphireAmpFast PCR master mix (Takara Bio Inc., Shiga, Japan). These samples were also subjected to screening for other infectious viral causes of gastroenteritis viz. calicivirus, astrovirus and picobirnaviruses following the primers and procedures outlined in the published reports [50] [51] [52] .

Characterization of the selected RVC strain (HuRVC/H28/2013/India) was done by amplification of the structural capsid genes (VP6, VP4 and VP7) and non-structural enterotoxin gene (NSP4). Primers used for the amplification are mentioned in the Supplementary Data 1. The published PCR conditions for the amplification of individual genes were followed viz. for the VP6 gene [53] , VP4 and VP7 genes [43] and NPS4 gene [14] . RT-PCR amplicons were visualized under a UV trans-illuminator after resolving in 1.5% agarose gel electrophoresis. Desired amplified products were cloned into the pDRIVE (Qiagen GmbH, Hilden, Germany, Figure 1 ) cloning vector and transformed in Escherichia coli DH5α competent cells [54] . The recombinant plasmids were obtained using the GeneJET plasmid Miniprep kit (Thermo Fisher Scientific, Waltham, MA, Vilnius, Lithuania). Positive recombinant clones were sequenced by the Bigdye terminator Sanger sequencing method in an ABI 3730 x l sequencer (Eurofins Genomic India Ltd., Bangalore, India). The sequences were deposited in GenBank with Accession Numbers MG553198 (NSP4 gene), MG553199 (VP4 gene) and MG553200 (VP6 gene).

GeneJET plasmid Miniprep kit (Thermo Fisher Scientific, Waltham, MA, Vilnius, Lithuania). Positive recombinant clones were sequenced by the Bigdye terminator Sanger sequencing method in an ABI 3730 x l sequencer (Eurofins Genomic India Ltd., Bangalore, India). The sequences were deposited in GenBank with Accession Numbers MG553198 (NSP4 gene), MG553199 (VP4 gene) and MG553200 (VP6 gene). 

The pair-wise similarity among the nucleotide/amino acid sequences were calculated for the VP6, VP4 and NSP4 genes after aligning the sequences by the Clustal V program in MegAlign software of the DNASTAR software package. For the genetic relatedness study, representative fulllength genes (VP6: n = 47; NSP4: n = 26; VP4: n = 30) of human, cattle, pig, dog and ferret RVC strains were retrieved from the NCBI database (Supplementary Data 2) [55] . Phylogenetic analysis was performed using the maximum likelihood method (1000 bootstrap replicates) in MEGA 6 software (v 6.06) [56] . The suitable dendrogram analysis model was identified as described earlier [13] using the find best DNA/protein model tool available in MEGA 6 (v 6.06), confirmed with the FindModel online tool [57] .

The N-glycosylation sites on the VP6, VP4 and NSP4 proteins of the HuRVC/H28/2013/India strain were predicted by using a NetNGlyc 1.0 server [58] . These sites were compared with available human RVC sequences (Supplementary Data 2).

Among the total 119 samples, 60 were found positive for RVA (unpublished data), and one sample out of 60 was found with picobirnavirus genogroup I (PBV GG-I), RVB and RVA co-infection [41] . Of the, 119 samples, 38 (31.93%) were found positive for RVC in the VP6 gene-based RT-PCR ( Figure 2 ). Preceding reports from different Asian countries showed both sporadic and outbreak cases with comparatively less prevalence of RVC, which may be to the lack of sensitive detection assay availability [43] [44] [45] [46] [47] [48] . In the previously mentioned reports from India, serological prevalence of 0.43% and 25.32% was reported by Brown et al. [39] and Mukhopadhya et al. [40] , respectively, from southern India (Vellore). The percent molecular detection of RVC by RT-PCR was 8.6% (in the outbreak cases) and 0.7% (sporadic cases) from western India [1] and 3.33% from Delhi [37] . Moreover, the molecular detection of RVCs throughout the world shows a range between 0.3% and 23.7% [1] .The present study foresees a high prevalence (31.93%) and supports a diverse epidemiology of RVC in north India among the children of less than one year of age in contrast to the earlier studies, which could emphasize the occurrence of RVC in adults [36] . Notably, none of the sample detected as positive for RVC showed co-infection with screened viruses (calicivirus, astrovirus and picobirnavirus). 

The pair-wise similarity among the nucleotide/amino acid sequences were calculated for the VP6, VP4 and NSP4 genes after aligning the sequences by the Clustal V program in MegAlign software of the DNASTAR software package. For the genetic relatedness study, representative full-length genes (VP6: n = 47; NSP4: n = 26; VP4: n = 30) of human, cattle, pig, dog and ferret RVC strains were retrieved from the NCBI database (Supplementary Data 2) [55] . Phylogenetic analysis was performed using the maximum likelihood method (1000 bootstrap replicates) in MEGA 6 software (v 6.06) [56] . The suitable dendrogram analysis model was identified as described earlier [13] using the find best DNA/protein model tool available in MEGA 6 (v 6.06), confirmed with the FindModel online tool [57] .

The N-glycosylation sites on the VP6, VP4 and NSP4 proteins of the HuRVC/H28/2013/India strain were predicted by using a NetNGlyc 1.0 server [58] . These sites were compared with available human RVC sequences (Supplementary Data 2).

Among the total 119 samples, 60 were found positive for RVA (unpublished data), and one sample out of 60 was found with picobirnavirus genogroup I (PBV GG-I), RVB and RVA co-infection [41] . Of the, 119 samples, 38 (31.93%) were found positive for RVC in the VP6 gene-based RT-PCR ( Figure 2 ). Preceding reports from different Asian countries showed both sporadic and outbreak cases with comparatively less prevalence of RVC, which may be to the lack of sensitive detection assay availability [43] [44] [45] [46] [47] [48] . In the previously mentioned reports from India, serological prevalence of 0.43% and 25.32% was reported by Brown et al. [39] and Mukhopadhya et al. [40] , respectively, from southern India (Vellore). The percent molecular detection of RVC by RT-PCR was 8.6% (in the outbreak cases) and 0.7% (sporadic cases) from western India [1] and 3.33% from Delhi [37] . Moreover, the molecular detection of RVCs throughout the world shows a range between 0.3% and 23.7% [1] .The present study foresees a high prevalence (31.93%) and supports a diverse epidemiology of RVC in north India among the children of less than one year of age in contrast to the earlier studies, which could emphasize the occurrence of RVC in adults [36] . Notably, none of the sample detected as positive for RVC showed co-infection with screened viruses (calicivirus, astrovirus and picobirnavirus).

One sample (designated as HuRVC/H28/2013/India) out of the 38 RVC positive samples got amplified partially for the VP6 and VP4 genes and the complete CDS (coding DNA sequence) region for the NSP4 gene using the reported primers (Supplementary Data 1). The absence of intact virus particles or degraded viral RNA (repeated freeze thawing of biological sample or inherent nucleases present in fecal matter) may be a cause for the inability to amplify full-length VP6 and VP4 gene of the virus with otherwise well-established primers. As observed in other studies, cross-priming or primer-binding failures at the 3 end of the primer may also cause the inability to obtain the full-length genome sequence [59, 60] . 

One sample (designated as HuRVC/H28/2013/India) out of the 38 RVC positive samples got amplified partially for the VP6 and VP4 genes and the complete CDS (coding DNA sequence) region for the NSP4 gene using the reported primers (Supplementary Data 1). The absence of intact virus particles or degraded viral RNA (repeated freeze thawing of biological sample or inherent nucleases present in fecal matter) may be a cause for the inability to amplify full-length VP6 and VP4 gene of the virus with otherwise well-established primers. As observed in other studies, cross-priming or primer-binding failures at the 3' end of the primer may also cause the inability to obtain the fulllength genome sequence [59, 60] .

The amplified VP6 gene of the HuRVC/H28/2013/India strain has 831 nucleotides (521st-1352nd position), encoding a protein of 228 amino acids in length. We failed to get the sequence stretch from the 5'-end of the VP6 gene. The phylogenetic tree (n = 47) of the VP6 gene consisted of two major clades ( Figure 3 ). One contained bovine, ferret and canine RVCs and the other consisted of porcine, human and human-like porcine RVCs. As per the cut-off value established for genotype determination, the HuRVC/H28/2013/India strain clustered with human RVCs (I2) along with all other Indian human RVCs, showing their common ancestry [32] . The study isolate clustered along with Indian human origin (AY795898, AY770980, AY786571, AY786570; Barman and Naik [61] ), and KT900227, KT900231, KT900236; [1] ) and human-like porcine RVCs (KX374486, KX374489, KX374492, KT932962, KT932963; [14] ) in the I2 lineage. Although coming in the same I2 lineage, some of the Indian human RVC sequences reported recently from the Delhi region, Northern India (KY886474-79; [37] ) showed comparatively distant phylogenetic relatedness with all other previously-reported Indian RVC sequences. Further, our results corroborated the findings of Tiku and co-workers [37] where an independent evolution of human RVCs of Indian origin at the nucleotide level based on the VP6 gene was noted. The HuRVC/H28/2013/India strain showed the highest sequence similarity with recently-isolated Indian human-like porcine RVC from Assam, PoRVC/ASM140/2013/India (nt: 99.9%, aa: 99.1%) and PoRVC/ASM132/2013/India (nt: 98.7%, aa: 95.7%) (Supplementary Data 3) . The presence of a similar VP6 gene in porcine and human RVC strains could suggest that the gene may have been reasserted and could indicate zoonotic or anthroponotic transmission. Interspecies transmission of RVC from pig to children has been reported in Brazil [42, 58] . In concordance with a previous Indian study [1] , the current study isolate clustered with Indian and Bangladeshi RVC isolates, showing a percent of identity of 97.7-98.5% (nt), and 99.1% (aa), respectively. However, in 

The amplified VP6 gene of the HuRVC/H28/2013/India strain has 831 nucleotides (521st-1352nd position), encoding a protein of 228 amino acids in length. We failed to get the sequence stretch from the 5'-end of the VP6 gene. The phylogenetic tree (n = 47) of the VP6 gene consisted of two major clades (Figure 3 ). One contained bovine, ferret and canine RVCs and the other consisted of porcine, human and human-like porcine RVCs. As per the cut-off value established for genotype determination, the HuRVC/H28/2013/India strain clustered with human RVCs (I2) along with all other Indian human RVCs, showing their common ancestry [32] . The study isolate clustered along with Indian human origin (AY795898, AY770980, AY786571, AY786570; Barman and Naik [61] ), and KT900227, KT900231, KT900236; [1] ) and human-like porcine RVCs (KX374486, KX374489, KX374492, KT932962, KT932963; [14] ) in the I2 lineage. Although coming in the same I2 lineage, some of the Indian human RVC sequences reported recently from the Delhi region, Northern India (KY886474-79; [37] ) showed comparatively distant phylogenetic relatedness with all other previously-reported Indian RVC sequences. Further, our results corroborated the findings of Tiku and co-workers [37] where an independent evolution of human RVCs of Indian origin at the nucleotide level based on the VP6 gene was noted. The HuRVC/H28/2013/India strain showed the highest sequence similarity with recently-isolated Indian human-like porcine RVC from Assam, PoRVC/ASM140/2013/India (nt: 99.9%, aa: 99.1%) and PoRVC/ASM132/2013/India (nt: 98.7%, aa: 95.7%) (Supplementary Data 3) . The presence of a similar VP6 gene in porcine and human RVC strains could suggest that the gene may have been reasserted and could indicate zoonotic or anthroponotic transmission. Interspecies transmission of RVC from pig to children has been reported in Brazil [42, 58] . In concordance with a previous Indian study [1] , the current study isolate clustered with Indian and Bangladeshi RVC isolates, showing a percent of identity of 97.7-98.5% (nt), and 99.1% (aa), respectively. However, in contrast to another study from India [37] where the isolates were more similar to Nigerian isolates, our isolate (HuRVC/H28/2013/India) showed lower similarity to Nigerian isolates (nt: 95.5-97.1%, aa: 98.7%) (Supplementary Data 3) . Sequence identity with other animal species is mentioned in Table 1 , where the current study isolate shows high sequence similarity with other human RVCs at the nucleotide and amino acid level (>97%). Of note, the percent similarity was low with animal-origin VP6 sequences (bovine, ferret and canine, <90%), except for porcine, wherein a wide range of similarity was seen contrast to another study from India [37] where the isolates were more similar to Nigerian isolates, our isolate (HuRVC/H28/2013/India) showed lower similarity to Nigerian isolates (nt: 95.5-97.1%, aa: 98.7%) (Supplementary Data 3) . Sequence identity with other animal species is mentioned in Table 1, where the current study isolate shows high sequence similarity with other human RVCs at the nucleotide and amino acid level (>97%). Of note, the percent similarity was low with animal-origin VP6 sequences (bovine, ferret and canine, <90%), except for porcine, wherein a wide range of similarity was seen (84-99%). This shows close relatedness of the current study's HuRVC strain to the porcine-origin RVC strain. Further, no significant change in the available amino acid sequence of the VP6 gene from HuRVC/H28/2013/India supports the conserved nature of this species-specific protein. Figure 3 . Phylogenetic analysis of rotavirus C based on VP6 genes (831 nucleotide) of human, swine, cattle, dog and ferret origin at the nucleotide level. The Tamura-3+G algorithm was identified using the find best DNA/protein model tool available in MEGA 6, which was confirmed with the FindModel online tool [57] . Numbers on branches indicate the percentages of bootstrap support from 1000 Figure 3 . Phylogenetic analysis of rotavirus C based on VP6 genes (831 nucleotide) of human, swine, cattle, dog and ferret origin at the nucleotide level. The Tamura-3+G algorithm was identified using the find best DNA/protein model tool available in MEGA 6, which was confirmed with the FindModel online tool [57] . Numbers on branches indicate the percentages of bootstrap support from 1000 replicates. VP6 gene-based I typing of RVC is denoted along with clusters. Host species depicted are human (Hu); cattle (Bo); pig (Por)); dog (Can); and ferret. Strains/isolates are represented according to their host species, accession number, gene, strain, year of isolation and country of origin. The isolate of the current study is denoted by the solid dot. Genotype designations are shown on the right. The scale bar indicates nucleotide substitutions per site. The amplified VP4 gene of the HuRVC/H28/2013/India isolate comprised of a 1224-nucleotide (1st-1224th position) long stretch encoded a polypeptide of 408 amino acids. Upon phylogenetic analysis with representative sequences of the RVC VP4 gene (n = 30) from different host species, the dendrogram scattered into two major clades, among which one clade contained bovine, canine and porcine RVCs, whereas the another contained all human RVCs ( Figure 4) . As per the cut-off values stated for genotype determination, the HuRVC/H28/2013/India isolate clustered with human RVCs of the P2 type. Sequence similarity revealed that the study sequence showed the least sequence identity (nt: 94.3%, aa: 95.6%) with the Japanese isolate (LC129057) and highest (nt: 97%, aa: 97.6%) with the Bangladeshi isolate (HQ185635) as determined through the sequence distance in MegAlign software (Supplementary Data 4). Other Indian RVC VP4 strains reported in the study of Tiku et al. [37] were not taken into consideration for analyzing the percentage of identity due to the small nucleotide length. Analyzing the current study strain with the earlier Indian study sequence may lead to misinterpretation of the obtained data. Even though the current study strain shared the same clade as the sequences reported by Tiku and co-workers [37] , the isolate of the present study showed a higher sequence identity (96.4-96.9%) (Supplementary Data 4) with other RVC VP4 gene sequences from India, which is not in accordance with the observation by the former report. According to the report, their isolates were more similar to Korean RVCs; however, the isolate HuRVC/H28/2013/India showed more identity to the Bangladeshi and Indian isolates. Though Indian isolates possess a common ancestor, by virtue of the divergence in nucleotide, they have been placed in different sub-clades ( Figure 4) . Additionally, upon analysis, six point mutations were observed in the amino acid sequence of the current study isolate, i.e., I31T, V73I, A124T, N141D, I146L and N192S, which may have a role in shifting the antigenicity of the circulating RVCs ( Figure 5 ). However, further work is needed to elucidate their position on the protein surface. Sequence identity of the VP4 gene with other animal species is mentioned in Table 1 , where the current study strain shows quite high similarity (95.6-97.6%) at the amino acid level with other human RVCs as compared to animal RVCs (60.4-68.7%). This further proves the human origin of our study strain. Further, more divergence in the outer capsid protein VP4 gene in comparison to inner capsid protein VP6 among Indian RVC isolates proves the rapid evolving nature to escape selection pressure.

The RVC strains NSP4 gene sequenced in this study has a length of 472 nucleotides, encoding a protein of 150 amino acids in length (complete ORF). The NSP4 gene of human RVC is smaller than that of RVA (175 amino acid) [62] and larger than that of bovine RVC (146 amino acid) [63] The sequence identity at the nucleotide and amino acid level revealed that the study isolate had 91.4-96% and 90.2-95.1% identity, respectively, with the available human RVC sequences. The highest sequence identity (nt: 96%, aa: 95.1%) was observed with Russian isolates (KP735977 and JN969079) and the least identity with Korean isolate Chungnam (nt: 91.4%, aa: 90.2) (Supplementary Data 5). Upon phylogenetic analysis (n = 26), all RVC isolates clustered into two major clades, i.e., one which contains human RVCs exclusively and another containing RVC isolated from animals (bovine, porcine and canine) ( Figure 6 ). As per the cut-off value established for the E genotype, HuRVC/H28/2013/India clustered with human RVCs of the E2 genotype. In comparison to the RVC isolates in the animal RVC clade, the isolates in the human RVC clade were more similar. In animal RVC clades, further sub-clades were seen (bovine, canine and porcine subclades). The study strain showed higher genetic relatedness to Eurasian RVC strains (strains from Hungary, Russia, the U.K. and Bangladesh) than to the more geographically-close Chinese RVC strains ( Figure 6 ). The amino acid sequence analysis of the NSP4 gene from HuRVC/H28/2013/India showed point mutations like D2E, H75N, H100Q, I120M and H133Y (Figure 7) . Although, the NSP4 gene's functional role has not been deciphered yet for RVC, the changes at position 120 and 133 might affect the functional property of the NSP4 gene as the amino acid sequence from 114 to 135 is considered responsible for enterotoxin activity and the binding capability of virus with caveolin, which help in the budding of the virus particles in RVA infection [62] . Nevertheless, this has to be proven through site-directed mutagenesis studies. These two point mutations (I120M and H133Y) lie in the well-described VP4 binding domain (described for RVA strains) [64] . This report describes the full-length NSP4 gene from an Indian human RVC. Sequence identity of the VP4 gene with other animal species is mentioned in Table 1 , as for the VP4 gene, the NSP4 gene of the current study also shows higher amino acid identity (90.2-95.4%) to other human RVCs as compared to animal RVCs (54.9-62.1%). 

The RVC strains NSP4 gene sequenced in this study has a length of 472 nucleotides, encoding a protein of 150 amino acids in length (complete ORF). The NSP4 gene of human RVC is smaller than that of RVA (175 amino acid) [62] and larger than that of bovine RVC (146 amino acid) [63] The sequence identity at the nucleotide and amino acid level revealed that the study isolate had 91.4-96% and 90.2-95.1% identity, respectively, with the available human RVC sequences. The highest sequence identity (nt: 96%, aa: 95.1%) was observed with Russian isolates (KP735977 and JN969079) and the least identity with Korean isolate Chungnam (nt: 91.4%, aa: 90.2) (Supplementary Data 5). Upon phylogenetic analysis (n = 26), all RVC isolates clustered into two major clades, i.e., one which contains human RVCs exclusively and another containing RVC isolated from animals (bovine, porcine and canine) ( Figure 6 ). As per the cut-off value established for the E genotype, HuRVC/H28/2013/India clustered with human RVCs of the E2 genotype. In comparison to the RVC isolates in the animal RVC clade, the isolates in the human RVC clade were more similar. In animal RVC clades, further sub-clades were seen (bovine, canine and porcine subclades). The study strain showed higher genetic relatedness to Eurasian RVC strains (strains from Hungary, Russia, the U.K. and Bangladesh) than to the more geographically-close Chinese RVC strains ( Figure 6 ). The amino acid sequence analysis of the NSP4 gene from HuRVC/H28/2013/India showed point mutations like D2E, H75N, H100Q, I120M and H133Y (Figure 7) . Although, the NSP4 gene's functional role has not been deciphered yet for RVC, the changes at position 120 and 133 might affect the functional property of the NSP4 gene as the amino acid sequence from 114 to 135 is considered responsible for enterotoxin activity and the binding capability of virus with caveolin, which help in the budding of the virus particles in RVA infection [62] . Nevertheless, this has to be proven through site-directed mutagenesis studies. These two point mutations (I120M and H133Y) lie in the well-described VP4 binding domain (described for RVA strains) [64] . This report describes the full-length NSP4 gene from an Indian human RVC. Sequence identity of the VP4 gene with other animal species is mentioned in Table 1 , as for the VP4 gene, the NSP4 gene of the current study also shows higher amino acid identity (90.2-95.4%) to other human RVCs as compared to animal RVCs (54.9-62.1%). Tamura 

In the current study, we failed to amplify the VP7 gene from any of the RVC-positive samples. A similar observation was noticed by Tiku et al. [37] , where they could only amplify 50% of the isolate's VP7 and VP4 genes. This may be due to the sequence variations observed in the circulating Indian RVC strains; however, newer approaches viz. adaptor ligation amplification could solve this problem.

The prediction and analysis of N-linked glycosylation have gained intense interest in recently due to the association of its relevance with the overall understanding of the virus biology, its ability to impart various advantages to virus survival, folding of the proteins, replication, maturation processes and also virulence. Here, we intended to know the possible glycosylation sites by analyzing the human RVC strains using the NetNGlyc 1.0 Server [58] . In the VP6 protein analysis, we predicted N-glycosylation at positions 36, 64, 76, 89 and 205, and these remained unaltered in the VP6 gene of the current study isolate as well; while in the VP4 protein, there was a deletion of the N-glycosylation site (N192S) ( Figure 5 ). Being an outer capsid protein, the predicted deletion in glycosylation (N192S) may alter the virus attachment and entry into the host cell, but needs further studies in this context. The NSP4 in RVs is a heavily glycosylated viral protein [62] , and that was also observed during the current analysis of available human RVC NPS4 sequences, where conserved N-glycosylation sites were seen in the NSP4 protein at 5-7 (NQT), 32-35 (NGS) and 65-67 (NNS) [64] . This study provides a preliminary addition to the available knowledge, which in particular is lacking with respect to Indian RVC isolates.

A number of studies on RVC from all over the world have shown growing concern about RVC as one of the causes of severe acute gastroenteritis in human and animals [43, 46, 65, 66] . RVCs have been confirmed in cow [26] , ferrets [27] and dogs [28] as the cause of both sporadic and outbreak disease. The prevalence data on RVC in humans from the Indian subcontinent is limited. In this study, we report a high presence (31.93%) of RVC in contrast to the earlier reports, which ascertains RVC as an important diarrheagenic pathogen. The genetic analysis of inner capsid protein VP6 and outer capsid VP4 revealed divergence from earlier Indian human RVC isolates. Overall, the findings of this study establish the presence of RVC in the human population of northern India along with insights Figure 7 . Deduced amino acid sequence alignment of the NSP4 genes of the current study HuRVC strain with that of previously-described HuRVC strains. The dots represent the identical amino acids, and five changes were noticed at positions 2, 75, 100, 120, 133 in the current study's RVC isolate, indicated in red boxes.

In the current study, we failed to amplify the VP7 gene from any of the RVC-positive samples. A similar observation was noticed by Tiku et al. [37] , where they could only amplify 50% of the isolate's VP7 and VP4 genes. This may be due to the sequence variations observed in the circulating Indian RVC strains; however, newer approaches viz. adaptor ligation amplification could solve this problem.

The prediction and analysis of N-linked glycosylation have gained intense interest in recently due to the association of its relevance with the overall understanding of the virus biology, its ability to impart various advantages to virus survival, folding of the proteins, replication, maturation processes and also virulence. Here, we intended to know the possible glycosylation sites by analyzing the human RVC strains using the NetNGlyc 1.0 Server [58] . In the VP6 protein analysis, we predicted N-glycosylation at positions 36, 64, 76, 89 and 205, and these remained unaltered in the VP6 gene of the current study isolate as well; while in the VP4 protein, there was a deletion of the N-glycosylation site (N192S) ( Figure 5 ). Being an outer capsid protein, the predicted deletion in glycosylation (N192S) may alter the virus attachment and entry into the host cell, but needs further studies in this context. The NSP4 in RVs is a heavily glycosylated viral protein [62] , and that was also observed during the current analysis of available human RVC NPS4 sequences, where conserved N-glycosylation sites were seen in the NSP4 protein at 5-7 (NQT), 32-35 (NGS) and 65-67 (NNS) [64] . This study provides a preliminary addition to the available knowledge, which in particular is lacking with respect to Indian RVC isolates.

A number of studies on RVC from all over the world have shown growing concern about RVC as one of the causes of severe acute gastroenteritis in human and animals [43, 46, 65, 66] . RVCs have been confirmed in cow [26] , ferrets [27] and dogs [28] as the cause of both sporadic and outbreak disease. The prevalence data on RVC in humans from the Indian subcontinent is limited. In this study, we report a high presence (31.93%) of RVC in contrast to the earlier reports, which ascertains RVC as an important diarrheagenic pathogen. The genetic analysis of inner capsid protein VP6 and outer capsid VP4 revealed divergence from earlier Indian human RVC isolates. Overall, the findings of this study establish the presence of RVC in the human population of northern India along with insights about the similarity with porcine RVC (based on the VP6 gene), indicating a potential source and an interspecies transmission event. In view of the emerging importance of RVC infection all around the world, there is a need to conclusively define the classification system for RVC as well. A few of the studies have used the classification criterion adopted for RVA, but need approval of the RCWG before their global acceptance. In the future, it is essential to look into RVC prevalence through extensive surveillance program both in human and animal populations to understand the geographical distribution of RVC infections in India.

Supplementary Materials: The following are available online at http://www.mdpi.com/2076-0817/7/1/23/s1.

Group C rotavirus infection in patients with acute gastroenteritis in outbreaks in western India between

Epidemiology of rotavirus in the Iranian children: A systematic review and meta-analysis

Molecular epidemiology of enteric adenovirus gastroenteritis in under-five-year-old children in Iran

Molecular epidemiology of noroviruses detected in Vietnamese children with acute gastroenteritis from 2012 to 2015

Fields Virology

VP6-sequence-based cutoff values as a criterion for rotavirus species demarcation

Virus Taxonomy: Eight Report of the International Committee on Taxonomy of Viruses; Fauquet

Candidate new rotavirus species in Schreiber's bats, Serbia

Whole genomic characterization of a human rotavirus strain B219 belonging to a novel group of the genus Rotavirus

Porcine pararotavirus: Detection, differentiation from rotavirus, and pathogenesis in gnotobiotic pigs

Avian rotavirus enteritis-An updated review

Unexpected detection of porcine rotavirus C strains carrying human origin VP6 gene

Molecular evidence of group D rotavirus in commercial broiler chicks in India

Evidence for Occurrence of Human group B rotavirus in Central India Based on Characterization of NSP2 Gene

Group D Rotaviruses: Structure, Epidemiology, Diagnosis, and Perspectives on Future Research Challenges

Candidate new rotavirus species in sheltered dogs

Rotavirus I in feces of a cat with diarrhea

First detection of group C rotavirus in children with acute diarrhea in Spain

Review of group A rotavirus strains reported in swine and cattle

Genomic diversity among group A rotaviruses from diarrheic children, piglets, buffalo and cow calves of Madhya Pradesh

Incidence and genetic diversity of group C rotavirus among adults

Rotavirus-like, calicivirus-like, and 23-nm virus-like particles associated with diarrhea in young pigs

Detection of a rotavirus-like agent associated with diarrhea in an infant

Detection of a bovine group C rotavirus from adult cows with diarrhea and reduced milk production

Detection of Group C Rotavirus in Juvenile Ferrets (Mustela putorius furo) with Diarrhea by Reverse Transcription Polymerase Chain Reaction: Sequencing and Analysis of the Complete Coding Region of the VP6 Gene

Canine rotavirus C strain detected in Hungary shows marked genotype diversity

Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains

Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments

RotaC: A web-based tool for the complete genome classification of group A rotaviruses

Genetic diversity of the VP7, VP4 and VP6 genes of Korean porcine group C rotaviruses

Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG)

Identification, phylogenetic analysis and classification of porcine group C rotavirus VP7 sequences from the United States and Canada

Whole-genome characterization of human group C rotaviruses: Identification of two lineages in the VP3 gene

An occurrence of diarrheal cases associated with group C rotavirus in adults

First study conducted in Northern India that identifies group C rotavirus as the etiological agent of severe diarrhea in children in Delhi

First detection of group C rotavirus in fecal specimens of children with diarrhea in the United States

Rotavirus epidemiology in Vellore, south India: Group, subgroup, serotype, and electrophoretype

Anti-VP6 IgG antibodies against group A and group C rotaviruses in South India

Molecular evidences of mixed viral infection in an infant hospitalized with acute gastroenteritis in northern India

Evidence for zoonotic transmission of group C rotaviruses among children in Belém, Brazil

Detection and characterization of human group C rotaviruses in Bangladesh

First genetic characterization of rotavirus C in Russia

Identification of a human group C rotavirus in Malaysia

First detection of group C rotavirus in children with acute gastroenteritis in South Korea

Group C rotavirus infections in patients with diarrhea in Thailand, Nepal, and England

Human group C rotaviruses identified in Kenya

Multispecies reassortant bovine rotavirus strain carries a novel simian G3-like VP7 genotype

Designand evaluation of a primer pair that detects both Norwalk-and Sapporo-like caliciviruses by RT-PCR

Human stool contains a previously unrecognized diversity of novel astroviruses

Cloning of human picobirnavirus genomic segments and development of an RT-PCR detection assay

The correct sequence of the porcine group C/Cowden rotavirus major inner capsid protein shows close homology with human isolates from Brazil and the

Role for nonstructural protein 1 of severe acute respiratory syndrome coronavirus in chemokine dysregulation

National Center for Biotechnology Information

Molecular Evolutionary Genetics Analysis Version 6.0

Testing the model of DNA substitution

Prediction of glycosylation across the human proteome and the correlation to protein function

Rotavirus genotyping: Keeping up with an evolving population of human rotaviruses

Designing and validation of new primers for efficient genotyping of animal G3 rotaviruses

Human rotavirus C strain V508 major inner capsid protein VP6 (VP6) gene, complete cds

Rotavirus Nonstructural Protein 4 (NSP4)-Viral Enterotoxin with Multiple roles in Pathogenesis of Diarrhoea in Children

Whole-genome analysis of two bovine rotavirus C strains:Shintoku and Toyama

Phylogenetic analysis of human group C rotavirus circulating in Brazil reveals a potential unique NSP4 genetic variant and high similarity with Asian strains

Detection and Quantification of Group C Rotaviruses in Communal Sewage

Phylogenetic characterization of VP6 gene (inner capsid) of porcine rotavirus C collected in Japan

Acknowledgments: All the authors acknowledge and thank their respective institutes, hospital and universities. The help of Sharad Saurabh and Kuldeep Sharma is acknowledged.Author Contributions: All the authors substantially contributed to the conception, design, analysis and interpretation of data, drafting and revising, updating, checking and approving the final version of manuscript. We all agree to be accountability for the content of this manuscript.

All authors declare that there exist no commercial or financial relationships that could in any way lead to a potential conflict of interest.