Outstanding diversity and microendemism in a clade of rare Atlantic Forest montane frogs
Graphical abstract
Introduction
Mountains affect the origin and maintenance of biodiversity directly and indirectly, and are often high biodiverse areas themselves (Hoorn et al., 2013). Montane systems in the tropics sustain even larger number of species and striking degrees of endemism at higher elevations (Cadena et al., 2012; Merckx et al., 2015). Mountain taxa diversification is often explained by neutral and allopatric models: populations become isolated in mountaintops and differentiate after enduring intense genetic drift (e.g. Firkowski et al., 2016), albeit non-neutral and sympatric diversification across elevation gradients have also been explored (e.g. Patton and Smith, 1992).
The Brazilian Atlantic Forest is a hyperdiverse biome distributed across the tropical and subtropical Atlantic margin of the South American continent. In eastern Brazil, topography is especially complex and landscapes vary in elevation from sea level to nearly 3,000 m above sea level (Pico da Bandeira, 2892 m a.s.l.). Diversification hypotheses invoke isolation of populations following Plio-Pleistocene climate-related habitat fragmentation (Carnaval et al., 2009) or populations being bounded by geographic barriers (Thome et al., 2014). In this context, whereas the complex relief is permeated by neotectonic faults and lineaments (Saadi et al., 2002), a preference for the cold would cause montane taxa to respond to past climate change in a particular way, becoming more widespread during glacial times and restrained to mountaintops at interglacials (Amaro et al., 2012, Carnaval et al., 2014, Rodrigues et al., 2009).
Atlantic Forest amphibians constitute one of the most diverse vertebrate groups globally; a recent account listing an impressive number of 625 species likely represents an underestimation as new taxa are frequently described (Rossa-Feres et al., 2017). Because of particular life history traits that enhance phylogeographic signal retention, such as low vagility and habitat specificity (Avise, 2000), Atlantic Forest amphibians have been used as models to elucidate the origins of the biome's biodiversity, attributed to a combination of distinct evolutionary mechanisms and deep histories (Turchetto-Zolet et al., 2013). Particularly for amphibians, diversity is sustained by a combination of large altitudinal and latitudinal ranges with a complex relief creating moisture gradients and climatic conditions that vary at local scale (da Silva et al., 2012, Haddad and Prado, 2005).
One of the most emblematic anuran families in the Atlantic Forest is the Brachycephalidae, comprising ground-dwelling, direct-developing frogs in the genera Brachycephalus and Ischnocnema (Hedges et al., 2008). Both show high species diversity and extreme endemism levels. About half of the Brachycephalidae species were described in the last decade (Frost, 2019), and the availability of molecular data has played a fundamental role in solving phylogenetic uncertainty within complexes of cryptic species. A relevant aspect of the ecology of the family is the preference for montane habitats and very restricted ranges; often, its species span a single mountaintop. This trait is prevalent in Brachycephalus but less frequently reported for Ischnocnema, although recent molecular studies revealed that several Ischnocnema species correspond, in fact, to complexes of multiple, highly divergent lineages with restricted distributions (e.g. Gehara et al., 2017, Taucce et al., 2018a), including cases of microendemism (I. guenteri, Gehara et al., 2013).
Ischnocnema currently comprises 38 described species distributed in five species series, as supported by recent phylogenetic studies; they are the I. guentheri, I. lactea, I. parva, I. venancioi, and I. verrucosa species series (Canedo and Haddad, 2012, Taucce et al., 2018b). The Ischnocnema lactea species series is at the base of the Ischnocnema tree and comprises 11 species (Canedo and Haddad, 2012, Taucce et al., 2019). Within this series, a particular clade containing I. lactea (Miranda-Ribeiro, 1923) and I. holti (Cochran, 1948) was highlighted as a possible species complex, given that I. holti was recovered paraphyletic in relation to I. lactea (Canedo and Haddad, 2012, Taucce et al., 2018a, Taucce et al., 2018b). The two species are associated to montane Atlantic Forest physiognomies in Southeastern Brazil, being notably rare in the field and poorly represented in museum collections. Ischnocnema lactea inhabits the Serra do Mar mountain range at moderate elevations (700–1120 m a.s.l., Silva-Soares et al., 2018), and was recently redescribed by Silva-Soares et al. (2018). Ischnocnema holti was redescribed by Targino and de Carvalho-e-Silva (2008), who restricted its occurrence to a single high elevation site in the Serra da Mantiqueira mountain range (Itamonte, 2000–2400 m a.s.l.). Costa et al. (2008) expanded this geographic distribution to Serra dos Órgãos, and specimens of dubious identification (morphologically similar to I. holti) were also found in other high elevation areas of southeastern Brazil (see Table 1, Fig. 1). Little is known about natural history of the populations assigned to the two species names, with individuals of I. lactea being reported to reproduce during cold fronts in the winter season, under relatively low temperatures (Silva-Soares et al., 2018).
In this study, we prospect tissue samples from major museum collections in southeastern Brazil to access the genetic diversity of the rare I. holti – I. lactea species complex. We delimit lineages by combining traditional frequency and tree-based methods with coalescent methods under a multilocus approach. We use a molecular clock to define a time window for diversification in the complex, and explore aspects of its evolutionary history by contrasting lineage distributions to maps of geographic barriers and past putative ranges obtained with ecological niche modeling (ENM). Finally, we project the I. holti – I. lactea species complex climatic niche to the future and discuss its conservation under a global climate change scenario.
Section snippets
Taxon and gene sampling
We gathered tissue samples from vouchers previously identified as Ischnocnema holti, I. cf. holti, I. aff. holti, I. lactea, I. cf. lactea, I. aff. lactea, I. gr. lactea, Ischnocnema sp. (gr. lactea), and Ischnocnema sp., deposited in six main tissue collections in Brazil. This strategy resulted in 37 tissue samples from 18 localities (Fig. 1, Table S1). We also sampled five species to be used as outgroups in phylogenetic analyses, following results in Canedo and Haddad, 2012, Taucce et al.,
Genetic structure and phylogenetic inferences
Structure analyses suggest the existence of seven unique genetic clusters within the I. holti – I. lactea species clade. The plot with mean LnP(K) across K values shows an inflexion point at K = 7, with similar mean LnP(K) values above that. The Evanno criterion ΔK peaks twice, at K = 2 and K = 7 (Fig. 2A), indicating hierarchical genetic structure: population limits at K = 2 are also recovered at K = 7 (Fig. 2, Table 1). At K = 2 all samples are attributed with high q values, being clustered
Discussion
In this study we combined assignment analyses with phylogenetic inference to prospect lineages within the I. holti – I. lactea species complex, validated those lineages under the multispecies coalescent model, and used ENM to explore its past and future putative range. We detected nine lineages that were fully supported as species, but their phylogenetic relationships were not completely resolved. Extensive suitable areas at the LGM, along with phylogenetic uncertainty, are consistent with
CRediT authorship contribution statement
Maria Tereza C. Thomé: Writing - original draft, Formal analysis. Mariana L. Lyra: Data curation, Writing - review & editing. Priscila Lemes: Formal analysis, Writing - review & editing. Laryssa S. Teixeira: Data curation. Ana Carolina Carnaval: Writing - review & editing. Célio F.B. Haddad: Resources, Funding acquisition, Data curation, Writing - review & editing. Clarissa Canedo: Conceptualization, Supervision, Writing - review & editing.
Acknowledgments
We thank José P. Pombal Jr. (Museu Nacional, Rio de Janeiro - MNRJ), Paulo Garcia (Coleção de Anfíbios do Centro de Coleções Taxonômicas da Universidade Federal de Minas Gerais - UFMG), Miguel Rodrigues (Coleção MTR do Instituto de Biociências da Universidade de São Paulo - MTR), Marcia dos Reis Gomes (Coleção de Anfíbios do Laboratório de Anfíbios e Répteis do Departamento de Zoologia, Instituto de Biologia, Universidade Federal do Rio de Janeiro - ZUFRJ) and Taran Grant (Coleção de Anfíbios
References (91)
- et al.
Demographic processes in the montane Atlantic rainforest: molecular and cytogenetic evidence from the endemic frog Proceratophrys boiei
Mol. Phylogenet. Evol.
(2012) - et al.
Ensemble forecasting of species distributions
Trends Ecol. Evol.
(2007) - et al.
Phylogenetic relationships within anuran clade Terrarana, with emphasis on the placement of Brazilian Atlantic rainforest frogs genus Ischnocnema (Anura: Brachycephalidae)
Mol. Phylogenet. Evol.
(2012) - et al.
Species delimitation, phylogeny and evolutionary demography of co-distributed, montane frogs in the southern Brazilian Atlantic Forest
Mol. Phylogenet. Evol.
(2016) - et al.
Coalescent-based species delimitation in an integrative taxonomy
TREE
(2012) - et al.
Model-based analyses reveal insular population diversification and cryptic frog species in the Ischnocnema parva complex in the Atlantic forest of Brazil
Mol. Phylogenet. Evol.
(2017) - et al.
Support vector machines for predicting distribution of Sudden Oak Death in California
Ecol. Model.
(2005) - et al.
Maximum entropy modeling of species geographic distributions
Ecol. Model.
(2006) - et al.
The Brazilian Atlantic Forest: how much is left, and how is the remaining forest distributed? Implications for conservation
Biol. Conserv.
(2009) - et al.
Neotectonic activity in the Serra do Mar rift system (southeastern Brazil)
J. South Am. Earth Sci.
(1989)
Molecular phylogeny of Neotropical rock frogs reveals a long history of vicariant diversification in the Atlantic Forest
Mol. Phylogen. Evol.
A new statistical method for haplotype reconstruction from population data
Am. J. Hum. Genet.
Accounting for decay of linkage disequilibrium in haplotype inference and missing data imputation
Am. J. Hum. Genet.
Molecular phylogeny of Ischnocnema (Anura: Brachycephalidae) with the redefinition of its series and the description of two new species
Mol. Phylogenet. Evol.
spThin: an R package for spatial thinning of species occurrence records for use in ecological niche models
Ecography
Phylogeography: The History and Formation of Species
Genealogical perspectives on the species problem
Monog. Syst. Botan.
bModelTest: Bayesian phylogenetic site model averaging and model comparison
BMC Evol. Biol.
BEAST 2: a software platform for Bayesian evolutionary analysis
PLoS Comput. Biol.
A new species of Ischnocnema parva species series (Anura, Brachycephalidae) from northern state of Rio de Janeiro, Brazil
Herpetologica
BIOCLIM – a bioclimate analysis and prediction system
Latitude, elevational climatic zonation and speciation in New World vertebrates
Proc. R. Soc. B Biol. Sci.
Stability predicts genetic diversity in the Brazilian Atlantic Forest hotspot
Science
Prediction of phylogeographic endemism in an environmentally complex biome
Proc. R. Soc. B Biol. Sci.
DOMAIN: a flexible modelling procedure for mapping potential distributions of plants and animals
Biodivers. Conserv.
A new subspecies of frog from Itatiaya
Brazil. Am. Museum Novit.
Amphibia, Anura, Brachycephalidae, Ischnocnema holti: distribution extension
Check List
Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change
Ecography
STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method
Conserv. Genet. Resour.
Detecting the number of clusters of individuals using the software structure: a simulation study
Mol. Ecol.
From high to microendemic: molecular and acoustic analyses show that Ischnocnema guentheri (Amphibia: Brachycephalidae) is endemic to Rio de Janeiro, Brazil
Conserv. Genet.
Response to increasing Southern Hemisphere winds in CCSM4
J. Clim.
Speciation in Amazonian forest birds
Science
Alternative models of vertebrate speciation in Amazonia: an overview
Biodivers. Conserv.
Reproductive modes in frogs and their unexpected diversity in the Atlantic Forest of Brazil
Bioscience
New World direct-developing frogs (Anura: Terrarana): molecular phylogeny, classification, biogeography, and con- servation
Zootaxa
The confounding effect of population structure on Bayesian skyline plot inferences of demographic history
PLoS ONE
The genetic legacy of the Quaternary ice ages
Nature
Very high resolution interpolated climate surfaces for global land areas
Int. J. Climatol.
Cited by (8)
Geographic restriction, genetic divergence, and morphological disparity in the Brazilian Atlantic Forests: Insights from Leposoma lizards (Gymnophthalmidae, Squamata)
2021, Molecular Phylogenetics and EvolutionCitation Excerpt :Moreover, there are several divergent intraspecific lineages that are geographically circumscribed, especially within L. scincoides (Figs. S4, S5). Geographically restricted species are numerous in the AF, as demonstrated by its high levels of species (micro)endemism (e.g. Fouquet et al., 2012a; Firkowski et al., 2016; Pie et al., 2018; Thomé et al., 2020). However, small-ranged and divergent lineages are less documented in the coastal forests in Bahia and northern Espírito Santo, where most Leposoma species are distributed.
Mountain passes are higher at low latitudes for madicolous insect communities of the Neotropical region
2023, Diversity and DistributionsEffects of Amazonian flying rivers on frog biodiversity and populations in the Atlantic rainforest
2023, Conservation Biology