key: cord-322738-nvyknvc7
authors: Delangue, Julie; Roca Sanchez, Yelin; Piorkowski, Géraldine; Bessaud, Maël; Baronti, Cécile; Thirion-Perrier, Laurence; Mafayle, Roxana Loayza; Ardaya, Cinthia Avila; Aguilera, Gabriela Añez; Guzman, Jimmy Revollo; Riera, Javier Lora; de Lamballerie, Xavier
title: Viral aetiology influenza like illnesses in Santa Cruz, Bolivia (2010–2012)
date: 2014-02-24
journal: Virol J
DOI: 10.1186/1743-422x-11-35
sha: 
doc_id: 322738
cord_uid: nvyknvc7

BACKGROUND: Acute respiratory infections represent a serious public health issue worldwide but virological aetiologies of Influenza Like Illnesses (ILIs) remain largely unknown in developing countries. This study represents the first attempt to characterise viral aetiologies of ILIs in Bolivia. METHODS: It was performed in Santa Cruz city from January 2010 to September 2012, based on 564 naso-pharyngeal swabs collected in a National Reference Laboratory and real-time PCR techniques, viral cultures and phylogenetic analyses. RESULTS: 50.2% of samples were positive for at least one virus with influenza viruses (Flu A: ~15%; Flu B: ~9%), rhinoviruses (~8%), coronaviruses (~5%) and hRSV (~4%) being the most frequently identified. The pattern of viral infections varied according to age groups. The elucidation rate was the highest (>60%) amongst patients under 10 yo and the lowest (<40%) amongst patients ≥60 yo. Nearly 3% of samples showed dual viral infections. Epidemiological peaks were associated with a predominant virus but generally included 30-50% of infections by different viruses. Unexpectedly, the frequency of influenza in the 0–4 yo population was very low and a complete hRSV eclipse occurred in 2011. Genetic analyses indicated that distinct evolutionary lineages of Flu A(H1N1)pdm2009, Flu A/H3N2 and Flu B have co-circulated in Bolivia in the study period, originating from Central and North America, Europe, Asia and Australia. CONCLUSION: Our results emphasise the requirement for a reinforced epidemiological and genetic follow-up of influenza and other ILIs in Bolivia to further inform the preparation of vaccines used in the region, guide vaccination campaigns and improve the medical management of patients.

Acute respiratory tract infections (ARIs) represent a leading cause of human acute illnesses worldwide and an important contributing factor of childhood morbidity and mortality, especially under the age of 5 [1] [2] [3] . Besides the issue of bacterial pneumonia, which represents a frequent cause of mortality in paediatric populations of developing countries, there is an increasing interest in the epidemiology and characterization of Influenza Like Illnesses (ILIs), which are predominantly of viral aetiology [4] . First, ILIs are responsible for a large number of cases in all age groups (with a specifically high frequency in children) and despite their frequent association with a mild clinical presentation, their medical and socio-economical impact is enormous [5] . Second, a proportion of ILIs is associated with severe clinical evolution [6] . In adults, viruses are the putative causative agents in a third of cases of community-acquired pneumonia, in particular influenza viruses, rhinoviruses, and coronaviruses. Dual viral infections are common, and a large proportion of children have evidence of viral-bacterial co-infection [4] . Deciphering the epidemiological and pathophysiological links between viral primary infection and secondary bacterial super-infection and pneumonia currently represents a major challenge for improving the prevention and medical management of pneumonia.

Despite their importance in terms of morbidity as well as infant mortality, microbiological aetiologies of ARIs represent a complex, which has yet to be fully characterized in developed countries and remains largely unknown in developing countries. Despite recent improvements due to prospective studies launched during or following the influenza A(H1N1)pdm09 in Peru and Ecuador [7, 8] , the information about ILIs in the tropical Hispanic regions of South America is scarce and the epidemiology remains insufficiently characterised.

