Next Article in Journal
Pathogens Manipulating Tick Behavior—Through a Glass, Darkly
Next Article in Special Issue
Molecular Epidemiology of Rotavirus A Strains Pre- and Post-Vaccine (Rotarix®) Introduction in Mozambique, 2012–2019: Emergence of Genotypes G3P[4] and G3P[8]
Previous Article in Journal
Comparative Pathology of Domestic Pigs and Wild Boar Infected with the Moderately Virulent African Swine Fever Virus Strain “Estonia 2014”
Previous Article in Special Issue
Epidemiological Trends of Five Common Diarrhea-Associated Enteric Viruses Pre- and Post-Rotavirus Vaccine Introduction in Coastal Kenya
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Characterisation of a Rare Reassortant Porcine-Like G5P[6] Rotavirus Strain Detected in an Unvaccinated Child in Kasama, Zambia

1
Next Generation Sequencing Unit, Division of Virology, Faculty of Health Sciences, University of the Free State, Bloemfontein 9300, South Africa
2
Virology Laboratory, Department of Pathology & Microbiology, University Teaching Hospital, Adult and Emergency Hospital, Lusaka 10101, Zambia
3
Department of Paediatrics & Child Health, School of Medicine, University of Zambia, Ridgeway, Lusaka RW50000, Zambia
4
World Health Organization, Regional Office for Africa, Brazzaville P.O. Box 06, Congo
5
Diarrhoeal Pathogens Research Unit, Faculty of Health Sciences, Sefako Makgatho Health Sciences University, Medunsa, Pretoria 0204, South Africa
6
South African Medical Research Council, 1 Soutpansberg Road, Pretoria 0001, South Africa
*
Author to whom correspondence should be addressed.
Submission received: 24 July 2020 / Revised: 14 August 2020 / Accepted: 14 August 2020 / Published: 17 August 2020
(This article belongs to the Special Issue Rotaviruses and Rotavirus Vaccines)

Abstract

:
A human-porcine reassortant strain, RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6], was identified in a sample collected in 2014 from an unvaccinated 12 month old male hospitalised for gastroenteritis in Zambia. We sequenced and characterised the complete genome of this strain which presented the constellation: G5-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1. The genotype A8 is often observed in porcine strains. Phylogenetic analyses showed that VP6, VP7, NSP2, NSP4, and NSP5 genes were closely related to cognate gene sequences of porcine strains (e.g., RVA/Pig-wt/CHN/DZ-2/2013/G5P[X] for VP7) from the NCBI database, while VP1, VP3, VP4, and NSP3 were closely related to porcine-like human strains (e.g., RVA/Human-wt/CHN/E931/2008/G4P[6] for VP1, and VP3). On the other hand, the origin of the VP2 was not clear from our analyses, as it was not only close to both porcine (e.g., RVA/Pig-tc/CHN/SWU-1C/2018/G9P[13]) and porcine-like human strains (e.g., RVA/Human-wt/LKA/R1207/2009/G4P[6]) but also to three human strains (e.g., RVA/Human-wt/USA/1476/1974/G1P[8]). The VP7 gene was located in lineage II that comprised only porcine strains, which suggests the occurrence of independent porcine-to-human reassortment events. The study strain may have collectively been derived through interspecies transmission, or through reassortment event(s) involving strains of porcine and porcine-like human origin. The results of this study underline the importance of whole-genome characterisation of rotavirus strains and provide insights into interspecies transmissions from porcine to humans.

1. Introduction

Group A rotaviruses (RVA), of the family Reoviridae, are the number one viral pathogens causing severe diarrhoea in children below five years of age [1]. In 2016, an estimated 128,000 deaths in children below five years were due to RVA infections, 90% of which occurred in developing countries [2,3]. Similarly, RVA are the primary cause of acute gastroenteritis in new-born piglets [4].
Rotaviruses have a distinctive morphology which comprises a nonenveloped, three-layered icosahedral protein shell. The rotavirus genome within the protein shell comprises 11 segments of double-stranded (dsRNA) that encode six structural viral proteins (VP1 to VP4, VP6, and VP7) and five or six nonstructural proteins (NSP1 to NSP5/6) [1]. A binary classification system is used to distinguish RVA based on the antigenic properties of the outer shell proteins, VP7 and VP4, that determine the G-genotype and P-genotype, respectively [1]. Furthermore, RVA can be separated into two main genogroups and one minor genogroup according to a whole-genome classification system, whereby a specific genotype is assigned to the 11 gene segments. These genogroups represent the genotype constellations that are present in most human strains globally [5,6]. Genogroup 1 (Wa-like) bears the constellation I1-R1-C1-M1-A1-N1-T1-E1-H1 and is often associated with the G genotypes G1, G3, G4, G9, and G12 and P genotype P[8]. Genogroup 2 (DS-1-like) includes G2P[4] strains and bears the constellation I2-R2-C2-M2-A2-N2-T2-E2-H2. Lastly, the minor genogroup 3 (AU-1-like) bears the I3-R3-C3-M3-A3-N3-T3-E3-H3 constellation and includes G3P[9] strains [7]. As of 5th May 2020, the Rotavirus Classification Working Group had identified at least 36 G, 51 P, 26 I, 22 R, 20 C, 20 M, 31 A, 22 N, 22 T, 27 E, and 22 H genotypes [8]. The whole-genome classification system has made it possible to analyse and understand the origin of various strains, interspecies transmission, and animal–human reassortment events [9]. Human Wa-like strains and porcine rotavirus strains share a common origin, whereas DS-1-like and AU-1-like strains have a common origin with bovine and feline strains, respectively [5].
In humans, G1-G4, G9, and G12 along with P[4], P[6], and P[8] are the most frequently detected, globally [10,11,12,13]. On the contrary, in porcine, predominant genotypes are G3-G5, G9, and G11 along with P[6], P[7], and P[13] [4,14]. Porcine rotaviruses bear the constellation I5-R1-C1-M1-A8-N1-T1/T7-E1-H1 [5,15,16,17,18,19,20]. While human Wa-like RVA differ from porcine rotaviruses in some gene segments (VP4, VP6, VP7, and NSP1), they both appear to have genotype 1 in the VP1, VP2, VP3, NSP2, NSP3, NSP4, and NSP5 gene segments. Hence, the suggestion that human Wa-like and porcine RVAs have arisen from a common ancestor [5].
The findings that show animals can serve as potential reservoirs for genetically diverse rotavirus strains that can be passed on to humans have elicited a large amount of interest and topics for further research [21]. Several novel and rare animal-like or animal–human reassortant rotavirus strains have been identified globally [22,23,24,25,26,27,28]. The detection of animal strains in humans is presumed to be as a result of zoonotic transmission, along with reassortment, which contributes to the diversity of circulating RVA [4,29,30]. Inter- and intragenogroup reassortment may occur when multiple RVA simultaneously infect a host. This is attributed to the segmented nature of the rotavirus genome [1,31]. It is, therefore, necessary to continuously carry out the monitoring of animal RVA and the role they play in contributing to the diversity of circulating RVA in humans.
The G5, one of the most common porcine genotypes, has sporadically been identified in human populations in Brazil (G5P[X]), Cameroon (G5P[7] and G5P[8]), Argentina (G5P[8]), and the United Kingdom(G5P[X]) [32,33,34,35,36]. The P[6] is presumed to be of porcine origin. They have also been identified in human populations [37,38,39,40]. The first human G5P[6] strain, LL36755, was detected in a child who had acute gastroenteritis in China in 2007 [41]. Other G5P[6] strains were detected in Vietnam, Taiwan, Bulgaria, Japan, and Thailand [37,42,43,44,45]. To date, the whole-genome of only two human G5P[6] strains—Bulgarian BG620 (nt sequences unavailable in the DDBJ, EMBL, and GenBank data libraries as of 13 August 2020) and Japanese Ryukyu-1120 (full open reading frame, available in GenBank)—have been analysed [45,46].
Diarrhoea is a burden for the Zambian healthcare system, with about 33% of the extreme cases being attributable to RVA [47,48,49]. In an attempt to generate disease burden attributable to rotavirus diarrhoea in children, the Zambian Ministry of Health, with support from WHO, launched rotavirus surveillance at the University Teaching Hospital (UTH) in 2006 [50,51]. Surveillance data generated provided evidence of the burden of rotavirus diarrhoea that supported the introduction of the rotavirus vaccine, Rotarix®, as a pilot project in Lusaka, Zambia in 2012, and was later rolled out nationwide in November 2013 [50]. According to the estimates reported by the World Health Organization (WHO) and the United Nations International Children’s Emergency Fund (WHO/UNICEF), rotavirus vaccine coverage in Zambia has been consistently high for the last six years, increasing from 73% in 2014 to 90% in 2019 [52]. Over this period, a sustained and significant reduction in rotavirus-associated hospitalisations and mortality was observed in children under 5 years [51].
The African Rotavirus Surveillance Network, coordinated by the World Health Organization Regional Office for Africa (WHO/AFRO), is actively monitoring the diversity and distribution of RVA genotypes in children hospitalised with acute diarrhoea [53]. Initially, the network was established with four countries in 2006, and expanded to 29 countries by the end of 2016 [54,55]. The Diarrhoeal Pathogens Research Unit at Sefako Makgatho University in Pretoria (South Africa) and the Noguchi Memorial Institute for Medical Research in Accra (Ghana) are the two WHO Rotavirus Regional Reference Laboratories (RRLs) for the network that conducts monitoring of rotavirus epidemiology in Africa [55]. The WHO/AFRO is currently supporting the University of the Free State-Next Generation Sequencing (UFS-NGS) unit to undertake rotavirus surveillance of rotavirus strains that circulated in Zambia between 2013 and 2016 at the whole-genome level. A G5P[6] strain, UFS-NGS-MRC-DPRU4723, was identified among these strains and was analysed so as to elucidate its origin and evolution. The sample was collected in 2014 from an unvaccinated 12 month old male hospitalised for gastroenteritis at Arthur Davison Children’s Hospital in Ndola, Zambia.

2. Results

2.1. Nucleotide Sequencing and Identity of the Strain

Illumina® MiSeq sequencing exhibited a phred score of Q30 and collectively yielded 98.8 Mbs of data for this specific sample. The whole genome of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] was 18272 bps in size. The length and ORF of the 11 gene segments as determined by nucleotide sequencing are shown in Table 1. A BLASTn search was performed, and it appeared to exhibit maximum sequence identities of 95.7%–98.0% with porcine and human porcine-like strains (Table 1). Based on the whole genome classification system, RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] exhibited a G5-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1 genotype constellation (Table 2). The genetic constellation of the study strain was compared to those of other G5 and non-G5 strains retrieved from the GenBank (Table 2).

2.2. Sequence and Phylogenetic Analysis

To investigate the potential origin of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6], phylogenetic trees were constructed for each of the 11 gene segments along with cognate gene sequences of RVA strains obtained from the GenBank.

2.2.1. Sequence and Phylogenetic Analysis of the VP7 Gene

Phylogenetically, there are three known VP7 G5 lineages (I-III) [63]. The VP7 genes of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered into lineage II, which consisted only of porcine G5 strains from mainly Asia and the Americas (Figure 1). The VP7 gene showed the highest nucleotide (nt) and amino acid (aa) identities with the Chinese porcine strains RVA/Pig-wt/CHN/DZ-2/2013/G5P[X] nt (aa), 98.6% (99.0%), and RVA/Pig-wt/CHN/JN-2/2014/G5P[X] 98.5% (99.0%) and was distantly related to the strains within lineage III with lower sequence identities (nt, 83.4%–86.5%; aa, 90.4%–94.5%) (Figure 1; Supplementary data 1). Overall, strains within lineage II exhibited sequence identities that were in the range nt, 89.6%–98.6%; aa, 92.4%–99.0% (Supplementary data 1).
The comparison of the amino acid sequence of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] to reference G5 strains e.g., RVA/Pig-wt/THA/CMP-001-12/2012/G5P[13] (lineage I), RVA/Pig-wt/BRA/ROTA24/2013/G5P [6] (lineage II) and RVA/Human-wt/JPN/Ryukyu-1120/2011/G5P[6] (lineage III) within each of the three lineages revealed a high identity (range 90.0%–94.9% (Supplementary data 1; Supplementary data 2a). Numerous substitutions were identified in the nine VP7 variable regions, VR-1 to VR-9 [64]: VR-1 (I9V and I19V), VR-2 (V27T and V29T), VR-3 M/F39L, I40V, V41I, L/I43V, I/L/V47F, R49K, and A50T), VR-4 (K/A65T, V/M68A, M/A72T, and M/Q75T), VR-5/antigenic site A (N/S/D/T96A), VR-6 (I129V and D130E), VR-7/antigenic site B (N145D and A/V/E146G), VR-8/antigenic site C (L/S208T, A210T, T/V212I, S/A213I, I/M217T, V218I, and S220N), and VR-9/antigenic site F (A/M241T and S242N).

2.2.2. Sequence and Phylogenetic Analysis of the VP4 Gene

The VP4 gene of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] was phylogenetically compared to the already established five lineages (I-V) of genotype P[6] [65] (Figure 2). The P[6] gene of the study strain clustered into lineage V, which consisted of porcine and putative human porcine-like strains detected in parts of Europe and one African strain. A similarity analysis of the P[6] gene of the study strain with strains obtained from GenBank showed that the Zambian G5P[6] exhibited the highest sequence identity of 98.1% (98.3%) with a porcine-like human strain RVA/Human-wt/COD/KisB332/2008/G4P[6] from the Democratic Republic of Congo (Supplementary data 1). All the African strains clustered into a separate lineage, lineage I, with sequence identities of 85.7%–86.8% (92.5%–93.9%) (Supplementary data 1).
The deduced amino acid sequences of the VP4 gene of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] along with the reference P[6] strain from each of the five lineages was compared (Supplementary data 2b). The reference strains shared high amino acid identities ranging from 91.0% to 98.3% (Supplementary data 1). Several amino acid changes were identified throughout the VP4 protein, and most of the substitutions were concentrated in the hypervariable region (amino acid 71-208) which houses the VR-3 (92%–192) and includes a neutralization site at amino acid 135 [66,67]. Several amino acid substitutions were observed among the P[6] lineage I strains [65] at the VR-3 (L105I, V108I and T134S) and VR-8 (D602N) variable regions. Other amino acid substitutions were identified among the P[6] lineages at VR-1 (S30N), VR-2 (I61V), VR-3 (V112I, N114S, V130I, H182N and T189S), VR-4 (I280V), and VR-9 (E698K). The potential trypsin cleavage sites at residues 241 and 247 [68] were highly conserved in all the strains with three substitutions at positions 242 (I to V), 243 (A to T), and 244 (H to Y).

2.2.3. Phylogenetic Analysis of the VP6 Gene

The VP6 gene of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered closely with divergent African porcine strains from Uganda (RVA/Pig-wt/UGA/BUW-14-A003/2014/G3P[13], RVA/Pig-wt/UGA/KYE-14-A048/2014/G3P[13], and RVA/Pig-wt/UGA/KYE-14-A047/2014/G3P[13]) and a human porcine-like strain from the Democratic Republic of Congo (RVA/Human-wt/COD/KisB332/2008/G4P[6]) which displayed nt(aa) sequence identities ranging from 98.6% to 98.9% (98.9%–99.7%) (Figure 3, Supplementary data 1). Porcine-like Asian strains such as RVA/Human-wt/CHN/GX54/2010/G4P[6] and RVA/Human-wt/CHN/E931/2008/G4P[6] clustered separately, displaying identities of 88.7%–90.2% (97.5%–98.7%) (Supplementary data 1).

2.2.4. Phylogenetic Analysis of VP1 Gene

The VP1 gene of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered only with porcine and porcine-like human strains from Asia (China and Vietnam) (Supplementary data 3a). The VP1 gene exhibited a maximum nt (aa) sequence identity of 96.8% (98.9%) with the Chinese human porcine-like reassortant strains RVA/Human-wt/CHN/GX82/2010/G4P[6], RVA/Human-wt/CHN/GX78/2010/G4P[6], RVA/Human-wt/CHN/GX77/2010/G4P[6], and RVA/Human-wt/CHN/GX54/2010/G4P[6] (Supplementary data 1). Overall, the Asian strains within the cluster showed sequence identities of 94.1%–96.8% (97.9%–98.9%). Human non-porcine African strains clustered separately, with lower identities of 88.2%–88.8% (96.3%–97.3%) (Supplementary data 1).

2.2.5. Phylogenetic Analysis of VP2 Gene

The VP2 gene of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] fell into a distinct cluster predominantly composed of porcine and porcine-like human strains from Asia (China, India, Vietnam, South Korea, and Sri Lanka) (Supplementary data 3b). The VP2 gene of the study strain showed a maximum nt (aa) sequence identity of 96.6% (90.9%) with a Sri Lankan porcine-like human strain RVA/Human-wt/LKA/R1207/2009/G4P[6] (Supplementary data 1).

2.2.6. Phylogenetic Analysis of VP3 Gene

The VP3 gene of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered in a lineage composed mainly of Asian (Asia and Thailand) porcine and porcine-like human strains (Supplementary data 3c), and exhibited the highest nt (aa) sequence identity with the Chinese porcine-like human strains—RVA/Human-wt/CHN/R946/2006/G3P[6], 95.8% (97.8%) and RVA/Human-wt/CHN/E931/2008/G4P[6], 95.7% (98.0%) (Supplementary data 1). The overall similarities of the Asian strains within the lineage ranged from 84.8% to 95.8% (92.7%–97.8%) (Supplementary data 1). Non-porcine African strains clustered separately and showed lower sequence identities of 84.1%–84.5% (92.1%–92.7%) (Supplementary data 1).

2.2.7. Phylogenetic Analysis of NSP1 Gene

The NSP1 gene of strain RVA/Human-wt/ZMB/UFS-NSG-MRC-DPRU4723/2014/G5P[6] was assigned to a porcine genotype A8 and clustered among Asian (Vietnam, China, and Bangladesh) porcine and porcine-like human strains and an African (Ghana) porcine strain (Supplementary data 3d). The NSP1 gene of the study strain was closest to strain RVA/Human-tc/VNM/NT0042/2007/G4P[6] displaying a nt(aa) sequence identity of 98.2% (97.9%) (Supplementary data 1). The porcine and porcine-like human strains from Europe and the Americas clustered separately showing sequence identities of 84.2%–85.9% (85.4%–88.2%) and 84.1%–85.9% (83.7%–88.3%), respectively (Supplementary data 1).

2.2.8. Phylogenetic Analysis of NSP2 Gene

The NSP2 gene of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered with Asian and European porcine and porcine-like human strains (Supplementary data 3e). The Nt(aa) similarity analysis showed that the NSP2 gene of the study strain was most similar to the Chinese porcine strains RVA/Pig-wt/CHN/YN/2012/GXP[X] and RVA/Pig-tc/CHN/SCMY-A3/2017/G9P[23]—96.8% (97.8%) (Supplementary data 1). Two African porcine strains, RVA/Pig-wt/ZAF/MRC-DPRU1487/2007/G3G5P[23] and RVA/Pig-wt/ZAF/MRC-DPRU1557/2008/G4G5P[23], were seen to cluster within the same lineage with sequence identities of 93.6%–93.7% (97.5%–97.8%) (Supplementary data 1).

