The majority of next-generation sequencing short-reads can be properly aligned by leading aligners at high speed. However, the alignment quality can still be further improved, since usually not all ...reads can be correctly aligned to large genomes, such as the human genome, even for simulated data. Moreover, even slight improvements in this area are important but challenging, and usually require significantly more computational endeavor. In this paper, we present CUSHAW3, an open-source parallelized, sensitive and accurate short-read aligner for both base-space and color-space sequences. In this aligner, we have investigated a hybrid seeding approach to improve alignment quality, which incorporates three different seed types, i.e. maximal exact match seeds, exact-match k-mer seeds and variable-length seeds, into the alignment pipeline. Furthermore, three techniques: weighted seed-pairing heuristic, paired-end alignment pair ranking and read mate rescuing have been conceived to facilitate accurate paired-end alignment. For base-space alignment, we have compared CUSHAW3 to Novoalign, CUSHAW2, BWA-MEM, Bowtie2 and GEM, by aligning both simulated and real reads to the human genome. The results show that CUSHAW3 consistently outperforms CUSHAW2, BWA-MEM, Bowtie2 and GEM in terms of single-end and paired-end alignment. Furthermore, our aligner has demonstrated better paired-end alignment performance than Novoalign for short-reads with high error rates. For color-space alignment, CUSHAW3 is consistently one of the best aligners compared to SHRiMP2 and BFAST. The source code of CUSHAW3 and all simulated data are available at http://cushaw3.sourceforge.net.
Uterine leiomyomas (ULs) constitute a considerable health burden in the general female population. The fumarate hydratase (FH) deficient subtype is found in up to 1.6% and can occur in hereditary ...leiomyomatosis and renal cell carcinoma (HLRCC) syndrome. We sequenced 13 FH deficient ULs from a previous immunohistochemical screen using a targeted panel and identified biallelic FH variants in all. In eight, we found an FH point mutation (two truncating, six missense) with evidence for loss of the second allele. Variant allele-frequencies in all cases with a point mutation pointed to somatic variants. Spatial clustering of the identified missense variants in the lyase domain indicated altered fumarase oligomerization with subsequent degradation as explanation for the observed FH deficiency. Biallelic FH deletions in five tumors confirm the importance of copy number loss as mutational mechanism. By curating all pathogenic FH variants and calculating their population frequency, we estimate a carrier frequency of up to 1/2,563. Comparing with the prevalence of FH deficient ULs, we conclude that most are sporadic and estimate 2.7–13.9% of females with an FH deficient UL to carry a germline FH variant. Further prospective tumor/normal sequencing studies are needed to develop a reliable screening strategy for HLRCC in women with ULs.
Intellectual disability (ID) is a clinically and genetically heterogeneous common condition that remains etiologically unresolved in the majority of cases. Although several hundred diseased genes ...have been identified in X-linked, autosomal-recessive, or syndromic types of ID, the establishment of an etiological basis remains a difficult task in unspecific, sporadic cases. Just recently, de novo mutations in SYNGAP1, STXBP1, MEF2C, and GRIN2B were reported as relatively common causes of ID in such individuals. On the basis of a patient with severe ID and a 2.5 Mb microdeletion including ARID1B in chromosomal region 6q25, we performed mutational analysis in 887 unselected patients with unexplained ID. In this cohort, we found eight (0.9%) additional de novo nonsense or frameshift mutations predicted to cause haploinsufficiency. Our findings indicate that haploinsufficiency of ARID1B, a member of the SWI/SNF-A chromatin-remodeling complex, is a common cause of ID, and they add to the growing evidence that chromatin-remodeling defects are an important contributor to neurodevelopmental disorders.
N-methyl-D-aspartate (NMDA) receptors mediate excitatory neurotransmission in the mammalian brain. Two glycine-binding NR1 subunits and two glutamate-binding NR2 subunits each form highly ...Ca2+-permeable cation channels which are blocked by extracellular Mg2+ in a voltage-dependent manner. Either GRIN2B or GRIN2A, encoding the NMDA receptor subunits NR2B and NR2A, was found to be disrupted by chromosome translocation breakpoints in individuals with mental retardation and/or epilepsy. Sequencing of GRIN2B in 468 individuals with mental retardation revealed four de novo mutations: a frameshift, a missense and two splice-site mutations. In another cohort of 127 individuals with idiopathic epilepsy and/or mental retardation, we discovered a GRIN2A nonsense mutation in a three-generation family. In a girl with early-onset epileptic encephalopathy, we identified the de novo GRIN2A mutation c.1845C>A predicting the amino acid substitution p.N615K. Analysis of NR1-NR2AN615K (NR2A subunit with the p.N615K alteration) receptor currents revealed a loss of the Mg2+ block and a decrease in Ca2+ permeability. Our findings suggest that disturbances in the neuronal electrophysiological balance during development result in variable neurological phenotypes depending on which NR2 subunit of NMDA receptors is affected.
Autosomal dominant polycystic kidney disease (ADPKD) represents the most common hereditary nephropathy. Despite growing evidence for genetic heterogeneity, ADPKD diagnosis is still primarily based ...upon clinical imaging criteria established before discovery of additional PKD genes. This study aimed at assessing the diagnostic value of genetic verification in clinical ADPKD.
In this prospective, diagnostic trial, 100 families with clinically diagnosed ADPKD were analyzed by PKD gene panel and multiplex ligation-dependent probe amplification (MLPA); exome sequencing (ES) was performed in panel/MLPA-negative families.
Diagnostic PKD1/2 variants were identified in 81 families (81%), 70 of which in PKD1 and 11 in PKD2. PKD1 variants of unknown significance were detected in another 9 families (9%). Renal survival was significantly worse upon PKD1 truncation versus nontruncation and PKD2 alteration. Ten percent of the cohort were PKD1/2-negative, revealing alternative genetic diagnoses such as autosomal recessive PKD, Birt-Hogg-Dubé syndrome, and ALG9-associated PKD. In addition, among unsolved cases, ES yielded potential novel PKD candidates.
By illustrating vast genetic heterogeneity, this study demonstrates the value of genetic testing in a real-world PKD cohort by diagnostic verification, falsification, and disease prediction. In the era of specific treatment for fast progressive ADPKD, genetic confirmation should form the basis of personalized patient care.
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 diagnostics of neurodevelopmental disorders with epilepsy (NDDE) are predominantly applied in children, thus limited information is available regarding adults or elderly.
We investigated 150 ...adult/elderly individuals with NDDE by conventional karyotyping, FMR1 testing, chromosomal microarray, panel sequencing, and for unresolved cases, also by exome sequencing (nsingle = 71, ntrios = 24).
We identified (likely) pathogenic variants in 71 cases (47.3%) comprising fragile X syndrome (n = 1), disease-causing copy number (n = 23), and single-nucleotide variants (n = 49). Seven individuals displayed multiple independent genetic diagnoses. The diagnostic yield correlated with the severity of intellectual disability. Individuals with anecdotal evidence of exogenic early-life events (e.g., nuchal cord, complications at delivery) with alleged/unproven association to the disorder had a particularly high yield of 58.3%. Screening for disease-specific comorbidities was indicated in 45.1% and direct treatment consequences arose in 11.8% of diagnosed individuals.
Panel/exome sequencing displayed the highest yield and should be considered as first-tier diagnostics in NDDE. This high yield and the numerous indications for additional screening or treatment modifications arising from genetic diagnoses indicate a current medical undersupply of genetically undiagnosed adult/elderly individuals with NDDE. Moreover, knowledge of the course of elderly individuals will ultimately help in counseling newly diagnosed individuals with NDDE.
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.
What's new?
About 2% of Wilms tumors run in families, and some of the mutations remain unknown. These authors identified a new Wilms tumor mutation, a truncation on the TRIM28 gene. They started by performing exome sequencing on tumors in pairs of affected children from 2 families. In these 4 patients, they found mutations in TRIM28, which encodes a scaffold protein involved in DNA repair and genome stability. They then screened a cohort of 269 cases and found 8 more patients bearing TRIM28 loss‐of‐function mutations. The gene appears to function as a tumor suppressor with loss of heterozygosity in the tumor cells.
The risk of epilepsy among individuals with intellectual disability (ID) is approximately ten times that of the general population. From a cohort of >5,000 families affected by neurodevelopmental ...disorders, we identified six consanguineous families harboring homozygous inactivating variants in MBOAT7, encoding lysophosphatidylinositol acyltransferase (LPIAT1). Subjects presented with ID frequently accompanied by epilepsy and autistic features. LPIAT1 is a membrane-bound phospholipid-remodeling enzyme that transfers arachidonic acid (AA) to lysophosphatidylinositol to produce AA-containing phosphatidylinositol. This study suggests a role for AA-containing phosphatidylinositols in the development of ID accompanied by epilepsy and autistic features.
A few de novo missense variants in the cytoplasmic FMRP-interacting protein 2 (CYFIP2) gene have recently been described as a novel cause of severe intellectual disability, seizures, and hypotonia in ...18 individuals, with p.Arg87 substitutions in the majority.
We assembled data from 19 newly identified and all 18 previously published individuals with CYFIP2 variants. By structural modeling and investigation of WAVE-regulatory complex (WRC)-mediated actin polymerization in six patient fibroblast lines we assessed the impact of CYFIP2 variants on the WRC.
Sixteen of 19 individuals harbor two previously described and 11 novel (likely) disease-associated missense variants. We report p.Asp724 as second mutational hotspot (4/19 cases). Genotype-phenotype correlation confirms a consistently severe phenotype in p.Arg87 patients but a more variable phenotype in p.Asp724 and other substitutions. Three individuals with milder phenotypes carry putative loss-of-function variants, which remain of unclear pathogenicity. Structural modeling predicted missense variants to disturb interactions within the WRC or impair CYFIP2 stability. Consistent with its role in WRC-mediated actin polymerization we substantiate aberrant regulation of the actin cytoskeleton in patient fibroblasts.
Our study expands the clinical and molecular spectrum of CYFIP2-related neurodevelopmental disorder and provides evidence for aberrant WRC-mediated actin dynamics as contributing cellular pathomechanism.