Article Text

Original article
Respiratory viruses associated with community-acquired pneumonia in children: matched case–control study
Free
  1. Samuel Rhedin1,
  2. Ann Lindstrand2,3,
  3. Annie Hjelmgren1,
  4. Malin Ryd-Rinder4,
  5. Lars Öhrmalm1,
  6. Thomas Tolfvenstam1,2,
  7. Åke Örtqvist1,5,
  8. Maria Rotzén-Östlund6,7,
  9. Benita Zweygberg-Wirgart6,7,
  10. Birgitta Henriques-Normark6,7,
  11. Kristina Broliden1,
  12. Pontus Naucler1
  1. 1Department of Medicine Solna, Infectious Diseases Unit, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden
  2. 2The Public Health Agency of Sweden, Stockholm, Sweden
  3. 3Department of Public Health Sciences, Karolinska Institutet, Stockholm, Sweden
  4. 4Department of Clinical Science and Education, Karolinska Institutet, Sachs’ Children and Youth Hospital, South General Hospital, Stockholm, Sweden
  5. 5Department of Communicable Disease and Control, Stockholm County Council, Stockholm, Sweden
  6. 6Department of Clinical Microbiology, Karolinska University Hospital, Stockholm, Sweden
  7. 7Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
  1. Correspondence to Dr Samuel Rhedin, Infectious Diseases Unit, Department of Medicine Solna, Karolinska Institutet and Karolinska University Hospital, Stockholm 171 76, Sweden; samuel.rhedin{at}ki.se

Abstract

Background Community-acquired pneumonia (CAP) is the leading cause of death in children worldwide and a substantial proportion of childhood CAP is caused by viruses. A better understanding of the role of virus infections in this condition is needed to improve clinical management and preventive measures. The aim of the study was therefore to assess the association between specific respiratory viruses and childhood CAP.

Methods A case–control study was conducted during 3 years in Stockholm, Sweden. Cases were children aged ≤5 years with radiological CAP. Healthy controls were consecutively enrolled at child health units during routine visits and matched to cases on age and calendar time. Nasopharyngeal aspirates were obtained and analysed by real-time PCR for 15 viruses. Multivariate conditional logistic regression was used to account for coinfections with other viruses and baseline characteristics.

Results A total of 121 cases, of which 93 cases met the WHO criteria for radiological pneumonia, and 240 controls were included in the study. Viruses were detected in 81% of the cases (n=98) and 56% of the controls (n=134). Influenza virus, metapneumovirus and respiratory syncytial virus were detected in 60% of cases and were significantly associated with CAP with ORs >10. There was no association with parainfluenza virus, human enterovirus or rhinovirus and coronavirus and bocavirus were negatively associated with CAP.

Conclusions Our study indicates viral CAP is an underestimated disease and points out hMPV as a new important target for the prevention of childhood CAP.

  • Pneumonia
  • Viral infection
  • Clinical Epidemiology
  • Paediatric Physician

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Key messages

What is the key question?

  • What is the role of virus infections in community-acquired pneumonia in children?

What is the bottom line?

  • Infection with respiratory syncytial virus, metapneumovirus and influenza virus infections were highly associated with community-acquired pneumonia in children.

Why read on?

  • We performed a carefully designed case–control study that points out metapneumovirus as a new important target for the prevention of childhood pneumonia.

Introduction

Community-acquired pneumonia (CAP) is the leading cause of death in children aged ≤5 years worldwide.1 A substantial proportion of childhood CAP is estimated to be caused by viruses, yet viral CAP is a condition that only recently is gaining increased attention.2 ,3 Successful vaccination programmes against Haemophilus influenzae and Streptococcus pneumoniae have decreased the mortality and incidence of bacterial CAP, and it is likely that the relative importance of viral CAP will increase.4 A better understanding of the contribution of specific respiratory viruses in childhood CAP is needed to improve clinical management and preventive measures.

Human adenovirus (HAdV), human bocavirus (HBoV), human coronavirus (HCoV) enterovirus (EV), influenza virus (Flu), metapneumovirus (hMPV), parainfluenza virus (PIV) and respiratory syncytial virus (RSV) infections have all been associated with childhood CAP.2 However, the majority of studies have been descriptive and case series of hospitalised patients and hence only provide a low-grade level of evidence.2 Few case–control studies have investigated the association between viruses and childhood CAP and most of these studies have been done in low-income countries that might not represent the epidemiological pattern of high-income countries, or they used hospital-based controls that might bias the results.5–9 Moreover, childhood CAP is difficult to diagnose. The Integrated Management of Childhood Illness (IMCI) diagnostic guidelines for the clinical diagnosis of CAP, developed by the WHO, have been important to reduce mortality, yet they lack sensitivity and specificity and might introduce misclassification bias in epidemiological studies.10 The gold standard today is diagnosis based on radiographic findings, which is the preferred method of case definition in epidemiological studies.11

Advances in molecular diagnostic methods have increased the sensitivity of viral detection and led to discovery of novel viruses that cannot be detected by traditional virus isolation.12 Real-time PCR is a sensitive method that has been widely introduced in the routine clinical management in high-income countries. However, the clinical significance of a real-time PCR detection is obscured by extensive asymptomatic detection in children.13 This underscores the importance of representative controls as a reference when assessing the impact of respiratory viruses on CAP.

The primary aim of the study was to assess the association between specific respiratory viruses and childhood CAP. A secondary aim was to investigate the clinical presentation of different respiratory virus infections in children.

Methods

Study design

A matched case–control study was conducted in Stockholm, Sweden, between September 2011 and March 2014, thus covering three winter seasons. No sampling was performed during summers. Notably, vaccination with pneumococcal conjugate vaccines (PCVs) in the child immunisation programme had been in place in the study setting since October 2007 (Prevenar 7 was replaced by Prevenar 13 in January 2010). Cases were children aged ≤5 years with radiographic evidence of CAP (defined by a primary radiologist) enrolled at the paediatric emergency room or inpatient wards (within 48 h following admission) at Sachs’ Children and Youth Hospital, South General Hospital and Astrid Lindgren Children's Hospital, Huddinge in Stockholm, Sweden. To reduce misclassification bias, all chest radiographs were retrospectively reviewed according to the WHO guidelines (endpoint alveolar consolidation or pleural effusion) by an independent senior paediatric radiologist.11

Controls were randomly selected from an ongoing study on carriage of S. pneumoniae in children. In this study, healthy children are consecutively sampled at child health units in Stockholm County during routine visits for regular check-ups or vaccination within the national childhood immunisation programme (covering >99% of the population).14 Reported respiratory symptoms within the last 7 days were recorded, but were not an exclusion criteria. Up to two controls were matched to each case on calendar time (±14 days) and age (±6 months). If eligible controls were not found, intervals were expanded to ±30 days and ±12 months.

A standardised questionnaire was used to collect information regarding baseline characteristics in both study groups and cases also underwent a clinical examination at the time of enrolment. Tachycardia was defined as a heart rate of ≥160 beats/min (bpm) in children aged <1 years and ≥120 bpm in children aged ≥1 years. Hypoxaemia was defined as peripheral oxygen saturation <90% by pulse oximetry. Tachypnoea was defined as a respiratory rate of ≥60 breaths/min in children aged <2 months, ≥50 breaths/min in children aged 2–12 months and as ≥40 breaths/min in children aged >1 year. Fever was defined as a body temperature ≥38.0°C. Data on management and outcome of the cases were collected from the medical records. Severity of disease was defined according to the IMCI guidelines developed by the WHO15 and the respiratory index of severity in children (RISC) score was calculated based on Swedish growth standards.16 ,17 Informed consents were collected from the guardians before enrolment.

Sampling and microbiological analyses

Nasopharyngeal aspirates were collected with an identical technique in both study groups.13 The samples were stored at 6°C–8°C during the day of sampling and subsequently frozen at −80°C. After thawing, a standardised real-time PCR assay was performed for the following 15 viruses at the accredited (ISO 15189:2012) Microbiological Laboratory at Karolinska University Hospital: HAdV, HBoV, HCoV (species HKU1, NL63, OC43 and 229E), EV, Flu (A, A(H1N1)pdm09 and B), hMPV, PIV 1–3, rhinovirus (RV), RSV (A and B).18 Single agent real-time PCRs were used for most viruses, nevertheless duplex real-time PCRs were used for PIV 1 and 3, for PIV 2 and HCoV 229E and for RSV A and B. Due to cross-reactivity between RV and EV, samples positive for both viruses were analysed with an in-house real-time PCR for EV.19 Cycle threshold (CT) values were recorded as an inverse estimate of viral load, yet all samples were also dichotomously coded as positive or negative.18

