Polar bears exhibit genome-wide signatures of bioenergetic adaptation to life in the arctic environment

Where the microbes aren't

Although a large fraction of Earth's volume and most places beyond the planet lack life because physical and chemical conditions are too extreme, intriguing scientific questions are raised in many environments within or at the edges of life's niche space in which active life is absent. This review explores the environments in which active microorganisms do not occur. Within the known niche space for life, uninhabited, but habitable physical spaces potentially offer opportunities for hypothesis testing, such as using them as negative control environments to investigate the influence of life on planetary processes. At the physico-chemical limits of life, questions such as whether spaces devoid of actively metabolizing or reproducing life constitute uninhabitable space or space containing vacant niches that could be occupied with appropriate adaptation are raised. We do not know the extent to which evolution has allowed life to occupy all niche space within its biochemical potential. The case of habitable extraterrestrial environments and the scientific and ethical questions that they raise is discussed.

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Insights into bear evolution from a Pleistocene polar bear genome

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear’s lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

Long-read isoform sequencing reveals tissue-specific isoform expression between active and hibernating brown bears (Ursus arctos)

Understanding hibernation in brown bears (Ursus arctos) can provide insight into some human diseases. During hibernation, brown bears experience periods of insulin resistance, physical inactivity, extreme bradycardia, obesity, and the absence of urine production. These states closely mimic aspects of human diseases such as type 2 diabetes, muscle atrophy, as well as renal and heart failure. The reversibility of these states from hibernation to active season enables the identification of mediators with possible therapeutic value for humans. Recent studies have identified genes and pathways that are differentially expressed between active and hibernation seasons in bears. However, little is known about the role of differential expression of gene isoforms on hibernation physiology. To identify both distinct and novel mRNA isoforms, full-length RNA-sequencing (Iso-Seq) was performed on adipose, skeletal muscle, and liver from three individual bears sampled during both active and hibernation seasons. The existing reference genome annotation was improved by combining it with the Iso-Seq data. Short-read RNA-sequencing data from six individuals were mapped to the new reference annotation to quantify differential isoform usage (DIU) between tissues and seasons. We identified differentially expressed isoforms in all three tissues, to varying degrees. Adipose had a high level of DIU with isoform switching, regardless of whether the genes were differentially expressed. Our analyses revealed that DIU, 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.

The Antarctic Weddell seal genome reveals evidence of selection on cardiovascular phenotype and lipid handling

The Weddell seal (Leptonychotes weddellii) thrives in its extreme Antarctic environment. We generated the Weddell seal genome assembly and a high-quality annotation to investigate genome-wide evolutionary pressures that underlie its phenotype and to study genes implicated in hypoxia tolerance and a lipid-based metabolism. Genome-wide analyses included gene family expansion/contraction, positive selection, and diverged sequence (acceleration) compared to other placental mammals, identifying selection in coding and non-coding sequence in five pathways that may shape cardiovascular phenotype. Lipid metabolism as well as hypoxia genes contained more accelerated regions in the Weddell seal compared to genomic background. Top-significant genes were SUMO2 and EP300; both regulate hypoxia inducible factor signaling. Liver expression of four genes with the strongest acceleration signals differ between Weddell seals and a terrestrial mammal, sheep. We also report a high-density lipoprotein-like particle in Weddell seal serum not present in other mammals, including the shallow-diving harbor seal.

Fatty acid profiles of feeding and fasting bears: estimating calibration coefficients, the timeframe of diet estimates, and selective mobilization during hibernation

Accurate information on diet composition is central to understanding and conserving carnivore populations. Quantitative fatty acid signature analysis (QFASA) has emerged as a powerful tool for estimating the diets of predators, but ambiguities remain about the timeframe of QFASA estimates and the need to account for species-specific patterns of metabolism. We conducted a series of feeding experiments with four juvenile male brown bears (Ursus arctos) to (1) track the timing of changes in adipose tissue composition and QFASA diet estimates in response to a change in diet and (2) quantify the relationship between consumer and diet FA composition (i.e., determine "calibration coefficients"). Bears were fed three compositionally distinct diets for 90-120 days each. Two marine-based diets were intended to approximate the lipid content and composition of the wild diet of polar bears (U. maritimus). Bear adipose tissue composition changed quickly in the direction of the diet and showed evidence of stabilization after 60 days. During hibernation, FA profiles were initially stable but diet estimates after 10 weeks were sensitive to calibration coefficients. Calibration coefficients derived from the marine-based diets were broadly similar to each other and to published values from marine-fed mink (Mustela vison), which have been used as a model for free-ranging polar bears. For growing bears on a high-fat diet, the temporal window for QFASA estimates was 30-90 days. Although our results reinforce the importance of accurate calibration, the similarities across taxa and diets suggest it may be feasible to develop a generalized QFASA approach for mammalian carnivores.

Whole genome sequencing reveals a complex introgression history and the basis of adaptation to subarctic climate in wild sheep

To predict species responses to anthropogenic disturbances and climate change, it is reasonable to use species with high sensitivity to such factors. Snow sheep (Ovis nivicola) could represent a good candidate for this; as the only large herbivore species adapted to the cold and alpine habitats of northeastern Siberia, it plays a crucial role in its ecosystem. Despite having an extensive geographical distribution among all ovine species, it is one of the least studied. In this study, we sequenced and analysed six genomes of snow sheep in combination with all other wild sheep species to infer key aspects of their evolutionary history and unveil the genetic basis of their adaptation to subarctic environments. Despite their large census population size, snow sheep genomes showed remarkably low heterozygosity, which could reflect the effect of isolation and historical bottlenecks that we inferred using the pairwise sequential Markovian coalescent and runs of homozygosity. F₄ -statistics indicated instances of introgression involving snow sheep with argali (Ovis ammon) and Dall (Ovis dalli) sheep, suggesting that these species might have been more widespread during the Pleistocene. Furthermore, the introgressed segments, which were identified using mainly minimum relative node depth, covered genes associated with immunity, adipogenesis and morphology-related traits, representing potential targets of adaptive introgression. Genes related to mitochondrial functions and thermogenesis associated with adipose tissue were identified to be under selection. Overall, our data suggest introgression as a mechanism facilitating adaptation in wild sheep species and provide insights into the genetic mechanisms underlying cold adaptation in snow sheep.

Patterns of Natural Selection on Mitochondrial Protein-Coding Genes in Lungless Salamanders: Relaxed Purifying Selection and Presence of Positively Selected Codon Sites in the Family Plethodontidae

There are two distinct lungless groups in caudate amphibians (salamanders and newts) (the family Plethodontidae and the genus Onychodactylus, from the family Hynobiidae). Lunglessness is considered to have evolved in response to environmental and/or ecological adaptation with respect to oxygen requirements. We performed selection analyses on lungless salamanders to elucidate the selective patterns of mitochondrial protein-coding genes associated with lunglessness. The branch model and RELAX analyses revealed the occurrence of relaxed selection (an increase of the dN/dS ratio = ω value) in most mitochondrial protein-coding genes of plethodontid salamander branches but not in those of Onychodactylus. Additional branch model and RELAX analyses indicated that direct-developing plethodontids showed the relaxed pattern for most mitochondrial genes, although metamorphosing plethodontids had fewer relaxed genes. Furthermore, aBSREL analysis detected positively selected codons in three plethodontid branches but not in Onychodactylus. One of these three branches corresponded to the most recent common ancestor, and the others corresponded with the most recent common ancestors of direct-developing branches within Hemidactyliinae. The positive selection of mitochondrial protein-coding genes in Plethodontidae is probably associated with the evolution of direct development.

Pervasive Genomic Signatures of Local Adaptation to Altitude Across Highland Specialist Andean Hummingbird Populations

Populations along steep environmental gradients are subject to differentiating selection that can result in local adaptation, despite countervailing gene flow, and genetic drift. In montane systems, where species are often restricted to narrow ranges of elevation, it is unclear whether the selection is strong enough to influence functional differentiation of subpopulations differing by a few hundred meters in elevation. We used targeted capture of 12 501 exons from across the genome, including 271 genes previously implicated in altitude adaptation, to test for adaptation to local elevations for 2 highland hummingbird species, Coeligena violifer (n = 62) and Colibri coruscans (n = 101). For each species, we described population genetic structure across the complex geography of the Peruvian Andes and, while accounting for this structure, we tested whether elevational allele frequency clines in single nucleotide polymorphisms (SNPs) showed evidence for local adaptation to elevation. Although the 2 species exhibited contrasting population genetic structures, we found signatures of clinal genetic variation with shifts in elevation in both. The genes with SNP-elevation associations included candidate genes previously discovered for high-elevation adaptation as well as others not previously identified, with cellular functions related to hypoxia response, energy metabolism, and immune function, among others. Despite the homogenizing effects of gene flow and genetic drift, natural selection on parts of the genome evidently optimizes elevation-specific cellular function even within elevation range-restricted montane populations. Consequently, our results suggest local adaptation occurring in narrow elevation bands in tropical mountains, such as the Andes, may effectively make them "taller" biogeographic barriers.

Diving deep: understanding the genetic components of hypoxia tolerance in marine mammals

Marine mammals have highly specialized physiology, exhibited in many species by extreme breath-holding capabilities that allow deep dives and extended submergence. Cardiovascular control and cell-level hypoxia tolerance are key features of this phenotype. Identifying genomic signatures tied to physiology will be valuable in understanding these natural model species, which may generate translational opportunities to human diseases arising from hypoxic stress or tissue injury. Genomic analyses have now been conducted in dolphins, river dolphins, minke whales, bowhead whales, and polar bears, with multispecies studies exploring evolutionary signals across marine mammal lineages, encompassing extinct and extant divers. Single-species genome studies for sirenians do not yet exist. Extant marine mammals arose in three lineages from separate aquatic recolonizations. Their physiological specializations, along with these independent origins create an interesting case to examine convergent evolution. Although molecular mechanisms of hypoxia tolerance are not universally apparent across marine mammal genomic studies, altered evolutionary rates have been identified for genes linked to oxygen binding and transport (e.g., MB, HBA, and HBB), blood pressure control (e.g., endothelin pathway genes), and cell protection in multiple species. Despite convergent phenotypes across clades, instances of identical molecular convergence have been uncommon. Given the inherent logistical and regulatory difficulties associated with functional genetic experiments in marine mammals, several avenues of further investigation are suggested to enable validation of candidate genes for hypoxia tolerance: leveraging phylogeny to better understand convergent phenotypes; ontogenic studies to identify regulation of key genes underlying the elite, adult, hypoxia-tolerant physiology; and cell culture manipulations to understand gene function.

Assessment of Risks of Dioxins for Aryl Hydrocarbon Receptor-Mediated Effects in Polar Bear (Ursus maritimus) by in Vitro and in Silico Approaches

Polar bear (Ursus maritimus) populations accumulate dioxins and related compounds (DRCs) at levels that are of health concern. The toxicities of DRCs are primarily mediated via aryl hydrocarbon receptor (AHR) signaling pathway. To evaluate the sensitivity and responses to DRCs in polar bears, we assessed the activation potencies of polar bear-specific AHR (pbAHR) by DRCs through in vitro and in silico approaches. In vitro assays showed that the pbAHR was as sensitive to DRCs as C3H/lpr mouse AHR, which is well-known to be highly sensitive to DRCs. Comparison of pbAHR transactivation potencies indicated that TCDF, 2,3,4,7,8-PeCDF, and BaP exhibited high induction equivalency factors (IEFs). Considering the accumulation levels of DRCs in polar bears, PCB126 was found to be the most active inducer of pbAHR. The in vitro transactivation potencies of ligands of pbAHR showed a significant relationship with in silico ligand docking energies in a pbAHR homology model. The protein ligand interaction fingerprint (PLIF) analysis showed different interaction patterns depending on the ligands. Several amino acids which are highly conserved among mammals may be involved in species-specific responses via backbone interactions with neighboring amino acid residues which are specific to pbAHR. We document high susceptibility of polar bears to DRCs, through a mechanistic approach, for the first time.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

Divergence history and hydrothermal vent adaptation of decapod crustaceans: A mitogenomic perspective

Decapod crustaceans, such as alvinocaridid shrimps, bythograeid crabs and galatheid squat lobsters are important fauna in the hydrothermal vents and have well adapted to hydrothermal vent environments. In this study, eighteen mitochondrial genomes (mitogenomes) of hydrothermal vent decapods were used to explore the evolutionary history and their adaptation to the hydrothermal vent habitats. BI and ML algorithms produced consistent phylogeny for Decapoda. The phylogenetic relationship revealed more evolved positions for all the hydrothermal vent groups, indicating they migrated from non-vent environments, instead of the remnants of ancient hydrothermal vent species, which support the extinction/repopulation hypothesis. The divergence time estimation on the Alvinocarididae, Bythograeidae and Galatheoidea nodes are located at 75.20, 56.44 and 47.41–50.43 Ma, respectively, which refers to the Late Cretaceous origin of alvinocaridid shrimps and the Early Tertiary origin of bythograeid crabs and galatheid squat lobsters. These origin stories are thought to associate with the global deep-water anoxic/dysoxic events. Total eleven positively selected sites were detected in the mitochondrial OXPHOS genes of three lineages of hydrothermal vent decapods, suggesting a link between hydrothermal vent adaption and OXPHOS molecular biology in decapods. This study adds to the understanding of the link between mitogenome evolution and ecological adaptation to hydrothermal vent habitats in decapods.

31° South: The physiology of adaptation to arid conditions in a passerine bird

Arid environments provide ideal ground for investigating the mechanisms of adaptive evolution. High temperatures and low water availability are relentless stressors for many endotherms, including birds; yet birds persist in deserts. While physiological adaptation probably involves metabolic phenotypes, the underlying mechanisms (plasticity, genetics) are largely uncharacterized. To explore this, we took an intraspecific approach that focused on a species that is resident over a mesic to arid gradient, the Karoo scrub-robin (Cercotrichas coryphaeus). Specifically, we integrated environmental (climatic and primary productivity), physiological (metabolic rates: a measure of energy expenditure), genotypic (genetic variation underlying the machinery of energy production) and microbiome (involved in processing food from where energy is retrieved) data, to infer the mechanism of physiological adaptation. We that found the variation in energetic physiology phenotypes and gut microbiome composition are associated with environmental features as well as with variation in genes underlying energy metabolic pathways. Specifically, we identified a small list of candidate adaptive genes, some of them with known ties to relevant physiology phenotypes. Together our results suggest that selective pressures on energetic physiology mediated by genes related to energy homeostasis and possibly microbiota composition may facilitate adaptation to local conditions and provide an explanation to the high avian intraspecific divergence observed in harsh environments.

Polar bear evolution is marked by rapid changes in gene copy number in response to dietary shift

Polar bear (Ursus maritimus) and brown bear (Ursus arctos) are recently diverged species that inhabit vastly differing habitats. Thus, analysis of the polar bear and brown bear genomes represents a unique opportunity to investigate the evolutionary mechanisms and genetic underpinnings of rapid ecological adaptation in mammals. Copy number (CN) differences in genomic regions between closely related species can underlie adaptive phenotypes and this form of genetic variation has not been explored in the context of polar bear evolution. Here, we analyzed the CN profiles of 17 polar bears, 9 brown bears, and 2 black bears (Ursus americanus). We identified an average of 318 genes per individual that showed evidence of CN variation (CNV). Nearly 200 genes displayed species-specific CN differences between polar bear and brown bear species. Principal component analysis of gene CN provides strong evidence that CNV evolved rapidly in the polar bear lineage and mainly resulted in CN loss. Olfactory receptors composed 47% of CN differentiated genes, with the majority of these genes being at lower CN in the polar bear. Additionally, we found significantly fewer copies of several genes involved in fatty acid metabolism as well as AMY1B, the salivary amylase-encoding gene in the polar bear. These results suggest that natural selection shaped patterns of CNV in response to the transition from an omnivorous to primarily carnivorous diet during polar bear evolution. Our analyses of CNV shed light on the genomic underpinnings of ecological adaptation during polar bear evolution.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

State of knowledge on current exposure, fate and potential health effects of contaminants in polar bears from the circumpolar Arctic

The polar bear (Ursus maritimus) is among the Arctic species exposed to the highest concentrations of long-range transported bioaccumulative contaminants, such as halogenated organic compounds and mercury. Contaminant exposure is considered to be one of the largest threats to polar bears after the loss of their Arctic sea ice habitat due to climate change. The aim of this review is to provide a comprehensive summary of current exposure, fate, and potential health effects of contaminants in polar bears from the circumpolar Arctic required by the Circumpolar Action Plan for polar bear conservation. Overall results suggest that legacy persistent organic pollutants (POPs) including polychlorinated biphenyls, chlordanes and perfluorooctane sulfonic acid (PFOS), followed by other perfluoroalkyl compounds (e.g. carboxylic acids, PFCAs) and brominated flame retardants, are still the main compounds in polar bears. Concentrations of several legacy POPs that have been banned for decades in most parts of the world have generally declined in polar bears. Current spatial trends of contaminants vary widely between compounds and recent studies suggest increased concentrations of both POPs and PFCAs in certain subpopulations. Correlative field studies, supported by in vitro studies, suggest that contaminant exposure disrupts circulating levels of thyroid hormones and lipid metabolism, and alters neurochemistry in polar bears. Additionally, field and in vitro studies and risk assessments indicate the potential for adverse impacts to polar bear immune functions from exposure to certain contaminants.

Environmental contaminants modulate the transcriptional activity of polar bear (Ursus maritimus) and human peroxisome proliferator-activated receptor alpha (PPARA)

Peroxisome proliferator-activated receptor alfa (PPARA/NR1C1) is a ligand activated nuclear receptor that is a key regulator of lipid metabolism in tissues with high fatty acid catabolism such as the liver. Here, we cloned PPARA from polar bear liver tissue and studied in vitro transactivation of polar bear and human PPARA by environmental contaminants using a luciferase reporter assay. Six hinge and ligand-binding domain amino acids have been substituted in polar bear PPARA compared to human PPARA. Perfluorocarboxylic acids (PFCA) and perfluorosulfonic acids induced the transcriptional activity of both human and polar bear PPARA. The most abundant PFCA in polar bear tissue, perfluorononanoate, increased polar bear PPARA-mediated luciferase activity to a level comparable to that of the potent PPARA agonist WY-14643 (~8-fold, 25 μM). Several brominated flame retardants were weak agonists of human and polar bear PPARA. While single exposures to polychlorinated biphenyls did not, or only slightly, increase the transcriptional activity of PPARA, a technical mixture of PCBs (Aroclor 1254) strongly induced the transcriptional activity of human (~8-fold) and polar bear PPARA (~22-fold). Polar bear PPARA was both quantitatively and qualitatively more susceptible than human PPARA to transactivation by less lipophilic compounds.

Landscape genomics: natural selection drives the evolution of mitogenome in penguins

Background

Mitochondria play a key role in the balance of energy and heat production, and therefore the mitochondrial genome is under natural selection by environmental temperature and food availability, since starvation can generate more efficient coupling of energy production. However, selection over mitochondrial DNA (mtDNA) genes has usually been evaluated at the population level. We sequenced by NGS 12 mitogenomes and with four published genomes, assessed genetic variation in ten penguin species distributed from the equator to Antarctica. Signatures of selection of 13 mitochondrial protein-coding genes were evaluated by comparing among species within and among genera (Spheniscus, Pygoscelis, Eudyptula, Eudyptes and Aptenodytes). The genetic data were correlated with environmental data obtained through remote sensing (sea surface temperature [SST], chlorophyll levels [Chl] and a combination of SST and Chl [COM]) through the distribution of these species.

Results

We identified the complete mtDNA genomes of several penguin species, including ND6 and 8 tRNAs on the light strand and 12 protein coding genes, 14 tRNAs and two rRNAs positioned on the heavy strand. The highest diversity was found in NADH dehydrogenase genes and the lowest in COX genes. The lowest evolutionary divergence among species was between Humboldt (Spheniscus humboldti) and Galapagos (S. mendiculus) penguins (0.004), while the highest was observed between little penguin (Eudyptula minor) and Adélie penguin (Pygoscelis adeliae) (0.097). We identified a signature of purifying selection (Ka/Ks < 1) across the mitochondrial genome, which is consistent with the hypothesis that purifying selection is constraining mitogenome evolution to maintain Oxidative phosphorylation (OXPHOS) proteins and functionality. Pairwise species maximum-likelihood analyses of selection at codon sites suggest positive selection has occurred on ATP8 (Fixed-Effects Likelihood, FEL) and ND4 (Single Likelihood Ancestral Counting, SLAC) in all penguins. In contrast, COX1 had a signature of strong negative selection. ND4 Ka/Ks ratios were highly correlated with SST (Mantel, p-value: 0.0001; GLM, p-value: 0.00001) and thus may be related to climate adaptation throughout penguin speciation.

Conclusions

These results identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment. Such knowledge may be particularly useful for developing predictive models of how these species may respond to severe climatic changes in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4424-9) contains supplementary material, which is available to authorized users.

Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes

Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.

De novo assembling and primary analysis of genome and transcriptome of gray whale Eschrichtius robustus

Background

Gray whale, Eschrichtius robustus (E. robustus), is a single member of the family Eschrichtiidae, which is considered to be the most primitive in the class Cetacea. Gray whale is often described as a “living fossil”. It is adapted to extreme marine conditions and has a high life expectancy (77 years). The assembly of a gray whale genome and transcriptome will allow to carry out further studies of whale evolution, longevity, and resistance to extreme environment.

Results

In this work, we report the first de novo assembly and primary analysis of the E. robustus genome and transcriptome based on kidney and liver samples. The presented draft genome assembly is complete by 55% in terms of a total genome length, but only by 24% in terms of the BUSCO complete gene groups, although 10,895 genes were identified. Transcriptome annotation and comparison with other whale species revealed robust expression of DNA repair and hypoxia-response genes, which is expected for whales.

Conclusions

This preliminary study of the gray whale genome and transcriptome provides new data to better understand the whale evolution and the mechanisms of their adaptation to the hypoxic conditions.

Electronic supplementary material

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

Adaptive Patterns of Mitogenome Evolution Are Associated with the Loss of Shell Scutes in Turtles

The mitochondrial genome encodes several protein components of the oxidative phosphorylation (OXPHOS) pathway and is critical for aerobic respiration. These proteins have evolved adaptively in many taxa, but linking molecular-level patterns with higher-level attributes (e.g., morphology, physiology) remains a challenge. Turtles are a promising system for exploring mitochondrial genome evolution as different species face distinct respiratory challenges and employ multiple strategies for ensuring efficient respiration. One prominent adaptation to a highly aquatic lifestyle in turtles is the secondary loss of keratenized shell scutes (i.e., soft-shells), which is associated with enhanced swimming ability and, in some species, cutaneous respiration. We used codon models to examine patterns of selection on mitochondrial protein-coding genes along the three turtle lineages that independently evolved soft-shells. We found strong evidence for positive selection along the branches leading to the pig-nosed turtle (Carettochelys insculpta) and the softshells clade (Trionychidae), but only weak evidence for the leatherback (Dermochelys coriacea) branch. Positively selected sites were found to be particularly prevalent in OXPHOS Complex I proteins, especially subunit ND2, along both positively selected lineages, consistent with convergent adaptive evolution. Structural analysis showed that many of the identified sites are within key regions or near residues involved in proton transport, indicating that positive selection may have precipitated substantial changes in mitochondrial function. Overall, our study provides evidence that physiological challenges associated with adaptation to a highly aquatic lifestyle have shaped the evolution of the turtle mitochondrial genome in a lineage-specific manner.

Multiple-stressor effects in an apex predator: combined influence of pollutants and sea ice decline on lipid metabolism in polar bears

There is growing evidence from experimental and human epidemiological studies that many pollutants can disrupt lipid metabolism. In Arctic wildlife, the occurrence of such compounds could have serious consequences for seasonal feeders. We set out to study whether organohalogenated compounds (OHCs) could cause disruption of energy metabolism in female polar bears (Ursus maritimus) from Svalbard, Norway (n = 112). We analyzed biomarkers of energy metabolism including the abundance profiles of nine lipid-related genes, fatty acid (FA) synthesis and elongation indices in adipose tissue, and concentrations of lipid-related variables in plasma (cholesterol, high-density lipoprotein, triglycerides). Furthermore, the plasma metabolome and lipidome were characterized by low molecular weight metabolites and lipid fingerprinting, respectively. Polychlorinated biphenyls, chlordanes, brominated diphenyl ethers and perfluoroalkyl substances were significantly related to biomarkers involved in lipid accumulation, FA metabolism, insulin utilization, and cholesterol homeostasis. Moreover, the effects of pollutants were measurable at the metabolome and lipidome levels. Our results indicate that several OHCs affect lipid biosynthesis and catabolism in female polar bears. Furthermore, these effects were more pronounced when combined with reduced sea ice extent and thickness, suggesting that climate-driven sea ice decline and OHCs have synergistic negative effects on polar bears.

Life in the fat lane: seasonal regulation of insulin sensitivity, food intake, and adipose biology in brown bears

Grizzly bears (Ursus arctos horribilis) have evolved remarkable metabolic adaptations including enormous fat accumulation during the active season followed by fasting during hibernation. However, these fluctuations in body mass do not cause the same harmful effects associated with obesity in humans. To better understand these seasonal transitions, we performed insulin and glucose tolerance tests in captive grizzly bears, characterized the annual profiles of circulating adipokines, and tested the anorectic effects of centrally administered leptin at different times of the year. We also used bear gluteal adipocyte cultures to test insulin and beta-adrenergic sensitivity in vitro. Bears were insulin resistant during hibernation but were sensitive during the spring and fall active periods. Hibernating bears remained euglycemic, possibly due to hyperinsulinemia and hyperglucagonemia. Adipokine concentrations were relatively low throughout the active season but peaked in mid-October prior to hibernation when fat content was greatest. Serum glycerol was highest during hibernation, indicating ongoing lipolysis. Centrally administered leptin reduced food intake in October, but not in August, revealing seasonal variation in the brain's sensitivity to its anorectic effects. This was supported by strong phosphorylated signal transducer and activator of transcription 3 labeling within the hypothalamus of hibernating bears; labeling virtually disappeared in active bears. Adipocytes collected during hibernation were insulin resistant when cultured with hibernation serum but became sensitive when cultured with active season serum. Heat treatment of active serum blocked much of this action. Clarifying the cellular mechanisms responsible for the physiology of hibernating bears may inform new treatments for metabolic disorders.

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea)

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni. Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

