INTRODUCTION
Group A rotavirus is the most common cause of viral gastroenteritis in young children and animals worldwide [Reference Gunn1–Reference Badaracco3]. The virus particle has an 11-segment double-stranded RNA genome encoding six structural viral proteins (VP1–4, VP6, VP7) and six non-structural proteins (NSP 1–5/6) [Reference Estes, Kapikian, Knipe, Howley, Griffin, Lamb, Martin, Roizman and Straus4]. The two outer capsid proteins, VP4 and VP7 are used to classify rotavirus strains into P (protease-sensitive) and G (glycoprotein) serotypes, respectively, based on neutralization determinants located on the surface of VP4 and VP7 [Reference Estes5, Reference Gentsch6]. To date, 37 P types and 27 G types have been reported, of these, 15 P types and 12 G types have been identified in humans and animals [Reference Cui7, Reference Kang8].
There have been persistent efforts made towards the development of effective rotavirus vaccines. However, success of rotavirus immunization mostly depends on incorporation of prevalent genotype(s) in the vaccine. Several researchers have reported that G1–G4 are the predominant genotypes circulating globally [Reference Degiuseppe9–Reference O'Ryan11]. The two rotavirus vaccines, a live human monovalent strain and a polyvalent human-bovine reassortant strain [Reference Taneja and Malik12, Reference Cowley13] include only the predominant four G types. However, circulation of unusual rotavirus strains have been reported repeatedly in different parts of the world, e.g. the prevalence of [P6]G9 strains in India [Reference Ramachandran14], [P6]G8 strains in Malawi [Reference Cunliffe15], and [P8]G5 strains in Brazil [Reference Alfieri16]. Moreover, lack of host restriction in rotaviruses facilitates the emergence of reassortant strains [Reference Heatona and Ciarlet17]. There have been several reports regarding human rotavirus genotypes commonly found in animals [Reference Desselberger, Iturriza-Goýmara and Gray18, Reference Griffin19]. In India, the northeastern (NE) states is an agrarian region where animal rearing is an essential part of farming. This provides ample opportunity for humans and animals to share a common habitat. Such conditions generate the scope for rotavirus co-circulation during mixed infections in young children and animals with diarrhoea. Moreover, in some settings, close contact of humans with farm animals could generate an opportunity for inter-species transmission and gene reassortment [Reference Gouvea20, Reference Laird21]. Genetic diversity between individual rotavirus genotypes is now suspected to be more extensive than previously anticipated. This may potentially alter the immune response to vaccination. Several studies have evidenced the predominance of subtypes or genetic lineages of an individual genotype. Jin et al. [Reference Jin22] identified four distinct lineages of G1 rotavirus (G1–1, G1–2, G1–3, G1–4) and Piec & Palombo [Reference Piec and Palombo23] reported the prevalence of genetic subtypes of G2 and G4 rotavirus. These genetic variants among the individual types have the potential to develop genetically distinct strains with similar genotypes. From this perspective, it is important to study the genetic diversity occurring in rotaviruses globally as well as in NE India. Previously, use of a restriction fragment length polymorphism (RFLP) assay by Gouvea et al. [Reference Gouvea24] suggested successful examination of the degree of genetic diversity in rotaviruses. In the present study, the RFLP technique was adopted to identify and differentiate G types of NE Indian rotavirus isolates together with global isolates of known genotype from three different host species, i.e. human, porcine and bovine, sharing similar and diverse settings. This investigation will prove valuable in establishing genetic relationships between the Indian and global isolates from this region for the first time. Moreover, the comparison of strains from human and animal origins may further provide insights into the inter-species transmission of this virus.
METHODS
Virus samples
A total of 200 randomly selected faecal specimens were collected from February 2012 until November 2013. The faecal specimens were of children (n = 50), piglets (n = 80) and calves (n = 70) collected from hospitals and field outbreaks in Assam (26·1400° N, 91·7700° E). The local environment of the children, calves and piglets was also explored for collection of stool specimens. The samples were screened for rotavirus using a commercial monoclonal antibody-based enzyme immunoassay (mAb-EIA) (Premier Rotaclone; Meridian Bioscience Inc., USA) for detection of VP6 antigen. Samples were also processed for RNA–PAGE (polyacrylamide gel electrophoresis) for detection of viral nucleic acid as described previously [Reference Sharma25].
RNA purification and reverse transcriptase–polymerase chain reaction (RT–PCR)
Double-stranded rotavirus RNA genome was extracted using a QIAamp Viral RNA Mini kit (Qiagen, The Netherlands) according to the manufacturer's protocol. Full-length (1062-bp) VP7 gene segments encoding the major neutralization protein, VP7, was amplified by one-step RT–PCR with sets of generic primers, Beg9 and End9 [Reference Gouvea26]. The amplified products were resolved by conventional agarose gel (1%) electrophoresis and were visualized after ethidium bromide staining (0·1 mg/ml).
RFLP analysis
The amplified cDNA product was analysed by individual and sequential direct digestion with 1 U VspI (AT|TAAT), HaeIII (GG|CC) and NlaIV (GGN^NCC) restriction enzymes. The enzymes were selected since they produced distinct digestion patterns on the basis of published sequence data by empirical investigation in Sequence Manipulation Suite (SMS), v. 2 (www.bioinformatics.org/sms/). Following 2 h of incubation at 37 °C, the digested products were analysed by conventional agarose gel (2%) electrophoresis.
Database rotavirus sequences
Thirty-eight rotavirus VP7 gene sequences with known G types (G1–G4, G6, G8–G10, G12) were retrieved from the GenBank database, representing strains from various geographical locations (Supplementary Table S1). The sequences were analysed for the corresponding VspI, HaeII and NlaIV restriction sites using in silico restriction digestion in SMS v. 2 and NebCutter v. 2·0 (New England Biolabs Inc., USA). The fragment arrays, thus generated were correlated with the banding patterns of study samples (previously noted) and the characteristic G types of Indian isolates were determined.
RESULTS
Rotavirus was detected by mAb-EIA in 42/200 samples, accounting for 21% of samples tested positive. RNA–PAGE showed a characteristic 4:2:3:2 pattern, thus confirming all isolates as group A rotavirus. Of these, 16% (n = 32/200) of samples were confirmed positive by VP7 RT–PCR. Restriction fragments were successfully developed for 29 samples with each of the restriction enzymes, i.e. VspI, HaeIII and NlaIV. The restriction digestion patterns of the studied samples are presented in Figure 1. The different band patterns produced by each of the enzymes, VspI, HaeIII and NlaIV are defined as profiles V, H and N, respectively. The properties associated in respect of each designated profile along with the concordant restriction pattern produced by database sequences of known genotype with the studied samples are presented in Table 1.
A total of 29 studied samples were assessed with the enzymes. VspI, HaeIII and NlaIV digested the samples in 7, 6 and 6 patterns, respectively. Database VP7 sequences of known genotype were assessed with the same enzymes. Representative database sequences were shown from a set of sequences of the same genotype showing consensus restriction pattern with the studied samples.
a Each of the restriction enzymes generated profiles for total sets (n = 29) of Indian isolates used in the study.
b Represents the band size in agarose gel.
c Numbers in parentheses indicate number of isolates associated with the corresponding profile.
d Represents the host origin along with sample code.
e Represents the fragment size generated of global strains empirically with corresponding enzymes used for Indian strains.
f Restriction pattern of known G types of global strain compared to enzyme profiles of Indian strains.
VspI enzyme digested profile
Seven restriction profiles were noted after VP7 gene digestion of 29 Indian isolates with VspI and were designated as V profiles (V1–V7) (Table 1). Briefly, 2/29 rotavirus isolates, from bovine (IA-212) and porcine (IA-18) samples showed the V4 restriction profile found associated with genotype G6 (Table 1). Eighty-six percent (n = 25/29) of isolates were observed to be more diverse, generating a number of V profiles including profiles V3 (Fig. 1 a, lanes 4, 7, 8) and V6 (Fig. 1 a, lanes 13, 14) for G1 viruses, and profiles V1 (Fig. 1 a, lane 1), V2 (Fig. 1 a, lanes 2, 3–6, 9, 12) and V7 (Fig. 1 a, lane 15; Fig. 1 b; lanes 16–24, 26, 27, 29) for G8 viruses. Of these, 60% (n = 3/5, V3 profile) of G1 viruses and 95% (n = 19/20, V2 and V7 profile) of G8 viruses demonstrated multiple host association between each of the three hosts, i.e. human, bovine and porcine by sharing similar restriction profiles (Table 1). Profile V5 demonstrated a VspI non-cleavage site for two bovine (IA-71, IA-98) isolates; however, empirical analysis showed a consistent profile for database sequences (data not shown).
HaeIII enzyme digested profile
Six restriction profiles were obtained after digestion of VP7 cDNAs of 29 Indian isolates with HaeIII (profiles H1–H6; Table 1). Each strain revealed consistent and in some cases, additional cleavage or no cleavage patterns representing variation in the enzyme restriction sites. A brief analysis of Indian G2 viruses after gel electrophoresis demonstrated two bovine (IA-68, IA-132) isolates were associated with a single enzyme profile, H6 (Fig. 1 c, lanes 12, 13). A single porcine isolate, IA-21 designated as G6 genotype demonstrated a unique H profile, H5 (Fig. 1 c, lane 8) not previously recognized in any of our isolates. HaeIII digestion also suggested diversity of G1 viruses through association of 83% (n = 24/29) of isolates in three different H profiles, i.e. H1 (Fig. 1 c, lanes 1, 2, 4, 14, 15; Fig. 1 d, lanes 16, 17, 19–22, 24–29), H2 (Fig. 1 c, lanes 3, 7, 9–11; Fig. 1 d, lane 23) and H3 (Fig 1 c, lane 5). Of these, 96% (n = 23/24) of isolates demonstrated multiple host association by exhibiting similar restriction band patterns (Table 1, H1 and H2 profiles). The majority of isolates were from diverse host origin that were part of mixed infections (21%, n = 6/29) as revealed by the study and demonstrated a single H profile, H2, which suggested that co-infecting G type strains were dominant (Table 1). However, isolates IA-224 and IA-92 of human and porcine origin, respectively, demonstrated an HaeIII non-cleavage site (Fig. 1 c, lane 6; Fig. 1 d, lane 18), although in-silico analysis barely showed the pattern for database sequences (data not shown).
NlaIV enzyme digested profile
The six N profiles (N1–N6) generated by NlaIV were informative (Table 1). Similarly to the VspI and HaeIII enzyme profiles, NlaIV also demonstrated G1 as the most diverse genotype, associated with 66% (n = 9/29) of isolates with a number of N profiles including profile N1 (Fig. 1 e, lanes 1, 8, 9), N2 (Fig. 1 e, lane 2) and N3 (Fig. 1 e, lanes 3–5; Fig. 1 f, lanes 16–25, 27, 29). The majority (79%, n = 23/29) of the viruses showed multiple host alliance sharing an identical band pattern (Table 1, N1, N2 and N4 profiles). Isolates (17%, n = 5/29) that were a part of mixed infections, again had a single N profile, N4, indicative of one of the strains responsible for co-infection (Table 1, N4 profile). Three isolates [IA-122, IA-109 (porcine), and IA-212 (bovine)] designated as profile N5 could not be classified into any known G type due to demonstration of unique and additional restriction sites (Fig. 1 e, lanes 7, 15; Fig. 1 f, lane 26). Isolates IA-18 (bovine) and IA-212 (porcine), were associated with profile N6; however, they represented a NlaIV non-cleavage site (Fig. 1 e, lane 10; Fig. 1 f, lane 28).
Computational analysis of 38 randomly selected global sequences from GenBank database allowed the construction of the corresponding V, H and N profiles. An additional 12 V profiles were noted after VspI digestion, and these were designated V8–V19. The profiles were distinct from the associated Indian strain V profiles. Sixteen distinct H profiles (profiles H7–H22) were related to the database strains which were obtained after computational digestion of the retrieved global strains by HaeIII and similarly, 10 additional N profiles (N7–N16) were obtained from NlaIV digestion, producing a total of 16 profiles for the complete collection (see Supplementary Table S1).
Combined RFLP profile
The enzyme profiles whose G serotypes were proposed previously were combined together to produce a characteristic RFLP pattern. The 29 Indian rotavirus VP7 genes were classified into a total of 20 RFLP patterns, which are represented in Figure 2. The majority of the RFLP combinations were unique, except for two patterns, i.e. V7H1N3 and V7H1N5 which were shared between 38% (n = 11) of the strains. The V7H1N3 pattern was shared commonly among each host of diverse locality, i.e. out of nine isolates sharing the single pattern, four were human, three were porcine, and one was of bovine host origin. However, pattern V7H1N5, was shared between two bovine isolates from a single epidemic. Most of the RFLP patterns (40%, n = 8) were associated with G1/G8 virus infections; however, two patterns demonstrated G1/G6 G type association. The combined RFLP could also demonstrate a considerably higher number (34%, n = 10) of mixed infections by G8/G1 + G9 types throughout the study period. Few of the isolates that were part of mixed infections corresponded to multiple RFLP patterns and were suggested as variants by polyacrylamide gel electrophoresis (data not shown).
The evaluation of combined RFLP profiles of the global strains showed marked variations in the Indian isolates, and none of the patterns were shared by the Indian and global strains (see Supplementary Table S1).
DISCUSSION
An RFLP assay was designed for identification and differentiation of 29 group A rotavirus G types from three different hosts, i.e. human, porcine and bovine, sharing similar or diverse settings. Examination of VP7 restriction profiles obtained after digestion with VspI, HaeIII and NlaIV revealed several interesting features of rotavirus diversity in NE India. The banding patterns showed isolates demonstrating a single enzyme profile, while others had a combination of enzyme profiles or a unique RFLP pattern. The enzyme profiles generated suggested that when a single G type infection occurred there was an obvious similarity between the enzyme profiles. Here, the study of restriction profiles of the isolates demonstrated relatedness of rotavirus VP7 genes belonging to the same G type. Major G types G1, G2, G6 and G8 along with a few mixed G types signify rotavirus diversity within the studied population. Detection of G1, G2 and G9 rotavirus G types from India by Kang et al. [Reference Kang8] supports the present finding which highlights a substantial rotavirus disease burden during these periods of time. Inter-species transmission or sharing of genotypes within different hosts is reported regularly by various workers [Reference Mukherjee27, Reference Steyer28]. Our finding is consistent with those reports where human rotavirus G type G1 is circulating significantly in the human population as well as in porcine and bovine populations. Although the G1 genotype is highly prevalent in human populations, infection with human G type in bovine and porcine populations suggests the sharing of specific rotavirus G types in these species [Reference Kobayashi and Mendez-Vilas10]. G6 and G8, being the major bovine G types [Reference Cui7], were observed in infections in both human and porcine neonates. Serotype G9, where pigs are suspected as a potential host reservoir, was also associated with infections in humans [Reference Matthijnssens, Palombo and Kirkwood29]. Hence, the co-circulation or sharing of intrinsic G types/RFLP patterns in diverse host ranges suggests the occurrence of inter-species transmission particularly during mixed infections and more generally in the setting of close contact between humans and farm animals. This phenomenon can be evidenced from the study area where closer contacts between animals and their human handlers are recognized. Moreover, backyard pig/cattle rearing is a customary activity of NE India where animals are housed indoors in group-housing or straw-lined sheds or pens, particularly for pigs which allows easy contact with waste matter for the human handlers and vice versa. The socioeconomic status of the studied population was notably low with poor hygiene practised. In such circumstances, the faecal–oral route is the best route for viral transmission within the associated hosts. Thus, the animals along with their associated human handlers, are both at risk of contracting infection, which is maintained in the environment. Reassortment of the rotavirus genome occurring after co-infection of a host, has been shown to be an important mechanism to generate diversity on many occasions [Reference Martella30] and, moreover, other less important mechanisms such as inter- or intragenic recombination are believed to occur less frequently [Reference Cao, Barro and Hoshino31]. The affinity of the Indian rotavirus isolates towards sharing of G types could significantly boost evolution of the viruses and subsequent emergence of atypical or novel strains.
Comparison of RFLP data for Indian and global isolates suggests the Indian rotavirus population is distinguishable from global strains as none of their associated RFLP patterns were shared. The global isolates produced a consistent profile through empirical analysis but the Indian isolates exhibited a greater diversity by providing additional or no cleavage sites. It could be assumed that there were some point mutations in the gene during replication, and existing restriction enzyme sites disappeared or new sites were generated. Restriction enzymes were also found to have the potential to identify the presence of mixed infections in the circulating rotavirus strains by Halloran et al. [Reference O'Halloran32]. The combined analysis of the RFLP patterns produced considerable numbers of mixed rotavirus infection within the hosts. Thus, it may be reasonable to suggest that one of the co-existing G types was dominant towards infection.
Recent studies have reported the alliance of diverse strains from different host origins, in a few cases this may be originated by the two major mechanisms that are believed to be responsible for the production of genetic/genetic variants, i.e. nucleotide substitution and gene reassortment [Reference Suzuki, Gojobori and Nakagomi33, Reference Iturriza-Gomara34]. With the series of rotavirus strains reported from India [Reference Wani35–Reference Dubal37], it appears likely that unique animal–human mixing patterns add to the potential inter-species transmission of rotaviruses. In this study, direct evidence of inter-species spread has been obtained by characterizing human and animal isolates that are linked by similar VP7 gene RFLP patterns along with epidemic, time and place. Such circumstances may lead to the realization that the diversity of co-circulating rotavirus is much greater than previously believed, and has shown that introduction of novel rotavirus types in the population are likely to be frequent events in NE India. Moreover, the detection of porcine and bovine G types along with differentiation of rotavirus types associated with three different hosts are reported for the first time from NE India.
Since rapid evolution of these viruses by generation of reassortment in multiple infections is evident, molecular epidemiological surveillance of the rotavirus types co-circulating in the population is indispensable. Therefore, the present study could differentiate rotavirus G types by RFLP into major G genotypes (G1, G2, G6, G8) along with a few mixed types (G1 + G9). Positive inter-species transmission and significant variations in the genotypes were also anticipated but the extent of such mechanisms are ambiguous and will need to be established by analysing more genes and their respective sequence data. Thus, RFLP assays could be a valuable tool for monitoring the emergence of uncommon strains and assessing the genetic diversity.
SUPPLEMENTARY MATERIAL
For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S0950268814003343.
ACKNOWLEDGEMENTS
The authors gratefully acknowledge Department of Microbiology, College of Veterinary Science, Assam Agricultural University for providing cell culture adapted standard reference strains of rotavirus. The authors are also grateful to the patients and animal handlers for their willingness to participate and to provide information about themselves and their animals. The authors also thank Dr Maloyjo Joyraj Bhattacharjee, Department of Biotechnology, Assam University for the detailed and critical editing of the paper.
This research received no specific grant from any funding agency, commercial or not-for-profit sectors.
DECLARATION OF INTEREST
None.