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.

Group A rotavirus (RVA) infections form a major public health problem, especially in low-income countries like the Democratic Republic of the Congo (COD). However, limited data on RVA diversity is available from sub-Saharan Africa in general and the COD in particular. Therefore, the first aim of this study was to determine the genetic diversity of 99 RVAs detected during 2007-2010 in Kisangani, COD. The predominant G-type was G1 (39%) and the most predominant P-type was P[6] (53%). A total of eight different G/P-combinations were found: G1P[8] (28%), G8P[6] (26%), G2P[4] (14%), G12P[6] (13%), G1P[6] (11%), G9P[8] (4%), G4P[6] (2%) and G8P[4] (1%). The second aim of this study was to gain insight into the diversity of P[6] RVA strains in the COD. Therefore, we selected five P[6] RVA strains in combination with the G1, G4, G8 (2x) or G12 genotype for complete genome analysis. Complete genome analysis showed that the genetic background of the G1P[6] and G12P[6] strains was entirely composed of genotype 1 (Wa-like), while the segments of the two G8P[6] strains were identified as genotype 2 (DS-1-like). Interestingly, all four strains possessed a NSP4 gene of animal origin. The analyzed G4P[6] RVA strain was found to possess the unusual G4-P[6]-I1-R1-C1-M1-A1-N1-T7-E1-H1 constellation. Although the majority of its genes (if not all), were presumably of porcine origin, this strain was able to cause gastro-enteritis in humans. The high prevalence of unusual RVA strains in the COD highlights the need for continued surveillance of RVA diversity in the COD. These results also underline the importance of complete genetic characterization of RVA strains and indicate that reassortments and interspecies transmission among human and animal RVAs strains occur regularly. Based on these data, RVA vaccines will be challenged with a wide variety of different RVA strain types in the COD.

Introduction

Group A rotaviruses (RVAs) are the world’s leading cause of severe diarrhea in children <5 years of age and form a major public health problem, with an estimated 453,000 deaths per year worldwide. More than 50% of these deaths occur in Africa (232,000) and 7% of all RVA related deaths occur in the Democratic Republic of the Congo (COD) (32,653; 95% CI 27,804–37,699) [1]. In fact, the COD, together with India, Nigeria, Pakistan and Ethiopia, is one of the countries that bear the highest RVA mortality rates worldwide.

Despite this high burden of disease, only limited regional and country specific data on RVA genotypes diversity is available from sub-Saharan Africa in general and the COD in particular. The last report of the African Rotavirus Surveillance Network reported the distribution of predominant RVA genotypes among hospitalized children aged <5 years between 2007 and 2011 from 16 African countries [2]. The major G- and P-genotype combinations detected in this study were G1P[8] (18.4%), G9P[8] (11.7%), G2P[4] (8.6%), G2P[6] (6.2%), G1P[6] (4.9%), G3P[6] (4.3%), G8P[6] (3.8%) and G12P[8] (3.1%). Worldwide the majority of RVA strains infecting humans belong to the G1-G4, G9 and G12 G-genotypes, in combination with the P[4] (for G2) and P[8] (for G1, G3, G4, G9 and G12) genotypes. When focusing on the differences in the genetic setup of RVA strains circulating in Africa, compared to the ones circulating in industrialized countries, multiple studies have reported the high detection rate of RVA strains containing the P[6] and G8 genotypes. Both genotypes are believed to be of animal origin: P[6] RVA strains (most likely) originated from pigs, while G8 is a genotype often detected in cows or other members of the mammalian order of the Artiodactyla [3]–[6] rotaviruses: evidence that the different human P[6] lineages have originated from multiple interspecies transmission events. Virology. 2006 ">10]. The majority of known human G8 RVAs strains have been described all over the African continent in combination with a large number of VP4 specificities and with different genetic backgrounds [11]–[6]. Infect Genet Evol. 2012 ">23].

Two human RVA vaccines, Rotarix (GlaxoSmithKline) and RotaTeq (Merck) have been licensed and proven safe in large clinical trials in the Americas and Europe [24], [25]. Data from clinical trials conducted in five African countries have shown a lower vaccine immunogenicity and efficacy compared to the ones observed in America and Europe [26]–[28]. Despite this observed lower vaccine efficacy, vaccine introduction in Africa will prevent a greater number of RVA related deaths due to the substantially larger mortality rates in Africa compared to the Americas and Europe. Between August 2009 and May 2014, 21 African countries (in chronological order: South-Africa, Botswana, Morocco, Sudan, Ghana, Rwanda, Malawi, Tanzania, The Gambia, Burkina Faso, Ethiopia, Libya, Zambia, Burundi, Mali, Cameroon, Sierra Leone, The Republic of the Congo, Angola, Madagascar and Zimbabwe) have implemented RVA vaccination into their national immunization program [29], [30].

