|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 40, 42211-42220, October 1, 2004
Altered Expression of Genes of the Bmp/Smad and Wnt/Calcium Signaling Pathways in the Cone-only Nrl-/- Mouse Retina, Revealed by Gene Profiling Using Custom cDNA Microarrays*From the Departments of aOphthalmology and Visual Sciences, bBiomedical Engineering, eBiostatistics, and hHuman Genetics, University of Michigan, Ann Arbor, Michigan 48015, the dTranslational Research Center, Kyoto University Hospital, Sakyo-ku, Kyoto 606-8507, Japan, and the gLaboratoire de Physiopathologie Cellulaire et Moléculaire de la Rétine, INSERM U.592, 67091 Strasbourg Cedex, France
Received for publication, July 20, 2004
Many mammalian retinas are rod-dominant, and hence our knowledge of cone photoreceptor biology is relatively limited. To gain insights into the molecular differences between rods and cones, we compared the gene expression profile of the rod-dominated retina of wild type mouse with that of the cone-only retina of Nrl-/- (Neural retina leucine zipper knockout) mouse. Our analysis, using custom microarrays of eye-expressed genes, provided equivalent data using either direct or reference-based experimental designs, confirmed differential expression of rod- and cone-specific genes in the Nrl-/- retina and identified novel genes that could serve as candidates for retinopathies or for functional studies. In addition, we detected altered expression of several genes that encode cell signaling or structural proteins. Prompted by these findings, additional real-time PCR analysis revealed that genes belonging to the Bmp/Smad and Wnt/Ca2+ signaling pathways are expressed in the mature wild type retina and that their expression is significantly altered in the Nrl-/- retina. Chromatin immunoprecipitation analysis of adult retina identified Bmp4 and Smad4, which are down-regulated in the Nrl-/- retina, as possible direct transcriptional targets of Nrl. Consistent with these studies, Bmp4 and Smad4 are expressed in the mature rod photoreceptors of mouse retina. Modulation of Bmp4 and/or Smad4 by Nrl may provide a mechanism for integrating diverse cell signaling networks in rods. We hypothesize that Bmp/Smad and Wnt/Ca2+ pathways participate in cell-cell communication in the mature retina, and expression changes observed in the Nrl-/- retina reflect their biased utilization in rod versus cone homeostasis.
In mammals, vision is initiated in the retina, which is a highly structured part of the brain consisting of over 50 types of neurons that are organized in three distinct layers (1, 2). The outer nuclear layer consists exclusively of two types of photoreceptors, rods and cones, responsible for detection and transduction of light energy. Rods function under low ambient light and form the major photoreceptor population of many mammals, including humans (95%) and mice (97%). Cones are responsible for phototransduction in bright light, providing high acuity and color vision (1). Cones are needed for maintaining central vision, and their survival is vital for preserving visual function in retinal and macular diseases (3, 4). Rods and cones possess distinct subsets of proteins involved in the phototransduction cascade (5, 6), but despite significant neuroanatomical and physiological advances (2, 7, 8), little progress has been made toward delineating the molecular mechanisms that underlie functional distinctions between the two photoreceptor types, their communication with other neurons, and their maintenance, survival, or remodeling in response to extrinsic or intrinsic insults. One approach to systematically dissect the regulatory networks and molecules associated with rod or cone photoreceptor function is to take advantage of animal models that exhibit preferential utilization of one or the other photoreceptor sub-type. Unfortunately, many species with cone-rich retinas (e.g. ground squirrel and chick) present difficulties with respect to experimentation and are less amenable to genetic manipulations. Because of a large number of naturally occurring and experimentally generated mutants (available at jaxmice.jax.org/info/index.html), the mouse offers a unique opportunity to examine this complex question. The Nrl1-knockout (Nrl-/-) mouse, recently generated in our laboratory, exhibits a unique rod-less and cone-only retinal phenotype (9). Nrl was originally identified from a subtracted human retina library and shown to be expressed, by Northern analysis, specifically in the retina and retinoblastoma cell lines (10). In situ hybridization analysis identified Nrl transcripts in developing mouse brain and lens, although the expression became restricted to the retina after birth (11). Later studies, however, demonstrated that Nrl is specifically and highly expressed in the rod (and not cone) photoreceptors (9, 12) and pineal gland2 and that the transcripts in developing brain and lens probably represented cross-hybridization with p45, L-Maf, or another homologous sequence in the mouse genome (12, 13). Nrl is shown to interact with other transcription factors, such as the homeodomain protein Crx, zinc finger protein Fiz1, and orphan nuclear receptor Nr2e3, and regulate (either alone or synergistically) the expression of several rod-specific genes (1420). Missense mutations in the human NRL gene are associated with autosomal dominant retinitis pigmentosa (21, 22). Consistent with these findings, the targeted deletion of Nrl (Nrl-/-) in mouse resulted in a retina with no rod photoreceptors; instead, a concomitant increase in functional S-opsin expressing cones was observed (9). This apparent functional switching of photoreceptor sub-types (from rods to cones) has been validated by histology, electrophysiology, and biochemical and molecular analysis (9, 23, 24).3
Microarray-based global profiling of gene expression, in combination with bioinformatic tools, can yield valuable insights into cell- or tissue-specific functions. Expression profiling of tissues from mice deficient in a transcription factor gene can point to downstream regulatory targets, provide candidates for functional studies, and facilitate positional cloning of human disease loci (2528). Analysis of the retinal transcriptome during development and aging and in mouse models of retinal dysfunction has been the subject of intense investigation (2831). We have used the Nrl-/- mouse model to identify molecular differences between rods and cones (24). However, commercially available microarrays do not have adequate representation of genes transcribed in developing and mature eye/retina; hence, several groups have produced custom slide microarrays of eye/retina-expressed genes (28, 30, 32). For gene profiling, we have isolated and sequenced cDNAs from mouse eye/retina libraries, annotated over 10,000 expressed sequence tags (ESTs), and produced cDNA microarrays (called I-gene microarrays) (33, 34). Here, we report the expression profile of the mature retina from the rod-less (and cone-only) Nrl-/- mice using I-gene microarrays and compare it to the gene profile of the rod-dominated wild type mouse retina. We demonstrate differential expression of several genes, encoding phototransduction, structural, and signaling proteins, in the Nrl-/- retina. Our data reveal novel differentially expressed genes for future functional studies. Of particular importance are the findings that genes encoding components of Bmp/Smad and Wnt/Ca2+ signaling pathways are expressed in the mature retina and their expression is altered in the Nrl-/- retina. Expression analysis by real-time PCR and chromatin immunoprecipitation studies suggest that the activity of the Smad-mediated Bmp signaling pathway is modulated by Nrl in the mature retina. In support of this, we show that Bmp4 and Smad4 are expressed in the rod photoreceptors of the mature mouse retina. We also propose that rods and cones exhibit selective bias in the utilization of different signaling pathways for cell-cell communication and controlling intracellular functions.
Tissue Preparation, RNA Isolation, and Northern AnalysisAll procedures involving mice were approved by the University Committee on Use and Care of Animals of the University of Michigan. Retinas were dissected from the wild type and Nrl-/- mice at postnatal day (P) 21 and snap-frozen on dry ice. Total RNA was isolated using TRIzol reagent (Invitrogen) and purified by using an RNeasy kit (Qiagen). RNA integrity was verified by denaturing formaldehyde-agarose gels. Total RNA samples with a 260- to 280-nm absorbance ratio of greater than 1.9 were used for studies. Northern analysis was performed as described previously (10). Reference RNA for Retinal Gene ProfilingTo generate the reference RNA for microarray hybridizations, total RNA was pooled from the following tissues and cell lines: mouse eye or retina at embryonic day (E) 1416, P23, P1012, and adult (5.5 mg); mouse adult brain (1.5 mg); P19 embryonic carcinoma cells (35) (3 mg); retinoic acid-induced P19 cells that differentiated into neuronal and glial-like cells (3 mg); and neuroblastoma cell line N1E-115 (36) (4 mg). The pooled RNA was divided into aliquots (at a concentration of 1.28 µg/µl) and stored at -80 °C until use. Target Labeling and Microarray HybridizationMouse I-gene cDNA microarray slides (34) contained PCR-amplified products from over 6,500 cDNAs, printed in duplicate; cDNAs were isolated from libraries constructed from E15.5 eyes, P2 eyes, and adult retinas, sequenced, and annotated (33). Target RNA (10 µg of total RNA) was labeled using a 3DNA Submicro Expression Array Detection kit (Genisphere, Hatfield, PA) and hybridized to microarray slides, as described (37). The slides were scanned using an Affymetrix 428 scanner (Affymetrix, Santa Clara, CA) to obtain the highest intensity of signal, without reaching saturation for a maximum of 10 out of the 13,440 spots. Image Processing and Data AnalysisScanned array images without major defects (such as scratches or blobs) were analyzed using AnalyzerDG (MolecularWare Inc., Cambridge, MA) in a batch mode. A "contour shape" was utilized to detect spots for intensity calculations, whereas a "cell method" was set to calculate background on an individual basis in a local square region centered on the spot. A data file containing spots' intensities and annotations was exported for each array in the tab-delimited text format and then imported into the statistical package, R (available at www.r-project.org/). Intensities of Cy3 and Cy5 channels for each spot were calculated separately by subtracting corresponding median background from the mean signal intensity. Genes with negative background-corrected intensities in either channel were filtered out. To reduce systematic variation caused by experimental procedures (such as dye effects), a data-driven normalization was applied to individual datasets using a cluster of least-altered genes on the array identified by a rank-based algorithm. Briefly, let (Rj, Gj) denote the measurements for the jth gene in the red (Cy5) and green (Cy3) channels, respectively, j = 1,..., m. We calculate the ranks in each channel separately, take the difference in ranks between the two channels, and fit a three-component normal mixture model to the difference in ranks for the m genes. There will be three classes of genes to consider: those for which the rank in the red channel is significantly higher than that in the green channel, those for which the rank in the green channel is significantly higher than that in the red channel, and those for which the ranks do not substantially change between the two channels. The genes whose ranks did not change between the two channels were used to perform a slide-dependent normalization based on a locally weighted linear squares procedure (38). This allowed us to normalize genes and redefine a new horizontal zero axis that was used to compute log ratios of intensities in two channels. For indirect comparisons, log ratios of the wild type slide were subtracted from the Nrl-/- slide to obtain reference-corrected values for each gene. An Empirical Bayes method was then applied to the replicated arrays to obtain a B statistic for each gene (39); the B statistic estimates the posterior log odds of differential expression. Quantitative Reverse Transcription-PCRTo validate gene expression changes, qRT-PCR analysis was performed as described (33). Briefly, total RNA (2.5 µg) treated with RQ1 RNase-free DNase (Promega, Madison, WI) was subjected to reverse transcription using oligod(T) primers and with (+RT sample) or without (-RT) SuperScript II (Invitrogen). For each gene, real-time PCR reactions from wild type and Nrl-/- samples were performed in triplicate with SybrGreen I (Molecular Probes, Eugene, OR) and analyzed with the iCycler IQ real-time PCR detection system (Bio-Rad). The average threshold cycle (Ct) differences between the two samples were normalized against Hprt (control) in the corresponding cDNA preparation.
ImmunohistochemistryAnti-Gnb3 polyclonal antibodies were generated in rabbit against the peptide (ADITLAELVSGLEVV) and affinity-purified (Invitrogen). Anti-Smad4 rabbit polyclonal antibody (H-552) was purchased from Santa Cruz Biotechnology Inc. (Santa Cruz, CA). Due to the low number of cones in wild type mouse retina (3% of photoreceptors), we also performed immunohistochemical analyses on frozen sections of adult pig retina, which contains Chromatin ImmunoprecipitationA commercially available assay kit (Upstate Biotechnologies, Charlottesville, VA) was used for ChIP studies. Briefly, four snap-frozen retinas from wild type mice were cross-linked for 15 min at 37 °C with 1% formaldehyde in PBS containing proteinase inhibitors. The retinas were washed four times in ice-cold PBS with proteinase inhibitors and then incubated on ice for 15 min. The tissue was then sonicated on ice eight times using 20-s pulses. The remaining steps were essentially performed as described by the manufacturer, using anti-NRL polyclonal antibody (12). Putative promoter regions (5' upstream of the transcription start site) for Rho, Bmp4, Smad4, and Bmpr1a were determined in silico (www.ncbi.nlm.nih.gov/mapview). Each DNA sequence was analyzed using MatInspector (www.genomatix.de/index.html). PCR primers were designed to flank the putative Nrl binding sites (Nrl response element, NRE) (18) either predicted manually or by MatInspector (Table I). The sequence closest to the transcription start site was chosen. ChIP DNAs from two independent experiments were used for PCR using equal amounts for input, with antibody, and no antibody reactions.
In Situ HybridizationThe 35S-labeled antisense and sense cRNA transcripts were synthesized from the mouse Bmp4 and Smad4 cDNAs cloned in pSPORT1 vector, according to the manufacturer's instructions (Ambion). The labeled transcripts were subjected to alkali hydrolysis (to give an average size of 70 nucleotides) and then hybridized overnight to 5- to 7-µm serial sections of adult mouse retina using standard protocols (Phylogeny Inc., Columbus, OH). After stringent washing, the slides were dipped in Kodak NTB-2 nuclear track emulsion and exposed for 12 weeks. Slides were developed in Kodak D-19, counterstained with toluidine blue, and analyzed by light and dark-field optics. Expression Analysis Using Purified Rod PhotoreceptorsTo obtain a highly enriched population of rod photoreceptors, we used Nrl-GFP transgenic mice that express enhanced green fluorescent protein (EGFP) under the control of Nrl promoter specifically in rods.4 Retinas of Nrl-GFP mice at P28 were dissected and dissociated using the Papain Dissociation System (Worthington). Cells that were positive and negative for EGFP were flow-sorted by FACSAria (BD Biosciences). Total RNA was used for RT-PCR using gene-specific primers, as follows. Bmpr1a: forward, AAGGAATGGGTGGGATTAGC; reverse, TGGCGATTGCCAACTAGATA; Grm6: forward, CAAGTAGCAAGGTTGAGTGT; reverse, GGAAGAATGCTGGAAGCAAG; Hprt: forward, CAAACTTTGCTTTCCCTGGT; reverse, CAAGGGCATATCCAACAACA; Nrl: forward, GCTGCATTTTCACCGAATCT; reverse, GGTGGTTTGGGTTGTGGTAG; Rho: forward, CTTCCTGATCTGCTGGCTTC; reverse, ACAGTCTCTGGCCAGGCTTA; Smad4: forward, ACCCGCGTATGCCGCCCCATCC; reverse, ACAGCGTCGCCAGGTGCTCGGC; Thy1: forward, AACTCTTGGCACCATGAACC; reverse, AGGCTGAACTCATGCTGGAT.
We chose to generate gene expression profiles of P21 retinas from rod-dominated wild type and rodless (cone-only) Nrl-/- mice. At this stage, the differentiation and laminar organization of retinal neurons are complete, and phototransduction pathways in the retina (from photoreceptors to ganglion cells) are fully functional. High Concordance between Direct and Indirect Microarray ComparisonsWe examined two different experimental designs, direct and indirect comparison (42), for their ability to identify differentially expressed genes. Two direct comparison experiments (two slides), in which Cy3-labeled wild type and Cy5-labeled Nrl-/- retinal RNA targets were hybridized simultaneously to the same slide, were performed. Five indirect comparisons were carried out with the reference RNA labeled by Cy3 and hybridized in conjunction with Cy5-labeled either wild type or Nrl-/- retinal RNA (total of 10 slides). Of the 13,440 spots on the I-gene microarray slides, 97.4% (13,092 spots) showed higher spot intensity than the background (i.e. were considered detected) in both direct comparison slides, whereas 91.6% (12,307 spots) spots were detected in all 10 slides with indirect comparisons. Scatter-plot analysis of signal intensities of the wild type and Nrl-/- retinas revealed high similarity between direct and indirect comparisons (Fig. 1). A majority of genes have log2 ratios centered at 0, indicating no change between the two tested samples. Only a few spots displayed over 4-fold change in both methods, although direct comparison showed a tighter scatter within the 4-fold lines. Both methods were able to successfully identify the duplicate spots of S-opsin and rhodopsin as the most up- or down-regulated genes in the Nrl-/- retina, respectively.
To compare the power of direct versus indirect methods in identifying differential gene expression, similar statistical analysis protocols were applied to both datasets to rank the genes based on their odds of being altered. The top 40 genes identified by direct comparison included 31 genes that were also ranked in top 40 by indirect comparison. Furthermore, 37 of the top 40 genes revealed by indirect experiments were present in the top 75 genes with direct comparison. This indicates a possible false-positive rate of 3/40 in the indirect comparison, assuming that only genes ranked in the top 75 by both methods are truly positive. Differential Expression of Genes in the Nrl-/- RetinaBased on the greatest odds of difference identified by both direct and indirect comparisons, a non-redundant set of 74 highly ranked genes was selected; of these, 50 genes are known (Table II), whereas 24 cDNAs only show homology to sequences in the EST or genomic databases (Table III). Of the 50 known differentially expressed genes, only 8 demonstrate higher expression in the Nrl-/- retina as compared with the wild type. The majority of these genes can be divided into the following functional groups: phototransduction, transcriptional regulation, signaling pathways, and structural or membrane-associated proteins. One-third of the altered genes (24/74) represent proteins of unknown function (Table III); such genes can be used as candidates for mutation screening of patients with rod and cone dystrophies or for biological studies.
We first wanted to determine whether our microarray profiling could validate changes in the expression of 16 phototransduction genes that had been analyzed previously by Northern analysis of P10 Nrl-/- retina (see Fig. 5A of Ref. 9). Our microarray data corroborated the expression changes of all genes (data not shown), except cone arrestin 3 (Arr3) and Pde6b, which did not exhibit any alteration. To further establish the fidelity of the microarray data, we examined the overlaid microarray images of spots corresponding to three genes: Gnb3 (cone transducin subunit) and its counterpart in rods, Gnb1, and one random EST (AC021049
[GenBank]
) (corresponding to the clone MRA-1648, with partial sequence homology to cold-shock domain protein A). The overlaid images of these three genes were in agreement with their microarray data (Fig. 2A). Northern blot and qRT-PCR analysis further confirmed the expression differences detected by microarray experiments (Fig. 2, B and C). Consistent with these observations, Gnb3 immunoreactivity was specific to cones: intense within outer segments, and moderate in cone cell bodies and synaptic pedicles (Fig. 2D, panels a and b). The staining within the outer segments was distinct from that observed for rod opsin (Fig. 2D, panel c), but co-localized precisely with PNA (Fig. 2D, panels d and e), a commonly used marker of cone outer matrix sheaths, demonstrating the cone-specific distribution of Gnb3 immunolabeling.
QRT-PCR Validation of Differentially Expressed GenesDue to the variability inherent in microarray technology, qRT-PCR analysis was applied to validate expression data (Tables II and III). For 33 out of 35 genes examined (95%), qRT-PCR data were in good agreement with the microarray results with respect to the direction of observed changes (i.e. qualitative). For 21 genes, microarray data underestimated the expression change. Two genes, Rxrg and Mtap6 (also called Stop), showed a discrepancy between the two methods. Microarray data indicated slight reduction, whereas qRT-PCR demonstrated clear augmentation of their expression in the mutant retina. This discrepancy may be the result of differential expression of various isoforms (generated by alternative splicing) encoded by the two genes. Although microarray method is hybridization-based and may reflect cumulative effects, qRT-PCR using gene-specific primers generates expression data for specific isoform(s). Repression of Bmp/Smad Signaling Pathway in the Nrl-/- RetinaMicroarray analysis showed that the expression of four genes belonging to the Bmp/Smad pathway (Bmp2, Bmp4, Bmpr1a, and Tob1) was reduced in the Nrl-/- retina. Because the expression and role of Bmp/Smad proteins have not been documented in mature retina, we decided to explore this further by qRT-PCR analysis of other pathway components that were not represented on I-gene microarray. All genes of the Smad-mediated Bmp signaling pathway exhibited lower expression in the cone-only Nrl-/- retina (Fig. 3A). The expression levels of Bmp4, Bmp2, and Bmpr1a were reduced by 12-, 4-, and 32-fold, respectively, in the Nrl-/- retina, whereas R-Smads (Smad 1 and Smad 5) and co-Smad (Smad 4) transcripts were marginally decreased. The expression of anti-Smads (Smad 6 and Smad 7) and Tob1 was also reduced in the Nrl-/- retina.
Altered Expression of the Components of Wnt/Ca2+ Signaling Pathway in the Nrl-/- RetinaAltered expression of several Ca2+/calmodulin-regulated genes, including Myo5a, Myo7a, and Mtap6, prompted us to examine the transcript levels of calmodulin (Calm2) and Camk2b by qRT-PCR. Both genes demonstrated altered expression in the Nrl-/- retinal RNA as compared with the wild type (Fig. 3B). We also detected higher expression of the Wnt5a ligand by qRT-PCR, indicating the overall increased activity and utilization of this pathway in the Nrl-/- retina. ChIP AnalysisTo investigate if the reduced expression of Bmp/Smad genes in the Nrl-/- retina is a direct result of the lack of Nrl, we performed ChIP assays that permit in vivo analysis of the binding of a transcription factor to its cognate cis-regulatory sequence elements (43, 44). DNA fragments bound to Nrl were isolated from mouse retina using an anti-Nrl antibody (12). The ChIP DNA was examined for the enrichment of putative promoter regions of the three Bmp/Smad pathway genes (Bmp4, Smad4, and Bmpr1a) (Fig. 4). PCR amplification of the NRE region in the rhodopsin promoter (14, 18) was used as positive control for ChIP assay. We were able to amplify the genomic DNA fragments, including the putative Nrl binding sites upstream of transcription start sites of Bmp4 and Smad4; however, no enrichment of the promoter region of Bmpr1a was observed.
Expression of Bmp4 and Smad4 in Mature Rod PhotoreceptorsTo further validate the possible role of the Bmp/Smad pathway and its regulation by Nrl, we decided to determine the expression of Bmp4 and Smad4 in mature rod photoreceptors. In situ hybridization to RNA in adult mouse retina sections revealed a positive signal for Bmp4 and Smad4 in all three nuclear layers, including the outer nuclear layer (ONL) of photoreceptors that primarily contains rods (Fig. 5A). RT-PCR analysis of the FACS-purified rod photoreceptors (GFP-positive and rhodopsin-positive) from the Nrl-GFP retina also demonstrated high enrichment of Smad4 transcripts in GFP(+) versus GFP(-) cells (Fig. 5B). GFP-positive cells also expressed Bmpr1a; however, the data for Bmp4 was inconsistent with the primers used for RT-PCR.
We then examined whether Bmp4 or Smad4 protein is detectable in mature photoreceptors. Various commercially available Bmp4 antibodies gave either high background or did not work well on retinal sections; hence, we were unable to observe specific Bmp4 immunostaining. Nevertheless, one Smad4 polyclonal antibody displayed a distinct pattern in accordance with its diminution in the Nrl-/- mouse retina (Fig. 5C). Within the adult wild type mouse retina, Smad4 immunostaining was observed throughout the inner segments (IS) and reaction product percolated down through the ONL as discrete fiber-like profiles (Fig. 5C, panels b and e). Immunostaining was more intense toward the scleral surface of the ONL but was far more widespread than the few cone cell bodies located in this position. Smad4 immunolabeling partially overlapped with that of rhodopsin (panel c) and did not correspond to PNA staining (panel f). Taken together, these results suggest Smad4 is within a subset of rod cells.
Our knowledge of signaling molecules and cell-cell communication pathways that specify physiological divergence between rods and cones and their homeostasis is currently limited. The availability of the Nrl-/- mouse, which exhibits a complete lack of rods and significantly enhanced cone function, permits dissection of molecular differences between the two types of photoreceptors and associated pathways. We took advantage of the advances in microarray technology and utilized custom I-gene microarrays to compare the gene expression profiles of mature P21 retinas of the rod-dominant wild type mouse with that of the Nrl-/- mouse. In addition to a large number of differentially expressed genes encoding known proteins, we identified 24 novel ESTs that might serve as candidate genes for human retinal diseases but would require additional molecular and functional analyses. The major findings from expression profiling presented in this report are that: (i) components of Bmp/Smad and Wnt/Ca2+ signaling pathways are expressed in the mature retina, suggesting their physiological relevance in regulating cell-cell communication and/or homeostasis; (ii) Bmp/Smad and Wnt/Ca2+ pathways are altered in the cone-only Nrl-/- retina, indicating their differential utilization in rod- and cone-associated inter- and intra-cellular signal transduction; (iii) at least two Bmp/Smad pathway genes (Bmp4 and Smad4) may be direct targets of regulation by Nrl; and (iv) Bmp4 and Smad4 are expressed in the rod photoreceptors of mature retina.
Global gene expression profiling using microarrays is a powerful methodology to generate molecular signatures of specific cell types, yielding valuable data for investigating complex biological processes. However, a careful experimental design is critical; for cDNA microarrays, issues related to direct versus indirect comparisons have been repeatedly emphasized (42, 45). Indirect comparison is more advantageous, because it allows researchers to cross-compare studies performed under different conditions. Recent advances have greatly minimized systematic sources of variations by improved technique (37), appropriate normalization (4648), and sufficient replication (49, 50). In this study, we used I-gene microarrays and evaluated the performance of indirect method using a reference RNA target relative to the direct comparison of wild type and Nrl-/- samples. Using the reference, we were able to obtain signals for >97% of spots on a single slide and about 91% of spots on all 12 slides studied. This high representation of retinal genes in the reference sample provides a basis for future indirect comparisons using I-gene microarrays and meta-analysis of the data. The percentage of positive spots on I-gene microarrays is substantially higher than that obtained with commercial microarray slides (52%) (31) or Affymetrix GeneChips ( Gene profiling of whole retina from wild type and Nrl-/- mice would yield expression differences that reflect the morphological and physiological changes associated with the loss of Nrl function. As predicted, the lack of rods and their replacement by cones resulted in higher expression of cone phototransduction genes with a concomitant reduction in rod gene transcripts. A majority of the differentially expressed genes showed lower expression in the Nrl-/- retina, reflecting the bias toward rod genes in the I-gene microarrays (containing cDNAs derived from rod-dominant wild type retinal tissue). Gene profiling suggested a remodeling of the Nrl-/- retina, as evident by altered expression of many genes encoding structural proteins (e.g. Prph2, Myo5a, Myo7a, Krt118, Cct4, and Mtap6) (see Table II). This is consistent with the collapse of the sub-retinal space, shortening of photoreceptor outer segments, and rewiring of neurons in the outer plexiform layer (9, 51). It should be noted that many of these expression changes probably reflect secondary rather than primary effects of the loss of Nrl.
Proteins of Wnt (wingless) and transforming growth factor-
Bone morphogenetic proteins (Bmps) are members of transforming growth factor- superfamily of secreted proteins that primarily transduce signals through Smad proteins (5356). Bmp signaling is essential for lens development and is involved in retinal differentiation (6770). Although a role for Bmps has been suggested in adult retinal pigment epithelium (71), little is known about Bmp and Smad pathways in the mature retinal photoreceptors. We postulate that reduced expression of genes of the Bmp/Smad pathway has relevance to the lack of rod function in the Nrl-/- retina. ChIP and expression data strongly suggest that Bmp4 and Smad4 are expressed in rod photoreceptors of adult mouse retina and that their expression is directly regulated by Nrl (see Figs. 4 and 5). The absence of enrichment of Bmpr1a (although it is expressed in rods) by ChIP suggests that Nrl may not directly act on its promoter (at least in the region tested). Modulation of Bmp/Smad signaling by Nrl argues in favor of its role in rod-mediated cell-cell communication. Activated Smad complexes can positively or negatively regulate transcription of numerous target genes and modify cellular response. We hypothesize that in addition to directly controlling the expression of specific phototransduction genes in mature rods (1518), Nrl influences rod function and homeostasis by integrating discrete signaling events.
* This research was supported in part by National Institutes of Health Grants EY11115 (including administrative supplements), EY07003, and GM72007, by The Foundation Fighting Blindness (Owings Mills, MD), by Research to Prevent Blindness (RPB, New York, NY), by the Macula Vision Research Foundation (West Conshohocken, PA), and by the British Retinitis Pigmentosa Society. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
c A recipient of Foundation Fighting Blindness-Canada Post-Doctoral Fellowship.
f Current address: University of Ottawa Eye Institute, Ottawa Health Research Institute, Ottawa, Ontario K1H 8L6, Canada. i To whom correspondence should be addressed: W. K. Kellogg Eye Center, University of Michigan, 1000 Wall St., Ann Arbor, MI 48105. Tel.: 734-615-2246; Fax: 734-647-0228; E-mail: swaroop{at}umich.edu.
1 The abbreviations used are: Nrl, neural retina leucine zipper; ChIP, chromatin immunoprecipitation; NRE, Nrl response element; Bmp, bone morphogenetic protein; EST, expressed sequence tag; P, postnatal day; E, embryonic day; qRT-PCR, quantitative reverse transcription-PCR; Hprt, hypoxanthine guanine phosphoribosyl transferase; PBS, phosphate-buffered saline; PNA, peanut agglutinin; GFP, green fluorescence protein; EGFP, enhanced GFP; ONL, outer nuclear layer; FACS, fluorescence-activated cell sorting; IS, inner segment.
2 A. J. Mears and A. Swaroop, unpublished data.
3 S. S. Nikonov, L. Daniele, C. Lillo, A. J. Mears, A. Swaroop, D. Williams, and E. N. Pugh, Jr., submitted for publication.
4 M. Akimoto and A. Swaroop, unpublished data.
We thank R. Farjo, M. I. Othman, and S. P. MacNee for generation of microarrays, Alfred Hero for discussions on data analysis, P. Hitchcock, T. Glaser, and A. J. Hunt for constructive suggestions, Jérôme Mutterer for confocal microscopy, Jennifer Masterson for technical, and Sharyn Ferrara for administrative assistance.
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||