A protocol for the isolation and cultivation of brown bear (Ursus arctos) adipocytes

Brown bears (Ursus arctos) exhibit hyperphagia each fall and can become obese in preparation for hibernation. They do this without displaying the physiological problems typically seen in obese humans, such as Type 2 diabetes and heart disease. The study of brown bear hibernation biology could therefore aid in the development of novel methods for combating metabolic diseases. To this end, we isolated mesenchymal stem cells from subcutaneous fat biopsies, and culture methods were developed to differentiate these into the adipogenic lineage. Biopsies were taken from 8 captive male (N = 6) and female (N = 2) brown bears, ages 2-12 years. Plastic adherent, fibroblast-like cells were proliferated and subsequently cryopreserved or differentiated. Differentiation conditions were optimized with respect to fetal bovine serum content and time spent in differentiation medium. Cultures were characterized through immunostaining, RT-qPCR, and Oil red O staining to quantify lipid accumulation. Adiponectin, leptin, and glycerol medium concentrations were also determined over the course of differentiation. The culturing protocol succeeded in generating hormone-sensitive lipase-expressing, lipid-producing white-type adipocytes (UCP1 negative). Serum concentration and time of exposure to differentiation medium were both positively related to lipid production. Cells cultured to low passage numbers retained similar lipid production and expression of lipid markers PLIN2 and FABP4. Ultimately, the protocols described here may be useful to biologists in the field investigating the health of wild bear populations and could potentially increase our understanding of metabolic disorders in humans.

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population-specific mitonuclear coadaptations create between-population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.

Do relaxed selection and habitat temperature facilitate biased mitogenomic introgression in a narrowly endemic fish?

Introgression might be exceptionally common during the evolution of narrowly endemic species. For instance, in the springs of the small and isolated Cuatro Ciénegas Valley, the mitogenome of the cichlid fish Herichthys cyanoguttatus could be rapidly introgressing into populations of the trophically polymorphic H. minckleyi. We used a combination of genetic and environmental data to examine the factors associated with this mitochondrial introgression. A reduced representation library of over 6220 single nucleotide polymorphisms (SNPs) from the nuclear genome showed that mitochondrial introgression into H. minckleyi is biased relative to the amount of nuclear introgression. SNP assignment probabilities also indicated that cichlids with more hybrid ancestry are not more commonly female providing no support for asymmetric backcrossing or hybrid‐induced sex‐ratio distortion in generating the bias in mitochondrial introgression. Smaller effective population size in H. minckleyi inferred from the SNPs coupled with sequences of all 13 mitochondrial proteins suggests that relaxed selection on the mitogenome could be facilitating the introgression of “H. cyanoguttatus” haplotypes. Additionally, we showed that springs with colder temperatures had greater amounts of mitochondrial introgression from H. cyanoguttatus. Relaxed selection in H. minckleyi coupled with temperature‐related molecular adaptation could be facilitating mitogenomic introgression into H. minckleyi.

Divergent ecological histories of two sister Antarctic krill species led to contrasted patterns of genetic diversity in their heat-shock protein (hsp70) arsenal

The Arctic and the Antarctic Peninsula are currently experiencing some of the most rapid rates of ocean warming on the planet. This raises the question of how the initial adaptation to extreme cold temperatures was put in place and whether or not directional selection has led to the loss of genetic variation at key adaptive systems, and thus polar species’ (re)adaptability to higher temperatures. In the Southern Ocean, krill represents the most abundant fauna and is a critical member at the base of the Antarctic food web. To better understand the role of selection in shaping current patterns of polymorphisms, we examined genetic diversity of the cox‐1 and hsp70 genes by comparing two closely related species of Euphausiid that differ in ecology. Results on mtcox‐1 agreed with previous studies, indicating high and similar effective population sizes. However, a coalescent‐based approach on hsp70 genes highlighted the role of positive selection and past demographic changes in their recent evolution. Firstly, some form of balancing selection was acting on the inducible isoform C, which reflected the maintenance of an ancestral adaptive polymorphism in both species. Secondly, E. crystallorophias seems to have lost most of its hsp70 diversity because of a population crash and/or directional selection to cold. Nonsynonymous diversities were always greater in E. superba, suggesting that it might have evolved under more heterogeneous conditions. This can be linked to species’ ecology with E. superba living in more variable pelagic conditions, while E. crystallorophias is strictly associated with continental shelves and sea ice.

Heterogeneous natural selection on oxidative phosphorylation genes among fishes with extreme high and low aerobic performance

Background

Oxidative phosphorylation (OXPHOS) is the primary source of ATP in eukaryotes and serves as a mechanistic link between variation in genotypes and energetic phenotypes. While several physiological and anatomical factors may lead to increased aerobic capacity, variation in OXPHOS proteins may influence OXPHOS efficiency and facilitate adaptation in organisms with varied energy demands. Although there is evidence that natural selection acts on OXPHOS genes, the focus has been on detection of directional (positive) selection on specific phylogenetic branches where traits that increase energetic demands appear to have evolved. We examined patterns of selection in a broader evolutionary context, i.e., on multiple lineages of fishes with extreme high and low aerobic performance.

Results

We found that patterns of natural selection on mitochondrial OXPHOS genes are complex among fishes with different swimming performance. Positive selection is not consistently associated with high performance taxa and appears to be strongest on lineages containing low performance taxa. In contrast, within high performance lineages, purifying (negative) selection appears to predominate.

Conclusions

We provide evidence that selection on OXPHOS varies in both form and intensity within and among lineages through evolutionary time. These results provide evidence for fluctuating selection on OXPHOS associated with divergence in aerobic performance. However, in contrast to previous studies, positive selection was strongest on low performance taxa suggesting that adaptation of OXPHOS involves many factors beyond enhancing ATP production in high performance taxa. The broader pattern indicates a complex interplay between organismal adaptations, ATP demand, and OXPHOS function.

Electronic supplementary material

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

Considerations for the use of transcriptomics in identifying the ‘genes that matter’ for environmental adaptation

Transcriptomics has emerged as a powerful approach for exploring physiological responses to the environment. However, like any other experimental approach, transcriptomics has its limitations. Transcriptomics has been criticized as an inappropriate method to identify genes with large impacts on adaptive responses to the environment because: (1) genes with large impacts on fitness are rare; (2) a large change in gene expression does not necessarily equate to a large effect on fitness; and (3) protein activity is most relevant to fitness, and mRNA abundance is an unreliable indicator of protein activity. In this review, these criticisms are re-evaluated in the context of recent systems-level experiments that provide new insight into the relationship between gene expression and fitness during environmental stress. In general, these criticisms remain valid today, and indicate that exclusively using transcriptomics to screen for genes that underlie environmental adaptation will overlook constitutively expressed regulatory genes that play major roles in setting tolerance limits. Standard practices in transcriptomic data analysis pipelines may also be limiting insight by prioritizing highly differentially expressed and conserved genes over those genes that undergo moderate fold-changes and cannot be annotated. While these data certainly do not undermine the continued and widespread use of transcriptomics within environmental physiology, they do highlight the types of research questions for which transcriptomics is best suited and the need for more gene functional analyses. Such information is pertinent at a time when transcriptomics has become increasingly tractable and many researchers may be contemplating integrating transcriptomics into their research programs.

Introgressive hybridization: brown bears as vectors for polar bear alleles

The dynamics and consequences of introgression can inform about numerous evolutionary processes. Biologists have therefore long been interested in hybridization. One challenge, however, lies in the identification of nonadmixed genotypes that can serve as a baseline for accurate quantification of admixture. In this issue of Molecular Ecology, Cahill et al. (2015) analyse a genomic data set of 28 polar bears, eight brown bears and one American black bear. Polar bear alleles are found to be introgressed into brown bears not only near a previously identified admixture zone on the Alaskan Admiralty, Baranof and Chichagof (ABC) Islands, but also far into the North American mainland. Elegantly contrasting admixture levels at autosomal and X chromosomal markers, Cahill and colleagues infer that male-biased dispersal has spread these introgressed alleles away from the Late Pleistocene contact zone. Compared to a previous study on the ABC Island population in which an Alaskan brown bear served as a putatively admixture-free reference, Cahill et al. (2015) utilize a newly sequenced Swedish brown bear as admixture baseline. This approach reveals that brown bears have been impacted by introgression from polar bears to a larger extent (up to 8.8% of their genome), than previously known, including the bear that had previously served as admixture baseline. No evidence for introgression of brown bear into polar bear is found, which the authors argue could be a consequence of selection. Besides adding new exciting pieces to the puzzle of polar/brown bear evolutionary history, the study by Cahill and colleagues highlights that wildlife genomics is moving from analysing single genomes towards a landscape genomics approach.