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.

Genitopatellar syndrome (GPS) is a rare disorder in which patellar aplasia or hypoplasia is associated with external genital anomalies and severe intellectual disability. Using an exome-sequencing approach, we identified de novo mutations of KAT6B in five individuals with GPS; a single nonsense variant and three frameshift indels, including a 4 bp deletion observed in two cases. All identified mutations are located within the terminal exon of the gene and are predicted to generate a truncated protein product lacking evolutionarily conserved domains. KAT6B encodes a member of the MYST family of histone acetyltranferases. We demonstrate a reduced level of both histone H3 and H4 acetylation in patient-derived cells suggesting that dysregulation of histone acetylation is a direct functional consequence of GPS alleles. These findings define the genetic basis of GPS and illustrate the complex role of the regulation of histone acetylation during development.

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.

Purpose Stargardt disease (STGD1) is caused by biallelic mutations in ABCA4, but many patients are genetically unsolved due to insensitive mutation‐scanning methods. We aimed to develop a cost‐effective sequencing method for ABCA4 exons and regions carrying known causal deep‐intronic variants. Methods Fifty exons and 12 regions containing 14 deep‐intronic variants of ABCA4 were sequenced using double‐tiled single molecule Molecular Inversion Probe (smMIP)‐based next‐generation sequencing. DNAs of 16 STGD1 cases carrying 29 ABCA4 alleles and of four healthy persons were sequenced using 483 smMIPs. Thereafter, DNAs of 411 STGD1 cases with one or no ABCA4 variant were sequenced. The effect of novel noncoding variants on splicing was analyzed using in vitro splice assays. Results Thirty‐four ABCA4 variants previously identified in 16 STGD1 cases were reliably identified. In 155/411 probands (38%), two causal variants were identified. We identified 11 deep‐intronic variants present in 62 alleles. Two known and two new noncanonical splice site variants showed splice defects, and one novel deep‐intronic variant (c.4539+2065C>G) resulted in a 170‐nt mRNA pseudoexon insertion (p.[Arg1514Lysfs*35,=]). Conclusions smMIPs‐based sequence analysis of coding and selected noncoding regions of ABCA4 enabled cost‐effective mutation detection in STGD1 cases in previously unsolved cases.

CHARGE (coloboma, heart defects, atresia of choanae, retardation of growth and development, genital hypoplasia, and ear abnormalities) and 22q11.2 deletion syndromes are variable, congenital malformation syndromes that show considerable phenotypic overlap. We further explored this clinical overlap and proposed recommendations for the genetic diagnosis of both syndromes. We described 2 patients clinically diagnosed with CHARGE syndrome, who were found to carry a 22q11.2 deletion, and searched the literature for more cases. In addition, we screened our cohort of CHD7 mutation carriers (n = 802) for typical 22q11.2 deletion features and studied CHD7 in 20 patients with phenotypically 22q11.2 deletion syndrome but without haploinsufficiency of TBX1. In total, we identified 5 patients with a clinical diagnosis of CHARGE syndrome and a proven 22q11.2 deletion. Typical 22q11.2 deletion features were found in 30 patients (30/802, 3.7%) of our CHD7 mutation-positive cohort. We found truncating CHD7 mutations in 5/20 patients with phenotypically 22q11.2 deletion syndrome. Differentiating between CHARGE and 22q11.2 deletion syndromes can be challenging. CHD7 and TBX1 probably share a molecular pathway or have common target genes in affected organs. We strongly recommend performing CHD7 analysis in patients with a 22q11.2 deletion phenotype without TBX1 haploinsufficiency and conversely, performing a genome-wide array in CHARGE syndrome patients without a CHD7 mutation.

Approximately 25–30% of colorectal cancer (CRC) cases are expected to result from a genetic predisposition, but in only 5–10% of these cases highly penetrant germline mutations are found. The remaining CRC heritability is still unexplained, and may be caused by a hitherto-undefined set of rare variants with a moderately penetrant risk. Here we aimed to identify novel risk factors for early-onset CRC using whole-exome sequencing, which was performed on a cohort of CRC individuals (n = 55) with a disease onset before 45 years of age. We searched for genes that were recurrently affected by rare variants (minor allele frequency ≤0.001) with potentially damaging effects and, subsequently, re-sequenced the candidate genes in a replication cohort of 174 early-onset or familial CRC individuals. Two functionally relevant genes with low frequency variants with potentially damaging effects, PTPN12 and LRP6, were found in at least three individuals. The protein tyrosine phosphatase PTP-PEST, encoded by PTPN12, is a regulator of cell motility and LRP6 is a component of the WNT-FZD-LRP5-LRP6 complex that triggers WNT signaling. All variants in LRP6 were identified in individuals with an extremely early-onset of the disease (≤30 years of age), and two of the three variants showed increased WNT signaling activity in vitro. In conclusion, we present PTPN12 and LRP6 as novel candidates contributing to the heterogeneous susceptibility to CRC.

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.

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.