Genomic signatures of thermal adaptation are associated with clinal shifts of life history in a broadly distributed frog

Frankenstein matrices: Among-population life history variation affects the reliability and predictions of demographic models

Population matrix models are routinely used to study the demography of wild populations and to guide management choices. When vital rates are unknown for a specific population or life history stage, researchers often replace them with estimates from other populations of the same species. Such 'hybrid' matrices might ignore among-population life history variation and lead to incorrect inferences. In this study, we examined the real-world effect of using hybrid matrices on demographic inference and management decisions, using a large dataset on yellow-bellied toad (Bombina variegata) populations, an amphibian species whose life history depends on human land use. We estimated stage-specific survival and recruitment for 18 populations across different habitat types. We then assessed how estimated population growth rates and elasticities changed when population-specific vital rates were replaced by estimates from other populations, chosen randomly or based on habitat, demographic or geographic proximity. The use of hybrid matrices mixing demographic estimates from different populations and habitats biased predictions. The mean bias was relatively minor even when sampling randomly across all populations, because our large dataset represented the whole range of life histories and errors cancelled out on average. However, borrowing estimates from geographically close or demographically similar populations substantially reduced the risk of extreme errors. Borrowing from populations from similar habitat types could also reduce bias, but results varied depending on the exact habitat types concerned. Our study illustrates how habitat-specific among-population variation in life history affects the reliability of population matrices commonly used in evolutionary demography, ecology and conservation. When the use of hybrid population matrices cannot be avoided, their creation can be informed by additional information about ecological or demographic patterns, helping reduce bias. When such information is not available, we recommend that studies should consider the whole space of parameter estimates (the complete range of estimates available), thus transparently describing the true uncertainty surrounding demographic estimates.

Increased Temperature and Exposure to Ammonium Alter the Life Cycle of an Anuran Species

Amphibian populations are undergoing a major recession worldwide, likely triggered by global change components such as the global warming and pollutants, among which agrochemicals, in general, and fertilizers, in particular, play a central role given their relevance in agriculture. Potential synergies among these stressors could maximize their individual effects. In this work, we investigated the consequences of a controlled chronic exposure to increased temperature and a sublethal dose of ammonium during the larval stage of Pelophylax perezi frogs on the growth, development, and locomotor performance of tadpoles and the metamorphs they gave rise to. To that end, tadpoles were reared either in heated or nonheated tanks, with or without ammonium added. The parents of these tadpoles came from either a pine grove or an agrosystem. Survival was reduced in agrosystem tadpoles reared with ammonium. Increased temperature potentiated tadpole growth while giving way to smaller metamorphs. Faster growth could be a consequence of increased metabolism, whereas the smaller size could follow an accelerated development and metamorphosis, which reduced the growth period. Also, swimming speed was greater in tadpoles reared in heated tanks, while jumping distance was greater in metamorphs reared in nonheated tanks. The effects of temperature were more marked in agrosystem than in pine grove individuals, which could mirror reduced adaptability. Thus, the ability to withstand the effects of these stressors was lower in agrosystem tadpoles.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Polygenic adaptation of a cosmopolitan pest to a novel thermal environment

The fluctuation in temperature poses a significant challenge for poikilothermic organisms, notably insects, particularly in the context of changing climatic conditions. In insects, temperature adaptation has been driven by polygenes. In addition to genes that directly affect traits (core genes), other genes (peripheral genes) may also play a role in insect temperature adaptation. This study focuses on two peripheral genes, the GRIP and coiled-coil domain containing 2 (GCC2) and karyopherin subunit beta 1 (KPNB1). These genes are differentially expressed at different temperatures in the cosmopolitan pest, Plutella xylostella. GCC2 and KPNB1 in P. xylostella were cloned, and their relative expression patterns were identified. Reduced capacity for thermal adaptation (development, reproduction and response to temperature extremes) in the GCC2-deficient and KPNB1-deficient P. xylostella strains, which were constructed by CRISPR/Cas9 technique. Deletion of the PxGCC2 or PxKPNB1 genes in P. xylostella also had a differential effect on gene expression for many traits including stress resistance, resistance to pesticides, involved in immunity, trehalose metabolism, fatty acid metabolism and so forth. The ability of the moth to adapt to temperature via different pathways is likely to be key to its ability to remain an important pest species under predicted climate change conditions.

Genetic-environment associations explain genetic differentiation and variation between western and eastern North Pacific rhinoceros auklet (Cerorhinca monocerata) breeding colonies

Animals are strongly connected to the environments they live in and may become adapted to local environments. Examining genetic-environment associations of key indicator species, like seabirds, provides greater insights into the forces that drive evolution in marine systems. Here we examined a RADseq dataset of 19,213 SNPs for 99 rhinoceros auklets (Cerorhinca monocerata) from five western Pacific and 10 eastern Pacific breeding colonies. We used partial redundancy analyses to identify candidate adaptive loci and to quantify the effects of environmental variation on population genetic structure. We identified 262 candidate adaptive loci, which accounted for 3.0% of the observed genetic variation among western Pacific and eastern Pacific breeding colonies. Genetic variation was more strongly associated with pH and maximum current velocity, than maximum sea surface temperature. Genetic-environment associations explain genetic differences between western and eastern Pacific populations; however, genetic variation within the western and eastern Pacific Ocean populations appears to follow a pattern of isolation-by-distance. This study represents a first to quantify the relationship between environmental and genetic variation for this widely distributed marine species and provides greater insights into the evolutionary forces that act on marine species.

Putative climate adaptation in American pikas (Ochotona princeps) is associated with copy number variation across environmental gradients

Improved understanding of the genetic basis of adaptation to climate change is necessary for maintaining global biodiversity moving forward. Studies to date have largely focused on sequence variation, yet there is growing evidence that suggests that changes in genome structure may be an even more significant source of adaptive potential. The American pika (Ochotona princeps) is an alpine specialist that shows some evidence of adaptation to climate along elevational gradients, but previous work has been limited to single nucleotide polymorphism based analyses within a fraction of the species range. Here, we investigated the role of copy number variation underlying patterns of local adaptation in the American pika using genome-wide data previously collected across the entire species range. We identified 37-193 putative copy number variants (CNVs) associated with environmental variation (temperature, precipitation, solar radiation) within each of the six major American pika lineages, with patterns of divergence largely following elevational and latitudinal gradients. Genes associated (n = 158) with independent annotations across lineages, variables, and/or CNVs had functions related to mitochondrial structure/function, immune response, hypoxia, olfaction, and DNA repair. Some of these genes have been previously linked to putative high elevation and/or climate adaptation in other species, suggesting they may serve as important targets in future studies.

Climate change and its effects on body size and shape: the role of endocrine mechanisms

In many organisms, rapidly changing environmental conditions are inducing dramatic shifts in diverse phenotypic traits with consequences for fitness and population viability. However, the mechanisms that underlie these responses remain poorly understood. Endocrine signalling systems often influence suites of traits and are sensitive to changes in environmental conditions; they are thus ideal candidates for uncovering both plastic and evolved consequences of climate change. Here, we use body size and shape, a set of integrated traits predicted to shift in response to rising temperatures with effects on fitness, and insulin-like growth factor-1 as a case study to explore these ideas. We review what is known about changes in body size and shape in response to rising temperatures and then illustrate why endocrine signalling systems are likely to be critical in mediating these effects. Lastly, we discuss research approaches that will advance understanding of the processes that underlie rapid responses to climate change and the role endocrine systems will have. Knowledge of the mechanisms involved in phenotypic responses to climate change will be essential for predicting both the ecological and the long-term evolutionary consequences of a warming climate. This article is part of the theme issue 'Endocrine responses to environmental variation: conceptual approaches and recent developments'.

Genomic vulnerability of a freshwater salmonid under climate change

Understanding the adaptive potential of populations and species is pivotal for minimizing the loss of biodiversity in this era of rapid climate change. Adaptive potential has been estimated in various ways, including based on levels of standing genetic variation, presence of potentially beneficial alleles, and/or the severity of environmental change. Kokanee salmon, the non-migratory ecotype of sockeye salmon (Oncorhynchus nerka), is culturally and economically important and has already been impacted by the effects of climate change. To assess its climate vulnerability moving forward, we integrated analyses of standing genetic variation, genotype-environment associations, and climate modeling based on sequence and structural genomic variation from 224 whole genomes sampled from 22 lakes in British Columbia and Yukon (Canada). We found that variables for extreme temperatures, particularly warmer temperatures, had the most pervasive signature of selection in the genome and were the strongest predictors of levels of standing variation and of putatively adaptive genomic variation, both sequence and structural. Genomic offset estimates, a measure of climate vulnerability, were significantly correlated with higher increases in extreme warm temperatures, further highlighting the risk of summer heat waves that are predicted to increase in frequency in the future. Levels of standing genetic variation, an important metric for population viability and resilience, were not correlated with genomic offset. Nonetheless, our combined approach highlights the importance of integrating different sources of information and genomic data to formulate more comprehensive and accurate predictions on the vulnerability of populations and species to future climate change.

