Evolutionary diversification of cryophilic Grylloblatta species (Grylloblattodea: Grylloblattidae) in alpine habitats of California

Environmental and Biogeographic Drivers behind Alpine Plant Thermal Tolerance and Genetic Variation

In alpine ecosystems, elevation broadly functions as a steep thermal gradient, with plant communities exposed to regular fluctuations in hot and cold temperatures. These conditions lead to selective filtering, potentially contributing to species-level variation in thermal tolerance and population-level genetic divergence. Few studies have explored the breadth of alpine plant thermal tolerances across a thermal gradient or the underlying genetic variation thereof. We measured photosystem heat (Tcrit-hot) and cold (Tcrit-cold) thresholds of ten Australian alpine species across elevation gradients and characterised their neutral genetic variation. To reveal the biogeographical drivers of present-day genetic signatures, we also reconstructed temporal changes in habitat suitability across potential distributional ranges. We found intraspecific variation in thermal thresholds, but this was not associated with elevation, nor underpinned by genetic differentiation on a local scale. Instead, regional population differentiation and considerable homozygosity within populations may, in part, be driven by distributional contractions, long-term persistence, and migrations following habitat suitability. Our habitat suitability models suggest that cool-climate-distributed alpine plants may be threatened by a warming climate. Yet, the observed wide thermal tolerances did not reflect this vulnerability. Conservation efforts should seek to understand variations in species-level thermal tolerance across alpine microclimates.

Phylogenomics of American pika (Ochotona princeps) lineage diversification

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Phylogeography of an endemic California silkmoth genus suggests the importance of an unheralded central California province in generating regional endemic biodiversity

The California Floristic province is a biodiversity hotspot. Understanding the phylogeographic patterns that exist in this diverse region is essential to understanding its evolution and for guiding conservation efforts. Calosaturnia is a charismatic silkmoth genus endemic to large portions of the region with three described species, C. mendocino, C. walterorum, and C. albofasciata. We sampled all three species from across their ranges, sequenced 1463 bp of mitochondrial COI and 1941 bp of nuclear DNA from three genes, and reconstructed phylogenetic relationships and estimated divergence times within the lineages. All three species show pronounced evidence of isolation and, in two cases, secondary reconnection. An unexpected monophyletic mtDNA lineage was found in the Central Coast region, in a region thought to represent an intergrade between C. mendocino and C. walterorum. Our genetic data also significantly revise previous hypotheses as to which species occur in which regions, suggesting that historical ecological changes around four Ma ago isolated some lineages, and a secondary isolation event two Ma ago led to isolation of populations both in the Central Coast region and the southern Sierra Nevada. Our results add to a currently under-appreciated pattern suggesting that coastal Central California is not a transition zone between Northern and Southern California Floristic Province faunas but rather its own unique, periodically isolated, biogeographic region. They also suggest cryptic diversity may be present in many other currently unrecognized groups. Additional research should focus on this central California region because many species may be highly restricted in range and in need of conservation attention.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.

Rapid speciation and ecological divergence into North American alpine habitats: the Nippononebria (Coleoptera: Carabidae) species complex

The climate-driven species pump hypothesis has been supported in a number of phylogeographic studies of alpine species. Climate-driven shifts in distribution, coupled with rapid demographic change, have led to strong genetic drift and lineage diversification. Although the species pump has been linked to rapid speciation in a number of studies, few studies have demonstrated that ecological divergence accompanies rapid speciation. Here we examine genetic, morphological and physiological variation in members of the ground beetle taxon Nippononebria, to test three competing hypotheses of evolutionary diversification: isolation and incomplete lineage sorting (no speciation), recent speciation without ecological divergence, or recent speciation with ecological divergence into alpine habitats. Genetic data are consistent with recent divergence, with major lineages forming in the last million years. A species tree analysis, in conjunction with morphological divergence in male reproductive traits, support the formation of three recognized Nippononebria taxa. Furthermore, both morphological and physiological traits demonstrate ecological divergence in alpine lineages, with convergent shifts in body shape and thermal tolerance breadth. This provides strong evidence that the climate-driven species pump can generate ecological novelty, though it is argued that spatial scale may be a key determinant of broader patterns of macroevolution in alpine communities.

Grylloblattodea of Canada

The enigmatic insect order Grylloblattodea comprises two described species in Canada, which are limited to the Montane Cordillera and Pacific Maritime ecozones. One of the described species has three Canadian subspecies of uncertain taxonomic ranking, and there are two additional undescribed or unreported species known in close proximity to the Canadian border in western Alberta and British Columbia that likely also occur in Canada. Thus, as much as 50% of the total taxonomic diversity of Grylloblattodea in Canada is still undocumented. Targeted surveys and taxonomic work, as well as studies that describe the ecology and conservation status of Grylloblattodea are important goals for future research.

