# Targeted Next Generation Sequencing Identifies Novel Mutations in RP1 as a Relatively Common Cause of Autosomal Recessive Rod-Cone Dystrophy

Said El Shamieh, Elise Boulanger-Scemama, Marie-Elise Lancelot, Aline Antonio, Vanessa Démontant, Christel Condroyer, Mélanie Letexier, Jean-Paul Saraiva, Saddek Mohand-Saïd, José-Alain Sahel, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Said El Shamieh, Elise Boulanger-Scemama, Marie-Elise Lancelot, Aline Antonio, Vanessa Démontant, et al. Targeted Next Generation Sequencing Identifies Novel Mutations in RP1 as a Relatively Common Cause of Autosomal Recessive Rod-Cone Dystrophy. BioMed Research International, 2015, 2015, pp.485624. 10.1155/2015/485624. hal-01211446

HAL Id: hal-01211446

https://hal.sorbonne-universite.fr/hal-01211446

Submitted on 5 Oct 2015

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
Research Article

# Targeted Next Generation Sequencing Identifies Novel Mutations in *RP1* as a Relatively Common Cause of Autosomal Recessive Rod-Cone Dystrophy

Said El Shamieh,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> Elise Boulanger-Scemama,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> Marie-Elise Lancelot,<fnref n="1" />,<fnref n="2" />,<fnref n="3" />

Aline Antonio,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> Vanessa Démontant,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> Christel Condroyer,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> Mélanie Letexier,<fnref n="4" />

Jean-Paul Saraiva,<fnref n="4" /> Saddek Mohand-Saïd,<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="5" /> José-Alain Sahel,<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="5" />,<fnref n="6" />,<fnref n="7" />,<fnref n="8" />

Isabelle Audo,<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="5" />,<fnref n="8" /> and Christina Zeitz<fnref n="1" />,<fnref n="2" />,<fnref n="3" />

Correspondence should be addressed to Isabelle Audo; isabelle.audo@inserm.fr and Christina Zeitz; christina.zeitz@inserm.fr

Received 27 June 2014; Accepted 10 July 2014

Academic Editor: Calvin Yu-Chian Chen

Copyright © 2015 Said El Shamieh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We report ophthalmic and genetic findings in families with autosomal recessive rod-cone dystrophy (arRCD) and *RP1* mutations. Detailed ophthalmic examination was performed in 242 sporadic and arRCD subjects. Genomic DNA was investigated using our customized next generation sequencing panel targeting up to 123 genes implicated in inherited retinal disorders. Stringent filtering coupled with Sanger sequencing and followed by cosegregation analysis was performed to confirm biallelism and the implication of the most likely disease causing variants. Sequencing identified 9 *RP1* mutations in 7 index cases. Eight of the mutations were novel, and all cosegregated with severe arRCD phenotype, found associated with additional macular changes. Among the identified mutations, 4 belong to a region, previously associated with arRCD, and 5 others in a region previously associated with adRCD. Our prevalence studies showed that *RP1* mutations account for up to 2.5% of arRCD. These results point out for the necessity of sequencing *RP1* when genetically investigating sporadic and arRCD. It further highlights the interest of unbiased sequencing technique, which allows investigating the implication of the same gene in different modes of inheritance. Finally, it reports that different regions of *RP1* can also lead to arRCD.

## 1. Introduction

Rod-cone dystrophy (RCD), also known as retinitis pigmentosa, is a heterogeneous group of inherited disorders affecting primary rod photoreceptors in the majority of cases with secondary cone degeneration [1, 2]. Population-based studies showed that 1 in 4,000 individuals is affected around the

world [1]. Patients diagnosed with RCD initially complain of night blindness due to rod dysfunction followed by progressive visual field constriction, abnormal color vision, and eventually loss of central vision due to cone photoreceptor involvement [1].

RCD is inherited as a Mendelian trait in most cases [3]. On the basis of its mode of inheritance and prevalence,

<footnote n="1">¹INSERM, U968, 75012 Paris, France</footnote>

<footnote n="2">²Sorbonne Universités, UPMC University, Paris 06, UMR_S 968, Institut de la Vision, 75012 Paris, France</footnote>

<footnote n="3">³CNRS, UMR_7210, 75012 Paris, France</footnote>

<footnote n="4">⁴IntegraGen SA, Genopole Campus 1, Building G8, 91030 Evry, France</footnote>

<footnote n="5">⁵Centre Hospitalier National d'Ophtalmologie des Quinze-Vingt's, DHU ViewMaintain, INSERM-DHOS CIC 1423, 75012 Paris, France</footnote>

<footnote n="6">⁶Fondation Ophtalmologique Adolphe de Rothschild, 75019 Paris, France</footnote>

<footnote n="7">⁷Académie des Sciences-Institut de France, 75006 Paris, France</footnote>

<footnote n="8">⁸University College London Institute of Ophthalmology, 11-43 Bath Street, London EC1V 9EL, UK</footnote>
RCD can be divided into 3 groups: autosomal dominant (ad) (30-40%), autosomal recessive (ar) (50-60%), and X-linked (xl) (5-15%) [3]. To date, mutations in at least 53 genes were reported to cause nonsyndromic RCD (till 25 June 2014, https://sph.uth.edu/retnet/). Prevalence studies revealed *rhodopsin* (*RHO*), *retinitis pigmentosa* *GTPase regulator* (*RPGR*), and *usherin* (*USH2A*) as being the most frequently mutated genes in adRCD [4, 5], xlRCD [4], and arRCD, respectively [6]. Of note is that many other genes with lower prevalence are also implicated in the genetic etiology of RCD [7, 8]. Mutations in *RP1* were first shown to cause adRCD [9-11]; however, since 2005, articles have shed light on its implication in arRCD etiology [12-20]. *RP1* mutations were shown to account for ≈5.5% and ≈1% of adRCD and arRCD cases, respectively [8-20]. Interestingly, Avila-Fernandez et al. [12] reported that a founder nonsense mutation in the Spanish population p.Ser542* is responsible for 4.5% of arRCD cases suggesting that *RP1* mutations are more prevalent in arRCD than previously thought [12].

Retinitis pigmentosa 1 (*RP1*) is a photoreceptor-specific gene encoding a protein regulated by oxygen [10]. *RP1* protein is required for correct orientation and higher-order stacking of outer segment disks [21] and was shown to be part of the photoreceptor axoneme [22]. *RP1* localizes to the connecting cilia of photoreceptors and may assist in maintenance of ciliary structure or transport down the photoreceptor [22]. Like many retinal degeneration genes, the mechanism by which mutations in *RP1* lead to photoreceptor cell death is still unclear.

We developed an unbiased and time-efficient retinal gene next generation sequencing array (NGS), which was further revised and improved to target more than 120 genes implicated in inherited retinal diseases (IRDs) (list available upon request) [23]. Using this NGS panel, we screened a total of 242 subjects with sporadic and recessive RCD in order to detect disease causing mutations and to report the prevalence of pathogenic mutations in *RP1* causing arRCD.

## 2. Methods

### 2.1. Ethics Statement and Clinical Diagnosis of Rod-Cone Dystrophy.
The study protocol adhered to the tenets of the Declaration of Helsinki and was approved by the local Ethics Committee (CPP, Ile de France V). Informed written consent was obtained from each study participant. Index patients underwent full ophthalmic examination as previously described [23].

### 2.2. Targeted Next Generation Sequencing.
A cohort of 242 subjects affected with sporadic and arRCD was investigated in the present study. Prior to NGS screening, molecular genetic analysis with microarray (Asper Ophthalmics, Tartu, Estonia), followed by direct Sanger sequencing of *EYS* and *C2orf71* (major and minor genes implicated in RCD, newly discovered at the beginning of our study), was performed in 201 index subjects (82%) [2, 24]. As *RPGR* exon *ORF15* (MIM 312610) is not targeted by existing NGS panels, we excluded mutations in this "hot spot" by Sanger sequencing.

Although our NGS panel was selected from the SureSelect Human All Exon Kits Version 4 (Agilent, Massy, France), this design was improved after analyzing the first 83 subjects with sporadic and arRCD. More precisely, a total of ≈300 Kb regions were added in order to cover all the previously nontargeted regions. Thus, whereas the first design covered the exons and the flanking intronic regions of 120 genes implicated in IRDs, the second covered 123 genes in total. The eArray web-based probe design tool was used for this purpose (https://earray.chem.agilent.com/earray). All probes were designed and synthesized by Agilent Technologies (Santa Clara, CA, USA). Sequence capture, enrichment, and elution were performed according to Agilent's instructions. The complete details were described elsewhere [23].

### 2.3. Assembly and Variant Calling.
Sequence reads were aligned to the reference human genome (UCSC hg19) using CASAVA1.7 software (Illumina) and the ELANDv2 alignment algorithm. Sequence variation annotation was performed using the IntegraGen in-house pipeline, which consisted of gene annotation (RefSeq), detection of known single nucleotide polymorphisms (dbSNP 135) followed by mutation characterization (missense, intronic, synonymous, nonsense, splice site, and insertions/deletions).

### 2.4. Quality Control and Coverage Assessment.
The first NGS retinal panel harbored 120 IRDs genes, encompassing 321,240 kb length per sample. However, after improvement, the same panel contained ≈600 Kb and covered 123 IRD genes. The depth of coverage was calculated by counting the number of sequenced bases mapping to the target regions. Mean depth of coverage was calculated per base pair for all samples; however, only the results of subjects having *RP1* mutations were shown.

### 2.5. Discrete Filtering of Annotated Variants.
In order to identify disease causing mutations among nonpathogenic single nucleotide polymorphisms, we used a filtering approach against a set of polymorphisms that are available in the public databases: dbSNP 137, 1000 Genomes Project [25], HapMap [26], and Exome Variant Server [27] with removal of variants with a minor allele frequency (MAF) ≥ 0.005 in case of presumed autosomal recessive mode of inheritance.

### 2.6. Pathogenicity Assessment.
We stratified candidate mutations based on their functional class by giving a priority to frameshifts, stop codons, and disruptions of canonical splice sites variants [28]. For missense changes, amino acid conservation across 46 different species was studied using the UCSC Genome Browser [29]. If no amino acid change was found, then the residue was considered as "highly conserved." If a different change was seen in less than four species and not in the primates, then it was considered as "moderately conserved" and if a change was present in 5-7, it was considered as "marginally conserved"; otherwise, the amino acid residue was considered as "not conserved." Pathogenic prediction was performed using two software programs: Polyphen2 [30] and SIFT [31], based on species/homologue conservation,
putative structural domains, and 3D structures (if available). Analysis of potential splice site variant consequences when relevant was done using human splicing finder [32].

2.7. *Known Genotype-Phenotype Correlations.* The search for previous genotype-phenotype associations was done by searching numerous literature databases, including Online Mendelian Inheritance in Man (http://omim.org/), Human Gene Mutation Database [33], Leiden Open Variation Database [34], and RetNet (https://sph.uth.edu/retnet/).

2.8. *Validation of the Identified Variants and Cosegregation Analyses.* Sanger sequencing was performed to validate disease causing mutations in *RPI*. The respective primer information can be communicated on request. In addition, blood samples were collected from additional family members when possible and cosegregation analyses on extracted DNA were performed as previously described [35, 36].

## 3. Results

3.1. *Clinical Data.* Clinical data are summarized in Table 1. Among identified patients, 5 were females, 2 were male, and ages at time of examination ranged from 25 to 42. All patients were diagnosed before age 20 mostly based on night blindness from early childhood and secondary central vision loss. They all showed severe RCD with constricted visual fields, no detectable responses on full field electroretinogram, and both peripheral involvement and macular involvement (Figure 1 presents fundus pictures of patient II.1 (CIC01245) in family F752 as an example). Comparing visual acuity and visual fields for these arRCD patients with those of adRCD cases published by Audo and coworkers [8], we noticed a more severe phenotype in recessive cases. However, more cases with *RPI* mutations would be needed to draw statistical conclusion.

3.2. *Sequencing Statistics.* In index patients, the overall sequencing coverage of the target regions was ≥88% for a 25X depth of coverage in each of the chromosomes (Figure 2(a)), resulting a mean sequencing depth of 299 times per base. Mean sequencing results per base in each target chromosome gene regions were shown in Figure 2(b). It is of importance to mention that <1% of target regions were not covered at all. These were fragments of 120 bp belonging in 66% of the cases only to a fraction of an exon. The remaining uncovered targets corresponded each to an entire exon in genes such as *CHM*, *PDZD7*, *RP9*, *CC2D2A*, *IMPDH1*, *CNGA1*, and *EYS*.

3.3. *Detection of Disease Causing Mutations in RPI Gene.* After data filtering, the total number of putative disease causing variants was reduced by 99.3%. Thus, in total, filtering enriched the percentage of putative disease causing mutations from 0.7% (25/339 variants) to 33.3% (9/25 variants) in the 7 subjects presented here (Table 2). These subjects exhibit *RPI* mutations in the last exon 4 that are predicted to lead to a premature stop codon. We found 9 pathogenic mutations in *RPI* among which one (p.Ser542* in CIC00445) was already

reported by Avila-Fernandez et al. [12] as a founder nonsense mutation in the Spanish population, responsible for 4.5% of arRCD. Although F303 is from French origin, we cannot exclude the possibility of a founder effect of p.Ser542* in our cohort.

Patient family F303: II.1 (CIC00445) was found to carry compound heterozygous variants: a nonsense mutation c.1625C>G, leading to a predicted premature stop (p.Ser542*) and a deletion c.4587_4590delTAAG leading to a frameshift and a premature termination codon (p.Ser1529Argfs*9) (Table 2, Figure 3). Patient family F752: II.1 (CIC01245) was also found to carry compound heterozygous variants: a 1bp duplication c.2025dupA leading to p.Ser676Ilefs*22 and a 1bp deletion c.2377delA leading to p.Arg793Glufs*55 (Table 2).

Patients from family F335: III.1 (CIC00491), family F674: III.6 (CIC01106), family F782: II.5 (CIC01300), family F1941: III.1 (CIC04130), and family F3110: III.5 (CIC05941) were found to carry homozygous deletions c.4089_4092delAAGA leading to p.Arg1364Valfs*8; c.1205delG leading to p.Gly402Alafs*7; c.1719_1723delCTCAA leading to p.Ser574Cysfs*7; c.1329delG leading to p.Lys443Asnfs*12; and c.2391_2392delAA leading to p.Asp799*, resp.) (Table 2 and Figure 3). It is important to note that consanguinity was reported in families F335, F674, F782 and F1941.

All *RPI* mutations detected by NGS were further validated by Sanger sequencing. All variants cosegregated with the phenotype in available family members. Based on the previous findings, the measured prevalence of *RPI*-associated arRCD in this cohort is ≈2.5%.

## 4. Discussion

The current study further demonstrates the usefulness of NGS as a comprehensive genetic diagnostic tool for IRDs with further impact on patients counseling and participation for potential therapeutic trials. Our study applied to a large cohort of sporadic and autosomal recessive cases of RCD identifies 8 novel mutations in a gene not classically screened in arRCD by other methods such as Sanger sequencing or microarray analysis, outlining the interest of this massive parallel sequencing method. Consequently, a prevalence of *RPI* mutation in 2.5% of sporadic or arRCD cases in the European population is herein reported.

*RPI* is a 15 kb single copy gene clustering the small arm of the chromosome 8 (8q12.1). It encodes a 2506 amino acid protein having a molecular weight of 241 kDa containing a *Drosophila melanogaster* bifocal (BIF) (amino acid 486–635) and two doublecortin domains. Whereas the BIF domain helps to maintain the photoreceptor morphogenesis, doublecortin domains bind microtubules and regulate their polymerization [22]. Along with *RP1L1* (Retinitis Pigmentosa 1-like 1, another retinal-specific protein), *RP1* plays essential and synergistic roles in outer segment morphogenesis of rod photoreceptors [22].

To date, at least 50 mutations in *RPI* were identified in RCD [8, 12–20], the majority of which are located in its
<table><caption>TABLE 1: Clinical data of the 7 index patients with RPI recessive mutations.</caption><thead><tr><th>Patient</th><th>Age at time of testing</th><th>Age of onset</th><th>Sex</th><th>Family history</th><th>Symptoms at time of diagnosis</th><th>BCVA OD/OS With refraction</th><th>Color vision (15 desaturated Hue)</th><th>Binocular kinetic visual field (III4e stimulus)</th><th>FF and mfERG</th><th>Fundus examination</th><th>FAF</th><th>Sd-OCT</th></tr></thead><tbody><tr><td>F303: II.1 (CIC00445)</td><td>42</td><td>6</td><td>F</td><td>No other affected FM, from France.</td><td>Night blindness</td><td>Hand motion in both eyes</td><td>Impossible due to low vision</td><td>Reduced to peripheral islands of perception</td><td>Both undetectable</td><td>Pale optic disc narrowed blood vessels, macular atrophic changes, and optic nerve drusen</td><td>Hypoautofluorescence in the macular region</td><td>Thinning of outer retina in the macular region</td></tr><tr><td>F335: III.1 (CIC00491)</td><td>36</td><td>3</td><td>M</td><td>Two other brothers affected; parents first cousins</td><td>Night blindness and rapid decreased vision</td><td>LP in both eyes</td><td>Impossible due to low vision</td><td>Impossible due to low vision</td><td>Both undetectable</td><td>Widespread RPE changes and retinal atrophy in both the periphery and the macular area</td><td>Widespread loss of FAF</td><td>Widespread thinning of outer retina</td></tr><tr><td>F674: III.6 (CIC01106)</td><td>25</td><td>19</td><td>F</td><td>Parents first cousins from Turkey, one female and male cousins affected also from a con- sanguineous union</td><td>Night blindness and decreased vision</td><td>HM -3 (-1) 20/160 -3 (-0.50) 0°</td><td>Dyschromatopsia with no specific axis</td><td>Reduced to 5 central degrees</td><td>Both undetectable</td><td>Well-colored optic disc and no narrowing of retinal vessels; RPE changes in the periphery and macular atrophic changes</td><td>Hypoautofluorescence in the macular region and outside the vascular arcades</td><td>Thinning of outer retina in the macular region</td></tr><tr><td>F752: II.1 (CIC01245)</td><td>31</td><td>Early teens</td><td>F</td><td>Two sisters affected</td><td>Night blindness</td><td>20/63 plano (-3) 180° 20/50 plano (-1.75) 180°</td><td>Deutan defect on both eyes</td><td>Reduced to 10° × 20°</td><td>Both undetectable</td><td>Pale optic disc head, narrowed retinal vessels, and RPE changes in the periphery with some macular atrophic changes</td><td>Hypoautofluorescence in the macular region and outside the vascular arcades</td><td>Thinning of outer retina in the macular region</td></tr><tr><td>F782 II.5 (CIC01300)</td><td>27</td><td>9</td><td>M</td><td>Parents from Algeria, first cousins</td><td>Night blindness and decreased vision</td><td>20/50 -9.25 (-2.50) 15° 20/50 -9 (-1.75) 100°</td><td>—</td><td>Reduced to the 10° central degrees</td><td>Both undetectable</td><td>Mild optic disc pallor, atrophic macular changes, and peripheral pigment deposits</td><td>Hypoautofluorescence in the macular region</td><td>Thinning of outer retina in the macular region</td></tr></tbody></table>
<table><caption>TABLE 1: Continued.</caption><thead><tr><th>Patient</th><th>Age at time of testing</th><th>Age of onset</th><th>Sex</th><th>Family history</th><th>Symptoms at time of diagnosis</th><th>BCVA OD/OS With refraction</th><th>Color vision (15 desaturated Hue)</th><th>Binocular kinetic visual field (III4e stimulus)</th><th>FF and mfERG</th><th>Fundus examination</th><th>FAF</th><th>Sd-OCT</th></tr></thead><tbody><tr><td>F1941: III.1 (CIC04130)</td><td>30</td><td>childhood</td><td>F</td><td>Parents from Algeria, first cousins</td><td>Night blindness</td><td>20/100 –4.25 (-1.25) 150°</td><td>Normal at the saturated test</td><td>Reduced to the 10 central degrees</td><td>Both undetectable</td><td>Well-colored optic disc but narrowed retinal vessels; RPE changes in the periphery and macular atrophic changes</td><td>Hypoautofluorescence in the macular region and outside the vascular arcades</td><td>Thinning of outer retina in the macular region</td></tr><tr><td>F3110: III.5 (CIC05941)</td><td>27</td><td>5</td><td>F</td><td>One cousin on mother side may have decreased RCD vision</td><td>Night blindness and decreased vision</td><td>20/125 +2 (-2) 95°</td><td>Dyschromatopsia with no specific axis</td><td>Reduced to the 10 central degrees</td><td>Both undetectable</td><td>Pale optic disc, narrowed retinal vessels, and RPE changes in the periphery with some macular atrophic changes</td><td>Hypoautofluorescence in the macular region and outside the vascular arcades</td><td>Thinning of outer retina in the macular region</td></tr></tbody></table>

F: female, FM: family member, M: male, BCVA: best corrected visual acuity; OD: ocula dextra (right eye); OS: ocula sinistra (left eye); FF and mfERG: full-field and multifocal ERG; FAF: fundus autofluorescence; Sd-OCT: spectral domain optical coherence tomography; RPE: retinal pigment epithelium; LP: light perception; HM: hand motion.
<table><caption>TABLE 2: List of mutations detected by next generation sequencing after applying relevant filters.</caption><thead><tr><th>Patient</th><th>Gene</th><th>Exon</th><th>Allele state</th><th>Nucleotide exchange</th><th>Protein effect</th><th>rs ID</th><th>Conservation</th><th>Polyphen 2</th><th>SIFT</th><th>Pathogenicity</th><th>Note</th></tr></thead><tbody><tr><td rowspan="2">F303: II.1 (CIC00445)</td><td>NPHP4</td><td>12</td><td>HTZ</td><td>A&gt;G</td><td>p.Ser481Asn</td><td>no</td><td>N C</td><td>B</td><td>T</td><td>Neutral</td><td>Disease causing R M [12]</td></tr><tr><td>RPI</td><td>4</td><td>HTZ</td><td>c.1625C&gt;G</td><td>p.Ser542*</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td rowspan="4">F335: III.1 (CIC00491)</td><td>RPI</td><td>4</td><td>HTZ</td><td>c.4587_4590delTAAG</td><td>p.Ser1529Argfs*9</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>N M</td></tr><tr><td>PROM1</td><td>4</td><td>HTZ</td><td>T&gt;C</td><td>p.Ile178Val</td><td>—</td><td>N C</td><td>B</td><td>T</td><td>Neutral</td><td>Disease causing</td></tr><tr><td>GPR98</td><td>29</td><td>HTZ</td><td>G&gt;A</td><td>p.Arg2128Gln</td><td>rs149390094</td><td>N C</td><td>B</td><td>T</td><td>Neutral</td><td>N M</td></tr><tr><td>RPI</td><td>4</td><td>HMZ</td><td>c.4089_4092delAAGA</td><td>p.Arg1364Valfs*8</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td>F674: III.6 (CIC01106)</td><td>USH2A</td><td>39</td><td>HTZ</td><td>T&gt;G</td><td>p.Ser2450Arg</td><td>No</td><td>H C</td><td>P D</td><td>D</td><td>Probably disease causing</td><td>N M</td></tr><tr><td rowspan="2">F752: II.1 (CIC01245)</td><td>RPI</td><td>4</td><td>HMZ</td><td>c.1205delG</td><td>p.Gly402Alafs*7</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td>USHIC</td><td>17</td><td>HTZ</td><td>G&gt;A</td><td>p.Arg477Trp</td><td>TMP_ESP_II_17532053</td><td>H C</td><td>P D</td><td>D</td><td>Probably disease causing</td><td>N M</td></tr><tr><td rowspan="4">F782: II.5 (CIC01300)</td><td>PDE6B</td><td>10</td><td>HTZ</td><td>T&gt;C</td><td>p.Thr432Ile</td><td>—</td><td>H C</td><td>B</td><td>T</td><td>Uncertain pathogenicity</td><td>N M</td></tr><tr><td>RPI</td><td>4</td><td>HTZ</td><td>c.2025dupA</td><td>p.Ser671lefs*22</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td>RPI</td><td>4</td><td>HTZ</td><td>c.2377delA</td><td>p.Arg793Glufs*55</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td>RPI</td><td>4</td><td>HMZ</td><td>c.1719_1723delCTCAA</td><td>p.Ser574Cysfs*7</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td rowspan="4">F194i: III.1 (CIC04130)</td><td>TULP1</td><td>5</td><td>HMZ</td><td>c.395_418dup</td><td>p.Asp124_Glu131del</td><td>rs63749128</td><td>—</td><td>—</td><td>—</td><td>—</td><td>N M</td></tr><tr><td>PCDH15</td><td>33</td><td>HTZ</td><td>C&gt;T</td><td>p.Arg188His</td><td>rs145851144</td><td>N C</td><td>B</td><td>T</td><td>Neutral</td><td>N M</td></tr><tr><td>C2orf71</td><td>1</td><td>HTZ</td><td>C&gt;A</td><td>p.Arg656Ser</td><td>rs201980758</td><td>N C</td><td>B</td><td>T</td><td>Neutral</td><td>N M</td></tr><tr><td>CACNA2D4</td><td>Exon37- Intron 37</td><td>HTZ</td><td>C&gt;T</td><td>—</td><td>rs80092457</td><td>N C</td><td>—</td><td>—</td><td>Neutral</td><td>N M</td></tr><tr><td rowspan="4">F3110: III.5 (CIC05941)</td><td>RPI</td><td>4</td><td>HMZ</td><td>c.1329delG</td><td>p.Lys443Asnfs*12</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr><tr><td>EYS</td><td>6</td><td>HTZ</td><td>C&gt;T</td><td>p.Ser326Asn</td><td>rs112822256</td><td>N C</td><td>B</td><td>T</td><td>Neutral</td><td>N M</td></tr><tr><td>MERTK</td><td>8</td><td>HTZ</td><td>C&gt;G</td><td>p.Arg421Trp</td><td>rs138908058</td><td>N C</td><td>B</td><td>D</td><td>Neutral</td><td>N M</td></tr><tr><td>PRPF6</td><td>21</td><td>HTZ</td><td>A&gt;G</td><td>p.Val915Met</td><td>rs139778757</td><td>M C</td><td>P D</td><td>D</td><td>Uncertain pathogenicity</td><td>N M</td></tr><tr><td rowspan="4">F3110: III.5 (CIC05941)</td><td>TULP1</td><td>14</td><td>HTZ</td><td>G&gt;A</td><td>p.Ala496Thr</td><td>rs141980901</td><td>M C</td><td>B</td><td>D</td><td>Neutral</td><td>N M</td></tr><tr><td>EYS</td><td>26</td><td>HTZ</td><td>G&gt;A</td><td>p.Lys1365Glu</td><td>rs16895519</td><td>N C</td><td>B</td><td>D</td><td>Neutral</td><td>N M</td></tr><tr><td>MERTK</td><td>18</td><td>HTZ</td><td>G&gt;C</td><td>p.Glu823Gln</td><td>rs55924349</td><td>M C</td><td>B</td><td>D</td><td>Neutral</td><td>N M</td></tr><tr><td>RPI</td><td>4</td><td>HMZ</td><td>c.2391_2392delAA</td><td>p.Asp799*</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Disease causing</td><td>N M</td></tr></tbody></table>

Probably disease causing mutations are highlighted in bold.

B: benign, HMZ: homozygous, HTZ: heterozygous, M C: marginally conserved, N C: not conserved, N M: novel mutation, R M: recurrent mutation, T: tolerated, P D: possibly damaging.
<figure><img src="image_3.png" /><figcaption>FIGURE 1: Ophthalmic features of family F752: II.1 (CIC01245): fundus color photographs ((a) and (d) for right and left eye resp.), autofluorescence ((b) and (e) for right and left eye resp.), and spectral domain optical coherence tomography horizontal macula scans ((c) and (f) for right and left eye resp.), showing severe rod-cone dystrophy signs with macular involvement.</figcaption></figure>

last exon (exon 4) and shown to be transmitted in an auto-
somal dominant mode of inheritance. Most of *RP1* disease
causing variants represent nonsense mutations, deletions, or
insertions. In mammalian genes, nonsense mutations lead to
unstable mRNAs that are degraded by nonsense-mediated

decay (NMD). However, exceptions might arise when prema-
ture stop codons occur in the last exon [37]. These variants are
thought to abolish RP1 function by resulting in a truncated
protein lacking important functional domains although still
able to interact with some of its protein partner(s) [21]. The
<figure><img src="image_4.png" /><figcaption>FIGURE 2: Sequencing statistics in index patients. (a) The overall sequencing coverage of the target regions at 25X depth of coverage is shown in each of the chromosomes. No values were indicated for chromosomes 13, 18, 21, and 22 as they were not targeted. The term chromosome 23 was used to designate the X chromosome. F1941: III.1 (CIC04130) and F3110: III.5 (CIC05941) showed the lowest coverage results. (b) The average mean depth per base pair is shown for each of the chromosomes. Most targets showed coverage around 300 times.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>FIGURE 3: Pedigrees of seven families with *RPI* mutations underlying autosomal recessive rod-cone dystrophy. Affected and unaffected individuals are represented by shapes filled with black and white colors, respectively. Men and women are indicated by squares and circles, respectively. Index subjects are marked by ^. Consanguinity is marked by a double horizontal line.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>FIGURE 4: Schematic presentation of *RP1* disease causing mutations. Disease causing mutations were represented based on the classification by Chen and coworkers [13]. Mutations responsible for recessive arRCD were shown in the upper half, whereas mutations causing adRCD were shown in the lower half. p.Gly402Alafs*7, p.Lys443Asnfs*12, p.Arg1364Valfs*8, and p.Ser1529Argfs*9 belong to class III. Although p.Ser574Cysfs*7, p.Ser676Ilefs*22, p.Arg793Glufs*55, and p.Asp799* are class II mutations, these variants do not cause adRCD but arRCD instead. Amino acid modifications shown in red and blue represent novel frameshift or nonsense mutations and the recurrent p.Ser542* mutation respectively. Protein localization of p.Ser542* was highlighted in blue as it marked a recurrent mutation. adRCD: autosomal dominant: rod-cone dystrophy, arRCD: autosomal recessive rod-cone dystrophy, BIF: drosophila melanogaster bifocal.</figcaption></figure>

latter observation is supported by finding that *RP1* mutant mRNA is expressed in a human cell line carrying a homozygous p.Arg677* mutation [21].

Based on Chen et al. [13], *RP1* truncating mutations leading to arRCD or adRCD can be divided into four distinct groups. Class I is composed of truncating mutations located in exons 2 and 3. These variants are sensitive to NMD and thus are considered as true loss-of-function alleles (Figure 4) [13]. Class II involves truncating mutations that are located in a spot between codons 500 and 1053 in exon 4 [13], the so called “*RP1* hot spot.” The “hot spot” variants tend to be insensitive to NMD process and thus result in a protein with a potential dominant negative effect leading to adRCD (Figure 4) [13]. Class III includes truncating mutations insensible to NMD located between codons 264 and 499 and between codons 1054 to 1751 in exon 4. These truncating proteins result in a loss of function leading to arRCD (Figure 4) [13]. Finally, class IV includes protein-truncating mutations near the 3' end of the fourth exon (Figure 4) [13]. Most likely, the resulting proteins display only a minor loss of their C-terminal portion, preserving the majority of functional domains and keeping a residual activity. According to the classification of Chen et al. [13], p.Gly402Alafs*7, p.Lys443Asnfs*12, p.Arg1364Valfs*8, and p.Ser1529Argfs*9 belong to class III (Figure 4).

The predicted physiopathology for p.Ser542*, p.Ser574Cysfs*7, p.Ser676Ilefs*22, p.Arg793Glufs*55, and p.Asp799* is more complex. According to Chen's classification, these frameshift deletions and nonsense mutations should belong to class II, previously only associated with adRCD. However, herein, they are causing presumably arRCD (Figure 4). To further confirm these findings, clinical and genetic testing of the reported unaffected parents should be done.

Based on the previous findings, we speculate that the classification by Chen and coworkers does not hold true for all mutations. Supporting this statement, Avila-Fernandez et al. [12] reported the same nonsense mutation (p.Ser542*) found in (F303: II.1 (CIC00445)) and located at the 5' end of the “hot spot” to cause *arRP* in the Spanish population [12]. These observations are of interest as they point out for an implication of hot spot region for adRCD-*RP1* mutation also in case of arRCD. Future studies will need to clarify why some class II mutations lead to adRCD and others to arRCD.

Patients with arRCD and *RP1* mutations show a more severe disease than adRCD-RP1 mutant patients with macular atrophy in all our cases. This was first outlined by Lafont et al. [17]. When patients are presenting with late, severe disease, the diagnostic distinction between RCD, with
initial rod involvement, and cone-rod dystrophy (CRD) with initial cone involvement is difficult. Of note is that one of the patients (CIC01300) in the present study was initially classified as possibly having severe CRD and his diagnosis was actually revisited after NGS results. This also outlines the interest of unbiased massive parallel sequencing methods for a more precise clinical diagnostic in case of end stage disease. This point will most likely become even more critical with the perspective of therapeutic trials.

**4.1. Strength and Limitations.** We estimate that 1% of our target regions were not covered. Partially uncovered exons are a real common issue when capturing the DNA sequences using commercially available probes; this bias might imply a loss of some candidate variants. However, we found that rate of 1% is very reasonable when compared with other NGS panels. In addition, in order to exclude the possibility of finding other candidate variants, we have sequenced by Sanger method the majority of these regions. Five of our patients carried homozygous *RP1* mutations. For four of the subjects carrying homozygous variants, namely CIC00491, F335; CIC01106, F674; CIC01300, F782 and CIC04130, F1941; co-segregation analysis needs to be done to confirm autosomal recessive inheritance but we do not have access to parent's DNA. CIC05941 was the only one not to report clear consanguinity in the family, and we cannot exclude the possibility of a large deletion on the second allele of *RP1* gene. Again, DNA of the father, not available for us, would be helpful to prove autosomal recessive inheritance and the homozygous state of the mutation.

In conclusion, we have reported 9 mutations in *RP1* of which 8 were novel causing *arRCD* [8, 12–20]. Interestingly, a prevalence of ≈2.5% points out for the necessity of sequencing *RP1* in sporadic and recessive cases of *RCD*. Further functional studies are needed to understand the impact of *RP1* structure on its function at the molecular level; such a step would strengthen our knowledge in the physiology of retinal photoreceptors.

**Conflict of Interests**

The authors declare that there is no conflict of interests regarding the publication of this paper.

**Authors’ Contribution**

Isabelle Audo and Christina Zeitz contributed equally to this work.

**Acknowledgments**

The authors express their sincere gratitude to the families who participated in this study and to the clinical staff for their help in clinical data and DNA collection. This work was supported by Foundation Voir et Entendre (CZ), Prix Dalloz for "la recherche en ophtalmologie" (CZ), Foundation Fighting Blindness (FFB) [CD-CL-0808-0466-CHNO] (IA) and FFB center (FFB grant C-GE-0912-0601-INSERM02),

Prix de la Fondation de l'Œil (IA), Ville de Paris and Region Ile de France, and the French State program "Investissements d'Avenir" managed by the Agence Nationale de la Recherche [LIFESENSES: ANR-10-LABX-65].

**References**

[1] D. T. Hartong, E. L. Berson, and T. P. Dryja, "Retinitis pigmentosa," *The Lancet*, vol. 368, no. 9549, pp. 1795-1809, 2006.

[2] I. Audo, M. Lancelot, S. Mohand-Saïd et al., "Novel C2orf71 mutations account for approximately ~1% of cases in a large French arRP cohort," *Human Mutation*, vol. 32, no. 4, pp. E2091-E2103, 2011.

[3] A. Anasagasti, C. Irigoyen, O. Barandika, A. López de Munain, and J. Ruiz-Ederra, "Current mutation discovery approaches in Retinitis Pigmentosa," *Vision Research*, vol. 75, pp. 117-129, 2012.

[4] L. S. Sullivan, S. J. Bowne, D. G. Birch et al., "Prevalence of disease-causing mutations in families with autosomal dominant retinitis pigmentosa: a screen of known genes in 200 families," *Investigative Ophthalmology and Visual Science*, vol. 47, no. 7, pp. 3052-3064, 2006.

[5] B. J. Seyedahmadi, C. Rivolta, J. A. Keene, E. L. Berson, and T. P. Dryja, "Comprehensive screening of the USH2A gene in Usher syndrome type II and non-syndromic recessive retinitis pigmentosa," *Experimental Eye Research*, vol. 79, no. 2, pp. 167-173, 2004.

[6] A. Ávila-Fernández, D. Cantalapiedra, E. Aller et al., "Mutation analysis of 272 Spanish families affected by autosomal recessive retinitis pigmentosa using a genotyping microarray," *Molecular Vision*, vol. 16, pp. 2550-2558, 2010.

[7] S. El Shamieh, M. Neuille, A. Terray et al., "Whole-exome sequencing identifies KIZ as a ciliary gene associated with autosomal-recessive rod-cone dystrophy," *American Journal of Human Genetics*, vol. 94, no. 4, pp. 625-633, 2014.

[8] I. Audo, S. Mohand-Saïd, C. M. Dhaenens et al., "RP1 and autosomal dominant rod-cone dystrophy: novel mutations, a review of published variants, and genotype-phenotype correlation," *Human Mutation*, vol. 33, no. 1, pp. 73-80, 2012.

[9] X. Guillonneau, N. I. Piriev, M. Danciger et al., "A nonsense mutation in a novel gene is associated with retinitis pigmentosa in a family linked to the RP1 locus," *Human Molecular Genetics*, vol. 8, no. 8, pp. 1541-1546, 1999.

[10] E. A. Pierce, T. Quinn, T. Meehan, T. L. McGee, E. L. Berson, and T. P. Dryja, "Mutations in a gene encoding a new oxygen-regulated photoreceptor protein cause dominant retinitis pigmentosa," *Nature Genetics*, vol. 22, no. 3, pp. 248-254, 1999.

[11] L. S. Sullivan, J. R. Heckenlively, S. J. Bowne et al., "Mutations in a novel retina-specific gene cause autosomal dominant retinitis pigmentosa," *Nature Genetics*, vol. 22, no. 3, pp. 255-259, 1999.

[12] A. Avila-Fernandez, M. Corton, K. M. Nishiguchi et al., "Identification of an RP1 prevalent founder mutation and related phenotype in Spanish patients with early-onset autosomal recessive retinitis," *Ophthalmology*, vol. 119, no. 12, pp. 2616-2621, 2012.

[13] L. J. Chen, T. Y. Y. Lai, P. O. S. Tam et al., "Compound heterozygosity of two novel truncation mutations in RP1 causing autosomal recessive retinitis pigmentosa," *Investigative Ophthalmology and Visual Science*, vol. 51, no. 4, pp. 2236-2242, 2010.

[14] S. Khaliq, A. Abid, M. Ismail et al., "Novel association of RP1 gene mutations with autosomal recessive retinitis pigmentosa," *Journal of Medical Genetics*, vol. 42, no. 5, pp. 436-438, 2005.
[15] A. M. Siemiatkowska, G. D. N. Astuti, K. Arimadyo et al., “Identification of a novel nonsense mutation in RPI that causes autosomal recessive retinitis pigmentosa in an Indonesian family,” *Molecular Vision*, vol. 18, pp. 2411–2419, 2012.

[16] B. Bocquet, N. A. Marzouka, M. Hebrard et al., “Homozygosity mapping in autosomal recessive retinitis pigmentosa families detects novel mutations,” *Molecular Vision*, vol. 19, pp. 2487–2500, 2013.

[17] E. M. Lafont, G. Manes, G. Sénéchal et al., “Patients with retinitis pigmentosa due to RPI mutations show greater severity in recessive than in dominant cases,” *Journal of Clinical & Experimental Ophthalmology*, vol. 2, article 194, 2011.

[18] M. Al-Rashed, L. Abu Safieh, H. Alkuraya et al., “RPI and retinitis pigmentosa: report of novel mutations and insight into mutational mechanism,” *British Journal of Ophthalmology*, vol. 96, no. 7, pp. 1018–1022, 2012.

[19] S. A. Riazuddin, A. Shahzadi, C. Zeitz et al., “A mutation in SLC24A1 implicated in autosomal-recessive congenital stationary night blindness,” *The American Journal of Human Genetics*, vol. 87, no. 4, pp. 523–531, 2010.

[20] H. P. Singh, S. Jalali, R. Narayanan, and C. Kannabiran, “Genetic analysis of indian families with autosomal recessive retinitis pigmentosa by homozygosity screening,” *Investigative Ophthalmology and Visual Science*, vol. 50, no. 9, pp. 4065–4071, 2009.

[21] Q. Liu, A. Lyubarsky, J. H. Skalet, E. N. Pugh Jr., and E. A. Pierce, “RPI is required for the correct stacking of outer segment discs,” *Investigative Ophthalmology & Visual Science*, vol. 44, no. 10, pp. 4171–4183, 2003.

[22] Q. Liu, J. Zuo, and E. A. Pierce, “The retinitis pigmentosa 1 protein is a photoreceptor microtubule-associated protein,” *Journal of Neuroscience*, vol. 24, no. 29, pp. 6427–6436, 2004.

[23] I. Audo, K. M. Bujakowska, T. Léveillard et al., “Development and application of a next-generation-sequencing (NGS) approach to detect known and novel gene defects underlying retinal diseases,” *Orphanet Journal of Rare Diseases*, vol. 7, no. 1, article 8, 2012.

[24] I. Audo, J. A. Sahel, S. Mohand-Saïd et al., “EYS is a major gene for rod-cone dystrophies in France,” *Human Mutation*, vol. 31, no. 5, pp. E1406–E1435, 2010.

[25] T. G. P. Consortium, “A map of human genome variation from population-scale sequencing,” *Nature*, vol. 467, no. 7319, pp. 1061–1073, 2010.

[26] D. M. Altshuler, R. A. Gibbs, L. Peltonen et al., “Integrating common and rare genetic variation in diverse human populations,” *Nature*, vol. 467, no. 7311, pp. 52–58, 2010.

[27] J. A. Tennessen, A. W. Bigham, T. D. O'Connor et al., “Evolution and functional impact of rare coding variation from deep sequencing of human exomes,” *Science*, vol. 337, no. 6090, pp. 64–69, 2012.

[28] M. J. Bamshad, S. B. Ng, A. W. Bigham et al., “Exome sequencing as a tool for Mendelian disease gene discovery,” *Nature Reviews Genetics*, vol. 12, no. 11, pp. 745–755, 2011.

[29] W. J. Kent, C. W. Sugnet, T. S. Furey et al., “The human genome browser at UCSC,” *Genome Research*, vol. 12, no. 6, pp. 996–1006, 2002.

[30] I. A. Adzhubei, S. Schmidt, L. Peshkin et al., “A method and server for predicting damaging missense mutations,” *Nature Methods*, vol. 7, no. 4, pp. 248–249, 2010.

[31] P. Kumar, S. Henikoff, and P. C. Ng, “Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm,” *Nature Protocols*, vol. 4, no. 7, pp. 1073–1082, 2009.

[32] F. O. Desmet, D. Hamroun, M. Lalande, G. Collod-Béroud, M. Claustres, and C. Béroud, “Human splicing finder: an online bioinformatics tool to predict splicing signals,” *Nucleic Acids Research*, vol. 37, no. 9, article no. e67, 2009.

[33] P. D. Stenson, E. V. Ball, M. Mort, A. D. Phillips, K. Shaw, and D. N. Cooper, “The Human Gene Mutation Database (HGMD) and its exploitation in the fields of personalized genomics and molecular evolution,” in *Current Protocols in Bioinformatics*, chapter 1–13, 2012.

[34] I. F. A. C. Fokkema, J. T. Den Dunnen, and P. E. M. Taschner, “LOVD: easy creation of a locus-specific sequence variation database using an “LSDB-in-a-Box” approach,” *Human Mutation*, vol. 26, no. 2, pp. 63–68, 2005.

[35] S. A. Miller, D. D. Dykes, and H. F. Polesky, “A simple salting out procedure for extracting DNA from human nucleated cells,” *Nucleic Acids Research*, vol. 16, no. 3, p. 1215, 1988.

[36] C. Zeitz, B. Kloeckener-Gruissem, U. Forster et al., “Mutations in *CABP4*, the gene encoding the Ca²⁺-binding protein 4, cause autosomal recessive night blindness,” *American Journal of Human Genetics*, vol. 79, no. 4, pp. 657–667, 2006.

[37] S. A. Riazuddin, F. Zulfiqar, Q. Zhang et al., “Autosomal recessive retinitis pigmentosa is associated with mutations in RPI in three consanguineous Pakistani families,” *Investigative Ophthalmology and Visual Science*, vol. 46, no. 7, pp. 2264–2270, 2005.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

Hindawi

Submit your manuscripts at
http://www.hindawi.com

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /></figure>