Transcriptome analysis of the response of Burmese python to digestion

Molecular regulation of reversible cardiac remodeling: lessons from species with extreme physiological adaptations

Some vertebrates evolved to have a remarkable capacity for anatomical and physiological plasticity in response to environmental challenges. One example of such plasticity can be found in the ambush-hunting snakes of the genus Python, which exhibit reversible cardiac growth with feeding. The predation strategy employed by pythons is associated with months-long fasts that are arrested by ingestion of large prey. Consequently, digestion compels a dramatic increase in metabolic rate and hypertrophy of multiple organs, including the heart. In this Review, we summarize the post-prandial cardiac adaptations in pythons at the whole-heart, cellular and molecular scales. We highlight circulating factors and cellular signaling pathways that are altered during digestion to affect cardiac form and function and propose possible mechanisms that may drive the post-digestion regression of cardiac mass. Adaptive physiological cardiac hypertrophy has also been observed in other vertebrates, including in fish acclimated to cold water, birds flying at high altitudes and exercising mammals. To reveal potential evolutionarily conserved features, we summarize the molecular signatures of reversible cardiac remodeling identified in these species and compare them with those of pythons. Finally, we offer a perspective on the potential of biomimetics targeting the natural biology of pythons as therapeutics for human heart disease.

Postprandial cardiac hypertrophy is sustained by mechanics, epigenetic, and metabolic reprogramming in pythons

Constricting pythons, known for their ability to consume infrequent, massive meals, exhibit rapid and reversible cardiac hypertrophy following feeding. Our primary goal was to investigate how python hearts achieve this adaptive response after feeding. Isolated myofibrils increased force after feeding without changes in sarcomere ultrastructure and without increasing energy cost. Ca2+ transients were prolonged after feeding with no changes in myofibril Ca2+ sensitivity. Feeding reduced titin-based tension, resulting in decreased cardiac tissue stiffness. Feeding also reduced the activity of sirtuins, a metabolically linked class of histone deacetylases, and increased chromatin accessibility. Transcription factor enrichment analysis on transposase-accessible chromatin with sequencing revealed the prominent role of transcription factors Yin Yang1 and NRF1 in postfeeding cardiac adaptation. Gene expression also changed with the enrichment of translation and metabolism. Finally, metabolomics analysis and adenosine triphosphate production demonstrated that cardiac adaptation after feeding not only increased energy demand but also energy production. These findings have broad implications for our understanding of cardiac adaptation across species and hold promise for the development of innovative approaches to address cardiovascular diseases.

Utility of the burmese Python as a model for studying plasticity of extreme physiological systems

Non-traditional animal models present an opportunity to discover novel biology that has evolved to allow such animals to survive in extreme environments. One striking example is the Burmese python (Python molurus bivittatus), which exhibits extreme physiological adaptation in various metabolic organs after consuming a large meal following long periods of fasting. The response to such a large meal in pythons involves a dramatic surge in metabolic rate, lipid overload in plasma, and massive but reversible organ growth through the course of digestion. Multiple studies have reported the physiological responses in post-prandial pythons, while the specific molecular control of these processes is less well-studied. Investigating the mechanisms that coordinate organ growth and adaptive responses offers the opportunity to gain novel insight that may be able to treat various pathologies in humans. Here, we summarize past research on the post-prandial physiological changes in the Burmese python with a focus on the gastrointestinal tract, heart, and liver. Specifically, we address our recent molecular discoveries in the post-prandial python liver which demonstrate transient adaptations that may reveal new therapeutic targets. Lastly, we explore new biology of the aquaporin 7 gene that is potently upregulated in mammalian cardiac myocytes by circulating factors in post-prandial python plasma.

Animal Models of Exercise From Rodents to Pythons

Acute and chronic animal models of exercise are commonly used in research. Acute exercise testing is used, often in combination with genetic, pharmacological, or other manipulations, to study the impact of these manipulations on the cardiovascular response to exercise and to detect impairments or improvements in cardiovascular function that may not be evident at rest. Chronic exercise conditioning models are used to study the cardiac phenotypic response to regular exercise training and as a platform for discovery of novel pathways mediating cardiovascular benefits conferred by exercise conditioning that could be exploited therapeutically. The cardiovascular benefits of exercise are well established, and, frequently, molecular manipulations that mimic the pathway changes induced by exercise recapitulate at least some of its benefits. This review discusses approaches for assessing cardiovascular function during an acute exercise challenge in rodents, as well as practical and conceptual considerations in the use of common rodent exercise conditioning models. The case for studying feeding in the Burmese python as a model for exercise-like physiological adaptation is also explored.

Convergent evolution of toxin resistance in animals

Convergence is the phenomenon whereby similar phenotypes evolve independently in different lineages. One example is resistance to toxins in animals. Toxins have evolved many times throughout the tree of life. They disrupt molecular and physiological pathways in target species, thereby incapacitating prey or deterring a predator. In response, molecular resistance has evolved in many species exposed to toxins to counteract their harmful effects. Here, we review current knowledge on the convergence of toxin resistance using examples from a wide range of toxin families. We explore the evolutionary processes and molecular adaptations driving toxin resistance. However, resistance adaptations may carry a fitness cost if they disrupt the normal physiology of the resistant animal. Therefore, there is a trade‐off between maintaining a functional molecular target and reducing toxin susceptibility. There are relatively few solutions that satisfy this trade‐off. As a result, we see a small set of molecular adaptations appearing repeatedly in diverse animal lineages, a phenomenon that is consistent with models of deterministic evolution. Convergence may also explain what has been called ‘autoresistance’. This is often thought to have evolved for self‐protection, but we argue instead that it may be a consequence of poisonous animals feeding on toxic prey. Toxin resistance provides a unique and compelling model system for studying the interplay between trophic interactions, selection pressures and the molecular mechanisms underlying evolutionary novelties.

Physiological effects of capture and short-term captivity in an invasive snake species, the Burmese python (Python bivittatus) in Florida

It is important to evaluate the role of captivity as a potential stressor. An understanding of stress responses to capture and transition to captivity may inform the limitations of laboratory studies on wild animals, aid in understanding the consequences of introducing animals into captive environments, and help predict which species may be successful invasives. We investigated physiological effects of captivity by comparing at-capture blood variables in wild Burmese pythons (Python bivittatus) in Florida to pythons recently brought into captivity (1-109 days). We conducted an acute restraint test by collecting samples at baseline (immediately at handling) and one hour post-restraint across wild field-sampled (n = 19) and recently-captive (n = 33) pythons to evaluate fluctuations in plasma corticosterone, bacterial killing ability, antibody response, leukogram, and serpentovirus infection. We observed higher baseline plasma corticosterone and monocytes in recently captive compared to wild snakes, which both subsided in snakes held for a longer time in captivity, and a mild decrease in lymphocytes in the middle of the captivity period. Functional immunity and viral infection were not affected by captivity, and pythons maintained restraint-induced responses in corticosterone, heterophil to lymphocyte ratio, and monocyte counts throughout captivity. Prevalence for serpentovirus was 50%, though infection status was related to sampling date rather than captivity, indicating that viral infection may be seasonal. The history of Burmese python as a common captive animal for research and pet trade, as well as its general resilience to effects of capture and short-term captivity, may contribute to its invasion success in Florida.

Burmese pythons exhibit a transient adaptation to nutrient overload that prevents liver damage

As an opportunistic predator, the Burmese python (Python molurus bivittatus) consumes large and infrequent meals, fasting for up to a year. Upon consuming a large meal, the Burmese python exhibits extreme metabolic responses. To define the pathways that regulate these postprandial metabolic responses, we performed a comprehensive profile of plasma metabolites throughout the digestive process. Following ingestion of a meal equivalent to 25% of its body mass, plasma lipoproteins and metabolites, such as chylomicra and bile acids, reach levels observed only in mammalian models of extreme dyslipidemia. Here, we provide evidence for an adaptive response to postprandial nutrient overload by the python liver, a critical site of metabolic homeostasis. The python liver undergoes a substantial increase in mass through proliferative processes, exhibits hepatic steatosis, hyperlipidemia-induced insulin resistance indicated by PEPCK activation and pAKT deactivation, and de novo fatty acid synthesis via FASN activation. This postprandial state is completely reversible. We posit that Burmese pythons evade the permanent hepatic damage associated with these metabolic states in mammals using evolved protective measures to inactivate these pathways. These include a transient activation of hepatic nuclear receptors induced by fatty acids and bile acids, including PPAR and FXR, respectively. The stress-induced p38 MAPK pathway is also transiently activated during the early stages of digestion. Taken together, these data identify a reversible metabolic response to hyperlipidemia by the python liver, only achieved in mammals by pharmacologic intervention. The factors involved in these processes may be relevant to or leveraged for remediating human hepatic pathology.

Catecholamines are key modulators of ventricular repolarization patterns in the ball python (Python regius)

Ectothermic vertebrates experience daily changes in body temperature, and anecdotal observations suggest these changes affect ventricular repolarization such that the T-wave in the ECG changes polarity. Mammals, in contrast, can maintain stable body temperatures, and their ventricular repolarization is strongly modulated by changes in heart rate and by sympathetic nervous system activity. The aim of this study was to assess the role of body temperature, heart rate, and circulating catecholamines on local repolarization gradients in the ectothermic ball python (Python regius). We recorded body-surface electrocardiograms and performed open-chest high-resolution epicardial mapping while increasing body temperature in five pythons, in all of which there was a change in T-wave polarity. However, the vector of repolarization differed between individuals, and only a subset of leads revealed T-wave polarity change. RNA sequencing revealed regional differences related to adrenergic signaling. In one denervated and Ringer's solution-perfused heart, heating and elevated heart rates did not induce change in T-wave polarity, whereas noradrenaline did. Accordingly, electrocardiograms in eight awake pythons receiving intra-arterial infusion of the β-adrenergic receptor agonists adrenaline and isoproterenol revealed T-wave inversion in most individuals. Conversely, blocking the β-adrenergic receptors using propranolol prevented T-wave change during heating. Our findings indicate that changes in ventricular repolarization in ball pythons are caused by increased tone of the sympathetic nervous system, not by changes in temperature. Therefore, ventricular repolarization in both pythons and mammals is modulated by evolutionary conserved mechanisms involving catecholaminergic stimulation.

Feeding effects on liver mitochondrial bioenergetics of Boa constrictor (Serpentes: Boidae)

Snakes are interesting examples of taxa that can overcome energy metabolism challenges, as many species can endure long periods without feeding, and their eventual meals are of reasonably large sizes, thus exhibiting dual extreme adaptations. Consequently, metabolic rate increases considerably to attend to the energetic demand of digestion, absorption and protein synthesis. These animals should be adapted to transition from these two opposite states of energy fairly quickly, and therefore we investigated mitochondrial function plasticity in these states. Herein, we compared liver mitochondrial bioenergetics of the boid snake Boa constrictor during fasting and after meal intake. We fasted the snakes for 60 days, and then we fed a subgroup with 30% of their body size and evaluated their maximum postprandial response. We measured liver respiration rates from permeabilized tissue and isolated mitochondria. From isolated mitochondria, we also measured Ca2+ retention capacity and redox status. Mitochondrial respiration rates were maximized after feeding, reaching an approximately 60% increase from fasting levels when energized with complex I-linked substrates. Interestingly, fasting and fed snakes exhibited similar respiratory control ratios and citrate synthase activity. Furthermore, we found no differences in Ca2+ retention capacity, indicating no increase in susceptibility to mitochondrial permeability transition, and no changes in mitochondrial redox state, although fed animals exhibited increases in the release of H2O2. Thus, we conclude that liver mitochondria from B. constrictor snakes increase respiration rates during the postprandial period and quickly improve the bioenergetic capacity without compromising redox balance.

The gut–brain axis in vertebrates: implications for food intake regulation

The gut and brain are constantly communicating and influencing each other through neural, endocrine and immune signals in an interaction referred to as the gut–brain axis. Within this communication system, the gastrointestinal tract, including the gut microbiota, sends information on energy status to the brain, which, after integrating these and other inputs, transmits feedback to the gastrointestinal tract. This allows the regulation of food intake and other physiological processes occurring in the gastrointestinal tract, including motility, secretion, digestion and absorption. Although extensive literature is available on the mechanisms governing the communication between the gut and the brain in mammals, studies on this axis in other vertebrates are scarce and often limited to a single species, which may not be representative for obtaining conclusions for an entire group. This Review aims to compile the available information on the gut–brain axis in birds, reptiles, amphibians and fish, with a special focus on its involvement in food intake regulation and, to a lesser extent, in digestive processes. Additionally, we will identify gaps of knowledge that need to be filled in order to better understand the functioning and physiological significance of such an axis in non-mammalian vertebrates.

The vasopressor action of angiotensin II (ANG II) in ball pythons (Python regius)

Angiotensin II (ANG II) is part of the renin-angiotensin system (RAS) in vertebrates and exert vasoconstriction in all species studied. The present study examines the vasopressor effect of ANG II in the ball python (Python regius), and examines whether ANG II exert its effect through direct angiotensin receptors or through an activation of α-adrenergic receptors. The studies were conducted in snakes with chronic arterial catheters that had recovered from anesthesia. In addition to demonstrating a clear and pronounced dose-dependent rise in arterial blood pressure upon repeated injections of boluses with ANG II (0.001-1 μg/kg), we demonstrate that the pressor response persisted following α-adrenergic blockade using the α-adrenergic antagonist phentolamine (2.5 mg/kg). Unfortunately, it proved impossible to block the ANG receptors using losartan (1, 3 or even 10 mg/kg). The pressor response to ANG II was associated with a significant rise in heart rate at the higher dosages, pointing to a resetting of the barostatic mechanism for heart rate regulation. The responses were similar in fasting and digesting pythons despite the expected rise in baseline values for blood pressure and heart rate of the digesting snakes.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

Transcriptome sequencing reveals signatures of positive selection in the Spot-Tailed Earless Lizard

The continual loss of threatened biodiversity is occurring at an accelerated pace. High-throughput sequencing technologies are now providing opportunities to address this issue by aiding in the generation of molecular data for many understudied species of high conservation interest. Our overall goal of this study was to begin building the genomic resources to continue investigations and conservation of the Spot-Tailed Earless lizard. Here we leverage the power of high-throughput sequencing to generate the liver transcriptome for the Northern Spot-Tailed Earless Lizard (Holbrookia lacerata) and Southern Spot-Tailed Earless Lizard (Holbrookia subcaudalis), which have declined in abundance in the past decades, and their sister species, the Common Lesser Earless Lizard (Holbrookia maculata). Our efforts produced high quality and robust transcriptome assemblies validated by 1) quantifying the number of processed reads represented in the transcriptome assembly and 2) quantifying the number of highly conserved single-copy orthologs that are present in our transcript set using the BUSCO pipeline. We found 1,361 1-to-1 orthologs among the three Holbrookia species, Anolis carolinensis, and Sceloporus undulatus. We carried out dN/dS selection tests using a branch-sites model and identified a dozen genes that experienced positive selection in the Holbrookia lineage with functions in development, immunity, and metabolism. Our single-copy orthologous sequences additionally revealed significant pairwise sequence divergence (~.73%) between the Northern H. lacerata and Southern H. subcaudalis that further supports the recent elevation of the Southern Spot-Tailed Earless Lizard to full species.

Effects of feeding and digestion on myocardial contractility and expression of calcium-handling proteins in Burmese pythons (Python molurus)

Pythons are important models of studies on postprandial metabolism because their physiological responses are exacerbated when digesting large prey. Prior studies of these animals have shown hypertrophy of the cardiac tissue 2 to 3 days after feeding, coinciding with the peak of the specific dynamic action (SDA), but the consequences of this remodeling in myocardial contractility have not been studied, which is the purpose of this work. Specimens of Python molurus were divided into two groups: a Digesting group (2 days after feeding, at the peak of SDA), and a Fasting group (28 days after feeding). When compared to the Fasting group, the Digesting group showed higher relative ventricular mass and calcium-handling protein expression such as sarcoplasmic reticulum Ca²⁺-ATPase (SERCA), phospholamban (PLB), and the Na⁺/Ca²⁺ exchanger (NCX). Digesting pythons also exhibited significant increases in the cardiac contraction force (Fc), rates of force development and relaxation, and cardiac pumping capacity. Therefore, the higher SERCA, PLB and NCX expression levels increased cytosolic Ca²⁺ transient amplitude, improving myofilament force. These changes are crucial to maintain cardiac output and a relatively high and continuous blood flow required by metabolic expenditure that occurs in postprandial animals.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Food consumption increases cell proliferation in the python brain

Pythons are model organisms for investigating physiological responses to food intake. While systemic growth in response to food consumption is well documented, what occurs in the brain is currently unexplored. In this study, male ball pythons (Python regius) were used to test the hypothesis that food consumption stimulates cell proliferation in the brain. We used 5-bromo-12'-deoxyuridine (BrdU) as a cell-birth marker to quantify and compare cell proliferation in the brain of fasted snakes and those at 2 and 6 days after a meal. Throughout the telencephalon, cell proliferation was significantly increased in the 6 day group, with no difference between the 2 day group and controls. Systemic postprandial plasticity occurs quickly after a meal is ingested, during the period of active digestion; however, the brain displays a surge of cell proliferation after most digestion and absorption is complete.

Transcriptome analysis of the response of Burmese python to digestion

Exceptional and extreme feeding behaviour makes the Burmese python (Python bivittatus) an interesting model to study physiological remodelling and metabolic adaptation in response to refeeding after prolonged starvation. In this study, we used transcriptome sequencing of 5 visceral organs during fasting as well as 24 hours and 48 hours after ingestion of a large meal to unravel the postprandial changes in Burmese pythons. We first used the pooled data to perform a de novo assembly of the transcriptome and supplemented this with a proteomic survey of enzymes in the plasma and gastric fluid. We constructed a high-quality transcriptome with 34 423 transcripts, of which 19 713 (57%) were annotated. Among highly expressed genes (fragments per kilo base per million sequenced reads > 100 in 1 tissue), we found that the transition from fasting to digestion was associated with differential expression of 43 genes in the heart, 206 genes in the liver, 114 genes in the stomach, 89 genes in the pancreas, and 158 genes in the intestine. We interrogated the function of these genes to test previous hypotheses on the response to feeding. We also used the transcriptome to identify 314 secreted proteins in the gastric fluid of the python. Digestion was associated with an upregulation of genes related to metabolic processes, and translational changes therefore appear to support the postprandial rise in metabolism. We identify stomach-related proteins from a digesting individual and demonstrate that the sensitivity of modern liquid chromatography/tandem mass spectrometry equipment allows the identification of gastric juice proteins that are present during digestion.