Transcriptomic Responses of the Heart and Brain to Anoxia in the Western Painted Turtle

Establishment and characterization of turtle liver organoids provides a potential model to decode their unique adaptations

Painted turtles are remarkable for their freeze tolerance and supercooling ability along with their associated resilience to hypoxia/anoxia and oxidative stress, rendering them an ideal biomedical model for hypoxia-induced injuries (including strokes), tissue cooling during surgeries, and organ cryopreservation. Yet, such research is hindered by their seasonal reproduction and slow maturation. Here we developed and characterized adult stem cell-derived turtle liver organoids (3D self-assembled in vitro structures) from painted, snapping, and spiny softshell turtles spanning ~175My of evolution, with a subset cryopreserved. This development is, to the best of our knowledge, a first for this vertebrate Order, and complements the only other non-avian reptile organoids from snake venom glands. Preliminary characterization, including morphological, transcriptomic, and proteomic analyses, revealed organoids enriched in cholangiocytes. Deriving organoids from distant turtles and life stages demonstrates that our techniques are broadly applicable to chelonians, permitting the development of functional genomic tools currently lacking in herpetological research. Such platform could potentially support studies including genome-to-phenome mapping, gene function, genome architecture, and adaptive responses to climate change, with implications for ecological, evolutionary, and biomedical research.

The vascular gene Apold1 is dispensable for normal development but controls angiogenesis under pathological conditions

The molecular mechanisms of angiogenesis have been intensely studied, but many genes that control endothelial behavior and fate still need to be described. Here, we characterize the role of Apold1 (Apolipoprotein L domain containing 1) in angiogenesis in vivo and in vitro. Single-cell analyses reveal that - across tissues - the expression of Apold1 is restricted to the vasculature and that Apold1 expression in endothelial cells (ECs) is highly sensitive to environmental factors. Using Apold1-/- mice, we find that Apold1 is dispensable for development and does not affect postnatal retinal angiogenesis nor alters the vascular network in adult brain and muscle. However, when exposed to ischemic conditions following photothrombotic stroke as well as femoral artery ligation, Apold1-/- mice display dramatic impairments in recovery and revascularization. We also find that human tumor endothelial cells express strikingly higher levels of Apold1 and that Apold1 deletion in mice stunts the growth of subcutaneous B16 melanoma tumors, which have smaller and poorly perfused vessels. Mechanistically, Apold1 is activated in ECs upon growth factor stimulation as well as in hypoxia, and Apold1 intrinsically controls EC proliferation but not migration. Our data demonstrate that Apold1 is a key regulator of angiogenesis in pathological settings, whereas it does not affect developmental angiogenesis, thus making it a promising candidate for clinical investigation.

Gene expression of hypoxia-inducible factor (HIF), HIF regulators, and putative HIF targets in ventricle and telencephalon of Trachemys scripta acclimated to 21 °C or 5 °C and exposed to normoxia, anoxia or reoxygenation

In anoxia-sensitive mammals, hypoxia inducible factor (HIF) promotes cellular survival in hypoxia, but also tumorigenesis. By comparison, anoxia-tolerant vertebrates likely need to circumvent a prolonged upregulation of HIF to survive long-term anoxia, making them attractive biomedical models for investigating HIF regulation. To lend insight into the role of HIF in anoxic Trachemys scripta ventricle and telencephalon, 21 °C- and 5 °C-acclimated turtles were exposed to normoxia, anoxia (24 h at 21 °C; 24 h or 14 d at 5 °C) or anoxia + reoxygenation and the gene expression of HIF-1α (hif1a) and HIF-2α (hif2a), two regulators of HIF, and eleven putative downstream targets of HIF quantified by qPCR. Changes in gene expression with anoxia at 21 °C differentially aligned with a circumvention of HIF activity. Whereas hif1a and hif2a expression was unaffected in ventricle and telencephalon, and BCL2 interacting protein 3 gene expression reduced by 30% in telencephalon, gene expression of vascular endothelial growth factor-A increased in ventricle (4.5-fold) and telencephalon (1.5-fold), and hexokinase 1 (2-fold) and hexokinase 2 (3-fold) gene expression increased in ventricle. At 5 °C, the pattern of gene expression in ventricle or telencephalon was unaltered with oxygenation state. However, cold acclimation in normoxia induced downregulation of HIF-1α, HIF-2α, and HIF target gene expression in telencephalon. Overall, the findings lend support to the postulation that prolonged activation of HIF is counterproductive for long-term anoxia survival. Nevertheless, quantification of the effect of anoxia and acclimation temperature on HIF binding activity and regulation at the protein level are needed to provide a strong scientific framework whereby new strategies for oxygen related pathologies can be developed.

Hypoxia Tolerant Species: The Wisdom of Nature Translated into Targets for Stroke Therapy

Human neurons rapidly die after ischemia and current therapies for stroke management are limited to restoration of blood flow to prevent further brain damage. Thrombolytics and mechanical thrombectomy are the available reperfusion treatments, but most of the patients remain untreated. Neuroprotective therapies focused on treating the pathogenic cascade of the disease have widely failed. However, many animal species demonstrate that neurons can survive the lack of oxygen for extended periods of time. Here, we reviewed the physiological and molecular pathways inherent to tolerant species that have been described to contribute to hypoxia tolerance. Among them, Foxo3 and Eif5A were reported to mediate anoxic survival in Drosophila and Caenorhabditis elegans, respectively, and those results were confirmed in experimental models of stroke. In humans however, the multiple mechanisms involved in brain cell death after a stroke causes translation difficulties to arise making necessary a timely and coordinated control of the pathological changes. We propose here that, if we were able to plagiarize such natural hypoxia tolerance through drugs combined in a pharmacological cocktail it would open new therapeutic opportunities for stroke and likely, for other hypoxic conditions.

Cytoskeletal Arrest: An Anoxia Tolerance Mechanism

Polymerization of actin filaments and microtubules constitutes a ubiquitous demand for cellular adenosine-5′-triphosphate (ATP) and guanosine-5′-triphosphate (GTP). In anoxia-tolerant animals, ATP consumption is minimized during overwintering conditions, but little is known about the role of cell structure in anoxia tolerance. Studies of overwintering mammals have revealed that microtubule stability in neurites is reduced at low temperature, resulting in withdrawal of neurites and reduced abundance of excitatory synapses. Literature for turtles is consistent with a similar downregulation of peripheral cytoskeletal activity in brain and liver during anoxic overwintering. Downregulation of actin dynamics, as well as modification to microtubule organization, may play vital roles in facilitating anoxia tolerance. Mitochondrial calcium release occurs during anoxia in turtle neurons, and subsequent activation of calcium-binding proteins likely regulates cytoskeletal stability. Production of reactive oxygen species (ROS) formation can lead to catastrophic cytoskeletal damage during overwintering and ROS production can be regulated by the dynamics of mitochondrial interconnectivity. Therefore, suppression of ROS formation is likely an important aspect of cytoskeletal arrest. Furthermore, gasotransmitters can regulate ROS levels, as well as cytoskeletal contractility and rearrangement. In this review we will explore the energetic costs of cytoskeletal activity, the cellular mechanisms regulating it, and the potential for cytoskeletal arrest being an important mechanism permitting long-term anoxia survival in anoxia-tolerant species, such as the western painted turtle and goldfish.

Cold acclimation induces life stage-specific responses in the cardiac proteome of western painted turtles (Chrysemys picta bellii): implications for anoxia tolerance

Western painted turtles (Chrysemys picta bellii) are the most anoxia-tolerant tetrapod. Survival time improves at low temperature and during ontogeny, such that adults acclimated to 3°C survive far longer without oxygen than either warm-acclimated adults or cold-acclimated hatchlings. As protein synthesis is rapidly suppressed to save energy at the onset of anoxia exposure, this study tested the hypothesis that cold acclimation would evoke preparatory changes in protein expression to support enhanced anoxia survival in adult but not hatchling turtles. To test this, adult and hatchling turtles were acclimated to either 20°C (warm) or 3°C (cold) for 5 weeks, and then the heart ventricles were collected for quantitative proteomic analysis. The relative abundance of 1316 identified proteins was compared between temperatures and developmental stages. The effect of cold acclimation on the cardiac proteome was only evident in the context of an interaction with life stage, suggesting that ontogenic differences in anoxia tolerance may be predicated on successful maturation of the heart. The main differences between the hatchling and adult cardiac proteomes reflect an increase in metabolic scope with age that included more myoglobin and increased investment in both aerobic and anaerobic energy pathways. Mitochondrial structure and function were key targets of the life stage- and temperature-induced changes to the cardiac proteome, including reduced Complex II proteins in cold-acclimated adults that may help down-regulate the electron transport system and avoid succinate accumulation during anoxia. Therefore, targeted cold-induced changes to the cardiac proteome may be a contributing mechanism for stage-specific anoxia tolerance in turtles.

Brain transcriptomic responses of Yarrow's spiny lizard, Sceloporus jarrovii, to conspecific visual or chemical signals

Species with multimodal communication integrate information from social cues in different modalities into behavioral responses that are mediated by changes in gene expression in the brain. Differences in patterns of gene expression between signal modalities may shed light on the neuromolecular mechanisms underlying multisensory processing. Here, we use RNA-Seq to analyze brain transcriptome responses to either chemical or visual social signals in a territorial lizard with multimodal communication. Using an intruder challenge paradigm, we exposed 18 wild-caught, adult, male Sceloporus jarrovii to either male conspecific scents (femoral gland secretions placed on a small pebble), the species-specific push-up display (a programmed robotic model), or a control (an unscented pebble). We conducted differential expression analysis with both a de novo S. jarrovii transcriptome assembly and the reference genome of a closely related species, Sceloporus undulatus. Despite some inter-individual variation, we found significant differences in gene expression in the brain across signal modalities and the control in both analyses. The most notable differences occurred between chemical and visual stimulus treatments, closely followed by visual stimulus versus the control. Altered expression profiles could explain documented aggression differences in the immediate behavioral response to conspecific signals from different sensory modalities. Shared differentially expressed genes between visually- or chemically-stimulated males are involved in neural activity and neurodevelopment and several other differentially expressed genes in stimulus-challenged males are involved in conserved signal-transduction pathways associated with the social stress response, aggression and the response to territory intruders across vertebrates.

Hepatic transcriptome of the freeze-tolerant Cope's gray treefrog, Dryophytes chrysoscelis: responses to cold acclimation and freezing

Background

Cope’s gray treefrog, Dryophytes chrysoscelis, withstands the physiological challenges of corporeal freezing, partly by accumulating cryoprotective compounds of hepatic origin, including glycerol, urea, and glucose. We hypothesized that expression of genes related to cryoprotectant mobilization and stress tolerance would be differentially regulated in response to cold. Using high-throughput RNA sequencing (RNA-Seq), a hepatic transcriptome was generated for D. chrysoscelis, and gene expression was compared among frogs that were warm-acclimated, cold-acclimated, and frozen.

Results

A total of 159,556 transcripts were generated; 39% showed homology with known transcripts, and 34% of all transcripts were annotated. Gene-level analyses identified 34,936 genes, 85% of which were annotated. Cold acclimation induced differential expression both of genes and non-coding transcripts; freezing induced few additional changes. Transcript-level analysis followed by gene-level aggregation revealed 3582 differentially expressed genes, whereas analysis at the gene level revealed 1324 differentially regulated genes. Approximately 3.6% of differentially expressed sequences were non-coding and of no identifiable homology. Expression of several genes associated with cryoprotectant accumulation was altered during cold acclimation. Of note, glycerol kinase expression decreased with cold exposure, possibly promoting accumulation of glycerol, whereas glucose export was transcriptionally promoted by upregulation of glucose-6-phosphatase and downregulation of genes of various glycolytic enzymes. Several genes related to heat shock protein response, DNA repair, and the ubiquitin proteasome pathway were upregulated in cold and frozen frogs, whereas genes involved in responses to oxidative stress and anoxia, both potential sources of cellular damage during freezing, were downregulated or unchanged.

Conclusion

Our study is the first to report transcriptomic responses to low temperature exposure in a freeze-tolerant vertebrate. The hepatic transcriptome of Dryophytes chrysoscelis is responsive to cold and freezing. Transcriptomic regulation of genes related to particular pathways, such as glycerol biosynthesis, were not all regulated in parallel. The physiological demands associated with cold and freezing, as well as the transcriptomic responses observed in this study, are shared with several organisms that face similar ecophysiological challenges, suggesting common regulatory mechanisms. The role of transcriptional regulation relative to other cellular processes, and of non-coding transcripts as elements of those responses, deserve further study.

Protein lysine methylation in the regulation of anoxia tolerance in the red eared slider turtle, Trachemys scripta elegans

The red eared slider turtle (Trachemys scripta elegans) is a champion vertebrate facultative anaerobe, capable of surviving for several months under conditions of exceptionally low oxygen availability. The ability of the turtle to facilitate this impressive tolerance to oxygen restriction is accomplished through a dramatic reduction in non-essential cellular processes. This is done in an attempt to conserve limited ATP stores and match demand in the anoxic state, with ATP supplied primarily through anaerobic glycolysis. Determining both the non-essential and the essential cellular processes that are deemed to be anoxia-responsive in the turtle has been an intense area of study over the past few decades. As a result, recent advancements have established the influence of global metabolic controls, such as post-transcriptional and post-translational regulation of gene expression in anoxia adaptation. A remaining question is whether or not epigenetic-level regulatory mechanisms are also utilized to allow for local control over gene expression. Recently, research has begun to document lysine methylation as an anoxia-responsive post-translational histone modification, as the activities of a number of methyl-lysine regulatory enzymes are extraordinarily sensitive to oxygen availability. As a result, oxygen-dependent methyl-lysine regulatory enzymes have been of particular interest to several recent studies of animal oxygen sensitivity, including the freshwater turtle. This review will introduce the concept of lysine methylation as an oxygen-sensitive protein modification as well as a prospectus on how this modification may contribute to anoxia tolerance in the turtle.

Navigating oxygen deprivation: liver transcriptomic responses of the red eared slider turtle to environmental anoxia

The best facultative anaerobes among vertebrates are members of the genera Trachemys (pond slider turtles) and Chrysemys (painted turtles), and are able to survive without oxygen for up to 12 to 18 weeks at ∼3 °C. In this study, we utilized RNAseq to profile the transcriptomic changes that take place in response to 20 hrs of anoxia at 5 °C in the liver of the red eared slide turtle (Trachemys scripta elegans). Sequencing reads were obtained from at least 18,169 different genes and represented a minimum 49x coverage of the C. picta bellii exome. A total of 3,105 genes showed statistically significant changes in gene expression between the two animal groups, of which 971 also exhibited a fold change equal to or greater than 50% of control normoxic values. This study also highlights a number of anoxia-responsive molecular pathways that are may be important to navigating anoxia survival. These pathways were enriched in mRNA found to significantly increase in response to anoxia and included molecular processes such as DNA damage repair and metabolic reprogramming. For example, our results indicate that the anoxic turtle may utilize succinate metabolism to yield a molecule of GTP in addition to the two molecules that results from lactate production, and agrees with other established models of anoxia tolerance. Collectively, our analysis provides a snapshot of the molecular landscape of the anoxic turtle and may provide hints into the how this animal is capable of surviving this extreme environmental stress.

Development-specific transcriptomic profiling suggests new mechanisms for anoxic survival in the ventricle of overwintering turtles

Oxygen deprivation swiftly damages tissues in most animals, yet some species show remarkable abilities to tolerate little or even no oxygen. Painted turtles exhibit a development-dependent tolerance that allows adults to survive anoxia ∼4x longer than hatchlings: adults survive ∼170 days and hatchlings survive ∼40 days at 3°C. We hypothesized this difference is related to development-dependent differences in ventricular gene expression. Using a comparative ontogenetic approach, we examined whole transcriptomic changes before, during, and five days after a 20-day bout of anoxic submergence at 3°C. Ontogeny accounted for more gene expression differences than treatment (anoxia or recovery): 1,175 vs. 237 genes, respectively. Of the 237 differences, 93 could confer protection against anoxia and reperfusion injury, 68 could be injurious, and 20 may be constitutively protective. Especially striking during anoxia was the expression pattern of all 76 annotated ribosomal protein (R-protein) mRNAs, which decreased in anoxia-tolerant adults, but increased in anoxia-sensitive hatchlings, suggesting adult-specific regulation of translational suppression. These genes, along with 60 others that decreased their levels in adults and either increased or remained unchanged in hatchlings, implicate antagonistic pleiotropy as a mechanism to resolve the long-standing question about why hatchling painted turtles overwinter in terrestrial nests, rather than emerge and overwinter in water during their first year. In sum, developmental differences in the transcriptome of the turtle ventricle revealed potentially protective mechanisms that contribute to extraordinary adult-specific anoxia tolerance, and provide a unique perspective on differences between the anoxia-induced molecular responses of anoxia-tolerant or anoxia-sensitive phenotypes within a species.

Small Non-coding RNA Expression and Vertebrate Anoxia Tolerance

Background: Extreme anoxia tolerance requires a metabolic depression whose modulation could involve small non-coding RNAs (small ncRNAs), which are specific, rapid, and reversible regulators of gene expression. A previous study of small ncRNA expression in embryos of the annual killifish Austrofundulus limnaeus, the most anoxia-tolerant vertebrate known, revealed a specific expression pattern of small ncRNAs that could play important roles in anoxia tolerance. Here, we conduct a comparative study on the presence and expression of small ncRNAs in the most anoxia-tolerant representatives of several major vertebrate lineages, to investigate the evolution of and mechanisms supporting extreme anoxia tolerance. The epaulette shark (Hemiscyllium ocellatum), crucian carp (Carassius carassius), western painted turtle (Chrysemys picta bellii), and leopard frog (Rana pipiens) were exposed to anoxia and recovery, and small ncRNAs were sequenced from the brain (one of the most anoxia-sensitive tissues) prior to, during, and following exposure to anoxia. Results: Small ncRNA profiles were broadly conserved among species under normoxic conditions, and these expression patterns were largely conserved during exposure to anoxia. In contrast, differentially expressed genes are mostly unique to each species, suggesting that each species may have evolved distinct small ncRNA expression patterns in response to anoxia. Mitochondria-derived small ncRNAs (mitosRNAs) which have a robust response to anoxia in A. limnaeus embryos, were identified in the other anoxia tolerant vertebrates here but did not display a similarly robust response to anoxia. Conclusion: These findings support an overall stabilization of the small ncRNA transcriptome during exposure to anoxic insults, but also suggest that multiple small ncRNA expression pathways may support anoxia tolerance, as no conserved small ncRNA response was identified among the anoxia-tolerant vertebrates studied. This may reflect divergent strategies to achieve the same endpoint: anoxia tolerance. However, it may also indicate that there are multiple cellular pathways that can trigger the same cellular and physiological survival processes, including hypometabolism.

Small noncoding RNA expression during extreme anoxia tolerance of annual killifish (Austrofundulus limnaeus) embryos

Small noncoding RNAs (sncRNA) have recently emerged as specific and rapid regulators of gene expression, involved in a myriad of cellular and organismal processes. MicroRNAs, a class of sncRNAs, are differentially expressed in diverse taxa in response to environmental stress, including anoxia. In most vertebrates, a brief period of oxygen deprivation results in severe tissue damage or death. Studies on sncRNA and anoxia have focused on these anoxia-sensitive species. Studying sncRNAs in anoxia-tolerant organisms may provide insight into adaptive mechanisms supporting anoxia tolerance. Embryos of the annual killifish Austrofundulus limnaeus are the most anoxia-tolerant vertebrates known, surviving over 100 days at their peak tolerance at 25°C. Their anoxia tolerance and physiology vary over development, such that both anoxia-tolerant and anoxia-sensitive phenotypes comprise the species. This allows for a robust comparison to identify sncRNAs essential to anoxia-tolerance. For this study, RNA sequencing was used to identify and quantify expression of sncRNAs in four embryonic stages of A. limnaeus in response to an exposure to anoxia and subsequent aerobic recovery. Unique stage-specific patterns of expression were identified that correlate with anoxia tolerance. In addition, embryos of A. limnaeus appear to constitutively express stress-responsive miRNAs. Most differentially expressed sncRNAs were expressed at higher levels during recovery. Many novel groups of sncRNAs with expression profiles suggesting a key role in anoxia tolerance were identified, including sncRNAs derived from mitochondrial tRNAs. This global analysis has revealed groups of candidate sncRNAs that we hypothesize support anoxia tolerance.

Molecular Physiology of Freeze Tolerance in Vertebrates

Freeze tolerance is an amazing winter survival strategy used by various amphibians and reptiles living in seasonally cold environments. These animals may spend weeks or months with up to ∼65% of their total body water frozen as extracellular ice and no physiological vital signs, and yet after thawing they return to normal life within a few hours. Two main principles of animal freeze tolerance have received much attention: the production of high concentrations of organic osmolytes (glucose, glycerol, urea among amphibians) that protect the intracellular environment, and the control of ice within the body (the first putative ice-binding protein in a frog was recently identified), but many other strategies of biochemical adaptation also contribute to freezing survival. Discussed herein are recent advances in our understanding of amphibian and reptile freeze tolerance with a focus on cell preservation strategies (chaperones, antioxidants, damage defense mechanisms), membrane transporters for water and cryoprotectants, energy metabolism, gene/protein adaptations, and the regulatory control of freeze-responsive hypometabolism at multiple levels (epigenetic regulation of DNA, microRNA action, cell signaling and transcription factor regulation, cell cycle control, and anti-apoptosis). All are providing a much more complete picture of life in the frozen state.

Transcriptomic responses to environmental temperature by turtles with temperature-dependent and genotypic sex determination assessed by RNAseq inform the genetic architecture of embryonic gonadal development

Vertebrate sexual fate is decided primarily by the individual's genotype (GSD), by the environmental temperature during development (TSD), or both. Turtles exhibit TSD and GSD, making them ideal to study the evolution of sex determination. Here we analyze temperature-specific gonadal transcriptomes (RNA-sequencing validated by qPCR) of painted turtles (Chrysemys picta TSD) before and during the thermosensitive period, and at equivalent stages in soft-shell turtles (Apalone spinifera-GSD), to test whether TSD's and GSD's transcriptional circuitry is identical but deployed differently between mechanisms. Our data show that most elements of the mammalian urogenital network are active during turtle gonadogenesis, but their transcription is generally more thermoresponsive in TSD than GSD, and concordant with their sex-specific function in mammals [e.g., upregulation of Amh, Ar, Esr1, Fog2, Gata4, Igf1r, Insr, and Lhx9 at male-producing temperature, and of β-catenin, Foxl2, Aromatase (Cyp19a1), Fst, Nf-kb, Crabp2 at female-producing temperature in Chrysemys]. Notably, antagonistic elements in gonadogenesis (e.g., β-catenin and Insr) were thermosensitive only in TSD early-embryos. Cirbp showed warm-temperature upregulation in both turtles disputing its purported key TSD role. Genes that may convert thermal inputs into sex-specific development (e.g., signaling and hormonal pathways, RNA-binding and heat-shock) were differentially regulated. Jak-Stat, Nf-κB, retinoic-acid, Wnt, and Mapk-signaling (not Akt and Ras-signaling) potentially mediate TSD thermosensitivity. Numerous species-specific ncRNAs (including Xist) were differentially-expressed, mostly upregulated at colder temperatures, as were unannotated loci that constitute novel TSD candidates. Cirbp showed warm-temperature upregulation in both turtles. Consistent transcription between turtles and alligator revealed putatively-critical reptilian TSD elements for male (Sf1, Amh, Amhr2) and female (Crabp2 and Hspb1) gonadogenesis. In conclusion, while preliminary, our data helps illuminate the regulation and evolution of vertebrate sex determination, and contribute genomic resources to guide further research into this fundamental biological process.