|
|
||||||||
1From the Department of Ophthalmology, Vitreoretinal Research Lab, University of California at Davis, School of Medicine, Davis, California; and the 2Beckman Vision Center, University of California, San Francisco School of Medicine, San Francisco, California.
| Abstract |
|---|
|
|
|---|
METHODS. C57BL/6J, A/J, or B.A-Chr6 mice were exposed to 75% oxygen (hyperoxia) or room air for 14 days. After death, one eye was fixed and processed for outer nuclear layer (ONL) thickness measurements. The retina and RPE/choroid were separately dissected from the fellow eye and processed for microarray analysis. Single nucleotide polymorphism (SNP) analysis for transcribed sequences from the C57BL/6J and A/J genomes was conducted using the NIH genome site.
RESULTS. C57BL/6J mice developed a significant retinal degeneration in the inferior hemisphere after 14 days of hyperoxia. Under identical conditions, A/J mice exhibited only minor changes. A significant genetic effect was located on chromosome 6. SNP analysis of known transcribed sequences on chromosome 6 combined with microarray expression analysis yielded 33 candidate genes.
CONCLUSIONS. A significant genetic effect of susceptibility to HRRD is located on chromosome 6. In silico analysis of transcribed sequences results in a fairly small number of candidate genes.
Although the complex retinal degenerations that are currently being studied in the mouse2 3 may not precisely mimic human diseases, genetic studies in the mouse are made easier by exact environmental controls, the availability of many inbred strains with various physiological properties, and specialized recombinant inbred and chromosome substitution strains.4 Recently developed informatics tools also take advantage of the published genome sequences of many of these strains to identify single nucleotide polymorphisms (SNPs), which may be the source of allelic variation leading to variable phenotypes. The rationale for studying complex retinal degenerations in the mouse is the expectation that at least some of the genes to be identified have human orthologues, which are also important in complex human retinal degenerations.
Yamada et al.5 6 first described hyperoxia-related retinal degeneration (HRRD). This degeneration was originally observed in the C57BL/6J mouse after a 14-day exposure to 75% oxygen (hyperoxia).5 At this time point, animals exhibited a significant thinning of the outer nuclear layer (ONL) in the posterior retina of the eye. Recently, Walsh et al.7 have identified BALB/cJ as a mouse strain that is resistant to HRRD when compared with C57BL/6J.
The observation of a strain difference in susceptibility to HRRD suggests that this retinal degeneration can be studied by quantitative genetic approaches to identifying underlying genes. Our long-term goal is to identify one or more genes with alleles that are responsible for susceptibility to HRRD. We report the initial measurement of susceptibility to HRRD as a quantitative trait. We then present a new strain difference between C57BL/6J and A/J and the mapping of a significant genetic effect to chromosome 6 by the use of B.A chromosome substitution strains (CSSs). Finally, sequence and microarray expression analyses comparing transcribed sequences on chromosome 6 in the C57BL/6J and A/J mouse that yield an initial set of candidate genes are presented.
| Materials and Methods |
|---|
|
|
|---|
Hyperoxia Exposure
Mice were exposed to hyperoxia (75% oxygen) or room air (21% oxygen) for up to 14 days. The hyperoxia chamber was a standard 7-L induction chamber used for small animal anesthesia (VetEquip Inc., Pleasanton, CA). A mixture of medical grade oxygen and air was bubbled through water to humidify the air mixture entering the hyperoxia chamber. The hyperoxia chamber had 17 volume changes of exposure gas per hour. In addition, the percentage of oxygen exiting the chamber was continuously monitored (MiniOX I oxygen analyzer; MSA Instrument Division, Pittsburgh, PA). The chamber was opened daily to assess animal health and as necessary to feed and water the mice and clean the cages.
Tissue Processing and Histology
The mice were killed with CO2. Eyes were fixed and processed according to published procedures.8 9 One eye from each mouse was fixed for 24 hours in 2% paraformaldehyde/2.5% glutaraldehyde in phosphate buffer. The eye was then bisected from the superior to inferior region through the optic disc, the lens was removed, and the tissue was embedded in a PolyBed 812/Araldite mixture. Single sections (1 µm) cut along the vertical meridian of the eye were then stained with 1% toluidine blue containing 1% sodium borate.
The remaining eye was processed for RNA isolation as described previously with minor modifications.10 The eye cup was immobilized to the dissecting Petri dish by a drop of Superglue (QuickTite; Loctite Corp. for Manco, Inc., Avon, OH), which allowed rapid dissection of the eye. Retinal tissue and the RPE/choroid were taken out sequentially, immediately immersed into 300 µL of RNA stabilizer (RNALater; Ambion, Austin, TX), and stored at 20°C.
Quantification of ONL Thickness
The thickness of the ONL was quantified by an image processing method with a microscope and 40x objective (Nikon, Tokyo, Japan), in conjunction with image-analysis software (ImageProPlus; Media Cybernetics, Silver Spring, MD). ONL thickness was measured in nine superior and nine inferior fields, with 3 measurements taken in each field, for a total of 27 measurements in each of the superior and inferior hemispheres of the eye.8 9 One section through the optic disc in a superior-to-inferior orientation was used for quantification of each eye.
RNA Isolation
Samples for microarray analysis were collected from three animals from each of the treatments (normoxia and hyperoxia), from two eye tissues (retina and RPE/choroid), and from two strains (C57BL/6J and B.A-Chr6). Before RNA isolation, tissue was homogenized with a 22-gauge needle and a spin column (QIAshredder; Qiagen, Valencia, CA). Total RNA was isolated (RNeasy kit; Qiagen) according to a previously published protocol.10 The purity of total RNA was initially determined from the A260/280 ratio with a spectrophotometer (Beckman Coulter Inc., Fullerton, CA). Total RNA quality and concentration were determined using a microfluidics-chipbased assay that can detect nanogram quantities of RNA (RNA 6000 Nano Assay LabChip; Agilent Technologies, Inc., Palo Alto, CA) and run on microfluidics-based platform instrument (2100 BioAnalyzer; Agilent Technologies, Inc.). A typical yield from processing a single retina was 4 µg and 0.4 µg from a single RPE/choroid.
Microarray Procedures and Data Analysis
Labeled RNA samples were prepared according to the manufacturers manual (Affymetrix, Santa Clara, CA).11 A labeled cRNA probe was prepared from the 16 highest-quality total RNA samples (two replicates per treatment/tissue/strain), and run on a total of 16 mouse genome microarrays (Mouse Genome 430A_v2.0 GeneChips; Affymetrix). We followed the manufacturers protocol for small sample which required only 0.100 µg of RNA/sample (GeneChip Eukaryotic Small Sample Target Labeling Protocol; Affymetrix). Hybridization, washing, and staining were performed on the fluidics station (GeneChip Fluidics Station 400; Affymetrix) according to the manufacturers standard protocols. Hybridized microarrays were then scanned using a microarray scanner (Agilent GeneArray Scanner; Agilent Technologies, Inc.).
The acquisition, processing, and basic analysis of microarray-generated data were performed using the manufacturers gene expression analysis software (Affymetrix Gene Chip Operating Software [GCOS] and Gene Chip RNA Expression Analysis Software [GREX]), and raw data were normalized based on 100 housekeeping genes present on the microarray, and scaled to an average median signal value of 500. The GCOS algorithm was used to evaluate the abundance of each transcript represented on the array and labeled it as present (P), absent (A), or marginal (M). Comparison of normoxic and hyperoxic data was performed with GREX software (Affymetrix) generating signal log ratios (SigLogRatios) and significance estimates. SigLogRatios for biological replicates were averaged and SDs were calculated. Generally, SigLogRatio values greater or equal to 0.5 were considered to show upregulation, values smaller or equal to 0.5 downregulated, and values in between 0.5 and 0.5 as no change (NC), although GREX data analysis algorithm would make the final call taking into account the given probes signal/noise ratio.
SNP Analysis
Single nucleotide polymorphism (SNP) analysis in the expressed sequences was performed using the NCBI SNP database, Build 125 (www.ncbi.nlm.nih.gov/SNP/MouseSNP/ provided by the National Center for Biotechnology Information, Bethesda, MD). We compared C57BL/6J and A/J strains, limiting retrieval to chromosome 6. Only those genes affected by SNPs and expressed (as confirmed by microarray analysis) in the RPE/choroid or retina were retained in the analysis.
| Results |
|---|
|
|
|---|
|
Figures 2A and 2B show photomicrographs of the C57BL/6J retina at 14 days after hyperoxia exposure as well as normoxic controls. These photomicrographs focus on the inferior hemisphere of the eye and illustrate a significant loss of ONL thickness in this region for hyperoxia versus control exposed eyes. At 10 days of hyperoxia, little or no thinning of the ONL was apparent in the inferior hemisphere (data not shown).
|
Figure 2 shows a comparison of the C57BL/6J (Fig. 2A 2B) and A/J (Figs. 2C 2D) strains after 14 days of hyperoxia versus control mice. It is clear from these micrographs that the A/J strain is much less susceptible to HRRD at 14 days when compared with the C57BL/6J strain.
Quantification of ONL Thickness in the C57BL/6J and A/J Strains after HRRD
To put the observed degeneration of the inferior hemisphere of the retina on a quantitative basis, we used the standard procedures first described by LaVail et al.8 9 Figure 3 illustrates ONL thickness measurements for C57BL/6J (Fig. 3A) and A/J (Fig. 3B) mice under hyperoxic and normoxic conditions. This figure clearly demonstrates a region in the inferior hemisphere (enclosed by the box in both figures) in which a large degeneration occurred in the C57BL/6J mouse relative to the A/J mouse. Note that the data represent six individuals for each point, and that the variance of the data for these experiments gives the experiments sufficient power for this determination.
|
We initially planned to survey the entire panel of B.A CSSs. To maximize the speed with which we might find a significant QTL, we prioritized CSSs by the chromosomal assignment of QTLs for oxidative stressrelated phenotypes in other systems. Chromosome 2 was chosen first, as it contains the gene for Nfe2l2 (the mouse homolog of NRF2), a quantitative trait gene for other hyperoxia-related traits.12 An initial evaluation of C57BL/6J versus B.A-Chr2 indicated no significant difference between these strains with respect to HRRD susceptibility (data not shown). Our second choice was chromosome 6. Several QTLs for oxidative stressinduced phenotypes related to ozone and nickel treatment, as well as several QTLs being investigated in photoreceptor degeneration in the mouse are located on chromosome 6.2 3 13 14
Our evaluation of B.A-Chr6 was striking (Figs. 2E 2F) . This CSSs phenotype appeared to be much more like the A/J strain than the C57BL/6J. There was little to no apparent HRRD after 14 days. Figure 2 shows photomicrographs for a standard vertical section through the optic disc for the C57BL/6J, A/J, and B.A-Chr6 strains. A large degeneration is apparent for the C57BL/6J strain but not the B.A-Chr6 CSS. Quantification of these results is given in Figure 3 .
Finally, in Figure 4 , the quantitative HRRD data are presented in a normalized fashion for comparison among strains and F1 animals from reciprocal C57BL/6J X B.A-Chr6 crosses.
|
SNPs on Chromosome 6 between the C57BL/6J and A/J Strains for Transcribed Sequences
We analyzed the C57BL/6J and A/J sequences using the NIH site Entrez SNP for RefSNPs on chromosome 6 (http://www.ncbi.nlm.nih.gov/SNP/MouseSNP.cgi/ Provided by the National Center for Biotechnology Information, Bethesda, MD) for the transcribed sequences for all known genes and found a total of 911 in db SNP Mouse Build 125 (Table 1) .
|
Expression Analysis for Genes Containing SNPs
Individual genes can have more than one SNP. These SNPs can be grouped together within identified genes. Not all identified genes will be expressed in any given tissue, however. To discover the total number of genes with SNPs that are actually expressed, either in the retina or the RPE/choroid, microarray analysis was performed. The results of these two analyses are presented in Table 2 . This table gives information only on genes that were found to be expressed and for which in silico analysis indicated the presence of SNPs in the coding region.
|
|
|
| Discussion |
|---|
|
|
|---|
The use of hyperoxia in retinal research is not new. This type of exposure for newborn mouse pups produces what is now considered to be a standard model for retinal neovascularization.20 Yamada et al.5 first described the progressive degeneration of the posterior retina in the C57BL/6J mouse after 14 days of hyperoxia. Recently, Walsh et al.7 have also described HRRD in the C57BL/6J mouse and identified BALB/cJ as a resistant strain when compared with the C57BL/6J strain.
Yu et al.21 previously measured the intraretinal PO2 in the rat retina. These experiments demonstrated that the PO2 at the level of the RPE and photoreceptors was a linear function of inspired PO2. This is undoubtedly related to the high PO2 in the choroid and the volume of the choroidal circulation. The inner retina, however, displays significant autoregulation as does much of the brain, making the relationship between inspired PO2 and tissue PO2 a nonlinear function.
The mechanism by which hyperoxia stimulates oxidative stress has been discussed in the literature.22 Oxygen is presumed to diffuse into the cell and from the cytoplasm into the mitochondria. At higher than normal PO2 inside the mitochondria, various elements of the electron transport chain can directly reduce molecular oxygen to the superoxide anion. This species can then be dismutated into hydrogen peroxide and water among other reactions leading to reactive oxygen intermediates.
We submitted animals after 14 days of hyperoxia for routine pathology. We were surprised to learn that no other tissues in the body exhibited significant pathology. In particular, the lungs appeared normal. This tissue is clearly massively affected when mice are exposed to 95% oxygen for up to 72 hours.23
Our original intent was to survey the entire B.A CSS panel. We chose a specific order of initial strains hoping to find a significant QTL early in our survey. Chromosome 2 was selected first, as Cho et al.12 have demonstrated that a polymorphism in the NRF2 gene is responsible for different susceptibilities to pulmonary damage induced by hyperoxia. NRF2 is a gene that codes for a key transcription factor in the oxidative stress response. We found that B.A-Chr2 had the same phenotype as the C57BL/6J (sensitive). Chromosome 6 was our second choice for a variety of reasons. Previous QTLs on chromosome 6 have been identified for retinal degenerations in the mouse as well as nickel, ozone, and hyperoxia-induced pulmonary damage.3 14 24 25 Our findings on the phenotype of B.A-Chr6 indicate that a significant genetic effect in HRRD susceptibility is accounted for by genes on this chromosome.
The identification of a single chromosome which has one or more QTLs allows us to focus our attention on only 5% of the genome. Chromosome 6 in the mouse has 150 Mb containing 1822 identified genes according to the current NCBI Mus musculus chromosome 6 map Build 35.1 (NIH Map Viewer http://www.ncbi.nlm.nih.gov/mapview/). It is likely that this represents only 50% of the transcribed sequences, when undiscovered novel genes, new exons, antisense RNAs, and various types of noncoding but functional RNA species are considered.26
Our analysis was limited to transcribed sequences, and we therefore did not examine polymorphisms in promoter sequences. This is a serious limitation of our study, and a SNP between the C57BL/6J and A/J found in any nontranscribed part of the genome could be responsible for the difference in phenotype.
Given this note of caution, however, mouse geneticists now typically use in silico sequence analysis combined with expression analysis to shorten the list of candidate genes that may be studied in advance of the completion of fine mapping.27 We performed our SNP analysis using mouse data on the NIH site mentioned earlier.
Individual candidates can be selected from our list of candidate genes for which enough information exists to make direct follow-up studies appropriate. The gene for paraoxonase 2 (PON2), for example, is identified as a candidate that has a SNP in the 5' UTR region of the mRNA at the SNP flank position 208 (in the NT_039340.4). This gene is a member of the paraoxonases (PON) gene cluster, which includes Pon1, -2, and -3. Polymorphisms in the human PON1 gene have previously been associated with susceptibility to age-related macular degeneration (AMD).28 29 The product of PON1 is a protein bound to lipoprotein complexes in the serum; its function is to reduce oxidized lipid species. PON1 polymorphisms have been associated with other diseases such as atherosclerosis,30 and the common role of oxidative stress in all these phenotypes suggests that this may be a susceptibility gene with influence on a variety of diseases.
Our next studies will therefore concentrate on the genetic fine mapping of chromosome 6. Ultimately, genetic proof from knock-in or complementation testing by transgenesis with appropriate BAC clones will be necessary to prove that any set of SNPs or genes is responsible for phenotypic variation in HRRD.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication July 24, 2006; accepted November 27, 2006.
Disclosure: Z. Smit-McBride, None; S.L. Oltjen, None; M.M. LaVail, None; L.M. Hjelmeland, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Leonard M. Hjelmeland, Department of Ophthalmology, Vitreoretinal Research Lab, University of California, One Shields Ave., Davis CA 95616; lmhjelmeland{at}ucdavis.edu.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |