key: cord-0713031-blizrzt9
authors: Fukutani, Kiyoshi F.; Nascimento-Carvalho, Cristiana M.; Van der Gucht, Winke; Wollants, Elke; Khouri, Ricardo; Dierckx, Tim; Van Ranst, Marc; Houspie, Lieselot; Bouzas, Maiara L.; Oliveira, Juliana R.; Barral, Aldina; Van Weyenbergh, Johan; de Oliveira, Camila I.
title: Pathogen transcriptional profile in nasopharyngeal aspirates of children with acute respiratory tract infection
date: 2015-06-10
journal: J Clin Virol
DOI: 10.1016/j.jcv.2015.06.005
sha: ae161c98329556419d76f0a33175e1ceae998a20
doc_id: 713031
cord_uid: blizrzt9

BACKGROUND: Acute respiratory tract infections (ARI) present a significant morbidity and pose a global health burden. Patients are frequently treated with antibiotics although ARI are most commonly caused by virus, strengthening the need for improved diagnostic methods. OBJECTIVES: Detect viral and bacterial RNA in nasopharyngeal aspirates (NPA) from children aged 6–23 months with ARI using nCounter. STUDY DESIGN: A custom-designed nCounter probeset containing viral and bacterial targets was tested in NPA of ARI patients. RESULTS: Initially, spiked control viral RNAs were detectable in ≥6.25 ng input RNA, indicating absence of inhibitors in NPA. nCounter applied to a larger NPA sample (n = 61) enabled the multiplex detection of different pathogens: RNA viruses Parainfluenza virus (PIV 1–3) and RSV A-B in 21%, Human metapneumovirus (hMPV) in 5%, Bocavirus (BoV), CoV, Influenza virus (IV) A in 3% and, Rhinovirus (RV) in 2% of samples, respectively. RSV A-B was confirmed by Real Time PCR (86.2–96.9% agreement). DNA virus (AV) was detected at RNA level, reflecting viral replication, in 10% of samples. Bacterial transcripts from Staphylococcus aureus, Haemophilus influenzae, Streptococcus pneumoniae, Moraxella catarrhalis, Mycoplasma pneumoniae and Chlamydophila pneumoniae were detected in 77, 69, 26, 8, 3 and 2% of samples, respectively. CONCLUSION: nCounter is robust and sensitive for the simultaneous detection of viral (both RNA and DNA) and bacterial transcripts in NPA with low RNA input (<10 ng). This medium-throughput technique will increase our understanding of ARI pathogenesis and may provide an evidence-based approach for the targeted and rational use of antibiotics in pediatric ARI.

Acute respiratory tract infection (ARI) accounts for both high mortality and morbidity in children <5 years old [1] . In developing countries, 75% of all acute morbidities were assigned as ARI but identification of the etiological agent in these cases is not common [2] . Upper respiratory tract infections are most commonly caused by virus and only 10% of cases are attributed to bacteria [3] ; nonetheless, most patients are treated with antibiotics [1] .

Respiratory viruses can be diagnosed by culture, antibodybased and nucleic acid-based detection methods. Polymerase chain reaction (PCR) is the most widespread technique for pathogen identification [4] ; however, simultaneous detection of different pathogens by multiplex PCR demands reaction standardization. To this end, nCounter (NanoString Technologies), a digital mRNA expression profiling method, can overcome these drawbacks: it identifies genetic material using two probes, one anchored to the genetic material and the other anchored to a singular bar code that identifies the target by specific hybridization. Upon bar code http://dx.doi.org/10.1016/j.jcv.2015.06.005 1386-6532/© 2015 Elsevier B.V. All rights reserved. hybridization, presence and frequency of the target is revealed and multiple barcodes in the same reaction allow multiplex analysis [5] .

We designed a codeset for quantification of viral (RV, RSV A and B, IV A and B, PIV (types 1 and 3), AV (types 2 and 5), hMPV, BoV, CoV) and bacterial (Streptococcus pneumoniae, Haemophilus influenzae, Moraxella catarrhalis, Staphylococcus aureus, Mycoplasma pneumoniae and Chlamydia pneumoniae) transcripts in one multiplex reaction employing nasopharingeal aspirates (NPA) from children presenting ARI.

Our study complied with the guidelines of the Declaration of Helsinki and was approved by the institutional review board of the Universidade Federal da Bahia (UFBA) (#067/2009). Written informed consent was obtained from all legal guardians. This was a prospective cohort study involving a sample of children aged 6-23 months, seen at UFBA's Emergency Unit between 2009 and 2011 [6] . Community-dwelling patients fulfilling inclusion criteria (fever, cough, sneeze or nasal obstruction for a period up to seven days and without previous episode of wheezing) were enrolled and NPA were collected. Control samples (nasopharyngeal wash) were obtained from individuals (n = 7) undergoing routine outpatient visits. After collection, samples were place in Nuclisens lysis buffer (Biomerieux) and frozen at −70 • C.

Viral and bacterial detection in NPA was performed using nCounter transcriptomic analysis (NanoString Technologies), allowing us to simultaneously assay for viral and bacterial targets. Total RNA (10-50 ng), obtained from a subset of randomly selected NPAs, was extracted using Qiagenís RNEasy, following manufacturerís instructions, and was subsequently hybridized against probes targeting RV, RSVA-B, IV A-B, PIV (types 1 and 3), AV (types 2 and 5), hMPV, BoV and CoV (229A and OC43), S. pneumoniae, H. influenzae, M. catarrhalis, S. aureus, C. pneumoniae and M. Table 1 ), synthesized by NanoString Technologies. Hybridization was also performed with control viral and bacterial RNA pools [kindly donated by K. Bruyninckx, UZ Leuven], in three different concentrations (6.25, 12.5 and 25 ng for viral RNA and 1.5, 3 and 6 ng for bacterial RNA).

RSV A and B were independently quantified by Real-Time PCR, as described [7] .

Raw data were pre-processed using both nSolver 2.0 software (Nanostring Technologies) and the NanoStringNorm R package [8] . Preprocessing sequentially corrects for three factors: technical variation, background and sample content. First, using a set of exogenous positive control RNAs present in each sample (see also Fig. 1 ), technical variation is accounted for by adjusting the counts in each sample to the geometric mean of counts for positive control counts in all samples. Subsequently, a background correction is performed by subtracting the maximum count value of the negative control probes in a sample from each probe output within the sample. Transcripts counts that are negative after background correction are set to 1. Each sample is then normalized for RNA content by adjusting the counts to the geometric mean of 15 housekeeping genes (Supplementary Table 2 ). Finally, the data were log 2 transformed. Full details of each step can be found in the supplementary information of Waggott et al. [8] . Network data were designed by association matrix and the graph file was created using Gephi. Nodes were aligned by Force Atlas algorithm and the centrality was calculated by the betweenness. Statistical analyses were performed using Graphpad Prism version 5.0.

The digital nCounter profiling technology captures and counts specific nucleic acid molecules in a complex mixture [5] . Briefly, unique pairs of capture and detection probes are designed against each gene of interest. The capture probe consists of a sequence complimentary to a target pathogen mRNA and a short common sequence coupled to an affinity tag (biotin). The reporter probe also contains the complementary sequence but is coupled to a colorcoded tag, providing the detection signal. The color code in each tag is unique, allowing the accurate discrimination of the different targets, in a complex mixture. Capture and reporter probes are mixed with total RNA in a single hybridization reaction, resulting in the formation of tripartite structures (target bound to the specific reporter and capture probes). Unhybridized probes are washed away and remaining complexes are exposed to streptavidin. Tar-gets bound to the specific probes are immobilized on a solid surface that is scanned; the unique reporter probe then identifies each target and the level of expression is measured by the digital readout of transcript counts. The reaction occurs in the absence of nucleic acid amplification and/or enzymatic reactions. To estimate target concentration, nCounter uses a series of positive controls covering a range of known concentrations (0.1-100 fM, comparable to 0.2-200 copies per cell). Based on the counts obtained for these known concentrations, the concentration of mRNA transcripts can be estimated. This estimate is done by plotting the counts obtained for the positive controls versus the concentrations tested and by fitting a regression line through the average (Fig. 1) . The regression line allows the quantification of mRNA transcript concentration. The correlation between the counts for positive controls and the concentration was 0.99 and 0.98 for the lower concentration range.

Next, we assessed the potential of nCounter as a diagnostic technique for ARI. In pilot experiments, we used RNA obtained from NPA of four patients and the set of probes custom-designed Table 1 ). We detected the presence of AV, CoV and RSVB ( Fig. 2A) , whereas BoV, PIV, RSVA, IV A-B, hMPV and RV were detected in none of the four samples initially tested. Detection was also positive in NPA spiked with known viral RNA (6.25-12.5-25 ng), indicating absence of inhibitors in RNA isolated from NPA (Fig. 2B ). In addition, use of tRNA (an RNA carrier) to ARI samples did not improve viral detection ( Supplementary Fig.  1 ), confirming that even NPA containing low input RNA (<10 ng) is suitable. Bacterial RNA (H. influenzae, M. catarrhalis, M. pneumoniae, S. aureus, and S. pneumoniae) was also detected in both patientsí and control NPA (Fig. 2C) .

