PLoS ONE
Home Colonization of different biomes drove the diversification of the Neotropical Eidmanacris crickets (Insecta: Orthoptera: Grylloidea: Phalangopsidae)
Colonization of different biomes drove the diversification of the Neotropical <i>Eidmanacris</i> crickets (Insecta: Orthoptera: Grylloidea: Phalangopsidae)
Colonization of different biomes drove the diversification of the Neotropical Eidmanacris crickets (Insecta: Orthoptera: Grylloidea: Phalangopsidae)

Competing Interests: The authors have declared that no competing interests exist.

Article Type: research-article Article History
Abstract

The phylogeny of the cricket genus Eidmanacris is used to analyse its historical distribution and diversification in three South American biomes: Atlantic Forest, Cerrado, and Chiquitano Dry Forest. A morphological phylogeny with all the 29 species of Eidmanacris and the Geographically explicit Event Model (GEM) is used to explain their colonization and diversification through three different biomes and their ancestral habitats and distributional areas. We analysed ecologically-significant characters, such as body size and metanotal characters, to test whether if morphology, habitat, or behaviour are connected. The relations of these features with the colonisation of wetter or drier biomes based on the distributional area, phylogeny and diversity of the genus were also tested. The results show that the ancestral distribution of the genus was the Atlantic Forest, and that biome occupancy, habitat, size, and mating behaviour evolved congruently through the phylogeny, drawing a coherent pattern of changes through Eidmanacris evolution toward the colonisation of drier biomes. Our results indicate that gallery forests could play a key role in the distribution and diversification of Eidmanacris species.

de Campos,de Souza-Dias,Desutter-Grandcolas,Shigueo Nihei,and ZHANG: Colonization of different biomes drove the diversification of the Neotropical Eidmanacris crickets (Insecta: Orthoptera: Grylloidea: Phalangopsidae)

Introduction

The Neotropical Region is the most species-rich region of the World, and the origins, diversification, and evolution of Neotropical organisms have been main subjects of interest for the last 150 years [1]. Until recently, most studies have focused on the Amazon Rainforest [2, 3]. However, other biomes have played a central role in the diversification of Neotropical biotas, like the Atlantic Forest and the drier vegetations located south and west of the Amazon Rainforest, i.e. the dry forests, Cerrado, and Chaco [4].

With more 8000 endemic species, the Atlantic Forest is the second-largest tropical forest in South America after the Amazon and is considered the most threatened diversity hotspot in South America [5]. The Atlantic Forest covers a wide range of vegetation formations, from tropical to subtropical and with heterogeneous composition, and covers elevations from the sea level to 2900 m [6]. It is separated from the Amazon Rainforest by the diagonal open formations, composed by Caatinga, Chaco and Cerrado [7], all equally threatened. The deciduous and semideciduous Chiquitano Dry Forest is a transition zone between the moist Amazon Rainforest and the dry forest of the Chaco [8], located in the east of Santa Cruz, Bolivia, and extending into western of Mato Grosso, Brazil. That forest is characterized by a pronounced dry season lasting 5–6 months per year, with a sparse and lower vegetation than tropical rain forests [7, 9]. Furthermore, it is one of the most diverse dry forests in the world [10]. The Cerrado hotspot covers more than 2 million km2 of South America [5, 11], between the Amazon, the southern Atlantic Forest, Caatinga and Chaco [1214]. It is a savannah vegetation with nutrient-poor, deep and well-drained soils [15], which physiognomies vary from an open field (“Campo limpo”) to a relatively tall and closed forest (“Cerradão”) [16]. It is crossed by gallery forests mainly formed by headwaters of rivers located in the plateaus of the Cerrado region [17].

The diversity of landscapes and environments of these biomes and their ecotone transitions may have represented evolutionary opportunities or an obstacle for the diversifications and/or dispersion of forest biotas. Several studies support the hypothesis that the Atlantic Forest was connected to the Amazon Rainforest in the past through gallery forests and/or by forests expansions [1823], which would have allowed the colonisation of the Cerrado by forest-dependent organisms, e.g. [2427], as well as faunal exchange between the two large South American forests through mesic corridors. Some groups of organisms are common to both forests, i.e. wood trees [7, 23], small mammals like rodents and marsupials [28, 29] and butterflies [21]. However, some vertebrate groups are exclusively present in the Amazon Rainforest, e.g. [30, 31], or in the Atlantic Forest, e.g. [32, 33], due to specific sensitivity to different climates, vicariance (geographical or ecological), specific niches and/or many other biotic or abiotic variables [1].

The colonisation of different biomes by a given evolutionary lineage, or its limitation to one biome, may be related to the biotic properties of the lineage, which must cope with environmental variables [34, 35]. This not only relates to climatic adaptation, but also resource use, and behavioural differences, in particular, reproduction.

Crickets are highly diversified in the Neotropical Region [3640], but they have been little used to test distributional hypotheses, due to the limited taxonomic knowledge, with many taxa still awaiting description, and the lack of phylogenetic analyses of well-supported monophyletic taxa [41]. In other hotspot areas, crickets are important model groups for studying modifications of insect communities through forest succession or the influence of invasive species [42, 43].

Here we analyse the diversification of the cricket genus Eidmanacris Chopard (Insecta, Orthoptera, Grylloidea, Phalangopsidae), which is distributed in the Atlantic Forest, Cerrado, and Chiquitano Dry Forest, but is absent from the Amazon Rainforest. Being flightless, the species of Eidmanacris, like other micropterous or apterous orthopterans, are potential models for studying migration corridors [44, 45]. Eidmanacris was recently revised, based on a large sampling effort through forested and open biomes [4648]. In the present paper, we perform an event-based biogeographical analysis (Geographically explicit Event Model–GEM) [49] using a morphology-based phylogenetic analysis including all the 29 species of Eidmanacris, to answer how did Eidmanacris colonise the three different biomes where it occurs today and what are the ancestral habitats and distributional areas. We also performed an ancestral state reconstruction analysis to examine ecologically-significant characters for the diversity of the genus, such as body size (resource use, climatic adaptation) and genital characters (mating behaviour), in order to verify if morphology, habitat, or behaviour are connected, and if these features are related to the colonisation of wetter or drier biomes. Furthermore, we discuss the conservation of the biomes in which Eidmanacris is present.

Materials and methods

Taxon sampling

The ingroup includes all the 29 species of Eidmanacris described to date (S1 Table) [48]. For E. paramarmorata Desutter-Grandcolas and E. longa Gorochov, we used data from the original descriptions and previous studies [50, 51]. High-resolution images of the holotypes of E. marmorata, E. paramarmorata and E. speluncae (Mello-Leitão) were accessed from the online catalog Orthoptera Species File [52].

The specimens examined belong to the following institutions: Laboratório de Insetos do Departamento de Zoologia da UNESP de Botucatu, São Paulo (UBTU); Laboratório de Orthopterologia da Universidade de Viçosa, Minas Gerais (UFV); Muséum national d’Histoire naturelle, Paris (MNHN); Museu de Zoologia da Universidade de São Paulo, São Paulo (MZSP); Museu de Zoologia da Universidade Estadual de Feira de Santana, Bahia (MZUEFS); The Academy of Natural Sciences of Drexel University, Philadelphia (ANSP). Additional material was collected in the field by LDC and PGBS and deposited in MZSP. Natural history observations were made during fieldworks. Outgroups were used to root the tree and polarize the characters [53]. Outgroups were chosen based on the phylogeny of Souza-Dias [54], which included four species of Eidmanacris. Nine outgroup taxa were selected (S2 Table): Adenopygus heikoi Bolfarini & de Mello; Bambuina bambui de Mello, Horta & Bolfarini; Guabamima lordelloi de Mello; Guabamima saiva de Mello; Melanotes ornata Desutter-Grandcolas; Ottedana cercalis de Mello & de Andrade; Strinatia brevipennis Chopard; Strinatia teresopolis Mesa, and an undescribed species of Modestozara Gorochov. Melanotes ornata was used to root the trees.

Characters and data matrix

