The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.

RNA sequencing studies have identified hundreds of non‐coding RNAs in bacteria, including regulatory small RNA (sRNA). However, our understanding of sRNA function has lagged behind their identification due to a lack of tools for the high‐throughput analysis of RNA–RNA interactions in bacteria. Here we demonstrate that in vivo sRNA–mRNA duplexes can be recovered using UV‐crosslinking, ligation and sequencing of hybrids (CLASH). Many sRNAs recruit the endoribonuclease, RNase E, to facilitate processing of mRNAs. We were able to recover base‐paired sRNA–mRNA duplexes in association with RNase E, allowing proximity‐dependent ligation and sequencing of cognate sRNA–mRNA pairs as chimeric reads. We verified that this approach captures bona fide sRNA–mRNA interactions. Clustering analyses identified novel sRNA seed regions and sets of potentially co‐regulated target mRNAs. We identified multiple mRNA targets for the pathotype‐specific sRNA Esr41, which was shown to regulate colicin sensitivity and iron transport in E. coli. Numerous sRNA interactions were also identified with non‐coding RNAs, including sRNAs and tRNAs, demonstrating the high complexity of the sRNA interactome.

Cystic fibrosis (CF) is an inherited disorder where individual disease etiology and response to therapeutic intervention is impacted by CF transmembrane regulator (CFTR) mutations and other genetic modifiers. CFTR regulates multiple mechanisms in a diverse range of epithelial tissues. In this Review, we consolidate the latest updates in the development of primary epithelial cellular model systems relevant for CF. We discuss conventional two-dimensional (2-D) airway epithelial cell cultures, the backbone of in vitro cellular models to date, as well as improved expansion protocols to overcome finite supply of the cellular source. We highlight a range of strategies for establishment of three dimensional (3-D) airway and intestinal organoid models and evaluate the limitations and potential improvements in each system, focusing on their application in CF. The in vitro CFTR functional assays in patient-derived organoids allow for preclinical pharmacotherapy screening to identify responsive patients. It is likely that organoids will be an invaluable preclinical tool to unravel disease mechanisms, design novel treatments, and enable clinicians to provide personalized management for patients with CF.

X chromosome inactivation in eutherian mammals has been thought to be tightly controlled, as expected from a mechanism that compensates for the different dosage of X-borne genes in XX females and XY males. However, many X genes escape inactivation in humans, inactivation of the X in marsupials is partial, and the unrelated sex chromosomes of monotreme mammals have incomplete and gene-specific inactivation of X-linked genes. The bird ZW sex chromosome system represents a third independently evolved amniote sex chromosome system with dosage compensation, albeit partial and gene-specific, via an unknown mechanism (i.e. upregulation of the single Z in females, down regulation of one or both Zs in males, or a combination). We used RNA-fluorescent in situ hybridization (RNA-FISH) to demonstrate, on individual fibroblast cells, inactivation of 11 genes on the chicken Z and 28 genes on the X chromosomes of platypus. Each gene displayed a reproducible frequency of 1Z/1X-active and 2Z/2X-active cells in the homogametic sex. Our results indicate that the probability of inactivation is controlled on a gene-by-gene basis (or small domains) on the chicken Z and platypus X chromosomes. This regulatory mechanism must have been exapted independently to the non-homologous sex chromosomes in birds and mammals in response to an over-expressed Z or X in the homogametic sex, highlighting the universal importance that (at least partial) silencing plays in the evolution on amniote dosage compensation and, therefore, the differentiation of sex chromosomes.

The dinucleotide CpG is highly underrepresented in the genome of human immunodeficiency virus type 1 (HIV-1). To identify the source of CpG depletion in the HIV-1 genome, we investigated two biological mechanisms: (1) CpG methylation-induced transcriptional silencing and (2) CpG recognition by Toll-like receptors (TLRs). We hypothesized that HIV-1 has been under selective evolutionary pressure by these mechanisms leading to the reduction of CpG in its genome. A CpG depleted genome would enable HIV-1 to avoid methylation-induced transcriptional silencing and/or to avoid recognition by TLRs that identify foreign CpG sequences. We investigated these two hypotheses by determining the sequence context dependency of CpG depletion and comparing it with that of CpG methylation and TLR recognition. We found that in both human and HIV-1 genomes the CpG motifs flanked by T/A were depleted most and those flanked by C/G were depleted least. Similarly, our analyses of human methylome data revealed that the CpG motifs flanked by T/A were methylated most and those flanked by C/G were methylated least. Given that a similar CpG depletion pattern was observed for the human genome within which CpGs are not likely to be recognized by TLRs, we argue that the main source of CpG depletion in HIV-1 is likely host-induced methylation. Analyses of CpG motifs in over 100 viruses revealed that this unique CpG representation pattern is specific to the human and simian immunodeficiency viruses.

The marsupial inactive X chromosome expresses a long noncoding RNA (lncRNA) called Rsx that has been proposed to be the functional analog of eutherian Xist. Despite the possibility that Xist and Rsx encode related functions, the two lncRNAs harbor no linear sequence similarity. However, both lncRNAs harbor domains of tandemly repeated sequence. In Xist, these repeat domains are known to be critical for function. Using k-mer based comparison, we show that the repeat domains of Xist and Rsx unexpectedly partition into two major clusters that each harbor substantial levels of nonlinear sequence similarity. Xist Repeats B, C, and D were most similar to each other and to Rsx Repeat 1, whereas Xist Repeats A and E were most similar to each other and to Rsx Repeats 2, 3, and 4. Similarities at the level of k-mers corresponded to domain-specific enrichment of protein-binding motifs. Within individual domains, protein-binding motifs were often enriched to extreme levels. Our data support the hypothesis that Xist and Rsx encode similar functions through different spatial arrangements of functionally analogous protein-binding domains. We propose that the two clusters of repeat domains in Xist and Rsx function in part to cooperatively recruit PRC1 and PRC2 to chromatin. The physical manner in which these domains engage with protein cofactors may be just as critical to the function of the domains as the protein cofactors themselves. The general approaches we outline in this report should prove useful in the study of any set of RNAs.

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. Electronic supplementary material The online version of this article (10.1186/s12864-019-5750-x) contains supplementary material, which is available to authorized users.

DNA methylation plays a key role in maintaining transcriptional silence on the inactive X chromosome of eutherian mammals. Beyond eutherians, there are limited genome wide data on DNA methylation from other vertebrates. Previous studies of X borne genes in various marsupial models revealed no differential DNA methylation of promoters between the sexes, leading to the conclusion that CpG methylation plays no role in marsupial X-inactivation. Using reduced representation bisulfite sequencing, we generated male and female CpG methylation profiles in four representative vertebrates (mouse, gray short-tailed opossum, platypus, and chicken). A variety of DNA methylation patterns were observed. Platypus and chicken displayed no large-scale differential DNA methylation between the sexes on the autosomes or the sex chromosomes. As expected, a metagene analysis revealed hypermethylation at transcription start sites (TSS) of genes subject to X-inactivation in female mice. This contrasted with the opossum, in which metagene analysis did not detect differential DNA methylation between the sexes at TSSs of genes subject to X-inactivation. However, regions flanking TSSs of these genes were hypomethylated. Our data are the first to demonstrate that, for genes subject to X-inactivation in both eutherian and marsupial mammals, there is a consistent difference between DNA methylation levels at TSSs and immediate flanking regions, which we propose has a silencing effect in both groups.