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.

Clinical exome sequencing-Mistakes and caveats.

Massive parallel sequencing technology has become the predominant technique for genetic diagnostics and research. Many genetic laboratories have wrestled with the challenges of setting up genetic testing workflows based on a completely new technology. The learning curve we went through as a laboratory was accompanied by growing pains while we gained new knowledge and expertise. Here we discuss some important mistakes that have been made in our laboratory through 10 years of clinical exome sequencing but that have given us important new insights on how to adapt our working methods. We provide these examples and the lessons that we learned to help other laboratories avoid to make the same mistakes.

Novel BRCA1 and BRCA2 Tumor Test as Basis for Treatment Decisions and Referral for Genetic Counselling of Patients with Ovarian Carcinomas.

With the recent introduction of Poly(ADP-ribose) polymerase inhibitors, a promising novel therapy has become available for ovarian carcinoma (OC) patients with inactivating BRCA1 or BRCA2 mutations in their tumor. To select patients who may benefit from these treatments, assessment of the mutation status of BRCA1 and BRCA2 in the tumor is required. For reliable evaluation of germline and somatic mutations in these genes in DNA derived from formalin-fixed, paraffin-embedded (FFPE) tissue, we have developed a single-molecule molecular inversion probe (smMIP)-based targeted next-generation sequencing (NGS) approach. Our smMIP-based NGS approach provides analysis of both strands of the open reading frame of BRCA1 and BRCA2, enabling the discrimination between real variants and formalin-induced artefacts. The single molecule tag enables compilation of unique reads leading to a high analytical sensitivity and enabling assessment of the reliability of mutation-negative results. Multiplex ligation-dependent probe amplification (MLPA) and Methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) were used to detect exon deletions of BRCA1 and methylation of the BRCA1 promoter, respectively. Here, we show that this combined approach allows the rapid and reliable detection of both germline and somatic aberrations affecting BRCA1 and BRCA2 in DNA derived from FFPE OCs, enabling improved hereditary cancer risk assessment and clinical treatment of ovarian cancer patients.

Detection of clinically relevant copy-number variants by exome sequencing in a large cohort of genetic disorders.

PURPOSE: Copy-number variation is a common source of genomic variation and an important genetic cause of disease. Microarray-based analysis of copy-number variants (CNVs) has become a first-tier diagnostic test for patients with neurodevelopmental disorders, with a diagnostic yield of 10-20%. However, for most other genetic disorders, the role of CNVs is less clear and most diagnostic genetic studies are generally limited to the study of single-nucleotide variants (SNVs) and other small variants. With the introduction of exome and genome sequencing, it is now possible to detect both SNVs and CNVs using an exome- or genome-wide approach with a single test. METHODS: We performed exome-based read-depth CNV screening on data from 2,603 patients affected by a range of genetic disorders for which exome sequencing was performed in a diagnostic setting. RESULTS: In total, 123 clinically relevant CNVs ranging in size from 727 bp to 15.3 Mb were detected, which resulted in 51 conclusive diagnoses and an overall increase in diagnostic yield of ~2% (ranging from 0 to -5.8% per disorder). CONCLUSIONS: This study shows that CNVs play an important role in a broad range of genetic disorders and that detection via exome-based CNV profiling results in an increase in the diagnostic yield without additional testing, bringing us closer to single-test genomics.Genet Med advance online publication 27 October 2016.

BRCA Testing by Single-Molecule Molecular Inversion Probes.

