Ancient RNA from Late Pleistocene permafrost and historical canids shows tissue-specific transcriptome survival

While sequencing ancient DNA (aDNA) from archaeological material is now commonplace, very few attempts to sequence ancient transcriptomes have been made, even from typically stable deposition environments such as permafrost. This is presumably due to assumptions that RNA completely degrades relatively quickly, particularly when dealing with autolytic, nuclease-rich mammalian tissues. However, given the recent successes in sequencing ancient RNA (aRNA) from various sources including plants and animals, we suspect that these assumptions may be incorrect or exaggerated. To challenge the underlying dogma, we generated shotgun RNA data from sources that might normally be dismissed for such study. Here, we present aRNA data generated from two historical wolf skins, and permafrost-preserved liver tissue of a 14,300-year-old Pleistocene canid. Not only is the latter the oldest RNA ever to be sequenced, but it also shows evidence of biologically relevant tissue specificity and close similarity to equivalent data derived from modern-day control tissue. Other hallmarks of RNA sequencing (RNA-seq) data such as exon-exon junction presence and high endogenous ribosomal RNA (rRNA) content confirms our data’s authenticity. By performing independent technical library replicates using two high-throughput sequencing platforms, we show not only that aRNA can survive for extended periods in mammalian tissues but also that it has potential for tissue identification. aRNA also has possible further potential, such as identifying in vivo genome activity and adaptation, when sequenced using this technology.

[1]  Anastacia M. Garcia,et al.  Optimization of phenol-chloroform RNA extraction , 2018, MethodsX.

[2]  B. Oliver,et al.  Comparison of normalization and differential expression analyses using RNA-Seq data from 726 individual Drosophila melanogaster , 2016, BMC Genomics.

[3]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[4]  S. Pääbo,et al.  Reconstructing the DNA Methylation Maps of the Neandertal and the Denisovan , 2014, Science.

[5]  M. Meyer,et al.  Quantifying and reducing spurious alignments for the analysis of ultra-short ancient DNA sequences , 2018, BMC Biology.

[6]  M. Gilbert,et al.  Long-term RNA persistence in postmortem contexts , 2013, Investigative Genetics.

[7]  Alexander F. Auch,et al.  MEGAN analysis of metagenomic data. , 2007, Genome research.

[8]  D. Shoham,et al.  Evidence of Influenza A Virus RNA in Siberian Lake Ice , 2006, Journal of Virology.

[9]  Christopher Bronk Ramsey,et al.  Recent and Planned Developments of the Program OxCal , 2013 .

[10]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[11]  D. Graham,et al.  The effect of sample type, temperature and RNAlater on the stability of avian influenza virus RNA. , 2008, Journal of virological methods.

[12]  Patricia P. Chan,et al.  GtRNAdb: a database of transfer RNA genes detected in genomic sequence , 2008, Nucleic Acids Res..

[13]  Sarah A. Palmer,et al.  Small RNA Activity in Archeological Barley Shows Novel Germination Inhibition in Response to Environment , 2017, Molecular biology and evolution.

[14]  T. Korneliussen,et al.  Ancient genomics , 2015, Philosophical Transactions of the Royal Society B: Biological Sciences.

[15]  E. Willerslev,et al.  Isolation of nucleic acids and cultures from fossil ice and permafrost. , 2004, Trends in ecology & evolution.

[16]  Peter A. Jones,et al.  Genome-wide nucleosome map and cytosine methylation levels of an ancient human genome , 2014, Genome research.

[17]  Sarah A. Palmer,et al.  Archaeogenetic Evidence of Ancient Nubian Barley Evolution from Six to Two-Row Indicates Local Adaptation , 2009, PloS one.

[18]  Eran Elinav,et al.  Use of Metatranscriptomics in Microbiome Research , 2016, Bioinformatics and biology insights.

[19]  S. Tirrell,et al.  The beagle dog MicroRNA tissue atlas: identifying translatable biomarkers of organ toxicity , 2016, BMC Genomics.

[20]  A. Westermann,et al.  Host-Pathogen Transcriptomics by Dual RNA-Seq. , 2018, Methods in molecular biology.

[21]  H. Ginsberg,et al.  Hepatic Synthesis of Lipoproteins and Apolipoproteins , 1992, Seminars in liver disease.

[22]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[23]  P. Reynier,et al.  Impact of various handling and storage conditions on quantitative detection of hepatitis C virus RNA. , 1996, Journal of hepatology.

[24]  J. Olsen,et al.  Proteomic profiling of archaeological human bone , 2017, Royal Society Open Science.

[25]  J. Khan,et al.  A Compendium of Canine Normal Tissue Gene Expression , 2011, PloS one.

[26]  Simon Kasif,et al.  GC/AT-content spikes as genomic punctuation marks. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[27]  C. Ramsey,et al.  Deep Sequencing of RNA from Ancient Maize Kernels , 2013, PloS one.

[28]  Paul Francis Harrison Varistran: Anscombe's variance stabilizing transformation for RNA-seq gene expression data , 2017, J. Open Source Softw..

[29]  M. Thomas P. Gilbert,et al.  mapDamage: testing for damage patterns in ancient DNA sequences , 2011, Bioinform..

[30]  Oliver Smith,et al.  A complete ancient RNA genome: identification, reconstruction and evolutionary history of archaeological Barley Stripe Mosaic Virus , 2014, Scientific Reports.

[31]  Stewart Shuman,et al.  Structure-Function Analysis of T4 RNA Ligase 2* , 2003, The Journal of Biological Chemistry.

[32]  Jesper V Olsen,et al.  Pathogens and host immunity in the ancient human oral cavity , 2014, Nature Genetics.

[33]  Gerald Zon,et al.  Small RNA Library Preparation Method for Next-Generation Sequencing Using Chemical Modifications to Prevent Adapter Dimer Formation , 2016, PloS one.

[34]  S. Clarke,et al.  Exceptional seed longevity and robust growth: ancient Sacred Lotus from China , 1995 .

[35]  D. Carraro,et al.  Biobanking practice: RNA storage at low concentration affects integrity. , 2014, Biopreservation and biobanking.

[36]  Xinhua Xiao,et al.  Downregulation of zinc-α2-glycoprotein in adipose tissue and liver of obese ob/ob mice and by tumour necrosis factor-α in adipocytes , 2009, The Journal of endocrinology.

[37]  W. Jahnen-Dechent,et al.  Tissue distribution and activity testing suggest a similar but not identical function of fetuin-B and fetuin-A. , 2003, The Biochemical journal.

[38]  Jun Wang,et al.  Genomic methylation patterns in archaeological barley show de-methylation as a time-dependent diagenetic process , 2014, Scientific Reports.

[39]  Michael Hahsler,et al.  Getting Things in Order: An Introduction to the R Package seriation , 2008 .

[40]  J. McLachlan,et al.  Ancient DNA from lake sediments: Bridging the gap between paleoecology and genetics , 2011, BMC Evolutionary Biology.

[41]  Marcel Martin Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[42]  C. Richards,et al.  The utility of aged seeds in DNA banks , 2006, Seed Science Research.

[43]  M. Esumi,et al.  Pitfalls of DNA Quantification Using DNA-Binding Fluorescent Dyes and Suggested Solutions , 2016, PloS one.

[44]  M. Pepys,et al.  Human plasma fibrinogen is synthesized in the liver. , 2007, Blood.

[45]  M. Collins,et al.  Preservation of the metaproteome: variability of protein preservation in ancient dental calculus , 2017, Science and technology of archaeological research.

[46]  P. Chomczyński,et al.  The single-step method of RNA isolation by acid guanidinium thiocyanate–phenol–chloroform extraction: twenty-something years on , 2006, Nature Protocols.

[47]  T. Sicheritz-Pontén,et al.  Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing , 2017, GigaScience.

[48]  F. J. Anscombe,et al.  THE TRANSFORMATION OF POISSON, BINOMIAL AND NEGATIVE-BINOMIAL DATA , 1948 .

[49]  T. Ng,et al.  Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch , 2014, Proceedings of the National Academy of Sciences.

[50]  R. Nielsen,et al.  Ancient DNA chronology within sediment deposits: are paleobiological reconstructions possible and is DNA leaching a factor? , 2007, Molecular biology and evolution.

[51]  Pelin Yilmaz,et al.  The SILVA ribosomal RNA gene database project: improved data processing and web-based tools , 2012, Nucleic Acids Res..

[52]  Vikas Bansal,et al.  A computational method for estimating the PCR duplication rate in DNA and RNA-seq experiments , 2017, BMC Bioinformatics.

[53]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[54]  A. Franke,et al.  miRNAs in Ancient Tissue Specimens of the Tyrolean Iceman , 2016, Molecular biology and evolution.