This study was based on morphological characters of adult males and females, including genitalia. The specimens were examined under stereomicroscope Leica EZ4. The morphological terminology adopted, including the male phallic complex, follows Desutter [36, 55], Desutter-Grandcolas [56], and Souza-Dias [54]. Male phallic complexes were dissected and cleared in an aqueous solution 10% KOH for 24 hours and stored in vial with 80% ethanol together with the respective specimen. The female copulatory papillae were also dissected, but not cleared in KOH, and stored in small vials with ethanol 80% together with the respective specimen. Drawings of the genitalia were made under a Leica MZ9.5 stereomicroscope coupled with camera lucida. The photographs were taken through a Leica MZ16 stereomicroscope with a Leica DFC-420 camera, using the Leica Application Suite LAS V4.0 software, with specimens, male genitalia, and female copulatory papilla immersed in ethanol 80%. Fresh samples were not available in order to complement the phylogeny with molecular data.

The metanotal and forewing glands are important sources of characters in Eidmanacris [47, 48, 50, 57]. We used scanning electron microscopy (SEM) to examine the metanotum of 21 species of Eidmanacris and six outgroup species. Due to the reduced number of specimens for SEM analysis, the glands of five additional species of Eidmanacris and three of outgroup were examined using light microscopy to avoid damaging the specimens. The species studied using SEM are indicated in S1 and S2 Tables. For the SEM analysis, we dissected a male specimen and removed its forewing and thorax. The samples were dehydrated in a graded ethanol series until 100%, critical point dried using CO2 as intermediate, mounted on stubs, and coated with gold. The samples were examined using a Scanning Electron Microscope Zeiss SIGMA VP at the Instituto de Biociências da Universidade de São Paulo (Biosciences Institute of the University of São Paulo).

The description of the characters follows Sereno [58] and the data matrix (S3 Table) was constructed in Mesquite 3.6 [59]. Inapplicable characters were coded as “-” and unobserved characters as “?”.

Ninety-eight (98) characters (83 binary and 15 multistate) from external morphology were scored, six from the head, 22 from the thorax (including forewings and metanotal gland), six from the abdomen, four from the legs, 54 from male genitalia, and six female exclusive characters (including copulatory papilla). Some of the proposed characters were based on previous studies [54, 60, 61]. However, most characters are proposed here for the first time after a detailed morphological study of Eidmanacris following the taxonomic revision of the genus [48]. The list of characters with states, indexes (ci and ri), comments, optimisations, and illustrations are presented in S1-S18 Figs in the S1 File.

Phylogenetic analysis

Phylogenetic analyses were carried out in TNT 1.5 [62], using TBR as branch swapping with 1000 replications, holding 100 trees per replicate. All characters were treated as unordered with equal weights. A strict consensus tree was constructed in TNT whenever the analyses found more than one most parsimonious tree. The resulting trees and patterns of character evolution were visualised in Winclada [63], where consistency index (ci) [64] and retention index (ri) [65] were calculated. Ambiguous characters were optimised using both ACCTRAN [66] and DELTRAN [67]. These optimisations are used to describe the preferences of secondary gains or losses, even though ACCTRAN does not necessarily result in secondary losses and DELTRAN in secondary gains [68, 69]. Sometimes, ACCTRAN may result in illogical optimisations for inapplicable states [68], which happened in some characters in our analysis. Therefore, in this study, both optimisations were applied individually for each character and are indicated when applied. Bremer support values [70] were calculated using the TNT script “BREMER.RUN” to assess branch support.

Two analyses were performed: one with all the 38 terminals studied (analysis 1) and one with 37 terminals (analysis 2), excluding E. paramarmorata. Known only from females and containing 94% of missing data, this species acts as a “wildcard”, potentially generating big polytomies in consensus trees [71].

Biogeographical analysis

A total of 61 localities data were obtained from literature [4648, 50, 51, 7275], labels of the studied specimens, the catalog Orthoptera Species File [52], and from additional material deposited at MZSP and MNHN. All collected data were plotted on a South America map and edited in QGIS 3.4 [76]. Biomes (sensu Dinerstein et al. [14]) were used for discussion of species distribution patterns and construction of the distribution map. The biomes are: Atlantic Forest (Tropical and Subtropical Moist Broadleaf Forests); Cerrado (Tropical and Subtropical Grasslands, Savannas, and Shrublands); and Chiquitano Dry Forest (Tropical and Subtropical Dry Broadleaf Forests).

Characterization of habitats was based on the nocturnal activity of the species and modified from Desutter-Grandcolas [38]: 1-straminicolous (str), species active in the leaf litter; 2-straminicolous/cavicolous (str/cav), species active in the leaf litter and cavities, including the entrance of caves; 3-cavicolous (cav), species foraging only in cavities at ground level and in the entrance of caves. No troglobitic species (i.e. species that inhabit only the interior of caves) is known so far in Eidmanacris. The data were obtained from the same literature of localities [4648, 50, 51, 7275], labels of specimens and personal observation in the field (Table 1).

Table 1
Distributional, ecological, and morphological characteristics of Eidmanacris species.
CladeTaxaBiomeHabitatSizeMetanotum bristles
AE. fuscaAtlstr3Anterior region
E. melloiAtlstr1Anterior region
E. endophallicaAtlstr2Anterior region
E. minutaAtlstr0Anterior region
CE. bidentataAtlstr1Anterior region
EE. putuhraAtlstr1Anterior region
GE. tridentataAtlstr1Anterior region
E. simoesiAtlstr1Anterior region
E. eliethaeAtlstr2Anterior region
IE. papaveroiAtlstr2Anterior region
KE. fontanettiaeAtlstr1Anterior region
ME. larvaeformisAtlstr/cav4Anterior region
E. septentrionalisAtlstr/cav2Anterior region
E. speluncaeAtlcav3Anterior region
E. multispinosaAtlstr3Anterior region
N1E. dissimilisAtl/Cerstr/cav2Entire metanotum
E. meridionalisAtlstr/cav2Entire metanotum
E. alboannulataAtl/Cerstr/cav3Entire metanotum
E. suassunaiAtlstr/cav3Entire metanotum
N2E. scopulaCercav3Entire metanotum
E. desutteraeCerstr/cav3Entire metanotum
E. corumbataiCerstr/cav3Entire metanotum
E. gigasCerstr/cav3Entire metanotum
E. caipiraAtl/Cerstr/cav5Entire metanotum
E. bernadiiCerstr/cav4Entire metanotum
E. neomarmorataCercav3Entire metanotum
E. marmorataChistr/cav4Entire metanotum
E. longaChistr/cav4Entire metanotum

(Atl) Atlantic Forest; (Cer) Cerrado; (Chi) Chiquitano Dry Forest; (str) straminicolous; (cav) cavicolous; (0) less than 12mm; (1) 12-15mm; (2) 15-18mm; (3) 18–2.

An event-based analysis was performed using the Geographically-explicit Event Model (GEM) method [49]. This method attributes different kinds of events (vicariance, sympatry, founder event, and point sympatry) to internal nodes of the phylogenetic tree applying most of the spatial diversification models implemented in BioGeoBEARS [77], such as DIVA, DEC and DEC+J. However, instead of using pre-defined areas as operational units, GEM considers the geographical coordinates directly, as advocated by Hovenkamp [78], therefore, spurious biogeographical results due to inadequately defined areas may be avoided. The reconstruction cost is calculated as the event cost plus the number of distribution changes along the branches, and the distributions are based on the presence of taxa in distribution grids. The best reconstruction is the reconstruction with the lowest cost. The method was implemented using the software package EVS [49].

EVS was set according to the author’s recommendations: a raster grid with pixels of 1x1 degree was used, with a filling of 1. The cost of all four events (vicariance, sympatry, founder event, and point sympatry) was set to 1. In order to penalize ancestral distributions and avoid widespread ancestors, we used Z = 10 (Z defines the size of the ancestral ranges) as proposed by Arias [49]. The searches were performed using a flipping algorithm with 10 independent runs, each one with 10,000 replicates. The input files are provided in S2 and S3 Files. Permutation tests (10) were done to check the results, as recommended by Arias [49]. After the analysis, 3x3 squares were used in the cladogram to represent hypothetical distributions of the nodes and terminals in order to illustrate and interpret the sequence of events (Fig 1): vicariance events fragment the former ancestral wide distribution into two non-overlapping descendant branches; in sympatry events the two descendant branches have both the same distribution and equal to ancestral node; and in founder event one descendant branch is distributed equal to ancestral node while the other descendant branch has a restricted and non-overlapping distribution (grey sub-squares, distribution present; white sub-squares, distribution absent).

