Human rhinovirus spatial-temporal epidemiology in rural coastal Kenya, 2015-2016, observed through outpatient surveillance

Seasonality and Genotype Diversity of Human Rhinoviruses during an Eight-Year Period in Slovenia

Due to the high socioeconomic burden of rhinoviruses, the development of prevention and treatment strategies is of high importance. Understanding the epidemiological and clinical features of rhinoviruses is essential in order to address these issues. Our study aimed to define the seasonality and molecular epidemiology of rhinoviruses in Slovenia. Over a period of eight years, a total of 20,425 patients from sentinel primary healthcare settings and sentinel hospitals were examined for a panel of respiratory viruses in the national programme for the surveillance of influenza-like illnesses and acute respiratory infections. The patients were from all age groups and had respiratory infections of various severity. Infection with a rhinovirus was confirmed using an RT-rPCR in 1834 patients, and 1480 rhinoviruses were genotyped. The molecular analysis was linked to demographical and meteorological data. We confirmed the year-round circulation of rhinoviruses with clear seasonal cycles, resulting in two seasonal waves with peaks in spring and autumn. High levels of genotype variability and co-circulation were confirmed between and within seasons and were analysed in terms of patient age, the patient source reflecting disease severity, and meteorological factors. Our study provides missing scientific information on the genotype diversity of rhinoviruses in Slovenia. As most previous investigations focused on exclusive segments of the population, such as children or hospitalised patients, and for shorter study periods, our study, with its design, size and length, contributes complementary aspects and new evidence-based knowledge to the regional and global understanding of rhinovirus seasonality and molecular epidemiology.

Rhinovirus dynamics across different social structures

Rhinoviruses (RV), common human respiratory viruses, exhibit significant antigenic diversity, yet their dynamics across distinct social structures remain poorly understood. Our study delves into RV dynamics within Kenya by analysing VP4/2 sequences across four different social structures: households, a public primary school, outpatient clinics in the Kilifi Health and Demographics Surveillance System (HDSS), and countrywide hospital admissions and outpatients. The study revealed the greatest diversity of RV infections at the countrywide level (114 types), followed by the Kilifi HDSS (78 types), the school (47 types), and households (40 types), cumulatively representing >90% of all known RV types. Notably, RV diversity correlated directly with the size of the population under observation, and several RV type variants occasionally fuelled RV infection waves. Our findings highlight the critical role of social structures in shaping RV dynamics, information that can be leveraged to enhance public health strategies. Future research should incorporate whole-genome analysis to understand fine-scale evolution across various social structures.

Identifying socio-ecological drivers of common cold in Bhutan: a national surveillance data analysis

The common cold is a leading cause of morbidity and contributes significantly to the health costs in Bhutan. The study utilized multivariate Zero-inflated Poisson regression in a Bayesian framework to identify climatic variability and spatial and temporal patterns of the common cold in Bhutan. There were 2,480,509 notifications of common cold between 2010 and 2018. Children aged < 15 years were twice (95% credible interval [CrI] 2.2, 2.5) as likely to get common cold than adults, and males were 12.4% (95 CrI 5.5%, 18.7%) less likely to get common cold than females. A 10 mm increase in rainfall lagged one month, and each 1 °C increase of maximum temperature was associated with a 5.1% (95% CrI 4.2%, 6.1%) and 2.6% (95% CrI 2.3%, 2.8%) increase in the risk of cold respectively. An increase in elevation of 100 m and 1% increase in relative humidity lagged three months were associated with a decrease in risk of common cold by 0.1% (95% CrI 0.1%, 0.2%) and 0.3% (95% CrI 0.2%, 0.3%) respectively. Seasonality and spatial heterogeneity can partly be explained by the association of common cold to climatic variables. There was statistically significant residual clustering after accounting for covariates. The finding highlights the influence of climatic variables on common cold and suggests that prioritizing control strategies for acute respiratory infection program to subdistricts and times of the year when climatic variables are associated with common cold may be an effective strategy.

Understanding Rhinovirus Circulation and Impact on Illness

Rhinoviruses (RVs) have been reported as one of the main viral causes for severe respiratory illnesses that may require hospitalization, competing with the burden of other respiratory viruses such as influenza and RSV in terms of severity, economic cost, and resource utilization. With three species and 169 subtypes, RV presents the greatest diversity within the Enterovirus genus, and despite the efforts of the research community to identify clinically relevant subtypes to target therapeutic strategies, the role of species and subtype in the clinical outcomes of RV infection remains unclear. This review aims to collect and organize data relevant to RV illness in order to find patterns and links with species and/or subtype, with a specific focus on species and subtype diversity in clinical studies typing of respiratory samples.

Trends and Intensity of Rhinovirus Invasions in Kilifi, Coastal Kenya, Over a 12-Year Period, 2007-2018

Background

Rhinoviruses (RVs) are ubiquitous pathogens and the principal etiological agents of common cold. Despite the high frequency of RV infections, data describing their long-term epidemiological patterns in a defined population remain limited.

Methods

Here, we analyzed 1070 VP4/VP2 genomic region sequences sampled at Kilifi County Hospital on the Kenya coast. The samples were collected between 2007 and 2018 from hospitalized pediatric patients (<60 months of age) with acute respiratory illness.

Results

Of 7231 children enrolled, RV was detected in 1497 (20.7%) and VP4/VP2 sequences were recovered from 1070 samples (71.5%). A total of 144 different RV types were identified (67 Rhinovirus A, 18 Rhinovirus B, and 59 Rhinovirus C) and at any month, several types co-circulated with alternating predominance. Within types, multiple genetically divergent variants were observed. Ongoing RV infections through time appeared to be a combination of (1) persistent types (observed up to 7 consecutive months), (2) reintroduced genetically distinct variants, and (3) new invasions (average of 8 new types annually).

Conclusions

Sustained RV presence in the Kilifi community is mainly due to frequent invasion by new types and variants rather than continuous transmission of locally established types/variants.

Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya, 2016-2018

Background: Virus genome sequencing is increasingly utilized in epidemiological surveillance. Genomic data allows comprehensive evaluation of underlying viral diversity and epidemiology to inform control. For human rhinovirus (HRV), genomic amplification and sequencing is challenging due to numerous types, high genetic diversity and inadequate reference sequences. Methods: We developed a tiled amplicon type-specific protocol for genome amplification and sequencing on the Illumina MiSeq platform of two HRV types, A15 and A101. We then assessed added value in analyzing whole genomes relative to the VP4/2 region only in the investigation of HRV molecular epidemiology within the community in Kilifi, coastal Kenya. Results: We processed 73 nasopharyngeal swabs collected between 2016-2018, and 48 yielded at least 70% HRV genome coverage. These included all A101 samples (n=10) and 38 (60.3%) A15 samples.  Phylogenetic analysis revealed that the Kilifi A101 sequences interspersed with global A101 genomes available in GenBank collected between 1999-2016. On the other hand, our A15 sequences formed a monophyletic group separate from the global genomes collected in 2008 and 2019. An improved phylogenetic resolution was observed with the genome phylogenies compared to the VP4/2 phylogenies. Conclusions: We present a type-specific full genome sequencing approach for obtaining HRV genomic data and characterizing infections.

Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya, 2016-2018

Background: Virus genome sequencing is increasingly utilized in epidemiological surveillance. Genomic data allows comprehensive evaluation of underlying viral diversity and epidemiology to inform control. For human rhinovirus (HRV), genomic amplification and sequencing is challenging due to numerous types, high genetic diversity and inadequate reference sequences. Methods: We developed a tiled amplicon type-specific protocol for genome amplification and sequencing on the Illumina MiSeq platform of two HRV types, A15 and A101. We then assessed added value in analyzing whole genomes relative to the VP4/2 region only in the investigation of HRV molecular epidemiology within the community in Kilifi, coastal Kenya. Results: We processed 73 samples collected between 2016-2018, and 48 yielded at least 70% HRV genome coverage. These included all A101 samples (n=10) and 38 (60.3%) A15 samples.  Phylogenetic analysis revealed that the Kilifi A101 sequences interspersed with global A101 genomes available in GenBank collected between 1999-2016. On the other hand, our A15 sequences formed a monophyletic group separate from the global genomes collected in 2008 and 2019. Improved phylogenetic resolution was observed with the genome phylogenies compared to the VP4/2 phylogenies. Conclusions: We present a type-specific full genome sequencing approach for obtaining HRV genomic data and characterizing infections.