Here, we have performed the first study of ILIs in Bolivia, which covers 3 years from January 2010 to September 2012. The viral aetiology of infections was investigated based on naso-pharyngeal swabs collected in a National Reference Laboratory in the city of Santa Cruz de la Sierra and real time PCR molecular biology techniques, viral cultures and phylogenetic analyses.

The number of samples that could be analysed each year at Centro Nacional de Enfermedades Tropicales, CENETROP was determined according to logistical considerations (257 in 2010, 163 in 2011 and 144 in 2012; total = 564). We did not attempt to build a representative sample of specimens received at CENETROP (supposedly poorly representative of ILIs in the general population of Santa Cruz), but rather examined a similar number of samples in each age group to allow a comprehensive aetiological analysis in all age classes. Our final sample has equivalent numbers of swabs in the 0-9, 10-19, 20-29 and 30-39 age classes (90 to 98 per class) but slightly lower numbers in the older classes (40-49: 70; 50-59: 59; ≥60: 61) due to the limited number of samples received at CENETROP for these age groups. As expected, it under-represents the 0-9 years old (yo) age group (which includes~35% of eligible samples received), and over-represents the 30-39, 40-49, 50-59 and ≥60 yo age groups (Additional file 1). The general characteristics and distribution of cases are presented in Table 1 . The gender ratio (M/F) is 0.84 and the history of influenza vaccination clearly shows that a limited proportion of the population received influenza vaccines, in particular in the elderly (0.53% have received at least one dose of influenza vaccine in their life in the ≥60 yo age group).

Of the 564 samples tested, 283 (50.2%) were positive for at least one virus (Table 1) with influenza viruses (Flu A:~15%; Flu B:~9%), rhinoviruses (hRV) (~8%), coronaviruses E229&OC43 (hCoV E229&OC43) (~5%) and respiratory syncytial virus hRSV (~4%) being the most frequently identified. The pattern of viral infections varied according to age groups ( Figure 1A ). The elucidation rate was the highest (>60%) in samples from patients under the age of 10 yo and the lowest (<40%) in patients ≥60 yo (Figures 1A and B) . Since all samples originated from symptomatic patients, this may be due either to the higher prevalence of the viruses tested or to the higher level of viral shedding in children.

The 0-4 yo age group displayed specific epidemiological characteristics with the majority of hRSV infections (infection rate: 23.19%) and a remarkably low proportion of influenza infections: <7.5% compared with~38% in the 5-9 yo age group,~26.5% in the 10-59 yo age classes and~16% in those ≥60 yo.

Adenoviruses were most frequently identified in children under the age of 10, bocaviruses in patients under the age of 20 and hCoV E229 & OC43, although present in all age groups, were more prevalent in patients ≥40 yo. Interestingly, metapneumoviruses (hMPV) were detected in both children and adults with the most elevated rate of infection in the ≥60 yo age group (~5%) and enteroviruses/parechoviruses were also found in both children and adults. Rhinoviruses were identified in all age groups (up to~10% in patient over 60 years of age).

Nearly 3% of samples showed dual viral infections ( Table 1 ): 8 (9.4%) of the 85 influenza A viruses, 5 (11.4%) of the 44 rhinoviruses and 5 (17.8%) of the 28 hCoV E229&OC43 infections detected were diagnosed in a situation of co-infection. Co-infections were specifically frequent in the case of adenoviruses (38.9%). Despite a high number of influenza B cases (50) no co-infection was noted. The rate of dual infections was high for bocaviruses and enteroviruses/parechoviruses (≥40%) and low for metapneumoviruses and parainfluenzaviruses (<10%) but this was based on a limited number of cases. The median age was not different in patients with single or dual infections.

