Go west: Population genomics reveals unexpected population fluctuations and little gene flow in Western hemisphere populations of the predatory lady beetle, Hippodamia convergens.

Hippodamia convergens-the convergent lady beetle, has been used extensively in augmentative biological control of aphids, thrips, and whiteflies across its native range in North America, and was introduced into South America in the 1950s. Overwintering H. convergens populations from its native western range in the United States are commercially collected and released across its current range in the eastern USA, with little knowledge of the effectiveness of its augmentative biological control. Here we use a novel ddRADseq-based SNP/haplotype discovery approach to estimate its range-wide population diversity, differentiation, and recent evolutionary history. Our results indicate (1) significant population differentiation among eastern USA, western USA, and South American populations of H. convergens, with (2) little to no detectable recent admixture between them, despite repeated population augmentation, and (3) continued recent population size expansion across its range. These results contradict previous findings using microsatellite markers. In light of these new findings, the implications for the effectiveness of augmentative biological control using H. convergens are discussed. Additionally, because quantifying the non-target effects of augmentative biological control is a difficult problem in migratory beetles, our results could serve as a cornerstone in improving and predicting the efficacy of future releases of H. convergens across its range.

The phylogenomic and biogeographic history of the gartersnakes, watersnakes, and allies (Natricidae: Thamnophiini).

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using âˆ¼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

A Pervasive History of Gene Flow in Madagascar's True Lemurs (Genus Eulemur).

In recent years, it has become widely accepted that interspecific gene flow is common across the Tree of Life. Questions remain about how species boundaries can be maintained in the face of high levels of gene flow and how phylogeneticists should account for reticulation in their analyses. The true lemurs of Madagascar (genus Eulemur, 12 species) provide a unique opportunity to explore these questions, as they form a recent radiation with at least five active hybrid zones. Here, we present new analyses of a mitochondrial dataset with hundreds of individuals in the genus Eulemur, as well as a nuclear dataset containing hundreds of genetic loci for a small number of individuals. Traditional coalescent-based phylogenetic analyses of both datasets reveal that not all recognized species are monophyletic. Using network-based approaches, we also find that a species tree containing between one and three ancient reticulations is supported by strong evidence. Together, these results suggest that hybridization has been a prominent feature of the genus Eulemur in both the past and present. We also recommend that greater taxonomic attention should be paid to this group so that geographic boundaries and conservation priorities can be better established.

Lemur Gut Microeukaryotic Community Variation Is Not Associated with Host Phylogeny, Diet, or Habitat.

Identifying the major forces driving variation in gut microbiomes enhances our understanding of how and why symbioses between hosts and microbes evolved. Gut prokaryotic community variation is often closely associated with host evolutionary and ecological variables. Whether these same factors drive variation in other microbial taxa occupying the animal gut remains largely untested. Here, we present a one-to-one comparison of gut prokaryotic (16S rRNA metabarcoding) and microeukaryotic (18S rRNA metabarcoding) community patterning among 12 species of wild lemurs. Lemurs were sampled from dry forests and rainforests of southeastern Madagascar and display a range of phylogenetic and ecological niche diversity. We found that while lemur gut prokaryotic community diversity and composition vary with host taxonomy, diet, and habitat, gut microeukaryotic communities have no detectable association with any of these factors. We conclude that gut microeukaryotic community composition is largely random, while gut prokaryotic communities are conserved among host species. It is likely that a greater proportion of gut microeukaryotic communities comprise taxa with commensal, transient, and/or parasitic symbioses compared with gut prokaryotes, many of which form long-term relationships with the host and perform important biological functions. Our study highlights the importance of greater specificity in microbiome research; the gut microbiome contains many "omes" (e.g., prokaryome, eukaryome), each comprising different microbial taxa shaped by unique selective pressures.

Significant effects of host dietary guild and phylogeny in wild lemur gut microbiomes.