What factors determine the lower vaccine efficacy in Africa is not yet completely understood but both host (malnutrition, competing enteric pathogens, differences in histo-blood group antigens, the level of anti-rotavirus antibodies in breast milk and a lower immunological response) and viral (a higher RVA strain diversity) characteristics have been proposed to (partially) explain this observation [31]–[33].

Specific data on RVA strain diversity from the Democratic Republic of the Congo are scarce with only two currently available studies dealing with the RVA strain diversity in the COD. The first study analyzed the complete genomes of two G8 strains (RVA/Human-wt/COD/DRC86/2003/G8P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] and RVA/Human-wt/COD/DRC88/2003/G8P[8]). These strains both possess a DS-1-like genotype constellation and are the only complete genomes of Congolese RVA strains available today [13]. The second study described the epidemiology of RVA associated gastroenteritis in Kinshasa, COD, from July to October in 2003, 2004, and 2005. The predominant genotypes in 2003 were G8P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] and G8P[8] strains while in 2004 and 2005 G1P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] strains were predominant [14].

The number of studies conducted in sub-Saharan Africa involving the analyses of complete genomes of human RVA strains isolated in Africa is also limited. To our knowledge the complete genotype constellation of 87 human RVA strains isolated in sub-Saharan Africa are retrievable from Genbank: 28 of the 87 RVA strains were isolated in Malawi, 22 in Cameroon, 20 in South-Africa, 4 in Kenya, 3 in Nigeria, 3 in Cote d’Ivoire, 3 in Zambia, 2 in the COD, 1 in Zimbabwe and 1 in Ethiopia [11], [13], [34]–[38]. Compared to the hundreds of complete genomes characterized from RVAs isolated in high-income countries, and considering the large number of co-circulating G/P genotype combinations, this number is rather low.

The present study aimed to determine the genetic diversity of RVAs detected during 2007–2010 from Kisangani, COD. In addition we selected five P[6] RVA strains in combination with the different observed G-genotypes (G1, G4, 2xG8 and G12) for complete genome analyses to gain insight into the diversity, origin and evolution of P[6] RVA strains in the COD.

Methods

Stool sample collection and detection of RVA

Fecal specimens were collected between May 2007 and December 2010 from children <5 years old admitted with gastroenteritis (GE) to the pediatric wards of three hospitals in Kisangani (University Hospital, General Referral Hospital Makisao of Kisangani – both public hospitals – and Village de Pédiatrie, a private hospital center) and screened by an immunochromatographic antigen test (Rota-CIT, BioConcept, Belgium) for the presence of RVA antigen. Fecal samples containing RVA were collected and shipped to the Rega Institute for Medical Research using chromatography paper strips as previously described [39].

Rotavirus genotyping

The paper strips were inserted into an Eppendorf tube with 500 µl of universal transport medium (Copan Diagnostics, Corona, CA, USA) and thoroughly squeezed using sterile forceps. An aliquot of 140 µl of the squeezed eluate was used for RNA extraction using the QIAamp Viral RNA mini kit (QIAGEN/Westburg, Leusden, The Netherlands) according to the manufacturer’s instructions. The G- and P-genotypes of 99 antigen-positive samples were characterized by reverse-transcription polymerase chain reaction (RT-PCR) using the QIAGEN OneStep RT-PCR kit (QIAGEN/Westburg) using primers Beg9 and End9 (VP7), and VP4_1-17F and Con2 (VP4). Primer sequences are shown in table S1. PCR products were run on a polyacrylamide gel, stained with EtBr and visualized under UV-light. The positive samples were purified with ExoSAP-IT (Affymetrix, Santa Clara, CA, USA), and Sanger sequenced using the forward primers (Beg9 or VP4_1-17F) with the ABI PRISM BigDye Terminator cycle sequencing reaction kit (Applied Biosystems, Foster City, CA, USA).

Sequencing complete genomes

For the complete genomes, RT-PCR was carried out with an initial reverse transcription step at 50°C for 30 min followed by a PCR activation step at 95°C for 15 min, 35 cycles of amplification, and a final extension step for 10 min at 70°C in a Biometra T3000 thermocycler (Biometra, Westburg BV, Netherlands). For the gene segments encoding VP6, VP7, NSP1, NSP2, NSP3, NSP4 and NSP5, the amplification cycle conditions were as follows: 30 s at 94°C, 30 s at 45°C, and 2 min at 72°C. For the larger segments encoding VP1, VP2, VP3 and VP4 the cycle conditions were 30 s at 94°C, 30 s at 48°C, and 6 min at 72°C. Primers used to amplify the VP1, VP2, VP3, VP6, NSP1, NSP2, NSP3, NSP4 and NSP5 gene segments are described in table S1. Primer walking was used for full-length sequencing of the longest genes (VP1–VP4). The 5′ and 3′ terminal sequences of the 11 genomic segments were determined for all five strains using a modified RACE technique as previously described [40].

Sequence alignment and phylogenetic analysis

The Mega 5.10 software was used for phylogenetic and molecular evolutionary analyses [41]. To compute evolutionary distances we used the p-distance nucleotide substitution model. In phylogenetic analyses we opted for maximum likelihood phylogenetic analyses using the complete deletion function in order to include partially sequenced genes. Bootstrap resampling analysis (500 replicates) was performed to measure the reliability of the tree topologies.

Nucleotide sequence accession numbers

Accession numbers from sequences of this study: KJ870690-KJ870932.

Ethics Statement

Pending the installation of an ethical committee in the Oriental Province of the COD and at the Kisangani hospitals, the study was assessed and approved by the Provincial Health Officer (the highest ranked medical authority in the Oriental Province) and the Director’s Board of UH Kisangani, General Reference Hospital of Makiso and Village de Pédiatrie, respectively. All samples were collected in a standard of care setting for therapeutic or diagnostic reasons. There has been no influence on the treatment of the patients based on the results of the RVA detection test, and analyses for this study were performed after installment of the patients’ treatment. Data were fully anonymized at collection.

Results

1. Rotavirus incidence from May 2007 to December 2010

RVA detection rate and study data

A total of 524 fecal specimens were collected over the entire study period of which 154 samples (29%) tested positive with the immunochromatographic antigen test. Figure 1 shows the number of children hospitalized with gastroenteritis (GE) and the number of RVA positive cases over the entire study period. There were clear fluctuations in the number of RVA positive cases and in the number of hospitalized GE cases per month with a three-month long period (April–June 2008) during which no samples were collected because of stock rupture. Duration of hospitalizations associated with RVA varied between 1 and 10 days and occurred among children aged 1 week-2 years, with a mean age of 10 months.

Figure 1

Numbers of children hospitalized with GE (blue bars) and children hospitalized with GE tested RVA positive (red bars) per month from May 2007 until December 2010.

To determine a possible seasonal variability in RVA disease in Kisangani, all data shown in figure 1 were pooled per month and plotted together with the percentage of RVA positive cases in figure 2. The proportion of RVA positive cases among children hospitalized with GE ranged between 18.8% and 37.1%, with an average of 29.0%. RVA infections occurred year-round with three periods of increased RVA detection rates: February, April–May and September–November.

Figure 2

Number of hospitalized GE cases (blue bars) and number of hospitalized GE cases tested RVA positive (red bars) pooled per month over the entire study period.

The green line indicates the percentage RVA positive cases tested in each month.

2. Genotype distribution

G- and P- genotyping of rotavirus strains

Ninety-nine samples out of 154 samples were selected, based on the availability of sufficient stool, for further G and P-genotyping. In the analyzed samples the predominant G-type was G1 (detected in 39% of specimens) and the most predominant P-type was P[6] (53%). A total of eight different G/P-combinations were found: G1P[8] (N = 28), G8P[6] (N = 26), G2P[4] (N = 14), G12P[6] (N = 13), G1P[6] (N = 11), G9P[8] (N = 4), G4P[6] (N = 2) and G8P[4] (N = 1) (figure 3). Only 46% of the isolates were found with the most typical human RVA genotype combinations (G1P[8], G2P[4] and G9P[8]). In contrast, 53% of the strains were found with the P[6] genotype, in combination with the VP7 genotypes G8, G1, G12 or G4.

Figure 3

Geographic location of the Democratic Republic of the Congo on the African continent; location of the city of Kisangani in the COD together with the genotype distribution of the 99 analyzed RVA strains.

Percentages of the different genotypes are indicated on the left. The number of strains per genotype is shown on the right hand side.

To investigate the genetic diversity within the VP4 and VP7 genes isolated in this study phylogenetic trees based on the partial nucleotide sequences of both genes were constructed (figures 4–5). For VP4, the majority (53%) of the characterized RVA strains clustered within the P[6] genotype, more specifically in two previously defined lineages, P[6]-I and P[6]-V, both of them containing human and porcine strains [10]. The Congolese RVA strains belonging to P[6]-I lineage analyzed in this study can be divided in three different clusters, two large and one smaller cluster (figure 4). The first large cluster contains the majority (24/26) of G8P[6] strains while the second big cluster contains all G12P[6] strains. G1P[6] strains were found in both the G8P[6] as the G12P[6] clusters (figure 4). The third, smaller cluster contains only 2 G8P[6] strains: RVA/Human-wt/COD/KisB101/2007/G8P[6] and RVA/Human-wt/COD/KisB605/2008/G8P[6], the only two G8P[6] strains isolated before 2009. Both G4P[6] strains detected in this study (KisB332 and KisB503) belonged to the P[6]-V lineage, a lineage containing mostly porcine or human strains believed to be of (partial) porcine RVA origin. Both samples were isolated in 2008, only one month apart from each other in the same hospital and shared 99.7% nucleotide (nt) identity in their VP4 genes.

Figure 4

Maximum likelihood tree based on 496 nucleotides of the VP4 encoding gene of the 99 RVA strains characterized in this study.

Tree was rooted using RVA/Human-tc/JPN/AU-1/1982/G3P[9] as an outgroup. Bootstrap values (500 replicates) above 70 are shown. The strains completely analyzed in this study are marked with a black triangle.

Figure 5

Maximum likelihood tree based on 363 nucleotides sequence of the VP7 encoding gene of the 99 RVA strains characterized in this study.

Tree was rooted using RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17] as an outgroup. Bootstrap values (500 replicates) above 70 are shown. The strains completely analyzed in this study are marked with a black triangle.

All P[4] strains were found in combination with the G2 genotype, except for a single G8 strain (KisB566). All G2P[4] strains were very closely related to each other despite the fact that they were isolated over the entire study period. The fact that the VP4 gene of KisB566 was almost identical to the G2P[4] strains suggests a rather recent reassortment event involving the VP7 gene. The P[8] strains from Kisangani showed limited diversity, with almost identical G1P[8] strains circulating in the Kisangani area between 2007 and 2009. The only characterized G1P[8] strain isolated in 2010 (RVA/Human-wt/COD/KisB563/2010/G1P[8]) was found to cluster closer to the four detected G9P[8] strains than to the other G1P[8] strains characterized in this study, again suggesting a recent reassortment event.

For VP7, we found limited diversity within each G-genotype, especially within G2, G4, G9 and G12 were we found at least 97.9%, 99.2%, 98.7% and 98.2% of nt similarity, respectively. Within the G1 genotype, two clusters could be defined, one of which contains only a single Kisangani strain: RVA/Human-wt/COD/KisB563/2010/G1P[8] (figure 5). This strain is the only G1P[8] RVA strain isolated in 2010, and shows 94.7% –95.7% nt similarity with G1P[8] strains isolated between 2007–2009. The large G1 cluster showed strains with P[6] and P[8] intermingled, further indicating the frequent occurrence of reassortments. A clustering according to year of isolation can also be found within the G8 RVA strains. More specifically, strains isolated in 2007 or 2008 and strains isolated in 2009 or 2010 differ 3.6% to 4.7% on the nucleotide level.

Both for VP4 and VP7, all strains isolated in Kisangani were closely related to RVAs present in the global RVA data collection, in particular to strains isolated from other African patients. For more detail, extended trees for both VP4 and VP7 were added as supporting information (Figure S1 for VP4 and Figure S2 for VP7). Both trees contain twelve representative strains from this study together with RVA VP4 or VP7 sequences available in Genbank.

3. Complete genome analysis

Genotype constellations of whole genomes

To gain insight into the diversity, origin and evolution of P[6] RVA strains in the COD we selected P[6] strains, based on the phylogenetic clustering in the VP7 tree and on the viral load of the samples. In total, the nucleotide sequences of the eleven genes of five P[6] RVA strains with the different G-genotypes detected in this study (G1, G4, 2xG8 and G12) were completely determined, except from the first 33 nucleotides of the VP2 gene of strain KisB332 which could not be completed due to the lack of sample available to complete the analysis. Complete genotype assignments of the five fully sequenced P[6] strains together with a number of other representative strains are shown in table 1.

Table 1

