The Atlantic Forest is one of the most disturbed Brazilian biomes, with 183 out of 298 species of mammals occurring in the state of Rio de Janeiro. In this study, we aimed to characterize the diversity, community structure, and habitat use of small mammals in the FIOCRUZ Atlantic Forest Campus (CFMA), including areas of Pedra Branca State Park (PBSP, subunit Pau da Fome), state of Rio de Janeiro. We also compared species diversity and composition between two moments 15 years apart (2001 and 2012–2015) and analyzed the population dynamics of the marsupial Didelphis aurita (Wied-Neuwied, 1826). Small mammal captures were made in different habitats: sylvatic-urban interface areas near human dwellings, disturbed forest, and preserved forest areas. Five marsupial species and four rodent species were captured in both periods. There was a reduction in species richness and β diversity between the two periods, indicating that disturbances in the environment over the years may have affected the small mammal community structure. The most altered environment showed the greatest species richness and abundance, while the forest areas showed the smallest values, which may be explained by the loss of mammal species, mainly specialist species in forested areas. We identified three groups of species according to habitat preferences: one related to environments with a higher density of vegetation in upper strata – Marmosa paraguayana (Tate, 1931) and Monodelphis americana (Müller, 1776), another related to a higher density in lower forest strata – Akodon cursor (Winge, 1887), and another with no association with the investigated habitat variables – D. aurita and Oligoryzomys nigripes (Olfers, 1818). The small mammal community structure showed a low level of nestedness in both sampling periods. This study is the first report to evaluate the community structure of small mammals in the sylvatic-urban interface area of Pedra Branca State Park, the largest forest reserve within an urban area in Brazil. The surveys indicate that the small mammal diversity was low in both sampling periods and in both areas, and a species loss in the Pau da Fome locality was observed, despite it is a conservation unit. The greater species abundance and richness in the most disturbed areas suggest an increase of factors favoring the occurrence of synanthropic and opportunistic species.

Atlantic Forest, diversity, habitat disturbance, marsupials, rodents

The continuous process of urbanization and expansion of human activities into forest areas leads to changes in the natural habitats, altering species diversity and abundance patterns (Soulé and Wilcox 1980, Ewers and Didham 2006, Devictor et al. 2008). In disturbed ecosystems, generalist and opportunistic species can be favored, increasing their densities and dispersion to rural or peridomicile areas. These patterns may have important consequences for conservation (Umetsu and Pardini 2007) by excluding more specialist species from the communities; for economics by increasing the number of agricultural pests (Stenseth et al. 2003); and for public health by changing the natural prevalence of zoonotic pathogens in wild reservoirs (Thompson 2013).

The fragmentation and degradation of natural habitats have been occurring quite rapidly for decades in tropical regions. This process has several effects on the demography, genetic structure of populations, and geographic distribution of a species, leading to changes in its community structure. As a result, the process of species extinction may be accelerated, biodiversity may be reduced, or community species composition may be altered by the introduction of invasive species. Local extinction would occur not only by the total elimination of a species’ habitat, but, mainly, by the accumulated effects on small populations.

The Atlantic Forest is considered one of the hotspots of biodiversity on the planet (Myers et al. 2000) due to its species richness, number of endemic species, and degree of threat (Mittermeier et al. 1998). It is one of the most disturbed Brazilian tropical biomes, with the largest number of endangered species by area until the most updated available information (Bergallo et al. 2000, Myers et al. 2000). Habitat fragmentation and degradation related to human activities have occurred for hundreds of years, reducing the biome’s original area by about 93% (Fundação SOS Mata Atlântica and INPE 2016). In the state of Rio de Janeiro, the removal of forested areas reached approximately 81% of the original cover, with only about 814,563 ha remaining (Fundação SOS Mata Atlântica and INPE 2016). The state of Rio de Janeiro constitutes a strategic portion along the continuum of remnants of the Atlantic Forest because it incorporates high species richness and endemism for different groups of animals (Rocha et al. 2003).