Data analysis

Data were analysed in STATA V.12.0. Mann–Whitney U test was used for continuous data and Fisher's exact test or χ2 test was used for categorical data as appropriate. Conditional logistic regression was used to compare differences in baseline characteristics and the presence of viruses between cases and controls. To account for confounding by coinfections and differences in baseline characteristics between cases and controls, a multivariate conditional logistic regression model was constructed that included variables that had p values <0.1 in initial univariate analyses. A backward elimination approach was then performed based on likelihood ratio tests. Data were presented as ORs and 95% CI. Due to high frequency of coinfections the initial matching strategy generated an unexpectedly large number of strata when adjusting for coinfections and resulted in large drift of the ORs (see online supplementary table S1). To account for that, individual matching sets were broken and instead study subjects were group-matched according to season (spring, fall, winter), year (2011, 2012, 2013) and age (<1 and ≥1 year), generating 18 strata. Separate analyses were performed for cases with pneumonia verified by a primary radiologist and for children with chest radiographs meeting the criteria by the WHO for epidemiological studies. The results of the first analysis represent the clinical reality of paediatricians and the second is according to gold standard.

Results

Baseline characteristics

A total of 123 cases were enrolled during the study period. Two cases were excluded due to incomplete sampling. Data were not collected on non-participating study subjects. The research staff involved in enrolment approximate that 10% of parents to cases and 50% of parents to controls declined to participate in the study, mainly due to fear of discomfort for the child or lack of time. The expanded matching criteria had to be applied in 52 controls and for one case no eligible controls were found, hence 121 cases and 240 matched controls were included in the analyses. Ninety-three cases had presence of alveolar consolidation or pleural effusion on chest radiographs, consistent with the WHO criteria. Underlying diseases were more common among cases than controls (24% compared with 4% of controls, p<0.001) and 18% of the cases suffered from asthma or reactive airway disease (table 1).

Table 1

Baseline characteristics of study subjects

A larger proportion of the cases was on current treatment with antibiotics (13% compared with 1% of controls, p<0.001) or had received antibiotics during the last year (46% and 28%, respectively, p<0.01). Moreover, smoking among parents was more common in cases, whereas controls were more likely to have parents with a university degree and also more likely to have been travelling abroad within the last 3 months (p=0.01, p=0.03 and p=0.02, respectively).

Viruses associated with CAP

One or more viruses were detected in 98 (81%) of the cases and 134 of the controls (56%; p<0.001; table 2). RSV (n=39, 32%), hMPV (n=28, 23%) and RV (n=28, 23%) were the most common viruses among cases, whereas RV (n=64, 27%), HBoV (n=50, 21%) and HCoV (n=28, 12%) were the most commonly detected viruses in the controls. In univariate analyses, HAdV (p=0.05), Flu (p=0.03), hMPV (p<0.001) and RSV (p<0.001) were significantly associated with CAP (table 2 and figure 1). In contrast, negative associations with CAP were seen for HCoV (p=0.03) and HBoV (p=0.04). A multivariate regression model was constructed, accounting for coinfection with multiple viruses and differences in baseline characteristics. RSV (OR=21.3, 95% CI 8.3 to 54.4), hMPV (OR=14.9, 95% CI 5.4 to 41.2) and Flu (OR=13.8, 95% CI 2.5 to 76.6) had the highest ORs for CAP (table 2), combined being present in 60% of the cases. The significant association between HAdV and CAP disappeared in adjusted analyses (OR=1.4, 95% CI 0.5 to 3.9), whereas a negative association with CAP for HCoV (OR=0.2, 95% CI 0.0 to 0.6) and HBoV (OR=0.3, 95% CI 0.1 to 0.8) was still present. In the analyses restricted to cases with radiological CAP consistent with WHO criteria (n=93), the same significant associations were seen as in the analyses of all patients and there were only minor differences in the OR point estimates (table 2). To assess potential effect modification by age, multivariate conditional logistic regression analyses were performed for the most frequently detected viruses with the virus and age (<1 or ≥1 year) included as binary variables together with an interaction term (age category multiplied with the virus). No significant effect modification by age was observed (p=0.68, p=0.55 and p=0.62 for hMPV, RSV and RV, respectively). The results were similar when age was categorised as ≤2 and >2 years.

