Transcription factor interactions explain the context-dependent activity of CRX binding sites.

The effects of transcription factor binding sites (TFBSs) on the activity of a cis-regulatory element (CRE) depend on the local sequence context. In rod photoreceptors, binding sites for the transcription factor (TF) Cone-rod homeobox (CRX) occur in both enhancers and silencers, but the sequence context that determines whether CRX binding sites contribute to activation or repression of transcription is not understood. To investigate the context-dependent activity of CRX sites, we fit neural network-based models to the activities of synthetic CREs composed of photoreceptor TFBSs. The models revealed that CRX binding sites consistently make positive, independent contributions to CRE activity, while negative homotypic interactions between sites cause CREs composed of multiple CRX sites to function as silencers. The effects of negative homotypic interactions can be overcome by the presence of other TFBSs that either interact cooperatively with CRX sites or make independent positive contributions to activity. The context-dependent activity of CRX sites is thus determined by the balance between positive heterotypic interactions, independent contributions of TFBSs, and negative homotypic interactions. Our findings explain observed patterns of activity among genomic CRX-bound enhancers and silencers, and suggest that enhancers may require diverse TFBSs to overcome negative homotypic interactions between TFBSs.

Active learning of enhancer and silencer regulatory grammar in photoreceptors.

Cis-regulatory elements (CREs) direct gene expression in health and disease, and models that can accurately predict their activities from DNA sequences are crucial for biomedicine. Deep learning represents one emerging strategy to model the regulatory grammar that relates CRE sequence to function. However, these models require training data on a scale that exceeds the number of CREs in the genome. We address this problem using active machine learning to iteratively train models on multiple rounds of synthetic DNA sequences assayed in live mammalian retinas. During each round of training the model actively selects sequence perturbations to assay, thereby efficiently generating informative training data. We iteratively trained a model that predicts the activities of sequences containing binding motifs for the photoreceptor transcription factor Cone-rod homeobox (CRX) using an order of magnitude less training data than current approaches. The model's internal confidence estimates of its predictions are reliable guides for designing sequences with high activity. The model correctly identified critical sequence differences between active and inactive sequences with nearly identical transcription factor binding sites, and revealed order and spacing preferences for combinations of motifs. Our results establish active learning as an effective method to train accurate deep learning models of cis-regulatory function after exhausting naturally occurring training examples in the genome.

A single-cell massively parallel reporter assay detects cell-type-specific gene regulation.

Massively parallel reporter gene assays are key tools in regulatory genomics but cannot be used to identify cell-type-specific regulatory elements without performing assays serially across different cell types. To address this problem, we developed a single-cell massively parallel reporter assay (scMPRA) to measure the activity of libraries of cis-regulatory sequences (CRSs) across multiple cell types simultaneously. We assayed a library of core promoters in a mixture of HEK293 and K562 cells and showed that scMPRA is a reproducible, highly parallel, single-cell reporter gene assay that detects cell-type-specific cis-regulatory activity. We then measured a library of promoter variants across multiple cell types in live mouse retinas and showed that subtle genetic variants can produce cell-type-specific effects on cis-regulatory activity. We anticipate that scMPRA will be widely applicable for studying the role of CRSs across diverse cell types.

A mechanism for red coloration in vertebrates.

Red coloration is a salient feature of the natural world. Many vertebrates produce red color by converting dietary yellow carotenoids into red ketocarotenoids via an unknown mechanism. Here, we show that two enzymes, cytochrome P450 2J19 (CYP2J19) and 3-hydroxybutyrate dehydrogenase 1-like (BDH1L), are sufficient to catalyze this conversion. In birds, both enzymes are expressed at the sites of ketocarotenoid biosynthesis (feather follicles and red cone photoreceptors), and genetic evidence implicates these enzymes in yellow/red color variation in feathers. In fish, the homologs of CYP2J19 and BDH1L are required for ketocarotenoid production, and we show that these enzymes are sufficient to produce ketocarotenoids in cell culture and when ectopically expressed in fish skin. Finally, we demonstrate that the red-cone-enriched tetratricopeptide repeat protein 39B (TTC39B) enhances ketocarotenoid production when co-expressed with CYP2J19 and BDH1L. The discovery of this mechanism of ketocarotenoid biosynthesis has major implications for understanding the evolution of color diversity in vertebrates.

Information content differentiates enhancers from silencers in mouse photoreceptors.

