Single transcript-level metabolic responsive landscape of human liver transcriptome

Direct knowledge of gene regulation in human liver by metabolic stimuli could fundamentally advance our understanding of metabolic physiology. This information, however, is largely unknown, and this void is deeply rooted in a paradox that is caused by the inaccessibility of human liver to treatments and the insufficient annotation of liver transcriptome. Recent advances have uncovered immense complexity of transcriptome, i.e., multiple transcripts produced from one gene and extensive RNA modifications, which are all highly regulated and often condition-dependent. Establishing an inclusive annotation to study liver transcriptome dynamics thus requires human liver samples of diverse conditions, which are paradoxically unavailable due to the inaccessibility of human liver to treatments. In this work, we addressed these challenges by coupling an isogenic humanized mouse model with Nanopore single-molecule direct RNA sequencing (DRS). We first generated mice that carried humanized livers of identical genetic background, which were equivalent to clones of a single human liver, and then subjected the mice to representative metabolic treatments. We then analyzed the humanized livers with Nanopore DRS, which directly reads full-length native RNAs to determine the expression level, m6A modification and poly(A) tail length of all RNA transcript isoforms. Thus, our system allows for constructing a de novo annotation of human liver transcriptomes reflecting metabolic responses and studying transcriptome dynamics in conjunction. Our analysis uncovered a vast number of novel genes and transcripts that have not been previously reported. Our transcript-level analysis of human liver transcriptomes also identified a multitude of regulated metabolic pathways that were otherwise invisible using conventional short read RNA-seq. We also revealed for the first time the dynamic changes in m6A and poly(A) tail length of human liver transcripts many of which are transcribed from key metabolic genes. Furthermore, we performed comparative analyses of gene regulation between human and mouse and between two individuals using the liver-specific humanized mice. This revealed that transcriptome dynamics are highly species- and genetic background-dependent, which may only be faithfully studied in a humanized system that entails clones of the same human liver. Hence our work revealed a complex metabolic responsive landscape of human liver transcriptome and also provided a framework to understand transcriptome dynamics of human liver in response to physiologically relevant metabolic stimuli.

[1]  A. Dziembowski,et al.  Measuring the tail: Methods for poly(A) tail profiling , 2022, Wiley interdisciplinary reviews. RNA.

[2]  J. Zuber,et al.  Dysregulated RNA polyadenylation contributes to metabolic impairment in non-alcoholic fatty liver disease , 2022, Nucleic acids research.

[3]  Yan Lu,et al.  Analysis of N6-Methyladenosine Methylation Modification in Fructose-Induced Non-Alcoholic Fatty Liver Disease , 2021, Frontiers in Endocrinology.

[4]  Paul J. Harrison,et al.  Accurate expression quantification from nanopore direct RNA sequencing with NanoCount , 2021, Nucleic acids research.

[5]  Gloria M. Sheynkman,et al.  Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing , 2021, Cell reports.

[6]  Oriol Vinyals,et al.  Highly accurate protein structure prediction with AlphaFold , 2021, Nature.

[7]  P. Tontonoz,et al.  Hepatic transcriptional responses to fasting and feeding , 2021, Genes & development.

[8]  Cuiping Yang,et al.  The role of m6A modification in the biological functions and diseases , 2021, Signal Transduction and Targeted Therapy.

[9]  Beryl B. Cummings,et al.  Transcriptome variation in human tissues revealed by long-read sequencing , 2021, Nature.

[10]  Maite Huarte,et al.  Gene regulation by long non-coding RNAs and its biological functions , 2020, Nature reviews. Molecular cell biology.

[11]  H. Suemizu,et al.  Comparative Transcriptomics Analyses in Livers of Mice, Humans, and Humanized Mice Define Human-Specific Gene Networks , 2020, Cells.

[12]  J. Rinn,et al.  Identification of human long non-coding RNAs associated with nonalcoholic fatty liver disease and metabolic homeostasis. , 2020, The Journal of clinical investigation.

[13]  J. Simpson,et al.  New Twists in Detecting mRNA Modification Dynamics , 2020, Trends in Biotechnology.

[14]  Y. Gilad,et al.  Gene expression variability in human and chimpanzee populations share common determinants , 2020, bioRxiv.

[15]  Bao-Fa Sun,et al.  m6A Regulates Liver Metabolic Disorders and Hepatogenous Diabetes , 2020, bioRxiv.

[16]  M. Pirooznia,et al.  In vivo functional analysis of non-conserved human lncRNAs associated with cardiometabolic traits , 2020, Nature Communications.

[17]  Alexa B. R. McIntyre,et al.  Direct RNA sequencing reveals m6A modifications on adenovirus RNA are necessary for efficient splicing , 2019, Nature Communications.

[18]  Anton J. Enright,et al.  RNA modifications detection by comparative Nanopore direct RNA sequencing , 2019, Nature Communications.

[19]  Daniel A. Lorenz,et al.  Direct RNA sequencing enables m6A detection in endogenous transcript isoforms at base-specific resolution , 2019, RNA.

[20]  J. Hadfield,et al.  RNA sequencing: the teenage years , 2019, Nature Reviews Genetics.

[21]  Xiaonan Zhang,et al.  Long-term functional maintenance of primary human hepatocytes in vitro , 2019, Science.

[22]  J. Jelsing,et al.  Hepatic transcriptome signatures in patients with varying degrees of nonalcoholic fatty liver disease compared with healthy normal-weight individuals. , 2019, American journal of physiology. Gastrointestinal and liver physiology.

[23]  M. Robinson,et al.  A comprehensive examination of Nanopore native RNA sequencing for characterization of complex transcriptomes , 2019, Nature Communications.

[24]  Ryan R. Wick,et al.  Performance of neural network basecalling tools for Oxford Nanopore sequencing , 2019, Genome Biology.

[25]  Angela N. Brooks,et al.  Nanopore native RNA sequencing of a human poly(A) transcriptome , 2018, bioRxiv.

[26]  Mark Gerstein,et al.  GENCODE reference annotation for the human and mouse genomes , 2018, Nucleic Acids Res..

[27]  Angela N. Brooks,et al.  Full-length transcript characterization of SF3B1 mutation in chronic lymphocytic leukemia reveals downregulation of retained introns , 2018, Nature Communications.

[28]  S. Blankenberg,et al.  Non-coding RNAs in cardiovascular diseases: diagnostic and therapeutic perspectives. , 2018, European heart journal.

[29]  D. Metzger,et al.  Short-Term Fasting Reveals Amino Acid Metabolism as a Major Sex-Discriminating Factor in the Liver , 2018, Cell metabolism.

[30]  Qiang Liu,et al.  Direct , 2018, Definitions.

[31]  Paul D Hutchins,et al.  Caloric Restriction Engages Hepatic RNA Processing Mechanisms in Rhesus Monkeys. , 2018, Cell metabolism.

[32]  Wolfgang Huber,et al.  Alternative start and termination sites of transcription drive most transcript isoform differences across human tissues , 2017, Nucleic acids research.

[33]  Nicola J. Rinaldi,et al.  Genetic effects on gene expression across human tissues , 2017, Nature.

[34]  Nicholas A. Sinnott-Armstrong,et al.  An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues , 2017, Nature Methods.

[35]  Ge Gao,et al.  CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features , 2017, Nucleic Acids Res..

[36]  Anthony O. Olarerin-George,et al.  MetaPlotR: a Perl/R pipeline for plotting metagenes of nucleotide modifications and other transcriptomic sites , 2017, Bioinform..

[37]  R. Guigó,et al.  Comparative transcriptomics in human and mouse , 2017, Nature Reviews Genetics.

[38]  K. Lindblad-Toh,et al.  FEELnc: a tool for long non-coding RNA annotation and its application to the dog transcriptome , 2017, Nucleic acids research.

[39]  Chengqi Yi,et al.  Epitranscriptome sequencing technologies: decoding RNA modifications , 2016, Nature Methods.

[40]  Jun Zhu,et al.  Integrative Transcriptome Analyses of Metabolic Responses in Mice Define Pivotal LncRNA Metabolic Regulators. , 2016, Cell metabolism.

[41]  Daniel R. Garalde,et al.  Highly parallel direct RNA sequencing on an array of nanopores , 2016, Nature Methods.

[42]  Fei Wang,et al.  Transcriptome-wide distribution and function of RNA hydroxymethylcytosine , 2016, Science.

[43]  J. Schug,et al.  Integration of ATAC-seq and RNA-seq identifies human alpha cell and beta cell signature genes , 2016, Molecular metabolism.

[44]  Wen J. Li,et al.  Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation , 2015, Nucleic Acids Res..

[45]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[46]  S. J. Coleman,et al.  Determinants and implications of mRNA poly(A) tail size--does this protein make my tail look big? , 2014, Seminars in cell & developmental biology.

[47]  Aimin Li,et al.  PLEK: a tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme , 2014, BMC Bioinformatics.

[48]  Yi Zhao,et al.  Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts , 2013, Nucleic acids research.

[49]  J. Harrow,et al.  Transcriptome analysis of human tissues and cell lines reveals one dominant transcript per gene , 2013, Genome Biology.

[50]  J. Kocher,et al.  CPAT: Coding-Potential Assessment Tool using an alignment-free logistic regression model , 2013, Nucleic acids research.

[51]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[52]  O. Elemento,et al.  Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and near Stop Codons , 2012, Cell.

[53]  M. Kupiec,et al.  Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq , 2012, Nature.

[54]  P. Tontonoz,et al.  Transcriptional integration of metabolism by the nuclear sterol-activated receptors LXR and FXR , 2012, Nature Reviews Molecular Cell Biology.

[55]  M. Suematsu,et al.  The reconstituted 'humanized liver' in TK-NOG mice is mature and functional. , 2011, Biochemical and biophysical research communications.

[56]  C. Sapienza,et al.  Explaining inter-individual variability in phenotype: Is epigenetics up to the challenge? , 2010, Epigenetics.

[57]  T. Nilsen,et al.  Expansion of the eukaryotic proteome by alternative splicing , 2010, Nature.

[58]  Johan Auwerx,et al.  Targeting bile-acid signalling for metabolic diseases , 2008, Nature Reviews Drug Discovery.

[59]  S. Stamm,et al.  Function of Alternative Splicing , 2004 .

[60]  T. Leonardi,et al.  Nanopore RNA Sequencing Analysis. , 2021, Methods in molecular biology.

[61]  G. Pertea,et al.  GFF Utilities: GffRead and GffCompare. , 2020, F1000Research.

[62]  Michael Lehrke,et al.  The many faces of PPARgamma. , 2005, Cell.