Comparison of the ABC and ACMG systems for variant classification

The ABC and ACMG variant classification systems were compared by asking mainly European clinical laboratories to classify variants in 10 challenging cases using both systems, and to state if the variant in question would be reported as a relevant result or not as a measure of clinical utility. In contrast to the ABC system, the ACMG system was not made to guide variant reporting but to determine the likelihood of pathogenicity. Nevertheless, this comparison is justified since the ACMG class determines variant reporting in many laboratories. Forty-three laboratories participated in the survey. In seven cases, the classification system used did not influence the reporting likelihood when variants labeled as "maybe report" after ACMG-based classification were included. In three cases of population frequent but disease-associated variants, there was a difference in favor of reporting after ABC classification. A possible reason is that ABC step C (standard variant comments) allows a variant to be reported in one clinical setting but not another, e.g., based on Bayesian-based likelihood calculation of clinical relevance. Finally, the selection of ACMG criteria was compared between 36 laboratories. When excluding criteria used by less than four laboratories (<10%), the average concordance rate was 46%. Taken together, ABC-based classification is more clear-cut than ACMG-based classification since molecular and clinical information is handled separately, and variant reporting can be adapted to the clinical question and phenotype. Furthermore, variants do not get a clinically inappropriate label, like pathogenic when not pathogenic in a clinical context, or variant of unknown significance when the significance is known.

The natural history and genotype-phenotype correlations of TMPRSS3 hearing loss: an international, multi-center, cohort analysis

TMPRSS3-related hearing loss presents challenges in correlating genotypic variants with clinical phenotypes due to the small sample sizes of previous studies. We conducted a cross-sectional genomics study coupled with retrospective clinical phenotype analysis on 127 individuals. These individuals were from 16 academic medical centers across 6 countries. Key findings revealed 47 unique TMPRSS3 variants with significant differences in hearing thresholds between those with missense variants versus those with loss-of-function genotypes. The hearing loss progression rate for the DFNB8 subtype was 0.3 dB/year. Post-cochlear implantation, an average word recognition score of 76% was observed. Of the 51 individuals with two missense variants, 10 had DFNB10 with profound hearing loss. These 10 all had at least one of 4 TMPRSS3 variants predicted by computational modeling to be damaging to TMPRSS3 structure and function. To our knowledge, this is the largest study of TMPRSS3 genotype-phenotype correlations. We find significant differences in hearing thresholds, hearing loss progression, and age of presentation, by TMPRSS3 genotype and protein domain affected. Most individuals with TMPRSS3 variants perform well on speech recognition tests after cochlear implant, however increased age at implant is associated with worse outcomes. These findings provide insight for genetic counseling and the on-going design of novel therapeutic approaches.

Stable long-term outcomes after cochlear implantation in subjects with TMPRSS3 associated hearing loss: a retrospective multicentre study

BACKGROUND

The spiral ganglion hypothesis suggests that pathogenic variants in genes preferentially expressed in the spiral ganglion nerves (SGN), may lead to poor cochlear implant (CI) performance. It was long thought that TMPRSS3 was particularly expressed in the SGNs. However, this is not in line with recent reviews evaluating CI performance in subjects with TMPRSS3-associated sensorineural hearing loss (SNHL) reporting overall beneficial outcomes. These outcomes are, however, based on variable follow-up times of, in general, 1 year or less. Therefore, we aimed to 1. evaluate long-term outcomes after CI implantation of speech recognition in quiet in subjects with TMPRSS3-associated SNHL, and 2. test the spiral ganglion hypothesis using the TMPRSS3-group.

METHODS

This retrospective, multicentre study evaluated long-term CI performance in a Dutch population with TMPRSS3-associated SNHL. The phoneme scores at 70 dB with CI in the TMPRSS3-group were compared to a control group of fully genotyped cochlear implant users with post-lingual SNHL without genes affecting the SGN, or severe anatomical inner ear malformations. CI-recipients with a phoneme score ≤ 70% at least 1-year post-implantation were considered poor performers and were evaluated in more detail.

RESULTS

The TMPRSS3 group consisted of 29 subjects (N = 33 ears), and the control group of 62 subjects (N = 67 ears). For the TMPRSS3-group, we found an average phoneme score of 89% after 5 years, which remained stable up to 10 years post-implantation. At both 5 and 10-year follow-up, no difference was found in speech recognition in quiet between both groups (p = 0.830 and p = 0.987, respectively). Despite these overall adequate CI outcomes, six CI recipients had a phoneme score of ≤ 70% and were considered poor performers. The latter was observed in subjects with residual hearing post-implantation or older age at implantation.

CONCLUSION

Subjects with TMPRSS3-associated SNHL have adequate and stable long-term outcomes after cochlear implantation, equal to the performance of genotyped patient with affected genes not expressed in the SGN. These findings are not in line with the spiral ganglion hypothesis. However, more recent studies showed that TMPRSS3 is mainly expressed in the hair cells with only limited SGN expression. Therefore, we cannot confirm nor refute the spiral ganglion hypothesis.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.

Optical genome mapping and revisiting short-read genome sequencing data reveal previously overlooked structural variants disrupting retinal disease-associated genes

PURPOSE

Structural variants (SVs) play an important role in inherited retinal diseases (IRD). Although the identification of SVs significantly improved upon the availability of genome sequencing, it is expected that involvement of SVs in IRDs is higher than anticipated. We revisited short-read genome sequencing data to enhance the identification of gene-disruptive SVs.

METHODS

Optical genome mapping was performed to improve SV detection in short-read genome sequencing-negative cases. In addition, reanalysis of short-read genome sequencing data was performed to improve the interpretation of SVs and to re-establish SV prioritization criteria.

RESULTS

