Next Article in Journal
Isolation, Identification, and Pathogenicity of a Goose Astrovirus Genotype 1 Strain in Goslings in China
Previous Article in Journal
PfAgo-Based Zika Virus Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Epidemiology and Clinical Characteristics of an Outbreak on Respiratory Virus Coinfection in Gansu, China

1
Jingyuan County Center for Disease Control and Prevention, Jingyuan 730699, China
2
Gansu Provincial Center for Disease Control and Prevention, Key Laboratory of Pathogenic Microbiology and Immunology (Gansu Province), Lanzhou 730050, China
3
Department of Health Inspection, Vaccine Research and Evaluation Center, School of Public Health, Peking University, Beijing 101100, China
4
Scientific Research Training Department, Baiyin City Center for Disease Control and Prevention, Baiyin 730900, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2024, 16(4), 540; https://doi.org/10.3390/v16040540
Submission received: 27 February 2024 / Revised: 26 March 2024 / Accepted: 26 March 2024 / Published: 30 March 2024
(This article belongs to the Section Human Virology and Viral Diseases)

Abstract

:
This study aims to analyze the epidemiological and pathogenic characteristics of an outbreak primarily caused by respiratory syncytial virus (RSV), human rhinovirus (HRV), and human metapneumovirus (HMPV) in a kindergarten and primary school. The outbreak was investigated by field epidemiological investigation, and the common respiratory pathogens were screened by RT-PCR detection technology. The attack rate of this outbreak was 63.95% (110/172). Main symptoms included cough (85.45%), sore throat (60.91%), and sneezing (60.00%). Multifactorial logistic regression analysis revealed that continuous handwashing and mouth and nose covering when sneezing were protective factors. All 15 collected throat swab specimens tested positive for viruses, with HMPV as the predominant pathogen (80.00%), followed by HRV (53.33%), and two cases of positive respiratory syncytial virus (13.33%). Among them, six samples showed coinfections of HMPV and HRV, and one had coinfections of HMPV and RSV, resulting in a coinfection rate of 46.67%. Genetic sequencing indicated that the HMPV genotype in this outbreak was A2c, and the HRV genotype was type A, resulting in a coinfection outbreak of HMPV, HRV, and RSV in schools and kindergartens, suggesting that multi-pathogen surveillance of respiratory tract infections should be strengthened.

1. Introduction

Respiratory tract infections are prevalent infectious diseases, with a significant impact on infants, the elderly, individuals with chronic illnesses, hematopoietic stem cell transplant recipients, and immunocompromised individuals, as evidenced by their ranking among the top five causes of global mortality in 2015 [1,2]. Viruses are the predominant pathogens responsible for acute respiratory infections in children, primarily transmitted through air, droplets, and close contact [3]. Given that schools serve as significant gathering places for children, they are susceptible to becoming hotspots for the transmission of respiratory infections.
Severe outcomes such as hospitalization and mortality in cases of infection are predominantly associated with viral respiratory infections. Notably, respiratory syncytial virus (RSV), human rhinovirus (HRV), and human metapneumovirus (HMPV) are significant pathogens responsible for acute respiratory infections across all age groups globally [4,5]. RSV, first isolated from children in 1957, is considered the most critical virus causing acute respiratory infections in children [6]. HRV, discovered in the 1950s, is a causative agent for over half of common cold-like illnesses [7,8]. HRV is commonly associated with the common cold and mild upper respiratory tract infections, but it can also cause severe respiratory infections in the immunocompromised or those with chronic lung diseases [9,10]. These three viruses share high similarities in transmission characteristics and clinical symptoms, leading to a spectrum of diseases ranging from mild upper respiratory infections to severe lower respiratory infections in children and adults alike [4]. The majority of infants and children have been infected with RSV, HMPV, and HRV, leading to prevalent diseases such as bronchiolitis and pneumonia. This trend persists in adults, where almost all individuals show serological evidence of prior infections with RSV, HMPV, and HRV [11,12,13]. While healthy adults typically experience self-limiting symptoms following infection, individuals with compromised immune systems are at an elevated risk of developing severe illnesses after being infected with RSV, HMPV, and HRV [9,14,15]. Due to insufficient attention given to this population, they contribute significantly to the global burden of hidden diseases each year.
Gansu Province in the northwest of China displays a wide range of climatic conditions, encompassing high-altitude mountainous areas, desert grasslands, and temperate semi-humid regions. Its geographical location, bordering multiple provinces, may impact population movements from other areas and consequently influence the speed and extent of epidemic spread. In contrast to the eastern coastal provinces of China, Gansu Province exhibits a relatively lower population density. Lower population density implies reduced inter-community contact, but it may also lead to the dispersed distribution of medical resources, presenting challenges for epidemic control. The respiratory infection epidemic simultaneously occurred in adjacent schools and kindergartens. Long-term monitoring data from the Gansu Provincial Center for Disease Control and Prevention indicate a different incidence pattern compared to previous reports, prompting special investigations. This research aims to comprehensively analyze the epidemiological features and pathogen characteristics of the epidemic.
Through an in-depth examination of this epidemic, we aim to better understand the transmission mechanisms of respiratory infections in closed environments, providing scientific evidence for developing more effective prevention and control strategies in the future. The outcomes of this study will contribute to a deeper understanding of the transmission mechanisms of respiratory infections, offer robust support for responding to similar events in the future, and serve as a reference for the prevention and control of respiratory outbreaks in collective settings.

2. Materials and Methods

2.1. Overview of the Epidemic

