Gregory, Brian D
Email Address
ORCID
Disciplines
37 results
Search Results
Now showing 1 - 10 of 37
Publication CoRAL: Predicting Non-Coding RNAs from Small RNA-Sequencing Data(2013-08-01) Leung, Yuk Y; Ryvkin, Paul; Ungar, Lyle H; Gregory, Brian D; Wang, Li-SanThe surprising observation that virtually the entire human genome is transcribed means we know little about the function of many emerging classes of RNAs, except their astounding diversities. Traditional RNA function prediction methods rely on sequence or alignment information, which are limited in their abilities to classify the various collections of non-coding RNAs (ncRNAs). To address this, we developed Classification of RNAs by Analysis of Length (CoRAL), a machine learning-based approach for classification of RNA molecules. CoRAL uses biologically interpretable features including fragment length and cleavage specificity to distinguish between different ncRNA populations. We evaluated CoRAL using genome-wide small RNA sequencing data sets from four human tissue types and were able to classify six different types of RNAs with ∼80% cross-validation accuracy. Analysis by CoRAL revealed that microRNAs, small nucleolar and transposon-derived RNAs are highly discernible and consistent across all human tissue types assessed, whereas long intergenic ncRNAs, small cytoplasmic RNAs and small nuclear RNAs show less consistent patterns. The ability to reliably annotate loci across tissue types demonstrates the potential of CoRAL to characterize ncRNAs using small RNA sequencing data in less well-characterized organisms.Publication The Functions of RNA-Dependent RNA Polymerases in Arabidopsis(2011-01-01) Willmann, Matthew R; Endres, Matthew W; Cook, Rebecca T; Gregory, Brian DOne recently identified mechanism that regulates mRNA abundance is RNA silencing, and pioneering work in Arabidopsis thaliana and other genetic model organisms helped define this process. RNA silencing pathways are triggered by either self-complementary fold-back structures or the production of double-stranded RNA (dsRNA) that gives rise to small RNAs (smRNAs) known as microRNAs (miRNAs) or small-interfering RNAs (siRNAs). These smRNAs direct sequence-specific regulation of various gene transcripts, repetitive sequences, viruses, and mobile elements via RNA cleavage, translational inhibition, or transcriptional silencing through DNA methylation and heterochromatin formation. Early genetic screens in Arabidopsis were instrumental in uncovering numerous proteins required for these important regulatory pathways. Among the factors identified by these studies were RNA-dependent RNA polymerases (RDRs), which are proteins that synthesize siRNA-producing dsRNA molecules using a single-stranded RNA (ssRNA) molecule as a template. Recently, a growing body of evidence has implicated RDR-dependent RNA silencing in many different aspects of plant biology ranging from reproductive development to pathogen resistance. Here, we focus on the specific functions of the six Arabidopsis RDRs in RNA silencing, their ssRNA substrates and resulting RDR-dependent smRNAs, and the numerous biological functions of these proteins in plant development and stress responses.Publication Lamin B1 Depletion in Senescent Cells Triggers Large-Scale Changes in Gene Expression and the Chromatin Landscape(2013-08-15) Shah, Parisha P; Donahue, Greg; Capell, Brian C; Aggarwala, Varun; Otte, Gabriel; Nelson, David M; Gregory, Brian D; Cao, Kajia; Berger, Shelley L; Cruickshanks, Hazel A; Rai, Taranjit Singh; McBryan, Tony; Adams, Peter DSenescence is a stable proliferation arrest, associated with an altered secretory pathway, thought to promote tumor suppression and tissue aging. While chromatin regulation and lamin B1 down-regulation have been implicated as senescence effectors, functional interactions between them are poorly understood. We compared genome-wide Lys4 trimethylation on histone H3 (H3K4me3) and H3K27me3 distributions between proliferating and senescent human cells and found dramatic differences in senescence, including large-scale domains of H3K4me3- and H3K27me3-enriched “mesas” and H3K27me3-depleted “canyons.” Mesas form at lamin B1-associated domains (LADs) in replicative senescence and oncogene-induced senescence and overlap DNA hypomethylation regions in cancer, suggesting that pre-malignant senescent chromatin changes foreshadow epigenetic cancer changes. Hutchinson-Gilford progeria syndrome fibroblasts (mutant lamin A) also show evidence of H3K4me3 mesas, suggesting a link between premature chromatin changes and accelerated cell senescence. Canyons mostly form between LADs and are enriched in genes and enhancers. H3K27me3 loss is correlated with up-regulation of key senescence genes, indicating a link between global chromatin changes and local gene expression regulation. Lamin B1 reduction in proliferating cells triggers senescence and formation of mesas and canyons. Our data illustrate profound chromatin reorganization during senescence and suggest that lamin B1 down-regulation in senescence is a key trigger of global and local chromatin changes that impact gene expression, aging, and cancer.Publication Transcriptomic Changes Due to Cytoplasmic TDP-43 Expression Revel Dysregulation of Histone Transcripts and Nuclear Chromatin(2015-01-01) Amlie-Wolf, Alexandre; Ryvkin, Paul; Tong, Rui; Suh, EunRan; Xu, Yan; Van Deerlin, Vivianna M; Gregory, Brian D; Trojanowski, John Q; Lee, Virginia Man-Yee; Wang, Li-San; Lee, Edward B; Dragomir, Isabelle; Kwong, Linda KAR DNA-binding protein 43 (TDP-43) is normally a nuclear RNA-binding protein that exhibits a range of functions including regulation of alternative splicing, RNA trafficking, and RNA stability. However, in amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration with TDP-43 inclusions (FTLD-TDP), TDP-43 is abnormally phosphorylated, ubiquitinated, and cleaved, and is mislocalized to the cytoplasm where it forms distinctive aggregates. We previously developed a mouse model expressing human TDP-43 with a mutation in its nuclear localization signal (ΔNLS-hTDP-43) so that the protein preferentially localizes to the cytoplasm. These mice did not exhibit a significant number of cytoplasmic aggregates, but did display dramatic changes in gene expression as measured by microarray, suggesting that cytoplasmic TDP-43 may be associated with a toxic gain-of-function. Here, we analyze new RNA-sequencing data from the ΔNLS-hTDP-43 mouse model, together with published RNA-sequencing data obtained previously from TDP-43 antisense oligonucleotide (ASO) knockdown mice to investigate further the dysregulation of gene expression in the ΔNLS model. This analysis reveals that the transcriptomic effects of the overexpression of the ΔNLS-hTDP-43 transgene are likely due to a gain of cytoplasmic function. Moreover, cytoplasmic TDP-43 expression alters transcripts that regulate chromatin assembly, the nucleolus, lysosomal function, and histone 3’ untranslated region (UTR) processing. These transcriptomic alterations correlate with observed histologic abnormalities in heterochromatin structure and nuclear size in transgenic mouse and human brains.Publication Regulatory Impact of RNA Secondary Structure across the Arabidopsis Transcriptome(2012-11-01) Li, Fan; Vandivier, Lee E; Gregory, Brian D; Willmann, Matthew R; Chen, YingThe secondary structure of an RNA molecule plays an integral role in its maturation, regulation, and function. However, the global influence of this feature on plant gene expression is still largely unclear. Here, we use a high-throughput, sequencing-based, structure-mapping approach in conjunction with transcriptome-wide sequencing of rRNA-depleted (RNA sequencing), small RNA, and ribosome-bound RNA populations to investigate the impact of RNA secondary structure on gene expression regulation in Arabidopsis thaliana. From this analysis, we find that highly unpaired and paired RNAs are strongly correlated with euchromatic and heterochromatic epigenetic histone modifications, respectively, providing evidence that secondary structure is necessary for these RNA-mediated posttranscriptional regulatory pathways. Additionally, we uncover key structural patterns across protein-coding transcripts that indicate RNA folding demarcates regions of protein translation and likely affects microRNA-mediated regulation of mRNAs in this model plant. We further reveal that RNA folding is significantly anticorrelated with overall transcript abundance, which is often due to the increased propensity of highly structured mRNAs to be degraded and/or processed into small RNAs. Finally, we find that secondary structure affects mRNA translation, suggesting that this feature regulates plant gene expression at multiple levels. These findings provide a global assessment of RNA folding and its significant regulatory effects in a plant transcriptome.Publication The Macronuclear Genome of Stentor coeruleus Reveals Tiny Introns in a Giant Cell(2017-02-20) Slabodnick, Mark M; Ruby, J. G; Reiff, Sarah B; Swart, Estienne C; Gosai, Sager J; Prabakaran, Sudhakaran; Witkowska, Ewa; Larue, Graham E; Gregory, Brian D; Nowacki, Mariusz; Derisi, Joseph; Roy, Scott W; Marshall, Wallace F; Sood, PranidhiThe giant, single-celled organism Stentor coeruleus has a long history as a model system for studying pattern formation and regeneration in single cells. Stentor [1, 2] is a heterotrichous ciliate distantly related to familiar ciliate models, such as Tetrahymena or Paramecium. The primary distinguishing feature of Stentor is its incredible size: a single cell is 1 mm long. Early developmental biologists, including T.H. Morgan [3], were attracted to the system because of its regenerative abilities—if large portions of a cell are surgically removed, the remnant reorganizes into a normal-looking but smaller cell with correct proportionality [2, 3]. These biologists were also drawn to Stentor because it exhibits a rich repertoire of behaviors, including light avoidance, mechanosensitive contraction, food selection, and even the ability to habituate to touch, a simple form of learning usually seen in higher organisms [4]. While early microsurgical approaches demonstrated a startling array of regenerative and morphogenetic processes in this single-celled organism, Stentor was never developed as a molecular model system. We report the sequencing of the Stentor coeruleus macronuclear genome and reveal key features of the genome. First, we find that Stentor uses the standard genetic code, suggesting that ciliate-specific genetic codes arose after Stentor branched from other ciliates. We also discover that ploidy correlates with Stentor’s cell size. Finally, in the Stentor genome, we discover the smallest spliceosomal introns reported for any species. The sequenced genome opens the door to molecular analysis of single-cell regeneration in Stentor.Publication HAMR: High-Throughput Annotation of Modified Ribonucleotides(2013-12-01) Ryvkin, Paul; Leung, Yuk Y; Childress, Micah; Valladares, Otto; Gregory, Brian D; Wang, Li-San; Silverman, Ian M; Dragomir, IsabelleRNA is often altered post-transcriptionally by the covalent modification of particular nucleotides; these modifications are known to modulate the structure and activity of their host RNAs. The recent discovery that an RNA methyl-6 adenosine demethylase (FTO) is a risk gene in obesity has brought to light the significance of RNA modifications to human biology. These noncanonical nucleotides, when converted to cDNA in the course of RNA sequencing, can produce sequence patterns that are distinguishable from simple base-calling errors. To determine whether these modifications can be detected in RNA sequencing data, we developed a method that can not only locate these modifications transcriptome-wide with single nucleotide resolution, but can also differentiate between different classes of modifications. Using small RNA-seq data we were able to detect 92% of all known human tRNA modification sites that are predicted to affect RT activity. We also found that different modifications produce distinct patterns of cDNA sequence, allowing us to differentiate between two classes of adenosine and two classes of guanine modifications with 98% and 79% accuracy, respectively. To show the robustness of this method to sample preparation and sequencing methods, as well as to organismal diversity, we applied it to a publicly available yeast data set and achieved similar levels of accuracy. We also experimentally validated two novel and one known 3-methylcytosine (3mC) sites predicted by HAMR in human tRNAs. Researchers can now use our method to identify and characterize RNA modifications using only RNA-seq data, both retrospectively and when asking questions specifically about modified RNA.Publication Analysis of the Host Transcriptome from Demyelinating Spinal Cord of Murine Coronavirus-Infected Mice(2013-09-18) Elliott, Ruth; Li, Fan; Dragomir, Isabelle; Gregory, Brian D; Chua, Ming Ming W; Weiss, Susan RPersistent infection of the mouse central nervous system (CNS) with mouse hepatitis virus (MHV) induces a demyelinating disease pathologically similar to multiple sclerosis and is therefore used as a model system. There is little information regarding the host factors that correlate with and contribute to MHV-induced demyelination. Here, we detail the genes and pathways associated with MHV-induced demyelinating disease in the spinal cord. High-throughput sequencing of the host transcriptome revealed that demyelination is accompanied by numerous transcriptional changes indicative of immune infiltration as well as changes in the cytokine milieu and lipid metabolism. We found evidence that a Th1-biased cytokine/chemokine response and eicosanoid-derived inflammation accompany persistent MHV infection and that antigen presentation is ongoing. Interestingly, increased expression of genes involved in lipid transport, processing, and catabolism, including some with known roles in neurodegenerative diseases, coincided with demyelination. Lastly, expression of several genes involved in osteoclast or bone-resident macrophage function, most notably TREM2 and DAP12, was upregulated in persistently infected mouse spinal cord. This study highlights the complexity of the host antiviral response, which accompany MHV-induced demyelination, and further supports previous findings that MHV-induced demyelination is immune-mediated. Interestingly, these data suggest a parallel between bone reabsorption by osteoclasts and myelin debris clearance by microglia in the bone and the CNS, respectively. To our knowledge, this is the first report of using an RNA-seq approach to study the host CNS response to persistent viral infection.Publication The Arabidopsis dwf/ste1 Mutant is Defective in the Δ7 Sterol C-5 Desaturation Step Leading to Brassinosteroid Biosynthesis(1999-02-01) Choe, Sunghwa; Noguchi, Takahiro; Fujioka, Shozo; Takatsuto, Suguru; Tissier, Christophe P; Gregory, Brian D; Ross, Amanda S; Tanaka, Atsushi; Yoshida, Shigeo; Tax, Frans E; Feldmann, Kenneth ALesions in brassinosteroid (BR) biosynthetic genes result in characteristic dwarf phenotypes in plants. Understanding the regulation of BR biosynthesis demands continued isolation and characterization of mutants corresponding to the genes involved in BR biosynthesis. Here, we present analysis of a novel BR biosynthetic locus, dwarf7 (dwf7). Feeding studies with BR biosynthetic intermediates and analysis of endogenous levels of BR and sterol biosynthetic intermediates indicate that the defective step in dwf7-1 resides before the production of 24-methylenecholesterol in the sterol biosynthetic pathway. Furthermore, results from feeding studies with 13C-labeled mevalonic acid and compactin show that the defective step is specifically the Δ7 sterol C-5 desaturation, suggesting that dwf7 is an allele of the previously cloned STEROL1 (STE1) gene. Sequencing of the STE1 locus in two dwf7 mutants revealed premature stop codons in the first (dwf7-2) and the third (dwf7-1) exons. Thus, the reduction of BRs in dwf7 is due to a shortage of substrate sterols and is the direct cause of the dwarf phenotype in dwf7.Publication Spatial and Functional Relationships Among Pol V-Associated Loci, Pol IV-Dependent siRNAs, and Cytosine Methylation in the Arabidopsis Epigenome(2012-08-15) Wierzbicki, Andrzej T; Cocklin, Ross; Mayampurath, Anoop; Lister, Ryan; Rowley, M. J; Gregory, Brian D; Ecker, Joseph R; Tang, Haixu; Pikaard, Craig SMultisubunit RNA polymerases IV and V (Pols IV and V) mediate RNA-directed DNA methylation and transcriptional silencing of retrotransposons and heterochromatic repeats in plants. We identified genomic sites of Pol V occupancy in parallel with siRNA deep sequencing and methylcytosine mapping, comparing wild-type plants with mutants defective for Pol IV, Pol V, or both Pols IV and V. Approximately 60% of Pol V-associated regions encompass regions of 24-nucleotide (nt) siRNA complementarity and cytosine methylation, consistent with cytosine methylation being guided by base-pairing of Pol IV-dependent siRNAs with Pol V transcripts. However, 27% of Pol V peaks do not overlap sites of 24-nt siRNA biogenesis or cytosine methylation, indicating that Pol V alone does not specify sites of cytosine methylation. Surprisingly, the number of methylated CHH motifs, a hallmark of RNA-directed de novo methylation, is similar in wild-type plants and Pol IV or Pol V mutants. In the mutants, methylation is lost at 50%–60% of the CHH sites that are methylated in the wild type but is gained at new CHH positions, primarily in pericentromeric regions. These results indicate that Pol IV and Pol V are not required for cytosine methyltransferase activity but shape the epigenome by guiding CHH methylation to specific genomic sites.