Tropical specialist vs. climate generalist: Diversification and demographic history of sister species of Carlia skinks from northwestern Australia

Dragons in the tropics - Phylogeography and speciation in Diporiphora lizards and common geographic breaks in co-distributed taxa

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species' ecology, and landscape features interact to shape regional diversity and endemism.

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.

Predictors of phylogeographic structure among codistributed taxa across the complex Australian monsoonal tropics

Differences in the geographic scale and depth of phylogeographic structure across codistributed taxa can reveal how microevolutionary processes such as population isolation and persistence drive diversification. In turn, environmental heterogeneity, species' traits, and historical biogeographic barriers may influence the potential for isolation and persistence. Using extensive SNP data and a combination of population genetic summary statistics and landscape genomic analyses, we explored predictors of the scale and depth of phylogeographic structure in codistributed lizard taxa from the topographically and climatically complex monsoonal tropics (AMT) of Australia. We first resolved intraspecific lineages and then tested whether genetic divergence across space within lineages is related to isolation by distance, resistance and/or environment and whether these factors differ across genera or between rock-related versus habitat generalist taxa. We then tested whether microevolutionary processes within lineages explain differences in the geographic scale and depth of intraspecific phylogeographic lineages. The results indicated that landscape predictors of phylogeographic structure differ between taxa. Within lineages, there was prevalent isolation by distance, but the strength of isolation by distance is independent of the taxonomic family, habitat specialization, and climate. Isolation by environment is the strongest predictor of landscape-scale genetic divergence for all taxa, with both temperature and precipitation acting as limiting factors. The strength of isolation by distance does not predict the geographic scale of the phylogeographic structure. However, more localized lineages had higher mean individual heterozygosity and less negative Tajima's D. This result implies that finer-scale phylogeographic structuring within species is associated with larger and more stable populations and, hence, persistence.

Current and Forthcoming Approaches for Benchmarking Genetic and Genomic Diversity

The current attrition of biodiversity extends beyond loss of species and unique populations to steady loss of a vast genomic diversity that remains largely undescribed. Yet the accelerating development of new techniques allows us to survey entire genomes ever faster and cheaper, to obtain robust samples from a diversity of sources including degraded DNA and residual DNA in the environment, and to address conservation efforts in new and innovative ways. Here we review recent studies that highlight the importance of carefully considering where to prioritize collection of genetic samples (e.g., organisms in rapidly changing landscapes or along edges of geographic ranges) and what samples to collect and archive (e.g., from individuals of little-known subspecies or populations, even of species not currently considered endangered). Those decisions will provide the sample infrastructure to detect the disappearance of certain genotypes or gene complexes, increases in inbreeding levels, and loss of genomic diversity as environmental conditions change. Obtaining samples from currently endangered, protected, and rare species can be particularly difficult, thus we also focus on studies that use new, non-invasive ways of obtaining genomic samples and analyzing them in these cases where other sampling options are highly constrained. Finally, biological collections archiving such samples face an inherent contradiction: their main goal is to preserve biological material in good shape so it can be used for scientific research for centuries to come, yet the technologies that can make use of such materials are advancing faster than collections can change their standardized practices. Thus, we also discuss current and potential new practices in biological collections that might bolster their usefulness for future biodiversity conservation research.

Testing which axes of species differentiation underlie covariance of phylogeographic similarity among montane sedge species

Co-distributed species may exhibit similar phylogeographic patterns due to shared environmental factors or discordant patterns attributed to the influence of species-specific traits. Although either concordant or discordant patterns could occur due to chance, stark differences in key traits (e.g., dispersal ability) may readily explain differences between species. Multiple species' attributes may affect genetic patterns, and it is difficult to isolate the contribution of each. Here we compare the relative importance of two attributes, range size, and niche breadth, in shaping the spatial structure of genetic variation in four sedge species (genus Carex) from the Rocky Mountains. Within two pairs of co-distributed species, one species exhibits narrow niche breadth, while the other species has broad niche breadth. Furthermore, one pair of co-distributed species has a large geographical distribution, while the other has a small distribution. The four species represent a natural experiment to tease apart how these attributes (i.e., range size and niche breadth) affect phylogeographic patterns. Investigations of genetic variation and structure revealed that range size, but not niche breadth, is related to spatial genetic covariation across species of montane sedges. Our study highlights how isolating key attributes across multiple species can inform their impact on processes driving intraspecific differentiation.

