Search

Hibernation in brown bears is an annual process involving multiple physiologically distinct seasons—hibernation, active, and hyperphagia. While recent studies have characterized broad patterns of differential gene regulation and isoform usage between hibernation and active seasons, patterns of gene and isoform expression during hyperphagia remain relatively poorly understood. The hyperphagia stage occurs between active and hibernation seasons and involves the accumulation of large fat reserves in preparation for hibernation. Here, we use time-series analyses of gene expression and isoform usage to interrogate transcriptomic regulation associated with all three seasons. We identify a large number of genes with significant differential isoform usage (DIU) across seasons and show that these patterns of isoform usage are largely tissue-specific. We also show that DIU and differential gene-level expression responses are generally non-overlapping, with only a small subset of multi-isoform genes showing evidence of both gene-level expression changes and changes in isoform usage across seasons. Additionally, we investigate nuanced regulation of candidate genes involved in the insulin signaling pathway and find evidence of hyperphagia-specific gene expression and isoform regulation that may enhance fat accumulation during hyperphagia. Our findings highlight the value of using temporal analyses of both gene- and isoform-level gene expression when interrogating complex physiological phenotypes and provide new insight into the mechanisms underlying seasonal changes in bear physiology.

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

A male mutation bias is observed across vertebrates, and, where data are available, this bias is accompanied by increased per-generation mutation rates with parental age. While continuing mitotic cell division in the male germline post puberty has been proposed as the major cellular mechanism underlying both patterns, little direct evidence for this role has been found. Understanding the evolution of the per-generation mutation rate among species requires that we identify the molecular mechanisms that change between species. Here, we study the per-generation mutation rate in an extended pedigree of the brown (grizzly) bear, Ursus arctos horribilis. Brown bears hibernate for one-third of the year, a period during which spermatogenesis slows or stops altogether. The reduction of spermatogenesis is predicted to lessen the male mutation bias and to lower the per-generation mutation rate in this species. However, using whole-genome sequencing, we find that both male bias and per-generation mutation rates are highly similar to that expected for a non-hibernating species. We also carry out a phylogenetic comparison of substitution rates along the lineage leading to brown bear and panda (a non-hibernating species) and find no slowing of the substitution rate in the hibernator. Our results contribute to accumulating evidence that suggests that male germline cell division is not the major determinant of mutation rates and mutation biases. The results also provide a quantitative basis for improved estimates of the timing of carnivore evolution.

The brown bear (Ursus arctos) is the second largest and most widespread extant terrestrial carnivore on Earth and has recently emerged as a medical model for human metabolic diseases. Here, we report a fully-phased chromosome-level assembly of a male North American brown bear built by combining Pacific Biosciences (PacBio) HiFi data and publicly available Hi-C data. The final genome size is 2.47 Gigabases (Gb) with a scaffold and contig N50 length of 70.08 and 43.94 Mb, respectively. BUSCO analysis revealed that 94.5% of single-copy orthologs from mammalia were present in the genome (the highest of any ursid genome to date). Repetitive elements accounted for 44.48% of the genome and a total of 20,480 protein coding genes were identified. Based on whole genome alignment, the brown bear is highly syntenic with the polar bear, and our phylogenetic analysis of 7,246 single-copy BUSCOs supports the currently proposed species tree for Ursidae. This highly contiguous genome assembly will support future research on both the evolutionary history of the bear family and the physiological mechanisms behind hibernation, the latter of which has broad medical implications.SignificanceBrown bears (Ursus arctos) are the most widespread, large terrestrial carnivore on the planet and represent an interesting example of speciation through hybridization, as well as a medical model for sedentary lifestyle-related disease. Although a previous genome for a brown bear has been published, the reported contig N50 was low (only ∼530 kb), despite being scaffolded into putative chromosomes. Genomes of this quality limit the accuracy of analyses which rely on long contiguous stretches of the genome to be assembled (such as with many demographic analyses) as well as attempts at connecting genotype to phenotype (such as in association analyses). In order to support studies on both the complex hybridization history of the brown bear and investigations into medically-relevant phenotypes, we generated a fully-phased, chromosome-level assembly from a male grizzly bear. The genome has a total size of 2.47 Gb and 90% of the genome is contained in 36 scaffolds, roughly corresponding to one autosome per scaffold. This high-quality genome will enable studies across a variety of disciplines, including conservation, evolution, and medicine.

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes. Yet, RE dynamics remain understudied in many taxonomic groups, preventing holistic understanding of how genomes and species evolve. Here, we investigated REs across 601 insect species (20 orders) to better understand the RE landscape of insects and to evaluate automated RE annotation methods in the era of biodiversity genomics. We identified wide variation in the types and frequency of REs across insect groups. We quantified associations between REs and protein-coding genes and found an elevated frequency of associations in insects with abundant long interspersed nuclear elements (LINEs). Sequencing technology impacts RE detection; ∼36% more REs could be identified in long-read versus short-read assemblies. Long terminal repeats (LTRs) showed markedly improved detection in long-read assemblies (162% more), while DNA transposons and LINEs showed less respective technology-related bias. We illustrate fundamental challenges to efficient study of REs in diverse groups, showing that in most insect lineages, 25–85% of repetitive sequences were “unclassified” compared to only ∼13% of unclassified repeats in Drosophila species. Our findings suggest this RE-annotation bottleneck, driven largely by uneven taxonomic representation in RE reference databases, is worsening. Although the diversity of available insect genomes has rapidly expanded, the rate of community contributions to RE databases (essential for RE annotation) has not kept pace, preventing high resolution study of REs in most groups. We highlight the tremendous opportunity and need for the field of biodiversity genomics to embrace REs and suggest collective steps for making progress towards this goal.

Natural history museums hold vast collections of biomaterials. The collections in museums, often painstakingly sampled, are largely unexplored treasures that may help us better understand biodiversity on the planet. Museum collections can provide a unique window into the past of species long gone or currently declining due to human activity. From a molecular perspective, however, many museum samples are stored under conditions that hasten the damage of DNA, RNA and proteins. For example, samples in wet collections are those stored in liquid preservatives, typically ethanol. These ethanol-preserved tissues are often, although not always, formalin-fixed prior to storage, which may damage DNA. In this and recent issues of Molecular Ecology Resources, Straube et al (2021), O'Connell et al (2021) and Hahn et al (2022) explore different types of specimens from museum wet collections as new sources of DNA for scientific studies. All three articles found that for wet museum collections, overall specimen condition mattered most for recovering high-quality genomic DNA., María de Maeztu, Grant/Award Number: Mobility_Fellowship.

The quality of genome assemblies has improved rapidly in recent years due to continual advances in sequencing technology, assembly approaches, and quality control. In the field of molecular ecology, this has led to the development of exceptional quality genome assemblies that will be important long-term resources for broader studies into ecological, conservation, evolutionary, and population genomics of naturally occurring species. Moreover, the extent to which a single reference genome represents the diversity within a species varies: pan-genomes will become increasingly important ecological genomics resources, particularly in systems found to have considerable presence-absence variation in their functional content. Here, we highlight advances in technology that have raised the bar for genome assembly and provide guidance on standards to achieve exceptional quality reference genomes. Key recommendations include the following: (a) Genome assemblies should include long-read sequencing except in rare cases where it is effectively impossible to acquire adequately preserved samples needed for high molecular weight DNA standards. (b) At least one scaffolding approach should be included with genome assembly such as Hi-C or optical mapping. (c) Genome assemblies should be carefully evaluated, this may involve utilising short read data for genome polishing, error correction, k-mer analyses, and estimating the percent of reads that map back to an assembly. Finally, a genome assembly is most valuable if all data and methods are made publicly available and the utility of a genome for further studies is verified through examples. While these recommendations are based on current technology, we anticipate that future advances will push the field further and the molecular ecology community should continue to adopt new approaches that attain the highest quality genome assemblies.

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

A male mutation bias is observed across vertebrates, and, where data are available, this bias is accompanied by increased per-generation mutation rates with parental age. While continuing mitotic cell division in the male germline post-puberty has been proposed as the major cellular mechanism underlying both patterns, little direct evidence for this role has been found. Understanding the evolution of the per-generation mutation rate among species requires that we identify the molecular mechanisms that change between species. Here, we study the per-generation mutation rate in an extended pedigree of the brown (grizzly) bear, Ursus arctos horribilis. Brown bears hibernate for one-third of the year, a period during which spermatogenesis slows or stops altogether. The cessation of spermatogenesis is predicted to lessen the male mutation bias and to lower the per-generation mutation rate in this species. However, using whole-genome sequencing, we find that both male bias and per-generation mutation rates are the same as expected for a non-hibernating species. We also carry out a phylogenetic comparison of substitution rates along the lineage leading to brown bear and panda (a non-hibernating species) and find no slowing of the substitution rate in the hibernator. Our results contribute to accumulating evidence that suggests that male germline cell division is not the major determinant of mutation rates and mutation biases. The results also provide a quantitative basis for improved estimates of the timing of carnivore evolution.

Understanding how metabolic reprogramming happens in cells will aid the progress in the treatment of a variety of metabolic disorders. Brown bears undergo seasonal shifts in insulin sensitivity, including reversible insulin resistance in hibernation. We performed RNA-sequencing on brown bear adipocytes and proteomics on serum to identify changes possibly responsible for reversible insulin resistance. We observed dramatic transcriptional changes, which depended on both the cell and serum season of origin. Despite large changes in adipocyte gene expression, only changes in eight circulating proteins were identified as related to the seasonal shifts in insulin sensitivity, including some that have not previously been associated with glucose homeostasis. The identified serum proteins may be sufficient for shifting hibernation adipocytes to an active-like state.

Changing climate and land‐use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host–plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome‐wide SNP data revealed two major genetic clusters—psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non‐crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co‐occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host–plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non‐crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Convergent evolution, where independent lineages evolve similar traits when adapting to similar habitats, is a common phenomenon and testament to the repeatability of evolutionary processes. Still, non-convergence is also common, and a major question is whether apparently idiosyncratic, lineage-specific evolutionary changes are reflective of chance events inherent to evolutionary processes, or whether they are also influenced by deterministic genetic or ecological factors. To address this question, we quantified the degree of convergence in genome-wide patterns of gene expression across lineages of livebearing fishes (family Poeciliidae) that span 40 million years of evolution and have colonized extreme environments in the form of toxic, hydrogen-sulfide-rich springs. We specifically asked whether the degree of convergence across lineage pairs was related to their phylogenetic relatedness or the ecological similarity of the habitats they inhabit. Using phylogenetic comparative analyses, we showed that the degree of convergence was highly variable across lineage pairs residing in sulfide springs. While closely related lineages did not exhibit higher degrees of convergence than distantly related ones, we uncovered a strong relationship between degree of convergence and ecological similarity. Our results indicate that variation in the degree of convergence is not merely noise associated with evolutionary contingency. Rather, cryptic environmental variation that is frequently ignored when we employ reductionist approaches can significantly contribute to adaptive evolution. This study highlights the importance of multivariate approaches that capture the complexities of both selective regimes and organismal design when assessing the roles of determinism and contingency in evolution.Significance StatementWhen different species adapt to similar environmental conditions, we frequently observe a mix between shared (convergent) and lineage-specific (nonconvergent) evolutionary changes. Shared changes provide evidence for the repeatability and predictability of evolution. However, it remains unclear whether lineage-specific changes are caused by random forces that limit the predictability of evolution, or whether they reflect deterministic processes shaped by unidentified genetic and environmental factors. By analyzing patterns of gene expression across fishes in extreme environments, we show that the degree of convergence between lineages is related to ecology, indicating that lineage-specific evolutionary changes are not just noise caused by random processes. Thus, acknowledging the complexity of nature in empirical research is critical if we want to predict evolution.

In less than 25 y, the field of animal genome science has transformed from a discipline seeking its first glimpses into genome sequences across the Tree of Life to a global enterprise with ambitions to sequence genomes for all of Earth’s eukaryotic diversity [H. A. Lewin et al. , Proc. Natl. Acad. Sci. U.S.A. 115, 4325–4333 (2018)]. As the field rapidly moves forward, it is important to take stock of the progress that has been made to best inform the discipline’s future. In this Perspective, we provide a contemporary, quantitative overview of animal genome sequencing. We identified the best available genome assemblies in GenBank, the world’s most extensive genetic database, for 3,278 unique animal species across 24 phyla. We assessed taxonomic representation, assembly quality, and annotation status for major clades. We show that while tremendous taxonomic progress has occurred, stark disparities in genomic representation exist, highlighted by a systemic overrepresentation of vertebrates and underrepresentation of arthropods. In terms of assembly quality, long-read sequencing has dramatically improved contiguity, whereas gene annotations are available for just 34.3% of taxa. Furthermore, we show that animal genome science has diversified in recent years with an ever-expanding pool of researchers participating. However, the field still appears to be dominated by institutions in the Global North, which have been listed as the submitting institution for 77% of all assemblies. We conclude by offering recommendations for improving genomic resource availability and research value while also broadening global representation.

SummaryUnderstanding hibernation in brown bears (Ursus arctos) can provide insight into many human diseases. During hibernation, brown bears experience states of insulin resistance, physical inactivity, extreme bradycardia, obesity, and the absence of urine production. These states closely mimic human diseases such as type 2 diabetes, muscle atrophy, renal and heart failure, cachexia, and obesity. The reversibility of these states from hibernation to active season allows for the identification of novel mediators with possible therapeutic value for humans. Recent studies have identified genes and pathways that are differentially expressed between active and hibernation seasons. However, little is known about the role of differential expression of gene isoforms on hibernation physiology. To identify both distinct and novel mRNA isoforms, we performed full-length RNA-sequencing (Iso-Seq) on three tissue types from three individuals sampled during both active and hibernation seasons. We combined the long-read data with the reference annotation for an improved transcriptome and mapped RNA-seq data from six individuals to the improved transcriptome to quantify differential isoform usage between tissues and seasons. We identified differentially expressed isoforms in all study tissues and showed that adipose has a high level of differential isoform usage with isoform switching, regardless of whether the genes were differentially expressed. Our analyses provide a comprehensive evaluation of isoform usage between active and hibernation states, revealing that differential isoform usage, even in the absence of differential gene expression, is an important mechanism for modulating genes during hibernation. These findings demonstrate the value of isoform expression studies and will serve as the basis for deeper exploration into hibernation biology.

Long-read sequencing is driving a new reality for genome science where highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we leveraged long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Opthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance, and vision have evidence of selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: hemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the hemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

