Whole Genome Analysis of Selected Human Group A Rotavirus Strains Revealed Evolution of DS-1-Like Single- and Double-Gene Reassortant Rotavirus Strains in Pakistan During 2015-2016

Characterization of Rotavirus Strains Responsible for Breakthrough Diarrheal Diseases among Zambian Children Using Whole Genome Sequencing

The occurrence of rotavirus (RV) infection among vaccinated children in high-burden settings poses a threat to further disease burden reduction. Genetically altered viruses have the potential to evade both natural infection and vaccine-induced immune responses, leading to diarrheal diseases among vaccinated children. Studies characterizing RV strains responsible for breakthrough infections in resource-limited countries where RV-associated diarrheal diseases are endemic are limited. We aimed to characterize RV strains detected in fully vaccinated children residing in Zambia using next-generation sequencing. We conducted whole genome sequencing on Illumina MiSeq. Whole genome assembly was performed using Geneious Prime 2023.1.2. A total of 76 diarrheal stool specimens were screened for RV, and 4/76 (5.2%) were RV-positive. Whole genome analysis revealed RVA/Human-wt/ZMB/CIDRZ-RV2088/2020/G1P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2 and RVA/Human-wt/ZMB/CIDRZ-RV2106/2020/G12P[4]-I1-R2-C2-M2-A2-N1-T2-E1-H2 strains were mono and multiple reassortant (exchanged genes in bold) respectively, whilst RVA/Human-wt/ZMB/CIDRZ-RV2150/2020/G12P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 was a typical Wa-like strain. Comparison of VP7 and VP4 antigenic epitope of breakthrough strains and Rotarix strain revealed several amino acid differences. Variations in amino acids in antigenic epitope suggested they played a role in immune evasion of neutralizing antibodies elicited by vaccination. Findings from this study have the potential to inform national RV vaccination strategies and the design of highly efficacious universal RV vaccines.

Comparative whole genome analysis reveals re-emergence of human Wa-like and DS-1-like G3 rotaviruses after Rotarix vaccine introduction in Malawi

G3 rotaviruses rank among the most common rotavirus strains worldwide in humans and animals. However, despite a robust long-term rotavirus surveillance system from 1997 at Queen Elizabeth Central Hospital in Blantyre, Malawi, these strains were only detected from 1997 to 1999 and then disappeared and re-emerged in 2017, 5 years after the introduction of the Rotarix rotavirus vaccine. Here, we analysed representative twenty-seven whole genome sequences (G3P[4], n = 20; G3P[6], n = 1; and G3P[8], n = 6) randomly selected each month between November 2017 and August 2019 to understand how G3 strains re-emerged in Malawi. We found four genotype constellations that were associated with the emergent G3 strains and co-circulated in Malawi post-Rotarix vaccine introduction: G3P[4] and G3P[6] strains with the DS-1-like genetic backbone genes (G3-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2 and G3-P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2), G3P[8] strains with the Wa-like genetic backbone genes (G3-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1), and reassortant G3P[4] strains consisting of the DS-1-like genetic backbone genes and a Wa-like NSP2 (N1) gene (G3-P[4]-I2-R2-C2-M2-A2-N1-T2-E2-H2). Time-resolved phylogenetic trees demonstrated that the most recent common ancestor for each ribonucleic acid (RNA) segment of the emergent G3 strains was between 1996 and 2012, possibly through introductions from outside the country due to the limited genetic similarity with G3 strains which circulated before their disappearance in the late 1990s. Further genomic analysis revealed that the reassortant DS-1-like G3P[4] strains acquired a Wa-like NSP2 genome segment (N1 genotype) through intergenogroup reassortment; an artiodactyl-like VP3 through intergenogroup interspecies reassortment; and VP6, NSP1, and NSP4 segments through intragenogroup reassortment likely before importation into Malawi. Additionally, the emergent G3 strains contain amino acid substitutions within the antigenic regions of the VP4 proteins which could potentially impact the binding of rotavirus vaccine-induced antibodies. Altogether, our findings show that multiple strains with either Wa-like or DS-1-like genotype constellations have driven the re-emergence of G3 strains. The findings also highlight the role of human mobility and genome reassortment events in the cross-border dissemination and evolution of rotavirus strains in Malawi necessitating the need for long-term genomic surveillance of rotavirus in high disease-burden settings to inform disease prevention and control.