In a monoallelic USH2A case, optical genome mapping identified a pericentric inversion (173 megabase), with 1 breakpoint disrupting USH2A. Retrospectively, the variant could be observed in genome sequencing data but was previously deemed false positive. Reanalysis of short-read genome sequencing data (427 IRD cases) was performed which yielded 30 pathogenic SVs affecting, among other genes, USH2A (n = 15), PRPF31 (n = 3), and EYS (n = 2). Eight of these (>25%) were overlooked during previous analyses.

CONCLUSION

Critical evaluation of our findings allowed us to re-establish and improve our SV prioritization and interpretation guidelines, which will prevent missing pathogenic events in future analyses. Our data suggest that more attention should be paid to SV interpretation and the current contribution of SVs in IRDs is still underestimated.

Genotype and Phenotype Analyses of a Novel WFS1 Variant (c.2512C>T p.(Pro838Ser)) Associated with DFNA6/14/38

The aim of this study is to contribute to a better description of the genotypic and phenotypic spectrum of DFNA6/14/38 and aid in counseling future patients identified with this variant. Therefore, we describe the genotype and phenotype in a large Dutch-German family (W21-1472) with autosomal dominant non-syndromic, low-frequency sensorineural hearing loss (LFSNHL). Exome sequencing and targeted analysis of a hearing impairment gene panel were used to genetically screen the proband. Co-segregation of the identified variant with hearing loss was assessed by Sanger sequencing. The phenotypic evaluation consisted of anamnesis, clinical questionnaires, physical examination and examination of audiovestibular function. A novel likely pathogenic WFS1 variant (NM_006005.3:c.2512C>T p.(Pro838Ser)) was identified in the proband and found to co-segregate with LFSNHL, characteristic of DFNA6/14/38, in this family. The self-reported age of onset of hearing loss (HL) ranged from congenital to 50 years of age. In the young subjects, HL was demonstrated in early childhood. At all ages, an LFSNHL (0.25-2 kHz) of about 50-60 decibel hearing level (dB HL) was observed. HL in the higher frequencies showed inter-individual variability. The dizziness handicap inventory (DHI) was completed by eight affected subjects and indicated a moderate handicap in two of them (aged 77 and 70). Vestibular examinations (n = 4) showed abnormalities, particularly in otolith function. In conclusion, we identified a novel WFS1 variant that co-segregates with DFNA6/14/38 in this family. We found indications of mild vestibular dysfunction, although it is uncertain whether this is related to the identified WFS1 variant or is an incidental finding. We would like to emphasize that conventional neonatal hearing screening programs are not sensitive to HL in DFNA6/14/38 patients, because high-frequency hearing thresholds are initially preserved. Therefore, we suggest screening newborns in DFNA6/14/38 families with more frequency-specific methods.

Diagnostic analysis of the highly complex OPN1LW/OPN1MW gene cluster using long-read sequencing and MLPA

Pathogenic variants in the OPN1LW/OPN1MW gene cluster are causal for a range of mild to severe visual impairments with color deficiencies. The widely utilized short-read next-generation sequencing (NGS) is inappropriate for the analysis of the OPN1LW/OPN1MW gene cluster and many patients with pathogenic variants stay underdiagnosed. A diagnostic genetic assay was developed for the OPN1LW/OPN1MW gene cluster, consisting of copy number analysis via multiplex ligation-dependent probe amplification and sequence analysis via long-read circular consensus sequencing. Performance was determined on 50 clinical samples referred for genetic confirmation of the clinical diagnosis (n = 43) or carrier status analysis (n = 7). A broad range of pathogenic haplotypes were detected, including deletions, hybrid genes, single variants and combinations of variants. The developed genetic assay for the OPN1LW/OPN1MW gene cluster is a diagnostic test that can detect both structural and nucleotide variants with a straightforward analysis, improving diagnostic care of patients with visual impairment.

Identification of a Complex Allele in IMPG2 as a Cause of Adult-Onset Vitelliform Macular Dystrophy

Purpose

Inherited retinal diseases are a group of clinically and genetically heterogeneous disorders with approximately 270 genes involved. IMPG2 is associated with adult-onset vitelliform macular dystrophy. Here, we investigated two unrelated patients with vitelliform macular dystrophy to identify the underlying genetic cause.

Methods

Whole-exome sequencing identified a putative causal complex allele consisting of c.3023-15T>A and c.3023G>A (p.(Gly1008Asp)) in IMPG2 in both individuals. To assess its effect, in vitro splice assays in HEK293T and further characterization in patient-derived photoreceptor precursor cells (PPCs) were conducted.

Results

The results of the midigene splice assays in HEK293T showed that the complex allele causes a variety of splicing defects ranging from a small deletion to (multiple-)exon skipping. This finding was further validated using patient-derived PPCs that showed a significant increase of out-of-frame transcripts lacking one or multiple exons compared to control-derived PPCs. Overall, control PPCs consistently showed low levels of aberrantly spliced IMPG2 transcripts that were highly elevated in patient-derived PPCs. These differences were even more obvious upon inhibition of nonsense-mediated decay with cycloheximide.

Conclusions

We report a heterozygous complex allele in IMPG2 causative for adult-onset vitelliform macular dystrophy in two unrelated individuals with mild visual loss and bilateral vitelliform lesions. The predicted causal missense mutation c.3023G>A, located in the consensus splice acceptor site, enhances the splicing effect of the upstream variant c.3023-15T>A, leading to the generation of aberrant transcripts that decrease the full-length IMPG2 levels. These results suggest a haploinsufficiency mechanism of action and highlight the complementarity of using different models to functionally assesses splicing defects.

Personalized genetic counseling for Stargardt disease: Offspring risk estimates based on variant severity

Recurrence risk calculations in autosomal recessive diseases are complicated when the effect of genetic variants and their population frequencies and penetrances are unknown. An example of this is Stargardt disease (STGD1), a frequent recessive retinal disease caused by bi-allelic pathogenic variants in ABCA4. In this cross-sectional study, 1,619 ABCA4 variants from 5,579 individuals with STGD1 were collected and categorized by (1) severity based on statistical comparisons of their frequencies in STGD1-affected individuals versus the general population, (2) their observed versus expected homozygous occurrence in STGD1-affected individuals, (3) their occurrence in combination with established mild alleles in STGD1-affected individuals, and (4) previous functional and clinical studies. We used the sum allele frequencies of these severity categories to estimate recurrence risks for offspring of STGD1-affected individuals and carriers of pathogenic ABCA4 variants. The risk for offspring of an STGD1-affected individual with the "severe|severe" genotype or a "severe|mild with complete penetrance" genotype to develop STGD1 at some moment in life was estimated at 2.8%-3.1% (1 in 36-32 individuals) and 1.6%-1.8% (1 in 62-57 individuals), respectively. The risk to develop STGD1 in childhood was estimated to be 2- to 4-fold lower: 0.68%-0.79% (1 in 148-126) and 0.34%-0.39% (1 in 296-252), respectively. In conclusion, we established personalized recurrence risk calculations for STGD1-affected individuals with different combinations of variants. We thus propose an expanded genotype-based personalized counseling to appreciate the variable recurrence risks for STGD1-affected individuals. This represents a conceptual breakthrough because risk calculations for STGD1 may be exemplary for many other inherited diseases.

Usher syndrome type IV: clinically and molecularly confirmed by novel ARSG variants

Usher syndrome (USH) is an autosomal recessively inherited disease characterized by sensorineural hearing loss (SNHL) and retinitis pigmentosa (RP) with or without vestibular dysfunction. It is highly heterogeneous both clinically and genetically. Recently, variants in the arylsulfatase G (ARSG) gene have been reported to underlie USH type IV. This distinct type of USH is characterized by late-onset RP with predominantly pericentral and macular changes, and late onset SNHL without vestibular dysfunction. In this study, we describe the USH type IV phenotype in three unrelated subjects. We identified three novel pathogenic variants, two novel likely pathogenic variants, and one previously described pathogenic variant in ARSG. Functional experiments indicated a loss of sulfatase activity of the mutant proteins. Our findings confirm that ARSG variants cause the newly defined USH type IV and support the proposed extension of the phenotypic USH classification.

Whole genome sequencing and in vitro splice assays reveal genetic causes for inherited retinal diseases

Inherited retinal diseases (IRDs) are a major cause of visual impairment. These clinically heterogeneous disorders are caused by pathogenic variants in more than 270 genes. As 30-40% of cases remain genetically unexplained following conventional genetic testing, we aimed to obtain a genetic diagnosis in an IRD cohort in which the genetic cause was not found using whole-exome sequencing or targeted capture sequencing. We performed whole-genome sequencing (WGS) to identify causative variants in 100 unresolved cases. After initial prioritization, we performed an in-depth interrogation of all noncoding and structural variants in genes when one candidate variant was detected. In addition, functional analysis of putative splice-altering variants was performed using in vitro splice assays. We identified the genetic cause of the disease in 24 patients. Causative coding variants were observed in genes such as ATXN7, CEP78, EYS, FAM161A, and HGSNAT. Gene disrupting structural variants were also detected in ATXN7, PRPF31, and RPGRIP1. In 14 monoallelic cases, we prioritized candidate noncanonical splice sites or deep-intronic variants that were predicted to disrupt the splicing process based on in silico analyses. Of these, seven cases were resolved as they carried pathogenic splice defects. WGS is a powerful tool to identify causative variants residing outside coding regions or heterozygous structural variants. This approach was most efficient in cases with a distinct clinical diagnosis. In addition, in vitro splice assays provide important evidence of the pathogenicity of rare variants.

PRPH2 mutation update: In silico assessment of 245 reported and 7 novel variants in patients with retinal disease

Mutations in PRPH2, encoding peripherin-2, are associated with the development of a wide variety of inherited retinal diseases (IRDs). To determine the causality of the many PRPH2 variants that have been discovered over the last decades, we surveyed all published PRPH2 variants up to July 2020, describing 720 index patients that in total carried 245 unique variants. In addition, we identified seven novel PRPH2 variants in eight additional index patients. The pathogenicity of all variants was determined using the ACMG guidelines. With this, 107 variants were classified as pathogenic, 92 as likely pathogenic, one as benign, and two as likely benign. The remaining 50 variants were classified as variants of uncertain significance. Interestingly, of the total 252 PRPH2 variants, more than half (n = 137) were missense variants. All variants were uploaded into the Leiden Open source Variation and ClinVar databases. Our study underscores the need for experimental assays for variants of unknown significance to improve pathogenicity classification, which would allow us to better understand genotype-phenotype correlations, and in the long-term, hopefully also support the development of therapeutic strategies for patients with PRPH2-associated IRD.

Development of a Genotype Assay for Age-Related Macular Degeneration: The EYE-RISK Consortium

PURPOSE

To develop a genotype assay to assess associations with common and rare age-related macular degeneration (AMD) risk variants, to calculate an overall genetic risk score (GRS), and to identify potential misdiagnoses with inherited macular dystrophies that mimic AMD.

DESIGN

Case-control study.

PARTICIPANTS

Individuals (n = 4740) from 5 European cohorts.

METHODS

We designed single-molecule molecular inversion probes for target selection and used next generation sequencing to sequence 87 single nucleotide polymorphisms (SNPs), coding and splice-site regions of 10 AMD-(related) genes (ARMS2, C3, C9, CD46, CFB, CFH, CFI, HTRA1, TIMP3, and SLC16A8), and 3 genes that cause inherited macular dystrophies (ABCA4, CTNNA1, and PRPH2). Genetic risk scores for common AMD risk variants were calculated based on effect size and genotype of 52 AMD-associated variants. Frequency of rare variants was compared between late AMD patients and control individuals with logistic regression analysis.

MAIN OUTCOME MEASURES

Genetic risk score, association of genetic variants with AMD, and genotype-phenotype correlations.

RESULTS

We observed high concordance rates between our platform and other genotyping platforms for the 69 successfully genotyped SNPs (>96%) and for the rare variants (>99%). We observed a higher GRS for patients with late AMD compared with patients with early/intermediate AMD (P < 0.001) and individuals without AMD (P < 0.001). A higher proportion of pathogenic variants in the CFH (odds ratio [OR] = 2.88; P = 0.006), CFI (OR = 4.45; P = 0.005), and C3 (OR = 6.56; P = 0.0003) genes was observed in late AMD patients compared with control individuals. In 9 patients, we identified pathogenic variants in the PRPH2, ABCA4, and CTNNA1 genes, which allowed reclassification of these patients as having inherited macular dystrophy.

CONCLUSIONS

This study reports a genotype assay for common and rare AMD genetic variants, which can identify individuals at intermediate to high genetic risk of late AMD and enables differential diagnosis of AMD-mimicking dystrophies. Our study supports sequencing of CFH, CFI, and C3 genes because they harbor rare high-risk variants. Carriers of these variants could be amendable for new treatments for AMD that currently are under development.

Extending the Spectrum of EYS-Associated Retinal Disease to Macular Dystrophy

Purpose

To assess the phenotypic variability and natural course of inherited retinal diseases (IRDs) caused by EYS mutations.

Methods

Multiethnic cohort study (N = 30) with biallelic EYS variants from a clinical IRD database (retinitis pigmentosa [RP], N = 27; cone-rod dystrophy [CRD], N = 1; and macular dystrophy, N = 2). In vitro minigene splice assay was performed to determine the effect on EYS pre-mRNA splicing of the c.1299+5_1299+8del variant in macular dystrophy patients.

Results

We found 27 different EYS variants in RP patients and 7 were novel. The rate of visual field loss of the V4e isopter area was -0.84 ± 0.44 ln(deg2) per year, and the rate of visual acuity loss was 0.75 Early Treatment Diabetic Retinopathy Study letters per year. Ellipsoid zone width was correlated with area of the hyperautofluorescent ring, with rs = 0.78 and P < 0.001. Rate of decline in ellipsoid zone width was -57 ± 17 μm per year (P < 0.01) (n = 14) or -3.69% ± 0.51% from baseline per year (P < 0.001). An isolated CRD patient carried a homozygous EYS variant (c.9405T>A), previously identified in RP patients. Two siblings with macular dystrophy carried compound heterozygous EYS variants: c.1299+5_1299+8del and c.6050G>T. The former was novel and shown to result in skipping of exon 8, and the latter was a known RP variant.

Conclusions

We report on EYS-associated macular dystrophy, extending the spectrum of EYS-associated IRDs. We observed heterogeneity between RP patients in age of onset and disease progression. Identical EYS variants were found in cases with RP, CRD, and macular dystrophy. Screening for EYS variants in CRD and macular dystrophy patients might increase the diagnostic yield in previously unsolved cases.

Robust and accurate detection and sizing of repeats within the DMPK gene using a novel TP-PCR test

Myotonic dystrophy type 1 is a multisystem disorder caused by the expansion of a trinucleotide repeat in the DMPK gene. In this study we evaluated the performance of the FastDM1^TMDMPK sizing kit in myotonic dystrophy type 1 testing. This commercially available triplet repeat-primed PCR based kit was validated using reference and clinical samples. Based on testing with 19 reference samples, the assay yielded repeat sizes within three repeats from the consensus reference length, demonstrating an accuracy of 100%. Additionally, the assay generated consistent repeat size information with a concentration range of template-DNA, and upon repetition and reproduction (CV 0.36% to 0.41%). Clinical performance was established with 235 archived prenatal and postnatal clinical samples, yielding results of 100% sensitivity (95% CI, 97.29% to 100%) and 100% specificity (95% CI, 96.19% to 100%) in classifying the samples into the respective genotype groups of 5–35 (normal), 36–50 (non-pathogenic pre-expansion), 51–150 (unstable intermediate-sized pathogenic) or >150 (unstable pathogenic) CTG repeats, respectively. Furthermore, the assay identified interrupted repeat expansions in all samples known to have interruptions, and also identified interruptions in a subset of the clinical samples.

Mutations in the polyglutamylase gene TTLL5, expressed in photoreceptor cells and spermatozoa, are associated with cone-rod degeneration and reduced male fertility

Hereditary retinal degenerations encompass a group of genetic diseases characterized by extreme clinical variability. Following next-generation sequencing and autozygome-based screening of patients presenting with a peculiar, recessive form of cone-dominated retinopathy, we identified five homozygous variants [p.(Asp594fs), p.(Gln117*), p.(Met712fs), p.(Ile756Phe), and p.(Glu543Lys)] in the polyglutamylase-encoding gene TTLL5, in eight patients from six families. The two male patients carrying truncating TTLL5 variants also displayed a substantial reduction in sperm motility and infertility, whereas those carrying missense changes were fertile. Defects in this polyglutamylase in humans have recently been associated with cone photoreceptor dystrophy, while mouse models carrying truncating mutations in the same gene also display reduced fertility in male animals. We examined the expression levels of TTLL5 in various human tissues and determined that this gene has multiple viable isoforms, being highly expressed in testis and retina. In addition, antibodies against TTLL5 stained the basal body of photoreceptor cells in rat and the centrosome of the spermatozoon flagellum in humans, suggesting a common mechanism of action in these two cell types. Taken together, our data indicate that mutations in TTLL5 delineate a novel, allele-specific syndrome causing defects in two as yet pathogenically unrelated functions, reproduction and vision.

Identification of Inherited Retinal Disease-Associated Genetic Variants in 11 Candidate Genes

Inherited retinal diseases (IRDs) display an enormous genetic heterogeneity. Whole exome sequencing (WES) recently identified genes that were mutated in a small proportion of IRD cases. Consequently, finding a second case or family carrying pathogenic variants in the same candidate gene often is challenging. In this study, we searched for novel candidate IRD gene-associated variants in isolated IRD families, assessed their causality, and searched for novel genotype-phenotype correlations. Whole exome sequencing was performed in 11 probands affected with IRDs. Homozygosity mapping data was available for five cases. Variants with minor allele frequencies ≤ 0.5% in public databases were selected as candidate disease-causing variants. These variants were ranked based on their: (a) presence in a gene that was previously implicated in IRD; (b) minor allele frequency in the Exome Aggregation Consortium database (ExAC); (c) in silico pathogenicity assessment using the combined annotation dependent depletion (CADD) score; and (d) interaction of the corresponding protein with known IRD-associated proteins. Twelve unique variants were found in 11 different genes in 11 IRD probands. Novel autosomal recessive and dominant inheritance patterns were found for variants in Small Nuclear Ribonucleoprotein U5 Subunit 200 (SNRNP200) and Zinc Finger Protein 513 (ZNF513), respectively. Using our pathogenicity assessment, a variant in DEAH-Box Helicase 32 (DHX32) was the top ranked novel candidate gene to be associated with IRDs, followed by eight medium and lower ranked candidate genes. The identification of candidate disease-associated sequence variants in 11 single families underscores the notion that the previously identified IRD-associated genes collectively carry > 90% of the defects implicated in IRDs. To identify multiple patients or families with variants in the same gene and thereby provide extra proof for pathogenicity, worldwide data sharing is needed.

EYS mutation update: In silico assessment of 271 reported and 26 novel variants in patients with retinitis pigmentosa

Mutations in Eyes shut homolog (EYS) are one of the most common causes of autosomal recessive (ar) retinitis pigmentosa (RP), a progressive blinding disorder. The exact function of the EYS protein and the pathogenic mechanisms underlying EYS-associated RP are still poorly understood, which hampers the interpretation of the causality of many EYS variants discovered to date. We collected all reported EYS variants present in 377 arRP index cases published before June 2017, and uploaded them in the Leiden Open Variation Database (www.LOVD.nl/EYS). We also describe 36 additional index cases, carrying 26 novel variants. Of the 297 unique EYS variants identified, almost half (n = 130) are predicted to result in premature truncation of the EYS protein. Classification of all variants using the American College of Medical Genetics and Genomics guidelines revealed that the predicted pathogenicity of these variants cover the complete spectrum ranging from likely benign to pathogenic, although especially missense variants largely fall in the category of uncertain significance. Besides the identification of likely benign alleles previously reported as being probably pathogenic, our comprehensive analysis underscores the need of functional assays to assess the causality of EYS variants, in order to improve molecular diagnostics and counseling of patients with EYS-associated RP.

Missense mutations in the WD40 domain of AHI1 cause non-syndromic retinitis pigmentosa

Background

Recent findings suggesting that Abelson helper integration site 1 (AHI1) is involved in non-syndromic retinal disease have been debated, as the functional significance of identified missense variants was uncertain. We assessed whether AHI1 variants cause non-syndromic retinitis pigmentosa (RP).

Methods

Exome sequencing was performed in three probands with RP. The effects of the identified missense variants in AHI1 were predicted by three-dimensional structure homology modelling. Ciliary parameters were evaluated in patient’s fibroblasts, and recombinant mutant proteins were expressed in ciliated retinal pigmented epithelium cells.

Results

In the three patients with RP, three sets of compound heterozygous variants were detected in AHI1 (c.2174G>A; p.Trp725* and c.2258A>T; p.Asp753Val, c.660delC; p.Ser221Glnfs*10 and c.2090C>T; p.Pro697Leu, c.2087A>G; p.His696Arg and c.2429C>T; p.Pro810Leu). All four missense variants were present in the conserved WD40 domain of Jouberin, the ciliary protein encoded by AHI1, with variable predicted implications for the domain structure. No significant changes in the percentage of ciliated cells, nor in cilium length or intraflagellar transport were detected. However, expression of mutant recombinant Jouberin in ciliated cells showed a significantly decreased enrichment at the ciliary base.

Conclusions

This report confirms that mutations in AHI1 can underlie autosomal recessive RP. Moreover, it structurally and functionally validates the effect of the RP-associated AHI1 variants on protein function, thus proposing a new genotype–phenotype correlation for AHI1 mutation associated retinal ciliopathies.

Whole-Exome Sequencing Identifies Biallelic IDH3A Variants as a Cause of Retinitis Pigmentosa Accompanied by Pseudocoloboma

PURPOSE

To identify the genetic cause of and describe the phenotype in 4 families with autosomal recessive retinitis pigmentosa (arRP) that can be associated with pseudocoloboma.

DESIGN

Case series.

PARTICIPANTS

Seven patients from 4 unrelated families with arRP, among whom 3 patients had bilateral early-onset macular pseudocoloboma.

METHODS

We performed homozygosity mapping and whole-exome sequencing in 5 probands and 2 unaffected family members from 4 unrelated families. Subsequently, Sanger sequencing and segregation analysis were performed in additional family members. We reviewed the medical history of individuals carrying IDH3A variants and performed additional ophthalmic examinations, including full-field electroretinography, fundus photography, fundus autofluorescence imaging, and optical coherence tomography.

MAIN OUTCOME MEASURES

IDH3A variants, age at diagnosis, visual acuity, fundus appearance, visual field, and full-field electroretinography, fundus autofluorescence, and optical coherence tomography findings.

RESULTS

We identified 7 different variants in IDH3A in 4 unrelated families, that is, 5 missense, 1 nonsense, and 1 frameshift variant. All participants showed symptoms early in life, ranging from night blindness to decreased visual acuity, and were diagnosed between the ages of 1 and 11 years. Four participants with biallelic IDH3A variants displayed a typical arRP phenotype and 3 participants were diagnosed with arRP and pseudocoloboma of the macula.

CONCLUSIONS

IDH3A variants were identified as a novel cause of typical arRP in some individuals associated with macular pseudocoloboma. We observed both phenotypes in 2 siblings carrying the same compound heterozygous variants, which could be explained by variable disease expression and warrants caution when making assertions about genotype-phenotype correlations.

Diagnostic exome sequencing in 266 Dutch patients with visual impairment

Inherited eye disorders have a large clinical and genetic heterogeneity, which makes genetic diagnosis cumbersome. An exome-sequencing approach was developed in which data analysis was divided into two steps: the vision gene panel and exome analysis. In the vision gene panel analysis, variants in genes known to cause inherited eye disorders were assessed for pathogenicity. If no causative variants were detected and when the patient consented, the entire exome data was analyzed. A total of 266 Dutch patients with different types of inherited eye disorders, including inherited retinal dystrophies, cataract, developmental eye disorders and optic atrophy, were investigated. In the vision gene panel analysis (likely), causative variants were detected in 49% and in the exome analysis in an additional 2% of the patients. The highest detection rate of (likely) causative variants was in patients with inherited retinal dystrophies, for instance a yield of 63% in patients with retinitis pigmentosa. In patients with developmental eye defects, cataract and optic atrophy, the detection rate was 50, 33 and 17%, respectively. An exome-sequencing approach enables a genetic diagnosis in patients with different types of inherited eye disorders using one test. The exome approach has the same detection rate as targeted panel sequencing tests, but offers a number of advantages. For instance, the vision gene panel can be frequently and easily updated with additional (novel) eye disorder genes. Determination of the genetic diagnosis improved the clinical diagnosis, regarding the assessment of the inheritance pattern as well as future disease perspective.

Frequency and characterization of known and novel RHD variant alleles in 37 782 Dutch D-negative pregnant women

To guide anti-D prophylaxis, Dutch D- pregnant women are offered a quantitative fetal-RHD-genotyping assay to determine the RHD status of their fetus. This allowed us to determine the frequency of different maternal RHD variants in 37 782 serologically D- pregnant women. A variant allele is present in at least 0·96% of Dutch D- pregnant women The D- serology could be confirmed after further serological testing in only 54% of these women, which emphasizes the potential relevance of genotyping of blood donors. 43 different RHD variant alleles were detected, including 15 novel alleles (11 null-, 2 partial D- and 2 DEL-alleles). Of those novel null alleles, one allele contained a single missense mutation (RHD*443C>G) and one allele had a single amino acid deletion (RHD*424_426del). The D- phenotype was confirmed by transduction of human D- erythroblasts, consolidating that, for the first time, a single amino acid change or deletion causes the D- phenotype. Transduction also confirmed the phenotypes for the two new variant DEL-alleles (RHD*721A>C and RHD*884T>C) and the novel partial RHD*492C>A allele. Notably, in three additional cases the DEL phenotype was observed but sequencing of the coding sequence, flanking introns and promoter region revealed an apparently wild-type RHD allele without mutations.

Non-syndromic retinitis pigmentosa due to mutations in the mucopolysaccharidosis type IIIC gene, heparan-alpha-glucosaminide N-acetyltransferase (HGSNAT)

Retinitis pigmentosa (RP), the most common form of inherited retinal degeneration, is clinically and genetically heterogeneous and can appear as syndromic or non-syndromic. Mucopolysaccharidosis type IIIC (MPS IIIC) is a lethal disorder, caused by mutations in the heparan-alpha-glucosaminide N-acetyltransferase (HGSNAT) gene and characterized by progressive neurological deterioration, with retinal degeneration as a prominent feature. We identified HGSNAT mutations in six patients with non-syndromic RP. Whole exome sequencing (WES) in an Ashkenazi Jewish Israeli RP patient revealed a novel homozygous HGSNAT variant, c.370A>T, which leads to partial skipping of exon 3. Screening of 66 Ashkenazi RP index cases revealed an additional family with two siblings homozygous for c.370A>T. WES in three Dutch siblings with RP revealed a complex HGSNAT variant, c.[398G>C; 1843G>A] on one allele, and c.1843G>A on the other allele. HGSNAT activity levels in blood leukocytes of patients were reduced compared with healthy controls, but usually higher than those in MPS IIIC patients. All patients were diagnosed with non-syndromic RP and did not exhibit neurological deterioration, or any phenotypic features consistent with MPS IIIC. Furthermore, four of the patients were over 60 years old, exceeding by far the life expectancy of MPS IIIC patients. HGSNAT is highly expressed in the mouse retina, and we hypothesize that the retina requires higher HGSNAT activity to maintain proper function, compared with other tissues associated with MPS IIIC, such as the brain. This report broadens the spectrum of phenotypes associated with HGSNAT mutations and highlights the critical function of HGSNAT in the human retina.

Homozygosity mapping and targeted sanger sequencing reveal genetic defects underlying inherited retinal disease in families from pakistan

Background

Homozygosity mapping has facilitated the identification of the genetic causes underlying inherited diseases, particularly in consanguineous families with multiple affected individuals. This knowledge has also resulted in a mutation dataset that can be used in a cost and time effective manner to screen frequent population-specific genetic variations associated with diseases such as inherited retinal disease (IRD).

Methods

We genetically screened 13 families from a cohort of 81 Pakistani IRD families diagnosed with Leber congenital amaurosis (LCA), retinitis pigmentosa (RP), congenital stationary night blindness (CSNB), or cone dystrophy (CD). We employed genome-wide single nucleotide polymorphism (SNP) array analysis to identify homozygous regions shared by affected individuals and performed Sanger sequencing of IRD-associated genes located in the sizeable homozygous regions. In addition, based on population specific mutation data we performed targeted Sanger sequencing (TSS) of frequent variants in AIPL1, CEP290, CRB1, GUCY2D, LCA5, RPGRIP1 and TULP1, in probands from 28 LCA families.

Results

Homozygosity mapping and Sanger sequencing of IRD-associated genes revealed the underlying mutations in 10 families. TSS revealed causative variants in three families. In these 13 families four novel mutations were identified in CNGA1, CNGB1, GUCY2D, and RPGRIP1.

Conclusions

Homozygosity mapping and TSS revealed the underlying genetic cause in 13 IRD families, which is useful for genetic counseling as well as therapeutic interventions that are likely to become available in the near future.

Impact of genetic variation in the SMIM1 gene on Vel expression levels

BACKGROUND

Serologic determination of the Vel- phenotype is challenging due to variable Vel expression levels. In this study we investigated the genetic basis for weak Vel expression levels and developed a high-throughput genotyping assay to detect Vel- donors.

STUDY DESIGN AND METHODS

In 548 random Caucasian and 107 Vel+(w) donors genetic variation in the SMIM1 gene was studied and correlated to Vel expression levels. A total of 3366 Caucasian, 621 black, and 333 Chinese donors were screened with a high-throughput genotyping assay targeting the SMIM1*64_80del allele.

RESULTS

The Vel+(w) phenotype is in most cases caused by the presence of one SMIM1 allele carrying the major allele of the rs1175550 SNP in combination with a SMIM1*64_80del allele or in few cases caused by the presence of the SMIM1*152T>A or SMIM1*152T>G allele. In approximately 6% of Vel+(w) donors genetic factors in SMIM1 could not explain the weak expression. We excluded the possibility that lack of expression of another blood group system was correlated with weak Vel expression levels. Furthermore, using a high-throughput Vel genotyping assay we detected two Caucasian Vel- donors.

CONCLUSION

Weak Vel expression levels are caused by multiple genetic factors in SMIM1 and probably also by other genetic or environmental factors. Due to the variation in Vel expression levels, serologic determination of the Vel- phenotype is difficult and a genotyping assay targeting the c.64_80del deletion in SMIM1 should be used to screen donors for the Vel- phenotype.

Novel alleles at the Kell blood group locus that lead to Kell variant phenotype in the Dutch population

BACKGROUND

Alloantibodies directed against antigens of the Kell blood group system are clinically significant. In the Netherlands, the KEL1 antigen is determined in all blood donors. In this study, after phenotyping of KEL:1-positive donors, genotyping analysis was conducted in KEL:1,-2 donors to identify possible KEL*02 variant alleles.

STUDY DESIGN AND METHODS

A total of 407 donors with the KEL:1,-2 phenotype were genotyped for the KEL*01/02 polymorphism, followed by direct sequencing of the KEL gene if the KEL*02 allele was detected. Two K0 patients were also included. Transcript analysis was conducted in two probands with the KEL*02. M05 allele defined by a synonymous mutation (G573G). Flow cytometry analysis to determine the expression of Kell antigen was performed.

RESULTS

Thirty KEL:1,-2 individuals (30/407, 7.4%) with discrepant KEL*01/02 genotype were identified. Seven novel alleles were identified: KEL*02(R86Q, R281W)mod, KEL*02(L133P)null, KEL*02(436delG)null, KEL*02(F418S)null, KEL*02(R492X)null, KEL*02(L611R)null, and KEL*02(R700X)null. Nine variant alleles described before were detected: KEL*02N.06, KEL*02N.15, KEL*02N.17, KEL*02N.19, KEL*02N.21, KEL*02M.02, KEL*02M.04, KEL*02M.05, and KEL*02(Q362K)mod. A transcript lacking Exon 16 was identified in two probands with the KEL*02M.05 allele as described before. Finally, flow cytometry analysis showed a decreased total Kell expression and a relatively increased KEL1 expression in individuals with the KEL:1,2null or KEL:1,2mod phenotype, compared to KEL:1,2 controls.

CONCLUSION

In 7.4% of a group of tested KEL:1,-2 Dutch donors, a KEL*02null or KEL*02mod allele was found. A relatively increased KEL1 antigen expression in KEL:1,2null and KEL:1,2mod individuals suggest that the expression of Kell-XK complexes depends on the availability of the XK protein.

Molecular analysis of immunized Jr(a-) or Lan- patients and validation of a high-throughput genotyping assay to screen blood donors for Jr(a-) and Lan- phenotypes

BACKGROUND

Individuals with anti-Jr(a) or anti-Lan are ideally transfused with rare Jr(a-) or Lan- red blood cells. We characterized mutations in Dutch Jr(a-) and Lan- individuals and developed a high-throughput genotyping assay to detect Jr(a-) and Lan- donors.

STUDY DESIGN AND METHODS

Six Jr(a-) and seven Lan- persons, who all made anti-Jr(a) or anti-Lan, were sequenced for ABCG2 or ABCB6 and the copy number of ABCG2 and ABCB6 was determined. A total of 3366 Caucasian, 621 black, and 333 Chinese donors were screened with a high-throughput screening assay targeting frequently occurring mutations causing the Jr(a-) or Lan- phenotype.

RESULTS

In the six tested Jr(a-) individuals previously described, c.376C > T, c.706C > T, and c.736C > T nonsense mutations in ABCG2 were detected. In the seven Lan- individuals 12 different mutations, of which 10 underlie the Lan- phenotype, were detected. No copy number variation was detected for ABCG2 and ABCB6. The high-throughput screening assay detected five Caucasian donors heterozygous for the c.706C > T or 736C > T mutation in ABCG2 and nine Caucasian donors heterozygous for the 574C > T mutation in ABCB6. No black or Chinese donors were found positive for a mutation.

CONCLUSION

We describe eight new mutations in ABCB6 of which seven, including three missense mutations, underlie the Lan- phenotype and determine that a complete gene deletion of ABCG2 or ABCB6 is not responsible for the Jr(a-) or Lan- phenotype, respectively. The extended heterogeneity of mutations causing the Jr(a-) or Lan- phenotype in most populations makes genetic screening for the Jr(a-) and Lan- phenotype inefficient in those populations.

Comprehensive genotyping for 18 blood group systems using a multiplex ligation-dependent probe amplification assay shows a high degree of accuracy

BACKGROUND

In recent years genotyping methods have been implemented in blood banks as alternative to comprehensive serologic typing. We evaluated a newly developed assay for convenient and comprehensive genotyping of blood group alleles based on multiplex ligation-dependent probe amplification (MLPA) technology.

STUDY DESIGN AND METHODS

We analyzed 103 random and 150 selected samples to validate the specificity of the blood-MLPA assay that is able to determine the presence, absence, and copy number of 48 blood group and 112 variant alleles of 18 blood group systems. A total of 4038 serologic typing results, including 52 different antigens, were available for these samples.

RESULTS

In 4018 (99.5%) of the 4038 serologic typing results the predicted phenotypes by the blood-MLPA were in concordance with serologic typing. Twenty discordant results were due to false-positive serologic results (n = 2), false-negative serologic results (n = 1), inability of routine serologic typing to detect variant antigens (n = 14), or false-positive prediction from the blood-MLPA due to the presence of a null allele (n = 3).

CONCLUSION

The blood-MLPA reliably predicts the presence or absence of blood group antigens, including almost all clinically relevant blood group antigens, except ABO, in patients and donors. Furthermore, it is the first assay that determines copy numbers of blood group alleles in the same test. It even provides more detailed and accurate information than serologic typing, because most variant alleles are immediately recognized. Since only standard laboratory equipment is needed, this assay finally offers the possibility to comprehensively type recipients and makes extensive matching for selected patients groups more feasible.

Serologic and molecular investigations of DAR1 (weak D Type 4.2), DAR1.2, DAR1.3, DAR2 (DARE), and DARA

BACKGROUND

The partial D variant DAR1 (weak D Type 4.2) is caused by three single-point mutations, 602C>G, 667T>G, and 1025T>C. Here we report a molecular study on different D variants belonging to the DAR category (DAR1, DAR1.2, DAR1.3, and DAR2) and their serologic data.

STUDY DESIGN AND METHODS

A total of 42 samples belonging to the DAR category were screened for the presence of the silent mutations 744C>T and 957G>A. The samples were phenotyped for RhD and RhCE, characterized for RhD epitope expression, and sequenced for RHD exons. Flow cytometry was performed to determine RhD antigen density.

RESULTS

The silent mutation 744C>T was found in all six samples previously typed as RHD*DAR2 (602C>G, 667T>G, 957G>A, 1025T>C). In addition to the three nucleotide changes originally reported for the RHD*DAR1 allele, the silent mutations 744C>T and 957G>A were found in 14 of 16 samples previously typed as RHD*DAR1. In the remaining two samples one additional silent mutation, 744C>T, was found. Serologically the DAR1.2 and DAR1.3 samples analyzed in this study showed no distinct difference in their anti-D reaction pattern compared to each other. The anti-D reaction pattern of DARA/DAR2 showed some distinct differences compared to those of DAR1.2 and DAR1.3.

CONCLUSION

RHD*DARA and RHD*DAR2 are the same allele. Furthermore, the alleles RHD*DAR1.2 and RHD*DAR1.3 both exist; however, the silent mutation 957G>A (V319) showed no influence on the RhD phenotype.

SMIM1 underlies the Vel blood group and influences red blood cell traits

The blood group Vel was discovered 60 years ago, but the underlying gene is unknown. Individuals negative for the Vel antigen are rare and are required for the safe transfusion of patients with antibodies to Vel. To identify the responsible gene, we sequenced the exomes of five individuals negative for the Vel antigen and found that four were homozygous and one was heterozygous for a low-frequency 17-nucleotide frameshift deletion in the gene encoding the 78-amino-acid transmembrane protein SMIM1. A follow-up study showing that 59 of 64 Vel-negative individuals were homozygous for the same deletion and expression of the Vel antigen on SMIM1-transfected cells confirm SMIM1 as the gene underlying the Vel blood group. An expression quantitative trait locus (eQTL), the common SNP rs1175550 contributes to variable expression of the Vel antigen (P = 0.003) and influences the mean hemoglobin concentration of red blood cells (RBCs; P = 8.6 × 10(-15)). In vivo, zebrafish with smim1 knockdown showed a mild reduction in the number of RBCs, identifying SMIM1 as a new regulator of RBC formation. Our findings are of immediate relevance, as the homozygous presence of the deletion allows the unequivocal identification of Vel-negative blood donors.