Next Article in Journal
Capsaicin-Sensitive Peptidergic Sensory Nerves Are Anti-Inflammatory Gatekeepers in the Hyperacute Phase of a Mouse Rheumatoid Arthritis Model
Next Article in Special Issue
Paternal Uniparental Isodisomy of Chromosome 2 in a Patient with CNGA3-Associated Autosomal Recessive Achromatopsia
Previous Article in Journal
Diabetes, Antidiabetic Medications and Cancer Risk in Type 2 Diabetes: Focus on SGLT-2 Inhibitors
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 4
10.3390/ijms22041681
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clinical and Molecular Characterization of Achromatopsia Patients: A Longitudinal Study

by
Raffaella Brunetti-Pierri
1,†,
Marianthi Karali
1,2,†,
Paolo Melillo
1,
Valentina Di Iorio
1,
Antonella De Benedictis
1,
Gennarfrancesco Iaccarino
1,
Francesco Testa
1,
Sandro Banfi
2,3 and
Francesca Simonelli
1,*
1
Eye Clinic, Multidisciplinary Department of Medical, Surgical and Dental Sciences, Università degli Studi della Campania “Luigi Vanvitelli”, via Pansini 5, 80131 Naples, Italy
2
Telethon Institute of Genetics and Medicine, via Campi Flegrei 34, 80078 Pozzuoli, Italy
3
Medical Genetics, Department of Precision Medicine, Università degli Studi della Campania “Luigi Vanvitelli”, via Luigi De Crecchio 7, 80138 Naples, Italy
*
Author to whom correspondence should be addressed.
†
Contributed equally to the work.
Int. J. Mol. Sci. 2021, 22(4), 1681; https://doi.org/10.3390/ijms22041681
Submission received: 30 December 2020 / Revised: 28 January 2021 / Accepted: 3 February 2021 / Published: 7 February 2021
(This article belongs to the Special Issue Achromatopsia: From Genetics to Therapy)
Browse Figures
Versions Notes

Abstract

:
Achromatopsia (ACHM) is a rare genetic disorder of infantile onset affecting cone photoreceptors. To determine the extent of progressive retinal changes in achromatopsia, we performed a detailed longitudinal phenotyping and genetic characterization of an Italian cohort comprising 21 ACHM patients (17 unrelated families). Molecular genetic testing identified biallelic pathogenic mutations in known ACHM genes, including four novel variants. At baseline, the patients presented a reduced best corrected visual acuity (BCVA), reduced macular sensitivity (MS), normal dark-adapted electroretinogram (ERG) responses and undetectable or severely reduced light-adapted ERG. The longitudinal analysis of 16 patients (mean follow-up: 5.4 ± 1.0 years) showed a significant decline of BCVA (0.012 logMAR/year) and MS (−0.16 dB/year). Light-adapted and flicker ERG responses decreased below noise level in three and two patients, respectively. Only two patients (12.5%) progressed to a worst OCT grading during the follow-up. Our findings corroborate the notion that ACHM is a progressive disease in terms of BCVA, MS and ERG responses, and affects slowly the structural integrity of the retina. These observations can serve towards the development of guidelines for patient selection and intervention timing in forthcoming gene replacement therapies.

1. Introduction

Complete achromatopsia (ACHM) (also termed total color blindness or rod monochromacy) is a rare inherited autosomal recessive disease with an estimated prevalence of approximately 1 in 30,000–50,000 worldwide [1]. It has been traditionally classified as a cone dysfunction syndrome rather than a cone dystrophy [2] to reflect its predominantly stationary nature. Nevertheless, there is also growing evidence on the occurrence of progressive morphological and functional retinal changes in this condition [3,4,5]. ACHM has an early age of onset and generally manifests itself with pendular nystagmus at birth or early infancy as a first sign. Affected individuals present photophobia, eccentric fixation, reduced best corrected visual acuity (BCVA) since birth, reduction or complete loss of color perception along the three axes of color vision and small central scotomas [6]. A broad range of refractive errors has been reported in ACHM patients, with hyperopic refraction being the most frequent [7]. The fundus appearance is often normal, but abnormal foveal reflex, pigmentary mottling and atrophic changes may be observed in the macula [5]. In terms of retinal activity, patients show markedly decreased or absent photopic responses (including the 30 Hz flicker response) in full-field electroretinograms (ERG), whereas scotopic responses are normal or only slightly abnormal [8]. The structure of the central retina can be compromised with variable degrees of foveal hypoplasia, disruption and/or loss of the photoreceptors’ ellipsoid zone (EZ band) as well as an attenuation of the retinal pigment epithelium (RPE) at the macula as revealed by optical coherence tomography (OCT) [9]. Fundus autofluorescence imaging (FAF) shows different extents of foveal hypofluorescence. In some cases, a hyper-autofluorescent ring with a central region devoid of autofluorescence (AF) was also reported [3].
There are two clinical forms of ACHM: the complete and incomplete subtypes. Unlike complete achromatopsia, in which cone function is totally abolished leading to profound color blindness, in the incomplete forms the retina retains some residual cone photoreceptor activity as not all three cone spectral types are impaired. These incomplete forms, also termed atypical achromatopsia or dyschromatopsia, are characterized by detectable cone-mediated responses in ERG, a higher BCVA (up to 20/80) and a better-preserved color vision [10].
Pathogenic variants in six genes encoding components of the cone phototransduction cascade (CNGA3, CNGB3, PDE6C, PDE6H, GNAT2) and of the unfolded protein response (ATF6) have been implicated in ACHM [11,12,13,14,15,16,17,18]. Collectively, these genes account for over 90% of ACHM cases, while CNGA3 and CNGB3 alone are responsible for 69% of cases [19]. Gene therapy trials are ongoing to correct the defects of CNGA3 and CNGB3 that encode, respectively, the α- and β-subunits of the cone photoreceptor cyclic nucleotide-gated ion channel (NCT02599922, NCT03001310, NCT02610582, NCT02935517, NCT03758404) [20,21].
ACHM serves as the ideal target for gene therapies directed to cone photoreceptors, because cones are primarily affected in this disease. Defining the disease course can have important implications for therapies, hence the need to conduct longitudinal studies. Although several studies [4,5,9,22,23,24,25] looked into the clinical course of ACHM patients to assess disease progression, there is still an open debate on whether ACHM is a progressive or stable condition. In particular, some reports showed that at least a subset of patients present progressive vision loss with impaired cone photoreceptor function and poorly preserved foveal integrity [4,5,9,26,27]. In this study, we performed a detailed longitudinal phenotyping and genetic characterization of an Italian cohort of ACHM patients. We describe in depth the clinical features and disease progression in ACHM and investigate possible genotype–phenotype correlations. Our findings, together with similar studies, can serve as a baseline to identify key elements of retinal structure and function in ACHM that can optimize patient selection and intervention timing in forthcoming gene replacement therapies for this condition.

2. Results

2.1. Patient Selection

The study cohort consisted of 21 patients (mean age at the study baseline: 18.0 ± 3.2 years) (Table S1) from 17 unrelated families. The cohort included three sibling pairs and an uncle–nephew pair. All patients had a clinical diagnosis of ACHM. The main baseline clinical findings are summarized in Table 1. Specifically, 14 patients (66.7%) presented complete ACHM and seven patients (33.3%) had an incomplete form. All patients reported photophobia and 19 patients (90.5%) had congenital nystagmus.

2.2. Ophthalmological Assessment of the ACHM Patients

The mean spherical equivalent was +0.16 ± 0.63 D (ranging from −9.0 D to +5 D) (Table S1). The most common refractive error was hyperopia, observed in nine patients (42.9%), whereas myopia was found only in three patients (14.3%). In terms of visual acuity, the mean BCVA at the study baseline was 0.88 ± 0.06 logMAR in both eyes (Table S1) and was generally symmetrical between the two eyes (mean inter-eye asymmetry: 0.04 ± 0.02 logMAR). Only two patients (9.5%) showed values of inter-eye asymmetry equal or higher than the established threshold for clinically significant changes (i.e., 0.3 logMAR). There was no significant correlation between BCVA and age (p = 0.102).
The longitudinal analysis of 16 patients over a mean follow-up of 5.4 ± 1.0 years showed a statistically significant decline of BCVA at a mean annual rate of 0.012 logMAR/year (95% confidence interval (CI): 0.007 to 0.018; standard error (SE): 0.003; p < 0.001). During the follow-up, BCVA decreased in 10 out of 16 patients (62.5%). In particular, BCVA worsened bilaterally in four patients (25.0%), whereas six patients (37.5%) presented a unilateral worsening (Figure 1). At the last observation (mean age: 23.5 ± 2.7 years), the mean BCVA was 0.97 ± 0.05 logMAR in both eyes.
All patients performed the Ishihara Pseudoisochromatic Plate test, while 14 patients were also administered the Farnsworth D-15 color test. The Ishihara color vision test ascertained the inability of all patients to read any pseudoisochromatic plate apart from the test plate. The Farnsworth D-15 color test revealed a complete loss of color perception along all three axes of color vision in 12 patients, whereas one patient (P19) was classified as having deuteranopia and tritanopia. One patient (P18) failed to perform the Farnsworth D-15 test due to low BCVA and poor collaboration (Table 1).
All patients had clear optic media (cornea and lens) and a normal appearing fundus (Figure 2), except for one case (P10) who had RPE macular dystrophy in both eyes (Figure 2m).