Testing models of refugial isolation, colonization and population connectivity in two species of montane salamanders

Pleistocene glaciations have profoundly affected patterns of genetic diversity within many species. Temperate alpine organisms likely experienced dramatic range shifts, given that much of their habitat was glaciated during this time. While the effects of glaciations are relatively well understood, the spatial locations of refugia and processes that gave rise to current patterns of diversity are less well known. We use a microsatellite data set to test hypotheses of population connectivity and refugial isolation in the web-toed salamanders (Hydromantes) of the Sierra Nevada. We reject models of refugia with subsequent expansion into either the high southern Sierra or low-elevation Owens Valley, in favor of a simple isolation model with no migration between current populations. We find no evidence of migration at even moderate spatial scales using a variety of analyses in the southern Sierra, and limited migration in the northern Sierra. These results suggest that divergence in isolation following fragmentation is the dominant process structuring genetic variation in these salamander species. In the context of anthropogenic climate change and habitat degradation, these results imply that salamanders and other low-vagility alpine organisms are at risk of decline as they are unlikely to migrate across large distances.

The Jurassic Bajanzhargalanidae (Insecta: Grylloblattida?): New genera and species, and data on postabdominal morphology

The presumed phylogenetic link between extant ice-crawlers (Grylloblattidae = 'crown-Grylloblattida') and fossil species of the taxon concept Grylloblattida sensu Storozhenko (2002) is essentially based on postabdominal morphology. However, the fossil data are limited, and the interpretation is open to debate. Here we investigate a sample of a poorly known fossil 'grylloblattidan' family, the Bajanzhargalanidae, collected from the Daohugou locality (Middle Jurassic, China). We describe Sinonele fangi gen. nov., sp. nov., Sinonele hei gen. nov., sp. nov., Sinonele phasmoides gen. nov., sp. nov., and Sinonele mini gen. nov., sp. nov. Thanks to the abundance and exceptional preservation of the material, we could document wing venation intra-specific variability, provide cues to identify male and female individuals, describe and tentatively interpret various body structures of both sexes, and discuss them with a broad pterygotan phylogenetic perspective. The Bajanzhargalanidae exhibit a puzzling combination of postabdominal characters leaving us inconclusive on their affinities, or lack thereof, with crown-Grylloblattida. Our contribution suggests that a substantial effort will be needed to further investigate postabdominal structures from comparatively ancient fossil insects preserved as rock imprints, because of their broad morphological disparity.

Conserved and narrow temperature limits in alpine insects: Thermal tolerance and supercooling points of the ice-crawlers, Grylloblatta (Insecta: Grylloblattodea: Grylloblattidae)

For many terrestrial species, habitat associations and range size are dependent on physiological limits, which in turn may influence large-scale patterns of species diversity. The temperature range experienced by individuals is considered to shape the breadth of the thermal niche, with species occupying temporally and/or geographically stable climates tolerating a narrow temperature range. High-elevation environments experience large temperature fluctuations, with frequent periods below 0 °C, but Grylloblatta (Grylloblattodea: Grylloblattidae) occupy climatically stable microhabitats within this region. Here we test critical thermal limits and supercooling points for five Grylloblatta populations from across a large geographic area, to examine whether the stable microhabitats of this group are associated with a narrow thermal niche and assess their capacity to tolerate cold conditions. Thermal limits are highly conserved in Grylloblatta, despite substantial genetic divergence among populations spanning 1500 m elevation and being separated by over 500 km. Further, Grylloblatta show exceptionally narrow thermal limits compared to other insect taxa with little capacity to improve cold tolerance via plasticity. In contrast, upper thermal limits were significantly depressed by cold acclimation. Grylloblatta maintain coordinated movement until they freeze, and they die upon freezing. Convergence of the critical thermal minima, supercooling point and lower lethal limits point to adaptation to a cold but, importantly, constant thermal environment. These physiological data provide an explanation for the high endemism and patchy distribution of Grylloblatta, which relies on subterranean retreats to accommodate narrow thermal limits. These retreats are currently buffered from temperature fluctuations by snow cover, and a declining snowpack thus places Grylloblatta at risk of exposure to temperatures beyond its tolerance capacity.

Sky island diversification meets the multispecies coalescent - divergence in the spruce-fir moss spider (Microhexura montivaga, Araneae, Mygalomorphae) on the highest peaks of southern Appalachia

