Molecular subtypes of respiratory Adenovirus infection outbreak in children in Northern Vietnam and risk factors of more severe cases

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.


49
Human Adenoviruses (HAdV) can cause many common diseases such as upper respiratory infection, 50 conjunctivitis, gastrointestinal disorders, etc. HAdV consists of many serotypes, which are responsible for 51 different disease contexts (1). 52 There had been some sporadic HAdV outbreaks with respiratory illnesses, conjunctivitis, and 53 gastrointestinal disorders worldwide before the COVID-19 era. Respiratory illnesses were reported in 54 most of these studies, and HAdV infections with severe respiratory illnesses in children had been reported 55 in some studies (3-26). The most prevalent respiratory HAdV reported in children in Taiwan were HAdV-56 B genotype 3 and HAdV-C (HAdV-2), and 7 (15,24,26), in China were HAdV-B3 and HAdV-C2, and 7 57 (9,10,12,17,27), in Japan were HAdV type 3, 2, and 1 (21), in Argentina were HAdV-B7 and HAdV-B3 58 (16,20), in Spain were HAdV-3, HAdV-6, and HAdV-5 (25), in Egypt were type 3 and 7 (28), in Palestine 59 were HAdV-C1, HAdV-C2, HAdV-B3 and HAdV-C5 (6). HAdV type 7, 3, and 2 were the most common types 60 which were associated with severe respiratory illnesses in children, especially in those with underlying 61 diseases or immunocompromised conditions (10-12, [16][17][18][19]. During the early stages of the COVID-19 62 pandemic with social distancing until 2021, respiratory tract infections in children including HAdV 63 infections exhibited a remarkable decrease in occurrence (14,(29)(30)(31). 64 In Vietnam, there have been some reports of adenoviral conjunctivitis (32,33), and gastroenteritis 65 (13,34,35). HAdV is among the less common viral causes of seasonal respiratory infections as compared 66 to many other more common viruses such as influenza, parainfluenza, respiratory syncytial virus (RSVs), 67 enteroviruses, rhinovirus, etc (36-43). HAdV is among the co-infection factors in patients with severe 68 pneumonia (39). HAdV pneumonia especially type 7 has been reported to be among the causes of death 69 in patients following measles infection (44,45). There has not been any information regarding respiratory 70 HAdV outbreaks in children in Vietnam until the first half of 2022. 71 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.   Mycoplasma pneumoniae antibody test. 121 Any co-infection was defined if a patient had either a viral or bacterial co-infection as described above. 122 HAdV molecular typing 123 The hyper-variable regions 1-6 (HRV1-6) of HadV exon are known to contain type-specific epitopes, 124 therefore, these sequences are commonly used to identify HAdV subtypes responsible for the diseases. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 23, 2023. ; https://doi.org/10.1101/2023.04.18.23288722 doi: medRxiv preprint size of PCR products were checked by DNA electrophoresis analysis on 1% agarose gels stained with 139 RedSafe DNA Stain (20,000 X) (Chembio, USA). Only samples with a detectable appropriate PCR product 140 from nested-PCR were used for Sanger Sequencing of the HVR1-6 of hexon gene. 141 Each nested-PCR product was sequenced using both forward AdhexF2 and reverse AdhexR2 primers. ClustalW. The robustness of the phylogenetic tree was constructed by the neighbor-joining method with 154 bootstrap analysis (n = 1000) by iTOL v6 and Adobe Illustrator software. The cut-off value is <80%. 155

Analysis of HAdV types and clinical contexts 156
HAdV types and other clinical features were compared between disease contexts such as lower 157 respiratory illnesses vs. no lower respiratory illnesses or severe pneumonia vs no severe pneumonia. The 158 purpose was to examine if HAdV type or any other clinical features were associated with the severity of 159 HAdV infection. Clinical features were also compared between HAdV types to examine if any clinical 160 characteristic was associated with HAdV types. 161 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 23, 2023. ; https://doi.org/10.1101/2023.04.18.23288722 doi: medRxiv preprint For comparison between groups for initial association exploration, Kruskall-Wallis was used for 162 continuous variables and Fisher's exact test was used for categorical variables. 163 Potential risk factors for lower respiratory illnesses or severe pneumonia which were clinically relevant or 164 statistically significant from the above exploration were further examined using univariate and 165 multivariate logistic regression models. Unadjusted and adjusted odds ratio (OR) and 95% confidence 166 interval (95%CI) by profile likelihood and p-values from Wald test were reported. 167 Data analysis was done using R statistical software version 4.1 (https://www.r-project.org). All statistical 168 tests were 2-sided with a significant level alpha of 0.05. 169

Patient characteristics and clinical features 172
From October 13th to November 9th, 2022, 138 patients tested positive for HAdV using a Q-PCR assay 173 from nasal swabs. Among these, 97 patients with Ct for HAdV <30 were selected for molecular typing. Of 174 these 97 samples, the median (IQR) Ct of HAdV was 23.6 (21.9, 25.6). 175 Patient characteristics overall and comparison between those with and without lower respiratory illnesses 176 and those with and without severe pneumonia are summarized in Table 1. About 68% of patients resided 177 in Hanoi and 32% resided in surrounding provinces. About 29% reported history of prior COVID-19 178 infection. About 4% reported history of allergy and about 7% reported history of prior severe respiratory 179 diseases. There were 73 (75%) patients with only upper respiratory symptoms whilst there were 24 (25%) 180 with lower respiratory manifestations of which 5 (5.1%) had severe pneumonia. There were no differences 181 in demographic and previous history of respiratory disease including COVID19 infection between those 182 with and without lower respiratory illnesses in these characteristics (Table 1, Figure 1, Figure 2). 183 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 23, 2023.  The mean (95%CI) age of all patients was 2.7 (2.4, 3) years and those with lower respiratory illnesses were 218 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 219 with severe pneumonia were the youngest (mean (95%CI) = 1.8 (0, 3.8)). Males accounted for 55% of all 220 patients and there were no gender bias in those with and without lower respiratory illnesses whereas 4/5 221 (80%) of those with severe pneumonia were males. About 43% of patients were born by Caesaren (C)-222 section and the percentage was similar between those with and without lower respiratory illnesses while 223 4/5 (80%) of those with severe pneumonia were born by C_section (p=0.051). There was no difference in 224 BMI, gestational age, birth weight and duration of breastfeeding between those with and without lower 225 respiratory illnesses or between those with and without severe pneumonia (Table 1). 226 Day of illness at hospital visit of those with lower respiratory illnesses was later than those without (mean 227 (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 228 cough (84%) at disease onset. Most of the patients (79%) had high fever (>39 C) and long fever duration 229 (mean (95%CI) = 6.4 (5.7, 7.1) days). Poor eating (as compared to usual) as reported by caregivers was 230 found more in patients with lower respiratory illnesses than those without. There was no significant 231 difference in symptoms at disease onset such as cough, red eyes, diarrhea, fever maximal temperature, 232 fever duration, ear pain, etc as well as in hematology or biochemistry indexes between those with and 233 without lower respiratory illnesses or between those with and without severe pneumonia. Almost all 234 patients with lower respiratory illnesses (95%) had detectable abnormalities in chest X-rays as well as 16% 235 of those without lower respiratory symptoms (Table 1, Figure 1 and Figure 2). 236 About 25% of the patients with lower respiratory illnesses received oxygen supplements, about 10% 237 received IVIG, 44% received nebulizers, 24% received general steroids while none of the patients without 238 lower respiratory illnesses received these interventions. Duration of hospitalization was longer for those 239 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 23, 2023. ; https://doi.org/10.1101/2023.04.18.23288722 doi: medRxiv preprint with lower respiratory illnesses than those without (mean (95%CI) = 6.1 (4.7, 7.4) vs. 4.5 (3.8, 5.1) days) 240 and longest in those with severe pneumonia (mean (95%CI) = 9.5 (4.9, 14.1) days). All patients except one 241 were discharged without complication (Table 1, Figure 1 and 2). 242 The five patients with severe pneumonia were of a younger age and were predominantly male (80%), of 243 which 80% were born by C-section and were all co-infected with at least one other respiratory viruses or 244 bacteria (p<0.05 as compared to other patients). There was no other remarkable difference in other 245 demographic, medical history or host characteristics as compared to other patients. 246

Co-infection of HAdV and other respiratory viruses or bacteria 247
Regarding viral co-infection, there were 19 (approx. 20%) samples tested positive for at least one other 248 virus from our respiratory virus PCR panel of which 5(5.1%) were co-infected with at least 2 viruses. The 249 number of samples co-infected with HEV, MPV, PIV1, PIV2, PIV3, PIV4 in the PCR panel were 7, 9, 4, 2, 0, 250 2 respectively. There were two samples co-infected with Influenza virus detected by rapid antigen test. In 251 overall, there were 20 samples co-infected with at least 1 virus detected by either PCR or rapid antigen 252 test. Regarding bacterial co-infection, there were 3, 11, 5 samples with positive culture for Streptococcus 253 pneumoniae, Haemophilus influenzae, Moraxella catarrhalis respectively. There were 2 patients with 254 blood antibody test positive for Mycoplasma pneumoniae. In overall, there were 20 patients co-infected 255 with at least one bacterium detected by either nasopharyngeal fluid culture or blood antibody test. There 256 were 5 patients co-infected with both at least one virus and at least one bacterium. In total, there were 257 35 patients co-infected with either virus or bacteria (Table 2). 258

264
The percentage of samples co-infected with at least one other respiratory virus was significantly higher in 265 patients with lower respiratory illnesses vs. those without (12 (50%) vs. 4 (9%), p< 0.001) and especially 266 higher in those with severe pneumonia (4/5 (80%), p=0.006). Co-infection with at least one respiratory 267 bacterium was detected in 20 (20.6%) patients, with 7 of these patients (29.2%) experiencing lower 268 respiratory illnesses (p=0.253) and in 3 (60.0%) with severe pneumonia (p=0.058). Any co-infection with 269 respiratory viruses or bacteria was found in 35 (36.1%) of all patients, significantly higher in those with 270 lower respiratory illnesses (17 (70.8%), p< 0.001) and especially higher in those with severe pneumonia (5 271 (100.0%), p=0.005) ( Table 1, Table 2). 272 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Risk factors for lower respiratory illnesses and severe pneumonia 273
Potential risk factors for lower respiratory illnesses or severe pneumonia which were clinically relevant 274 and statistically significant from the above exploration were viral co-infection of HAdV with ≥1 respiratory 275 viruses or bacterial co-infection of HAdV with ≥1 respiratory bacteria or any co-infection of HAdV with ≥1 276 respiratory viruses or bacteria, age at hospital visit, sex, birth delivery method, day of disease at hospital 277 visit. The results from univariate logistic regression and multivariate logistic regression adjusting for these 278 factors in the models for the comparison between those with lower respiratory illnesses or severe 279 pneumonia vs. those without are summarized in Table 3. Co-infection of HAdV with ≥1 other respiratory 280 viruses was remarkably significantly associated with either lower respiratory illnesses or severe 281 pneumonia (unadjusted OR (95%CI); (p-value) for those with lower respiratory illnesses or severe 282 pneumonia vs. those without were 11.14 (3.77, 36.05); <0.0001 and 19.00 (2.60, 385.24); 0.0105 283 respectively). After adjusting for age at hospital visit, sex, birth delivery method, day of disease at hospital 284 visit, the association between co-infection and lower respiratory illnesses or severe pneumonia remained 285 statistically significant (adjusted OR (95%CI); (p-value) for those with lower respiratory illnesses or severe 286 pneumonia vs. those without were 7.4 (1.4, 52.1); 0.03 and 16.9 (1.8, 421.9); 0.03 respectively). Co-287 infection of HAdV with any respiratory viruses or bacteria was also remarkably significantly associated 288 with lower respiratory illnesses (unadjusted OR (95%CI); (p-value) = 7.42 (2.75, 21.98); 0.0001 and 289 adjusted OR (95%CI); (p-value) = 5.21 (1.60, 19.36); 0.0084). 290 is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 23, 2023.  is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

317
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 23, 2023. 7, 2 in this study except a significant difference in hematology test mentioned above (Figure 4). Seoul Korea and this 11-year study also suggested a genome type shift between successive outbreaks 345 (11). Therefore, further studies examining the genome types or variants of HAdV may be useful. 346 There was no significant association between the 3 HAdV types identified and clinical features in this 347 study, except that WBC and neutrophil numbers were higher in those infected with HAdV type 3 than with 348 type 7. The proportion of severe cases seemed to have similar distribution between the more 349 predominant HAdV type 3 and type 7. This is consistent with other previously published studies in which 350 HAdV type 3 and 7 were the most common types among more severe respiratory HAdV infections (10-351 12,15-19). 352 Co-infections of HAdV with other respiratory viruses or bacteria stood out to be highly significantly 353 associated with increased severity of respiratory diseases. About 70.8% of those with lower respiratory 354 illnesses and 100% of those with severe pneumonia were co-infected with at least one respiratory viruses 355 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 23, 2023. ; https://doi.org/10.1101/2023.04.18.23288722 doi: medRxiv preprint or bacteria while the percentage was about 24.7% in those without lower respiratory illnesses. This 356 association remained statistically significant in a multivariate analysis adjusting for other potential risk 357 factors which were either clinically relevant or also statistically associated with lower respiratory illnesses 358 or severe pneumonia. This indicates that the association is both clinically meaningful and statistically 359 robust. This remarkable association suggests that HAdV might be the cause of the outbreak but might not 360 be the main cause of severe cases. Instead, the co-infection of other respiratory viruses or bacteria 361 together with HAdV might increase the severity of the disease. This advocates the advantage of multi-362 factor respiratory PCR panel(s) for common respiratory viruses and probably also respiratory bacteria 363 together with other microbial tests such as culture for the diagnosis and prognosis of children with 364 respiratory infections. These results may be helpful for the diagnosis, treatment and prognosis of patients 365 in future outbreaks of respiratory HAdV infection. 366 Although the sample size of this study is limited, it is similar to those of many other previously published 367 studies analysing molecular typing for HAdV outbreaks (5,9,16,49-51). The study samples were from 368 pediatric patients of a private general hospital in Hanoi, Vietnam and molecular typing was performed 369 using samples with Ct for HAdV <30. As such, the study samples may not fully represent the general child 370 population of the outbreak. Additionally, it should be noted that clinical data were collected 371 retrospectively for samples with HadV-positive PCR and thus some clinical data features maybe 372 incomplete, especially for outpatients. Moreover, molecular typing data may not be sufficient to 373 understand the characteristics of HAdV causing the outbreak and thus further study to examine HAdV 374 genome types or variants may be useful. 375 In summary, this study revealed that HAdV type 3 and type 7 were predominant in the current ongoing 376 outbreak of respiratory HAdV infection in children in northern Vietnam and HAdV type might not be 377 associated with the severity of the diseases. Instead, co-infection of HAdV together with other respiratory 378 viruses or bacteria appeared to be a significant risk factor for lower respiratory tract illnesses and severe 379 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 23, 2023. ; https://doi.org/10.1101/2023.04.18.23288722 doi: medRxiv preprint pneumonia. These findings advocate the advantages of multi-factor microbial panels for the diagnosis and 380 prognosis of respiratory infections in children. These results may be helpful for the diagnosis, treatment, 381 and prognosis of patients in future outbreaks of respiratory HAdV infection. 382 383 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 23, 2023. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 23, 2023. ; https://doi.org/10.1101/2023.04.18.23288722 doi: medRxiv preprint