Divergent selection along ecological gradients can lead to speciation, and replicated speciation events occur when populations of multiple lineages undergo divergence following colonization of similar environments. In such instances, it remains unclear to what extent reproductive isolation evolves via convergent mechanisms at the genomic level due to biases associated with shared ancestry and variation in selection associated with geographic replicates typically studied in this context. We used a unique system in which three species of poeciliid fishes occur in sympatry in a toxic, hydrogen sulfide (H2S)-rich spring and an adjacent nonsulfidic stream to examine shared patterns of adaptive divergence across multiple levels of biological organization. Despite small spatial scales, we found strong genetic differentiation between populations in sulfidic and nonsulfidic habitats mediated by strong selection against migrants between habitat types. High levels of reproductive isolation were accompanied by convergent patterns of adaptation in morphological and physiological traits, as well as genome-wide patterns of gene expression across all three species. Furthermore, the mitochondrial genomes of each species exhibit shared signatures of selection on key genes involved in H2S toxicity. However, contrary to predictions of speciation theory, analyses of divergence across the nuclear genomes neither revealed evidence for clear genomic islands of speciation nor substantial congruence of outlier regions across population pairs. Instead, heterogenous regions of divergence spread across the genome suggest that selection for polygenic physiological adaptations likely facilitated the rapid evolution of high levels of reproductive isolation. Overall, we demonstrate that substantial convergence across multiple levels of biological organization can be mediated by non-convergent modifications at the genomic level. By disentangling variation in natural selection from lineage-specific evolution in this system of highly divergent, yet sympatric lineages, our results emphasize the outsized role of the genomic substrate upon which selection acts in driving convergent evolution at the phenotypic level.Significance StatementDivergent lineages that coexist in sympatry and are exposed to the same sources of natural selection provide a unique opportunity to study convergent evolution across levels of organization because confounding factors associated with geographic replications are eliminated. Using three sympatric lineages of livebearing fishes inhabiting toxic and adjacent nontoxic habitats, we show that the convergent evolution of phenotypic adaptation and reproductive isolation can evolve in the absence of substantial convergence at the genomic level. Our results suggest that the nature of selection is less important in producing genomic convergence than the initial genomic substrates that are available for selection to act on.

Significance Environmental factors can promote phenotypic variation through alterations in the epigenome and mediate adaptation of an organism to the environment. Observations suggest the adaptation of Poecilia mexicana fish to toxic, hydrogen sulfide–rich environments in southern Mexico may, in part, be promoted through epigenetic DNA methylation alterations that became generationally stable and are inherited to subsequent generations independent of the environment. Environmental epigenetics may provide an important mechanism mediating adaptation in this species. This is an observation that the epigenome is stably inherited generationally through the germline after the removal of an environmental stressor (i.e., hydrogen sulfide) from a wild population., Environmental factors can promote phenotypic variation through alterations in the epigenome and facilitate adaptation of an organism to the environment. Although hydrogen sulfide is toxic to most organisms, the fish Poecilia mexicana has adapted to survive in environments with high levels that exceed toxicity thresholds by orders of magnitude. Epigenetic changes in response to this environmental stressor were examined by assessing DNA methylation alterations in red blood cells, which are nucleated in fish. Males and females were sampled from sulfidic and nonsulfidic natural environments; individuals were also propagated for two generations in a nonsulfidic laboratory environment. We compared epimutations between the sexes as well as field and laboratory populations. For both the wild-caught (F0) and the laboratory-reared (F2) fish, comparing the sulfidic and nonsulfidic populations revealed evidence for significant differential DNA methylation regions (DMRs). More importantly, there was over 80% overlap in DMRs across generations, suggesting that the DMRs have stable generational inheritance in the absence of the sulfidic environment. This is an example of epigenetic generational stability after the removal of an environmental stressor. The DMR-associated genes were related to sulfur toxicity and metabolic processes. These findings suggest that adaptation of P. mexicana to sulfidic environments in southern Mexico may, in part, be promoted through epigenetic DNA methylation alterations that become stable and are inherited by subsequent generations independent of the environment.

The first insect genome assembly (Drosophila melanogaster) was published two decades ago. Today, nuclear genome assemblies are available for a staggering 601 insect species representing 20 orders. In this study, we analyzed the most-contiguous assembly for each species and provide a “state-of-the-field” perspective, emphasizing taxonomic representation, assembly quality, gene completeness, and sequencing technologies. Relative to species richness, genomic efforts have been biased toward four orders (Diptera, Hymenoptera, Collembola, and Phasmatodea), Coleoptera are underrepresented, and 11 orders still lack a publicly available genome assembly. The average insect genome assembly is 439.2 Mb in length with 87.5% of single-copy benchmarking genes intact. Most notable has been the impact of long-read sequencing; assemblies that incorporate long reads are ∼48× more contiguous than those that do not. We offer four recommendations as we collectively continue building insect genome resources: 1) seek better integration between independent research groups and consortia, 2) balance future sampling between filling taxonomic gaps and generating data for targeted questions, 3) take advantage of long-read sequencing technologies, and 4) expand and improve gene annotations.

With more than 3,700 described species, stoneflies (Order Plecoptera) are an important component of global aquatic biodiversity. The meltwater stonefly Lednia tumana (Ricker, 1952) is endemic to al...

Hibernation is characterized by depression of many physiological processes. To determine if this state is reversible in a non-food caching species, we fed hibernating grizzly bears (Ursus arctos horribilis) dextrose for 10 days to replace 53% or 100% of the estimated minimum daily energetic cost of hibernation. Feeding caused serum concentrations of glycerol and ketones (β-hydroxybutyrate) to return to active season levels irrespective of the amount of glucose fed. By contrast, free fatty acids (FFAs) and indices of metabolic rate, such as general activity, heart rate, strength of heart rate circadian rhythm, and insulin sensitivity were restored to approximately 50% of active season levels. Body temperature was unaffected by feeding. To determine the contribution of adipose to the metabolic effects observed after glucose feeding, we cultured bear adipocytes collected at the beginning and end of the feeding and performed metabolic flux analysis. We found a ∼33% increase in energy metabolism after feeding. Moreover, basal metabolism before feeding was 40% lower in hibernation cells compared with fed cells or active cells cultured at 37°C, thereby confirming the temperature independence of metabolic rate. The partial depression of circulating FFAs with feeding likely explains the incomplete restoration of insulin sensitivity and other metabolic parameters in hibernating bears. Further depression of metabolic function is likely to be an active process. Together, the results provide a highly controlled model to examine the relationship between nutrient availability and metabolism on the hibernation phenotype in bears.

Hibernation is characterized by suppression of many physiological processes. To determine if this state is reversible in a non-food caching species, we fed hibernating grizzly bears (Ursus arctos horribilis) glucose for 10 days to replace 53% or 100% of the estimated minimum daily energetic cost of hibernation. Feeding caused serum concentrations of glycerol and ketones (ß-hydroxybutyrate) to return to active season levels irrespective of the amount of glucose fed. By contrast, free-fatty acids and indices of metabolic rate, such as general activity, heart rate, and strength of the daily heart rate rhythm and insulin sensitivity were restored to roughly 50% of active season levels. Body temperature was unaffected by feeding. To determine the contribution of adipose to these metabolic effects of glucose feeding we cultured bear adipocytes collected at the beginning and end of the feeding and performed metabolic flux analysis. We found a roughly 33% increase in energy metabolism after feeding. Moreover, basal metabolism before feeding was 40% lower in hibernation cells compared to fed cells or active cells cultured at 37°C, thereby confirming the temperature independence of metabolic rate. The partial suppression of circulating FFA with feeding likely explains the incomplete restoration of insulin sensitivity and other metabolic parameters in hibernating bears. Further suppression of metabolic function is likely an active process. Together, the results provide a highly controlled model to examine the relationship between nutrient availability and metabolism on the hibernation phenotype in bears.

Modern genetic data sets present unprecedented opportunities to understand the evolutionary origins of taxonomic groups comprising hundreds to thousands of species. When the timing of key events are known, it is also possible to investigate biogeographic history in the context of major phenomena (e.g., continental drift). In this study, we investigated the biogeographic history of the suborder Zoarcoidei, a globally distributed fish group that includes species inhabiting both poles and multiple taxa that produce antifreeze proteins to survive chronic subfreezing temperatures. We first generated a multi-locus, time-calibrated phylogeny for the group. We then used biogeographic modeling to reconstruct ancestral ranges across the tree and quantify the type and frequency of biogeographic events (e.g., founder, dispersal). With these results, we considered how the cooling of the Southern and Arctic Oceans, which reached their present-day subfreezing temperatures 10-15 million years ago (Mya) and 2-3 Mya, respectively, may have shaped the evolutionary history of Zoarcoidei, with an emphasis on the most speciose and widely distributed family, eelpouts (family Zoarcidae). Our phylogenetic results clarified standing issues in the Zoarcoidei taxonomy and showed that the group began to diversify in the Oligocene ∼31-32 Mya, with the center of origin for all families in north temperate waters. Within-area speciation was the most common biogeographic event in the group’s history (80% of all events) followed by dispersal (20%). Finally, we found mixed evidence for polar ocean cooling underpinning Zoarcoidei diversification, with support limited to eelpout speciation in the Southern Ocean over the last 10 million years.

The global cryosphere, Earth's frozen water, is in precipitous decline. The ongoing and predicted impacts of cryosphere loss are diverse, ranging from disappearance of entire biomes to crises of water availability. Covering approximately one-fifth of the planet, mass loss from the terrestrial cryosphere is driven primarily by a warming atmosphere but reductions in albedo (the proportion of reflected light) also contribute by increasing absorption of solar radiation. In addition to dust and other abiotic impurities, biological communities substantially reduce albedo worldwide. In this review, we provide a global synthesis of biological albedo reduction (BAR) in terrestrial snow and ice ecosystems. We first focus on known drivers—algal blooms and cryoconite (granular sediment on the ice that includes both mineral and biological material)—as they account for much of the biological albedo variability in snow and ice habitats. We then consider an array of potential drivers of BAR whose impacts may be overlooked, such as arthropod deposition, resident organisms (e.g., dark-bodied glacier ice worms), and larger vertebrates, including humans, that transiently visit the cryosphere. We consider both primary (e.g., BAR due to the presence of pigmented algal cells) and indirect (e.g., nutrient addition from arthropod deposition) effects, as well as interactions among biological groups (e.g., birds feeding on ice worms). Collectively, we highlight that in many cases, overlooked drivers and interactions among factors have considerable potential to alter BAR, perhaps rivaling the direct effects of algal blooms and cryoconite. We conclude by highlighting knowledge gaps for the field with an emphasis on the underrepresentation of genomic tools, understudied areas (particularly high-elevation glaciers at tropical latitudes), and a dearth of temporal sampling in current efforts. We detail a global framework for long-term BAR monitoring that, if implemented, would yield a tremendous amount of insight for BAR and would be particularly valuable in light of the rapid ecological and physical changes occurring in the contemporary cryosphere.

Dam construction and longitudinal river habitat fragmentation disrupt important life histories and movement of aquatic species. This is especially true for Oncorhynchus mykiss that exhibits both migratory (steelhead) and non-migratory (resident rainbow) forms. While the negative effects of dams on salmonids have been extensively documented, few studies have had the opportunity to compare population genetic diversity and structure prior to and following dam removal. Here we examine the impacts of the removal of two dams on the Elwha River on the population genetics of O. mykiss. Genetic data were produced from &gt, 1200 samples collected prior to dam removal from both life history forms, and post-dam removal from steelhead. We identified three genetic clusters prior to dam removal primarily explained by isolation due to dams and natural barriers. Following dam removal, genetic structure decreased and admixture increased. Despite large O. mykiss population declines after dam construction, we did not detect shifts in population genetic diversity or allele frequencies of loci putatively involved in migratory phenotypic variation. Steelhead descendants from formerly below and above dammed populations recolonized the river rapidly after dam removal, suggesting that dam construction did not significantly reduce genetic diversity underlying O. mykiss life history strategies. These results have significant evolutionary implications for the conservation of migratory adaptive potential in O. mykiss populations above current anthropogenic barriers.

The most remarkable feature of our planet is the diversity of its life forms, ranging from viruses and nanobacteria to blue whales and giant sequoias to satanic leaf‐tailed geckos and leafy seadragons (look them up!). Life is found in essentially all environments on earth, and the number of species living on our planet is many times greater than we could have imagined a century ago. A well‐regarded estimate pegs the number of eukaryotic species on earth at 8.7 million (±1.3 million), of which fewer than 15% are currently described (Mora et al., 2011). The diversity of prokaryotes is less clear (and highly controversial), but an analysis of 1.6 billion 16S ribosomal RNA sequences estimated that 0.8–1.6 million prokaryotic operational taxonomic units exist globally (Louca et al., 2019). While we do not know how many species are currently extant, or have existed in the past, we do know that this biodiversity is valuable, providing food, fibre and medicine, furnishing ecosystem services such as water and air purification, nutrient cycling, pollination and carbon uptake, and contributing to technological innovations ranging from biotechnology to robotics to material science. Moreover, biodiversity underlies the cultural identity of human populations and is important to human health and well‐being. Geographically, species richness increases from the Polar Regions to the tropics in terrestrial and surface marine ecosystems. Thus, some countries, especially those in tropical and subtropical regions, are endowed with much greater biodiversity than others. Unfortunately, benefits arising from the access and utilization of this biodiversity have been unequally shared, with (paradoxically) biodiversity‐poor countries often accruing the lion's share of economic gains. There can be imbalances within countries as well, wherein some segments of the population obtain greater economic benefits from biodiversity and associated traditional knowledge than indigenous peoples. The “Nagoya Protocol on Access to Genetic Resources and the Fair and Equitable Sharing of Benefits Arising from their Utilization,” which came into force in 2014, is an international agreement designed to ensure that the benefits arising from biodiversity are shared equitably (https://www.cbd.int/abs/). However, few scientific journals require compliance with the Nagoya Protocol or the reporting of benefits from biodiversity research. In this editorial, we (the editors of Molecular Ecology and Molecular Ecology Resources) express support for the Nagoya Protocol and the principle of benefit sharing. We believe that scientific journals publishing research on biodiversity can play an important role in implementing the Nagoya Protocol and in reporting on benefits generated from such research. Below, we provide background on the Nagoya Protocol, discuss the kinds of benefits that may arise from biodiversity research, describe the rationale for reporting on these benefits and introduce changes to the journals’ Data Accessibility Statements to incorporate the requirements and goals of the Nagoya Protocol.

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H2 S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H2 S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H2 S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H2 S adaptation by promoting functions needed for survival and reducing functions affected by H2 S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H2 S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

Simple Summary The genome is the basic evolutionary unit underpinning life on Earth. Knowing its sequence, including the many thousands of genes coding for proteins in an organism, empowers scientific discovery for both the focal organism and related species. Aquatic insects represent 10% of all insect diversity, can be found on every continent except Antarctica, and are key components of freshwater ecosystems. However, aquatic insect genome biology lags dramatically behind that of terrestrial insects. If genomic effort was spread evenly, one aquatic insect genome would be sequenced for every ~9 terrestrial insect genomes. Instead, ~24 terrestrial insect genomes have been sequenced for every aquatic insect genome. A lack of aquatic genomes is limiting research progress in the field at both fundamental and applied scales. We argue that the limited availability of aquatic insect genomes is not due to practical limitations—small body sizes or overly complex genomes—but instead reflects a lack of research interest. We call for targeted efforts to expand the availability of aquatic insect genomic resources to empower future research. Abstract Aquatic insects comprise 10% of all insect diversity, can be found on every continent except Antarctica, and are key components of freshwater ecosystems. However, aquatic insect genome biology lags dramatically behind that of terrestrial insects. If genomic effort was spread evenly, one aquatic insect genome would be sequenced for every ~9 terrestrial insect genomes. Instead, ~24 terrestrial insect genomes have been sequenced for every aquatic insect genome. This discrepancy is even more dramatic if the quality of genomic resources is considered; for instance, while no aquatic insect genome has been assembled to the chromosome level, 29 terrestrial insect genomes spanning four orders have. We argue that a lack of aquatic insect genomes is not due to any underlying difficulty (e.g., small body sizes or unusually large genomes), yet it is severely hampering aquatic insect research at both fundamental and applied scales. By expanding the availability of aquatic insect genomes, we will gain key insight into insect diversification and empower future research for a globally important taxonomic group.

How aquatic insects cope with cold temperatures is poorly understood. This is particularly true for high-elevation species that often experience a seasonal risk of freezing. In the Rocky Mountains, nemourid stoneflies (Plecoptera: Nemouridae) are a major component of mountain stream biodiversity and are typically found in streams fed by glaciers and snowfields, which due to climate change, are rapidly receding. Predicting the effects of climate change on mountain stoneflies is difficult because their thermal physiology is largely unknown. We investigated cold tolerance of several alpine stoneflies (Lednia tumana, Lednia tetonica, and Zapada spp.) from the Rocky Mountains, USA. We measured the supercooling point (SCP) and tolerance to ice enclosure of late-instar nymphs collected from a range of thermal regimes. SCPs varied among species and populations, with the lowest SCP measured for nymphs from an alpine pond, which is much more likely to freeze solid in winter than flowing streams. We also show that L. tumana cannot survive being enclosed in ice, even for short periods of time (less than three hours) at relatively mild temperatures (−0.5 °C). Our results indicate that high-elevation stoneflies at greater risk of freezing may have correspondingly lower SCPs, and despite their common association with glacial meltwater, they appear to be living near their lower thermal limits.

Revealing the mechanisms underlying the reversible physiology of hibernation could have applications to both human and animal health as hibernation is often associated with disease-like states. The present study uses RNA-sequencing to reveal the tissue and seasonal transcriptional changes occurring in grizzly bears (

Changing climate and land‐use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host–plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome‐wide SNP data revealed two major genetic clusters—psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non‐crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co‐occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host–plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non‐crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Rapid glacier recession is altering the physical conditions of headwater streams. Stream temperatures are predicted to rise and become increasingly variable, putting entire meltwater-associated biological communities at risk of extinction. Thus, there is a pressing need to understand how thermal stress affects mountain stream insects, particularly where glaciers are likely to vanish on contemporary timescales. In this study, we tested the critical thermal maximum (CTMAX) of stonefly nymphs representing multiple species and a range of thermal regimes in the high Rocky Mountains, USA. We then collected RNA-sequencing data to assess how organismal thermal stress translated to the cellular level. Our focal species included the meltwater stonefly, Lednia tumana, which was recently listed under the U.S. Endangered Species Act due to climate-induced habitat loss. For all study species, critical thermal maxima (CTMAX > 20°C) far exceeded the stream temperatures mountain stoneflies experience (< 10°C). Moreover, while evidence for a cellular stress response was present, we also observed constitutive expression of genes encoding proteins known to underlie thermal stress (i.e., heat shock proteins) even at low temperatures that reflected natural conditions. We show that high-elevation aquatic insects may not be physiologically threatened by short-term exposure to warm temperatures and that longer term physiological responses or biotic factors (e.g., competition) may better explain their extreme distributions.

Eye regression occurs across cave-dwelling populations of many species and is often coupled with a decrease or loss in eye function. Teleost fishes are among the few vertebrates to undergo widespread colonization of caves and often exhibit eye regression with blindness. Cave populations of the poeciliid fish Poecilia mexicana (cave molly) exhibit reduced—albeit functional—eyes, offering the opportunity to investigate partial eye regression. We sequenced eye transcriptomes of cave and surface populations of P. mexicana to identify differentially expressed genes that potentially underlie eye regression in cave mollies. We identified 28 significantly differentially expressed genes, 20 of which were directly related to light sensitivity, eye structure and visual signaling. Twenty-six of these genes were downregulated in cave compared to surface populations. Functional enrichment analysis revealed eye-related gene ontologies that were under-represented in cave mollies. In addition, a set of co-expressed genes related to vision and circadian rhythm was correlated with habitat type (cave versus surface). Our study suggests that differential gene expression plays a key role in the beginning evolutionary stages of eye regression in P. mexicana , shedding further light on regressive evolution in cavefish.

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

The mangrove rivulus, Kryptolebias marmoratus, is one of only two self-fertilizing hermaphroditic fish species and inhabits mangrove forests. While selfing can be advantageous, it reduces heterozygosity and decreases genetic diversity. Studies using microsatellites found that there are variable levels of selfing among populations of K. marmoratus, but overall, there is a low rate of outcrossing and, therefore, low heterozygosity. In this study, we used whole-genome data to assess the levels of heterozygosity in different lineages of the mangrove rivulus and infer the phylogenetic relationships among those lineages. We sequenced whole genomes from 15 lineages that were completely homozygous at microsatellite loci and used single nucleotide polymorphisms (SNPs) to determine heterozygosity levels. More variation was uncovered than in studies using microsatellite data because of the resolution of full genome sequencing data. Moreover, missense polymorphisms were found most often in genes associated with immune function and reproduction. Inferred phylogenetic relationships suggest that lineages largely group by their geographic distribution. The use of whole-genome data provided further insight into genetic diversity in this unique species. Although this study was limited by the number of lineages that were available, these data suggest that there is previously undescribed variation within lineages of K. marmoratus that could have functional consequences and (or) inform us about the limits to selfing (e.g., genetic load, accumulation of deleterious mutations) and selection that might favor the maintenance of heterozygosity. These results highlight the need to sequence additional individuals within and among lineages.

Modern genetic data sets present unprecedented opportunities to understand the evolutionary origins of diverse taxonomic groups. When the timing of key events is known, it is possible to investigate biogeographic history in the context of major phenomena (e.g., cooling of a major ocean). In this study, we investigated the biogeographic history of the suborder Zoarcoidei, a globally distributed fish group that includes species inhabiting both poles that produce antifreeze proteins to survive chronic subfreezing temperatures. We first generated a multi-locus, time-calibrated phylogeny for the group. We then used biogeographic modeling to reconstruct ancestral ranges across the tree and to quantify the type and frequency of biogeographic events (e.g., founder, dispersal). With these results, we considered how the cooling of the Southern and Arctic Oceans, which reached their present-day subfreezing temperatures 10–15 million years ago (Mya) and 2–3 Mya, respectively, may have shaped the group’s evolutionary history, with an emphasis on the most speciose and widely distributed family, eelpouts (family Zoarcidae). Our phylogenetic results clarified the Zoarcoidei taxonomy and showed that the group began to diversify in the Oligocene ~31–32 Mya, with the center of origin for all families in north temperate waters. Within-area speciation was the most common biogeographic event in the group’s history (80% of all events) followed by dispersal (20%). Finally, we only found evidence, albeit limited, for ocean cooling underpinning diversification of eelpouts living in the high Antarctic over the last 10 million years.

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

Comparative studies struggle with the challenge of balancing technical properties with the requirement of obtaining samples from multiple species. Guidelines are needed to minimize the effect of confounders and increase the robustness of the inference. We argue for the paramount importance of extensive record keeping and reporting. Comparative studies struggle to balance technical properties with the need to obtain samples from multiple species. The authors argue for extensive record keeping and reporting of metadata to minimize the effect of confounders and increase the robustness of inferences from these studies.

Hydrogen sulfide (H2S) is a respiratory toxicant that creates extreme environments tolerated by few organisms. H2S is also produced endogenously by metazoans and plays a role in cell signaling. The mechanisms of H2S toxicity and its physiological functions serve as a basis to discuss the multifarious strategies that allow animals to survive in H2S-rich environments. Despite their toxicity, H2S-rich environments also provide ecological opportunities, and complex selective regimes of covarying abiotic and biotic factors drive trait evolution in organisms inhabiting H2S-rich environments. Furthermore, adaptation to H2S-rich environments can drive speciation, giving rise to biodiversity hot spots with high levels of endemism in deep-sea hydrothermal vents, cold seeps, and freshwater sulfide springs. The diversity of H2S-rich environments and their inhabitants provides ideal systems for comparative studies of the effects of a clear-cut source of selection across vast geographic and phylogenetic scales, ultimately informing our understanding of how environmental stressors affect ecological and evolutionary processes.