Mammals harbor diverse gut microbiomes (GMs) that perform critical functions for host health and fitness. Identifying factors associated with GM variation can help illuminate the role of microbial symbionts in mediating host ecological interactions and evolutionary processes, including diversification and adaptation. Many mammals demonstrate phylosymbiosis-a pattern in which more closely-related species harbor more similar GMs-while others show overwhelming influences of diet and habitat. Here, we generated 16S rRNA sequence data from fecal samples of 15 species of wild lemurs across southern Madagascar to (1) test a hypothesis of phylosymbiosis, and (2) test trait correlations between dietary guild, habitat, and GM diversity. Our results provide strong evidence of phylosymbiosis, though some closely-related species with substantial ecological niche overlap exhibited greater GM similarity than expected under Brownian motion. Phylogenetic regressions also showed a significant correlation between dietary guild and UniFrac diversity, but not Bray-Curtis or Jaccard. This discrepancy between beta diversity metrics suggests that older microbial clades have stronger associations with diet than younger clades, as UniFrac weights older clades more heavily. We conclude that GM diversity is predominantly shaped by host phylogeny, and that microbes associated with diet were likely acquired before evolutionary radiations within the lemur families examined.

Geography is more important than life history in the recent diversification of the tiger salamander complex.

The North American tiger salamander species complex, including its best-known species, the Mexican axolotl, has long been a source of biological fascination. The complex exhibits a wide range of variation in developmental life history strategies, including populations and individuals that undergo metamorphosis; those able to forego metamorphosis and retain a larval, aquatic lifestyle (i.e., paedomorphosis); and those that do both. The evolution of a paedomorphic life history state is thought to lead to increased population genetic differentiation and ultimately reproductive isolation and speciation, but the degree to which it has shaped population- and species-level divergence is poorly understood. Using a large multilocus dataset from hundreds of samples across North America, we identified genetic clusters across the geographic range of the tiger salamander complex. These clusters often contain a mixture of paedomorphic and metamorphic taxa, indicating that geographic isolation has played a larger role in lineage divergence than paedomorphosis in this system. This conclusion is bolstered by geography-informed analyses indicating no effect of life history strategy on population genetic differentiation and by model-based population genetic analyses demonstrating gene flow between adjacent metamorphic and paedomorphic populations. This fine-scale genetic perspective on life history variation establishes a framework for understanding how plasticity, local adaptation, and gene flow contribute to lineage divergence. Many members of the tiger salamander complex are endangered, and the Mexican axolotl is an important model system in regenerative and biomedical research. Our results chart a course for more informed use of these taxa in experimental, ecological, and conservation research.

Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life.

Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent.

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Insights from Population Genomics to Enhance and Sustain Biological Control of Insect Pests.

Biological control-the use of organisms (e.g., nematodes, arthropods, bacteria, fungi, viruses) for the suppression of insect pest species-is a well-established, ecologically sound and economically profitable tactic for crop protection. This approach has served as a sustainable solution for many insect pest problems for over a century in North America. However, all pest management tactics have associated risks. Specifically, the ecological non-target effects of biological control have been examined in numerous systems. In contrast, the need to understand the short- and long-term evolutionary consequences of human-mediated manipulation of biological control organisms for importation, augmentation and conservation biological control has only recently been acknowledged. Particularly, population genomics presents exceptional opportunities to study adaptive evolution and invasiveness of pests and biological control organisms. Population genomics also provides insights into (1) long-term biological consequences of releases, (2) the ecological success and sustainability of this pest management tactic and (3) non-target effects on native species, populations and ecosystems. Recent advances in genomic sequencing technology and model-based statistical methods to analyze population-scale genomic data provide a much needed impetus for biological control programs to benefit by incorporating a consideration of evolutionary consequences. Here, we review current technology and methods in population genomics and their applications to biological control and include basic guidelines for biological control researchers for implementing genomic technology and statistical modeling.

A spatial genomic approach identifies time lags and historical barriers to gene flow in a rapidly fragmenting Appalachian landscape.

The resolution offered by genomic data sets coupled with recently developed spatially informed analyses are allowing researchers to quantify population structure at increasingly fine temporal and spatial scales. However, both empirical research and conservation measures have been limited by questions regarding the impacts of data set size, data quality thresholds and the timescale at which barriers to gene flow become detectable. Here, we used restriction site associated DNA sequencing to generate a 2,140 single nucleotide polymorphism (SNP) data set for the copperhead snake (Agkistrodon contortrix) and address the population genomic impacts of recent and widespread landscape modification across an ~1,000-km2 region of eastern Kentucky, USA. Nonspatial population-based assignment and clustering methods supported little to no population structure. However, using individual-based spatial autocorrelation approaches we found evidence for genetic structuring which closely follows the path of a historically important highway which experienced high traffic volumes from c. 1920 to 1970 before losing most traffic to a newly constructed alternative route. We found no similar spatial genomic signatures associated with more recently constructed highways or surface mining activity, although a time lag effect may be responsible for the lack of any emergent spatial genetic patterns. Subsampling of our SNP data set suggested that similar results could be obtained with as few as 250 SNPs, and a range of thresholds for missing data exhibited limited impacts on the spatial patterns we detected. While we were not able to estimate relative effects of land uses or precise time lags, our findings highlight the importance of temporal factors in landscape genetics approaches, and suggest the potential advantages of genomic data sets and fine-scale, spatially informed approaches for quantifying subtle genetic patterns in temporally complex landscapes.

Comparative Genomic Analysis of the Pheromone Receptor Class 1 Family (V1R) Reveals Extreme Complexity in Mouse Lemurs (Genus, Microcebus) and a Chromosomal Hotspot across Mammals.

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R "hotspot."

Extensive variability in the gut microbiome of a highly-specialized and critically endangered lemur species across sites.

Deforestation continues to jeopardize Malagasy primates as viable habitats become smaller, more fragmented, and more disturbed. This deforestation can lead to changes in diet, microhabitat, and gene flow between populations of endangered species, and it remains unclear how these changes may affect gut microbiome (GM) characteristics. The black-and-white ruffed lemur (Varecia variegata), which is among Madagascar's most threatened lemur species, provides a critical model for understanding the relationships between historical and on-going deforestation (habitat disturbance), feeding ecology, and GM composition and diversity. We studied four populations inhabiting two rainforests (relatively pristine vs. highly disturbed) in southeastern Madagascar. We conducted full-day focal animal behavioral follows and collected fecal samples opportunistically across a three-month period. Our results indicate that lemurs inhabiting sites characterized by habitat disturbance and low dietary diversity exhibited reduced gut microbial alpha diversity. We also show that these same factors were associated with high community dissimilarity using weighted and unweighted UniFrac metrics. Finally, an indicator species analysis showed that the most pristine site was characterized by an abundance of methanogenic archaea. While it is impossible to disentangle the relative contributions of each confounding variable presented by our sampling design, these results provide crucial information about GM variability, thereby underscoring the importance of monitoring endangered species at the population-level.

Genomic Data Reveal Conserved Female Heterogamety in Giant Salamanders with Gigantic Nuclear Genomes.