Microhexura montivaga is a miniature tarantula-like spider endemic to the highest peaks of the southern Appalachian mountains and is known only from six allopatric, highly disjunct montane populations. Because of severe declines in spruce-fir forest in the late 20th century, M. montivaga was formally listed as a US federally endangered species in 1995. Using DNA sequence data from one mitochondrial and seven nuclear genes, patterns of multigenic genetic divergence were assessed for six montane populations. Independent mitochondrial and nuclear discovery analyses reveal obvious genetic fragmentation both within and among montane populations, with five to seven primary genetic lineages recovered. Multispecies coalescent validation analyses [guide tree and unguided Bayesian Phylogenetics and Phylogeography (BPP), Bayes factor delimitation (BFD)] using nuclear-only data congruently recover six or seven distinct lineages; BFD analyses using combined nuclear plus mitochondrial data favour seven or eight lineages. In stark contrast to this clear genetic fragmentation, a survey of secondary sexual features for available males indicates morphological conservatism across montane populations. While it is certainly possible that morphologically cryptic speciation has occurred in this taxon, this system may alternatively represent a case where extreme population genetic structuring (but not speciation) leads to an oversplitting of lineage diversity by multispecies coalescent methods. Our results have clear conservation implications for this federally endangered taxon and illustrate a methodological issue expected to become more common as genomic-scale data sets are gathered for taxa found in naturally fragmented habitats.

Strong genetic structure corresponds to small-scale geographic breaks in the Australian alpine grasshopper Kosciuscola tristis

Background

Mountain landscapes are topographically complex, creating discontinuous ‘islands’ of alpine and sub-alpine habitat with a dynamic history. Changing climatic conditions drive their expansion and contraction, leaving signatures on the genetic structure of their flora and fauna. Australia’s high country covers a small, highly fragmented area. Although the area is thought to have experienced periods of relative continuity during Pleistocene glacial periods, small-scale studies suggest deep lineage divergence across low-elevation gaps. Using both DNA sequence data and microsatellite markers, we tested the hypothesis that genetic partitioning reflects observable geographic structuring across Australia’s mainland high country, in the widespread alpine grasshopper Kosciuscola tristis (Sjösted).

Results

We found broadly congruent patterns of regional structure between the DNA sequence and microsatellite datasets, corresponding to strong divergence among isolated mountain regions. Small and isolated mountains in the south of the range were particularly distinct, with well-supported divergence corresponding to climate cycles during the late Pliocene and Pleistocene. We found mixed support, however, for divergence among other mountain regions. Interestingly, within areas of largely contiguous alpine and sub-alpine habitat around Mt Kosciuszko, microsatellite data suggested significant population structure, accompanied by a strong signature of isolation-by-distance.

Conclusions

Consistent patterns of strong lineage divergence among different molecular datasets indicate genetic breaks between populations inhabiting geographically distinct mountain regions. Three primary phylogeographic groups were evident in the highly fragmented Victorian high country, while within-region structure detected with microsatellites may reflect more recent population isolation. Despite the small area of Australia’s alpine and sub-alpine habitats, their low topographic relief and lack of extensive glaciation, divergence among populations was on the same scale as that detected in much more extensive Northern hemisphere mountain systems. The processes driving divergence in the Australian mountains might therefore differ from their Northern hemisphere counterparts.

Electronic supplementary material

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

Australia's arid-adapted butcherbirds experienced range expansions during Pleistocene glacial maxima

A model of range expansions during glacial maxima (GM) for cold-adapted species is generally accepted for the Northern Hemisphere. Given that GM in Australia largely resulted in the expansion of arid zones, rather than glaciation, it could be expected that arid-adapted species might have had expanded ranges at GM, as cold-adapted species did in the Northern Hemisphere. For Australian biota, however, it remains paradigmatic that arid-adapted species contracted to refugia at GM. Here we use multilocus data and ecological niche models (ENMs) to test alternative GM models for butcherbirds. ENMs, mtDNA and estimates of nuclear introgression and past population sizes support a model of GM expansion in the arid-tolerant Grey Butcherbird that resulted in secondary contact with its close relative--the savanna-inhabiting Silver-backed Butcherbird--whose contemporary distribution is widely separated. Together, these data reject the universal use of a GM contraction model for Australia's dry woodland and arid biota.

GenGIS 2: geospatial analysis of traditional and genetic biodiversity, with new gradient algorithms and an extensible plugin framework

