Sub-genotype phylogeny of the non-G, non-P genes of genotype 2 Rotavirus A strains

Full genome sequence analysis of the predominant and uncommon G9P[4] rotavirus strains circulating in Tehran, Iran, 2021-2022: Evidence for inter and intra-genotype recombination

Group A rotaviruses (RVAs) are a major cause of acute gastroenteritis in children under 5 years of age worldwide. Herein, the genetic sequences of 11 RNA segments from three uncommon G9P[4] RVA strains found in the stool samples of children under 5 years of age in Iran were analyzed using next-generation sequencing (NGS) technology. The genomic constellations of these three uncommon G9P[4] strains indicated the presence of the double and quadruple reassortants of two G9P[4] strains, containing the VP7/NSP2 and VP7/VP2/NSP2/NSP4 genes on a DS-1-like genetic background, respectively. The genome of one strain indicated a Wa-like genetic backbone in a single-reassortant with the VP4 of the DS1-like human strains. With the exception of VP1, VP2, VP7, NSP2, NSP3, and NSP4 genes, which clustered with RVA of human origins belonging to cognate gene sequences of genogroup 1/2 genotypes/lineages, the remaining five genes (VP8/VP4, VP3, VP6, NSP1, NSP5) displayed direct evidence of recombination. It is presumed that the presence of uncommon G9P[4] strains in Iran is not linked to vaccination pressure, but rather to the high prevalence of RVA co-infection or the direct import of these uncommon RVA reassortants strains from other countries (especially those that have implemented RV vaccination).

Whole-Genome Characterization of Rotavirus G9P[6] and G9P[4] Strains That Emerged after Rotavirus Vaccine Introduction in Mozambique

Mozambique introduced the Rotarix® vaccine into the National Immunization Program in September 2015. Following vaccine introduction, rotavirus A (RVA) genotypes, G9P[4] and G9P[6], were detected for the first time since rotavirus surveillance programs were implemented in the country. To understand the emergence of these strains, the whole genomes of 47 ELISA RVA positive strains detected between 2015 and 2018 were characterized using an Illumina MiSeq-based sequencing pipeline. Of the 29 G9 strains characterized, 14 exhibited a typical Wa-like genome constellation and 15 a DS-1-like genome constellation. Mostly, the G9P[4] and G9P[6] strains clustered consistently for most of the genome segments, except the G- and P-genotypes. For the G9 genotype, the strains formed three different conserved clades, separated by the P type (P[4], P[6] and P[8]), suggesting different origins for this genotype. Analysis of the VP6-encoding gene revealed that seven G9P[6] strains clustered close to antelope and bovine strains. A rare E6 NSP4 genotype was detected for strain RVA/Human-wt/MOZ/HCN1595/2017/G9P[4] and a genetically distinct lineage IV or OP354-like P[8] was identified for RVA/Human-wt/MOZ/HGJM0644/2015/G9P[8] strain. These results highlight the need for genomic surveillance of RVA strains detected in Mozambique and the importance of following a One Health approach to identify and characterize potential zoonotic strains causing acute gastroenteritis in Mozambican children.

Evolution of DS-1-like G8P[8] rotavirus A strains from Vietnamese children with acute gastroenteritis (2014-21): Adaptation and loss of animal rotavirus-derived genes during human-to-human spread

Animal rotaviruses A (RVAs) are considered the source of emerging, novel RVA strains that have the potential to cause global spread in humans. A case in point was the emergence of G8 bovine RVA consisting of the P[8] VP4 gene and the DS-1-like backbone genes that appeared to have jumped into humans recently. However, it was not well documented what evolutionary changes occurred on the animal RVA-derived genes during circulation in humans. Rotavirus surveillance in Vietnam found that DS-1-like G8P[8] strains emerged in 2014, circulated in two prevalent waves, and disappeared in 2021. This surveillance provided us with a unique opportunity to investigate the whole process of evolutionary changes, which occurred in an animal RVA that had jumped the host species barrier. Of the 843 G8P[8] samples collected from children with acute diarrhoea in Vietnam between 2014 and 2021, fifty-eight strains were selected based on their distinctive electropherotypes of the genomic RNA identified using polyacrylamide gel electrophoresis. Whole-genome sequence analysis of those fifty-eight strains showed that the strains dominant during the first wave of prevalence (2014-17) carried animal RVA-derived VP1, NSP2, and NSP4 genes. However, the strains from the second wave of prevalence (2018-21) lost these genes, which were replaced with cognate human RVA-derived genes, thus creating strain with G8P[8] on a fully DS-1-like human RVA gene backbone. The G8 VP7 and P[8] VP4 genes underwent some point mutations but the phylogenetic lineages to which they belonged remained unchanged. We, therefore, propose a hypothesis regarding the tendency for the animal RVA-derived genes to be expelled from the backbone genes of the progeny strains after crossing the host species barrier. This study underlines the importance of long-term surveillance of circulating wild-type strains in order to better understand the adaptation process and the fate of newly emerging, animal-derived RVA among the human population. Further studies are warranted to disclose the molecular mechanisms by which spillover animal RVAs become readily transmissible among humans, and the roles played by the expulsion of animal-derived genes and herd immunity formed in the local population.

Genetic diversity and evolution of G12P[6] DS-1-like and G12P[9] AU-1-like Rotavirus strains in Brazil

In the early 2000s, the global emergence of rotavirus (RVA) G12P[8] genotype was noted, while G12P[6] and G12P[9] combinations remained rare in humans. This study aimed to characterize and phylogenetically analyze three Brazilian G12P[9] and four G12P[6] RVA strains from 2011 to 2020, through RT-PCR and sequencing, in order to enhance our understanding of the genetic relationship between human and animal-origin RVA strains. G12P[6] strains displayed a DS-1-like backbone, showing a distinct genetic clustering. G12P[6] IAL-R52/2020, IAL-R95/2020 and IAL-R465/2019 strains clustered with 2019 Northeastern G12P[6] Brazilian strains and a 2018 Benin strain, whereas IAL-R86/2011 strain grouped with 2010 Northern G12P[6] Brazilian strains and G2P[4] strains from the United States and Belgium. These findings suggest an African genetic ancestry and reassortments with co-circulating American strains sharing the same DS-1-like constellation. No recent zoonotic reassortment was observed, and the DS-1-like constellation detected in Brazilian G12P[6] strains does not seem to be genetically linked to globally reported intergenogroup G1/G3/G9/G8P[8] DS-1-like human strains. G12P[9] strains exhibited an AU-1-like backbone with two different genotype-lineage constellations: IAL-R566/2011 and IAL-R1151/2012 belonged to a VP3/M3.V Lineage, and IAL-R870/2013 to a VP3/M3.II Lineage, suggesting two co-circulating strains in Brazil. This genetic diversity is not observed elsewhere, and the VP3/M3.II Lineage in G12P[9] strains seems to be exclusive to Brazil, indicating its evolution within the country. All three G12P[9] AU-1-like strains were closely relate to G12P[9] strains from Paraguay (2006-2007) and Brazil (2010). Phylogenetic analysis also highlighted that all South American G12P[9] AU-1-like strains had a common origin and supports the hypothesis of their importation from Asia, with no recent introduction from globally circulating G12P[9] strains or reassortments with local G12 strains P[8] or P[6]. Notably, certain genes in the Brazilian G12P[9] AU-1-like strains share ancestry with feline/canine RVAs (VP3/M3.II, NSP4/E3.IV and NSP2/N3.II), whereas NSP1/A3.VI likely originated from artiodactyls, suggesting a history of zoonotic transmission with human strains. This genomic data adds understanding to the molecular epidemiology of G12P[6] and G12P[9] RVA strains in Brazil, offering insights into their genetic diversity and evolution.

Emergence of Intergenogroup Reassortant G9P[4] Strains Following Rotavirus Vaccine Introduction in Ghana