Enhancers and silencers often depend on the same transcription factors (TFs) and are conflated in genomic assays of TF binding or chromatin state. To identify sequence features that distinguish enhancers and silencers, we assayed massively parallel reporter libraries of genomic sequences targeted by the photoreceptor TF cone-rod homeobox (CRX) in mouse retinas. Both enhancers and silencers contain more TF motifs than inactive sequences, but relative to silencers, enhancers contain motifs from a more diverse collection of TFs. We developed a measure of information content that describes the number and diversity of motifs in a sequence and found that, while both enhancers and silencers depend on CRX motifs, enhancers have higher information content. The ability of information content to distinguish enhancers and silencers targeted by the same TF illustrates how motif context determines the activity of cis-regulatory sequences.

Different cell types are established by activating and repressing the activity of specific sets of genes, a process controlled by proteins called transcription factors. Transcription factors work by recognizing and binding short stretches of DNA in parts of the genome called cis-regulatory sequences. A cis-regulatory sequence that increases the activity of a gene when bound by transcription factors is called an enhancer, while a sequence that causes a decrease in gene activity is called a silencer. To establish a cell type, a particular transcription factor will act on both enhancers and silencers that control the activity of different genes. For example, the transcription factor cone-rod homeobox (CRX) is critical for specifying different types of cells in the retina, and it acts on both enhancers and silencers. In rod photoreceptors, CRX activates rod genes by binding their enhancers, while repressing cone photoreceptor genes by binding their silencers. However, CRX always recognizes and binds to the same DNA sequence, known as its binding site, making it unclear why some cis-regulatory sequences bound to CRX act as silencers, while others act as enhancers. Friedman et al. sought to understand how enhancers and silencers, both bound by CRX, can have different effects on the genes they control. Since both enhancers and silencers contain CRX binding sites, the difference between the two must lie in the sequence of the DNA surrounding these binding sites. Using retinas that have been explanted from mice and kept alive in the laboratory, Friedman et al. tested the activity of thousands of CRX-binding sequences from the mouse genome. This showed that both enhancers and silencers have more copies of CRX-binding sites than sequences of the genome that are inactive. Additionally, the results revealed that enhancers have a diverse collection of binding sites for other transcription factors, while silencers do not. Friedman et al. developed a new metric they called information content, which captures the diverse combinations of different transcription binding sites that cis-regulatory sequences can have. Using this metric, Friedman et al. showed that it is possible to distinguish enhancers from silencers based on their information content. It is critical to understand how the DNA sequences of cis-regulatory regions determine their activity, because mutations in these regions of the genome can cause disease. However, since every person has thousands of benign mutations in cis-regulatory sequences, it is a challenge to identify specific disease-causing mutations, which are relatively rare. One long-term goal of models of enhancers and silencers, such as Friedman et al.'s information content model, is to understand how mutations can affect cis-regulatory sequences, and, in some cases, lead to disease.

Genetic Basis of De Novo Appearance of Carotenoid Ornamentation in Bare Parts of Canaries.

Unlike wild and domestic canaries (Serinus canaria), or any of the three dozen species of finches in genus Serinus, the domestic urucum breed of canaries exhibits bright red bills and legs. This novel trait offers a unique opportunity to understand the mechanisms of bare-part coloration in birds. To identify the mutation producing the colorful phenotype, we resequenced the genome of urucum canaries and performed a range of analyses to search for genotype-to-phenotype associations across the genome. We identified a nonsynonymous mutation in the gene BCO2 (beta-carotene oxygenase 2, also known as BCDO2), an enzyme involved in the cleavage and breakdown of full-length carotenoids into short apocarotenoids. Protein structural models and in vitro functional assays indicate that the urucum mutation abrogates the carotenoid-cleavage activity of BCO2. Consistent with the predicted loss of carotenoid-cleavage activity, urucum canaries tended to have increased levels of full-length carotenoid pigments in bill tissue and reduced levels of carotenoid-cleavage products (apocarotenoids) in retinal tissue compared with other breeds of canaries. We hypothesize that carotenoid-based bare-part coloration might be readily gained, modified, or lost through simple switches in the enzymatic activity or regulation of BCO2 and this gene may be an important mediator in the evolution of bare-part coloration among bird species.

Cis-regulatory basis of sister cell type divergence in the vertebrate retina.

Multicellular organisms evolved via repeated functional divergence of transcriptionally related sister cell types, but the mechanisms underlying sister cell type divergence are not well understood. Here, we study a canonical pair of sister cell types, retinal photoreceptors and bipolar cells, to identify the key cis-regulatory features that distinguish them. By comparing open chromatin maps and transcriptomic profiles, we found that while photoreceptor and bipolar cells have divergent transcriptomes, they share remarkably similar cis-regulatory grammars, marked by enrichment of K50 homeodomain binding sites. However, cell class-specific enhancers are distinguished by enrichment of E-box motifs in bipolar cells, and Q50 homeodomain motifs in photoreceptors. We show that converting K50 motifs to Q50 motifs represses reporter expression in bipolar cells, while photoreceptor expression is maintained. These findings suggest that partitioning of Q50 motifs within cell type-specific cis-regulatory elements was a critical step in the evolutionary divergence of the bipolar transcriptome from that of photoreceptors.

A massively parallel reporter assay reveals context-dependent activity of homeodomain binding sites in vivo.

Cone-rod homeobox (CRX) is a paired-like homeodomain transcription factor (TF) and a master regulator of photoreceptor development in vertebrates. The in vitro DNA binding preferences of CRX have been described in detail, but the degree to which in vitro binding affinity is correlated with in vivo enhancer activity is not known. In addition, paired-class homeodomain TFs can bind DNA cooperatively as both homodimers and heterodimers at inverted TAAT half-sites separated by 2 or 3 nucleotides. This dimeric configuration is thought to mediate target specificity, but whether monomeric and dimeric sites encode distinct levels of activity is not known. Here, we used a massively parallel reporter assay to determine how local sequence context shapes the regulatory activity of CRX binding sites in mouse photoreceptors. We assayed inactivating mutations in more than 1700 TF binding sites and found that dimeric CRX binding sites act as stronger enhancers than monomeric CRX binding sites. Furthermore, the activity of dimeric half-sites is cooperative, dependent on a strict 3-bp spacing, and tuned by the identity of the spacer nucleotides. Saturating single-nucleotide mutagenesis of 195 CRX binding sites showed that, on average, changes in TF binding site affinity are correlated with changes in regulatory activity, but this relationship is obscured when considering mutations across multiple cis-regulatory elements (CREs). Taken together, these results demonstrate that the activity of CRX binding sites is highly dependent on sequence context, providing insight into photoreceptor gene regulation and illustrating functional principles of homeodomain binding sites that may be conserved in other cell types.

Cell Type-Specific Epigenomic Analysis Reveals a Uniquely Closed Chromatin Architecture in Mouse Rod Photoreceptors.

Rod photoreceptors are specialized neurons that mediate vision in dim light and are the predominant photoreceptor type in nocturnal mammals. The rods of nocturnal mammals are unique among vertebrate cell types in having an 'inverted' nuclear architecture, with a dense mass of heterochromatin in the center of the nucleus rather than dispersed clumps at the periphery. To test if this unique nuclear architecture is correlated with a unique epigenomic landscape, we performed ATAC-seq on mouse rods and their most closely related cell type, cone photoreceptors. We find that thousands of loci are selectively closed in rods relative to cones as well as >60 additional cell types. Furthermore, we find that the open chromatin profile of photoreceptors lacking the rod master regulator Nrl is nearly indistinguishable from that of native cones, indicating that Nrl is required for selective chromatin closure in rods. Finally, we identified distinct enrichments of transcription factor binding sites in rods and cones, revealing key differences in the cis-regulatory grammar of these cell types. Taken together, these data provide insight into the development and maintenance of photoreceptor identity, and highlight rods as an attractive system for studying the relationship between nuclear organization and local changes in gene regulation.

A Simple Grammar Defines Activating and Repressing cis-Regulatory Elements in Photoreceptors.

Transcription factors often activate and repress different target genes in the same cell. How activation and repression are encoded by different arrangements of transcription factor binding sites in cis-regulatory elements is poorly understood. We investigated how sites for the transcription factor CRX encode both activation and repression in photoreceptors by assaying thousands of genomic and synthetic cis-regulatory elements in wild-type and Crx-/- retinas. We found that sequences with high affinity for CRX repress transcription, whereas sequences with lower affinity activate. This rule is modified by a cooperative interaction between CRX sites and sites for the transcription factor NRL, which overrides the repressive effect of high affinity for CRX. Our results show how simple rearrangements of transcription factor binding sites encode qualitatively different responses to a single transcription factor and explain how CRX plays multiple cis-regulatory roles in the same cell.

Massively parallel cis-regulatory analysis in the mammalian central nervous system.

Cis-regulatory elements (CREs, e.g., promoters and enhancers) regulate gene expression, and variants within CREs can modulate disease risk. Next-generation sequencing has enabled the rapid generation of genomic data that predict the locations of CREs, but a bottleneck lies in functionally interpreting these data. To address this issue, massively parallel reporter assays (MPRAs) have emerged, in which barcoded reporter libraries are introduced into cells, and the resulting barcoded transcripts are quantified by next-generation sequencing. Thus far, MPRAs have been largely restricted to assaying short CREs in a limited repertoire of cultured cell types. Here, we present two advances that extend the biological relevance and applicability of MPRAs. First, we adapt exome capture technology to instead capture candidate CREs, thereby tiling across the targeted regions and markedly increasing the length of CREs that can be readily assayed. Second, we package the library into adeno-associated virus (AAV), thereby allowing delivery to target organs in vivo. As a proof of concept, we introduce a capture library of about 46,000 constructs, corresponding to roughly 3500 DNase I hypersensitive (DHS) sites, into the mouse retina by ex vivo plasmid electroporation and into the mouse cerebral cortex by in vivo AAV injection. We demonstrate tissue-specific cis-regulatory activity of DHSs and provide examples of high-resolution truncation mutation analysis for multiplex parsing of CREs. Our approach should enable massively parallel functional analysis of a wide range of CREs in any organ or species that can be infected by AAV, such as nonhuman primates and human stem cell-derived organoids.

Redesign of the monomer-monomer interface of Cre recombinase yields an obligate heterotetrameric complex.

Cre recombinase catalyzes the cleavage and religation of DNA at loxP sites. The enzyme is a homotetramer in its functional state, and the symmetry of the protein complex enforces a pseudo-palindromic symmetry upon the loxP sequence. The Cre-lox system is a powerful tool for many researchers. However, broader application of the system is limited by the fixed sequence preferences of Cre, which are determined by both the direct DNA contacts and the homotetrameric arrangement of the Cre monomers. As a first step toward achieving recombination at arbitrary asymmetric target sites, we have broken the symmetry of the Cre tetramer assembly. Using a combination of computational and rational protein design, we have engineered an alternative interface between Cre monomers that is functional yet incompatible with the wild-type interface. Wild-type and engineered interface halves can be mixed to create two distinct Cre mutants, neither of which are functional in isolation, but which can form an active heterotetramer when combined. When these distinct mutants possess different DNA specificities, control over complex assembly directly discourages recombination at unwanted half-site combinations, enhancing the specificity of asymmetric site recombination. The engineered Cre mutants exhibit this assembly pattern in a variety of contexts, including mammalian cells.

Massively parallel in vivo enhancer assay reveals that highly local features determine the cis-regulatory function of ChIP-seq peaks.

Transcription factors (TFs) recognize short sequence motifs that are present in millions of copies in large eukaryotic genomes. TFsmust distinguish their target binding sites from a vast genomic excess of spurious motif occurrences; however, it is unclear whether functional sites are distinguished from nonfunctional motifs by local primary sequence features or by the larger genomic context in which motifs reside. We used a massively parallel enhancer assay in living mouse retinas to compare 1,300 sequences bound in the genome by the photoreceptor transcription factor Cone-rod homeobox (Crx), to 3,000 control sequences. We found that very short sequences bound in the genome by Crx activated transcription at high levels, whereas unbound genomic regions with equal numbers of Crx motifs did not activate above background levels, even when liberated from their larger genomic context. High local GC content strongly distinguishes bound motifs from unbound motifs across the entire genome. Our results show that the cis-regulatory potential of TF-bound DNA is determined largely by highly local sequence features and not by genomic context.

Reprogramming of adult rod photoreceptors prevents retinal degeneration.

A prime goal of regenerative medicine is to direct cell fates in a therapeutically useful manner. Retinitis pigmentosa is one of the most common degenerative diseases of the eye and is associated with early rod photoreceptor death followed by secondary cone degeneration. We hypothesized that converting adult rods into cones, via knockdown of the rod photoreceptor determinant Nrl, could make the cells resistant to the effects of mutations in rod-specific genes, thereby preventing secondary cone loss. To test this idea, we engineered a tamoxifen-inducible allele of Nrl to acutely inactivate the gene in adult rods. This manipulation resulted in reprogramming of rods into cells with a variety of cone-like molecular, histologic, and functional properties. Moreover, reprogramming of adult rods achieved cellular and functional rescue of retinal degeneration in a mouse model of retinitis pigmentosa. These findings suggest that elimination of Nrl in adult rods may represent a unique therapy for retinal degeneration.

Quantifying the activity of cis-regulatory elements in the mouse retina by explant electroporation.

Transcription factors control gene expression by binding to noncoding regions of DNA known as -cis-regulatory elements (CREs; i.e., enhancer/promoters). Traditionally, cis-regulatory analysis has been carried out via mouse transgenesis which is time-consuming and nonquantitative. Electroporation of DNA reporter constructs into living mouse tissue is a rapid and effective alternative to transgenesis which permits quantitative assessment of cis-regulatory activity. Here, we present a simple technique for quantifying the activity of photoreceptor-specific CREs in living explanted mouse retinas.

Complex effects of nucleotide variants in a mammalian cis-regulatory element.

Cis-regulatory elements (CREs) control gene expression by recruiting transcription factors (TFs) and other DNA binding proteins. We aim to understand how individual nucleotides contribute to the function of CREs. Here we introduce CRE analysis by sequencing (CRE-seq), a high-throughput method for producing and testing large numbers of reporter genes in mammalian cells. We used CRE-seq to assay >1,000 single and double nucleotide mutations in a 52-bp CRE in the Rhodopsin promoter that drives strong and specific expression in mammalian photoreceptors. We find that this particular CRE is remarkably complex. The majority (86%) of single nucleotide substitutions in this sequence exert significant effects on regulatory activity. Although changes in the affinity of known TF binding sites explain some of these expression changes, we present evidence for complex phenomena, including binding site turnover and TF competition. Analysis of double mutants revealed complex, nucleotide-specific interactions between residues in different TF binding sites. We conclude that some mammalian CREs are finely tuned by evolution and function through complex, nonadditive interactions between bound TFs. CRE-seq will be an important tool to uncover the rules that govern these interactions.

Exome sequencing and cis-regulatory mapping identify mutations in MAK, a gene encoding a regulator of ciliary length, as a cause of retinitis pigmentosa.

A fundamental challenge in analyzing exome-sequence data is distinguishing pathogenic mutations from background polymorphisms. To address this problem in the context of a genetically heterogeneous disease, retinitis pigmentosa (RP), we devised a candidate-gene prioritization strategy called cis-regulatory mapping that utilizes ChIP-seq data for the photoreceptor transcription factor CRX to rank candidate genes. Exome sequencing combined with this approach identified a homozygous nonsense mutation in male germ cell-associated kinase (MAK) in the single affected member of a consanguineous Turkish family with RP. MAK encodes a cilium-associated mitogen-activated protein kinase whose function is conserved from the ciliated alga, Chlamydomonas reinhardtii, to humans. Mutations in MAK orthologs in mice and other model organisms result in abnormally long cilia and, in mice, rapid photoreceptor degeneration. Subsequent sequence analyses of additional individuals with RP identified five probands with missense mutations in MAK. Two of these mutations alter amino acids that are conserved in all known kinases, and an in vitro kinase assay indicates that these mutations result in a loss of kinase activity. Thus, kinase activity appears to be critical for MAK function in humans. This study highlights a previously underappreciated role for CRX as a direct transcriptional regulator of ciliary genes in photoreceptors. In addition, it demonstrates the effectiveness of CRX-based cis-regulatory mapping in prioritizing candidate genes from exome data and suggests that this strategy should be generally applicable to a range of retinal diseases.

Quantifying the activity of cis-regulatory elements in the mouse retina by explant electroporation.

Transcription factors within cellular gene networks control the spatiotemporal pattern and levels of expression of their target genes by binding to cis-regulatory elements (CREs), short (˜300-600 bp) stretches of genomic DNA which can lie upstream, downstream, or within the introns of the genes they control. CREs (i.e., enhancers/promoters) typically consist of multiple clustered binding sites for both transcriptional activators and repressors(1-3). They serve as logical integrators of transcriptional input giving a unitary output in the form of spatiotemporally precise and quantitatively exact promoter activity. Most studies of mammalian cis-regulation to date have relied on mouse transgenesis as a means of assaying the enhancer function of CREs(4-5). This technique is time-consuming, costly and, on account of insertion site effects, largely non-quantitative. On the other hand, quantitative assays for mammalian CRE function have been developed in tissue culture systems (e.g., dual luciferase assays), but the in vivo relevance of these results is often uncertain. Electroporation offers an excellent alternative to traditional mouse transgenesis in that it permits both spatiotemporal and quantitative assessment of cis-regulatory activity in living mammalian tissue. This technique has been particularly useful in the analysis of cis-regulation in the central nervous system, especially in the cerebral cortex and the retina(6-8). While mouse retinal electroporation, both in vivo and ex vivo, has been developed and extensively described by Matsuda and Cepko(6-7,9), we have recently developed a simple approach to quantify the activity of photoreceptor-specific CREs in electroporated mouse retinas(10). Given that the amount of DNA that is introduced into the retina by electroporation can vary from experiment to experiment, it is necessary to include a co-electroporated 'loading control' in all experiments. In this respect, the technique is very similar to the dual luciferase assay used to quantify promoter activity in cultured cells. When assaying photoreceptor cis-regulatory activity, electroporation is usually performed in newborn mice (postnatal day 0, P0) which is the time of peak rod production(11-12). Once retinal cell types become post-mitotic, electroporation is much less efficient. Given the high rate of rod birth in newborn mice and the fact that rods constitute more than 70% of the cells in the adult mouse retina, the majority of cells that are electroporated at P0 are rods. For this reason, rod photoreceptors are the easiest retinal cell type to study via electroporation. The technique we describe here is primarily useful for quantifying the activity of photoreceptor CREs.

Risk- and non-risk-associated variants at the 10q26 AMD locus influence ARMS2 mRNA expression but exclude pathogenic effects due to protein deficiency.

Fifteen variants in 10q26 are in strong linkage disequilibrium and are associated with an increased risk for age-related macular degeneration (AMD), a frequent cause of blindness in developed countries. These variants tag a single-risk haplotype encompassing the genes ARMS2 (age-related maculopathy susceptibility 2) and part of HTRA1 (HtrA serine peptidase 1). To define the true AMD susceptibility gene in 10q26, several studies have focused on the influence of risk alleles on the expression of ARMS2 and/or HTRA1, but the results have been inconsistent. By heterologous expression of genomic ARMS2 variants, we now show that ARMS2 mRNA levels transcribed from the risk haplotype are significantly reduced compared with non-risk mRNA isoforms. Analyzing variant ARMS2 constructs, this effect could specifically be assigned to the known insertion/deletion polymorphism (c.(*)372_815del443ins54) in the 3'-untranslated region of ARMS2. Reporter gene assays with HTRA1 promoter sequences demonstrated the presence of a Müller glia-specific cis-regulatory region further upstream of the transcription start site. However, AMD risk alleles had little or no effect on HTRA1 promoter activity in the retina. Analysis of a large series of human post-mortem retina/retinal pigment epithelial samples heterozygous for the risk haplotype confirmed the in vitro/ex vivo results and demonstrated that the risk haplotype affects ARMS2 but not HTRA1 mRNA expression. Furthermore, we provide in vivo evidence that a common non-risk-associated non-synonymous variant (rs2736911) also leads to decreased ARMS2 transcript levels. Consequently, our data suggest that pathogenic effects due to ARMS2 protein deficiency are unlikely to account for AMD pathology.

CRX ChIP-seq reveals the cis-regulatory architecture of mouse photoreceptors.

Approximately 98% of mammalian DNA is noncoding, yet we understand relatively little about the function of this enigmatic portion of the genome. The cis-regulatory elements that control gene expression reside in noncoding regions and can be identified by mapping the binding sites of tissue-specific transcription factors. Cone-rod homeobox (CRX) is a key transcription factor in photoreceptor differentiation and survival, but its in vivo targets are largely unknown. Here, we used chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) on CRX to identify thousands of cis-regulatory regions around photoreceptor genes in adult mouse retina. CRX directly regulates downstream photoreceptor transcription factors and their target genes via a network of spatially distributed regulatory elements around each locus. CRX-bound regions act in a synergistic fashion to activate transcription and contain multiple CRX binding sites which interact in a spacing- and orientation-dependent manner to fine-tune transcript levels. CRX ChIP-seq was also performed on Nrl(-/-) retinas, which represent an enriched source of cone photoreceptors. Comparison with the wild-type ChIP-seq data set identified numerous rod- and cone-specific CRX-bound regions as well as many shared elements. Thus, CRX combinatorially orchestrates the transcriptional networks of both rods and cones by coordinating the expression of photoreceptor genes including most retinal disease genes. In addition, this study pinpoints thousands of noncoding regions of relevance to both Mendelian and complex retinal disease.