Closely related species show species-specific environmental responses and different spatial conservation needs: Prionailurus cats in the Indian subcontinent

Phylogenetically closely related species are often assumed to have similar responses to environmental conditions, but species-specific responses have also been described. These two scenarios may have different conservation implications. We tested these two hypotheses for Prionailurus cats (P. rubiginosus, P. bengalensis, P. viverrinus) in the Indian subcontinent and show its implications on species current protected area coverage and climatic suitability trends through time. We fitted ecological niche models with current environmental conditions and calculated niche overlap. In addition, we developed a model for the Jungle Cat Felis chaus to compare species responses and niche overlap estimates within Prionailurus with those for a related sympatric small cat species. Then we estimated the proportion of current suitable environment covered by protected area and projected climatic models from past (last interglacial) to future (2070; RCP4.5 and RCP8.5) conditions to show implications on population management and conservation. The hypothesis of a similar response and niche overlap among closely related species is not supported. Protected area coverage was lowest for P. viverrinus (mean = 0.071, SD = 0.012) and highest for P. bengalensis (mean = 0.088, SD = 0.006). In addition, the proportion of the subcontinent with suitable climate varied through time and was species-specific. For P. bengalensis, climatic suitability shrunk since at least the mid-Holocene, a trend that can be intensified by human-induced climate warming. Concerning P. viverrinus, most predictions show stable future climatic suitability, but a few indicated potential loss. Climatic suitability for P. rubiginous was predicted to remain stable but the species exhibited a negative association with intensive agriculture. Similar responses to environmental change by phylogenetically closely related species should not be assumed and have implications on protected area coverage and natural trends of species climatic suitability over time. This should be taken into account during conservation and management actions.

Trophic niche segregation allows range-extending coral reef fishes to co-exist with temperate species under climate change

Changing climate is forcing many terrestrial and marine species to extend their ranges poleward to stay within the bounds of their thermal tolerances. However, when such species enter higher latitude ecosystems, they engage in novel interactions with local species, such as altered predator-prey dynamics and competition for food. Here, we evaluate the trophic overlap between range-extending and local fish species along the east coast of temperate Australia, a hotspot for ocean warming and species range extensions. Stable isotope ratios (δ¹⁵ N and δ¹³ C) of muscle tissue and stomach content analysis were used to quantify overlap of trophic niche space between vagrant tropical and local temperate fish communities along a 730 km (6°) latitudinal gradient. Our study shows that in recipient temperate ecosystems, sympatric tropical and temperate species do not overlap significantly in their diet-even though they forage on broadly similar prey groups-and are therefore unlikely to compete for trophic niche space. The tropical and temperate species we studied, which are commonly found in shallow-water coastal environments, exhibited moderately broad niche breadths and local-scale dietary plasticity, indicating trophic generalism. We posit that because these species are generalists, they can co-exist under current climate change, facilitating the existence of novel community structures.

Multiple trans-Torres Strait colonisations by tree frogs in the Litoria caerulea group, with the description of a new species from New Guinea

