Mobster: Accurate detection of mobile element insertions in next generation sequencing data

Mobile elements are major drivers in changing genomic architecture and can cause disease. The detection of mobile elements is hindered due to the low mappability of their highly repetitive sequences. We have developed an algorithm, called Mobster, to detect non-reference mobile element insertions in next generation sequencing data from both whole genome and whole exome studies. Mobster uses discordant read pairs and clipped reads in combination with consensus sequences of known active mobile elements. Mobster has a low false discovery rate and high recall rate for both L1 and Alu elements. Mobster is available at http://sourceforge.net/projects/mobster. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13059-014-0488-x) contains supplementary material, which is available to authorized users

Oto-facial syndrome and esophageal atresia, intellectual disability and zygomatic anomalies: expanding the phenotypes associated with EFTUD2 mutations

Background: Mutations in EFTUD2 were proven to cause a very distinct mandibulofacial dysostosis type Guion-Almeida (MFDGA, OMIM #610536). Recently, gross deletions and mutations in EFTUD2 were determined to cause syndromic esophageal atresia (EA), as well. We set forth to find further conditions caused by mutations in the EFTUD2 gene (OMIM *603892). Methods and results: We performed exome sequencing in two familial cases with clinical features overlapping with MFDGA and EA, but which were previously assumed to represent distinct entities, a syndrome with esophageal atresia, hypoplasia of zygomatic complex, microcephaly, cup-shaped ears, congenital heart defect, and intellectual disability in a mother and her two children [AJMG 143A(11):1135-1142, 2007] and a supposedly autosomal recessive oto-facial syndrome with midline malformations in two sisters [AJMG 132(4):398-401, 2005]. While the analysis of our exome data was in progress, a recent publication made EFTUD2 mutations highly likely in these families. This hypothesis could be confirmed with exome as well as with Sanger sequencing. Also, in three further sporadic patients, clinically overlapping to these two families, de novo mutations within EFTUD2 were identified by Sanger sequencing. Our clinical and molecular workup of the patients discloses a broad phenotypic spectrum, and describes for the first time an instance of germline mosaicism for an EFTUD2 mutation. Conclusions: The clinical features of the eight patients described here further broaden the phenotypic spectrum caused by EFTUD2 mutations or deletions. We here show, that it not only includes mandibulofacial dysostosis type Guion-Almeida, which should be reclassified as an acrofacial dysostosis because of thumb anomalies (present in 12/35 or 34% of patients) and syndromic esophageal atresia [JMG 49(12). 737-746, 2012], but also the two new syndromes, namely oto-facial syndrome with midline malformations published by Megarbane et al. [AJMG 132(4): 398-401, 2005] and the syndrome published by Wieczorek et al. [AJMG 143A(11):1135-1142, 2007] The finding of mild phenotypic features in the mother of one family that could have been overlooked and the possibility of germline mosaicism in apparently healthy parents in the other family should be taken into account when counseling such families

Overlapping SETBP1 gain-of-function mutations in Schinzel-Giedion syndrome and hematologic malignancies

Schinzel-Giedion syndrome (SGS) is a rare developmental disorder characterized by multiple malformations, severe neurological alterations and increased risk of malignancy. SGS is caused by de novo germline mutations clustering to a 12bp hotspot in exon 4 of SETBP1. Mutations in this hotspot disrupt a degron, a signal for the regulation of protein degradation, and lead to the accumulation of SETBP1 protein. Overlapping SETBP1 hotspot mutations have been observed recurrently as somatic events in leukemia. We collected clinical information of 47 SGS patients (including 26 novel cases) with germline SETBP1 mutations and of four individuals with a milder phenotype caused by de novo germline mutations adjacent to the SETBP1 hotspot. Different mutations within and around the SETBP1 hotspot have varying effects on SETBP1 stability and protein levels in vitro and in in silico modeling. Substitutions in SETBP1 residue I871 result in a weak increase in protein levels and mutations affecting this residue are significantly more frequent in SGS than in leukemia. On the other hand, substitutions in residue D868 lead to the largest increase in protein levels. Individuals with germline mutations affecting D868 have enhanced cell proliferation in vitro and higher incidence of cancer compared to patients with other germline SETBP1 mutations. Our findings substantiate that, despite their overlap, somatic SETBP1 mutations driving malignancy are more disruptive to the degron than germline SETBP1 mutations causing SGS. Additionally, this suggests that the functional threshold for the development of cancer driven by the disruption of the SETBP1 degron is higher than for the alteration in prenatal development in SGS. Drawing on previous studies of somatic SETBP1 mutations in leukemia, our results reveal a genotype-phenotype correlation in germline SETBP1 mutations spanning a molecular, cellular and clinical phenotype

Solve-RD: systematic pan-European data sharing and collaborative analysis to solve rare diseases.

For the first time in Europe hundreds of rare disease (RD) experts team up to actively share and jointly analyse existing patient's data. Solve-RD is a Horizon 2020-supported EU flagship project bringing together >300 clinicians, scientists, and patient representatives of 51 sites from 15 countries. Solve-RD is built upon a core group of four European Reference Networks (ERNs; ERN-ITHACA, ERN-RND, ERN-Euro NMD, ERN-GENTURIS) which annually see more than 270,000 RD patients with respective pathologies. The main ambition is to solve unsolved rare diseases for which a molecular cause is not yet known. This is achieved through an innovative clinical research environment that introduces novel ways to organise expertise and data. Two major approaches are being pursued (i) massive data re-analysis of >19,000 unsolved rare disease patients and (ii) novel combined -omics approaches. The minimum requirement to be eligible for the analysis activities is an inconclusive exome that can be shared with controlled access. The first preliminary data re-analysis has already diagnosed 255 cases form 8393 exomes/genome datasets. This unprecedented degree of collaboration focused on sharing of data and expertise shall identify many new disease genes and enable diagnosis of many so far undiagnosed patients from all over Europe

Solving unsolved rare neurological diseases-a Solve-RD viewpoint.

Funder: Durch Princess Beatrix Muscle Fund Durch Speeren voor Spieren Muscle FundFunder: University of Tübingen Medical Faculty PATE programFunder: European Reference Network for Rare Neurological Diseases | 739510Funder: European Joint Program on Rare Diseases (EJP-RD COFUND-EJP) | 44140962

Solving patients with rare diseases through programmatic reanalysis of genome-phenome data.

Funder: EC | EC Seventh Framework Programm | FP7 Health (FP7-HEALTH - Specific Programme "Cooperation": Health); doi: https://doi.org/10.13039/100011272; Grant(s): 305444, 305444Funder: Ministerio de Economía y Competitividad (Ministry of Economy and Competitiveness); doi: https://doi.org/10.13039/501100003329Funder: Generalitat de Catalunya (Government of Catalonia); doi: https://doi.org/10.13039/501100002809Funder: EC | European Regional Development Fund (Europski Fond za Regionalni Razvoj); doi: https://doi.org/10.13039/501100008530Funder: Instituto Nacional de Bioinformática ELIXIR Implementation Studies Centro de Excelencia Severo OchoaFunder: EC | EC Seventh Framework Programm | FP7 Health (FP7-HEALTH - Specific Programme "Cooperation": Health)Reanalysis of inconclusive exome/genome sequencing data increases the diagnosis yield of patients with rare diseases. However, the cost and efforts required for reanalysis prevent its routine implementation in research and clinical environments. The Solve-RD project aims to reveal the molecular causes underlying undiagnosed rare diseases. One of the goals is to implement innovative approaches to reanalyse the exomes and genomes from thousands of well-studied undiagnosed cases. The raw genomic data is submitted to Solve-RD through the RD-Connect Genome-Phenome Analysis Platform (GPAP) together with standardised phenotypic and pedigree data. We have developed a programmatic workflow to reanalyse genome-phenome data. It uses the RD-Connect GPAP's Application Programming Interface (API) and relies on the big-data technologies upon which the system is built. We have applied the workflow to prioritise rare known pathogenic variants from 4411 undiagnosed cases. The queries returned an average of 1.45 variants per case, which first were evaluated in bulk by a panel of disease experts and afterwards specifically by the submitter of each case. A total of 120 index cases (21.2% of prioritised cases, 2.7% of all exome/genome-negative samples) have already been solved, with others being under investigation. The implementation of solutions as the one described here provide the technical framework to enable periodic case-level data re-evaluation in clinical settings, as recommended by the American College of Medical Genetics

Twist exome capture allows for lower average sequence coverage in clinical exome sequencing

Background Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques. Results We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection. Conclusion We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques

A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing

Purpose Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the “ClinVar low-hanging fruit” reanalysis, reasons for the failure of previous analyses, and lessons learned. Methods Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. Results We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). Conclusion The “ClinVar low-hanging fruit” analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock

Chromosomal abnormalities in hepatic cysts point to novel polycystic liver disease genes

Autosomal dominant polycystic liver disease (ADPLD) is caused by variants in PRKCSH, SEC63, and LRP5, whereas autosomal dominant polycystic kidney disease is caused by variants in PKD1 and PKD2. Liver cyst development in these disorders is explained by somatic loss-of-heterozygosity (LOH) of the wild-type allele in the developing cyst. We hypothesize that we can use this mechanism to identify novel disease genes that reside in LOH regions. In this study, we aim to map abnormal genomic regions using high-density SNP microarrays to find novel PLD genes. We collected 46 cysts from 23 patients with polycystic or sporadic hepatic cysts, and analyzed DNA from those cysts using high-resolution microarray (n=24) or Sanger sequencing (n=22). We here focused on regions of homozygosity on the autosomes (>3.0 Mb) and large CNVs (>1.0 Mb). We found frequent LOH in PRKCSH (22/29) and PKD1/PKD2 (2/3) cysts of patients with known heterozygous germline variants in the respective genes. In the total cohort, 12/23 patients harbored abnormalities outside of familiar areas. In individual ADPLD cases, we identified germline events: a 2q13 complex rearrangement resulting in BUB1 haploinsufficiency, a 47XXX karyotype, chromosome 9q copy-number loss, and LOH on chromosome 3p. The latter region was overlapping with an LOH region identified in two other cysts. Unique germline and somatic abnormalities occur frequently in and outside of known genes underlying cysts. Each liver cyst has a unique genetic makeup. LOH driver gene BUB1 may imply germline causes of genetic instability in PLD