Population scale mapping of transposable element diversity reveals links to gene regulation and epigenomic variation
Abstract
Variation in the presence or absence of transposable elements (TEs) is a major source of genetic variation between individuals. Here, we identified 23,095 TE presence/absence variants between 216 Arabidopsis accessions. Most TE variants were rare, and we find these rare variants associated with local extremes of gene expression and DNA methylation levels within the population. Of the common alleles identified, two thirds were not in linkage disequilibrium with nearby SNPs, implicating these variants as a source of novel genetic diversity. Many common TE variants were associated with significantly altered expression of nearby genes, and a major fraction of inter-accession DNA methylation differences were associated with nearby TE insertions. Overall, this demonstrates that TE variants are a rich source of genetic diversity that likely plays an important role in facilitating epigenomic and transcriptional differences between individuals, and indicates a strong genetic basis for epigenetic variation.
Data availability
-
Population scale mapping of transposable element diversity reveals links to gene regulation and epigenomic variationAvailable at Dryad Digital Repository under a CC0 Public Domain Dedication.
-
Patterns of population epigenomic diversityPublicly available at the NCBI Short Read Archive (accession no: SRA012474).
-
Patterns of population epigenomic diversityPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE43858).
-
Patterns of population epigenomic diversityPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE43857).
-
Reference-guided assembly of four diverse Arabidopsis thaliana genomesAvailable at the 1001 Genomes Project.
-
Environmentally responsive genome-wide accumulation of de novo Arabidopsis thaliana mutations and epimutationsPublicly available at the NCBI Short Read Archive (accession no: SRP045804).
-
Mapping and Dynamics of Regulatory DNA and Transcription Factor Networks in A. thalianaPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE53322).
Article and author information
Author details
Funding
Centre of Excellence in Plant Energy Biology, Australian Research Council (CE140100008)
- Tim Stuart
- Steven Eichten
- Jonathan Cahn
- Yuliya Karpievitch
- Justin Borevitz
- Ryan Lister
Australian Research Council
- Tim Stuart
- Steven Eichten
- Jonathan Cahn
- Yuliya Karpievitch
- Justin Borevitz
- Ryan Lister
Sylvia and Charles Viertel Charitable Foundation
- Ryan Lister
Australian Research Council (FT120100862)
- Ryan Lister
Australian Research Council (DE150101206)
- Steven Eichten
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Daniel Zilberman, University of California, Berkeley, United States
Version history
- Received: August 18, 2016
- Accepted: December 1, 2016
- Accepted Manuscript published: December 2, 2016 (version 1)
- Version of Record published: December 19, 2016 (version 2)
Copyright
© 2016, Stuart et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,714
- views
-
- 1,200
- downloads
-
- 178
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Chromosomes and Gene Expression
Recent findings indicate that the translation elongation rate influences mRNA stability. One of the factors that has been implicated in this link between mRNA decay and translation speed is the yeast DEAD-box helicase Dhh1p. Here, we demonstrated that the human ortholog of Dhh1p, DDX6, triggers the deadenylation-dependent decay of inefficiently translated mRNAs in human cells. DDX6 interacts with the ribosome through the Phe-Asp-Phe (FDF) motif in its RecA2 domain. Furthermore, RecA2-mediated interactions and ATPase activity are both required for DDX6 to destabilize inefficiently translated mRNAs. Using ribosome profiling and RNA sequencing, we identified two classes of endogenous mRNAs that are regulated in a DDX6-dependent manner. The identified targets are either translationally regulated or regulated at the steady-state-level and either exhibit signatures of poor overall translation or of locally reduced ribosome translocation rates. Transferring the identified sequence stretches into a reporter mRNA caused translation- and DDX6-dependent degradation of the reporter mRNA. In summary, these results identify DDX6 as a crucial regulator of mRNA translation and decay triggered by slow ribosome movement and provide insights into the mechanism by which DDX6 destabilizes inefficiently translated mRNAs.
-
- Chromosomes and Gene Expression
Tardigrades are microscopic animals renowned for their ability to withstand extreme conditions, including high doses of ionizing radiation (IR). To better understand their radio-resistance, we first characterized induction and repair of DNA double- and single-strand breaks after exposure to IR in the model species Hypsibius exemplaris. Importantly, we found that the rate of single-strand breaks induced was roughly equivalent to that in human cells, suggesting that DNA repair plays a predominant role in tardigrades’ radio-resistance. To identify novel tardigrade-specific genes involved, we next conducted a comparative transcriptomics analysis across three different species. In all three species, many DNA repair genes were among the most strongly overexpressed genes alongside a novel tardigrade-specific gene, which we named Tardigrade DNA damage Response 1 (TDR1). We found that TDR1 protein interacts with DNA and forms aggregates at high concentration suggesting it may condensate DNA and preserve chromosome organization until DNA repair is accomplished. Remarkably, when expressed in human cells, TDR1 improved resistance to Bleomycin, a radiomimetic drug. Based on these findings, we propose that TDR1 is a novel tardigrade-specific gene conferring resistance to IR. Our study sheds light on mechanisms of DNA repair helping cope with high levels of DNA damage inflicted by IR.