2.2.9. Phylogenetic Analysis of NSP3 Gene

The NSP3 gene of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered closely with porcine and porcine-like human strains mainly from Asia (Thailand and Vietnam) and exhibited a maximum nt(aa) sequence identities of 96.5%–97.0% (98.4%–98.7%) with the strains RVA/Human-wt/VNM/30378/2009/G26P[19], RVA/Pig-wt/VNM/12070-4/2012/GXP[X], RVA/Human-wt/VNM/NT0205/2007/G4P[6], and RVA/Human-wt/VNM/NT0621/2008/G4P[6] (Supplementary data 1; Supplementary data 3f).

2.2.10. Phylogenetic Analysis of NSP4 Gene

The NSP4 gene of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered with porcine and porcine-like human strains identified in Asia (China and Vietnam) and a porcine-like human strain from the Americas (Brazil) (Supplementary data 3g). In this cluster, the closest strains to UFS-NGS-MRC-DPRU4723 were the wild pig strains (RVA/WildBoar-wt/CZE/P828/2015/G9P[23] and RVA/WildBoar-wt/CZE/P830/2015/G9P[23]) from the Czech Republic, with nt(aa) sequence identities of 97.5% (98.3%) (Supplementary data 1). The Asian strains within the cluster showed nt(aa) similarities of 96.2%–97.3% (97.7%–98.9%). Porcine and porcine-like human strains from the Americas clustered separately and exhibited identities of 87.2%–96.4% (94.3%–98.9%) (Supplementary data 1).

2.2.11. Phylogenetic Analysis of the NSP5 Gene

The NSP5 gene of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] clustered with porcine strains from Asia and showed the highest nt(aa) sequence identity of 98.6% (100%) with the porcine strains RVA/Pig-wt/CHN/TM-a/2009/G3P[8] and RVA/Pig-tc/CHN/TM-a-P20/2018/G9P[23] identified in China (Supplementary data 1; Supplementary data 3h). Overall, the porcine and porcine-like human strains from Asia and the Americas displayed nt(aa) identities of in the range 94.8%–98.6% (98.0%–100%) and 93.9%–96.1% (95.9%–99.0%), respectively (Supplementary data 1).

2.3. Reassortment Analysis

The concatenated whole genome alignment of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6], together with the Japanese G5P[6] strain and selected Chinese porcine-like human P[6] strains, was visualised (Figure 4). The whole genome of the Zambian G5P[6] strain demonstrated a relatively high degree of conservation with the Japanese G5P[6] strain and the two Chinese G4P[6] strains. With the exception of VP7 and VP4, the genome of the Chinese strain E931 exhibited the overall highest genomic conservation to the study strain. With the exception of VP7, VP3, and NSP1 genes, the Chinese strain GX54 shared a highly conserved genome with the study strain. The Japanese strain Ryukyu-1120 demonstrated a highly similar genome to the study strain for seven of the 11 genes, the exceptions being VP1, VP3, VP6, and VP7. The results of this analysis confirmed the genetic similarity between RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and Asian (Chinese) porcine-like human strains, hence suggesting that the Zambian G5P[6] strain may have been derived via reassortment events.

3. Discussion

The detection of genotype G5 in humans, which is typical for pigs, is possibly due to interspecies transmission [35,45]. In Zambia, as with many countries in Africa, humans and farm animals live in proximity. The interaction between humans and animals could be the primary cause for zoonotic transmission, which could result in genetic reassortments and perhaps other mechanisms of genetic diversity, ultimately leading to the introduction and spread of animal genotypes into human populations [69].
In this study, an analysis was conducted on a sample collected from a child admitted to a paediatric ward presenting with clinical symptoms (vomiting, diarrhoea, and fever) that are usually present during typical rotavirus infection. This raises the question whether such animal-derived strains are capable of mutating and effectively spreading within/across human populations as in the case of established typical Wa-like and DS-1-like genotype constellations, with the same magnitude of rotavirus disease severity. Furthermore, taking into consideration that the G5 and P[6] genotypes are not included in the currently available vaccines, the probability for such strains to have the potential to spread more swiftly from human to human may have implications for the effectiveness of current rotavirus vaccine candidates that are in use in African countries.
This study identified the complete genome of a reassortant porcine-like human strain, G5P[6], that showed the genotype constellation G5-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1, which is commonly found in porcine and porcine-like human rotavirus strains [19]. RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] was found to share the same constellation (I1-R1-C1-M1-A8-N1-T1-E1-H1) with the archival porcine strain, Gottfried, and porcine-like human strains—BG260, E931, and GX54 [5,46,56,58]. In addition, porcine strains 12R002, 12R005, and 12R006, as well as porcine-like human strains Ryukyu-1120, mani-97, 30378, rj24598, and BE2001 shared the same constellation with strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] with the exception of VP6 (I5 instead of I1) and NSP3 (T7 instead of T1 gene segments) [20,25,26,45,70].
A phylogenetic analysis of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] showed that this strain was a possible reassortant, as it was closely related to both porcine and porcine-like human strains, predominantly from Asia, than to typical human RVA strains. The VP6, VP7, NSP2, NSP4, and NSP5 segments of this strain showed a close similarity to porcine strains. Although the remaining gene segments (VP1, VP3, VP4, and NSP3) were closely related to human strains, all of these were porcine-like human strains [26,56,58,59,60,70]. With a genotype 1 (Wa-like) backbone, this finding is consistent with the hypothesis that human Wa-like strains and porcine strains have a common ancestor [5]. However, the origin of the VP2 gene of the study strain was not very definitive, as it was not only close to porcine and porcine-like human strains but also to three human strains (DC1476, DC582, and DC1127). Phylogenetically, the clusters of these three strains were shown to be distinctive from the genes of contemporary, wild-type human strains [71]. Notably, the VP7 gene of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] was located in lineage II, which comprised only porcine strains, hence implying the possibility of porcine-to-human interspecies transmission [63]. Phylogenetic analysis of porcine and human P[6] strains indicated that both porcine and human P[6] strains were present in P[6] lineages I, III, and V, hence showing that human P[6] strains might have separately emerged from at least three porcine-to-human transmissions [65]. This finding supports the Zambian G5P[6] strain, as the VP4 gene clustered and shared high nucleotide and amino acid identities with lineage V of P[6] porcine and porcine-like human strains. The NSP1 gene was most similar to porcine-like human strains. However, it was revealed to have the porcine genotype A8. Taking this together, it is likely that RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] originated by zoonotic transmission, coupled with reassortment events.
Several amino acid changes were identified in the nine variable regions when the VP7 gene of the study strain was compared to other G5 strains within each of the three lineages [64]. Additionally, the previously described conserved N-glycosylation site at residues 69–71 within the variable region 4 (VR-4) was found to be conserved in all the G5 strains used in this analysis [64,72]. Four major antigenic regions have been described for the VP7 protein in rotaviruses (A, B, C and F) [73,74]. Marked differences in the antigenic regions of RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] were seen when it was compared to other globally circulating G5 strains. Usually, antigenic regions A and C are said to be conserved within serotypes [75]. However, multiple substitutions were observed in these regions when comparing the Zambian G5 strain to other G5 strains globally.
The amino acid sequence for the VP4 gene was 775 amino acids long and displayed amino acid identity values ranging from 91.0% to 98.3% with the reference P[6] strains. Considering it has been established that strains with amino acid identities greater than 89% belong to the same P genotype [76], our findings show that RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] belongs to the genotype P[6]. The analysis of the amino acid sequences showed that the hypervariable region (amino acid 71-208) which houses the variable region 3 (VR-3) contained most of the substitutions. Furthermore, the potential trypsin cleavage sites [68] were conserved in all the P[6] strains. Several amino acid substitutions were observed among the lineage I P[6] strains. The presence of several amino acid changes in the VP4 gene of this strain compared to other circulating P[6] strains globally is in agreement with the hypothesis that the P[6] gene has been introduced to humans via independent reassortment events [40,65,77].
Rotaviruses are genetically diverse in nature and are host-species specific, suggesting that host species barriers and restrictions exist. However, rotaviruses of animal origin may cross the host species barrier and may acquire human rotavirus gene segments, which enables the viruses to efficiently spread across human populations [4]. In this regard, G5 rotavirus strains have sporadically been documented in Latin America, Asia, Europe, and Africa [33,34,35,36,37,41,45,46]. Porcine P[6] strains seem to pose a lesser species barrier to humans [20]. Even though the relationship between porcine and human rotaviruses has already been established [5], whole genome analysis in this study presented the possible occurrence of interspecies transmission and reassortment between human and porcine rotaviruses.

4. Materials and Methods

4.1. Ethics Statement

This is a subset of a major project which involved the whole genome characterisation of 133 specimens collected in Zambia from 2013 through 2016 as part of the surveillance supported by the WHO/AFRO (reference 2017/757922-0) in collaboration with the University of the Free State (UFS-NGS). Ethical clearance for the main project was obtained under ethics number HSREC130/2016(UFS-HSD2016/1082) from the Health Science Research Ethics Committee (HSREC), University of the Free State, Bloemfontein, South Africa. Furthermore, this specific study was approved by the HSREC under ethics number UFS-HSD2020/0277/2104.

4.2. Sample Collection