RVA genotype constellations of the P[6] strains characterized in the present study (in bold) together with selected human and porcine P[6] RVA strains.

	VP7	VP4	VP6			VP1	VP2	VP3	NSP1	NSP2	NSP3	NSP4	NSP5
G1P[6] strains
RVA/Human-wt/COD/KisB504/2009/G1P [6]	G1	P [6]			I1	R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-wt/BRA/BA17290/2009/G1P[6]	G1	P[6]			I1	R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-tc/JPN/AU19/1997/G1P[6]	G1	P [6]			I5	R1	C1	M1	A8	N1	T1	E1	H2
G12P[6] strains
RVA/Human-wt/COD/KisB521/2008/G12P [6]	G12	P [6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-wt/ZAF/3176WC/2009/G12P[6]	G12	P[6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6]	G12	P[6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-wt/BGD/Dhaka12/2003/G12P[6]	G12	P[6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-wt/DEU/GER172/2008/G12P[6]	G12	P[6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
G4P[6] strains
RVA/Human-wt/COD/KisB332/2008/G4P [6]	G4	P [6]		I1		R1	C1	M1	A1	N1	T7	E1	H1
Hungarian strains BP190, BP1792 and BP1490	G4	P[6]		I1		R1	C1	M1	A1	N1	T7	E1	H1
Hungarian strains BP271, BP1125 and BP1547	G4	P[6]		I1		R1	C1	M1	A8	N1	T7	E1	H1
Hungarian strains BP1227 and BP1231	G4	P[6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-tc/GBR/ST3/1975/G4P[6]	G4	P[6]		I1		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Pig-tc/USA/Gottfried/1983/G4P[6]	G4	P[6]		I1		R1	C1	M1	A8	N1	T1	E1	H1
RVA/Human-wt/ARG/Arg4671/2006/G4P[6]	G4	P[6]		I1		R1	C1	M1	A8	N1	T1	E1	H1
RVA/Human-wt/ARG/Arg4605/2006/G4P[6]	G4	P [6]		I1		R1	C1	M1	A8	N1	T7	E1	H1
RVA/Human-wt/BRA/HSP180/1999/G4P[6]	G4	P [6]		I1		R1	C1	M1	A8	N1	T7	E1	H1
RVA/Human-wt/BRA/NB150/1997/G1G4P[6]	G1G4	P[6]		I1		R1	C1	M1	A1	Nx	T1	E1	H1
RVA/Human-tc/CHN/R479/2004/G4P[6]	G4	P[6]		I5		R1	C1	M1	A1	N1	T7	E1	H1
RVA/Human-wt/CHN/GX54/2010/G4P[6]	G4	P[6]		I1		R1	C1	M1	A8	N1	T1	E1	H1
RVA/Human-wt/IND/mani-362/2007/G4P[6]	G4	P[6]		I1		R1	C1	M1	A8	N1	T1	E1	H1
RVA/Human-wt/BRA/HST327/1999/G4P[6]	G4	P[6]		Ix		R1	C1	M1	A1	N1	T1	E1	H1
G8P[6] strains
RVA/Human-wt/COD/KisB554/2010/G8P [6]	G8	P [6]		I2		R2	C2	M2	A2	N2	T2	E2	H2
RVA/Human-wt/COD/KisB565/2010/G8P [6]	G8	P [6]		I2		R2	C2	M2	A2	N2	T2	E2	H2
RVA/Human-wt/COD/DRC86/2003/G8P[6]	G8	P[6]		I2		R2	C2	M2	A2	N2	T2	E2	H2
RVA/Human-tc/MWI/QEC29/2005/G8P[6]	G8	P [6]		I2		R2	C2	M2	A2	N2	T2	E2	H2
RVA/Human-tc/MWI/MW1-860/1999/G8P[6]	G8	P [6]		I2		R2	C2	M2	A2	N1	T2	E2	H1
RVA/Human-wt/CMR/6809/2000/G8P[6]	G8	P [6]		I2		R1	C1	M1	A1	N1	T1	E1	H1
RVA/Human-wt/CMR/6782/2000/G8P[6]	G8	P [6]		I1		R1	C1	M1	A1	N1	T1	E1	H1

Gene segments of which no complete open reading frame was available in Genbank are indicated in italic.

The genetic background of the G1P[6] and G12P[6] strains was entirely composed of genotype 1 (Wa-like). To our knowledge, no strain with the same genotype constellation as RVA/Human-wt/COD/KisB504/2009/G1P[6] has previously been reported from Africa. However, the same constellation was reported in a study investigating 6 Brazilian G1P[6] RVA strains, isolated from vaccinated children in 2009 and 2010 [42]. The G12-P[6]-I1-R1-C1-M1-A1-N1-T1-E1-H1 configuration has been previously reported in two African countries: South-Africa (RVA/Human-wt/ZAF/3176WC/2009/G12P[6]) and Zambia (RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6]) [35]. Both G8P[6] strains possessed a conserved genomic configuration: G8-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. This constellation was exactly the same as that of RVA/Human-wt/COD/DRC86/2003/G8P[6], also isolated in the COD, and RVA/Human-tc/MWI/QEC29/2005/G8P[6], isolated in Malawi [11], [13].

The G4P[6] strain was found to possess the, for humans unusual, G4-P[6]-I1-R1-C1-M1-A1-N1-T7-E1-H1 genotype constellation. Although unusual in humans, a recent study describing G4P[6] strains in Hungarian children, reported three strains with the exact same genotype constellation as KisB332 [43]. The genotype constellation I1/I5-R1-C1-M1-A1/A8-N1-T1/T7-E1-H1 is typically found in pigs with a variety of different G- and P-genotypes [9], [44]–[46]. A number of recent studies have identified human G4P[6] strains detected in different parts of the world, which are believed to be examples of direct interspecies transmissions, or may have undergone reassortments with human RVA strains [43], [47]–[51].

Phylogenetic analyses

Despite the different G-genotype, a very close genetic relationship exists between the G1P[6] (KisB504) and G12P[6] (KisB521) strain within their segments VP4 (99.6%), VP6 (99.9%), NSP1 (99.6%), NSP4 (99.8%) and NSP5 (99.8%). For the remaining genome segments encoding VP1, VP2, VP3, NSP2 and NSP3 a slightly lower nucleotide sequence similarity was detected, ranging from 94.4% (for NSP3) to 98.2% (for NSP2). The genetic similarity between the VP1-4, VP6, VP7, NSP1-3 and NSP5 of RVA strains KisB521, KisB504 and other recently isolated human Wa-like RVA strains from all over the world was very high (figure 4–8). The majority of gene segments of strain KisB521 (VP7, VP4, VP1-3, NSP1, NSP4) and five gene segments of strain KisB504 (VP4 and NSP1-4) are closely related to a recently isolated strain from Zambia (RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6]).

Figure 8

Maximum likelihood trees inferred from the nucleotide sequences of RVA NSP3, NSP4 and NSP5 genes.

The gene sequence lengths compared were: NSP3 (841 bp), NSP4 (501 bp), NSP5 (542 bp). The bootstrap values (500 replicates) are shown at the branch nodes (values <70% not shown). The scale bar is proportional to the genetic distance. An outgroup sequence (RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17]) for each gene was included. Completely sequenced RVA strains from Kisangani are indicated with black triangles.

Both G8P[6] strains clustered very closely together in the phylogenetic trees of all gene segments, with a nucleotide sequence similarity ranging between 99.5% and 100% (figure 4–8). These results indicate a clonal origin of both strains, which were isolated within a 2-month interval, from patients admitted to the same hospital. For VP7, VP3, NSP3 and NSP5 these strains shared the highest nucleotide similarity with RVA/Human-wt/COD/DRC86/2003/G8P[6]. Human RVA strain DRC86 is the only completely characterized Congolese G8P[6] strain available in Genbank until now [13]. Maximal nucleotide sequence identity of KisB554 and KisB565 was found with RVA strain RVA/Human-wt/USA/06-242/2006/G2P[6] (VP4, VP6, VP2, NSP1 and NSP2) and RVA/Human-wt/BEL/B1711/2002/G6P[6] (VP1). Strain 06–242 was isolated in an American child in 2006 and was the first G2P[6] strain characterized with a DS-1-like genotype constellation [52]. Strain B1711 also possessed a complete DS-1-like genotype constellation and was isolated from a 13-month-old child admitted to the university hospital in Leuven, Belgium, with severe gastroenteritis after returning from a trip to Mali [53], [6] rotavirus isolate are of bovine origin. J Gen Virol. 2008 ">54].

Overall, most gene segments of the G8P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6], G1P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] and G12P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] strains appeared to cluster closely to gene segments of recently isolated human RVA strains from all around the world, with the exception of NSP4. This gene segment showed more similarity with porcine/human porcine-like (for the G1P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6], G12P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] and G4P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] strains) or bovine/human bovine-like (for the G8P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] strains) RVAs (figure 8). More specific, the three Congolese strains (and the Zimbabwean strain RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6]) belonging to genotype E1 clustered together with RVA strains RVA/Human-wt/BRA/NB150/1997/G1G4P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] and RVA/Human-wt/ECU/EC2184/200x/G11P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6], both previously classified as porcine-human reassortants [51], [14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6], from Ecuador. Arch Virol. 2009 ">55]. Furthermore, these six strains also felt into a larger cluster containing porcine RVA strains, or strains believed to be of porcine origin. On the other hand, the NSP4 encoding gene segments of the G8P[14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6] strains belonged to an E2 genotype cluster containing both animal and human P[14] RVA strains, which are believed to be the result of direct interspecies transmissions from sheep or other members of the mammalian order Artiodactyla [14] rotavirus strains the result of interspecies transmissions from sheep or other ungulates that belong to the mammalian order Artiodactyla?. J Virol. 2009 ">6].

For the G4P[6] strain (KisB332) it appeared that all gene segments were more closely related to porcine RVA strains or human porcine-like RVA strains than to typical human RVA strains. A possible exception is VP7, because it is difficult to determine the exact origin of the VP7 gene, since there are no porcine G4 RVA sequences from Africa known today. In general, the gene segments of KisB332 showed lower nucleotide sequence identities to known contemporary and historic human or animal RVAs, especially in the VP6, VP1, VP2 and NSP1 genes (figures 6 and [7]). Although these gene segments of KisB332 were only relatively distantly related to known RVA strains, they showed the highest relatedness to RVA/Human-wt/NCA/24J/2010/G1P[8] (for VP6), RVA/Pig-tc/USA/OSU/1977/G5P9[7] (VP1), RVA/Pig-wt/KOR/K71/2006/G5P7 (VP2) and human reference strain RVA/Human-tc/USA/WI61/1983/G9P1A[8] (NSP1) [56], [57]. The VP3 of KisB332 was most closely related to porcine RVA strain OSU (figure 7). The VP4, NSP2, NSP3, NSP4 and NSP5 genes clustered together with putative porcine interspecies transmitted RVA strains or human-porcine reassortant strains (BE2001, Arg4671, NB150, EC2184) [47], [51], [6], from Ecuador. Arch Virol. 2009 ">55], [58]. Interestingly, the NSP3 gene segment of KisB332 belongs to the rarely described T7 genotype (figure 8). The T7 genotype was first detected in a bovine RVA strain (RVA/Cow-tc/GBR/UK/1973/G6P7[5]) and was more recently described in a bovine-like human strain (RVA/Human-xx/IND/Mani-265/2007/G10P[6]) [48],[59]. However, an increasing number of studies have reported the T7 genotype in porcine-like human strains (R479, BE2001, BP1901, BP1490, BP271, BP1792, BP1125, BP1547 and Arg4605) and porcine RVA strains (CMP45, AB82, F6–4, F8–4-A and TM-a) derived from different geographical locations [43], [46], [47], [50], [58]. The detection of the T7 genotype in strain RVA/Human-wt/COD/KisB332/2008/G4P[6], a strain with a human porcine-like gene constellation, makes it more likely that T7 genotypes have their origins in pigs as was previously suggested in other studies [46], [58].

Figure 6

Maximum likelihood trees inferred from the nucleotide sequences of RVA VP6, VP1 and VP2 genes.

The gene sequence lengths compared were: VP6 (1064 bp), VP1 (3168 bp), VP2 (2534 bp). The bootstrap values (500 replicates) are shown at the branch nodes (values <70% not shown). The scale bar is proportional to the genetic distance. An outgroup sequence (RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17]) for each gene was included. Completely sequenced RVA strains from Kisangani are indicated with black triangles.

Figure 7

Maximum likelihood trees inferred from the nucleotide sequences of RVA VP3, NSP1 and NSP2 genes.

The gene sequence lengths compared were: VP3 (2442 bp), NSP1 (1139 bp), NSP2 (938 bp). The bootstrap values (500 replicates) are shown at the branch nodes (values <70% not shown). The scale bar is proportional to the genetic distance. An outgroup sequence (RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17] for VP3 and NSP2 and RVA/Horse-wt/ARG/E403/2006/G14P[12] for NSP1) was included. Completely sequenced RVA strains from Kisangani are indicated with black triangles.

Discussion

Data shown in figure 2, suggest three periods of increased RVA detection rates over the course of this study. However, these fluctuations can be partially explained by sampling bias inherent to the circumstances in which these samples have been collected. Kisangani is known for its complex political and economical climate, making it difficult to collect sufficient samples over a large period of time without any period of lower or no sampling. Although the data available in this study do have their limitations, the available data confirm previous studies that reported year-round RVA infections in tropical settings, with ups and downs but no clear period with minimal or no RVA hospitalizations [60], [61].

The six most common G/P-genotype combinations in countries with a more moderate climate and by extension worldwide are G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8]. In this study a high proportion of P[6] RVA strains was observed together with a large number of unusual strains (G8P[6], G12P[6], G1P[6], G4P[6] and G8P[4]). These unusual genotype combinations may represent potential reassortants between human and animal RVA strains as a result of zoonotic infection. The only way to investigate this is to perform phylogenetic analysis of complete genomes. Because of the high number of P[6] RVA strains detected in this study and the fact that the P[6] genotype is much less prevalent in high-income countries this study focused on the four different P[6] genotype combinations detected in Kisangani between May 2007 and December 2010. The majority of the gene segments from the five P[6] strains characterized in this study appeared to have a genetic background most closely related to human RVA strains circulating all over the world. However, all five strains possessed at least one gene segment (NSP4) that was most closely related to animal RVAs. One strain, RVA/Human-wt/COD/KisB332/2008/G4P[6], possessed a majority of gene segments (if not all) most closely related to porcine, or human porcine-like RVA strains. For most genome segments KisB332 had still a considerable genetic distance to known porcine RVA strains. This can partly be explained by sampling bias since surveillance studies for animal rotaviruses are much less common than human rotavirus surveillance studies. In addition, porcine RVA sequences obtained from African pigs are not available in Genbank thus far. The lack of availability of complete genome sequences from porcine RVA strains in general and from African pig RVAs in particular makes it difficult to differentiate between the porcine or human origin of certain gene segments.

The child from which strain KisB332 was isolated lived in close proximity of a pig farm, providing a plausible route of transmission. Unfortunately it was not possible to take samples from the pigs living in close contact with the child. These findings suggest that KisB332 is of porcine origin and that this strain, after an interspecies transmission and potential reassortment with human RVAs, was able to cause gastro-enteritis in a human child. Surprisingly, a second G4P[6] strain (KisB603) was found in this study with very similar VP4 and VP7 gene sequences as those of KisB332. In theory, both strains can be the result of independent interspecies transmission events but a more plausible theory would be that the G4P[6] strain was able to spread from one human to another, although this hypothesis needs further investigation.

Although interspecies transmission events mostly result in so-called death-end infections, interspecies transmissions from pigs to humans seem to occur rather frequently, probably due to the fact that porcine RVA strains share a similar genetic background with human Wa-like RVAs [62]. Overall, these results underline the importance of complete genetic characterization of RVA strains and indicate that reassortments and interspecies transmission between human RVAs and porcine RVA strains occur frequently in the COD. Interspecies transmission events can have a potential negative influence on the efficacy of RVA vaccines, highlighting the need for continued surveillance of RVA diversity in the COD.

Based on the data from this study, RVA vaccines will be challenged with a wider variety of different G- and P- genotypes in the COD compared to those found in high-income countries. Of particular importance for the potential implementation of the currently available RVA vaccines in the COD is the predominance of P[6] detected in this study (52%) and the high proportion of strains with the G8 VP7 specificity (27%). Since these genotypes are not included in the currently available RVA vaccines and the genotype specific vaccine efficacy against these genotypes is not yet known, the success of vaccination programs could be lower. Nevertheless, implementation of rotavirus vaccines in the COD would be highly recommended to prevent rotavirus related morbidity and mortality.

Phylogenetic dendrogram based on 595 nucleotides of the VP4 encoding gene of the 12 representative RVA strains characterized in the study (black triangle) and RVA strains isolated all over the world. Tree was rooted using RVA/Human-tc/JPN/AU-1/1982/G3P9] as an outgroup. Bootstrap values (500 replicates) above 70 are shown.

(TIF)

Click here for additional data file.

Phylogenetic dendrogram based on 503 nucleotides of the VP7 encoding gene of the 12 representative RVA strains characterized in this study (black triangle) and RVA strains isolated all over the world. Tree was rooted using RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17] as an outgroup. Bootstrap values (500 replicates) above 70 are shown.

(TIF)

Click here for additional data file.

Primers used to amplify the VP1, VP2, VP3, VP4, VP6, VP7, NSP1, NSP2, NSP3, NSP4 and NSP5 gene segments described in this study.

(DOC)

Click here for additional data file.