Rotavirus (RVA) is a leading cause of childhood gastroenteritis. RVA vaccines have reduced the global disease burden; however, the emergence of intergenogroup reassortant strains is a growing concern. During surveillance in Ghana, we observed the emergence of G9P[4] RVA strains in the fourth year after RVA vaccine introduction. To investigate whether Ghanaian G9P[4] strains also exhibited the DS-1-like backbone, as seen in reassortant G1/G3/G8/G9 strains found in other countries in recent years, this study determined the whole genome sequences of fifteen G9P[4] and two G2P[4] RVA strains detected during 2015-2016. The results reveal that the Ghanaian G9P[4] strains exhibited a double-reassortant genotype, with G9-VP7 and E6-NSP4 genes on a DS-1-like backbone (G9-P[4]-I2-R2-C2-M2-A2-N2-T2-E6-H2). Although they shared a common ancestor with G9P[4] DS-1-like strains from other countries, further intra-reassortment events were observed among the original G9P[4] and co-circulating strains in Ghana. In the post-vaccine era, there were significant changes in the distribution of RVA genotype constellations, with unique strains emerging, indicating an impact beyond natural cyclical fluctuations. However, reassortant strains may exhibit instability and have a limited duration of appearance. Current vaccines have shown efficacy against DS-1-like strains; however, ongoing surveillance in fully vaccinated children is crucial for addressing concerns about long-term effectiveness.

Genomic Analysis of Rwandan G9P[8] Rotavirus Strains Pre- and Post-RotaTeq® Vaccine Reveals Significant Distinct Sub-Clustering in a Post-Vaccination Cohort

Although the introduction of rotavirus vaccines has substantially contributed to the reduction in rotavirus morbidity and mortality, concerns persist about the re-emergence of variant strains that might alter vaccine effectiveness in the long term. The G9 strains re-emerged in Africa during the mid-1990s and have more recently become predominant in some countries, such as Ghana and Zambia. In Rwanda, during the 2011 to 2015 routine surveillance period, G9P[8] persisted during both the pre- and post-vaccine periods. The pre-vaccination cohort was based on the surveillance period of 2011 to 2012, and the post-vaccination cohort was based on the period of 2013 to 2015, excluding 2014. The RotaTeq® vaccine that was first introduced in Rwanda in 2012 is genotypically heterologous to Viral Protein 7 (VP7) G9. This study elucidated the whole genome of Rwandan G9P[8] rotavirus strains pre- and post-RotaTeq® vaccine introduction. Fecal samples from Rwandan children under the age of five years (pre-vaccine n = 23; post-vaccine n = 7), conventionally genotyped and identified as G9P[8], were included. Whole-genome sequencing was then performed using the Illumina® MiSeq platform. Phylogenetic analysis and pair-wise sequence analysis were performed using MEGA6 software. Distinct clustering of three post-vaccination study strains was observed in all 11 gene segments, compared to the other Rwandan G9P[8] study strains. Specific amino acid differences were identified across the gene segments of these three 2015 post-vaccine strains. Important amino acid differences were identified at position N242S in the VP7 genome segment of the three post-vaccine G9 strains compared to the other G9 strains. This substitution occurs at a neutralization epitope site and may slightly affect protein interaction at that position. These findings indicate that the Rwandan G9P[8] strains revealed a distinct sub-clustering pattern among post-vaccination study strains circulating in Rwanda, with changes at neutralization epitopes, which may play a role in neutralization escape from vaccine candidates. This emphasizes the need for continuous whole-genome surveillance to better understand the evolution and epidemiology of the G9P[8] strains post-vaccination.

Detection and full-genotype determination of rare and reassortant rotavirus A strains in Nizhny Novgorod in the European part of Russia

Reassortant DS-1-like rotavirus A strains have been shown to circulate widely in many countries around the world. In Russia, the prevalence of such strains remains unclear due to the preferred use of the traditional binary classification system. In this work, we obtained partial sequence data from all 11 genome segments and determined the full-genotype constellations of rare and reassortant rotaviruses circulating in Nizhny Novgorod in 2016-2019. DS-1-like G3P[8] and G8P[8] strains were found, reflecting the global trend. Most likely, these strains were introduced into the territory of Russia from other countries but subsequently underwent further evolutionary changes locally. G3P[8], G9P[8], and G12P[8] Wa-like strains of subgenotypic lineages that are unusual for the territory of Russia were also identified. Reassortant G2P[8], G4P[4], and G9P[4] strains with one Wa-like gene (VP4 or VP7) on a DS-1-like backbone were found, and these apparently had a local origin. Feline-like G3P[9] and G6P[9] strains were found to be phylogenetically close to BA222 isolated from a cat in Italy but carried some traces of reassortment with human strains from Russia and other countries. Thus, full-genotype determination of rotavirus A strains in Nizhny Novgorod has clarified some questions related to their origin and evolution.