The sample was collected in 2014 from an unvaccinated 12 month old male at Arthur Davidson Children’s Hospital (ADCH) in Ndola, a rotavirus surveillance sentinel site. The child had travelled with parents from Kasama, a town in the Northern Province of Zambia which is approximately 760 km away from Ndola, Zambia. This child was admitted to a paediatric ward at ADCH, with gastroenteritis of four days duration and a history of fever. Frequency of vomiting and diarrhoea was three episodes and two episodes, respectively, in the previous 24 h. The level of dehydration was assessed as mild and the child received an oral rehydration solution and was discharged after a few days. The stool sample was screened using the enzyme immunoassay (EIA) technique for the presence of RVA antigen in the Virology laboratory in Lusaka. It was randomly picked and sent to the Diarrhoeal Pathogens Research Unit (DPRU), a World Health Organization Rotavirus Regional Reference Laboratory (WHO-RRL) in Pretoria, South Africa, as part of the WHO/AFRO annual rotavirus surveillance. Conventional genotyping was carried out at DPRU. Thereafter, the sample was shipped to the UFS-NGS unit for sequencing and whole-genome analysis.

4.3. Viral dsRNA Extraction

The viral double-stranded RNA (dsRNA) was extracted from human stool suspensions using a previously described method with modifications [78]. Approximately 100 mg stool was suspended in 200 µL phosphate-buffered saline (PBS) solution (Sigma-Aldrich®, St Louis, MO, United States). The faecal suspension was mixed with 900 µL TRI Reagent® LS (Molecular Research Centre, Cincinnati, OH, United States) and homogenized for five minutes. A 300 µL volume of chloroform (Sigma-Aldrich®, St Louis, MO, United States) was used to achieve phase separation, which was followed by centrifugation (Eppendorf microcentrifuge 5427 R, Germany) at 17,319× g for 20 min at 4 °C. The supernatant was precipitated using 700 µL ice-cold isopropanol (Sigma-Aldrich®, United States) and centrifuged (Eppendorf microcentrifuge 5427 R, Germany) at 17,319× g for 30 min at 4 °C. The supernatant was discarded, and the tubes were air-dried for 5 min, followed by the precipitation of single-stranded RNA (ssRNA) using 30 µL 8 M lithium chloride (Sigma, St Louis, MO, United States) at 4 °C for 16 h. The dsRNA was purified using the MinElute gel extraction kit (Qiagen, Hilden, Germany). RNA integrity was determined by electrophoresis on 1% TBE agarose gel stained with ethidium bromide (Sigma-Aldrich®, St Louis, MO, United States), which was visualised on a G: Box UV transilluminator (Syngene, Cambridge, United Kingdom).

4.4. cDNA Synthesis and Purification

cDNA synthesis was carried out using the Maxima H Minus Double-stranded cDNA kit (Thermo Fisher Scientific, Waltham, MA, United States) according to the manufacturer’s instructions with minor modifications captured at the UFS-NGS SOP, whereby the dsRNA was denatured at 95 °C for 5 min. First strand synthesis was carried out for two hours at 50 °C. Random hexamer primer was employed for cDNA synthesis. The cDNA was purified using the MSB® Spin PCRapace purification kit (Stratec, Invitek Molecular, Berlin, Germany).

4.5. DNA Library Preparation and Illumina® MiSeq Sequencing

DNA libraries for Illumina® sequencing were prepared using the Nextera® XT DNA library preparation kit (Illumina, San Diego, CA, United States) according to the manufacturer’s instructions. Briefly, DNA was tagmented at 55 °C for five minutes followed by ligation to Illumina® sequencing index 1 and index 2 adapters by PCR amplification. Size selection and clean-up of the DNA libraries was performed using Agencourt AMPure XP beads (Beckman Coulter, South Kraemer Boulevard Brea, CA, United States). The quantity of DNA was determined on the Qubit 2.0 fluorimeter (Invitrogen, Carlsbad, CA, United States), and a quality check of the libraries was performed on a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, United States). After this, sequencing was performed on an Illumina® MiSeq sequencer (Illumina, San Diego, CA, United States) using a MiSeq reagent kit v3 for 600 cycles (2 × 300 bp paired reads) with a 10% PhiX DNA control spike-in.

4.6. Genome Assembly

The raw reads obtained in FASTQ format were assembled using Geneious Prime® 2019.2.1 (https://www.geneious.com/; [79]). Briefly, the paired-end reads were merged into single reads and trimmed to remove low quality and short reads. The reads were mapped to reference sequences obtained from GenBank. Consensus sequences covering the complete open reading frame (ORF) were submitted to the National Centre for Biotechnology Information (NCBI) GenBank and assigned accession numbers MT271025–MT271035. The ORF lengths were 3267 (VP1), 2673 (VP2), 2508 (VP3), 2328 (VP4), 1194 (VP6), 981 (VP7), 1482 (NSP1), 954 (NSP2), 942 (NSP3), 528 (NSP4), and 594 (NSP5).

4.7. Assignment of Genotypes

The genotypes of each of the 11 rotavirus genome segments were determined using the online Virus Pathogen Resource (ViPR).

4.8. Phylogenetic Analysis

Gene-specific multiple sequence alignments were made using the MAFFT plugin implemented in Geneious Prime® 2019.2.1 and the MUSCLE algorithm embedded in MEGA 6.06 (for the VP2 and NSP1 segments) [80,81]. Once aligned, the DNA Model Test program in MEGA 6.06 was used to identify the optimal evolutionary model for each genome segment [82]. Using an Akaike information criterion (corrected) (AICc), the following models were found to best fit the data: HKY+G+I (VP1), GTR+G+I (VP2, VP3, and VP4), T92+G (VP6, NSP1, NSP2, NSP3, NSP4, and NSP5), and T92+G+I (VP7). Maximum likelihood trees were constructed using the optimal models in MEGA version 6.06 [82,83] with 1000 bootstrap replicates to estimate branch support [84]. The shared nucleotide and amino acid sequence identities among strains were calculated for each gene using the p-distance algorithm in MEGA 6.06. Analysis and visualization of the aligned concatenated whole genomes was performed on the mVISTA online platform [85].

5. Conclusions

In summary, RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] was a reassortant possessing gene segment of porcine and porcine-like human origin, and was closest to Asian strains. It is presumed that pigs play a crucial part as a source for new or newly-evolved emerging human rotaviruses. This highlights the need for continuous large-scale surveillance and whole genome analysis of circulating porcine and human rotaviruses. Furthermore, it was imperative to examine the prevalence of G5P[6] strains in Zambia. Eventually, this should result in a greater understanding of the genes that determine the transmission between hosts successfully as well as to gain insights on complex reassortment patterns between porcine and human rotaviruses.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2076-0817/9/8/663/s1. Supplementary data 1 (S1): Identity matrices for the VP1, VP2, VP3, VP4, VP6, VP7, NSP1, NSP2, NSP3, NSP4, and NSP5 nucleotide and deduced amino acid identities among strains calculated using the p-distance algorithm in MEGA 6.06. Supplementary data 2: Comparison of amino acid sequences. Supplementary data 3: Additional phylograms of the VP1, VP2, VP3, NSP1, NSP2, NSP3, NSP4 and NSP5.

Author Contributions

M.M.N., M.J.M., M.L.S., and J.M.M. conceptualised the main project. W.M.M., P.N.M., J.S., E.M.M., and M.M.N. performed the laboratory experiments. M.M.N., J.S., E.M.M., M.J.M., I.P., M.L.S., and J.M.M. facilitated the obtaining of the sample. Formal analysis was performed by W.M.M., M.D.E., and M.M.N. Data curation was performed by W.M.M., P.N.M., and M.M.N. Writing (original draft preparation) was performed by W.M.M. Review of the drafts was performed by all co-authors. Supervision, project administration, and funding acquisition was conducted by M.M.N. and J.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was principally funded by a grant awarded to M.M.N. by the World Health Organization (agreement number: UFS-AGR17-000378). Other sources of funding that supported the study were grants awarded to M.M.N.: the Bill and Melinda Gates Foundation (BMGF-OPP1180423_2017), South African Medical Research Foundation through the Self-Initiated Research grant (SAMRC-SIR), National Research Foundation (NRF-120814), and a grant from the Poliomyelitis Research Foundation (PRF-19/95).

Acknowledgments

We greatly thank the Zambia team who assisted with sample collection and ELISA testing. We would like to thank Khutso Mothapo, Kebareng Rakau, and Nonkululeko Magagula at the WHO-RRL (Pretoria, South Africa) who assisted with the retrieval of samples. We would also like to thank Sebotsana Rasebotsa, Milton Mogotsi, Emmanuel Ayodeji, Teboho Mooko, and Gilmore Pambuka for their assistance with laboratory work. We are grateful to Stephanus Riekert for technical ICT support.

Conflicts of Interest

The authors declare no conflict of interest.

Disclaimer

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the World Health Organization.

