Sequence composition similarities with the 7SL RNA are highly predictive of functional genomic features

Transposable elements derived from the 7SL RNA gene, such as Alu elements in primates, have had remarkable success in several mammalian lineages. The results presented here show a broad spectrum of functions for genomic segments that display sequence composition similarities with the 7SL RNA gene. Using thoroughly documented loci, we report that DNaseI-hypersensitive sites can be singled out in large genomic sequences by an assessment of sequence composition similarities with the 7SL RNA gene. We apply a root word frequency approach to illustrate a distinctive relationship between the sequence of the 7SL RNA gene and several classes of functional genomic features that are not presumed to be of transposable origin. Transposable elements that show noticeable similarities with the 7SL sequence include Alu sequences, as expected, but also long terminal repeats and the 5′-untranslated regions of long interspersed repetitive elements. In sequences masked for repeated elements, we find, when using the 7SL RNA gene as query sequence, distinctive similarities with promoters, exons and distal gene regulatory regions. The latter being the most notoriously difficult to detect, this approach may be useful for finding genomic segments that have regulatory functions and that may have escaped detection by existing methods.     

Potential for genomic instability associated with retrotranspositionally-incompetent L1 loci

Expression of the L1 retrotransposon can damage the genome through insertional mutagenesis and the generation of DNA double-strand breaks (DSBs). The majority of L1 loci in the human genome are 5′-truncated and therefore incapable of retrotransposition. While thousands of full-length L1 loci remain, most are retrotranspositionally-incompetent due to inactivating mutations. However, mutations leading to premature stop codons within the L1 ORF2 sequence may yield truncated proteins that retain a functional endonuclease domain. We demonstrate that some truncated ORF2 proteins cause varying levels of toxicity and DNA damage when chronically overexpressed in mammalian cells. Furthermore, transfection of some ORF2 constructs containing premature stop codons supported low levels of Alu retrotransposition, demonstrating the potential for select retrotranspositionally-incompetent L1 loci to generate genomic instability. This result suggests yet another plausible explanation for the relative success of Alu elements in populating the human genome. Our data suggest that a subset of retrotranspositionally-incompetent L1s, previously considered to be harmless to genomic integrity, may have the potential to cause chronic DNA damage by introducing DSBs and mobilizing Alu. These results imply that the number of known L1 loci in the human genome that potentially threaten its stability may not be limited to the retrotranspositionally active loci.     

Expression of human IFN-gamma genomic DNA in transgenic mice

We have introduced an 8.6-kb fragment of human genomic DNA containing the full length IFN-gamma gene into the mouse germline. The transgenic animals had no biologic or developmental defects as human IFN-gamma does not bind to the mouse IFN R. Regulation of the transgene paralleled that of the endogenous murine IFN-gamma gene in that: 1) it is not expressed constitutively in any tissue examined thus far, 2) it can be induced in thymus and spleen cells by T cell mitogens, 3) it is not expressed in B cells stimulated by LPS, and 4) it produces normal mRNA and biologically active IFN protein. Whereas expression of the transgene is likely restricted to T cells, we had observed that both fibroblasts and B cell lines could express the same DNA when transfected in vitro; this indicates that in vivo, developmental factors restrict expression of the IFN-gamma gene to T cells. These findings also indicate that the 8.6-kb fragment contains the regulatory elements necessary for normal tissue specific expression in vivo. Moreover, they indicate that the regulatory elements for this gene are completely preserved over the phylogenetic distance separating mouse and man, even though substantial drift has occurred in the structural gene, and probably in the IFN-gamma R as well.     

Cell-type-specific long-range looping interactions identify distant regulatory elements of the CFTR gene

Identification of regulatory elements and their target genes is complicated by the fact that regulatory elements can act over large genomic distances. Identification of long-range acting elements is particularly important in the case of disease genes as mutations in these elements can result in human disease. It is becoming increasingly clear that long-range control of gene expression is facilitated by chromatin looping interactions. These interactions can be detected by chromosome conformation capture (3C). Here, we employed 3C as a discovery tool for identification of long-range regulatory elements that control the cystic fibrosis transmembrane conductance regulator gene, CFTR. We identified four elements in a 460-kb region around the locus that loop specifically to the CFTR promoter exclusively in CFTR expressing cells. The elements are located 20 and 80 kb upstream; and 109 and 203 kb downstream of the CFTR promoter. These elements contain DNase I hypersensitive sites and histone modification patterns characteristic of enhancers. The elements also interact with each other and the latter two activate the CFTR promoter synergistically in reporter assays. Our results reveal novel long-range acting elements that control expression of CFTR and suggest that 3C-based approaches can be used for discovery of novel regulatory elements.     

Functional validation of mouse tyrosinase non-coding regulatory DNA elements by CRISPR–Cas9-mediated mutagenesis

Newly developed genome-editing tools, such as the clustered regularly interspaced short palindromic repeat (CRISPR)–Cas9 system, allow simple and rapid genetic modification in most model organisms and human cell lines. Here, we report the production and analysis of mice carrying the inactivation via deletion of a genomic insulator, a key non-coding regulatory DNA element found 5′ upstream of the mouse tyrosinase (Tyr) gene. Targeting sequences flanking this boundary in mouse fertilized eggs resulted in the efficient deletion or inversion of large intervening DNA fragments delineated by the RNA guides. The resulting genome-edited mice showed a dramatic decrease in Tyr gene expression as inferred from the evident decrease of coat pigmentation, thus supporting the functionality of this boundary sequence in vivo, at the endogenous locus. Several potential off-targets bearing sequence similarity with each of the two RNA guides used were analyzed and found to be largely intact. This study reports how non-coding DNA elements, even if located in repeat-rich genomic sequences, can be efficiently and functionally evaluated in vivo and, furthermore, it illustrates how the regulatory elements described by the ENCODE and EPIGENOME projects, in the mouse and human genomes, can be systematically validated.     

QuIN: A Web Server for Querying and Visualizing Chromatin Interaction Networks

Recent studies of the human genome have indicated that regulatory elements (e.g. promoters and enhancers) at distal genomic locations can interact with each other via chromatin folding and affect gene expression levels. Genomic technologies for mapping interactions between DNA regions, e.g., ChIA-PET and HiC, can generate genome-wide maps of interactions between regulatory elements. These interaction datasets are important resources to infer distal gene targets of non-coding regulatory elements and to facilitate prioritization of critical loci for important cellular functions. With the increasing diversity and complexity of genomic information and public ontologies, making sense of these datasets demands integrative and easy-to-use software tools. Moreover, network representation of chromatin interaction maps enables effective data visualization, integration, and mining. Currently, there is no software that can take full advantage of network theory approaches for the analysis of chromatin interaction datasets. To fill this gap, we developed a web-based application, QuIN, which enables: 1) building and visualizing chromatin interaction networks, 2) annotating networks with user-provided private and publicly available functional genomics and interaction datasets, 3) querying network components based on gene name or chromosome location, and 4) utilizing network based measures to identify and prioritize critical regulatory targets and their direct and indirect interactions. AVAILABILITY: QuIN’s web server is available at http://quin.jax.org QuIN is developed in Java and JavaScript, utilizing an Apache Tomcat web server and MySQL database and the source code is available under the GPLV3 license available on GitHub: https://github.com/UcarLab/QuIN/.

This is a PLOS Computational Biology Software paper.

     

Annotating regulatory DNA based on man-mouse genomic comparison

Non-coding DNA segments that are conserved between the human and mouse genomic sequence are good indicators of possible regulatory sequences. Here we report on a systematic approach to delineate such conserved elements from upstream regions of orthologous gene pairs from man and mouse. We focus on orthologous genes in order to maximize our chances to find functionally similar regulatory elements. The identification of conserved elements is effected using the Waterman-Eggert local suboptimal alignment algorithm. We have modified an implementation of this algorithm such that it integrates the determination of statistical significance for the local suboptimal alignments. This has the effect of outputting a dynamically determined number of suboptimal alignments that are deemed statistically significant. Comparison with experimentally determined annotation shows a striking enrichement of regulatory sites among the conserved regions. Furthermore, the conserved regions tend to cover the promotor region described in the EPD database.     

Rapid genomic DNA changes in allotetraploid fish hybrids

Rapid genomic change has been demonstrated in several allopolyploid plant systems; however, few studies focused on animals. We addressed this issue using an allotetraploid lineage (4nAT) of freshwater fish originally derived from the interspecific hybridization of red crucian carp (Carassius auratus red var., ♀, 2n=100) × common carp (Cyprinus carpio L., ♂, 2n=100). We constructed a bacterial artificial chromosome (BAC) library from allotetraploid hybrids in the 20th generation (F₂₀) and sequenced 14 BAC clones representing a total of 592.126 kb, identified 11 functional genes and estimated the guanine–cytosine content (37.10%) and the proportion of repetitive elements (17.46%). The analysis of intron evolution using nine orthologous genes across a number of selected fish species detected a gain of 39 introns and a loss of 30 introns in the 4nAT lineage. A comparative study based on seven functional genes among 4nAT, diploid F₁ hybrids (2nF₁) (first generation of hybrids) and their original parents revealed that both hybrid types (2nF₁ and 4nAT) not only inherited genomic DNA from their parents, but also demonstrated rapid genomic DNA changes (homoeologous recombination, parental DNA fragments loss and formation of novel genes). However, 4nAT presented more genomic variations compared with their parents than 2nF₁. Interestingly, novel gene fragments were found for the iqca1 gene in both hybrid types. This study provided a preliminary genomic characterization of allotetraploid F₂₀ hybrids and revealed evolutionary and functional genomic significance of allopolyploid animals.     

Identifying co‐opted transposable elements using comparative epigenomics

The human genome gives rise to different epigenomic landscapes that define each cell type and can be deregulated in disease. Recent efforts by ENCODE, the NIH Roadmap and the International Human Epigenome Consortium (IHEC) have made significant advances towards assembling reference epigenomic maps of various tissues. Notably, these projects have found that approximately 80% of human DNA was biochemically active in at least one epigenomic assay while only approximately 10% of the sequence displayed signs of purifying selection. Given that transposable elements (TEs) make up at least 50% of the human genome and can be actively transcribed or act as regulatory elements either for their own purposes or be co‐opted for the benefit of their host; we are interested in exploring their overall contribution to the “functional” genome. Traditional methods used to identify functional DNA have relied on comparative genomics, conservation analysis and low throughput validation assays. To discover co‐opted TEs, and distinguish them from noisy genomic elements, we argue that comparative epigenomic methods will also be important.