Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia.

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

State-transition modeling of blood transcriptome predicts disease evolution and treatment response in chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention, before phenotypic changes become detectable.

State-transition Modeling of Blood Transcriptome Predicts Disease Evolution and Treatment Response in Chronic Myeloid Leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention even before phenotypic changes become detectable.

Integration of single-cell transcriptomes and biological function reveals distinct behavioral patterns in bone marrow endothelium.

Heterogeneity of endothelial cell (EC) populations reflects their diverse functions in maintaining tissue's homeostasis. However, their phenotypic, molecular, and functional properties are not entirely mapped. We use the Tie2-CreERT2;Rosa26-tdTomato reporter mouse to trace, profile, and cultivate primary ECs from different organs. As paradigm platform, we use this strategy to study bone marrow endothelial cells (BMECs). Single-cell mRNA sequencing of primary BMECs reveals that their diversity and native molecular signatures is transitorily preserved in an ex vivo culture that conserves key cell-to-cell microenvironment interactions. Macrophages sustain BMEC cellular diversity and expansion and preserve sinusoidal-like BMECs ex vivo. Endomucin expression discriminates BMECs in populations exhibiting mutually exclusive properties and distinct sinusoidal/arterial and tip/stalk signatures. In contrast to arterial-like, sinusoidal-like BMECs are short-lived, form 2D-networks, contribute to in vivo angiogenesis, and support hematopoietic stem/progenitor cells in vitro. This platform can be extended to other organs' ECs to decode mechanistic information and explore therapeutics.

Dynamic patterns of microRNA expression during acute myeloid leukemia state-transition.

MicroRNAs (miRNAs) have been shown to hold prognostic value in acute myeloid leukemia (AML); however, the temporal dynamics of miRNA expression in AML are poorly understood. Using serial samples from a mouse model of AML to generate time-series miRNA sequencing data, we are the first to show that the miRNA transcriptome undergoes state-transition during AML initiation and progression. We modeled AML state-transition as a particle undergoing Brownian motion in a quasi-potential and validated the AML state-space and state-transition model to accurately predict time to AML in an independent cohort of mice. The critical points of the model provided a framework to align samples from mice that developed AML at different rates. Our mathematical approach allowed discovery of dynamic processes involved during AML development and, if translated to humans, has the potential to predict an individual's disease trajectory.

Pharmacometabonomic association of cyclophosphamide 4-hydroxylation in hematopoietic cell transplant recipients.

The widely used alkylating agent cyclophosphamide (CY) has substantive interpatient variability in the area under the curve (AUC) of it and its metabolites. Numerous factors may influence the drug-metabolizing enzymes that metabolize CY to 4-hydroxycyclophosphamide (4HCY), the principal precursor to CY's cytotoxic metabolite. We sought to identify endogenous metabolomics compounds (EMCs) associated with 4HCY formation clearance (ratio of 4HCY/CY AUC) using global metabolomics. Patients who undergo hematopoietic cell transplantation receiving post-transplant CY (PT-CY) were enrolled, cohort 1 (n = 26) and cohort 2 (n = 25) donating longitudinal blood samples before they started HCT (pre-HCT), before infusion of the donor allograft (pre-graft), before the first dose of PT-CY (pre-CY), and 24 h after the first dose of PT-CY (24-h post-CY), which is also immediately before the second dose of CY. A total of 512 and 498 EMCs were quantitated in two cohorts, respectively. Both univariate linear regression with false discovery rate (FDR), and pathway enrichment analyses using a global association test were performed. At the pre-CY time point, no EMCs were associated at FDR less than 0.1. At pre-HCT, cohort 1 had one EMC (levoglucosan) survive the FDR threshold. At pre-graft, cohort 1 and cohort 2 had 20 and 13 EMCs, respectively, exhibiting unadjusted p values less than 0.05, with the only EMCs having an FDR less than 0.1 being two unknown EMCs. At 24-h post-CY, there were three EMCs, two ketones, and threitol, at FDR less than 0.1 in cohort 2. These results demonstrate the potential of pharmacometabonomics, but future studies in larger samples are needed to optimize CY.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation.

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