Systems of genetic sex determination and the homology of sex chromosomes in different taxa vary greatly across vertebrates. Much progress remains to be made in understanding systems of genetic sex determination in non-model organisms, especially those with homomorphic sex chromosomes and/or large genomes. We used reduced representation genome sequencing to investigate genetic sex determination systems in the salamander family Cryptobranchidae (genera Cryptobranchus and Andrias), which typifies both of these inherent difficulties. We tested hypotheses of male- or female-heterogamety by sequencing hundreds of thousands of anonymous genomic regions in a panel of known-sex cryptobranchids and characterized patterns of presence/absence, inferred zygosity, and depth of coverage to identify sex-linked regions of these 56 gigabase genomes. Our results strongly support the hypothesis that all cryptobranchid species possess homologous systems of female heterogamety, despite maintenance of homomorphic sex chromosomes over nearly 60 million years. Additionally, we report a robust, non-invasive genetic assay for sex diagnosis in Cryptobranchus and Andrias which may have great utility for conservation efforts with these endangered salamanders. Co-amplification of these W-linked markers in both cryptobranchid genera provides evidence for long-term sex chromosome stasis in one of the most divergent salamander lineages. These findings inform hypotheses about the ancestral mode of sex determination in salamanders, but suggest that comparative data from other salamander families are needed. Our results further demonstrate that massive genomes are not necessarily a barrier to effective genome-wide sequencing and that the resulting data can be highly informative about sex determination systems in taxa with homomorphic sex chromosomes.

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet.

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep-sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep-sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back-arc basins, currently of great interest for deep-sea mineral extraction. A total of 42 loci were sequenced from each individual using high-throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.

Long-distance dispersal, ice sheet dynamics and mountaintop isolation underlie the genetic structure of glacier ice worms.

Disentangling the contemporary and historical factors underlying the spatial distributions of species is a central goal of biogeography. For species with broad distributions but little capacity to actively disperse, disconnected geographical distributions highlight the potential influence of passive, long-distance dispersal (LDD) on their evolutionary histories. However, dispersal alone cannot completely account for the biogeography of any species, and other factors-e.g. habitat suitability, life history-must also be considered. North American ice worms ( Mesenchytraeus solifugus) are ice-obligate annelids that inhabit coastal glaciers from Oregon to Alaska. Previous studies identified a complex biogeographic history for ice worms, with evidence for genetic isolation, unexpectedly close relationships among geographically disjunct lineages, and contemporary migration across large (e.g. greater than 1500 km) areas of unsuitable habitat. In this study, we analysed genome-scale sequence data for individuals from most of the known ice worm range. We found clear support for divergence between populations along the Pacific Coast and the inland flanks of the Coast Mountains (mean FST = 0.60), likely precipitated by episodic ice sheet expansion and contraction during the Pleistocene. We also found support for LDD of ice worms from Alaska to Vancouver Island, perhaps mediated by migrating birds. Our results highlight the power of genomic data for disentangling complex biogeographic patterns, including the presence of LDD.

Microbial assemblages reflect environmental heterogeneity in alpine streams.

Alpine streams are dynamic habitats harboring substantial biodiversity across small spatial extents. The diversity of alpine stream biota is largely reflective of environmental heterogeneity stemming from varying hydrological sources. Globally, alpine stream diversity is under threat as meltwater sources recede and stream conditions become increasingly homogeneous. Much attention has been devoted to macroinvertebrate diversity in alpine headwaters, yet to fully understand the breadth of climate change threats, a more thorough accounting of microbial diversity is needed. We characterized microbial diversity (specifically Bacteria and Archaea) of 13 streams in two disjunct Rocky Mountain subranges through 16S rRNA gene sequencing. Our study encompassed the spectrum of alpine stream sources (glaciers, snowfields, subterranean ice, and groundwater) and three microhabitats (ice, biofilms, and streamwater). We observed no difference in regional (γ) diversity between subranges but substantial differences in diversity among (ß) stream types and microhabitats. Within-stream (α) diversity was highest in groundwater-fed springs, lowest in glacier-fed streams, and positively correlated with water temperature for both streamwater and biofilm assemblages. We identified an underappreciated alpine stream type-the icy seep-that are fed by subterranean ice, exhibit cold temperatures (summer mean <2°C), moderate bed stability, and relatively high conductivity. Icy seeps will likely be important for combatting biodiversity losses as they contain similar microbial assemblages to streams fed by surface ice yet may be buffered against climate change by insulating debris cover. Our results show that the patterns of microbial diversity support an ominous trend for alpine stream biodiversity; as meltwater sources decline, stream communities will become more diverse locally, but regional diversity will be lost. Icy seeps, however, represent a source of optimism for the future of biodiversity in these imperiled ecosystems.

