Variant Profiling of a Large Cohort of 138 Chinese Families With Autosomal Dominant Retinitis Pigmentosa

Retinitis pigmentosa (RP) is the most common form of inherited retinal dystrophy, and 15–25% of RP is transmitted as an autosomal dominant (ad) trait. The objectives of this study were to establish the variant profile in a large cohort of adRP families and to elucidate the variant spectrum of each adRP gene in Chinese patients. A total of 138 probands clinically diagnosed with RP as a presumed autosomal dominant trait were recruited. All probands underwent ophthalmic examinations by specialists. A combination of molecular screening methods, including targeted next-generation sequencing, Sanger DNA sequencing, and multiplex ligation probe amplification assay, was used to detect variants. We identified heterozygous variants of 11 adRP genes in 73 probands, hemizygous, or heterozygous variants of X-linked RP genes in six patients, compound heterozygous variants of autosomal recessive RP genes in three pseudodominant families, and one heterozygous variant of one ad cone and rod dystrophy gene in one proband. One proband was found carrying both variants in RPGR and FAM161A. The overall detection rate was 59.4% (82/138). We detected 72 distinct disease-causing variants involving 16 RP genes and one cone-rod dystrophy gene; 33 of these variants have not been reported previously. Disease-causing variants were identified in the adRP genes in 52.9% of the families, followed by 4.3% in the X-linked RP genes, and 2.2% in the autosomal recessive genes. The most frequent mutant genes were RHO, PRPF31, RP1, SNRNP200, and PRPF8, which explained up to 78.0% of the genetically diagnosed families. Most of the variants identified in adRP genes were missense, and copy number variations were common (7/20) in the PRPF31 gene. We established the profile of the mutated genes and the variant spectrum of adRP genes in a large cohort of Chinese patients, providing essential information for genetic counseling and future development of therapeutics for retinal dystrophy inherited as a dominant trait.

Introduction

Retinitis pigmentosa (RP) is the most common form of inherited retinal dystrophy (IRD), with a prevalence of about 1 in 4,000 (Ayuso and Millan, 2010). RP is a progressive disorder characterized by initial degeneration of rod photoreceptors, followed by degeneration of cone cells (Ayuso and Millan, 2010). Clinical features include night blindness (usually occurring in adolescence), progressive defects of the peripheral visual field, and ultimately, severe damage to central visual acuity (Ayuso and Millan, 2010; Daiger et al., 2014a). The typical fundus appearance consists of black bone-spicule pigmentation in the midperipheral retina, attenuated retinal arterioles, and a pale optic disc. Electroretinograms (ERGs) present reduced or non-recordable signals (Ayuso and Millan, 2010; Daiger et al., 2014a). Most RP cases (about 70–80%) are non-syndromic, which means that patients display only ocular dysfunction. However, some patients may present highly variable clinical symptoms and progression, even if they come from the same family (Ayuso and Millan, 2010).

RP is a highly genetically heterogeneous disorder and can be transmitted in an autosomal recessive (ar), an autosomal dominant (ad), or an X-linked recessive pattern (xl) (Ayuso and Millan, 2010; Daiger et al., 2014a). Digenic inheritance has been reported in some rare cases (Ayuso and Millan, 2010). Autosomal dominant RP (adRP) accounts for about 15–25% of the total RP cases (Daiger et al., 2014a). At present, 30 causative genes have been identified for adRP (RetNet: https://sph.uth.edu/retnet/home.htm), and more than 1,000 different kinds of variants have been described in those genes (Daiger et al., 2014a; Dias et al., 2018). The most commonly mutated gene is RHO, which is responsible for ~20–30% of adRP cases (Daiger et al., 2014a; Dias et al., 2018). Most of the detected variants in the adRP genes are private variants, and a small fraction of the common variants are ethnicity specific; for example, the RHO variant p.(P23H) is a founder variant almost exclusively described in Americans of European origin (Sullivan et al., 2013; Daiger et al., 2014a,b; Dias et al., 2018). Copy number variants (CNVs) have also been reported in some adRP genes, and incomplete penetrance has been observed in some specific genes (Sullivan et al., 2013; Daiger et al., 2014a,b; Dias et al., 2018). Taken together, these findings for adRP increase the complexity of its genetic diagnosis, which is very crucial for genetic consulting and gene therapy in patients with this disorder.

In recent years, the application of next-generation sequencing (NGS), mostly as targeted exome sequencing (TES) and whole exome sequencing (WES), has greatly increased the genetic diagnosis rates of different form of IRD (Daiger et al., 2014b; Xu et al., 2014; Costa et al., 2017; Van Cauwenbergh et al., 2017; Martin-Merida et al., 2018; Gao et al., 2019). The profile of mutated genes and the variant spectrum of adRP genes have been established in European and American patients with variant detection rates between 50 and 75% (Daiger et al., 2014b; Costa et al., 2017; Van Cauwenbergh et al., 2017; Martin-Merida et al., 2018). One very recent study that reported genetic analysis using TES in a large cohort of Chinese patients has shown that 72% of the patients could receive a molecular diagnosis (Gao et al., 2019). However, the large cohort included several other IRD cases, such as Bietti crystalline retinopathy, Leber congenital amaurosis, and retinitis punctata albescens, and it did not report the gene profiles for adRP patients (Gao et al., 2019). Several previous studies have reported the gene profiles in the relatively small Chinese adRP cohort (not more than 78 families) (Li et al., 2010; Xu et al., 2014; Huang et al., 2017).

In the current study, we have described the outcomes of a comprehensive molecular analysis of 138 pedigrees with possible adRP by a combination of methods, including TES, Sanger sequencing, and real-time quantitative polymerase chain reaction (q-PCR) analysis or multiplex ligation-dependent probe amplifications (MLPAs).

Subjects and Methods

Patients

In total, 138 unrelated families with a clinical diagnosis of RP were enrolled at the Genetics Laboratory of the Beijing Institute of Ophthalmology, Beijing Tongren Ophthalmic Center. Dominant inheritance was presumed when each family had affected members in at least two consecutive generations. This cohort included four previously reported families (Pan et al., 2012; Dong et al., 2013). Patients were diagnosed with RP based on the following criteria: a history of night blindness, progressive visual field defects, fundus displaying bone spicule-like pigment clumping in the midperipheral or peripheral retina and attenuation of retinal vessels, and severe rod-cone dysfunction or non-recordable ERG recording (Ayuso and Millan, 2010). The molecular testing processes were prospectively evaluated and approved by the ethics committee of Beijing Tongren Hospital, and all tests were implemented under the official guidelines of the Beijing Tongren Hospital Joint Committee on Clinical Investigation in compliance with the Declaration of Helsinki. Informed consent was obtained from probands after a detailed explanation of the processes.

Each proband and available family members underwent a regular ophthalmic examination that included best-corrected visual acuity, slit-lamp biomicroscopy, and a fundus examination. Most participants also underwent optical coherence tomography, visual field, and ERG evaluations. Peripheral blood samples were collected from patients and their family members, and genomic DNA was extracted from the leukocytes with a genomic DNA extraction and purification kit (vigorous whole-blood genomic DNA extraction kit; Vigorous Beijing, China), according to the producer's protocol.

PCR-Based Sequencing of the RHO Gene

Five exons and flanking splicing sites of the RHO gene were first sequenced for all probands diagnosed with adRP, except for the four families that we have previously reported (Pan et al., 2012; Dong et al., 2013). Those four previously reported families were analyzed by linkage mapping following Sanger sequencing (Pan et al., 2012; Dong et al., 2013). The PCR amplifications were performed with regular reaction mixtures, and the purified amplicons were sequenced on an ABI Prism 373A DNA sequencer (Applied Biosystems, Foster City, CA, USA). The sequencing outcomes were matched to the available cDNA sequence of RHO (GenBank NM_000539).

TES and Bioinformatics Analysis

We performed TES in 113 patients who did not carry any RHO variants using a capture panel developed and evaluated by our group (Sun et al., 2018). This panel comprised 188 known IRD genes, and 26 of them were adRP genes. Details of the procedures, which included the Illumina library preparation, capture experiments, and the enrichment libraries sequencing, have been described previously (Sun et al., 2018). The raw sequencing data processing, calling, and evaluation were carried out as previously reported (Sun et al., 2018). The bioinformatics programs PolyPhen2 (http://genetics.bwh.harvard.edu/pph/), Mutation Taster (http://www.mutationtaster.org/), and SIFT (http://sift.jcvi.org/) were used to predict the pathogenicity of each variant. The programs NetGene2 Server (http://www.cbs.dtu.dk/services/NetGene2/), Human Splice Finder (http://www.umd.be/HSF3/), and Berkeley Drosophila Genome Project (http://www.fruitfly.org/seq_tools/splice.html) were used to analyze any variants involving a splicing effect. We further determined the pathogenicity of each variant by searching the reported pathogenic variants in the HGMD database (http://www.hgmd.cf.ac.uk/ac/index.php) and the LOVD database (https://www.lovd.nl/). We ultimately classified the variants into pathogenic or likely pathogenic variants, variants of uncertain significance (VUS), and benign or likely benign variants according to the standards described by the American College of Medical Genetics and Genomics (ACMG) (Richards et al., 2015). Sanger sequencing was conducted to verify the supposed disease-causing variants and VUS. Segregation analysis was done for the probands and their family members.

CNV Analysis and Validation

We employed the CNV kit software (https://github.com/etal/cnvkit) to identify CNVs from variations in the read depth for the patients who had TES data (Talevich et al., 2016). Real-time q-PCR was then performed to confirm the presence of the presumed CNVs of PRPF31 in five families and of FAM161A in one family, as we previously reported (Dong et al., 2013). MLPA assays were conducted for three probands and their family members using the SALSA MLPA Kit P235 (Amsterdam, the Netherlands), following the producer's protocols.

Supplementary PCR-Based Sequencing

We performed Sanger sequencing of exon 15 [open reading frame 15 (ORF15)] of RPGR in all male patients whose disease-causing variants were not found after TES and whose pedigrees could not exclude X-linked transmission.

Statistical Analysis

All statistical analyses were performed using SPSS statistical software (version 25.0, SPSS Inc., Chicago, IL). The Kolmogorov–Smirnov test was used to evaluate whether the onset age of a single group conformed to a normal distribution. The Wilcoxon rank sum test was used to evaluate any difference in the onset age between two groups. P ≤ 0.05 was considered statistically significant.

Results

Variant Detection Rate and Variant Spectrum

We identified heterozygous autosomal dominant variants that were pathogenic or likely pathogenic in 74 probands, hemizygous, or heterozygous in six probands, and compound heterozygous in three patients, for a general variant detection rate of 59.4% (82/138) (Table 1, Figure 1A). We found a total of 72 different variants in 17 genes, including 11 adRP genes, two xlRP genes, three arRP genes, and one autosomal dominant cone and rod dystrophy (adCORD) gene (Figure 1). Four of the 72 variants were identified three times or more, and the remaining 68 variants were detected either once (88.9%) or twice (5.6%). The most common variants were p.(R135W) and p.(P347L) in RHO, with a gene-specific allele frequency of 17.4% (4/23), followed by p.(R677^*) in the RP1 (37.5%, 3/8) and a whole PRPF31 deletion (15.0%, 3/20) (Supplementary Table 1).

TABLE 1

Table 1. Summary of Family Composition and Variant Screening Results in 138 Chinese adRP Families.

FIGURE 1

Figure 1. Summary of the proportion of patients in this study with variants in the involved genes. (A) Proportion of all probands with variants in each involved gene. (B) Proportion in all genetically diagnosed families with variants in each involved gene. The shaded areas indicate a proband carrying both variants of RPGR and FAM161A.

Of the 72 putative disease-causing variants, 33 variants were first identified in the current study (Supplementary Table 1). These 33 novel variants comprised 16 missense, five frameshift indel, four non-sense, three splicing effect, three CNV, one synonymous, and one stop-codon-lost variant. None of these novel variants were recorded in our in-house and any public databases, such as the Exome Variant Server and 1000 Genomes Database (Supplementary Table 1). The 16 novel missense variants were defined as pathogenic or likely pathogenic variants according to the ACMG guidelines and standards (Supplementary Table 1). One synonymous variant c.1146G>A, p.(E382E), located in the last base of the exon 11 of PRPF31, was predicted to alter the downstream splice effect by NetGene2, HSF, and Mutation Taster. The remaining variants (frameshift small indel, non-sense, splicing effect, and run-on variants) or CNVs were considered to be obviously pathogenic variants. Five variants of PRPF8 were identified in this cohort, and they all were novel. These five variants contained four missense variants and one frameshift small deletion that escaped non-sense-mediated decay. Three of these variants were located in the C-terminal Jab1 domain, and the other two were in the Linker and RNAseH-like domain (Figure 2A). We also detected an unreported missense variant, p.(S216G), which is not located in any domains of PRPF8 (Figure 2A). As p.(S216G) was predicted to be benign by both Polyphen2 and SIFT, we defined this variant as a VUS according to the ACMG guidelines and standards.

FIGURE 2

Figure 2. The distribution of 30 distinct variants of PRPF31, SNRNP200, PRPF8, and PRPF6 identified in our study. (A) The distribution of 21 missense variants in the corresponding domain of PRPF31, SNRNP200, PRPF8, and PRPF6 and the distribution of five splicing effect variants in the exons and introns of PRPF31 and SNRNP200. Numbers under domains indicate amino acid location. Red asterisks indicate missense variants detected in the current cohort. Blue asterisks indicate missense variants reported in HGMD. Yellow diamond indicates the VUS of PRPF8 detected in the current cohort. Red triangles indicate splicing effect variants detected in the current cohort. Blue triangles indicate splicing effect variants reported in HGMD. Prp1, pre-mRNA processing domain; TPR, tetratricopeptide repeat domain; NT, PRO8NT domain; RT, reverse transcriptase homology domain; Endo, restriction endonuclease homology domain; RH, RNase H homology domain; Jab1, Jun kinase activation binding protein; NOP, nucleolar protein domain; RecA1, repeat helicase ATP-binding domain; RecA2, repeat helicase C-terminal domain; Sec63, secretory-63 domain. (B) Lengths and positions of the four gross deletions involving PRPF31 on chromosome 19q13.42.

Variant Profile in the adRP Genes and Related Clinical Features

The 59 putative disease-causing variants of the 11 adRP genes were identified in 73 pedigrees (52.9%, 73/138) and accounted for 89.0% (73/82) of the genetically diagnosed cases in this cohort (Figure 1B). Cosegregation analyses were performed in 65 pedigrees of the 73 families (Supplementary Table 2). The most frequently mutated gene was RHO, identified in 23 of the 138 probands (16.7%), followed by PRPF31 in 20 unrelated patients (14.5%), RP1 and SNRNP200 each in 8 patients (5.8%), and PRPF8 in 5 patients (3.6%) (Figure 1A). Variants in the pre-mRNA splicing factor genes accounted for 24.6% of the families and the variant locations in each gene are displayed in Figure 2A. Most of the variants in RHO, SNRNP200, PRPF8, and IMPDH1 were missense variants, whereas loss-of-function variants, including non-sense, frameshift small indels, and CNVs, were more frequently observed in PRPF31, RP1, CRX, and TOPORS (Table 2). In the current cohort, CNVs were more prevalent in PRPF31, with a gene-specific frequency of 35.0% (7/20) (Figure 2B).

TABLE 2

Table 2. Prevalence of Disease-Causing Variants in 138 Chinese adRP families and mutation type for each gene.

All probands with the adRP gene variants had typical RP symptoms, and their representative fundus appearances of the patients with variants in the first five common genes are shown in Figure 3. The median onset age of all patients with adRP genes was 5 years (range, 1–45 years) (Table 3). The median onset age was significantly older for the probands carrying the variants of RP1 than for all probands with adRP (Table 3). Incomplete penetrance was observed in five pedigrees—three carrying variants in PRPF31, one in PRPF8, and one in TOPORS.

FIGURE 3

Figure 3. Colored fundus (CF) photographs of patients with variants of RHO, PRPF31, SNRNP200, RP1, and PRPF8. (A) CF image of proband 019875 with the missense variant of RHO. (B) CF photograph of proband 0191443 with exons 2–3 deletion of PRPF31. (C1,C2) CF photograph of proband 019417 with the non-sense variant of PRRF31 shows dense black pigmentation in the midperipheral retina, but CF image of his 45-year-old father with the same variant displaying a normal fundus appearance (D1,D2) CF photographs of two affected sisters of pedigree 019917 with the missense variants of PRPF8 located in the Jab1 domain. (E,F) CF photographs of two probands of 010215 and 019444 with the missense variants of PRPF8 located in the Linker and RH region or domain. (G,H) CF photographs of probands 019545 and 0191212 with the variants of RP1. (I1,I2) CF images of two patients from pedigree 019756 with the missense variant of SNRNP200.

TABLE 3

Table 3. Correlations between onset age of patients with different gene variants.

Variant Profiles in the xlRP and arRP Genes

Six hemizygous or heterozygous variants of the two xlRP genes were identified in six probands (4.3%, 6/138) and three compound heterozygous variants of three arRP genes in three probands (2.2%) (Figure 4). All six families with xlRP gene variants had at least two generations of inheritance; however, none of them showed male-to-male transmission that would exclude X-linked inheritance (Figure 4). Proband 0191318 was found to carry one heterozygous variant, p.(Q565Rfs^*17), of RPGR and a compound heterozygous variant, M1:p.(Q548^*) and M2:exon 1-5del, of FAM161A. Cosegregation analysis showed that his affected father carried a hemizygous variant of RPGR and a heterozygous variant, p.(Q548^*), of FAMI161A. Pedigree 019207 was a small two-generation family with male-to-male transmission; however, the proband was found to harbor a compound heterozygous variant, p.[(G268R)]; [(P2811T)], of the USH2A gene. Cosegregation analysis showed that his affected father only carried a heterozygous variant, p.(G268R). We then conducted TES analysis on his affected father, but we did not detect any other disease-causing variant. Therefore, these two families were considered to be pseudodominant pedigrees (Figure 4). Proband 019645 was found to carry a compound heterozygous variant of EYS. During his family history review, he stated that his father and sister also had night blindness and visual acuity defects, but his father did not suffer from RP, according to a later ophthalmologic examination.

FIGURE 4

Figure 4. Pedigrees of eight families with xlRP and arRP gene variants and their cosegregation results.

Variants in Other IRD Genes

Proband 0191323 was found to carry a heterozygous variant, p.(S679T), of RIMS1, which is a disease-causing gene for adCORD. Cosegregation analysis showed that his affected mother also harbored this variant. The 51-year-old proband stated that he had experienced night-blindness since the age of 8 years and had developed obvious visual acuity loss and photophobia at around age 40 years. His fundus examination showed macular atrophy and black spicule-like pigment in the posterior pole and peripheral retina, and his ERG recording was extinguished. This patient might be in the late stage of CORD, so his fundus examination presented an RP-like appearance.

Discussion

In the current study, we conducted a comprehensive molecular analysis in a large Chinese adRP cohort. We obtained an overall variant detection rate of 59.4% by means of several molecular methods, including linkage mapping, Sanger sequencing, targeted-exon sequencing, MLPA, and real-time q-PCR analysis. This detection rate is close to the rates recently reported in a large Spanish adRP cohort comprising 258 families (60%), as well as in a previously described small Belgian adRP cohort comprising 86 families (56%) (Van Cauwenbergh et al., 2017; Martin-Merida et al., 2018). However, it is much lower than the rate (70%) described previously in a large American adRP cohort that included 253 families (Daiger et al., 2014b). The variant detection rate is related to the accuracy of the probands' clinical diagnoses. In the current study, the solving rate was much higher for the families with male-to-male transmission (69.3%) than for the families without male-to-male transmission (51.3%). The variant detection rate is also related to molecular screening methods. For the large American adRP cohort, application of next-generation sequencing increased the detection rate from 65% (by only Sanger sequencing) to 70% (Daiger et al., 2014b).

Consistent with several previous studies, RHO was the most frequently mutated gene in the current adRP cohort, and its proportion (16.7%) was higher than the proportions reported in Chinese (8.9%) and Japanese (11.5%) adRP cohorts (Xu et al., 2014; Koyanagi et al., 2019) and close to those reported in French (18.3%), Italian (16.0%), and Belgian (14.0%) adRP cohorts (Ziviello et al., 2005; Coussa et al., 2015; Van Cauwenbergh et al., 2017), but much lower than those observed in American (26.9%) adRP cohorts (Daiger et al., 2014a,b). This difference might reflect, in part, the absence of the P23H variant, which is a founder variant and accounts for about 13.2% of the adRP in American families (Daiger et al., 2014a). The second most common mutated gene was PRPF31, which showed a 14.5% variant frequency in the current cohort. This frequency was higher than that observed in several adRP cohorts (6.8–10.5%) from Japan, Belgium, America, and Spain (Daiger et al., 2014a; Van Cauwenbergh et al., 2017; Martin-Merida et al., 2018; Koyanagi et al., 2019). The variant frequency of SNRNP200 was also much higher in the current cohort than in the Spanish and American cohorts (5.8% vs. 2.3 or 1.5%) (Daiger et al., 2014a; Martin-Merida et al., 2018). Therefore, the findings that the variants in the pre-mRNA splicing factor genes were the major cause of the current Chinese adRP cohort and were responsible for 24.6% of the probands were not surprising (Figure 1B). By contrast, the variant frequency was much lower for PRPH2 (0.7%) in our cohort than in the American, Belgian, and Spanish cohorts (7.0, 4.7, and 3.9%, respectively) (Daiger et al., 2014a; Van Cauwenbergh et al., 2017; Martin-Merida et al., 2018). The rate was also lower than the rate previously reported (3.7%, 3/79) in a small Chinese adRP cohort (Xu et al., 2014). Several previous studies have indicated that the variant frequency of PRPH2 in adRP varies extensively with ethnicity, ranging from 0% (Mexican cohort) up to 10.3% (French cohort) (Sullivan et al., 2006; Matias-Florentino et al., 2009; Manes et al., 2015).

Our cohort showed a greater allelic diversity, and the majority of the observed variants (almost 90%) were private variants responsible for their own respective families. Only four variants in RHO, PRPF31, and RP1 were detected three or more times. The most frequent variants, p.(R135W) and p.(P347L) in RHO, were variant hotspots that have been reported in several previous studies (Xu et al., 2014; Van Cauwenbergh et al., 2017; Martin-Merida et al., 2018). The common variant p.(R677^*) in RP1 was also a variant hotspot that has been described in many adRP families with different ethnicities (Bowne et al., 1999; Jacobson et al., 2000; Audo et al., 2012; Daiger et al., 2014a,b; Martin-Merida et al., 2018).

The proportion of novel variants identified in the current cohort was relatively high. In the current study, two novel missense variants, p.(R1384W) and p.(T1931M) of PRPF8, were not located in the C-terminal Jab1 domain (Figure 3A), where almost all reported variants are clustered (RuŽičková and Staněk, 2017). These two variants were situated in highly conserved regions and were predicted to be disease-causing by three in silico analysis programs and to cosegregate with their phenotype in the families. They were not found in any public databases, and we defined them as disease-causing variants according to the ACMG guidelines and standards. Additional functional studies are needed to verify their pathogenicity in the future.

Five distinct and large genomic DNA deletions of two genes were identified in eight probands. PRPF31 showed a high prevalence (5.1%, 7/138) of CNVs in our adRP cohort; this prevalence was much higher than the prevalence of 1.9% recently reported in a large Spanish cohort (Martin-Merida et al., 2018). All five CNVs were in a heterozygous state; therefore, they were undetected by Sanger sequencing. Next-generation sequencing has the capability to detect CNVs, but this capability is related to the coverage depth. A previous study found that CNV analysis could generate ambiguous results in the target regions when the coverage was <250 × (Aparisi et al., 2014). In the current cohort, we used comprehensive screening methods, including linkage analysis in the early period, q-PCR, MLPA, and TES read count analysis, to detect the CNVs of PRPF31. This comprehensive analysis might be one of the reasons for the high CNV frequency for PRPF31 observed in the present study.

Six probands (4.3%, 6/138) in the current cohort were found to harbor variants of the xlRP genes (RPGR and RP2). Four of the six probands were female and did not present a milder phenotype than was observed for the affected individuals in their families. Variant p.(E1031fs^*58) (c.3092delA) located in the ORF15 of RPGR was identified in the supplementary PCR-based sequencing. The ORF15 where most RPGR variants are clustered is poorly covered in the NGS because of its highly repetitive sequence (Huang et al., 2015). Some variants in the ORF15 might be missed by NGS assays; therefore, alternative Sanger sequencing is necessary in patients from any adRP families without male-to-male transmission. A previous study using Sanger sequencing identified variants of RPGR and RP2 in 8.5% (22/258) of adRP families (Churchill et al., 2013).

Obtaining a clear genotype–phenotype correlation is difficult because most probands harbor their own private variants in different genes. In the current cohort, we observed that patients with RP1 variants had a late onset age and a mild phenotype when compared with patients with other adRP genes variants. This observation was consistent with previous descriptions (Jacobson et al., 2000; Audo et al., 2012). Incomplete penetrance was relatively common in the families with variants of PRPF31; however, its prevalence (15%, 3/20) was much lower than the 66.7% observed in Spanish patients (Martin-Merida et al., 2017). Three of our probands carried biallelic variants in three known arRP genes. Cosegregation analysis revealed pseudodominant inheritance in two families, with the remaining family unclear because of an incorrect family history report; this is one of the limitations—we did not perform ophthalmic evaluation in all affected family members of the current study.

Of the 56 probands whose disease-causing genes were not identified, 37 pedigrees did not show male-to-male transmission, and three families presented incomplete penetrance. In addition, 21 of the probands did not undergo ERG recording, which is another limitation of the current study as it means that the probability of a misdiagnosis could not be excluded. Their causal variants might be in areas not covered by our TES panel, such as the promoter or deep intronic regions of the targeted genes, or they may reside in newly identified genes that were not included in our TES panel. In the future, we will perform WES or WGS analyses for these patients.

In conclusion, we have established the profile of the mutated genes and the variant spectrum of adRP genes in a large cohort of Chinese patients, thereby providing essential information for genetic counseling and future therapeutic development for retinal dystrophy inherited as a dominant trait.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics committee of Beijing Tongren Hospital. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

YL and Z-BJ designed and supervised the whole study. TX, YX, XinZ, KX, and XiaZ performed the experiments and interpret the results. TX wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Key R&D Program of China, 2016YFC0905200 and the High-level Talents training plan of the health system of Beijing (Grant No. 2013-2-021). The funding organization had no role in designing or conducting this research.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2020.629994/full#supplementary-material

References

Aparisi, M. J., Aller, E., Fuster-García, C., García-García, G., Rodrigo, R., Vázquez-Manrique, R. P., et al. (2014). Targeted next generation sequencing for molecular diagnosis of usher syndrome. Orphanet J. Rare Dis. 9:168. doi: 10.1186/s13023-014-0168-7

CrossRef Full Text | Google Scholar

Audo, I., Mohand-Saïd, S., Dhaenens, C. M., Germain, A., Orhan, E., Antonio, A., et al. (2012). RP1 and autosomal dominant rod-cone dystrophy: novel variants, a review of published variants, and genotype-phenotype correlation. Hum. Mutat. 33, 73–80. doi: 10.1002/humu.21640

CrossRef Full Text | Google Scholar

Ayuso, C., and Millan, J. M. (2010). Retinitis pigmentosa and allied conditions today: a paradigm of translational research. Genome Med 2:34. doi: 10.1186/gm155

CrossRef Full Text | Google Scholar

Bowne, S. J., Daiger, S. P., Hims, M. M., Sohocki, M. M., Malone, K. A., McKie, A. B., et al. (1999). Variants in the RP1 gene causing autosomal dominant retinitis pigmentosa. Hum. Mol. Genet. 8, 2121–2128. doi: 10.1093/hmg/8.11.2121

PubMed Abstract | CrossRef Full Text | Google Scholar

Churchill, J. D., Bowne, S. J., Sullivan, L. S., Lewis, R. A., Wheaton, D. K., Birch, D. G., et al. (2013). Variants in the X-linked retinitis pigmentosa genes RPGR and RP2 found in 8.5% of families with a provisional diagnosis of autosomal dominant retinitis pigmentosa. Invest. Ophthalmol. Vis. Sci. 54, 1411–1416. doi: 10.1167/iovs.12-11541

CrossRef Full Text | Google Scholar

Costa, K. A., Salles, M. V., Whitebirch, C., Chiang, J., and Sallum, J. M. F. (2017). Gene panel sequencing in Brazilian patients with retinitis pigmentosa. Int. J. Retina Vitreous 3:33. doi: 10.1186/s40942-017-0087-6

CrossRef Full Text | Google Scholar

Coussa, R. G., Chakarova, C., Ajlan, R., Taha, M., Kavalec, C., Gomolin, J., et al. (2015). Genotype and phenotype studies in autosomal dominant retinitis pigmentosa (adRP) of the French Canadian founder population. Invest. Ophthalmol. Vis. Sci. 56, 8297–8305. doi: 10.1167/iovs.15-17104

CrossRef Full Text | Google Scholar

Daiger, S. P., Bowne, S. J., and Sullivan, L. S. (2014a). Genes and variants causing autosomal dominant retinitis pigmentosa. Cold Spring Harb. Perspect. Med. 5:a017129. doi: 10.1101/cshperspect.a017129

PubMed Abstract | CrossRef Full Text | Google Scholar

Daiger, S. P., Bowne, S. J., Sullivan, L. S., Blanton, S. H., Weinstock, G. M., Koboldt, D. C., et al. (2014b). Application of next-generation sequencing to identify genes and variants causing autosomal dominant retinitis pigmentosa (adRP). Adv. Exp. Med. Biol. 801, 123–129. doi: 10.1007/978-1-4614-3209-8_16

CrossRef Full Text | Google Scholar

Dias, M. F., Joo, K., Kemp, J. A., Fialho, S. L., da Silva Cunha, A. Jr., Woo, S. J., et al. (2018). Molecular genetics and emerging therapies for retinitis pigmentosa: basic research and clinical perspectives. Prog. Retin. Eye Res. 63, 107–131. doi: 10.1016/j.preteyeres.2017.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, B., Chen, J., Zhang, X., Pan, Z., Bai, F., and Li, Y. (2013). Two novel PRP31 premessenger ribonucleic acid processing factor 31 homolog variants including a complex insertion-deletion identified in Chinese families with retinitis pigmentosa. Mol. Vis. 19, 2426–2435.

Google Scholar

Gao, F. J., Li, J. K., Chen, H., Hu, F. Y., Zhang, S. H., Qi, Y. H., et al. (2019). Genetic and clinical findings in a large cohort of Chinese patients with suspected retinitis pigmentosa. Ophthalmology 126, 1549–1556. doi: 10.1016/j.ophtha.2019.04.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, L., Zhang, Q., Huang, X., Qu, C., Ma, S., Mao, Y., et al. (2017). Variant screening in genes known to be responsible for retinitis pigmentosa in 98 small han Chinese families. Sci. Rep. 7:1948. doi: 10.1038/s41598-017-00963-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X. F., Wu, J., Lv, J. N., Zhang, X., and Jin, Z. B. (2015). Identification of false-negative variants missed by next-generation sequencing in retinitis pigmentosa patients: a complementary approach to clinical genetic diagnostic testing. Genet. Med. 17, 307–311. doi: 10.1038/gim.2014.193

CrossRef Full Text | Google Scholar

Jacobson, S. G., Cideciyan, A. V., Iannaccone, A., Weleber, R. G., Fishman, G. A., Maguire, A. M., et al. (2000). Disease expression of RP1 variants causing autosomal dominant retinitis pigmentosa. Invest. Ophthalmol. Vis. Sci. 41, 1898–1908.

Google Scholar

Koyanagi, Y., Akiyama, M., Nishiguchi, K. M., Momozawa, Y., Kamatani, Y., Takata, S., et al. (2019). Genetic characteristics of retinitis pigmentosa in 1204 Japanese patients. J. Med. Genet. 56, 662–670. doi: 10.1136/jmedgenet-2018-105691

CrossRef Full Text | Google Scholar

Li, S., Xiao, X., Wang, P., Guo, X., and Zhang, Q. (2010). Variant spectrum and frequency of the RHO gene in 248 Chinese families with retinitis pigmentosa. Biochem. Biophys. Res. Commun. 401, 42–47. doi: 10.1016/j.bbrc.2010.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Manes, G., Guillaumie, T., Vos, W. L., Devos, A., Audo, I., Zeitz, C., et al. (2015). High prevalence of PRPH2 in autosomal dominant retinitis pigmentosa in france and characterization of biochemical and clinical features. Am. J. Ophthalmol. 159, 302–314. doi: 10.1016/j.ajo.2014.10.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin-Merida, I., Aguilera-Garcia, D., Fernandez-San Jose, P., Blanco-Kelly, F., Zurita, O., Almoguera, B., et al. (2018). Toward the variantal landscape of autosomal dominant retinitis pigmentosa: a comprehensive analysis of 258 Spanish families. Invest. Ophthalmol. Vis. Sci.59, 2345–2354. doi: 10.1167/iovs.18-23854

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin-Merida, I., Sanchez-Alcudia, R., Fernandez-San Jose, P., Blanco-Kelly, F., Perez-Carro, R., Rodriguez-Jacy da Silva, L., et al. (2017). Analysis of the PRPF31 gene in Spanish autosomal dominant retinitis pigmentosa patients: a novel genomic rearrangement. Invest. Ophthalmol. Vis. Sci. 58, 1045–1053. doi: 10.1167/iovs.16-20515

PubMed Abstract | CrossRef Full Text | Google Scholar

Matias-Florentino, M., Ayala-Ramirez, R., Graue-Wiechers, F., and Zenteno, J. C. (2009). Molecular screening of rhodopsin and peripherin/RDS genes in Mexican families with autosomal dominant retinitis pigmentosa. Curr. Eye Res. 34, 1050–1056. doi: 10.3109/02713680903283169

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, Z., Lu, T., Zhang, X., Dai, H., Yan, W., Bai, F., et al. (2012). Identification of two variants of the RHO gene in two Chinese families with retinitis pigmentosa: correlation between genotype and phenotype. Mol. Vis. 18, 3013–3020.

PubMed Abstract | Google Scholar

Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American college of medical genetics and genomics and the association for molecular pathology. Genet. Med. 17, 405–424. doi: 10.1038/gim.2015.30

PubMed Abstract | CrossRef Full Text | Google Scholar

RuŽičková, Š., and Staněk, D. (2017). Variants in spliceosomal proteins and retina degeneration. RNA Biol. 14, 544–552. doi: 10.1080/15476286.2016.1191735

PubMed Abstract | CrossRef Full Text | Google Scholar

Sullivan, L. S., Bowne, S. J., Birch, D. G., Hughbanks-Wheaton, D., Heckenlively, J. R., Lewis, R. A., et al. (2006). Prevalence of disease-causing variants in families with autosomal dominant retinitis pigmentosa: a screen of known genes in 200 families. Invest. Ophthalmol. Vis. Sci. 47, 3052–3064. doi: 10.1167/iovs.05-1443

PubMed Abstract | CrossRef Full Text | Google Scholar

Sullivan, L. S., Bowne, S. J., Reeves, M. J., Blain, D., Goetz, K., Ndifor, V., et al. (2013). Prevalence of variants in eyeGENE probands with a diagnosis of autosomal dominant retinitis pigmentosa. Invest. Ophthalmol. Vis. Sci. 54, 6255–6261. doi: 10.1167/iovs.13-12605

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, T., Xu, K., Ren, Y., Xie, Y., Zhang, X., Tian, L., et al. (2018). Comprehensive molecular screening in Chinese usher syndrome patients. Invest. Ophthalmol. Vis. Sci. 59, 1229–1237. doi: 10.1167/iovs.17-23312

PubMed Abstract | CrossRef Full Text | Google Scholar

Talevich, E., Shain, A. H., Botton, T., and Bastian, B. C. (2016). CNVkit: genome-wide copy number detection and visualization from targeted DNA sequencing. PLoS Comput. Biol. 12:e1004873. doi: 10.1371/journal.pcbi.1004873

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Cauwenbergh, C., Coppieters, F., Roels, D., De Jaegere, S., Flipts, H., De Zaeytijd, J., et al. (2017). Variants in splicing factor genes are a major cause of autosomal dominant retinitis pigmentosa in belgian families. PLoS ONE 12:e0170038. doi: 10.1371/journal.pone.0170038

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Guan, L., Shen, T., Zhang, J., Xiao, X., Jiang, H., et al. (2014). Variants of 60 known causative genes in 157 families with retinitis pigmentosa based on exome sequencing. Hum. Genet. 133, 1255–1271. doi: 10.1007/s00439-014-1460-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ziviello, C., Simonelli, F., Testa, F., Anastasi, M., Marzoli, S. B., Falsini, B., et al. (2005). Molecular genetics of autosomal dominant retinitis pigmentosa (ADRP): a comprehensive study of 43 Italian families. J. Med. Genet. 42:e47. doi: 10.1136/jmg.2005.031682

PubMed Abstract | CrossRef Full Text | Google Scholar