Biology:STARR-seq

STARR-seq Methodology

STARR-seq (short for self-transcribing active regulatory region sequencing) is a method to assay enhancer activity for millions of candidates from arbitrary sources of DNA. It is used to identify the sequences that act as transcriptional enhancers in a direct, quantitative, and genome-wide manner.

In eukaryotes, transcription is regulated by sequence-specific DNA-binding proteins (transcription factors) associated with a gene’s promoter and also by distant control sequences including enhancers. Enhancers are non-coding DNA sequences, containing several binding sites for a variety of transcription factors. They typically recruit transcriptional factors that modulate chromatin structure and directly interact with the transcription machinery placed at the promoter of a gene. Enhancers are able to regulate transcription of target genes in a cell type-specific manner, independent of their location or distance from the promoter of genes. In certain contexts (see Transvection (genetics)), they can even regulate transcription of genes located in a different chromosome. However, the knowledge about enhancers so far has been limited to studies of a small number of enhancers, as they have been difficult to identify accurately at a genome-wide scale. Moreover, many regulatory elements function only in certain cell types and specific conditions.

Enhancer detection

Enhancer detection in Drosophila is an original methodology using random insertion of transposon-derived vector that encodes a reporter protein downstream of a minimal promoter. This approach allows to observe the expression of reporter in transgenic animals and provides information about nearby genes that are regulated by these sequences. The discovery and characterization of cell types along with genes involved in their determination have been significantly improved by the discovery of this technique.

During the past few years, post-genomic technologies, have displayed specific features of poised and active enhancers that have improved enhancer discovery. Development of new methods such as deep sequencing of DNase I hypersensitive sites (DNase-Seq), formaldehyde-assisted isolation of regulatory elements sequencing (FAIRE-Seq), chromatin immunoprecipitation followed by deep sequencing (ChIP-sequencing), and MNase-defined cistrome-Occupancy Analysis (MOA-seq), provide genome-wide enhancer predictions by enhancer-associated chromatin features.

Application

However, DNase-seq and FAIRE-seq alone fail to provide a direct functional or quantitative readout of enhancer activity, so reporter assays that can deduce enhancer strength from the quantitative enrichment of reporter transcripts are needed to assess enhancer activity quantitatively. Yet, these assays are not high-throughput (High throughput biology), as it is impossible to conduct millions of tests required for identification of enhancers in a genome-wide manner. The development of STARR-seq attempts to circumvent this analytical barrier. Taking advantage of the knowledge that enhancers can work independently of their relative locations, candidate sequences are placed downstream of a minimal promoter, allowing the active enhancers to transcribe themselves. The strength of each enhancer is then reflected by its relative enrichment among cellular RNAs. Such a direct coupling of candidate sequences to enhancer activity enables the parallel evaluation of millions of DNA fragments from arbitrary sources.

Methodology

Genomic DNA is randomly sheared and broken down to small fragments. Adaptors are ligated to size-selected DNA fragments. Next, adaptor-linked fragments are amplified and the PCR products are purified followed by placing candidate sequences downstream of a minimal promoter of screening vectors, giving them an opportunity to transcribe themselves. Candidate cells are then transfected with reporter library and cultured. Thereafter, total RNAs are extracted and poly-A RNAs isolated. Using reverse transcription method, cDNAs are produced, amplified and then candidate fragments are used for high-throughput paired end sequencing. Sequence reads are mapped to the reference genome and computational processing of data is carried out.

Identification of enhancers

Applying this technology to Drosophila genome, Arnold et al. found 96% of the non-repetitive genome with at least 10-fold coverage. Authors discovered that most identified enhancers (55.6%) were placed within introns, particularly in the first intron and intergenic regions. 4.5% of enhancers were located at transcription start sites (TSS), suggesting that these enhancers can start transcription and also improve transcription from a distant TSS. The strongest enhancers were near housekeeping genes such as enzymes or component of the cytoskeleton and developmental regulators such as the transcription factors. The strongest enhancer was located within the intron of the transcription factor zfh1. This transcription factor regulates neuropeptide expression and growth of larval neuromuscular junctions in Drosophila. The ribosomal protein genes were the only class of genes with poor enhancers ranking. Moreover, authors demonstrated that many genes are regulated by several independent active enhancers even in a single cell type. Furthermore, gene expression levels on average were correlated with the sum of the enhancer strengths per gene, supporting direct link between gene expression and enhancer activity.

Characterization of variant alleles

Applying this technology to the characterization and discovery of regulatory variant alleles, Vockley et al. characterized the effects of human genetic variation on non-coding regulatory element function, measuring the activity of 100 putative enhancers captured directly from the genomes of 95 members of a study cohort. This approach enables the functional fine-mapping of causal regulatory variants in regions of high linkage disequilibrium identified by eQTL analyses. This approach provides a general path forward to identify perturbations in gene regulatory elements that contribute to complex phenotypes.

Quantifying enhancer activity

STARR-seq has been used to measure the regulatory activity of DNA fragments that have been enriched for sites occupied by specific transcription factors. Cloning ChIP DNA libraries generated from chromatin immunoprecipitation of the glucocorticoid receptor into STARR-seq enabled genome-scale quantification of glucocorticoid-induced enhancer activity. This approach is useful for measuring the differences in enhancer activity between sites that are bound by the same transcription factor.

Advantages

• A quantitative genome-wide assay for enhancer detection.
• Applicable technique for screening arbitrary sources of DNA in any cell type or tissue that allows adequate introduction of reporter constructs.
• A method with high detection rate (>99%) by employing pair-end sequencing, even for sequences that contain transcript-destabilizing elements.
• Technique to evaluate the strength of enhancers quantitatively, and identify endogenously silenced enhancers by integrating them into a chromosomal context.

Future directions

By combining traditional approach with high-throughput sequencing technology and highly specialized bio-computing methods, STARR-seq is able to detect enhancers in a quantitative and genome-wide manner. The study of gene regulation and their responsible pathways in the genome during normal development and also in disease can be very demanding. Therefore, applying STARR-seq to many cell types across organisms supports identifying cell type-specific gene regulatory elements and practically assesses non-coding mutations causing disease. Recently, a related approach coupling capture of regions of interest to STARR-seq technique have been developed and extensively validated in mammalian cell lines.