Genome assembly and gene expression in the American black bear provides new insights into the renal response to hibernation

Cognitive Impairment in Chronic Kidney Disease Is Associated with Glymphatic System Dysfunction

Introduction

This study was designed to explore the associations between impaired cognition in chronic kidney disease (CKD) patients and the dysfunction of the glymphatic system.

Method

Data were obtained from 77 CKD patients and 50 age-matched healthy control individuals from the First Affiliated Hospital of Zhengzhou University. CKD patients were stratified into with and without impaired cognitive function. T2-weighted magnetic resonance imaging results were used to assess area ratios for the perivascular space and ventricles in participants, while the Montreal Cognitive Assessment and the Mini-Mental State Examination were employed to measure cognitive function. Correlations between the perivascular space or ventricle area ratios and cognitive impairment were assessed in CKD patients.

Results

Significant increases in the burden of enlarged perivascular spaces in the frontal cortex and basal ganglia were observed in CKD patients with cognitive impairment relative to those without such impairment, with a concomitant increase in analyzed ventricle area ratios. Enlarged perivascular spaces in the frontal cortex, basal ganglia and increased area ratios of lateral ventricles and 4th ventricle exhibited relatively high sensitivity and specificity as means of differing between the CKD patients with and without cognitive impairment.

Conclusion

These results indicate that the burden of enlarged perivascular spaces in the frontal cortex and basal ganglia and increases in ventricle area ratio values may offer utility as biomarkers that can aid in detection of even mild cognitive decline in individuals with CKD. The dysfunction of the glymphatic system may play a key role in the pathogenesis of CKD-related cognitive impairment.

Ancient bears provide insights into Pleistocene ice age refugia in Southeast Alaska

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

Genetic architecture and evolution of color variation in American black bears

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Sequencing and assembling bear genomes: the bare necessities

Unique genetic adaptations are present in bears of every species across the world. From (nearly) shutting down important organs during hibernation to preventing harm from lifestyles that could easily cause metabolic diseases in humans, bears may hold the answer to various human ailments. However, only a few of these unique traits are currently being investigated at the molecular level, partly because of the lack of necessary tools. One of these tools is well-annotated genome assemblies from the different, extant bear species. These reference genomes are needed to allow us to identify differences in genetic variants, isoforms, gene expression, and genomic features such as transposons and identify those that are associated with biomedical-relevant traits. In this review we assess the current state of the genome assemblies of the eight different bear species, discuss current gaps, and the future benefits these reference genomes may have in informing human biomedical applications, while at the same time improving bear conservation efforts.

A beary good genome: Haplotype-resolved, chromosome-level assembly of the brown bear (Ursus arctos)

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 Megabases (Mb), respectively. Benchmarking Universal Single-Copy Ortholog (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 to the polar bear, the brown bear is highly syntenic with the polar bear, and our phylogenetic analysis of 7,246 single-copy orthologs 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.

Uncovering the enigmatic evolution of bears in greater depth: The hybrid origin of the Asiatic black bear

Bears are fascinating mammals because of their complex pattern of speciation and rapid evolution of distinct phenotypes. Interspecific hybridization has been common and has shaped the complex evolutionary history of bears. In this study, based on the largest population-level genomic dataset to date involving all Ursinae species and recently developed methods for detecting hybrid speciation, we provide explicit evidence for the hybrid origin of Asiatic black bears, which arose through historical hybridization between the ancestor of polar bear/brown bear/American black bears and the ancestor of sun bear/sloth bears. This was inferred to have occurred soon after the divergence of the two parental lineages in Eurasia due to climate-driven population expansion and dispersal. In addition, we found that the intermediate body size of this hybrid species arose from its combination of relevant genes derived from two parental lineages of contrasting sizes. This and alternate fixation of numerous other loci that had diverged between parental lineages may have initiated the reproductive isolation of the Asiatic black bear from its two parents. Our study sheds further light on the evolutionary history of bears and documents the importance of hybridization in new species formation and phenotypic evolution in mammals.

Transformation of Mitochondrial Architecture and Dynamics in the Chinese Soft-Shelled Turtle (Pelodiscus sinensis) During Hibernation

Hibernation is a biological status during which hibernating animals acclimatize themselves to reduced energy consumption through extreme but governed decline in self-metabolism. The role of mitochondria (Mt) in metabolic suppression during hibernation has already been elaborated in different organs and species. Nonetheless, the concretely changing process of mitochondrial architecture and the mechanism underlying this transformation during hibernation remains unclear. Herein, the present study was aimed at clarifying the detailed alteration of mitochondrial morphology and its potential role in the Chinese soft-shelled turtle (Pelodiscus sinensis) during different stages of hibernation. Compared with the nonhibernation period, the mitochondrial architecture was changing from round to crescent, and lipid droplet (LD)/Mt interaction was enhanced during hibernation, as observed by transmission electron microscopy (TEM). Further ultrastructural analysis uncovered that mitochondrial fusion was promptly accelerated in the early stage of hibernation, followed by mitochondrial fission in the middle stage, and mitophagy was boosted in the late stage. Moreover, gene and protein expression related to mitochondrial fusion, fission, and mitophagy accorded closely with the mitochondrial ultrastructural changes in different stages of hibernation. Taken together, our results clarified that the transformation of mitochondrial architecture and mitochondrial dynamics are of vital importance in maintaining internal environment homeostasis of Pelodiscus sinensis.

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.

Conservation of chromatin conformation in carnivores

We found the three-dimensional (3D) structure of chromatin at the chromosome level to be highly conserved for more than 50 million y of carnivore evolution. Intrachromosomal contacts were maintained even after chromosome rearrangements within carnivore lineages, demonstrating that the maintenance of 3D chromatin architecture is essential for conserved genome functions. These discoveries enabled the identification of orthologous chromosomal DNA segments among related species, a method we call 3D comparative scaffotyping. The method has application for putative chromosome assignment of chromosome-scale DNA sequence scaffolds produced by de novo genome sequencing. Broadly applied to biodiversity genome sequencing efforts, the approach can reduce costs associated with karyotyping and the physical mapping of DNA segments to chromosomes.

High throughput chromatin conformation capture (Hi-C) of leukocyte DNA was used to investigate the evolutionary stability of chromatin conformation at the chromosomal level in 11 species from three carnivore families: Felidae, Canidae, and Ursidae. Chromosome-scale scaffolds (C-scaffolds) of each species were initially used for whole-genome alignment to a reference genome within each family. This approach established putative orthologous relationships between C-scaffolds among the different species. Hi-C contact maps for all C-scaffolds were then visually compared and found to be distinct for a given reference chromosome or C-scaffold within a species and indistinguishable for orthologous C-scaffolds having a 1:1 relationship within a family. The visual patterns within families were strongly supported by eigenvectors from the Hi-C contact maps. Analysis of Hi-C contact maps and eigenvectors across the three carnivore families revealed that most cross-family orthologous subchromosomal fragments have a conserved three-dimensional (3D) chromatin structure and thus have been under strong evolutionary constraint for ∼54 My of carnivore evolution. The most pronounced differences in chromatin conformation were observed for the X chromosome and the red fox genome, whose chromosomes have undergone extensive rearrangements relative to other canids. We also demonstrate that Hi-C contact map pattern analysis can be used to accurately identify orthologous relationships between C-scaffolds and chromosomes, a method we termed “3D comparative scaffotyping.” This method provides a powerful means for estimating karyotypes in de novo sequenced species that have unknown karyotype and no physical mapping information.

An Annotated Draft Genome for the Andean Bear, Tremarctos ornatus

The Andean bear is the only extant member of the Tremarctine subfamily and the only extant ursid species to inhabit South America. Here, we present an annotated de novo assembly of a nuclear genome from a captive-born female Andean bear, Mischief, generated using a combination of short and long DNA and RNA reads. Our final assembly has a length of 2.23 Gb, and a scaffold N50 of 21.12 Mb, contig N50 of 23.5 kb, and BUSCO score of 88%. The Andean bear genome will be a useful resource for exploring the complex phylogenetic history of extinct and extant bear species and for future population genetics studies of Andean bears.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

A Translational Model for Venous Thromboembolism: MicroRNA Expression in Hibernating Black Bears

BACKGROUND

Hibernating American black bears have significantly different clotting parameters than their summer active counterparts, affording them protection against venous thromboembolism during prolonged periods of immobility. We sought to evaluate if significant differences exist between the expression of microRNAs in the plasma of hibernating black bears compared with their summer active counterparts, potentially contributing to differences in hemostasis during hibernation.

MATERIALS AND METHODS

MicroRNA sequencing was assessed in plasma from 21 American black bears in summer active (n = 11) and hibernating states (n = 10), and microRNA signatures during hibernating and active state were established using both bear and human genome. MicroRNA targets were predicted using messenger RNA (mRNA) transcripts from black bear kidney cells. In vitro studies were performed to confirm the relationship between identified microRNAs and mRNA expression, using artificial microRNA and human liver cells.

RESULTS

Using the bear genome, we identified 15 microRNAs differentially expressed in the plasma of hibernating black bears. Of these microRNAs, three were significantly downregulated (miR-141-3p, miR-200a-3p, and miR-200c-3p), were predicted to target SERPINC1, the gene for antithrombin, and demonstrated regulatory control of the gene mRNA expression in cell studies.

CONCLUSIONS

Our findings suggest that the hibernating black bears' ability to maintain hemostasis and achieve protection from venous thromboembolism during prolonged periods of immobility may be due to changes in microRNA signatures and possible upregulation of antithrombin expression.

Insights in the regulation of trimetylamine N-oxide production using a comparative biomimetic approach suggest a metabolic switch in hibernating bears

Experimental studies suggest involvement of trimethylamine N-oxide (TMAO) in the aetiology of cardiometabolic diseases and chronic kidney disease (CKD), in part via metabolism of ingested food. Using a comparative biomimetic approach, we have investigated circulating levels of the gut metabolites betaine, choline, and TMAO in human CKD, across animal species as well as during hibernation in two animal species. Betaine, choline, and TMAO levels were associated with renal function in humans and differed significantly across animal species. Free-ranging brown bears showed a distinct regulation pattern with an increase in betaine (422%) and choline (18%) levels during hibernation, but exhibited undetectable levels of TMAO. Free-ranging brown bears had higher betaine, lower choline, and undetectable TMAO levels compared to captive brown bears. Endogenously produced betaine may protect bears and garden dormice during the vulnerable hibernating period. Carnivorous eating habits are linked to TMAO levels in the animal kingdom. Captivity may alter the microbiota and cause a subsequent increase of TMAO production. Since free-ranging bears seems to turn on a metabolic switch that shunts choline to generate betaine instead of TMAO, characterisation and understanding of such an adaptive switch could hold clues for novel treatment options in burden of lifestyle diseases, such as CKD.

Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role

BACKGROUND

Polar bears are uniquely adapted to an Arctic existence. Since their relatively recent divergence from their closest living relative, brown bears, less than 500,000 years ago, the species has evolved an array of novel traits suited to its Arctic lifestyle. Previous studies sought to uncover the genomic underpinnings of these unique characteristics, and disclosed the genes showing the strongest signal of positive selection in the polar bear lineage. Here, we survey a comprehensive dataset of 109 polar bear and 33 brown bear genomes to investigate the genomic variants within these top genes present in each species. Specifically, we investigate whether fixed homozygous variants in polar bears derived from selection on standing variation in the ancestral gene pool or on de novo mutation in the polar bear lineage.

RESULTS

We find that a large number of sites fixed in polar bears are biallelic in brown bears, suggesting selection on standing variation. Moreover, we uncover sites in which polar bears are fixed for a derived allele while brown bears are fixed for the ancestral allele, which we suggest may be a signal of de novo mutation in the polar bear lineage.

CONCLUSIONS

Our findings suggest that, among other mechanisms, natural selection acting on changes in genes derived from a combination of variation already in the ancestral gene pool, and from de novo missense mutations in the polar bear lineage, may have enabled the rapid adaptation of polar bears to their new Arctic environment.

Using Genetic and Species Diversity to Tackle Kidney Disease

Progress in the identification of causal genes and understanding of the mechanism underlying kidney disease is hindered by the almost exclusive use of a few animal models with restrictive monogenic backgrounds that may be more resistant to kidney disease compared with humans and, therefore, poor models. Exploring the large genetic diversity in classical animal models, such as mice and rats, and leveraging species diversity will allow us to use the genetic advantages of zebrafish, Drosophila, and other species, to develop both new animal models that are more relevant to the study of human kidney disease and potential therapies.

Advances and applications of environmental stress adaptation research

Evolution has produced animals that survive extreme fluctuations in environmental conditions including freezing temperatures, anoxia, desiccating conditions, and prolonged periods without food. For example, the wood frog survives whole-body freezing every winter, arresting all gross physiological functions, but recovers functions upon thawing in the spring. Likewise, many small mammals hibernate for months at a time with minimal metabolic activity, organ perfusion, and movement, yet do not suffer significant muscle atrophy upon arousal. These conditions and the biochemical adaptations employed to deal with them can be viewed as Nature's answer to problems that humans wish to answer, particularly in a biomedical context. This review focuses on recent advances in the field of animal environmental stress adaptation, starting with an emphasis on new areas of research such as epigenetics and microRNA. We then examine new and emerging technologies such as genome editing, novel sequencing applications, and single cell analysis and how these can push us closer to a deeper understanding of biochemical adaptation. Next, evaluate the potential contributions of new high-throughput technologies (e.g. next-generation sequencing, mass spectrometry proteomics) to better understanding the adaptations that support these extreme phenotypes. Concluding, we examine some of the human applications that can be gained from understanding the principles of biochemical adaptation including organ preservation and treatments for conditions such as ischemic stroke and muscle disuse atrophy.

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.

Molecular interactions underpinning the phenotype of hibernation in mammals

Mammals maintain a constant warm body temperature, facilitating a wide variety of metabolic reactions. Mammals that hibernate have the ability to slow their metabolism, which in turn reduces their body temperature and leads to a state of hypothermic torpor. For this metabolic rate reduction to occur on a whole-body scale, molecular interactions that change the physiology of cells, tissues and organs are required, resulting in a major departure from normal mammalian homeostasis. The aim of this Review is to cover recent advances in the molecular biology of mammalian hibernation, including the role of small molecules, seasonal changes in gene expression, cold-inducible RNA-binding proteins, the somatosensory system and emerging information on hibernating primates. To underscore the importance of differential gene expression across the hibernation cycle, mRNA levels for 14,261 ground squirrel genes during periods of activity and torpor are made available for several tissues via an interactive transcriptome browser. This Review also addresses recent findings on molecular interactions responsible for multi-day survival of near-freezing body temperatures, single-digit heart rates and a slowed metabolism that greatly reduces oxygen consumption. A better understanding of how natural hibernators survive these physiological extremes is beginning to lead to innovations in human medicine.