Full genotype characterization of Brazilian canine G3P[3] strains during a 10-year survey (2012-2021) of rotavirus infection in domestic dogs and cats

There is a dearth of information on the molecular epidemiology of rotaviruses in pets in Brazil. The aim of this study was to monitor rotavirus infections in household dogs and cats, determine full-genotype constellations, and obtain data on evolutionary relationships. Between 2012 and 2021, 600 fecal samples from dogs and cats (516 and 84, respectively) were collected at small animal clinics in São Paulo state, Brazil. Rotavirus screening was conducted using ELISA, PAGE, RT-PCR, sequencing, and phylogenetic analysis. Rotavirus type A (RVA) was detected in 0.5% (3/600) of the animals. No non-RVA types were detected. The three canine RVA strains were found to have a novel genetic constellation, G3-P[3] -I2-R3-C2-M3-A9-N2-T3-E3-H6, which has never been reported in dogs. As expected, all of the viral genes, except those encoding NSP2 and VP7, were closely related to the corresponding genes from canine, feline, and canine-like-human RVA strains. A novel N2 (NSP2) lineage was identified, grouping together Brazilian canine, human, rat and bovine strains, suggesting that genetic reassortment had occurred. Uruguayan G3 strains obtained from sewage contained VP7 genes that were phylogenetically close to those of the Brazilian canine strains, which suggests that these strains are widely distributed in pet populations in South American countries. For the NSP2 (I2), NSP3 (T3), NSP4 (E3), NSP5 (H6), VP1 (R3), VP3 (M3), and VP6 (I2) segments, phylogenetic analysis revealed possibly new lineages. The epidemiological and genetic data presented here point out the necessity for collaborative efforts to implement the One Health strategy in the field of RVA research and to provide an updated understanding of RVA strains circulating canines in Brazil.

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.

The Evolution of Post-Vaccine G8P[4] Group a Rotavirus Strains in Rwanda; Notable Variance at the Neutralization Epitope Sites

Africa has a high level of genetic diversity of rotavirus strains, which is suggested to be a possible reason contributing to the suboptimal effectiveness of rotavirus vaccines in this region. One strain that contributes to this rotavirus diversity in Africa is the G8P[4]. This study aimed to elucidate the entire genome and evolution of Rwandan G8P[4] strains. Illumina sequencing was performed for twenty-one Rwandan G8P[4] rotavirus strains. Twenty of the Rwandan G8P[4] strains had a pure DS-1-like genotype constellation, and one strain had a reassortant genotype constellation. Notable radical amino acid differences were observed at the neutralization sites when compared with cognate regions in vaccine strains potentially playing a role in neutralization escape. Phylogenetic analysis revealed that the closest relationship was with East African human group A rotavirus (RVA) strains for five of the genome segments. Two genome sequences of the NSP4 genome segment were closely related to bovine members of the DS-1-like family. Fourteen VP1 and eleven VP3 sequences had the closest relationships with the RotaTeq™ vaccine WC3 bovine genes. These findings suggest that the evolution of VP1 and VP3 might have resulted from reassortment events with RotaTeq™ vaccine WC3 bovine genes. The close phylogenetic relationship with East African G8P[4] strains from Kenya and Uganda suggests co-circulation in these countries. These findings highlight the need for continued whole-genomic surveillance to elucidate the evolution of G8P[4] strains, especially after the introduction of rotavirus vaccination.

Genomic Constellation of Human Rotavirus G8 Strains in Brazil over a 13-Year Period: Detection of the Novel Bovine-like G8P[8] Strains with the DS-1-like Backbone

Rotavirus (RVA) G8 is frequently detected in animals, but only occasionally in humans. G8 strains, however, are frequently documented in nations in Africa. Recently, an increase in G8 detection was observed outside Africa. The aims of the study were to monitor G8 infections in the Brazilian human population between 2007 and 2020, undertake the full-genotype characterization of the four G8P[4], six G8P[6] and two G8P[8] RVA strains and conduct phylogenetic analysis in order to understand their genetic diversity and evolution. A total of 12,978 specimens were screened for RVA using ELISA, PAGE, RT-PCR and Sanger sequencing. G8 genotype represented 0.6% (15/2434) of the entirely RVA-positive samples. G8P[4] comprised 33.3% (5/15), G8P[6] 46.7% (7/15) and G8P[8] 20% (3/15). All G8 strains showed a short RNA pattern. All twelve selected G8 strains displayed a DS-1-like genetic backbone. The whole-genotype analysis on a DS-1-like backbone identified four different genotype-linage constellations. According to VP7 analysis, the Brazilian G8P[8] strains with the DS-1-like backbone strains were derived from cattle and clustered with newly DS-1-like G1/G3/G9/G8P[8] strains and G2P[4] strains. Brazilian IAL-R193/2017/G8P[8] belonged to a VP1/R2.XI lineage and were grouped with bovine-like G8P[8] strains with the DS-1-like backbone strains detected in Asia. Otherwise, the Brazilian IAL-R558/2017/G8P[8] possess a "Distinct" VP1/R2 lineage never previously described and grouped apart from any of the DS-1-like reference strains. Collectively, our findings suggest that the Brazilian bovine-like G8P[8] strains with the DS-1-like backbone strains are continuously evolving and likely reassorting with local RVA strains rather than directly relating to imports from Asia. The Brazilian G8P[6]-DS-1-like strains have been reassorted with nearby co-circulating American strains of the same DS-1 genotype constellation. However, phylogenetic analyses revealed that these strains have some genetic origin from Africa. Finally, rather than being African-born, Brazilian G8P[4]-DS-1-like strains were likely imported from Europe. None of the Brazilian G8 strains examined here exhibited signs of recent zoonotic reassortment. G8 strains continued to be found in Brazil according to their intermittent and localized pattern, thus, does not suggest that a potential emergence is taking place in the country. Our research demonstrates the diversity of G8 RVA strains in Brazil and adds to the understanding of G8P[4]/P[6]/P[8] RVA genetic diversity and evolution on a global scale.

Genomic Analysis of G2P[4] Group A Rotaviruses in Zambia Reveals Positive Selection in Amino Acid Site 7 of Viral Protein 3

The G2P[4] genotype is among the rotavirus strains that circulate commonly in humans. Several countries have reported its immediate upsurge after the introduction of rotavirus vaccination, raising concern about sub-optimal vaccine effectiveness against this genotype in the long term. This study aimed to gain insight into the evolution of post-vaccine Zambian G2P[4] group A rotavirus (RVA) strains and their overall genetic make-up by analysis of sequence alignments at the amino acid (AA) level. Twenty-nine Zambian G2P[4] rotavirus strains were subjected to whole-genome sequencing using the Illumina MiSeq® platform. All the strains exhibited the typical DS-1-like genotype constellation, and the nucleotide sequences of the 11 genome segments showed high nucleotide similarities (>97%). Phylogenetic analyses together with representative global G2P[4] RVA showed that Zambian strains clustered into human lineages IV (for VP2, VP4, VP7, NSP1, and NSP5), V (for VP1, VP3, VP6, NSP2, and NSP3), and XXIII (for NSP4). The AA differences between the lineages where the study strains clustered and lineages of global reference strains were identified and analyzed. Selection pressure analysis revealed that AA site seven in the Viral Protein 3 (VP3) genome segment was under positive selection. This site occurs in the region of intrinsic disorder in the VP3 protein, and Zambian G2P[4] strains could potentially be utilizing this intrinsically disordered region to survive immune pressure. The Zambian G2P[4] strains from 2012 to 2016 comprised the G2P[4] strains that have been circulating globally since the early 2000s, highlighting the epidemiological fitness of these contemporary G2P[4] strains. Continuous whole-genome surveillance of G2P[4] strains remains imperative to understand their evolution during the post-vaccination period.

Evolutionary changes between pre- and post-vaccine South African group A G2P[4] rotavirus strains, 2003-2017

The transient upsurge of G2P[4] group A rotavirus (RVA) after Rotarix vaccine introduction in several countries has been a matter of concern. To gain insight into the diversity and evolution of G2P[4] strains in South Africa pre- and post-RVA vaccination introduction, whole-genome sequencing was performed for RVA positive faecal specimens collected between 2003 and 2017 and samples previously sequenced were obtained from GenBank (n=103; 56 pre- and 47 post-vaccine). Pre-vaccine G2 sequences predominantly clustered within sub-lineage IVa-1. In contrast, post-vaccine G2 sequences clustered mainly within sub-lineage IVa-3, whereby a radical amino acid (AA) substitution, S15F, was observed between the two sub-lineages. Pre-vaccine P[4] sequences predominantly segregated within sub-lineage IVa while post-vaccine sequences clustered mostly within sub-lineage IVb, with a radical AA substitution R162G. Both S15F and R162G occurred outside recognised antigenic sites. The AA residue at position 15 is found within the signal sequence domain of Viral Protein 7 (VP7) involved in translocation of VP7 into endoplasmic reticulum during infection process. The 162 AA residue lies within the hemagglutination domain of Viral Protein 4 (VP4) engaged in interaction with sialic acid-containing structure during attachment to the target cell. Free energy change analysis on VP7 indicated accumulation of stable point mutations in both antigenic and non-antigenic regions. The segregation of South African G2P[4] strains into pre- and post-vaccination sub-lineages is likely due to erstwhile hypothesized stepwise lineage/sub-lineage evolution of G2P[4] strains rather than RVA vaccine introduction. Our findings reinforce the need for continuous whole-genome RVA surveillance and investigation of contribution of AA substitutions in understanding the dynamic G2P[4] epidemiology.

The Common Occurrence of Theileria ovis in Tibetan Sheep and the First Report of Theileria sinensis in Yaks from Southern Qinghai, China

PURPOSE

The tick-borne protozoa piroplasms, including Theileria and Babesia, are the cause of substantial economic losses to the livestock industry. However, in southern Qinghai province, China, there are limited information on the molecular characteristics of piroplasms. This study therefore aimed at determining the prevalence and genetic diversity of piroplasms.

METHODS

In order to detect piroplasms, we examined 330 yaks and 236 Tibetan sheep blood samples by nested PCR. The differences in piroplasms prevalence in relation to different risk factors was analyzed using SPSS 26. Phylogenetic analysis based on 18S rRNA sequences was inferred using MEGA 7.

RESULTS

During this study, Theileria spp. were detected in 33.6% (111/330) of yaks and 94.1% (222/236) of Tibetan sheep, but no Babesia was identified. Importantly, a comparison study revealed that T. ovis infection was highly prevalent in sheep (94.1%) but infrequent in yaks (6.1%), while T. sinensis was host-specific to yaks with an infection rate of 27.6%. In addition, male animals were more likely to be infected by T. sinensis and female animals were more likely to be infected by T. ovis. And animals below 4000 m areas reported a higher infection rate with T. sinensis (26.1% vs. 2.9%, p < 0.001). Alongside these differences in prevalence, we found a significantly higher T. sinensis infection rate in separated-grazing livestock (22.2% vs. 3.7%, p < 0.001), while mixed-grazing ruminants had a higher T. ovis infection rate (50.0% vs. 39.0%, p = 0.014). Furthermore, sequence analysis revealed that the 18S rRNA sequences obtained in this study shared 86.9-100.0% identities with each other and they were clustered into T. sinensis or T. ovis.

CONCLUSION

To our knowledge, this is the first report of T. sinensis in Qinghai region. In addition, high prevalence of the generally sub-clinical T. ovis in sheep indicates extensive exposure to ticks and transmission of tick-borne pathogens with a significant economic impact. This study provides insights into the distribution and genetic diversity of Theileria in China.

Genetic Characterisation of South African and Mozambican Bovine Rotaviruses Reveals a Typical Bovine-like Artiodactyl Constellation Derived through Multiple Reassortment Events

