1. INTRODUCTION
Dengue fever, dengue hemorrhagic fever and dengue shock syndrome are considered to be the most important arthropod-borne viral diseases due to its high prevalence, and significant morbidity and mortality. The global incidences of dengue has now increased enormously and about 50 - 100 million dengue infections are now reported annually from more than 100 tropical and subtropical countries of the world [1]. Infection with any of the dengue viruses can be manifested in a range of symptoms from a relatively mild flu-like syndrome with rash (dengue fever, DF) to a severe and potentially fatal disease characterized by capillary leakage, thrombocytopenia, and sometimes hypovolemic shock (dengue hemorrhagic fever, DHF/dengue shock syndrome, DSS). Immune enhancement mechanisms following a second heterologous dengue virus infection have been associated with a higher risk of developing DHF/DSS [2], however, cases of DHF/DSS in confirmed primary infections have been observed [3]. Different studies have shown that changes in the viral genetic structure may contribute to shifts in the epidemic potential or pathogenicity of dengue virus serotypes [4,5]. These observations point to the importance of viral factors in determining disease severity. In addition, some strains have been associated with increased pathogenic potential [6].
DENV-3 re-appeared in Latin American in 1994 after an absence of 17 years [7]. Before, serotypes DENV-1 and DENV-2 were circulating for more than a decade in Brazil. In 2002, the introduction of DENV-3 in the state of Bahia produced the first cases of DHF. This event transformed Bahia into a hiperendemic area with three serotypes circulating at the same time [8]. DENV-3, after its introduction, was the serotype responsible for the majority of infections in both Brazil and in Bahia; isolates from different parts of Brazil have been analyzed with respect to their genetic variability. Molecular techniques including sequencing of the viral genomic RNA and direct sequencing of reverse transcriptase/polymerase chain reaction (RT/PCR)-amplified genes have provided many information on geographic distribution, origin and evolution among dengue viruses. Many genes have been used as 5′3′UTR regions, E gene and E/NS1 junction [9,10].
Here we report the analysis of a genomic interval spanning the E/NS1 junction of the dengue genome used for molecular surveillance of this flavivirus and to determine which strains of dengue virus for DENV-2 and DENV-3 were circulating in Bahia state during the 2006- 2007 transmission seasons.
2. MATERIALS AND METHODS
2.1. Sampling
For administrative and epidemiologic purposes, the state is divided into municipalities, which are the basic reporting unit for the State Secretariat of Health (SSH). Dengue is a reportable disease in Brazil. The secretary of health in each municipality is responsible for collecting serum from each suspected case of dengue within five days of onset of symptoms, and sending this on to the Bahia State Central Laboratory (LACEN) for viral isolation and serology. Officially, all cases of dengue in the state should send samples the LACEN-Bahia for viral serotyping. Between 2006 and 2007, 4018 samples were sent to LACEN to be serotyped. 108 samples were positive for isolating in C6/36 Aedes albopictus cells, and 31 positive samples were further analyzed by RT-PCR and sequencing.
2.2. Viral Isolation and Serotyping
Viral isolation was performed in all samples (in 4018 samples) by inoculation of an Aedes albopictus C6/36 cell line with 25 µL of undiluted acute serum. The cell line was then incubated for seven days at 32˚C in Leibovitz’s medium supplemented with 5% fetal calf serum. After fixation to slides in quadruplicate with cold buffered acetone, the cells were incubated with serotypespecific monoclonal antibodies for each of the 4 serotypes followed by fluorescein conjugated anti-mouse antibodies. The slides were then scored as positive or negative by fluorescent microscopy [11].
2.3. RNA Extraction, Amplification and Sequencing
Viral RNA was extracted from serum of infected patients of all positive samples for cell culture using QIAamp Viral RNA Mini kit (Qiagen, CA) following the manufacturer’s instructions. A Mastercycle Eppendorf (model 5331) was used throughout. Five µl of viral RNA solution were added to 45 µl of a medium compatible with both the reverse transcription and PCR amplification steps of the QIAGENOneStep RT-PCR kit. The RTPCR reactions were carried out using an initial reverse transcription step at 41˚C for 45 minutes followed by a denaturation and Hot Star Taq polymerase activation step (94˚C, 15 minutes) and 40 cycles of denaturation (94˚C, 30 seconds), primer annealing (55˚C, 1 minute), and primer extension (72˚C, 30 seconds). A final incubation was carried out at 72˚C for 5 minutes. A second amplification reaction (nested PCR) was seeded with 1 µl of the initial amplification product. The reaction mixture contained 1 buffer B (60 mM Tris-HCl pH 8.5, 2 mM MgCl2, 15 mM (NH4)2SO4, 40 pmol of each sense and antisense primer [9] and 2.5 U of DNA Taq Polymerase (Invitrogen). The samples were subjected to a denaturation step (94˚C, 2 minutes) followed by 40 cycles of denaturation (94˚C, 30 seconds), primer annealing (57˚C, 4 minutes), and primer extension (72˚C, 30 seconds) and a further extension step at 72˚C for 5 minutes.
2.4. Sequence Analysis
Sequencing reactions were performed as recommended in the BigDye Dideoxy Terminator sequencing kit (ABI Prism, Foster City, CA) and the products were analyzed using an automated 3100 DNA sequencer (ABI Prism, Applied Biosystems). Original sequence data were first analyzed by the CHROMAS software (version 1.3, McCarthy 1996; School of Biomolecular and Biomedical Science, Faculty of Science and Technology, Griffith University, Brisbane, Queensland, Australia); forward and reverse sequence data of each sample were aligned using the EDITSEQ program (DNASTAR Inc. Software, Madison, Wisconsin, USA). The consensus sequence was compared and aligned to other samples and DNA database sequences using the software CLUSTAL X, version 1.83 [12]. Program from MEGA package [13] were used to produce phylogenetic trees using neighbor joining as the method to reconstruct the phylogeny and Kimura-2p as nucleotide substitution calculation method. The statistical significance of a particular tree topology was evaluated by bootstrap resampling of the sequences 1000 times. Published sequences used in the comparisons were obtained from the GenBank databases (75 DENV-2 and 73 DENV-3 sequences).
3. RESULTS
3.1. Outbreak Molecular Surveillance
Dengue in Brazil is a notified disease, in 2002 (when DENV-3 was introduced) almost 80,000 cases were notified in all Bahia State (data not showed), and being the first time that DHF cases occurred in this region. After four years, between 2006 and 2007, about 25,000 cases were notified. The most important variable influencing this action is the number of cases [8], this agree that during outbreaks the system became more sensitive to detect disease than in inter-epidemic periods. For example, in 2006 Buritirama, a small city in Bahia with 20,495 inhabitants, between September and December notified about 300 cases of dengue, comparing with Salvador with almost 3,000,000 people that notified for all year just 664 cases. Despite that the surveillance is weak, the outbreaks could be detected, and even in municipalities far from the State Central Laboratory (LACEN) the number of cases remains as the more important aspect to a sample be collected and sent to LACEN.
3.2. Isolation and Typing of the Viruses
Between 2006 and 2007 were analyzed 4018 samples by fluorescent microscopy. Only 108 samples had been the serotype determined. RT-PCR was attempted directly in the serum of the patients of the only 31 viruses isolated in C6/36 cells, the positive samples and the RNA isolation was confirmed. The results indicated that 30 samples were positive for DENV-3, and 1 sample was positive for DENV-2 specific RNA.
3.3. Nucleotide Sequence Analysis
We have determined the nucleotide sequences of the E-NS1 gene junction (220 bp) of 31 dengue viruses in the present study. Pair-wise sequence analysis using Needleman-Wunsch global alignment was carried out on the 220 bp sequence for comparison. As expected, comparisons between serotypes showed a low sequence similarity and could be easily grouped (Figure 1). The results were consistent between distance and characterbased methods, with minimal differences in topologies. Samples from different years, for DENV-3, were grouped in the same genotype (INDIA genotype III) with 2006 samples showing more strong clustering (Figure 2). This genotype has been implicated in DHF outbreaks in American and Asian continents. DENV-2 segregated into South East Asia genotype, as Brazilian samples in the last ten years. In both analyzes the bootstraps values were satisfactory, with higher values showed by DENV- 2 sample.
4. DISCUSSION
The regions, in the dengue virus genome, encoding the envelope protein (E) and the nonstructural protein 1 (NS1) seem to be the most appropriate to characterize genotypes within serotypes [14,15]. In this way, for some flaviviruses, as yellow fever viruses, most molecular markers for pathogenicity have been found in the E gene. However for dengue, several studies of cloned or cell culture-passaged wild and attenuated viruses have been performed, but each report described different sites of nucleotide and amino acid changes possibly associated with virulence [16,17].
DENV-3 introduction was responsible, in 2002, for almost 80,000 cases notified in all Bahia State [8], and being the first time that DHF cases occurred in this region. Between 2006 and 2007, about 25,000 cases were notified. Dengue circulation in Bahia (Brazil) from 2006 to 2007 has consisted mainly of a mixture of DENV serotypes. Between two serotypes (DENV-2 and DENV-3), the former has been consistently detected in outbreaks/ epidemics since its introduction 2002. Before this study, dengue virus genotypes in Bahia state were not available or few samples have been analyzed. In this report, we analyzed 31 local DENV isolates collected in a span of two years to further understand the ecology and dynamics of this virus in the region. As known, the potential for causing severe disease has been described for all four serotypes of DENV, and the main factors considered to explain its pathogenicity are host genetic susceptibility [18], antibody dependent enhancement and differences in virulence among strains [19], in this way determinants of virulence have been located in three genomic regions [5] and have been tested in vitro [20], so the genotypes have been more or less related with potential to cause DHF. Different studies concluded that the spread of genotype III, of DENV-3, from the Indian subcontinent to Africa, and then consequently to Americas, was correlated with an increase in severe cases of dengue disease in this region [21,22].
The routine laboratory diagnosis of dengue virus infection is primarily achieved by the isolation of virus, detection of IgM/IgG antibodies by serodiagnosis and/or molecular detection by the demonstration of viral RNA by RT-PCR [23,24]. RT-PCR is one of the most important confirmatory tests to confirm dengue infection. But, its positivity is highly dependent of the time in the clinical diagnosis and conditions in the sample maintenance (better when is collected within first five days of onset of the symptoms). The lower rate of the virus isolation may be attributed to the absence of live virus in the sample. This must be due to misclassification of the cases by epidemiological staff and others factors as improper maintenance or storage system. All positive cell culture isolations were confirmed and identified by RT-PCR and nested PCR (31 available samples). Still, we have carried out the molecular epidemiology and genotyping study of the detected DENV-2 and DENV-3 viruses, during 2006 and 2007 circulation seasons. We studied the sequence of these dengue viruses, directly from patient serum sample, as recently advocated by several researchers [5,25]. Various genomic regions of dengue viruses have been selected
Figure 1. Phylogenetic tree constructed with the E/NS1 fragment which detected the samples in two different dengue serotypes.
(a) DENV-2,(b) DENV-3
Figure 2. (a), (b)—Maximum likelihood tree depicting the phylogenetic relationships of dengue samples based on the E/NS1 junction (DENV-2 and DENV-3). Numbers on branches represent bootstrap support for each branch. Genotypes are also indicated.
for molecular phylogenetic analysis [26]. However, we have selected the E/NS1 gene junction as it also harbours epidemiologically important sequence information and sensitive for surveillance application [9].
The results obtained for DENV-3 in this study bring some interesting findings how can samples from two years, 2006 and 2007, display a high identity with isolated from Peru in 2002? And in the same time samples from 2007 showed a high similarity with samples differrent countries in all American Continent. Other important finding was that almost all grouping were with samples from the same area in Bahia state. In some way, the microregion where the virus circulate in any level determine how they become grouped. During 2006 and 2007, the main serotype circulating in Bahia state was DENV-3, in agreement with the displacement of serotypes 1 and 2 of DENV after its introduction in Brazil in 2002 [27]. A study developed in Rondônia and Acre found both genotype I and V [10], while we found just genotype III after two years of surveillance. This confirm the epidemic potential of this genotype by the displacement of serotypes 1 and 2 of DENV when a highly virulent genotype of DENV-3 was introduced in Brazil in 2002 [27]. Since its introduction this genotype have been related with severe disease in Bahia and Brazil [8,28].
In conclusion, the effective worldwide control of dengue virus into Bahia and consequently in all Brazil, requires the definition of the epidemical viral sources, and the precise identification of virus genotypes. A key objective of DF and DHF surveillance programmes is early detection of outbreaks to permit the implementation of control measures by the public health authorities. DHF outbreaks can be anticipated by monitoring the emergence of new genotypes in a region like Bahia State.
5. ACKNOWLEDGEMENTS
This project was funded by FAPESB PPSUS, CNPq and FIOCRUZ.
NOTES