This epidemic occurred in a local kindergarten and primary school, consisting of 9 grades and nine classes. The educational institutions had a combined population of 140 students and children, and 32 staff members. Among these were 84 primary school students, 23 primary school staff, 56 kindergarten children, and 9 kindergarten staff. Both schools operated as day schools, occupying spacious buildings with classrooms. During school hours, doors and windows could be opened for ventilation, and no healthcare physicians were on-site. The schools did not provide meals or accommodation, and they were approximately 400 m from residential buildings. Students commuted daily between school and home on foot, escorted by parents. On 25 May 2023, the kindergarten reported that more than 20 students were absent due to illness, and the primary school reported that more than 10 students in a class were absent due to illness, with the main symptoms being cough, fever, sore throat, etc. The county CDC immediately entered the scene to investigate, including the basic information about the patients, the incidence, epidemiological characteristics, respiratory symptoms, signs, and so on. In running a case search on suspicious patients between 1 May and 15 June, 101 cases meeting the case definition were identified in the two institutions, affecting children, students, and staff, resulting in an overall attack rate of 72.14% (101/140). The highest attack rate was observed in the kindergarten’s small class, reaching 92.86% (13/14), while the staff’s attack rate was 28.13% (9/32).

2.2. Study Subjects

All kindergarten children, primary school students, and staff with or without respiratory symptoms between 1 May and 15 June 2023, in kindergartens and primary schools within a township of Gansu Province were included as subjects. The case definitions were as follows: (a) Suspected cases: Individuals who exhibited coughing, phlegm production, sore throat, runny nose, or other respiratory clinical symptoms lasting more than one day. They may also experience fatigue, fever (axillary temperature ≥ 37.3 °C), headache, limb pain, or other systemic symptoms lasting more than one day. (b) Confirmed cases: Suspected cases whose throat swabs tested positive for respiratory virus nucleic acid through Reverse Transcription Polymerase Chain Reaction (RT-PCR). (c) Control group: Kindergarten children, primary school students, and staff who had not shown any signs of respiratory tract infection in the past 45 days within kindergartens and primary schools in a township of Gansu Province were selected as the healthy control group through investigation.

2.3. Outbreak Investigation

The investigation involved searching for cases through local medical records, interviews with kindergarten caregivers and teachers, and on-site and telephone follow-ups with parents. A standardized questionnaire covering basic information, illness details, clinical data, influenza vaccination status, contact history, and risk factors was designed and administered to each class. Both the kindergartens and primary schools implemented a daily epidemic reporting system.
A comprehensive survey was conducted on the overall layout of kindergartens and primary schools, children’s home living conditions, environmental hygiene, hygiene facilities in kindergarten bathrooms, handwashing facilities, and disinfection status of two locations. Parents’ and children’s travel history and participation in extracurricular activities in the last seven days were also investigated.
On-site questionnaires were conducted for parents of affected kindergarten children, caregivers, and class teachers. Controls were randomly selected within the class, meeting the criteria of no corresponding respiratory clinical symptoms from 1 May to 15 June 2023. The investigation of case and control groups in kindergartens and primary schools was conducted based on the questionnaire content.
Local Centers for Disease Control (CDC) epidemiological investigators underwent unified training, adhering to standardized questionnaire understanding. The questionnaires were cross-checked by designated personnel, and any missing or incomplete items were supplemented through telephone follow-ups. Double-entry of questionnaire data was performed, and consistency checks were conducted on the entered results.

2.4. Sample Collection and Pathogen Detection

Throat swab samples of suspected cases were collected by qualified medical personnel, and case data were collected and stored in a non-inactivated sampling tube containing 3–4 mL sterile sampling solution (Yocon Biotechnology (Beijing) Co., Ltd., Beijing, China, Batch number: 01230319). Nasal swab samples were collected, and the novel coronavirus (2019-nCoV) antigen detection was performed using the colloidal gold method on site (Shanghai Biogerm Medical Technology Co., Ltd., Shanghai, China, batch number 20221109A). The sampling and detection procedures were conducted per the kit instructions. Throat swab samples were collected and stored at 4 °C within 24 h and transported to the Virus Laboratory of Gansu Provincial Center for Disease Control and Prevention. Nucleic acid was extracted from pharyngeal swab samples of 15 cases by the magnetic bead method (Guangzhou Daan Gene Co., Ltd., Guangzhou, China, article No: DA1003), and the extraction process was carried out strictly following the kit instructions. RT-PCR method was used to detect six respiratory viruses (Shanghai Berger Medical Technology Co., Ltd., Shanghai, China, article No.: SJ-HX-807-2): respiratory syncytial virus (RSV), parainfluenza virus I (PIV-I), parainfluenza virus III (PIV-III), human rhinovirus (HRV), adenovirus (ADV), and human metapneumovirus (HMPV). The experimental methods and results were interpreted per the kit instructions.

2.5. Amplification and Genotype Identification of Target Genes

The pathogen detection results showed mainly HMPV and HRV infection, so the HMPV- and HRV-positive samples were further identified. RT-PCR was used to amplify the G gene in HMPV-positive samples, and the VP1 region was identified in HRV-positive samples. The Measles Room provided the identification primer, the National Institute of Viral Diseases, Chinese Center for Disease Control and Prevention, and the amplification reagent was used with One Step RT-PCR kit (TAKARA, Tokyo, Japan, article number: DRR057A). Sanger sequencing was used on PCR amplification products. BLAST online through https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 25 May 2023). was used to obtain the nucleotide sequence alignment sequence and determine the most similar sequences. Regarding the literature from Genbank (https://www.ncbi.nlm.nih.gov/) (accessed on 25 May 2023), we downloaded RSV for each genotype reference sequence through MEGA7.0 software, together with the above sequence-building related evolutionary tree.

2.6. Statistical Analysis

EpiDate3.1 established the database, and SPSS 19.0 was used for statistical analysis. Descriptive epidemiological analysis methods were employed to illustrate the temporal, spatial and population distribution of the outbreak, and logistic regression models were applied to analyze factors influencing the outbreak for the case–control studies (comparison between 110 cases and 64 controls). A significance level of α = 0.05 (two-tailed) was used, with p < 0.05 indicating statistical significance.

3. Results

The affected kindergarten and primary school relocated to the current residential building in 2018. The kindergarten is adjacent to the primary school, with classrooms covering an area of 90 m2. The average space per person is 2 m2, including coatrooms and bathrooms. The bathroom facilities are complete, and no windows are directly connected to the outdoors, ensuring a well-maintained indoor environment. The primary school classrooms have an area of 72 m2, offering good ventilation. Each floor is equipped with two separate bathrooms for males and females. Although the bathrooms have windows, they are kept closed. The handwashing facilities are complete, but hand soap is not provided. Most faucets are not functioning correctly. The primary school offers teachers bottled, purified water while students bring their own water. The kindergarten, on the other hand, centrally supplies purified water for drinking.
In the outbreak, 110 cases of acute upper respiratory infections were reported, resulting in a morbidity rate of 63.95% (110/172). All cases were mild, with no reports of severe or fatal cases. The main symptoms included coughing (85.45%), sore throat (60.91%), sneezing (60.00%), nasal congestion (54.55%), runny nose (54.55%), fever (53.64%), and fatigue (53.64%), among others, as detailed in Table 1.

3.1. Epidemiological Characteristics

3.1.1. Temporal Distribution

The first case in the kindergarten occurred on 10 May, affecting a preschooler in the senior class, and the last case appeared on 14 June, spanning 36 days. The epidemic exhibited a bimodal distribution, with peaks on 17 May and 8 June, consistent with the characteristics of intergenerational transmission in respiratory infectious diseases. The value-added pattern revealed five generations of cases (Figure 1A).
The first case emerged on 1 May; subsequent generations occurred every 7 to 8 days. A similar pattern was observed in the primary school, with the first case on 1 May and a 46-day duration. The epidemic showed a sawtooth distribution, and five generations appeared, following a similar pattern of intervals (Figure 1B).

3.1.2. Spatial Distribution

Cases were distributed across various classes in both the primary school and kindergarten. The primary school had 57 cases with a morbidity rate of 53.27% (57/107), while the kindergarten had 53 cases with a morbidity rate of 81.54% (53/65). The difference in morbidity rates between the kindergarten and primary school was statistically significant (χ2 = 14.016, p < 0.05). Morbidity rates varied significantly among different grades in the primary school, ranging from 18.75% (3/16) in the fifth grade to 86.67% (13/15) in the first grade (χ2 = 19.095, p < 0.05). The kindergarten showed no significant difference in morbidity rates between classes (χ2 = 1.131, p > 0.05). Table 2 provides further details. The analysis indicated clustering in the distribution of cases in both kindergarten and primary school classes (χ2 = 13.83, p < 0.01; χ2 = 268.38, p < 0.01). Intra-family relationships, especially among siblings, contributed to a significant contact history between different classes within the same school. Out of 44 cases from 20 families, family clustering constituted 40% (44/110) of the total cases.

3.1.3. Population Distribution

Out of the 110 cases, students accounted for 101, resulting in a morbidity rate of 72.14% (101/140). Faculty members experienced nine cases, with a morbidity rate of 28.13% (9/32). The difference in morbidity rates between students and faculty was statistically significant (χ2 = 57.607, p = 0.000). Kindergarten children had a morbidity rate of 87.50% (49/56), while faculty had a rate of 44.44% (4/9), with a significant difference (adjusted χ2 = 6.903, p = 0.009). In the primary school, student morbidity was 61.90% (52/84), and faculty morbidity was 21.74% (5/23), also showing a significant difference (χ2 = 11.702, p = 0.001). Gender-wise analysis revealed no significant difference in morbidity rates between males and females in both the kindergarten and primary school (Table 2).

3.2. Case–Control Study

A total of 174 questionnaires were collected, with a response rate of 100.00% (110 cases, 64 controls) after excluding two duplicate questionnaires. Variables such as infection status, population category, history of COVID-19 infection, contact with flu-like cases in the last seven days, location of contact, occupation, age, handwashing duration, covering mouth and nose when sneezing, immediate handwashing after sneezing, immediate handwashing after class, touching lips/fingers, eye-rubbing habits, cleanliness of hands and nails, use of transportation to school, time spent in public places (excluding school), family members covering mouth and nose when sneezing, home ventilation status, preventive medication, sleep time, frequency of colds in the past six months, and participation in extracurricular training were considered as independent variables. After adjusting for confounding factors such as gender, multiple logistic regression analysis showed that handwashing for 10–30 s or >30 s and family members covering mouth and nose when sneezing were protective factors (OR = 0.06, 95.0% CI = 0.01–0.51, p = 0.01; OR = 0.06, 95.0% CI = 0.01–0.34, p = 0.00; OR = 0.23, 95.0% CI = 0.06–0.95, p = 0.04). Factors such as contact with flu-like cases in the last seven days, walking to and from school, window ventilation (<2 h), frequency of colds (≥3 times or = 2 times), and participation in extracurricular training were identified as risk factors for the outbreak, as shown in Table 3 (OR = 1.18, 95% CI = 1.05–1.60, p = 0.01; OR = 13.47, 95.0% CI = 1.92–94.75, p = 0.01; OR = 27.99, 95.0% CI = 1.33–589.54, p = 0.03; OR = 424.73, 95.0% CI = 22.36–8066.18, p < 0.01; OR = 18.04, 95.0% CI = 4.80–67.74, p < 0.01; OR = 32.55, 95.0% CI = 3.16–335.29, p < 0.01).

3.3. Laboratory Examination

The results of 15 nasal swab samples tested in the field were negative for 2019-nCoV antigen. The results of 15 nasopharyngeal swab samples sent to the laboratory of the Gansu Provincial Center for Disease Control and Prevention were positive for HMPV (80.00%, 12/15), followed by HRV (53.33%, 8/15) and RSV (13.33%, 2/15) in 2 cases. Among these, six samples exhibited coinfections of HMPV and HRV, and one sample showed a combination of HMPV and RSV, resulting in a coinfection rate of 46.67% (Figure 2).
Sanger sequencing revealed three complete HMPV G gene sequences, designated as Gasu-2023-36, Gasu-2023-37, and Gasu-2023-41. Blast comparisons on GenBank classified them within the A2c sub-lineage. The three HMPV sequences were used to construct a phylogenetic tree with 59 HMPV reference strains from around the globe, as shown in Figure 3A, indicating the A2c genotype for HMPV in this outbreak. Additionally, five complete HRV VP1 gene sequences were obtained, named Gasu-2023-A3, Gasu-2023-A5, Gasu-2023-A6, Gasu-2023-A7, and Gasu-2023-A8. Blast comparisons on GenBank confirmed these as HRV-A type. A phylogenetic tree in Figure 3B illustrates the HRV infection in this outbreak as type A. Due to low viral load, gene typecasting could not be performed because an RSV sequence was not obtained in this study.

4. Discussion

The common viruses that can cause human respiratory tract infection are influenza virus (flu), respiratory syncytial virus (RSV), parainfluenza virus (PIV), rhinovirus (HRV), human metapneumovirus (HMPV), adenovirus (ADV), etc., among which there are also mixed infections [16,17]. Due to the challenge of establishing lasting immunity post-respiratory virus infection, populations are generally susceptible, with transmission intensifying in densely populated areas, resulting in outbreaks. Numerous studies report school outbreaks caused by HMPV, HRV, or RSV [18,19,20,21,22,23]. In this paper, the epidemiological investigation and etiology results of an upper respiratory tract infection outbreak in a kindergarten and primary school in Jingyuan County, Baiyin City, Gansu Province, were summarized and analyzed. We used the RT-PCT technique to type HRV genotypes, which has been proven effective when differentiating HRV from enteroviruses and human parechoviruses with close genetic relatedness [24,25,26]. The outbreak was a mixed infection of HMPV + HRV, HMPV + RSV, and an outbreak caused by RSV infection. In Gansu Province and Yinchuan City of Ningxia Hui Autonomous Region, cases of HMPV + HRV coinfection were found in the detection of viruses in hospitalized children. Still, cases of HMPV and RSV coinfection were rarely reported [27,28]. Studies have shown that HMPV and RSV are the primary pathogens of lower respiratory tract infections in children, often affecting infants under two years of age; in particular, children under six months of age frequently present with bronchopneumonia, the main clinical manifestations of which are fever, cough, wheezing, shortness of breath, vomiting, runny nose, sore throat, headache, etc. [18,19,20,21,22,23]. Conversely, HMPV infections in adults are often asymptomatic [19]. Monitoring data from Gansu suggests that HRV infections primarily affect young children [27], resulting in symptoms of lower respiratory tract infections, with HRV-positive children exhibiting more pronounced fever and vomiting symptoms [28]. HRV infections are associated with a higher risk of other diseases, often leading to respiratory symptoms such as fever and cough [28]. This study reveals that coinfections of HMPV, HRV, and RSV not only cause illnesses in children but also pose a risk to adult staff. Clinical presentations of cases in this outbreak primarily featured fever accompanied by respiratory symptoms, similar to cases caused by single HMPV, RSV, or HRV infections [29,30,31,32,33,34]. Additionally, the outbreak exhibited flu-like symptoms. The duration of the outbreak in the two schools was 36 and 46 days, longer than acute respiratory infection outbreaks caused by single infections of HMPV, RSV, or HRV [29,30,31,32,33,34]. HMPV and RSV infections typically show seasonal patterns, with higher prevalence in winter and spring [20,32]. In contrast, HRV infections exhibit a more noticeable seasonal peak in autumn [29], being the most common virus causing common colds in the summer and autumn. This outbreak, involving a mixture of HMPV, RSV, and HRV, occurred in May, with a deviation in the seasonal prevalence of HMPV and RSV and an early trend for HRV infections. The prolonged and unclear seasonal patterns suggest a multi-source exposure for this outbreak. Due to sibling relationships, the outbreak occurred in a primary school with significant contact history among different classes and schools. The cases demonstrated apparent class clustering, with higher attack rates in younger classes, particularly in the first grade and kindergarten, where students and preschoolers had higher attack rates than teaching staff.
Multiple logistic regression analysis revealed that, when compared to the control group and after adjusting for factors that may affect the continued spread of the epidemic, contact with flu-like cases within seven days, walking to and from school, window ventilation (<2 h), frequency of colds (≥3 times or =2 times), and participation in extracurricular training were identified as significant risk factors for the continued spread of the outbreak. These risk factors played an essential role in the sustained transmission of the epidemic. Contact with flu-like cases within seven days may lead to inhaling virus aerosols from infected individuals. Walking to and from school on a single road increased contact opportunity, and participation in extracurricular training may contribute to the spread between schools and grades. A higher frequency of colds in the past six months (≥3 times or =2 times) may indicate a weakened immune system due to repeated infections with COVID-19 or HRV, making individuals more susceptible to viral infections. Conversely, handwashing for 10–30 s or >30 s and family members covering their mouth and nose when sneezing were identified as protective factors. Handwashing and covering the mouth and nose when sneezing are effective measures to prevent respiratory infections and reduce the spread of viruses. Therefore, kindergartens, primary schools, and other crowded places should establish a working mechanism for the prevention and control of infectious diseases, the implementation of morning and afternoon inspection systems, and the follow-up of children absent due to illness through the morning and afternoon inspection so that children with abnormal symptoms are found in time and isolated for observation and treatment as soon as possible, or taken to the hospital for treatment. At the same time, it is essential to publicize the health knowledge of respiratory infectious diseases in schools and kindergartens, improve the personal protection awareness of students and staff, improve behavioral habits, and often ventilate classrooms and families, which can effectively reduce the occurrence of clustered epidemics. Timely reporting of outbreaks and continuous monitoring and research by disease control and medical institutions on the variation and evolution of infectious disease pathogens provide crucial data support for developing preventive and control strategies for respiratory pathogens such as HMPV, HRV, and RSV in the local area.
This outbreak identified 12 HMPV-positive samples, obtaining four complete sequences. Evolutionary tree analysis indicated that the HMPV in this outbreak belonged to the A2c genotype; it was most similar to the circulating strains in Beijing, Qingdao, and Henan (95.75–96.64%), among which the full gene sequence homology with the sequencing strains in Beijing 2018 and Qingdao 2019 reached 100%. The first strain of the HMPV A2c genotype was discovered in Okinawa, Japan in 2011 [35], with similar strains reported in China [36]. This outbreak identified eight HRV-positive samples, obtaining five complete VP1 region sequences. VP1 region sequencing revealed that the HRV in this outbreak belonged to species A, and evolutionary tree analysis suggested that the strains were most similar to those prevalent in Fujian (97.03–97.15%). Evolutionary tree analysis closely resembles prevalent strains in Fujian and other regions. This outbreak occurred in schools and kindergartens in a newly established immigrant community, with residents engaged in work from other areas. Based on epidemiological investigations and laboratory genetic sequencing results, the outbreak was caused by mixed infection of HMPV, HRV, and RSV viruses in local schools and kindergartens. The high genetic homology of the HMPV A2c strain suggests possible importation from Beijing, Qingdao, Henan, and other areas. Similarly, the high genetic homology of the HRV A strain suggests possible importation from Fujian and other areas. Both HMPV A2c and HRV A strains are currently dominant circulating strains. Further research is needed to investigate whether the combined infection with RSV increases pathogenicity and invasiveness compared to single infections, as well as to elucidate the mechanisms of triple infection. However, there are still some limitations to this study. Firstly, the outbreak occurred in northwest China and therefore the results may not be generalizable to other regions due to variations in climate, temperature, humidity, and human immune response following COVID-19 infection. Secondly, samples were not collected and tested for each case, which hindered accurate determination of pathogen type and corresponding clinical manifestations. Additionally, the long time interval between infection and the survey as well as inaccurate recall by some participants may have introduced bias into the analysis results. Finally, while RT-PCR was employed for virus identification in this study, critical values for experimental quantification were not utilized to confirm active infections. In conclusion, this paper presents preliminary evidence of mixed HMPV, HRV, and RSV infections in schools and kindergartens, but further research is warranted.

Author Contributions

Conceptualization, W.L. and H.Z.; methodology, H.Z.; formal analysis, W.L. and T.Z. (Tianshuo Zhao); investigation, L.Y., G.G., X.C. (Xiaoyan Che) and Z.Z.; writing—original draft preparation, W.L. and H.Z.; writing—review and editing, F.C. and T.Z. (Tongxia Zeng); visualization, T.Z. (Tianshuo Zhao) and X.C. (Xianming Cai); funding acquisition, W.L. and H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Science and Technology Major Project of China (2017ZX10103006), Sponsor (Health Commission of China); Natural Science Foundation of Gansu Province (20JR5RA142), Sponsor (Department of Science and Technology of Gansu Province); Baiyin Science and Technology Project (2022-2-46Y), Sponsor (Baiyin Science and Technology Bureau).

Institutional Review Board Statement

This study was approved by the ethics committee of Jingyuan County Center for Disease Control and Prevention, Gansu, China (approval number: 2023-05) in May 2023. Surveys and sampling among children in kindergartens and primary schools are permitted. All procedures were performed in compliance with relevant laws and institutional guidelines and have been approved by the appropriate institutional committee.

Informed Consent Statement

Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

All data and statistical code to reproduce the results in the manuscript are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. GBD 2015 Mortality and Causes of Death Collaborators. Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980–2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet 2016, 388, 1459–1544. [Google Scholar] [CrossRef] [PubMed]
  2. Madi, N.; Al-Nakib, W.; Mustafa, A.S.; Habibi, N. Metagenomic analysis of viral diversity in respiratory samples from patients with respiratory tract infections in Kuwait. J. Med. Virol. 2018, 90, 412–420. [Google Scholar] [CrossRef] [PubMed]
  3. Li, Z.J.; Zhang, H.Y.; Ren, L.L.; Lu, Q.B.; Ren, X.; Zhang, C.H.; Wang, Y.F.; Lin, S.H.; Zhang, X.A.; Li, J.; et al. Etiological and epidemiological features of acute respiratory infections in China. Nat. Commun. 2021, 12, 5026. [Google Scholar] [CrossRef] [PubMed]
  4. Tsukagoshi, H.; Ishioka, T.; Noda, M.; Kozawa, K.; Kimura, H. Molecular epidemiology of respiratory viruses in virus-induced asthma. Front. Microbiol. 2013, 4, 278. [Google Scholar] [CrossRef] [PubMed]
  5. Windsor, W.J.; Lamb, M.M.; Dominguez, S.R.; Mistry, R.D.; Rao, S. Clinical characteristics and illness course based on pathogen among children with respiratory illness presenting to an emergency department. J. Med. Virol. 2022, 94, 6103–6110. [Google Scholar] [CrossRef] [PubMed]
  6. Huang, Y.; Hua, J.; Wang, D.; Chen, L.; Zhang, J.; Zhu, H.; Tian, J.; Zhang, T.; Zhao, G. Risk factors of respiratory syncytial virus infection among pediatric influenza-like illness and severe acute respiratory infections in Suzhou, China. J. Med. Virol. 2018, 90, 397–404. [Google Scholar] [CrossRef] [PubMed]
  7. Bouzid, D.; Hadad, O.; Bertine, M.; Houhou-Fidouh, N.; Mirand, A.; Duval, X.; Bunel, V.; Borie, R.; Lucet, J.C.; Descamps, D.; et al. Rhinoviruses: Molecular diversity and clinical characteristics. Int. J. Infect. Dis. 2022, 118, 144–149. [Google Scholar] [CrossRef] [PubMed]
  8. Nichol, K.L.; D’Heilly, S.; Ehlinger, E. Colds and influenza-like illnesses in university students: Impact on health, academic and work performance, and health care use. Clin. Infect. Dis. 2005, 40, 1263–1270. [Google Scholar] [CrossRef] [PubMed]
  9. Jacobs, S.E.; Lamson, D.M.; St George, K.; Walsh, T.J. Human rhinoviruses. Clin. Microbiol. Rev. 2013, 26, 135–162. [Google Scholar] [CrossRef]
  10. The Common Cold, Again. JAMA 2023, 330, 2400. [CrossRef]
  11. Shapiro, D.; Bodinayake, C.K.; Nagahawatte, A.; Devasiri, V.; Kurukulasooriya, R.; Hsiang, J.; Nicholson, B.; De Silva, A.D.; Østbye, T.; Reller, M.E.; et al. Burden and Seasonality of Viral Acute Respiratory Tract Infections among Outpatients in Southern Sri Lanka. Am. J. Trop. Med. Hyg. 2017, 97, 88–96. [Google Scholar] [CrossRef] [PubMed]
  12. Li, Y.; Zhai, Y.; Lin, Y.; Lu, C.; He, Z.; Wu, S.; Yang, C.; Zhou, R.; Liu, W.; Chen, D. Epidemiology of respiratory syncytial virus in hospitalized children with community-acquired pneumonia in Guangzhou: A 10-year study. J. Thorac. Dis. 2023, 15, 967–976. [Google Scholar] [CrossRef] [PubMed]
  13. Chadha, M.; Prabhakaran, A.O.; Choudhary, M.L.; Biswas, D.; Koul, P.; Kaveri, K.; Dar, L.; Mamta, C.S.; Jadhav, S.; Bhardwaj, S.D.; et al. Multisite surveillance for influenza and other respiratory viruses in India: 2016–2018. PLOS Glob. Public Health 2022, 2, e0001001. [Google Scholar] [CrossRef] [PubMed]
  14. Schuster, J.E.; Williams, J.V. Human Metapneumovirus. Microbiol. Spectr. 2014, 2, 5. [Google Scholar] [CrossRef] [PubMed]
  15. Wildenbeest, J.G.; Billard, M.N.; Zuurbier, R.P.; Korsten, K.; Langedijk, A.C.; van de Ven, P.M.; Snape, M.D.; Drysdale, S.B.; Pollard, A.J.; Robinson, H.; et al. The burden of respiratory syncytial virus in healthy term-born infants in Europe: A prospective birth cohort study. Lancet Respir. Med. 2023, 11, 341–353. [Google Scholar] [CrossRef] [PubMed]
  16. Hodinka, R.L. Respiratory RNA Viruses. Microbiol. Spectr. 2016, 4, 4. [Google Scholar] [CrossRef] [PubMed]
  17. Yu, D.S.; Ren, L.L.; Chen, J.H.; Wang, P.; Qiao, R.J.; Kang, Q.; Zhang, R.X.; Ni, F.A. Etiology of acute lower respiratory tract infections among hospitalized children in Baiyincity, Gansu province of China. Chin. J. Viral. Dis. 2017, 7, 6. [Google Scholar]
  18. Chen, Q.X.; Luo, J.W.; Xu, S.J.; Zhou, S.Q.; Peng, W.J.; Lai, X.L. An outbreak of acute respiratory disease caused by human metapneumovirus. J. Tmp. Med. 2020, 20, 257–260. [Google Scholar] [CrossRef]
  19. Xiao, X.C.; Wu, L.P.; Wei, L.Y.; Liu, H.; Chen, Z.Q.; Lu, J.Y. An outbreak of respiratory tract infection caused by human metapneumovirus among adults in Guangzhou. Chin. J. Epidemiol. 2013, 34, 943–944. [Google Scholar] [CrossRef]
  20. Liu, X.L.; Luo, X.J.; Feng, Z.G.; Liu, C.Y.; Bao, Y.M.; Zheng, Y.J. Clinical analysis of 558 hospitalized children with human metapneumovirus pneumonia. J. Clin. Pediatr. 2023, 41, 459–463. [Google Scholar] [CrossRef]
  21. Wang, H.B.; Lei, N.; Ma, J.X.; Sun, L.L. Pathogen spectrum of influenza, like illness and infection status of human metapneumovirus in Chaoyang district of Beijing. Cap. J. Public Health 2020, 14, 259–261. [Google Scholar] [CrossRef]
  22. Liu, L.Z.; Chen, H.B.; Liang, J.; Luo, J.W.; Zhou, S.Q. Epidemiological investigation and risk factors of an outbreak of human metapneumovirus in a primary school. J. Med. Pest. Control 2021, 37, 265–268. [Google Scholar]
  23. Liu, B.H.; Peng, E.L.; Peng, C.W. Investigation on an outbreak of human metapneumovirus respiratory tract infection in a primary school of Huzhou City. Chin. J. Sch. Health 2022, 43, 120–122+128. [Google Scholar] [CrossRef]
  24. Liu, B. Universal PCR Primers Are Critical for Direct Sequencing-Based Enterovirus Genotyping. J. Clin. Microbiol. 2017, 55, 339–340. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, B.M.; Mulkey, S.B.; Campos, J.M.; DeBiasi, R.L. Laboratory diagnosis of CNS infections in children due to emerging and re-emerging neurotropic viruses. Pediatr. Res. 2024, 95, 543–550. [Google Scholar] [CrossRef] [PubMed]
  26. Liu, B.; Forman, M.; Valsamakis, A. Optimization and evaluation of a novel real-time RT-PCR test for detection of parechovirus in cerebrospinal fluid. J. Virol. Methods 2019, 272, 113690. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, S.; Mao, N.Y.; Yu, D.S.; Huang, G.H.; Li, X.X.; Li, H.Y.; Li, B.D. Characterization of Human Rhinovirus in Children with Acute Respiratory Infections inGansu Province during 2011. Chin. J. Virol. 2013, 29, 273–279. [Google Scholar] [CrossRef]
  28. Hei, Z.P.; Guo, J.H.; Ji, K.; Yan, Y.; Wang, C.L.; Li, F.; Sun, Y.N. Epidemiological and clinical characteristics of human rhinovirus infection in children hospitalized due to lower respiratory tract infections in Yinchuan, China. J. Chongqing Med. Univ. 2023, 48, 1337–1343. [Google Scholar] [CrossRef]
  29. Yang, X.H.; Yu, X.F.; Zhang, C.Y.; Wang, F.; Zhu, L.; Jiang, L.F.; Wang, J. Investigation on outbreaks of acute respiratory tract infection caused by respiratory syncytial virus in kindergartens in Hangzhou. Chin. J. Sch. Health 2022, 43, 142–145. [Google Scholar] [CrossRef]
  30. Yu, P.; Cheng, H.J.; Xiong, Y.; Fu, W.J.; Shi, Y.; Yuan, H. An outbreak of respiratory syncytial virus B infection in a private school. Chin. J. Epidemiol. 2013, 34, 2. [Google Scholar]
  31. Zhao, G.B.; Yi, G.G.; Yi, X.M.; Yang, G.; Hu, S.H.; Mei, Y.F.; Dong, F. Investigation of a fever outbreak with respiratory syncytial virus. Chin. J. School Doctor 2007, 21, 1. [Google Scholar]
  32. Yan, T.; Qin, M.; Yang, X.X.; Feng, H.Y.; Yu, X.X.; Yang, Y.; Yu, H. An investigation and pathogen identification of respiratory syncytial virus outbreak in a school. Mod. Prev. Med. 2021, 48, 4. [Google Scholar]
  33. Zhou, P.; Qin, X.; Tan, Z.Y.; Chen, H.B. Epidemiological features and clinical manifestations of children with respiratory syncytial virus of subtype A and B. China Med. Her. 2016, 13, 4. [Google Scholar]
  34. Wang, L.L.; Liu, Z.; Peng, H.Y.; Wu, S.T.; Yang, Q.L.; Li, L.P. Epidemic characteristics of respiratory syncytial virus infection in children before and after the COVlD-19 outbreak. Pract. Prev. Med. 2021, 28, 1487–1489. [Google Scholar]
  35. Nidaira, M.; Taira, K.; Hamabata, H.; Kawaki, T.; Gushi, K.; Mahoe, Y.; Maeshiro, N.; Azama, Y.; Okano, S.; Kyan, H.; et al. Molecular epidemiology of human metapneumovirus from 2009 to 2011 in Okinawa, Japan. Jpn. J. Infect. Dis. 2012, 65, 337–340. [Google Scholar] [CrossRef]
  36. Yi, L.; Zou, L.; Peng, J.; Yu, J.; Song, Y.; Liang, L.; Guo, Q.; Kang, M.; Ke, C.; Song, T.; et al. Epidemiology, evolution and transmission of human metapneumovirus in Guangzhou China, 2013–2017. Sci. Rep. 2019, 9, 14022. [Google Scholar] [CrossRef]
Figure 1. Time distribution of febrile respiratory syndrome cases in kindergartens/primary schools in Gansu Province in 2023. (A) The cases in kindergartens; (B) the cases in primary schools.
Figure 1. Time distribution of febrile respiratory syndrome cases in kindergartens/primary schools in Gansu Province in 2023. (A) The cases in kindergartens; (B) the cases in primary schools.
Viruses 16 00540 g001
Figure 2. Coinfection and detection of respiratory viruses. RSV infection is indicated by blue, HRV by orange, and HMPV by green. The overall detection of each pathogen is displayed in the lower left corner, while the combined and single pathogen detections are shown in the upper right corner.
Figure 2. Coinfection and detection of respiratory viruses. RSV infection is indicated by blue, HRV by orange, and HMPV by green. The overall detection of each pathogen is displayed in the lower left corner, while the combined and single pathogen detections are shown in the upper right corner.
Viruses 16 00540 g002
Figure 3. The phylogenetic tree of human metapneumovirus (HMPV) strains and human rhinovirus (HRV). (A) The phylogenetic tree of HMPV strains is based on G genes (B). The phylogenetic tree of HRV strains is based on VP1 genes.
Figure 3. The phylogenetic tree of human metapneumovirus (HMPV) strains and human rhinovirus (HRV). (A) The phylogenetic tree of HMPV strains is based on G genes (B). The phylogenetic tree of HRV strains is based on VP1 genes.
Viruses 16 00540 g003
Table 1. Main clinical manifestations of febrile respiratory syndrome in a township of Gansu Province in 2023.
Table 1. Main clinical manifestations of febrile respiratory syndrome in a township of Gansu Province in 2023.
SymptomsNo. Cases = 110Proportion (%)
Cough9485.45
Sore throat6760.91
Sneezing6660
Nasal congestion6054.55
Runny nose6054.55
Fever (°C)5953.64
Low-grade fever (37.5–38.0)87.27
Moderate fever (38.1–39.0)4440
High fever (≥39.1)76.37
Fatigue5953.64
Loss of appetite5751.82
Wheezing4641.82
Enlarged tonsils3531.82
Abdominal pain3430.91
Headache3027.27
Nausea3027.27
Chills2724.55
Phlegm production2220
Vomiting2018.18
Muscle aches1311.82
Diarrhea98.18
Table 2. Grade distribution of febrile respiratory syndrome cases in a township of Gansu Province in 2023.
Table 2. Grade distribution of febrile respiratory syndrome cases in a township of Gansu Province in 2023.
SchoolClassClass SizeNo. CasesIncidence (%)
MaleFemaleMaleFemale
Primary SchoolGrade 1878586.67
Grade 281241175.00
Grade 3853561.54
Grade 4342142.86
Grade 5791218.75
Grade 6674569.23
Teaching Staff1675126.09
KindergartenJunior Class959492.86
Middle Class1199990.00
Senior Class71571181.82
Teaching Staff451344.44
Table 3. Multivariate logistic regression analysis of respiratory virus coinfection.
Table 3. Multivariate logistic regression analysis of respiratory virus coinfection.
VariableβSEWaldpOR95% CI
InferiorSuperior
COVID-19 infection history (2)0.670.680.97 0.32 1.950.527.33
Contact with flu-like cases in the past seven days (2)1.710.627.71 0.01 1.181.051.60
Location of contact with flu-like cases (2)−0.74 0.98 0.57 0.45 0.480.07 3.24
Occupation--0.10 0.95 ---
Student−1.563.06 0.26 0.61 0.21 0.00 85.35
Preschool children2.30 2.71 0.72 0.40 10.00 0.05 818.88
Age - 0.63 0.43 - - -
7–14 years old (2)1.74 1.91 0.83 0.36 5.71 0.13 242.28
3–6 years old (3)2.33 2.76 0.79 0.40 9.97 0.06 818.88
Handwashing time--10.47 0.01 ---
Handwashing time (1)−2.88 1.13 6.57 0.01 0.06 0.01 0.51
Handwashing time (2)−2.89 0.92 9.92 0.00 0.06 0.01 0.34
Covering mouth and nose when sneezing0.29 0.67 0.19 0.66 1.34 0.36 5.02
Washing hands immediately after sneezing−0.94 0.77 1.52 0.22 0.39 0.09 1.75
Handwashing after class0.74 0.73 1.02 0.31 2.10 0.50 8.86
Thumb-sucking--0.25 0.88 ---
Thumb-sucking (1)−0.32 1.04 0.09 0.76 0.73 0.10 5.59
Thumb-sucking (2)0.34 1.08 0.10 0.75 1.41 0.17 11.62
Rubbing eyes--0.42 0.81 ---
Rubbing eyes (1)−0.65 1.39 0.22 0.64 0.52 0.03 7.94
Rubbing eyes (2)−0.37 0.67 0.32 0.57 0.69 0.19 2.55
Clean hands and nails (1)−0.14 0.92 0.02 0.88 0.87 0.14 5.29
Mode of transportation--10.74 0.01 ---
Mode of transportation (1)−1.53 1.77 0.75 0.39 0.22 0.01 6.96
Mode of transportation (2)−0.56 1.80 0.10 0.76 0.57 0.02 19.34
Mode of transportation (3)2.60 1.00 6.83 0.01 13.47 1.92 94.75
Duration of stay--1.10 0.58 ---
Duration of stay (1)−2.86 3.49 0.67 0.41 0.06 0.00 54.11
Duration of stay (2)−0.47 0.64 0.54 0.46 0.62 0.18 2.18
Wearing masks−0.25 0.64 0.15 0.70 0.78 0.22 2.75
Covering mouth and nose when sneezing−1.45 0.72 4.11 0.04 0.23 0.06 0.95
Ventilating by opening windows (<2 h)3.33 1.55 4.59 0.03 27.99 1.33 589.54
Prophylactic medication0.50 0.71 0.49 0.49 1.65 0.41 6.68
Sleep duration−0.68 1.57 0.19 0.67 0.51 0.02 11.01
Frequency of colds--26.65 0.00 ---
Frequency of colds (1)6.05 1.50 16.23 0.00 424.73 22.36 8066.18
Frequency of colds (2)2.89 0.68 18.36 0.00 18.04 4.80 67.74
Extracurricular training3.48 1.19 8.57 0.00 32.55 3.16 335.29
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, W.; Zhang, H.; Zhao, T.; Cai, X.; Yang, L.; Gao, G.; Che, X.; Zhu, Z.; Zeng, T.; Cui, F. Molecular Epidemiology and Clinical Characteristics of an Outbreak on Respiratory Virus Coinfection in Gansu, China. Viruses 2024, 16, 540. https://doi.org/10.3390/v16040540

AMA Style

Liu W, Zhang H, Zhao T, Cai X, Yang L, Gao G, Che X, Zhu Z, Zeng T, Cui F. Molecular Epidemiology and Clinical Characteristics of an Outbreak on Respiratory Virus Coinfection in Gansu, China. Viruses. 2024; 16(4):540. https://doi.org/10.3390/v16040540

Chicago/Turabian Style

Liu, Wu, Hui Zhang, Tianshuo Zhao, Xianming Cai, Liguo Yang, Genxia Gao, Xiaoyan Che, Zhenhong Zhu, Tongxia Zeng, and Fuqiang Cui. 2024. "Molecular Epidemiology and Clinical Characteristics of an Outbreak on Respiratory Virus Coinfection in Gansu, China" Viruses 16, no. 4: 540. https://doi.org/10.3390/v16040540

APA Style

Liu, W., Zhang, H., Zhao, T., Cai, X., Yang, L., Gao, G., Che, X., Zhu, Z., Zeng, T., & Cui, F. (2024). Molecular Epidemiology and Clinical Characteristics of an Outbreak on Respiratory Virus Coinfection in Gansu, China. Viruses, 16(4), 540. https://doi.org/10.3390/v16040540

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop