BACKGROUND
one of the approaches that has shown to be effective in the analysis of multi-species assemblages is to consider resource utilization pattern data  <cit> .
the examination of resource utilization patterns reveals aspects of species coexistence in an ecological community  <cit> .
the degree of spatial overlap in resource use is an important parameter to explain community organization and the pattern of species coexistence  <cit> .
interspecific resource competition within larval containers , all of which are food-limiting environments, results in the simultaneous presence of multiple mosquito species.
for the bromeliad-dwelling mosquito community, gilbert et al.
<cit>  proposed two coexistence mechanisms among mosquitoes.
one provides co-occurrence of species within a given larval container via specialization of feeding behavior, while the other performs spatial segregation of species with overlapping use of resource by specializing them at different bromeliad sizes.
the third mechanism of coexistence was assessed from the best studied pair of competitors, aedes albopictus and aedes aegypti, regarding the asymmetrical competition, which is significantly associated with the diversity of resources in a landscape scale  <cit> .
species coexistence may occur in heterogeneous environments by spatial storage effect: covariance between the environment and competition  <cit> .
for container-dwelling mosquitoes, different container types may represent a tradeoff for mosquito larvae between the risk of dry out and presence of large predators  <cit> .
the differences in the species responses to this kind of conflicting pressures affect the competitive outcomes at different bromeliads.
as the distribution and abundances of bromeliads and other larval sites vary in space and time, spatial heterogeneity may promote coexistence by changing the scale in which competing species coexists  <cit> .
here we focus on host-seeking adult females and larvae of mosquito species across a gradient of vegetation physiognomies in a natural landscape of atlantic rainforest.
previous research on mosquito communities has provided evidence for coexistence mechanisms within the scope of container-dwelling mosquitoes.
herein, coexistence mechanisms in the whole mosquito community, from species occupying ground waters to other that occupy natural containers, are considered at micro-habitat, habitat and landscape scales.
moreover, most of the mosquito surveys performed in areas of tropical rainforest in southeastern brazil was focused on testing hypotheses concerning particular mosquito fauna along a gradient of anthropogenic modifications  <cit> .
this approach has led to baseline knowledge about mosquito species’ responses under the influences of human-caused changes in ecosystems, reinforcing the importance of synanthropy .
interpreting results through the lens of ecosystem changes, these studies left gaps in the knowledge about properties of natural and undisturbed mosquito communities in tropical rainforest.
our approach is to fill some of these gaps, by testing hypotheses about mosquito species coexistence.
more specifically we were interested in investigating co-occurrences of potential competitors at the micro-habitat scale   <cit> , habitat scale   <cit>  and landscape scale   <cit> .figure  <dig>
coexistence mechanisms.
a) larva of the species  <dig>  and larva of the species  <dig>  can co-occur at the micro-habitat  if they can vertically partition food and space .
b) the coexistence among larvae of competitor species is mediated by vegetation heterogeneity that allows occurrence of the sp <dig> in the ground  and the sp <dig> in the canopy .
c) along a physiognomic gradient of vegetation, interactions among potential competitors  are determined by the landscape context, causing one species to be superior competitor in a physiological stressful scenario  and the other as being an efficient competitor in the forest.
methods
study system
the brazilian atlantic forest is a tropical biome with mountainous forests and seasonally flooded coastal ecosystems  <cit> .
within the atlantic forest domain, regions of the coastal plain contain forest remnants  that include dense ombrophilous forests intermixed with mangroves, dune ecosystems and coastal islands .
bernardi et al.
<cit>  measured abiotic and biotic variables and applied remote sensing techniques to produce a thematic map of the vegetation types in the parque estadual da ilha do cardoso  .
these authors also indicated that physiological stress factors tend to decrease along the physiognomic gradient, with inland forests containing high levels of above ground biomass and soil nutrients.
studying bromeliad food webs in an understory of mixed arboreal vegetation in the peic, romero and srivastava  <cit>  found a diverse arthropod fauna composed of aquatic insect larvae of several functional groups.
the functional group of filter feeders is represented by the mosquito community in which the larvae assemblage, composed of different species, may experience interspecific resource competition.figure  <dig>
study area.
a) the atlantic forest domain and its forest remnants.
b) atlantic forest remnants in são paulo and paraná states.
c) location of the parque estadual da ilha do cardoso .
d) collection sites for measurements of vegetation complexity, the number of ground and phytotelmata oviposition sites and the number and identification of immature and adult mosquitoes.
source: conservation international , critical ecosystem partnership funding , instituto nacional de pesquisas espaciais .
sampling design
mosquito assemblages were assessed from  <dig> collection sites that spanned physiognomies varying from scrub vegetation to dense forest in the peic .
each collection site was composed of a cross consisting of four 15-m long transects .
previous long-term field collections conducted in the southeastern atlantic forest at the end of a hot and rainy summer showed that mosquito densities were high and that microclimatic factors were favorable for maintaining high levels of diversity in mosquito assemblages  <cit> .
given those findings, we sampled mosquitoes over 20 days in both march  <dig> and april  <dig>  the data were aggregated across years in subsequent analyses.
adults were collected at the middle  in figure 2d using automatic cdc-co <dig> as described by laporta and sallum  <cit>   and cdc-light traps  over a 24-hour sampling period.
each cdc trap was set  <dig>  m above ground level to attract both canopy and ground mosquitoes  <cit> .
measures of vegetation complexity and its heterogeneity were assessed at the middle of each transect .
the numbers of phytotelmata and ground pools were counted in each transect .
immature specimens were sampled in phytotelmata and ground pools .figure  <dig>
mosquito larvae collections across the gradient of forest physiognomies.
this drawing represents the vegetation complexity and mosquito larvae natural containers and ground waters which varied accordingly as follows: a) ephemeral ground pool; b) terrestrial bromeliad plant ; c) terrestrial bromeliad plant ; d) ephemeral ground pool; e) epiphytic bromeliad plant ; f) ephemeral ground pool; g) epiphytic bromeliad plant ; h) leaf axils in palm tree ; i) leaf axils in phytotelmata ; j) hollow internode of bamboo plant ; k) tree hole in long-lived tree; l) ephemeral ground pool; m) flower bracts in phytotelmata  and n) margins of permanent lotic waters.
adult mosquitoes collected in the field were maintained inside collection vials with silica to prevent fungal growth and damage to the specimens.
mosquito larvae were maintained in the laboratory to obtain males and females associated with both larval and pupal exuviae when possible.
species were identified based on external morphological characters using identification keys  <cit> .
mosquito nomenclature followed knight and stone  <cit> .
larvae identifications were undertaken using procedures adopted in marques et al.
<cit> .
the complexity of vegetation gradient was assessed by  <dig> measurements taken in each collection site .
data from vegetation complexity was obtained following the protocol described by pardini et al.
<cit> .
accordingly, foliage density  was measured for plants of each of the five forest strata 0- <dig>  1- <dig>  5- <dig>  10- <dig>  > <dig> meters above the ground.
the length of plant leaves was measured, along a vertical column, with a laser rangefinder bosch™ dle <dig>  by considering a vertical line from the bottom of the plant to the top of the leaves when straightened.
a 4-meter pole was also utilized to help to establish a vertical line.
foliage density values were estimated for each stratum and for each transect in all collection sites.
a principal component analysis  on the values of foliage density for all strata and for all transects was conducted.
the inter-transect mean of eigenvalues of the first axis was estimated for each collection site to represent the vegetation complexity.
the coefficient of variation of vegetation complexity was estimated and then it was utilized as indicative of heterogeneity of vegetation complexity.
the present research was certificated with the authorization by the instituto brasileiro do meio ambiente e dos recursos naturais renováveis  of the brazilian ministry of natural environment .
data analysis
response patterns for the distribution of mosquito species as adults were analyzed as a function of the vegetation complexity.
the response patterns of potential competitors were used as exploratory analysis of species co-occurrences or spatial segregation.
a model selection approach was applied using huisman-olff-fresco  models and extensions as described by jansen and oksanen  <cit>  as follow: i  - there was no response pattern along the vegetation complexity; ii  - the response pattern was linearly correlated with the vegetation complexity; iii  - the response pattern was partially linear and partially asymptotic; iv  - the response pattern was similar to a bell-shaped response curve; v  - the response pattern was a skewed non-monotonic response curve; vi  - the response pattern had a symmetric bimodal shape; vii  - the response pattern had a bimodal shape with two skewed maxima.
five model parameters, a-d and f, were optimised during the process of model fitting.
the number of parameters decreased with increasing model simplicity: model vii had five parameters , models vi and v had four parameters , models iv and iii had three parameters , model ii had two parameters  and model i had one parameter .
the process of model fitting for the hof models and their extensions was based on the maximisation of a log-likelihood function using iterative non-linear methods to estimate parameters a-d and f. the process of model fitting was applied with the ehof package in the r computational environment  <dig> .
<dig>  <cit> .
the model selection approach was based on the akaike information criteria corrected for small samples .
the best model was selected as that associated with the smallest aic value, and a value of evidence of ∆aicc > <dig> was adopted to differentiate plausible models.
a value of evidence of ∆aicc > <dig> is equivalent to three consecutive successes in a bernoulli independent trial.
because candidate models can adjust similarly to the data, we also considered functional redundancy, defined as the presence of both simple and complex models with similar response patterns  <cit> .
co-occurrences or spatial segregation of mosquito species as larvae were evaluated in micro-habitat, habitat and landscape scales in order to determine evidences of species coexistence.
thus, logistic binary regression models were employed in order to fit biological data within r  <dig> .
<dig> programming environment with the aid of epicalc package  <cit> .
the logistic binary regression model was applied such that species presence-absence was fitted with another species presence-absence in a specific set of bromeliads or ephemeral ground pools.
odds ratio  was estimated in order to assess associations between the two potential competitor species, under the null hypothesis or = <dig> .
an or statistically higher than  <dig>  provides evidence for species coexistence in micro-habitat scale.
for an or lower than  <dig> , species coexistence may occur either at habitat or landscape scales.
in order to see overall responses from a community perspective, hof model selection was employed to determine whether species richness response patterns varied monotonically or non-monotonically against vegetation gradient.
furthermore, the coefficient of variance of vegetation structure, i.e., the square root of inter-transect variance divided by inter-transect mean of eigenvalues of the first axis of pca, was estimated and then it was utilized as indicative of heterogeneity of vegetation structure, according to chaves et al.
<cit> .
hof model selection was again employed to determine whether heterogeneity of vegetation structure can be the underlying mechanism that influences patterns of diversity  of the mosquito fauna.
finally, 95% confidence intervals  were calculated by applying gam  with the employ of mgcv package in the r  <dig> .
<dig> program  <cit> .
RESULTS
after 720 hours of automatic trap sampling effort,  <dig>  adult specimens  were collected from the mosquito assemblages.
coquillettidia chrysonotum was the most abundant species .
other species were, for example, aedes scapularis , anopheles cruzii , wyeomyia quasilongirostris , aedes serratus , wyeomyia muehlensi , culex imitator  and anopheles bellator .
the sample also contained rare species, a total of  <dig> singletons.
no specimens of anopheles  or synanthropic species  were recorded.
mosquito larvae were collected from natural containers and ground waters according to figure  <dig>  although it was possible to collect wyeomyia aporonoma and sabethes intermedius from a tree hole in ficus spp.
in the forest, most specimens were obtained from bromeliads and ground pools.
an average of three bromeliads and one ground pool per site were searched for the presence of mosquito larvae, totaling  <dig> bromeliads and  <dig> ground pools.
mosquito fauna from permanent ground pools with emergent vegetation was not sampled because it is supposed to be less controlled by interspecific resource competition than that in small pools, thus larvae of cq.
chrysonotum and uranotaenia were not recorded.
in the herein study site, permanent ground pools were present along margins of the river that encounters the sea on the plain coast .
instead, temporary ground pools  were selected and provided mostly larvae of aedes , psorophora and culex .
both terrestrial and epiphytic bromeliads occurred along the gradient, terrestrial more common in scrublands  and epiphytic bromeliads in the forest .
in the middle of the gradient where mixed arboreal and scrub vegetation occurred, both terrestrial  and epiphytic  bromeliads were found, mainly in the genera quesnelia, vriesea and nidularium.
common mosquito inhabitants were anopheles , culex  and wyeomyia .
in total,  <dig> specimens were collected as larvae.
species that could be identified to species level were ae.
scapularis , ae.
serratus , an.
bellator , an.
cruzii , cx.
imitator , wy.
muhelensi  and wy.
quasilongirostris .
results of the pca on the values of foliage density showed that the first component explained approximately  <dig> % of the variation, with the highest importance associated with foliage densities in the strata that were >15 m, 10-15 m and 1-5 m above ground level.
this component was used as the main component of vegetation complexity.
on scrublands foliage density in the lower strata increased, whereas foliage density in the lower strata decreased in multi-layered forest physiognomies .
vegetation complexities in scrublands and dense forests did not change very much, the former being composed of scrub and the latter by a multi-layer forest.
in the middle of the gradient, however, vegetation complexities were a heterogeneous mix of arboreal and scrub vegetation presenting an irregular and spatially discontinuous mosaic of trees and bushes .
the relationships among mosquito potential competitor species as adults and vegetation gradient varied substantially, showing that some species could co-occur at habitat scale  and others that co-occurred at landscape scale  .
co-occurrences and spatial segregation of mosquito potential competitor species as larvae showed evidences of coexistence mechanisms at three scales: micro-habitat, habitat and landscape .
abundance distributions of an.
bellator and cx.
imitator were associated with scrublands .
in the micro-habitat on scrublands, the co-occurrence of an.
bellator and cx.
imitator in  <dig> bromeliads  was statistically significant .
abundance distributions of wy.
muhelensis and wy.
quasilongirostris were associated with the arboreal scrub vegetation .
although wy.
quasilongirostris and wy.
muhelensis could occur in the same habitat , they did not at micro-habitat scale and occurred together only in  <dig> bromeliads, with the former species mostly found in epiphytic bromeliads and the latter species in terrestrial bromeliads.
the two other pairs of species, an.
bellator-an.
cruzii and ae.
scapularis-ae.
serratus, could only co-occur on a landscape scale .
abundance distribution of an.
bellator was complementary to the abundance distribution of an.
cruzii  and larvae of an.
bellator and an.
cruzii did not frequently co-occur .
abundance distribution of ae.
scapularis was associated with scrubland, whereas that of ae.
serratus had two peaks, one at mid-gradient and the other in forest .
of the  <dig> ground pools examined, ae.
scapularis and ae.
serratus as larvae co-occurred only in  <dig>  and did not co-occur in  <dig> .figure  <dig>
observed pattern of co-occurrences or spatial segregation among potential competitors.
a)
an.
bellator and cx.
imitator co-occurred on scrublands at habitat scale.
b)
wy.
muehlensi and wy.
quasilongirostris co-occurred in the transition zone at habitat scale.
c)
an.
bellator and an.
cruzii and d)
ae.
scapularis and ae.
serratus co-occurred at landscape scale and spatially segregated at habitat scale.
gam for non-monotonic hof models and 95%ci are in additional file 7: figure s <dig>  additional file 9: figure s <dig>
co-occurrences or spatial segregation of mosquito potential competitor species as larvae, parque estadual da ilha do cardoso, 2009-2010
logistic regression output
pairs of species
type of containers/ground waters
or
p
pattern
support
process
an.
bellator - cx.
imitator
wy.
muehlensi - wy.
quasilongirostris
an.
bellator - an.
cruzii
ae.
scapularis - ae.
serratus
a: this is a significant result under the null hypothesis or = <dig>  thus there is co-occurrence between the two species.
b: this is a significant result under the null hypothesis or = <dig>  thus the two species do not co-occur.
the association between mosquito species richness and vegetation gradient had a bimodal hump-shaped pattern .
the vegetation heterogeneity also showed a bimodal hump-shaped pattern related to vegetation gradient .
mosquito species richness had a plateau correlation with vegetation gradient, i.e., it linearly increased up to a plateau .
the increase of vegetation heterogeneity was also associated with the rise of evenness .figure  <dig>
relationships between mosquito diversity and the gradient of forest physiognomies.
a) mosquito species richness has a bimodal, skewed response pattern.
the strong peak is related to the heterogeneity of vegetation complexity.
b) the heterogeneity of vegetation complexity has also a bimodal skewed response pattern, similar to mosquito species richness.
c) mosquito species richness increased linearly with heterogeneity of vegetation complexity until it reached an asymptote where a mix of arboreal and scrub vegetation occurred.
d) the mosquito assemblage evenness-dominance profile as a function of heterogeneity of vegetation complexity.
high levels of heterogeneity of vegetation complexity resulted in assemblage evenness, whereas medium to low levels of heterogeneity of vegetation complexity resulted in the presence of dominant assemblages.
gam for non-monotonic hof models and 95%ci are in additional file 8: figure s <dig>
analyses performed with gam, a form-free method, detected the same overall pattern yielded with the employ of hof models, i.e., a form-imposed function .
discussion
co-occurrences and spatial segregation of mosquito species were assessed at micro-habitat, habitat and landscape scales along vegetation gradient in a natural and undisturbed landscape of tropical rainforest, where ecologically and taxonomically similar mosquito species coexist.
these similar species usually exploit limiting resources in different ways, such as feeding at different heights inside the bromeliad water column at micro-habitat scale, colonizing either terrestrial or epiphytic bromeliads at the habitat scale and selecting more or less physically stressful habitats at landscape scale.
thus, resource partitioning may be an important structuring process in these assemblages.
mosquito species diversity was at least partly caused by the several ways in which potential competitors have evolved specialized traits.
anopheles bellator and cx.
imitator significantly co-occurred in the same bromeliad, especially on scrublands.
both species are frequent in quesnelia bromeliads that are highly exposed to sunlight.
in these bromeliads, larvae of both species collect and filter detritus and microorganisms  <cit> .
because of the lack of a respiratory siphon which characterizes all anopheles species, an.
bellator feed at the air-water interface, whilst cx.
imitator feeds mainly in the water column, as do other culex species  <cit> .
similarly, marques et al.
<cit>  observed statistically significant associations between other species of anopheles  and culex  in vriesea and nidularium bromeliads in lowlands and hillsides at another study site in the atlantic rainforest.
the non-random co-occurrence of these species at micro-habitat scale is indicative of the first mechanism of species coexistence in mosquito assemblages .
prior to gilbert et al.
<cit> , yee et al.
<cit>  proposed this mechanism after video recording the feeding behavior of cx.
pipiens, ae.
albopictus and ae.
triseriatus specimens in different food environments within an experimental study design.
interestingly, this is a consistent explanation for the co-occurrence of an.
bellator with cx.
imitator in the present work and an.
cruzii with other cx.
species and an.
homunculus with cx.
ocellatus in marques et al.
<cit> .
mosquito species richness varied non-linearly across the gradient of vegetation gradient.
more species occurred in forested habitats than on scrublands, showing that only some mosquito species can endure the physiological stress of harsh environments.
more species occurred in the forest than was expected, attributable to diverse phytotelmata-inhabiting mosquito species of the tribe sabethini  <cit> .
however, the maximum species richness was observed in the middle where mixed vegetation occurs, with both arboreal and scrub elements.
this represented a transition zone .
as expected, this is part of the gradient with highest heterogeneity of vegetation complexity.
the resulting mosaic of trees and bushes caused variation in light incidence and thus the occurrence of both terrestrial and epiphytic bromeliads.
wyeomyia  mosquito species occurred in this physiognomy.
the feeding behavior of these species should be very similar at the micro-habitat scale and because of that these species did not co-occur in the same type of bromeliad.
wyeomyia muhelensis occurred in terrestrial quesnelia bromeliads, whereas wy.
quasilongirostris occurred in epiphytic vriesea bromeliads.
hence, coexistence probably is possible because of the presence of these terrestrial and epiphytic bromeliad species and provides evidence in support of the second mechanism of species coexistence in mosquito assemblages .
the pairs of potential competitor species anopheles bellator-an.
cruzii and ae.
scapularis-ae.
serratus did not co-occur at micro-habitat and habitat.
however, all species did co-occur at the landscape scale.
because both pairs of species are taxonomically and ecologically similar, they cannot possibly allow their persistence at local scales which may suggest a mechanism of coexistence at the landscape scale.
this observational evidence may be related to the asymmetrical competition mechanism which results in a greater likelihood of local competitive exclusion, but it can provide a regional coexistence depending on the context  <cit> .
this finding represents evidence for the third mechanism of species coexistence in mosquito assemblages .
CONCLUSIONS
the present work did not experimentally test the effects of competition among species across the gradient of vegetation physiognomies to ensure that the observed patterns match the predicted mechanism.
thus, it was found herein empirical evidence for three mechanisms of resource-niche partitioning that provided greater biodiversity in the species-rich mosquito community of atlantic rainforest.
the mechanisms are: 1) spatial niche segregation in vertical strata in bromeliads, 2) niche partitioning mediated by spatial heterogeneity which promotes species coexistence at habitat scale and 3) asymmetrical competition with context dependence in groups of taxonomically similar species in bromeliads and ephemeral pools.
additional files
additional file 1: table s <dig> distribution of the foliage vertical length  per forest strata, transect, and collection site, parque estadual da ilha do cardoso, cananéia, são paulo state, brazil,  <dig>
additional file 2: table s <dig> results of the aic model selection approach with mosquito  individual abundance in function of vegetation gradient, parque estadual da ilha do cardoso, 2009- <dig>
additional file 3: table s <dig> distribution of the number of collected adult mosquitoes in cdc traps  per species and collection site, parque estadual da ilha do cardoso, cananéia, são paulo state, brazil, 2009- <dig>
additional file 4: figure s1-s <dig> the abundance of a given species as adults per collection site as a function of vegetation gradient.
additional file 5: table s <dig> distribution of numbers of larvae of mosquito species per bromeliad in a gradient of vegetation, parque estadual da ilha do cardoso, southeastern atlantic forest, brazil, 2009- <dig>
additional file 6: table s <dig> results of the aic model selection approach with mosquito species richness in function of vegetation gradient and its heterogeneity, parque estadual da ilha do cardoso, 2009- <dig>
additional file 7: figure s <dig> solid lines represent gam fitted to the data using cubic regression splines.
dashed lines define the 95% ci.
a) figure 4b-wy.
muhelensi.
b) figure 4b-wy.
quasilongirostris.
c) figure 4d-ae.
scapularis.
d) figure 4d-ae.
serratus.
additional file 8: figure s <dig> solid lines represent gam fitted to the data using cubic regression splines.
dashed lines define the 95% ci.
a) figure 5a.
b) figure 5b.
additional file 9: figure s <dig> solid lines represent gam fitted to the data using cubic regression splines.
dashed lines define the 95% ci.
a) figure 4a-an.
bellator.
b) figure 4a-cx.
imitator.
c) figure 4c-an.
bellator.
d) figure 4c-an.
cruzii.
abbreviations
aicakaike information criteria
cdccenters for disease control
ciconfidence interval
gamgeneralized additive models
hofhuisman-olff-fresco
pcaprincipal component analysis
peicparque estadual da ilha do cardoso
orodds ratio
competing interests
the authors declare that they have no competing interests.
authors’ contributions
gzl and mams designed the study.
gzl conducted the field and laboratory research.
mams provided equipment and materials.
gzl analysed the data.
gzl and mams wrote the manuscript.
both authors read and approved the final manuscript.
authors’ information
gzl is a biologist with ph.d. in science.
his main interests are landscape epidemiology, disease ecology and eco-epidemiology.
he holds a postdoctoral research position at university of são paulo.
mams is a professor with ph.d. in public health.
her main interests are culicidae systematic, evolution and ecology.
she holds an associate professor position at university of são paulo.

 one of the approaches that has shown to be effective in the analysis of multi-species assemblages is to consider resource utilization pattern data. the examination of resource utilization patterns reveals aspects of species coexistence in an ecological community. the degree of spatial overlap in resource use is an important parameter to explain community organization and the pattern of species coexistence. interspecific resource competition within larval containers, all of which are food-limiting environments, results in the simultaneous presence of multiple mosquito species. for the bromeliad-dwelling mosquito community, gilbert et al. proposed two coexistence mechanisms among mosquitoes. one provides co-occurrence of species within a given larval container via specialization of feeding behavior, while the other performs spatial segregation of species with overlapping use of resource by specializing them at different bromeliad sizes. the third mechanism of coexistence was assessed from the best studied pair of competitors, aedes albopictus and aedes aegypti, regarding the asymmetrical competition, which is significantly associated with the diversity of resources in a landscape scale. species coexistence may occur in heterogeneous environments by spatial storage effect: covariance between the environment and competition. for container-dwelling mosquitoes, different container types may represent a tradeoff for mosquito larvae between the risk of dry out and presence of large predators. the differences in the species responses to this kind of conflicting pressures affect the competitive outcomes at different bromeliads. as the distribution and abundances of bromeliads and other larval sites vary in space and time, spatial heterogeneity may promote coexistence by changing the scale in which competing species coexists. here we focus on host-seeking adult females and larvae of mosquito species across a gradient of vegetation physiognomies in a natural landscape of atlantic rainforest. previous research on mosquito communities has provided evidence for coexistence mechanisms within the scope of container-dwelling mosquitoes. herein, coexistence mechanisms in the whole mosquito community, from species occupying ground waters to other that occupy natural containers, are considered at micro-habitat, habitat and landscape scales. moreover, most of the mosquito surveys performed in areas of tropical rainforest in southeastern brazil was focused on testing hypotheses concerning particular mosquito fauna along a gradient of anthropogenic modifications. this approach has led to baseline knowledge about mosquito species’ responses under the influences of human-caused changes in ecosystems, reinforcing the importance of synanthropy. interpreting results through the lens of ecosystem changes, these studies left gaps in the knowledge about properties of natural and undisturbed mosquito communities in tropical rainforest. our approach is to fill some of these gaps, by testing hypotheses about mosquito species coexistence. more specifically we were interested in investigating co-occurrences of potential competitors at the micro-habitat scale, habitat scale and landscape scale.figure coexistence mechanisms. a) larva of the species and larva of the species can co-occur at the micro-habitat if they can vertically partition food and space. b) the coexistence among larvae of competitor species is mediated by vegetation heterogeneity that allows occurrence of the sp in the ground and the sp in the canopy. c) along a physiognomic gradient of vegetation, interactions among potential competitors are determined by the landscape context, causing one species to be superior competitor in a physiological stressful scenario and the other as being an efficient competitor in the forest. methods study system the brazilian atlantic forest is a tropical biome with mountainous forests and seasonally flooded coastal ecosystems. within the atlantic forest domain, regions of the coastal plain contain forest remnants that include dense ombrophilous forests intermixed with mangroves, dune ecosystems and coastal islands. bernardi et al. measured abiotic and biotic variables and applied remote sensing techniques to produce a thematic map of the vegetation types in the parque estadual da ilha do cardoso. these authors also indicated that physiological stress factors tend to decrease along the physiognomic gradient, with inland forests containing high levels of above ground biomass and soil nutrients. studying bromeliad food webs in an understory of mixed arboreal vegetation in the peic, romero and srivastava found a diverse arthropod fauna composed of aquatic insect larvae of several functional groups. the functional group of filter feeders is represented by the mosquito community in which the larvae assemblage, composed of different species, may experience interspecific resource competition.figure study area. a) the atlantic forest domain and its forest remnants. b) atlantic forest remnants in são paulo and paraná states. c) location of the parque estadual da ilha do cardoso. d) collection sites for measurements of vegetation complexity, the number of ground and phytotelmata oviposition sites and the number and identification of immature and adult mosquitoes. source: conservation international, critical ecosystem partnership funding, instituto nacional de pesquisas espaciais. sampling design mosquito assemblages were assessed from collection sites that spanned physiognomies varying from scrub vegetation to dense forest in the peic. each collection site was composed of a cross consisting of four 15-m long transects. previous long-term field collections conducted in the southeastern atlantic forest at the end of a hot and rainy summer showed that mosquito densities were high and that microclimatic factors were favorable for maintaining high levels of diversity in mosquito assemblages. given those findings, we sampled mosquitoes over 20 days in both march and april the data were aggregated across years in subsequent analyses. adults were collected at the middle in figure 2d using automatic cdc-co as described by laporta and sallum and cdc-light traps over a 24-hour sampling period. each cdc trap was set  m above ground level to attract both canopy and ground mosquitoes. measures of vegetation complexity and its heterogeneity were assessed at the middle of each transect. the numbers of phytotelmata and ground pools were counted in each transect. immature specimens were sampled in phytotelmata and ground pools.figure mosquito larvae collections across the gradient of forest physiognomies. this drawing represents the vegetation complexity and mosquito larvae natural containers and ground waters which varied accordingly as follows: a) ephemeral ground pool; b) terrestrial bromeliad plant; c) terrestrial bromeliad plant; d) ephemeral ground pool; e) epiphytic bromeliad plant; f) ephemeral ground pool; g) epiphytic bromeliad plant; h) leaf axils in palm tree; i) leaf axils in phytotelmata; j) hollow internode of bamboo plant; k) tree hole in long-lived tree; l) ephemeral ground pool; m) flower bracts in phytotelmata and n) margins of permanent lotic waters. adult mosquitoes collected in the field were maintained inside collection vials with silica to prevent fungal growth and damage to the specimens. mosquito larvae were maintained in the laboratory to obtain males and females associated with both larval and pupal exuviae when possible. species were identified based on external morphological characters using identification keys. mosquito nomenclature followed knight and stone. larvae identifications were undertaken using procedures adopted in marques et al.. the complexity of vegetation gradient was assessed by measurements taken in each collection site. data from vegetation complexity was obtained following the protocol described by pardini et al.. accordingly, foliage density was measured for plants of each of the five forest strata 0- 1- 5- 10- > meters above the ground. the length of plant leaves was measured, along a vertical column, with a laser rangefinder bosch™ dle by considering a vertical line from the bottom of the plant to the top of the leaves when straightened. a 4-meter pole was also utilized to help to establish a vertical line. foliage density values were estimated for each stratum and for each transect in all collection sites. a principal component analysis on the values of foliage density for all strata and for all transects was conducted. the inter-transect mean of eigenvalues of the first axis was estimated for each collection site to represent the vegetation complexity. the coefficient of variation of vegetation complexity was estimated and then it was utilized as indicative of heterogeneity of vegetation complexity. the present research was certificated with the authorization by the instituto brasileiro do meio ambiente e dos recursos naturais renováveis of the brazilian ministry of natural environment. data analysis response patterns for the distribution of mosquito species as adults were analyzed as a function of the vegetation complexity. the response patterns of potential competitors were used as exploratory analysis of species co-occurrences or spatial segregation. a model selection approach was applied using huisman-olff-fresco models and extensions as described by jansen and oksanen as follow: i - there was no response pattern along the vegetation complexity; ii - the response pattern was linearly correlated with the vegetation complexity; iii - the response pattern was partially linear and partially asymptotic; iv - the response pattern was similar to a bell-shaped response curve; v - the response pattern was a skewed non-monotonic response curve; vi - the response pattern had a symmetric bimodal shape; vii - the response pattern had a bimodal shape with two skewed maxima. five model parameters, a-d and f, were optimised during the process of model fitting. the number of parameters decreased with increasing model simplicity: model vii had five parameters, models vi and v had four parameters, models iv and iii had three parameters, model ii had two parameters and model i had one parameter. the process of model fitting for the hof models and their extensions was based on the maximisation of a log-likelihood function using iterative non-linear methods to estimate parameters a-d and f. the process of model fitting was applied with the ehof package in the r computational environment.. the model selection approach was based on the akaike information criteria corrected for small samples. the best model was selected as that associated with the smallest aic value, and a value of evidence of ∆aicc > was adopted to differentiate plausible models. a value of evidence of ∆aicc > is equivalent to three consecutive successes in a bernoulli independent trial. because candidate models can adjust similarly to the data, we also considered functional redundancy, defined as the presence of both simple and complex models with similar response patterns. co-occurrences or spatial segregation of mosquito species as larvae were evaluated in micro-habitat, habitat and landscape scales in order to determine evidences of species coexistence. thus, logistic binary regression models were employed in order to fit biological data within r. programming environment with the aid of epicalc package. the logistic binary regression model was applied such that species presence-absence was fitted with another species presence-absence in a specific set of bromeliads or ephemeral ground pools. odds ratio was estimated in order to assess associations between the two potential competitor species, under the null hypothesis or =. an or statistically higher than provides evidence for species coexistence in micro-habitat scale. for an or lower than, species coexistence may occur either at habitat or landscape scales. in order to see overall responses from a community perspective, hof model selection was employed to determine whether species richness response patterns varied monotonically or non-monotonically against vegetation gradient. furthermore, the coefficient of variance of vegetation structure, i.e., the square root of inter-transect variance divided by inter-transect mean of eigenvalues of the first axis of pca, was estimated and then it was utilized as indicative of heterogeneity of vegetation structure, according to chaves et al.. hof model selection was again employed to determine whether heterogeneity of vegetation structure can be the underlying mechanism that influences patterns of diversity of the mosquito fauna. finally, 95% confidence intervals were calculated by applying gam with the employ of mgcv package in the r. program.  after 720 hours of automatic trap sampling effort, adult specimens were collected from the mosquito assemblages. coquillettidia chrysonotum was the most abundant species. other species were, for example, aedes scapularis, anopheles cruzii, wyeomyia quasilongirostris, aedes serratus, wyeomyia muehlensi, culex imitator and anopheles bellator. the sample also contained rare species, a total of singletons. no specimens of anopheles or synanthropic species were recorded. mosquito larvae were collected from natural containers and ground waters according to figure  although it was possible to collect wyeomyia aporonoma and sabethes intermedius from a tree hole in ficus spp. in the forest, most specimens were obtained from bromeliads and ground pools. an average of three bromeliads and one ground pool per site were searched for the presence of mosquito larvae, totaling bromeliads and ground pools. mosquito fauna from permanent ground pools with emergent vegetation was not sampled because it is supposed to be less controlled by interspecific resource competition than that in small pools, thus larvae of cq. chrysonotum and uranotaenia were not recorded. in the herein study site, permanent ground pools were present along margins of the river that encounters the sea on the plain coast. instead, temporary ground pools were selected and provided mostly larvae of aedes, psorophora and culex. both terrestrial and epiphytic bromeliads occurred along the gradient, terrestrial more common in scrublands and epiphytic bromeliads in the forest. in the middle of the gradient where mixed arboreal and scrub vegetation occurred, both terrestrial and epiphytic bromeliads were found, mainly in the genera quesnelia, vriesea and nidularium. common mosquito inhabitants were anopheles, culex and wyeomyia. in total, specimens were collected as larvae. species that could be identified to species level were ae. scapularis, ae. serratus, an. bellator, an. cruzii, cx. imitator, wy. muhelensi and wy. quasilongirostris. results of the pca on the values of foliage density showed that the first component explained approximately % of the variation, with the highest importance associated with foliage densities in the strata that were >15 m, 10-15 m and 1-5 m above ground level. this component was used as the main component of vegetation complexity. on scrublands foliage density in the lower strata increased, whereas foliage density in the lower strata decreased in multi-layered forest physiognomies. vegetation complexities in scrublands and dense forests did not change very much, the former being composed of scrub and the latter by a multi-layer forest. in the middle of the gradient, however, vegetation complexities were a heterogeneous mix of arboreal and scrub vegetation presenting an irregular and spatially discontinuous mosaic of trees and bushes. the relationships among mosquito potential competitor species as adults and vegetation gradient varied substantially, showing that some species could co-occur at habitat scale and others that co-occurred at landscape scale. co-occurrences and spatial segregation of mosquito potential competitor species as larvae showed evidences of coexistence mechanisms at three scales: micro-habitat, habitat and landscape. abundance distributions of an. bellator and cx. imitator were associated with scrublands. in the micro-habitat on scrublands, the co-occurrence of an. bellator and cx. imitator in bromeliads was statistically significant. abundance distributions of wy. muhelensis and wy. quasilongirostris were associated with the arboreal scrub vegetation. although wy. quasilongirostris and wy. muhelensis could occur in the same habitat, they did not at micro-habitat scale and occurred together only in bromeliads, with the former species mostly found in epiphytic bromeliads and the latter species in terrestrial bromeliads. the two other pairs of species, an. bellator-an. cruzii and ae. scapularis-ae. serratus, could only co-occur on a landscape scale. abundance distribution of an. bellator was complementary to the abundance distribution of an. cruzii and larvae of an. bellator and an. cruzii did not frequently co-occur. abundance distribution of ae. scapularis was associated with scrubland, whereas that of ae. serratus had two peaks, one at mid-gradient and the other in forest. of the ground pools examined, ae. scapularis and ae. serratus as larvae co-occurred only in and did not co-occur in.figure observed pattern of co-occurrences or spatial segregation among potential competitors. a) an. bellator and cx. imitator co-occurred on scrublands at habitat scale. b) wy. muehlensi and wy. quasilongirostris co-occurred in the transition zone at habitat scale. c) an. bellator and an. cruzii and d) ae. scapularis and ae. serratus co-occurred at landscape scale and spatially segregated at habitat scale. gam for non-monotonic hof models and 95%ci are in additional file 7: figure s additional file 9: figure s co-occurrences or spatial segregation of mosquito potential competitor species as larvae, parque estadual da ilha do cardoso, 2009-2010 logistic regression output pairs of species type of containers/ground waters or p pattern support process an. bellator - cx. imitator wy. muehlensi - wy. quasilongirostris an. bellator - an. cruzii ae. scapularis - ae. serratus a: this is a significant result under the null hypothesis or = thus there is co-occurrence between the two species. b: this is a significant result under the null hypothesis or = thus the two species do not co-occur. the association between mosquito species richness and vegetation gradient had a bimodal hump-shaped pattern. the vegetation heterogeneity also showed a bimodal hump-shaped pattern related to vegetation gradient. mosquito species richness had a plateau correlation with vegetation gradient, i.e., it linearly increased up to a plateau. the increase of vegetation heterogeneity was also associated with the rise of evenness.figure relationships between mosquito diversity and the gradient of forest physiognomies. a) mosquito species richness has a bimodal, skewed response pattern. the strong peak is related to the heterogeneity of vegetation complexity. b) the heterogeneity of vegetation complexity has also a bimodal skewed response pattern, similar to mosquito species richness. c) mosquito species richness increased linearly with heterogeneity of vegetation complexity until it reached an asymptote where a mix of arboreal and scrub vegetation occurred. d) the mosquito assemblage evenness-dominance profile as a function of heterogeneity of vegetation complexity. high levels of heterogeneity of vegetation complexity resulted in assemblage evenness, whereas medium to low levels of heterogeneity of vegetation complexity resulted in the presence of dominant assemblages. gam for non-monotonic hof models and 95%ci are in additional file 8: figure s analyses performed with gam, a form-free method, detected the same overall pattern yielded with the employ of hof models, i.e., a form-imposed function. discussion co-occurrences and spatial segregation of mosquito species were assessed at micro-habitat, habitat and landscape scales along vegetation gradient in a natural and undisturbed landscape of tropical rainforest, where ecologically and taxonomically similar mosquito species coexist. these similar species usually exploit limiting resources in different ways, such as feeding at different heights inside the bromeliad water column at micro-habitat scale, colonizing either terrestrial or epiphytic bromeliads at the habitat scale and selecting more or less physically stressful habitats at landscape scale. thus, resource partitioning may be an important structuring process in these assemblages. mosquito species diversity was at least partly caused by the several ways in which potential competitors have evolved specialized traits. anopheles bellator and cx. imitator significantly co-occurred in the same bromeliad, especially on scrublands. both species are frequent in quesnelia bromeliads that are highly exposed to sunlight. in these bromeliads, larvae of both species collect and filter detritus and microorganisms. because of the lack of a respiratory siphon which characterizes all anopheles species, an. bellator feed at the air-water interface, whilst cx. imitator feeds mainly in the water column, as do other culex species. similarly, marques et al. observed statistically significant associations between other species of anopheles and culex in vriesea and nidularium bromeliads in lowlands and hillsides at another study site in the atlantic rainforest. the non-random co-occurrence of these species at micro-habitat scale is indicative of the first mechanism of species coexistence in mosquito assemblages. prior to gilbert et al., yee et al. proposed this mechanism after video recording the feeding behavior of cx. pipiens, ae. albopictus and ae. triseriatus specimens in different food environments within an experimental study design. interestingly, this is a consistent explanation for the co-occurrence of an. bellator with cx. imitator in the present work and an. cruzii with other cx. species and an. homunculus with cx. ocellatus in marques et al.. mosquito species richness varied non-linearly across the gradient of vegetation gradient. more species occurred in forested habitats than on scrublands, showing that only some mosquito species can endure the physiological stress of harsh environments. more species occurred in the forest than was expected, attributable to diverse phytotelmata-inhabiting mosquito species of the tribe sabethini. however, the maximum species richness was observed in the middle where mixed vegetation occurs, with both arboreal and scrub elements. this represented a transition zone. as expected, this is part of the gradient with highest heterogeneity of vegetation complexity. the resulting mosaic of trees and bushes caused variation in light incidence and thus the occurrence of both terrestrial and epiphytic bromeliads. wyeomyia mosquito species occurred in this physiognomy. the feeding behavior of these species should be very similar at the micro-habitat scale and because of that these species did not co-occur in the same type of bromeliad. wyeomyia muhelensis occurred in terrestrial quesnelia bromeliads, whereas wy. quasilongirostris occurred in epiphytic vriesea bromeliads. hence, coexistence probably is possible because of the presence of these terrestrial and epiphytic bromeliad species and provides evidence in support of the second mechanism of species coexistence in mosquito assemblages. the pairs of potential competitor species anopheles bellator-an. cruzii and ae. scapularis-ae. serratus did not co-occur at micro-habitat and habitat. however, all species did co-occur at the landscape scale. because both pairs of species are taxonomically and ecologically similar, they cannot possibly allow their persistence at local scales which may suggest a mechanism of coexistence at the landscape scale. this observational evidence may be related to the asymmetrical competition mechanism which results in a greater likelihood of local competitive exclusion, but it can provide a regional coexistence depending on the context. this finding represents evidence for the third mechanism of species coexistence in mosquito assemblages.  the present work did not experimentally test the effects of competition among species across the gradient of vegetation physiognomies to ensure that the observed patterns match the predicted mechanism. thus, it was found herein empirical evidence for three mechanisms of resource-niche partitioning that provided greater biodiversity in the species-rich mosquito community of atlantic rainforest. the mechanisms are: 1) spatial niche segregation in vertical strata in bromeliads, 2) niche partitioning mediated by spatial heterogeneity which promotes species coexistence at habitat scale and 3) asymmetrical competition with context dependence in groups of taxonomically similar species in bromeliads and ephemeral pools. additional files additional file 1: table s distribution of the foliage vertical length per forest strata, transect, and collection site, parque estadual da ilha do cardoso, cananéia, são paulo state, brazil, additional file 2: table s results of the aic model selection approach with mosquito individual abundance in function of vegetation gradient, parque estadual da ilha do cardoso, 2009- additional file 3: table s distribution of the number of collected adult mosquitoes in cdc traps per species and collection site, parque estadual da ilha do cardoso, cananéia, são paulo state, brazil, 2009- additional file 4: figure s1-s the abundance of a given species as adults per collection site as a function of vegetation gradient. additional file 5: table s distribution of numbers of larvae of mosquito species per bromeliad in a gradient of vegetation, parque estadual da ilha do cardoso, southeastern atlantic forest, brazil, 2009- additional file 6: table s results of the aic model selection approach with mosquito species richness in function of vegetation gradient and its heterogeneity, parque estadual da ilha do cardoso, 2009- additional file 7: figure s solid lines represent gam fitted to the data using cubic regression splines. dashed lines define the 95% ci. a) figure 4b-wy. muhelensi. b) figure 4b-wy. quasilongirostris. c) figure 4d-ae. scapularis. d) figure 4d-ae. serratus. additional file 8: figure s solid lines represent gam fitted to the data using cubic regression splines. dashed lines define the 95% ci. a) figure 5a. b) figure 5b. additional file 9: figure s solid lines represent gam fitted to the data using cubic regression splines. dashed lines define the 95% ci. a) figure 4a-an. bellator. b) figure 4a-cx. imitator. c) figure 4c-an. bellator. d) figure 4c-an. cruzii. abbreviations aicakaike information criteria cdccenters for disease control ciconfidence interval gamgeneralized additive models hofhuisman-olff-fresco pcaprincipal component analysis peicparque estadual da ilha do cardoso orodds ratio competing interests the authors declare that they have no competing interests. authors’ contributions gzl and mams designed the study. gzl conducted the field and laboratory research. mams provided equipment and materials. gzl analysed the data. gzl and mams wrote the manuscript. both authors read and approved the final manuscript. authors’ information gzl is a biologist with ph.d. in science. his main interests are landscape epidemiology, disease ecology and eco-epidemiology. he holds a postdoctoral research position at university of são paulo. mams is a professor with ph.d. in public health. her main interests are culicidae systematic, evolution and ecology. she holds an associate professor position at university of são paulo. 
environmental characteristics, mosquito fauna as adults and larvae were assessed along vegetation gradient in a natural landscape of tropical rainforest. biodiversitybiotic interactionscoexistencecommunityecologymosquitoesresource partitioningtropical rainforestissue-copyright-statement© the author 2014 species coexistence in mosquito assemblages may depend on mechanisms related to interspecific resource partitioning occurring at multiple scales. in the present work we investigated co-occurrence or spatial segregation in mosquito assemblages sharing resources at micro-habitat, habitat and landscape scales. huisman-olff-fresco and generalized additive models were employed to explore relationships between abundances of potential competitors in mosquito assemblages and vegetation gradient. co-occurrences and spatial segregation of mosquito species showed evidences of three scales of coexistence mechanisms: 1) micro-habitat - scale 1: different behaviors in response to food availability in specific vertical strata within larval container; 2) habitat - scale 2: specialized strategies related to heterogeneity of resource availability among larval containers and 3) landscape - scale 3: asymmetrical competition dependent upon the context of abiotic and biotic variables. results of the present work suggest that coexistence mechanisms can concomitantly work at multiple scales. we tested hypotheses concerning mosquito species co-occurrence or spatial segregation employing binomial logistic regression models.
 species coexistence in mosquito assemblages may depend on mechanisms related to interspecific resource partitioning occurring at multiple scales. in the present work we investigated co-occurrence or spatial segregation in mosquito assemblages sharing resources at micro-habitat, habitat and landscape scales. environmental characteristics, mosquito fauna as adults and larvae were assessed along vegetation gradient in a natural landscape of tropical rainforest. huisman-olff-fresco and generalized additive models were employed to explore relationships between abundances of potential competitors in mosquito assemblages and vegetation gradient. we tested hypotheses concerning mosquito species co-occurrence or spatial segregation employing binomial logistic regression models.  co-occurrences and spatial segregation of mosquito species showed evidences of three scales of coexistence mechanisms: 1) micro-habitat - scale 1: different behaviors in response to food availability in specific vertical strata within larval container; 2) habitat - scale 2: specialized strategies related to heterogeneity of resource availability among larval containers and 3) landscape - scale 3: asymmetrical competition dependent upon the context of abiotic and biotic variables.  results of the present work suggest that coexistence mechanisms can concomitantly work at multiple scales. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords biodiversitybiotic interactionscoexistencecommunityecologymosquitoesresource partitioningtropical rainforestissue-copyright-statement© the author 2014 