Table 2

Viral detections in cases and controls

Figure 1

Viral PCR positivity in cases and controls. Relative number of real-time PCR detections for cases (grey, n=121) and controls (white, n=240). Number above bars indicates percentage of positive detections within group. p Values calculated with conditional logistic regression. *p<0.05, **p<0.01 and ***p<0.001.

CT values for each virus were compared between cases and controls as an inverse estimate of viral load (figure 2). Cases positive for hMPV, RV and RSV had significantly lower CT values as compared with controls (p<0.01, p=0.02 and p<0.01, respectively).

Figure 2

PCR cycle threshold (CT) values in cases and controls. Box plot of CT values for all PCR-positive cases (grey) and controls (white) sorted by virus. p Values calculated with Mann–Whitney U test. *p<0.05, **p<0.01, ***p<0.001.

Of the 121 cases, 60 (50%) had a single infection and 38 (31%) had coinfections. In 31 cases (26%) two viruses were detected, five cases (4%) had three detected viruses, one case (1%) had four detected viruses and in one case (1%) five viruses were detected. Coinfections were seen in 50 (21%) of the controls, with two viruses detected in 36 (15%), three viruses detected in 11 (5%), four viruses detected in one and five viruses detected in two (1%) controls. In fact, all viruses except for RSV and Flu were predominantly detected as a coinfection in the cases (see online supplementary table S2). The most common virus pairs among cases were HAdV/RV (n=7), HAdV/hMPV (n=5), HBoV/RSV (n=5) and hMPV/RV (n=5). Flu, hMPV and RSV appeared to follow a seasonal distribution with peaks during winter/spring (see online supplementary figure S1). Flu tended to be more common during the winter season 2012/2013 (p=0.07), whereas no difference in detection rates between the seasons was observed for RSV or hMPV.

Clinical presentation in relation to viral findings

Tachycardia, increased respiratory rate and fever were the most commonly observed clinical signs seen in 85%, 72% and 72% of the cases, respectively (table 3). Wheezing was noted in 41% of the cases and tended to be associated with detection of HBoV (p=0.08) and to be negatively associated with HAdV and HCoV (p=0.06 and p=0.08, respectively). Ninety-two (76%) cases were admitted to inpatient wards for a median of 3 days, three patients were transferred to the paediatric intensive care unit and none of the children died during their hospitalisation. Antibiotic treatment was started on 115 (95%) cases, of these 45 (37%) received parenteral antibiotics. Seventy-seven (64%) cases received inhalation therapy. As part of routine management nasopharyngeal bacterial cultures were collected from 10 cases, all were positive for one or more bacteria: H. influenzae (n=6), S. pneumoniae (n=5), Moraxella catharrhalis (n=5). Mycoplasma pneumoniae was detected by PCR in two of 26 tested cases. Five cases had radiographic evidence of pleural effusion. One out of 29 tested cases had a positive blood culture (positive for S. pneumoniae). This patient also had S. pneumoniae in the nasopharynx, tested positive for RSV and HCoV and had a massive infiltrate in the lower left lobe. He had received two doses of Prevenar 13 and serotyping revealed a non-vaccine type (serotype 38) of S. pneumoniae.

Table 3

Clinical signs and management of cases

C reactive protein (CRP) analyses were ordered in 118 of 121 cases with a median of 73 mg/L. HBoV, EV and hMPV were associated with low CRP (p=0.04, p=0.04 and p=0.02 respectively). The median RISC score was 1 and a RISC score of ≥3 was negatively associated with detection of RV (p=0.02). Severe pneumonia according to the clinical IMCI guidelines was observed in 48 (40%) cases and was negatively associated with detection of RV (p=0.04). Seventeen cases did not meet the IMCI criteria for CAP. Yet, 13 of these were classified as radiological CAP according to the WHO guidelines in the independent radiological revision.

Discussion