The curve of positive samples followed that of samples tested ( Figure 2B ) and that of samples received at CENETROP (Figure 2A ). Peaks of detection matched peaks of reception of samples, which globally followed the seasonal subtropical pattern of ILIs observed in Bolivia, with one peak during the austral autumn and a second during the austral winter season. In 2010, a great number of samples could be studied (257) allowing to precisely follow several peaks of activity corresponding to a mosaic of viruses including Flu B, hRSV, hRV, Flu A A/H3N2 and A (H1N1)pdm09 and hCoV E229&OC43. In 2011 and 2012, the number of samples studied annually was lower (163 and 144, respectively), but still sufficient to elucidate the main peaks of activity with a major incidence of Flu A (H1N1)pdm09 in 2011 and Flu B in 2012, co-circulating with a variety of other viruses (see Figure 2C ). Epidemiological peaks were generally associated with a predomin- This is summarised in Figure 1C that shows a clear epidemiological evolution. In 2010, following the 2009 Flu A(H1N1)pdm09 pandemic wave, influenza A infections were for the most part due to Flu A/H3N2 and associated to a large number of Flu B and rhinovirus infections and a variety of less prevalent viruses. By contrast, year 2011 was characterised by the massive come back of Flu A(H1N1)pdm09, the low number of Flu B, the emergence of a large number of hCoV E229&OC43 infections and the absence of hRSV infections. In 2012, the situation changed again with a large number of Flu B and hRSV infections and the persistence of a significant proportion of rhinovirus and Flu A(H1N1)pdm09 infections. Altogether, this provides a striking picture of a moving infection pattern, which is presumably modelled by the herd immunity in the most epidemiologically susceptible populations. The complete disappearance of hRSV infections in 2011 represents a specific case that was never encountered by the authors in Europe.

Cell culture isolates of influenza A/H3N2 (N = 4), influenza A(H1N1)pdm09 (N = 19) and influenza B2 (N = 11) viruses were obtained. Eight-segment genomic sequences were produced for Flu A/H3N2 isolates (GenBank accession numbers KF612026-KF612180) and Flu A/H3N2 isolates (GenBank accession numbers KF612181-KF612211). Partial sequences (segments HA, NA, NS, M) were produced for Flu B isolates (GenBank accession numbers KF612212-KF612267). The corresponding sequences were used for further phylogenetic analyses.

Influenza B viruses are divided into two antigenic lineages (Victoria/87 and Yamagata/88) since 1983. Both of them frequently co-circulate during epidemic periods [9] and segment reassortment between strains of both lineages is frequent, making it difficult to produce robust complete-genome phylogenetic analyses. Here, evolutionary trees were reconstructed from the alignment of 485 HA sequences and showed that ten out 

These data suggest that variants belonging to distinct evolutionary clusters have circulated simultaneously in South America during the 2011 and 2012 influenza seasons. This information is relevant for optimising the formula of the seasonal influenza vaccine in order to provide efficient protection [10] .

Influenza A(H1N1)pdm09 Figure 5) .

Altogether, the main information arising is that several evolutionary lineages co-circulated in Bolivia during the 2009 H1N1 pandemic wave.

This study represents the first attempt to characterise the viral aetiology of ILIs in Bolivia and the first in South America in the post-influenza A(H1N1)pdm09 period. It was not designed for providing a representative picture of the general Bolivian population: (i) samples originated from the region of Santa Cruz only (this item is further discusses below); (ii) samples sent to CENETROP were not expected to faithfully represent ILIs in the region of Santa Cruz, with obvious recruitment biases such as clinical severity or access to public medical care; (iii) sampling was performed with the objective of obtaining equivalent numbers of cases in all age classes.