GenGIS is free and open source software designed to integrate biodiversity data with a digital map and information about geography and habitat. While originally developed with microbial community analyses and phylogeography in mind, GenGIS has been applied to a wide range of datasets. A key feature of GenGIS is the ability to test geographic axes that can correspond to routes of migration or gradients that influence community similarity. Here we introduce GenGIS version 2, which extends the linear gradient tests introduced in the first version to allow comprehensive testing of all possible linear geographic axes. GenGIS v2 also includes a new plugin framework that supports the development and use of graphically driven analysis packages: initial plugins include implementations of linear regression and the Mantel test, calculations of alpha-diversity (e.g., Shannon Index) for all samples, and geographic visualizations of dissimilarity matrices. We have also implemented a recently published method for biomonitoring reference condition analysis (RCA), which compares observed species richness and diversity to predicted values to determine whether a given site has been impacted. The newest version of GenGIS supports vector data in addition to raster files. We demonstrate the new features of GenGIS by performing a full gradient analysis of an Australian kangaroo apple data set, by using plugins and embedded statistical commands to analyze human microbiome sample data, and by applying RCA to a set of samples from Atlantic Canada. GenGIS release versions, tutorials and documentation are freely available at http://kiwi.cs.dal.ca/GenGIS, and source code is available at https://github.com/beiko-lab/gengis.

Conceptual models of the processes and patterns of the ecological, evolutionary and bio-geographical consequences of global climate changes

A unifying conceptual model is constructed for the major effects of alternating periods of global warming and cooling and sea-level changes on the geographical distributions and the ecological and genetic characteristics of species and ecological communities.

The main results found are:

Colliding fragment islands transport independent lineages of endemic rock-crawlers (Grylloblattodea: Grylloblattidae) in the Japanese archipelago

Fragment islands, viewed from the paradigm of island biogeographic theory, depend on continual immigration from continental sources to maintain levels of species diversity, or otherwise undergo a period of relaxation where species diversity declines to a lower equilibrium. Japan is a recently derived fragment island with a rich endemic flora and fauna. These endemic species have been described as paleoendemics, and conversely as recently derived Pleistocene colonists. Geological events in the Miocene period, notably the fragmentation and collision of islands, and the subsequent uplift of mountains in central Japan, provided opportunities for genetic isolation. More recently, cyclical climatic change during the Pliocene and Pleistocene periods led to intermittent land bridge connections to continental Asia. Here we investigate the pattern and timing of diversification in a diverse endemic lineage in order to test whether ongoing migration has sustained species diversity, whether there is evidence of relaxation, and how geological and climatic events are associated with lineage diversification. Using multi-locus genetic data, we test these hypotheses in a poorly dispersing, cold-adapted terrestrial insect lineage (Grylloblattodea: Grylloblattidae) sampled from Japan, Korea, and Russia. In phylogenetic analyses of concatenated data and a species tree approach, we find evidence of three deeply divergent lineages of rock-crawlers in Japan consistent with the pattern of island fragmentation from continental Asia. Tests of lineage diversification rates suggest that relaxation has not occurred and instead endemism has increased in the Japanese Grylloblattidae following mountain-building events in the Miocene. Although the importance of climate change in generating species diversity is a commonly held paradigm in Japanese biogeography, our analyses, including analyses of demographic change and phylogeographic range shifts in putative species, suggests that Pleistocene climatic change has had a limited effect on the diversification of rock-crawlers.

Hybridization in East African swarm-raiding army ants

BACKGROUND

Hybridization can have complex effects on evolutionary dynamics in ants because of the combination of haplodiploid sex-determination and eusociality. While hybrid non-reproductive workers have been found in a range of species, examples of gene-flow via hybrid queens and males are rare. We studied hybridization in East African army ants (Dorylus subgenus Anomma) using morphology, mitochondrial DNA sequences, and nuclear microsatellites.

RESULTS

While the mitochondrial phylogeny had a strong geographic signal, different species were not recovered as monophyletic. At our main study site at Kakamega Forest, a mitochondrial haplotype was shared between a "Dorylus molestus-like" and a "Dorylus wilverthi-like" form. This pattern is best explained by introgression following hybridization between D. molestus and D. wilverthi. Microsatellite data from workers showed that the two morphological forms correspond to two distinct genetic clusters, with a significant proportion of individuals being classified as hybrids.

CONCLUSIONS

We conclude that hybridization and gene-flow between the two army ant species D. molestus and D. wilverthi has occurred, and that mating between the two forms continues to regularly produce hybrid workers. Hybridization is particularly surprising in army ants because workers have control over which males are allowed to mate with a young virgin queen inside the colony.