This study presents whole genomes of seven bovine rotavirus strains from South Africa and Mozambique. Double-stranded RNA, extracted from stool samples without prior adaptation to cell culture, was used to synthesise cDNA using a self-annealing anchor primer ligated to dsRNA and random hexamers. The cDNA was subsequently sequenced using an Illumina MiSeq platform without prior genome amplification. All strains exhibited bovine-like artiodactyl genome constellations (G10/G6-P[11]/P[5]-I2-R2-C2-M2-A3/A11/A13-N2-T6-E2-H3). Phylogenetic analysis revealed relatively homogenous strains, which were mostly related to other South African animal strains or to each other. It appears that these study strains represent a specific bovine rotavirus population endemic to Southern Africa that was derived through multiple reassortment events. While one Mozambican strain, MPT307, was similar to the South African strains, the second strain, MPT93, was divergent from the other study strains, exhibiting evidence of interspecies transmission of the VP1 and NSP2 genes. The data presented in this study not only contribute to the knowledge of circulating African bovine rotavirus strains, but also emphasise the need for expanded surveillance of animal rotaviruses in African countries in order to improve our understanding of rotavirus strain diversity.

Whole Genome Analysis of Human Rotaviruses Reveals Single Gene Reassortant Rotavirus Strains in Zambia

Rotarix® vaccine was implemented nationwide in Zambia in 2013. In this study, four unusual strains collected in the post-vaccine period were subjected to whole genome sequencing and analysis. The four strains possessed atypical genotype constellations, with at least one reassortant genome segment within the constellation. One of the strains (UFS-NGS-MRC-DPRU4749) was genetically and phylogenetically distinct in the VP4 and VP1 gene segments. Pairwise analyses demonstrated several amino acid disparities in the VP4 antigenic sites of this strain compared to that of Rotarix®. Although the impact of these amino acid disparities remains to be determined, this study adds to our understanding of the whole genomes of reassortant strains circulating in Zambia following Rotarix® vaccine introduction.

Nosocomial acute gastroenteritis outbreak caused by an equine-like G3P[8] DS-1-like rotavirus and GII.4 Sydney[P16] norovirus at a pediatric hospital in Rio de Janeiro, Brazil, 2019

Worldwide, rotavirus (RVA) and norovirus are considered major etiological agents of acute gastroenteritis (AGE) in pediatric population admitted to hospitals. This study describes the investigation of nosocomial infections caused by emergent RVA and norovirus strains reported at a pediatric hospital in southern Brazil in May 2019. This outbreak affected 30 people among children and adults. Nine stool samples (eight children and one nurse) were obtained and analyzed by RT-qPCR to detect and quantify RVA and norovirus. Positive samples were genotyped by sequencing and subjected to phylogenetic analysis. We detected RVA in 44.4% (4/9) and norovirus in 55.5% (5/9) at high viral loads, ranging from 3.5 × 10⁷ to 6.1 × 10⁷ and 3.2 × 10² to 3.2 × 10⁹ genome copies/g of stool, respectively. Co-infections were not observed. RVA VP4 and VP7 gene sequencing in combination with polyacrylamide gel electrophoresis identified the circulation of equine-like G3P[8] DS-1-like, and the partial sequencing of the other nine genes revealed that strains possessed I2-R2-C2-M2-A2-N1-T2-E2-H2 genotype background. The emergent recombinant norovirus variant, GII.4 Sydney[P16], was identified by ORF1-2 sequencing. Active surveillance and effective prevention measures should be constantly reinforced to avoid the spread of nosocomial viral infections into hospitals, which could severely affect pediatric patients admitted with underlying health conditions.

Molecular analysis of G3P[6] rotavirus in the Amazon region of Brazil: evidence of reassortment with equine-like strains

Aim: To perform a molecular analysis of rotavirus A (RVA) G3P[6] strains detected in 2012 and 2017 in the Amazon region of Brazil. Materials & methods: Eighteen RVA G3P[6] strains were collected from children aged under 10 years hospitalized with acute gastroenteritis, and partial sequencing of each segment genome was performed using Sanger sequencing. Results: Phylogenetic analysis showed that all G3P[6] strains had a DS-1-like genotype constellation. Two strains had the highest nucleotide identities with equine-like G3P[6]/G3P[8] genotypes. Several amino acid alterations in VP4 and VP7 neutralizing epitopes of equine-like RVA G3P[6] strains were observed in comparison with vaccine strains. Conclusion: These findings suggest that equine-like RVA G3P[6] strains have been circulating in the Amazon region of Brazil as a result of direct importation, and support natural RVA evolutionary mechanisms.

