Abstract:BackgroundInstead of testing predefined hypotheses, the goal of exploratory data analysis (EDA) is to find what data can tell us. Following this strategy, we re-analyzed a large body of genomic data to study the complex gene regulation in mouse pre-implantation development (PD).ResultsStarting with a single-cell RNA-seq dataset consisting of 259 mouse embryonic cells derived from zygote to blastocyst stages, we reconstructed the temporal and spatial gene expression pattern during PD. The dynamics of gene expression can be partially explained by the enrichment of transposable elements in gene promoters and the similarity of expression profiles with those of corresponding transposons. Long Terminal Repeats (LTRs) are associated with transient, strong induction of many nearby genes at the 2-4 cell stages, probably by providing binding sites for Obox and other homeobox factors. B1 and B2 SINEs (Short Interspersed Nuclear Elements) are correlated with the upregulation of thousands of nearby genes during zygotic genome activation. Such enhancer-like effects are also found for human Alu and bovine tRNA SINEs. SINEs also seem to be predictive of gene expression in embryonic stem cells (ESCs), raising the possibility that they may also be involved in regulating pluripotency. We also identified many potential transcription factors underlying PD and discussed the evolutionary necessity of transposons in enhancing genetic diversity, especially for species with longer generation time. ConclusionsTogether with other recent studies, our results provide further evidence that many transposable elements may play a role in establishing the expression landscape in early embryos. It also demonstrates that exploratory bioinformatics investigation can pinpoint developmental pathways for further study, and serve as a strategy to generate novel insights from big genomic data.

Biological function of unannotated transcription during the early development of Drosophila melanogaster

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

Transcription start site profiling uncovers divergent transcription and enhancer-associated RNAs in Drosophila melanogaster

Diversity and dynamics of the Drosophila transcriptome

Identification of Functional Elements and Regulatory Circuits by Drosophila Modencode

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq.

Comprehensive Identification and Alternative Splicing of Microexons in Drosophila

Genome-wide identification and developmental expression profiling of long noncoding RNAs during Drosophila metamorphosis

Organogenetic transcriptomes of the Drosophila embryo at single cell resolution

Origin of New Genes and Source for N-terminal Domain of the Chimerical Gene, Jingwei, in Drosophila.

Detecting Novel Low-Abundant Transcripts In Drosophila

Exploratory bioinformatics investigation reveals importance of “junk” DNA in early embryo development

Evolution of Alternative Splicing in Newly Evolved Genes of Drosophila.

Ultra-deep Sequencing of Ribosome-Associated Poly-Adenylated RNA in Early Drosophila Embryos Reveals Hundreds of Conserved Translated Sorfs.

Extensive Divergence of Transcription Factor Binding in Drosophila Embryos with Highly Conserved Gene Expression

The Transcriptional Diversity of 25 Drosophila Cell Lines

Precise Annotation of Drosophila Mitochondrial Genomes Leads to Insights into AT-rich Regions

The origin and structural evolution of de novo genes in Drosophila

Identification and quantification of transposable element transcripts using Long-Read RNA-seq in Drosophila germline tissues

On the Origin of New Genes in Drosophila.

Variant functional assessment in Drosophila by overexpression: what can we learn?