Widespread Deviant Patterns of Heterozygosity in Whole-Genome Sequencing Due to Autopolyploidy, Repeated Elements, and Duplication

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths.

Importance of Genetic-Fitness Correlations for the Conservation of Amphibians

Endangered animals suffer from isolation of their habitats. Isolation leads to a reduction in population size as well as a decrease in genetic diversity and a concomitant increase in the risk of extinction. Amphibians are the most endangered vertebrate class. Besides habitat loss, fragmentation and isolation, amphibians are threatened by emerging diseases e.g., chytrid fungus or Ranavirus. By employing experiments, researchers investigate whether changes in genetic diversity within or among isolated populations affect amphibian fitness. While genetic diversity estimates are based on molecular markers, typically microsatellites, fitness is mostly measured as tadpole performance in rearing experiments often under varying environmental conditions. Tadpole performances (e.g., body mass, growth rate and survival) have been found to be negatively affected by low genetic diversity, as several studies have found a positive association between genetic diversity and these fitness traits. Moreover, infection with pathogens also seems to be more likely in individuals or populations with lower genetic diversity. Overall, these genetic-fitness correlations seem to be more pronounced or detectable in smaller, declining populations but not in larger populations. Genomic studies, which sample a larger fraction of the genome, are still scarce in the conservation genetic literature on amphibians. These are likely to increase in upcoming years and may reveal adaptive variants that protect against dangerous pathogens or environmental changes. Altogether, genetic-fitness correlation studies should be a priority in order to develop effective management plans for the genetic rescue of isolated, imperilled amphibian populations.

Climate Adaptation and Drift Shape the Genomes of Two Eel-Goby Sister Species Endemic to Contrasting Latitude

Deciphering the role of climate adaptation in generating genetic divergence and hence speciation is a central question in evolution. Comparisons of genomes of closely related species spanning selective climate gradients are particularly informative in discerning the signatures of selection and thereby providing valuable information concerning the role of climate adaptation in speciation. Here we re-sequenced 99 genomes of the two sister eel-goby species Odontamblyopus lacepedii and O. rebecca, which are endemic to tidal mudflats spanning contrasting latitude gradients, to estimate the influence of divergent climate selection on shaping genome-wide patterns of divergence. The results indicated that genome-wide differentiation between the two species was evident (genome-wide FST = 0.313). Against a background of high baseline genomic divergence, 588 and 1202 elevated divergent loci were detected to be widespread throughout their genomes, as opposed to focused within small islands of genomic regions. These patterns of divergence may arise from divergent climate selection in addition to genetic drift acting through past glacial segregation (1.46 million years ago). We identified several candidate genes that exhibited elevated divergence between the two species, including genes associated with substance metabolism, energy production, and response to environmental cues, all putative candidates closely linked to thermal adaptation expected from the latitude gradient. Interestingly, several candidates related to gamete recognition and time of puberty, and also exhibited elevated divergence, indicating their possible role in pre-zygote isolation and speciation of the two species. Our results would expand our knowledge on the roles of latitude climate adaptation and genetic drift in generating and maintaining biodiversity in marine teleosts.

Slow and negligible senescence among testudines challenges evolutionary theories of senescence

Is senescence inevitable and universal for all living organisms, as evolutionary theories predict? Although evidence generally supports this hypothesis, it has been proposed that certain species, such as turtles and tortoises, may exhibit slow or even negligible senescence-i.e., avoiding the increasing risk of death from gradual deterioration with age. In an extensive comparative study of turtles and tortoises living in zoos and aquariums, we show that ~75% of 52 species exhibit slow or negligible senescence. For ~80% of species, aging rates are lower than those in modern humans. We find that body weight positively relates to adult life expectancy in both sexes, and sexual size dimorphism explains sex differences in longevity. Unlike humans and other species, we show that turtles and tortoises may reduce senescence in response to improvements in environmental conditions.

Diverse aging rates in ectothermic tetrapods provide insights for the evolution of aging and longevity

Comparative studies of mortality in the wild are necessary to understand the evolution of aging; yet, ectothermic tetrapods are underrepresented in this comparative landscape, despite their suitability for testing evolutionary hypotheses. We present a study of aging rates and longevity across wild tetrapod ectotherms, using data from 107 populations (77 species) of nonavian reptiles and amphibians. We test hypotheses of how thermoregulatory mode, environmental temperature, protective phenotypes, and pace of life history contribute to demographic aging. Controlling for phylogeny and body size, ectotherms display a higher diversity of aging rates compared with endotherms and include phylogenetically widespread evidence of negligible aging. Protective phenotypes and life-history strategies further explain macroevolutionary patterns of aging. Analyzing ectothermic tetrapods in a comparative context enhances our understanding of the evolution of aging.

Over a decade of field physiology reveals life-history specific strategies to drought in garter snakes (Thamnophis legans)

Changing climates and severe weather events can affect population viability. Individuals need to buffer such negative fitness consequences through physiological plasticity. Whether certain life-history strategies are more conducive to surviving changing climates is unknown, but theory predicts that strategies prioritizing maintenance and survival over current reproduction should be better able to withstand such change. We tested this hypothesis in a meta-population of garter snakes having naturally occurring variation in life-history strategies. We tested whether slow pace-of-life (POL) animals, that prioritize survival over reproduction, are more resilient than fast POL animals as measured by several physiological biomarkers. From 2006 to 2019, which included two multi-year droughts, baseline and stress-induced reactivity of plasma corticosterone and glucose varied annually with directionalities consistent with life-history theory. Slow POL animals exhibited higher baseline corticosterone and lower baseline glucose, relative to fast POL animals. These patterns were also observed in stress-induced measures; thus, reactivity was equivalent between ecotypes. However, in drought years, measures of corticosterone did not differ between different life histories. Immune cell distribution showed annual variation independent of drought or life history. These persistent physiological patterns form a backdrop to several extirpations of fast POL populations, suggesting a limited physiological toolkit to surviving periods of extreme drought.

Interpopulation differences in male reproductive effort drive the population dynamics of a host exposed to an emerging fungal pathogen

Compensatory recruitment is a key demographic mechanism that has allowed the coexistence of populations of susceptible amphibians with Batrachochytrium dendrobatidis (Bd), a fungus causing one of the most devastating emerging infectious disease ever recorded among vertebrates. However, the underlying processes (e.g. density-dependent increase in survival at early life stages, change in reproductive traits) as well as the level of interpopulation variation in this response are poorly known. We explore potential mechanisms of compensatory recruitment in response to Bd infection by taking advantage of an amphibian system where male reproductive traits are easy to quantify in free-living populations. The Southern Darwin's frog Rhinoderma darwinii is a vocal sac-brooding species that exhibits a high susceptibility to lethal Bd infection. Using a 7-year capture-recapture study at four populations with contrasting Bd infection status (one high prevalence, one low prevalence and two Bd-free populations), we evaluated whether Bd-positive populations exhibited a higher adult recruitment and a higher male reproductive effort than Bd-negative populations. We also estimated population growth rates to explore whether recruitment compensated for the negative impacts of Bd on the survival of adults. In addition, we evaluated a potential demographic signal of compensatory recruitment (i.e. positive relationship between the proportion of juveniles and Bd prevalence) in response to Bd infection using raw count data from 13 R. darwinii populations. The high Bd prevalence population exhibited the highest male reproductive effort and the highest recruitment among the four monitored populations. This led to a growing population during the study period despite high mortality of adult hosts. In contrast, males from the population with low Bd prevalence had a low reproductive effort and this population, which had the lowest adult recruitment, was declining during the study period despite adults having a higher survival in comparison to the high Bd prevalence population. We also found a demographic signal of compensatory recruitment in response to Bd infection in our broader analysis of 13 R. darwinii populations. Our study underlines the importance of interpopulation variation in life-history strategies on the fate of host populations after infectious disease emergence. Our results also suggest that an increase in reproductive effort can be one of the processes underlying compensatory recruitment in populations of Bd-susceptible amphibians.