Concerning the mammals, 298 species are recognized in the Atlantic Forest (Fundação SOS Mata Atlântica and INPE 2016) and 42 are considered as threatened (Chiarello et al. 2008), with 183 occurring in the state of Rio de Janeiro, at least three of which are endemic to the state (Rocha et al. 2003). These data indicate that this region is an area of ​​high richness of mammal species. Medium and large mammals are extremely threatened by hunting and habitat disturbance (Fonseca et al. 1994, Cullen et al. 2000, 2001) because of their position in the trophic chain and their widespread habitat needs. In some localities of the Atlantic Forest, populations of these species disappeared completely, leading to changes in ecological interactions and likely increasing the abundance of generalist species of small mammals, such as the common opossum Didelphis aurita (Wied-Neuwied, 1826) (Fonseca and Robinson 1990, Gentile et al. 2004).

These factors make the Atlantic Forest remnants extremely important for biological conservation in the State of Rio de Janeiro. In this study, we aimed to characterize the diversity and community structure of small mammals, also evaluating their habitat use, in a sylvatic-urban interface area of the FIOCRUZ Atlantic Forest Campus (CFMA) of the Oswaldo Cruz Foundation (Campus FIOCRUZ Mata Atlântica – CFMA) and its surroundings, including areas of Pedra Branca State Park (PBSP), state of Rio de Janeiro. Herein, we also compared species diversity and composition between two moments, 15 years apart (2001 and 2012–2015), and analyzed the population dynamics of D. aurita, the most abundant and constant species throughout the study. To understand the community structure of small mammals and how the species are related to different habitat types, we addressed five inter-related questions: (1) Does diversity or species composition vary among habitat types? (2) What is the pattern of the community structure of the small mammals? (3) Did diversity, species composition and community structure varied between 2001 and 2012–2015 samplings? (4) Are small mammal species more related to a specific area or habitat variable? (5) Is there any influence of habitat heterogeneity on the mammal assemblage structure?

The study was carried out in areas of the CFMA and surroundings, encompassing areas of the FIOCRUZ Atlantic Forest Biological Station (EFMA) and PBSP in the municipality of Rio de Janeiro. The CFMA is part of the buffer zone of PBSP, including preserved areas of the park, and is also adjacent to an expanding urban region in the western area of Rio de Janeiro. The area that today comprises the CFMA has been previously used for sugarcane mills and coffee farms, with construction of large residences and roads. In the 20^th century, a psychiatric colony with various urban facilities was established. In 2003, the area was incorporated to FIOCRUZ, and from then on, non-edification areas were defined and a more effective process of environmental protection and recovery was established in order to mitigate the effects of anthropic impacts. Thus, the 2001 samplings preceded this area management and conservation program, and the 2012–2015 samplings were carried out a decade after the implementation of this program. PBSP is the largest forested reserve located within an urban environment in Brazil. It covers an area of 12,492 ha, and the predominant phytophysiognomy is Ombrophilous Dense Atlantic Forest vegetation. PBSP was created in 1974. Although the region was also occupied by sugarcane and coffee cycles in the colonial period, part of this forest was preserved in 1908 as a protection for the hydric resources. In the last decade, with the population growth of the surrounding area, expansion of low-sanitation communities, implementation of major structuring works and increase of the road network, the area began to suffer great anthropic pressure from the surrounding area. These factors have impacted on the quality and structure of the forest borders and boundary areas of the PBSP, including the areas of our study. Currently, those border areas have a disturbed forest with little understory and sparse midstory.

Small mammal captures were made in transects inside the CFMA and PBSP (subunit Pau da Fome), in areas with different habitat characteristics: (1) sylvatic-urban interface areas near human dwellings (peridomicile) in the CFMA and Pau da Fome (22°56’18”S 43°24’11”W and 22°55’57”S 43°26’34”W), (2) disturbed forest areas in the CFMA and Pau da Fome (22°56’28”S, 43°24’34”W; 22°56’45.5”S 43°24’45.6”W and 22°55’57”S, 43°26’36”W), and (3) preserved forest areas in the CFMA (22°56’47”S, 43°25’07”W and 22°56’41”S, 43°24’52”W) (Fig. 1).

Figure 1. 

Map of the study area indicating the sampling localities in CFMA and Pau da Fome region, State of Rio de Janeiro, Brazil. 1) CFMA – Peridomicile, 2–3) CFMA – Disturbed Forest, 4) CFMA – Preserved Forest, 5) Pau da Fome – Peridomicile, 6) Pau da Fome – Disturbed Forest.

The preserved forest areas presented a vegetation canopy ranging from 10 to 40 m and an irregular slope varying from flat to steep. The disturbed forest areas were in the regeneration process, with vegetation in different ecological succession stages located between areas of human occupation and preserved areas. They had a medium vegetation height ranging from 6 to 20 m and a flat slope. The peridomicile areas were located in the backyards of houses and had a lower canopy than the other areas, ranging from 6 to 15 m, a predominance of bushes and small trees, some flooded areas, and a flat to moderate slope. In all areas, the understory was sparse to semi-open in most transects. The climate of the region is humid mesothermal, with hot, rainy summers and mild winters.

The study was done in two periods: in 2001, with two samplings in each area; and between 2012 and 2015, with two samplings in Pau da Fome and six in the CFMA in all habitat types. Small mammal captures were carried out in 20-point transects established in different habitat types. Capture points were equidistant, 15 m apart, with two traps per point on the ground: a Sherman (7.62 cm × 9.53 cm × 30.48 cm) and a Tomahawk (40.64 × 12.70 × 12.70 cm), both suitable for the capture of live animals up to 3 kg. Each trapping session lasted five consecutive nights, except for one sampling in 2001 in the CFMA, which lasted 10 nights. The bait used was composed of a mixture of bacon, oats, bananas, and peanut butter. Trapping effort for each area and sampling period are better described in Table 1.

The captured animals were anesthetized, euthanized, and had their bionomic data collected (age, sex, reproductive activity, body mass, and external measurements, such as body, tail, foot and ear lenghts) in a laboratory field base. Young animals, lactating females and specimens exceeding the limit of the capture license (only D. aurita) were marked with ear-tags and released at their trapping points. All the other animals were taxidermized and deposited as voucher specimens in the scientific collection of the Department of Vertebrates of the National Museum of Rio de Janeiro. Specimens were identified by external and cranial morphology. Rodent species of the genus Akodon Meyen, 1833 were identified by their diploid number after karyotyping.

Animals were captured under authorization of the Brazilian Government’s Chico Mendes Institute for Biodiversity and Conservation (ICMBIO, license number 13373) and the Environmental Institute of Rio de Janeiro State (INEA, license number 020/2011). All procedures followed the guidelines for the capture, handling, and care of animals of the Ethical Committee on Animal Use of the Oswaldo Cruz Foundation (CEUA, licenses number P-0083-01, L-049/08, LW81/12, and LW-39/14). Biosafety procedures and personal safety equipment were used during all procedures involving animal handling and biological sampling.

Microhabitat variables were measured at each trapping point, considering an area of 9 m² around the traps only for the 2012–2015 period, and only for the CFMA transects. The variables were measured at two points, 1.5 m apart, within each trap station according to a modified version of Freitas et al. (2002) method. Nine quantitative variables were measured: percentage of vegetation cover on the soil (VCS), percentage of litter (LITT), percentage of exposed soil (SOIL), percentage of rocks (ROCK), percentage of canopy cover (CANO), canopy height (CAHE), number of trees with diameter at breast height ≥ 5 (TREE), number of fallen trunks (TRUN), and vertical vegetation obstruction at 50 cm to 1 m (VVO). The variables were measured using a wire square of 0.5 × 0.5 m divided into 100 smaller squares, as proposed by Freitas et al. (2002). We quantified the frequency of squares occupied by each variable, except for CAHE, which was measured in meters, and TREE and TRUN, which were directly counted in numbers.

Population dynamics were analyzed in the 2012–2015 period only for the opossum D. aurita, which was the most abundant species. The sex ratio was tested by chi-square test. The reproductive activity was analyzed, considering the presence of lactating females. The age structure analysis was performed by dividing the animals into age classes according to dental development classes and reproductive condition according to Gentile et al. (1995).

Trapping success was calculated as the number of captured animals divided by the number of trap-nights multiplied by 100. The species richness was considered as the number of species observed during the study. In addition, a species richness index was estimated using the Jackknife 1 method (Palmer 1990). Estimators were calculated for the 2001 and 2012–2015 data using the “vegan” package (Oksanen et al. 2016) of RStudio software, version 0.99.902. The relationship between abundance and species richness was investigated for the two periods using Spearman’s correlation. In this case, abundance values were standardized by considering the total trapping effort of each transect.

