From genotype to phenotype: genetics of mammalian long non-coding RNAs in vivo.

Genome-wide sequencing has led to the discovery of thousands of long non-coding RNA (lncRNA) loci in the human genome, but evidence of functional significance has remained controversial for many lncRNAs. Genetically engineered model organisms are considered the gold standard for linking genotype to phenotype. Recent advances in CRISPR-Cas genome editing have led to a rapid increase in the use of mouse models to more readily survey lncRNAs for functional significance. Here, we review strategies to investigate the physiological relevance of lncRNA loci by highlighting studies that have used genetic mouse models to reveal key in vivo roles for lncRNAs, from fertility to brain development. We illustrate how an investigative approach, starting with whole-gene deletion followed by transcription termination and/or transgene rescue strategies, can provide definitive evidence for the in vivo function of mammalian lncRNAs.

Diverse epigenetic mechanisms maintain parental imprints within the embryonic and extraembryonic lineages.

Genomic imprinting and X chromosome inactivation (XCI) require epigenetic mechanisms to encode allele-specific expression, but how these specific tasks are accomplished at single loci or across chromosomal scales remains incompletely understood. Here, we systematically disrupt essential epigenetic pathways within polymorphic embryos in order to examine canonical and non-canonical genomic imprinting as well as XCI. We find that DNA methylation and Polycomb group repressors are indispensable for autosomal imprinting, albeit at distinct gene sets. Moreover, the extraembryonic ectoderm relies on a broader spectrum of imprinting mechanisms, including non-canonical targeting of maternal endogenous retrovirus (ERV)-driven promoters by the H3K9 methyltransferase G9a. We further identify Polycomb-dependent and -independent gene clusters on the imprinted X chromosome, which appear to reflect distinct domains of Xist-mediated suppression. From our data, we assemble a comprehensive inventory of the epigenetic pathways that maintain parent-specific imprinting in eutherian mammals, including an expanded view of the placental lineage.

Correction: The Airn lncRNA does not require any DNA elements within its locus to silence distant imprinted genes.

[This corrects the article DOI: 10.1371/journal.pgen.1008268.].

Cis and trans effects differentially contribute to the evolution of promoters and enhancers.

BACKGROUND: Gene expression differences between species are driven by both cis and trans effects. Whereas cis effects are caused by genetic variants located on the same DNA molecule as the target gene, trans effects are due to genetic variants that affect diffusible elements. Previous studies have mostly assessed the impact of cis and trans effects at the gene level. However, how cis and trans effects differentially impact regulatory elements such as enhancers and promoters remains poorly understood. Here, we use massively parallel reporter assays to directly measure the transcriptional outputs of thousands of individual regulatory elements in embryonic stem cells and measure cis and trans effects between human and mouse. RESULTS: Our approach reveals that cis effects are widespread across transcribed regulatory elements, and the strongest cis effects are associated with the disruption of motifs recognized by strong transcriptional activators. Conversely, we find that trans effects are rare but stronger in enhancers than promoters and are associated with a subset of transcription factors that are differentially expressed between human and mouse. While we find that cis-trans compensation is common within promoters, we do not see evidence of widespread cis-trans compensation at enhancers. Cis-trans compensation is inversely correlated with enhancer redundancy, suggesting that such compensation may often occur across multiple enhancers. CONCLUSIONS: Our results highlight differences in the mode of evolution between promoters and enhancers in complex mammalian genomes and indicate that studying the evolution of individual regulatory elements is pivotal to understand the tempo and mode of gene expression evolution.

In vivo Firre and Dxz4 deletion elucidates roles for autosomal gene regulation.

Recent evidence has determined that the conserved X chromosome mega-structures controlled by the Firre and Dxz4 loci are not required for X chromosome inactivation (XCI) in cell lines. Here, we examined the in vivo contribution of these loci by generating mice carrying a single or double deletion of Firre and Dxz4. We found that these mutants are viable, fertile and show no defect in random or imprinted XCI. However, the lack of these elements results in many dysregulated genes on autosomes in an organ-specific manner. By comparing the dysregulated genes between the single and double deletion, we identified superloop, megadomain, and Firre locus-dependent gene sets. The largest transcriptional effect was observed in all strains lacking the Firre locus, indicating that this locus is the main driver for these autosomal expression signatures. Collectively, these findings suggest that these X-linked loci are involved in autosomal gene regulation rather than XCI biology.

The Airn lncRNA does not require any DNA elements within its locus to silence distant imprinted genes.

Long non-coding (lnc) RNAs are numerous and found throughout the mammalian genome, and many are thought to be involved in the regulation of gene expression. However, the majority remain relatively uncharacterised and of uncertain function making the use of model systems to uncover their mode of action valuable. Imprinted lncRNAs target and recruit epigenetic silencing factors to a cluster of imprinted genes on the same chromosome, making them one of the best characterized lncRNAs for silencing distant genes in cis. In this study we examined silencing of the distant imprinted gene Slc22a3 by the lncRNA Airn in the Igf2r imprinted cluster in mouse. Previously we proposed that imprinted lncRNAs may silence distant imprinted genes by disrupting promoter-enhancer interactions by being transcribed through the enhancer, which we called the enhancer interference hypothesis. Here we tested this hypothesis by first using allele-specific chromosome conformation capture (3C) to detect interactions between the Slc22a3 promoter and the locus of the Airn lncRNA that silences it on the paternal chromosome. In agreement with the model, we found interactions enriched on the maternal allele across the entire Airn gene consistent with multiple enhancer-promoter interactions. Therefore, to test the enhancer interference hypothesis we devised an approach to delete the entire Airn gene. However, the deletion showed that there are no essential enhancers for Slc22a2, Pde10a and Slc22a3 within the Airn gene, strongly indicating that the Airn RNA rather than its transcription is responsible for silencing distant imprinted genes. Furthermore, we found that silent imprinted genes were covered with large blocks of H3K27me3 on the repressed paternal allele. Therefore we propose an alternative hypothesis whereby the chromosome interactions may initially guide the lncRNA to target imprinted promoters and recruit repressive chromatin, and that these interactions are lost once silencing is established.

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.

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.

Entry of human rhinovirus 89 via ICAM-1 into HeLa epithelial cells is inhibited by actin skeleton disruption and by bafilomycin.

HRV89, a major-group rhinovirus, uses intercellular adhesion molecule 1 (ICAM-1) for cell entry, while minor-group HRV2 uses the LDL receptor for clathrin-mediated endocytosis. Entry of HRV89 into HeLa epithelial cells was found to be inefficient, and infectious virus was still detected on the plasma membrane after 3 h of incubation with the cells. Endocytosis, and consequently infection, of HRV89 but not of HRV2, was almost completely blocked by the actin-polymerization inhibitor cytochalasin D, while the phosphatidylinositol 3-kinase inhibitor LY294002 had no effect on infection with either virus. Cytochalasin D also inhibited major-group HRV infection of rhabdomyosarcoma cells expressing ICAM-1 when the time available for uncoating was limited to 30 min. Although cholesterol depletion strongly inhibited HRV89 infection of HeLa cells, it only slightly affected HRV89 endocytosis, indicating that a lipid raft environment was not essential for virus uptake. The sodium-proton exchange inhibitor 5-(N-ethyl-N-isopropyl) amiloride (EIPA) significantly reduced cell entry and infection by HRV89 only at a concentration that also inhibited HRV2 infection and Alexa 488-transferrin entry. These data rule out classical macropinocytosis as an infectious entry pathway of HRV89 in HeLa cells. Notably, the proton ATPase inhibitor bafilomycin strongly affected cell entry of both viruses, suggesting a role for submembraneous pH in rhinovirus endocytosis.