A deletion in nephronophthisis 4 (NPHP4) is associated with recessive cone-rod dystrophy in standard wire-haired dachshund

  1. Anne Caroline Wiik1,
  2. Claire Wade2,3,
  3. Tara Biagi2,
  4. Ernst-Otto Ropstad4,
  5. Ellen Bjerkås4,
  6. Kerstin Lindblad-Toh2,5,6, and
  7. Frode Lingaas1,6
  1. 1 Department of Basic Sciences and Aquatic Medicine, Division of Genetics, Norwegian School of Veterinary Science, 0033 Oslo, Norway;
  2. 2 Broad Institute of Harvard and MIT, Cambridge, Massachusetts 02141, USA;
  3. 3 Center for Human Genetic Research, Massachusetts General Hospital, Boston, Massachusetts 02114, USA;
  4. 4 Department of Companion Animal Clinical Sciences, Norwegian School of Veterinary Science, 0033 Oslo, Norway;
  5. 5 Department of Medical Biochemistry and Microbiology, Uppsala University, 751 24 Uppsala, Sweden

Abstract

Cone-rod dystrophy is a retinal degenerative disorder occurring naturally in man and dog. Here we identify a novel gene for early-onset cone-rod dystrophy in the wire-haired dachshund. For the first time, we use genome-wide association-based Sibling Transmission Disequilibrium Test (sibTDT) analysis of only 13 discordant sib-pairs to identify a single significantly associated 6.5-Mb region (PrawTDT = 4.8 × 10−5, PgenomeTDT = 6 × 10−4) on canine chromosome 5, containing more than 70 genes. Segregation studies using microsatellites in the candidate region including additional meiosis supported the sibTDT analysis but could not further reduce the area. Candidate gene resequencing identified a 180-bp deletion in exon/intron 5 of NPHP4 (nephronophthisis 4, also known as nephroretinin). RT-PCR analysis of NPHP4 in cases and controls showed exon skipping of exon 5, resulting in a truncated protein that retains the binding domain interacting with nephronophthisis 1 (also known as nephrocystin-1) in the kidney but lacks the domain interacting with RPGRIP1 in retina. We suggest that this deletion in the canine NPHP4 gene is the cause of cone-rod dystrophy in the standard wire-haired dachshund. In humans, mutations in NPHP4 have been associated with simultaneous eye and kidney disease. Here we describe the first naturally occurring mutation in NPHP4 without additional kidney disease. Further studies will permit elucidation of the complex molecular mechanism of this retinopathy and the development of potential therapies.

The cone-rod dystrophies (crds) are a heterogeneous group of disorders occurring naturally in man and dog and are characterized by a predominant loss of cone function, with relative preservation of rod function (Udar et al. 2003; Michaelides et al. 2005). A few genes have been associated with autosomal recessively inherited crd in man, including ABCA4 (Allikmets et al. 1997a, b) and RPGRIP1 (Hameed et al. 2003) as well as the loci CORD8 (Khaliq et al. 2000; Ismail et al. 2006) and CORD9 (Danciger et al. 2001). This number of identified loci associated with crd is low compared with the large number of ∼130 genes identified for eye diseases in general (RetNet).

Cone-rod dystrophy is classified as a retinal ciliopathy, together with many other diseases including retinitis pigmentosa, macular degeneration, cone dystrophy, Leber congenital amaurosis (LCA), and Joubert syndrome (JBTS) (Adams et al. 2007). Several groups of genes have been implicated in these studies, including retinitis pigmentosa-1 (RP1), retinitis pigmentosa GTPase regulator (RPGR), retinitis pigmentosa GTPase regulator interacting protein (RPGRIP1), as well as the Usher, Bardet-Biedl, and nephronophthisis genes (Adams et al. 2007). These genes are expressed in the outer segment of the photoreceptors, which is a specialized and modified cilium, where they have a structural role and are crucial for transport function. They also are involved in protein–protein interaction with transport-like proteins (Adams et al. 2007).