2.3. Retinal Structure in ACHM Patients

At the baseline OCT, nine patients (42.9%) had a continuous EZ band in both eyes (foveal structure of grade 1) (Figure 2c), five (23.8%) had an EZ band disruption (grade 2) (Figure 2f), one patient (4.8%) had an absent EZ band (grade 3) (Figure 2i), three patients (14.3%) had a hyporeflective zone (HRZ) (grade 4) (Figure 2l), and one patient (4.8%) showed outer retinal atrophy including RPE loss (grade 5) in both eyes (Figure 2o). Representative examples of each OCT grade are shown in Figure 2. Moreover, two patients had a HRZ in one eye and either an EZ band disruption (4.8%) or an absent EZ band in the contralateral eye (4.8%) (Table 1). Significant differences (ANOVA, p = 0.02) were observed in the age of the patients stratified according to the OCT classification, when considering the groups that comprised at least three cases (i.e., grades 1, 2, 4). In particular, grade 4 patients were significantly older (35.3 ± 2.6 years) than grade 1 (11.0 ± 8.6 years; p = 0.006) and grade 2 (18.0 ± 7.8 years; p = 0.048).
Mean central retinal thickness (CRT) was 224.4 ± 9.8 µm in the right eyes and 213.9 ± 8.2 µm in the left eyes (Table S1). Only three patients (14.3%) showed foveal hypoplasia bilaterally (Table 1). Cross-sectional analysis showed that CRT significantly decreased with patient’s age at baseline (β: −1.3; SE: 0.6; p = 0.042), whereas the longitudinal data did not show significant changes in CRT during the follow-up period (p = 0.993). Furthermore, in terms of OCT grading, only two patients (12.5%) progressed to a worse OCT stage over the follow-up: one patient from grade 2 to 3 in one eye only (P8; left eye), the other from grade 1 to 2 in both eyes (P2) (Figure S1).
FAF imaging (available for 14 patients) revealed significant abnormalities in 11 subjects: one patient (7.1%) showed a central region of absent AF with a surrounding hyper-autofluorescent ring, three patients (21.4%) presented a reduced AF with subtle hyper-autofluorescence around the fovea, three patients (21.4%) had widespread reduced AF, and four patients (28.6%) showed foveal and parafoveal hyper-autofluorescence (Table 1). We also observed that the patient with a grade 5 OCT had a central region devoid of AF surrounded by a hyper-autofluorescent ring (Figure 2n). OCT grades 3 and 4 were associated with foveal and parafoveal hyper-autofluorescence (three out of five cases) (Figure 2h,k). On the other hand, patients belonging to OCT grades 1 and 2 (eight cases) had variable distributions of fundus autofluorescence that included a reduced AF with subtle hyper-autofluorescence around the fovea (two out of eight cases) (Figure 2e), a foveal and parafoveal hyperautofluorescence (one out of eight cases), a widespread reduced AF (two out of eight cases) or a normal AF (three out of eight cases) (Figure 2b, Table 1).

2.4. Retinal Function in the ACHM Patients

Microperimetric assessment (available for 18 patients) showed a reduced mean macular sensitivity (MS) in both eyes (14.6 ± 0.9 dB in right eyes, 15.2 ± 0.9 dB in left eyes) at the baseline (Table S1). MS decreased significantly with patient’s age (β = −0.312 dB/year; SE: 0.09; p = 0.001). It declined further during the follow-up period (5.4 ± 1.0 years) at a mean annual rate of −0.16 dB/year (SE: 0.05; p = 0.003). In terms of fixation stability, one patient (5.6%) showed stable fixation in both eyes, seven (38.9%) had a relatively stable fixation and two (11.1%) patients had unstable fixation in both eyes. The remaining patients (44.4%) showed inter-ocular asymmetry in the fixation stability: stable vs. relatively stable fixation (one patient) and relatively stable vs. unstable fixation (seven patients) (Table 1).
At the baseline ERG assessment which was performed for the entire cohort (n = 21), all patients showed normal dark-adapted 0.01 ERG responses (mean b-wave amplitude: 147.2 ± 11.7 µV in the right eyes; 163.0 ± 10.9 µV in the left eyes) (Table S1). On the contrary, light-adapted 3.0 ERG responses were undetectable in most patients (15/21; 71.4%) (Table 1). Only six patients (28.6%) had detectable light-adapted 3.0 ERG (mean age: 13.1 ± 2.7 years) with a mean b-wave amplitude of 12.2 ± 3.5 µV in the right eyes and 14.7 ± 7.4 µV in the left eyes (Table S1). Similarly, flicker ERG responses were extremely reduced, yet detectable, in the majority of patients (14/21; 66.7%) (Table 1), with a mean N1-P1 amplitude of 5.1 ± 1.2 µV in the right eyes and of 6.0 ± 1.4 µV in the left eyes (Table S1). Moreover, detectable light-adapted 3.0 responses associated with a significantly higher CRT (β = 18.9 µm; p = 0.031), higher MS (β = 6.6 dB; p < 0.001), and with higher values of b-wave amplitudes in the dark-adapted response (β = 61.9 µV; p = 0.004). Finally, detectable flicker responses correlated with a better BCVA (β = −0.46 logMAR; p = 0.002) and higher MS (β = 4.6 dB; p = 0.007).
Longitudinal analysis (mean follow-up: 5.4 ± 0.9 years) did not show a significant reduction neither of the b-wave amplitude in dark- (−0.6 µV per year; p = 0.624) and light-adapted 3.0 ERG responses (−0.3 µV per year; p = 0.318) nor of the N1-P1 amplitude in flicker ERG responses (−0.4 µV per year; p = 0.084). However, we observed a reduction of light-adapted 3.0 ERG below noise level in three (P1, P9, P12) out of the six patients, and reduction below noise level of flicker ERG in only two patients out of 14 (P8, P9).

2.5. Genetic Analysis of ACHM Patients

To determine the genetic defect underlying their condition, we performed high-throughput sequencing analysis of patients who gave consent to genetic testing (n = 18; P19–P21 did not consent). These samples were analyzed either by panel-based approaches or by clinical exome sequencing. One case (P18) that remained unsolved after a panel-based analysis of 130 retinopathy genes (RETplex) was subsequently analyzed by whole-exome sequencing (WES).
Overall, we identified causative variants that could explain the clinical phenotype in 16 ACHM patients (from 13 unrelated families) (Table 2). These patients harbored biallelic variants in one of the six genes that are more commonly associated with ACHM. Specifically, seven patients had biallelic mutations in CNGA3, five patients in CNGB3, three in GNAT2, and one in PDE6C (Table 2). Two patients with typical ACHM remained unsolved. They included a female patient (P17) with a monoallelic missense variant of unknown significance (c.1055G>A; p.Arg352Lys) in CNGB3 that alone cannot explain the clinical phenotype, and a male patient (P18) analyzed by WES following a negative first-tier screening using the RETplex panel.
A third of the causative variants identified in the cohort (i.e., 4 out of 12) were novel and absent from reference databases (e.g., HGMD, LOVD). One variant (GNAT2, c.832dup) was predicted to generate a frameshift leading to a premature stop codon. The remaining mutations were missense variants of uncertain significance according to the classification of the American College of Medical Genetics, predicted in silico to be deleterious for the protein function (Table 3).

2.6. Genotype-Phenotype Correlation Analysis

To explore possible genotype–phenotype correlations, we sought for differences in the clinical presentation between the patients with mutations in the GNAT2, CNGB3 and CNGA3 genes. There were no significant differences in the selected clinical parameters among the three genotypic groups (Table 4). Nonetheless, we observed a higher frequency of detectable flicker ERG responses in CNGB3-ACHM (80%) and GNAT2-ACHM patients (100%) compared to CNGA3-ACHM patients (28.6%) (Table 4).

3. Discussion