Habitats ancestral state reconstruction of Eidmanacris species based on the analysis without E. paramarmorata.
Fig 1

Habitats ancestral state reconstruction of Eidmanacris species based on the analysis without E. paramarmorata.

White branches, straminicolous; black branches, straminicolous/cavicolous; grey branches, cavicolous. Phylogeny reconstructed with 98 morphological characters and 28 species of Eidmanacris. Symbols on the nodes indicate the events of GEM method: black square, vicariance; white square, sympatry; white triangle, founder event. Squares 3x3 represent the total distribution of Eidmanacris in order to illustrate and interpret the sequence of events: grey sub-squares represent the presence and whit sub-squares the absence of species distribution. Colored terminal squares indicate the biomes where species came from: (1) Atlantic Forest, Parque Nacional de Itatiaia; (2) Cerrado, Parque Nacional Grande Sertão Veredas; (3) Chiquitano Dry Forest, Concepción, department of Santa Cruz (Reprinted from Mendivelso et al., 2013 under a CC BY license, with permission from PLOS One, original copyright 2013); (a) ancestral distribution of the genus Eidmanacris.

According to Arias [49], vicariance assignment is based on the sum of overlap of each descendant plus the cost of the vicariance, penalizing the overlap to indicate the difficulty of crossing barriers. The cost of assignment of sympatry is the cost of gained and lost pixels of both descendants in relation to the ancestral distribution plus the cost of the event. Point of sympatry assignment is based on the size of the range of overlapped pixels of descendants minus one, plus the event cost, the pixels outside the ancestral range being considered an extra cost. Founder event is based on the size of the range of overlapped pixels of descendants minus one, plus the cost of the event, the overlapped pixels with ancestral ranges being an extra cost (differently of point of sympatry). For detailed information about the method and its application, see Arias [49] and Ruiz et al. [79].

Morphology and ancestral state reconstruction

Two characters were selected to test the influence of habitats and biomes on Eidmanacris biology: body size and the metanotal glands in males. The body size is measured as the length between the fastigium and supra anal plate tip in dorsal view [46, 48]. Both characters were optimised on the phylogeny without E. paramarmorata.

The size is indicative of resource use [80]. Six non overlapping size categories separated by 3mm each were considered here. The arbitrary division has been chosen to better show size trends onto the topology (Table 1).

The male metanotal glands produce secretions and the females bite the metanotal structure during copulation. These secretions are considered as a male nuptial gift, but they also allow the male to increase mating duration, thus delaying the removal and chewing of the spermatophore by the female [81]. In Eidmanacris, the male glandular surface structure is delimited by an anterior median crest, two lateral projections, and a more or less extended area covered with bristles, where the secretions accumulate [48]. To take into account variation in male gland structure (= male investment in mating, or nuptial gift), we consider the surface of the bristle area, with two states: anterior portion of area covered by bristles, and entire area covered by bristles.

Ancestral state reconstructions were performed in Mesquite 3.6 using the parsimony criteria. The reconstructed characters were habitat (3 states), size (6 states) and bristle area in metanotal gland (2 states). The states of reconstructions are in Table 1.

Results

Phylogenetic analysis

The cladograms resulting from the analyses with character states and support values are presented in S19 and S20 Figs. The analysis with 37 terminals (including all the species of Eidmanacris except E. paramarmorata, plus all the outgroups) resulted in one most parsimonious tree (S19 Fig), with 240 steps, ci = 0.47 and ri = 0.78. The analysis with all 38 terminals resulted in 11 most parsimonious trees, with 240 steps, ci = 0.47 and ri = 0.78. The consensus tree is similar to the topology with 37 terminals, except for a polytomy in clade N2 (S4 File and S20 Fig) as expected with the “wildcard” E. paramarmorata. The discussion is based on analysis without E. paramarmorata due the absence of the polytomy of clade N2.

Eidmanacris is monophyletic, including the combinations Eidmanacris endophallica (de Mello) and E. minuta (de Mello), as proposed by Campos et al. [48]. All the clades are well supported, including clade N2 (S19 and S20 Figs).

Eidmanacris is supported by seven apomorphic and nine homoplastic changes (Table 2 and S4 Table) and the clade receives a high support value (Bremer support = 11). It is divided into two main clades: clade A with four species (E. fusca, E. melloi, E. endophallica, and E. minuta) and clade B with all the 25 remaining species. Clade B includes eight clades, labelled C to N (composed by N1 and N2) (S19 and S20 Figs), which are all characterised by distinct apomorphies (S5 Table). The clades C, E, I and K each include only one taxon, respectively E. bidentata, E. putuhra, E. papaveroi and E. fontanettiae.

Table 2
Exclusive synapomorphies for the clade Eidmanacris.
Character/StateState of character description
32(1)Presence of latero-posterior projections of the supra anal plate of male
64(1)Presence of anterior projection of the pseudepiphallic sclerite
67(1)No sclerotized connection between the PsP1 and PsP2
84(1)Ectophallic fold membranous
90(1)Endophallic apodeme crest-shaped
95(1)Medio-posterior projection of the endophallic sclerite elongate, through all the ectophallic fold
98(1)Presence of latero-posterior lobes of endophallic sclerite

Ancestral state reconstructions

Our analyses show that, with 5 steps, the ancestral habitat of Eidmanacris is straminicolous, with the species only foraging in the leaf litter. This pattern was modified in clade L becoming straminicolous/cavicolous, i.e. in addition to being active in the litter, species are also active in cavities and in the entrance of caves. The species E. speluncae (clade M), E. scopula and E. neomarmorata (clade N2) became cavicolous and a subsequent reversal to straminicoly occurs in clade M (E. multispinosa) (Table 1 and Fig 1).

The reconstruction of size, with 12 steps, displayed a general tendency of size increase along the phylogeny. The ancestral size of the genus is 12-15mm. The size increased in clade L (18-21mm) and the largest species are found in the clade N2 (larger than 21mm) following the same pattern of habitat modification (Table 1 and Fig 2).

Ancestral state reconstruction of Eidmanacris species for bristles on the metanotal gland (left) and size (right).
Fig 2

Ancestral state reconstruction of Eidmanacris species for bristles on the metanotal gland (left) and size (right).

A metanotal gland entirely covered by bristles is present in all the species in clade N (Fig 2B), while in the ancestral state of the genus bristles cover only the anterior region of metanotal gland (1 step) (Fig 2A). Species that are found in Cerrado are bigger, and their metanotum is entirely covered by bristles, while species that are only present in the Atlantic Forest, are smaller and have fewer bristles.

Biogeographical analysis

Applying the GEM method, four reconstructions were returned, two with a cost of 44 and two with a cost of 45. The 45 cost-reconstructions are similar to the 44 cost-reconstructions, except for the direction of a founder event between E. dissimilis and E. meridionalis (clade N1). We will discuss here only the 44 cost-reconstructions, which differ only in the node of clade G (a different direction in the founder event). The biogeographic scenario obtained implies three vicariances, five sympatries, and 17 founder events, but no point of sympatry (Fig 1). From the point of view of biome occupancy, the ancestral distribution of Eidmanacris was found to be the Atlantic Forest (Fig 1A). Species in clade N2 dispersed to Cerrado but two species, E. marmorata and E. longa, colonised the Bolivian Dry forest, and one species, E. caipira, returned to the Atlantic Forest. E. alboannulata and E. dissimilis (clade N1) are found both in Cerrado and Atlantic Forest. The optimisations of size and metanotal glands parallel the changes in the biome occupancy, with species becoming larger and the males having a more extended glandular structure as they colonise drier biomes (Fig 1).

The first vicariance event in the cladogram, in the basal node of Eidmanacris, divides the ancestral distribution of the genus (Fig 1A), between clade A and all the remaining species (clade B). A second vicariance event separates E. fusca of E. melloi + E. endophallica + E. minuta in clade A, with the sister species E. endophallica and E. minuta being sympatric. According to the event reconstruction, clade A, together with the species E. alboannulata and E. suassunai in clade N1, are present in the ancestral distribution area of the genus (Fig 1A).

As mentioned above, the two 44 cost-reconstructions only differ in the base of clade G. Clade G shows a founder event at its base, between E. tridentata and two sympatric species, E. simoesi + E. eliethae, but either E. tridentata or the sympatric E. simoesi and E. eliethae could occur in the ancestral distribution of the clade.

According to our results, the ancestral distribution of clade M is different from that of clade B because of the direction of the founder event that occurred in clade L. Clade M and N are separated by a founder event which defines the ancestral distribution of clade N as similar to that of clade B. Clade N is divided by a sympatry, i.e. clades N1 and N2 have the same ancestral distribution.

Clade N2 requires five founder events in seven nodes, and a vicariance event in the node separating E. corumbatai + E. gigas and E. caipira + E. bernardii + E. marmorata + E. longa + E. neomarmorata (Fig 3). All the extant distributions of clade N2 species differ from the ancestral distribution of the genus, but E. caipira is distributed in Atlantic Forest and Cerrado, like species in the clade N1. The other taxa of clade N2 have a western distribution, either in the Cerrado of central Brazil or in the Chiquitano Dry Forest from Bolivia (E. longa and E. marmorata).

Distribution map of Eidmanacris clades and species.
Fig 3

Distribution map of Eidmanacris clades and species.

Through the whole clade, sympatries are observed, either between species belonging to different clades (e.g., E. tridentata, E. papaveroi, and E. multispinosa) or between sister species (e.g., E. minuta and E. endophallica).

Discussion

Eidmanacris, with 29 species, is one of the most species rich genera of crickets from the Neotropical region [52], and the first in which a well-supported phylogeny has been constructed for all known species, following intensive field surveys. Our phylogenetic study confirms the monophyly of the genus, provided it includes Endophallusia de Mello, as proposed earlier [48]. The genus is subdivided into nine clades with strong support for all the nodes in the resultant cladogram (S19 and S20 Figs).

Ancestral distribution and diversification

The distribution of the Eidmanacris species covers the south, southeast, and center-west Brazil, as well as eastern Paraguay and eastern Bolivia, which corresponds to the Atlantic Forest, Cerrado, and the Chiquitano Dry Forest regions. Further cricket surveys in South America show that Eidmanacris is absent in the entire Amazonian region, as well as the southeastern and southern Atlantic Forest and southern South America. Using GEM, the ancestral distribution of the genus is hypothesised to have been the Atlantic Forest region, where Eidmanacris had its ancient diversification (Fig 1A). Only the more posterior clades have expanded the distribution of the genus further west, within the Brazilian central region to Bolivia.

This distribution pattern documents several biome changes, from the Atlantic Forest to the dryer biomes of Cerrado and the Chiquitano Dry Forest. Within Cerrado, Eidmanacris has been found only in gallery forests and never in other physiognomies, e.g. Campo limpo, Campo sujo, and Cerradão, although several collecting trips have been conducted to find them (LDC, pers. obs.). As mentioned above, gallery forests are known to act as remnant humid forests in dry and open regions and as dispersal corridors, allowing forest-dwelling species to reach the savanna areas [82]. This pattern of distribution is also documented for plants [18, 20, 21], butterflies [21, 83, 84], anurans [85], lizards [86], mammals [24, 87] and birds [26, 88].

Studies on mammals from Cerrado also suggest that gallery forests provide habitats for distinct communities throughout the Cerrado, and that these environments favor Cerrado colonisation by both Atlantic and Amazon forests species [24, 89].

According to our phylogenetic and biogeographic results, Eidmanacris originated in the Atlantic Forest, dispersed to Cerrado gallery forests, and from there colonised the Chiquitano Dry Forests in eastern Bolivia; some lineages returned to the Atlantic Forest, again through Cerrado gallery forests. However, Eidmanacris did not reach the Amazon Rainforest through Cerrado gallery forests. Even though Eidmanacris species from clade N are found in regions with environmental conditions different from Atlantic Forest, they are active in places that retain humidity, such as caves, rock crevices, and cavities near rivers that are mainly associated with gallery forests [47, 48, 50].

Two non-excludent hypotheses may explain how these species may have colonised the Cerrado and Bolivian Dry Forests:

    Eidmanacris was present in Cerrado and forced into a refuge in gallery forests by climatic variations: The climatic changes that favored vegetations with dryer characteristics during the Miocene and Pliocene could have contributed to a distributional retraction of Eidmanacris species to the gallery forests (and any forest refuges) of the Cerrado and Chiquitano Dry Forest [9092].

    Eidmanacris species were not present in the region of Cerrado and Chiquitano Dry Forest prior to the expansion of the xeric vegetation in the Quaternary period but may have colonised these biomes by dispersal through the gallery forests in more recent times. Palynological records showed that after the humidity increased again, an expansion of humid forests occurred in the Cerrado [93] allowing species of Atlantic Forest to expand into the Cerrado through gallery forests and stational forests [94].

Without a dated phylogeny, it is not possible to support any alternative. However, after xeric vegetations were established, the species may have used the gallery forests as dispersal corridors as well as refuges, which supports the idea that speciation and diversification in the Neotropics cannot be explained by only a single model of vicariance or climatic changes [29, 95]. Thus, the species would have colonised other regions in Cerrado and had the possibility to return to the Atlantic Forest, as did E. caipira. This is supported by the observations that two species in clade N1, E. dissimilis and E. alboannulata, are present both in the Atlantic Forest and in Cerrado, in geographically very approximate localities (Fig 3). In Cerrado, however, both species were collected only in cave entrances. The microclimate inside cavities and caves could maintain conditions of humidity and temperature close to that found in humid forests, and so allow the permanent presence of Eidmanacris in these dryer biomes.

The non-occurrence of Eidmanacris species in the Amazon Rainforest (LDC and LDG, pers. obs.; [52]) could fit the hypothesis of altitudinal segregation of fauna and flora species between Amazon and Atlantic Forests [26]. Amazonian elements occur at lower altitudes than the Atlantic elements, and more Atlantic Forest species occur in higher altitudes of Cerrado than in the Amazon area [26]. According to Oliveira-Filho & Fontes [96], Cerrado is more influenced by Atlantic Forest than by the Amazon Rainforest, mainly because the Atlantic Forest climate is more similar to that of Cerrado, with lower temperatures and pronounced dry seasons in the winter [97]. Furthermore, in the Amazon we found Phalangopsis Serville, which occupies the same habitat as Eidmanacris and does not occur in the Atlantic Forest, Cerrado, and Bolivian Dry Forest [52]. Similarly, Strinatia Chopard occurs in southern and southeastern Atlantic Forest areas where Eidmanacris does not occur, while both occupy the same habitat. These mutually exclusive distributional areas between closely-related phalangopsid genera with similar habits, Eidmanacris and Phalangopsis on one hand, Eidmanacris and Strinatia on the other, may represent cases of ecological vicariance [98], a question which will have to be further investigated to increase our understanding of cricket diversification patterns in the Neotropics.

A founder event occurs when a new population is established from a few individuals of an ancestral population [99, 100], leading to a founder effect, often documented in population genetics on speciation and island colonization, e.g. [101103]. Here, the number of founder events in the phylogeny suggests that Eidmanacris colonised new regions over time based on dispersals especially in clade N in which seven founder events occur in ten nodes. This scenario may have been supported by gallery forests within the Cerrado.

Among the three vicariance events recorded on Eidmanacris phylogeny, two are related to the presence of clade A species in mountain ranges and one is related to gallery forests. The species in clade A occur in forests at the slopes of Serra do Mar, Serra dos Orgãos and Serra Geral. According to Haffer [104], refuges in the Atlantic Forest region were probably located on mountain range slopes during dry periods maintained by orographic rainfalls. A basal vicariance event occurred between E. fusca, which is restricted to Serra Geral mountain range (state of Santa Catarina and northeastern of state of Rio Grande do Sul, Brazil) and its sister group E. melloi + E. endophallica + E. minuta that occurs in Serra do Mar, northeastern of São Paulo State and Serra dos Orgãos, in state of Rio de Janeiro state respectively.

The sympatric distributions observed in Eidmanacris occur either between species belonging to different clades or between sister species. Sympatry is usually accompanied by differences between species so that interspecific hybridization or overlap in resource use are avoided. Among the characters we studied, metanotal structures (related to mating behaviour) and size (related to resource use) can influence species interactions. The sister species E. endophallica and E. minuta (clade A), which occur together in the litter at the Parque Nacional da Serra dos Orgãos, state of Rio de Janeiro, are distinct morphologically: E. endophallica has no metanotal structures and a mean size of 16mm, while E. minuta is only 11mm in length (the smallest species of the genus) and has glandular structures in the metanotum. The morphological divergences of both species likely interfere directly on their mating behavior and resource use.

A similar situation exists between the species E. tridentata, E. papaveroi and E. multispinosa, which belong to clades G, I and M of Eidmanacris (Fig 1) and are sympatric in an Atlantic Forest fragment in Santa Teresa, state of Espírito Santo: the three species differ in size (see measurements in [46, 48]; Table 1) and in metanotal morphology (glandular structures absent in E. tridentata; vestigial in E. papaveroi; complete in E. multispinosa).

Habitat, resources, and behaviour changes

As a whole, the phylogeny illustrates a very coherent pattern of evolution for biome occupancy, habitat, resource use and mating behaviour in Eidmanacris. During the diversification of the genus, species became more straminicolous/cavicolous rather than only straminicolous, thus shifting toward more humid habitats, and species size increased significantly as they conquered dryer biomes; male nuptial gifts also increased, consequently allowing longer mating duration.

The changes in habitat and biome occupancy co-occur. Straminicolous species are present in the Atlantic Forest, while straminicolous/cavicolous or cavicolous species are found in the Atlantic Forest, but also Cerrado and the Chiquitano Dry Forest. For example, E. endophallica (clade A) is straminicolous in the Atlantic Forest; E. larvaeformis (clade M), E. bernardii (clade N2), and E. marmorata (clade N2) are straminicolous/cavicolous in the Atlantic Forest, Cerrado and the Chiquitano Dry Forest respectively; E. speluncae (clade M) and E. neomarmorata (clade N2) are cavicolous in the Atlantic Forest and Cerrado (Fig 1).

The changes in preferred habitats also follows the tendency of size increase within the genus (Table 1, Fig 2). Different resources are available in plain forest versus gallery forest, and in leaf litter versus cavities. Larger size could mean different/more available resources for the straminicolous/cavicolous species inside gallery forests or caves in Cerrado and the Chiquitano Dry Forest, compared to straminicolous species in humid forests. No data on Eidmanacris diet is available to explain these observations. Furthermore, larger size could also mean better resistance to desiccation [105, 106], through a greater capacity to store water [107] and reduced water loss with a smaller surface area/volume ratio [108].

With respect to mating behaviour, the modification in the metanotal gland is the most evident change in clade N, with bristles occupying the whole metanotum surface (Fig 2). While the female feeds on the nuptial gift, the male couples the genitalia in the female copulatory papilla and attaches the spermatophore tube [109]. More secretions in the bristles entertain the female for a longer time, increasing the male success to transfer its sperm [47, 110]. Clade N2 is also the most divergent in the genus: it is defined by seven apomorphies, six of which are related to reproduction, i.e. five related to male genitalia, and one related to aperture of the female copulatory papilla (see S5 Table). These modifications could interfere with the attachment of the male phallic complex in the female copulatory papilla, also modifying how the spermatophore tube is coupled in the female genitalia. The modifications of the metanotal gland, and male and female genitalia indicate that the mating behaviours of clade N2 species probably differ from the other species of the genus present in the Atlantic Forest.

Thus, Eidmanacris managed to occupy drier places by shifting habitats, becoming larger and changing mating behaviour.

Conservation issues

Conservation issues in the Atlantic Forests are becoming of crucial importance, as stated by several authors [33, 96, 111, 112]. The Atlantic Forest is one of 34 hotspots for conservation globally, and one of the most threatened [113]. According to Ribeiro et al. [114], only 11.73% of the original vegetation in the Atlantic Forest remains preserved, while more than 50% of animals and plants are endemic [113].

For harvestmen (Opiliones), twelve areas of endemism have been identified in the Atlantic Forest [115]. Eidmanacris, which has similar affinities to moist forests and humid microhabitats, occurs in seven of these areas, namely Bahia, Espírito Santo, Serra da Mantiqueira, Serra dos Orgãos, Southern Rio de Janeiro coast, Serra do Mar of São Paulo and Santa Catarina. It is absent in five endemism areas, three areas poorly sampled for crickets (Pernambuco, Serra do Espinhaço, Serra da Bocaina), and two areas further south, which are occupied by the genus Strinatia (Southern of São Paulo, Parana, see above).

The similar patterns of occurrence of Eidmanacris and harvestmen in the Atlantic Forest, highlight the need for biogeographical studies using non model organisms and consequently further enhancing the conservation policies of this biome. The pattern of distribution of Eidmanacris species, although documented by few localities, increases the originality of the fauna in the Atlantic Forest and supports the same core areas as for harvestmen.

Similarly, gallery forests play an important role for species from the Atlantic Forest in Cerrado, as illustrated by Eidmanacris species in clade N, but also for other species of invertebrates [21, 83], vertebrates [24, 26, 29] and plants [20, 23, 116]. These forests increase the biodiversity of the Cerrado biome, acting as a secondary habitat for organisms [29, 82]. Nevertheless, Cerrado and its gallery forests are threatened by human activities, such as cattle raising and agricultural expansion in the central region of Brazil [117119]. In addition to their function as dispersal corridors and their general ecological role, the gallery forests play a key role in the preservation of the environment, stabilising soils and protecting rivers and its headwaters [120].

Conclusion

The phylogenetic patterns of Eidmanacris recovered here, with congruent changes in distribution, habitat, resource use and mating behaviour, shows the fragile balance that may allow species survival and dispersal. This case study clearly highlights the need to protect what remains of these Neotropical biomes, and their interactions in terms of taxa survival and evolution. In the future, additional geographical and ecological data per species will allow niche reconstructions [121, 122], complementing direct observations on habitats and nychthemeral activity, in order to further document the criteria of distribution of Eidmanacris species. Also, additional behavioural observations, the relations between biome occupancy, habitat, resource use to further analyses could reinforce the hypotheses on the motor of the evolution of Eidmanacris and future conservation issues.

Acknowledgements

We thank Dr. Eliana Marques Cancello (MZSP), Dr. Freddy Bravo (MZUEFS), Dr. Francisco de Assis Ganeo de Mello (UBTU) for the loan the material. LDC thanks to the entomology collection manager of the ANSP, Jason D. Weintraub, for receiving him as a visitant researcher. Part of this paper was developed during LDC internship in the MNHN in 2016 and 2019/2020. LDG and SSN are both co-last authors of this article.

References

AAntonelli, AZizka, FAntunes, RScharn, CDBacon, DSilvestro, et al Amazonia is the primary source of Neotropical biodiversity. Proc Natl Acad Sci USA. 2018;115: 60346039. 10.1073/pnas.1713819115

J.Haffer Speciation in Amazonian Forest Birds. Science. 1969;165: 131135. 10.1126/science.165.3889.131

GTPrance, TELovejoy. Amazon. Oxford: Pergamon Press; 1985.

AAntonelli, ISanmartín. Why are there so many plant species in the Neotropics? Taxon. 2011;60: 403414. 10.1002/tax.602010

NMyers, RAMittermeier, CGMittermeier, GABDa Fonseca, JKent. Biodiversity hotspots for conservation priorities. Nature. 2000;403: 853 10.1038/35002501

W.Mantovani Patrimônio Ambiental Brasileiro. São Paulo: EDUSP 2003.

DEPrado, PEGibbs. Patterns of Species Distributions in the Dry Seasonal Forests of South America. Ann Missouri Bot Gard. 1993;80: 902927. 10.2307/2399937

TAParkerIII, AHGentry, RBFoster, LHEmmons, JVRemsenJr. The lowland dry forests of Santa Cruz, Bolivia: a global conservation priority. Cornell: Conservation International; 1993.

PGMurphy, AELugo. Ecology of tropical dry forest. Annu Rev Ecol Evol Syst. 1986;17: 6788.

10 

EDinerstein, DMOlson, DJGraham, ALWebster, SAPrimm, MPBookbinder, et al A conservation assessment of the terrestrial ecoregions of Latin America and the Caribbean. The World Bank; 1995.

11 

ALKlein. Eugen Warming e o cerrado brasileiro: um século depois. São Paulo: Editora da UNESP; 2002.

12 

ANAb’Saber. Os domínios morfoclimáticos na América do Sul. Geomofol. 1977;52: 121.

13 

DMOlson, EDinerstein, EDWikramanayake, NDBurgess, GVNPowell, ECUnderwood, et al Terrestrial Ecoregions of the World: A New Map of Life on Earth. Bioscience. 2001;51: 933938.

14 

EDinerstein, DOlson, AJoshi, CVynne, NDBurgess, EWikramanayake, et al An Ecoregion-Based Approach to Protecting Half the Terrestrial Realm. Bioscience. 2017;67: 534545. 10.1093/biosci/bix014

15 

PAFurley, JARatter. Soil resources and plant communities of the central Brazilian cerrado and their development. Journal of Biogeography. 1988; 97108.

16 

JFRibeiro, BMTWalter. Fitofisionomias do bioma Cerrado. Embrapa Cerrados-Capítulo em livro científico (ALICE) 1998.

17 

PSOliveira, RJMarquis, editors. The cerrados of Brazil: ecology and natural history of a neotropical savanna New York: Columbia University Press; 2002.

18 

JARatter, PWRichards, GArgent, DRGifford. Observations on the vegetation of northeastern Mato Grosso: I. The woody vegetation types of the Xavantina-Cachimbo Expedition area. Philos Trans R Soc Lond. 1973;266: 449492.

19 

CTRizzini. Tratado de fitogeografia do Brasil: aspectos florísticos e estruturais. São Paulo: HUCITEC/EDUSP; 1979.

20 

SAMori, BMBoom, GTPrance. Distribution Patterns and Conservation of Eastern Brazilian Coastal Forest Tree Species. Brittonia. 1981;33: 233 10.2307/2806330

21 

KSBrownJr. Biogeography and evolution of Neotropical butterflies In: TCWhitmore, GTPrance, editors. Biogeography and Quaternary history in tropical America. Oxford; 1987 pp. 66104.

22 

AFouquet, DLoebmann, SCastroviejo-Fisher, JMPadial, VGOrrico, MLLyra, et al From Amazonia to the Atlantic forest: Molecular phylogeny of Phyzelaphryninae frogs reveals unexpected diversity and a striking biogeographic pattern emphasizing conservation challenges. Mol Phylogenet Evol. 2012;65: 547561. 10.1016/j.ympev.2012.07.012

23 

TSobral-Souza, MSLima-Ribeiro, VNSolferini. Biogeography of Neotropical Rainforests: past connections between Amazon and Atlantic Forest detected by ecological niche modeling. Evol Ecol. 2015;29: 643655. 10.1007/s10682-015-9780-9

24 

KHRedford, GABda Fonseca. The Role of Gallery Forests in th Zoogeography of the Cerrado’s Non-volant Mammalian Fauna. Biotropica. 1986;18: 126135. 10.2307/2388755

25 

ATOliveira-Filho, JARatter. A study of the origin of central Brazilian forests by the analysis of plant species distribution patterns. Edinb J Bot. 1995;52: 141194.

26 

JMC daSilva Distribution of Amazonian and Atlantic birds in gallery forests of the Cerrado region, South America. Ornitol Neotrop. 1996;7: 118.

27 

IRDiniz, KKitayama. Seasonality of vespid species (Hymenoptera: Vespidae) in a central Brazilian cerrado. Rev Biol Trop. 1998;46: 109114.

28 

P.Müller Dispersal centres of terrestrial vertebrates in the neotropic realm: A study in the evolution of the neotropical biota and its native landscapes. 1st ed. Haugue: Springer; 1973.

29 

LPCosta. The historical bridge between the Amazon and the Atlantic Forest of Brazil: a study of molecular phylogeography with small mammals. J Biogeogr. 2003;30: 7186. 10.1046/j.1365-2699.2003.00792.x

30 

DFStotz, JWFitzpatrick, TAParkerIII, DKMoskovits. Neotropical Birds: Ecology and Conservation. Chicago: University of Chicago Press; 1996.

31 

WEDuellman. Patterns of Distribution of Amphibians: A Global Perspective. Baltimore: JHU Press; 1999.

32 

WLROWillian, IBSantos. Threatened endemic mammals of the Atlantic forest region of south-east Brazil. Michigan: Jersey Wildlife Preservation Trust; 1991.

33 

RSRidgely, JAGwynne, GTudor, MArgel. Wildlife Conservation Society Birds of Brazil: The Atlantic Forest of Southeast Brazil, including São Paulo and Rio de Janeiro. Cornell University Press; 2016.

34 

CERLehmann, SAArchibald, WAHoffmann, WJBond. Deciphering the distribution of the savanna biome. New Phytol. 2011;191: 197209. 10.1111/j.1469-8137.2011.03689.x

35 

MJDonoghue, EJEdwards. Biome Shifts and Niche Evolution in Plants. Annu Rev Ecol Evol Syst. 2014;45: 547572. 10.1146/annurev-ecolsys-120213-091905

36 

L.Desutter Structure et évolution du complexe phallique des Gryllidea (Orthoptères) et classification des genres Néotropicaux de Grylloidea. Primière Partie. Ann Soc entomol Fr. 1987;23: 213239.

37 

L.Desutter Structure et évolution du complexe phallique des Gryllidea (Orthoptères) et classification des genres Nèotropicaux de Grylloidea.-Deuxième Partie. Ann Soc entomol Fr. 1988;24: 342373.

38 

L.Desutter-Grandcolas Toward the knowledge of the evolutionary biology of Phalangopsid crickets (Orthoptera: Grylloidea: Phalangopsidae): Data, questions and Evolutionary scenarios. J Orthoptera Res. 1995;4: 163175.

39 

TRobillard, LDesutter-Grandcolas. A revision of Neotropical Eneopterinae crickets (Orthoptera, Grylloidea, Eneopteridae) with phylogenetic discussion. Insect Syst Evol. 2005;35: 411435.

40 

LD deCampos, PGBSouza-Dias, SSNihei New Brazilian Tafaliscina increase the diversity of this Neotropical cricket clade (Orthoptera: Grylloidea: Gryllidae: Oecanthinae: Paroecanthini). Zoosystema. 2020;42: 331 10.5252/zoosystema2020v42a19

41 

LDesutter-Grandcolas, TRobillard. Acoustic evolution in crickets: Need for phylogenetic study and a reappraisal of signal effectiveness. An Acad Bras Ciênc. 2004;76: 301315. 10.1590/s0001-37652004000200019

42 

JAnso, LBarrabé, LDesutter-Grandcolas, HJourdan, PGrandcolas, JDong, et al Old Lineage on an Old Island: Pixibinthus, a New Cricket Genus Endemic to New Caledonia Shed Light on Gryllid Diversification in a Hotspot of Biodiversity. DFontaneto, editor. PLOS ONE. 2016;11: 134. 10.1371/journal.pone.0150920

43 

AGasc, JAnso, JSueur, HJourdan, LDesutter-Grandcolas. Cricket calling communities as an indicator of the invasive ant Wasmannia auropunctata in an insular biodiversity hotspot. Biol Invasions. 2018;20: 10991111. 10.1007/s10530-017-1612-0

44 

