High throughput sequencing has greatly advanced disease gene identification, especially in heterogeneous entities. Despite falling costs this is still an expensive and laborious technique, particularly when studying large cohorts. To address this problem we applied Exome Pool-Seq as an economic and fast screening technology in neurodevelopmental disorders (NDDs). Sequencing of 96 individuals can be performed in eight pools of 12 samples on less than one Illumina sequencer lane. In a pilot study with 96 cases we identified 27 variants, likely or possibly affecting function. Twenty five of these were identified in 923 established NDD genes (based on SysID database, status November 2016) (ACTB, AHDC1, ANKRD11, ATP6V1B2, ATRX, CASK, CHD8, GNAS, IFIH1, KCNQ2, KMT2A, KRAS, MAOA, MED12, MED13L, RIT1, SETD5, SIN3A, TCF4, TRAPPC11, TUBA1A, WAC, ZBTB18, ZMYND11), two in 543 (SysID) candidate genes (ZNF292, BPTF), and additionally a de novo loss-of-function variant in LRRC7, not previously implicated in NDDs. Most of them were confirmed to be de novo, but we also identified X-linked or autosomal-dominantly or autosomal-recessively inherited variants. With a detection rate of 28%, Exome Pool-Seq achieves comparable results to individual exome analyses but reduces costs by >85%. Compared with other large scale approaches using Molecular Inversion Probes (MIP) or gene panels, it allows flexible re-analysis of data. Exome Pool-Seq is thus well suited for large-scale, cost-efficient and flexible screening in characterized but heterogeneous entities like NDDs.

The TUBA1A-associated tubulinopathy is clinically heterogeneous with brain malformations, microcephaly, developmental delay and epilepsy being the main clinical features. It is an autosomal dominant disorder mostly caused by de novo variants in TUBA1A. In three individuals with developmental delay we identified heterozygous de novo missense variants in TUBA1A using exome sequencing. While the c.1307G > A, p.(Gly436Asp) variant was novel, the two variants c.518C > T, p.(Pro173Leu) and c.641G > A, p.(Arg214His) were previously described. We compared the variable phenotype observed in these individuals with a carefully conducted review of the current literature and identified 166 individuals, 146 born and 20 fetuses with a TUBA1A variant. In 107 cases with available clinical information we standardized the reported phenotypes according to the Human Phenotype Ontology. The most commonly reported features were developmental delay (98%), anomalies of the corpus callosum (96%), microcephaly (76%) and lissencephaly (agyria-pachygyria) (70%), although reporting was incomplete in the different studies. We identified a total of 121 specific variants, including 15 recurrent ones. Missense variants cluster in the C-terminal region around the most commonly affected amino acid position Arg402 (13.3%). In a three-dimensional protein model, 38.6% of all disease-causing variants including those in the C-terminal region are predicted to affect the binding of microtubule-associated proteins or motor proteins. Genotype-phenotype analysis for recurrent variants showed an overrepresentation of certain clinical features. However, individuals with these variants are often reported in the same publication. With 166 individuals, we present the most comprehensive phenotypic and genotypic standardized synopsis for clinical interpretation of TUBA1A variants. Despite this considerable number, a detailed genotype-phenotype characterization is limited by large inter-study variability in reporting.

Genetic integrity of induced pluripotent stem cells (iPSCs) is essential for their validity as disease models and for potential therapeutic use. We describe the comprehensive analysis in the ForIPS consortium: an iPSC collection from donors with neurological diseases and healthy controls. Characterization included pluripotency confirmation, fingerprinting, conventional and molecular karyotyping in all lines. In the majority, somatic copy number variants (CNVs) were identified. A subset with available matched donor DNA was selected for comparative exome sequencing. We identified single nucleotide variants (SNVs) at different allelic frequencies in each clone with high variability in mutational load. Low frequencies of variants in parental fibroblasts highlight the importance of germline samples. Somatic variant number was independent from reprogramming, cell type and passage. Comparison with disease genes and prediction scores suggest biological relevance for some variants. We show that high-throughput sequencing has value beyond SNV detection and the requirement to individually evaluate each clone.

Recently, the Mucin-1 (MUC1) gene has been identified as a causal gene of autosomal dominant tubulointerstitial kidney disease (ADTKD). Most causative mutations are buried within a GC-rich 60 basepair variable number of tandem repeat (VNTR), which escapes identification by massive parallel sequencing methods due to the complexity of the VNTR. We established long read single molecule real time sequencing (SMRT) targeted to the MUC1-VNTR as an alternative strategy to the snapshot assay. Our approach allows complete VNTR assembly, thereby enabling the detection of all variants residing within the VNTR and simultaneous determination of VNTR length. We present high resolution data on the VNTR architecture for a cohort of snapshot positive (n = 9) and negative (n = 7) ADTKD families. By SMRT sequencing we could confirm the diagnosis in all previously tested cases, reconstruct both VNTR alleles and determine the exact position of the causative variant in eight of nine families. This study demonstrates that precise positioning of the causative mutation(s) and identification of other coding and noncoding sequence variants in ADTKD-MUC1 is feasible. SMRT sequencing could provide a powerful tool to uncover potential factors encoded within the VNTR that associate with intra- and interfamilial phenotype variability of MUC1 related kidney disease.

Two percent of patients with Wilms tumors have a positive family history. In many of these cases the genetic cause remains unresolved. By applying germline exome sequencing in two families with two affected individuals with Wilms tumors, we identified truncating mutations in TRIM28. Subsequent mutational screening of germline and tumor DNA of 269 children affected by Wilms tumor was performed, and revealed seven additional individuals with germline truncating mutations, and one individual with a somatic truncating mutation in TRIM28. TRIM28 encodes a complex scaffold protein involved in many different processes, including gene silencing, DNA repair and maintenance of genomic integrity. Expression studies on mRNA and protein level showed reduction of TRIM28, confirming a loss‐of‐function effect of the mutations identified. The tumors showed an epithelial‐type histology that stained negative for TRIM28 by immunohistochemistry. The tumors were bilateral in six patients, and 10/11 tumors are accompanied by perilobar nephrogenic rests. Exome sequencing on eight tumor DNA samples from six individuals showed loss‐of‐heterozygosity (LOH) of the TRIM28‐locus by mitotic recombination in seven tumors, suggesting that TRIM28 functions as a tumor suppressor gene in Wilms tumor development. Additionally, the tumors showed very few mutations in known Wilms tumor driver genes, suggesting that loss of TRIM28 is the main driver of tumorigenesis. In conclusion, we identified heterozygous germline truncating mutations in TRIM28 in 11 children with mainly epithelial‐type Wilms tumors, which become homozygous in tumor tissue. These data establish TRIM28 as a novel Wilms tumor predisposition gene, acting as a tumor suppressor gene by LOH.

Intellectual disability (ID) is often associated with behavioral problems or disorders. Mutations in the GRIN2B gene (MRD6, MIM613970) have been identified as a common cause of ID (prevalence of 0.5 – 1% in individuals with ID) associated with EEG and behavioral problems. We assessed five GRIN2B mutation carriers aged between 3 and 14 years clinically and via standardized questionnaires to delineate a detailed behavioral phenotype. Parents and teachers rated problem behavior of their affected children by completing the Developmental Behavior Checklist (DBC) and the Conners’ Rating Scales Revised (CRS-R:L). All individuals had mild to severe ID and needed guidance in daily routine. They showed characteristic behavior problems with prominent hyperactivity, impulsivity, distractibility and a short attention span. Stereotypies, sleeping problems and a friendly but boundless social behavior were commonly reported. Our observations provide an initial delineation of the behavioral phenotype of GRIN2B mutation carriers.