The role of respiratory viruses in childhood CAP is not fully understood. To address this issue, we set up a matched case–control study with strict definitions of CAP based upon radiological findings. Controls were selected to represent the source population, sampling and collection of baseline data were standardised and hence we believe that the results have high validity.

Viruses were detected in >80% of the cases. This is in line with Jain et al3 who recently reported detection of viruses in 73% of children hospitalised with CAP in a large prospective study in the USA. RSV, hMPV and Flu were present in 60% of the cases and were highly associated with CAP with adjusted OR >10. Importantly, no positive association was observed between CAP and infections with RV, HBoV, EV nor with any of the HCoV. This needs to be considered when interpreting viral findings with similar real-time PCR detection methods in children with CAP.

Recent publications suggest that hMPV accounts for extensive morbidity in children with an estimated hospitalisation rate of 1 in 1000 in children aged ≤5 years in the USA.20 The virus is a member of the paramyxoviridae family (subfamily pneumovirinae), together with RSV, hence it shares many similarities with RSV both in terms of virus structure and clinical presentation.21 The epidemiology is characterised by a seasonal distribution largely varying in incidence from year to year; however in this study, hMPV was detected at high rates during all three respiratory seasons.22 In the northern hemisphere, the highest incidence is seen during late winter and early spring.22 Infection with hMPV is associated with both upper and lower respiratory tract infections and the virus has been shown to successfully replicate in pneumocytes of macaques, but the role in CAP is not fully understood.22 ,23 Being detected in more than a fifth of all children with CAP with a significant over-representation among cases, our study suggests the virus to play an important role in childhood CAP. Moreover, hMPV was associated with lower CRP levels, which is indicative of viral CAP.24

Flu and RSV were also associated with CAP and only seldom detected in controls with detection rates concordant with Jain et al3 who detected RSV and Flu in 28% and 7%, respectively, in hospitalised children with CAP in the USA. This indicates that a substantial burden of CAP can be attributed to RSV, whereas the role of Flu seems to be less important in this condition.

Viruses were detected in more than half of the control subjects, the most common were RV, HBoV and HCoV. The high detection rate can partly be explained by our decision not to exclude controls with respiratory symptoms reported within the last 7 days. However, exclusion of controls with respiratory tract infection would make the control group less representative of the source population and would likely overestimate the contribution of viruses in CAP.25 Prolonged shedding of viruses such as HBoV that has been reported for up to 10 weeks is likely to reduce the clinical specificity of real-time PCR.26

Interestingly, HBoV and HCoV were negatively associated with CAP. Negative interaction between viruses, such as between Flu and RV, has been reported.27 Many infections are followed by an excessive proinflammatory response including accumulation of immune cells. It can thus be speculated whether the immune responses elicited by the more pathogenic microbes in the symptomatic children were eliminated or protected against infection by these presumably less pathogenic viruses.

HAdV was not significantly associated with CAP in adjusted analyses and was rather frequently detected in controls. HAdV is a pathogen that for long has been known to cause CAP,2 but the virus has a stable configuration and has been reported to persist in lymphatic tissue for weeks.28 The prevalence of PIV and EV was low in both cases and controls. PIV has previously been associated with childhood CAP, in contrary to our results.3 ,5 ,7–9 ,29 The real-time PCRs for PIV are monitored by yearly participation in quality assessment panels. Nevertheless, as with all PCRs there is always a risk of escape by unknown mutations.

High viral load of certain viruses has been associated with severe respiratory disease in children.12 In our data, CT values for hMPV, RSV and RV were significantly lower in cases than in controls indicating a higher viral load. However, the estimation of viral load is associated with uncertainty in nasopharyngeal sampling and we did not have an internal control to account for sample quantity.

A difference in educational level and smoking among parents was observed between cases and controls. To adjust for these differences, parental smoking and educational level were included in the initial multivariate regression models but did not affect the point estimates. To further investigate socioeconomic differences between the groups, we obtained data on purchase power from the child healthcare services in Stockholm County based on postal codes of the study subjects.14 The proportion of families with low purchase power in the whole of Stockholm County was 20.9% (index 1.0). The residential areas of cases had a lower purchase power compared with controls (index 1.16 vs 0.75, p<0.001). Since we did not have individual data on purchase power we did not include it in the multivariate model. Hence, we cannot rule out some residual socioeconomic confounding in the analyses.

To minimise misclassification bias, all chest radiographs were revised by an independent senior radiologist and a subanalysis of cases meeting the WHO criteria for radiological CAP in children was done generating similar results as the full-scale analysis.

It is difficult to obtain valid samples to define whether CAP is caused by bacteria in children; blood cultures are rarely positive and representative respiratory samples are difficult to obtain.30 We did not systematically obtain bacterial sputum samples or blood cultures, hence it is not clear to what extent cases were coinfected with bacteria. We did perform pneumococcal bacterial cultures on all nasopharyngeal samples where 32% of the cases and 33% of the controls were positive (p=0.747), which as previous studies show that pneumococcal isolation in the nasopharynx in children with CAP is of limited use to determine aetiology.31 Several studies both in mice and in humans indicate that certain viral infections predispose for bacterial superinfection.32–34 A recent study of lung aspirates in children with CAP in the Gambia showed S. pneumoniae in 90% of the samples, yet this was prior to introduction of PCV and most detected pneumococci were serotypes included in the PCV.35 They did not find any hMPV which might indicate that the aetiology of CAP in high- and low-income countries differ.10 This is further supported by the low detection rate of hMPV (3%) in the study by Berkley et al5 on Kenyan children with severe CAP. Yet, it is also possible that viral detection in the nasopharynx might not be fully representative of findings in the lung.

Our study indicates viral CAP is an underestimated disease and points out hMPV as a new important target for the prevention of childhood CAP. Better diagnostic methods are warranted to distinguish viral from bacterial CAP, which will clarify if these viral infections cause CAP on their own or merely predispose for bacterial coinfection. This will be pivotal for decision making regarding antibiotic treatment.

Acknowledgments

We would like to acknowledge paediatric resident Dr Sara Leeb, pediatric radiologist Dr Helen Stålemark, medical students Karolin Petersson, Henrik Hurtigh and Hritul Karim, research nurses Kerstin Jämtberg, Carita Krokstrand, Marie Olander-Carlson, Eva Sjögren and Karin Grunditz-Melin, as well as the staff at Sachs’ Children and Youth Hospital, at Karolinska Children's Hospital Huddinge and at the Karolinska Microbiological Laboratory in Stockholm, above all Pia Andersson.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Twitter Follow Samuel Rhedin at @rhedinsamuel

  • Contributors SR coordinated enrolment of case patients at the study sites, collected and interpreted data from patient history files of the case patients, performed statistical analyses and wrote the initial draft of the manuscript; PN designed and conceptualised the study and supervised collection of samples, statistical analyses and writing of the manuscript; AH supervised data collection at one of the study sites, critically re-evaluated medical histories, performed a pilot analysis of the initial data and critically reviewed and revised the manuscript; AL, BH-N and ÅÖ designed the data collection instruments, managed and supervised data collection of the control subjects, interpreted sociodemographic and clinical data of the control subjects and reviewed and revised the manuscript; LÖ, MR-Ö, TT and BZ-W designed the standardised sampling protocol, conducted, supervised and interpreted the microbiological analyses and reviewed and revised the manuscript from a clinical microbiological point of view and MR-R supervised data collection at one of the study sites, critically re-evaluated medical histories and interpreted clinical signs of infection of the case patients and reviewed and revised the manuscript from a paediatric point of view. KB designed and supervised the study, interpreted data analyses and critically reviewed and revised the manuscript; all authors approved the final version of the manuscript.

  • Funding The study was supported by grants from the Stockholm County Council, the Shizu Matsumura Donation, the Swedish Society of Medicine and the Swedish Childhood Cancer Foundation.

  • Competing interests MR-R has been clinical investigator in an epidemiological study on rotavirus with Glaxo Smith Kline and Sanofi Pasteur MSD 2008–2009 and a clinical study on RotaTeq vaccine for Sanofi Pasteur MSD 2013–2014.

  • Patient consent Obtained.

  • Ethics approval Regional Ethical Review Board in Stockholm (ref 2011/879-31/1).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement As to date, we do not hold any unpublished data from the study. Nevertheless, we have plans on performing additional analyses on the samples collected.

Linked Articles