Transcriptome analysis of ciliary-dependent MCH signaling in differentiating 3T3-L1 pre-adipocytes.

An understanding of adipocyte responsiveness to G-protein-coupled receptor-(GPCR) derived signals must take into consideration the role of membrane microenvironments; that individual sub-populations of proteins may vary significantly across different regions of the cell, and that cell differentiation alters those microenvironments. 3T3-L1 pre-adipocytes undergo a dramatic phenotypic transformation during differentiation into adipocytes, requiring the development of a transient primary cilium. We demonstrate that melanin-concentrating hormone (MCH) receptor 1, a GPCR that stimulates appetite, translocates to the transient primary cilium during early 3T3-L1 cell adipogenesis. Furthermore, we used RNA-Seq to investigate whether MCH signaling is influenced by its receptor localization and whether MCH can influence the transcriptome of early adipocyte development. We found that MCH signaling is sensitive to receptor localization to cilia, and this alters the adipogenic transcriptional program. Also, novel MCH signaling pathways in 3T3-L1 cells are identified, including those for circadian rhythm, the inflammatory response, and ciliary biogenesis. The presence of active MCH-signaling pathways in pre-adipocytes and the discovery that these pathways intersect with the early adipogenic program, among other newly-identified signaling pathways, suggests that the use of MCH receptor 1 antagonists for clinical interventions may have unintended consequences on adipose tissue development.

AML suppresses hematopoiesis by releasing exosomes that contain microRNAs targeting c-MYB.

Exosomes are paracrine regulators of the tumor microenvironment and contain complex cargo. We previously reported that exosomes released from acute myeloid leukemia (AML) cells can suppress residual hematopoietic stem and progenitor cell (HSPC) function indirectly through stromal reprogramming of niche retention factors. We found that the systemic loss of hematopoietic function is also in part a consequence of AML exosome-directed microRNA (miRNA) trafficking to HSPCs. Exosomes isolated from cultured AML or the plasma from mice bearing AML xenografts exhibited enrichment of miR-150 and miR-155. HSPCs cocultured with either of these exosomes exhibited impaired clonogenicity, through the miR-150- and miR-155-mediated suppression of the translation of transcripts encoding c-MYB, a transcription factor involved in HSPC differentiation and proliferation. To discover additional miRNA targets, we captured miR-155 and its target transcripts by coimmunoprecipitation with an attenuated RNA-induced silencing complex (RISC)-trap, followed by high-throughput sequencing. This approach identified known and previously unknown miR-155 target transcripts. Integration of the miR-155 targets with information from the protein interaction database STRING revealed proteins indirectly affected by AML exosome-derived miRNA. Our findings indicate a direct effect of AML exosomes on HSPCs that, through a stroma-independent mechanism, compromises hematopoiesis. Furthermore, combining miRNA target data with protein-protein interaction data may be a broadly applicable strategy to define the effects of exosome-mediated trafficking of regulatory molecules within the tumor microenvironment.

Transcriptional Profiling of Newly Generated Dentate Granule Cells Using TU Tagging Reveals Pattern Shifts in Gene Expression during Circuit Integration.

Despite representing only a small fraction of hippocampal granule cells, adult-generated newborn granule cells have been implicated in learning and memory (Aimone et al., 2011). Newborn granule cells undergo functional maturation and circuit integration over a period of weeks. However, it is difficult to assess the accompanying gene expression profiles in vivo with high spatial and temporal resolution using traditional methods. Here we used a novel method ["thiouracil (TU) tagging"] to map the profiles of nascent mRNAs in mouse immature newborn granule cells compared with mature granule cells. We targeted a nonmammalian uracil salvage enzyme, uracil phosphoribosyltransferase, to newborn neurons and mature granule cells using retroviral and lentiviral constructs, respectively. Subsequent injection of 4-TU tagged nascent RNAs for analysis by RNA sequencing. Several hundred genes were significantly enhanced in the retroviral dataset compared with the lentiviral dataset. We compared a selection of the enriched genes with steady-state levels of mRNAs using quantitative PCR. Ontology analysis revealed distinct patterns of nascent mRNA expression, with newly generated immature neurons showing enhanced expression for genes involved in synaptic function, and neural differentiation and development, as well as genes not previously associated with granule cell maturation. Surprisingly, the nascent mRNAs enriched in mature cells were related to energy homeostasis and metabolism, presumably indicative of the increased energy demands of synaptic transmission and their complex dendritic architecture. The high spatial and temporal resolution of our modified TU-tagging method provides a foundation for comparison with steady-state RNA analyses by traditional transcriptomic approaches in defining the functional roles of newborn neurons.

miR-218 is essential to establish motor neuron fate as a downstream effector of Isl1-Lhx3.

While microRNAs have emerged as an important component of gene regulatory networks, it remains unclear how microRNAs collaborate with transcription factors in the gene networks that determines neuronal cell fate. Here we show that in the developing spinal cord, the expression of miR-218 is directly upregulated by the Isl1-Lhx3 complex, which drives motor neuron fate. Inhibition of miR-218 suppresses the generation of motor neurons in both chick neural tube and mouse embryonic stem cells, suggesting that miR-218 plays a crucial role in motor neuron differentiation. Results from unbiased RISC-trap screens, in vivo reporter assays and overexpression studies indicated that miR-218 directly represses transcripts that promote developmental programs for interneurons. In addition, we found that miR-218 activity is required for Isl1-Lhx3 to effectively induce motor neurons and suppress interneuron fates. Together our results reveal an essential role of miR-218 as a downstream effector of the Isl1-Lhx3 complex in establishing motor neuron identity.

Crucial roles of mixed-lineage leukemia 3 and 4 as epigenetic switches of the hepatic circadian clock controlling bile acid homeostasis in mice.

UNLABELLED: The histone H3-lysine-4 methyltransferase mixed-lineage leukemia 3 (MLL3) and its closest homolog, MLL4 (aka KMT2D), belong to two homologous transcriptional coactivator complexes, named MLL3 and MLL4 complexes, respectively. MLL3 plays crucial roles in multiple metabolic processes. However, the physiological roles of MLL4 in metabolism and the relationship between MLL3 and MLL4 in metabolic gene regulation are unclear. To address these issues, we analyzed the phenotypes of newly generated MLL4 mutant mice, along with MLL3 mutant and MLL3;MLL4 compound mutant mice. We also performed comparative genome-wide transcriptome analyses in livers of MLL3, MLL4, and MLL3;MLL4 mutant mice. These analyses revealed that MLL3 and MLL4 complexes are key epigenetic regulators of common metabolic processes and the hepatic circadian clock. Subsequent mechanistic analyses uncovered that MLL3/4 complexes function as pivotal coactivators of the circadian transcription factors (TFs), retinoid-related orphan receptor (ROR)-α and -γ, in the hepatic circadian clock. Consistent with disturbed hepatic clock gene expression in MLL4 mutant mice, we found that rhythmic fluctuation of hepatic and serum bile acid (BA) levels over the circadian cycle is abolished in MLL4 mutant mice. Our analyses also demonstrate that MLL4 primarily impinges on hepatic BA production among several regulatory pathways to control BA homeostasis. Together, our results provide strong in vivo support for important roles of both MLL3 and MLL4 in similar metabolic pathways. CONCLUSION: Both MLL3 and MLL4 complexes act as major epigenetic regulators of diverse metabolic processes (including circadian control of bile acid homeostasis) and as critical transcriptional coactivators of the circadian TFs, RORs.

Isl1 directly controls a cholinergic neuronal identity in the developing forebrain and spinal cord by forming cell type-specific complexes.

The establishment of correct neurotransmitter characteristics is an essential step of neuronal fate specification in CNS development. However, very little is known about how a battery of genes involved in the determination of a specific type of chemical-driven neurotransmission is coordinately regulated during vertebrate development. Here, we investigated the gene regulatory networks that specify the cholinergic neuronal fates in the spinal cord and forebrain, specifically, spinal motor neurons (MNs) and forebrain cholinergic neurons (FCNs). Conditional inactivation of Isl1, a LIM homeodomain factor expressed in both differentiating MNs and FCNs, led to a drastic loss of cholinergic neurons in the developing spinal cord and forebrain. We found that Isl1 forms two related, but distinct types of complexes, the Isl1-Lhx3-hexamer in MNs and the Isl1-Lhx8-hexamer in FCNs. Interestingly, our genome-wide ChIP-seq analysis revealed that the Isl1-Lhx3-hexamer binds to a suite of cholinergic pathway genes encoding the core constituents of the cholinergic neurotransmission system, such as acetylcholine synthesizing enzymes and transporters. Consistently, the Isl1-Lhx3-hexamer directly coordinated upregulation of cholinergic pathways genes in embryonic spinal cord. Similarly, in the developing forebrain, the Isl1-Lhx8-hexamer was recruited to the cholinergic gene battery and promoted cholinergic gene expression. Furthermore, the expression of the Isl1-Lhx8-complex enabled the acquisition of cholinergic fate in embryonic stem cell-derived neurons. Together, our studies show a shared molecular mechanism that determines the cholinergic neuronal fate in the spinal cord and forebrain, and uncover an important gene regulatory mechanism that directs a specific neurotransmitter identity in vertebrate CNS development.

STAT3 promotes motor neuron differentiation by collaborating with motor neuron-specific LIM complex.

The motor neuron (MN)-hexamer complex consisting of LIM homeobox 3, Islet-1, and nuclear LIM interactor is a key determinant of motor neuron specification and differentiation. To gain insights into the transcriptional network in motor neuron development, we performed a genome-wide ChIP-sequencing analysis and found that the MN-hexamer directly regulates a wide array of motor neuron genes by binding to the HxRE (hexamer response element) shared among the target genes. Interestingly, STAT3-binding motif is highly enriched in the MN-hexamer-bound peaks in addition to the HxRE. We also found that a transcriptionally active form of STAT3 is expressed in embryonic motor neurons and that STAT3 associates with the MN-hexamer, enhancing the transcriptional activity of the MN-hexamer in an upstream signal-dependent manner. Correspondingly, STAT3 was needed for motor neuron differentiation in the developing spinal cord. Together, our studies uncover crucial gene regulatory mechanisms that couple MN-hexamer and STAT-activating extracellular signals to promote motor neuron differentiation in vertebrate spinal cord.

Capturing microRNA targets using an RNA-induced silencing complex (RISC)-trap approach.

Identifying targets is critical for understanding the biological effects of microRNA (miRNA) expression. The challenge lies in characterizing the cohort of targets for a specific miRNA, especially when targets are being actively down-regulated in miRNA- RNA-induced silencing complex (RISC)-messengerRNA (mRNA) complexes. We have developed a robust and versatile strategy called RISCtrap to stabilize and purify targets from this transient interaction. Its utility was demonstrated by determining specific high-confidence target datasets for miR-124, miR-132, and miR-181 that contained known and previously unknown transcripts. Two previously unknown miR-132 targets identified with RISCtrap, adaptor protein CT10 regulator of kinase 1 (CRK1) and tight junction-associated protein 1 (TJAP1), were shown to be endogenously regulated by miR-132 in adult mouse forebrain. The datasets, moreover, differed in the number of targets and in the types and frequency of microRNA recognition element (MRE) motifs, thus revealing a previously underappreciated level of specificity in the target sets regulated by individual miRNAs.

Fusion protein Isl1-Lhx3 specifies motor neuron fate by inducing motor neuron genes and concomitantly suppressing the interneuron programs.

Combinatorial transcription codes generate the myriad of cell types during development and thus likely provide crucial insights into directed differentiation of stem cells to a specific cell type. The LIM complex composed of Isl1 and Lhx3 directs the specification of spinal motor neurons (MNs) in embryos. Here, we report that Isl1-Lhx3, a LIM-complex mimicking fusion, induces a signature of MN transcriptome and concomitantly suppresses interneuron differentiation programs, thereby serving as a potent and specific inducer of MNs in stem cells. We show that an equimolar ratio of Isl1 and Lhx3 and the LIM domain of Lhx3 are crucial for generating MNs without up-regulating interneuron genes. These led us to design Isl1-Lhx3, which maintains the desirable 1:1 ratio of Isl1 and Lhx3 and the LIM domain of Lhx3. Isl1-Lhx3 drives MN differentiation with high specificity and efficiency in the spinal cord and embryonic stem cells, bypassing the need for sonic hedgehog (Shh). RNA-seq analysis revealed that Isl1-Lhx3 induces the expression of a battery of MN genes that control various functional aspects of MNs, while suppressing key interneuron genes. Our studies uncover a highly efficient method for directed MN generation and MN gene networks. Our results also demonstrate a general strategy of using embryonic transcription complexes for producing specific cell types from stem cells.

Supersplat--spliced RNA-seq alignment.

MOTIVATION: High-throughput sequencing technologies have recently made deep interrogation of expressed transcript sequences practical, both economically and temporally. Identification of intron/exon boundaries is an essential part of genome annotation, yet remains a challenge. Here, we present supersplat, a method for unbiased splice-junction discovery through empirical RNA-seq data. RESULTS: Using a genomic reference and RNA-seq high-throughput sequencing datasets, supersplat empirically identifies potential splice junctions at a rate of approximately 11.4 million reads per hour. We further benchmark the performance of the algorithm by mapping Illumina RNA-seq reads to identify introns in the genome of the reference dicot plant Arabidopsis thaliana and we demonstrate the utility of supersplat for de novo empirical annotation of splice junctions using the reference monocot plant Brachypodium distachyon. AVAILABILITY: Implemented in C++, supersplat source code and binaries are freely available on the web at http://mocklerlab-tools.cgrb.oregonstate.edu/.

Genome-wide mapping of alternative splicing in Arabidopsis thaliana.

Alternative splicing can enhance transcriptome plasticity and proteome diversity. In plants, alternative splicing can be manifested at different developmental stages, and is frequently associated with specific tissue types or environmental conditions such as abiotic stress. We mapped the Arabidopsis transcriptome at single-base resolution using the Illumina platform for ultrahigh-throughput RNA sequencing (RNA-seq). Deep transcriptome sequencing confirmed a majority of annotated introns and identified thousands of novel alternatively spliced mRNA isoforms. Our analysis suggests that at least approximately 42% of intron-containing genes in Arabidopsis are alternatively spliced; this is significantly higher than previous estimates based on cDNA/expressed sequence tag sequencing. Random validation confirmed that novel splice isoforms empirically predicted by RNA-seq can be detected in vivo. Novel introns detected by RNA-seq were substantially enriched in nonconsensus terminal dinucleotide splice signals. Alternative isoforms with premature termination codons (PTCs) comprised the majority of alternatively spliced transcripts. Using an example of an essential circadian clock gene, we show that intron retention can generate relatively abundant PTC(+) isoforms and that this specific event is highly conserved among diverse plant species. Alternatively spliced PTC(+) isoforms can be potentially targeted for degradation by the nonsense mediated mRNA decay (NMD) surveillance machinery or regulate the level of functional transcripts by the mechanism of regulated unproductive splicing and translation (RUST). We demonstrate that the relative ratios of the PTC(+) and reference isoforms for several key regulatory genes can be considerably shifted under abiotic stress treatments. Taken together, our results suggest that like in animals, NMD and RUST may be widespread in plants and may play important roles in regulating gene expression.

Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology.

Organellar DNA sequences are widely used in evolutionary and population genetic studies, however, the conservative nature of chloroplast gene and genome evolution often limits phylogenetic resolution and statistical power. To gain maximal access to the historical record contained within chloroplast genomes, we have adapted multiplex sequencing-by-synthesis (MSBS) to simultaneously sequence multiple genomes using the Illumina Genome Analyzer. We PCR-amplified approximately 120 kb plastomes from eight species (seven Pinus, one Picea) in 35 reactions. Pooled products were ligated to modified adapters that included 3 bp indexing tags and samples were multiplexed at four genomes per lane. Tagged microreads were assembled by de novo and reference-guided assembly methods, using previously published Pinus plastomes as surrogate references. Assemblies for these eight genomes are estimated at 88-94% complete, with an average sequence depth of 55x to 186x. Mononucleotide repeats interrupt contig assembly with increasing repeat length, and we estimate that the limit for their assembly is 16 bp. Comparisons to 37 kb of Sanger sequence show a validated error rate of 0.056%, and conspicuous errors are evident from the assembly process. This efficient sequencing approach yields high-quality draft genomes and should have immediate applicability to genomes with comparable complexity.

Network discovery pipeline elucidates conserved time-of-day-specific cis-regulatory modules.

Correct daily phasing of transcription confers an adaptive advantage to almost all organisms, including higher plants. In this study, we describe a hypothesis-driven network discovery pipeline that identifies biologically relevant patterns in genome-scale data. To demonstrate its utility, we analyzed a comprehensive matrix of time courses interrogating the nuclear transcriptome of Arabidopsis thaliana plants grown under different thermocycles, photocycles, and circadian conditions. We show that 89% of Arabidopsis transcripts cycle in at least one condition and that most genes have peak expression at a particular time of day, which shifts depending on the environment. Thermocycles alone can drive at least half of all transcripts critical for synchronizing internal processes such as cell cycle and protein synthesis. We identified at least three distinct transcription modules controlling phase-specific expression, including a new midnight specific module, PBX/TBX/SBX. We validated the network discovery pipeline, as well as the midnight specific module, by demonstrating that the PBX element was sufficient to drive diurnal and circadian condition-dependent expression. Moreover, we show that the three transcription modules are conserved across Arabidopsis, poplar, and rice. These results confirm the complex interplay between thermocycles, photocycles, and the circadian clock on the daily transcription program, and provide a comprehensive view of the conserved genomic targets for a transcriptional network key to successful adaptation.

Molecular interactions involving Escherichia coli nucleoside diphosphate kinase.

Nucleoside diphosphate kinase plays a distinctive metabolic role as the enzyme poised between the last reaction of deoxyribonucleoside triphosphate (dNTP) biosynthesis and the DNA polymerization apparatus. In bacteriophage T4 infection, NDP kinase is one of very few enzymes of host cell origin to participate in either dNTP synthesis or DNA replication. Yet NDP kinase forms specific contacts with phage-coded proteins of dNTP and DNA synthesis. This article summarizes work from our laboratory that identifies and characterizes these interactions. Despite these specific interactions, the enzyme appears to be dispensable, both for T4 replication and for growth of the host, Escherichia coli, because site-specific disruption of ndk, the structural gene for NDP kinase, does not interfere with growth of the host cell and only partly inhibits phage replication. However, ndk disruption unbalances the dNTP pools and stimulates mutagenesis. We discuss our attempts to understand the basis for this enhanced mutagenesis.

Adenylate kinase of Escherichia coli, a component of the phage T4 dNTP synthetase complex.

Adenylate kinase, which catalyzes the reversible ATP-dependent phosphorylation of AMP to ADP and dAMP to dADP, can also catalyze the conversion of nucleoside diphosphates to the corresponding triphosphates. Lu and Inouye (Lu, Q., and Inouye, M. (1996) Proc. Natl. Acad. Sci. U. S. A. 93, 5720-5725) showed that an Escherichia coli ndk mutant, lacking nucleoside diphosphate kinase, can use adenylate kinase as an alternative source of nucleoside triphosphates. Bacteriophage T4 can reproduce in an Escherichia coli ndk mutant, implying that adenylate kinase can meet a demand for deoxyribonucleoside triphosphates that increases by up to 10-fold as a result of T4 infection. In terms of kinetic linkage and specific protein-protein associations, NDP kinase is an integral component of T4 dNTP synthetase, a multienzyme complex containing phage-coded enzymes, which facilitates the synthesis of dNTPs and their flow into DNA. Here we asked whether, by similar criteria, adenylate kinase of the host cell is also a specific component of the complex. Experiments involving protein affinity chromatography, immunoprecipitation, optical biosensor measurements, and glutathione S-transferase pulldowns demonstrated direct interactions between adenylate kinase and several phage-coded enzymes, as well as E. coli nucleoside diphosphate kinase. These results identify adenylate kinase as a specific component of the complex. The rate of DNA synthesis after infection of an ndk mutant was found to be about 40% of the rate seen in wild-type infection, implying that complementation of the missing NDP kinase function by adenylate kinase is fairly efficient, but that adenylate kinase becomes rate-limiting for DNA synthesis when it is the sole source of dNTPs.

Protein-DNA interactions in the T4 dNTP synthetase complex dependent on gene 32 single-stranded DNA-binding protein.

Our laboratory has reported data suggesting a role for T4 phage gene 32 single-stranded DNA-binding protein in organizing a complex of deoxyribonucleotide-synthesizing enzymes at the replication fork. In this article we examined the effects of gene 32 ablation on the association of these enzymes with DNA-protein complexes. These experiments showed several deoxyribonucleotide-synthesizing enzymes to be present in DNA-protein complexes, with some of these associations being dependent on gene 32 protein. To further understand the role of gp32, we created amber mutations at codons 24 and 204 of gene 32, which encodes a 301-residue protein. We used the newly created mutants along with several experimental approaches--DNA-cellulose chromatography, immunoprecipitation, optical biosensor analysis and glutathione-S-transferase pulldowns--to identify relevant protein-protein and protein-DNA interactions. These experiments identified several proteins whose interactions with DNA depend on the presence of intact gp32, notably thymidylate synthase, dihydrofolate (DHF) reductase, ribonucleotide reductase (RNR) and Escherichia coli nucleoside diphosphate (NDP) kinase, and they also demonstrated direct associations between gp32 and RNR and NDP kinase, but not dCMP hydroxymethylase, deoxyribonucleoside monophosphate kinase, or DHF reductase. Taken together, the results support the hypothesis that the gene 32 protein helps to recruit enzymes of deoxyribonucleoside triphosphates synthesis to DNA replication sites.

Escherichia coli nucleoside diphosphate kinase interactions with T4 phage proteins of deoxyribonucleotide synthesis and possible regulatory functions.

In both prokaryotic and eukaryotic organisms, nucleoside diphosphate kinase is a multifunctional protein, with well defined functions in ribo- and deoxyribonucleoside triphosphate biosynthesis and more recently described functions in genetic and metabolic regulation, signal transduction, and DNA repair. This paper concerns two unusual properties of nucleoside diphosphate (NDP) kinase from Escherichia coli: 1) its ability to interact specifically with enzymes encoded by the virulent bacteriophage T4 and 2) its roles in regulating metabolism of the host cell. By means of optical biosensor analysis, fluorescence spectroscopy, immunoprecipitation, and glutathione S-transferase pull-down assays, we have shown that E. coli NDP kinase interacts directly with T4 thymidylate synthase, aerobic ribonucleotide reductase, dCTPase-dUTPase, gene 32 single-strand DNA-binding protein, and deoxycytidylate hydroxymethylase. The interactions with ribonucleotide reductase and with gp32 are enhanced by nucleoside triphosphates, suggesting that the integrity of the T4 dNTP synthetase complex in vivo is influenced by the composition of the nucleotide pool. The other investigations in this work stem from the unexpected finding that E. coli NDP kinase is dispensable for successful T4 phage infection, and they deal with two observations suggesting that the NDP kinase protein plays a genetic role in regulating metabolism of the host cell: 1) the elevation of CTP synthetase activity in an ndk mutant, in which the structural gene for NDP kinase is disrupted, and 2) the apparent ability of NDP kinase to suppress anaerobic growth in a pyruvate kinase-negative E. coli mutant. Our data indicate that the regulatory roles are metabolic, not genetic, in nature.