Genome-wide Occupancy of NRL revealed by ChIP–Seq using Illumina and ABI/SOLiD sequencing platforms.

(A) Analysis workflow. Raw sequence reads from Illumina or ABI/SOLiD were mapped to the mouse genome (NCBI build 37) using the Genomatix Mining Station (GMS) and the reads mapped to unique genomic locations (uniquely mapped reads) were used for further analyses. ChIP–Seq peaks were called using NGS Analyzer (Genomatix) or MACS (Zhang et al., 2008), and the common peaks were used for further analyses. The NRL ChIP–Seq peaks were compared to the CRX ChIP–Seq peaks for overlapping using GenomeInspector (Genomatix) software. The ChIP–Seq peaks were assigned to the nearest gene. Transcription profile analyses of flow-sorted photoreceptors of WT and Nrl−/− were performed using ChipInsepector program (Genomatix) and 1.5 fold expression change was used as a criterion for NRL target genes. TF motif enrichment analyses were performed on the NRL ChIP–Seq peak regions that were associated with NRL target genes. Comparison was made between CRX-overlapping and non CRX-overlapping NRL ChIP–Seq peaks. Gene regulatory network was constructed based on TF enrichment analysis. (B) Correlation of ChIP–Seq peaks by Illumina and ABI. The number of correlations (y-axis) was plotted to the distance of ABI ChIP–Seq peaks to Illumina ChIP–Seq peaks (x-axis). The Venn diagram (inset) calculated the percentage of ABI and Illumina peaks within 500 bp of each other: 88% of Illumina peaks are within 500 bp of ABI peaks and 49% of ABI peaks are within 500 bp of Illumina peaks. (C) Correlation of ChIP–Seq peaks to promoters. The number of correlation (y-axis) was plotted to the distance of ABI ChIP–Seq peaks (green graph) or Illumina ChIP–Seq peaks (blue graph) to the transcription start site (TSS) (x-axis). The Venn diagram (inset) calculated the percentage of ABI (75%) or Illumina (72%) peaks within 10,000 bp from the TSS. (D) Genomic distribution of NRL ChIP–Seq peaks relative to the nearest annotated genes. Promoters and exons account for 2.3% and 5.4% of the mouse genome, respectively.

(A) Identification of direct NRL transcriptional target genes by ChIP–Seq and transcription profiling. ABI ChIP–Seq peaks and Illumina ChIP–Seq peaks were assigned to the nearest genes (ABI and Illumina). Transcriptional profiling of flow-sorted photoreceptors of WT or Nrl−/− mice was generated using microarrays. Up: genes up-regulated in Nrl−/− photoreceptors. Down: genes down-regulated in Nrl−/− photoreceptors. (B) Correlation of NRL ChIP–Seq peaks to the promoters of its target genes. The number of correlation (y-axis) was plotted to the distance (x-axis) of ABI ChIP–Seq peaks (green graph) or Illumina ChIP–Seq peaks (blue graph) to the promoters of the genes that were down-regulated (top) or up-regulated (bottom) in the Nrl−/− photoreceptors. (C) TF enrichment and TF binding site (TFBS) positional bias analysis. NRL ChIP–Seq peak regions were analyzed for TF enrichment using Genomatix RegionMiner. The positional bias of TFBS (P) was calculated and plotted as −log(P) (y-axis) to the distance of TFBS to the peak center (x-axis). Positions where TFBS are overrepresented appear as peaks in these plots. * significantly enriched. (D) Correlation of NRL ChIP–Seq peaks with CRX ChIP–Seq peaks. The number of correlation (y-axis) was plotted to the distance of NRL ChIP–Seq peaks to CRX ChIP–Seq peaks (x-axis). The Venn diagram (inset) calculated the percentage of NRL ChIP–Seq peaks (Illumina and ABI) within 500 bp of CRX peaks: 65% of Illumina peaks and 48% of ABI peaks are within 500 bp of CRX peaks.

ChIP-qPCR was performed to validate NRL binding to 26 ChIP–Seq peak regions (left panel), and 5 non-peak regions (right panel) served as negative controls. The amount of ChIP DNA was measured by qPCR in triplicates using primers flanking the regions of interest. Normal IgG served as an antibody control when ChIP was performed using WT retinas (white bars). White bars (NRL Ab/IgG) represent fold change (FC) of qPCR signals comparing NRL ChIP DNA to the IgG control ChIP DNA. A separate set of ChIP assays was performed using NRL antibody to compare signals from WT retina to signals from Nrl−/− retina (tissue control). Black bars (WT/Nrl−/−) represent fold increase (Fc) of qPCR signals comparing NRL ChIP DNA from wild type C57BL/6 mouse retina to NRL ChIP DNA from Nrl−/− mouse retina. The ChIP-qPCR assays were performed twice. The representative results were shown as mean ± SD. P<0.01 for all by Student's t test.

Twenty-six NRL peak regions were cloned into pGL3-promoter vector in front of a SV40 basal promoter and a luciferase reporter. The constructs were transfected in HEK293T cells together with mouse Nrl (mNrl) expression plasmid (in pC4C vector) or empty pC4C vector. The y-axis is fold change (Fc) of normalized luciferase readings. Control: enhancer constructs co-transfected with empty pC4C vector. Five non-peak regions served as additional negative controls. The experiments were performed three times. Representative results are shown as mean ± SD. P<0.05 for all by Student's t test.

Microarray analysis of flow-sorted photoreceptors after in vivo Kdm5b knockdown.

(A) Experimental workflow. CD-1 mouse retinas were transfected at P0 with Ub-GFP and shRNA against Gapdh or Kdm5b by sub-retinal injection and in vivo electroporation. Retinas were dissociated at P20 and shRNA-transfected retinal cells were isolated by flow-sorting. The effect of Kdm5b shRNA on transcriptional profile was measured by microarray analysis. (B) Ontology analysis of common targets of KDM5b and NRL. Kdm5b shRNA up: genes that are up-regulated in retinal cells electroporated with Kdm5b shRNA. Kdm5b shRNA down: genes that are down-regulated in retinal cells electroporated with Kdm5b shRNA. Nrl-ko up: genes that are up-regulated in Nrl−/− photoreceptor cells. Nrl-ko down: genes that are down-regulated in Nrl−/− photorceptor cells.