Nephronophthisis (NPHP) comprises a group of autosomal recessive cystic kidney diseases, which constitute the most frequent genetic cause for end-stage renal failure in children and young adults. Six genes, NPHP1, INVS (formerly known as NPHP2), NPHP3, NPHP4, IQCB1 (NPHP5), CEP290 (NPHP6), have been identified as mutated in human NPHP. In JBTS, retinal degeneration is associated with mental retardation, cerebellar vermis aplasia/hypoplasia, and NPHP (Parisi et al. 2004). Mutations in CEP290 and RPGRIP1L, homolog of RPGRIP1 and nephronophthisis 4 (NPHP4, also known as nephroretinin and nephrocystin-4) interactor, are involved in JBTS (Valente et al. 2006; Arts et al. 2007). In Senior-Loken syndrome (SLSN), NPHP are associated with retinitis pigmentosa, tapetoretinal degeneration, or retinal dysplasia (Hildebrandt et al. 2007). SLSN is present in about 10% of people with mutations in NPHP1, INVS, NPHP3, or NPHP4 (Caridi et al. 1998; Otto et al. 2002; Olbrich et al. 2003; O'Toole et al. 2007) and >90% of people with mutations in IQCB1 or CEP290 (Otto et al. 2005; Sayer et al. 2006). RPGRIP1 and NPHP4 interact strongly in vitro and in vivo, and they colocalize in the retina (Roepman et al. 2005). Their interaction is disrupted either by mutations in RPGRIP1, found in patients with LCA, or by mutations in NPHP4, found in patients with SLSN (Roepman et al. 2005).

The standard wire-haired dachshund (SWHD), the miniature long-haired dachshund (MLHD), and the pit bull terrier (PBT) are the only dog breeds to date known to be affected by crd (Kijas et al. 2004; Mellersh et al. 2006; Ropstad et al. 2007a). A mutation in canine RPGRIP1 has been demonstrated to cause autosomal recessive crd in MLHDs (Mellersh et al. 2006), while the genetic basis of the crd found in the PBT is still unknown (Kijas et al. 2004).

A family of SWHDs segregating recessively inherited crd was utilized to characterize the phenotype (Ropstad et al. 2007a, b) and, initially, to rule out candidate genes (Wiik et al. 2008). In this study, a whole-genome association of discordant sib-pairs genotyped with the canine single nucleotide polymorphism (SNP) array (Karlsson et al. 2007) was followed by fine-mapping and mutation detection to identify a deletion in NPHP4 that causes crd in these dogs. This mutation is the first nephroretinin mutation to cause retinal disease without kidney involvement.

Results

Genome-wide association analysis

A significant association was detected to chromosome 5 (CFA5) using sibTDT analysis of ∼50,000 SNPs in 13 cases and 13 controls selected as discordant sib-pairs from a pedigree of SWHD (Fig. 1). All samples were genotyped on the Affymetrix v2 canine SNP array. One control was excluded from further analysis based on low call rate (<75%). On the basis of the close relationship of individuals within this pedigree, a Sibling Transmission Disequilibrium Test (sibTDT) analysis as part of the PLINK software (Purcell et al. 2007) was performed, and genome-wide significance (PrawTDT = 4.8 × 10−5 and PgenomeTDT = 0.00060) was achieved over the region CFA5: 61,122,335–67,576,960. Homozygosity was detected in affected cases from 61.10 to 69.07 Mb (Fig. 2A,B). Case-control association identified the same region but reported the highest association to a spurious SNP in the “shoulder” of the linkage region (PrawASSOC = 1.3 × 10−5, PgenomeASSOC = 5.2 × 10−3; Fig. 2B). In this region, allelic frequencies differed maximally between cases and controls; however, one control dog was homozygous for the “affected” allele at the point of association (Chr5:69351798) and the foundation sire (201) was heterozygous in the near vicinity (although not at this SNP).

Figure 1.

Pedigree of the dachshund family segregating crd that was used for genome-wide association. Discordant sibs of half sibs were genotyped (*) on the canine SNP-array, and sibTDT genome-wide association was performed. Microsatellite markers in the associated CFA5 region were used for identification of recombinant breakpoints and haplotypes are shown for each dog.


Figure 2.

Genome-wide association mapping identified a locus on CFA5. (A) The best association was identified by sibTDT analysis (PgenomeTDT = 0.00060), but case-control association was also significant (PgenomeASSOC = 5.2 × 10−3) using PLINK (Purcell et al. 2007). (B) Comparison of sibTDT and single SNP association statistics across the CFA5 region.


Fine-mapping and mutation detection

Fine-mapping and haplotype analysis of microsatellites of the region between 52.8 Mb and 77.0 Mb, in 36 informative offspring, revealed crossovers at 61.4 Mb and 72.0 Mb (Fig. 1), thus it could not further reduce the area of association, despite including additional meiosis. The associated region contains >70 genes that were examined to find possible candidate genes expressed in the retina. DAB1 (55.3 Mb, outside region), AJAP1 (62.0 Mb), VAMP3 (64.4 Mb), and NPHP4 (62.9 Mb) were identified as potential candidate genes. Exon sequencing, including exon–intron boundaries, showed no mutations in the first three genes. PCR products from exon 5 of NPHP4 showed a different size in unaffected (519 bp) and affected (339 bp) individuals and with both bands seen in carriers (Fig. 3A). Exon sequencing showed a deletion of 180 bp in exon 5/intron 5 (location 62,913,591–62,913,770) in the affected animal, where only the first 18 bases of exon 5 were retained (Figs. 3C, 4). Screening of exon 5 in 84 unrelated SWHD revealed 8 carriers, giving a gene frequency of the mutated gene of ∼9.5% in the SWHD population. The exon screening also revealed five SNPs within exons, and 4 SNPs were found in intronic sequence close to exons (Supplemental Table S4). All of the expressed SNPs give rise to amino acid changes, but all are present as polymorphisms in multiple breeds and are not associated to disease within the dachshunds (Supplemental Table S4). No mutations were found in the exon/intron junctions. The putative mutation was not found in any of 32 examined dogs from eight different breeds (Supplemental Table S5).

Figure 3.

A deletion in exon 5 is visualized on the genomic and transcript level. (A) PCR products of exon 5 in affected dogs (339 bp, lanes 1,2), carrier dogs (heterozygotes, lanes 35), and unaffected dog (519 bp, lane 6). (B) PCR amplification of cDNA of one affected (329 bp, lane 1) and one healthy dog (392 bp, lane 2) using primers localized in exon 3 and exon 6. (C) Alignment of sequenced PCR product of mutant and wild-type DNA, where exonic sequence is represented in uppercase letters and intronic sequence in lowercase letters. (D) cDNA.


Figure 4.

Illustration of the 180-bp deletion in exon 5 of the NPHP4 gene at DNA, RNA, and protein level. The mutation leads to exon-skipping of exon 5 and a stop codon in exon 6, which truncates the affected protein at 155 aa. The truncated protein contains the nephrocystin interacting domain but not the RPGRIP1 domain.


Transcript analysis

Comparison of cDNA from affected and unaffected dogs using RT-PCR revealed complete exon-skipping of exon 5 in the affected animals (Fig. 3B,D). Translation of the 4819-bp cDNA sequence of the healthy dogs gives a 1429-amino-acid (aa) protein, while translation of the cDNA of the affected dog results in an early stop codon in exon 6 resulting in a short protein of 155 aa (Fig. 4). The RACE procedure of the 3′- and 5′-ends gave a complete mRNA sequence, “EU659993,” consisting of 4819 bp distributed on 30 exons. Exon 1 (73 bp) and the first 35 bp of exon 2 are untranslated (5′-UTR), while the last 420 bp of exon 30 are 3′-UTR. The NPHP4 exons are located at the minus strand of CFA5 from 62,935,638 to 62,819,750 bp. The amplification of the transcript ends revealed no differences between the affected and the unaffected dogs.

Discussion

Thirteen affected and 13 sex- and age-matched controls were selected for genome-wide association and were genotyped using the Affymetrix version 2 SNP array. After association to CFA5, the haplotypes were used to identify a region of homozygosity among case individuals in the associated region. Both case-control association analysis and sibTDT identified the same general region of association to the disease trait with genome-wide significance. As has been previously demonstrated (Karlsson et al. 2007), a Mendelian recessive trait can be successfully mapped within a single dog breed with fewer than 15 cases and 15 control dogs. However, the SNP with the most significant case-control association probability was not in the region of the best sibTDT association. In this instance, where a relatively small number of individuals were analyzed and where these individuals were predominantly mismatched sibs, the employment of allelic case-control analysis was not the optimal analytical method. This study shows us that in such cases, supportive analyses, such as autozygosity mapping and TDT analysis, may be required to identify the best possible association. Since 6% of the genome is homozygous at the 0.5-Mb level within a breed (Karlsson et al. 2007), simple homozygosity mapping will identify too many loci, and the power offered by information about transmission in a sibTDT test is valuable. The analysis demonstrates that lack of relationship among included individuals is not a prerequisite of successful association mapping and that knowledge of the pedigree structure can be used to better inform statistical results from the analysis. However, it is important to consider that the identified region will be considerably larger than for a case-control study where the typical associated region is 0.5–1 Mb in size (Karlsson et al. 2007; Salmon Hillbertz et al. 2007)

Resequencing of candidate genes within the region of homozygosity of affected dogs was conducted in parallel with fine-mapping. DAB1 coding for disabled homolog 1 (55.3 Mb), AJAP1-adherens junction associated protein (62.0 Mb), and VAMP3-vesicle associated protein (64.4 Mb) were excluded, but resequencing of the NPHP4 gene identified a 180-bp deletion in exon/intron 5 that causes exon skipping and early truncation of the protein in exon 6. This 155-aa, shortened form of NPHP4 contains the binding site for nephronophthisis 1 (NPHP1, also known as nephrocystin 1), which is known to interact with NPHP4 in kidney, but has lost the binding site for RPGRIP1, the protein known to interact with NPHP4 in retina. This is particularly interesting since mutations in NPHP4 have not yet been described as a cause of an eye disease without affecting the kidneys simultaneously. Mutations in RPGRIP1, however, have been reported to be a cause of LCA (Dryja et al. 2001; Gerber et al. 2001), while mutations causing residual RPGRIP1 activity may lead to phenotypes such as RP or crd, which are less severe than LCA (Cremers et al. 2002). Some RPGRIP1 gene mutations are also associated with recessive crd in human (Hameed et al. 2003) and dog (Mellersh et al. 2006).

More specifically, the RPGRIP-interacting domain of NPHP4 has previously been pinpointed to amino acids 591–960 of GenBank entry NP_055917, encoded by human exons 15–21 (Roepman et al. 2005), and is clearly not present in the 155-aa NPHP4 protein of affected dachshunds. The size of the deletion in NPHP4 and the colocalization and interaction with RPGRIP1 may explain why this mutation is causing pure crd and why this disease found in the SWHDs is relatively similar to the one found in MLHD, which has a RPGRIP1 mutation. Time of onset of the initial crd in the affected SWHDs varied from 10 mo to 3 yr of age, with a complete retinal atrophy evident within 6 yr (Ropstad et al. 2007a). The most characteristic opthalmoscopic changes were initial cellophane-like increase in tapetal sheen and marked pigment migration in the nontapetal fundus. The MLHD showed a granular appearance of the tapetal fundus at the age of 26 wk. Electrophysiologically, crd-affected SWHDs showed reduced cone-derived responses already at 5 wk, while the rod-derived responses were less affected. At 52 wk, most affected animals still had recordable electroretinography (ERG) responses. In the MLHDs, the ERG responses were nonrecordable already at the age of 40 wk (Mellersh et al. 2006).

The N-teminal region of NPHP4 is crucial for the interaction with NPHP1 (Mollet et al. 2005). Our findings show that the first 155 amino acids are intact in the mutated protein of the SWHD, and this may provide another explanation why the mutation does not affect the kidneys. The binding region with NPHP1 was known to be located in the first 176 amino acids of NPHP4 (Mollet et al. 2005), but this study refines the region further to the first 155 amino acids—assuming that NPHP4 has the same function in humans and dogs—as no clinical or pathological symptoms of kidney disease were observed in the crd-affected dachshunds (oldest dog followed until 6 yr of age). It will now be interesting to search for mutations in NPHP4 in human crd patients without kidney disease.

This paper demonstrates clearly that small canine pedigrees can be utilized for association mapping to identify mutations in genes of human relevance and that such data can be used to shed light on protein function as well as to generate new candidate genes for human disease.

Methods

Animals and diagnosis

Day blindness was diagnosed in a wire-haired standard dachshund, and a purpose-bred family of dogs confirmed that this most likely was a simple recessively inherited phenotype. Detailed clinical descriptions of the affected family have been reported previously (Ropstad et al. 2007a, b), establishing the disease as an early-onset crd. All procedures used in this study were in accordance with the Norwegian Animal Research Authority (“Forsøksdyrutvalget”) and adhered to the Association for Research in Vision and Ophthalmology Statement for the Use of Animals in Opthalmic Vision and Research.

Samples

Blood samples were collected in EDTA-coated tubes from all the families segregating the disease and from two unaffected, unrelated dogs. Genomic DNA was isolated from the blood samples using DNeasy Tissue kit (Qiagen). Thirteen affected and 13 sex- and age-matched controls (preferably full sibs) were selected for whole-genome SNP analysis (Fig. 1). As some of the dogs were euthanized, a number of samples from various tissues (including retina and muscle) were collected immediately post-mortem and frozen on liquid nitrogen for RNA analysis.

Whole-genome association analysis

From the SWHD family, 13 pairs of cases and controls were chosen. Where possible, mismatched full sibs (n = 10) were used in the analysis to equally distribute population stratification effects among cases and controls. Where a full sib-pair was not available, mismatched half-sibs were used (n = 2 pairs). The founding sire (case 201) and dam (control 208) were included in the analysis.

Genome-wide association analysis was undertaken using the Affymetrix version 2 Dog SNP Array “Platinum Panel” containing 49,663 SNP markers (Affymetrix). SNP genotypes were obtained by following the Affymetrix protocol for the human 500K array protocol but with a smaller hybridization volume to allow for the smaller surface area of the canine array as described in Karlsson et al. (2007). Detailed information on the arrays is available at http://www.broad.mit.edu/mammals/dog/caninearray/.

The sibTDT was carried out on all 49,663 SNPs using the “Dfam” option in the analysis package PLINK (Purcell et al. 2007). Data were filtered to exclude individual calls with BRLMM calling probabilities lower than 0.01, loci with low minor allele frequency (<0.05), loci with low genotyping success rate (failed calls >0.1), and individual dogs with greater than 25% of genotypes missing. After applying these filters, 24,798 SNPs and 25 dogs (13 cases and 12 controls) remained in the analysis. The mean genotyping rate in the individual SWHDs included in the analysis was 96%. Genome-wide significance was ascertained through permutation of phenotypes (n = 10,000 permutations) over all 24,798 SNPs in the analysis (EMP2 in PLINK).

Case-control association mapping using PLINK was applied to the same data with the same filtering parameters to test the efficacy of case-control association strategies in related samples. A region of association spanning from 57.6 to 72.8 Mb was identified using the array analysis.

Fine-mapping

Primers for microsatellite analysis in the region of CFA5 (52.8 and 77.0 Mb) were identified by searching selected sequence windows of ∼200,000 bp for microsatellite motifs. Tetra- and dinucleotide repeat microsatellites with more than 10 and 19 repeats, respectively, were selected. Primers for amplification of selected microsatellites were designed using Primer3 (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) (Supplemental Table S1). One primer in each pair was labeled with fluorescein. For genotyping, PCR was carried out as described below with 28 cycles. The sizes of the alleles were estimated with an automated genotyper (ABI PRISM 3100 Genetic Analyzer, Applied Biosystems) with software for fragment analysis.

Mutation screening

Chromosomal locations of the predicted canine orthologs of the candidate genes were determined by using the http://genome.ucsc.edu website. Primers were designed in introns, using primer3 (Supplemental Table S2), at a distance of ∼100–200 bp from the intron/exon borders, to provide optimal PCR products for bidirectional sequencing. The primers for amplification on each of the canine NPHP4-exons were positioned in introns according to the described human sequence NM_015102 (Supplemental Table S3). The exon-specific primers for amplification of cDNA were based on the predicted canine sequence XM_546745.

PCR amplification reactions were performed using 83 μM of each primer in a 15-μL reaction containing 1.5 μL of DNA prepared as described above, 1× Ammonium PCR buffer containing 1.5 mM MgCl2, 83 μM each of dATP, dCTP, dGTP, and dTTP, and 0.4 units Taq DNA polymerase (Ambion). After an initial denaturation at 95°C for 150 sec, samples were amplified for 34 cycles at 95°C for 30 sec, 58°C for 40 sec, and 72°C for 50 sec, followed by a final extension of 72°C for 330 sec. All PCR products were bidirectionally sequenced using the PCR primers and the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems). Amplification was done in a 10-μL reaction mixture containing 2 μL of BigDye sequencing buffer, 1.5 μL of BigDye sequencing enzyme, 1 μL of PCR product, 0.5 μL of sequencing primer (5 pmol/μL), and 5 μL of dH2O. The sequencing reaction was performed by cycle sequencing with the following protocol: initially 95°C for 1 min, then 28 cycles each of 95°C for 10 sec, 50°C for 15 sec, and 60°C for 2 min. Sequence reaction products were purified by ammonium acetate precipitation and sequenced on the ABI3100.

