Deep Phenotyping of Sleep in Drosophila.

Sleep is an evolutionarily conserved behavior, whose function is unknown. Here, we present a method for deep phenotyping of sleep in Drosophila, consisting of a high-resolution video imaging system, coupled with closed-loop laser perturbation to measure arousal threshold. To quantify sleep-associated microbehaviors, we trained a deep-learning network to annotate body parts in freely moving flies and developed a semi-supervised computational pipeline to classify behaviors. Quiescent flies exhibit a rich repertoire of microbehaviors, including proboscis pumping (PP) and haltere switches, which vary dynamically across the night. Using this system, we characterized the effects of optogenetically activating two putative sleep circuits. These data reveal that activating dFB neurons produces micromovements, inconsistent with sleep, while activating R5 neurons triggers PP followed by behavioral quiescence. Our findings suggest that sleep in Drosophila is polyphasic with different stages and set the stage for a rigorous analysis of sleep and other behaviors in this species.

Secondary bile acids function through the vitamin D receptor in myeloid progenitors to promote myelopoiesis.

Metabolic products of the microbiota can alter hematopoiesis. However, the contribution and site of action of bile acids is poorly understood. Here, we demonstrate that the secondary bile acids, deoxycholic acid (DCA) and lithocholic acid (LCA), increase bone marrow myelopoiesis. Treatment of bone marrow cells with DCA and LCA preferentially expanded immunophenotypic and functional colony-forming unit-granulocyte and macrophage (CFU-GM) granulocyte-monocyte progenitors (GMPs). DCA treatment of sorted hematopoietic stem and progenitor cells (HSPCs) increased CFU-GMs, indicating that direct exposure of HSPCs to DCA sufficed to increase GMPs. The vitamin D receptor (VDR) was required for the DCA-induced increase in CFU-GMs and GMPs. Single-cell RNA sequencing revealed that DCA significantly upregulated genes associated with myeloid differentiation and proliferation in GMPs. The action of DCA on HSPCs to expand GMPs in a VDR-dependent manner suggests microbiome-host interactions could directly affect bone marrow hematopoiesis and potentially the severity of infectious and inflammatory disease.

Pathogenic human variant that dislocates GATA2 zinc fingers disrupts hematopoietic gene expression and signaling networks.

Although certain human genetic variants are conspicuously loss of function, decoding the impact of many variants is challenging. Previously, we described a patient with leukemia predisposition syndrome (GATA2 deficiency) with a germline GATA2 variant that inserts 9 amino acids between the 2 zinc fingers (9aa-Ins). Here, we conducted mechanistic analyses using genomic technologies and a genetic rescue system with Gata2 enhancer-mutant hematopoietic progenitor cells to compare how GATA2 and 9aa-Ins function genome-wide. Despite nuclear localization, 9aa-Ins was severely defective in occupying and remodeling chromatin and regulating transcription. Variation of the inter-zinc finger spacer length revealed that insertions were more deleterious to activation than repression. GATA2 deficiency generated a lineage-diverting gene expression program and a hematopoiesis-disrupting signaling network in progenitors with reduced granulocyte-macrophage colony-stimulating factor (GM-CSF) and elevated IL-6 signaling. As insufficient GM-CSF signaling caused pulmonary alveolar proteinosis and excessive IL-6 signaling promoted bone marrow failure and GATA2 deficiency patient phenotypes, these results provide insight into mechanisms underlying GATA2-linked pathologies.

PKM2-TMEM33 axis regulates lipid homeostasis in cancer cells by controlling SCAP stability.

The pyruvate kinase M2 isoform (PKM2) is preferentially expressed in cancer cells to regulate anabolic metabolism. Although PKM2 was recently reported to regulate lipid homeostasis, the molecular mechanism remains unclear. Herein, we discovered an ER transmembrane protein 33 (TMEM33) as a downstream effector of PKM2 that regulates activation of SREBPs and lipid metabolism. Loss of PKM2 leads to up-regulation of TMEM33, which recruits RNF5, an E3 ligase, to promote SREBP-cleavage activating protein (SCAP) degradation. TMEM33 is transcriptionally regulated by nuclear factor erythroid 2-like 1 (NRF1), whose cleavage and activation are controlled by PKM2 levels. Total plasma cholesterol levels are elevated by either treatment with PKM2 tetramer-promoting agent TEPP-46 or by global PKM2 knockout in mice, highlighting the essential function of PKM2 in lipid metabolism. Although depletion of PKM2 decreases cancer cell growth, global PKM2 knockout accelerates allografted tumor growth. Together, our findings reveal the cell-autonomous and systemic effects of PKM2 in lipid homeostasis and carcinogenesis, as well as TMEM33 as a bona fide regulator of lipid metabolism.

MLG: multilayer graph clustering for multi-condition scRNA-seq data.

Single-cell transcriptome sequencing (scRNA-seq) enabled investigations of cellular heterogeneity at exceedingly higher resolutions. Identification of novel cell types or transient developmental stages across multiple experimental conditions is one of its key applications. Linear and non-linear dimensionality reduction for data integration became a foundational tool in inference from scRNA-seq data. We present multilayer graph clustering (MLG) as an integrative approach for combining multiple dimensionality reduction of multi-condition scRNA-seq data. MLG generates a multilayer shared nearest neighbor cell graph with higher signal-to-noise ratio and outperforms current best practices in terms of clustering accuracy across large-scale benchmarking experiments. Application of MLG to a wide variety of datasets from multiple conditions highlights how MLG boosts signal-to-noise ratio for fine-grained sub-population identification. MLG is widely applicable to settings with single cell data integration via dimension reduction.

The SWI/SNF-Related, Matrix Associated, Actin-Dependent Regulator of Chromatin A4 Core Complex Represses Respiratory Syncytial Virus-Induced Syncytia Formation and Subepithelial Myofibroblast Transition.

Epigenetics plays an important role in the priming the dynamic response of airway epithelial cells to infectious and environmental stressors. Here, we examine the epigenetic role of the SWI/SNF Related, Matrix Associated, Actin Dependent Regulator of Chromatin A4 (SMARCA4) in the epithelial response to RSV infection. Depletion of SMARCA4 destabilized the abundance of the SMARCE1/ARID1A SWI/SNF subunits, disrupting the innate response and triggering a hybrid epithelial/mesenchymal (E/M) state. Assaying SMARCA4 complex-regulated open chromatin domains by transposase cleavage -next generation sequencing (ATAC-Seq), we observed that the majority of cleavage sites in uninfected cells have reduced chromatin accessibility. Paradoxically, SMARCA4 complex-depleted cells showed enhanced RSV-inducible chromatin opening and gene expression in the EMT pathway genes, MMP9, SNAI1/2, VIM, and CDH2. Focusing on the key MMP9, we observed that SMARCA4 complex depletion reduced basal BRD4 and RNA Polymerase II binding, but enhanced BRD4/Pol II binding in response to RSV infection. In addition, we observed that MMP9 secretion in SMARCA4 complex deficient cells contributes to mesenchymal transition, cellular fusion (syncytia) and subepithelial myofibroblast transition. We conclude the SMARCA4 complex is a transcriptional repressor of epithelial plasticity, whose depletion triggers a hybrid E/M state that affects the dynamic response of the small airway epithelial cell in mucosal remodeling via paracrine MMP9 activity.

GATA factor-regulated solute carrier ensemble reveals a nucleoside transporter-dependent differentiation mechanism.

Developmental-regulatory networks often include large gene families encoding mechanistically-related proteins like G-protein-coupled receptors, zinc finger transcription factors and solute carrier (SLC) transporters. In principle, a common mechanism may confer expression of multiple members integral to a developmental process, or diverse mechanisms may be deployed. Using genetic complementation and enhancer-mutant systems, we analyzed the 456 member SLC family that establishes the small molecule constitution of cells. This analysis identified SLC gene cohorts regulated by GATA1 and/or GATA2 during erythroid differentiation. As >50 SLC genes shared GATA factor regulation, a common mechanism established multiple members of this family. These genes included Slc29a1 encoding an equilibrative nucleoside transporter (Slc29a1/ENT1) that utilizes adenosine as a preferred substrate. Slc29a1 promoted erythroblast survival and differentiation ex vivo. Targeted ablation of murine Slc29a1 in erythroblasts attenuated erythropoiesis and erythrocyte regeneration in response to acute anemia. Our results reveal a GATA factor-regulated SLC ensemble, with a nucleoside transporter component that promotes erythropoiesis and prevents anemia, and establish a mechanistic link between GATA factor and adenosine mechanisms. We propose that integration of the GATA factor-adenosine circuit with other components of the GATA factor-regulated SLC ensemble establishes the small molecule repertoire required for progenitor cells to efficiently generate erythrocytes.

PRAM: a novel pooling approach for discovering intergenic transcripts from large-scale RNA sequencing experiments.

Publicly available RNA-seq data is routinely used for retrospective analysis to elucidate new biology. Novel transcript discovery enabled by joint analysis of large collections of RNA-seq data sets has emerged as one such analysis. Current methods for transcript discovery rely on a '2-Step' approach where the first step encompasses building transcripts from individual data sets, followed by the second step that merges predicted transcripts across data sets. To increase the power of transcript discovery from large collections of RNA-seq data sets, we developed a novel '1-Step' approach named Pooling RNA-seq and Assembling Models (PRAM) that builds transcript models from pooled RNA-seq data sets. We demonstrate in a computational benchmark that 1-Step outperforms 2-Step approaches in predicting overall transcript structures and individual splice junctions, while performing competitively in detecting exonic nucleotides. Applying PRAM to 30 human ENCODE RNA-seq data sets identified unannotated transcripts with epigenetic and RAMPAGE signatures similar to those of recently annotated transcripts. In a case study, we discovered and experimentally validated new transcripts through the application of PRAM to mouse hematopoietic RNA-seq data sets. We uncovered new transcripts that share a differential expression pattern with a neighboring gene Pik3cg implicated in human hematopoietic phenotypes, and we provided evidence for the conservation of this relationship in human. PRAM is implemented as an R/Bioconductor package.

Constructing and deconstructing GATA2-regulated cell fate programs to establish developmental trajectories.

Stem and progenitor cell fate transitions constitute key decision points in organismal development that enable access to a developmental path or actively preclude others. Using the hematopoietic system, we analyzed the relative importance of cell fate-promoting mechanisms versus negating fate-suppressing mechanisms to engineer progenitor cells with multilineage differentiation potential. Deletion of the murine Gata2-77 enhancer, with a human equivalent that causes leukemia, downregulates the transcription factor GATA2 and blocks progenitor differentiation into erythrocytes, megakaryocytes, basophils, and granulocytes, but not macrophages. Using multiomics and single-cell analyses, we demonstrated that the enhancer orchestrates a balance between pro- and anti-fate circuitry in single cells. By increasing GATA2 expression, the enhancer instigates a fate-promoting mechanism while abrogating an innate immunity-linked, fate-suppressing mechanism. During embryogenesis, the suppressing mechanism dominated in enhancer mutant progenitors, thus yielding progenitors with a predominant monocytic differentiation potential. Coordinating fate-promoting and -suppressing circuits therefore averts deconstruction of a multifate system into a monopotent system and maintains critical progenitor heterogeneity and functionality.

CTCF loss mediates unique DNA hypermethylation landscapes in human cancers.

BACKGROUND: The chromatin insulator CCCTC-binding factor (CTCF) displays tissue-specific DNA binding sites that regulate transcription and chromatin organization. Despite evidence linking CTCF to the protection of epigenetic states through barrier insulation, the impact of CTCF loss on genome-wide DNA methylation sites in human cancer remains undefined. RESULTS: Here, we demonstrate that prostate and breast cancers within The Cancer Genome Atlas (TCGA) exhibit frequent copy number loss of CTCF and that this loss is associated with increased DNA methylation events that occur preferentially at CTCF binding sites. CTCF sites differ among tumor types and result in tissue-specific methylation patterns with little overlap between breast and prostate cancers. DNA methylation and transcriptome profiling in vitro establish that forced downregulation of CTCF leads to spatially distinct DNA hypermethylation surrounding CTCF binding sites, loss of CTCF binding, and decreased gene expression that is also seen in human tumors. DNA methylation inhibition reverses loss of expression at these CTCF-regulated genes. CONCLUSION: These findings establish CTCF loss as a major mediator in directing localized DNA hypermethylation events in a tissue-specific fashion and further support its role as a driver of the cancer phenotype.

SURF: integrative analysis of a compendium of RNA-seq and CLIP-seq datasets highlights complex governing of alternative transcriptional regulation by RNA-binding proteins.

Advances in high-throughput profiling of RNA-binding proteins (RBPs) have resulted inCLIP-seq datasets coupled with transcriptome profiling by RNA-seq. However, analysis methods that integrate both types of data are lacking. We describe SURF, Statistical Utility for RBP Functions, for integrative analysis of large collections of CLIP-seq and RNA-seq data. We demonstrate SURF's ability to accurately detect differential alternative transcriptional regulation events and associate them to local protein-RNA interactions. We apply SURF to ENCODE RBP compendium and carry out downstream analysis with additional reference datasets. The results of this application are browsable at http://www.statlab.wisc.edu/shiny/surf/.

Single-nucleotide human disease mutation inactivates a blood-regenerative GATA2 enhancer.

The development and function of stem and progenitor cells that produce blood cells are vital in physiology. GATA-binding protein 2 (GATA2) mutations cause GATA-2 deficiency syndrome involving immunodeficiency, myelodysplastic syndrome, and acute myeloid leukemia. GATA-2 physiological activities necessitate that it be strictly regulated, and cell type-specific enhancers fulfill this role. The +9.5 intronic enhancer harbors multiple conserved cis-elements, and germline mutations of these cis-elements are pathogenic in humans. Since mechanisms underlying how GATA2 enhancer disease mutations impact hematopoiesis and pathology are unclear, we generated mouse models of the enhancer mutations. While a multi-motif mutant was embryonically lethal, a single-nucleotide Ets motif mutant was viable, and steady-state hematopoiesis was normal. However, the Ets motif mutation abrogated stem/progenitor cell regeneration following stress. These results reveal a new mechanism in human genetics, in which a disease predisposition mutation inactivates enhancer regenerative activity, while sparing developmental activity. Mutational sensitization to stress that instigates hematopoietic failure constitutes a paradigm for GATA-2 deficiency syndrome and other contexts of GATA-2-dependent pathogenesis.

Mechanisms of erythrocyte development and regeneration: implications for regenerative medicine and beyond.

Hemoglobin-expressing erythrocytes (red blood cells) act as fundamental metabolic regulators by providing oxygen to cells and tissues throughout the body. Whereas the vital requirement for oxygen to support metabolically active cells and tissues is well established, almost nothing is known regarding how erythrocyte development and function impact regeneration. Furthermore, many questions remain unanswered relating to how insults to hematopoietic stem/progenitor cells and erythrocytes can trigger a massive regenerative process termed 'stress erythropoiesis' to produce billions of erythrocytes. Here, we review the cellular and molecular mechanisms governing erythrocyte development and regeneration, and discuss the potential links between these events and other regenerative processes.

O2 availability impacts iron homeostasis in Escherichia coli.

The ferric-uptake regulator (Fur) is an Fe2+-responsive transcription factor that coordinates iron homeostasis in many bacteria. Recently, we reported that expression of the Escherichia coli Fur regulon is also impacted by O2 tension. Here, we show that for most of the Fur regulon, Fur binding and transcriptional repression increase under anaerobic conditions, suggesting that Fur is controlled by O2 availability. We found that the intracellular, labile Fe2+ pool was higher under anaerobic conditions compared with aerobic conditions, suggesting that higher Fe2+ availability drove the formation of more Fe2+-Fur and, accordingly, more DNA binding. O2 regulation of Fur activity required the anaerobically induced FeoABC Fe2+ uptake system, linking increased Fur activity to ferrous import under iron-sufficient conditions. The increased activity of Fur under anaerobic conditions led to a decrease in expression of ferric import systems. However, the combined positive regulation of the feoABC operon by ArcA and FNR partially antagonized Fur-mediated repression of feoABC under anaerobic conditions, allowing ferrous transport to increase even though Fur is more active. This design feature promotes a switch from ferric import to the more physiological relevant ferrous iron under anaerobic conditions. Taken together, we propose that the influence of O2 availability on the levels of active Fur adds a previously undescribed layer of regulation in maintaining cellular iron homeostasis.

GATA Factor-Regulated Samd14 Enhancer Confers Red Blood Cell Regeneration and Survival in Severe Anemia.

An enhancer with amalgamated E-box and GATA motifs (+9.5) controls expression of the regulator of hematopoiesis GATA-2. While similar GATA-2-occupied elements are common in the genome, occupancy does not predict function, and GATA-2-dependent genetic networks are incompletely defined. A "+9.5-like" element resides in an intron of Samd14 (Samd14-Enh) encoding a sterile alpha motif (SAM) domain protein. Deletion of Samd14-Enh in mice strongly decreased Samd14 expression in bone marrow and spleen. Although steady-state hematopoiesis was normal, Samd14-Enh-/- mice died in response to severe anemia. Samd14-Enh stimulated stem cell factor/c-Kit signaling, which promotes erythrocyte regeneration. Anemia activated Samd14-Enh by inducing enhancer components and enhancer chromatin accessibility. Thus, a GATA-2/anemia-regulated enhancer controls expression of an SAM domain protein that confers survival in anemia. We propose that Samd14-Enh and an ensemble of anemia-responsive enhancers are essential for erythrocyte regeneration in stress erythropoiesis, a vital process in pathologies, including ß-thalassemia, myelodysplastic syndrome, and viral infection.

Regulation of WNT Signaling by VSX2 During Optic Vesicle Patterning in Human Induced Pluripotent Stem Cells.

Few gene targets of Visual System Homeobox 2 (VSX2) have been identified despite its broad and critical role in the maintenance of neural retina (NR) fate during early retinogenesis. We performed VSX2 ChIP-seq and ChIP-PCR assays on early stage optic vesicle-like structures (OVs) derived from human iPS cells (hiPSCs), which highlighted WNT pathway genes as direct regulatory targets of VSX2. Examination of early NR patterning in hiPSC-OVs from a patient with a functional null mutation in VSX2 revealed mis-expression and upregulation of WNT pathway components and retinal pigmented epithelium (RPE) markers in comparison to control hiPSC-OVs. Furthermore, pharmacological inhibition of WNT signaling rescued the early mutant phenotype, whereas augmentation of WNT signaling in control hiPSC-OVs phenocopied the mutant. These findings reveal an important role for VSX2 as a regulator of WNT signaling and suggest that VSX2 may act to maintain NR identity at the expense of RPE in part by direct repression of WNT pathway constituents. Stem Cells 2016;34:2625-2634.

Impact of Anaerobiosis on Expression of the Iron-Responsive Fur and RyhB Regulons.

UNLABELLED: Iron, a major protein cofactor, is essential for most organisms. Despite the well-known effects of O2 on the oxidation state and solubility of iron, the impact of O2 on cellular iron homeostasis is not well understood. Here we report that in Escherichia coli K-12, the lack of O2 dramatically changes expression of genes controlled by the global regulators of iron homeostasis, the transcription factor Fur and the small RNA RyhB. Using chromatin immunoprecipitation sequencing (ChIP-seq), we found anaerobic conditions promote Fur binding to more locations across the genome. However, by expression profiling, we discovered that the major effect of anaerobiosis was to increase the magnitude of Fur regulation, leading to increased expression of iron storage proteins and decreased expression of most iron uptake pathways and several Mn-binding proteins. This change in the pattern of gene expression also correlated with an unanticipated decrease in Mn in anaerobic cells. Changes in the genes posttranscriptionally regulated by RyhB under aerobic and anaerobic conditions could be attributed to O2-dependent changes in transcription of the target genes: aerobic RyhB targets were enriched in iron-containing proteins associated with aerobic energy metabolism, whereas anaerobic RyhB targets were enriched in iron-containing anaerobic respiratory functions. Overall, these studies showed that anaerobiosis has a larger impact on iron homeostasis than previously anticipated, both by expanding the number of direct Fur target genes and the magnitude of their regulation and by altering the expression of genes predicted to be posttranscriptionally regulated by the small RNA RyhB under iron-limiting conditions. IMPORTANCE: Microbes and host cells engage in an "arms race" for iron, an essential nutrient that is often scarce in the environment. Studies of iron homeostasis have been key to understanding the control of iron acquisition and the downstream pathways that enable microbes to compete for this valuable resource. Here we report that O2 availability affects the gene expression programs of two Escherichia coli master regulators that function in iron homeostasis: the transcription factor Fur and the small RNA regulator RyhB. Fur appeared to be more active under anaerobic conditions, suggesting a change in the set point for iron homeostasis. RyhB preferentially targeted iron-containing proteins of respiration-linked pathways, which are differentially expressed under aerobic and anaerobic conditions. Such findings may be relevant to the success of bacteria within their hosts since zones of reduced O2 may actually reduce bacterial iron demands, making it easier to win the arms race for iron.

Multivariate Boosting for Integrative Analysis of High-Dimensional Cancer Genomic Data.

In this paper, we propose a novel multivariate component-wise boosting method for fitting multivariate response regression models under the high-dimension, low sample size setting. Our method is motivated by modeling the association among different biological molecules based on multiple types of high-dimensional genomic data. Particularly, we are interested in two applications: studying the influence of DNA copy number alterations on RNA transcript levels and investigating the association between DNA methylation and gene expression. For this purpose, we model the dependence of the RNA expression levels on DNA copy number alterations and the dependence of gene expression on DNA methylation through multivariate regression models and utilize boosting-type method to handle the high dimensionality as well as model the possible nonlinear associations. The performance of the proposed method is demonstrated through simulation studies. Finally, our multivariate boosting method is applied to two breast cancer studies.