Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
[Preprint]. 2023 Dec 2:2023.05.27.542576.
doi: 10.1101/2023.05.27.542576.

Pathogenic variants in Crx have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

James L Shepherdson  1   2 , Ryan Z Friedman  1   2 , Yiqiao Zheng  3 , Chi Sun  3 , Inez Y Oh  3 , David M Granas  1   2 , Barak A Cohen  1   2 , Shiming Chen  3   4 , Michael A White  1   2
Affiliations

Pathogenic variants in Crx have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

James L Shepherdson et al. bioRxiv. .

Update in

Abstract

Dozens of variants in the photoreceptor-specific transcription factor (TF) CRX are linked with human blinding diseases that vary in their severity and age of onset. It is unclear how different variants in this single TF alter its function in ways that lead to a range of phenotypes. We examined the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in live mouse retinas carrying knock-ins of two variants, one in the DNA binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation caused by the variants corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, while p.E168d2 has stronger effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are de-repressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci were partially predictive of episomal MPRA activity, and silencers were notably enriched among distal elements whose accessibility increases later in retinal development. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers, while having a qualitatively different impact on silencers.

PubMed Disclaimer

Conflict of interest statement

Competing Interest Statement B.A.C is on the scientific advisory board of Patch Biosciences. The authors declare no other competing interests.

Figures

Figure 1:
Figure 1:
(A) Schematics of wild-type CRX protein and the p.R90W and p.E168d2 pathogenic variants, containing a point DNA binding domain mutant or resulting in truncation of the transcriptional effector domain, respectively. (B) Outline of the experimental procedure for constructing and testing the CRE libraries. A library of CRE sequences is cloned upstream of a rhodopsin (“Rho”) driving expression of the dsRed fluorescent protein, with each CRE marked by a unique sequence barcode (“BC”) in the 3’ UTR. Each library was electroporated into retinal explants with the six indicated Crx genotypes. RNA was collected from the retinas, and transcript counts for each element were measured and normalized to abundance in the input DNA library to calculate transcriptional activity.
Figure 2:
Figure 2:
(A,B) Correlations of the transcriptional activity of each CRE with intact CRX motifs (A) or CRE with killed CRX motifs (B) in the indicated genotypes. Superimposed black cross symbols indicate the transcriptional activity of the 100 scrambled control sequences, and the yellow circle shows the activity of the basal rhodopsin promoter alone. Using the scrambled control sequences to establish an empirical null distribution of CRE activity, we performed robust regression using a trimmed mean M-estimator (fit plotted as a dashed line) for each genotype comparison. The outer dotted lines indicate the 80th percentile of the residuals of the library CREs relative to the regression of the scrambled sequences. The outer vertical and horizontal lines indicate the thresholds for the “strong silencer” and “strong enhancer” classes, while the inner horizontal and vertical lines correspond to the basal rhodopsin promoter alone. (C) Activity of the basal promoter construct in each of the indicated genotypes. Each black dot indicates a different unique sequence barcode; the mean across barcodes is shown as a yellow dot. In all panels, transcriptional activity was adjusted for visualization purposes using a fixed pseudocount of 2e−5; the out-of-bound region is indicated by the light gray stripes.
Figure 3:
Figure 3:
(A) Distribution of activity classes in +/+ retina for 66 and 117 CREs near cone- or rod-enriched genes, respectively (Ruzycki et al. 2015). (B) Fold change relative to +/+ of 66 CREs near cone-enriched genes. (C) Fold change relative to +/+ of 117 CREs near rod-enriched genes. In (B) and (C), rows are sorted by fold change in −/−.
Figure 4:
Figure 4:
(A) Heatmap comparing the transcriptional activity of each CRE across the six CRX genotypes. In the top right (upper arrow), note a group of CREs with generally low activity in all but the E168d2/d2 and −/− genotypes. Along the bottom (lower arrow) note the stepwise reduction in activity with increasing phenotypic severity. CREs are sorted by the ratio of transcriptional activity in +/+ and −/− genotypes. (B) Classifications across genotypes for CREs classified as strong silencers (left column of each heatmap) or weak silencers (right column) in wildtype retinas. (C) Transcriptional activity of wildtype or motif mutant CREs across genotypes for 234 silencers in +/+ that convert to enhancers in E168d2/d2 (“Converted”), or the remaining 251 silencers in +/+ that do not convert to enhancers in E168d2/d2 (“Unconverted”). The yellow dot indicates the transcriptional activity of the basal promoter alone in each genotype. (D) Predicted CRX occupancy of the same silencer subgroups as in panel C and CREs that act as enhancers in +/+. (E) Information content of the same subgroups as panel D. *: p < 0.05, **: p < 0.01, ****: p < 0.0001 via two-sided Mann-Whitney-Wilcoxon test. In panel C, transcriptional activity was adjusted for visualization purposes using a fixed pseudocount of 2e-5; the out-of-bound region is indicated by the light gray stripes.
Figure 5:
Figure 5:
(A) Distribution of CRE activity classes in +/+ retina separated into groups based on chromatin annotations of the site of genomic origin and CRX-dependent ATAC accessibility (Ruzycki et al. 2018). Group A: H3K4me3+ H3K27Ac+, Group C: H3K4me3- H3K27Ac+, Group D: H3K4me3- H3K27Ac-. (B) Transcriptional activity of all CREs, separated by chromatin and ATAC group, across all genotypes. (C) Chromatin classification breakdown of E168d2/d2 converted silencers (elements classified as weak or strong silencers in +/+ that become weak or strong enhancers in E168d2/d2). In panel B, transcriptional activity was adjusted for visualization purposes using a fixed pseudocount of 2e−5; the out-of-bound region is indicated by the light gray stripes.
Figure 6:
Figure 6:
A model of differential enhancer and silencer sensitivity to CRX mutants. Silencers, with higher predicted CRX occupancy, engage in effector-domain mediated homotypic interactions that permit CRX to remain bound even when the DBD is weakened by the p.R90W variant. The p.E168d2 variant lacks the majority of the effector domain and cannot engage in these interactions. Enhancers have lower CRX predicted occupancy and engage in fewer homotypic interactions.
See this image and copyright information in PMC

Similar articles

See all similar articles

References

    1. Bailey TL, Grant CE. 2021. SEA: Simple Enrichment Analysis of motifs. bioRxiv.
    1. Berger W, Kloeckener-Gruissem B, Neidhardt J. 2010. The molecular basis of human retinal and vitreoretinal diseases. Progress in Retinal and Eye Research 29: 335–375. - PubMed
    1. Chau K-Y, Chen S, Zack DJ, Ono SJ. 2000. Functional Domains of the Cone-Rod Homeobox (CRX) Transcription Factor. Journal of Biological Chemistry 275: 37264–37270. - PubMed
    1. Chen S, Wang Q-L, Nie Z, Sun H, Lennon G, Copeland NG, Gilbert DJ, Jenkins NA, Zack DJ. 1997. Crx, a Novel Otx-like Paired-Homeodomain Protein, Binds to and Transactivates Photoreceptor Cell-Specific Genes. Neuron 19: 1017–1030. - PubMed
    1. Chen S, Wang Q-L, Xu S, Liu I, Li LY, Wang Y, Zack DJ. 2002. Functional analysis of conerod homeobox (CRX) mutations associated with retinal dystrophy. Human Molecular Genetics 11: 873–884. - PubMed

Publication types