Emergence in 2017-2019 of novel reassortant equine-like G3 rotavirus strains in Palermo, Sicily

Rotavirus A (RVA) is a major etiologic agent of gastroenteritis in children worldwide. Hospital-based surveillance of viral gastroenteritis in paediatric population in Palermo (Italy) from 2017 onwards revealed a sharp increase in G3P[8] RVAs, accounting for 71% of all the RVAs detected in 2019. This pattern had not been observed before in Italy, with G3 RVA usually being detected at rates lower than 3%. In order to investigate this unique epidemiological pattern, the genetic diversity of G3 RVAs identified during a 16-year long surveillance (2004-2019) was explored by systematic sequencing of the VP7 and VP4 genes and by whole genome sequencing of selected G3 strains, representative of the various RVA seasons. Sequence and phylogenetic analyses of the VP7 and VP4 genes revealed the emergence, in 2017 of reassortant equine-like G3P[8], which gradually replaced former G3P[8] strains. The G3P[8] circulating before 2017 showed a Wa-like constellation of genome segments while the G3P[8] that emerged in 2017 had a DS-1-like backbone. On direct inspection of the VP7 and VP4 antigenic epitopes, the equine-like G3P[8] strains possessed several amino acid variations in neutralizing regions compared with vaccine strains. The equine-like G3P[8] RVAs are a further example of the zoonotic impact of animal viruses on human health.

Genotype constellations of the rotavirus A strains circulating in Nizhny Novgorod, Russia, 2017-2018

Currently, the full-genome-based classification is widely used to investigate rotavirus A (RVA) strains found in different countries around the world. However, the information on the full genotypes of rotaviruses circulating in Russia is limited. Using partial sequencing, this study determined the full genotype constellations of 15 RVA strains in total commonly detected in Nizhny Novgorod (European part of Russia) in 2017-2018, three from each of the following genotypes G1P[8], G4P[8], and G9P[8] and six from G2P[4]. There were two intergenogroup mono-reassortants possessing an identical genotype constellation of G4-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1 with the DS-1-like NSP4 gene of probably local origin. A variety of subgenotype lineages and their combinations of Wa-like rotaviruses and genetic heterogeneity among G9P[8] and G1P[8] strains were shown on the basis of phylogenetic analysis of each gene. Moreover, two distinct co-circulating variants that differed in all 11 genome segments were found among DS-1-like rotaviruses.

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.

High prevalence of equine-like G3P[8] rotavirus in children and adults with acute gastroenteritis in Thailand

Group A rotavirus (RVA) is a major cause of acute gastroenteritis in infants and young children worldwide. This study aims to clarify the distribution of G/P types and genetic characteristics of RVAs circulating in Thailand. Between January 2014 and September 2016, 1867 stool specimens were collected from children and adults with acute gastroenteritis in six provinces in Thailand. RVAs were detected in 514/1867 (27.5%) stool specimens. G1P[8] (44.7%) was the most predominant genotype, followed by G3P[8] (33.7%), G2P[4] (11.5%), G8P[8] (7.0%), and G9P[8] (1.3%). Unusual G3P[9] (0.8%), G3P[10] (0.4%), G4P[6] (0.4%), and G10P[14] (0.2%) were also detected at low frequencies. The predominant genotype, G1P[8] (64.4%), in 2014 decreased to 6.1% in 2016. In contrast, the frequency of G3P[8] markedly increased from 5.5% in 2014 to 65.3% in 2015 and 89.8% in 2016. On polyacrylamide gel electrophoresis, most (135/140; 96.4%) of the G3P[8] strains exhibited a short RNA profile. Successful determination of the nucleotide sequences of the VP7 genes of 98 G3P[8] strains with a short RNA profile showed that they are all equine-like G3P[8] strains. On phylogenetic analysis of genome segments of two representative Thai equine-like G3P[8] strains, it was noteworthy that they possessed distinct NSP4 genes, one bovine-like and the other human-like. Thus, we found that characteristic equine-like G3P[8] strains with a short RNA electropherotype are becoming highly prevalent in children and adults in Thailand.