Genomic data reject the hypothesis of sympatric ecological speciation in a clade of Desmognathus salamanders.

Closely related taxa with dissimilar morphologies are often considered to have diverged via natural selection favoring different phenotypes. However, some studies have found these scenarios to be paired with limited or no genetic differentiation. Desmognathus quadramaculatus and D. marmoratus are sympatric salamander species thought to represent a case of ecological speciation based on distinct morphologies, but the results of previous studies have not resolved corresponding patterns of lineage divergence. Here, we use genome-wide data to test this hypothesis of ecological speciation. Population structure analyses partitioned individuals geographically, but not morphologically, into two adjacent regions of western North Carolina: Pisgah and Nantahala. Phylogenetic analyses confirmed the nominal species are nonmonophyletic and resolved deep divergence between the two geographic clusters. Model-testing overwhelmingly supported the hypothesis that lineage divergence followed geography. Finally, ecological niche modeling showed that Pisgah and Nantahala individuals occupy different climatic niches, and geographic boundaries for the two lineages correspond to differences in precipitation regimes across southern Appalachia. Overall, we reject the previous hypothesis of ecological speciation based on microhabitat partitioning. Instead, our results suggest that there are two cryptic lineages, each containing the same pair of morphotypes.

Estimation of contemporary effective population size and population declines using RAD sequence data.

Large genomic data sets generated with restriction site-associated DNA sequencing (RADseq), in combination with demographic inference methods, are improving our ability to gain insights into the population history of species. We used a simulation approach to examine the potential for RADseq data sets to accurately estimate effective population size (N e) over the course of stable and declining population trends, and we compare the ability of two methods of analysis to accurately distinguish stable from steadily declining populations over a contemporary time scale (20 generations). Using a linkage disequilibrium-based analysis, individual sampling (i.e., n ≥ 30) had the greatest effect on N e estimation and the detection of population size declines, with declines reliably detected across scenarios ~10 generations after they began. Coalescent-based inference required fewer sampled individuals (i.e., n = 15), and instead was most influenced by the size of the SNP data set, with 25,000-50,000 SNPs required for accurate detection of population trends and at least 20 generations after decline began. The number of samples available and targeted number of RADseq loci are important criteria when choosing between these methods. Neither method suffered any apparent bias due to the effects of allele dropout typical of RAD data. With an understanding of the limitations and biases of these approaches, researchers can make more informed decisions when designing their sampling and analyses. Overall, our results reveal that demographic inference using RADseq data can be successfully applied to infer recent population size change and may be an important tool for population monitoring and conservation biology.

Climate change and alpine stream biology: progress, challenges, and opportunities for the future.

In alpine regions worldwide, climate change is dramatically altering ecosystems and affecting biodiversity in many ways. For streams, receding alpine glaciers and snowfields, paired with altered precipitation regimes, are driving shifts in hydrology, species distributions, basal resources, and threatening the very existence of some habitats and biota. Alpine streams harbour substantial species and genetic diversity due to significant habitat insularity and environmental heterogeneity. Climate change is expected to affect alpine stream biodiversity across many levels of biological resolution from micro- to macroscopic organisms and genes to communities. Herein, we describe the current state of alpine stream biology from an organism-focused perspective. We begin by reviewing seven standard and emerging approaches that combine to form the current state of the discipline. We follow with a call for increased synthesis across existing approaches to improve understanding of how these imperiled ecosystems are responding to rapid environmental change. We then take a forward-looking viewpoint on how alpine stream biologists can make better use of existing data sets through temporal comparisons, integrate remote sensing and geographic information system (GIS) technologies, and apply genomic tools to refine knowledge of underlying evolutionary processes. We conclude with comments about the future of biodiversity conservation in alpine streams to confront the daunting challenge of mitigating the effects of rapid environmental change in these sentinel ecosystems.