Epidemic and Evolutionary Characteristics of Swine Enteric Viruses in South-Central China from 2018 to 2021

Swine enteric viruses are a major cause of piglet diarrhea, causing a devastating impact on the pork industry. To further understand the molecular epidemiology and evolutionary diversity of swine enteric viruses, we carried out a molecular epidemiological investigation of swine enteric viruses (PEDV, PDCoV, PoRVA, and TGEV) on 7107 samples collected from pig farms in south-central China. The results demonstrated that PEDV is the predominant pathogen causing piglet diarrhea, and its infection occurs mainly in relatively cold winter and spring in Hunan and Hubei provinces. The positive rate of PEDV showed an abnormal increase from 2020 to 2021, and that of PoRVA and PDCoV exhibited gradual increases from 2018 to 2021. PEDV-PoRVA and PEDV-PDCoV were the dominant co-infection modes. A genetic evolution analysis based on the PEDV S1 gene and ORF3 gene revealed that the PEDV GII-a is currently epidemic genotype, and the ORF3 gene of DY2020 belongs to a different clade relative to other GII-a strains isolated in this study. Overall, our results indicated that the variant PEDV GII-a is the main pathogen of piglet diarrhea with a trend of outbreak. G9 is the dominant PoRVA genotype and has the possibility of outbreak as well. It is therefore critical to strengthen the surveillance of PEDV and PoRVA, and to provide technical reserves for the prevention and control of piglet diarrhea.

Effect of rotavirus genetic diversity on vaccine impact

Group A rotaviruses (RVAs) are the leading cause of gastroenteritis, causing 0.2 million deaths and several million hospitalisations globally each year. Four rotavirus vaccines (Rotarix^TM , RotaTeq^TM , Rotavac^® and ROTASIIL^® ) have been pre-qualified by the World Health Organization (WHO), but the two newly pre-qualified vaccines (Rotavac^® and ROTASIIL^® ) are currently only in use in Palestine and India, respectively. In 2009, WHO strongly proposed that rotavirus vaccines be included in the routine vaccination schedule of all countries around the world. By the end of 2019, a total of 108 countries had administered rotavirus vaccines, and 10 countries have currently been approved by Gavi for the introduction of rotavirus vaccine in the near future. With 39% of global coverage, rotavirus vaccines have had a substantial effect on diarrhoeal morbidity and mortality in different geographical areas, although efficacy appears to be higher in high income settings. Due to the segmented RNA genome, the pattern of RVA genotypes in the human population is evolving through interspecies transmission and/or reassortment events for which the vaccine might be less effective in the future. However, despite the relative increase in some particular genotypes after rotavirus vaccine use, the overall efficacy of rotavirus mass vaccination worldwide has not been affected. Some of the challenges to improve the effect of current rotavirus vaccines can be solved in the future by new rotavirus vaccines and by vaccines currently in progress.

Whole genome analysis of Aichivirus A, isolated from a child, suffering from gastroenteritis, in Pakistan

Viruses are the primary cause of acute gastroenteritis in children all over the world. Understanding the emergence and genetic variation of these viruses may help to prevent infections. Aichivirus (AiV) is a member of the Kobuvirus genus, which currently contains six officially recognized species: Aichivirus A-F. The species AiV A contains six types including Aichivirus 1 (AiV 1) and eventually, three genotypes have been identified in the human AiV 1 (named A to C). The present study describes the identification and sequencing of the polyprotein gene of a human AiV 1 strain PAK419 via NGS in Pakistani children with acute gastroenteritis. Our study strain PAK419 was classified as AiV 1 genotype A, most commonly found in Japan and Europe, and closely related to non-Japanese and European strains on the phylogenetic tree. PAK419 showed 95-98 % nucleotide sequence identity with strains isolated from Ethiopia (ETH/2016/P4), Australia (FSS693) and China (Chshc7). On phylogenetic observation PAK419 formed a distinct cluster in the AiV 1 genotype A with the above mentioned and other human AiV strains detected around the world (Germany, Brazil, Japan, Thailand, Korea and Vietnam). The data clearly showed that Pakistani AiV strains and human strains identified from all over the world are distinct from Aichivirus strains found in bovine, swine, canine, feline, caprine, ferret, bat, and environmental samples. The distinguishing characteristics of the AiV genome showed a lower probability of inter-genotypic recombination events, which may support the lack of AiV serotypes. PAK419 also had a high content of C nucleotide (37.4 %), as found in previous studies, which could also restrict the possible genetic variation of AiV. This study demonstrate the power of NGS in uncovering unknown gastroenteric etiological agents circulating in the population.

Molecular characterization of the complete genome sequence of human Parechovirus 1 in Pakistan

Human parechoviruses (HPeVs) are highly common pathogens in children under 2 years of age. Of the 19 distinct HPeV genotypes identified worldwide, HPeV1 is still the most prevalent type associated with respiratory and gastrointestinal symptoms in infants and young children. Pakistan's previous studies have focused only on the detection and partial sequencing of HPeV genotypes. In the present study, we have obtained the complete genomes of 2 HPeV1 strains (PAK419 and PAK663) from children using NGS method on Illumina Hiseq Platform. These samples were collected from children suffering from acute gastroenteritis in Rawalpindi, Pakistan during 2016. The near complete genome sequences obtained for two HPeV1 strains (PAK419 and PAK663) consist of total 6877 nucleotides with a single, large open reading frame (ORF) encoding a polyprotein gene. Phylogenetic analysis showed that both HPeV1 strains exhibited maximum amino acid similarity (97 %) to HPeV1 strains from The Nederlands (2007-863, GQ183034) and clustered closely with this and with other HPeV1 strains isolated from other countries in the world (Ethiopia, Taiwan, Russia and Brazil). A motif of arginine-glycine-aspartic acid (RGD) in the VP1 (Outer capsid protein) C-terminus region that is suggested to help virus entry into the host cell also identified in PAK419 and PAK663. SimPlot analysis revealed that intergenotypic recombination events may have take place in the non-structural region between both HPeV1 strains (PAK419, PAK663), two major strains of HPeV1 (GQ183034 and MG873157) and four minor strains of HPeV4 (AM235750), HPeV7 (EU556224), HPeV15 (MN265386) and HPeV18 (KT879915). The full genome of HPeV1 strains characterized in the current study will provide complete information on these newly isolated strains for further preventive or treatment measures.

Full genome characterization of novel DS-1-like G9P[8] rotavirus strains that have emerged in Thailand

The emergence and rapid spread of unusual DS-1-like intergenogroup reassortant rotaviruses having G1/3/8 genotypes have been recently reported from major parts of the world (Africa, Asia, Australia, Europe, and the Americas). During rotavirus surveillance in Thailand, three novel intergenogroup reassortant strains possessing the G9P[8] genotype (DBM2017-016, DBM2017-203, and DBM2018-291) were identified in three stool specimens from diarrheic children. In the present study, we determined and analyzed the full genomes of these three strains. On full-genomic analysis, all three strains were found to share a unique genotype constellation comprising both genogroup 1 and 2 genes: G9-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Phylogenetic analysis demonstrated that each of the 11 genes of the three strains was closely related to that of emerging DS-1-like intergenogroup reassortant, human, and/or locally circulating human strains. Thus, the three strains were suggested to be multiple reassortants that had acquired the G9-VP7 genes from co-circulating Wa-like G9P[8] rotaviruses in the genetic background of DS-1-like intergenogroup reassortant (likely equine-like G3P[8]) strains. To our knowledge, this is the first description of emerging DS-1-like intergenogroup reassortant strains having the G9P[8] genotype. Our observations will add to the growing insights into the dynamic evolution of emerging DS-1-like intergenogroup reassortant rotaviruses through reassortment.