Demographic processes in the montane Atlantic rainforest: Molecular and cytogenetic evidence from the endemic frog Proceratophrys boiei

Mitonuclear and phenotypic discordance in an Atlantic Forest frog hybrid zone

Discordance between mitochondrial and nuclear DNA is common among animals and can be the result of a number of evolutionary processes, including incomplete lineage sorting and introgression. Particularly relevant in contact zones, mitonuclear discordance is expected because the mitochondrial genome is haploid and primarily uniparentally inherited, whereas nuclear loci are evolving at slower rates. In addition, when closely related taxa come together in hybrid zones, the distribution of diagnostic phenotypic characters and their concordance with the mitochondrial or nuclear lineages can also inform on historical and ongoing dynamics within hybrid zones. Overall, genetic and phenotypic discordances provide evidence for evolutionary divergence and processes that maintain boundaries among sister species or lineages. In this study, we characterized patterns of genetic and phenotypic variation in a contact zone between Cycloramphus dubius and Cycloramphus boraceiensis, two sister species of frogs endemic to the Atlantic Coastal Forest of Brazil. We examined genomic-scale nuclear diversification across 19 populations, encompassing the two parental forms and a contact zone between them. We compared the distribution of genomic DNA variability with that of a mitochondrial locus (16S) and two morphological traits (dorsal tubercles and body size). Our results reveal multiple divergent lineages with ongoing admixture. We detected discordance in patterns of introgression across the three data types. Cycloramphus dubius males are significantly larger than C. boraceiensis males, and we posit that competition among males in the hybrid zone, coupled with mate choice by females, may be one mechanism leading to patterns of introgression observed between the species.

Montane diversification as a mechanism of speciation in neotropical butterflies

The mountains in the Atlantic Forest domain are environments that harbor a high biodiversity, including species adapted to colder climates that were probably influenced by the climatic variations of the Pleistocene. To understand the phylogeographic pattern and assess the taxonomic boundaries between two sister montane species, a genomic study of the butterflies Actinote mantiqueira and A. alalia (Nymphalidae: Acraeini) was conducted. Analyses based on partial sequences of the mitochondrial gene COI (barcode region) failed to recover any phylogenetic or genetic structure discriminating the two species or sampling localities. However, single nucleotide polymorphisms gathered using Genotyping-by-Sequencing provided a strong isolation pattern in all analyses (genetic distance, phylogenetic hypothesis, clustering analyses, and F ST statistics) which is consistent with morphology, separating all individuals of A. alalia from all populations of A. mantiqueira. The three sampled mountain ranges where A. mantiqueira populations occur-Serra do Mar, Serra da Mantiqueira, and Poços de Caldas Plateau-were identified as three isolated clusters. Paleoclimate simulations indicate that both species' distributions changed according to climatic oscillations in the Pleistocene period, with the two species potentially occurring in areas of lower altitude during glacial periods when compared to the interglacial periods (as the present). Besides, a potential path between their distribution through the Serra do Mar Mountain range was inferred. Therefore, the Pleistocene climatic fluctuation had a significant impact on the speciation process between A. alalia and A. mantiqueira, which was brought on by isolation at different mountain summits during interglacial periods, as shown by the modeled historical distribution and the observed genetic structure.

Interspecific cytogenomic comparison reveals a potential chromosomal centromeric marker in Proceratophrys frog species

Among the repetitive elements, satellite DNA (SatDNA) emerges as extensive arrays of highly similar tandemly repeated units, spanning megabases in length. Given that the satDNA PboSat01-176, previously characterized in P. boiei, prompted our interest for having a high abundance in P. boiei and potential for centromeric satellite, here, we employed various approaches, including low coverage genome sequencing, followed by computational analysis and chromosomal localization techniques in four Proceratophrys species and, investigating the genomic presence and sharing, as well as its potential for chromosomal centromere marker in Proceratophrys frog species. Our findings demonstrate that PboSat01-176 exhibits high abundance across all four Proceratophrys species, displaying distinct characteristics that establish it as the predominant repetitive DNA element in these species. The satellite DNA is prominently clustered in the peri/centromeric region of the chromosomes, particularly in the heterochromatic regions. The widespread presence of PboSat01-176 in closely related Proceratophrys species reinforces the validity of the library hypothesis for repetitive sequences. Thus, this study highlighted the utility of the satDNA family PboSat01-176 as a reliable centromeric marker in Proceratophrys species, with potential to be applied in other species of anuran amphibians. The observed sharing and maintenance of this sequence within the genus suggest possibilities for future research, particularly through expanded sampling to elucidate parameters that underlie the library hypothesis and the evolutionary dynamics of satDNA sequences.

Analysis in Proceratophrys boiei genome illuminates the satellite DNA content in a frog from the Brazilian Atlantic forest

Satellite DNAs (satDNAs) are one of the most abundant elements in genomes. Characterized as tandemly organized sequences that can be amplified into multiple copies, mainly in heterochromatic regions. The frog P. boiei (2n = 22, ZZ♂/ZW♀) is found in the Brazilian Atlantic forest and has an atypical pattern of heterochromatin distribution when compared to other anuran amphibians, with large pericentromeric blocks on all chromosomes. In addition, females of Proceratophrys boiei have a metacentric sex chromosome W showing heterochromatin in all chromosomal extension. In this work, we performed high-throughput genomic, bioinformatic, and cytogenetic analyses to characterize the satellite DNA content (satellitome) in P. boiei, mainly due to high amount of C-positive heterochromatin and the highly heterochromatic W sex chromosome. After all the analyses, it is remarkable that the satellitome of P. boiei is composed of a high number of satDNA families (226), making P. boiei the frog species with the highest number of satellites described so far. Consistent with the observation of large centromeric C-positive heterochromatin blocks, the genome of P. boiei is enriched with high copy number of repetitive DNAs, with total satDNA abundance comprising 16.87% of the genome. We successfully mapped via Fluorescence in situ hybridization the two most abundant repeats in the genome, PboSat01-176 and PboSat02-192, highlighting the presence of certain satDNAs sequences in strategic chromosomal regions (e.g., centromere and pericentromeric region), which leads to their participation in crucial processes for genomic organization and maintenance. Our study reveals a great diversity of satellite repeats that are driving genomic organization in this frog species. The characterization and approaches regarding satDNAs in this species of frog allowed the confirmation of some insights from satellite biology and a possible relationship with the evolution of sex chromosomes, especially in anuran amphibians, including P. boiei, for which data were not available.

Assessing diversity of King Crab Lithodes spp. in the south-eastern pacific using phylogeny and molecular species delimitation methods

The purpose of this study was to test the hypothesis that the genetic diversity of commercially significant species of King Crabs (Lithodes spp.) along the south-eastern Pacific (SEP) comprises different independent evolutionary units (IEUs) with spatially isolated distribution. Nine localities from inner and open waters along the SEP Chilean coast (39°S-55°S) were sampled. We analyzed sequences from 173 individuals for the mitochondrial gene Cytochrome oxidase I (COX-I), 151 individuals for the Internal Transcribed Spacer 1 (ITS) and 135 for the structural ribosomal RNA (28S). Genetic delimitation was performed through three analytical methods: ABGD, GMYC, and its Bayesian implementation, bGMYC. Bayesian phylogenetic analyses and haplotype networks were also performed. Divergence time between clades was assessed for the COX-I marker and estimated from known evolutionary rates for this marker in other crustacean species and fossil calibration from other Anomuran species. Delimitation analyses, phylogenetic analyses, and mitochondrial haplotype networks suggested the presence of two deeply divergent mitochondrial lineages of Lithodes in the SEP, referred to as Clade1 and Clade 2. Nuclear markers showed low phylogenetic resolution and therefore were unsuitable for molecular species delimitation. Divergence time analysis of the mitochondrial lineages suggests a separation between Clades of approximately 2.3 Mya. The divergence time obtained suggested that Pliocene glaciations and deglaciations cycles could be involved in hybridization events between Lithodes IEUs at southern tip of South American coasts. The different frequencies of Lithodes haplotypes in inner and open water environments along SEP coasts could be explained by events such as the last glacial maximum or by differences in the adaptation of each clade to different environments. These findings support the necessity of evaluating the taxonomic status of Lithodes individuals found along SEP coasts under an integrative taxonomy approach or through markers with other evolution rates than those already used.

The dispersal between Amazonia and Atlantic Forest during the Early Neogene revealed by the biogeography of the treefrog tribe Sphaenorhynchini (Anura, Hylidae)

The Amazonia and the Atlantic Forest, separated by the diagonal of open formations, are two ecoregions that comprise the most diverse tropical forests in the world. The Sphaenorhynchini tribe is among the few tribes of anurans that occur in both rainforests, and their historical biogeographic have never been proposed. In this study, we infer a dated phylogeny for the species of the Sphaenorhynchini and we reconstructed the biogeographic history describing the diversification chronology, and possible patterns of dispersion and vicariance, providing information about how orogeny, forest dynamics and allopatric speciation affected their evolution in South America. We provided a dated phylogeny and biogeography study for the Sphaenorhynchini tribe using mitochondrial and nuclear genes. We analyzed 41 samples to estimate the ancestral areas using biogeographical analysis based on the estimated divergence times and the current geographical ranges of the species of Sphaenorhynchini. We recovered three characteristic clades that we recognize as groups of species (S. lacteus, S. planicola, and S. platycephalus groups), with S. carneus and G. pauloalvini being the sister taxa of all other species from the tribe. We found that the diversification of the tribe lineages coincided with the main climatic and geological factors that shaped the Neotropical landscape during the Cenozoic. The most recent common ancestor of the Sphaenorhynchini species emerged in the North of the Atlantic Forest and migrated to the Amazonia in different dispersion events that occurred during the connections between these ecoregions. This is the first large-scale study to include an almost complete calibrated phylogeny of Sphaenorhynchini, presenting important information about the evolution and diversification of the tribe. Overall, we suggest that biogeographic historical of Sphaenorhynchini have resulted from a combination of repeated range expansion and contraction cycles concurrent with climate fluctuations and dispersal events between the Atlantic Forest and Amazonia.

Polymorphism in a Neotropical toad species: ontogenetic, populational and geographic approaches to chromatic variation in Proceratophrys cristiceps (Müller, 1883) (Amphibia, Anura, Odontophrynidae)

Quantifying variability is important for understanding how evolution operates in polymorphic species such as those of the genus Proceratophrys Miranda-Ribeiro, 1920, which is widely distributed in South America. P. cristiceps distribution is limited to the Caatinga biome in Brazil. We examined its chromatic variation from a populational perspective, looking at different phenetic polymorphism levels and probable chromotypic association by applying statistical and GIS tools that could facilitate future taxonomic research regarding this and other species. We characterized P. cristiceps colour patterns and re-evaluated its geographic variation, highlighting potential consequences for the taxonomy of the genus. Our results revealed six principle chromotypes whose frequencies varied among sex and ontogenetic classes. Phenotypic expression appeared to respect defined proportions and evidenced selective value for the species. We conclude that individual variation, together with typological traditionalism may overestimate the polymorphic magnitude at the population level and cause taxonomic inflation. Our data support the usefulness of P. cristiceps as a model for microevolutionary studies.

Cytogenetic approaches provide evidence of a conserved diploid number and cytological differences between Proceratophrys species (Anura: Odontophrynidae)

Taxonomic and cytogenetic aspects of Proceratophrys have not been thoroughly clarified in the literature; thus, the objective of the present study was to provide unprecedent karyotype data regarding P. schirchi, P. laticeps and P. melanopogon. Additionally, the karyotype of P. boiei (2n = 22) and its ZZ / ZW sex chromosome system was analyzed for different populations of southeastern and southern Brazil. All Proceratophrys species have a diploid number of 2n = 22 chromosomes. In P. schirchi, a strong C-band was found in the long arm in one of the homologues of the pair 8 in female metaphasic cells, denoting that this chromosome pair could represent the heteromorphic sex chromosome in a ZZ / ZW sex system. Despite the conserved number of chromosomes, there are considerable chromosomal differences in P. melanopogon and P. boiei (southern Brazil), strongly differentiating them cytogenetically from other species of the genus. Moreover, with the confirmation of chromosomal heteromorphism related to sexual differentiation in P. boiei and the possible description of this system in P. schirchi, the Proceratophrys genus can be regarded as an interesting group for evolutionary studies and sex chromosome differentiation in anurans.

Systematics and phylogeography of the Brazilian Atlantic Forest endemic harvestmen Neosadocus Mello-Leitão, 1926 (Arachnida: Opiliones: Gonyleptidae)

Neosadocus harvestmen are endemic to the Southern Brazilian Atlantic Forest. Although they are conspicuous and display great morphological variation, their evolutionary history and the biogeographical events underlying their diversification and distribution are still unknown. This contribution about Neosadocus includes the following: a taxonomic revision; a molecular phylogenetic analysis using mitochondrial and nuclear markers; an investigation of the genetic structure and species' diversity in a phylogeographical framework. Our results show that Neosadocus is a monophyletic group and comprises four species: N. bufo, N. maximus, N. robustus and N. misandrus (which we did not find on fieldwork and only studied the female holotype). There is astonishing male polymorphism in N. robustus, mostly related to reproductive strategies. The following synonymies have resulted from this work: "Bunoweyhia" variabilis Mello-Leitão, 1935 = Neosadocus bufo (Mello-Leitão, 1926); and "Bunoweyhia" minor Mello-Leitão, 1935 = Neosadocus maximus (Giltay, 1928). Most divergences occurred during the Miocene, a geological epoch marked by intense orogenic and climatic events in the Brazilian Atlantic Forest. Intraspecific analyses indicate strong population structure, a pattern congruent with the general behavior and physiological constraints of Neotropical harvestmen.

North-south and climate-landscape-associated pattern of population structure for the Atlantic Forest White Morpho butterflies

Atlantic Forest White Morpho butterflies, currently classified as Morpho epistrophus and M. iphitus, are endemic to the Atlantic Forest, where they are widely distributed throughout heterogeneous environmental conditions. Studies with endemic butterflies allow to elucidate questions on both patterns of diversity distribution and current and past processes acting on insect groups in this biodiversity hotspot. In the present study, we characterized one mtDNA marker (COI sequences) and developed 11 polymorphic loci of microsatellite for 22 sampling locations distributed throughout the entire Atlantic Forest domain. We investigated both the taxonomic limits of taxa classified as White Morpho and the structure and distribution of the genetic diversity throughout their populations. Genetic markers and distribution data failed to identify species diversification, population structure, or isolation among subpopulations attributed to different taxa proposed for the White Morpho, suggesting that the current distinction between two species is unreasonable. The Bayesian coalescent tree based on COI sequences also failed to recover monophyletic clades for the putative species, and pointed instead to a north-south oriented pattern of genetic structure, with the northern clade coalescing later than the southern clade. Northern samples also showed more intragroup structure than southern samples based on mtDNA data. Clustering tests based on microsatellites indicated the existence of three genetic clusters, with turnover between the states of Paraná and São Paulo. The north-south pattern found for the White Morpho populations is showed for the first time to a endemic AF insect and coincides with the two different bioclimatic domains previously described for vertebrates and plants. Population structure observed for these butterflies is related to climate- and landscape-associated variables, mainly precipitation and elevation.

A review of the diagnosis and geographical distribution of the recently described flea toad Brachycephalus sulfuratus in relation to B. hermogenesi (Anura: Brachycephalidae)

Background

The flea toad Brachycephalus sulfuratus was recently described from southeastern and southern Brazil. In its description, the authors overlooked previous records of flea toads that had been identified as “Brachycephalus sp. nov.” and B. hermogenesi occurring in the same regions, which could suggest the possibility of up to three flea toads coexisting in southern Brazil. In addition, B. sulfuratus is characterized by substantial phenotypic variability, to an extent that compromises its current diagnosis with respect to its congener B. hermogenesi. Therefore, the current state-of-affairs regarding the geographical distribution of these two species and the identification of previously known populations is hitherto uncertain. Our goals are to reassess previous records of flea toads attributable to B. hermogenesi, B. sulfuratus and “Brachycephalus sp. nov.”, considering the description of B. sulfuratus, and to review the diagnosis of B. sulfuratus.

Methods

A critical analysis of the species identity of flea toad specimens attributable to B. hermogenesi, B. sulfuratus, or to a potentially undescribed species from southeastern and southern Brazil was based either on the analysis of morphology or on their advertisement calls. These analyses include our independent examinations of specimens and, when not possible, examinations of published descriptions. To allow for a consistent comparison of advertisement calls between B. hermogenesi and B. sulfuratus, we made recordings of both species, including in the type locality of the former.

Results

We found that morphological and call characters originally proposed as diagnostic for B. sulfuratus in relation to B. hermogenesi vary intraspecifically. Live individuals with ventral yellow spots correspond to B. sulfuratus; individuals without yellow spots can be either B. sulfuratus or B. hermogenesi. In preservative, they are indistinguishable. Previous records of Brachycephalus sp. nov. correspond to B. sulfuratus. We propose that the reduced number of notes per call and the presence of only isolated notes in the call of B. sulfuratus, as opposed to a high number of notes per call with isolated notes and note groups in the call of B. hermogenesi, as the only diagnostic characters between them. Regarding their distributions and based in our assessment, only B. sulfuratus occurs in southern Brazil, without any overlap with B. hermogenesi. There is a narrow gap between the distributions of these species around the southeast of the city of São Paulo. Our revision also revealed that some records previously attributed to B. hermogenesi in Rio de Janeiro and north São Paulo represent a distinct, unidentified flea toad that is not B. sulfuratus. Both species occur side by side in Corcovado, São Paulo, a locality from where five paratypes of B. hermogenesi were obtained. Biogeographic events that might have led to vicariance between B. hermogenesi and B. sulfuratus are discussed.

Quaternary climatic fluctuations influence the demographic history of two species of sky-island endemic amphibians in the Neotropics

We evaluated the role of Quaternary climatic fluctuations on the demographic history and population structure of amphibian species endemic to the 'campo rupestre' in the Neotropics, evaluating their distributional shifts, demographic changes, and lineage formation from the end of Pleistocene to present. We chose two anurans endemic to the high-elevation 'campo rupestre' in the Espinhaço Range (ER) in northeastern and southeastern Brazil (Bokermannohyla alvarengai and Bokermannohyla oxente), as models to test the role of Quaternary climatic fluctuations over their distribution range in this region. We collected tissue samples throughout their distribution range and used statistical phylogeography to examine processes of divergence and population demography. We generated spatial-temporal reconstructions using Bayesian inference in a coalescent framework in combination with hind-cast projections of species distribution models (SDMs). We also used the results and literature information to test alternative diversification scenarios via approximate Bayesian computation (ABC). Our results show that Quaternary climatic fluctuations influenced the geographic ranges of both species showing population expansion during the last glacial maximum (LGM) and range contraction during interglacial periods, as inferred from selected ABC models and from past projections of SDMs. We recovered Pleistocene diversification for both species occuring in distinctly unique periods for each taxon. An older and range-restricted lineage was recovered in a geographically isolated geological massif, deserving conservation and further taxonomic study. The diversification and distribution of these amphibian species endemic to the Neotropical 'campo rupestre' were influenced by Quaternary climatic fluctuations. The expansion of cold adapted species restricted to higher elevations during glacial periods and their concomitant retraction during interglacial periods may have been crucial for producing patterns of species richness and endemism along elevation gradients in tropical and subtropical domains. Such processes may influence the evolution of the biota distributed in heterogeneous landscapes with varied topography.

Geographic restriction, genetic divergence, and morphological disparity in the Brazilian Atlantic Forests: Insights from Leposoma lizards (Gymnophthalmidae, Squamata)

Lineage differentiation, long-term persistence, and range limitation promote high levels of phylogenetic and phylogeographic endemisms and likely underlie the abundant morphologically cryptic diversity observed in the Brazilian Atlantic Forests (AF). We explore lineage differentiation and range restriction in the AF and ask if genetic divergence and morphological disparity are correlated by integrating coalescent-based species delimitation, molecular phylogenetic, and morphological analyses in the lizard genus Leposoma. We present the first species tree for Leposoma and of their tribe, the Ecpleopodini. The analyses are based on the largest dataset ever assembled for Leposoma in terms of number of species (all represented), genetic markers (12 loci), and geographic coverage (~2,500 km). The exercise allows us to robustly delimit species within the genus and phylogeographic lineages within all species. We find support for the monophyly of the genus and for the recognition of a yet undescribed species around the Baía de Todos-os-Santos, in the state of Bahia; this form is distinct from all other congeners, both genetically and morphologically. We find that L. baturitensis, from the northeastern state of Ceará, is basal to the genus - and sister to a clade of six species restricted to the AF across the eastern coast of Brazil. Relationships within this coastal clade are ((((L. annectans, Leposoma sp.), L. scincoides), L. puk) (L. nanodactylus, L. sinepollex)). Phylogenetic and phylogeographic analyses, together with precise distribution data, allowed us to update the ranges of species and phylogeographic lineages. We reveal pervasive geographic restriction of divergent lineages in Leposoma at and below species level and discuss how forest refuges and rivers might have contributed to it. We find that morphological disparity lags behind genetic divergence in the genus because although they are correlated, the first accumulates at a much slower rate than the latter. We hope to encourage new studies in the area of AF north of the Doce river; phylogeographic sampling in that region has been much less common relative to southern sites, yet it may hold the key to several important processes defining biodiversity patterns in eastern Brazil. This appears to specially apply to processes underlying geographic restriction of morphologically cryptic, yet genetic divergent lineages, as the case of Leposoma.

Comparative phylogeographic and demographic analyses reveal a congruent pattern of sister relationships between bird populations of the northern and south-central Atlantic Forest

The Pernambuco Center of Endemism (PCE) is the northernmost strip of the Atlantic Forest (AF). Biogeographic affinities among avifaunas in the PCE, the southern-central Atlantic Forest (SCAF), and Amazonia (AM) have not been studied comprehensively, and current patterns of genetic diversity in the PCE remain unclear. The interplay between species' ecological attributes and historical processes, such as Pleistocene climate fluctuations or the appearance of rivers, may have affected population genetic structures in the PCE. Moreover, the role of past connections between the PCE and AM and the elevational distribution of species in assembling the PCE avifauna remain untested. Here, we investigated the biogeographic history of seven taxa endemic to the PCE within a comparative phylogeographic framework based on a mean of 3,618 independent single nucleotide polymorphisms (SNPs) extracted from flanking regions of ultraconserved elements (UCEs) and one mitochondrial gene. We found that PCE populations were more closely related to SCAF populations than they were to those in AM, regardless of their elevational range, with divergence times placed during the Mid-Pleistocene. These splits were consistent with a pattern of allopatric divergence with gene flow until the upper Pleistocene and no signal of rapid changes in population sizes. Our results support the existence of a Pleistocene refugium driving current genetic diversity in the PCE, thereby rejecting the role of the São Francisco River as a primary barrier for population divergence. Additionally, we found that connections with Amazonia also played a significant role in assembling the PCE avifauna through subsequent migration events.

Life history and ecology might explain incongruent population structure in two co-distributed montane bird species of the Atlantic Forest

Comparative phylogeography is a powerful approach to investigate the role of historical and environmental processes in the evolution of biodiversity within a region. In this regard, comparative studies of species with similar habitat preferences are valuable to reduce the confounding influence of habitat association when interpreting phylogeographic patterns. In the Atlantic Forest of South America, phylogeographic studies of highland and lowland species have shown distinct population structure patterns so far, suggesting that such species have responded differently to Pleistocene glacial cycles. Herein, we performed a comparative analysis using molecular data and paleodistribution models of two Montane Atlantic Forest (MAF) co-distributed passerine birds with similar habitat requirements but with distinct life-history traits and ecologies: the frugivore lek-breeding Blue Manakin (Chiroxiphia caudata) and the insectivore and socially monogamous Drab-Breasted Bamboo Tyrant (Hemitriccus diops). We aimed to shed light on the role of contrasting life histories and ecologies onto the demography and population structure of MAF species. We sampled both species throughout most of their distribution range, sequenced a mitochondrial and a nuclear molecular marker, and used standard phylogeographic methods to investigate population structure and ecological niche modeling (ENM) to infer the species' paleodistributions. Our analyses recovered a phylogeographic break in H. diops in the region of the Doce River, but no genetic structure in C. caudata. We also found higher differentiation among subpopulations within each lineage of H. diops than among subpopulations of C. caudata. We suggest that these discrepancies in population structure might be due to distinct life-history traits and their impact on gene flow and generation time. For example, while H. diops is an insectivore species, C. caudata is a frugivore and the latter ecological aspect likely selects for a higher dispersion distance. Additionally, because C. caudata is a lek-breeding species, it has a longer generation time than H. diops. These traits could hinder genetic differentiation when populations become geographically isolated. Nonetheless, both species showed some common biological features, such as signatures of synchronous population expansion and larger distribution ranges during the Last Glacial Maximum, possibly due to similar cold tolerance.

Areas of endemism of land planarians (Platyhelminthes: Tricladida) in the Southern Atlantic Forest

Areas of endemism (AoE) are the main study units in analytical biogeographic methods, and are often defined as an area with two or more endemic species living in them, presenting substantial congruence among their range limits. We explored the distribution of land planarians (Geoplanidae, Platyhelminthes) across the southern region of the Brazilian Atlantic forest (from the state of Rio de Janeiro, to the state of Rio Grande do Sul) utilizing DaSilva's et al. (2015) protocol. We used two methods, Endemicity Analysis (EA), and Geographical Interpolation of Endemism (GIE). We identified nine AoE of terrestrial flatworms in the Southern Atlantic forest. Performance of the methodologies is discussed. These AoE of land planarians can be explained through vicariance events combined with their physiological and ecological own limitations. Interestingly, these AoE are congruent with fine-scale approaches such as that with harvestmen. Most land planarians have revealed to present a very small distributional range evidencing their potential as a good model for fine-scale studies of AoE.

Vicariance, dispersal, extinction and hybridization underlie the evolutionary history of Atlantic forest fire-eye antbirds (Aves: Thamnophilidae)

In order to gain insights into the biogeographic processes underlying biotic diversification in the Atlantic Forest (AF), we used a multi-locus approach to examine the evolutionary history of the White-shouldered Fire-eye (Pyriglena leucoptera) and the Fringe-backed Fire-eye (Pyriglena atra), two parapatric sister species endemic to the AF. We sequenced one mitochondrial, three Z chromosome-linked and three anonymous markers of 556 individuals from 66 localities. We recovered four lineages throughout the AF: P. atra and three populations within P. leucoptera. All populations diverged during the late Pleistocene and presented varying levels of admixture. One Z-linked locus showed the highest level of differentiation between the two species. On the other hand, a mitochondrial haplotype was shared extensively between them. Our data supported vicariance driving speciation along with extinction and dispersal as processes underlying intraspecific diversification. Furthermore, signatures of demographic expansion in most populations and areas of genetic admixture were recovered throughout the AF, suggesting that forest fragmentation was also important in differentiation. Genetic admixture areas are located between large rivers suggesting that AF rivers may diminish gene flow. Our results indicated a complex and dynamic biogeographic history of Pyriglena in the AF, with vicariance, extinction, dispersal and secondary contact followed by introgression likely influencing the current patterns of genetic distribution.

Multilocus phylogeny of Paratelmatobiinae (Anura: Leptodactylidae) reveals strong spatial structure and previously unknown diversity in the Atlantic Forest hotspot

The Brazilian Atlantic Forest harbors high levels of anuran diversity and endemism, including several taxa restricted to small geographic ranges. Here, we provide a multilocus phylogeny for Paratelmatobiinae, a leptodactylid subfamily composed of small-ranged species distributed in the Brazilian Atlantic Forest and in the campo rupestre ecosystem. We performed Bayesian inference and maximum likelihood analyses using three mitochondrial and five nuclear markers, and a matrix comprising a broad taxonomic sampling. We then delimitated independently evolving lineages within the group. We recovered Paratelmatobiinae and each of its four genera as monophyletic and robustly supported. Five putatively new species included in our analyses were unambiguously supported in the phylogenetic trees and delimitation analyses. We also recovered other deeply divergent and geographically structured lineages within the four genera of Paratelmatobiinae. Our estimation of divergence times indicates that diversification in the subfamily began in the Eocene and continued until the Pleistocene. We discuss possible scenarios of diversification for the four genera of Paratelmatobiinae, and outline the implications of our findings for taxonomy and conservation.

Outstanding diversity and microendemism in a clade of rare Atlantic Forest montane frogs

Distributed across topographically complex landscapes that vary from lowland to high elevation, the Atlantic Forest harbors one of the richest biotas worldwide. Atlantic Forest amphibians are particularly speciose, taxonomic accounts are rising and the group is used as model for biogeographic inference. Past climate-related habitat fragmentation is often invoked to explain diversification, with montane taxa expected to become more widespread during glacial times and restrained at interglacials. In this study we investigate diversification in Ischnocnema lactea and I. holti (Anura: Brachycephalidae), two rare frog species inhabiting Atlantic Forest montane regions in Southeastern Brazil. Previous phylogenetic accounts have suggested uncertain limits between these two sister species. We assembled a multilocus DNA dataset, delimited lineages in this clade, and used ecological niche modeling to explore past and future putative ranges. Assignment analyses and traditional and coalescent phylogenetic methods confirmed the existence of a species complex of Miocene origin comprising nine lineages, most of which show very narrow ranges. Lineages were fully supported as species based in coalescent species delimitation, but the phylogenetic relationships among lineages in higher elevation were unresolved. Models of past ranges suggest extensive suitable areas at the last glacial maximum which, along with phylogenetic uncertainty, are consistent with a hypothesis that climate driven vicariance at higher elevation areas resulted in hard polytomies. Species distribution models under future climates suggest narrower ranges of the lineages relative to now, but no species are currently considered endangered. Overall, our results argue in favor for the reassessment of the taxonomic and conservation status of the I. holti - I. lactea species complex.

Climatic dynamics and topography control genetic variation in Atlantic Forest montane birds

Montane organisms responded to Quaternary climate change by tracking suitable habitat along elevational gradients. However, it is unclear whether these past climatic dynamics generated predictable patterns of genetic diversity in co-occurring montane taxa. To test if the genetic variation is associated with historical changes in the elevational distribution of montane habitats, we integrated paleoclimatic data and a model selection approach for testing the demographic history of five co-distributed bird species occurring in the southern Atlantic Forest sky islands. We found that changes in historical population sizes and current genetic diversity are attributable to habitat dynamics among time periods and the current elevational distribution of populations. Taxa with populations restricted to the more climatically dynamic southern mountain block (SMB) had, on average, a six-fold demographic expansion, whereas the populations from the northern mountain block (NMB) remained constant. In the current configuration of the southern Atlantic Forest montane habitats, populations in the SMB have more widespread elevational distributions, occur at lower elevations, and harbor higher levels of genetic diversity than NMB populations. Despite the apparent coupling of demographic and climatic oscillations, our data rejected simultaneous population structuring due to historical habitat fragmentation. Demographic modeling indicated that the species had different modes of differentiation, and varied in the timing of divergence and the degree of gene flow across mountain blocks. Our results suggest that the heterogeneous distribution of genetic variation in birds of the Atlantic Forest sky islands is associated with the interplay between topography and climate of distinct mountains, leading to predictable patterns of genetic diversity.

Cryptic diversity and ancient diversification in the northern Atlantic Forest Pristimantis (Amphibia, Anura, Craugastoridae)

The Neotropical region is known both for its megadiverse fauna and for the significant deficiency of our knowledge on species limits. The Amazon and Atlantic Forest are the two most diverse and large rainforests in South America, and they harbor many groups of sister species and sister genera. The frog genus Pristimantis is the most speciose genus of terrestrial vertebrates with 546 species, but only three of them occur in the Atlantic forest. Herein, we investigate the diversification history and phylogenetic relationship among the Atlantic Forest Pristimantis lineages in a spatial-temporal framework, using mitochondrial and nuclear genes. Our results reveal a remarkable hidden diversity, with nine highly structure lineages that may correspond to undescribed species, with many cases of sympatry among these divergent lineages. Atlantic forest Pristimantis form a monophyletic group that started to diversify over 40 million years ago. This ancient group shows diversification events that remount the early bursts of the Atlantic Forest diversification history, as well as lineage diversification likely resulting from recent Pleistocene climatic fluctuations. Future work must concentrate in comparing these lineages under an integrative framework including morphology, advertisement calls and other ecology traits to confidently delimit species of Pristimantis in the Atlantic Forest.

Great Abundance of Satellite DNA in Proceratophrys (Anura, Odontophrynidae) Revealed by Genome Sequencing

Most eukaryotic genomes contain substantial portions of repetitive DNA sequences. These are located primarily in highly compacted heterochromatin and, in many cases, are one of the most abundant components of the sex chromosomes. In this sense, the anuran Proceratophrys boiei represents an interesting model for analyses on repetitive sequences by means of cytogenetic techniques, since it has a karyotype with large blocks of heterochromatin and a ZZ/ZW sex chromosome system. The present study describes, for the first time, families of satellite DNA (satDNA) in the frog P. boiei. Its genome size was estimated at 1.6 Gb, of which 41% correspond to repetitive sequences, including satDNAs, rDNAs, transposable elements, and other elements characterized as non-repetitive. The satDNAs were mapped by FISH in the centromeric and pericentromeric regions of all chromosomes, suggesting a possible involvement of these sequences in centromere function. SatDNAs are also present in the W sex chromosome, occupying the entire heterochromatic area, indicating a probable contribution of this class of repetitive DNA to the differentiation of the sex chromosomes in this species. This study is a valuable contribution to the existing knowledge on repetitive sequences in amphibians. We show the presence of repetitive DNAs, especially satDNAs, in the genome of P. boiei that might be of relevance in genome organization and regulation, setting the stage for a deeper functional genome analysis of Proceratophrys.

Historical biogeography and multi-trait evolution in miniature toadlets of the genus Brachycephalus (Anura: Brachycephalidae)

Evolutionary changes towards a miniaturized body plan may directly affect other important phenotypic traits related to the physiology, behaviour and ecology of organisms. The frog genus Brachycephalus is an outstanding example of a radiation of miniaturized species endemic to the Brazilian Atlantic Forest. We inferred ancestral states and historical changes in body size, body colour and hyperossification to test hypotheses about diversification and selective environmental mechanisms leading to the evolution of these specialized traits. The ancestral distribution was associated with high-elevation regions in the northern Serra do Mar mountain range, and diversification in the genus was coincident with important geological and climatic events during the history of the Atlantic Forest. The dynamic historical changes provided an opportunity for multiple lowland lineages and for speciation via dispersal and vicariance in multiple invasions of the highlands. The ancestral Brachycephalus was reconstructed as miniaturized and dull coloured, without hyperossification in the skin, skull or postcranial skeleton. A parallel evolution of phenotypic traits has occurred in northern and southern Atlantic Forest lineages, beginning in the Miocene. Shifts in body size are not related to elevation range or latitude. However, we found a significant correlation between the evolution of hyperossification and aposematism with increasing body size.

Sufficient versus optimal climatic stability during the Late Quaternary: using environmental quality to guide phylogeographic inferences in a Neotropical montane system

Quaternary climatic oscillations affected species distributions worldwide, creating cycles of connectivity and isolation that impacted population demography and promoted lineage divergence. These effects have been well studied in temperate regions. Taxa inhabiting mesic montane habitats in tropical ecosystems show high levels of endemism and diversification in the distinct mountain ranges they inhabit; such a pattern has commonly been ascribed to past climatic oscillations, but few phylogeographic studies have tested this hypothesis. Here, we combine ecological niche models of species distributions with molecular data to study phylogeographic patterns in two rodents endemic to the highlands of Costa Rica and western Panama (Reithrodontomys creper and Nephelomys devius). In so doing, we apply a novel approach that incorporates a basic ecological principle: the expected positive relationship between environmental suitability and population abundance. Specifically, we use niche models to predict potential patterns of population connectivity and stability of different suitability levels during climatic extremes of the last glacial–interglacial cycle; we then test these predictions with population genetic analyses of a mitochondrial and a nuclear marker. The detailed predictions arising from the different levels of suitability were moderately to highly congruent with the molecular data depending on the species. Overall, results suggest that in these tropical montane ecosystems, cycles of population connectivity and isolation followed a pattern opposite to that typically described for temperate or lowland tropical ecosystems: namely, higher connectivity during the colder glacials, with isolation in montane refugia during the interglacials, including today. Nevertheless, the individualistic patterns for each species indicate a potentially wide gamut of phylogeographic histories reflecting particularities of their niches. Taken together, this study illustrates how phylogeographic inferences may benefit from niche model outputs that provide more detailed predictions of connectivity and finer characterizations of potential refugia through time.

Historical climate changes and hybridization shaped the evolution of Atlantic Forest spinetails (Aves: Furnariidae)

Combining phylogeographic approaches and hybrid zone inference in a single framework is a robust way to depict respectively the biogeographic history of lineages and the evolutionary processes responsible for speciation. Here, we studied the spatiotemporal patterns of diversification and characterize the hybrid zone between two Atlantic Forest spinetails (Synallaxis ruficapilla and Synallaxis cinerea) using mitochondrial DNA and nuclear (autosomal and Z-linked) genes. We consistently recovered divergence between and within the two species during the late Pliocene and Pleistocene using an isolation with migration model. Also, our results indicate distinct levels of introgression among lineages. Ecological niche models and demographic inferences, used to infer range distributions throughout the late Quaternary, were not consistent with the hypothesis of a large river as a primary barrier responsible for the divergence of the two species. Instead, a scenario of isolation and divergence followed by geographic expansion and admixture as a consequence of Quaternary climatic oscillations was supported. Paleomodels also were not consistent with the idea that the hybrid zone originated in primary differentiation and favor a secondary contact scenario. Model fitting indicated that clines of different loci spanning the hybrid zone are coincident and concordant. The narrow cline for one Z-linked locus could be indicative of some form of post-zygotic selection hindering genetic homogenization between the two species.

Phylogeography of Atlantic Forest glassfrogs (Vitreorana): when geography, climate dynamics and rivers matter

Recent advances in the field of landscape genetics provide ways to jointly analyze the role of present-day climate and landscape configuration in current biodiversity patterns. Expanding this framework into a phylogeographic study, we incorporate information on historical climatic shifts, tied to descriptions of the local topography and river configuration, to explore the processes that underlie genetic diversity patterns in the Atlantic Forest hotspot. We study two montane, stream-associated species of glassfrogs: Vitreorana eurygnatha and V. uranoscopa. By integrating species distribution modeling with geographic information systems and molecular data, we find that regional patterns of molecular diversity are jointly explained by geographic distance, historical (last 120 ky) climatic stability, and (in one species) river configuration. Mitochondrial DNA genealogies recover significant regional structure in both species, matching previous classifications of the northern and southern forests in the Atlantic Forest, and are consistent with patterns reported in other taxa. Yet, these spatial patterns of genetic diversity are only partially supported by nuclear data. Contrary to data from lowland taxa, historical climate projections suggest that these montane species were able to persist in the southern Atlantic Forest during glacial periods, particularly during the Last Glacial Maximum. These results support generally differential responses to climatic cycling by northern (lowland) and southern (montane) Atlantic Forest species, triggered by the joint impact of regional landscape configuration and climate change.

Neither barriers nor refugia explain genetic structure in a major biogeographic break: phylogeography of praying mantises in the Brazilian Atlantic Forest

The Atlantic Forest is one of the world's top biodiversity hotspots, but the diversification processes of its biota are still poorly known, with competing models attributing dominant roles to either Quaternary climatic changes or geographic barriers. Many studies identify the Doce river as a major phylogeographic break, but the reasons for this phenomenon are highly debated. Here we test the predictions of the refugial and barrier models for a common species of praying mantis, Miobantia fuscata, focusing in the areas immediately south and north of the Doce river. Our analyses show high intraspecific genetic diversity, deep coalescence times and no evidence for recent population expansion. Phylogeographic structure is inconsistent with a refugial hypothesis. Significant gene flow between northern and southern populations also conflicts with a strong role for geographic barriers. This study highlights the need for considering invertebrate taxa to infer recent landscape changes, and points towards a more complex picture of genetic diversification in the Atlantic Forest.

Integrated Analyses of Cuticular Hydrocarbons, Chromosome and mtDNA in the Neotropical Social Wasp Mischocyttarus consimilis Zikán (Hymenoptera, Vespidae)

In the present work, we explored multiple data from different biological levels such as cuticular hydrocarbons, chromosomal features, and mtDNA sequences in the Neotropical social wasp Mischocyttarus consimilis (J.F. Zikán). Particularly, we explored the genetic and chemical differentiation level within and between populations of this insect. Our dataset revealed shallow intraspecific differentiation in M. consimilis. The similarity among the analyzed samples can probably be due to the geographical proximity where the colonies were sampled, and we argue that Paraná River did not contribute effectively as a historical barrier to this wasp.

Historical relationships of areas of endemism of the Brazilian Atlantic rain forest: a cladistic biogeographic analysis of harvestman taxa (Arachnida: Opiliones)

Based on a cladistic biogeographic analysis of 6 species-level phylogenies of harvestman taxa, we searched for congruence in the historical relationships of 12 areas of endemism of the Brazilian Atlantic Rain Forest. We constructed general area cladograms using Primary Brooks Parsimony Analysis (BPA), BPA of nodes, and paralogy-free subtree analysis. These analyses resulted in 6 general area cladograms, that allow to infer a general pattern of the relationships among areas of endemism from the Brazilian Atlantic Rain Forest. Northern areas resulted related basally showing main disjunctions at the Doce River Valley and Todos os Santos Bay/São Francisco River Valley. The remaining areas of endemism were included in a southern and a southeastern block, separated by the Ribeira do Iguape Valley. Incongruence Length Differences tests showed no significant incongruence among the resulting cladograms and other matrix partitions. We concluded that tectonism and ancient marine transgressions were the probable processes responsible for the main disjunctions, whereas Neogene refugia seem to have caused the more recent disjunctions. The general pattern and redundancy in area relationships suggest a model of main reiterative barriers in diversification at multiple times for the evolution of the Atlantic Rain Forest. The renewal of cladistic biogeography and the search for common biogeographic patterns are discussed.

The roles of barriers, refugia, and chromosomal clines underlying diversification in Atlantic Forest social wasps

Phylogeographic studies have sought to explain the genetic imprints of historical climatic changes and geographic barriers within the Brazilian Atlantic Forest (AF) biota, and consequently two processes of diversification (refugia and barriers) have been proposed. Additionally, there is evidence that eustatic changes influenced the biogeographic history of the AF. Here we evaluate these contrasting diversification processes using two AF social wasp species – the mid-montane Synoeca cyanea and the lowland Synoeca aff. septentrionalis. We analyzed several sources of data including multilocus DNA sequence, climatic niche models and chromosomal features. We find support for idiosyncratic phylogeographic patterns between these wasps, involving different levels of population structure and genetic diversity, contrary suitable climatic conditions during the last glaciation, and contrasting historical movements along the AF. Our data indicate that neotectonics and refugia played distinct roles in shaping the genetic structure of these wasps. However, we argue that eustatic changes influenced the demographic expansion but not population structure in AF biota. Notably, these wasps exhibited chromosomal clines, involving chromosome number and decreasing of GC content, latitudinally oriented along the AF. Together, these results reinforce the need to consider individual organismal histories and indicate that barriers and refugia are significant factors in understanding AF evolution.

Environmental constraints on the compositional and phylogenetic beta-diversity of tropical forest snake assemblages

The ongoing biodiversity crisis increases the importance and urgency of studies addressing the role of environmental variation on the composition and evolutionary history of species assemblages, but especially the tropics and ectotherms remain understudied. In regions with rainy summers, coexistence of tropical ectothermic species may be determined by the partitioning of the climatic niche, as ectotherms can rely on water availability and thermoregulatory behaviour to buffer constraints along their climatic niche. Conversely, tropical ectotherms facing dry summers would have fewer opportunities to climatic niche partitioning and other processes rather than environmental filtering would mediate species coexistence. We used 218 snake assemblages to quantify the compositional (CBD) and phylogenetic (PBD) beta-diversity of snakes in the Atlantic Forest (AF) hotspot, South America. We identify two AF regions with distinct climatological regimes: dry summers in the northern-AF and rainy summers in the southern-AF. While accounting for the influence of multiscale spatial processes, we disentangle the relative contribution of thermal, water-related and topographic conditions in structuring the CBD and PBD of snake assemblages, and determine the extent in which snake assemblages under distinct climatological regimes are affected by environmental filtering. Thermal conditions best explain CBD and PBD of snakes for the whole AF, whereas water-related factors best explain the structure of snake assemblages within a same climatological regime. CBD and PBD patterns are similarly explained by spatial factors but snake assemblages facing dry summers are more affected by spatial processes operating at fine to intermediate spatial scale, whereas those assemblages in regions with rainy summers have a stronger signature of coarse-scale processes. As expected, environmental filtering plays a stronger role in southern-AF than northern-AF, and the synergism between thermal and water-related conditions is the key cause behind this difference. Differences in climatological regimes within the tropics affect processes mediating species coexistence. The influence of broad-scale gradients (e.g. temperature and precipitation) in structuring tropical ectothermic assemblages is greater in regions with rainy summers where climatic niche partitioning is more likely. Our findings highlight the potential stronger role of biotic interactions and neutral processes in structuring ectothermic assemblages facing changes towards warmer and dryer climates.

Genetic diversity of Morato's Digger Toad, Proceratophrys moratoi: spatial structure, gene flow, effective size and the need for differential management strategies of populations

The Morato's Digger Toad, Proceratophrys moratoi, is a critically endangered toad species with a marked population decline in southern Brazilian Cerrado. Despite this, new populations are being discovered, primarily in the northern part of the distribution range, which raises a number of questions with regard to the conservation status of the species. The present study analyzed the genetic diversity of the species based on microsatellite markers. Our findings permitted the identification of two distinct management units. We found profound genetic structuring between the southern populations, on the left margin of the Tietê River, and all other populations. A marked reduction was observed in the contemporary gene flow among the central populations that are most affected by anthropogenic impacts, such as extensive sugar cane plantations, which presumably decreases habitat connectivity. The results indicated reduced diversity in the southern populations which, combined with a smaller effective population size, may make these populations more susceptible to extinction. We recommend the reclassification of P. moratoi as vulnerable and the establishment of a special protection program for the southern populations. Our results provide important insights about the local extinction of southern populations of this toad.

Variability of NOR patterns in European water frogs of different genome composition and ploidy level

We studied water frogs from a complex composed of two species: Pelophylax lessonae (Camerano, 1882) (genome LL, 2n = 26) and P. ridibundus (Pallas, 1771) (RR, 2 = 26), and their natural hybrid P. esculentus (Fitzinger, 1843) of various ploidy and genome composition (RL, 2n = 26, and RRL or RLL, 3n = 39). Tetraploids RRLL were found (4n = 52) in juveniles. We applied cytogenetic techniques: AgNO3, chromomycin A3, PI and fluorescent in situ hybridization with a 28S rDNA probe. Results obtained by silver staining corresponded well with those stained with CMA3, PI and FISH. As a rule, NORs are situated on chromosomes 10. The number of Ag-NORs visible on metaphase plates was the same as the number of Ag-nucleoli present in interphase nuclei of the same individual. In all analyzed metaphases, NORs exhibited variations in size after AgNO3 and CMA3 stainings. Sixty-six individuals (out of 407 analyzed) were polymorphic for the localization and number of NORs. Fifty-one diploids had NORs only on one chromosome of pair 10. Three triploids (LLR and RRL) displayed two NORs, and two other triploid RRL individuals displayed one, instead of expected three NORs. In ten individuals extra NORs were detected on chromosomes other than 10 (chromosomes 2 and 9).

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0803-0) contains supplementary material, which is available to authorized users.

Geographical and altitudinal distribution of Brachycephalus (Anura: Brachycephalidae) endemic to the Brazilian Atlantic Rainforest

Mountains of the Brazilian Atlantic Forest can act as islands of cold and wet climate, leading to the isolation and speciation of species with low dispersal capacity, such as the toadlet species of the genus Brachycephalus. This genus is composed primarily by diurnal species, with miniaturized body sizes (<2.5 cm), inhabiting microhabitats in the leaf litter of montane forests. Still, little is known about the geographical distribution, altitudinal range, and ecological limits of most Brachycephalus species. In this study, we review the available data on the geographical and altitudinal distribution of Brachycephalus based on occurrence records compiled from literature and museums, both for the genus as a whole and separately for the three recently proposed groups of species (ephippium, didactylus, and pernix). The final ensemble dataset comprised 333 records, 120 localities, 28 described species, and six undescribed ones. Species were recorded in six relief units, the richest of which being the Serra do Mar, with 30 species. When the Serra do Mar is subdivided into three subunits, Northern, Central and Southern Serra do Mar, the number of species increase from north to the south, with records of six, nine, and 16 species, respectively. We were able to estimate the extent of occurrence of nearly half of the described species, and the resulting estimates indicate that many of them show remarkably small ranges, some of which less than 50 ha. Brachycephalus species are present from sea level to roughly 1,900 m a.s.l., with the highest richness being found between 751 and 1,000 m a.s.l. (21 spp.). The species with the broadest altitudinal range were B. didactylus (1,075 m) and Brachycephalus sp. 1 (1,035 m), both in the didactylus group, and B. ephippium (1,050 m), of the ephippium group. The broadest altitudinal amplitude for species of the pernix group was recorded for B. brunneus (535 m). The lowest altitudinal records for the pernix group were at 845 m a.s.l. in the state of Paraná and at 455 m a.s.l. in the state of Santa Catarina. The altitudinal occurrence in the pernixspecies group seems to decrease southward. Syntopy between species is also reviewed.