Elevational speciation in action? Restricted gene flow associated with adaptive divergence across an altitudinal gradient

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Population diversification in the frog Mantidactylus bellyi on an isolated massif in northern Madagascar based on genetic, morphological, bioacoustic and ecological evidence

In the processes that give rise to new species, changes first occur at the population level. But with the continuous nature of the divergence process, change in biological properties delimiting the shift from “individuals of divergent populations” towards “individuals of distinct species”, as well as abiotic factors driving the change, remain largely ambivalent. Here we study diversification processes at the population level in a semi-aquatic frog, Mantidactylus (Brygoomantis) bellyi, across the diverse vegetation types of Montagne d’Ambre National Park (MANP), Madagascar. Genetic diversity was assessed with seven newly developed microsatellite markers as well as mitochondrial DNA sequences and concordance with patterns of ecological, morphological, and bioacoustic divergence evaluated. We found M. bellyi lacking mitochondrial differentiation within MANP, while microsatellite datasets partitioned them into three highly differentiated, geographically separated subpopulations (with indications for up to five subpopulations). The molecular grouping–primarily clustering individuals by geographic proximity–was coincident with differences in mean depth and width of waters, suggesting a possible role of fluvial characteristics in genetic exchange in this stream-breeding species. Genetic clustering not consistent with differences in call properties, except for dominant call frequencies under the two-subpopulations model. Morphological divergence was mostly consistent with the genetic clustering; subpopulations strongly differed by their snout-vent length, with individuals from high-elevation subpopulations smaller than those from populations below 1000 m above sea level. These results exemplify how mountains and environmental conditions might primarily shape genetic and morphological divergence in frog populations, without strongly affecting their calls.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

Spatiotemporal maintenance of flora in the Himalaya biodiversity hotspot: Current knowledge and future perspectives

Mountain ecosystems support a significant one-third of all terrestrial biodiversity, but our understanding of the spatiotemporal maintenance of this high biodiversity remains poor, or at best controversial. The Himalaya hosts a complex mountain ecosystem with high topographic and climatic heterogeneity and harbors one of the world's richest floras. The high species endemism, together with increasing anthropogenic threats, has qualified the Himalaya as one of the most significant global biodiversity hotspots. The topographic and climatic complexity of the Himalaya makes it an ideal natural laboratory for studying the mechanisms of floral exchange, diversification, and spatiotemporal distributions. Here, we review literature pertaining to the Himalaya in order to generate a concise synthesis of the origin, distribution, and climate change responses of the Himalayan flora. We found that the Himalaya supports a rich biodiversity and that the Hengduan Mountains supplied the majority of the Himalayan floral elements, which subsequently diversified from the late Miocene onward, to create today's relatively high endemicity in the Himalaya. Further, we uncover links between this Miocene diversification and the joint effect of geological and climatic upheavals in the Himalaya. There is marked variance regarding species dispersal, elevational gradients, and impact of climate change among plant species in the Himalaya, and our review highlights some of the general trends and recent advances on these aspects. Finally, we provide some recommendations for conservation planning and future research. Our work could be useful in guiding future research in this important ecosystem and will also provide new insights into the maintenance mechanisms underpinning other mountain systems.

Amphibian Speciation Rates Support a General Role of Mountains as Biodiversity Pumps

AbstractContinental mountain areas cover <15% of global land surface, yet these regions concentrate >80% of global terrestrial diversity. One prominent hypothesis to explain this pattern proposes that high mountain diversities could be explained by higher diversification rates in regions of high topographic complexity (HTC). While high speciation in mountains has been detected for particular clades and regions, the global extent to which lineages experience faster speciation in mountains remains unknown. Here we addressed this issue using amphibians as a model system (>7,000 species), and we found that families showing high speciation rates contain a high proportion of species distributed in mountains. Moreover, we found that lineages inhabiting areas of HTC speciate faster than lineages occupying areas that are topographically less complex. When comparing across regions, we identified the same pattern in five biogeographical realms where higher speciation rates are associated with higher levels of complex topography. Low-magnitude differences in speciation rates between some low and high complex topographies suggest that high mountain diversity is also affected by low extinction and/or high colonization rates. Nevertheless, our results bolster the importance of mountains as engines of speciation at different geographical scales and highlight their importance for the conservation of global biodiversity.

Evidence of Elevational Speciation in Kerteszia cruzii (Diptera: Culicidae) in the Ribeira Valley, São Paulo, Brazil

Kerteszia cruzii [former Anopheles (Kerteszia) cruzii] is a bromeliad mosquito widespread in the Brazilian Atlantic rainforest. In South-eastern Brazil, it plays an important role in malaria transmission because it was infected with at least four Plasmodium species. There is robust evidence that Ke. cruzii is a species complex. We used single nucleotide polymorphisms (SNPs) from a nextRAD sequence (nextera-tagmented, reductively amplified DNA) to investigate the genetic structure of Ke. cruzii in the Ribeira Valley, South-eastern Brazil. Furthermore, we verified whether the genetic structure was associated with forest cover, elevation, slope, and vegetation physiognomy. Our results showed two distinct lineages in the studied region associated with elevation and isolation by distance. The first lineage included samples from coastal localities and the second comprised specimens from inland or mountain sites. At one sampling locality (Esteiro do Morro in Cananéia municipality), both lineages are sympatric. These results are in accordance with previously published data that showed elevated stratification in Ke. cruzii. However, Fst values did not indicate the existence of cryptic or sister species in Ke. cruzii in this region, we concluded that elevational speciation probably occurs, and we hypothesized that differences in population structure found might be associated with the distribution of bromeliad species.

Phenotypic and genomic diversification with isolation by environment along elevational gradients in a neotropical treefrog

Understanding how geographic and environmental heterogeneity drive local patterns of genetic variation is a major goal of ecological genomics and a key question in evolutionary biology. The tropical Andes and inter-Andean valleys are shaped by markedly heterogeneous landscapes, where species experience strong selective processes. We examined genome-wide SNP data together with behavioural and ecological traits (mating calls and body size) known to contribute to genetic isolation in anurans in the banana tree-dwelling frog, Boana platanera, distributed across an environmental gradient in Central Colombia (northern South America). Here, we analysed the relationships between environmentally (temperature and precipitation) associated genetic and phenotypic differentiation and the potential drivers of isolation by environment along an elevation gradient. We identified candidate SNPs associated with temperature and body size, which follow a clinal pattern of genome-wide differentiation tightly coupled with phenotypic variation: as elevation increases, B. platanera exhibits larger body size and longer call duration with more pulses but lower pulse rate and frequency. Thus, the environmental landscape has rendered a scenario where isolation by environment and candidate loci show concordance with phenotypic divergence in this tropical frog along an elevation gradient in the Colombian Andes. Our study sets the basis for evaluating the role of temperature in the genetic structure and local adaptation in tropical treefrogs and its putative effect on life cycle (embryos, tadpoles, adults) along elevation gradients.

Perspectives on studying molecular adaptations of amphibians in the genomic era

Understanding the genetic mechanisms underlying particular adaptations/phenotypes of organisms is one of the core issues of evolutionary biology. The use of genomic data has greatly advanced our understandings on this issue, as well as other aspects of evolutionary biology, including molecular adaptation, speciation, and even conservation of endangered species. Despite the well-recognized advantages, usages of genomic data are still limited to non-mammal vertebrate groups, partly due to the difficulties in assembling large or highly heterozygous genomes. Although this is particularly the case for amphibians, nonetheless, several comparative and population genomic analyses have shed lights into the speciation and adaptation processes of amphibians in a complex landscape, giving a promising hope for a wider application of genomics in the previously believed challenging groups of organisms. At the same time, these pioneer studies also allow us to realize numerous challenges in studying the molecular adaptations and/or phenotypic evolutionary mechanisms of amphibians. In this review, we first summarize the recent progresses in the study of adaptive evolution of amphibians based on genomic data, and then we give perspectives regarding how to effectively identify key pathways underlying the evolution of complex traits in the genomic era, as well as directions for future research.

Density Dependence and Adult Survival Drive Dynamics in Two High Elevation Amphibian Populations

Amphibian conservation has progressed from the identification of declines to mitigation, but efforts are hampered by the lack of nuanced information about the effects of environmental characteristics and stressors on mechanistic processes of population regulation. Challenges include a paucity of long-term data and scant information about the relative roles of extrinsic (e.g., weather) and intrinsic (e.g., density dependence) factors. We used a Bayesian formulation of an open population capture-recapture model and >30 years of data to examine intrinsic and extrinsic factors regulating two adult boreal chorus frogs (Pseudacris maculata) populations. We modelled population growth rate and apparent survival directly, assessed their temporal variability, and derived estimates of recruitment. Populations were relatively stable (geometric mean population growth rate >1) and regulated by negative density dependence (i.e., higher population sizes reduced population growth rate). In the smaller population, density dependence also acted on adult survival. In the larger population, higher population growth was associated with warmer autumns. Survival estimates ranged from 0.30–0.87, per-capita recruitment was <1 in most years, and mean seniority probability was >0.50, suggesting adult survival is more important to population growth than recruitment. Our analysis indicates density dependence is a primary driver of population dynamics for P. maculata adults.

New insights into the evolutionary history of Megacodon: Evidence from a newly discovered species

Megacodon is an ideal genus to study speciation and ecological adaptation in the Sino-Himalayan region. The genus contains two species distributed at different elevations and in two separate areas. However, studies of this genus have long been impeded by a lack of fieldwork on one of its species, Megacodon venosus. In this study, we collected specimens of two Megacodon species and found an extraordinary new species of Megacodon in Lushui county of north-west Yunnan province, which we have since named Megacodon lushuiensis. We propose new species based on both morphological and molecular evidence. The finding of this new species emphasized the importance of ecological divergence in the divergence of Megacodon stylophorus and its parapatric low-elevation Megacodon species. To identify genetic determinants that underlie adaptations to different elevations, we characterized transcriptomes of the new species M. lushuiensis, which is distributed at low elevations, and M. stylophorus, which is distributed at high elevations. Comparative transcriptome analysis identified 8926 orthogroups containing single-copy genes, and 370 orthogroups containing significantly positively selected genes. The set of positively selected genes was enriched into 25 Gene Ontology terms, including "response to water deprivation", "response to osmotic stress", and "cellular response to external stimulus". Our results provide new insights into how ecological adaptation and speciation occurred in Megacodon and highlight the role of heterogeneous habitats in the speciation of plants in the Sino-Himalayan region.

Testing theoretical metapopulation conditions with genotypic data from Boreal Chorus Frogs (Pseudacris maculata)

The metapopulation concept has far-reaching implications in ecology and conservation biology. Hanski’s criteria operationally define metapopulations, yet testing them is hindered by logistical and financial constraints inherent to the collection of long-term demographic data. Hence, ecologists and conservationists often assume metapopulation existence for dispersal-limited species that occupy patchy habitats. To advance understanding of metapopulation theory and improve conservation of metapopulations, we used population and landscape genetic tools to develop a methodological framework for evaluating Hanski’s criteria. We used genotypic data (11 microsatellite loci) from a purported metapopulation of Boreal Chorus Frogs (Pseudacris maculata (Agassiz, 1850)) in Colorado, U.S.A., to test Hanski’s four criteria. We found support for each criterion: (1) significant genetic differentiation between wetlands, suggesting distinct breeding populations; (2) wetlands had small effective population sizes and recent bottlenecks, suggesting populations do not experience long-term persistence; (3) population graphs provided evidence of gene flow between patches, indicating potential for recolonization; and (4) multiscale bottleneck analyses suggest asynchrony, indicating that simultaneous extinction of all populations was unlikely. Our methodological framework provides a logistically and financially feasible alternative to long-term demographic data for identifying amphibian metapopulations.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Comparative landscape genetics of pond-breeding amphibians in Mediterranean temporal wetlands: The positive role of structural heterogeneity in promoting gene flow

Comparative landscape genetics studies can provide key information to implement cost-effective conservation measures favouring a broad set of taxa. These studies are scarce, particularly in Mediterranean areas, which include diverse but threatened biological communities. Here, we focus on Mediterranean wetlands in central Iberia and perform a multi-level, comparative study of two endemic pond-breeding amphibians, a salamander (Pleurodeles waltl) and a toad (Pelobates cultripes). We genotyped 411 salamanders from 20 populations and 306 toads from 16 populations at 18 and 16 microsatellite loci, respectively, and identified major factors associated with population connectivity through the analysis of three sets of variables potentially affecting gene flow at increasingly finer levels of spatial resolution. Topographic, land use/cover, and remotely sensed vegetation/moisture indices were used to derive optimized resistance surfaces for the two species. We found contrasting patterns of genetic structure, with stronger, finer scale genetic differentiation in Pleurodeles waltl, and notable differences in the role of fine-scale patterns of heterogeneity in vegetation cover and water content in shaping patterns of regional genetic structure in the two species. Overall, our results suggest a positive role of structural heterogeneity in population connectivity in pond-breeding amphibians, with habitat patches of Mediterranean scrubland and open oak woodlands ("dehesas") facilitating gene flow. Our study highlights the usefulness of remotely sensed continuous variables of land cover, vegetation and water content (e.g., NDVI, NDMI) in conservation-oriented studies aimed at identifying major drivers of population connectivity.

First estimates of the probability of survival in a small-bodied, high-elevation frog (Boreal Chorus Frog, Pseudacris maculata), or how historical data can be useful

In an era of shrinking budgets yet increasing demands for conservation, the value of existing (i.e., historical) data are elevated. Lengthy time series on common, or previously common, species are particularly valuable and may be available only through the use of historical information. We provide first estimates of the probability of survival and longevity (0.67–0.79 and 5–7 years, respectively) for a subalpine population of a small-bodied, ostensibly common amphibian, the Boreal Chorus Frog (Pseudacris maculata (Agassiz, 1850)), using historical data and contemporary, hypothesis-driven information–theoretic analyses. We also test a priori hypotheses about the effects of color morph (as suggested by early reports) and of drought (as suggested by recent climate predictions) on survival. Using robust mark–recapture models, we find some support for early hypotheses regarding the effect of color on survival, but we find no effect of drought. The congruence between early findings and our analyses highlights the usefulness of historical information in providing raw data for contemporary analyses and context for conservation and management decisions.

Unravelling the role of host plant expansion in the diversification of a Neotropical butterfly genus

BACKGROUND

Understanding the processes underlying diversification is a central question in evolutionary biology. For butterflies, access to new host plants provides opportunities for adaptive speciation. On the one hand, locally abundant host species can generate ecologically significant selection pressure. But a diversity of host plant species within the geographic range of each population and/or species might also eliminate any advantage conferred by specialization. This paper focuses on four Melinaea species, which are oligophagous on the family Solanaceae: M. menophilus, M. satevis, M. marsaeus, and finally, M. mothone. We examined both female preference and larval performance on two host plant species that commonly occur in this butterfly's native range, Juanulloa parasitica and Trianaea speciosa, to determine whether the different Melinaea species show evidence of local adaptation.

RESULT

In choice experiments, M. mothone females used both host plants for oviposition, whereas all other species used J. parasitica almost exclusively. In no choice experiment, M. mothone was the only species that readily accepted T. speciosa as a larval host plant. Larval survival was highest on J. parasitica (82.0 % vs. 60.9 %) and development took longer on T. speciosa (14.12 days vs. 13.35 days), except for M. mothone, which did equally well on both host plants. For all species, average pupal weight was highest on J. parasitica (450.66 mg vs. 420.01 mg), although this difference was least apparent in M. mothone.

CONCLUSION

We did not find that coexisting species of Melinaea partition host plant resources as expected if speciation is primarily driven by host plant divergence. Although M. mothone shows evidence of local adaptation to a novel host plant, T. speciosa, which co-occurs, it does not preferentially lay more eggs on or perform better on this host plant than on host plants used by other Melinaea species and not present in its distributional range. It is likely that diversification in this genus is driven by co-occurring Müllerian mimics and the resulting predation pressure, although this is also likely made possible by greater niche diversity as a consequence of plasticity for potential hosts.

Evolution of a Planktonic Foraminifer during Environmental Changes in the Tropical Oceans

Ecological adaptation to environmental changes is a strong driver of evolution, enabling speciation of pelagic plankton in the open ocean without the presence of effective physical barriers to gene flow. The tropical ocean environment, which plays an important role in shaping marine biodiversity, has drastically and frequently changed since the Pliocene. Nevertheless, the evolutionary history of tropical pelagic plankton has been poorly understood, as phylogeographic investigations are still in the developing state and paleontological approaches are insufficient to obtain a sequential record from the deep-sea sediments. The planktonic foraminifer Pulleniatina obliquiloculata is widely distributed in the tropical area throughout the world’s oceans, and its phylogeography is well established. It is thus one of the best candidates to examine how past environmental changes may have shifted the spatial distribution and affected the diversification of tropical pelagic plankton. Such an examination requires the divergence history of the planktonic foraminifers, yet the gene marker (partial small subunit (SSU) rDNA) previously used for phylogeographic studies was not powerful enough to achieve a high accuracy in estimating the divergence times. The present study focuses on improving the precision of divergence time estimates for the splits between sibling species (genetic types) of planktonic foraminifers by increasing the number of genes as well as the number of nucleotide bases used for molecular clock estimates. We have amplified the entire coding regions of two ribosomal RNA genes (SSU rDNA and large subunit (LSU) rDNA) of three genetic types of P. obliquiloculata and two closely related species for the first time and applied them to the Bayesian relaxed clock method. The comparison of the credible intervals of the four datasets consisting either of sequences of the partial SSU rDNA, the complete SSU rDNA, LSU rDNA, or a combination of both genes (SSU+LSU) clearly demonstrated that the two-gene dataset improved the accuracy of divergence time estimates. The P. obliquiloculata lineage diverged twice, first at the end of the Pliocene (3.1 Ma) and again in the middle Pleistocene (1.4 Ma). Both timings coincided with the environmental changes, which indirectly involved geographic separation of populations. The habitat of P. obliquiloculata was expanded toward the higher latitudinal zones during the stable warm periods and subsequently placed on the steep environmental gradients following the global cooling. Different environmental conditions in the stable warm tropics and unstable higher latitudes may have triggered ecological divergence among the populations, leading to adaptive differentiation and eventually speciation. A comprehensive analysis of divergence time estimates combined with phylogeography enabled us to reveal the evolutionary history of the pelagic plankton and to find the potential paleoenvironmental events, which could have changed their biogeography and ecology.