For comparison purposes, we performed RSV A and B detection by Real-Time PCR. Agreement between nCounter and Real-Time PCR was 86.2% for RSVA (kappa = 0.231) and 96.9% for RSVB (kappa = 0.82) ( Table 1 ). RNA counts detected by nCounter were significantly correlated to the number of copies detected by Real-Time PCR for both RSVA (r = 0.25, Spearman; p = 0.049) and RSVB (r = 0.83; Spearman p = 0.001). 

We then expanded nCounter analysis to a set of 61 randomly selected NPA and included seven healthy control samples (Fig. 3) . Using this extended group, we detected RV (2% of samples), BoV (3%), CoV (229E and OC43, 3%), IV A (3%), hMPV (5%), PIV (genotypes 1 and 3, 21%) and RSV A and B (21%). In addition, detection of DNA virus AV (types 2 and 5) was confirmed at the RNA level in 10% of samples, indicating ongoing (intracellular) viral replication. Transcripts for C. pneumoniae, M. catarrhalis, M. pneumoniae and S. pneumoniae were detected in 2, 3, 8 and 26% of samples, respectively. H. influenzae and S. aureus transcripts were ubiquitous (69 and 77% of samples, respectively) and were also detected in healthy controls (28 and 100% of samples, respectively). Using this data set, we observed bacterial presence in the majority of children (55.1%), whereas in 43.5% of children, both viruses and bacteria were present (Fig. 4A) . A low percentage of children (1.5%) presented viruses only. Upon observation of the number of transcripts per type of pathogen (Fig. 4B) , AV 2 and RSB (A and B) showed Network analysis also revealed the complex interactions that occur between individual viruses and bacteria detected in NPA (Fig. 5) . For example, the strongest links were observed for bacteria: H. influenzae/S. aureus and, secondly, H. influenzae/S. aureus or H. influenzae/S. pneumoniae. Comparing virus and bacteria, PIV1/H. influenzae showed a strong interaction as did RSV (A and B)/S. aureus. Of note, the interactions of CoV 229E were different from those observed for CoV OC43.

Our results demonstrate that nCounter is a robust mediumthroughput technique for detection of both viral and bacterial RNA transcripts in ARI patients. The sensitivity and reproducibility of nCounter has already been compared with Openarray ® [9], Affymetrix microarray [10] and TaqMan PCR [5] . In the latter, nCounter was equally reproducible and sensitive regarding detection of differentially modulated genes; however, nCounter detected low abundant transcripts, not detected by microarray. nCounter was also successfully used for detection of pathogen RNA transcripts in clinical samples [11] , and an important advantage is transcript detection without the need for purification or amplification, yielding a rapid indicator of pathogen replication.

Transcriptomic profiling by nCounter allowed simultaneous detection of both viral and bacterial transcripts in ARI patients. We detected the presence of AV, BoV, CoV, IV, hMPV, PIV, RV and RSVA-B, which are associated with respiratory tract infections [1] . To our knowledge, this is the first description of detection of a pathogenic DNA virus (Adenovirus type 2 and 5) at the RNA level in clinical samples. In contrast to a standard DNA PCR, which quantifies DNA in (mostly extracellular) viral particles, nCounter enables demonstration of ongoing (intracellular) viral replication, which we hypothesize to be more relevant to ARI pathogenesis. We also, simultaneously detected transcripts for C. pneumoniae, M. catarrhalis and S. pneumoniae, which are commonly isolated from the respiratory tract [12] , as is H. influenzae. We observed a good agreement between Real-Time PCR and nCounter for RSV (A and B) detection, reinforcing the possibility of using the latter as an alternative method for pathogen detection. Of note, RSVA strains present a high variability in the target used for detection (G protein domain) [13] , also in Brazil [6] , which could explain the weaker correlation observed for nCounter and Real-Time PCR, as compared to RSVB.

