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.

Convergent and environmentally associated chromatic polymorphism in Bryconops Kner, 1858 (Ostariophysi: Characiformes: Iguanodectidae)

Bryconops Kner, 1858, includes two well defined subgenera based on morphological evidence, with each containing at least one species (B. (Bryconops) caudomaculatus and B. (Creatochanes) melanurus) with a very wide distribution, within which regional populations present color variations. To test if phenotypic variation is related to cladogenetic events, we performed tests for phylogenetic independence and determined the strength of convergence for color characters in relation to water type, as the variation between clear, black and white waters is considered to be one of the major driving forces in the evolution of Amazonian fishes. Color characters for fins above the median line of the body were generally found to be independent from phylogeny and the Wheatsheaf test strongly supports convergence of the dorsal fin color between populations of species in the same type of water, with a similar trend suggested for the color of the dorsal lobe of the caudal fin. This means that simple color characters cannot necessarily be relied upon for taxonomic revisions of the genus as local phenotypic variants may represent environmentally determined plasticity or convergent evolution. Further studies are required to determine the validity of these characters.

Phylogenomic and morphological data reveal hidden patterns of diversity in the national tree of Brazil, Paubrasilia echinata

PREMISE

Paubrasilia echinata (common names, pau brasil, brazilwood) is the national tree of Brazil and an endangered species endemic to the Brazilian Atlantic Forest. Over its wide distribution of 2000 km, its leaflets morphology exhibits extensive plasticity. Three morphotypes are commonly identified based on leaf size, but it is unclear if they represent distinct taxa or a single polymorphic species. This study aims to clarify the taxonomic position of the three morphotypes to inform conservation decisions.

METHODS

A morphometric study of leaf characters of herbarium specimens was coupled with genetic analyses using genotype-by-sequencing data. We used maximum-likelihood and coalescent methods to evaluate the phylogenetic and population structure of the species. We compared these with a morphological dendrogram built from hierarchical clustering.

RESULTS

Two of the three morphotypes formed separately evolving lineages, the third morphotype formed two geographically separate lineages, and northern trees with intermediate leaf morphology formed a separate fifth lineage. Leaflet size varied by over 35-fold, and although morphological clustering generally matched the genetic patterns, there were some overlaps, highlighting the cryptic diversity within this group.

CONCLUSIONS

Our genetic and morphological results provide some evidence that cultivated trees from different states in Brazil seem to have a limited genetic origin and do not reflect the broader genetic and geographical diversity of the species. As a result, more care is likely needed to preserve the overall genomic diversity of this endangered and iconic species.

Assessing model adequacy for Bayesian Skyline plots using posterior predictive simulation

Bayesian skyline plots (BSPs) are a useful tool for making inferences about demographic history. For example, researchers typically apply BSPs to test hypotheses regarding how climate changes have influenced intraspecific genetic diversity over time. Like any method, BSP has assumptions that may be violated in some empirical systems (e.g., the absence of population genetic structure), and the naïve analysis of data collected from these systems may lead to spurious results. To address these issues, we introduce P2C2M.Skyline, an R package designed to assess model adequacy for BSPs using posterior predictive simulation. P2C2M.Skyline uses a phylogenetic tree and the log file output from Bayesian Skyline analyses to simulate posterior predictive datasets and then compares this null distribution to statistics calculated from the empirical data to check for model violations. P2C2M.Skyline was able to correctly identify model violations when simulated datasets were generated assuming genetic structure, which is a clear violation of BSP model assumptions. Conversely, P2C2M.Skyline showed low rates of false positives when models were simulated under the BSP model. We also evaluate the P2C2M.Skyline performance in empirical systems, where we detected model violations when DNA sequences from multiple populations were lumped together. P2C2M.Skyline represents a user-friendly and computationally efficient resource for researchers aiming to make inferences from BSP.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

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.

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.

Strict allopatric speciation of sky island Pyrrhula erythaca species complex

Increasing evidence of post-divergence gene flow between taxa is shifting our understanding on the mode of speciation. A fundamental question arises concerning the circumstances under which strict allopatric speciation still holds true. Sky island populations might undergo reduced gene flow by niche conservatism to highland habitats and follow divergence in an allopatric manner. In this study, we tested this hypothesis in the sky island Grey-headed Bullfinch (Pyrrhula erythaca) species complex via statistical analyses of both genetic and ecological data. Results of coalescent-based analysis of multiple nuclear loci suggested that P. e. owstoni likely colonized Taiwan island during the severe mid-Pleistocene glacial climate followed by strictly allopatric divergence from P. e. erythaca distributed in Himalayas-Hengduan mountains and central North China. Results of ecological niche modeling suggested that their speciation may be attributed to the niche conservatism of these birds and the lack of a suitable ecological corridor during subsequent milder glacial episodes. In addition, we delimited the traditionally defined P. erythaca into two full species, P. erythaca in the Asian mainland and P. owstoni on the island of Taiwan, based on both genetic and behavioural evidences. These results suggest that ecology can have a dynamic role in allowing highland populations to expand their ranges and isolated by habitat barriers to diversify in a strictly allopatric manner.

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.

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.

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.

Signature of the Paleo-Course Changes in the São Francisco River as Source of Genetic Structure in Neotropical Pithecopus nordestinus (Phyllomedusinae, Anura) Treefrog

Historical processes that have interrupted gene flow between distinct evolutionary lineages have played a fundamental role in the evolution of the enormous diversity of species found in the Neotropical region. Numerous studies have discussed the role of geographic barriers and Pleistocene forest refugia in the diversification of the region's biodiversity. In the present study, we investigated the relative contribution of these different factors to the evolutionary history of Pithecopus nordestinus, a Neotropical tree frog, which is amply distributed in the Brazilian Atlantic Forest and adjacent areas of the Caatinga biome. We used an extensive sample and multilocus DNA sequences to provide an overview of the intraspecific genetic diversity of P. nordestinus, characterize historical diversification patterns, and identify possible phylogenetic splits. We tested different scenarios of diversification based on Pleistocene Refugia and river barrier models using approximate Bayesian computation (ABC) and ecological niche modeling (ENM). The phylogenetic approach indicate the occurrence of processes of phylogeographic divergence in both time and space, related to historical shifts in the course of the São Francisco River during Plio-Pleistocene period, resulting in two principal, highly divergent clades. The ABC model provided strong statistical support for this scenario, confirming the hypothesis that the São Francisco River acted as an effective geographical barrier during vicariant events in the evolutionary history of P. nordestinus. We believe that the climatic changes that occurred during the Pleistocene also played a secondary role in the genetic signatures identified, reinforcing the divergence of populations isolated by physical barriers. These findings reinforce the conclusion that the two models of diversification (geographic barriers and refugia) are not mutually exclusive in the Neotropical domain but may interact extensively during the diversification of species on a regional scale.

Symmetrical discrimination despite weak song differentiation in 2 suboscine bird sister species

Song mediates territorial competition and mate choice in birds and population divergence in this signal can have important evolutionary consequences. For example, divergent songs can act in specific recognition and limit gene flow and, hence, have a fundamental role on the origin and/or integrity of evolutionary lineages. Especially interesting systems to test the role of song in specific recognition are species pairs that present small structural differences in this signal. Here, we perform song play-back experiments on males of a long-diverged sister pair of Neotropical Suboscine species, the squamate antbird (Myrmoderus squamosus) and the white-bibbed antbird (Myrmoderus loricatus), which occur in parapatry in the Atlantic Forest and that overlap extensively in song variation. Previous evidence indicates that genetic introgression between these species is either absent or negligible, suggesting that vocal discrimination or other mechanisms function as effective barriers to gene flow. Our results show that responses to heterospecific songs were symmetrical and intermediary compared with responses to conspecific songs in both species. A stronger response to conspecific territorial songs suggests that conspecific individuals pose greater competitive threat than heterospecifics. An important implication of our study is that even small song differences can play an important role in specific recognition.

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.

What happened in the South American Gran Chaco? Diversification of the endemic frog genus Lepidobatrachus Budgett, 1899 (Anura: Ceratophryidae)

The Chaco is one the most neglected and least studied regions of the world. This highly-seasonal semiarid biome is an extensive continuous plain without any geographic barrier, and in spite of its high species diversity, the events and processes responsible have never been assessed. Miocene marine introgressions and Pleistocene glaciations have been mentioned as putative drivers of diversification for some groups of vertebrates in adjacent biomes of southern South America. Here we used multilocus data (one mitochondrial and six nuclear loci) from the three species of the endemic frog genus Lepidobatrachus (Lepidobatrachus asper, Lepidobatrachus laevis, and Lepidobatrachus llanensis) to determine if any of the historical events suggested as drivers of vertebrate diversification in southern South America are related to the diversification of the genus and if the Chaco is indeed a biome without barriers. Using fossil calibration in a coalescent framework we estimated that the genus diversified in the second half of the Miocene, coinciding with marine introgressions. Genetic patterns and historical demography suggest an important role of old archs and cratons as refuges during floods. In one species of the genus, L. llanensis, genetic structure reveals some breaks along the landscape, the main one of which corresponds to an area of the central Chaco that may act as a climatic barrier. Additionally, we found differential effects of the main Chacoan rivers on species of Lepidobatrachus that could be related to the time of persistence of populations in the areas influenced by these rivers.

Colonization and diversification of the white-browed shortwing (Aves: Muscicapidae: Brachypteryx montana) in the Philippines

Molecular phylogenetic approaches have greatly improved our knowledge of the pattern and process of biological diversification across the globe; however, many regions remain poorly documented, even for well-studied vertebrate taxa. The Philippine archipelago, one of the least-studied 'biodiversity hotspots', is an ideal natural laboratory for investigating the factors driving diversification in an insular and geologically dynamic setting. We investigated the history and geography of diversification of the Philippine populations of a widespread montane bird, the White-browed Shortwing (Brachypteryx montana). Leveraging dense archipelago-wide sampling, we generated a multi-locus genetic dataset (one nuclear and two mtDNA markers), which we analyzed using phylogenetic, population genetic, and coalescent-based methods. Our results demonstrate that Philippine shortwings (1) likely colonized the Philippines from the Sunda Shelf to Mindanao in the late Miocene or Pliocene, (2) diversified across inter-island barriers into three divergent lineages during the Pliocene and early Pleistocene, (3) have not diversified within the largest island, Luzon, contrary to patterns observed in other montane taxa, and (4) colonized Palawan from the oceanic Philippines rather than from Borneo, challenging the assumption of Palawan functioning exclusively as a biogeographic extension of the Sunda Shelf. Additionally, our finding that divergent (c. 4.0 mya) lineages are coexisting in secondary sympatry on Mindanao without apparent gene flow suggests that the speciation process is likely complete for these shortwing lineages. Overall, these investigations provide insight into how topography and island boundaries influence diversification within remote oceanic archipelagos and echo the results of many other studies in demonstrating that taxonomic diversity continues to be underestimated in the Philippines.

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.

Gene flow during glacial habitat shifts facilitates character displacement in a Neotropical flycatcher radiation

BACKGROUND

Pleistocene climatic fluctuations are known to be an engine of biotic diversification at higher latitudes, but their impact on highly diverse tropical areas such as the Andes remains less well-documented. Specifically, while periods of global cooling may have led to fragmentation and differentiation at colder latitudes, they may - at the same time - have led to connectivity among insular patches of montane tropical habitat with unknown consequences on diversification. In the present study we utilized ~5.5 kb of DNA sequence data from eight nuclear loci and one mitochondrial gene alongside diagnostic morphological and bioacoustic markers to test the effects of Pleistocene climatic fluctuations on diversification in a complex of Andean tyrant-flycatchers of the genus Elaenia.

RESULTS

Population genetic and phylogenetic approaches coupled with coalescent simulations demonstrated disparate levels of gene flow between the taxon chilensis and two parapatric Elaenia taxa predominantly during the last glacial period but not thereafter, possibly on account of downward shifts of montane forest habitat linking the populations of adjacent ridges. Additionally, morphological and bioacoustic analyses revealed a distinct pattern of character displacement in coloration and vocal traits between the two sympatric taxa albiceps and pallatangae, which were characterized by a lack of gene flow.

CONCLUSION

Our study demonstrates that global periods of cooling are likely to have facilitated gene flow among Andean montane Elaenia flycatchers that are more isolated from one another during warm interglacial periods such as the present era. We also identify a hitherto overlooked case of plumage and vocal character displacement, underpinning the complexities of gene flow patterns caused by Pleistocene climate change across the Andes.

A history of hybrids? Genomic patterns of introgression in the True Geese

Background

The impacts of hybridization on the process of speciation are manifold, leading to distinct patterns across the genome. Genetic differentiation accumulates in certain genomic regions, while divergence is hampered in other regions by homogenizing gene flow, resulting in a heterogeneous genomic landscape. A consequence of this heterogeneity is that genomes are mosaics of different gene histories that can be compared to unravel complex speciation and hybridization events. However, incomplete lineage sorting (often the outcome of rapid speciation) can result in similar patterns. New statistical techniques, such as the D-statistic and hybridization networks, can be applied to disentangle the contributions of hybridization and incomplete lineage sorting. We unravel patterns of hybridization and incomplete lineage sorting during and after the diversification of the True Geese (family Anatidae, tribe Anserini, genera Anser and Branta) using an exon-based hybridization network approach and taking advantage of discordant gene tree histories by re-sequencing all taxa of this clade. In addition, we determine the timing of introgression and reconstruct historical effective population sizes for all goose species to infer which demographic or biogeographic factors might explain the observed patterns of introgression.

Results

We find indications for ancient interspecific gene flow during the diversification of the True Geese and were able to pinpoint several putative hybridization events. Specifically, in the genus Branta, both the ancestor of the White-cheeked Geese (Hawaiian Goose, Canada Goose, Cackling Goose and Barnacle Goose) and the ancestor of the Brent Goose hybridized with Red-breasted Goose. One hybridization network suggests a hybrid origin for the Red-breasted Goose, but this scenario seems unlikely and it not supported by the D-statistic analysis. The complex, highly reticulated evolutionary history of the genus Anser hampered the estimation of ancient hybridization events by means of hybridization networks. The reconstruction of historical effective population sizes shows that most species showed a steady increase during the Pliocene and Pleistocene. These large effective population sizes might have facilitated contact between diverging goose species, resulting in the establishment of hybrid zones and consequent gene flow.

Conclusions

Our analyses suggest that the evolutionary history of the True Geese is influenced by introgressive hybridization. The approach that we have used, based on genome-wide phylogenetic incongruence and network analyses, will be a useful procedure to reconstruct the complex evolutionary histories of many naturally hybridizing species groups.

Electronic supplementary material

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

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.

Ultraconserved Elements Sequencing as a Low-Cost Source of Complete Mitochondrial Genomes and Microsatellite Markers in Non-Model Amniotes

Sequence capture of ultraconserved elements (UCEs) associated with massively parallel sequencing has become a common source of nuclear data for studies of animal systematics and phylogeography. However, mitochondrial and microsatellite variation are still commonly used in various kinds of molecular studies, and probably will complement genomic data in years to come. Here we show that besides providing abundant genomic data, UCE sequencing is an excellent source of both sequences for microsatellite loci design and complete mitochondrial genomes with high sequencing depth. Identification of dozens of microsatellite loci and assembly of complete mitogenomes is exemplified here using three species of Poospiza warbling finches from southern and southeastern Brazil. This strategy opens exciting opportunities to simultaneously analyze genome-wide nuclear datasets and traditionally used mtDNA and microsatellite markers in non-model amniotes at no additional cost.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Phenotypic evolution of an Atlantic Forest passerine (Xiphorhynchus fuscus): biogeographic and systematic implications

We studied the phenotypic variation of the Atlantic Forest passerine Xiphorhynchus fuscus (Aves: Dendrocolaptidae) with the broad aim of addressing whether the history and type of forest affected the evolution of endemic taxa. We also tested whether the different subspecies and genetic lineages of X. fuscus could be considered full species. We collected plumage and body size measurements and, in combination with genetic data, used multivariate tests to evaluate the working hypotheses. Our results, combined with previous biogeographic analyses, indicate that vicariant events have been important determinants in the evolution of phenotypic characters of X. fuscus, once genetic isolation was complete. Our analysis also suggests that forest heterogeneity and ecotones are important factors in the early evolution of Atlantic Forest taxa, perhaps via divergent selection. Forest instability during the Pleistocene was critical in the evolution of phenotypic traits. We confirm that the subspecies atlanticus should be considered a full species. Other lineages or populations are also phenotypically differentiated but we do not suggest considering them as full species. They share high levels of gene flow and are part of a continuous latitudinal cline of phenotypic variation. Our study suggests that not all the historic events in the Atlantic Forest that affected the evolution of genetic lineages also influenced the evolution of phenotypic characters in the same direction and intensity. Undoubtedly, natural selection played a major role in the evolution of Atlantic Forest organisms.