- Open Access
Rhinovirus is an important pathogen in upper and lower respiratory tract infections in Mexican children
Virology Journalvolume 12, Article number: 31 (2015)
Most of the studies characterizing the incidence of rhinovirus (RV) have been carried out in hospitalized children and in developed countries. In those studies, RV-C has been associated with more severe respiratory tract infections than RV species A and B. In this study we determined the frequency and diversity of RV strains associated with upper and lower respiratory tract infections (URTI, LRTI) in Mexico, and describe the clinical characteristics of the illness associated with different RV species.
A prospective surveillance of 526 and 250 children with URTI and LRTI was carried out. Respiratory samples were analyzed by RT-PCR for viruses. The 5′ untranslated region of the RV genome was amplified and sequenced.
In the case of URTI, 17.5% were positive for RV, while this virus was found in 24.8% of LRTI. The RV species was determined in 73 children with URTI: 61.6% were RV-A, 37% RV-C and, 1.4% RV-B; and in 43 children with LRTI: 51.2% were RV-A, 41.8% RV-C, and 7% RV-B. No significant differences in clinical characteristics were found in patients with RV-A or RV-C infections. A high genetic diversity of RV strains was found in both URTI and LRTI.
Both RV-A and RV-C species were frequently found in hospitalized as well as in outpatient children. This study underlines the high prevalence and genetic diversity of RV strains in Mexico and the potential severity of disease associated with RV-A and RV-C infections.
Pneumonia continues to be a major killer of young children in developing countries and elderly people in developed countries. Despite recent advances, further studies are needed to examine regional variation in its etiology, particularly in developing countries, where most of the deaths from serious respiratory diseases occur [1,2].
Rhinovirus (RV) is the most frequent cause of acute respiratory illness worldwide [3-5]. This virus has been typically associated with upper respiratory tract infections (URTI); however, with the development of molecular methods RV has been found to be also commonly associated with lower respiratory tract infections (LRTI). These findings are changing the long held view of RV as a minor pathogen, as it is now being involved in a wide variety of respiratory illnesses, ranging from mild common colds and asthma exacerbation to serious lower respiratory tract disease [3-5].
RVs are antigenically quite diverse; historically, these viruses were classified into 99 serotypes through neutralization assays. These different RV serotypes were later classified genotypically into two species, RV-A (74 serotypes/genotypes) and RV-B (25 serotypes/genotypes) . More recently, based only on genomic sequence information, a third virus species, RV-C, was described [6-9]. The current genotypification of these viruses is based on the nucleotide sequence encoding the VP1 and VP4/VP2 proteins . Nevertheless, a variable region in the 5′-untranslated region (UTR) of the viral genome can also accurately distinguish among virus species . Since its description in 2006, RV-C was suggested to be associated with more severe respiratory illness as compared to RV-A and RV-B, as well as with more frequent asthma exacerbations [11-16]. However, other studies have found no difference in illness severity among RV species [17-20], thus, more information is needed to clarify this issue.
Most studies have described the presence of RV genotypes in hospitalized patients with severe respiratory illness, and only a few studies have described the prevalence of virus genotypes in URTIs. In this work, we carried out a prospective multicenter study of two children populations having either URTI or LRTI.
Nasopharyngeal samples from 526 pediatric patients with URTI attending the private consult were tested for the presence of RV and other respiratory viruses. Ninety-two (17.5%) children (47 males and 45 females, age range 0 to 175 months) were RV positive. The species of 73 (80.2%) of these RVs was determined by sequencing the 5′-UTR region of the viral genome; 45 (61.6%) were RV-A, 27 (37%) RV-C, and 1 (1.4%) was RV-B (Table 1). The incidence of RV infection was highest in September and November 2011 as well as in April 2012, when it accounted for 24.4%, 27.1% and 30.3%, respectively, of the total samples collected (Figure 1).
On the other hand, 250 hospitalized children with clinical diagnosis suggestive of viral pneumonia were included in the study. Nasal washings were obtained and tested for the presence of respiratory viruses. Sixty-two children (41 males and 21 females, age range 1 to 76 months) were positive for RV (24.8%). The species of 43 (69.4%) of these viruses was determined by sequencing the 5'-UTR region of the viral genome. Of the typed viruses, 22 (51.2%) were RV-A, 18 (41.8%) RV-C, and 3 (7%) were RV-B (Table 1). The difference in the incidence of RV-A and RV-C was not statistically significant. RV-positive children were detected more frequently (p < 0.0001) in children between 25 and 36 months as compared to the other age groups (Table 1). There was no significant gender between children infected with RV or another respiratory virus or between RV-A and RV-C. Hospitalization for RV was significantly higher (p < 0.0001) in the spring as compared to other seasons.
Co-infection of RV with other respiratory viruses
An additional respiratory virus was found in 30% (28/92) and 37% (23/62) of RV-positive patients with URTI and LRTI, respectively (Table 1). The virus species could be determined in 19 (14 RV-A, 5 RV-C) and 14 (7 RV-A, 1 RV-B, 6 RV-C) of URTI and LRTI mixed infections, respectively. Neither RV-A nor RV-C types were significantly associated with a particular respiratory virus in mixed infections of URTI or LRTI. However, in URTI, RV-A was found more frequently in co-infections than RV-C (p = 0.025). In both children populations, when RV was present in mixed infections, it was more frequently found associated with a single virus (URTI, 20/28; LRTI, 21/23) than with 2 other viruses (URTI, 8/28; LRTI, 2/23) (URTI p = 0.014, LRTI p < 0.0001).
Clinical characteristics of RV-positive children
The major symptoms and signs observed in children with RV-associated URTI in single or mixed infections were cough, rhinorrhea, and sore throat. However, rhinorrhea was the only symptom significantly higher in RV-positive children as compared to patients with infections caused by other respiratory viruses (Table 2). Increased respiratory rate was found in more than one third of the patients with a single or mixed infection with RV, as well as in infections caused by other viruses. The clinical features of children having a co-infection with RV and other respiratory virus were no statistically different from those with single RV infections. Also, no significant difference was found between patients infected with RV species A or C. RV-B was excluded from the analysis given that a single child was infected with this virus.
In 92.7% of the RV-positive patients with LRTI the diagnosis of pneumonia was confirmed radiologically, as it was in 94.5% of children positive for other viruses (Table 3). Children with RV had fever >38°C and increased respiratory rate more frequently than children with any other respiratory virus. This was true for single and mixed RV infections. Other major signs detected upon examination of RV-positive patients were thoracoabdominal dissociation, and intercostal retraction. According to the Silverman-Anderson score there was no difference between children with RV or any other virus; 76.4% of the patients with RV co-infections and 82.3% of RV single infections had a score between 1 and 3. In the chest-X ray examination, 78.1% of the RV-positive patients had signs of interstitial pneumonia; these percentages were similar for RV-single infections. There was no difference between the clinical signs observed in children infected with either RV-A or RV-C.
Genetic diversity of RV
The phylogenetic trees in Figure 2 depict the wide distribution of 5′-UTR sequences of the RV-A and RV-C strains isolated from URTI and LRTI patients. The 112 characterized RV strains could be classified into 35 different 5′-UTR-based genotypes (5′-genotypes); in this work, two 5′-genotypes were defined as different when their 5′-UTR sequence identity was equal or less than 87%, similar to the cutoff determined for genotypes based on VP1 nucleotide sequences . Seventeen genotypes were associated with RV-A and 18 with RV-C. Each genotype included from 1 to 14 different virus strains. In the case of RV-A, 5′-genotypes were significantly more frequently found in either hospitalized (genotype 2, p < 0001) or outpatient (genotype 6, p < 0.001) children, while for RV-C some 5′-genotypes also associated more frequently with either URTIs (genotypes 2, p < 0.001 and 18, p = 0.03) or LRTIs (genotype 6, p < 0.01) (Figure 2). Importantly, the correlation of 5′-genotypes with reference genotypes defined by full genomic sequences  was remarkably high. Comparison of the phylogenetic trees constructed with either the 5′-genotypes or the full genomic sequences showed that 83% (12/72) of the terminal clades present in the 5′-phylogenetic tree were the same as compared with the terminal clades of the reference tree (Figure 2). These results suggest that 5′-genotypes could serve as a simpler and initial approach to genomically classify RVs. Further studies need to be conducted to determine if the association between some RV 5′-genotypes and the severity of infection observed in this work holds, or if it is the result of particular RV genomic types that circulated more frequently in the years that URTI and LRTI samples were collected.
Here, we report the frequency of RV in pediatric patients with URTI and LRTI in different regions of Mexico. The frequencies observed in this work are in agreement with those reported in other studies, in which RV has varied from 16% to 37% in children with either URTI or LRTI. Regarding seasonality, RV infections were associated most frequently to URTI in autumn and early spring, as reported , although it should de noted that the samples from patients with URTI miss the months from May to August. On the other hand, a clear peak of prevalence of the virus in LTRI was observed only during spring.
An additional respiratory virus was found in about one-third of the RV-positive patients with either URTI or LRTI. Detection of several viruses in a high proportion of cases has been a characteristic of respiratory infections in which a PCR-based diagnostic method was used. In particular, RV has been found in dual infections with another respiratory virus in 26% to 63% of cases [14,18,20]. However, the clinical relevance of detection of several viruses in pneumonia, and the association with severe illness is uncertain [2,23-25].
Despite technological advances, establishing the cause of pneumonia remains challenging . In this study diagnosis of viral pneumonia relied on nasopharyngeal specimens, what might be misleading since detection of a virus in the nasopharynx could represent a coincidental upper-respiratory infection, the asymptomatic presence of the virus, or a genuine pneumonia pathogen. Measurement of background prevalence of asymptomatic nasopharyngeal viral infections in a healthy control group and in patients with mild-to-moderate disease might help to clarify the relevance of this diagnostic issue at a population level . In most of this type of studies, but not in all, a higher RV detection among participants with illness was noted [11,13,19,27-32].
The overall identification of RV was similar but statistically different between the two patient populations studied. Children with URTI had a RV rate of 17% as compared to 24.5% in patients with LRTI (p < 0.017). The observation that the frequency of RV detection in serious respiratory disease is similar to that observed in mild-to-moderate illness, or even in asymptomatic cases [11,13,19,27-32], is intriguing. One would think that if the respiratory disease caused by RV was more severe than other pathogens, the percentage of hospitalized children positive for RV should be significantly higher than that found in outpatient or healthy children, as observed, for instance, in rotavirus infections . On the contrary, the incidence of RV should have to be significantly lower in LRTI as compared to URTI if the virus infection caused a mild disease. Then, why is there about the same incidence of RV found in respiratory infections with different degrees of severity? It is tempting to hypothesize that specific virus serotypes/genotypes are preferentially associated to severe or mild clinical symptoms. This potential association might be obscured by a wide diversity of circulating RV serotypes/genotypes with potentially different intrinsic virulence, such that in different children populations similar frequencies of RV infection are observed but with different clinical outcomes. It cannot be discarded, however, that mixed infections with bacteria (not searched in this work) or different host factors could influence the clinical severity of RV infections, as reported . Of interest, there have been suggestions since the earliest epidemiological studies of variation of virulence among the different RV serotypes [35,36] and, in this line, we noted in this work that some 5′-genotypes were detected only in either hospitalized or outpatient children. It would be important to characterize the diversity of RV genotypes in different parts of the world to determine if the prevalent viruses vary among regions and in time, and associate with different severity of the disease. Identification of this potential association would have a great impact in the development of prophylactic measures to control the infection of this very common pathogen.
In conclusion, this study underlines the high RV exposition and diversity of circulating strains in Mexico and the potential severity of disease associated with both RV-A and RV-C infections.
This study describes the frequency of detection of rhinovirus species in children with upper and lower respiratory tract infections in Mexico and their genetic diversity, determined by sequencing the 5′ UTR region of the viral genome. We report that the 5′UTR sequence can be useful for an initial approach to determine the rhinovirus genotypes. We also describe the clinical characteristics of illness associated with specific RV species. We found that both RV-A and RV-C species were very frequently found in both hospitalized and outpatient children and no statistically significant differences were found in the severity of disease associated with RV-A and RV-C infections in neither of the two children populations. This study underlines the high RV prevalence and genetic diversity of circulating strains in Mexico and the potential severity of disease associated with both RV-A and RV-C infections. This is of particular relevance, since the information about respiratory viruses in Mexico is very limited, and studies characterizing viruses circulating in the community level are even scarcer.
Study populations and clinical samples
Two patient populations were included in this study. The first was composed of pediatric patients that attended the private consult in five different cities of the state of Veracruz, Mexico (Poza Rica, Veracruz, Córdoba, Tierra Blanca, and Minatitlán). The children were enrolled in the study if they were younger than 16 years, clinically diagnosed as having an acute respiratory infection, with an onset of illness less than one week, and the parent or guardian signed the informed consent form. Since none of these children required hospitalization, they were all considered as having URTIs. From September 2011 to April 2012 nasopharyngeal swabs (rayon-tipped, BD BBL) were collected from 526 children (male:female ratio, 1.27; median age, 19 months; average age, 39 months; age range, 0–191 months). The second population consisted of children that required hospital admission with suspected pneumonia in four cities (Guadalajara, Mexico City, San Luis Potosí, and Durango) in different states of Mexico. These children were all considered to have LRTIs. Nasal washings were collected from 250 children (male:female ratio, 1.43; median age, 10 months; average age, 16 months; age range, 1–76 months, with only one child 143 months old) between March 2010 and April 2011. For both populations, demographic and clinical information was collected. The samples were placed in vials containing viral transport medium (Microtest M4-RT, Remel, Lenexa, KS) and sent either to the Institute of Biotechnology in Cuernavaca (URTI samples) or to the School of Medicine in Mexico City (LRTI samples) and stored at −70°C until analyzed. The study was approved by the institutional review boards of the School of Medicine and the Institute of Biotechnology of the National University of Mexico and from the institutional review board and ethics committee of each participant hospital. Written informed consent was obtained from each parent or guardian prior to enrollment.
Nucleic acid extraction and RT-PCR assay
Nucleic acids were extracted from 100 μl of respiratory specimens using the PureLink Viral RNA/DNA mini kit (Invitrogen, Carlsbad, CA). RV and other fourteen respiratory viruses were detected with the Seeplex®RV15 OneStep ACE Detection kit (Seegene, Seoul, South Korea) following the manufacturer’s instructions. This diagnostic test detects the following viruses: parainfluenza viruses 1, 2, 3, and 4, adenovirus A/B/C/D/E, human coronaviruses 229E/NL63 and HuCoV-OC43, rinovirus A/B/C, influenza A, influenza B, respiratory syncyctial virus-A and -B, human bocavirus 1/2/3/4, human metapneumovirus, and enterovirus. The RNA present in the samples positive for RV was reverse transcribed with random hexamers using standard protocols. PCR was then performed with previously described primers DK001 and DK004, which target a fragment of approximately 400 bp of the hypervariable region of the RV 5′-UTR virus genome . The amplified DNA fragment was purified with the High Pure PCR Product Purification kit (Roche).
The purified PCR products were sequenced in an Applied Biosystems (Foster City, CA), model 3130xl apparatus. The sequences obtained were edited using the Unipro UGENE and Seaview softwares. A database of RV complete genomes was created using all sequences available on GenBank until September 2013. BLAST was performed and the genotype of the RV strains was assigned based on the first five best scores. The 5′ sequences of the RV strains characterized in this work were deposited in GenBank with accession numbers KJ765008 to KJ765123.
Multiple sequence alignments were made using MUSCLE method; maximum likelihood trees were generated with 1000 repetitions bootstrap using the MEGA 5.2.2 program. To determine the 5′-genotypes, clustering was made using the cd-hit-test at 87% of identity.
All statistical data analyses were carried out using PASW statistical software version 18.0 (SPSS Inc.). Significant differences between groups were evaluated using chi-square tests when it was possible, or Fisher’s exact tests when values were smaller than 25. Associations of demographic and clinical features of children with a) RV detected versus other respiratory virus and, b) RV-A versus RV-C, were examined using logistic regression. To determine the 5′-genotypes that were significantly enriched for RV-A or RV-C, the probability of randomly selecting different viral strains in each 5′-genotype was calculated by a combinatorial approach. In all statistical analysis p-values <0.05 were considered statistically significant.
Written informed consent was obtained from each parent or guardian prior to enrollment in the study.
Black RE, Cousens S, Johnson HL, Lawn JE, Rudan I, Bassani DG, et al. Global, regional, and national causes of child mortality in 2008: a systematic analysis. Lancet. 2010;375:1969–87.
Ruuskanen O, Lahti E, Jennings LC, Murdoch DR. Viral pneumonia. Lancet. 2011;377:1264–75.
Gern JE PA. Rhinoviruses. In: David M, Knipe PMH, editors. Fields Virology, vol. 1. Philadelphia: Wolters Kluger/Lippincott Williams & Wilkins; 2013. p. 531–49.
Heikkinen T, Jarvinen A. The common cold. Lancet. 2003;361:51–9.
Pitkaranta A, Hayden FG. Rhinoviruses: important respiratory pathogens. Ann Med. 1998;30:529–37.
Palmenberg AC, Spiro D, Kuzmickas R, Wang S, Djikeng A, Rathe JA, et al. Sequencing and analyses of all known human rhinovirus genomes reveal structure and evolution. Science. 2009;324:55–9.
Lau SK, Yip CC, Tsoi HW, Lee RA, So LY, Lau YL, et al. Clinical features and complete genome characterization of a distinct human rhinovirus (HRV) genetic cluster, probably representing a previously undetected HRV species, HRV-C, associated with acute respiratory illness in children. J Clin Microbiol. 2007;45:3655–64.
Lee WM, Kiesner C, Pappas T, Lee I, Grindle K, Jartti T, et al. A diverse group of previously unrecognized human rhinoviruses are common causes of respiratory illnesses in infants. PLoS One. 2007;2:e966.
McErlean P, Shackelton LA, Lambert SB, Nissen MD, Sloots TP, Mackay IM. Characterisation of a newly identified human rhinovirus, HRV-QPM, discovered in infants with bronchiolitis. J Clin Virol. 2007;39:67–75.
Mubareka S, Louie L, Wong H, Granados A, Chong S, Luinstra K, et al. Co-circulation of multiple genotypes of human rhinovirus during a large outbreak of respiratory illness in a veterans’ long-term care home. J Clin Virol. 2013;58:455–60.
Calvo C, Casas I, Garcia-Garcia ML, Pozo F, Reyes N, Cruz N, et al. Role of rhinovirus C respiratory infections in sick and healthy children in Spain. Pediatr Infect Dis J. 2010;29:717–20.
Gern JE. The ABCs of rhinoviruses, wheezing, and asthma. J Virol. 2010;84:7418–26.
Khetsuriani N, Lu X, Teague WG, Kazerouni N, Anderson LJ, Erdman DD. Novel human rhinoviruses and exacerbation of asthma in children. Emerg Infect Dis. 2008;14:1793–6.
Lauinger IL, Bible JM, Halligan EP, Bangalore H, Tosas O, Aarons EJ, et al. Patient characteristics and severity of human rhinovirus infections in children. J Clin Virol. 2013;58:216–20.
Linder JE, Kraft DC, Mohamed Y, Lu Z, Heil L, Tollefson S, et al. Human rhinovirus C: Age, season, and lower respiratory illness over the past 3 decades. J Allergy Clin Immunol. 2013;131:69–77. e61-66.
Miller EK, Edwards KM, Weinberg GA, Iwane MK, Griffin MR, Hall CB, et al. A novel group of rhinoviruses is associated with asthma hospitalizations. J Allergy Clin Immunol. 2009;123:98–104.
Garcia J, Espejo V, Nelson M, Sovero M, Villaran MV, Gomez J, et al. Human rhinoviruses and enteroviruses in influenza-like illness in Latin America. Virol J. 2013;10:305.
Han TH, Chung JY, Hwang ES, Koo JW. Detection of human rhinovirus C in children with acute lower respiratory tract infections in South Korea. Arch Virol. 2009;154:987–91.
Iwane MK, Prill MM, Lu X, Miller EK, Edwards KM, Hall CB, et al. Human rhinovirus species associated with hospitalizations for acute respiratory illness in young US children. J Infect Dis. 2011;204:1702–10.
Xiang Z, Gonzalez R, Xie Z, Xiao Y, Liu J, Chen L, et al. Human rhinovirus C infections mirror those of human rhinovirus A in children with community-acquired pneumonia. J Clin Virol. 2010;49:94–9.
McIntyre CL, Knowles NJ, Simmonds P. Proposals for the classification of human rhinovirus species A, B and C into genotypically assigned types. J Gen Virol. 2013;94:1791–806.
Makela MJ, Puhakka T, Ruuskanen O, Leinonen M, Saikku P, Kimpimaki M, et al. Viruses and bacteria in the etiology of the common cold. J Clin Microbiol. 1998;36:539–42.
Calvo C, Garcia-Garcia ML, Blanco C, Vazquez MC, Frias ME, Perez-Brena P, et al. Multiple simultaneous viral infections in infants with acute respiratory tract infections in Spain. J Clin Virol. 2008;42:268–72.
Jennings LC, Anderson TP, Werno AM, Beynon KA, Murdoch DR. Viral etiology of acute respiratory tract infections in children presenting to hospital: role of polymerase chain reaction and demonstration of multiple infections. Pediatr Infect Dis J. 2004;23:1003–7.
Midulla F, Scagnolari C, Bonci E, Pierangeli A, Antonelli G, De Angelis D, et al. Respiratory syncytial virus, human bocavirus and rhinovirus bronchiolitis in infants. Arch Dis Child. 2010;95:35–41.
Murdoch DR, O’Brien KL, Scott JA, Karron RA, Bhat N, Driscoll AJ, et al. Breathing new life into pneumonia diagnostics. J Clin Microbiol. 2009;47:3405–8.
Fry AM, Lu X, Olsen SJ, Chittaganpitch M, Sawatwong P, Chantra S, et al. Human rhinovirus infections in rural Thailand: epidemiological evidence for rhinovirus as both pathogen and bystander. PLoS One. 2011;6:e17780.
Jartti T, Lehtinen P, Vuorinen T, Koskenvuo M, Ruuskanen O. Persistence of rhinovirus and enterovirus RNA after acute respiratory illness in children. J Med Virol. 2004;72:695–9.
Kusel MM, de Klerk NH, Holt PG, Kebadze T, Johnston SL, Sly PD. Role of respiratory viruses in acute upper and lower respiratory tract illness in the first year of life: a birth cohort study. Pediatr Infect Dis J. 2006;25:680–6.
Piotrowska Z, Vazquez M, Shapiro ED, Weibel C, Ferguson D, Landry ML, et al. Rhinoviruses are a major cause of wheezing and hospitalization in children less than 2 years of age. Pediatr Infect Dis J. 2009;28:25–9.
Singleton RJ, Bulkow LR, Miernyk K, DeByle C, Pruitt L, Hummel KB, et al. Viral respiratory infections in hospitalized and community control children in Alaska. J Med Virol. 2010;82:1282–90.
Wisdom A, Kutkowska AE, McWilliam Leitch EC, Gaunt E, Templeton K, Harvala H, et al. Genetics, recombination and clinical features of human rhinovirus species C (HRV-C) infections; interactions of HRV-C with other respiratory viruses. PLoS One. 2009;4:e8518.
Estes MK, Greenberg HB. Rotaviruses. In: Knipe DM, Howley PM, editors. Fields Virology, vol. 1. 6th ed. Philadelphia: Wolters Kluwer/Lippincott Williams & Wilkins; 2013. p. 1347–401.
Miller EK, Williams JV, Gebretsadik T, Carroll KN, Dupont WD, Mohamed YA, et al. Host and viral factors associated with severity of human rhinovirus-associated infant respiratory tract illness. J Allergy Clin Immunol. 2011;127:883–91.
Hayden FG, Turner RB. Rhinovirus genetics and virulence: looking for needles in a haystack. Am J Respir Crit Care Med. 2012;186:818–20.
Monto AS, Cavallaro JJ. The Tecumseh study of respiratory illness. IV. Prevalence of rhinovirus serotypes, 1966–1969. Am J Epidemiol. 1972;96:352–60.
Kiang D, Kalra I, Yagi S, Louie JK, Boushey H, Boothby J, et al. Assay for 5′ noncoding region analysis of all human rhinovirus prototype strains. J Clin Microbiol. 2008;46:3736–45.
We thank Arely García Botello for her valuable help in integrating the databases used in this work. We also thank Miguel L. García-León for his help in handling and organizing all pneumonia samples. This work was supported by grants 153639 (to J.I. Santos) and “Influenza 2009” (to C.F. Arias) from the National Council for Science and Technology—Mexico (CONACYT). F.E.A. is recipient of a scholarship from CONACYT.
The authors declare that they have no competing interests.
Conceived and designed the experiments: CFA, RMWC. Performed the experiments: FEA, BT, MAE. Analyzed the data: CFA, FEA, BT, RMWC, SL, JISP. Contributed materials and acquired data: MAAO, JMM, RRV, FDH, FZV, VFR, CNRA, MGB, DEN, LFPG Wrote the paper: CFA, FEA, BT, SL. All authors revised the manuscript critically for intellectual content and gave final approval of the final version.
Fernando E Aponte and Blanca Taboada contributed equally to this work.