Australia and New Guinea (together referred to as Sahul) were linked by land for much of the late Tertiary and share many biotic elements. However, New Guinea is dominated by rainforest, and northern Australia by savannah. Resolving patterns of biotic interchange between these two regions is critical to understanding the expansion and contraction of both habitat types. The green tree frog (Litoria caerulea) has a vast range across northern and eastern Australia and New Guinea. An assessment of mitochondrial and morphological diversity in this nominal taxon in New Guinea reveals two taxa. True Litoria caerulea occurs in disjunct savannahs of the Trans-Fly, Central Province and across northern Australia, with very low genetic divergence, implying late Pleistocene connectivity. A previously unrecognised taxon is endemic to New Guinea and widespread in lowland swampy rainforest. Date estimates for the divergence of the new species suggest Pliocene connectivity across lowland tropical habitats of northern Australia and New Guinea. In contrast, the new species shows shallow phylogeographic structuring across the central mountains of New Guinea, implying recent dispersal between the northern and southern lowlands. These results emphasise that the extent and connectivity of lowland rainforest and savannah environments across northern Australia and southern New Guinea have undergone profound shifts since the late Pliocene. http://zoobank.org/urn:lsid:zoobank.org:pub:A577A415-0B71-4663-B4C1-7271B97298CD

On and off the rocks: persistence and ecological diversification in a tropical Australian lizard radiation

Background

Congruent patterns in the distribution of biodiversity between regions or habitats suggest that key factors such as climatic and topographic variation may predictably shape evolutionary processes. In a number of tropical and arid biomes, genetic analyses are revealing deeper and more localised lineage diversity in rocky ranges than surrounding habitats. Two potential drivers of localised endemism in rocky areas are refugial persistence through climatic change, or ecological diversification and specialisation. Here we examine how patterns of lineage and phenotypic diversity differ across two broad habitat types (rocky ranges and open woodlands) in a small radiation of gecko lizards in the genus Gehyra (the australis group) from the Australian Monsoonal Tropics biome.

Results

Using a suite of approaches for delineating evolutionarily independent lineages, we find between 26 and 41 putative evolutionary units in the australis group (versus eight species currently recognised). Rocky ranges are home to a greater number of lineages that are also relatively more restricted in distribution, while lineages in open woodland habitats are fewer, more widely distributed, and, in one case, show evidence of range expansion. We infer at least two shifts out of rocky ranges and into surrounding woodlands. Phenotypic divergence between rocky ranges specialist and more generalist taxa is detected, but no convergent evolutionary regimes linked to ecology are inferred.

Conclusions

In climatically unstable biomes such as savannahs, rocky ranges have functioned as zones of persistence, generators of diversity and a source of colonists for surrounding areas. Phenotypic divergence can also be linked to the use of differing habitat types, however, the extent to which ecological specialisation is a primary driver or secondary outcome of localised diversification remains uncertain.

Electronic supplementary material

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

Unraveling the influences of climate change in Lepidosauria (Reptilia)

In recent decades, changes in climate have caused impacts on natural and human systems on all continents and across the oceans and many species have shifted their geographic ranges, seasonal activities, migration patterns, abundances and interactions in response to these changes. Projections of future climate change are uncertain, but the Earth's warming is likely to exceed 4.8 °C by the end of 21th century. The vulnerability of a population, species, group or system due to climate change is a function of impact of the changes on the evaluated system (exposure and sensitivity) and adaptive capacity as a response to this impact, and the relationship between these elements will determine the degree of species vulnerability. Predicting the potential future risks to biodiversity caused by climate change has become an extremely active field of research, and several studies in the last two decades had focused on determining possible impacts of climate change on Lepidosaurians, at a global, regional and local level. Here we conducted a systematic review of published studies in order to seek to what extent the accumulated knowledge currently allow us to identify potential trends or patterns regarding climate change effects on lizards, snakes, amphisbaenians and tuatara. We conducted a literature search among online literature databases/catalogues and recorded 255 studies addressing the influence of climate change on a total of 1918 species among 49 Lepidosaurian's families. The first study addressing this subject is dated 1999. Most of the studies focused on species distribution, followed by thermal biology, reproductive biology, behavior and genetics. We concluded that an integrative approach including most of these characteristics and also bioclimatic and environmental variables, may lead to consistent and truly effective strategies for species conservation, aiming to buffer the climate change effects on this group of reptiles.

Phylogenomics of a rapid radiation: the Australian rainbow skinks

Background

The application of target capture with next-generation sequencing now enables phylogenomic analyses of rapidly radiating clades of species. But such analyses are complicated by extensive incomplete lineage sorting, demanding the use of methods that consider this process explicitly, such as the multispecies coalescent (MSC) model. However, the MSC makes strong assumptions about divergence history and population structure, and when using the full Bayesian implementation, current computational limits mean that relatively few loci and samples can be analysed for even modest sized radiations. We explore these issues through analyses of an extensive (> 1000 loci) dataset for the Australian rainbow skinks. This group consists of 3 genera and 41 described species, which likely diversified rapidly in Australia during the mid-late Miocene to occupy rainforest, woodland, and rocky habitats with corresponding diversity of morphology and breeding colouration. Previous phylogenetic analyses of this group have revealed short inter-nodes and high discordance among loci, limiting the resolution of inferred trees. A further complication is that many species have deep phylogeographic structure – this poses the question of how to sample individuals within species for analyses using the MSC.

Results

Phylogenies obtained using concatenation and summary coalescent species tree approaches to the full dataset are well resolved with generally consistent topology, including for previously intractable relationships near the base of the clade. As expected, branch lengths at the tips are substantially overestimated using concatenation. Comparisons of different strategies for sampling haplotypes for full Bayesian MSC analyses (for one clade and using smaller sets of loci) revealed, unexpectedly, that combining haplotypes across divergent phylogeographic lineages yielded consistent species trees.

Conclusions

This study of more than 1000 loci provides a strongly-supported estimate of the phylogeny of the Australian rainbow skinks, which will inform future research on the evolution and taxonomy of this group. Our analyses suggest that species tree estimation with the MSC can be quite robust to violation of the assumption that the individuals representing a taxon are sampled from a panmictic population.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1130-4) contains supplementary material, which is available to authorized users.

Validation and description of two new north-western Australian Rainbow skinks with multispecies coalescent methods and morphology

While methods for genetic species delimitation have noticeably improved in the last decade, this remains a work in progress. Ideally, model based approaches should be applied and considered jointly with other lines of evidence, primarily morphology and geography, in an integrative taxonomy framework. Deep phylogeographic divergences have been reported for several species of Carlia skinks, but only for some eastern taxa have species boundaries been formally tested. The present study does this and revises the taxonomy for two species from northern Australia, Carlia johnstonei and C. triacantha. We introduce an approach that is based on the recently published method StarBEAST2, which uses multilocus data to explore the support for alternative species delimitation hypotheses using Bayes Factors (BFD). We apply this method, jointly with two other multispecies coalescent methods, using an extensive (from 2,163 exons) data set along with measures of 11 morphological characters. We use this integrated approach to evaluate two new candidate species previously revealed in phylogeographic analyses of rainbow skinks (genus Carlia) in Western Australia. The results based on BFD StarBEAST2, BFD* SNAPP and BPP genetic delimitation, together with morphology, support each of the four recently identified Carlia lineages as separate species. The BFD StarBEAST2 approach yielded results highly congruent with those from BFD* SNAPP and BPP. This supports use of the robust multilocus multispecies coalescent StarBEAST2 method for species delimitation, which does not require a priori resolved species or gene trees. Compared to the situation in C. triacantha, morphological divergence was greater between the two lineages within Kimberley endemic C. johnstonei, which also had deeper divergent histories. This congruence supports recognition of two species within C. johnstonei. Nevertheless, the combined evidence also supports recognition of two taxa within the more widespread C. triacantha. With this work, we describe two new species, Carlia insularis sp. nov and Carlia isostriacantha sp. nov. in the northwest of Australia. This contributes to increasing recognition that this region of tropical Australia has a rich and unique fauna.