BACKGROUND: Despite advances in next generation DNA sequencing (NGS), NGS-based single gene tests for diagnostic purposes require improvements in terms of completeness, quality, speed, and cost. Single-molecule molecular inversion probes (smMIPs) are a technology with unrealized potential in the area of clinical genetic testing. In this proof-of-concept study, we selected 2 frequently requested gene tests, those for the breast cancer genes BRCA1 and BRCA2, and developed an automated work flow based on smMIPs. METHODS: The BRCA1 and BRCA2 smMIPs were validated using 166 human genomic DNA samples with known variant status. A generic automated work flow was built to perform smMIP-based enrichment and sequencing for BRCA1, BRCA2, and the checkpoint kinase 2 (CHEK2) c.1100del variant. RESULTS: Pathogenic and benign variants were analyzed in a subset of 152 previously BRCA-genotyped samples, yielding an analytical sensitivity and specificity of 100%. Following automation, blind analysis of 65 in-house samples and 267 Norwegian samples correctly identified all true-positive variants (>3000), with no false positives. Consequent to process optimization, turnaround times were reduced by 60% to currently 10-15 days. Copy number variants were detected with an analytical sensitivity of 100% and an analytical specificity of 88%. CONCLUSIONS: smMIP-based genetic testing enables automated and reliable analysis of the coding sequences of BRCA1 and BRCA2. The use of single-molecule tags, double-tiled targeted enrichment, and capturing and sequencing in duplo, in combination with automated library preparation and data analysis, results in a robust process and reduces routine turnaround times. Furthermore, smMIP-based copy number variation analysis could make independent copy number variation tools like multiplex ligation-dependent probes amplification dispensable.

Is the $1000 Genome as Near as We Think? A Cost Analysis of Next-Generation Sequencing.

BACKGROUND: The substantial technological advancements in next-generation sequencing (NGS), combined with dropping costs, have allowed for a swift diffusion of NGS applications in clinical settings. Although several commercial parties report to have broken the $1000 barrier for sequencing an entire human genome, a valid cost overview for NGS is currently lacking. This study provides a complete, transparent and up-to-date overview of the total costs of different NGS applications. METHODS: Cost calculations for targeted gene panels (TGP), whole exome sequencing (WES) and whole genome sequencing (WGS) were based on the Illumina NextSeq500, HiSeq4000, and HiSeqX5 platforms, respectively. To anticipate future developments, sensitivity analyses are performed. RESULTS: Per-sample costs were €1669 for WGS, € 792 for WES and €333 for TGP. To reach the coveted $1000 genome, not only is the long-term and efficient use of the sequencing equipment needed, but also large reductions in capital costs and especially consumable costs are also required. CONCLUSIONS: WES and TGP are considerably lower-cost alternatives to WGS. However, this does not imply that these NGS approaches should be preferred in clinical practice, since this should be based on the tradeoff between costs and the expected clinical utility of the approach chosen. The results of the present study contribute to the evaluation of such tradeoffs.

Translating sanger-based routine DNA diagnostics into generic massive parallel ion semiconductor sequencing.

BACKGROUND: Dideoxy-based chain termination sequencing developed by Sanger is the gold standard sequencing approach and allows clinical diagnostics of disorders with relatively low genetic heterogeneity. Recently, new next generation sequencing (NGS) technologies have found their way into diagnostic laboratories, enabling the sequencing of large targeted gene panels or exomes. The development of benchtop NGS instruments now allows the analysis of single genes or small gene panels, making these platforms increasingly competitive with Sanger sequencing. METHODS: We developed a generic automated ion semiconductor sequencing work flow that can be used in a clinical setting and can serve as a substitute for Sanger sequencing. Standard amplicon-based enrichment remained identical to PCR for Sanger sequencing. A novel postenrichment pooling strategy was developed, limiting the number of library preparations and reducing sequencing costs up to 70% compared to Sanger sequencing. RESULTS: A total of 1224 known pathogenic variants were analyzed, yielding an analytical sensitivity of 99.92% and specificity of 99.99%. In a second experiment, a total of 100 patient-derived DNA samples were analyzed using a blind analysis. The results showed an analytical sensitivity of 99.60% and specificity of 99.98%, comparable to Sanger sequencing. CONCLUSIONS: Ion semiconductor sequencing can be a first choice mutation scanning technique, independent of the genes analyzed.

Best practice guidelines for the use of next-generation sequencing applications in genome diagnostics: a national collaborative study of Dutch genome diagnostic laboratories.

Next-generation sequencing (NGS) methods are being adopted by genome diagnostics laboratories worldwide. However, implementing NGS-based tests according to diagnostic standards is a challenge for individual laboratories. To facilitate the implementation of NGS in Dutch laboratories, the Dutch Society for Clinical Genetic Laboratory Diagnostics (VKGL) set up a working group in 2012. The results of their discussions are presented here. We provide best practice guidelines and criteria for implementing and validating NGS applications in a clinical setting. We introduce the concept of "diagnostic yield" as the main performance characteristic for evaluating diagnostic tests. We recommend that the laboratory procedures, including the tested genes, should be recorded in a publicly available document describing the complete "diagnostic routing." We also propose that laboratories should use a list of "core disease genes" for specific genetic diseases. This core list contains the essential genes for each disease, and they should all be included in a diagnostic test to establish a reliable and accurate molecular diagnosis. The guidelines will ensure a clear and standardized quality of care provided by genetic diagnostic laboratories. The best practice guidelines and criteria that are presented here were adopted by the VKGL in January 2013.

A post-hoc comparison of the utility of sanger sequencing and exome sequencing for the diagnosis of heterogeneous diseases.

The advent of massive parallel sequencing is rapidly changing the strategies employed for the genetic diagnosis and research of rare diseases that involve a large number of genes. So far it is not clear whether these approaches perform significantly better than conventional single gene testing as requested by clinicians. The current yield of this traditional diagnostic approach depends on a complex of factors that include gene-specific phenotype traits, and the relative frequency of the involvement of specific genes. To gauge the impact of the paradigm shift that is occurring in molecular diagnostics, we assessed traditional Sanger-based sequencing (in 2011) and exome sequencing followed by targeted bioinformatics analysis (in 2012) for five different conditions that are highly heterogeneous, and for which our center provides molecular diagnosis. We find that exome sequencing has a much higher diagnostic yield than Sanger sequencing for deafness, blindness, mitochondrial disease, and movement disorders. For microsatellite-stable colorectal cancer, this was low under both strategies. Even if all genes that could have been ordered by physicians had been tested, the larger number of genes captured by the exome would still have led to a clearly superior diagnostic yield at a fraction of the cost.

Newborn screening for Cerebrotendinous Xanthomatosis: A retrospective biomarker study using both flow-injection and UPLC-MS/MS analysis in 20,000 newborns.

BACKGROUND AND AIMS: Cerebrotendinous Xanthomatosis (CTX) is a treatable disorder of bile acid synthesis caused by deficiency of 27-sterol hydroxylase -encoded by CYP27A1- leading to gastrointestinal and progressive neuropsychiatric symptoms. Biochemically, CTX is characterized by accumulation of the bile alcohol cholestanetetrol glucuronide (GlcA-tetrol) and the deficiency of tauro-chenodeoxycholic acid (t-CDCA) and tauro-trihydroxycholestanoic acid (t-THCA). MATERIALS AND METHODS: To ascertain the feasibility of CTX newborn screening (NBS) we performed a study with deidentified Dutch dried blood spots using reagents and equipment that is frequently used in NBS laboratories. 20,076 deidentified newborn blood spots were subjected to flow-injection (FIA)-MS/MS and UPLC-MS/MS analysis to determine the concentration of GlcA-tetrol and calculate the GlcA-tetrol/t-CDCA and t-THCA/GlcA-tetrol ratios. RESULTS: Using UPLC-MS/MS analysis both GlcA-tetrol concentration and/or metabolite ratios GlcA-tetrol/t-CDCA proved to be informative biomarkers; newborn DBS results did not overlap with those of the CTX patients. For FIA-MS/MS, GlcA-tetrol also was an excellent marker but when using the combination of the GlcA-tetrol/t-CDCA and t-THCA/GlcA-tetrol ratios also did not yield any screen positives. CONCLUSION: Newborn screening for CTX using only metabolite ratios following the measurement of three CTX biomarkers is possible using either FIA-MS/MS or UPLC-MS/MS, which paves the way for introduction of CTX NBS.

The performance of genome sequencing as a first-tier test for neurodevelopmental disorders.

Genome sequencing (GS) can identify novel diagnoses for patients who remain undiagnosed after routine diagnostic procedures. We tested whether GS is a better first-tier genetic diagnostic test than current standard of care (SOC) by assessing the technical and clinical validity of GS for patients with neurodevelopmental disorders (NDD). We performed both GS and exome sequencing in 150 consecutive NDD patient-parent trios. The primary outcome was diagnostic yield, calculated from disease-causing variants affecting exonic sequence of known NDD genes. GS (30%, n = 45) and SOC (28.7%, n = 43) had similar diagnostic yield. All 43 conclusive diagnoses obtained with SOC testing were also identified by GS. SOC, however, required integration of multiple test results to obtain these diagnoses. GS yielded two more conclusive diagnoses, and four more possible diagnoses than ES-based SOC (35 vs. 31). Interestingly, these six variants detected only by GS were copy number variants (CNVs). Our data demonstrate the technical and clinical validity of GS to serve as routine first-tier genetic test for patients with NDD. Although the additional diagnostic yield from GS is limited, GS comprehensively identified all variants in a single experiment, suggesting that GS constitutes a more efficient genetic diagnostic workflow.

Towards Next-Generation Sequencing (NGS)-Based Newborn Screening: A Technical Study to Prepare for the Challenges Ahead.

Newborn screening (NBS) aims to identify neonates with severe conditions for whom immediate treatment is required. Currently, a biochemistry-first approach is used to identify these disorders, which are predominantly inherited meta1bolic disorders (IMD). Next-generation sequencing (NGS) is expected to have some advantages over the current approach, for example the ability to detect IMDs that meet all screening criteria but lack an identifiable biochemical footprint. We have now designed a technical study to explore the use of NGS techniques as a first-tier approach in NBS. Here, we describe the aim and set-up of the NGS-first for the NBS (NGSf4NBS) project, which will proceed in three steps. In Step 1, we will identify IMDs eligible for NGS-first testing, based on treatability. In Step 2, we will investigate the feasibility, limitations and comparability of different technical NGS approaches and analysis workflows for NBS, eventually aiming to develop a rapid NGS-based workflow. Finally, in Step 3, we will prepare for the incorporation of this workflow into the existing Dutch NBS program and propose a protocol for referral of a child after a positive NGS test result. The results of this study will be the basis for an additional analytical route within NBS that will be further studied for its applicability within the NBS program, e.g., regarding the ethical, legal, financial and social implications.

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.

Common breast cancer susceptibility alleles are associated with tumour subtypes in BRCA1 and BRCA2 mutation carriers: results from the Consortium of Investigators of Modifiers of BRCA1/2.

INTRODUCTION: Previous studies have demonstrated that common breast cancer susceptibility alleles are differentially associated with breast cancer risk for BRCA1 and/or BRCA2 mutation carriers. It is currently unknown how these alleles are associated with different breast cancer subtypes in BRCA1 and BRCA2 mutation carriers defined by estrogen (ER) or progesterone receptor (PR) status of the tumour. METHODS: We used genotype data on up to 11,421 BRCA1 and 7,080 BRCA2 carriers, of whom 4,310 had been affected with breast cancer and had information on either ER or PR status of the tumour, to assess the associations of 12 loci with breast cancer tumour characteristics. Associations were evaluated using a retrospective cohort approach. RESULTS: The results suggested stronger associations with ER-positive breast cancer than ER-negative for 11 loci in both BRCA1 and BRCA2 carriers. Among BRCA1 carriers, single nucleotide polymorphism (SNP) rs2981582 (FGFR2) exhibited the biggest difference based on ER status (per-allele hazard ratio (HR) for ER-positive = 1.35, 95% CI: 1.17 to 1.56 vs HR = 0.91, 95% CI: 0.85 to 0.98 for ER-negative, P-heterogeneity = 6.5 × 10-6). In contrast, SNP rs2046210 at 6q25.1 near ESR1 was primarily associated with ER-negative breast cancer risk for both BRCA1 and BRCA2 carriers. In BRCA2 carriers, SNPs in FGFR2, TOX3, LSP1, SLC4A7/NEK10, 5p12, 2q35, and 1p11.2 were significantly associated with ER-positive but not ER-negative disease. Similar results were observed when differentiating breast cancer cases by PR status. CONCLUSIONS: The associations of the 12 SNPs with risk for BRCA1 and BRCA2 carriers differ by ER-positive or ER-negative breast cancer status. The apparent differences in SNP associations between BRCA1 and BRCA2 carriers, and non-carriers, may be explicable by differences in the prevalence of tumour subtypes. As more risk modifying variants are identified, incorporating these associations into breast cancer subtype-specific risk models may improve clinical management for mutation carriers.

Long-read trio sequencing of individuals with unsolved intellectual disability.

Long-read sequencing (LRS) has the potential to comprehensively identify all medically relevant genome variation, including variation commonly missed by short-read sequencing (SRS) approaches. To determine this potential, we performed LRS around 15×-40× genome coverage using the Pacific Biosciences Sequel I System for five trios. The respective probands were diagnosed with intellectual disability (ID) whose etiology remained unresolved after SRS exomes and genomes. Systematic assessment of LRS coverage showed that ~35 Mb of the human reference genome was only accessible by LRS and not SRS. Genome-wide structural variant (SV) calling yielded on average 28,292 SV calls per individual, totaling 12.9 Mb of sequence. Trio-based analyses which allowed to study segregation, showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming that these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, ~100,000 were located in the ~35 Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known Mendelian disease genes. Our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

Evidence for SMAD3 as a modifier of breast cancer risk in BRCA2 mutation carriers.

INTRODUCTION: Current attempts to identify genetic modifiers of BRCA1 and BRCA2 associated risk have focused on a candidate gene approach, based on knowledge of gene functions, or the development of large genome-wide association studies. In this study, we evaluated 24 SNPs tagged to 14 candidate genes derived through a novel approach that analysed gene expression differences to prioritise candidate modifier genes for association studies. METHODS: We successfully genotyped 24 SNPs in a cohort of up to 4,724 BRCA1 and 2,693 BRCA2 female mutation carriers from 15 study groups and assessed whether these variants were associated with risk of breast cancer in BRCA1 and BRCA2 mutation carriers. RESULTS: SNPs in five of the 14 candidate genes showed evidence of association with breast cancer risk for BRCA1 or BRCA2 carriers (P < 0.05). Notably, the minor alleles of two SNPs (rs7166081 and rs3825977) in high linkage disequilibrium (r² = 0.77), located at the SMAD3 locus (15q22), were each associated with increased breast cancer risk for BRCA2 mutation carriers (relative risk = 1.25, 95% confidence interval = 1.07 to 1.45, P(trend) = 0.004; and relative risk = 1.20, 95% confidence interval = 1.03 to 1.40, P(trend) = 0.018). CONCLUSIONS: This study provides evidence that the SMAD3 gene, which encodes a key regulatory protein in the transforming growth factor beta signalling pathway and is known to interact directly with BRCA2, may contribute to increased risk of breast cancer in BRCA2 mutation carriers. This finding suggests that genes with expression associated with BRCA1 and BRCA2 mutation status are enriched for the presence of common genetic modifiers of breast cancer risk in these populations.

A co-segregating microduplication of chromosome 15q11.2 pinpoints two risk genes for autism spectrum disorder.

High resolution genomic copy-number analysis has shown that inherited and de novo copy-number variations contribute significantly to autism pathology, and that identification of small chromosomal aberrations related to autism will expedite the discovery of risk genes involved. Here, we report a microduplication of chromosome 15q11.2, spanning only four genes, co-segregating with autism in a Dutch pedigree, identified by SNP microarray analysis, and independently confirmed by FISH and MLPA analysis. Quantitative RT-PCR analysis revealed over 70% increase in peripheral blood mRNA levels for the four genes present in the duplicated region in patients, and RNA in situ hybridization on mouse embryonic and adult brain sections revealed that two of the four genes, CYFIP1 and NIPA1, were highly expressed in the developing mouse brain. These findings point towards a contribution of microduplications at chromosome 15q11.2 to autism, and highlight CYFIP1 and NIPA1 as autism risk genes functioning in axonogenesis and synaptogenesis. Thereby, these findings further implicate defects in dosage-sensitive molecular control of neuronal connectivity in autism. However, the prevalence of this microduplication in patient samples was statistically not significantly different from control samples (0.94% in patients vs. 0.42% controls, P = 0.247), which suggests that our findings should be interpreted with caution and indicates the need for studies that include large numbers of control subjects to ascertain the impact of these changes on a population scale.

Activation delay and VT parameters in arrhythmogenic right ventricular dysplasia/cardiomyopathy: toward improvement of diagnostic ECG criteria.

INTRODUCTION: Desmosomal changes, electrical uncoupling, and surviving myocardial bundles embedded in fibrofatty tissue are hallmarks of activation delay in arrhythmogenic right ventricular dysplasia/cardiomyopathy (ARVD/C). Currently, generally accepted task force criteria (TFC) are used for clinical diagnosis. We propose additional criteria based on activation delay and ventricular tachycardia (VT) to improve identification of affected individuals. METHODS AND RESULTS: Activation delay and VT-related 12-lead electrocardiographic (ECG) criteria were studied, while off drugs, in 42 patients with proven ARVD/C according to TFC, and 27 controls with idiopathic VT from the RV outflow tract. Two of three measured TFC could only be identified in a small minority of ARVD/C patients. Additional ECG criteria proposed in this study included (a) prolonged terminal activation duration, an indicator of activation delay; (b) VT with LBBB morphology and superior axis; and (c) multiple different VT morphologies. These criteria were met in 30 (71%), 28 (67%), and 37 (88%) ARVD/C patients, respectively, and in one control patient (P < 0.001). Electrophysiologic studies contributed importantly to yield different VT morphologies. Pathogenic plakophilin-2 mutations were identified in 25 (60%) of ARVD/C patients and in none of the controls. In ARVD/C patients, parameters measured were not significantly different between mutation carriers and noncarriers, except for negative T waves in V1-3, occurring more frequently in patients with mutation. CONCLUSIONS: The proposed additional criteria are specific for ARVD/C and more sensitive than the current TFC. Therefore, adding the newly proposed criteria to current TFC could improve ARVD/C diagnosis, independent of DNA analysis.

DNA analysis of AHI1, NPHP1 and CYCLIN D1 in Joubert syndrome patients from the Netherlands.

Joubert syndrome (JBS) is a clinically variable and genetically heterogeneous developmental brain disorder with autosomal recessive inheritance. Five genes, AHI1, NPHP1, CEP290, MKS3, and RPGRIP1L, and two additional loci on chromosome 9 and 11 have been identified so far. The relative contributions of AHI1 mutations and NPHP1 deletions have not yet been determined in a population-based JBS patient cohort. We therefore undertook a nationwide survey of JBS in the Netherlands and performed DNA analysis of the AHI1 and NPHP1 genes, as well as a new candidate gene CYCLIN D1. We obtained clinical data and DNA samples of 25 Dutch JBS patients. DNA analysis of AHI1 revealed pathogenic homozygous or compound heterozygous AHI1 mutations in four patients (16%). Based on the birth prevalence of about 1 in 100,000 for JBS in the Netherlands, we estimated a carrier frequency of AHI1 mutations of approximately 1 in 400. In another two patients, the AHI1 mutation Arg830Trp was identified (homozygously and heterozygously), a possible low penetrance allele. No deletions of NPHP1 or CYCLIN D1 mutations were detected in these 25 patients. In the four patients with AHI1 mutations, retinal disease (Leber congenital amaurosis or retinal dystrophy) was present in two, whereas none had renal disease. Pooling our data and data from the literature, retinal disease seems to occur in 75% of AHI1-associated JBS patients. Renal disease is present in 10% at most. We conclude that AHI1 mutations are an important cause of JBS in Dutch patients, and should always be looked for in patients suspected of JBS, especially when retinal dystrophy is present. Patients with AHI1 mutations should be regularly checked for retinal and renal disease up until adolescence.