References

  1. Estes, M.K.; Greenberg, H.B. Rotaviruses. In Fields Virology, 6th ed.; Knipe, D.M., Howley, P.M., Eds.; Wolters Kluwer Health/Lippincott, Williams and Wilkins: Philadelphia, PA, USA, 2013; pp. 1347–1401. [Google Scholar]
  2. Tate, J.E.; Burton, A.H.; Boschi-Pinto, C.; Parashar, U.D. World Health Organization—Coordinated Global Rotavirus Surveillance Network. Global, regional, and national estimates of rotavirus mortality in children <5 years of age, 2000–2013. Clin. Infect. Dis. 2016, 62 (Suppl. 2), S96–S105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Troeger, C.; Khalil, I.A.; Rao, P.C.; Cao, S.; Blacker, B.F.; Ahmed, T.; Armah, G.; Bines, J.E.; Brewer, T.G.; Colombara, D.V.; et al. Rotavirus vaccination and the global burden of rotavirus diarrhoea among children younger than 5 years. JAMA Pediatr. 2018, 172, 958–965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Martella, V.; Bányai, K.; Matthijnssens, J.; Buonavoglia, C.; Ciarlet, M. Zoonotic aspects of rotaviruses. Vet. Microbiol. 2010, 140, 246–255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Matthijnssens, J.; Ciarlet, M.; Heiman, E.; Arijs, I.; Delbeke, T.; McDonald, S.M.; Palombo, E.A.; Iturriza-Gomara, M.; Maes, P.; Patton, J.T.; et al. 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. J. Virol. 2008, 82, 3204–3219. [Google Scholar] [CrossRef] [Green Version]
  6. Matthijnssens, J.; Ciarlet, M.; McDonald, S.M.; Attoui, H.; Bányai, K.; Brister, J.R.; Buesa, J.; Esona, M.D.; Estes, M.K.; Gentsch, J.R.; et al. Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG). Arch. Virol. 2011, 156, 1397–1413. [Google Scholar] [CrossRef] [Green Version]
  7. Matthijnssens, J.; Van Ranst, M. Genotype constellation and evolution of group A rotaviruses infecting humans. Curr. Opin. Virol. 2012, 2, 426–433. [Google Scholar] [CrossRef]
  8. RCWG. Rotavirus Classification Working Group–Laboratory of Viral Metagenomics. 2018. Available online: https://rega.kuleuven.be/cev/viralmetagenomics/virus-classification/rcwg (accessed on 5 May 2020).
  9. Ghosh, S.; Kobayashi, N. Whole-genomic analysis of rotavirus strains: Current status and future prospects. Future Microbiol. 2011, 6, 1049–1065. [Google Scholar] [CrossRef]
  10. Iturriza-Gómara, M.; Dallman, T.; Bányai, K.; Böttiger, B.; Buesa, J.; Diedrich, S.; Fiore, L.; Johansen, K.; Koopmans, M.; Korsun, N.; et al. Rotavirus genotypes co-circulating in Europe between 2006 and 2009 as determined by EuroRotaNet, a pan-European collaborative strain surveillance network. Epidemiol. Infect. 2011, 139, 895–909. [Google Scholar] [CrossRef]
  11. Matthijnssens, J.; Heylen, E.; Zeller, M.; Rahman, M.; Lemey, P.; Van Ranst, M. Phylodynamic analyses of rotavirus genotypes G9 and G12 underscore their potential for swift global spread. Mol. Biol. Evol. 2010, 27, 2431–2436. [Google Scholar] [CrossRef] [Green Version]
  12. Patel, M.M.; Steele, D.; Gentsch, J.R.; Wecker, J.; Glass, R.I.; Parashar, U.D. Real-world impact of rotavirus vaccination. Pediatr. Infect. Dis. J. 2011, 30, 1–5. [Google Scholar] [CrossRef]
  13. Rahman, M.; Matthijnssens, J.; Yang, X.; Delbeke, T.; Arijs, I.; Taniguchi, K.; Iturriza-Gomara, M.; Iftekharuddin, N.; Azim, T.; Van Ranst, M. Evolutionary history and global spread of the emerging G12 human rotaviruses. J. Virol. 2007, 81, 2382–2390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Papp, H.; László, B.; Jakab, F.; Ganesh, B.; De Grazia, S.; Matthijnssens, J.; Ciarlet, M.; Martella, V.; Bányai, K. Review of group A rotavirus strains reported in swine and cattle. Vet. Microbiol. 2013, 165, 190–199. [Google Scholar] [CrossRef] [PubMed]
  15. Agbemabiese, C.A.; Nakagomi, T.; Gauchan, P.; Sherchand, J.B.; Pandey, B.D.; Cunliffe, N.A.; Nakagomi, O. Whole genome characterisation of a porcine-like human reassortant G26P[19] Rotavirus A strain detected in a child hospitalised for diarrhoea in Nepal, 2007. Infect. Genet. Evol. 2017, 54, 164–169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kim, H.H.; Matthijnssens, J.; Kim, H.J.; Kwon, H.J.; Park, J.G.; Son, K.Y.; Ryu, E.H.; Kim, D.S.; Lee, W.S.; Kang, M.I.; et al. Full-length genomic analysis of porcine G9P[23] and G9P[7] rotavirus strains isolated from pigs with diarrhoea in South Korea. Infect. Genet. Evol. 2012, 12, 1427–1435. [Google Scholar] [CrossRef] [PubMed]
  17. Martel-Paradis, O.; Laurin, M.A.; Martella, V.; Sohal, J.S.; L’Homme, Y. Full-length genome analysis of G2, G9 and G11 porcine group A rotaviruses. Vet. Microbiol. 2013, 162, 94–102. [Google Scholar] [CrossRef]
  18. Monini, M.; Zaccaria, G.; Ianiro, G.; Lavazza, A.; Vaccari, G.; Ruggeri, F.M. Full-length genomic analysis of porcine rotavirus strains isolated from pigs with diarrhoea in Northern Italy. Infect. Genet. Evol. 2014, 25, 4–13. [Google Scholar] [CrossRef]
  19. Silva, F.D.; Gregori, F.; McDonald, S.M. Distinguishing the genotype 1 genes and proteins of human Wa-like rotaviruses vs. porcine rotaviruses. Infect. Genet. Evol. 2016, 43, 6–14. [Google Scholar] [CrossRef] [Green Version]
  20. Theuns, S.; Heylen, E.; Zeller, M.; Roukaerts, I.D.; Desmarets, L.M.; Van Ranst, M.; Nauwynck, H.J.; Matthijnssens, J. Complete genome characterisation of recent and ancient Belgian pig Group A rotaviruses and assessment of their evolutionary relationship with human rotaviruses. J. Virol. 2015, 89, 1043–1057. [Google Scholar] [CrossRef] [Green Version]
  21. Doro, R.; Farkas, S.L.; Martella, V.; Banyai, K. Zoonotic transmission of rotavirus: Surveillance and control. Expert Rev. Anti. Infect. Ther. 2015, 13, 1337–1350. [Google Scholar] [CrossRef]
  22. Cowley, D.; Donato, C.M.; Roczo-Farkas, S.; Kirkwood, C.D. Novel G10P[14] Rotavirus Strain, Northern Territory, Australia. Emerg. Infect. Dis. 2013, 19, 1324–1327. [Google Scholar] [CrossRef]
  23. Komoto, S.; Tacharoenmuang, R.; Guntapong, R.; Ide, T.; Sinchai, P.; Upachai, S.; Fukuda, S.; Yoshikawa, T.; Tharmaphornpilas, P.; Sangkitporn, S.; et al. Identification and characterisation of a human G9P[23] rotavirus strain from a child with diarrhoea in Thailand: Evidence for porcine-to-human interspecies transmission. J. Gen. Virol. 2017, 98, 532–538. [Google Scholar] [CrossRef] [PubMed]
  24. Malasao, R.; Khamrin, P.; Kumthip, K.; Ushijima, H.; Maneekarn, N. Complete genome sequence analysis of rare G4P[6] rotavirus strains from human and pig reveals the evidence for interspecies transmission. Infect. Genet. Evol. 2018, 65, 357–368. [Google Scholar] [CrossRef] [PubMed]
  25. Mukherjee, A.; Ghosh, S.; Bagchi, P.; Dutta, D.; Chattopadhyay, S.; Kobayashi, N.; Chawla-Sarkar, M. Full genomic analyses of human rotavirus G4P[4], G4P[6], G9P[19] and G10P[6] strains from north-eastern India: Evidence for interspecies transmission and complex reassortment events. Clin. Microbiol. Infect. 2011, 17, 1343–1346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. My, P.V.; Rabaa, M.A.; Donato, C.; Cowley, D.; Phat, V.V.; Dung, T.T.; Anh, P.H.; Vinh, H.; Bryant, J.E.; Kellam, P.; et al. Novel porcine-like human G26P[19] rotavirus identified in hospitalised paediatric diarrhoea patients in Ho Chi Minh City, Vietnam. J. Gen. Virol. 2014, 95, 2727–2733. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Quaye, O.; Roy, S.; Rungsrisuriyachai, K.; Esona, M.D.; Xu, Z.; Tam, K.I.; Banegas, D.J.; Rey-Benito, G.; Bowen, M.D. Characterisation of a rare, reassortant human G10P[14] rotavirus strain detected in Honduras. Mem. Inst. Oswaldo Cruz. 2018, 113, 9–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Tacharoenmuang, R.; Komoto, S.; Guntapong, R.; Ide, T.; Singchai, P.; Upachai, S.; Fukuda, S.; Yoshida, Y.; Murata, T.; Yoshikawa, T.; et al. Characterisation of a G10P[14] rotavirus strain from a diarrhoeic child in Thailand: Evidence for bovine-to-human zoonotic transmission. Infect. Genet. Evol. 2018, 63, 43–57. [Google Scholar] [CrossRef]
  29. Bwogi, J.; Jere, K.C.; Karamagi, C.; Byarugaba, D.K.; Namuwulya, P.; Baliraine, F.N.; Desselberger, U.; Iturriza-Gomara, M. Whole genome analysis of selected human and animal rotaviruses identified in Uganda from 2012 to 2014 reveals complex genome reassortment events between human, bovine, caprine and porcine strains. PLoS ONE 2017, 12, e0178855. [Google Scholar] [CrossRef] [Green Version]
  30. Matthijnssens, J.; Bilcke, J.; Ciarlet, M.; Martella, V.; Bányai, K.; Rahman, M.; Zeller, M.; Beutels, P.; Van Damme, P.; Van Ranst, M. Rotavirus disease and vaccination: Impact on genotype diversity. Future Microbiol. 2009, 4, 1303–1316. [Google Scholar] [CrossRef]
  31. Nyaga, M.M.; Jere, K.C.; Esona, M.D.; Seheri, M.L.; Stucker, K.M.; Halpin, R.A.; Akopov, A.; Stockwell, T.B.; Peenze, I.; Diop, A.; et al. Whole genome detection of rotavirus mixed infections in human, porcine and bovine samples co-infected with various rotavirus strains collected from sub-Saharan Africa. Infect. Genet. Evol. 2015, 31, 321–334. [Google Scholar] [CrossRef] [Green Version]
  32. Beards, G.; Graham, C. Temporal distribution of rotavirus G-serotypes in the West Midlands region of the United Kingdom, 1983–1994. J. Diarrhoeal Dis. Res. 1995, 13, 235–237. [Google Scholar]
  33. Bok, K.; Castagnaro, N.; Borsa, A.; Nates, S.; Espul, C.; Fay, O.; Fabri, A.; Grinstein, S.; Miceli, I.; Matson, D.O.; et al. Surveillance for rotavirus in Argentina. J. Med. Virol. 2001, 65, 190–198. [Google Scholar] [CrossRef] [PubMed]
  34. Esona, M.D.; Armah, G.E.; Geyer, A.; Steele, A.D. Detection of an unusual human rotavirus strain with G5P[8] specificity in a Cameroonian child with diarrhoea. J. Clin. Microbiol. 2004, 42, 441–444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Esona, M.D.; Geyer, A.; Banyai, K.; Page, N.; Aminu, M.; Armah, G.E.; Hull, J.; Steele, D.A.; Glass, R.I.; Gentsch, J.R. Novel human rotavirus genotype G5P[7] from child with diarrhoea, Cameroon. Emerg. Infect. Dis. 2009, 15, 83–86. [Google Scholar] [CrossRef] [PubMed]
  36. Gouvea, V.; de Castro, L.; Timenetsky Mdo, C.; Greenberg, H.; Santos, N. Rotavirus serotype G5 associated with diarrhoea in Brazilian children. J. Clin. Microbiol. 1994, 32, 1408–1409, Erratum appears in J. Clin. Microbiol. 1994, 32, 1834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hwang, K.P.; Wu, F.T.; Bányai, K.; Wu, H.S.; Yang, D.C.; Huang, Y.C.; Lin, J.S.; Hsiung, C.A.; Huang, J.C.; Jiang, B.; et al. Identification of porcine rotavirus-like genotype P[6] strains in Taiwanese children. J. Med. Microbiol. 2012, 61, 990–997. [Google Scholar] [CrossRef]
  38. Lorenzetti, E.; Da Silva Medeiros, T.N.; Alfieri, A.F.; Alfieri, A.A. Genetic heterogeneity of wild-type G4P[6] porcine rotavirus strains detected in a diarrhoea outbreak in a regularly vaccinated pig herd. Vet. Microbiol. 2011, 154, 191–196. [Google Scholar] [CrossRef]
  39. Martella, V.; Ciarlet, M.; Bányai, K.; Lorusso, E.; Cavalli, A.; Corrente, M.; Elia, G.; Arista, S.; Camero, M.; Desario, C.; et al. Identification of a novel VP4 genotype carried by a serotype G5 porcine rotavirus strain. Virology 2006, 346, 301–311. [Google Scholar] [CrossRef] [Green Version]
  40. Nyaga, M.M.; Tan, Y.; Seheri, M.L.; Halpin, R.A.; Akopov, A.; Stucker, K.M.; Fedorova, N.B.; Shrivastava, S.; Duncan Steele, A.; Mwenda, J.M.; et al. Whole-genome sequencing and analyses identify high genetic heterogeneity, diversity and endemicity of rotavirus genotype P[6] strains circulating in Africa. Infect. Genet. Evol. 2018, 63, 79–88. [Google Scholar] [CrossRef]
  41. Li, D.D.; Duan, Z.J.; Zhang, Q.; Liu, N.; Xie, Z.P.; Jiang, B.; Steele, D.; Jiang, X.; Wang, Z.S.; Fang, Z.Y. Molecular characterisation of unusual human G5P[6] rotaviruses identified in China. J. Clin. Virol. 2008, 42, 141–148. [Google Scholar] [CrossRef]
  42. Ahmed, K.; Dang, D.A.; Nakagomi, O. Rotavirus G5P[6] in child with diarrhoea, Vietnam. Emerg. Infect. Dis. 2007, 13, 1232–1235. [Google Scholar] [CrossRef]
  43. Chieochansin, T.; Vutithanachot, V.; Phumpholsup, T.; Posuwan, N.; Theamboonlers, A.; Poovorawan, Y. The prevalence and genotype diversity of Human Rotavirus A circulating in Thailand, 2011–2014. Infect. Genet. Evol. 2016, 37, 129–136. [Google Scholar] [CrossRef]
  44. Duan, Z.J.; Li, D.D.; Zhang, Q.; Liu, N.; Huang, C.P.; Jiang, X.; Jiang, B.; Glass, R.; Steele, D.; Tang, J.Y.; et al. Novel human rotavirus of genotype G5P[6] identified in a stool specimen from a Chinese girl with diarrhoea. J. Clin. Microbiol. 2007, 45, 1614–1617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Komoto, S.; Maeno, Y.; Tomita, M.; Matsuoka, T.; Ohfu, M.; Yodoshi, T.; Akeda, H.; Taniguchi, K. Whole genomic analysis of a porcine-like human G5P[6] rotavirus strain isolated from a child with diarrhoea and encephalopathy in Japan. J. Gen. Virol. 2013, 94, 1568–1575. [Google Scholar] [CrossRef] [PubMed]
  46. Mladenova, Z.; Papp, H.; Lengyel, G.; Kisfali, P.; Steyer, A.; Steyer, A.F.; Esona, M.D.; Iturriza-Gómara, M.; Bányai, K. Detection of rare reassortant G5P[6] rotavirus, Bulgaria. Infect. Genet. Evol. 2012, 12, 1676–1684. [Google Scholar] [CrossRef] [PubMed]
  47. Chilengi, R.; Rudd, C.; Bolton, C.; Guffey, B.; Masumbu, P.K.; Stringer, J. Successes, challenges and lessons learned in accelerating introduction of rotavirus immunisation in Zambia. World J. Vaccines 2015, 5, 43–53. [Google Scholar] [CrossRef] [Green Version]
  48. Zambia Ministry of Health. The 2012 Annual Health Statistical Bulletin. 2014. Available online: https://www.moh.gov.zm/docs/reports/2012_Annual_Health_Statistical_Bulletin_Version_1.pdf (accessed on 18 March 2020).
  49. Mpabalwani, M.; Oshitani, H.; Kasolo, F.; Mizuta, K.; Luo, N.; Matsubayashi, N.; Bhat, G.; Suzuki, H.; Numazaki, Y. Rotavirus gastro-enteritis in hospitalised children with acute diarrhoea in Zambia. Ann. Trop. Paediatr. 1995, 15, 39–43. [Google Scholar] [CrossRef]
  50. Mpabalwani, E.M.; Simwaka, C.J.; Mwenda, J.M.; Mubanga, C.P.; Monze, M.; Matapo, B.; Parashar, U.D.; Tate, J.E. Impact of rotavirus vaccination on diarrhoeal Hospitalisations in children aged <5 Years in Lusaka, Zambia. Clin. Infect. Dis. 2016, 62, S183–S187. [Google Scholar] [CrossRef] [Green Version]
  51. Mpabalwani, E.M.; Simwaka, J.C.; Mwenda, J.M.; Matapo, B.; Parashar, U.D.; Tate, J.E. Sustained impact of rotavirus vaccine on rotavirus hospitalisations in Lusaka, Zambia, 2009–2016. Vaccine 2018, 36, 7165–7169. [Google Scholar] [CrossRef]
  52. WHO Vaccine-Preventable Diseases: Monitoring System. 2020 Global Summary. Available online: https://apps.who.int/immunization_monitoring/globalsummary/countries?countrycriteria%5Bcountry%5D%5B%5D=ZMB (accessed on 7 August 2020).
  53. Mwenda, J.M.; Tate, J.E.; Parashar, U.D.; Mihigo, R.; Agócs, M.; Serhan, F.; Nshimirimana, D. African Rotavirus Surveillance Network: A brief overview. Pediatr. Infect. Dis. J. 2014, 33, S6–S8. [Google Scholar] [CrossRef]
  54. Mwenda, J.M.; Burke, R.M.; Shaba, K.; Mihigo, R.; Tevi-Benissan, M.C.; Mumba, M.; Biey, J.N.; Cheikh, D.; Poy, A.; Zawaira, F.R.; et al. Implementation of rotavirus surveillance and vaccine introduction—World Health Organization African region, 2007–2016. Morb. Mortal. Wkly. Rep. 2017, 66, 1192–1196. [Google Scholar] [CrossRef] [Green Version]
  55. Mwenda, J.M.; Ntoto, K.M.; Abebe, A.; Enweronu-Laryea, C.; Amina, I.; Mchomvu, J.; Kisakye, A.; Mpabalwani, E.M.; Pazvakavambwa, I.; Armah, G.E.; et al. Burden and epidemiology of rotavirus diarrhoea in selected African countries: Preliminary results from the African Rotavirus Surveillance Network. J. Infect. Dis. 2010, 202, S5–S11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Dong, H.J.; Qian, Y.; Huang, T.; Zhu, R.N.; Zhao, L.Q.; Zhang, Y.; Li, R.C.; Li, Y.P. Identification of circulating porcine-human reassortant G4P[6] rotavirus from children with acute diarrhoea in China by whole genome analyses. Infect. Genet. Evol. 2013, 20, 155–162. [Google Scholar] [CrossRef] [PubMed]
  57. Yahiro, T.; Takaki, M.; Chandrasena, T.G.A.; Rajindrajith, S.; Iha, H.; Ahmed, K. Human-porcine reassortant rotavirus generated by multiple reassortment events in a Sri Lankan child with diarrhoea. Infect. Genet. Evol. 2018, 65, 170–186. [Google Scholar] [CrossRef] [PubMed]
  58. Zhou, X.; Wang, Y.H.; Ghosh, S.; Tang, W.F.; Pang, B.B.; Liu, M.Q.; Peng, J.S.; Zhou, D.J.; Kobayashi, N. Genomic characterisation of G3P[6], G4P[6] and G4P[8] human rotaviruses from Wuhan, China: Evidence for interspecies transmission and reassortment events. Infect. Genet. Evol. 2015, 33, 55–71. [Google Scholar] [CrossRef] [PubMed]
  59. Heylen, E.; Likele, B.B.; Zeller, M.; Stevens, S.; De Coster, S.; Conceição-Neto, N.; Van Geet, C.; Jacobs, J.; Ngbonda, D.; Van Ranst, M.; et al. Rotavirus surveillance in Kisangani, the Democratic Republic of the Congo, reveals a high number of unusual genotypes and gene segments of animal origin in non-vaccinated symptomatic children. PLoS ONE 2014, 9, e100953. [Google Scholar] [CrossRef] [PubMed]
  60. Kaneko, M.; Do, L.P.; Doan, Y.H.; Nakagomi, T.; Gauchan, P.; Agbemabiese, C.A.; Dang, A.D.; Nakagomi, O. Porcine-like G3P[6] and G4P[6] rotavirus A strains detected from children with diarrhoea in Vietnam. Arch. Virol. 2018, 163, 2261–2263. [Google Scholar] [CrossRef] [PubMed]
  61. Phan, M.V.; Anh, P.H.; Cuong, N.V.; Munnink, B.B.; Hoek, L.; My, P.T.; Tri, T.N.; Bryant, J.E.; Baker, S.; Thwaites, G.; et al. Unbiased whole-genome deep sequencing of human and porcine stool samples reveals circulation of multiple groups of rotaviruses and a putative zoonotic infection. Virus Evol. 2016, 2, 1–15. [Google Scholar] [CrossRef] [Green Version]
  62. Wang, Y.H.; Kobayashi, N.; Zhou, D.J.; Yang, Z.Q.; Zhou, X.; Peng, J.S.; Zhu, Z.R.; Zhao, D.F.; Liu, M.Q.; Gong, J. Molecular epidemiologic analysis of group A rotaviruses in adults and children with diarrhoea in Wuhan city, China, 2000–2006. Arch. Virol. 2007, 152, 669–685. [Google Scholar] [CrossRef]
  63. da Silva, M.F.; Tort, L.F.; Goméz, M.M.; Assis, R.M.; Volotão, E.d.M.; de Mendonça, M.C.; Bello, G.; Leite, J.P. VP7 Gene of human rotavirus A genotype G5: Phylogenetic analysis reveals the existence of three different lineages worldwide. J. Med. Virol. 2011, 83, 357–366. [Google Scholar] [CrossRef]
  64. Green, K.Y.; Hoshino, Y.; Ikegami, N. Sequence analysis of the gene encoding the serotype-specific glycoprotein (VP7) of two new human rotavirus serotypes. Virology 1989, 168, 429–433. [Google Scholar] [CrossRef]
  65. Martella, V.; Bányai, K.; Ciarlet, M.; Iturriza-Gómara, M.; Lorusso, E.; De Grazia, S.; Arista, S.; Decaro, N.; Elia, G.; Cavalli, A.; et al. Relationships among porcine and human P[6] rotaviruses: Evidence that the different human P[6] lineages have originated from multiple interspecies transmission events. Virology 2006, 344, 509–519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Burke, B.; Bridger, J.C.; Desselberger, U. Temporal correlation between a single amino acid change in the VP4 of a porcine rotavirus and a marked change in pathogenicity. Virology 1994, 202, 754–759. [Google Scholar] [CrossRef] [PubMed]
  67. Mackow, E.R.; Shaw, R.D.; Matsui, S.M.; Vo, P.T.; Dang, M.N.; Greenberg, H.B. The rhesus rotavirus gene encoding protein VP3: Location of amino acids involved in homologous and heterologous rotavirus neutralization and identification of a putative fusion region. Proc. Natl. Acad. Sci. USA 1988, 85, 645–649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Arias, C.F.; Romero, P.; Alvarez, V.; López, S. Trypsin activation pathway of rotavirus infectivity. J. Virol. 1996, 70, 5832–5839. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Steyer, A.; Poljšak-Prijatelj, M.; Barlič-Maganja, D.; Marin, J. Human, porcine and bovine rotaviruses in Slovenia: Evidence of interspecies transmission and genome reassortment. J. Gen. Virol. 2008, 89, 1690–1698. [Google Scholar] [CrossRef]
  70. Zeller, M.; Heylen, E.; De Coster, S.; Van Ranst, M.; Matthijnssens, J. Full genome characterisation of a porcine-like human G9P[6] rotavirus strain isolated from an infant in Belgium. Infect. Genet. Evol. 2012, 12, 1492–1500. [Google Scholar] [CrossRef]
  71. Zhang, S.; McDonald, P.W.; Thompson, T.A.; Dennis, A.F.; Akopov, A.; Kirkness, E.F.; Patton, J.T.; McDonald, S.M. Analysis of human rotaviruses from a single location over an 18-Year time span suggests that protein co-adaption influences gene constellations. J. Virol. 2014, 88, 9842–9863. [Google Scholar] [CrossRef] [Green Version]
  72. Ciarlet, M.; Ludert, J.E.; Liprandi, F. Comparative amino acid sequence analysis of the major outer capsid protein (VP7) of porcine rotaviruses with G3 and G5 serotype specificities isolated in Venezuela and Argentina. Arch. Virol. 1995, 140, 437–451. [Google Scholar] [CrossRef]
  73. Dyall-Smith, M.L.; Lazdins, I.; Tregear, G.W.; Holmes, I.H. Location of the major antigenic sites involved in rotavirus serotype-specific neutralisation. Proc. Natl. Acad. Sci. USA 1986, 83, 3465–3468. [Google Scholar] [CrossRef] [Green Version]
  74. Kobayashi, N.; Taniguchi, K.; Urasawa, S. Analysis of the newly identified neutralisation epitopes on VP7 of human rotavirus serotype 1. J. Gen. Virol. 1991, 72, 117–124. [Google Scholar] [CrossRef]
  75. Green, K.Y.; Sears, J.F.; Taniguchi, K.; Midthun, K.; Hoshino, Y.; Gorziglia, M.; Nishikawa, K.; Urasawa, S.; Kapikian, A.Z.; Chanock, R.M. Prediction of human rotavirus serotype by nucleotide sequence analysis of the VP7 protein gene. J. Virol. 1988, 62, 1819–1823. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Gorziglia, M.; Larralde, G.; Kapikian, A.Z.; Chanock, R.M. Antigenic relationships among human rotaviruses as determined by outer capsid protein VP4. Proc. Natl. Acad. Sci. USA 1990, 87, 7155–7159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Bányai, K.; Martella, V.; Jakab, F.; Melegh, B.; Szücs, G. Sequencing and phylogenetic analysis of human genotype P[6] rotavirus strains detected in Hungary provides evidence for genetic heterogeneity within the P[6] VP4 gene. J. Clin. Microbiol. 2004, 42, 4338–4343. [Google Scholar] [CrossRef] [Green Version]
  78. Potgieter, A.C.; Page, N.A.; Liebenberg, J.; Wright, I.M.; Landt, O.; van Dijk, A.A. Improved strategies for sequence-independent amplification and sequencing of viral double-stranded RNA genomes. J. Gen. Virol. 2009, 90, 1423–1432. [Google Scholar] [CrossRef]
  79. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  80. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 1–19. [Google Scholar] [CrossRef] [Green Version]
  81. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  82. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA 6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  83. Guindon, S.; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [Green Version]
  84. Felsenstein, J. Confidence Limits on Phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef]
  85. Frazer, K.A.; Pachter, L.; Poliakov, A.; Rubin, E.M.; Dubchak, I. VISTA: Computational tools for comparative genomics. Nucleic Acids Res. 2004, 32, 273–279. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogenetic tree constructed from the nucleotide sequences of the VP7 genes of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and representative strains. The position of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] is shown by the black square (▪). Reference strains obtained from GenBank are represented by accession number, strain name, country, and year of isolation. The three closest strains, as identified by BLASTn, are also included. Bootstrap values ≥70% are shown adjacent to each branch node. Scale bar: 0.05 substitutions per nucleotide.
