Divergence in cis-regulatory networks: taking the 'species' out of cross-species analysis
© BioMed Central Ltd 2008
Published: 04 November 2008
Skip to main content
© BioMed Central Ltd 2008
Published: 04 November 2008
Many essential transcription factors have conserved roles in regulating biological programs, yet their genomic occupancy can diverge significantly. A new study demonstrates that such variations are primarily due to cis-regulatory sequences, rather than differences between the regulators or nuclear environments.
Genetic studies in a range of organisms reveal that essential transcription factors (TFs) tend not only to be conserved in sequence but also in function. For example, the NKx2.5 TFs are essential for heart development in species as diverse as mice , zebrafish , Xenopus , humans  and Drosophila . At a structural level, the DNA-binding domains of many orthologous TFs are highly similar over large evolutionary distances, allowing them to bind to identical DNA motifs. In fact, cross-species experiments demonstrate that orthologous TFs can regulate the same target genes and even rescue some mutant phenotypes [6, 7]. It is thus reasonable to assume that conserved TFs, which lead to the development and maintenance of orthologous tissues , regulate conserved sets of downstream target genes as part of conserved gene-regulatory networks.
It therefore came as a surprise when recent studies on DNA binding of the TFs Zeste among Drosophila species  and Ste12 and Tec1 across yeast species  indicated that individual binding events turn over rapidly during evolution. A similar discovery has been made for liver-specific TFs among vertebrates . Mouse and human hepatocytes have a similar complement of gene expression  and are defined by a set of highly conserved TFs , yet the underlying cis-regulatory network appears to have diverged extensively. Odom et al.  showed that relatively few TF-binding events - perhaps even a small minority in some cases - are conserved between the two species. Their results indicate that the target genes of hepatocyte TFs differ significantly from mouse to human, and even when orthologous genes are targeted by the same TF, the exact pattern of binding events at the conserved DNA motifs is different. These results, together with those from Drosophila and yeast, argue that binding events are subject to less selective pressure than previously anticipated, which has important implications for the degree of divergence in cis-regulatory networks.
Despite the high conservation of the TFs assayed in the studies mentioned above, it is conceivable that the differences in binding signatures between species were due to differential interaction with cofactors (owing to differences in protein-protein interactions or cofactor availability), other species-specific nuclear conditions, or simply because of experimental variables. Alternatively, the genomic sequences themselves might be different enough to trigger species-specific TF-binding signatures. A new study by Wilson et al.  addresses precisely this question by using a mouse model for human trisomy 21. This partially mosaic 'Tc1' mouse line carries most of human chromosome 21 in addition to the entire murine chromosome complement . Assaying TF binding to both the mouse and human chromosomes in the same cells eliminates many technical variables, as well as variables pertaining to interspecies differences in nuclear environment. Importantly, all assayed TFs are derived from the mouse genome, as none of them, nor any known cofactors or other hepatocyte-specific factors, are encoded on human chromosome 21 . The authors were therefore able to ask: 'Does a human chromosome in the murine nuclear context exhibit human-like, mouse-like, or a mixture of TF binding signatures?' In other words, does the human genetic material direct where TFs bind, or do mouse TFs bind elsewhere - maybe even to sites orthologous to the cognate mouse chromosome sites?
The authors focus on the binding events exhibited by three hepatocyte-specific TFs (HNF1a, HNF4a, and HNF6) across the orthologous regions of human chromosome 21 (WT-HsChr21) in human liver tissue, human chromosome 21 in mice (Tc1-HsChr21) and mouse chromosome 16 (Tc1-MmChr16) . Only about a third to a half of identified bound regions are shared among all three chromosomes, confirming the stark differences in TF-binding events between mouse and human observed previously . Importantly, the vast majority of the remaining peaks on human chromosome 21 are not found on the mouse chromosome, but rather recapitulate peaks found on chromatin isolated from human liver tissue . The fact that mouse TFs, in the mouse nuclear environment, still recapitulate human-like binding signatures on a human-derived chromosome strongly indicate that it is the human chromosomal sequence that is primarily responsible for the placement of transcription factors (cis-directed), rather than changes in the regulators or the regulative environment (trans-directed). It is interesting to note that a small number of peaks (5 out of 173 non-shared peaks) appear to be trans-directed (Tc1-HsChr21 peaks align with Tc1-MmChr16 peaks), and may warrant further investigation in their own right.
Having established that the TFs are placed on the DNA in a species-specific sequence-dependent manner, the authors examined an event downstream of TF recruitment - the placement of the basal transcriptional machinery. They did this by chromatin immunoprecipitation followed by microarray analysis (ChIP-chip) against the trimethylated state of lysine 4 on histone H3 (H3K4me3) . Whereas the majority of the H3K4me3 peaks detected can be identified in equivalent positions on human chromosome 21 and the corresponding mouse regions, some of these methylation marks appear species-specific, as indicated previously .
In Tc1 mice, the authors report 78 alignable H3K4me3 marks, of which about two-thirds (53) are shared between mouse and human. Of the remaining 25 peaks, 18 Tc1-HsChr21 peaks were also found on the WT-HsChr21 (cis-directed, mostly not at transcriptional start sites (TSSs)), indicating that the human chromosomal sequence plays a significant (albeit not necessarily direct) role in the placement of at least some epigenetic marks . Curiously, the remaining seven H3K4me3 marks appear trans-directed (also found on Tc1-MmChr16, mostly at TSSs) and may represent cases where human chromosomal regions are recognized and treated by the mouse nuclear environment in a mouse-specific manner. Finally, the authors find that the transcriptional profile of human chromosome 21 genes in Tc1 mice resembles their transcription in the native human environment, rather than the transcriptional profile of their murine orthologs .
Studies of cis-evolution have largely focused on individual enhancers or cis-regulatory modules (CRMs) [16–19]; however, more recent studies venture to identify cis-regulatory differences on a global scale [10, 11, 20]. The use of the trans-chromosomic Tc1 mice  to address species-specific differences in transcriptional regulation is certainly elegant, and one wonders if, in principle, a similar system might be extendable to other chromosomes, transcription factors, tissues, developmental contexts and species.
The study by Wilson et al.  provides strong evidence that it is the genomic sequence, rather than differences in nuclear environment, which is primarily responsible for the differences in mouse versus human TF occupancy. This underlines the importance of measuring TF binding directly rather than inferring occupancy through sequence and phylogenomic analysis. The ability of murine hepatocyte TFs to 'read' the transcriptional program of a human chromosome, even when placed in the nuclear environment of the mouse, a species separated from humans by approximately 75-100 million years, adds to the growing evidence that cis-regulatory changes are a major (if not the) driving force of evolutionary change .
As with all interspecies comparisons, the conclusions that can be drawn from these studies are largely dependent on reliable alignment of the genomes and the faithful mapping of orthologous regions . For example, misalignment of ChIP peaks will skew data, as orthologous peaks could easily be misannotated as trans-, rather than cis-directed. The task of sequence alignment is relatively tractable when performing interspecies comparisons of coding regions, but the challenge is exponentially more difficult when comparing noncoding regions. Even with largely syntenic chromosomes (such as mouse chromosome 16 versus human chromosome 21), defining orthologous peaks is very difficult. Choosing the proper species for cross-species analyses is extremely important and depends on the precise question being asked (for example, ): whereas comparisons over large evolutionary distances might yield insights into gross changes in gene regulatory networks [10, 12], comparisons over smaller distances might be more fruitful when dissecting differences in the underlying cis-regulatory networks [9, 16].
One important remaining question from the hepatocyte studies [11, 12] concerns the functional activity of species-specific TF binding. Although the authors show by Solexa sequencing that most of the species-unique H3K4me3 marks are associated with transcription, a precise analysis of the overlap of TF-bound regions with regions of active transcription (deduced from either H3K4me3 marks or expression profiling) was not presented. Do the genomic regions bound in both human and mouse correspond to regulatory regions in the vicinity of active transcription (that is, in close proximity to shared H3K4me3 peaks), whereas uniquely bound regions do not? In other words, do conserved binding events represent the functional sites? If this is the case, it suggests that once 'functional' cis-binding events are distilled from non-functional ones, there may be significant conservation in cis-regulatory networks. Alternatively, although the general properties of gene regulatory networks are conserved, the underlying cis-regulatory networks may have undergone significant divergence. No doubt future cis-evolutionary studies, both at individual loci and genome-wide, will begin to unravel this question and provide exciting insights into the general principles underlying the changes in cis-regulatory networks during speciation.