A Large Animal Model for CNGB1 Autosomal Recessive Retinitis Pigmentosa

Abstract

Retinal dystrophies in dogs are invaluable models of human disease. Progressive retinal atrophy (PRA) is the canine equivalent of retinitis pigmentosa (RP). Similar to RP, PRA is a genetically heterogenous condition. We investigated PRA in the Papillon breed of dog using homozygosity mapping and haplotype construction of single nucleotide polymorphisms within a small family group to identify potential positional candidate genes. Based on the phenotypic similarities between the PRA-affected Papillons, mouse models and human patients, CNGB1 was selected as the most promising positional candidate gene. CNGB1 was sequenced and a complex mutation consisting of the combination of a one basepair deletion and a 6 basepair insertion was identified in exon 26 (c.2387delA;2389_2390insAGCTAC) leading to a frameshift and premature stop codon. Immunohistochemistry (IHC) of pre-degenerate retinal sections from a young affected dog showed absence of labeling using a C-terminal CNGB1 antibody. Whereas an antibody directed against the N-terminus of the protein, which also recognizes the glutamic acid rich proteins arising from alternative splicing of the CNGB1 transcript (upstream of the premature stop codon), labeled rod outer segments. CNGB1 combines with CNGA1 to form the rod cyclic nucleotide gated channel and previous studies have shown the requirement of CNGB1 for normal targeting of CNGA1 to the rod outer segment. In keeping with these previous observations, IHC showed a lack of detectable CNGA1 protein in the rod outer segments of the affected dog. A population study did not identify the CNGB1 mutation in PRA-affected dogs in other breeds and documented that the CNGB1 mutation accounts for ∼70% of cases of Papillon PRA in our PRA-affected canine DNA bank. CNGB1 mutations are one cause of autosomal recessive RP making the CNGB1 mutant dog a valuable large animal model of the condition.

About the Authors:

Paige A. Winkler

Affiliations Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, Michigan, United States of America, Genetics Program, Michigan State University, East Lansing, Michigan, United States of America

Kari J. Ekenstedt

Affiliation: Department of Animal and Food Sciences, University of Wisconsin-River Falls, River Falls, Wisconsin, United States of America

Laurence M. Occelli

Affiliation: Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, Michigan, United States of America

Anton V. Frattaroli

Affiliation: Health Information Technology, Michigan State University, East Lansing, Michigan, United States of America

Joshua T. Bartoe

Affiliation: Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, Michigan, United States of America

Patrick J. Venta

Affiliations Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, Michigan, United States of America, Genetics Program, Michigan State University, East Lansing, Michigan, United States of America, Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, Michigan, United States of America

Simon M. Petersen-Jones

* E-mail: [email protected]

Affiliations Department of Small Animal Clinical Sciences, College of Veterinary Medicine, Michigan State University, East Lansing, Michigan, United States of America, Genetics Program, Michigan State University, East Lansing, Michigan, United States of America

Introduction

Retinitis pigmentosa (RP) is the leading cause of inherited blindness in humans affecting about 1 in 4,000 people [1]. It can be inherited in a dominant, recessive or X-linked fashion and shows considerable locus heterogeneity, with mutations in over 40 genes identified as causing non-syndromic RP (RetNet: https://sph.uth.edu/retnet/sum-dis.htm). Proteins encoded by these genes are necessary for a variety of functions within photoreceptors and their supporting cells. The age at onset and rate of progression of RP vary such that some patients have a history of night blindness from childhood while others may not notice symptoms until they are adults. The variability depends on the gene involved and the effect of the mutation on gene function, but there is also variability between patients with the same mutation [2], [3]. Rod photoreceptors are affected initially, resulting in loss of night (rod-mediated) vision and constriction of the visual fields. Loss of cone-mediated (daytime and color) vision may occur secondarily to rod-loss, even when RP is caused by a mutation of a gene exclusively expressed in rods, and can lead to complete blindness.

Retinal dystrophies analogous to RP occur in dogs, with reports of such conditions in over 100 different breeds [4]. The canine RP equivalent is known as progressive retinal atrophy (PRA) [5], [6]. The gene mutations underlying several forms of PRA have been identified and many have proven to be in genes analogous to those known to cause RP [7], [8], [9], [10], [11], [12], [13] or in some instances have suggested new candidate genes for investigation in RP patients [14], [15], [16].

Spontaneously occurring retinal dystrophies in canine models are of particular interest because the canine eye is similar in size to the human eye. This morphological similarity allows for identical surgical approaches for intravitreal and subretinal injection of therapeutic agents and testing for approaches such as implantation of intravitreal sustained-release devices. An additional advantage of canine models over rodent models is that the canine eye has regions of higher photoreceptor density (of both rods and cones), namely the area centralis and the visual streak that are somewhat analogous to the human macula [17]. In contrast, the retina of laboratory rodents lacks an equivalent region having an even density of photoreceptors across the retina [18]. Dogs with spontaneous mutations resulting in retinal dystrophies have proven to be important in preclinical assessment of therapies destined for use in human patients. For example, dogs with a mutation in RPE65 as a model for Leber congenital amaurosis type II were crucial for preclinical proof-of-concept gene therapy trials [19] which led to phase 1/2 human clinical trials [20], [21], [22]. The RPE65 mutant dog and other dog retinal dystrophy models have subsequently been used in several other preclinical trials for retinal gene and drug therapy [19], [23], [24], [25], [26], [27], [28], [29], [30], [31]. Identification of the gene mutations underlying other forms of canine PRA may provide additional spontaneous canine models to allow study of disease mechanisms and proof-of-concept therapy trials.

The Papillon breed of dog was initially reported to have PRA in 1995 [32]. Studies of the phenotype of affected dogs suggested loss of rod electrophysiological responses but maintenance of cone-driven responses at least until late in the disease process [33], [34]. Our unpublished studies of PRA in Papillons show a wide range in age of onset. This phenotypic variability could either suggest within-breed locus heterogeneity or could merely be the result of background genetic or environmental influences.

In this study we report a frameshift mutation in CNGB1 that is the cause of one form of PRA in Papillon dogs providing a large-animal model of autosomal recessive RP (RP45) due to CNGB1 mutations.

Results

Phenotypic Description

DNA samples were collected from 23 PRA-affected Papillons and 119 unaffected Papillons. The dogs had all been examined by a veterinary ophthalmologist. For the affected dogs an ophthalmoscopic diagnosis was made between 10 months and 13 years of age (data not shown). A small breeding colony of PRA-affected Papillons was established, consisting of an affected female and two affected offspring. Electroretinography (ERG) showed markedly reduced or absent rod-mediated ERG responses from an early age with preservation of cone photoreceptor responses (Figure 1). Spectral Domain-Optical Coherence Tomography (SD-OCT) performed on PRA-affected Papillons (confirmed to have the CNGB1 mutation described in this paper) showed that at the time the canine retina reaches maturity (approximately 8 weeks of age [35]) retinal layer thicknesses were comparable to a normal control (Figure 2A and 2B) and that affected dogs have a progressive thinning of the outer nuclear layer with age (Figure 2C & D).

[Figure omitted. See PDF.]

Figure 1. Representative ERG tracings from a normal control Papillon and a PRA-affected Papillon, both 10 weeks of age.

A. Dark-adapted ERG recordings at −2.4, −1.2 and 0.4 log cdS/m2. B. Light-adapted flash and flicker (33 Hz) ERG tracings. Background white light of 30 cd/m2 and flash intensity of 0.4 log cdS/m2. The vertical bars on the flicker ERG indicates the flash timing.

https://doi.org/10.1371/journal.pone.0072229.g001

[Figure omitted. See PDF.]

Figure 2. SD-OCT cross sectional images of the central retina from CNGB1 mutant Papillons and wild-type controls.

Note the progressive thinning of the outer nuclear layer of the retina in the affected animals. NFL/GCL – nerve fiber layer/ganglion cell layer; IPL – inner plexiform layer; INL – inner nuclear layer; OPL – outer plexiform layer; ONL – outer nuclear layer; PR – photoreceptor inner and outer segments; C – choroid. (Note the retinal pigment epithelium is not labeled).

https://doi.org/10.1371/journal.pone.0072229.g002

Mapping of the Papillon PRA Locus

An initial genome-wide association analysis performed using PLINK software [36] and including the genotypes from 23 Papillons (9 cases, 4 obligate carriers and 10 controls) yielded no significant associations (data not shown). Because within breed locus heterogeneity for PRA in dogs is a common occurrence, and because we and others had noted a wide range of age of onset [32] between the affected dogs, we suspected that more than one form of PRA may be segregating in the Papillon breed. Therefore, we analyzed the genotyping data from a small family group of 3 affected dogs within our breeding colony that we felt were very likely to share the same gene mutation (pedigree in Figure 3). We compared the genotyping results of the 3 affected dogs with that of 2 obligate carriers from the family and 11 control dogs (an additional control was included in this analysis). Our pedigree analysis supported an autosomal recessive mode of inheritance (data not shown) so we performed homozygosity mapping using a custom written computer program. The program was set to identify regions of homozygosity containing runs of at least six SNPs in the cases and for which the control animals did not share homozygosity (see Materials and Methods). This revealed 13 such regions of homozygosity greater than 1.5 Mb but only 4 of these regions contained obvious positional autosomal recessive RP candidate genes; CNGB1 (CFA2), RBP3 and RGR (CFA4), RD3 and CRB1(CFA7) and TULP1 (CFA12) (Table 1). Figure 3 shows a section of the run of homozygosity surrounding CNGB1 and the p-values for each marker, resulting from a chi-square association test corrected for multiple testing (the full region of homozygosity is shown in Table S4). These four regions were then subjected to haplotype construction, and haplotypes were examined within the small family group (data not shown). Only in the CFA2 region did the affected dogs from this family have a unique extended haplotype which was not present in the homozygous state in control (non-obligate carrier) dogs. Obligate carrier dogs each possessed one copy of this haplotype. Furthermore, after comparing the phenotype of the PRA-affected dogs in our colony with that reported for human families and mouse models with CNGB1 mutations, CNGB1 was considered the strongest candidate. Based on the haplotype analysis and phenotypic information, CNGB1 was selected to screen first for mutations.

[Figure omitted. See PDF.]

Figure 3. SNP analysis in the CNGB1 region using a small family group of affected and unaffected Papillons.

The SNPs are located in a 1.84 Mb region of homozygosity among the affected dogs (Only SNPs near the center of the region are shown). CNGB1 is located at chr2:61,454,476–61,520,336 in the reference genome (canFam2.0). SNP genotypes are given in columns for each dog, with the family group data displayed to the left of the black back and all others to the right of the black bar. All affected dogs (1–3) share the same haplotype, obligate carriers (C1–2) each have one copy of the haplotype seen in the affected dogs and none of the unaffected dogs (U1–11) possess the haplotype seen in the affected dogs. A-Major allele, B-minor allele, as designated by Illumina (Illumin Inc, San Diego, CA).

https://doi.org/10.1371/journal.pone.0072229.g003

[Figure omitted. See PDF.]

Table 1. Regions of homozygosity above 1.5 MB from small Papillon family.

https://doi.org/10.1371/journal.pone.0072229.t001

Genomic Structure of Canine CNGB1

To establish the genomic structure of canine CNGB1 we sequenced cDNA from a control canine retinal library and genomic DNA from a control Papillon and compared these to the published canine genomic sequence for this region (CanFam2.0). From these comparisons we deduced the intron/exon boundaries for CNGB1 (Table 2) which differed from the predicted structure seen on the University of California, Santa Cruz (UCSC) Genome Browser [37] (http://genome.ucsc.edu/).

[Figure omitted. See PDF.]

Table 2. Intron and exon boundries for canine CNGB1 gene.

https://doi.org/10.1371/journal.pone.0072229.t002

Sequence Analysis of CNGB1

The full coding region of canine CNGB1 was sequenced from the cDNA of a control canine retinal library (cDNA sequence submitted to GenBank KC527595). Sequencing cDNA from the control dog revealed three single nucleotide variants (SNVs), at locations c.3378C>A, c.3440T>C and c.3534G>A, and one previously reported single nucleotide polymorphism (SNP) at location c.151G>A (rs22870569). CNGB1 exons and exon/intron boundaries were sequenced from genomic DNA from control and affected Papillons. This revealed the same SNVs and SNP as seen in the sequenced cDNA but also contained two additional variants; a previously described SNP (c.27G>A, rs22870567) and a SNV allele that segregates with the affected phenotype (c.215C>T; p.P72L). PolyPhen-2 predicts that this is a benign change in amino acid (HumVar 0.412) [38].

The affected Papillon had a frameshift mutation in exon 26. This consisted of a 1 bp deletion (chr2: 61,502,597; c.2387delA) and a 6 bp insertion (between chr2:61,502,599–61,502,600; c. 2389_2390insAGCTAC). This mutation (c.2387delA;2389_2390insAGCTAC, which for simplicity we will refer to as CNGB1-fs26) is predicted to result in a premature stop codon, 17 bp downstream and is present in affected Papillons but not in unaffected Papillons or in the canine reference genome (UCSC Genome Browser CanFam2.0) (Figure 4).

[Figure omitted. See PDF.]

Figure 4. Papillon mutation in CNGB1 exon 26.

Sanger dideoxy-sequencing traces for part of CNGB1 exon 26 are shown for an unaffected (A) and an affected (B) Papillon. Panel C shows the codon and amino acid alignment inferred from the traces in panels A and B, for the unaffected sequence (top) and affected Papillon mutation sequence (bottom). The complex Papillon mutation includes a 6 bp AGCTAC insertion between reference bases chr2:61,502,599–61,502,600 (yellow highlight in panel C and within yellow box in panel B) and an adenine deletion at chr2: 61,502,597 (red highlight in panel C and red triangle in panel B). The deletion causes a frameshift and premature stop codon within seven residues, including the two new, inserted codons.

https://doi.org/10.1371/journal.pone.0072229.g004

Table 3 shows the SNVs and SNPs detected and Figure S1 shows the predicted canine protein amino acid sequence and alignment with other species. Primers for sequencing canine gDNA and cDNA are supplied in Table S1 and Table S2, respectively.

[Figure omitted. See PDF.]

Table 3. Single nucleotide variant (SNV) locations in genomic Papillon CNGB1.

https://doi.org/10.1371/journal.pone.0072229.t003

A genotyping assay for the mutation was developed (methods S1). The assay was used to genotype 20 Papillons that had been diagnosed by a veterinary ophthalmologist to have PRA (not including purpose-bred colony dogs) and 119 Papillons whose owners reported no abnormal vision (Table 4). The mutation was identified in the homozygous state in 13 of 20 Papillons that had been diagnosed with PRA. Of the phenotypically normal Papillons, none were homozygous for the mutation and 20 were heterozygous for it. This indicates a 16.5% carrier rate, and a mutated allele frequency of ∼31% in the Papillon breed, however, ascertainment bias almost certainly falsely inflates these values. In addition, we genotyped 33 dogs from 8 different breeds that had been diagnosed with PRA and 66 dogs from 9 different breeds that were clinically normal, none of the dogs of non-Papillon breeds had the mutation (Table S3).

[Figure omitted. See PDF.]

Table 4. PRA type 1 genotypes and clinical status for 139 Papillons.

https://doi.org/10.1371/journal.pone.0072229.t004

Immunohistochemistry Shows Lack of Detectable Full-length CNGB1 Protein in Affected Retina

To confirm that the CNGB1-fs26 mutation does disrupt CNGB1 expression in the homozygous animal, we performed immunohistochemistry (IHC) on retinal sections from an 8 week-old PRA-affected Papillon from the breeding colony that was homozygous for the CNGB1-fs26 mutation and compared it to retinal sections from an 8 week-old normal dog that was confirmed not to have the CNGB1-fs26 mutation.

CNGB1 in other species codes for multiple transcripts via alternative splicing [39], [40]. The CNGB1 locus has been described to code for four sensory transcripts; three retinal transcripts and one olfactory sensory transcript, as well as other splice variants expressed in kidney, brain, testes and spermatozoa [41], [42], [43], [44]. The 5′ portion of the gene encodes two glutamic acid rich proteins (GARPs) while the full-length transcript encodes the CNGB1 protein. The position of the CNGB1-fs26 mutation is predicted to allow normal expression of the two GARPs but to disrupt production of the full-length CNGB1 protein. To test this prediction we used two CNGB1 antibodies: one that targets the amino terminal (GARP region) of CNGB1 and the GARPs and a second antibody that targets the carboxyl end of CNGB1 downstream of both the GARP region and the predicted premature stop codon in the mutant canine CNGB1 gene. The results from this study showed that while the rod outer segments of the wild-type retina were labeled by both antibodies (Figure 5 A, C), the rod outer segments of the PRA-affected Papillon were labeled with the amino terminal antibody but not the carboxyl terminal antibody (Figure 5 B, D). This provides strong evidence that the mutation disrupts production of full-length CNGB1 protein while still allowing expression of GARPs as predicted from the CNGB1-fs26 mutation.

[Figure omitted. See PDF.]

Figure 5. Immunohistochemistry on frozen retinal sections from age (8 wk) and sex matched (female) control and affected dogs.

Images on the left are from the control dog and images on the right are from the affected dog. Panels A and B are stained with GARP-CNGB1 N-terminal antibody (green) and DAPI (blue). GARP proteins are present in both the control and affected samples. Panels C and D were stained with CNGB1 C-terminal antibody (green) and DAPI. CNGB1 full length protein is not detected in the affected sample. Panels E and F are stained with CNGA1 antibody (green) and DAPI. CNGA1 is not detected in the affected sample, presumably due to necessity for CNGB1 to form viable channels and normal trafficking. Size bar: 20 µm. OS- photoreceptor outer segment, IS – photoreceptor inner segment, ONL – outer nuclear layer, IPL – inner plexiform layer, INL – inner nuclear layer.

https://doi.org/10.1371/journal.pone.0072229.g005

We also performed IHC with a CNGA1 antibody. The rod cGMP-gated channel consists of both CNGA1 and CNGB1 subunits. Studies in mouse models have shown that lack of CNGB1 also disrupts trafficking of CNGA1 to the outer segments resulting in very reduced or absent CNGA1 protein [45], [46]. As predicted by the mouse model, the retinal sections for the affected Papillon showed lack of detectable CNGA1 protein while in the retinal sections from the normal dog CNGA1 was appropriately expressed and correctly targeted to the rod outer segments (Figure 5 E, F).

Discussion

We used homozygosity mapping of SNP microarray genotyping data from PRA-affected Papillons to identify regions where the affected dogs had runs of homozygosity spanning greater than 1.5 Mb and for which control dogs showed allelic variability. Positional candidate genes were identified mapping to those regions and haplotype analysis revealed only one region (on CFA2) in which obligate-carrier dogs had one copy of the haplotype observed in affected dogs, and none of the unrelated control dogs were homozygous for the haplotype observed in affected dogs. CNGB1 was selected as the most promising positional candidate gene because of its location within the CFA2 region and of the similarity of the phenotype of CNGB1 retinal dystrophy in mice and humans to that of the PRA-affected Papillons in our breeding colony (early lack of rod function and yet a slow retinal degeneration). The canine genome assembly (CanFam2.0) on the UCSC Genome Browser had incorrect exon predictions for CNGB1. We established the normal gene structure by sequencing cDNA from a retinal library developed from a control dog (GenBank KC527595). Sequencing of the confirmed exons and nearby flanking intronic regions of CNGB1, in affected and phenotypically normal Papillons, revealed a frameshift mutation (CNGB1-fs26) that segregated with disease status in our breeding colony and was not present in the homozygous state in any unaffected Papillons. A missense variation was also detected in the affected dogs (p.P72L). This converts a proline to leucine but at a residue that in some species is a leucine and in others a proline, making it unlikely that this would have a major effect on the function of CNGB1 in the dog [47]. Furthermore, this amino acid change was predicted to be benign by the PolyPhen-2 program [38].

Cyclic nucleotide gated (CNG) channels are necessary for normal phototransduction. There are different CNG channels encoded by paralogous genes specific to the rod and cone photoreceptors [41]. The rod CNG channel is a heterotetrameric protein consisting of three CNGA1 subunits and one CNGB1 subunit [48], [49], [50]. In dark conditions, the rod has higher levels of cyclic guanosine monophosphate (cGMP) which act to open a proportion of the CNG channels allowing an influx of cations, resulting in depolarization of the rod photoreceptor. Light stimulation triggers the phototransduction cascade activating the cGMP phosphodiesterase that hydrolyzes cGMP lowering its concentration and leading to closure of the CNG channels. This halts the influx of cations through the channel which, coupled with the continued action of other ion pumps in the cell membrane, results in the light-induced hyperpolarization of the rod photoreceptor.

CNGA1 can form functional channels in vitro without CNGB1 [51], [52], [53], while CNGB1 on its own does not form a functional channel. However, in vivo CNGB1 is important for rod CNG formation and normal functionality of the channel [45], [46]. The CNGB1 gene codes for three splice variants in the retina that encode the soluble glutamic acid rich proteins 1 and 2 (GARP1 and GARP2) and a full length CNGB1 protein [39]. The full length transcript in the retina is known as CNGB1a. Full length CNGB1 protein in the retina consists of an N-terminal “GARP region” and a C-terminal “channel domain”. GARP1 is of low abundance whereas GARP2 is more highly expressed [46]. The functions of the GARP subunits are still being explored. GARP2 does bind to PDE6 and may act to reduce dark level noise [54]. It is also postulated to have a structural role in the rod outer segment [46] where it interacts with peripherin-2 at the rod outer segment disk rim [55]. As discussed further below GARP2 may also play a role in control of the opening of the CNG channel [56]. Other than these three splice variants expressed in the retina, CNGB1 also has splice variants expressed in the olfactory epithelium (CNGB1b –which encodes a shorter protein than expressed in the retina lacking the GARP region), kidney, brain and testes [41], [42], [43], [44]. Two mouse gene targeting models, CNGB1X-1 and CNGB1X-26, have been created to study the functions of these proteins in the retina.

The CNGB1X-1 mouse knockout model, which lacks all retinal CNGB1 products, has shown that GARPs are necessary in outer segment disk development and the structural integrity of the rod outer segments [46]. Interestingly, the CNGB1X-1 mice have a weak rod response that is detectable on single cell recording. This activity originates from a low level of homomeric CNGA1 channel formation in the rod outer segments. Despite the presence of the channels and the resulting weak light-induced rod response, these mice have severe rod photoreceptor degeneration. It was suggested that degeneration may be due to rod structural, rather than functional, failure [46]. Another mouse model, CNGB1X-26, has an engineered mutation that leads to a deletion of exon 26 resulting in a frameshift that introduces a premature stop codon in the first triplet in exon 27 of CNGB1. The CNGB1X-26 mouse has a complete absence of mRNA for full-length CNGB1 and thus a lack of the CNGB1 protein but still produces GARP proteins [45]. These mice have normal development of the rod outer segments and a slower retinal degeneration than the CNGB1X-1 mouse.

The premature stop codon in the CNGB1X-26 mouse is positioned at a site homologous to 10 codons upstream of the premature stop codon predicted in the CNGB1-fs26 mutation identified in this study. While we anticipate that the premature stop codon in CNGB1 in the PRA-affected Papillon will lead to nonsense-mediated decay of the abnormal CNGB1 mRNA and a complete absence of the CNGB1 protein, as was reported for the CNGB1X-26 mouse [45], rather than production of a truncated protein, further studies are required to confirm this.

IHC using an antibody that binds to the carboxyl end of CNGB1 confirmed the lack of full-length CNGB1 in the rod outer segments of a young affected Papillon prior to photoreceptor degeneration. The CNGB1-fs26 mutation is not predicted to affect expression of the GARPs, and IHC using an antibody that recognizes GARPs as well as the GARP region of the full-length CNGB1 labeled the rod outer segments of the affected Papillon suggesting that GARPs and/or a truncated CNGB1 product are expressed. If the truncated mRNA avoids degradation and allows the production of a truncated protein, it would be missing domains essential for channel gating and cyclic nucleotide binding and would therefore not be expected to function normally. IHC also showed a lack of detectable CNGA1 protein in the rod outer segments in the young affected Papillon. Similar IHC results were observed in the CNGB1X-26 mouse. It was suggested that the GARPs may have a dominant-negative effect on the transport of CNGA1 to the outer segment accounting for the lack of homomeric CNGA1 channels in the outer segments of CNGB1X-26 mice (which have GARPs) contrasted with the presence of low levels of CNGA1 in outer segments of CNGB1X-1 mice (which lack GARPs) [46]. The reduced level of CNGA1 in the outer segment of the CNGB1X-1 mouse is probably a reflection of the importance of CNGB1 in transport of CNGA1 to the outer segment [46]. A region in the N-terminal section of the CNGB1 subunit (bovine amino acids #677-764) has been described to interact and promote protein-protein interactions with the CNGA1 subunit C-terminal region; this region is upstream of the CNGB1-fs26 mutation and is not predicted to be effected by this mutation [57], [58]. As would be expected in animals with a lack of rod CNG channels, the ERG changes in the affected Papillons indicate a lack of rod function. Specifically, prior to a thinning of the outer nuclear layer, there was an elevated dark-adapted ERG threshold with only a recordable response to light intensities above cone response threshold, a marked reduction in a- and b-wave amplitudes for the dark-adapted brighter flashes that elicit a mixed rod cone response in normal dogs.

The CNGB1X-26 mutation also lacks expression of the CNGB1 splice variant expressed in the olfactory epithelium (CNGB1b). These mice show absence of proper channel localization and delayed olfaction [59]. More detailed phenotypic characterization of the dog model is required to ascertain how closely the canine phenotype, both retinal and olfactory, mimics the CNGB1X-26 mouse model.

Mutations in CNGB1 have been identified in human patients with autosomal recessive RP (RP45) and are reported to account for ∼4% RP of cases [1]. In 2001, Bariel et al. reported a consanguineous French family with RP that they mapped to an interval containing CNGB1 [60]. They then identified a missense mutation that converted an evolutionarily conserved glycine to valine (c.2978G>T; p.G993V) which they predicted would alter the cyclic nucleotide-binding domain (CNBD). The effect of this mutation was further elucidated in an elegant study reported by Michalakis et al [56]; they showed that the p.G993V mutation prevented binding of cGMP and that binding is required for removal of an inhibitory effect that the GARP domain and also GARP2 has on channel opening, such that if cGMP cannot bind to CNGB1 the CNG channel is silent [56]. More recently resequencing of candidate RP genes led to the identification of a simplex RP patient homozygous for a missense substitution in CNGB1 (c.2957A>T; p.N986I) resulting in substitution of a conserved amino acid 7 codons upstream of the mutation identified by Bariel et al. and also in the CNBD [61]. Kondo et al. used homozygosity mapping to screen known arRP genes in Japanese RP patients and in one patient identified a mutation at the donor site of exon 32 (c.3444+1G>A) of CNGB1 [62]. Subsequently, Becirovic et al. performed exon trapping experiments to investigate the effect of the mutation. Their studies suggest that the mutation leads to skipping of exon 32 and replacement of the last 170 amino acids by 68 unrelated amino acids [63]. The probands in the Bariel et al. and Kondo et al. studies had night blindness from a young age and were diagnosed with RP in their 30′s [60], [62].

The CNGB1-fs26 mutation identified in Papillon dogs was present in 13 of the 20 PRA-affected Papillons tested. This suggests that there is at least one additional PRA locus segregating within the breed. Within-dog breeds genetic heterogeneity for PRA is becoming more evident [13], [16]. Additional studies will be required to find the gene mutation(s) responsible for the other form(s) of PRA segregating in Papillons.

The early onset of loss of rod function in the CNGB1-fs26 mutant dog, coupled with a slow retinal degeneration that we have observed in our colony dogs, seems to accurately parallel the described disease course in human patients as well as the comparable mouse model (CNGB1X-26). Recently recombinant adeno-associated viral vector-mediated gene therapy to deliver a normal copy of CNGB1 to CNGB1−/− mice was reported to allow for CNG channel formation, restoration of rod function and retinal morphological preservation [64]. The early loss of rod function and yet slow rod photoreceptor loss in the animal models suggests that CNGB1 RP is a good target for gene augmentation therapy. The CNGB1-fs26 mutant dog promises to be a valuable model for preclinical trials of such therapy.

Materials and Methods

Ethics Statement

All procedures were in compliance with the ARVO statement for the Use of Animals in Ophthalmic and Vision Research and approved by the Michigan State University Institutional Animal Care and Use Committee (AUF number 05-11-106-00; Institutional NIH/PHS Animal Welfare Assurance number A3955-01).

Electroretinography

To assess rod and cone photoreceptor function, electroretinograms (ERGs) were recorded using a modification of a previously described technique [26]. Briefly, ERGs were recorded using an Espion E2 Electrophysiology system with ColorDome Ganzfeld (Diagnosys LLC, Lowell, MA) and bandpass set between 0.5 and 500 Hz. Dogs were dark-adapted for one hour and anesthetized with injectable propofol (10 mg/kg PropoFlo, Abbott Animal Health, North Chicago, IL), intubated and maintained on inhaled 1 to 2% isoflurane (Isoflo, Abbott Laboratories, North Chicago, IL) delivered in oxygen. The ERG assessment consisted of three dark-adapted flash intensities at −2.4, −1.2 and 0.4 log cdS/m2 to record rod and mixed rod-cone responses. This was followed by light adaptation at 30 cd/m2 for 10 minutes and recording of a light-adapted response to a 0.4 log cdS/m2 flash and then 33 Hz flicker responses to the same intensity both superimposed on the same background light.

Spectral Domain-Optical Coherence Tomography

Assessment of retinal morphology was performed by Spectral Domain-Optical Coherence Tomography (SD-OCT; Spectralis OCT+HRA Heidelberg Engineering Inc., Heidelberg, Germany). Dogs were anesthetized as described for ERG, the pupil dilated with 1% topical tropicamide (Mydriacyl, Alcon Laboratories, Honolulu, HI, USA), a lid speculum fitted and the eye positioned in primary gaze using a stay suture in the inferior perilimbal conjunctiva. High-resolution cross-section images obtained by line and volume scanning and images from the same region of the central retina of affected and control (wild-type) Papillons were assessed.

Animal Use and Sample Collection

A pregnant female Papillon dog that had been diagnosed with PRA and had been mated with a PRA-affected Papillon stud dog was donated to the Michigan State University Comparative Ophthalmology laboratory with the consent of the owner to allow the study of the phenotype of PRA in the breed. This female and her offspring were used to establish a small breeding colony of dogs. The colony was kept under standard laboratory housing with 12:12 hr light:dark cycles.

Blood samples from client-owned Papillon dogs were donated with owner consent. DNA was extracted from blood samples using a commercial DNA extraction kit with a modified protocol (Qiagen Sciences, Germantown, MD). Briefly, a red blood cell lysis buffer (0.32 M sucrose, 10 mM Tris, 5 mM MgCl2) was added in a 2 step fashion to whole blood (3X volume and then 2X volume, respectively). Cell lysis solution (Qiagen Sciences, Germantown, MD) was added to lyse the white blood cells followed by addition of protein precipitation solution (Qiagen Sciences, Germantown, MD), isopropanol DNA precipitation and a 70% ethanol wash step.

The retina from a mixed breed dog was dissected from an enucleated eye and placed in an RNA stabilization buffer (RNAlater, Qiagen Sciences, Germantown, MD) and stored in a −80°C freezer until RNA extraction. RNA was extracted using an RNEasy kit according to manufacturer’s protocol (Qiagen Sciences, Germantown, MD). cDNA was made from mRNA using a 3′ RACE kit according to manufacturer’s protocol (Invitrogen, Carlsbad, CA).

Genome-Wide Association Mapping

Twenty-four Papillons (9 cases, 15 controls) were genotyped for 173,662 single nucleotide polymorphisms (SNPs) using Illumina Canine HD BeadChips. Initial genome-wide association analysis was conducted using the genome analysis toolset PLINK [36]. SNPs with a minor allele frequency (MAF) of <5% and with missing genotype calls of >10% were removed from the analysis. The final data set consisted of 116,235 markers. All 24 individuals genotyped successfully for over 90% of the SNPs and were retained in the analysis. One of the control dogs was removed from the PLINK analysis but was used in the run of homozygosity analysis. The final genotyping rate was >99.8%. Chi-square association mapping was conducted in PLINK, and correction for multiple testing was achieved using the Max(T) permutation procedure (10,000 permutations) in PLINK.

Custom Sorting Program and Haplotype Construction

Homozygosity mapping was performed using a custom sorting program. An algorithm was written to search for blocks of SNPs where there was a difference in the calls between the affected and unaffected dogs. The high quality SNP data, generated from PLINK, was imported to a Microsoft SQL Server 2008R2 database (Microsoft Corporation, Redmond, WA) to take advantage of its efficient set theory based querying mechanism. The received data was formatted with dog identifiers as columns and SNPs as rows. Indexing the data provided a method to measure continuity. The query used scalar-valued functions to assess the criteria described below and attached a flag to each row identified. For each individual SNP, if the affected dogs shared the same homozygous genotype, then the unaffected dogs were compared and rows identified where 95% of the unaffected dogs did not share the same genotype as the affected dogs. Upon those criteria being met, the affected and unaffected groups were considered ‘different’ by the algorithm. Identified SNPs were then sorted into groups formed by having a level of adjacency of at most four SNPs apart. Only sections with 6 or more SNPs meeting the above criteria were marked for further analysis. These regions were then sorted by size of the region and all regions 1.5 MB or larger were inspected for arRP candidate genes using the University of California, Santa Cruz (UCSC) Genome Browser [37] (http://genome.ucsc.edu/). Four regions of particular interest (on CFA2, CFA4, CFA7, and CFA12) were then subjected to haplotype construction, using fastPHASE [65]. Haplotypes were manually examined for shared regions in related affected family members.

DNA Sequencing

The UCSC Genome Browser (http://genome.ucsc.edu/) CanFam2.0 was used in conjunction with the cDNA sequences to identify the exons for the CNGB1 gene. Primers (Table S1) were designed flanking the entire exon and the splice sites using Primer3 (http://frodo.wi.mit.edu/). Sanger dideoxy-sequencing was done by an ABI 3730 Genetic Analyzer (Applied Biosystems, Inc., Foster City, CA) at Michigan State University’s Research Technology Support Facility.

CNGB1 Genotyping Assay

A restriction enzyme digest was designed to quickly screen dogs for the mutation in CNGB1 (see methods S1). This assay was used to test for the presence of the mutation in 139 Papillon dogs, 33 PRA affected dogs from 8 different breeds and 66 unaffected dogs from 9 different breeds.

Immunohistochemistry

A PRA-affected female Papillon from the research colony and an unaffected mixed breed female control dog were humanely euthanized at 8 weeks of age. The eyes were enucleated and the right eye was fixed in paraformaldehyde following a previously describe protocol [66].

Frozen sections were immunolabeled with either a rabbit anti-mouse CNGB1 N-terminal (kindly provided by Dr. Stephen Pittler), rabbit anti-human CNGB1 C-terminal (Sigma-Aldrich, St Louis, MO) or mouse monoclonal CNGA1 ([67] kindly provided by Dr. Bob Molday) antibody. Sections were blocked with 10% horse serum (Sigma-Aldrich, St Louis, MO) for 2 hours at room temperature and labeled with primary antibodies (dilutions of: N-terminal CNGB1 1:100, C-terminal CNGB1 1:300 and CNGA1 1:10) overnight at 4°C. Secondary antibodies (anti-rabbit or anti-mouse Alexa Fluor 488, 1∶500) (Invitrogen Molecular Probes, Carlsbad, CA) were placed on the sections for 2 hours at room temperature. All sections were counterstained with the nuclear stain DAPI (Invitrogen Molecular Probes, Carlsbad, CA).

Sections were imaged using a fluorescent microscope (Nikon Eclipse 80i, Nikon Instruments Inc., Melville NY) using commercial image capture software (MetaVue, Molecular Devices, Sunnyvale CA).

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

CNGB1 amino acid alignments. Sequence alignments performed using muscle alignment in SeaView software ([68]. Single nucleotide polymorphism and variants found in Papillon gDNA sequencing are numbered (1: p.E69K, 2: p.P72L, 3: p.L1165P) and the changed amino acid is underlined. 2: p.P72L is the SNV that has only been seen in Papillons with the mutation and is marked with red text. The CNGB1X-26 mouse stop codon is highlighted in red. The Papillon mutation is highlighted in green. The human mutations are marked in yellow (highlighted N [61] G [60]and arrow which represents a splice mutation [62]). Epitopes for N-terminal antibody (mouse, highlighted in purple) and C-terminal antibody (human, highlighted in teal).

https://doi.org/10.1371/journal.pone.0072229.s001

(DOCX)

Table S1.

Primers for genomic DNA sequencing.

https://doi.org/10.1371/journal.pone.0072229.s002

(XLS)

Table S2.

Primers for cDNA sequencing.

https://doi.org/10.1371/journal.pone.0072229.s003

(XLS)

Table S3.

Non-Papillon breeds tested for CNGB1 mutation.

https://doi.org/10.1371/journal.pone.0072229.s004

(XLS)

Table S4.

The entire region of homozygosity surrounding the site of the CNGB1 mutation.

https://doi.org/10.1371/journal.pone.0072229.s005

(XLS)

Methods S1.

https://doi.org/10.1371/journal.pone.0072229.s006

(DOCX)

Acknowledgments

The authors of this paper would like to thank the Papillon Club of America for their continued support throughout this project. Their effort in acquiring the affected female dog that started the Michigan State University research colony was integral to finding the mutation. We would also like to thank the breeders and owners who donated samples from their dogs. The authors thank Dr. Bob Molday for the CNGA1 antibody and Dr. Stephen Pittler for the N-terminal CNGB1 antibody. The authors would also like to thank Janice Querubin for her help with management of the Papillon research colony and recording the ERGs.

Author Contributions

Conceived and designed the experiments: PAW SMP-J PJV JTB KJE. Performed the experiments: PAW LMO SMP-J. Analyzed the data: KJE AVF PAW PJV SMP-J. Contributed reagents/materials/analysis tools: KJE AVF. Wrote the paper: PAW SMP-J.

Citation: Winkler PA, Ekenstedt KJ, Occelli LM, Frattaroli AV, Bartoe JT, Venta PJ, et al. (2013) A Large Animal Model for CNGB1 Autosomal Recessive Retinitis Pigmentosa. PLoS ONE 8(8): e72229. https://doi.org/10.1371/journal.pone.0072229

References

1. Hartong DT, Berson EL, Dryja TP (2006) Retinitis pigmentosa. Lancet 368: 1795–1809.

2. Maubaret CG, Vaclavik V, Mukhopadhyay R, Waseem NH, Churchill A, et al. (2011) Autosomal dominant retinitis pigmentosa with intrafamilial variability and incomplete penetrance in two families carrying mutations in PRPF8. Invest Ophthalmol Vis Sci 52: 9304–9309.

3. Berson EL, Rosner B, Sandberg MA, Weigel-DiFranco C, Dryja TP (1991) Ocular findings in patients with autosomal dominant retinitis pigmentosa and rhodopsin, proline-347-leucine. Am J Ophthalmol 111: 614–623.

4. Whitley RD (1988) Focusing on eye disorders among purebred dogs. Veterinary Medicine 83: 50–63.

5. Petersen-Jones SM (1998) Animal models of human retinal dystrophies. Eye 12 (Pt 3b): 566–570.

6. Aguirre GD, Acland GM (2006) Models, mutants, and man: searching for unique phenotypes and genes in the dog model of inherited retinal degeneration. In: Ostrander EA, Giger U, Lindblad-Toh K, editors. The dog and its genome. Cold Spring Harbor: Cold Spring harbor Laboratory Press. 291–326.

7. Suber ML, Pittler SJ, Quin N, Wright GC, Holcombe N, et al. (1993) Irish setter dogs affected with rod-cone dysplasia contain a nonsense mutation in the rod cGMP phosphodiesterase beta-subunit gene. Proc Natl Acad Sci U S A 90: 3968–3972.

8. Clements PJ, Gregory CY, Peterson-Jones SM, Sargan DR, Bhattacharya SS (1993) Confirmation of the rod cGMP phosphodiesterase beta subunit (PDE beta) nonsense mutation in affected rcd-1 Irish setters in the UK and development of a diagnostic test. Curr Eye Res 12: 861–866.

9. Petersen-Jones SM, Entz DD, Sargan DR (1999) cGMP phosphodiesterase-α mutation causes progressive retinal atrophy in the Cardigan Welsh corgi dog. Invest Ophthalmol Vis Sci 40: 1637–1644.

10. Dekomien G, Runte M, Godde R, Epplen JT (2000) Generalized progressive retinal atrophy of Sloughi dogs is due to an 8-bp insertion in exon 21 of the PDE6B gene. Cytogenet Cell Genet 90: 261–267.

11. Zhu L, Jang GF, Jastrzebska B, Filipek S, Pearce-Kelling SE, et al. (2004) A naturally occurring mutation of the opsin gene (T4R) in dogs affects glycosylation and stability of the G protein-coupled receptor. J Biol Chem 279: 53828–53839.

12. Kukekova AV, Goldstein O, Johnson JL, Richardson MA, Pearce-Kelling SE, et al. (2009) Canine RD3 mutation establishes rod-cone dysplasia type 2 (rcd2) as ortholog of human and murine rd3. Mamm Genome 20: 109–123.

13. Downs LM, Bell JS, Freeman J, Hartley C, Hayward LJ, et al. (2012) Late-onset progressive retinal atrophy in the Gordon and Irish Setter breeds is associated with a frameshift mutation in C2orf71. Anim Genet 44: 169–177.

14. Goldstein O, Kukekova AV, Aguirre GD, Acland GM (2010) Exonic SINE insertion in STK38L causes canine early retinal degeneration (erd). Genomics 96: 362–368.

15. Zangerl B, Goldstein O, Philp AR, Lindauer SJ, Pearce-Kelling SE, et al. (2006) Identical mutation in a novel retinal gene causes progressive rod-cone degeneration in dogs and retinitis pigmentosa in humans. Genomics 88: 551–563.

16. Downs LM, Wallin-Hakansson B, Boursnell M, Marklund S, Hedhammar A, et al. (2011) A frameshift mutation in golden retriever dogs with progressive retinal atrophy endorses SLC4A3 as a candidate gene for human retinal degenerations. PLoS One 6: e21452.

17. Mowat FM, Petersen-Jones SM, Williamson H, Williams DL, Luthert PJ, et al. (2008) Topographical characterization of cone photoreceptors and the area centralis of the canine retina. Mol Vis 14: 2518–2527.

18. Carter-Dawson LD, LaVail MM (1979) Rods and cones in the mouse retina. I. Structural analysis using light and electron microscopy. J Comp Neurol 188: 245–262.

19. Acland GM, Aguirre GD, Ray J, Zhang Q, Aleman TS, et al. (2001) Gene therapy restores vision in a canine model of childhood blindness. Nat genet 28: 92–95.

20. Bainbridge JW, Smith AJ, Barker SS, Robbie S, Henderson R, et al. (2008) Effect of gene therapy on visual function in Leber's congenital amaurosis. N Engl J Med 358: 2231–2239.

21. Maguire AM, Simonelli F, Pierce EA, Pugh EN Jr, Mingozzi F, et al. (2008) Safety and efficacy of gene transfer for Leber's congenital amaurosis. N Engl J Med 358: 2240–2248.

22. Cideciyan AV, Hauswirth WW, Aleman TS, Kaushal S, Schwartz SB, et al. (2009) Human RPE65 Gene Therapy for Leber Congenital Amaurosis: Persistence of Early Visual Improvements and Safety at 1 Year. Hum Gene Ther 20: 999–1004.

23. Acland GM, Aguirre GD, Bennett J, Aleman TS, Cideciyan AV, et al. (2005) Long-term restoration of rod and cone vision by single dose rAAV-mediated gene transfer to the retina in a canine model of childhood blindness. Mol Ther 12: 1072–1082.

24. Narfström K, Katz ML, Bragadottir R, Seeliger M, Boulanger A, et al. (2003) Functional and structural recovery of the retina after gene therapy in the RPE65 null mutation dog. Invest Ophthalmol Vis Sci 44: 1663–1672.

25. Jacobson SG, Acland GM, Aguirre GD, Aleman TS, Schwartz SB, et al. (2006) Safety of recombinant adeno-associated virus type 2-RPE65 vector delivered by ocular subretinal injection. Mol Ther 13: 1074–1084.

26. Annear MJ, Bartoe JT, Barker SE, Smith AJ, Curran PG, et al. (2011) Gene therapy in the second eye of RPE65-deficient dogs improves retinal function. Gene Ther 18: 53–61.

27. Gearhart PM, Gearhart C, Thompson DA, Petersen-Jones SM (2010) Improvement of visual performance with intravitreal administration of 9-cis-retinal in Rpe65-mutant dogs. Arch Ophthalmol 128: 1442–1448.

28. Beltran WA, Cideciyan AV, Lewin AS, Iwabe S, Khanna H, et al. (2012) Gene therapy rescues photoreceptor blindness in dogs and paves the way for treating human X-linked retinitis pigmentosa. Proc Natl Acad Sci U S A 109: 2132–2137.

29. Komaromy AM, Alexander JJ, Rowlan JS, Garcia MM, Chiodo VA, et al. (2010) Gene therapy rescues cone function in congenital achromatopsia. Hum Mol Genet 19: 2581–2593.

30. Mowat FM, Bartoe JT, Bruewer A, Dinculescu A, Boye SL, et al. (2012) Evaluation Of Rod Photoreceptor Function And Preservation Following Retinal Gene Therapy In The PDE6A Mutant Dog. ARVO Meeting Abstracts 53: 1928.

31. Petit L, Lheriteau E, Weber M, Le Meur G, Deschamps JY, et al. (2012) Restoration of vision in the pde6beta-deficient dog, a large animal model of rod-cone dystrophy. Mol Ther 20: 2019–2030.

32. Håkanson N, Narfström K (1995) Progressive retinal atrophy in papillon dogs in Sweden: A clinical survey. Vet Comp Ophthalmol 5: 83–87.

33. Narfström K, Ekesten B (1998) Electroretinographic evaluation of Papillons with and without hereditary retinal degeneration. Am J Vet Res 59: 221–226.

34. Narfström K, Wrigstad A (1999) Clinical, electrophysiological and morphological changes in a case of hereditary retinal degeneration in the Papillon dog. Vet Ophthalmol 2: 67–74.

35. Gum GG, Gelatt KC, Samuelson DA (1984) Maturation of the retina of the canine neonate as determined by electroretinography and histology. Am J Vet Res 45: 1166–1171.

36. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.

37. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, et al. (2002) The human genome browser at UCSC. Genome Res 12: 996–1006.

38. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, et al. (2010) A method and server for predicting damaging missense mutations. Nat Methods 7: 248–249.

39. Colville CA, Molday RS (1996) Primary structure and expression of the human beta-subunit and related proteins of the rod photoreceptor cGMP-gated channel. J Biol Chem 271: 32968–32974.

40. Ardell MD, Bedsole DL, Schoborg RV, Pittler SJ (2000) Genomic organization of the human rod photoreceptor cGMP-gated cation channel beta-subunit gene. Gene 245: 311–318.

41. Kaupp UB, Seifert R (2002) Cyclic nucleotide-gated ion channels. Physiol Rev 82: 769–824.

42. Sautter A, Zong X, Hofmann F, Biel M (1998) An isoform of the rod photoreceptor cyclic nucleotide-gated channel beta subunit expressed in olfactory neurons. Proc Natl Acad Sci U S A 95: 4696–4701.

43. Chen TY, Peng YW, Dhallan RS, Ahamed B, Reed RR, et al. (1993) A new subunit of the cyclic nucleotide-gated cation channel in retinal rods. Nature 362: 764–767.

44. Wiesner B, Weiner J, Middendorff R, Hagen V, Kaupp UB, et al. (1998) Cyclic nucleotide-gated channels on the flagellum control Ca2+ entry into sperm. J Cell Biol 142: 473–484.

45. Huttl S, Michalakis S, Seeliger M, Luo DG, Acar N, et al. (2005) Impaired channel targeting and retinal degeneration in mice lacking the cyclic nucleotide-gated channel subunit CNGB1. J Neurosci 25: 130–138.

46. Zhang Y, Molday LL, Molday RS, Sarfare SS, Woodruff ML, et al. (2009) Knockout of GARPs and the beta-subunit of the rod cGMP-gated channel disrupts disk morphogenesis and rod outer segment structural integrity. J Cell Sci 122: 1192–1200.

47. Sugimoto Y, Yatsunami K, Tsujimoto M, Khorana HG, Ichikawa A (1991) The amino acid sequence of a glutamic acid-rich protein from bovine retina as deduced from the cDNA sequence. Proc Natl Acad Sci U S A 88: 3116–3119.

48. Zheng J, Trudeau MC, Zagotta WN (2002) Rod cyclic nucleotide-gated channels have a stoichiometry of three CNGA1 subunits and one CNGB1 subunit. Neuron 36: 891–896.

49. Zhong H, Molday LL, Molday RS, Yau KW (2002) The heteromeric cyclic nucleotide-gated channel adopts a 3A:1B stoichiometry. Nature 420: 193–198.

50. Weitz D, Ficek N, Kremmer E, Bauer PJ, Kaupp UB (2002) Subunit stoichiometry of the CNG channel of rod photoreceptors. Neuron 36: 881–889.

51. Shuart NG, Haitin Y, Camp SS, Black KD, Zagotta WN (2011) Molecular mechanism for 3:1 subunit stoichiometry of rod cyclic nucleotide-gated ion channels. Nat Commun 2: 457.

52. Chen TY, Illing M, Molday LL, Hsu YT, Yau KW, et al. (1994) Subunit 2 (or beta) of retinal rod cGMP-gated cation channel is a component of the 240-kDa channel-associated protein and mediates Ca(2+)-calmodulin modulation. Proc Natl Acad Sci U S A 91: 11757–11761.

53. Kaupp UB, Niidome T, Tanabe T, Terada S, Bönigk W, et al. (1989) Primary structure and functional expression from complementary DNA of the rod photoreceptor cyclic GMP-gated channel. Nature 342: 762–766.

54. Pentia DC, Hosier S, Cote RH (2006) The glutamic acid-rich protein-2 (GARP2) is a high affinity rod photoreceptor phosphodiesterase (PDE6)-binding protein that modulates its catalytic properties. J Biol Chem 281: 5500–5505.

55. Ritter LM, Khattree N, Tam B, Moritz OL, Schmitz F, et al. (2011) In situ visualization of protein interactions in sensory neurons: glutamic acid-rich proteins (GARPs) play differential roles for photoreceptor outer segment scaffolding. J Neurosci 31: 11231–11243.

56. Michalakis S, Zong X, Becirovic E, Hammelmann V, Wein T, et al. (2011) The glutamic acid-rich protein is a gating inhibitor of cyclic nucleotide-gated channels. J Neurosci 31: 133–141.

57. Trudeau MC, Zagotta WN (2002) An intersubunit interaction regulates trafficking of rod cyclic nucleotide-gated channels and is disrupted in an inherited form of blindness. Neuron 34: 197–207.

58. Zheng J, Zagotta WN (2004) Stoichiometry and assembly of olfactory cyclic nucleotide-gated channels. Neuron 42: 411–421.

59. Michalakis S, Reisert J, Geiger H, Wetzel C, Zong X, et al. (2006) Loss of CNGB1 protein leads to olfactory dysfunction and subciliary cyclic nucleotide-gated channel trapping. J Biol Chem 281: 35156–35166.

60. Bareil C, Hamel CP, Delague V, Arnaud B, Demaille J, et al. (2001) Segregation of a mutation in CNGB1 encoding the beta-subunit of the rod cGMP-gated channel in a family with autosomal recessive retinitis pigmentosa. Hum Genet 108: 328–334.

61. Simpson DA, Clark GR, Alexander S, Silvestri G, Willoughby CE (2011) Molecular diagnosis for heterogeneous genetic diseases with targeted high-throughput DNA sequencing applied to retinitis pigmentosa. J Med Genet 48: 145–151.

62. Kondo H, Qin M, Mizota A, Kondo M, Hayashi H, et al. (2004) A homozygosity-based search for mutations in patients with autosomal recessive retinitis pigmentosa, using microsatellite markers. Invest Ophthalmol Vis Sci 45: 4433–4439.

63. Becirovic E, Nakova K, Hammelmann V, Hennel R, Biel M, et al. (2010) The retinitis pigmentosa mutation c.3444+1G>A in CNGB1 results in skipping of exon 32. PLoS One 5: e8969.

64. Koch S, Sothilingam V, Garcia Garrido M, Tanimoto N, Becirovic E, et al. (2012) Gene therapy restores vision and delays degeneration in the CNGB1(−/−) mouse model of retinitis pigmentosa. Hum Mol Genet 21: 4486–4496.

65. Scheet P, Stephens M (2006) A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet 78: 629–644.

66. Mowat FM, Breuwer AR, Bartoe JT, Annear MJ, Zhang Z, et al.. (2012) RPE65 gene therapy slows cone loss in Rpe65-deficient dogs. Gene Ther DOI:10.1038/gt.2012.63.

67. Cook NJ, Molday LL, Reid D, Kaupp UB, Molday RS (1989) The cGMP-gated channel of bovine rod photoreceptors is localized exclusively in the plasma membrane. J Biol Chem 264: 6996–6999.

68. Gouy M, Guindon S, Gascuel O (2010) SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol 27: 221–224.