Transcript analysis

One affected family member was compared with an unaffected English cocker spaniel (ECS) at the mRNA level. Total RNA was isolated from 30 mg of muscle tissue using RNeasy fibrous tissue kit (Qiagen) following the manufacturer’s instructions. The RNA was precipitated using standard procedures with ethanol, sodium acetate, and glycogen. The RNA was dissolved in 10 μL of RNase free water and first-strand cDNA synthesis was performed using Superscript II RT following the manufacturer’s descriptions. The final cDNA product was diluted 1:1 with RNase free water. Exon-specific primers were designed based on the predicted canine NPHP4 sequence XM_546745. PCR with these primers were performed as described above with 1 μL of cDNA as template.

3′- and 5′-RACE was performed to amplify the ends of the NPHP4 transcript in both dogs. RNA was isolated from 30 mg of muscle tissue as described above, precipitated, and dissolved in RNase free water. mRNA for use in the Marathon RACE protocol was isolated from the total RNA using Dynabeads oligo-dT particles (Invitrogen). 5′- and 3′-end RACE was performed with GeneRacer Kit (Invitrogen) and Marathon cDNA Amplification Kit (Clontech), according to the manufacturers’ protocols. PCR and nested PCR was performed using Kit-specific 5′- and 3′-primers with Reverse end Forward gene-specific primers, and cloning of PCR products was performed using TOPO cloning kit. The inserts of the clones were amplified with M13 forward and reverse primer and sequenced with T3 and T7 primers.

Acknowledgments

We acknowledge the Norwegian School of Veterinary Science for funding the project and Casey Gates and the Broad Institute Genetic analysis platform for performing the genotyping. We thank Michael C. Zody for comments on the manuscript and Leslie Gaffney for practical help with the manuscript preparation.

Footnotes

  • 6 Corresponding authors.

    6 E-mail kersli{at}broad.mit.edu; fax (617) 258-0903.

    6 E-mail Frode.Lingaas{at}veths.no; fax +47-22597310.

  • [Supplemental material is available online at www.genome.org. The sequence data from this study have been submitted to GenBank under accession nos. EU707775–EU707790.]

  • Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.074302.107.

    • Received February 6, 2008.
    • Accepted May 19, 2008.

References