Species diversity between areas and between the 2001 and 2012–2015 samplings were compared using Nonmetric Multidimensional Scaling (NMDS), with the “stats” package of RStudio software, version 0.99.902 (R Core Team 2017). The Gower distance measure was used for the relative abundances of the species, considering a stress value of up to 0.2. Only for this analysis, relative abundances were calculated by dividing the abundance value of each species by the abundance of the most abundant species for each case based on Williams et al. (2002). We also calculated the total Whitaker β diversity and the β diversity between each transect pair for 2001 and 2012–2015, separately. These analyses were done in Past software, version 3.08 (Hammer et al. 2001).

The habitat variables were transformed into a standard area unit according to Cerqueira and Freitas (1999). Subsequently, habitat heterogeneity indexes for each CFMA transect were calculated according to Paglia et al. (1995). The variance of each habitat variable per transect was divided by the total variance of that variable, considering all transects. Then, these quotients were added for all the variables in each transect in order to obtain a heterogeneity index. The relationship between habitat heterogeneity and species richness was investigated using Spearman’s correlation.

In order to investigate the relationship among the abundances of small mammals, the habitat variables, and the areas, abundance data were ordered along the canonical axes of a Canonical Correspondence Analysis (CCA), according to its ecological optima (Borcard et al. 2011). The analysis was performed using the means of each habitat variable for each transect. Before ordering, correlated and multimodal variables were excluded. The CCA was made with the following variables: VCS, VVO, CANO, and TREE. The rodent Rattus rattus (Linnaeus, 1758) was not included in this analysis. The correlation analysis was done with the “Hmisc” package (Harrell 2016) and the CCA with “Vegan” package (Oksanen et al. 2016) of RStudio software, version 0.99.902, only for the CFMA-2012–2015 data.

The mammalian community structure was investigated for both the 2001 and 2012–2015 data, separately, using the abundance data standardized for trapping effort. The structure was analyzed based on the nesting measure using the weighted “NODF” method (Almeida-Neto and Ulrich 2011) in the “Bipartite” package (Dormann et al. 2008) of RStudio software, version 0.99.902.

The small mammal species richness observed was similar between 2001 (8) and 2012–2015 (7), although the estimated richness (10) was slightly higher for the 2001 period. However, there was a reduction in species richness in Pau da Fome due to the absence of captures of the marsupials G. microtarsus, M. parguayana, and M. nudicaudatus in the recent samplings, as well as a reduction of β diversity between the two periods. The lack of correlation between species abundance and richness for the most recent period may also indicate that disturbances in the environment are continuously occurring in spite of both areas are under conservation measures, affecting the small mammal community structure, as discussed hereafter.

Another small mammal survey, also carried out in several altitudes of PBSP by Oliveira et al. (2012), registered the occurrence of 10 species of wild small mammals from 2005 to 2007, with a similar dominance of D. aurita. The species richness was similar to our study; however, the species composition showed a few differences. Oliveira et al. (2012) reported Marmosops incanus (Lund, 1840), Philander frenatus (Olfers, 1818), and Guerlinguetus brasiliensis (Gmelin, 1788), which were not captured in the present study, although the last one was visually recorded in 2001. Moreover, G. microtarsus was captured in the present study, but not in Oliveira et al. (2012). Considering all the information gathered, we may consider a total of 11 wild non-volant small mammals observed in PBSP so far, excluding the two invasive Rattus species. When comparing the species richness in this study with other studies of small mammals in the Atlantic Forest, some of which, carried out in continuous forest areas, they related more than 20 species (23 in Pardini and Umetsu 2006, 26 in Machado et al. 2013, 21 in Maestri et al. 2014), meanwhile we observed a very low species richness, especially of rodents.

The results of the NMDS analyses indicated a greater difference in species diversity when compared the 2001 with the 2012–2015 sampling periods than among localities. Comparing the localities within each period, the greater discrepancies were observed between peridomicile areas in the CFMA and the other transects, which was expected since the peridomicile areas represent a different environment with more distinct characteristics when compared to the preserved and disturbed forest environments. However, the peridomicile areas, which are the most altered environments, showed the greatest species richness and abundance, while the forest areas showed the smallest values. We suggest two hypotheses. The direct or indirect impact caused by anthropic activities in the studied forested areas could result in a process of loss of mammal species, mainly the specialist ones, leading to a phenomenon known as empty forest (Redford 1992). Another possibility is that the occasional captures of some species in the forest environment is due to the limits of the dispersion of the individuals coming from more preserved areas in the interior of the forest. The latter hypothesis could explain the only capture of M. nudicaudatus during the study, because this species may roam large areas inside their habitats (Gentile and Cerqueira 1995).

The greater abundance and richness of species in the most disturbed areas suggest an increase of factors favoring the occurrence of generalist and synanthropic species, such as D. aurita and the introduced rodent R. rattus, which are frequently found in areas close to human dwellings, in addition to the opportunistic rodents A. cursor and O. nigripes. The highest trapping success observed for the most disturbed areas may also be associated with the conditions present in these areas, which favored the greater abundance and occurrence of these species. The occurrence of D. aurita and O. nigripes in the three types of environment indicates an opportunistic and plastic character in relation to the present environmental conditions of the areas. These results highlights the importance of a continuous monitoring in the areas, since these species are widely known as reservoir of zoonotic pathogens in the state of Rio de Janeiro (Vaz et al. 2007, Oliveira et al. 2009, 2017).

The marsupial D. aurita showed much greater abundance compared to the other species and was present throughout the study. Its greater abundance in the peridomicile areas supports other studies that have reported its presence as a dominant species of small mammal communities in disturbed areas (Cáceres and Monteiro-Filho 1998, Gentile et al. 2000). Other studies conducted in sylvatic-urban interface areas around Serra dos Órgãos National Park also recorded the occurrence of this marsupial in high abundance (Gentile et al. 2004, Kajin et al. 2008).

The population dynamics of the common opossum in this study were similar to those in other studies carried out in Atlantic Forest areas, where the greatest abundances of this marsupial were usually recorded between May and August, at the end of the rainy season, and during the dry season (Cáceres and Monteiro-Filho 1997, Gentile et al. 2000, Graipel and Santos-Filho 2006, Bonecker et al. 2009). This pattern could be related to the capture of young and subadults born in the previous reproductive season. In other studies of population dynamics of D. aurita in the state of Rio de Janeiro, the opossum also presented a pattern of seasonal reproduction (Cerqueira 2005), with a reproductive period usually beginning in July to August and ending in March to April (Gentile et al. 1995, Macedo et al. 2007, Kajin et al. 2008).

The habitat structure based on its heterogeneity did not show a relationship with the species richness, indicating that the degree of habitat heterogeneity may have varied within each environment. However, this analysis allowed us to identify three groups of species: one related to environments with a higher density of vegetation in upper strata, another related to a higher density in lower forest strata, and another with no association with the investigated habitat variables, as discussed bellow.

Although the occurrence of D. aurita is associated with the presence of certain variables in the environment (Gentile and Fernandez 1999, Gentile et al. 2004, Freitas et al. 1997), no association of this species was observed with the habitat variables in the present study. The higher abundance of the rodents O. nigripes and A. cursor in areas with greater anthropic influence (peridomicile), also confirms the opportunistic nature of these species in relation to habitat occupancy, as previously observed in disturbed regions and along forest borders (Olifiers et al. 2005, Vaz et al. 2007, Geise 2012, Cardoso et al. 2016, Gonçalves et al. 2016). However, A. cursor showed more association with habitat variables, such as a close vegetation near the ground, which can provide shelter, protection against predators, and food resources (Dalmagro and Vieira 2005, Lima et al. 2010, Melo et al. 2013), a higher vegetation density in the understory and an open canopy, which are characteristic of open areas (Cirelli and Penteado-Dias 2003).

The marsupial M. nudicaudatus was one of the least abundant species in the study area, although it was recorded in both sampling periods. Other studies also reported a low abundance of this species in the Atlantic Forest (Pardini 2004, Caldara-Junior and Leite 2007, Santos-Filho et al. 2012). The very low abundance of this species may have contributed to the lack of association between this species and the environmental variables. However, an association of M. nudicaudatus with areas of closed understory (Freitas et al. 1997), lower canopy cover, and a moderate degree of disturbance has already been reported (Moura et al. 2005).

Marmosa paraguayana, which is a scansorial marsupial (Paglia et al. 2012) and may occur in forest environments (Prevedello et al. 2008), fragments and is able to cross matrices of grassland (Pires et al. 2002), was also found in low abundance in the preserved forest areas, with only one individual in the peridomicile. The relationship between this species and environments with higher vegetation density in the upper strata was also found by Kajin and Grelle (2012), where the marsupial showed a relationship with areas with a more closed canopy.

The marsupial M. americana, also captured in low abundance and only in preserved forest areas, has previously been found in disturbed areas (D’Andrea et al. 1999), although it is more common in forest environments (Alho et al. 1986, Mares et al. 1986). This species exhibited the same habitat associations observed for M. paraguayana, which can be characterized by the presence of more trees and a closed canopy. These two species could represent visitor species in this interface area, presenting occasional captures of individuals coming from more preserved areas in the upper parts of the forest.

The absence of traps in the midstory and in the upper strata may have underestimated the occurrence of G. microtarsus, since it is an arboreal marsupial (Paglia et al. 2012) and it was captured in low abundance only in 2001 samplings. A study carried out in the Pantanal biome indicated that the congener species G. agillis had significantly more captures in trees than on the ground (Andreazzi et al. 2011).

Nestedness distribution patterns are related to certain community processes because, according to the model, species occurring in species-poor sites are subsets of species occurring in species-rich sites (Patterson and Atmar 1986). The understanding of community structure patterns allows us to evaluate niche specialization, and is an important feature in the study of ecological networks (Sugihara and Ye 2009) and in conservation programs (Patterson 1987, Southwood 1996, Horváth et al. 2011). The observed small mammal community structure presented a pattern of low level of nestedness in both periods, differing from a random pattern, which indicates that some species are dominant along the environmental gradient, such as D. aurita and the opportunist rodents. Others might have a stochastic fluctuation in their distributions, which is a characteristic of species poor communities, since less generalist species, such as M. americana and M. paraguayana, occurred in species-rich sites, especially in the 2012–2015 period.

The higher degree of nestedness recorded for the 2012–2015 data in relation to 2001 may reinforce the higher dominance of the synanthropic and generalist species in the most recent sampling. There is also the possibility that there may be other species not sampled in the study, mainly in the more preserved and upper areas of PBSP. Menezes and Fernandez (2013) also found low nestedness values in a meta-analysis review of mammal communities in fragmented landscapes in different regions of the world. Surveys together with ecological studies on populations and communities are urgent in those areas, not only because of biodiversity concerns but also to expand understanding of the mechanisms of the ecological patterns in multiple scales of spatial-temporal organization, taking into account the metacommunity scale and the patterns and processes that emerge from it (Leibold et al. 2004).

This study is the first report to evaluate the community structure of small mammals in the sylvatic-urban interface area of PBSP, the largest forest reserve within an urban area in Brazil. These surveys indicate that the small mammal diversity was low in both sampling periods in spite of the preservation and forest restoration program recently implemented in CFMA areas, and a species loss in the Pau da Fome locality, despite it is a conservation unit. The greater species abundance in the most disturbed areas suggests the increase of factors favoring the occurrence of synanthropic species, such as R. rattus, the generalist marsupial D. aurita and the opportunistic rodents, A. cursor and O. nigripes. These results have important consequences not only for conservation but also for public health, providing potential subsidies for future public policy actions.

We would like to thank the staff and students of Laboratório de Biologia e Parasitologia de Mamíferos Silvestres Reservatórios, Laboratório de Hantaviroses e Rickettsioses, Laboratório de Biologia de Tripanossomatíedeos and Laboratório de Pesquisa Clínica em Dermatozoonose from FIOCRUZ who helped in the field work; to the coordinator of FIOCRUZ Atlantic Forest Campus, Gilson Antunes, for providing local facilities and supporting the general project; to the staff of FIOCRUZ Atlantic Forest Campus, especially to Marta L. Brandão and José L.P. Cordeiro for local facilities and help in the field work; to Marcelo F. Freitas for helping in the field work in all campaigns; to Cibele R. Bonvicino for rodents identification, to Fabiana P. Caramaschi for marsupials identification, to Daniela Dias for the curator of the collected animals, to Natália A. Costa and Fernando O. Santos for helping in the figures edition, to Instituto Estadual do Meio Ambiente (INEA) for the collecting licence and for local facilities and to an anonymous reviewer for suggestions on the manuscript. This project was financially supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico-CNPq (485074/2012-5 and 306352/2014-1), Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro-FAPERJ (E-26/103.285/2011 and E-26/111.296/2014), and Instituto Oswaldo Cruz (FIOCRUZ) (POM and PAEF). TSC, SFCN and BRT received grants from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).