The healthy upper respiratory tract is colonized by commensals and potential pathogens kept in check by the immune system. There is evidence that viral respiratory infections can stimulate an increase in bacterial load and this association can lead to secondary complications [14] . Preterm infants with RSV and bacterial co-infection were hospitalized for longer periods and more frequently admitted to the intensive care unit, compared to children infected with RSV alone [15] . In our samples, the majority of ARI patients presented H. influenzae and S. aureus, which, together with S. pneumoniae and M. catarrhalis are the usual opportunistic pathogens co-detected with RSV [16] [17] [18] . Indeed, transcripts for all these agents were detected in our NPA samples from children with ARI. Network analysis showed strong links between PIV 1 and H. influenzae and between RSV B and S. pneumoniae, for example. Both S. pneumoniae [19] and RSV activate the inflammasome [20] ; S. pneumoniae triggers NF-kB activation leading to an inflammatory response whereas, RSV infection results in a strong interferon response. The possible overlap between the signaling pathways triggered by each pathogen may synergize and amplify the inflammatory signal during co-infection. This exacerbated immune response could accentuate damage in situ and increase morbidity. Given this scenario, knowing the effective pathogen load, as shown here by nCounter, can contribute towards understanding the contribution of each pathogen to disease severity. As a proof-of-principle, our group has already used nCounter to simultaneously quantify both human and retroviral transcriptomes; thus, revealing intact IFN signaling in HTLV-1 infection [21] .

We conclude that nCounter can be used as an alternative medium-throughput method to identify, in a single reaction, the presence of multiple pathogens in NPA from ARI patients. Given that respiratory infections are one of the most common reasons for visits to physicians, improved diagnostic methods such as nCounter can contribute to better patient management. Moreover, a transcriptomic approach involving pathogen detection, and as discussed above, immune response transcripts will enable the identification of biomarkers in ARI clinical samples, even those presenting low RNA yield (10-50 ng) such as NPA. This approach, combined with long-term clinical follow-up, will increase our understanding of ARI pathogenesis, disease progression (wheez- ing, asthma) and, more importantly, provide an evidence-based approach for the targeted and rational use of antibiotics in pediatric ARI.

This work was supported by Fundaç ão de Amparo à Pesquisa do Estado da Bahia (FAPESB; grant PNX0019/2009). KFF was supported by fellowships from Science without Borders (Conselho Nacional de Desenvolvimento Científico e Tecnológico; CNPQ) and FAPESB. RK was supported by a fellowship from CNPq. AB, CN-C and CIO are senior investigators from CNPq.

Our study complied with the guidelines of the Declaration of Helsinki and was approved by the institutional review board of the Universidade Federal da Bahia (UFBA) (#067/2009). Written informed consent was obtained from all legal guardians.

Viral infections of the lower respiratory tract: old viruses, new viruses, and the role of diagnosis

Detection of respiratory viruses by molecular methods

Systematic review of the treatment of upper respiratory tract infection

PCR-based diagnostics for infectious diseases: uses, limitations, and future applications in acute-care settings

Direct multiplexed measurement of gene expression with color-coded probe pairs

Validating report of first episode of wheezing with pediatrician-detected wheezing among children

Circulation of HRSV in Belgium: from multiple genotype circulation to prolonged circulation of predominant genotypes

NanoStringNormn: an extensible R package for the pre-processing of NanoString mRNA and miRNA data

Systematic evaluation of medium-throughput mRNA abundance platforms

High throughput digital quantification of mRNA abundance in primary human acute myeloid leukemia samples

RNA signatures allow rapid identification of pathogens and antibiotic susceptibilities

Colonization and infection of the respiratory tract: what do we know

Epidemiological aspects of influenza A related to climatic conditions during and after a pandemic period in the city of Salvador, Northeastern Brazil

Viral-bacterial interactions and risk of acute otitis media complicating upper respiratory tract infection

Risk of concurrent bacterial infection in preterm infants hospitalized due to respiratory syncytial virus infection

Incidence of bacterial coinfection with respiratory syncytial virus bronchopulmonary infection in pediatric inpatients

Pulmonary and systemic bacterial co-infections in severe RSV bronchiolitis

High incidence of pulmonary bacterial co-infection in children with severe respiratory syncytial virus (RSV) bronchiolitis

TLR2/MyD88/NF-kappaB pathway, reactive oxygen species, potassium efflux activates NLRP3/ASC inflammasome during respiratory syncytial virus infection

Critical roles of ASC inflammasomes in caspase-1 activation and host innate resistance to Streptococcus pneumoniae infection

Simultaneous RNA quantification of human and retroviral genomes reveals intact interferon signaling in HTLV-1-infected CD4+ T cell lines

Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.jcv.2015.06.005