key: cord-0025014-mpyd6e19 authors: Mwita Morobe, John; Kamau, Everlyn; Murunga, Nickson; Gatua, Winfred; Luka, Martha M; Lewa, Clement; Cheruiyot, Robinson; Mutunga, Martin; Odundo, Calleb; James Nokes, D; Agoti, Charles N title: Trends and Intensity of Rhinovirus Invasions in Kilifi, Coastal Kenya, Over a 12-Year Period, 2007–2018 date: 2021-11-16 journal: Open Forum Infect Dis DOI: 10.1093/ofid/ofab571 sha: babf9df630129dcb0a5b46f770d6eae95092a45e doc_id: 25014 cord_uid: mpyd6e19 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. Rhinoviruses (RVs) are a highly prevalent group of viruses and are the principal cause of common cold syndrome in humans globally [1, 2] . RV infections result in a wide range of clinical outcomes spanning from asymptomatic and mild illness in the upper airways to severe illness in the lower airways [3, 4] . The infections occur in all ages, with severe presentation more likely in children under the age of 5 years [5, 6] , the elderly [7] , and immunocompromised persons [8] . Despite the clinical significance of RV infections, there is little information on the long-term trends and diversity of circulating RV types. RV belongs to the genus Enterovirus of the family Picornaviridae. The viral single-stranded positive sense RNA genome consists of approximately 7200 nucleotides and encodes 4 structural proteins (VP4, VP2, VP3, and VP1) and 7 nonstructural proteins (2A, 2B, 2C, 3A, 3B, 3C, and 3D) [2] . The 3 surface-exposed capsid proteins (VP1, VP2, and VP3) carry the antigenically critical sites [9] [10] [11] . The high genetic variability in the VP4/VP2 and VP1 genomic regions of RVs have been instrumental in molecular typing [12, 13] and molecular epidemiological investigations of RV infections [14] [15] [16] . Currently, a total of 169 RV types have been described and classified into 3 species: Rhinovirus A, Rhinovirus B, and Rhinovirus C (https://www.picornaviridae.com/sg3_ensavirinae/enterovirus/enterovirus.htm). RV infections occur year-round in most geographical locations, although peaking in the early autumn and late spring in many temperate countries, and in the rainy season in tropical countries [2, 17] . Unclear seasonality and year-round transmission of RVs have been attributed to lack of intertype crossprotective immunity [18, 19] , coupled with the high genetic diversity within the 3 species, each with the ability to spread independently in a population [14, 16, 20] . A recent study in Kilifi County, located in coastal Kenya, that spanned over a 12-month period [14] (December 2015-November 2016) found that multiple RV types co-circulate over varied time periods ranging from 1 to 9 months and, in most cases, each displaying a typical epidemic curve at the local population level; transmission is presumably constrained by the decline in susceptibles to that type within the locality. Typespecific (homologous) immunity has been reported to wane approximately after a 1-year period [21] , and individuals who were previously immune to a particular type gradually become susceptible to the type again [21, 22] . Previous studies found that introduction of new RV types or sequential invasion by different genetic variants could be due to declining levels of population immunity as well as viral evolution [23, 24] . These assertions of perpetually changing RV types during year-round RV transmission have not been fully investigated in a longitudinal manner [16] . In this study, we analyzed VP4/VP2 sequences of samples collected from hospitalized children with acute respiratory illness between 2007 and 2018 on the Kenyan coast to evaluate the long-term incidence of the different RV types, their temporal patterns, and intensity of new invasions in a local population. The study was conducted at the Kilifi County Hospital (KCH) as part of long-term surveillance initially aimed at understanding the epidemiology and disease burden of respiratory syncytial virus-associated pneumonia cases [25] and expanded to a range of respiratory viruses from 2007 onward [6, [26] [27] [28] [29] [30] . KCH, located on the coast of Kenya, is a referral hospital serving the wider Kilifi County, which has a population of 1 453 787 and covers an area of approximately 12 254 km 2 . Details of study design, participant recruitment, and sampling procedures have been described elsewhere [25, 29] . In brief, upon presentation to the pediatric ward, a detailed medical review was undertaken by the clinician and the decision to admit was made. For this study, children (<60 months of age) admitted to the pediatric ward between January 2007 and December 2018 were eligible if they presented with symptoms of syndromic severe or very severe pneumonia. Clinical definitions include a history of cough or difficulty in breathing for <30 days, which if accompanied by lower chest wall indrawing was defined as severe pneumonia; or if accompanied by any 1 of prostration, coma, or hypoxemia was defined as very severe pneumonia (prostration included the inability to feed or drink, and hypoxemia defined by oxygen saturation [pO 2 ] <90%) [25] . Following a written informed consent from the parent or guardian, a nasopharyngeal flocked swab, nasal wash, or combination of nasopharyngeal swab and oropharyngeal swab was collected from each child and transferred into viral transport medium for laboratory screening. Ethical approval for the study protocol was obtained from the Scientific and Ethics Review Unit (SERU number 3443) ethics committee, Kenya Medical Research Institute, Nairobi, Kenya. Viral RNA was extracted from each sample using QIAamp Viral RNA kit (Qiagen, Valencia, California) and screened for respiratory viruses using a multiplex real-time reverse-transcription PCR (rRT-PCR) (Applied Biosystems, United Kingdom) as described elsewhere [31, 32] . A sample was considered RV positive if the rRT-PCR cycle threshold was <35.0 [30] . A section of VP4/VP2 viral genomic region (~420 nucleotides long) of positive samples was amplified and sequenced as previously described [14] . Consensus sequences were assembled using Sequencher software version 5.4.6 (Gene Codes Corporation, Ann Arbor, Michigan). We used the term "type" to refer to RV sequences classified by either cross-neutralization or genetic comparisons as distinct as described previously [13] . Based on this approach, sequences were assigned into the same RV type based on >90% nucleotide similarity to RV prototype sequences (also referred to as reference sequences, http://www.picornaviridae.com/ sequences/sequences.htm) and phylogenetic clustering with bootstrap support value >70% [13] . Distributions of pairwise genetic distances were assessed for evaluation of intertype and intratype divergence [13] . Intratype "variant" was defined on the basis of a divergence threshold value determined as the least frequent value between the first and second modes in a pairwise nucleotide difference distribution plot. Here we are implicitly assuming that sequences with pairwise nucleotide difference falling into the distribution with the low (first) mode are members of the same phylogenetic clade, whereas those with pairwise nucleotide difference within the second distribution with higher mode are members of different phylogenetic clades. A group of viruses within the first, lower distribution were classified as belonging to the same RV type variant. The definitions used to describe the temporal occurrence of RV types are summarized as follows: 1. Persistent: Continued detection, in consecutive or nonconsecutive years, of a group of viruses belonging to the same variant of a RV type. 2. Recurrent: Detection of a virus or group of viruses not observed in the preceding years (>1 year) that belong to a different variant of a previously observed RV type. 3. Invasion: Detection of a new RV type not previously locally documented. Multiple sequence alignments were generated using MAFFT version 7.220 [33] and maximum likelihood phylogenetic trees estimated using IQ-TREE version 1.6.12 [34] . Branch support was assessed by 1000 bootstrap iterations. Temporal signal in the data was examined using TempEst version 1.5.3 [35] . To infer time-scaled phylogenies, Bayesian phylogenetic analyses were undertaken in BEAST version 1.10.4 assuming an uncorrelated log-normal relaxed molecular model [36] . The Markov chain Monte Carlo convergence was assessed in Tracer version Between January 2007 and December 2018, a total of 7231 nasopharyngeal swab samples were collected from children (<60 months of age) admitted with severe or very severe pneumonia in KCH (Supplementary Table 1 ). RV was detected in 20.7% (1497/7231), with the proportion positive across the years ranging from 15.6% to 38.3% (Supplementary Table 1 ). The monthly frequency of detection of RV in the study population is shown in Figure 1 . RV infections were observed to occur year-round, frequently peaking between the months of May and September each year ( Figure 1 ). A total of 1070 (71.5%) VP4/VP2 sequences (~420 nucleotides, some previously reported [30] ) were available for this analysis. Of these, 520 ( (Table 1) . Twenty-four sequences could not be assigned to known RV types based on the criterion proposed by McIntyre et al [13] due to these sequences having p-distance of >10.5% with respect to their closest reference sequences (Supplementary Table 2 ). Other enteroviruses were also detected on sequencing the rRT-PCR RV detections: enterovirus D68 (EV-D68) (n = 5), coxsackievirus B3 (CVB3) (n = 1), coxsackievirus B2 (CVB2) (n = 1), and echovirus 19 (E19) (n = 1). We detected on average, 39 RV types annually (range, 35-47), a mean of 8 (range, 1-29) of which were new RV types identified for the first time in the population each year from 2008 as other previously detected types disappeared ( Figure 2B ). The cumulative number of new RV types detected annually increased rapidly since the beginning of the surveillance period and then saturated after approximately 9 years ( Figure 2B ). RV types commonly co-circulated and with varying frequency in the 12-year period ( Figure 2C, Supplementary File 1) . Several Abbreviation: RV, rhinovirus. types were present at high prevalence whereas others occurred once or sporadically. Some types circulated consecutively for months; for example, RV-A56 was detected in 7 consecutive months (May-November 2007); RV-C11 was present for 6 consecutive months (February-July 2016); and RV-C38, A40, and C2 types circulated consecutively for 5 months (November 2009 to Additionally, several types recurred after considerable periods of absence. For example, RV-A12, first seen in February 2007, was not detected again until February 2009, 23 months later, whereas C38 viruses were detected 4 years apart between 2012 and 2016 ( Figure 3, Supplementary File 1) . Temporally, several RV types exhibited synchronized co-circulation and recurrence, for example: (1) RV-C1, C11, C2, C38, C22, and C21; (2) RV-A75, A89, A12, A28, A96, A106, A80, and A10; (3) RV-A90, A55, A61, A45, A54, and A60; (4) RV-C14, C41, C45, C10, C16, C25, C32, and C47. The nucleotide sequence identity among Rhinovirus A, Rhinovirus B, and Rhinovirus C viruses was determined as 57.3%-100%, 66.0%-100%, and 45.1%-100%, respectively, and 59.8%-100%, 79.9%-100% and 53.3%-100% at the amino acid level, respectively. Intratype nucleotide variation was observed in the VP4/VP2 region of viruses sampled over the 12-year study period ( Figure 4A, Supplementary Figure 1) . Nonetheless, the substitutions were mostly synonymous, that is, not amino acid changing. The distribution of pairwise nucleotide distances showed multimodal peaks suggesting circulation of distinct variants within individual RV types ( Figure 4B, Supplementary Figure 2) . These observations were congruent with multiple within-type phylogenetic clusters. Several RV types were characterized by genetically distinct temporal clusters, for example, RV-A49, C38, and A101 ( Figure 5 ). RV-A49 was detected as 11 distinct variants circulating at different periods, 3 of which occurred as singletons (single sequences), suggesting undersampled genetic diversity ( Figure 5 , Table 2 ). Multiple genetic variants of type RV-C6 Table 2 ), which likely indicates separate virus introductions into the Kilifi population. Several RV types had variants that contained sequences from multiple years indicating variant persistence over an extended period or repeated reintroductions, for example, RV-A101 variant 5 comprised of viruses observed from 2010 to 2013 ( Figure 5 , Table 2 ). We describe the long-term pattern of co-circulation, persistence, and invasion of RV types in hospitalized children (<60 months old) with pneumonia in Kilifi, coastal Kenya, over a 12-year period (2007-2018) . Consistent with other studies, RV was ubiquitous and multiple types co-circulated even within a single month [14, 16, 20] . Among the RV cases detected, Rhinovirus B was least frequently detected. It is not clear why Rhinovirus B is less diverse and each type within it was on average less frequent. The observed annual proportions of RV species in Kilifi are consistent with recent similar epidemiological studies in Brazil, Nigeria, and Cameroon [37] [38] [39] . Although children <5 years of age are not a comprehensive representative of the community, this demographic gives insight into the pattern of RV transmission since RV burden is highest in children <5 years of age [14] . RV detection rates decrease with increasing age as adults have had multiple and widespread exposures to RV types [14] . Other social groups are vital in RV transmission [6, 40, 41] and it would be useful to evaluate RV transmission patterns and prevalence within these groups. The majority (99%) of our sequences were within the proposed divergence thresholds for RV typing and classification using the VP4/VP2 region (10.5% for Rhinovirus A, 9.5% for Rhinovirus B, and 10.5% for Rhinovirus C) [13] . This exemplifies significant sequence conservation in the VP4/VP2 region within a type allowing robust genotypic assignment. However, 24 sequences did not fit the classification system for VP4/VP2 region and require whole genome sequencing to check for variation in the VP1 region and determine if they correspond to new types [13] . Detection of other enteroviruses reflects PCR cross-reactivity due to nucleotide conservation at the 5ʹ-untranslated target region [42, 43] . EV-D68, CVB3, CVB2, and E19 have been associated with respiratory disease or detected in respiratory samples [44, 45] . The frequent invasions of new types could be explained by lack of preexisting immune memory or weak heterotypic immunity [46] . The number of new types decreased over time, levelling off in 2016, perhaps indicating the period a population takes to experience the maximum number of RV types. Recurrence of RV types could be promoted by antigenic variation on the other surface proteins (VP2, VP1, and VP3) allowing infection where prior exposure confers incomplete or short-lived immunity to future genetic variants. Recurrence, particularly where the recurring strains were genetically identical to older strains, may also be observed in a population not previously exposed to a RV type. Some RV types occurred sporadically and could be associated with mild disease or asymptomatic infections or have reduced transmission rates probably suppressed by preexisting immunity [47] . For some RV types, the sequenced VP4/VP2 region remained conserved after periods of quiescence, which probably ensures strain survival by maintaining low-level genetic variation. In a linear strain space, strains interact via cross-immunity to nearby strains with shared epitopes, and this interaction tails off with genetic distance [48] . Yet, the VP4/VP2 region might not be primarily antigenic [49] , and genetic changes could have occurred at immunogenic sites located in other capsid proteins (VP1 or VP3). Genome-wide sequence data would therefore be useful to confirm strain conservation and maintenance. The evident intratype genetic diversity with differential temporal distribution could suggest sequential virus introductions or diversification of locally circulating variants [50] . This study had 2 limitations. First, in some years (2010-2013 and 2015) we only sequenced a proportion of the positive cases (Supplementary Table 1 ), which might underestimate the circulating RV diversity. Samples selected were prioritized based on viral load and monthly distribution. Second, we only sequenced the VP4/VP2 coding region, but more reliable phylogenetic relationships would be defined from full-length genome analysis [51] . In conclusion, this study describes the nature of RV infections in hospitalized children <60 months old and enhances our understanding on RV transmission dynamics in a community. RV dynamics in Kilifi during 2007-2018 were characterized by repeated invasions by heterogeneous types rather than longterm continuity of the same RV types and continuous diversification of circulating variants. Improved understanding on the RV types circulating in a community may support better guidance of future therapeutic interventions in clinical practice. The high diversity and rates of invasion of RV as observed in this study, even within a short duration (week or month), underpins the application of molecular typing for surveillance and understanding virus epidemiological dynamics. study samples. We also thank all members of the Virus Epidemiology and Control (VEC) Research Group in Kilifi who were involved at various stages of this study. We also thank the parents and guardians of the children for agreeing to participate in this study. Data availability. All data generated and analysis script for this manuscript are available from the Virus Epidemiology and Control, Kenya Medical Research Institute (KEMRI)-Wellcome Trust Research Programme, data server: https://doi.org/10.7910/DVN/OL399P. Disclaimer. The funder had no role in other aspects of the study including its design, data collection, data analysis, data interpretation, or writing of this manuscript. This was submitted for publication with permission from Director of KEMRI. Financial support. This study was supported by the Wellcome Trust, United Kingdom (grant number 102975). Potential conflicts of interest. All authors: No reported conflicts of interest. Viral etiology of respiratory tract infections in children at the pediatric hospital in Ouagadougou Human rhinoviruses Human rhinovirus infections in rural Thailand: epidemiological evidence for rhinovirus as both pathogen and bystander Predominance of rhinovirus in the nose of symptomatic and asymptomatic infants New Vaccine Surveillance Network. Rhinovirus-associated hospitalizations in young children Molecular epidemiology of human rhinovirus from 1-year surveillance within a school setting in rural coastal Kenya Unexpectedly higher morbidity and mortality of hospitalized elderly patients associated with rhinovirus compared with influenza virus respiratory tract infection Persistent rhinovirus infection in pediatric hematopoietic stem cell transplant recipients with impaired cellular immunity Evidence for at least two dominant neutralization antigens on human rhinovirus 14 Use of monoclonal antibodies to identify four neutralization immunogens on a common cold picornavirus, human rhinovirus 14 Overview of human rhinovirus immunogenic epitopes for rational vaccine design Proposals for the classification of human rhinovirus species C into genotypically assigned types Proposals for the classification of human rhinovirus species A, B and C into genotypically assigned types Human rhinovirus spatial-temporal epidemiology in rural coastal Kenya, 2015-2016, observed through outpatient surveillance Molecular epidemiological study of human rhinovirus species A, B and C from patients with acute respiratory illnesses in Japan A molecular epidemiological perspective of rhinovirus types circulating in Amsterdam from Human rhinoviruses and enteroviruses in influenza-like illness in Latin America Challenges in developing a cross-serotype rhinovirus vaccine Cross-serotype immunity induced by immunization with a conserved rhinovirus capsid protein Rhinovirus infections in western Sweden: a four-year molecular epidemiology study comparing local and globally appearing types The time course of the humoral immune response to rhinovirus infection Evaluation of an enzyme-linked immunosorbent assay that measures rhinovirus-specific antibodies in human sera and nasal secretions Serial viral infections in infants with recurrent respiratory illnesses Highly frequent infections with human rhinovirus in healthy young children: a longitudinal cohort study Incidence and severity of respiratory syncytial virus pneumonia in rural Kenyan children identified through hospital surveillance Viral etiology of severe pneumonia among Kenyan infants and children A preliminary study of pneumonia etiology among hospitalized children in Kenya Human metapneumovirus prevalence and patterns of subgroup persistence identified through surveillance of pediatric pneumonia hospital admissions in coastal Kenya Surveillance of endemic human coronaviruses (HCoV-NL63, OC43 and 229E) associated with childhood pneumonia in Kilifi Molecular epidemiology of human rhinovirus infections in Kilifi, coastal Kenya Real-time RT-PCR detection of 12 respiratory viral infections in four triplex reactions Added value of an oropharyngeal swab in detection of viruses in children hospitalized with lower respiratory tract infection MAFFT multiple sequence alignment software version 7: improvements in performance and usability IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen) Model averaging and Bayes factor calculation of relaxed molecular clocks in Bayesian phylogenetics Genotypes of rhinovirus detected among children in two communities of south-west Nigeria Molecular epidemiology of enteroviruses and rhinoviruses in patients with acute respiratory infections in Yaounde Molecular and clinical characteristics related to rhinovirus infection in Brasília, Brazil Rhinovirus transmission within families with children: incidence of symptomatic and asymptomatic infections An intensive, active surveillance reveals continuous invasion and high diversity of rhinovirus in households New molecular detection tools adapted to emerging rhinoviruses and enteroviruses Improved molecular identification of enteroviruses by RT-PCR and amplicon sequencing Rhinoviruses and respiratory enteroviruses: not as simple as ABC Enterovirus D68 and human respiratory infections Full-genome sequence and analysis of a novel human rhinovirus strain within a divergent HRV-A clade Infectious disease dynamics: what characterizes a successful invader Dynamics and selection of many-strain pathogens VP1 sequencing of all human rhinovirus serotypes: insights into genus phylogeny and susceptibility to antiviral capsid-binding compounds Successive respiratory syncytial virus epidemics in local populations arise from multiple variant introductions, providing insights into virus persistence Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya Acknowledgments. We thank the study participants for providing the All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed. Supplementary materials are available at Open Forum Infectious Diseases online. Consisting of data provided by the authors to benefit the reader, the posted materials are not copyedited and are the sole responsibility of the authors, so questions or comments should be addressed to the corresponding author.