A SARS-CoV-2 targeted siRNA-nanoparticle therapy for COVID-19.

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in humans. Despite several emerging vaccines, there remains no verifiable therapeutic targeted specifically to the virus. Here we present a highly effective small interfering RNA (siRNA) therapeutic against SARS-CoV-2 infection using a novel lipid nanoparticle (LNP) delivery system. Multiple siRNAs targeting highly conserved regions of the SARS-CoV-2 virus were screened, and three candidate siRNAs emerged that effectively inhibit the virus by greater than 90% either alone or in combination with one another. We simultaneously developed and screened two novel LNP formulations for the delivery of these candidate siRNA therapeutics to the lungs, an organ that incurs immense damage during SARS-CoV-2 infection. Encapsulation of siRNAs in these LNPs followed by in vivo injection demonstrated robust repression of virus in the lungs and a pronounced survival advantage to the treated mice. Our LNP-siRNA approaches are scalable and can be administered upon the first sign of SARS-CoV-2 infection in humans. We suggest that an siRNA-LNP therapeutic approach could prove highly useful in treating COVID-19 disease as an adjunctive therapy to current vaccine strategies.

Koala cathelicidin PhciCath5 has antimicrobial activity, including against Chlamydia pecorum.

Devastating fires in Australia over 2019-20 decimated native fauna and flora, including koalas. The resulting population bottleneck, combined with significant loss of habitat, increases the vulnerability of remaining koala populations to threats which include disease. Chlamydia is one disease which causes significant morbidity and mortality in koalas. The predominant pathogenic species, Chlamydia pecorum, causes severe ocular, urogenital and reproductive tract disease. In marsupials, including the koala, gene expansions of an antimicrobial peptide family known as cathelicidins have enabled protection of immunologically naïve pouch young during early development. We propose that koala cathelicidins are active against Chlamydia and other bacteria and fungi. Here we describe ten koala cathelicidins, five of which contained full length coding sequences that were widely expressed in tissues throughout the body. Focusing on these five, we investigate their antimicrobial activity against two koala C. pecorum isolates from distinct serovars; MarsBar and IPTaLE, as well as other bacteria and fungi. One cathelicidin, PhciCath5, inactivated C. pecorum IPTaLE and MarsBar elementary bodies and significantly reduced the number of inclusions compared to the control (p<0.0001). Despite evidence of cathelicidin expression within tissues known to be infected by Chlamydia, natural PhciCath5 concentrations may be inadequate in vivo to prevent or control C. pecorum infections in koalas. PhciCath5 also displayed antimicrobial activity against fungi and Gram negative and positive bacteria, including methicillin-resistant Staphylococcus aureus (MRSA). Electrostatic interactions likely drive PhciCath5 adherence to the pathogen cell membrane, followed by membrane permeabilisation leading to cell death. Activity against E. coli was reduced in the presence of 10% serum and 20% whole blood. Future modification of the PhciCath5 peptide to enhance activity, including in the presence of serum/blood, may provide a novel solution to Chlamydia infection in koalas and other species.

Inhibition of CD38 and supplementation of nicotinamide riboside ameliorate lipopolysaccharide-induced microglial and astrocytic neuroinflammation by increasing NAD.

Neuroinflammation is initiated by activation of the brain's innate immune system in response to an inflammatory challenge. Insufficient control of neuroinflammation leads to enhanced or prolonged pathology in various neurological conditions including multiple sclerosis and Alzheimer's disease. Nicotinamide adenine dinucleotide (NAD+ ) plays critical roles in cellular energy metabolism and calcium homeostasis. Our previous study demonstrated that deletion of CD38, which consumes NAD+ , suppressed cuprizone-induced demyelination, neuroinflammation, and glial activation. However, it is still unknown whether CD38 directly affects neuroinflammation through regulating brain NAD+ level. In this study, we investigated the effect of CD38 deletion and inhibition and supplementation of NAD+ on lipopolysaccharide (LPS)-induced neuroinflammation in mice. Intracerebroventricular injection of LPS significantly increased CD38 expression especially in the hippocampus. Deletion of CD38 decreased LPS-induced inflammatory responses and glial activation. Pre-administration of apigenin, a flavonoid with CD38 inhibitory activity, or nicotinamide riboside (NR), an NAD+ precursor, increased NAD+ level, and significantly suppressed induction of cytokines and chemokines, glial activation and subsequent neurodegeneration after LPS administration. In cell culture, LPS-induced inflammatory responses were suppressed by treatment of primary astrocytes or microglia with apigenin, NAD+ , NR or 78c, the latter a specific CD38 inhibitor. Finally, all these compounds suppressed NF-κB signaling pathway in microglia. These results suggest that CD38-mediated neuroinflammation is linked to NAD+ consumption and that boosting NAD+ by CD38 inhibition and NR supplementation directly suppress neuroinflammation in the brain.

A SARS-CoV-2 targeted siRNA-nanoparticle therapy for COVID-19.

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in humans. Despite several emerging vaccines, there remains no verifiable therapeutic targeted specifically to the virus. Here we present a highly effective siRNA therapeutic against SARS-CoV-2 infection using a novel lipid nanoparticle delivery system. Multiple small-interfering RNAs (siRNAs) targeting highly conserved regions of the SARS-CoV-2 virus were screened and three candidate siRNAs emerged that effectively inhibit virus by greater than 90% either alone or in combination with one another. We simultaneously developed and screened two novel lipid nanoparticle formulations for the delivery of these candidate siRNA therapeutics to the lungs, an organ that incurs immense damage during SARS-CoV-2 infection. Encapsulation of siRNAs in these LNPs followed by in vivo injection demonstrated robust repression of virus in the lungs and a pronounced survival advantage to the treated mice. Our LNP-siRNA approaches are scalable and can be administered upon the first sign of SARS-CoV-2 infection in humans. We suggest that an siRNA-LNP therapeutic approach could prove highly useful in treating COVID-19 disease as an adjunctive therapy to current vaccine strategies.

RAMP2-AS1 Regulates Endothelial Homeostasis and Aging.

The homeostasis of vascular endothelium is crucial for cardiovascular health and endothelial cell (EC) aging and dysfunction could negatively impact vascular function. Leveraging transcriptome profiles from ECs subjected to various stimuli, including time-series data obtained from ECs under physiological pulsatile flow vs. pathophysiological oscillatory flow, we performed principal component analysis (PCA) to identify key genes contributing to divergent transcriptional states of ECs. Through bioinformatics analysis, we identified that a long non-coding RNA (lncRNA) RAMP2-AS1 encoded on the antisense of RAMP2, a determinant of endothelial homeostasis and vascular integrity, is a novel regulator essential for EC homeostasis and function. Knockdown of RAMP2-AS1 suppressed RAMP2 expression and caused EC functional changes promoting aging, including impaired angiogenesis and increased senescence. Our study demonstrates an integrative approach to quantifying EC aging based on transcriptome changes, which also identified a number of novel regulators, including protein-coding genes and many lncRNAs involved EC functional modulation, exemplified by RAMP2-AS1.

Broadly active zinc finger protein-guided transcriptional activation of HIV-1.

Human immunodeficiency virus type 1 (HIV-1) causes a persistent viral infection resulting in the demise of immune regulatory cells. Clearance of HIV-1 infection results in integration of proviral DNA into the genome of host cells, which provides a means for evasion and long-term persistence. A therapeutic compound that specifically targets and sustainably activates a latent HIV-1 provirus could be transformative and is the goal for the "shock-and-kill" approach to a functional cure for HIV-1. Substantial progress has been made toward the development of recombinant proteins that target specific genomic loci for gene activation, repression, or inactivation by directed mutations. However, most of these modalities are too large or too complex for efficient therapeutic application. We describe here the development and testing of a novel recombinant zinc finger protein transactivator, ZFP-362-VPR, which specifically and potently enhances proviral HIV-1 transcription both in established latency models and activity across different viral clades. Additionally, ZFP-362-VPR-activated HIV-1 reporter gene expression in a well-established primary human CD4+ T cell latency model and off-target pathways were determined by transcriptome analyses. This study provides clear proof of concept for the application of a novel, therapeutically relevant, protein transactivator to purge cellular reservoirs of HIV-1.

MicroRNA dynamics during hibernation of the Australian central bearded dragon (Pogona vitticeps).

Hibernation is a physiological state employed by many animals that are exposed to limited food and adverse winter conditions. Controlling tissue-specific and organism wide changes in metabolism and cellular function requires precise regulation of gene expression, including by microRNAs (miRNAs). Here we profile miRNA expression in the central bearded dragon (Pogona vitticeps) using small RNA sequencing of brain, heart, and skeletal muscle from individuals in late hibernation and four days post-arousal. A total of 1295 miRNAs were identified in the central bearded dragon genome; 664 of which were novel to central bearded dragon. We identified differentially expressed miRNAs (DEmiRs) in all tissues and correlated mRNA expression with known and predicted target mRNAs. Functional analysis of DEmiR targets revealed an enrichment of differentially expressed mRNA targets involved in metabolic processes. However, we failed to reveal biologically relevant tissue-specific processes subjected to miRNA-mediated regulation in heart and skeletal muscle. In brain, neuroprotective pathways were identified as potential targets regulated by miRNAs. Our data suggests that miRNAs are necessary for modulating the shift in cellular metabolism during hibernation and regulating neuroprotection in the brain. This study is the first of its kind in a hibernating reptile and provides key insight into this ephemeral phenotype.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction.

BACKGROUND: Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes. RESULTS: We generated 96.7 × 109 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 109 paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females. CONCLUSIONS: The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

Author Correction: Transcriptomic changes in the pre-implantation uterus highlight histotrophic nutrition of the developing marsupial embryo.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

State-Transition Analysis of Time-Sequential Gene Expression Identifies Critical Points That Predict Development of Acute Myeloid Leukemia.

Temporal dynamics of gene expression inform cellular and molecular perturbations associated with disease development and evolution. Given the complexity of high-dimensional temporal genomic data, an analytic framework guided by a robust theory is needed to interpret time-sequential changes and to predict system dynamics. Here we model temporal dynamics of the transcriptome of peripheral blood mononuclear cells in a two-dimensional state-space representing states of health and leukemia using time-sequential bulk RNA-seq data from a murine model of acute myeloid leukemia (AML). The state-transition model identified critical points that accurately predict AML development and identifies stepwise transcriptomic perturbations that drive leukemia progression. The geometry of the transcriptome state-space provided a biological interpretation of gene dynamics, aligned gene signals that are not synchronized in time across mice, and allowed quantification of gene and pathway contributions to leukemia development. Our state-transition model synthesizes information from multiple cell types in the peripheral blood and identifies critical points in the transition from health to leukemia to guide interpretation of changes in the transcriptome as a whole to predict disease progression. SIGNIFICANCE: These findings apply the theory of state transitions to model the initiation and development of acute myeloid leukemia, identifying transcriptomic perturbations that accurately predict time to disease development.See related commentary by Kuijjer, p. 3072 GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/15/3157/F1.large.jpg.

ZW Sex Chromosomes in Australian Dragon Lizards (Agamidae) Originated from a Combination of Duplication and Translocation in the Nucleolar Organising Region.

Sex chromosomes in some reptiles share synteny with distantly related amniotes in regions orthologous to squamate chromosome 2. The latter finding suggests that chromosome 2 was formerly part of a larger ancestral (amniote) super-sex chromosome and raises questions about how sex chromosomes are formed and modified in reptiles. Australian dragon lizards (Agamidae) are emerging as an excellent model for studying these processes. In particular, they exhibit both genotypic (GSD) and temperature-dependent (TSD) sex determination, show evidence of transitions between the two modes and have evolved non-homologous ZW sex microchromosomes even within the same evolutionary lineage. They therefore represent an excellent group to probe further the idea of a shared ancestral super-sex chromosome and to investigate mechanisms for transition between different sex chromosome forms. Here, we compare sex chromosome homology among eight dragon lizard species from five genera to identify key cytological differences and the mechanisms that may be driving sex chromosome evolution in this group. We performed fluorescence in situ hybridisation to physically map bacterial artificial chromosome (BAC) clones from the bearded dragon, Pogona vitticeps' ZW sex chromosomes and a nucleolar organising region (NOR) probe in males and females of eight Agamid species exhibiting either GSD or TSD. We show that the sex chromosome derived BAC clone hybridises near the telomere of chromosome 2q in all eight species examined. This clone also hybridises to the sex microchromosomes of three species (P vitticeps, P. barbata and Diporiphora nobbi) and a pair of microchromosomes in three others (Ctenophorus pictus, Amphibolurus norrisi and Amphibolurus muricatus). No other chromosomes are marked by the probe in two species from the closely related genus Physignathus. A probe bearing nucleolar organising region (NOR) sequences maps close to the telomere of chromosome 2q in all eight species, and to the ZW pair in P. vitticeps and P. barbata, the W microchromosome in D. nobbi, and several microchromosomes in P. cocincinus. Our findings provide evidence of sequence homology between chromosome 2 and the sex chromosomes of multiple agamids. These data support the hypothesis that there was an ancestral sex chromosome in amniotes that gave rise to squamate chromosome 2 and raises the prospect that some particular property of this chromosome has favoured its role as a sex chromosome in amniotes. It is likely that the amplification of repetitive sequences associated with this region has driven the high level of heterochromatinisation of the sex-specific chromosomes in three species of agamid. Our data suggest a possible mechanism for chromosome rearrangement, including inversion and duplication near the telomeric regions of the ancestral chromosome 2 and subsequent translocation to the ZW sex microchromosomes in three agamid species. It is plausible that these chromosome rearrangements involving sex chromosomes also drove speciation in this group.

Waking the sleeping dragon: gene expression profiling reveals adaptive strategies of the hibernating reptile Pogona vitticeps.

BACKGROUND: Hibernation is a physiological state exploited by many animals exposed to prolonged adverse environmental conditions associated with winter. Large changes in metabolism and cellular function occur, with many stress response pathways modulated to tolerate physiological challenges that might otherwise be lethal. Many studies have sought to elucidate the molecular mechanisms of mammalian hibernation, but detailed analyses are lacking in reptiles. Here we examine gene expression in the Australian central bearded dragon (Pogona vitticeps) using mRNA-seq and label-free quantitative mass spectrometry in matched brain, heart and skeletal muscle samples from animals at late hibernation, 2 days post-arousal and 2 months post-arousal. RESULTS: We identified differentially expressed genes in all tissues between hibernation and post-arousal time points; with 4264 differentially expressed genes in brain, 5340 differentially expressed genes in heart, and 5587 differentially expressed genes in skeletal muscle. Furthermore, we identified 2482 differentially expressed genes across all tissues. Proteomic analysis identified 743 proteins (58 differentially expressed) in brain, 535 (57 differentially expressed) in heart, and 337 (36 differentially expressed) in skeletal muscle. Tissue-specific analyses revealed enrichment of protective mechanisms in all tissues, including neuroprotective pathways in brain, cardiac hypertrophic processes in heart, and atrophy protective pathways in skeletal muscle. In all tissues stress response pathways were induced during hibernation, as well as evidence for gene expression regulation at transcription, translation and post-translation. CONCLUSIONS: These results reveal critical stress response pathways and protective mechanisms that allow for maintenance of both tissue-specific function, and survival during hibernation in the central bearded dragon. Furthermore, we provide evidence for multiple levels of gene expression regulation during hibernation, particularly enrichment of miRNA-mediated translational repression machinery; a process that would allow for rapid and energy efficient reactivation of translation from mature mRNA molecules at arousal. This study is the first molecular investigation of its kind in a hibernating reptile, and identifies strategies not yet observed in other hibernators to cope stress associated with this remarkable state of metabolic depression.

Draft genome assembly of the invasive cane toad, Rhinella marina.

Background: The cane toad (Rhinella marina formerly Bufo marinus) is a species native to Central and South America that has spread across many regions of the globe. Cane toads are known for their rapid adaptation and deleterious impacts on native fauna in invaded regions. However, despite an iconic status, there are major gaps in our understanding of cane toad genetics. The availability of a genome would help to close these gaps and accelerate cane toad research. Findings: We report a draft genome assembly for R. marina, the first of its kind for the Bufonidae family. We used a combination of long-read Pacific Biosciences RS II and short-read Illumina HiSeq X sequencing to generate 359.5 Gb of raw sequence data. The final hybrid assembly of 31,392 scaffolds was 2.55 Gb in length with a scaffold N50 of 168 kb. BUSCO analysis revealed that the assembly included full length or partial fragments of 90.6% of tetrapod universal single-copy orthologs (n = 3950), illustrating that the gene-containing regions have been well assembled. Annotation predicted 25,846 protein coding genes with similarity to known proteins in Swiss-Prot. Repeat sequences were estimated to account for 63.9% of the assembly. Conclusions: The R. marina draft genome assembly will be an invaluable resource that can be used to further probe the biology of this invasive species. Future analysis of the genome will provide insights into cane toad evolution and enrich our understanding of their interplay with the ecosystem at large.