Figure 1. Phylogenetic tree constructed from the nucleotide sequences of the VP7 genes of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and representative strains. The position of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] is shown by the black square (▪). Reference strains obtained from GenBank are represented by accession number, strain name, country, and year of isolation. The three closest strains, as identified by BLASTn, are also included. Bootstrap values ≥70% are shown adjacent to each branch node. Scale bar: 0.05 substitutions per nucleotide.
Pathogens 09 00663 g001
Figure 2. Phylogenetic tree constructed from the nucleotide sequences of the VP4 genes of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and representative strains. The position of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] is shown by the black square (▪). Reference strains obtained from GenBank are represented by accession number, strain name, country, and year of isolation. The three closest strains, as identified by BLASTn, are also included. Bootstrap values ≥70% are shown adjacent to each branch node. Scale bar: 0.05 substitutions per nucleotide.
Figure 2. Phylogenetic tree constructed from the nucleotide sequences of the VP4 genes of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and representative strains. The position of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] is shown by the black square (▪). Reference strains obtained from GenBank are represented by accession number, strain name, country, and year of isolation. The three closest strains, as identified by BLASTn, are also included. Bootstrap values ≥70% are shown adjacent to each branch node. Scale bar: 0.05 substitutions per nucleotide.
Pathogens 09 00663 g002
Figure 3. Phylogenetic tree constructed from the nucleotide sequences of the VP6 genes of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and representative strains. The position of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] is shown by the black square (▪). Reference strains obtained from GenBank are represented by accession number, strain name, country, and year of isolation. The three closest strains, as identified by BLASTn, are also included. Bootstrap values ≥70% are shown adjacent to each branch node. Scale bar: 0.05 substitutions per nucleotide.
Figure 3. Phylogenetic tree constructed from the nucleotide sequences of the VP6 genes of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] and representative strains. The position of strain RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6] is shown by the black square (▪). Reference strains obtained from GenBank are represented by accession number, strain name, country, and year of isolation. The three closest strains, as identified by BLASTn, are also included. Bootstrap values ≥70% are shown adjacent to each branch node. Scale bar: 0.05 substitutions per nucleotide.
Pathogens 09 00663 g003
Figure 4. mVISTA whole genome nucleotide alignment comparing the Zambian G5P[6] strain (RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014G5P[6]) with the G5P[6] strain from Japan (Ryukyu-1120), whose whole genome sequence had been determined, and with selected porcine-like human P[6] strains from China (GX54 and E931). Strain names are shown on the left, and the proteins VP1-VP4, VP6-VP7, and NSP1-NSP5 are indicated on the top. The bottom scale indicates distance in kb. Percentile values on the right indicate sequence-based similarity between the study strain and the respective reference strains. Shading indicates the level of conservation.
Figure 4. mVISTA whole genome nucleotide alignment comparing the Zambian G5P[6] strain (RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014G5P[6]) with the G5P[6] strain from Japan (Ryukyu-1120), whose whole genome sequence had been determined, and with selected porcine-like human P[6] strains from China (GX54 and E931). Strain names are shown on the left, and the proteins VP1-VP4, VP6-VP7, and NSP1-NSP5 are indicated on the top. The bottom scale indicates distance in kb. Percentile values on the right indicate sequence-based similarity between the study strain and the respective reference strains. Shading indicates the level of conservation.
Pathogens 09 00663 g004
Table 1. The segment and ORF lengths of strain UFS-NGS-MRC-DPRU4723 and the highest sequence identities obtained using the Basic Local Alignment Search Tool (BLAST).
Table 1. The segment and ORF lengths of strain UFS-NGS-MRC-DPRU4723 and the highest sequence identities obtained using the Basic Local Alignment Search Tool (BLAST).
GENOME SEGMENT EncodingGenBank Accession no.Segment LengthORF LengthResults of Blast Search
Most Similar StrainGenBank Accession no.Similarity (%)Reference
VP1MT27102533023267GX54KF04144196.7[56]
VP2MT27102626732673R1207LC38988696.5[57]
VP3MT27102725912508R946KF72606095.7[58]
VP4MT27102823592328KisB332KJ87090398.0[59]
NSP1MT27102915121482NT0042LC09589498.1[60]
VP6MT27103013561194KYE-14-A048KX98827998.7[29]
NSP3MT271031107694212070-4KX36328797.1[61]
NSP2MT271032954954YNKJ46698796.8[https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/nuccore/KJ466987]
VP7MT2710331054981JN-2KT82077798.0[https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/nuccore/KT820777]
NSP4MT27103475152814150-54KX36335497.7[61]
NSP5MT271035644594R479GU18955997.6[62]
Table 2. Genotype natures of the 11 gene segments of Zambian strain UFS-NGS-MRC-DPRU4723 compared with those of selected human and porcine strains.
Table 2. Genotype natures of the 11 gene segments of Zambian strain UFS-NGS-MRC-DPRU4723 compared with those of selected human and porcine strains.
StrainGenotype
VP7VP4VP6VP1VP2VP3NSP1NSP2NSP3NSP4NSP5
RVA/Human-wt/ZMB/UFS-NGS-MRC-DPRU4723/2014/G5P[6]G5P[6]I1R1C1M1A8N1T1E1H1
RVA/Human-wt/BGR/BG260/2008/G5P[6] *G5P[6]I1R1C1M1A8N1T1E1H1
RVA/Human-wt/JPN/Ryukyu-1120/2011/G5P[6]G5P[6]I5R1C1M1A8N1T1E1H1
RVA/Human-wt/CHN/LL3354/2000/G5P[6]G5P[6]I5------E1-
RVA/Human-wt/CHN/LL4260/2001/G5P[6]G5P[6]-------E1-
RVA/Human-wt/CHN/LL36755/2003/G5P[6]G5P[6]-------E1-
RVA/Human-wt/VNM/KH210/2004/G5P[6]G5P[6]-------E1-
RVA/Human-wt/TWN/03-98P50/2009/G5P[6] *G5P[6]I5------E1-
RVA/Human-wt/CMR/6784/ARN/2000/G5P[7]G5P[7]I5R1C1M1A1N1T1E1H1
RVA/Human-tc/BRA/IAL28/1992/G5P[8]G5P[8]I5R1C1M1A1N1T1E1H1
RVA/Pig-tc/USA/OSU/1975/G5P[7]G5P[7]I5R1C1M1A1N1T1E1H1
RVA/Pig-wt/BEL/12R002/2012/G5P[7]G5P[7]I5R1C1M1A8N1T7E1H1
RVA/Pig-wt/JPN/BU2/2014/G5P[7]G5P[7]I5R1C1M1A8N1T1E1H1
RVA/Human-tc/USA/Wa/1974/G1P[8]G1P[8]I1R1C1M1A1N1T1E1H1
RVA/Human-tc/USA/DS-1/1976/G2P[4]G2P[4]I2R2C2M2A2N2T2E2H2
RVA/Human-tc/JPN/AU-1/1982/G3P[9]G3P[9]I3R3C3M3A3N3T3E3H3
RVA/Pig-wt/BEL/12R006/2012/G3P[6]G3P[6]I5R1C1M1A8N1T1E1H1
RVA/Human-tc/GBR/ST3/1974/G4P[6]G4P[6]I1R1C1M1A1N1T1E1H1
RVA/Pig-tc/USA/Gottfried/1975/G4P[6]G4P[6]I1R1C1M1A8N1T1E1H1
RVA/Human-tc/CHN/R479/2004/G4P[6]G4P[6]I5R1C1M1A1N1T7E1H1
RVA/Human-wt/CHN/E931/2008/G4P[6]G4P[6]I1R1C1M1A8N1T1E1H1
RVA/Human-wt/COD/KisB332/2008/G4P[6]G4P[6]I1R1C1M1A1N1T7E1H1
RVA/Human-wt/CHN/GX54/2010/G4P[6]G4P[6]I1R1C1M1A8N1T1E1H1
RVA/Pig-wt/BEL/12R005/2012/G4P[7]G4P[7]I5R1C1M1A8N1T7E1H1
RVA/Human-wt/BEL/BE2001/2009/G9P[6]G9P[6]I5R1C1M1A8N1T7E1H1
RVA/Human-tc/USA/WI61/1983/G9P[8]G9P[8]I1R1C1M1A1N1T1E1H1
RVA/Human-wt/BEL/B3458/2003/G9P[8]G9P[8]I1R1C1M1A1N1T1E1H1
RVA/Human-tc/IND/mani-97/2006/G9P[19]G9P[19]I5R1C1M1A8N1T1E1H1
RVA/Human-wt/BGD/Dhaka6/2001/G11P[25]G11P[25]I1R1C1M1A1N1T1E1H1
RVA/Human-wt/VNM/30378/2009/G26P[19]G26P[19]I5R1C1M1A8N1T1E1H1
RVA/Human-wt/BRA/rj24598/2015/G26P[19]G26P[19]I5R1C1M1A8N1T1E1H1
Blue shading indicates the gene segments with genotypes identical to those of UFS-NGS-MRC-DPRU4723. Bold font indicates genotypes associated with porcine strains. “−” indicates that no sequence data were available in GenBank/EMBL/DDBJ data banks. * Genotype assignment based on reports by [37] (strain 03-98sP50) and (strain BG260) [46]. To date, the nucleotide accession numbers for the 11 gene segments of strains 03-98sP50 and BG260 are not available in the GenBank, EMBL, or DDBJ data banks.

Share and Cite

MDPI and ACS Style

Maringa, W.M.; Mwangi, P.N.; Simwaka, J.; Mpabalwani, E.M.; Mwenda, J.M.; Peenze, I.; Esona, M.D.; Mphahlele, M.J.; Seheri, M.L.; Nyaga, M.M. Molecular Characterisation of a Rare Reassortant Porcine-Like G5P[6] Rotavirus Strain Detected in an Unvaccinated Child in Kasama, Zambia. Pathogens 2020, 9, 663. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens9080663

AMA Style

Maringa WM, Mwangi PN, Simwaka J, Mpabalwani EM, Mwenda JM, Peenze I, Esona MD, Mphahlele MJ, Seheri ML, Nyaga MM. Molecular Characterisation of a Rare Reassortant Porcine-Like G5P[6] Rotavirus Strain Detected in an Unvaccinated Child in Kasama, Zambia. Pathogens. 2020; 9(8):663. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens9080663

Chicago/Turabian Style

Maringa, Wairimu M., Peter N. Mwangi, Julia Simwaka, Evans M. Mpabalwani, Jason M. Mwenda, Ina Peenze, Mathew D. Esona, M. Jeffrey Mphahlele, Mapaseka L. Seheri, and Martin M. Nyaga. 2020. "Molecular Characterisation of a Rare Reassortant Porcine-Like G5P[6] Rotavirus Strain Detected in an Unvaccinated Child in Kasama, Zambia" Pathogens 9, no. 8: 663. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens9080663

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop