Morphological differentiation correlates with ecological but not with genetic divergence in a Gehyra gecko

Evaluating the drivers and engines of morphological diversification in the invasive gecko Hemidactylus mabouia (Moreau de Jonnès, 1818) (Squamata: Gekkonidae)

Development determines the range of possible phenotypes that can be produced and exposed to selection and has a major role in the evolutionary trajectories of species. Nevertheless, development is itself subject to evolutionary forces. Here, we describe differences at the ontogenetic and population levels in head and limb proportions of the invasive gecko Hemidactylus mabouia, to assess the developmental mechanisms and extrinsic forces associated with morphological diversification during colonization of novel habitats. We have found that allometric trajectories of most skeletal traits remain constant throughout postnatal development. Linear morphometric analysis did not find multivariate differences between ontogenetic stages or sexes. When comparing populations, our results showed that the divergence of the corresponding external measures was explained by shifts in the intercept of static allometry curves, indicating that differences arose early in development. Populations aggregated into two morphological groups that did not correspond to the groups formed on the basis of genetic structure. Using two different approaches, we found support for an adaptive hypothesis when comparing observed patterns of morphological variation with that expected under neutral evolutionary models.

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.

Conspicuously concealed: revision of the arid clade of the Gehyra variegata (Gekkonidae) group in Western Australia using an integrative molecular and morphological approach, with the description of five cryptic species

The methods used to detect and describe morphologically cryptic species have advanced in recent years, owing to the integrative nature of molecular and morphological techniques required to elucidate them. Here we integrate recent phylogenomic work that sequenced many genes but few individuals, with new data from mtDNA and morphology from hundreds of gecko specimens of the Gehyra variegata group from the Australian arid zone. To better understand morphological and geographical boundaries among cryptic forms, we generated new sequences from 656 Gehyra individuals, largely assigned to G. variegata group members over a wide area in Western Australia, with especially dense sampling in the Pilbara region, and combined them with 566 Gehyra sequences from GenBank, resulting in a dataset of 1,222 specimens. Results indicated the existence of several cryptic species, from new species with diagnostic morphological characters, to cases when there were no useful characters to discriminate among genetically distinctive species. In addition, the cryptic species often showed counter-intuitive distributions, including broad sympatry among some forms and short range endemism in other cases. Two new species were on long branches in the phylogram and restricted to the northern Pilbara region: most records of the moderately sized G. incognita sp. nov. are near the coast with isolated inland records, whereas the small-bodied saxicoline G. unguiculata sp. nov. is only known from a small area in the extreme north of the Pilbara. Three new species were on shorter branches in the phylogram and allied to G. montium. The moderately sized G. crypta sp. nov. occurs in the western and southern Pilbara and extends south through the Murchison region; this species was distinctive genetically, but with wide overlap of characters with its sister species, G. montium. Accordingly, we provide a table of diagnostic nucleotides for this species as well as for all other species treated here. Two small-bodied species occur in isolated coastal regions: G. capensis sp. nov. is restricted to the North West Cape and G. ocellata sp. nov. occurs on Barrow Island and other neighbouring islands. The latter species showed evidence of introgression with the mtDNA of G. crypta sp. nov., possibly due to recent connectivity with the mainland owing to fluctuating sea levels. However, G. ocellata sp. nov. was more closely related to G. capensis sp. nov. in the phylogenomic data and in morphology. Our study illustrates the benefits of combining phylogenomic data with extensive screens of mtDNA to identify large numbers of individuals to the correct cryptic species. This approach was able to provide sufficient samples with which to assess morphological variation. Furthermore, determination of geographic distributions of the new cryptic species should greatly assist with identification in the field, demonstrating the utility of sampling large numbers of specimens across wide areas.

Diversification across biomes in a continental lizard radiation

Ecological opportunity is a powerful driver of evolutionary diversification, and predicts rapid lineage and phenotypic diversification following colonization of competitor-free habitats. Alternatively, topographic or environmental heterogeneity could be key to generating and sustaining diversity. We explore these hypotheses in a widespread lineage of Australian lizards: the Gehyra variegata group. This clade occurs across two biomes: the Australian monsoonal tropics (AMT), where it overlaps a separate, larger bodied clade of Gehyra and is largely restricted to rocks; and in the larger Australian arid zone (AAZ) where it has no congeners and occupies trees and rocks. New phylogenomic data and coalescent analyses of AAZ taxa resolve lineages and their relationships and reveal high diversity in the western AAZ (Pilbara region). The AMT and AAZ radiations represent separate radiations with no difference in speciation rates. Most taxa occur on rocks, with small geographic ranges relative to widespread generalist taxa across the vast central AAZ. Rock-dwelling and generalist taxa differ morphologically, but only the lineage-poor central AAZ taxa have accelerated evolution. This accords with increasing evidence that lineage and morphological diversity are poorly correlated, and suggests environmental heterogeneity and refugial dynamics have been more important than ecological release in elevating lineage diversity.

There's more than one way to climb a tree: Limb length and microhabitat use in lizards with toe pads

Ecomorphology links microhabitat and morphology. By comparing ecomorphological associations across clades, we can investigate the extent to which evolution can produce similar solutions in response to similar challenges. While Anolis lizards represent a well-studied example of repeated convergent evolution, very few studies have investigated the ecomorphology of geckos. Similar to anoles, gekkonid lizards have independently evolved adhesive toe pads and many species are scansorial. We quantified gecko and anole limb length and microhabitat use, finding that geckos tend to have shorter limbs than anoles. Combining these measurements with microhabitat observations of geckos in Queensland, Australia, we observed geckos using similar microhabitats as reported for anoles, but geckos with relatively longer limbs were using narrower perches, differing from patterns observed in anoles and other lizards. We also observed arboreal geckos with relatively shorter proximal limb segments as compared to rock-dwelling and terrestrial geckos, similar to patterns observed for other lizards. We conclude that although both geckos and anoles have adhesive pads and use similar microhabitats, their locomotor systems likely complement their adhesive pads in unique ways and result in different ecomorphological patterns, reinforcing the idea that species with convergent morphologies still have idiosyncratic characteristics due to their own separate evolutionary histories.

Why do different oceanic archipelagos harbour contrasting levels of species diversity? The macaronesian endemic genus Pericallis (Asteraceae) provides insight into explaining the 'Azores diversity Enigma'

Background

Oceanic archipelagos typically harbour extensive radiations of flowering plants and a high proportion of endemics, many of which are restricted to a single island (Single Island Endemics; SIEs). The Azores represents an anomaly as overall levels of endemism are low; there are few SIEs and few documented cases of intra-archipelago radiations. The distinctiveness of the flora was first recognized by Darwin and has been referred to as the ‘Azores Diversity Enigma’ (ADE). Diversity patterns in the Macaronesian endemic genus Pericallis (Asteraceae) exemplify the ADE. In this study we used morphometric, Amplified Length Polymorphisms, and bioclimatic data for herbaceous Pericallis lineages endemic to the Azores and the Canaries, to test two key hypotheses proposed to explain the ADE: i) that it is a taxonomic artefact or Linnean shortfall, ie. the under description of taxa in the Azores or the over-splitting of taxa in the Canaries and (ii) that it reflects the greater ecological homogeneity of the Azores, which results in limited opportunity for ecological diversification compared to the Canaries.

Results

In both the Azores and the Canaries, morphological patterns were generally consistent with current taxonomic classifications. However, the AFLP data showed no genetic differentiation between the two currently recognized Azorean subspecies that are ecologically differentiated. Instead, genetic diversity in the Azores was structured geographically across the archipelago. In contrast, in the Canaries genetic differentiation was mostly consistent with morphology and current taxonomic treatments. Both Azorean and Canarian lineages exhibited ecological differentiation between currently recognized taxa.

Conclusions

Neither a Linnean shortfall nor the perceived ecological homogeneity of the Azores fully explained the ADE-like pattern observed in Pericallis. Whilst variation in genetic data and morphological data in the Canaries were largely congruent, this was not the case in the Azores, where genetic patterns reflected inter-island geographical isolation, and morphology reflected intra-island bioclimatic variation. The combined effects of differences in (i) the extent of geographical isolation, (ii) population sizes and (iii) geographical occupancy of bioclimatic niche space, coupled with the morphological plasticity of Pericallis, may all have contributed to generating the contrasting patterns observed in the archipelagos.

Electronic supplementary material

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

Hierarchical genetic structure shaped by topography in a narrow-endemic montane grasshopper

BACKGROUND

Understanding the underlying processes shaping spatial patterns of genetic structure in free-ranging organisms is a central topic in evolutionary biology. Here, we aim to disentangle the relative importance of neutral (i.e. genetic drift) and local adaptation (i.e. ecological divergence) processes in the evolution of spatial genetic structure of the Morales grasshopper (Chorthippus saulcyi moralesi), a narrow-endemic taxon restricted to the Central Pyrenees. More specifically, we analysed range-wide patterns of genetic structure and tested whether they were shaped by geography (isolation-by-distance, IBD), topographic complexity and present and past habitat suitability models (isolation-by-resistance, IBR), and environmental dissimilarity (isolation-by-environment, IBE).

RESULTS

Different clustering analyses revealed a deep genetic structure that was best explained by IBR based on topographic complexity. Our analyses did not reveal a significant role of IBE, a fact that may be due to low environmental variation among populations and/or consequence of other ecological factors not considered in this study are involved in local adaptation processes. IBR scenarios informed by current and past climate distribution models did not show either a significant impact on genetic differentiation after controlling for the effects of topographic complexity, which may indicate that they are not capturing well microhabitat structure in the present or the genetic signal left by dispersal routes defined by habitat corridors in the past.

CONCLUSIONS

Overall, these results indicate that spatial patterns of genetic variation in our study system are primarily explained by neutral divergence and migration-drift equilibrium due to limited dispersal across abrupt reliefs, whereas environmental variation or spatial heterogeneity in habitat suitability associated with the complex topography of the region had no significant effect on genetic discontinuities after controlling for geography. Our study highlights the importance of considering a comprehensive suite of potential isolating mechanisms and analytical approaches in order to get robust inferences on the processes promoting genetic divergence of natural populations.

Has contemporary climate change played a role in population declines of the lizard Ctenophorus decresii from semi-arid Australia?

Whilst contemporary climatic changes are small in magnitude compared to those predicted for the coming decades, they have already been linked to species range shifts and local extinctions. Elucidating the drivers behind species' responses to contemporary climate change will better inform management strategies for vulnerable and pest species alike. A recent proposal to explain worldwide local extinctions in lizards is that increasing maximum temperatures have constrained lizard activity time in the breeding season beyond extinction thresholds. Here we document a significant population decline and potential local extinction at the warm (northern) range margin of the tawny dragon, Ctenophorus decresii, a rock-dwelling lizard from the Flinders Ranges in semi-arid Australia. We developed and tested a biophysical model of tawny dragon thermoregulatory behaviour and drove the model with daily weather data for the period 1990-2009 across the Flinders Ranges. Our results indicate that potential annual activity time has likely increased over this period throughout the historic range, with within-season declines only in the summer months at the northern range limit. However, populations that have declined since 2000 have also likely experienced higher active body temperatures and more stringent retreat-site requirements (deeper crevices) than have regions where the species remains common, during a period of declining rainfall. Our laboratory estimates of thermal preference in this species were insensitive to altered nutritional and hydric state. Thus it is possible that recent population declines are linked to desiccation stress driven by higher body temperatures and declining rainfall. Our study illustrates that simple indices of the impact of climate warming on animals, such as activity restriction, may in fact reflect a variety of potential mechanisms whose ultimate outcome will be contingent on other factors such as water and shelter availability.

Evaluating evolutionary history in the face of high gene tree discordance in Australian Gehyra (Reptilia: Gekkonidae)

Species tree methods have provided improvements for estimating species relationships and the timing of diversification in recent radiations by allowing for gene tree discordance. Although gene tree discordance is often observed, most discordance is attributed to incomplete lineage sorting rather than other biological phenomena, and the causes of discordance are rarely investigated. We use species trees from multi-locus data to estimate the species relationships, evolutionary history and timing of diversification among Australian Gehyra-a group renowned for taxonomic uncertainty and showing a large degree of gene tree discordance. We find support for a recent Asian origin and two major clades: a tropically adapted clade and an arid adapted clade, with some exceptions, but no support for allopatric speciation driven by chromosomal rearrangement in the group. Bayesian concordance analysis revealed high gene tree discordance and comparisons of Robinson-Foulds distances showed that discordance between gene trees was significantly higher than that generated by topological uncertainty within each gene. Analysis of gene tree discordance and incomplete taxon sampling revealed that gene tree discordance was high whether terminal taxon or gene sampling was maximized, indicating discordance is due to biological processes, which may be important in contributing to gene tree discordance in many recently diversified organisms.

Speciation on the rocks: integrated systematics of the Heteronotia spelea species complex (Gekkota; Reptilia) from Western and Central Australia

The isolated uplands of the Australian arid zone are known to provide mesic refuges in an otherwise xeric landscape, and divergent lineages of largely arid zone taxa have persisted in these regions following the onset of Miocene aridification. Geckos of the genus Heteronotia are one such group, and have been the subject of many genetic studies, including H. spelea, a strongly banded form that occurs in the uplands of the Pilbara and Central Ranges regions of the Australian arid zone. Here we assess the systematics of these geckos based on detailed examination of morphological and genetic variation. The H. spelea species complex is a monophyletic lineage to the exclusion of the H. binoei and H. planiceps species complexes. Within the H. spelea complex, our previous studies based on mtDNA and nine nDNA loci found populations from the Central Ranges to be genetically divergent from Pilbara populations. Here we supplement our published molecular data with additional data gathered from central Australian samples. In the spirit of integrative species delimitation, we combine multi-locus, coalescent-based lineage delimitation with extensive morphological analyses to test species boundaries, and we describe the central populations as a new species, H. fasciolatussp. nov. In addition, within the Pilbara there is strong genetic evidence for three lineages corresponding to northeastern (type), southern, and a large-bodied melanic population isolated in the northwest. Due to its genetic distinctiveness and extreme morphological divergence from all other Heteronotia, we describe the melanic form as a new species, H. atrasp. nov. The northeastern and southern Pilbara populations are morphologically indistinguishable with the exception of a morpho-type in the southeast that has a banding pattern resembling H. planiceps from the northern monsoonal tropics. Pending more extensive analyses, we therefore treat Pilbara H. spelea as a single species with phylogenetic structure and morphological heterogeneity.

Congruent morphological and genetic differentiation as a signature of range expansion in a fragmented landscape

Phenotypic differentiation is often interpreted as a result of local adaptation of individuals to their environment. Here, we investigated the skull morphological differentiation in 11 populations of the white-footed mouse (Peromyscus leucopus). These populations were sampled in an agricultural landscape in the Montérégie region (Québec, Canada), at the northern edge of the distribution of the white-footed mouse. We found a strong pattern of phenotypic differentiation matching the genetic structure across these populations. Landscape fragmentation and the presence of geographic barriers, in particular north-south oriented rivers, contribute to this differentiation and modulate the pattern of rapid ongoing northward range expansion of the white-footed mouse in response to climate warming. We conclude that while large rivers and postglacial recolonization routes have shaped the current pattern of distribution and differentiation of white-footed mouse populations, further local differentiation is occurring, at the scale of the landscape. We posit that the northern expansion of the white-footed mouse is achieved through successive independent founder events in a fragmented landscape at the northern range edge of the species. The phenotypic differentiation we observe is thus a result of a number of mechanisms operating at different spatial and temporal scales.