AHemp, CHemp. Broken bridges: The isolation of Kilimanjaro’s ecosystem. Glob Change Biol. 2018;24: 34993507. 10.1111/gcb.14078

45 

CHemp, SMKuechler, SKehl, JWWaegele, AHemp. Dispersal and speciation of East African mountain bushcrickets of the genus Phlesirtes (Orthoptera: Tettigoniidae: Conocephalinae, Conocephalini) related to climate fluctuations. Syst Entomol. 2019;44: 408417.

46 

Souza-Dias PGBCampos LD de, SSNihei, F de AGDe Mello. New species of Eidmanacris Chopard, 1956 from Brazil (Orthoptera: Phalangopsidae: Luzarinae). Zootaxa. 2015;4018: 228248. 10.11646/zootaxa.4018.2.4

47 

PGBSouza-Dias, LD deCampos, SSNihei. Two New Species of Eidmanacris (Orthoptera: Grylloidea: Phalangopsidae) from the Atlantic Forest of São Paulo State, Brazil. Fla Entomol. 2015;98: 547555. 10.1653/024.098.0223

48 

Souza-Dias PGBCampos LD de, SSNihei. Taxonomic review of Eidmanacris Chopard, 1956 (Orthoptera: Grylloidea: Phalangopsidae). Zootaxa. 2017;4321: 193. 10.11646/zootaxa.4321.1.1

49 

JSArias. An event model for phylogenetic biogeography using explicitly geographical ranges. J Biogeogr. 2017;44: 22252235. 10.1111/jbi.13024

50 

L.Desutter-Grandcolas Le genre Eidmanacris Chopard, 1956 (Orthoptera. Grylloidea, Phalangopsidae, Luzarinae): habitat, répartition et espèces nouvelles. Bull Mus Natl Hist Nat. 1995;16: 453474.

51 

AVGorochov. Classification of the Phalangopsinae subfamily group, and new taxa from the subfamilies Phalangopsinae and Phaloriinae (Orthoptera: Gryllidae). Zoosyst Ross. 2014;23: 788.

52 

Cigliano MM, Braun H, Eades DC, Otte D. Orthoptera Species File. Version 5.0/5.0 [WWW document]. 2020 [cited 12 Sep 2020]. Available: http://orthoptera.speciesfile.org

53 

KCNixon, JMCarpenter. On outgroups. Cladistics. 1993; 413426.

54 

PGBde Souza-Dias. Análise cladística e morfologia do complexo fálico de Phalangopside, com ênfase em Luzarinae (Orthoptera, Ensifera, Grylloidea). Instituto de Biociências da Universiadade de São Paulo. 2015;1: 72.

55 

Desutter L. Étude phylogénétique, biogéografique et écologique des Grylloidea néotropicaux (Insectes, Orthoptères). PhD Thesis, Université de Paris-Sud—Centre d’Orsay. 1990. p. 348.

56 

LDesutter-Grandcolas. Phylogeny and the evolution of acoustic communication in extant Ensifera (Insecta, Orthoptera). Zool Scr. 2003;32: 525561. 10.1046/j.1463-6409.2003.00142.x

57 

RAPrado, CSFontanetti. Metanotal gland of the genus Eidmanacris (Grylloidea, Phalangopsidae): taxonomic importance. Iheringia, Sér Zool. 2005;95: 8387. 10.1590/S0073-47212005000100012

58 

PCSereno. Logical basis for morphological characters in phylogenetics. Cladistics. 2007;23: 565587. 10.1111/j.1096-0031.2007.00161.x

59 

WPMaddison, DRMaddison. Mesquite: a modular system for evolutionary analysis. 2019 Available: http://www.mesquiteproject.org

60 

L.Desutter-Grandcolas Les Phalangopsidae de Guyane française (Orthoptères, Grylloidea): systématique, élements de phylogénie et de biologie. Bull Soc entomol Fr. 1992;14: 93177.

61 

L.Desutter-Grandcolas The cricket fauna of Chiapanecan caves (Mexico): systematics, phylogeny and the evolution of troglobitic life (Orthoptera, Grylloidea, Phalangopsidae, Luzarinae). Int J Speleol. 1993;22: 182. 10.5038/1827-806X.22.1.1

62 

PAGoloboff, SACatalano. TNT version 1.5, including a full implementation of phylogenetic morphometrics. Cladistics. 2016;32: 221238. 10.1111/cla.12160

63 

KCNixon. WinClada ver. 1.00. 08. Ithaca, New York: Published by the author; 2002.

64 

AGKluge, JSFarris. Quantitative Phyletics and the Evolution of Anurans. Syst Zool. 1969;18: 132. 10.2307/2412407

65 

JSFarris. The retention index and the rescaled consistency index. Cladistics. 1989;5: 417419. 10.1111/j.1096-0031.1989.tb00573.x

66 

JSFarris. Methods for Computing Wagner Trees. Syst Biol. 1970;19: 8392. 10.1093/sysbio/19.1.83

67 

DLSwofford, WPMaddison. Reconstructing ancestral character states under Wagner parsimony. Math Biosci. 1987;87: 199229. 10.1016/0025-5564(87)90074-5

68 

I.Agnarsson Morphological phylogeny of cobweb spiders and their relatives (Araneae, Araneoidea, Theridiidae). Zool J Linn Soc. 2004;141: 447626. 10.1111/j.1096-3642.2004.00120.x

69 

IAgnarsson, JAMiller. Is ACCTRAN better than DELTRAN? Cladistics. 2008;24: 10321038. 10.1111/j.1096-0031.2008.00229.x

70 

K.Bremer Branch support and tree stability. Cladistics. 1994;10: 295304. 10.1111/j.1096-0031.1994.tb00179.x

71 

KCNixon, QDWheeler. Extinction and the origin of species In: MJNovacek, QDWheeler, editors. Extinction and phylogeny. New York: Columbia University Press; 1992 pp. 119143.

72 

L.Chopard Description d’un Grylllide nouveau du Brésil. Bull Soc entomol Fr. 1938;43: 159160.

73 

AMesa, F de AGDe Mello. A new species of the genus Eidmanacris (Orthoptera, Grylloidea, Phalangopsidae). Rev Bras Entomol. 1985;29: 199204.

74 

AMesa, CFSperber, PCGarcia. Two New Species of the Cricket Genus Eidmanacris and a New Combination Name for a Third Species (Orthoptera, Grylloidea, Phalangopsidae). Trans Am Entomol Soc. 1998;124: 4357.

75 

MPBolfarini. A new species of Eidmanacris Chopard, 1956, with notes on its distribution in Brazilian caves (Phalangopsidae, Luzarinae). Zootaxa. 2016;4083: 526532. 10.11646/zootaxa.4083.4.4

76 

QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project 2020 Available: https://qgis.org

77 

NJMatzke. BioGeoBEARS: BioGeography with Bayesian (and likelihood) evolutionary analysis in R Scripts. 2013.

78 

P.Hovenkamp Vicariance Events, not Areas, Should be Used in Biogeographical Analysis. Cladistics. 1997;13: 6779. 10.1111/j.1096-0031.1997.tb00241.x

79 

MSRuiz, SNPortelli, TNHibbard, ASQuinteros. Discovering the biogeographic history using predefined areas and explicit geographical data in the South American Liolaemus elongatus group (Iguania: Liolaemidae). Herpetol J. 2020;30: 5368. 10.33256/hj30.2.5368

80 

D.Tilman Resource competition and community structure. Princenton: Princeton university press; 1982

81 

DTGwynne. A secondary copulatory structure in a female insect: A clasp for a nuptial meal? Naturwissenschaften. 2002;89: 125127. 10.1007/s00114-002-0298-y

82 

JFEisenberg, KHRedford. Mammals of the Neotropics, Volume 3: Ecuador, Bolivia, Brazil. Chicago: University of Chicago Press; 1999.

83 

KSBrownJr, OHHMielke. Lepidoptera of the Central Brazil Plateau. I. Preliminary list of Rhopalocera (continued): Lycaenidae, Pieridae, Papilionidae, Hesperiidae. J Lepid Soc. 1967;21: 145168.

84 

HSRCabette, JRSouza, YShimano, LJuen. Effects of changes in the riparian forest on the butterfly community (Insecta: Lepidoptera) in Cerrado areas. Rev Bras Entomol. 2017;61: 4350. 10.1016/j.rbe.2016.10.004

85 

TSVasconcelos, JLPDoro. Assessing how habitat loss restricts the geographic range of Neotropical anurans. Ecol Res. 2016;31: 913921. 10.1007/s11284-016-1401-8

86 

CNogueira, PHValdujo, FGRFrança. Habitat variation and lizard diversity in a Cerrado area of Central Brazil. Stud Neotrop Fauna Environ. 2005;40: 105112. 10.1080/01650520500129901

87 

IRBishop. An annotated list of caviomorph rodents collected in north-eastern Mato Grosso Brazil. Mammalia. 1974;38: 489502.

88 

EOWillis. Zoogeographical origins of eastern Brazilian birds. Ornitol Neotrop. 1992;3: 115.

89 

MAJohnson, PMSaraiva, DCoelho. The role of gallery forests in the distribution of cerrado mammals. Rev Bras Biol. 1999;59: 421427. 10.1590/S0034-71081999000300006

90 

MMCole. Cerrado, Caatinga and Pantanal: The Distribution and Origin of the Savanna Vegetation of Brazil. Geogr J. 1960;126: 168179. 10.2307/1793957

91 

AGFernandes. Fitogeografia brasileira. 2nd ed. Fortaleza: Multigraf; 2000.

92 

MHOPinheiro, RMonteiro. Contribution to the discussions on the origin of the cerrado biome: Brazilian savanna. Braz J Biol. 2010;70: 95102. 10.1590/s1519-69842010000100013

93 

M-PLedru. Late Quaternary Environmental and Climatic Changes in Central Brazil. Quat Res. 1993;39: 9098. 10.1006/qres.1993.1011

94 

BBMéio, CVFreitas, LJatobá, MEFSilva, JFRibeiro, RPBHenriques. Influência da flora das florestas Amazônica e Atlântica na vegetação do cerrado sensu stricto. Rev Bras Bot. 2003;26: 437444. 10.1590/S0100-84042003000400002

95 

SSNihei, CJBDCarvalho. Systematics and biogeography of Polietina Schnabl & Dziedzicki (Diptera, Muscidae): Neotropical area relationships and Amazonia as a composite area. Syst Entomol. 2007;32: 477501. 10.1111/j.1365-3113.2006.00376.x

96 

ATOliveira-Filho, MALFontes. Patterns of Floristic Differentiation among Atlantic Forests in Southeastern Brazil and the Influence of Climate. Biotropica. 2000;32: 793810. 10.1111/j.1744-7429.2000.tb00619.x

97 

E.Nimer Climatologia do Brasil. Rio de Janeiro: Fundação Instituto Brasileiro de Geografia e Estatística; 1979.

98 

JACoyne, HAOrr. Speciation Oxford, New York: Oxford University Press; 2004.

99 

E.Mayr Change of genetic environment and evolution. In: JHuxley, ACHardy, EBFord, editors. Evolution as a process. London: Macmilan; 1954 pp. 157180.

100 

ARTempleton. The reality and importance of founder speciation in evolution. BioEssays. 2008;30: 470479. 10.1002/bies.20745

101 

SMClegg, SMDegnan, JKikkawa, CMoritz, AEstoup, IPOwens. Genetic consequences of sequential founder events by an island-colonizing bird. Proc Natl Acad Sci USA. 2002;99: 81278132. 10.1073/pnas.102583399

102 

MJCardoso, MDBEldridge, MOakwood, BRankmore, WBSherwin, KBFirestone. Effects of founder events on the genetic variation of translocated island populations: Implications for conservation management of the northern quoll. Conserv Genet. 2009;10: 17191733. 10.1007/s10592-008-9774-z

103 

NJMatzke. Model selection in historical biogeography reveals that founder-event speciation is a crucial process in island clades. Syst Biol. 2014;63: 951970. 10.1093/sysbio/syu056

104 

J.Haffer Biogeography of neotropical birds In: TCWhitmore, GTPrance, editors. Biogeography and Quaternary history in tropical America. Oxford: Oxford Science Publications; 1987 pp. 105150.

105 

DWWhitman. The significance of body size in the Orthoptera: a review. J Orthoptera Res. 2008;17: 117134. 10.1665/1082-6467-17.2.117

106 

WEWinterhalter, TAMousseau. The strength of temperature-mediated selection on body size in a wild insect population. J Orthoptera Res. 2008;17: 347351. 10.1665/1082-6467-17.2.347

107 

SLChown. Desiccation resistance in six sub-Antarctic weevils (Coleoptera: Curculionidae): humidity as an abiotic factor influencing assemblage structure. Funct Ecol. 1993;7: 318325.

108 

SLChown, PPistorius, CHScholtz. Morphological correlates of flightlessness in southern African Scarabaeinae (Coleoptera: Scarabaeidae): Testing a condition of the water-conservation hypothesis. Can J Zool. 1998;76: 11231133. 10.1139/z98-036

109 

RDAlexander, DOtte. The Evolution of Genitalia and Mating Behavior in Crickets (Gryllidae) and other Orthoptera. Misc Publ Mus Zool Univ Mich 1967; 162.

110 

R.Prado Reproductive behavior of Eidmanacris corumbatai Garcia (Orthoptera: Phalangopsidae). Neotrop Entomol. 2006;35: 452457. 10.1590/s1519-566x2006000400005

111 

MÁlvarez-Presas, ASánchez-Gracia, FCarbayo, JRozas, MRiutort. Insights into the origin and distribution of biodiversity in the Brazilian Atlantic forest hot spot: A statistical phylogeographic study using a low-dispersal organism. Heredity. 2014;112: 656665. 10.1038/hdy.2014.3

112 

MBDa Silva, RPinto-da-Rocha, JJMorrone. Historical relationships of areas of endemism of the Brazilian Atlantic rain forest: a cladistic biogeographic analysis of harvestman taxa (Arachnida: Opiliones). Curr Zool. 2016;63: 525535. 10.1093/cz/zow092

113 

RAMittermeier, GABDa Fonseca, ABRylands, KBrandon. A Brief History of Biodiversity Conservation in Brazil. Conserv Biol. 2005;19: 601607. 10.1111/j.1523-1739.2005.00709.x

114 

MCRibeiro, JPMetzger, ACMartensen, FJPonzoni, MMHirota. The Brazilian Atlantic Forest: How much is left, and how is the remaining forest distributed? Implications for conservation. Biol Conserv. 2009;142: 11411153. 10.1016/j.biocon.2009.02.021

115 

MBDa Silva, RPinto-da-Rocha, AMDe Souza. A protocol for the delimitation of areas of endemism and the historical regionalization of the Brazilian Atlantic Rain Forest using harvestmen distribution data. Cladistics. 2015;31: 692705. 10.1111/cla.12121

116 

TCWhitmore, GTPrance. Biogeography and Quaternary History in Tropical America. Oxford: Clarendon Press; 1987.

117 

PAFurley. The nature and diversity of neotropical savanna vegetation with particular reference to the Brazilian cerrados. Glob Ecol Biogeogr. 1999;8: 223241. 10.1046/j.1466-822X.1999.00142.x

118 

JMCardoso Da Silva, JMBates. Biogeographic Patterns and Conservation in the South American Cerrado: A Tropical Savanna Hotspot. Bioscience. 2002;52: 225 10.1641/0006-3568(2002)052[0225:BPACIT]2.0.CO;2

119 

FMVCarvalho, PDe Marco, LGFerreira. The Cerrado into-pieces: Habitat fragmentation as a function of landscape use in the savannas of central Brazil. Biol Conserv. 2009;142: 13921403. 10.1016/j.biocon.2009.01.031

120 

HFVogel, CHZawadzki, RMetri. Florestas ripárias: importância e principais ameaças. SaBios. 2009;4: 2430.

121 

SJPhillips, RPAnderson, RESchapire. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190: 231259. 10.1016/j.ecolmodel.2005.03.026

122 

SJPhillips, MDudík. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008;31: 161175.