Abstract
Background Under the pressure of the outbreak of respiratory Human Adenovirus (HAdV) infections in children in Northern Vietnam in the end of 2022, this study was initiated to identify the HAdV subtype(s) responsible for the outbreak in relation to the clinical features of the patients and examine the risk factors of more severe cases.
Methods The study was conducted on pediatric patients tested positive with HAdV using multiplex real- time PCR between October and November 2022. Nasal swab samples were used for sequencing to identify HAdV subtypes and clinical data were collected retrospectively.
Results Among 97 successfully sequenced samples, the predominant subtypes were HAdV-B3 (83%), HAdV-B7 (16%) and HAdV-C2 (1%). Lower respiratory manifestations were found in 25% of patients (5% with severe pneumonia). There was no significant association between HAdV type and clinical features except that those infected with HAdV type 3 exhibited higher WBC and neutrophil % (p<0.001). Co- infection of HAdV with ≥1 other respiratory viruses or bacteria was found in 70.8% of those with lower respiratory illnesses (OR (95%CI); p-value vs. those without =5.21 (1.60, 19.36); 0.0084 after adjusting for age at hospital visit, sex, birth delivery method, day of disease at hospital visit), and in 100% of those with severe pneumonia vs. 33% of those without (p=0.005).
Conclusion HAdV-B3 and HAdV-B7 were predominant in the outbreak. Co-infection of HAdV together with other respiratory viruses or bacteria was a strong risk factor for lower respiratory tract illnesses and severe pneumonia. The findings advocate the advantages of multi-factor microbial panels for the diagnosis and prognosis of respiratory infections in children.
Introduction
Human Adenoviruses (HAdV) can cause many common diseases such as upper respiratory infection, conjunctivitis, gastrointestinal disorders, etc. HAdV consists of many serotypes, which are responsible for different disease contexts (1).
There had been some sporadic HAdV outbreaks with respiratory illnesses, conjunctivitis, and gastrointestinal disorders worldwide before the COVID-19 era. Respiratory illnesses were reported in most of these studies, and HAdV infections with severe respiratory illnesses in children had been reported in some studies (3–26). The most prevalent respiratory HAdV reported in children in Taiwan were HAdV- B genotype 3 and HAdV-C (HAdV-2), and 7 (15,24,26), in China were HAdV-B3 and HAdV-C2, and 7 (9,10,12,17,27), in Japan were HAdV type 3, 2, and 1 (21), in Argentina were HAdV-B7 and HAdV-B3 (16,20), in Spain were HAdV-3, HAdV-6, and HAdV-5 (25), in Egypt were type 3 and 7 (28), in Palestine were HAdV-C1, HAdV-C2, HAdV-B3 and HAdV-C5 (6). HAdV type 7, 3, and 2 were the most common types which were associated with severe respiratory illnesses in children, especially in those with underlying diseases or immunocompromised conditions (10–12,16–19). During the early stages of the COVID-19 pandemic with social distancing until 2021, respiratory tract infections in children including HAdV infections exhibited a remarkable decrease in occurrence (14,29–31).
In Vietnam, there have been some reports of adenoviral conjunctivitis (32,33), and gastroenteritis (13,34,35). HAdV is among the less common viral causes of seasonal respiratory infections as compared to many other more common viruses such as influenza, parainfluenza, respiratory syncytial virus (RSVs), enteroviruses, rhinovirus, etc (36–43). HAdV is among the co-infection factors in patients with severe pneumonia (39). HAdV pneumonia especially type 7 has been reported to be among the causes of death in patients following measles infection (44,45). There has not been any information regarding respiratory HAdV outbreaks in children in Vietnam until the first half of 2022.
Since the second half of 2022, HAdV had emerged as the cause of the outbreak of respiratory infections in children in Northern Vietnam. Whilst most of the cases were mild, there had been a considerable number of cases with severe complications such as pneumonia, acute respiratory distress. which required intensive treatment and prolonged hospitalization. In addition, quite a few deaths attributed to HAdV infection had been reported. While the COVID-19 pandemic had been somewhat under control, the outbreak of respiratory HAdV infection in children in Northern Vietnam with a large number of cases and an alarming number of severe cases and deaths was unprecedented and thus especially worrisome and required immediate actions.
Up till now, in Northern Vietnam, the diagnosis of HAdV infection was confirmed by positive PCR of nasal swabs performed by a few centers. However, molecular typing of HAdV had not been performed to identify the HAdV subtype(s) responsible for the outbreak, especially in severe or fatal cases. As such, this study aimed to perform molecular typing of HAdV and examine the relationship between HAdV types and clinical contexts of children with respiratory HAdV infections in Northern Vietnam to understand the characteristics of the pathogen in relation to host characteristics. The study also examined potential risk factors for the severity of the diseases to help orientate the treatment and prognosis of the patients.
Materials & Methods
Sample Collection
This study was done in Vinmec Times City International Hospital, a private general hospital in Hanoi Vietnam which receives patients mostly from Hanoi and surrounding provinces. The study was approved by the hospital’s ethical committee before being carried out (Reference number: 152/2022/CN-HDDD VMEC).
Nasopharyngeal wash or swab samples of pediatric inpatients and outpatients with respiratory symptoms suspected of Adenoviral infection were collected in a 15.0 ml centrifuge tube and transported to the Department of Microbiology. Specimens were stored at 4°C and partly used to extract RNA by QIAamp Viral RNA Kits (QIAGEN, German) within 24 hours after collection. Extracted RNA samples were tested for HadV and 6 other respiratory pathogens (human enterovirus (HEV), human metapneumovirus (hMPV) and parainfluenza virus 1-4 (PIV1-4)) using multiplex real-time PCR Allplex™ Respiratory Panel 2 (Avicalab Diagnostics, Kenya) according to the manufacturer’s instruction (46,47). Samples with quantification cycle (Ct) value < 38 were considered HAdV-positive.
The leftover nasopharyngeal samples (NS) of those with HAdV-positive PCR were stored and transported to Vinmec Medical Genetic Department for further processing to perform HAdV molecular typing. To ensure a sufficient DNA concentration for sequencing and to minimize possible artifacts due to contamination, only samples with Ct for HAdV <30 were selected for molecular typing.
For pediatric patients with HAdV-positive PCR, clinical data regarding medical history especially prior to COVID infection, clinical manifestations, lab tests (hematology, biochemistry and also nasal-pharyngeal swab culture, PCR panel for 7 respiratory bacteria, blood Mycoplasma pneumoniae Antibody test, influenza antigen rapid test if available), treatments, outcomes, etc were retrospectively collected by Vinmec Pediatric Center. These data were entered and managed in RedCap (https://www.project-redcap.org/).
All patients included in this study had respiratory symptoms such as cough, runny nose, etc. Patients were classified as having lower respiratory illnesses if having lower respiratory symptoms such as rales, rhonchi, wheezing, etc or having lesions in chest x-ray. Patients were classified as having severe pneumonia if having signs of pneumonia clinically and/or in chest x-ray and respiratory distress.
Viral co-infection was defined if a patient with HAdV-positive PCR was also positive with at least one virus in the PCR panel or in the influenza antigen rapid test.
Bacterial co-infection was defined if a patient with HAdV-positive PCR was also positive with at least one bacteria in nasopharyngeal fluid culture or in a PCR panel of 7 respiratory bacteria (H. influenzae, S. pneumonia, S. aereus, M. catarrhalis, M. pneumoniae, C. pneumoniae, L. pneumoniae) or in blood Mycoplasma pneumoniae antibody test.
Any co-infection was defined if a patient had either a viral or bacterial co-infection as described above.
HAdV molecular typing
The hyper-variable regions 1-6 (HRV1-6) of HadV exon are known to contain type-specific epitopes, therefore, these sequences are commonly used to identify HAdV subtypes responsible for the diseases. HAdV molecular typing was performed on 97 HAdV-positive samples. Viral nucleic acids were extracted from NS using QIAamp DNA Blood Mini Kit (QIAGEN, Germany). Following the protocol recommended by CDC, the hexon of HadV was amplified using nested PCR targeting the HVR1-6 and amplicons were expected to have the size between 764-896 base pairs (bp) (48). The outer primer pairs used for 1st round of nested-PCR were AdhexF1 (5′- TICTTTGACATICGIGGIGTICTIGA-3′) and AdhexR1 (5′- CTGTCIACIGCCTGRTTCCACA-3′). The inner primer pairs used for 2nd round of nested-PCR were AdhexF2 (5′- GGYCCYAGYTTYAARCCCTAYTC-3′) and AdhexR2 (5′- GGTTCTGTC ICCCAGAGARTCIAGCA-3′). Nested-PCR aws performed using a final concentration of 1x GoTaq Green Master Mix (Promega), 0.4 μM of each primer, 2.0 μl of extracted viral nucleic acid or 0.2 μl of 1st nested-PCR product, and made up to a final volume of 25 μl with double distilled water. The PCR thermal program for both rounds of nested-PCR reactions were comprised an initial denaturation step at 94°C for 2 min, followed by 11 touched-down cycles of 95°C for 30 sec, 55°C (decrease 1°C after each cycle) for 45 sec, 72°C for 1 min, followed by 20 thermal cycles of 95°C for 30 sec, 45°C for 45 sec, 72°C for 1 min, and a 5 min final extension at 72°C. The size of PCR products were checked by DNA electrophoresis analysis on 1% agarose gels stained with RedSafe DNA Stain (20,000 X) (Chembio, USA). Only samples with a detectable appropriate PCR product from nested-PCR were used for Sanger Sequencing of the HVR1-6 of hexon gene.
Each nested-PCR product was sequenced using both forward AdhexF2 and reverse AdhexR2 primers. Sequencing PCR reactions were performed in a final volume of 10 μl comprising 0.5 μl of BigDye Terminator v3.1 Cycle Sequencing (ABI), 2.0 μl of Sequencing Buffer 5X, 0.85 μl (10 μM) of primer, 0.2μl 2nd nested-PCR product, and distilled water (up-to-10 μl). The sequencing thermal cycle for sequencing reactions were comprised of an initial denaturation step at 96°C for 2 min, followed by 25 thermal cycle of 95°C for 10 sec, 50°C for 10 sec, and 60°C for 2 min. Cycle sequencing products were purified by Bigdye X Terminator purification kit before analysis on th Applied Biosystems 3500 Dx Genetic Analyzer.
Data analysis
Sequence data and phylogenetic analysis
The sequence data were used to identify the subtype of HAdV by Basic Local Alignment Search Tool (BLAST) search (NCBI). The sequences of HVR1-6 of hexon gene from all samples were aligned with related reference strains obtained from Genbank (reference strains are listed in Supplementary Table 1) using ClustalW. The robustness of the phylogenetic tree was constructed by the neighbor-joining method with bootstrap analysis (n = 1000) by iTOL v6 and Adobe Illustrator software. The cut-off value is <80%.
Analysis of HAdV types and clinical contexts
HAdV types and other clinical features were compared between disease contexts such as lower respiratory illnesses vs. no lower respiratory illnesses or severe pneumonia vs no severe pneumonia. The purpose was to examine if HAdV type or any other clinical features were associated with the severity of HAdV infection. Clinical features were also compared between HAdV types to examine if any clinical characteristic was associated with HAdV types.
For comparison between groups for initial association exploration, Kruskall-Wallis was used for continuous variables and Fisher’s exact test was used for categorical variables.
Potential risk factors for lower respiratory illnesses or severe pneumonia which were clinically relevant or statistically significant from the above exploration were further examined using univariate and multivariate logistic regression models. Unadjusted and adjusted odds ratio (OR) and 95% confidence interval (95%CI) by profile likelihood and p-values from Wald test were reported.
Data analysis was done using R statistical software version 4.1 (https://www.r-project.org). All statistical tests were 2-sided with a significant level alpha of 0.05.
Results
Patient characteristics and clinical features
From October 13th to November 9th, 2022, 138 patients tested positive for HAdV using a Q-PCR assay from nasal swabs. Among these, 97 patients with Ct for HAdV <30 were selected for molecular typing. Of these 97 samples, the median (IQR) Ct of HAdV was 23.6 (21.9, 25.6).
Patient characteristics overall and comparison between those with and without lower respiratory illnesses and those with and without severe pneumonia are summarized in Table 1. About 68% of patients resided in Hanoi and 32% resided in surrounding provinces. About 29% reported history of prior COVID-19 infection. About 4% reported history of allergy and about 7% reported history of prior severe respiratory diseases. There were 73 (75%) patients with only upper respiratory symptoms whilst there were 24 (25%) with lower respiratory manifestations of which 5 (5.1%) had severe pneumonia. There were no differences in demographic and previous history of respiratory disease including COVID19 infection between those with and without lower respiratory illnesses in these characteristics (Table 1, Figure 1, Figure 2).
The mean (95%CI) age of all patients was 2.7 (2.4, 3) years and those with lower respiratory illnesses were a bit younger than those without (mean (95%CI) = 2.3 (1.6, 2.9) vs. 3.1 (2.6, 3.6) respectively) and those with severe pneumonia were the youngest (mean (95%CI) = 1.8 (0, 3.8)). Males accounted for 55% of all patients and there were no gender bias in those with and without lower respiratory illnesses whereas 4/5 (80%) of those with severe pneumonia were males. About 43% of patients were born by Caesaren (C)- section and the percentage was similar between those with and without lower respiratory illnesses while 4/5 (80%) of those with severe pneumonia were born by C_section (p=0.051). There was no difference in BMI, gestational age, birth weight and duration of breastfeeding between those with and without lower respiratory illnesses or between those with and without severe pneumonia (Table 1).
Day of illness at hospital visit of those with lower respiratory illnesses was later than those without (mean (95%CI) = 4.7 (3.7, 5.6) days vs. 3.3 (2.8, 3.8)) days (p=0.002)). Most of the patients had fever (94%) and cough (84%) at disease onset. Most of the patients (79%) had high fever (>39 C) and long fever duration (mean (95%CI) = 6.4 (5.7, 7.1) days). Poor eating (as compared to usual) as reported by caregivers was found more in patients with lower respiratory illnesses than those without. There was no significant difference in symptoms at disease onset such as cough, red eyes, diarrhea, fever maximal temperature, fever duration, ear pain, etc as well as in hematology or biochemistry indexes between those with and without lower respiratory illnesses or between those with and without severe pneumonia. Almost all patients with lower respiratory illnesses (95%) had detectable abnormalities in chest X-rays as well as 16% of those without lower respiratory symptoms (Table 1, Figure 1 and Figure 2).
About 25% of the patients with lower respiratory illnesses received oxygen supplements, about 10% received IVIG, 44% received nebulizers, 24% received general steroids while none of the patients without lower respiratory illnesses received these interventions. Duration of hospitalization was longer for those with lower respiratory illnesses than those without (mean (95%CI) = 6.1 (4.7, 7.4) vs. 4.5 (3.8, 5.1) days) and longest in those with severe pneumonia (mean (95%CI) = 9.5 (4.9, 14.1) days). All patients except one were discharged without complication (Table 1, Figure 1 and 2).
The five patients with severe pneumonia were of a younger age and were predominantly male (80%), of which 80% were born by C-section and were all co-infected with at least one other respiratory viruses or bacteria (p<0.05 as compared to other patients). There was no other remarkable difference in other demographic, medical history or host characteristics as compared to other patients.
Co-infection of HAdV and other respiratory viruses or bacteria
Regarding viral co-infection, there were 19 (approx. 20%) samples tested positive for at least one other virus from our respiratory virus PCR panel of which 5(5.1%) were co-infected with at least 2 viruses. The number of samples co-infected with HEV, MPV, PIV1, PIV2, PIV3, PIV4 in the PCR panel were 7, 9, 4, 2, 0, 2 respectively. There were two samples co-infected with Influenza virus detected by rapid antigen test. In overall, there were 20 samples co-infected with at least 1 virus detected by either PCR or rapid antigen test. Regarding bacterial co-infection, there were 3, 11, 5 samples with positive culture for Streptococcus pneumoniae, Haemophilus influenzae, Moraxella catarrhalis respectively. There were 2 patients with blood antibody test positive for Mycoplasma pneumoniae. In overall, there were 20 patients co-infected with at least one bacterium detected by either nasopharyngeal fluid culture or blood antibody test. There were 5 patients co-infected with both at least one virus and at least one bacterium. In total, there were 35 patients co-infected with either virus or bacteria (Table 2).
The percentage of samples co-infected with at least one other respiratory virus was significantly higher in patients with lower respiratory illnesses vs. those without (12 (50%) vs. 4 (9%), p< 0.001) and especially higher in those with severe pneumonia (4/5 (80%), p=0.006). Co-infection with at least one respiratory bacterium was detected in 20 (20.6%) patients, with 7 of these patients (29.2%) experiencing lower respiratory illnesses (p=0.253) and in 3 (60.0%) with severe pneumonia (p=0.058). Any co-infection with respiratory viruses or bacteria was found in 35 (36.1%) of all patients, significantly higher in those with lower respiratory illnesses (17 (70.8%), p< 0.001) and especially higher in those with severe pneumonia (5 (100.0%), p=0.005) (Table 1, Table 2).
Risk factors for lower respiratory illnesses and severe pneumonia
Potential risk factors for lower respiratory illnesses or severe pneumonia which were clinically relevant and statistically significant from the above exploration were viral co-infection of HAdV with ≥1 respiratory viruses or bacterial co-infection of HAdV with ≥1 respiratory bacteria or any co-infection of HAdV with ≥1 respiratory viruses or bacteria, age at hospital visit, sex, birth delivery method, day of disease at hospital visit. The results from univariate logistic regression and multivariate logistic regression adjusting for these factors in the models for the comparison between those with lower respiratory illnesses or severe pneumonia vs. those without are summarized in Table 3. Co-infection of HAdV with ≥1 other respiratory viruses was remarkably significantly associated with either lower respiratory illnesses or severe pneumonia (unadjusted OR (95%CI); (p-value) for those with lower respiratory illnesses or severe pneumonia vs. those without were 11.14 (3.77, 36.05); <0.0001 and 19.00 (2.60, 385.24); 0.0105 respectively). After adjusting for age at hospital visit, sex, birth delivery method, day of disease at hospital visit, the association between co-infection and lower respiratory illnesses or severe pneumonia remained statistically significant (adjusted OR (95%CI); (p-value) for those with lower respiratory illnesses or severe pneumonia vs. those without were 7.4 (1.4, 52.1); 0.03 and 16.9 (1.8, 421.9); 0.03 respectively). Co- infection of HAdV with any respiratory viruses or bacteria was also remarkably significantly associated with lower respiratory illnesses (unadjusted OR (95%CI); (p-value) = 7.42 (2.75, 21.98); 0.0001 and adjusted OR (95%CI); (p-value) = 5.21 (1.60, 19.36); 0.0084).
HAdV types and the relationship between HAdV types and clinical features
Ninety-seven (98.98%) HAdV samples were genotyped successfully. BLAST results of 97 HAdV sequences indicated that 82 (84.5%) samples belonged to subtype B3, 14 (14.4%) samples belonged to subtype C7 and 1 (1.0%) samples belonged to sybtype C2. The phylogenetic analysis were applied to compare 97 hexon gene sequences to each other and 24 reference sequences downloaded from GenBank (Supplement Table 1). The phylogenetic tree of hexon gene confirmed the hexon sequences of HAdV samples belonged to three distinguishable clusters of HAdV-B3, -B7 and C2 (Figure 3).
Discussion
This study, to the best of our knowledge is probably the first providing a quick comprehensive snapshot of both clinical features and molecular typing of the respiratory HAdV infection outbreak in children in northern Vietnam in the end of 2022.
The predominant HAdV type in this study was species B type 3 (83%), followed by species B type 7 (16%) and species C type 2 (1%). This is similar to the findings of other studies on circulating or outbreaks of respiratory HAdV infections in children before the COVID-19 pandemic in Asia such as studies on circulating respiratory HAdV in Hong Kong in 2014 (49), Taiwan from 2002 to 2013 (15), Guangzhou China from 2017 to 2019 (27) or a study on febrile respiratory HAdV outbreaks in Hangzhou China in 2011 (17) or a study on successive respiratory HAdV outbreaks in Seoul, Korea from 1990 to 2000 (11). Novel genome types were identified during the outbreaks of lower respiratory tract HAdV infection in Seoul Korea and this 11-year study also suggested a genome type shift between successive outbreaks (11). Therefore, further studies examining the genome types or variants of HAdV may be useful.
There was no significant association between the 3 HAdV types identified and clinical features in this study, except that WBC and neutrophil numbers were higher in those infected with HAdV type 3 than with type 7. The proportion of severe cases seemed to have similar distribution between the more predominant HAdV type 3 and type 7. This is consistent with other previously published studies in which HAdV type 3 and 7 were the most common types among more severe respiratory HAdV infections (10–12,15–19).
Co-infections of HAdV with other respiratory viruses or bacteria stood out to be highly significantly associated with increased severity of respiratory diseases. About 70.8% of those with lower respiratory illnesses and 100% of those with severe pneumonia were co-infected with at least one respiratory viruses or bacteria while the percentage was about 24.7% in those without lower respiratory illnesses. This association remained statistically significant in a multivariate analysis adjusting for other potential risk factors which were either clinically relevant or also statistically associated with lower respiratory illnesses or severe pneumonia. This indicates that the association is both clinically meaningful and statistically robust. This remarkable association suggests that HAdV might be the cause of the outbreak but might not be the main cause of severe cases. Instead, the co-infection of other respiratory viruses or bacteria together with HAdV might increase the severity of the disease. This advocates the advantage of multi- factor respiratory PCR panel(s) for common respiratory viruses and probably also respiratory bacteria together with other microbial tests such as culture for the diagnosis and prognosis of children with respiratory infections. These results may be helpful for the diagnosis, treatment and prognosis of patients in future outbreaks of respiratory HAdV infection.
Although the sample size of this study is limited, it is similar to those of many other previously published studies analysing molecular typing for HAdV outbreaks (5,9,16,49–51). The study samples were from pediatric patients of a private general hospital in Hanoi, Vietnam and molecular typing was performed using samples with Ct for HAdV <30. As such, the study samples may not fully represent the general child population of the outbreak. Additionally, it should be noted that clinical data were collected retrospectively for samples with HadV-positive PCR and thus some clinical data features maybe incomplete, especially for outpatients. Moreover, molecular typing data may not be sufficient to understand the characteristics of HAdV causing the outbreak and thus further study to examine HAdV genome types or variants may be useful.
In summary, this study revealed that HAdV type 3 and type 7 were predominant in the current ongoing outbreak of respiratory HAdV infection in children in northern Vietnam and HAdV type might not be associated with the severity of the diseases. Instead, co-infection of HAdV together with other respiratory viruses or bacteria appeared to be a significant risk factor for lower respiratory tract illnesses and severe pneumonia. These findings advocate the advantages of multi-factor microbial panels for the diagnosis and prognosis of respiratory infections in children. These results may be helpful for the diagnosis, treatment, and prognosis of patients in future outbreaks of respiratory HAdV infection.
Data Availability
The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.
Notes
Author Contributions
D.D.N, L.T.T, H.T.T.T., and N.T.H. designed the study as well as drafted the initial manuscript. N.T.H. led the data management and statistical analysis. All authors provided (1) additional contributions to the conception or design of this work or the acquisition, analysis, or interpretation of data for the work; (2) drafting of the manuscript or revising it critically for important intellectual content; and (3) approval of the final version submitted. D.D.N., H.T.T.T., and N.T.H had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Financial support
This work was supported by the Research & Development fund of Vinmec High-tech Centre, Vinmec Healthcare System.