The information obtained should therefore be analysed with precautions. Firstly, the great year-to-year variation of infection epidemiological patterns is a reminder that the epidemiology of ILIs should be studied over long periods of time. Each of the three years studied here would have provided, if analysed separately, a partial view of the situation and would have lead to erroneous conclusions (e.g., the absence of hRSV infections in Bolivian children extrapolated from 2011 data). The analysis of a 3-year period is expected to smooth such epidemiological jolts but an extended observation period would have undoubtedly revealed more information about the viruses that circulated at low rates in the period studied (e.g., parechoviruses, parainfluenzaviruses). Secondly, Bolivia includes territories with diverse geographical, climatic and ecological characteristics. Santa Cruz is located in an eastern low-altitude area with a humid tropical climate whilst La Paz is in a western Andean high-altitude area with a colder and dryer climate [11] . This may explain important epidemiological differences observed between both cities, as exemplified by the case of influenza A(H1N1)pdm09: in 2010, we observed a very limited number of Flu A(H1N1)pdm09 infections in Santa Cruz whilst the PAHO surveillance system reported an intense circulation of the virus in La Paz during the first months of the year (PAHO surveillance website: http://ais.paho. org/phip/viz/ed_flu.asp); in 2011, a large outbreak occurred in Santa Cruz (see Figures 1 and 2) , whilst only a few cases were observed in La Paz; in 2012, an important outbreak was observed in La Paz, but not in Santa Cruz. This clearly indicates that epidemiological data collected in Santa Cruz cannot be extrapolated easily to the Andean regions of Bolivia. Divergent influenza epidemiological patterns according to geographical distribution were finely documented in a previous Peruvian study [7] .

Despite the fact that the local epidemiology of ARIs is modelled by herd immunity and environmental factors, the global picture of ARIs that arises from our study conforms with that observed in different regions of the world. In particular, the choice of the pathogens tested seems to be relevant since (i) the rate of elucidation is close to what is observed in France, using the same molecular assays (authors' personal data), (ii) peaks of detection match the ILI epidemiological peaks and (iii) in previous studies performed in Brazil, the different viruses identified in the current study were detected, confirming their circulation in the region [12] [13] [14] [15] [16] . The distribution of cases in age groups is also essentially similar with what was reported elsewhere, e.g., in Ecuador [8] . The most surprising items, which deserve investigations on an extended period of time, were the unusually low frequency of influenza infections in the 0-4 yo age group and the complete RSV eclipse in 2011. However, precise comparison with data collected in South America is most often difficult, e.g., with a molecular study in Ecuador [8] and a previous large study in Peru [7] (different periods investigated, viruses detected using immunofluorescence tests, parechoviruses, coronaviruses, metapneumoviruses and rhinoviruses not tested).

The genetic characterisation of influenza strains indicates that viruses belonging to distinct evolutionary lineages of Flu A(H1N1)pdm09, Flu A/H3N2 and Flu B have co-circulated in Bolivia. Interestingly, an increasing number of influenza sequences has been made available since 2008 and phylogenetic reconstructions revealed, on the one hand, the existence of South American evolutionary clusters of influenza strains, and on the other hand a complex epidemiological relationship with the rest of the world. Trees suggested the circulation in Bolivia during the study period of viruses originating from Central and North America, Europe, Asia and Australia.

This study provided the first systematic and broadspectrum information regarding the viral aetiology of ILIs in Bolivia. A continued effort is needed to accurately characterise the epidemiology of ILIs in the different regions of the country to guide vaccination campaigns for at-risk populations and improve the medical management of patients. In accordance with previous investigations performed in Peru [7] , our results emphasise the requirement for a reinforced genetic follow-up of influenza strains circulating in South America to further inform the preparation of vaccines used in the region.

The recruitment of cases was operated following the broad spectrum WHO case definition of ILIs: anybody with fever of 37.8°C or higher, associated with one respiratory symptom was identified as ILI infection and eligible for investigation. We used naso-pharyngeal swabs sampled in hospitals and health centres in the region of Santa Cruz and sent to the CENETROP (Centro Nacional de Enfermedades Tropicales, Santa Cruz, Bolivia) for diagnosis from January 2010 to September 2012. All specimens had been sampled between day 0 and day 4 following clinical onset, sent to the laboratory at room temperature using a universal transport media and kept at −20°C upon reception. All patients or their legal representative gave oral consent for the shipment of samples to CENETROP and the detection of respiratory pathogens. To provide a preliminary report of viral aetiologies in all age groups (i.e., information that may be translated easily for Public health purposes), a random selection of an equivalent number of samples was operated in age groups (0-9, 10-19, 20-29, 30-39, 40-49, 50-59 and ≥60 years of age) for each year of the study.

Aetiological identification was performed using TaqMan probe-based real-time PCR and RT-PCR protocols. Samples were spiked by MS2 and T4 phages, used as internal controls [17] . RNA [20] ), Respiratory syncytial viruses A and B (hRSV A, hRSV B, [21] ), Coronaviruses E229&OC43 (hCoV E229&OC43, [22] ), Rhinovirus (hRV, [23] , Adenoviruses (hAdV, [24] ), Enteroviruses & Parechoviruses (hEVs, hPeV, [25] ), human metapneumoviruses A and B (hMPV, [26] ), human Bocaviruses (hBoV, [27] ).

Samples providing a strong positive real-time PCR result (Ct < 30) for influenza A or B viruses were used for virus isolation using MDCK cells and a standard isolation protocol (Additional file 11).

In addition to strains isolated in the course of the current study, 6 influenza A(H1N1)pdm09 virus strains isolated at CENETROP in the June-August 2009 period were sequenced to reinforce phylogenetic analyses (see below).

Sequencing 200 μL of clarified supernatant from Influenza A or B positive cell cultures were extracted using the EZ1® Qiagen instrument and the EZ1® Virus Mini Kit v2.0 without RNA Carrier (Qiagen). All 8 segments of influenza A virus were amplified in a single reaction using the SuperScript® III One-Step RT-PCR System with the Plat-inum® Taq High Fidelity Kit and the universal primers univ12 and univ13 [28] .

To obtain influenza B virus sequences, we used a similar protocol and one set of primers (BM-NS-1 and BM-NS-2) to amplify the HA, NA and NS segments [29] and another set of newly designed primers (M-2eFor: AG AAGYASAGCATTTTCTTGTGA and M-2eRev: TAAA CACCCACATYCCAAACGT) to amplify the M segment.

PCR products were purified (Nucleofast PCR plate, Machrey Nagel) and sequenced using the NGS Ion Torrent technology: briefly, DNA concentrations were adjusted to 200nM (Qubit® dsDNA BR Assay Kit and Qubit® 2.0 Fluorometer) and a library of barcoded 100-200 nt fragments was built using the Diagenode Bioruptor UCD-200, the Xpress™ Barcode Adapters kit and the AB Library Builder™ (Life Technologies). Clonal amplification was performed using the Ion PGM™ Template OT2 200 Kit and the Ion OneTouch™ 2 Instrument (both from Life Technologies). Final sequencing was performed using the Personal Genome Machine® (PGM™) System, the Ion PGM™ Sequencing 200 Kit v2 and the Ion 316™ Chip Kit (Life Technologies).

Estimates of world-wide distribution of child deaths from acute respiratory infections

WHO estimates of the causes of death in children

Global burden of acute respiratory infections in children: implications for interventions

Viral pneumonia

Rhinoviruses: markers of, or causative for, recurrent wheeze and asthma?

Child Health Epidemiology Reference Group of WHO and UNICEF: Global, regional, and national causes of child mortality: an updated systematic analysis for 2010 with time trends since

Epidemiology of influenza-like illness in the Amazon Basin of Peru

Sentinel surveillance of influenza-like-illness in two cities of the tropical country of Ecuador

The evolutionary dynamics of human influenza B virus

Writing Committee of the World Health Organization Consultation on Southern Hemisphere Influenza Vaccine Composition for 2012: WHO recommendations for the viruses to be used in the 2012 Southern Hemisphere Influenza Vaccine: Epidemiology, antigenic and genetic characteristics of influenza A(H1N1)pdm09, A(H3N2) and B influenza viruses collected from

Global influenza seasonality: reconciling patterns across temperate and tropical regions

Acute respiratory infection and influenza-like illness viral etiologies in Brazilian adults

Sentinel surveillance of influenza and other respiratory viruses, Brazil

Epidemiological and clinical features of human coronavirus infections among different subsets of patients

Acute respiratory viral infections in children in Rio de Janeiro and Teresopolis

Full genome sequence analysis of parechoviruses from Brazil reveals geographical patterns in the evolution of non-structural genes and intratypic recombination in the capsid region

RNA and DNA bacteriophages as molecular diagnosis controls in clinical virology: a comprehensive study of more than 45,000 routine PCR tests

A simple method for molecular detection of Swine-origin and human-origin influenza a virus

Comparison of real-time PCR assays with fluorescent-antibody assays for diagnosis of respiratory virus infections in children

Applicability of a real-time quantitative PCR assay for diagnosis of respiratory syncytial virus infection in immunocompromised adults

Frequent detection of human coronaviruses in clinical specimens from patients with respiratory tract infection by use of a novel real-time reverse-transcriptase polymerase chain reaction

Lower respiratory viral illnesses: improved diagnosis by molecular methods and clinical impact

Pring-Akerblom P: Rapid and quantitative detection of human adenovirus DNA by real-time PCR

Rapid simultaneous detection of enterovirus and parechovirus RNAs in clinical samples by one-step real-time reverse transcription-PCR assay

Real-time reverse transcriptase PCR assay for detection of human metapneumoviruses from all known genetic lineages

Real-time PCR assays for detection of bocavirus in human specimens

Full genomic amplification and subtyping of influenza A virus using a single set of universal primers

Screening of the high yield influenza B virus on MDCK cell and cloning of its whole genome

MUSCLE: multiple sequence alignment with high accuracy and high throughput

MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods

Viral aetiology influenza like illnesses in

To produce consensus sequences for each isolated strains, reads obtained from using the Ion Torrent protocol were processed using the CLC Genomics Workbench software 6.04 (CLCbio, Aarhu, Denmark) and a home developed program and then mapped onto concatenated reference sequences.In the case of A/H1N1(2009) phylogenies, concatenated coding sequences of the eight segments from the complete genomic sequences produced as reported above from Bolivian strains were aligned with complete genomic sequences retrieved from databases (total: 3,412 genomic sequences). In the case of A/H3N2 and influenza B phylogenies we collected complete genomic sequences from databases but used only the HA gene sequences to produce alignments (A/H3N2: 2,643 sequences; influenza B: 485 sequences). Alignments were conducted using the MUSCLE programme [30] , implemented in the Influenza Research Database (Multiple Sequence Alignment: http:// www.fludb.org/brc/) and subsequently used to perform phylogenetic analyses using the MEGA 5.05 programme [31] , the Jukes Cantor algorithm for distance calculation and the Neighbour joining method for tree building.

Additional file 1: Distribution of the eligible and studied populations in age groups.Additional file 2: Influenza B phylogeny cluster 1. This is a magnification of Influenza B phylogeny cluster 1 ( Figure 3 ). Legend shape represents the year of strain isolation: square for 2012; circle for 2011; triangle with the tip up for 2010, triangle with the tip down for 2009; rhombus for <2009. The colour represents the geographical origin: light blue for Bolivian strains from this study; pink for Bolivian strains from databases; dark blue for South American strains from databases. 

The authors declare that they have no competing interests.Authors' contributions JD: Participate on sampling samples molecular biology, culture and sequencing parts and draft of paper; YRS: Participate to contributions to conception and design of project. GP: have made substantial contributions to conception and design of influenza sequencing participated in the sequence treatment alignment. MB: helps analysis and interpretation of data. CB; have made substantial contributions to conception and design of influenza sequencing. LTP; have been involved in drafting the manuscript. RLM: Participate on molecular biology and viral culture. CAA: Participate on sampling samples molecular biology. GAA: Participate on sampling samples molecular biology. JRG; Participate to contributions to conception and design of project. JLR: Participate to contributions to conception and design of project. XdL: have given correction and final approval of the version to be published. All authors read and approved the final manuscript.

Submit your next manuscript to BioMed Central and take full advantage of: