Mapping the mouse Allelome reveals tissue-specific regulation of allelic expression.

To determine the dynamics of allelic-specific expression during mouse development, we analyzed RNA-seq data from 23 F1 tissues from different developmental stages, including 19 female tissues allowing X chromosome inactivation (XCI) escapers to also be detected. We demonstrate that allelic expression arising from genetic or epigenetic differences is highly tissue-specific. We find that tissue-specific strain-biased gene expression may be regulated by tissue-specific enhancers or by post-transcriptional differences in stability between the alleles. We also find that escape from X-inactivation is tissue-specific, with leg muscle showing an unexpectedly high rate of XCI escapers. By surveying a range of tissues during development, and performing extensive validation, we are able to provide a high confidence list of mouse imprinted genes including 18 novel genes. This shows that cluster size varies dynamically during development and can be substantially larger than previously thought, with the Igf2r cluster extending over 10 Mb in placenta.

Long non-coding RNAs display higher natural expression variation than protein-coding genes in healthy humans.

BACKGROUND: Long non-coding RNAs (lncRNAs) are increasingly implicated as gene regulators and may ultimately be more numerous than protein-coding genes in the human genome. Despite large numbers of reported lncRNAs, reference annotations are likely incomplete due to their lower and tighter tissue-specific expression compared to mRNAs. An unexplored factor potentially confounding lncRNA identification is inter-individual expression variability. Here, we characterize lncRNA natural expression variability in human primary granulocytes. RESULTS: We annotate granulocyte lncRNAs and mRNAs in RNA-seq data from 10 healthy individuals, identifying multiple lncRNAs absent from reference annotations, and use this to investigate three known features (higher tissue-specificity, lower expression, and reduced splicing efficiency) of lncRNAs relative to mRNAs. Expression variability was examined in seven individuals sampled three times at 1- or more than 1-month intervals. We show that lncRNAs display significantly more inter-individual expression variability compared to mRNAs. We confirm this finding in two independent human datasets by analyzing multiple tissues from the GTEx project and lymphoblastoid cell lines from the GEUVADIS project. Using the latter dataset we also show that including more human donors into the transcriptome annotation pipeline allows identification of an increasing number of lncRNAs, but minimally affects mRNA gene number. CONCLUSIONS: A comprehensive annotation of lncRNAs is known to require an approach that is sensitive to low and tight tissue-specific expression. Here we show that increased inter-individual expression variability is an additional general lncRNA feature to consider when creating a comprehensive annotation of human lncRNAs or proposing their use as prognostic or disease markers.

A human haploid gene trap collection to study lncRNAs with unusual RNA biology.

Many thousand long non-coding (lnc) RNAs are mapped in the human genome. Time consuming studies using reverse genetic approaches by post-transcriptional knock-down or genetic modification of the locus demonstrated diverse biological functions for a few of these transcripts. The Human Gene Trap Mutant Collection in haploid KBM7 cells is a ready-to-use tool for studying protein-coding gene function. As lncRNAs show remarkable differences in RNA biology compared to protein-coding genes, it is unclear if this gene trap collection is useful for functional analysis of lncRNAs. Here we use the uncharacterized LOC100288798 lncRNA as a model to answer this question. Using public RNA-seq data we show that LOC100288798 is ubiquitously expressed, but inefficiently spliced. The minor spliced LOC100288798 isoforms are exported to the cytoplasm, whereas the major unspliced isoform is nuclear localized. This shows that LOC100288798 RNA biology differs markedly from typical mRNAs. De novo assembly from RNA-seq data suggests that LOC100288798 extends 289kb beyond its annotated 3' end and overlaps the downstream SLC38A4 gene. Three cell lines with independent gene trap insertions in LOC100288798 were available from the KBM7 gene trap collection. RT-qPCR and RNA-seq confirmed successful lncRNA truncation and its extended length. Expression analysis from RNA-seq data shows significant deregulation of 41 protein-coding genes upon LOC100288798 truncation. Our data shows that gene trap collections in human haploid cell lines are useful tools to study lncRNAs, and identifies the previously uncharacterized LOC100288798 as a potential gene regulator.

Allelome.PRO, a pipeline to define allele-specific genomic features from high-throughput sequencing data.

Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the 'allelome' by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Imprinted expression in cystic embryoid bodies shows an embryonic and not an extra-embryonic pattern.

A large subset of mammalian imprinted genes show extra-embryonic lineage (EXEL) specific imprinted expression that is restricted to placental trophectoderm lineages and to visceral yolk sac endoderm (ysE). Isolated ysE provides a homogenous in vivo model of a mid-gestation extra-embryonic tissue to examine the mechanism of EXEL-specific imprinted gene silencing, but an in vitro model of ysE to facilitate more rapid and cost-effective experiments is not available. Reports indicate that ES cells differentiated into cystic embryoid bodies (EBs) contain ysE, so here we investigate if cystic EBs model ysE imprinted expression. The imprinted expression pattern of cystic EBs is shown to resemble fetal liver and not ysE. To investigate the reason for this we characterized the methylome and transcriptome of cystic EBs in comparison to fetal liver and ysE, by whole genome bisulphite sequencing and RNA-seq. Cystic EBs show a fetal liver pattern of global hypermethylation and low expression of repeats, while ysE shows global hypomethylation and high expression of IAPEz retroviral repeats, as reported for placenta. Transcriptome analysis confirmed that cystic EBs are more similar to fetal liver than ysE and express markers of early embryonic endoderm. Genome-wide analysis shows that ysE shares epigenetic and repeat expression features with placenta. Contrary to previous reports, we show that cystic EBs do not contain ysE, but are more similar to the embryonic endoderm of fetal liver. This explains why cystic EBs reproduce the imprinted expression seen in the embryo but not that seen in the ysE.

Considerations when investigating lncRNA function in vivo.

Although a small number of the vast array of animal long non-coding RNAs (lncRNAs) have known effects on cellular processes examined in vitro, the extent of their contributions to normal cell processes throughout development, differentiation and disease for the most part remains less clear. Phenotypes arising from deletion of an entire genomic locus cannot be unequivocally attributed either to the loss of the lncRNA per se or to the associated loss of other overlapping DNA regulatory elements. The distinction between cis- or trans-effects is also often problematic. We discuss the advantages and challenges associated with the current techniques for studying the in vivo function of lncRNAs in the light of different models of lncRNA molecular mechanism, and reflect on the design of experiments to mutate lncRNA loci. These considerations should assist in the further investigation of these transcriptional products of the genome.

Genomic imprinting in mammals.

Genomic imprinting affects a subset of genes in mammals and results in a monoallelic, parental-specific expression pattern. Most of these genes are located in clusters that are regulated through the use of insulators or long noncoding RNAs (lncRNAs). To distinguish the parental alleles, imprinted genes are epigenetically marked in gametes at imprinting control elements through the use of DNA methylation at the very least. Imprinted gene expression is subsequently conferred through lncRNAs, histone modifications, insulators, and higher-order chromatin structure. Such imprints are maintained after fertilization through these mechanisms despite extensive reprogramming of the mammalian genome. Genomic imprinting is an excellent model for understanding mammalian epigenetic regulation.

Gene regulation by the act of long non-coding RNA transcription.

Long non-protein-coding RNAs (lncRNAs) are proposed to be the largest transcript class in the mouse and human transcriptomes. Two important questions are whether all lncRNAs are functional and how they could exert a function. Several lncRNAs have been shown to function through their product, but this is not the only possible mode of action. In this review we focus on a role for the process of lncRNA transcription, independent of the lncRNA product, in regulating protein-coding-gene activity in cis. We discuss examples where lncRNA transcription leads to gene silencing or activation, and describe strategies to determine if the lncRNA product or its transcription causes the regulatory effect.

Imprinted silencing is extended over broad chromosomal domains in mouse extra-embryonic lineages.

Gene silencing in imprinted gene clusters is established by an epigenetic initiator that is often a long non-coding (lnc) RNA. The clustered organization of known imprinted genes indicates that the initiator extends imprinted silencing over broader chromosomal domains in extra-embryonic lineages compared to the embryo. We propose that extension of imprinted gene clusters may result from known epigenetic differences between extra-embryonic and embryonic lineages that alter the behavior of epigenetic initiators. New RNA sequencing technology will enable the full extent of imprinted silencing in embryonic and extra-embryonic lineages to be defined, but appropriate analysis and cell systems are required, which we define here based on a review of recent studies.

Imprinted Igf2r silencing depends on continuous Airn lncRNA expression and is not restricted to a developmental window.

The imprinted Airn macro long non-coding (lnc) RNA is an established example of a cis-silencing lncRNA. Airn expression is necessary to initiate paternal-specific silencing of the Igf2r gene, which is followed by gain of a somatic DNA methylation imprint on the silent Igf2r promoter. However, the developmental requirements for Airn initiation of Igf2r silencing and the role of Airn or DNA methylation in maintaining stable Igf2r repression have not been investigated. Here, we use inducible systems to control Airn expression during mouse embryonic stem cell (ESC) differentiation. By turning Airn expression off during ESC differentiation, we show that continuous Airn expression is needed to maintain Igf2r silencing, but only until the paternal Igf2r promoter is methylated. By conditionally turning Airn expression on, we show that Airn initiation of Igf2r silencing is not limited to one developmental 'window of opportunity' and can be maintained in the absence of DNA methylation. Together, this study shows that Airn expression is both necessary and sufficient to silence Igf2r throughout ESC differentiation and that the somatic methylation imprint, although not required to initiate or maintain silencing, adds a secondary layer of repressive epigenetic information.

Airn transcriptional overlap, but not its lncRNA products, induces imprinted Igf2r silencing.

Mammalian imprinted genes often cluster with long noncoding (lnc) RNAs. Three lncRNAs that induce parental-specific silencing show hallmarks indicating that their transcription is more important than their product. To test whether Airn transcription or product silences the Igf2r gene, we shortened the endogenous lncRNA to different lengths. The results excluded a role for spliced and unspliced Airn lncRNA products and for Airn nuclear size and location in silencing Igf2r. Instead, silencing only required Airn transcriptional overlap of the Igf2r promoter, which interferes with RNA polymerase II recruitment in the absence of repressive chromatin. Such a repressor function for lncRNA transcriptional overlap reveals a gene silencing mechanism that may be widespread in the mammalian genome, given the abundance of lncRNA transcripts.

Macro lncRNAs: a new layer of cis-regulatory information in the mammalian genome.

In the past ten years, long non-protein-coding RNAs (lncRNAs) have been shown to comprise a major part of the mammalian transcriptome and are predicted from their highly specific expression patterns, to play a role in regulating protein-coding gene expression in development and disease. Many lncRNAs particularly those lying in imprinted clusters, are predominantly unspliced "macro" transcripts that can also show a low level of splicing, and both the unspliced and spliced transcript have the potential to be functional. Three known imprinted macro lncRNAs have been shown to silence from three to ten genes in cis in imprinted gene clusters. We review here the potential for functional macro lncRNAs, defined here as "inefficiently-spliced lncRNAs" to play a wider cis-regulatory role in the mammalian genome. This potential has been underestimated by the inability of current RNA-seq technology to annotate unspliced macro lncRNAs.

Mechanisms of long range silencing by imprinted macro non-coding RNAs.

Non-coding (nc) RNA silencing of imprinted genes in extra-embryonic tissues provides a good model for understanding an underexamined aspect of gene regulation by macro or long ncRNAs, that is their action as long-range cis-silencers. Numerous long intergenic ncRNAs (lincRNAs) have been recently discovered that are thought to regulate gene expression, some of which have been associated with disease. The few shown to regulate protein-coding genes are suggested to act by targeting repressive or active chromatin marks. Correlative evidence also indicates that imprinted macro ncRNAs cause long-range cis-silencing in placenta by targeting repressive histone modifications to imprinted promoters. It is timely, however, to consider alternative explanations consistent with the published data, whereby transcription alone could cause gene silencing at a distance.

A downstream CpG island controls transcript initiation and elongation and the methylation state of the imprinted Airn macro ncRNA promoter.

A CpG island (CGI) lies at the 5' end of the Airn macro non-protein-coding (nc) RNA that represses the flanking Igf2r promoter in cis on paternally inherited chromosomes. In addition to being modified on maternally inherited chromosomes by a DNA methylation imprint, the Airn CGI shows two unusual organization features: its position immediately downstream of the Airn promoter and transcription start site and a series of tandem direct repeats (TDRs) occupying its second half. The physical separation of the Airn promoter from the CGI provides a model to investigate if the CGI plays distinct transcriptional and epigenetic roles. We used homologous recombination to generate embryonic stem cells carrying deletions at the endogenous locus of the entire CGI or just the TDRs. The deleted Airn alleles were analyzed by using an ES cell imprinting model that recapitulates the onset of Igf2r imprinted expression in embryonic development or by using knock-out mice. The results show that the CGI is required for efficient Airn initiation and to maintain the unmethylated state of the Airn promoter, which are both necessary for Igf2r repression on the paternal chromosome. The TDRs occupying the second half of the CGI play a minor role in Airn transcriptional elongation or processivity, but are essential for methylation on the maternal Airn promoter that is necessary for Igf2r to be expressed from this chromosome. Together the data indicate the existence of a class of regulatory CGIs in the mammalian genome that act downstream of the promoter and transcription start.

An RNA-Seq strategy to detect the complete coding and non-coding transcriptome including full-length imprinted macro ncRNAs.

Imprinted macro non-protein-coding (nc) RNAs are cis-repressor transcripts that silence multiple genes in at least three imprinted gene clusters in the mouse genome. Similar macro or long ncRNAs are abundant in the mammalian genome. Here we present the full coding and non-coding transcriptome of two mouse tissues: differentiated ES cells and fetal head using an optimized RNA-Seq strategy. The data produced is highly reproducible in different sequencing locations and is able to detect the full length of imprinted macro ncRNAs such as Airn and Kcnq1ot1, whose length ranges between 80-118 kb. Transcripts show a more uniform read coverage when RNA is fragmented with RNA hydrolysis compared with cDNA fragmentation by shearing. Irrespective of the fragmentation method, all coding and non-coding transcripts longer than 8 kb show a gradual loss of sequencing tags towards the 3' end. Comparisons to published RNA-Seq datasets show that the strategy presented here is more efficient in detecting known functional imprinted macro ncRNAs and also indicate that standardization of RNA preparation protocols would increase the comparability of the transcriptome between different RNA-Seq datasets.

Genomic imprinting: a mammalian epigenetic discovery model.

Genomic imprinting is an epigenetic process leading to parental-specific expression of one to two percent of mammalian genes that offers one of the best model systems for a molecular analysis of epigenetic regulation in development and disease. In the twenty years since the first imprinted gene was identified, this model has had a significant impact on decoding epigenetic information in mammals. So far it has led to the discovery of long-range cis-acting control elements whose epigenetic state regulates small clusters of genes and of unusual macro noncoding RNAs (ncRNAs) that directly repress genes in cis, and critically, it has demonstrated that one biological role of DNA methylation is to allow expression of genes normally repressed by default. This review describes the progress in understanding how imprinted protein-coding genes are silenced; in particular, it focuses on the role of macro ncRNAs that have broad relevance as a potential new layer of regulatory information in the mammalian genome.

Extra-embryonic-specific imprinted expression is restricted to defined lineages in the post-implantation embryo.

A subset of imprinted genes in the mouse have been reported to show imprinted expression that is restricted to the placenta, a short-lived extra-embryonic organ. Notably, these so-called "placental-specific" imprinted genes are expressed from both parental alleles in embryo and adult tissues. The placenta is an embryonic-derived organ that is closely associated with maternal tissue, and as a consequence, maternal contamination can be mistaken for maternal-specific imprinted expression. The complexity of the placenta, which arises from multiple embryonic lineages, poses additional problems in accurately assessing allele-specific repressive epigenetic modifications in genes that also show lineage-specific silencing in this organ. These problems require that extra evidence be obtained to support the imprinted status of genes whose imprinted expression is restricted to the placenta. We show here that the extra-embryonic visceral yolk sac (VYS), a nutritive membrane surrounding the developing embryo, shows a similar "extra-embryonic-lineage-specific" pattern of imprinted expression. We present an improved enzymatic technique for separating the bilaminar VYS and show that this pattern of imprinted expression is restricted to the endoderm layer. Finally, we show that VYS "extra-embryonic-lineage-specific" imprinted expression is regulated by DNA methylation in a similar manner as shown for genes showing multi-lineage imprinted expression in extra-embryonic, embryonic, and adult tissues. These results show that the VYS is an improved model for studying the epigenetic mechanisms regulating extra-embryonic-lineage-specific imprinted expression.

Developmental control of imprinted expression by macro non-coding RNAs.

Genomic imprinting is a developmentally regulated epigenetic phenomenon. The majority of imprinted genes only show parent-of-origin specific expression in a subset of tissues or at defined developmental stages. In some cases, imprinted expression is controlled by an imprinted macro non-coding RNA (ncRNA) whose expression pattern and repressive activity does not necessarily correlate with that of the genes whose imprinted expression it controls. This suggests that developmentally regulated factors other than the macro ncRNA are involved in establishing or maintaining imprinted expression. Here, we review how macro ncRNAs control imprinted expression during development and differentiation and consider how this impacts on target choice in epigenetic therapy.

Genomic imprinting-an epigenetic gene-regulatory model.

Epigenetic mechanisms (Box 1) are considered to play major gene-regulatory roles in development, differentiation and disease. However, the relative importance of epigenetics in defining the mammalian transcriptome in normal and disease states is unknown. The mammalian genome contains only a few model systems where epigenetic gene regulation has been shown to play a major role in transcriptional control. These model systems are important not only to investigate the biological function of known epigenetic modifications but also to identify new and unexpected epigenetic mechanisms in the mammalian genome. Here we review recent progress in understanding how epigenetic mechanisms control imprinted gene expression.

The function of non-coding RNAs in genomic imprinting.

Non-coding RNAs (ncRNAs) that regulate gene expression in cis or in trans are a shared feature of prokaryotic and eukaryotic genomes. In mammals, cis-acting functions are associated with macro ncRNAs, which can be several hundred thousand nucleotides long. Imprinted ncRNAs are well-studied macro ncRNAs that have cis-regulatory effects on multiple flanking genes. Recent advances indicate that they employ different downstream mechanisms to regulate gene expression in embryonic and placental tissues. A better understanding of these downstream mechanisms will help to improve our general understanding of the function of ncRNAs throughout the genome.