Genomic identification of regulatory elements by evolutionary sequence comparison and functional analysis

SArKS: de novo discovery of gene expression regulatory motif sites and domains by suffix array kernel smoothing

MOTIVATION

We set out to develop an algorithm that can mine differential gene expression data to identify candidate cell type-specific DNA regulatory sequences. Differential expression is usually quantified as a continuous score-fold-change, test-statistic, P-value-comparing biological classes. Unlike existing approaches, our de novo strategy, termed SArKS, applies non-parametric kernel smoothing to uncover promoter motif sites that correlate with elevated differential expression scores. SArKS detects motif k-mers by smoothing sequence scores over sequence similarity. A second round of smoothing over spatial proximity reveals multi-motif domains (MMDs). Discovered motif sites can then be merged or extended based on adjacency within MMDs. False positive rates are estimated and controlled by permutation testing.

RESULTS

We applied SArKS to published gene expression data representing distinct neocortical neuron classes in Mus musculus and interneuron developmental states in Homo sapiens. When benchmarked against several existing algorithms using a cross-validation procedure, SArKS identified larger motif sets that formed the basis for regression models with higher correlative power.

AVAILABILITY AND IMPLEMENTATION

https://github.com/denniscwylie/sarks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Zebrafish: A Powerful Model for Understanding the Functional Relevance of Noncoding Region Mutations in Human Genetic Diseases

Determining aetiology of genetic disorders caused by damaging mutations in protein-coding genes is well established. However, understanding how mutations in the vast stretches of the noncoding genome contribute to genetic abnormalities remains a huge challenge. Cis-regulatory elements (CREs) or enhancers are an important class of noncoding elements. CREs function as the primary determinants of precise spatial and temporal regulation of their target genes during development by serving as docking sites for tissue-specific transcription factors. Although a large number of potential disease-associated CRE mutations are being identified in patients, lack of robust methods for mechanistically linking these mutations to disease phenotype is currently hampering the understanding of their roles in disease aetiology. Here, we have described the various systems available for testing the CRE potential of stretches of noncoding regions harbouring mutations implicated in human disease. We highlight advances in the field leading to the establishment of zebrafish as a powerful system for robust and cost-effective functional assays of CRE activity, enabling rapid identification of causal variants in regulatory regions and the validation of their role in disruption of appropriate gene expression.

Distribution of Initiation Times Reveals Mechanisms of Transcriptional Regulation in Single Cells

Transcription is the dominant point of control of gene expression. Biochemical studies have revealed key molecular components of transcription and their interactions, but the dynamics of transcription initiation in cells is still poorly understood. This state of affairs is being remedied with experiments that observe transcriptional dynamics in single cells using fluorescent reporters. Quantitative information about transcription initiation dynamics can also be extracted from experiments that use electron micrographs of RNA polymerases caught in the act of transcribing a gene (Miller spreads). Inspired by these data, we analyze a general stochastic model of transcription initiation and elongation and compute the distribution of transcription initiation times. We show that different mechanisms of initiation leave distinct signatures in the distribution of initiation times that can be compared to experiments. We analyze published data from micrographs of RNA polymerases transcribing ribosomal RNA genes in Escherichia coli and compare the observed distributions of interpolymerase distances with the predictions from previously hypothesized mechanisms for the regulation of these genes. Our analysis demonstrates the potential of measuring the distribution of time intervals between initiation events as a probe for dissecting mechanisms of transcription initiation in live cells.

Sequence based prediction of enhancer regions from DNA random walk

Regulatory elements play a critical role in development process of eukaryotic organisms by controlling the spatio-temporal pattern of gene expression. Enhancer is one of these elements which contributes to the regulation of gene expression through chromatin loop or eRNA expression. Experimental identification of a novel enhancer is a costly exercise, due to which there is an interest in computational approaches to predict enhancer regions in a genome. Existing computational approaches to achieve this goal have primarily been based on training of high-throughput data such as transcription factor binding sites (TFBS), DNA methylation, and histone modification marks etc. On the other hand, purely sequence based approaches to predict enhancer regions are promising as they are not biased by the complexity or context specificity of such datasets. In sequence based approaches, machine learning models are either directly trained on sequences or sequence features, to classify sequences as enhancers or non-enhancers. In this paper, we derived statistical and nonlinear dynamic features along with k-mer features from experimentally validated sequences taken from Vista Enhancer Browser through random walk model and applied different machine learning based methods to predict whether an input test sequence is enhancer or not. Experimental results demonstrate the success of proposed model based on Ensemble method with area under curve (AUC) 0.86, 0.89, and 0.87 in B cells, T cells, and Natural killer cells for histone marks dataset.

Sex chromosomes drive gene expression and regulatory dimorphisms in mouse embryonic stem cells

BACKGROUND

Pre-implantation embryos exhibit sexual dimorphisms in both primates and rodents. To determine whether these differences reflected sex-biased expression patterns, we generated transcriptome profiles for six 40,XX, six 40,XY, and two 39,X mouse embryonic stem (ES) cells by RNA sequencing.

RESULTS

We found hundreds of coding and non-coding RNAs that were differentially expressed between male and female cells. Surprisingly, the majority of these were autosomal and included RNA encoding transcription and epigenetic and chromatin remodeling factors. We showed differential Prdm14-responsive enhancer activity in male and female cells, correlating with the sex-specific levels of Prdm14 expression. This is the first time sex-specific enhancer activity in ES cells has been reported. Evaluation of X-linked gene expression patterns between our XX and XY lines revealed four distinct categories: (1) genes showing 2-fold greater expression in the female cells; (2) a set of genes with expression levels well above 2-fold in female cells; (3) genes with equivalent RNA levels in male and female cells; and strikingly, (4) a small number of genes with higher expression in the XY lines. Further evaluation of autosomal gene expression revealed differential expression of imprinted loci, despite appropriate parent-of-origin patterns. The 39,X lines aligned closely with the XY cells and provided insights into potential regulation of genes associated with Turner syndrome in humans. Moreover, inclusion of the 39,X lines permitted three-way comparisons, delineating X and Y chromosome-dependent patterns.

CONCLUSIONS

Overall, our results support the role of the sex chromosomes in establishing sex-specific networks early in embryonic development and provide insights into effects of sex chromosome aneuploidies originating at those stages.

The most common technologies and tools for functional genome analysis

Since the sequence of the human genome is complete, the main issue is how to understand the information written in the DNA sequence. Despite numerous genome-wide studies that have already been performed, the challenge to determine the function of genes, gene products, and also their interaction is still open. As changes in the human genome are highly likely to cause pathological conditions, functional analysis is vitally important for human health.

For many years there have been a variety of technologies and tools used in functional genome analysis. However, only in the past decade there has been rapid revolutionizing progress and improvement in high-throughput methods, which are ranging from traditional real-time polymerase chain reaction to more complex systems, such as next-generation sequencing or mass spectrometry. Furthermore, not only laboratory investigation, but also accurate bioinformatic analysis is required for reliable scientific results. These methods give an opportunity for accurate and comprehensive functional analysis that involves various fields of studies: genomics, epigenomics, proteomics, and interactomics. This is essential for filling the gaps in the knowledge about dynamic biological processes at both cellular and organismal level. However, each method has both advantages and limitations that should be taken into account before choosing the right method for particular research in order to ensure successful study. For this reason, the present review paper aims to describe the most frequent and widely-used methods for the comprehensive functional analysis.

In Vivo Hepatic Enhancer Elements in the Human ABCG2 Locus

ABCG2 encodes the mitoxantrone resistance protein (MXR; breast cancer resistance protein), an ATP-binding cassette (ABC) efflux membrane transporter. Computational analysis of the ∼300 kb region of DNA surrounding ABCG2 (chr4:88911376-89220011, hg19) identified 30 regions with potential cis-regulatory capabilities. These putative regulatory regions were tested for their enhancer and suppressor activity in a human liver cell line using luciferase reporter assays. The in vitro enhancer and suppressor assays identified four regions that decreased gene expression and five regions that increased expression >1.6-fold. Four of five human hepatic in vitro enhancers were confirmed as in vivo liver enhancers using the mouse hydrodynamic tail vein injection assay. Two of the in vivo liver enhancers (ABCG2RE1 and ABCG2RE9) responded to 17β-estradiol or rifampin in human cell lines, and ABCG2RE9 had ChIP-seq evidence to support the binding of several transcription factors and the transcriptional coactivator p300 in human hepatocytes. This study identified genomic regions surrounding human ABCG2 that can function as regulatory elements, some with the capacity to alter gene expression upon environmental stimulus. The results from this research will drive future investigations of interindividual variation in ABCG2 expression and function that contribute to differences in drug response.

Transcriptional Interference Promotes Rapid Expression Divergence of Drosophila Nested Genes

Nested genes are the most common form of protein-coding overlap in eukaryotic genomes. Previous studies have shown that nested genes accumulate rapidly over evolutionary time, typically via the insertion of short young duplicate genes into long introns. However, the evolutionary relationship between nested genes remains unclear. Here, I compare RNA-seq expression profiles of nested, proximal intra-chromosomal, intermediate intra-chromosomal, distant intra-chromosomal, and inter-chromosomal gene pairs in two Drosophila species. I find that expression profiles of nested genes are more divergent than those of any other class of genes, supporting the hypothesis that concurrent expression of nested genes is deleterious due to transcriptional interference. Further analysis reveals that expression profiles of derived nested genes are more divergent than those of their ancestral un-nested orthologs, which are more divergent than those of un-nested genes with similar genomic features. Thus, gene expression divergence between nested genes is likely caused by selection against nesting of genes with insufficiently divergent expression profiles, as well as by continued expression divergence after nesting. Moreover, expression divergence and sequence evolutionary rates are elevated in young nested genes and reduced in old nested genes, indicating that a burst of rapid evolution occurs after nesting. Together, these findings suggest that similarity between expression profiles of nested genes is deleterious due to transcriptional interference, and that natural selection addresses this problem both by eradicating highly deleterious nestings and by enabling rapid expression divergence of surviving nested genes, thereby quickly limiting or abolishing transcriptional interference.

Otx2 ChIP-seq reveals unique and redundant functions in the mature mouse retina

During mouse retinal development and into adulthood, the transcription factor Otx2 is expressed in pigment epithelium, photoreceptors and bipolar cells. In the mature retina, Otx2 ablation causes photoreceptor degeneration through a non-cell-autonomous mechanism involving Otx2 function in the supporting RPE. Surprisingly, photoreceptor survival does not require Otx2 expression in the neural retina, where the related Crx homeobox gene, a major regulator of photoreceptor development, is also expressed. To get a deeper view of mouse Otx2 activities in the neural retina, we performed chromatin-immunoprecipitation followed by massively parallel sequencing (ChIP-seq) on Otx2. Using two independent ChIP-seq assays, we identified consistent sets of Otx2-bound cis-regulatory elements. Comparison with our previous RPE-specific Otx2 ChIP-seq data shows that Otx2 occupies different functional domains of the genome in RPE cells and in neural retina cells and regulates mostly different sets of genes. To assess the potential redundancy of Otx2 and Crx, we compared our data with Crx ChIP-seq data. While Crx genome occupancy markedly differs from Otx2 genome occupancy in the RPE, it largely overlaps that of Otx2 in the neural retina. Thus, in accordance with its essential role in the RPE and its non-essential role in the neural retina, Otx2 regulates different gene sets in the RPE and the neural retina, and shares an important part of its repertoire with Crx in the neural retina. Overall, this study provides a better understanding of gene-regulatory networks controlling photoreceptor homeostasis and disease.

Generation and application of signaling pathway reporter lines in zebrafish

In the last years, we have seen the emergence of different tools that have changed the face of biology from a simple modeling level to a more systematic science. The transparent zebrafish embryo is one of the living models in which, after germline transformation with reporter protein-coding genes, specific fluorescent cell populations can be followed at single-cell resolution. The genetically modified embryos, larvae and adults, resulting from the transformation, are individuals in which time lapse analysis, digital imaging quantification, FACS sorting and next-generation sequencing can be performed in specific times and tissues. These multifaceted genetic and cellular approaches have permitted to dissect molecular interactions at the subcellular, intercellular, tissue and whole-animal level, thus allowing integration of cellular and developmental genetics with molecular imaging in the resulting frame of modern biology. In this review, we describe a new step in the zebrafish road to system biology, based on the use of transgenic biosensor animals expressing fluorescent proteins under the control of signaling pathway-responsive cis-elements. In particular, we provide here the rationale and details of this powerful tool, trying to focus on its huge potentialities in basic and applied research, while also discussing limits and potential technological evolutions of this approach.

Islet1 is a direct transcriptional target of the homeodomain transcription factor Shox2 and rescues the Shox2-mediated bradycardia

The heart's rhythm is initiated and regulated by a group of specialized cells in the sinoatrial node (SAN), the primary pacemaker of the heart. Abnormalities in the development of the SAN can result in irregular heart rates (arrhythmias). Although several of the critical genes important for SAN formation have been identified, our understanding of the transcriptional network controlling SAN development remains at a relatively early stage. The homeodomain transcription factor Shox2 is involved in the specification and patterning of the SAN. While the Shox2 knockout in mice results in embryonic lethality due to severe cardiac defects including improper SAN development, Shox2 knockdown in zebrafish causes a reduced heart rate (bradycardia). In order to gain deeper insight into molecular pathways involving Shox2, we compared gene expression levels in right atria of wildtype and Shox2 (-/-) hearts using microarray experiments and identified the LIM homeodomain transcription factor Islet1 (Isl1) as one of its putative target genes. The downregulation of Isl1 expression in Shox2 (-/-) hearts was confirmed and the affected region narrowed down to the SAN by whole-mount in situ hybridization. Using luciferase reporter assays and EMSA studies, we identified two specific SHOX2 binding sites within intron 2 of the ISL1 locus. We also provide functional evidence for Isl1 as a transcriptional target of Shox2 by rescuing the Shox2-mediated bradycardia phenotype with Isl1 using zebrafish as a model system. Our findings demonstrate a novel epistatic relationship between Shox2 and Isl1 in the heart with important developmental consequences for SAN formation and heart beat.

New advances in the genetic basis of atrial fibrillation

Over the past decade, compelling evidence has emerged from population-based studies to suggest that AF is a heritable disease. More recently, we have begun to elucidate the genetic substrate underlying AF. Genome-wide association studies (GWAS) have led to the identification of multiple risk loci that confer increased susceptibility to the arrhythmia. These loci harbor intriguing candidate genes including those encoding ion channels, transcription factors, and signaling molecules. Current efforts are ongoing to functionally validate the role of these genes in disease pathogenesis. In the future, novel genotyping technologies such as exome sequencing and whole-genome sequencing promise to uncover a greater proportion of the heritability underlying AF. In this article we review recent advances in AF genetics research and discuss future developments in the field.

Transcriptional enhancers in protein-coding exons of vertebrate developmental genes

Many conserved noncoding sequences function as transcriptional enhancers that regulate gene expression. Here, we report that protein-coding DNA also frequently contains enhancers functioning at the transcriptional level. We tested the enhancer activity of 31 protein-coding exons, which we chose based on strong sequence conservation between zebrafish and human, and occurrence in developmental genes, using a Tol2 transposable GFP reporter assay in zebrafish. For each exon we measured GFP expression in hundreds of embryos in 10 anatomies via a novel system that implements the voice-recognition capabilities of a cellular phone. We find that 24/31 (77%) exons drive GFP expression compared to a minimal promoter control, and 14/24 are anatomy-specific (expression in four anatomies or less). GFP expression driven by these coding enhancers frequently overlaps the anatomies where the host gene is expressed (60%), suggesting self-regulation. Highly conserved coding sequences and highly conserved noncoding sequences do not significantly differ in enhancer activity (coding: 24/31 vs. noncoding: 105/147) or tissue-specificity (coding: 14/24 vs. noncoding: 50/105). Furthermore, coding and noncoding enhancers display similar levels of the enhancer-related histone modification H3K4me1 (coding: 9/24 vs noncoding: 34/81). Meanwhile, coding enhancers are over three times as likely to contain an H3K4me1 mark as other exons of the host gene. Our work suggests that developmental transcriptional enhancers do not discriminate between coding and noncoding DNA and reveals widespread dual functions in protein-coding DNA.

An embryonic stem cell-based system for rapid analysis of transcriptional enhancers