Here we provide the first detailed clinical description and molecular characterization of an Italian cohort of ACHM patients. To the best of our knowledge, this is one of the longest longitudinal studies of ACHM patients with a mean follow-up period of more than 5 years. Delineating the clinical features and the natural history of ACHM is both relevant and timely, given that phase I/II gene supplementation trials aimed at patients with causative mutations in CNGA3 and CNGB3 are ongoing [20,21].
Our results suggest that ACHM is characterized by a progressive decline in visual function and a slow deterioration of the macular structure. First, we found that BCVA worsened significantly in 62.5% of our patients during a follow-up period of approximately 5 years. Previous studies reported a progressive decline in BCVA in a smaller percentage of cases (about 10%) [9,31], but these discrepancies may be related to differences in the follow-up length or in the frequency of mutated genes among the cohorts. Second, our findings indicate that the integrity of retinal structure is overall well preserved during the course of the disease within the time-frame studied. In particular, CRT did not change significantly during the follow-up period and only two patients (12.5%) progressed to a more severe OCT grade. Along the same line, other studies reported a progressive loss of foveal integrity in a subset of ACHM patients (about 5–10%) over a mean follow-up between 2 and 5 years [3,4,5,9,22]. However, we found that older patients had a lower CRT and a more advanced OCT grade (i.e., grade 4 patients were mainly in the fourth decade of life, whereas grade 1 and grade 2 patients were mainly in the first two decades) suggesting that changes in retinal structure occur slowly, for instance over decades of life. Therefore, monitoring patients for periods longer than 5 years may reveal progressive changes in a greater proportion of individuals. Moreover, extended follow-up periods may also enable to assess whether the aging of the retina progresses at a different speed in patients with achromatopsia.
As regards the retinal function, we found that MS negatively correlated with patient’s age and decreased significantly over the follow-up period. To date, only few studies assessed the MS of ACHM patients by MP [9,23,24,32]. Our findings are in agreement with the observations of Sundram et al. who reported a decline of MS with age in 40 patients [23]. Conversely, Aboshiha et al. and more recently Georgiou et al. did not observe significant changes in MS in 38 and 18 patients, respectively [9,24]. These discrepancies may be due to the shorter follow-up reported in the first case (i.e., 19.5 months) [9] and to the different distribution of mutated genes [24]. Finally, we observed unrecordable or markedly reduced light-adapted responses that did not change significantly during the follow-up. On the other hand, the amplitude of the photopic ERG decreased below noise level in half of patients with detectable photopic responses at baseline. Likewise, Khan et al. observed that some patients retain low levels of residual cone function into adulthood which is then progressively lost [33]. Taken together, the above functional read-outs suggest that patients with ACHM may have substantial residual cone function despite having a significantly abnormal ERG response and lower MS. In light of these findings, MP could provide useful insights into the retinal function of ACHM patients and guide therapeutic decision-making and patient selection for inclusion in clinical trials.
In our cohort, we identified biallelic pathogenic variants in four of the six genes that are typically implicated in ACHM, namely CNGA3, CNGB3, GNAT2, PDE6C. Mutations in CNGA3 and CNGB3 accounted for the majority of molecularly defined cases (75%, 12 out of 16 cases), consistent with the high prevalence of CNGA3 and CNGB3 mutations in autosomal recessive achromatopsia [28,34]. The most recurrent CNGB3 variant was the c.1148delC (representing 80% of the identified CNGB3 alleles), in line with the observed frequency of this allele in other CNGB3-ACHM cohorts [13,34,35]. We also identified two novel missense mutations in CNGA3 (c.1162G>A; p.(Gly388Ser)) and CNGB3 (c.1055G>A; p.(Arg352Lys)) that are absent from population databases and extend the list of known likely pathogenic variants. Additional novel disease-causing variants were identified in GNAT2, which was less frequently mutated in our cohort. The majority of solved cases (81%, 13 out of 16 cases) were homozygous for the identified variants, as reported also in other ACHM cohorts [36].
Two analyzed cases remained unsolved. One subject was heterozygous for a likely pathogenic variant in CNGB3. We believe that this finding is not incidental, given the low allele frequency of the identified variant (absent from gnomAD) and the patient’s clinical presentation of typical ACHM. In this case, the missing causative variant could lie in noncoding sequences, such as a deep-intronic splice variant that causes pseudoexon activation in CNGB3 [35]. Another hypothesis may be the presence of a copy number variation (CNV), which is expected to be undetected by our capture-based targeted NGS screening. CNGB3 was predicted to be prone to CNVs due to its genomic features and high Alu repeat content [37], and complementary targeted CNV screening proved decisive in the analysis of monoallelic CNGB3 patients [36]. These explanations are valid also for the second case that presented the hallmarks of ACHM but remained unsolved even after WES.
In an attempt to assess whether mutations in the three most prominent genes of the cohort (CNGA3, CNGB3 and GNAT2) associate with distinct clinical parameters, we performed a genotype–phenotype correlation analysis. We did not identify statistically significant differences in the disease course and clinical features among the three genotypic groups, which may be attributable to the limited size of the analyzed cohort. Since it was hypothesized that foveal hypoplasia may represent a gene-specific feature [38] (as it was observed in the majority of previously reported CNGA3- and CNGB3-ACHM subjects but not in patients carrying defects in GNAT2 or PDE6C [38]), we looked into this feature in relation to the patient’s genotype. In our cohort, foveal hypoplasia was not observed in GNAT2- and PDE6C-ACHM patients while it was observed only in a subgroup of patients harboring pathogenic variants in CNGB3 and CNGA3, although at a much lower frequency (about 17%).
This study has some limitations in its retrospective design. First, although the patients were invited to undergo annual follow-up visits, the follow-up length was not the same for all patients and only cross-sectional data were available for some patients or some tests. Moreover, microperimetric measurements and FAF imaging were available only for a subset of the cohort because some patients were not always willing or capable to perform the examination (i.e., because of poor BCVA, young age). Finally, we acknowledge that the relatively small sample size may hinder the detection of statistically significant associations in some comparisons (i.e., groups by genotype) due to the insufficient power of the test.
In conclusion, a better understanding of the disease course and progression are crucial for patient management, with important implications for therapeutic considerations. Our findings sustain that ACHM is a progressive disease in terms of BCVA, MS and ERG responses and affects slowly the structural integrity of the retina. In particular, MP could be useful in ACHM as an informative test to guide therapeutic decisions and establish appropriate timing for interventions. Ultimately, our study can contribute towards the development of guidelines for patient selection and intervention timing in forthcoming gene replacement therapies.

4. Materials and Methods

4.1. Ethics Statement

The retrospective observational study was performed at the Referral Centre for Inherited Retinopathies of the Eye Clinic, Multidisciplinary Department of Medical, Surgical and Dental Sciences, University of Campania “Luigi Vanvitelli”. All procedures of this study adhered to the tenets of the Declaration of Helsinki and were approved by the Ethics Board of the Università degli Studi della Campania “Luigi Vanvitelli” (for adult protocol n. 8189/2015, 9 April 2015; for pediatric subjects’ protocol n. 500/2017, 12 September 2017). Peripheral blood samples were collected upon written informed consent of the patients to sample collection and genetic analysis. For minors, informed consent was obtained by the parents or legal guardians.

4.2. Patients Inclusion Criteria

All patients evaluated at the Referral Centre for Inherited Retinopathies of the Eye Clinic of the University of Campania “Luigi Vanvitelli” with a clinical diagnosis of ACHM were included in the study. The clinical diagnosis of ACHM was performed according to previously established criteria [39]. Specifically, the diagnostic criteria included pendular nystagmus, photophobia, reduced BCVA (varies from 20/200 or less in complete ACHM and may be as high as 20/80 in incomplete ACHM), reduction (incomplete ACHM) or complete loss (complete ACHM) of color perception along all three axes of color vision, normal fundus appearance in many affected patients (although absence of the foveal reflex, pigment mottling or atrophy in the macular region may occur in older individuals), and absent (complete ACHM) or residual (incomplete ACHM) cone responses with normal rod responses in ERG.

4.3. Ophthalmological Examination

The patient’s medical history was obtained prior to examination. All patients underwent a full ophthalmological examination which included BCVA measurements, color vision testing, slit lamp anterior segment examination, fundus examination, OCT, FAF, MP, and ERG. For 16 subjects, follow-up measurements of BCVA, OCT, MP, and ERG were available for a mean period of 5.4 ± 1.0 years.
Color vision testing was carried out using the Pseudoisochromatic Plate Ishihara and/or the Farnsworth D-15 test.
Based on the fundus appearance, each patient was classified into one of the following categories: (i) normal (no RPE disturbance), (ii) RPE macular dystrophy, or (iii) RPE macular atrophy [9].
OCT was performed on both eyes following pupillary dilation (tropicamide 1% and phenylephrine 2.5% eye drops), using a spectral domain OCT (Spectralis OCT plus with blue peak, Heidelberg Engineering, Heidelberg, Germany) by experienced operators (RBP, VDI, FT). At follow-up visits, the device was switched to the follow-up mode in order to scan the same location both at baseline and in subsequent follow-ups. Qualitative assessment of foveal structure was performed by grading spectral domain OCT images into one of the five distinct categories, according to Aboshiha et al. [9]: continuous EZ band (grade 1), EZ band disruption (grade 2), EZ band absence (grade 3), presence of a HRZ (grade 4), outer retinal atrophy including RPE loss (grade 5). The presence or absence of foveal hypoplasia, defined as the persistence of one or more inner retinal layers (outer plexiform layer, inner nuclear layer, inner plexiform layer or ganglion cell layer) through the fovea, was also assessed [23]. Retinal thickness was calculated manually by identifying the thinnest point in the fovea.
FAF imaging was performed subsequently to the acquisition of the OCT images, due to the relatively intense lights used during FAF acquisition and the photophobic nature of ACHM patients. Pupillary dilatations were performed as described above. Consensus OCT grading and FAF pattern were established by two authors (RBP, VDI). When interpretations differed, a discussion was made to reach an agreement and a third author (FT) also gave his interpretation, if necessary.
MP was performed by an automatic fundus-related perimeter (MP1 Microperimeter, Nidek Technologies, Padova, Italy). The following parameters were used: a fixation target of 2 degrees in diameter consisting of a red ring; a white, monochromatic background with a luminance of 1.27 cd/m²; a Goldmann III-size stimulus with a projection time of 200 ms; a predefined automatic test pattern (Humphrey 10-2) covering 10 degrees centered onto the gravitational center of the fixation points with 68 stimuli [40].
ERG was recorded according to the International Guidelines of the International Society of Clinical Electrophysiology of Vision (ISCEV) [41,42].

4.4. Target Enrichment and Next-Generation Sequencing

Samples were analyzed either by panel-based analysis of known retinopathy genes or by clinical exome sequencing. One sample underwent whole-exome sequencing. Genomic DNA was extracted from peripheral blood using the DNeasy Blood and Tissue Kit (QIAGEN) according to standard protocols. Targeted NGS analysis was performed using the RETplex targeted sequencing panel as previously described [30]. Clinical exome and whole-exome sequencing libraries were prepared using the ClearSeq Inherited Disease Panel (Agilent Technologies, Santa Clara, CA, USA) and the SureSelect Human All Exon V7 (Agilent Technologies, Santa Clara, CA, USA), respectively. Targeted regions were enriched using the SureSelectQXT Target Enrichment system (Agilent Technologies, Santa Clara, CA, USA). Libraries were run on a NextSeq500 sequencing platform (Illumina inc., San Diego, CA, USA). Sequencing data were analyzed using a previously described pipeline [43].

4.5. Variant Interpretation and Validation

The single nucleotide variants (SNVs) and indel variants were annotated using ANNOVAR [44] with the relative position in genes using the RefSeq87 gene model, amino acid change, presence in dbSNP v137, frequency in the EXAC (http://exac.broadinstitute.org) and gnomAD databases (http://gnomad.broadinstitute.org), 1000 genomes project [45], presence in the Human Gene Mutation Database (HGMD) [46], Clinvar database [47], multiple cross-species conservation [48], and prediction scores of damage on protein activity [49,50,51,52]. Only variants with minor allele frequency (MAF) <0.01 were considered. The alignments at candidate positions were visually inspected using the Integrative Genomics Viewer (IGV).
The identified variants were validated by Sanger sequencing of the corresponding genomic fragments. Amplicons were obtained by PCR using Taq polymerase according to standard protocols and were Sanger sequenced. Mutation detection was performed using the CodonCode Aligner software. Segregation analysis of the identified variants was performed whenever samples from additional family members were available.

4.6. Statistical Analysis

Continuous variables are reported as mean ± standard error of mean (SEM) and categorical variables are reported as counts (percentage).
The natural history of disease was analyzed using previously applied methods [53,54]. In order to include all the patients (i.e., also those with no longitudinal data), we performed a cross-sectional analysis on the data of the baseline visit for each of the selected outcome measures (i.e., BCVA, CRT, MS). In particular, linear regression models, estimated by a generalized estimating equation (GEE), were performed with each outcome measure as the dependent variable and age as the independent variable in order to estimate a slope (mean rate of change per year of age) with its relative standard error (SE). Moreover, the GEE was fitted on longitudinal data, using baseline values as offset and follow-up length expressed in years as an independent variable, to estimate the annual change of each outcome over the follow-up period. Finally, regression models were fitted with each outcome measure as the dependent variable and mutated genes (i.e., CNGA3, CNGB3 and GNAT2) as independent variable, in order to investigate whether the genotype influenced the selected outcome.
GEE was applied because this method could accommodate the inter-eye correlation (i.e., between the two eyes of the same subject at a given visit) and longitudinal correlation (i.e., between values of the same eye followed over time) by adopting an appropriate covariance structure, as previously described [55]. Intercepts were included in all the models, except those regarding BCVA, since we assume that BCVA should start with the value of 0 logMAR (corresponding to healthy status). Moreover, asymmetry in BCVA between the two eyes was defined as a difference of 0.3 logMAR (15 Early Treatment Diabetic Retinopathy Study (ETDRS) letters), which is the threshold for clinical significance of BCVA changes [56].
Finally, analysis of variance (ANOVA) with Bonferroni correction for multiple comparisons performed with t-tests was adopted to evaluate differences in the age between the patients classified according to OCT grade.
p-values (p) lower than 0.05 were considered statistically significant. These statistical analyses were performed using the IBM SPSS Statistics platform (Version 21.0.0.0).

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/22/4/1681/s1.

Author Contributions

Conceptualization, F.S.; Formal analysis, R.B.-P., M.K., P.M. and F.T.; Funding acquisition, S.B. and F.S.; Investigation, R.B.-P., M.K., P.M., V.D.I., A.D.B., G.I. and F.T.; Methodology, R.B.-P., M.K., P.M., V.D.I. and F.T.; Project administration, F.T., S.B. and F.S.; Resources, S.B. and F.S.; Supervision, F.T., S.B. and F.S.; Visualization, R.B.-P., M.K. and F.T.; Writing—original draft, R.B.-P., M.K., P.M. and F.T.; Writing—review and editing, F.T., S.B. and F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University of Campania “Luigi Vanvitelli” under “VALERE: VAnviteLli pEr la RicErca” (to S.B. and F.S.). Some parts of this study were supported by grants from the Italian Fondazione Roma (to S.B. and F.S.) and from the Ministero dell “Istruzione dell” Università e della Ricerca (MIUR) under PRIN 2015 n 20158Y77NT (to S.B. and F.S.).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of the Università degli Studi della Campania “Luigi Vanvitelli” (for adult protocol n. 8189/2015; for pediatric subjects’ protocol n. 500/2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data is contained within the article or supplementary material.

Acknowledgments

We thank the participating patients and their families. We are grateful to Mariateresa Pizzo for technical support. We also wish to acknowledge the European Reference Network dedicated to Rare Eye Diseases (ERN-EYE) for useful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACHMAchromatopsia
AFAutofluorescence
BCVABest corrected visual acuity
CRTCentral retinal thickness
ERGElectroretinogram
FAFFundus autofluorescence
HGMDHuman Gene Mutation Database
LOVDLeiden Open Variation Database
MPMicroperimentry
MSMacular sensitivity
NGSNext generation sequencing
OCTOptical coherence tomography
RPERetinal pigment epithelium
SNVSingle nucleotide variants

References

  1. Sharpe, L.T.; Stockman, A.; Jagle, H. Opsin Genes, Cone Photopigments amd Colourblindness; Cambridge University Press: Cambridge, UK, 1999; pp. 3–52. [Google Scholar]
  2. Michaelides, M.; Hunt, D.M.; Moore, A.T. The cone dysfunction syndromes. Br. J. Ophthalmol. 2004, 88, 291–297. [Google Scholar] [CrossRef]
  3. Fahim, A.T.; Khan, N.W.; Zahid, S.; Schachar, I.H.; Branham, K.; Kohl, S.; Wissinger, B.; Elner, V.M.; Heckenlively, J.R.; Jayasundera, T. Diagnostic Fundus Autofluorescence Patterns in Achromatopsia. Am. J. Ophthalmol. 2013, 156, 1211–1219.e2. [Google Scholar] [CrossRef]
  4. Greenberg, J.P.; Sherman, J.; Zweifel, S.A.; Chen, R.W.S.; Duncker, T.; Kohl, S.; Baumann, B.; Wissinger, B.; Yannuzzi, L.A.; Tsang, S.H. Spectral-Domain Optical Coherence Tomography Staging and Autofluorescence Imaging in Achromatopsia. JAMA Ophthalmol. 2014, 132, 437–445. [Google Scholar] [CrossRef]
  5. Thiadens, A.; Somervuo, V.; Born, L.V.D.; Roosing, S.; Van Schooneveld, M.; Kuijpers, R.; Van Moll-Ramirez, N.; Cremers, F.; Hoyng, C.; Klaver, C.C.W. Progressive Loss of Cones in Achromatopsia: An Imaging Study Using Spectral-Domain Optical Coherence Tomography. Investig. Opthalmol. Vis. Sci. 2010, 51, 5952–5957. [Google Scholar] [CrossRef] [Green Version]
  6. Hirji, N.; Aboshiha, J.; Georgiou, M.; Bainbridge, J.; Michaelides, M. Achromatopsia: Clinical features, molecular genetics, animal models and therapeutic options. Ophthalmic Genet. 2018, 39, 149–157. [Google Scholar] [CrossRef]
  7. Simunovic, M.P.; Moore, A. The cone dystrophies. Eye 1998, 12, 553–565. [Google Scholar] [CrossRef] [Green Version]
  8. Andréasson, S.; Tornqvist, K. Electroretinograms in patients with achromatopsia. Acta Ophthalmol. 1991, 69, 711–716. [Google Scholar] [CrossRef] [PubMed]
  9. Aboshiha, J.; Dubis, A.M.A.; Cowing, J.; Fahy, R.T.A.; Sundaram, V.; Bainbridge, J.W.; Ali, R.R.; Dubra, A.; Nardini, M.; Webster, A.R.; et al. A Prospective Longitudinal Study of Retinal Structure and Function in Achromatopsia. Investig. Opthalmol. Vis. Sci. 2014, 55, 5733–5743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Pokorny, J.; Smith, V.C.; Pinckers, A.J.L.G.; Cozijnsen, M. Classification of complete and incomplete autosomal recessive achromatopsia. Graefe’s Arch. Clin. Exp. Ophthalmol. 1982, 219, 121–130. [Google Scholar] [CrossRef] [PubMed]
  11. Ansar, M.; University of Washington Center for Mendelian Genomics; Santos-Cortez, R.L.P.; Saqib, M.A.N.; Zulfiqar, F.; Lee, K.; Ashraf, N.M.; Ullah, E.; Wang, X.; Sajid, S.; et al. Mutation of ATF6 causes autosomal recessive achromatopsia. Qual. Life Res. 2015, 134, 941–950. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Chang, B.; Grau, T.; Dangel, S.; Hurd, R.; Jurklies, B.; Sener, E.C.; Andreasson, S.; Dollfus, H.; Baumann, B.; Bolz, S.; et al. A homologous genetic basis of the murine cpfl1 mutant and human achromatopsia linked to mutations in the PDE6C gene. Proc. Natl. Acad. Sci. USA 2009, 106, 19581–19586. [Google Scholar] [CrossRef] [Green Version]
  13. Kohl, S.; Baumann, B.; Broghammer, M.; Jägle, H.; Sieving, P.; Kellner, U.; Spegal, R.; Anastasi, M.; Zrenner, E.; Sharpe, L.T.; et al. Mutations in the CNGB3 gene encoding the beta-subunit of the cone photoreceptor cGMP-gated channel are responsible for achromatopsia (ACHM3) linked to chromosome 8q21. Hum. Mol. Genet. 2000, 9, 2107–2116. [Google Scholar] [CrossRef]
  14. Kohl, S.; Baumann, B.; Rosenberg, T.; Kellner, U.; Lorenz, B.; Vadalà, M.; Jacobson, S.G.; Wissinger, B. Mutations in the Cone Photoreceptor G-Protein α-Subunit Gene GNAT2 in Patients with Achromatopsia. Am. J. Hum. Genet. 2002, 71, 422–425. [Google Scholar] [CrossRef] [Green Version]
  15. Kohl, S.; Coppieters, F.; Meire, F.; Schaich, S.; Roosing, S.; Brennenstuhl, C.; Bolz, S.; Van Genderen, M.M.; Riemslag, F.C.; Lukowski, R.; et al. A Nonsense Mutation in PDE6H Causes Autosomal-Recessive Incomplete Achromatopsia. Am. J. Hum. Genet. 2012, 91, 527–532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kohl, S.; Marx, T.; Giddings, I.; Jägle, H.; Jacobson, S.G.; Apfelstedt-Sylla, E.; Zrenner, E.; Sharpe, L.T.; Wissinger, B. Total colourblindness is caused by mutations in the gene encoding the alpha-subunit of the cone photoreceptor cGMP-gated cation channel. Nat. Genet. 1998, 19, 257–259. [Google Scholar] [CrossRef] [PubMed]
  17. Kohl, S.; Zobor, D.; Chiang, W.-C.; Weisschuh, N.; Staller, J.; Menendez, I.G.; Chang, S.; Beck, S.C.; Garrido, M.G.; Sothilingam, V.; et al. Mutations in the unfolded protein response regulator ATF6 cause the cone dysfunction disorder achromatopsia. Nat. Genet. 2015, 47, 757–765. [Google Scholar] [CrossRef] [Green Version]
  18. Felden, J.; Baumann, B.; Ali, M.; Audo, I.; Ayuso, C.; Bocquet, B.; Casteels, I.; Garcia-Sandoval, B.; Jacobson, S.G.; Jurklies, B.; et al. Mutation spectrum and clinical investigation of achromatopsia patients with mutations in the GNAT2 gene. Hum. Mutat. 2019, 40, 1145–1155. [Google Scholar] [CrossRef] [PubMed]
  19. Weisschuh, N.; Stingl, K.; Audo, I.; Biskup, S.; Bocquet, B.; Branham, K.; Burstedt, M.S.; De Baere, E.; De Vries, M.J.; Golovleva, I.; et al. Mutations in the gene PDE6C encoding the catalytic subunit of the cone photoreceptor phosphodiesterase in patients with achromatopsia. Hum. Mutat. 2018, 39, 1366–1371. [Google Scholar] [CrossRef]
  20. Fischer, M.D.; Michalakis, S.; Wilhelm, B.; Zobor, D.; Muehlfriedel, R.; Kohl, S.; Weisschuh, N.; Ochakovski, A.; Klein, R.; Schoen, C.; et al. Safety and Vision Outcomes of Subretinal Gene Therapy Targeting Cone Photoreceptors in Achromatopsia: A Nonrandomized Controlled Trial. JAMA Ophthalmol. 2020, 138, 643–651. [Google Scholar] [CrossRef]
  21. Thompson, D.A.; Iannaccone, A.; Ali, R.R.; Arshavsky, V.Y.; Audo, I.; Bainbridge, J.W.B.; Besirli, C.G.; Birch, D.G.; Branham, K.E.; Cideciyan, A.V.; et al. Advancing Clinical Trials for Inherited Retinal Diseases: Recommendations from the Second Monaciano Symposium. Transl. Vis. Sci. Technol. 2020, 9, 2. [Google Scholar] [CrossRef]
  22. Hirji, N.; Georgiou, M.; Kalitzeos, A.; Bainbridge, J.W.; Kumaran, N.; Aboshiha, J.; Carroll, J.; Michaelides, M. Longitudinal Assessment of Retinal Structure in Achromatopsia Patients With Long-Term Follow-up. Investig. Opthalmol. Vis. Sci. 2018, 59, 5735–5744. [Google Scholar] [CrossRef] [Green Version]
  23. Sundaram, V.; Wilde, C.; Aboshiha, J.; Cowing, J.; Han, C.; Langlo, C.S.; Chana, R.; Davidson, A.E.; Sergouniotis, P.I.; Bainbridge, J.W.; et al. Retinal structure and function in achromatopsia: Implications for gene therapy. Ophthalmology 2014, 121, 234–245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Georgiou, M.; Singh, N.; Kane, T.; Zaman, S.; Hirji, N.; Aboshiha, J.; Kumaran, N.; Kalitzeos, A.; Carroll, J.; Weleber, R.G.; et al. Long-Term Investigation of Retinal Function in Patients with Achromatopsia. Investig. Opthalmol. Vis. Sci. 2020, 61, 38. [Google Scholar] [CrossRef]
  25. Langlo, C.S.; Patterson, E.J.; Higgins, B.P.; Summerfelt, P.; Razeen, M.M.; Erker, L.R.; Parker, M.; Collison, F.T.; Fishman, G.A.; Kay, C.N.; et al. Residual Foveal Cone Structure in CNGB3-Associated Achromatopsia. Investig. Opthalmol. Vis. Sci. 2016, 57, 3984–3995. [Google Scholar] [CrossRef] [Green Version]
  26. Georgiou, M.; Robson, A.G.; Singh, N.; Pontikos, N.; Kane, T.; Hirji, N.; Ripamonti, C.; Rotsos, T.; Dubra, A.; Kalitzeos, A.; et al. Deep Phenotyping of PDE6C-Associated Achromatopsia. Investig. Opthalmol. Vis. Sci. 2019, 60, 5112–5123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Thomas, M.G.; McLean, R.J.; Kohl, S.; Sheth, V.; Gottlob, I. Early signs of longitudinal progressive cone photoreceptor degeneration in achromatopsia. Br. J. Ophthalmol. 2012, 96, 1232–1236. [Google Scholar] [CrossRef]
  28. Wissinger, B.; Gamer, D.; Jägle, H.; Giorda, R.; Marx, T.; Mayer, S.; Tippmann, S.; Broghammer, M.; Jurklies, B.; Rosenberg, T.; et al. CNGA3 Mutations in Hereditary Cone Photoreceptor Disorders. Am. J. Hum. Genet. 2001, 69, 722–737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Sundin, O.H.; Yang, J.-M.; Li, Y.; Zhu, D.; Hurd, J.N.; Mitchell, T.N.; Silva, E.D.; Maumenee, I.H. Genetic basis of total colourblindness among the Pingelapese islanders. Nat. Genet. 2000, 25, 289–293. [Google Scholar] [CrossRef]
  30. Di Iorio, V.; Karali, M.; Brunetti-Pierri, R.; Filippelli, M.; Di Fruscio, G.; Pizzo, M.; Mutarelli, M.; Nigro, V.; Testa, F.; Banfi, S.; et al. Clinical and Genetic Evaluation of a Cohort of Pediatric Patients with Severe Inherited Retinal Dystrophies. Genes 2017, 8, 280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Thiadens, A.A.H.J.; Slingerland, N.W.R.; Roosing, S.; Van Schooneveld, M.J.; Van Lith-Verhoeven, J.J.C.; Van Moll-Ramirez, N.; Born, L.I.V.D.; Hoyng, C.B.; Cremers, F.P.; Klaver, C.C.W. Genetic Etiology and Clinical Consequences of Complete and Incomplete Achromatopsia. Ophthalmology 2009, 116, 1984–1989.e1. [Google Scholar] [CrossRef]
  32. Genead, M.A.; Fishman, G.A.; Rha, J.; Dubis, A.M.; Bonci, D.M.O.; Dubra, A.; Stone, E.M.; Neitz, M.; Carroll, J. Photoreceptor Structure and Function in Patients with Congenital Achromatopsia. Investig. Opthalmol. Vis. Sci. 2011, 52, 7298–7308. [Google Scholar] [CrossRef]
  33. Khan, N.W.; Wissinger, B.; Kohl, S.; Sieving, P.A. CNGB3Achromatopsia with Progressive Loss of Residual Cone Function and Impaired Rod-Mediated Function. Investig. Opthalmol. Vis. Sci. 2007, 48, 3864–3871. [Google Scholar] [CrossRef] [Green Version]
  34. Kohl, S.; Varsanyi, B.; Antunes, G.A.; Baumann, B.; Hoyng, C.B.; Jägle, H.; Rosenberg, T.; Kellner, U.; Lorenz, B.; Salati, R.; et al. CNGB3 mutations account for 50% of all cases with autosomal recessive achromatopsia. Eur. J. Hum. Genet. 2004, 13, 302–308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Weisschuh, N.; Sturm, M.; Baumann, B.; Audo, I.; Ayuso, C.; Bocquet, B.; Branham, K.; Brooks, B.P.; Catalá-Mora, J.; Giorda, R.; et al. Deep-intronic variants in CNGB3 cause achromatopsia by pseudoexon activation. Hum. Mutat. 2020, 41, 255–264. [Google Scholar] [CrossRef] [PubMed]
  36. Mayer, A.K.; Van Cauwenbergh, C.; Rother, C.; Baumann, B.; Reuter, P.; De Baere, E.; Wissinger, B.; Kohl, S.; ACHM Study Group. CNGB3 mutation spectrum including copy number variations in 552 achromatopsia patients. Hum. Mutat. 2017, 38, 1579–1591. [Google Scholar] [CrossRef]
  37. Van Schil, K.; CNV Study Group; Naessens, S.; Van De Sompele, S.; Carron, M.; Aslanidis, A.; Van Cauwenbergh, C.; Mayer, A.K.; Van Heetvelde, M.; Bauwens, M.; et al. Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations. Genet. Med. 2017, 20, 202–213. [Google Scholar] [CrossRef] [Green Version]
  38. Georgiou, M.; Singh, N.; Kane, T.; Robson, A.G.; Kalitzeos, A.; Hirji, N.; Webster, A.R.; Dubra, A.; Carroll, J.; Michaelides, M. Photoreceptor Structure in GNAT2-Associated Achromatopsia. Investig. Opthalmol. Vis. Sci. 2020, 61, 40. [Google Scholar] [CrossRef] [Green Version]
  39. Kohl, S.; Jägle, H.; Wissinger, B.; Zobor, D. Achromatopsia. In GeneReviews®; Pagon, R.A., Ardinger, H.H., Wallace, S.E., Eds.; University of Washington: Seattle, WA, USA, 1993. [Google Scholar]
  40. Di Iorio, V.; Karali, M.; Melillo, P.; Testa, F.; Brunetti-Pierri, R.; Musacchia, F.; Condroyer, C.; Neidhardt, J.; Audo, I.; Zeitz, C.; et al. Spectrum of Disease Severity in Patients With X-Linked Retinitis Pigmentosa Due to RPGR Mutations. Investig. Opthalmol. Vis. Sci. 2020, 61, 36. [Google Scholar] [CrossRef]
  41. Marmor, M.F.; Fulton, A.B.; Holder, G.E.; Miyake, Y.; Brigell, M.; Bach, M.; for the International Society for Clinical Electrophysiology of Vision. ISCEV Standard for full-field clinical electroretinography (2008 update). Doc. Ophthalmol. 2008, 118, 69–77. [Google Scholar] [CrossRef] [Green Version]
  42. McCulloch, D.L.; Marmor, M.F.; Brigell, M.G.; Hamilton, R.; Holder, G.E.; Tzekov, R.; Bach, M. ISCEV Standard for full-field clinical electroretinography (2015 update). Doc. Ophthalmol. 2015, 130, 1–12. [Google Scholar] [CrossRef] [PubMed]
  43. Musacchia, F.; Ciolfi, A.; Mutarelli, M.; Bruselles, A.; Castello, R.; Pinelli, M.; Basu, S.; Banfi, S.; Casari, G.; Tartaglia, M.; et al. VarGenius executes cohort-level DNA-seq variant calling and annotation and allows to manage the resulting data through a PostgreSQL database. BMC Bioinform. 2018, 19, 477. [Google Scholar] [CrossRef]
  44. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput se-quencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  45. 1000 Genomes Project Consortium. A map of human genome variation from population-scale sequencing. Nature 2010, 467, 1061–1073. [Google Scholar] [CrossRef] [Green Version]
  46. Stenson, P.D.; Mort, M.; Ball, E.V.; Shaw, K.; Phillips, A.D.; Cooper, D.N. The Human Gene Mutation Database: Building a comprehensive mutation repository for clinical and mo-lecular genetics, diagnostic testing and personalized genomic medicine. Hum. Genet. 2014, 133, 1–9. [Google Scholar] [CrossRef] [Green Version]
  47. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Hoover, J.; et al. ClinVar: Public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016, 44, D862–D868. [Google Scholar] [CrossRef] [Green Version]
  48. Goode, D.L.; Cooper, G.M.; Schmutz, J.; Dickson, M.; Gonzales, E.; Tsai, M.; Karra, K.; Davydov, E.; Batzoglou, S.; Myers, R.M.; et al. Evolutionary constraint facilitates interpretation of genetic variation in resequenced human genomes. Genome Res. 2010, 20, 301–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Adzhubei, I.; Jordan, D.M.; Sunyaev, S.R. Predicting Functional Effect of Human Missense Mutations Using PolyPhen. Curr. Protoc. Hum. Genet. 2013, 76, 7.20.1–7.20.41. [Google Scholar] [CrossRef] [Green Version]
  50. Kumar, P.; Henikoff, S.; Ng, P.C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 2009, 4, 1073–1081. [Google Scholar] [CrossRef]
  51. Liu, X.; Jian, X.; Boerwinkle, E. dbNSFP: A lightweight database of human nonsynonymous SNPs and their functional predictions. Hum. Mutat. 2011, 32, 894–899. [Google Scholar] [CrossRef]
  52. Schwarz, J.M.; Rödelsperger, C.; Schuelke, M.; Seelow, D. MutationTaster evaluates disease-causing potential of sequence alterations. Nat. Methods 2010, 7, 575–576. [Google Scholar] [CrossRef]
  53. Kong, X.; Strauss, R.W.; Michaelides, M.; Cideciyan, A.V.; Sahel, J.-A.; Muñoz, B.; West, S.; Scholl, H.P.; Wolfson, Y.; Bittencourt, M.; et al. Visual Acuity Loss and Associated Risk Factors in the Retrospective Progression of Stargardt Disease Study (ProgStar Report No. 2). Ophthalmology 2016, 123, 1887–1897. [Google Scholar] [CrossRef]
  54. Testa, F.; Melillo, P.; Bonnet, C.; Marcelli, V.; De Benedictis, A.; Colucci, R.; Gallo, B.; Kurtenbach, A.; Rossi, S.; Marciano, E.; et al. Clinical presentation and disease course of usher syndrome because of mutations in myo7a or ush2a. Retina 2017, 37, 1581–1590. [Google Scholar] [CrossRef]
  55. Glynn, R.J.; Rosner, B. Regression methods when the eye is the unit of analysis. Ophthalmic Epidemiol. 2012, 19, 159–165. [Google Scholar] [CrossRef] [PubMed]
  56. Csaky, K.G.; Richman, E.A.; Iii, F.L.F.; Ferris, F.L. Report from the NEI/FDA Ophthalmic Clinical Trial Design and Endpoints Symposium. Investig. Opthalmol. Vis. Sci. 2008, 49, 479–489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 22 01681 g001 550
Figure 1. Change in the best corrected visual acuity in individual ACHM patients showing a worsening over the follow-up period. Baseline and last time-point visual acuities are represented by circles and triangles, respectively. Data from the right eye (RE) are in red and left eye (LE) in blue. The duration of the follow-up period (in years) for each patient is shown on the right. Patient ID is shown on the left.
Figure 1. Change in the best corrected visual acuity in individual ACHM patients showing a worsening over the follow-up period. Baseline and last time-point visual acuities are represented by circles and triangles, respectively. Data from the right eye (RE) are in red and left eye (LE) in blue. The duration of the follow-up period (in years) for each patient is shown on the right. Patient ID is shown on the left.
Ijms 22 01681 g001
Ijms 22 01681 g002 550
Figure 2. Ophthalmological findings in five representative ACHM patients showing the OCT grading (c,f,i,l,o) with corresponding retinography (a,d,g,j,m) and fundus AF images (b,e,h,k,n). (a–c) Patient P15 (12 years old, homozygous for the c.619G>A variant in GNAT2) presented a normal fundus appearance, normal fundus autofluorescence and a grade 1 OCT with a continuous ellipsoid zone (EZ band). (d–f) Patient P21 (16 years old) had a normal appearing fundus, reduced autofluorescence with subtle hyper-autofluorescence around the nasal fovea, and a grade 2 OCT with a disrupted EZ band. (g–i) Patient P19 (19 years old) presented a normal fundus, foveal and parafoveal hyper-autofluorescence, and a grade 3 OCT characterized by an absent EZ band. (j–l) Patient P6 (32 years old, homozygous for the c.847C>T variant in CNGA3) had a normal fundus appearance, foveal and parafoveal hyper-autofluorescence, and a grade 4 OCT characterized by the presence of a hypo-reflective zone. (m–o) Patient P10 (46 years old, homozygous for the c.1148delC variant in CNGB3) presented retinal pigment epithelium (RPE) dystrophy, a central region devoid of autofluorescence surrounded by a hyper-autofluorescent ring, and a grade 5 OCT characterized by outer retinal atrophy including RPE loss.
Figure 2. Ophthalmological findings in five representative ACHM patients showing the OCT grading (c,f,i,l,o) with corresponding retinography (a,d,g,j,m) and fundus AF images (b,e,h,k,n). (a–c) Patient P15 (12 years old, homozygous for the c.619G>A variant in GNAT2) presented a normal fundus appearance, normal fundus autofluorescence and a grade 1 OCT with a continuous ellipsoid zone (EZ band). (d–f) Patient P21 (16 years old) had a normal appearing fundus, reduced autofluorescence with subtle hyper-autofluorescence around the nasal fovea, and a grade 2 OCT with a disrupted EZ band. (g–i) Patient P19 (19 years old) presented a normal fundus, foveal and parafoveal hyper-autofluorescence, and a grade 3 OCT characterized by an absent EZ band. (j–l) Patient P6 (32 years old, homozygous for the c.847C>T variant in CNGA3) had a normal fundus appearance, foveal and parafoveal hyper-autofluorescence, and a grade 4 OCT characterized by the presence of a hypo-reflective zone. (m–o) Patient P10 (46 years old, homozygous for the c.1148delC variant in CNGB3) presented retinal pigment epithelium (RPE) dystrophy, a central region devoid of autofluorescence surrounded by a hyper-autofluorescent ring, and a grade 5 OCT characterized by outer retinal atrophy including RPE loss.
Ijms 22 01681 g002
Table
Table 1. Main clinical findings in each achromatopsia patient at the study baseline.
Table 1. Main clinical findings in each achromatopsia patient at the study baseline.
IDAge (yrs)ACHM SubtypeCong. Nystag. *Color Vision TestBCVA † (RE-LE)OCT Grade (RE-LE)CRT $ (RE-LE)Foveal Hypo-PlasiaFAF PatternMS $ (RE-LE)Fixation Stability (RE-LE)Dark-Adapted 0.01 ERGLight-Adapted 3.0 ERGLight-Adapted 30 Hz ERG
P17INCNoCB (I)0.700.7011219232Yesn/a18.318.7RSRSNormalReducedReduced
P29INCYesCB (F)0.701.0011265274NoReduced AF, subtle perifoveal hyper-AF19.619.7UNRSNormalReducedReduced
P329COMYesCB (F)1.001.0022163171Non/a13.213.9RSRSNormalUndetectableUndetectable
P418COMYesCB (F)1.001.0034167189NoWidespread reduced AF9.817.6SSNormalUndetectableUndetectable
P519COMYesCB (I)0.700.7011267253Non/a12.615.4RSRSNormalUndetectableUndetectable
P632COMYesCB (F)1.001.0044239228NoFoveal and parafoveal hyper-AF109.7RSRSNormalUndetectableUndetectable
P759COMYesCB (F)1.001.0044187233NoReduced AF, subtle perifoveal hyper-AF1010UNRSNormalUndetectableUndetectable
P826INCYesCB (F)0.700.7022235221YesWidespread reduced AF16.117.8UNRSNormalReducedReduced
P98INCYesCB (I)1.001.0011271267NoWidespread reduced AF1717.8UNRSNormalReducedReduced
P1046COMYesCB (F)1.301.3055125115NoAbsent AF at center, hyper-AF ring9.79.3RSRSNormalUndetectableUndetectable
P1110COMYesCB (F)1.000.7022244238Non/a16.816.5UNRSNormalUndetectableReduced
P124INCYesCB (F)1.001.0022227217Non/a18.518.3RSUNNormalReducedReduced
P137COMYesCB (I)1.001.0011304196NoFoveal and parafoveal hyper-AF9.38.5RSUNNormalUndetectableReduced
P145COMYesCB (I)1.001.0011258230Non/a14.415UNUNNormalUndetectableReduced
P152INCYesCB (I)0.400.4011260166NoNormal19.619.7RSRSNormalReducedReduced
P1615COMYesCB (F)0.700.7044196218NoFoveal and parafoveal hyper-AF14.514UNUNNormalUndetectableReduced
P177COMYesCB (I)1.301.3011229184Non/an/an/an/an/aNormalUndetectableReduced
P1815COMYesCB (I)1.001.0022183206NoNormal13.513.2RSRSNormalUndetectableReduced
P1915INCNoDeu, Tri (F) 0.150.3033182187NoFoveal and parafoveal hyper-AF19.418.7SRSNormalUndetectableReduced
P2031COMYesCB (F)1.001.0011220207NoNormaln/an/an/an/aNormalUndetectableUndetectable
P2116COMYesCB (F)0.700.7024271259YesReduced AF, subtle perifoveal hyper-AFn/an/an/an/aNormalUndetectableReduced
AF, autofluorescence; BCVA, best corrected visual acuity; CB, color blindness; COM, complete form; CRT, central retinal thickness; Deu, deuteranopia; ERG, electroretinogram; F, Farnsworth D-15 color test; I, Pseudoisochromatic Plate Ishihara; INC, incomplete form; LE, left eye; MS, macular sensitivity; n/a, not available; RE, right eye; RS, relatively stable; S, stable; Tri, tritanopia; UN, unstable; yrs, years. * All patients had photophobia. † BCVA expressed in logMAR. $ CRT and MS expressed in µm.
Table
Table 2. Genetic findings in the achromatopsia patients.
Table 2. Genetic findings in the achromatopsia patients.
IDFamilyGeneZygosityRefSeqAllele 1Allele 2
NucleotideProteinReferenceNucleotideProteinReference
P1F1CNGA3Hom.NM_001298c.1114C>Tp.Pro372Ser[28]c.1114C>Tp.Pro372Ser[28]
P2F2CNGA3Hom.NM_001298c.1641C>Ap.Phe547Leu[16]c.1641C>Ap.Phe547Leu[16]
P3F3CNGA3Hom.NM_001298c.1641C>Ap.Phe547Leu[16]c.1641C>Ap.Phe547Leu[16]
P4F4CNGA3Hom.NM_001298c.1162G>Ap.(Gly388Ser)This studyc.1162G>Ap.(Gly388Ser)This study
P5F5CNGA3C. Het.NM_001298c.667C>Tp.Arg223Trp[28]c.1060T>Cp.Ser341Pro[28]
P6F6CNGA3Hom.NM_001298c.847C>Tp.Arg283Trp[16]c.847C>Tp.Arg283Trp[16]
P7F6CNGA3Hom.NM_001298c.847C>Tp.Arg283Trp[16]c.847C>Tp.Arg283Trp[16]
P8F7CNGB3Hom.NM_019098c.1148delCp.(Thr383Ilefs*13)[29]c.1148delCp.(Thr383Ilefs*13)[29]
P9F8CNGB3Hom.NM_019098c.1148delCp.(Thr383Ilefs*13)[29]c.1148delCp.(Thr383Ilefs*13)[29]
P10F9CNGB3Hom.NM_019098c.1148delCp.(Thr383Ilefs*13)[29]c.1148delCp.(Thr383Ilefs*13)[29]
P11F10CNGB3C. Het.NM_019098c.1148delCp.(Thr383Ilefs*13)[29]c.970A>Gp.(Arg324Gly)[30]
P12F10CNGB3C. Het.NM_019098c.1148delCp.(Thr383Ilefs*13)[29]c.970A>Gp.(Arg324Gly)[30]
P13F11GNAT2Hom.NM_005272c.832dupp.(Ile278Asnfs*14)This studyc.832dupp.(Ile278Asnfs*14)This study
P14F11GNAT2Hom.NM_005272c.832dupp.(Ile278Asnfs*14)This studyc.832dupp.(Ile278Asnfs*14)This study
P15F12GNAT2Hom.NM_005272c.619G>Ap.(Glu207Lys)This studyc.619G>Ap.(Glu207Lys)This study
P16F13PDE6CHom.NM_006204 c.2017G>Tp.(Asp673Tyr)[30]c.2017G>Tp.(Asp673Tyr)[30]
P17F14CNGB3 ?Het.NM_019098c.1055G>Ap.(Arg352Lys)This studyn.d.n.d.n.d.
P18F15n.d.n.d.n.d.n.d.n.d.n.d.n.d.n.d.n.d.
C. Het., compound heterozygous; Het., heterozygous; Hom., homozygous; n.d., not determined; ?, uncertain.
Table
Table 3. Pathogenicity predictions for the novel missense variants reported in this study.
Table 3. Pathogenicity predictions for the novel missense variants reported in this study.
GeneRefSeqNucleotideProteinIn Silico Pathogenicity Analysis
MutationTaster †PolyPhen-2 ‡Cadd13 #
CNGA3NM_001298c.1162G>Ap.(Gly388Ser)Disease causingProb. damaging26.1
CNGB3NM_019098c.1055G>Ap.(Arg352Lys)Disease causingProb. damaging37.0
GNAT2NM_005272c.619G>A *p.(Glu207Lys)Disease causingProb. damaging29.4
†http://www.mutationtaster.org/. ‡ Polymorphism Phenotyping v2; http://genetics.bwh.harvard.edu/pph2/. # Combined Annotation-Dependent Depletion; http://cadd.gs.washington.edu/. * present in gnomAD.
Table
Table 4. Comparison of ACHM patients stratified according to the genotype.
Table 4. Comparison of ACHM patients stratified according to the genotype.
ParametersCNGA3CNGB3GNAT2p
(Seven Patients)(Five Patients)(Three Patients)
Age (years)24.7 ± 6.718.7 ± 7.74.64 ± 1.60.238
Age at diagnosis (years)5.4 ± 2.811.2 ± 18.44.6 ± 1.50.609
Detectable light-adapted 3.0 ERG 2 (28.6%)3 (60.0%)1 (33.3%)0.530
(n, %)
Detectable flicker ERG (n, %)2 (28.6%)4 (80.0%)3 (100%)0.057
Right eyeLeft eyeRight eyeLeft eyeRight eyeLeft eyep
BCVA (logMAR)0.87 ± 0.060.91 ± 0.061.00 ± 0.100.93 ± 0.110.80 ± 0.200.80 ± 0.20.493
CRT (µm)215.3 ± 16.6225.7 ± 13.4220.4 ± 25211.6 ± 25.7274 ± 15197.3 ± 18.50.803
MS (dB)13.4 ± 1.515 ± 1.515.6 ± 1.515.9 ± 1.714.4 ± 314.4 ± 3.20.371
Dark-adapted 0.01 ERG136.9 ± 11.9159.9 ± 21.1151.3 ± 24.2140.4 ± 27.6145 ± 27143.5 ± 23.50.515
(b-wave amplitude, µV)
Light adapted 3.0 ERG5.7 ± 3.79.4 ± 8.312 ± 8.714.3 ± 11.322.624.00.793
(b-wave amplitude, µV)
30 Hz Flicker ERG 3.9 ± 2.83.3 ± 2.13.7 ± 3.24.5 ± 3.45 ± 0.95.1 ± 0.50.103
(N1-P1, µV)
BCVA, best corrected visual acuity; CRT, central retinal thickness; D, diopter; ERG, electroretinogram; MS, macular sensitivity. Data are expressed as mean ± standard error of the mean for continuous variables, and as number of patients (frequency) for categorical variables.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brunetti-Pierri, R.; Karali, M.; Melillo, P.; Di Iorio, V.; De Benedictis, A.; Iaccarino, G.; Testa, F.; Banfi, S.; Simonelli, F. Clinical and Molecular Characterization of Achromatopsia Patients: A Longitudinal Study. Int. J. Mol. Sci. 2021, 22, 1681. https://doi.org/10.3390/ijms22041681

AMA Style

Brunetti-Pierri R, Karali M, Melillo P, Di Iorio V, De Benedictis A, Iaccarino G, Testa F, Banfi S, Simonelli F. Clinical and Molecular Characterization of Achromatopsia Patients: A Longitudinal Study. International Journal of Molecular Sciences. 2021; 22(4):1681. https://doi.org/10.3390/ijms22041681

Chicago/Turabian Style

Brunetti-Pierri, Raffaella, Marianthi Karali, Paolo Melillo, Valentina Di Iorio, Antonella De Benedictis, Gennarfrancesco Iaccarino, Francesco Testa, Sandro Banfi, and Francesca Simonelli. 2021. "Clinical and Molecular Characterization of Achromatopsia Patients: A Longitudinal Study" International Journal of Molecular Sciences 22, no. 4: 1681. https://doi.org/10.3390/ijms22041681

APA Style

Brunetti-Pierri, R., Karali, M., Melillo, P., Di Iorio, V., De Benedictis, A., Iaccarino, G., Testa, F., Banfi, S., & Simonelli, F. (2021). Clinical and Molecular Characterization of Achromatopsia Patients: A Longitudinal Study. International Journal of Molecular Sciences, 22(4), 1681. https://doi.org/10.3390/ijms22041681

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop