INTRODUCTION
Hantavirus cardiopulmonary syndrome (HCPS) is a potentially fatal zoonosis caused by hantaviruses (family Bunyaviridae, genus Hantavirus) described in the Americas. These viruses are associated mainly with rodents of the family Cricetidae, more specifically the subfamilies Sigmodontinae and Neotominae [Reference Musser, Carleton, Wilson and Reeder1, Reference Nichol, Fauquet, Mayo and Maniloff2]. However, with the expansion of studies on wildlife diseases, other animals infected with hantavirus have been identified, such as bats, shrews and moles in Asia, Africa and Europe [Reference Oliveira3]. In South America, the first evidence was described in bats and opossums in the state of São Paulo, Brazil [Reference Araújo4]. More than 1780 HCPS cases have been confirmed in Brazil since 1993, with the second highest incidence in the state of Minas Gerais (258 cases), where the majority of HCPS cases are confirmed by means of serological assays that do not identify a specific aetiological agent [Reference Limongi5, Reference Araujo6].
Several rodent species in the Americas are reservoirs of different hantavirus genotypes, with high levels of specificity among them [Reference Jonsson, Figueiredo and Vapalahti7]. However, many hantaviruses have been identified in more than one host, and additional studies are required to clarify the roles of host specificity, spillover, host-jumping and co-speciation with respect to hantavirus–host association [Reference Oliveira3]. Seven hantavirus genotypes are often associated with HCPS in Brazil: Juquitiba, Araraquara, Paranoa, Laguna Negra, Castelo dos Sonhos, Rio Mamoré and Anajatuba. Three other genotypes, Rio Mearim, Jaborá, and Seoul-like, have been identified in rodents, but an association with human disease has not yet been established [Reference Oliveira3, Reference Jonsson, Figueiredo and Vapalahti7, Reference Oliveira8]. Araraquara virus (ARQV) is a genotype of clade III of Andes virus (ANDV), as is Paranoa virus, and both circulate in central Brazil [Reference Souza and Figueiredo9–Reference Melo-Silva11].
The rodent species Necromys lasiurus is considered the main reservoir of ARQV and is endemic throughout the Cerrado (savanna-like), Caatinga and along the Atlantic Forest–Cerrado ecotone [Reference Oliveira3, Reference Figueiredo12]. ARQV has been identified across the dry northern Cerrado region to the edge of the Atlantic rainforest in Brazil's southeastern region. Studies in this area have demonstrated that ARQV is considered the most virulent genotype and asymptomatic infections are relatively frequent. ARQV has a high fatality rate that may reach up to 50% [Reference Limongi5, Reference Figueiredo12, Reference Figueiredo13].
Juquitiba virus (JUQV), a genotype of clade IV of ANDV that occurs in the Atlantic rainforest and Araucaria subtropical forest environments in Brazil's southern states, was recently identified in the states of Mato Grosso do Sul, associated with Oligoryzomys mattogrossae (specimens found in the Cerrado and Caatinga domains, previously identified as O. fornesi, are here considered as O. mattogrossae) [Reference Weksler, Bonvicino, Patton, Pardiñas and D'Elía14] and Minas Gerais, from HCPS cases, both regions comprising Cerrado environment [Reference Araujo6, Reference Souza and Figueiredo9, Reference Guterres15]. There is a need to clarify which reservoir host is involved in human cases caused by JUQV in the Cerrado, since the main reservoirs of this virus are the rodent species Oligoryzomys nigripes, which is not abundant in this biome, and O. mattogrossae.
Rodents were captured in different phytophysiognomies of the Cerrado (Savanna-like) biome in order to identify the diversity of these mammal fauna, their population dynamics and their relationship with hantavirus. Additionally, an analysis was made of the partial S segment sequences of the hantaviruses obtained from serologically confirmed human HCPS cases and from rodent specimens in an endemic area in the state of Minas Gerais, Brazil.
METHODS
Study area
The study on humans was conducted from December 2011to November 2012 in the HCPS endemic area in the state of Minas Gerais, southeastern Brazil (Fig. 1). This entire area lies within the Cerrado biome, which is characterized by scrub vegetation adapted to long periods of drought, but which is undergoing varying degrees of degradation. The region is characterized by a subtropical climate with two well-defined seasons: a dry winter (May–September) and a rainy summer (October–April). The mean annual temperature and rainfall are 22 °C and 1650 mm, respectively. The rodent study was performed in the municipality of Uberlândia (18° 55′ 07′′ S, 48° 16′ 38′′ W, Greenwich meridian) in regions with different Cerrado physiognomies [Reference Ribeiro, Walter, Sano, Almeida and Ribeiro16].
Samples
Rodents
Rodents were captured and handled according to recommended safety procedures [Reference Mills17]. These mammals were trapped during four campaigns between December 2011 and November 2012, and the trapping effort amounted to a total of 4700 trap-nights. Sherman and Tomahawk traps were used in a 4:1 ratio. Trapping success was calculated as the number of captured animals per trapping effort multiplied by 100, where the trapping effort is the number of traps set multiplied by the number of nights of capture. At the moment of capture, the date and capture site, external measurements, body size, weight, sex, reproductive status and presence of scars of each rodent specimen were recorded. Each animal was anaesthetized and euthanized by intramuscular injection of ketamine hydrochloride (Ketalar® Laboratorio ELEA S.A.C.I.F. y A., Argentina), based on an anaesthetic protocol according to allometric scaling [Reference Pachaly, Brito, Fowler and Cubas18]. Blood was drawn by cardiac puncture from each anaesthetized rodent. Sera and tissue samples of liver, spleen, kidney, lung, and heart were collected and stored immediately in liquid nitrogen for further processing. All the collected specimens were scarified and fixed in 10% formalin, or prepared as skin and skull, and deposited as voucher specimens at the National Museum of the Federal University of Rio de Janeiro. Age structures of the most abundant species were based on body weight (g) and sex, and the individuals were divided into three groups [Reference Mills19]. Age groups of Necromys lasiurus were divided as follows: males [group I (younger) ⩽29 g, group II (adult) 30–50·5 g, group III (older) ⩾50·6 g], females (group I ⩽27 g, group II 28–45 g, group III ⩾46 g).
Most of the species were identified based on an external morphological analysis [Reference Bonvicino, Oliveira and D́andrea20]. Cryptic rodent species taxonomy was confirmed by karyotypic analysis [Reference Bonvicino21]. DNA sequence data of the mitochondrial gene cytochrome b were also analysed to confirm species identification of specimens without cell material for karyotypic analysis, following standard procedures [Reference Sambrook, Fritsch and Maniatis22]. Cytochrome b gene (~1140 bp) was amplified with L14724 [Reference Irwin, Kocher and Wilson23] and MVZ 14 primers [Reference Smith and Patton24] and sequenced with the same primers, in addition to Citb-Sot-in1 and Citb-Sot-in2 [Reference Cassens25] in an ABI PRISM 377 automatic DNA sequencer (Applied Biosystems, USA).
Humans
Serum samples, which were collected from suspected HCPS cases at hospitals after the onset of symptoms, were immediately placed in liquid nitrogen (−196 °C) and then stored in an ultrafreezer (−80 °C). All the samples were tested at the Regional Hantavirus Reference Laboratory of the Oswaldo Cruz Foundation in Rio de Janeiro, RJ, Brazil.
Serological tests
The human serum samples were analysed by ELISA, using an IgM kit from Hantec (Carlos Chagas Institute, ICC/Fiocruz-PR), according to the manufacturer's protocol [Reference Raboni26].
Rodent serum samples were examined for IgG antibodies against the recombinant Araraquara nucleocapsid protein used as the specific antigen by ELISA. The IgG ELISA was performed by coating a microtitre plate with an N-ARA antigen and a negative control antigen overnight at 37 °C. After washing and blocking the plate with 0·1% PBS-Tween 20 and 5% dry milk, serum samples were added. IgG antibodies was detected in rodent sera diluted 1:400 using peroxidase-labelled affinity-purified IgG antibody anti-Peromyscus leucopus and anti-rat conjugate for Rattus rattus. The chromogenic substrate used was 2·2’-azino-di-(3-ethylbenzthiazoline sulfonate), and absorbance was measured at 405 nm. A serum dilution was considered positive if its optical density (OD) was >0·3. The final OD of each dilution of serum sample was calculated as the difference between the OD values measured in wells coated with N-ARQV protein and those coated with negative control antigen. A titre of >1:400 was considered positive [Reference Figueiredo27, Reference Figueiredo28].
Characterization of hantavirus
Rodents
The lung tissues of antibody-reactive rodents were subjected to viral RNA extraction using a PureLink kit (Invitrogen, USA), according to the manufacturer's instructions.
Humans
The serum samples from antibody-reactive patients were subjected to viral RNA extraction using a QIAamp Viral RNA mini kit (Qiagen, USA), according to the manufacturer's instructions.
PCR
One-step RT–PCR was performed, followed by semi-nested PCR. Amplification and sequencing were performed using specific primers for the nucleoprotein coding region of the S segment of the viral genome [Reference Guterres29]. Amplicons of expected sizes (479 bp) were purified from agarose gel, using a Wizard® Genomic DNA Purification kit (Promega, USA). The PCR product purified from the gel was sequenced using a commercial BigDye Terminator v. 3.1 cycle sequencing kit (Applied Biosystems). The nucleotide sequences were obtained in an ABI PRISM 3130xl automated sequencer (Applied Biosystems). The nucleotide sequences and their sequencing electropherograms of the amplified fragments of genomic S segment were analysed using the MEGA 5.0 program [Reference Tamura30]. A consensus sequence was established, and the divergences between the nucleotide sequences were clarified by analysing the sequencing electropherograms. This sequence was compared with a reference panel of sequences that comprises most of the American hantaviruses, for which there are at least 390 nt sequences available in GenBank. Multiple sequence alignment and comparison of nucleotides were performed using muscle, in the SeaView v. 4 software program (http://pbil.univ-lyon1.fr/software/seaview). The phylogenetic relationships were estimated by the Markov Chain Monte Carlo (MCMC) method implemented in MrBayes v. 3·1·2 using the GTR + G nucleotide substitution model [Reference Ronquist and Huelsenbeck31]. The MCMC settings consisted of two simultaneous independent runs with four chains each, which were executed for 10 million generations and sampled every 100 generations, yielding 100 000 trees. After eliminating 25% of the samples as burn-in, a consensus tree was built. Statistical support of the clades was measured by Bayesian posterior probabilities. Hantaan virus (NC005218) and Seoul virus (AY027040) sequences were used as outgroups for the analyses.
Nucleotide similarity among hantavirus partial S segment amplimers was determined, along with other GenBank sequences representing different hantavirus strains, using MegAlign v. 5.05 (DNASTAR Inc., USA).
Statistical analysis
The χ 2 with Yates correction and Fisher's exact tests were applied for comparison between proportions, using Epi Info v. 3·5·4 (CDC, USA). Statistical significance was set at 5% (P⩽0·05).
Ethical standards
The authors state that all the procedures involved in this work are in compliance with the ethical standards of the pertinent national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008, and with institutional guidelines on the care and use of laboratory animals. The project was approved by the Ethics Committee for Human Research of the Oswaldo Cruz Institute in Rio de Janeiro. The patients' written consent was obtained. Permits for field collection were provided by IBAMA (Brazilian Institute of Environment and Renewable Natural Resources). Animal handling procedures were approved by the Animal Research Ethics Committee of the Federal University of Uberlândia.
RESULTS
Rodents
Two hundred and fifty-eight rodents were captured, 245 (95%) in Sherman traps and 13 (5%) in Tomahawk traps. Eleven taxa were identified to the species level and two to the generic level, all belonging to the subfamily Sigmodontinae (Tables 1 and 2). The most numerous species captured were N. lasiurus (181, 70·1%) and Calomys expulsus (30, 11·6%). The phytophysiognomy of semi-deciduous dry forest was the most biologically diverse environment for Sigmodontinae. The total number of captured animals was directly related to climate condition, and this number was significantly higher in the dry season than in the wet season (P < 0·0001) (Tables 1 and 2).
Nl, Necromys lasiurus.
A larger number of males (133, 51·6%) than females were captured. There was a higher prevalence of pregnancy in rainy summer periods than in the dry winter season (P < 0·0001) (Table 1).
N. lasiurus was the main rodent found in all types of vegetation (min-max: 50–83·7%). C. expulsus was found in four types of environment, but was absent from Brachiaria spp. (Table 2).
There was a predominance of individuals of age group III (older) during the rainy season (60%). During the dry season, the number of individuals of age groups I and II corresponded to 85·5%, indicating a higher percentage of young individuals in the population during this period.
Four rodents were seroreactive (overall prevalence 1·6%). Only N. lasiurus had IgG antibodies to hantavirus (N. lasiurus prevalence 2·2%) (Table 1). Three (75%) of the seroreactive animals were males and all seroreactive animals presented scars. Scarring was associated with the presence of antibodies to hantavirus (P = 0·01, Fisher's exact test). Fragments of S hantavirus genome segment were detected in three rodents (Supplementary Fig. S1). Two sequences from the rodents [LBCE 18 197 (Genbank: KT257178) and LBCE 17 849 (Genbank: KT257177)] were obtained and subjected to phylogenetic analysis.
Humans
All patients suffering from HCPS were male, mean age 34·6 ± 14·0 years. The fatality rate of the cases was 60% (Table 3). The classic symptoms of fever, headache, dyspnoea, myalgia, and dizziness were observed in all cases. The time between the onset of symptoms and hospital admission was 4·8 ± 2·1 days.
* All the patients were male and were reactive in ELISA and positive by PCR. Genotype ARQV was found in all the patients.
With regard to risk activities, 85·7% of the patients had engaged in occupational or occasional farming activities up to 60 days before the onset of symptoms.
Human serum samples were collected on average 4 days after the onset of symptoms of six suspected HCPS cases. All were serologically confirmed HCPS cases and five showed amplified PCR products of expected size (Table 3).
Phylogenetic analysis
Phylogenetic analysis of nucleotide sequences of human and rodent samples showed only the circulation of the ARQV genotype (Fig. 2).
The sequenced rodent samples, from rural areas in the municipality of Uberlândia, exhibited a high phylogenetic identity with human sample LH_206_12 (Genbank: KT257181), with a well-supported cluster [posterior probability (pp) = 0·98]. Samples of rodents, with human samples [LH_206_12 and LH_180_12 (Genbank: KT257179)], were more closely related to the Paranoa virus variant described in the federal district, located at a distance of about 463 km from the municipalities of Uberlândia and Patrocínio, than to the other cluster composed of three human samples of this study, located at a distance of about 261 km from these two municipalities. Pairwise sequence comparisons showed 91·6–93·7% nucleotide and 97·6–100% amino-acid identity between the samples of this study and Paranoa virus.
The LH-Nat (Genbank: KT257183), LH_212_12 (Genbank: KT257182) and LH_197_12 (Genbank: KT257180) samples were phylogenetically more closely related, forming a well-defined cluster in the three samples, with a well-supported cluster of pp = 1. The first two came from patients in the municipalities of Rio Paranaíba, MG and São Gotardo, MG, respectively, located only 30 km from each other. However, the LH_197_12 sample came from a patient in São Tomás de Aquino, a municipality located 384 km from the municipalities where the other samples occurred. Pairwise comparisons of these sequences showed 95·8–98·2% nucleotide and 100% amino-acid identity.
According to a similarity analysis, the greatest divergence in our sequences was 12·2% (nucleotide) and 2·4% (amino acid) between LH_Nat (Rio Paranaíba) and the rodent sample LBCE_17849, and the lowest was 1·6% (nucleotide) and 1·6% (amino acid) between rodent samples.
DISCUSSION
Climate is an important abiotic factor that influences the availability of food for rodents. The greater supply of food in the rainy season is directly associated with a higher reproductive rate of rodents. In this study, females showed a higher rate of pregnancy in the rainy season. This relationship leads to a higher proportion of juveniles in the dry season, resulting in an increased number of captures, and of adults in the rainy season, which have a higher prevalence of hantavirus. The population dynamics of these rodents is consistent with the results obtained in other studies [Reference Teixeira32, Reference Limongi33].
Eleven taxa were identified to the species level and two to the generic level, all belonging to the subfamily Sigmodontinae. This information reflects the territory in which the rodents were captured, in wild environments far away from peri-urban areas.
N. lasiurus was the most prevalent (70·1%) in the captured species. These rodents are aggressive and generally dominate the area that they infest [Reference Streilen34]. The general behaviour of N. lasiurus makes this species highly adaptable to man-made environments and to less complex ecosystems, especially grassy areas (Brachiaria spp.). N. lasiurus specimens represented >30% of the individuals captured in most of the areas sampled in the Cerrado biome [Reference Marinho-Filho35].
In our study, N. lasiurus represented at least 50% of the captured specimens (min-max: 50–83·7%). The higher prevalence of N. lasiurus in this study compared to earlier reports can be directly attributed to the intense process of degradation this biome has undergone in recent decades, with large areas planted with monoculture crops, which have favoured this species in interspecific competition with other rodents, mostly habitat specialists [Reference Ribeiro, Walter, Sano, Almeida and Ribeiro16, Reference Streilen34]. The ecological pattern of N. lasiurus associated with deep anthropogenic changes in this biome greatly increases the chances of contact between humans and infected rodents, mainly due to the potential increase in population density of these rodents and the narrowing of environments shared with humans [Reference Limongi36].
The rodent O. mattogrossae was captured in three phytophysiognomies. This rodent is associated with the hantavirus Anajatuba in northeastern Brazil and has also recently been considered a reservoir of JUQV in the Cerrado biome. In southern Brazil and in the state of São Paulo, the main reservoir of this virus is the species O. nigripes, which was also found in this study [Reference Guterres15, Reference Rosa37, Reference Suzuki38]. The species of the genus Calomys found in this study are not associated with any hantavirus, despite a report of ARQV that was sequenced from a lung sample from a C. tener specimen in the state of São Paulo. However, the authors attributed the infection to a likely spillover of the virus of N. lasiurus [Reference Figueiredo39].
The seroreactive rodents had scars on their ears and/or tail. The presence of scars in rodents may be associated with seroreactivity to hantavirus in reservoirs, and studies have shown that the prevalence of infection is higher in adult males than in females and juveniles [Reference Teixeira32, Reference Limongi33].
Since the first record of a case of HCPS in the state of Minas Gerais in 1998, only three molecular studies in Brazil have included human samples of HCPS cases from these regions, with 15 samples analysed phylogenetically, collected in the municipalities of Araxá (1), Patrocínio (4), Uberlândia (2), Uberaba (4), Patos de Minas (1), São Tomás de Aquino (1), São Sebastião do Paraíso (1) and Passos (1). In these studies, the cases of HCPS were associated primarily with ARQV [Reference Araujo6, Reference Figueiredo12, Reference Suzuki38]. However, evidence of the circulation of JUQV, as well as a highly divergent strain of ARQV and JUQV, has been reported in humans and rodents, respectively, in these areas [Reference Araujo6, Reference Figueiredo12]. In one of these studies, a virus related to JUQV was reported in human serum samples of from the municipalities of Uberaba, Patrocínio, São Sebastião do Paraíso and Passos [Reference Figueiredo12]. In our study, we found only the ARQV genotype in the samples from the municipalities of Patrocínio and São Tomás de Aquino, a municipality located 23 km from the municipality of São Sebastião do Paraíso. In the same study in which JUQV was reported in the state of Minas Gerais, ARQV was also reported in the municipalities of Passos, Teófilo Otoni, Uberlândia, Uberaba, Patos de Minas, São Tomás de Aquino, Sete Lagoas and Piumhi. These data corroborate the findings of our study in Uberlândia, where only the ARQV genotype was found.
JUQV is one of the most widely distributed hantaviruses found in South America and has been detected in O. mattogrossae, O. nigripes, Akodon montensis, Oxymycterus judex, Akodon paranaensis, O. nigripes, Oryzomys sp. and O. fornesi rodents in Brazil, Argentina, Paraguay and Uruguay [Reference Guterres15] JUQV is often related to HCPS cases occurring in southern Brazil and in the state of São Paulo [Reference Jonsson, Figueiredo and Vapalahti7, Reference Figueiredo12, Reference Suzuki38]. However, the rodents O. mattogrossae and O. nigripes also occur in the Cerrado biome in the state of Minas Gerais [Reference Bonvicino, Oliveira and D́andrea20]. Five specimens of O. nigripes were found in a study of hantavirus reservoirs in the municipality of Uberlândia, but they were not hantavirus seropositive [Reference Limongi33]. In the present study, one O. nigripes specimen and six O. mattogrossae specimens were captured, all of them seronegative for hantavirus. It should be noted that the JUQV genome was recently sequenced from specimens of O. mattogrossae from the Cerrado biome in the state of Mato Grosso do Sul, confirming the complexity of the hantavirus–wild rodent interaction [Reference Guterres15].
A highly divergent strain of the ARQV and JUQV genotypes was recently identified in N. lasiurus in the municipality São Gotardo, MG, and the authors sequenced ARQV genotype in two other N. lasiurus specimens captured in the same municipality [Reference Figueiredo12]. In our study, only the ARQV genotype was identified in a human sample from this municipality, as well as in one from Rio Paranaíba, MG, located 30 km from São Gotardo, MG.
In a phylogenetic study of hantavirus genotypes in samples from southern and southeastern Brazil, the viral sequence of a patient with HCPS from the municipality of Patrocínio, MG showed the highest divergence in the members of the Araraquara clade (16·7%) [Reference Suzuki38]. In this same study, but in the municipality of Uberlândia, another viral sequence extracted from the lungs of a N. lasiurus rodent also proved highly divergent from the other genotypes of the Andes clade. In the phylogenetic analysis of our study, the LH_180_12 sample of the patient from Patrocínio who died, along with the LH_206_12 sample of a patient with HCPS in Uberlândia, were also the most divergent between clade samples. These data indicate that genetic diversity of ARQV is greater than previously believed.
The samples analysed in this study came from an area whose natural ecosystems are in large part degraded as a result of intensive farming and cattle ranching. The environmental and demographic changes that have taken place in recent decades have probably affected the ecology of wild rodents, facilitating the occurrence of human hantavirus infections and the emergence of HCPS in this region, particularly ARQV [Reference Figueiredo12]. Samples identified with ARQV analysed in this study were distributed ~400 km from one another.
Higher case-fatality rates are observed in regions where the ARQV genotype is endemic in Brazil than in areas where JUQV occurs [Reference Figueiredo12]. In our study, the case-fatality rate of 60% of HCPS cases provides evidence that ARQV is the most prevalent genotype in the study area. Although other studies are necessary to determined which factors influence the outcome of HCPS cases caused by ARQV or JUQV, the occurrence of different genotypes/viral variants may partly explain the clinical spectrum of hantavirus infection, which ranges from asymptomatic infections to cases with high fatality rates [Reference Limongi5, Reference Limongi36].
This study reported the circulation of a single genotype, i.e. ARQV, in the most endemic area of the state of Minas Gerais, Brazil, where N. lasiurus is the predominant species. However, it is reasonable to assume that other reservoir hosts of hantavirus can be associated with human disease in this region, due to the presence of JUQV associated with HCPS cases and Oligoryzomys species. The environmental changes that have occurred in the Cerrado biome in recent decades have favoured N. lasiurus in interspecific competition of habitats, which, in turn, has considerably increased the risk of contact of these rodents with humans. All our data corroborate the definition of N. lasiurus as the main hantavirus reservoir in the Cerrado biome [Reference Suzuki38].
SUPPLEMENTARY MATERIAL
For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S095026881500237X.
ACKNOWLEDGEMENTS
The authors thank Matias Pablo Juan Szabó for assistance and technical support.
This work was supported by CNPq Papes FIOCRUZ V 403601/2008-8 and Papes FIOCRUZ VI 407664/2012-2 (Brazil).
DECLARATION OF INTEREST
None.