With the growing use of genome-wide screens for cis-regulatory elements, there is a pressing need for platforms that enable fast and cost-effective experimental validation of identified hits in relevant developmental and tissue contexts. Here, we describe a murine embryonic stem cell (ESC)-based system that facilitates rapid analysis of putative transcriptional enhancers. Candidate enhancers are targeted with high efficiency to a defined genomic locus via recombinase-mediated cassette exchange. Targeted ESCs are subsequently differentiated in vitro into desired cell types, where enhancer activity is monitored by reporter gene expression. As a proof of principle, we analyzed a previously characterized, Sonic hedgehog (Shh)-dependent, V3 interneuron progenitor (pV3)-specific enhancer for the Nkx2.2 gene, and observed highly specific enhancer activity. Given the broad potential of ESCs to generate a spectrum of cell types, this system can serve as an effective platform for the characterization of gene regulatory networks controlling cell fate specification and cell function.

Proliferation-associated POU4F2/Brn-3b transcription factor expression is regulated by oestrogen through ERα and growth factors via MAPK pathway

INTRODUCTION

In cancer cells, elevated transcription factor-related Brn-3a regulator isolated from brain cDNA (Brn-3b) transcription factor enhances proliferation in vitro and increases tumour growth in vivo whilst conferring drug resistance and migratory potential, whereas reducing Brn-3b slows growth both in vitro and in vivo. Brn-3b regulates distinct groups of key target genes that control cell growth and behaviour. Brn-3b is elevated in >65% of breast cancer biopsies, but mechanisms controlling its expression in these cells are not known.

METHODS

Bioinformatics analysis was used to identify the regulatory promoter region and map transcription start site as well as transcription factor binding sites. Polymerase chain reaction (PCR) cloning was used to generate promoter constructs for reporter assays. Chromatin immunoprecipitation and site-directed mutagenesis were used to confirm the transcription start site and autoregulation. MCF-7 and Cos-7 breast cancer cells were used. Cells grown in culture were transfected with Brn-3b promoter and treated with growth factors or estradiol to test for effects on promoter activity. Quantitative reverse transcriptase PCR assays and immunoblotting were used to confirm changes in gene and protein expression.

RESULTS

We cloned the Brn-3b promoter, mapped the transcription start site and showed stimulation by estradiol and growth factors, nerve growth factor and epidermal growth factor, which are implicated in breast cancer initiation and/or progression. The effects of growth factors are mediated through the mitogen-activated protein kinase pathway, whereas hormone effects act via oestrogen receptor α (ERα). Brn-3b also autoregulates its expression and cooperates with ERα to further enhance levels.

CONCLUSIONS

Key regulators of growth in cancer cells, for example, oestrogens and growth factors, can stimulate Brn-3b expression, and autoregulation also contributes to increasing Brn-3b in breast cancers. Since increasing Brn-3b profoundly enhances growth in these cells, understanding how Brn-3b is increased in breast cancers will help to identify strategies for reducing its expression and thus its effects on target genes, thereby reversing its effects in breast cancer cells.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

Efficient discovery of ASCL1 regulatory sequences through transgene pooling

Zebrafish transgenesis is a powerful and increasingly common strategy to assay vertebrate transcriptional regulatory control. Several challenges remain, however, to the broader application of this technique; they include increasing the rate with which transgenes can be analyzed and maximizing the informational value of the data generated. Presently, many rely on the injection of individual constructs and the analysis of resulting reporter expression in mosaic G0 embryos. Here, we contrast these approaches, examining whether injecting pooled transgene constructs can increase the efficiency with which regulatory sequences can be assayed, restricting analysis to the offspring of germ line transmitting transgenic zebrafish in an effort to reduce potential subjectivity. We selected a 64kb interval encompassing the human ASCL1 locus as our model interval and report the analysis of 9 highly conserved putative enhancers therein. We identified 32 transgene-positive zebrafish, transmitting one or more independent constructs displaying ASCL1-like regulatory control. Through examination of embryos harboring one or more transgenes, we demonstrate that five of the nine sequences account for the observed control and describe their likely roles in ASCL1 regulation. These data demonstrate the utility of this approach and its potential for further adaptation and higher throughput application.

Genetic modifiers of Hb E/beta0 thalassemia identified by a two-stage genome-wide association study

Background

Patients with Hb E/β⁰ thalassemia display remarkable variability in disease severity. To identify genetic modifiers influencing disease severity, we conducted a two-stage genome scan in groups of 207 mild and 305 severe unrelated patients from Thailand with Hb E/β⁰ thalassemia and normal α-globin genes.

Methods

First, we estimated and compared the allele frequencies of approximately 110,000 gene-based single nucleotide polymorphisms (SNPs) in pooled DNAs from different severity groups. The 756 SNPs that showed reproducible allelic differences at P < 0.02 by pooling were selected for individual genotyping.

Results

After adjustment for age, gender and geographic region, logistic regression models showed 50 SNPs significantly associated with disease severity (P < 0.05) after Bonferroni adjustment for multiple testing. Forty-one SNPs in a large LD block within the β-globin gene cluster had major alleles associated with severe disease. The most significant was bthal_bg200 (odds ratio (OR) = 5.56, P = 2.6 × 10^-13). Seven SNPs in two distinct LD blocks within a region centromeric to the β-globin gene cluster that contains many olfactory receptor genes were also associated with disease severity; rs3886223 had the strongest association (OR = 3.03, P = 3.7 × 10^-11). Several previously unreported SNPs were also significantly associated with disease severity.

Conclusions

These results suggest that there may be an additional regulatory region centromeric to the β-globin gene cluster that affects disease severity by modulating fetal hemoglobin expression.

Fetal hemoglobin in sickle cell anemia: genome-wide association studies suggest a regulatory region in the 5' olfactory receptor gene cluster

In a genome-wide association study of 848 blacks with sickle cell anemia, we identified single nucleotide polymorphisms (SNPs) associated with fetal hemoglobin concentration. The most significant SNPs in a discovery sample were tested in a replication set of 305 blacks with sickle cell anemia and in subjects with hemoglobin E or beta thalassemia trait from Thailand and Hong Kong. A novel region on chromosome 11 containing olfactory receptor genes OR51B5 and OR51B6 was identified by 6 SNPs (lowest P = 4.7E-08) and validated in the replication set. An additional olfactory receptor gene, OR51B2, was identified by a novel SNP set enrichment analysis. Genome-wide association studies also validated a previously identified SNP (rs766432) in BCL11A, a gene known to affect fetal hemoglobin levels (P = 2.6E-21) and in Thailand and Hong Kong subjects. Elements within the olfactory receptor gene cluster might play a regulatory role in gamma-globin gene expression.

COX-2 Is Associated with Periodontitis in Europeans

COX-2 plays an important role in periodontitis by mediating inflammatory reactions in periodontal tissues, and the COX-2 polymorphisms rs20417 and rs689466 have been reported to be associated with periodontitis in populations of Taiwanese and Chinese ethnicity. To test whether these variants were associated with periodontitis in populations of European ethnicity, we genotyped the single-nucleotide polymorphisms (SNPs) rs689466 and rs6681231, the latter a haplotype tagging SNP (htSNP) for rs20417 (r2>0.95), in our large-analysis population of individuals with aggressive (n = 532) and chronic periodontitis (n = 1052), and 2873 healthy control individuals. The rare G allele of htSNP rs6681231 was associated with aggressive periodontitis prior to and after adjustment for the covariates smoking, diabetes, and gender, with an odds ratio of 1.57 (95% confidence interval 1.18–2.08; p = 0.002). The validation of the association of rs20417 by the htSNP rs6681231 provides evidence for a general genetic risk of COX-2 variants in the pathogenesis of periodontitis.

Conserved expression without conserved regulatory sequence: the more things change, the more they stay the same

Regulatory regions with similar transcriptional output often have little overt sequence similarity, both within and between genomes. Although cis- and trans-regulatory changes can contribute to sequence divergence without dramatically altering gene expression outputs, heterologous DNA often functions similarly in organisms that share little regulatory sequence similarities (e.g. human DNA in fish), indicating that trans-regulatory mechanisms tend to diverge more slowly and can accommodate a variety of cis-regulatory configurations. This capacity to 'tinker' with regulatory DNA probably relates to the complexity, robustness and evolvability of regulatory systems, but cause-and-effect relationships among evolutionary processes and properties of regulatory systems remain a topic of debate. The challenge of understanding the concrete mechanisms underlying cis-regulatory evolution - including the conservation of function without the conservation of sequence - relates to the challenge of understanding the function of regulatory systems in general. Currently, we are largely unable to recognize functionally similar regulatory DNA.