Molecular Epidemiology of Human Rhinovirus From 1-Year Surveillance Within a School Setting in Rural Coastal Kenya

Background

Human rhinovirus (HRV) is the most common cause of the common cold but may also lead to more severe respiratory illness in vulnerable populations. The epidemiology and genetic diversity of HRV within a school setting have not been previously described. The objective of this study was to characterize HRV molecular epidemiology in a primary school in a rural location of Kenya.

Methods

Between May 2017 and April 2018, over 3 school terms, we collected 1859 nasopharyngeal swabs (NPS) from pupils and teachers with symptoms of acute respiratory infection in a public primary school in Kilifi County, coastal Kenya. The samples were tested for HRV using real-time reverse transcription polymerase chain reaction. HRV-positive samples were sequenced in the VP4/VP2 coding region for species and genotype classification.

Results

A total of 307 NPS (16.4%) from 164 individuals were HRV positive, and 253 (82.4%) were successfully sequenced. The proportion of HRV in the lower primary classes was higher (19.8%) than upper primary classes (12.2%; P < .001). HRV-A was the most common species (134/253; 53.0%), followed by HRV-C (73/253; 28.9%) and HRV-B (46/253; 18.2%). Phylogenetic analysis identified 47 HRV genotypes. The most common genotypes were A2 and B70. Numerous (up to 22 in 1 school term) genotypes circulated simultaneously, there was no individual re-infection with the same genotype, and no genotype was detected in all 3 school terms.

Conclusions

HRV was frequently detected among school-going children with mild acute respiratory illness symptoms, particularly in the younger age groups (<5-year-olds). Multiple HRV introductions were observed that were characterized by considerable genotype diversity.

Human rhinovirus spatial-temporal epidemiology in rural coastal Kenya, 2015-2016, observed through outpatient surveillance

Background: Human rhinovirus (HRV) is the predominant cause of upper respiratory tract infections, resulting in a significant public health burden. The virus circulates as many different types (168), each generating strong homologous, but weak heterotypic, immunity. The influence of these features on transmission patterns of HRV in the community is understudied. Methods: Nasopharyngeal swabs were collected from patients with symptoms of acute respiratory infection (ARI) at nine out-patient facilities across a Health and Demographic Surveillance System between December 2015 and November 2016. HRV was diagnosed by real-time RT-PCR, and the VP4/VP2 genomic region of the positive samples sequenced. Phylogenetic analysis was used to determine the HRV types. Classification models and G-test statistic were used to investigate HRV type spatial distribution. Demographic characteristics and clinical features of ARI were also compared. Results: Of 5,744 NPS samples collected, HRV was detected in 1057 (18.4%), of which 817 (77.3%) were successfully sequenced. HRV species A, B and C were identified in 360 (44.1%), 67 (8.2%) and 390 (47.7%) samples, respectively. In total, 87 types were determined: 39, 10 and 38 occurred within species A, B and C, respectively. HRV types presented heterogeneous temporal patterns of persistence. Spatially, identical types occurred over a wide distance at similar times, but there was statistically significant evidence for clustering of types between health facilities in close proximity or linked by major road networks. Conclusion: This study records a high prevalence of HRV in out-patient presentations exhibiting high type diversity. Patterns of occurrence suggest frequent and independent community invasion of different types. Temporal differences of persistence between types may reflect variation in type-specific population immunity. Spatial patterns suggest either rapid spread or multiple invasions of the same type, but evidence of similar types amongst close health facilities, or along road systems, indicate type partitioning structured by local spread.