Diminished TMEM100 Expression in a Newborn With Acinar Dysplasia and a Novel TBX4 Variant: A Case Report

Acinar dysplasia (AcDys) of the lung is a rare lethal developmental disorder in neonates characterized by severe respiratory failure and pulmonary arterial hypertension refractory to treatment. Recently, abnormalities of TBX4-FGF10-FGFR2-TMEM100 signaling regulating lung development have been reported in patients with AcDys due to heterozygous single-nucleotide variants or copy-number variant deletions involving TBX4, FGF10, or FGFR2. Here, we describe a female neonate who died at 4 hours of life due to severe respiratory distress related to AcDys diagnosed by postmortem histopathologic evaluation. Genomic analyses revealed a novel deleterious heterozygous missense variant c.728A>C (p.Asn243Thr) in TBX4 that arose de novo on paternal chromosome 17. We also identified 6 candidate hypomorphic rare variants in the TBX4 enhancer in trans to TBX4 coding variant. Gene expression analyses of proband's lung tissue showed a significant reduction of TMEM100 expression with near absence of TMEM100 within the endothelium of arteries and capillaries by immunohistochemistry. These results support the pathogenicity of the detected TBX4 variant and provide further evidence that disrupted signaling between TBX4 and TMEM100 may contribute to severe lung phenotypes in humans, including AcDys.

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

SNV/indel hypermutator phenotype in biallelic RAD51C variant: Fanconi anemia

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C. Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5 and 37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the hypothesis of a rare phenomenon of hypermutation. Eight-hundred-thirty apparent de novo SNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component having an average VAF of 25% (range: 18.7-37.41%) (n = 446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range 9.55-18.69%) (n = 384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA. We propose that other FA groups should be investigated for genome-wide de novo variants.

SNV/indel hypermutator phenotype in biallelic RAD51C variant - Fanconi anemia

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C . Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5%-37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the rare phenomenon of hypermutation. Eight-hundred-thirty apparent dnSNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component with an average VAF of 25% (range: 18.7-37.41%) (n=446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range: 9.55-18.69%) (n=384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA, and pregnancy loss in general. We propose that other FA groups should be investigated for genome-wide de novo variants.

Cloud-native distributed genomic pileup operations

MOTIVATION

Pileup analysis is a building block of many bioinformatics pipelines, including variant calling and genotyping. This step tends to become a bottleneck of the entire assay since the straightforward pileup implementations involve processing of all base calls from all alignments sequentially. On the other hand, a distributed version of the algorithm faces the intrinsic challenge of splitting reads-oriented file formats into self-contained partitions to avoid costly data exchange between computational nodes.

RESULTS

Here, we present a scalable, distributed and efficient implementation of a pileup algorithm that is suitable for deploying in cloud computing environments. In particular, we implemented: (i) our custom data-partitioning algorithm optimized to work with the alignment reads, (ii) a novel and unique approach to process alignment events from sequencing reads using the MD tags, (iii) the source code micro-optimizations for recurrent operations, and (iv) a modular structure of the algorithm. We have proven that our novel approach consistently and significantly outperforms other state-of-the-art distributed tools in terms of execution time (up to 6.5× faster) and memory usage (up to 2× less), resulting in a substantial cloud cost reduction. SeQuiLa is a cloud-native solution that can be easily deployed using any managed Kubernetes and Hadoop services available in public clouds, like Microsoft Azure Cloud, Google Cloud Platform, or Amazon Web Services. Together with the already implemented distributed range join and coverage calculations, our package provides end-users with a unified SQL interface for convenient analyses of population-scale genomic data in an interactive way.

AVAILABILITY AND IMPLEMENTATION

https://biodatageeks.github.io/sequila/.

Variable expressivity in a four-generation ACDMPV family with a non-coding hypermorphic SNV in trans to the frameshifting FOXF1 variant

Heterozygous single nucleotide variants (SNVs) or copy-number variant deletions involving FOXF1 or its distant lung-specific enhancer on chromosome 16q24.1 have been identified in 80-90% of patients with Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a lethal neonatal lung developmental disorder. We describe a four-generation family with a deceased ACDMPV neonate, her sibling from the electively terminated pregnancy, healthy mother with a history of pulmonary arterial hypertension (PAH), an unaffected aunt, an aunt deceased due to findings consistent with ACDMPV, and a reportedly unaffected grandmother, all with the frameshifting variant c.881_902dup (p.Gly302Profs*46) in FOXF1, and a deceased great-grandmother with a history of PAH. Genome sequencing analyses in the proband's unaffected mother revealed a non-coding putative regulatory SNV rs560517434-A within the lung-specific distant FOXF1 enhancer in trans to the FOXF1 frameshift mutation. Functional testing of this variant using an in vitro luciferase reporter assay showed that it increased FOXF1 promoter activity 10-fold. Our studies further demonstrate that non-coding SNVs in the FOXF1 enhancer region can rescue the lethal ACDMPV phenotype and support the compound inheritance gene dosage model.

High-level gonosomal mosaicism for a pathogenic non-coding CNV deletion of the lung-specific FOXF1 enhancer in an unaffected mother of an infant with ACDMPV

BACKGROUND

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) results from haploinsufficiency of the mesenchymal transcription factor FOXF1 gene. To date, only one case of an ACDMPV-causative CNV deletion inherited from a very-low level somatic mosaic mother has been reported.

METHODS

Clinical, histopathological, and molecular studies, including whole genome sequencing, chromosomal microarray analysis, qPCR, and Sanger sequencing, followed by in vitro fertilization (IVF) with preimplantation genetic testing (PGT) were used to study a family with a deceased neonate with ACDMPV.

RESULTS

A pathogenic CNV deletion of the lung-specific FOXF1 enhancer in the proband was found to be inherited from an unaffected mother, 36% mosaic for this deletion in her peripheral blood cells. The qPCR analyses of saliva, buccal cells, urine, nail, and hair samples revealed 19%, 18%, 15%, 19%, and 27% variant allele fraction, respectively, indicating a high recurrence risk. Grandparental studies revealed that the deletion arose on the mother's paternal chromosome 16. PGT studies revealed 44% embryos with the deletion, reflecting high-level germline mosaicism.

CONCLUSION

Our data further demonstrate the importance of parental testing in ACDMPV families and reproductive usefulness of IVF with PGT in families with high-level parental gonosomal mosaicism.

Loss-of-function variant in chymotrypsin like elastase 3B (CELA3B) is associated with non-alcoholic chronic pancreatitis

BACKGROUND

Genetic alterations in digestive enzymes have been associated with chronic pancreatitis (CP). Recently, chymotrypsin like elastase 3B (CELA3B) emerged as a novel risk gene. Thus, we evaluated CELA3B in two European cohorts with CP.

METHODS

We analyzed all 8 CELA3B exons in 550 German non-alcoholic CP (NACP) patients and in 241 German controls by targeted DNA sequencing. In addition, we analyzed exons 6 and 7 by Sanger sequencing and the c.129+1G>A variant by melting curve analysis in 1078 further German controls. As replication cohort, we investigated up to 243 non-German European NACP patients and up to 1665 controls originating from Poland, Hungary, and Sweden. We assessed the cellular secretion and the elastase activity of recombinant CELA3B variants.

RESULTS

In the German discovery cohort, we detected a splice-site variant in intron 2, c.129+1G>A, in 9/550 (1.64%) CP patients and in 5/1319 (0.38%) controls (P=0.007, OR=4.4, 95% CI=1.5-13.0). In the European replication cohort, this variant was also enriched in patients (9/178 [5.06%]) versus controls (13/1247 [1.04%]) (P=0.001, OR=5.1, 95% CI=2.1-12.0). We did not find the two previously reported codon 90 variants, p.R90C and p.R90L.

CONCLUSIONS

Our data indicate that CELA3B is a susceptibility gene for CP. In contrast to previous reports suggesting that increased CELA3B activity is associated with CP risk, the splice-site variant identified here is predicted to cause diminished CELA3B expression. How reduced CELA3B function predisposes to pancreatitis remains to be elucidated.

Postzygotic mutations and where to find them - Recent advances and future implications in the field of non-neoplastic somatic mosaicism

The technological progress of massively parallel sequencing (MPS) has triggered a remarkable development in the research on postzygotic mutations. Although the overwhelming majority of studies in the field focus on oncogenesis, non-neoplastic diseases are attracting more and more attention. The aim of this review was to summarize some of the most recent findings in the field of somatic mosaicism in diseases other than neoplastic events. We discuss the abundance and role of postzygotic mutations, with a special emphasis on disorders which occur only in a mosaic form (obligatory mosaic diseases; OMDs). Based on the list of OMDs compiled from the published literature and three databases (OMIM, Orphanet and MosaicBase), we demonstrate the prevalence of cancer-related genes across OMDs and suggest other sources to further explore OMDs and OMD-related genes. Additionally, we comment on some practical aspects related to mosaic diseases, such as approaches to tissue sampling, the MPS coverage required to detect variants at a very low frequency, as well as on bioinformatic and molecular tools dedicated to detect somatic mutations in MPS data.

Variants in the pancreatic CUB and zona pellucida-like domains 1 (CUZD1) gene in early-onset chronic pancreatitis - A possible new susceptibility gene

OBJECTIVE

Non-alcoholic chronic pancreatitis (NACP) frequently develops in the setting of genetic susceptibility associated with alterations in genes that are highly expressed in the pancreas. However, the genetic basis of NACP remains unresolved in a significant number of patients warranting a search for further risk genes.

DESIGN

We analyzed CUZD1, which encodes the CUB and zona pellucida-like domains 1 protein that is found in high levels in pancreatic acinar cells. We sequenced the coding region in 1163 European patients and 2018 European controls. In addition, we analyzed 297 patients and 1070 controls from Japan. We analyzed secretion of wild-type and mutant CUZD1 from transfected cells using Western blotting.

RESULTS

In the European cohort, we detected 30 non-synonymous variants. Using different prediction tools (SIFT, CADD, PROVEAN, PredictSNP) or the combination of these tools, we found accumulation of predicted deleterious variants in patients (p-value range 0.002-0.013; OR range 3.1-5.2). No association was found in the Japanese cohort, in which 13 non-synonymous variants were detected. Functional studies revealed >50% reduced secretion of 7 variants, however, these variants were not significantly enriched in European CP patients.

CONCLUSION

Our data indicate that CUZD1 might be a novel susceptibility gene for NACP. How these variants predispose to pancreatitis remains to be elucidated.

Do paternal deletions involving the FOXF1 locus on chromosome 16q24.1 manifest with more severe non-lung anomalies?

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disorder in neonates due to heterozygous loss-of-function of the mesenchymal transcription factor gene, FOXF1. Interestingly, unlike ACDMPV-causing point mutations in FOXF1 that can be inherited from the mother or father, causative copy-number variant (CNV) deletions arise de novo and almost exclusively on chromosome 16 inherited from the mother (n = 50 vs. n = 3). Here, we describe a fourth case of a de novo paternal CNV deletion encompassing FOXF1, its neighboring long non-coding RNA gene FENDRR, and their distant lung-specific enhancer, identified in a 21-week-old fetus with tetralogy of Fallot, gastrointestinal and genitourinary abnormalities, a single umbilical artery, and patchy histopathological findings of ACDMPV in autopsy lung. We also review the ACDMPV-causative CNV deletions detected prenatally and propose that the majority of paternal deletions manifest with more severe additional non-lung abnormalities.

Exacerbation of mild lung disorders to lethal pulmonary hypoplasia by a noncoding hypomorphic SNV in a lung-specific enhancer in trans to the frameshifting TBX4 variant

Variants involving TBX4 are associated with a wide variety of disorders, including pulmonary arterial hypertension, ischiocoxopodopatellar syndrome (ICPPS)/small patella syndrome (SPS), lethal lung developmental disorders (LLDDs) in neonates, heart defects, and prenatally lethal posterior amelia with pelvic and pulmonary hypoplasia syndrome. The objective of our study was to elucidate the wide variable phenotypic expressivity and incomplete penetrance in a three-generation family with a truncating variant in TBX4. In addition to exome and genome sequencing analyses, a candidate noncoding regulatory single nucleotide variant (SNV) within the lung-specific TBX4 enhancer was functionally tested using an in vitro luciferase reporter assay. A heterozygous frameshift variant c.1112dup (p.Pro372Serfs*14) in TBX4 was identified in patients with mild interstitial lung disease (1), bronchiolitis obliterans (1), recurrent pneumothorax (1), ICPPS/SPS (1), LLDD (2), and in unaffected individuals (4). In two deceased neonates with LLDD, we identified a noncoding SNV rs62069651-C located in trans to the mutated TBX4 allele that reduced the TBX4 promoter activity by 63% in the reporter assay. Our findings provide a functional evidence for the recently reported model of complex compound inheritance in which both TBX4 coding and in trans noncoding hypomorphic variants in the lung-specific enhancer of TBX4 contribute to LLDD.

Detection of low-level parental somatic mosaicism for clinically relevant SNVs and indels identified in a large exome sequencing dataset

BACKGROUND

Due to the limitations of the current routine diagnostic methods, low-level somatic mosaicism with variant allele fraction (VAF) < 10% is often undetected in clinical settings. To date, only a few studies have attempted to analyze tissue distribution of low-level parental mosaicism in a large clinical exome sequencing (ES) cohort.

METHODS

Using a customized bioinformatics pipeline, we analyzed apparent de novo single-nucleotide variants or indels identified in the affected probands in ES trio data at Baylor Genetics clinical laboratories. Clinically relevant variants with VAFs between 30 and 70% in probands and lower than 10% in one parent were studied. DNA samples extracted from saliva, buccal cells, redrawn peripheral blood, urine, hair follicles, and nail, representing all three germ layers, were tested using PCR amplicon next-generation sequencing (amplicon NGS) and droplet digital PCR (ddPCR).

RESULTS

In a cohort of 592 clinical ES trios, we found 61 trios, each with one parent suspected of low-level mosaicism. In 21 parents, the variants were validated using amplicon NGS and seven of them by ddPCR in peripheral blood DNA samples. The parental VAFs in blood samples varied between 0.08 and 9%. The distribution of VAFs in additional tissues ranged from 0.03% in hair follicles to 9% in re-drawn peripheral blood.

CONCLUSIONS

Our study illustrates the importance of analyzing ES data using sensitive computational and molecular methods for low-level parental somatic mosaicism for clinically relevant variants previously diagnosed in routine clinical diagnostics as apparent de novo.

Loss of function TRPV6 variants are associated with chronic pancreatitis in nonalcoholic early-onset Polish and German patients

PURPOSE

Loss of function variants of the transient receptor potential cation channel, subfamily V, member 6 (TRPV6) have been recently associated with chronic pancreatitis (CP) in Japanese, German and French patients. Here, we investigated the association of TRPV6 variants with CP in independent European cohorts of early-onset CP patients from Poland and Germany.

PATIENTS AND METHODS

We enrolled 152 pediatric CP patients (median age 8.6 yrs) with no history of alcohol/smoking abuse and 472 controls from Poland as well as 157 nonalcoholic young CP patients (median age 20 yrs) and 750 controls from Germany. Coding regions of TRPV6 were screened by Sanger and next generation sequencing. Selected, potentially pathogenic TRPV6 variants were expressed in HEK293T cells and TRPV6 activity was analyzed using ratiometric Ca²⁺ measurements.

RESULTS

Overall, we identified 10 novel (3 nonsense and 7 missenses) TRPV6 variants in CP patients. TRPV6 p.V239SfsX53 nonsense variant and the variants showing significant decrease in intracellular Ca²⁺ concentration in HEK293T cells (p.R174X, p.L576R, p.R342Q), were significantly overrepresented in Polish patients as compared to controls (6/152, 3.9% vs. 0/358, 0%; P = 0,0007). Nonsense TRPV6 variants predicted as loss of function (p.V239SfsX53 and p.R624X) were also significantly overrepresented in German patients (3/157; 2.0% vs 0/750; 0%, P = 0.005).

CONCLUSIONS

We showed that TRPV6 loss of function variants are associated with elevated CP risk in early-onset Polish and German patients confirming that TRPV6 is a novel CP susceptibility gene.

Perturbation of semaphorin and VEGF signaling in ACDMPV lungs due to FOXF1 deficiency

BACKGROUND

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal congenital lung disorder in neonates characterized by severe progressive respiratory failure and refractory pulmonary hypertension, resulting from underdevelopment of the peripheral pulmonary tree. Causative heterozygous single nucleotide variants (SNVs) or copy-number variant (CNV) deletions involving FOXF1 or its distant lung-specific enhancer on chromosome 16q24.1 have been identified in 80-90% of ACDMPV patients. FOXF1 maps closely to and regulates the oppositely oriented FENDRR, with which it also shares regulatory elements.

METHODS

To better understand the transcriptional networks downstream of FOXF1 that are relevant for lung organogenesis, using RNA-seq, we have examined lung transcriptomes in 12 histopathologically verified ACDMPV patients with or without pathogenic variants in the FOXF1 locus and analyzed gene expression profile in FENDRR-depleted fetal lung fibroblasts, IMR-90.

RESULTS

RNA-seq analyses in ACDMPV neonates revealed changes in the expression of several genes, including semaphorins (SEMAs), neuropilin 1 (NRP1), and plexins (PLXNs), essential for both epithelial branching and vascular patterning. In addition, we have found deregulation of the vascular endothelial growth factor (VEGF) signaling that also controls pulmonary vasculogenesis and a lung-specific endothelial gene TMEM100 known to be essential in vascular morphogenesis. Interestingly, we have observed a substantial difference in gene expression profiles between the ACDMPV samples with different types of FOXF1 defect. Moreover, partial overlap between transcriptome profiles of ACDMPV lungs with FOXF1 SNVs and FENDRR-depleted IMR-90 cells suggests contribution of FENDRR to ACDMPV etiology.

CONCLUSIONS

Our transcriptomic data imply potential crosstalk between several lung developmental pathways, including interactions between FOXF1-SHH and SEMA-NRP or VEGF/VEGFR2 signaling, and provide further insight into complexity of lung organogenesis in humans.

Floppy infant syndrome as a first manifestation of LMNA-related congenital muscular dystrophy

LMNA-related congenital muscular dystrophy (L-CMD) is the most severe phenotypic form of skeletal muscle laminopathies. This paper reports clinical presentation of the disease in 15 Polish patients from 13 families with genetically confirmed skeletal muscle laminopathy. In all these patients floppy infant syndrome was the first manifestation of the disease. The genetic diagnosis was established by next generation sequencing (targeted panel or exome; 11 patients) or classic Sanger sequencing (4 patients). In addition to known pathogenic LMNA variants: c.116A > G (p.Asn39Ser), c.745C > T (p.Arg249Trp), c.746G > A (p.Arg249Gln), c.1072G > A (p.Glu358Lys), c.1147G > A (p.Glu383Lys), c.1163G > C (p.Arg388Pro), c.1357C > T (p.Arg453Trp), c.1583C > G (p.Thr528Arg), we have identified three novel ones: c.121C > G (p.Arg41Gly), c.1127A > G (p.Tyr376Cys) and c.1160T > C (p.Leu387Pro). Eleven patients had de novo mutations, 4 - familial. In one family we observed intrafamilial variability of clinical course: severe L-CMD in the male proband, intermediate form in his sister and asymptomatic in their mother. One asymptomatic father had somatic mosaicism. L-CMD should be suspected in children with hypotonia in infancy and delayed motor development, who have poor head control, severe hyperlordosis and unstable and awkward gait. Serum creatine kinase may be high (~1000IU/l). Progression of muscle weakness is fast, leading to early immobilization. In some patients with L-CMD joint contractures can develop with time. MRI shows that the most frequently affected muscles are the serratus anterior, lumbar paraspinal, gluteus, vastus, adductor magnus, hamstrings, medial head of gastrocnemius and soleus. Ultra-rare laminopathies can be a relatively common cause of generalized hypotonia in children. Introduction of wide genome sequencing methods was a breakthrough in diagnostics of diseases with great clinical and genetic variability and allowed approach "from genotype do phenotype". However target sequencing of LMNA gene could be considered in selected patients with clinical picture suggestive for laminopathy.

Lung-specific distant enhancer cis regulates expression of FOXF1 and lncRNA FENDRR

The FOXF1 gene, causative for a neonatal lethal lung developmental disorder alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), maps 1.7 kb away from the long noncoding RNA gene FENDRR on the opposite strand, suggesting they may be coregulated. Using RNA sequencing in lung tissue from ACDMPV patients with heterozygous deletions of the FOXF1 distant enhancer located 286 kb upstream, leaving FOXF1 and FENDRR intact, we have found that the FENDRR and FOXF1 expressions were reduced by approximately 75% and 50%, respectively, and were monoallelic from the intact chromosome 16q24.1. In contrast, ACDMPV patients with FOXF1 SNVs had biallelic FENDRR expression reduced by 66%-82%. Corroboratively, depletion of FOXF1 by small interfering RNA in lung fibroblasts resulted in a 50% decrease of FENDRR expression. These data indicate that FENDRR expression in the lungs is regulated both in cis by the FOXF1 distant enhancer and in trans by FOXF1. Our findings are compatible with the involvement of FENDRR in FOXF1-related disorders, including ACDMPV.

Potential interactions between the TBX4-FGF10 and SHH-FOXF1 signaling during human lung development revealed using ChIP-seq

Background

The epithelial-mesenchymal signaling involving SHH-FOXF1, TBX4-FGF10, and TBX2 pathways is an essential transcriptional network operating during early lung organogenesis. However, precise regulatory interactions between different genes and proteins in this pathway are incompletely understood.

Methods

To identify TBX2 and TBX4 genome-wide binding sites, we performed chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) in human fetal lung fibroblasts IMR-90.

Results

We identified 14,322 and 1,862 sites strongly-enriched for binding of TBX2 and TBX4, respectively, 43.95% and 18.79% of which are located in the gene promoter regions. Gene Ontology, pathway enrichment, and DNA binding motif analyses revealed a number of overrepresented cues and transcription factor binding motifs relevant for lung branching that can be transcriptionally regulated by TBX2 and/or TBX4. In addition, TBX2 and TBX4 binding sites were found enriched around and within FOXF1 and its antisense long noncoding RNA FENDRR, indicating that the TBX4-FGF10 cascade may directly interact with the SHH-FOXF1 signaling.

Conclusions

We highlight the complexity of transcriptional network driven by TBX2 and TBX4 and show that disruption of this crosstalk during morphogenesis can play a substantial role in etiology of lung developmental disorders.

Integrated sequencing and array comparative genomic hybridization in familial Parkinson disease

Objective

To determine how single nucleotide variants (SNVs) and copy number variants (CNVs) contribute to molecular diagnosis in familial Parkinson disease (PD), we integrated exome sequencing (ES) and genome-wide array-based comparative genomic hybridization (aCGH) and further probed CNV structure to reveal mutational mechanisms.

Methods

We performed ES on 110 subjects with PD and a positive family history; 99 subjects were also evaluated using genome-wide aCGH. We interrogated ES and aCGH data for pathogenic SNVs and CNVs at Mendelian PD gene loci. We confirmed SNVs via Sanger sequencing and further characterized CNVs with custom-designed high-density aCGH, droplet digital PCR, and breakpoint sequencing.

Results

Using ES, we discovered individuals with known pathogenic SNVs in GBA (p.Glu365Lys, p.Thr408Met, p.Asn409Ser, and p.Leu483Pro) and LRRK2 (p.Arg1441Gly and p.Gly2019Ser). Two subjects were each double heterozygotes for variants in GBA and LRRK2. Based on aCGH, we additionally discovered cases with an SNCA duplication and heterozygous intragenic GBA deletion. Five additional subjects harbored both SNVs (p.Asn52Metfs*29, p.Thr240Met, p.Pro437Leu, and p.Trp453*) and likely disrupting CNVs at the PRKN locus, consistent with compound heterozygosity. In nearly all cases, breakpoint sequencing revealed microhomology, a mutational signature consistent with CNV formation due to DNA replication errors.

Conclusions

Integrated ES and aCGH yielded a genetic diagnosis in 19.3% of our familial PD cohort. Our analyses highlight potential mechanisms for SNCA and PRKN CNV formation, uncover multilocus pathogenic variation, and identify novel SNVs and CNVs for further investigation as potential PD risk alleles.

Low-level parental somatic mosaic SNVs in exomes from a large cohort of trios with diverse suspected Mendelian conditions

Purpose:

The goal of this study was to assess the scale of low-level parental mosaicism in exome sequencing (ES) databases.

Methods:

We analyzed approximately 2000 family trio ES datasets from the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) and Baylor Genetics (BG). Among apparent de novo single nucleotide variants (SNVs) identified in the affected probands, we selected rare unique variants with variant allele fraction (VAF) between 30-70% in the probands and lower than 10% in one of the parents.

Results:

Out of 102 candidate mosaic variants validated using amplicon-based NGS, droplet digital PCR, or blocker displacement amplification, 27 (26.4%) were confirmed to be low- (VAF between 1-10%) or very low- (VAF <1%) level mosaic. Detection precision in parental samples with two or more alternate reads was 63.6% (BHCMG) and 43.6% (BG). In nine investigated individuals, we observed variability of mosaic ratios among blood, saliva, fibroblast, buccal, hair, and urine samples.

Conclusion:

Our computational pipeline enables robust discrimination between true and false positive candidate mosaic variants and efficient detection of low-level mosaicism in ES samples. We confirm that the presence of two or more alternate reads in the parental sample is a reliable predictor of low-level parental somatic mosaicism.

SeQuiLa: an elastic, fast and scalable SQL-oriented solution for processing and querying genomic intervals

SUMMARY

Efficient processing of large-scale genomic datasets has recently become possible due to the application of 'big data' technologies in bioinformatics pipelines. We present SeQuiLa-a distributed, ANSI SQL-compliant solution for speedy querying and processing of genomic intervals that is available as an Apache Spark package. Proposed range join strategy is significantly (∼22×) faster than the default Apache Spark implementation and outperforms other state-of-the-art tools for genomic intervals processing.

AVAILABILITY AND IMPLEMENTATION

The project is available at http://biodatageeks.org/sequila/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Wolff-Parkinson-White syndrome: De novo variants and evidence for mutational burden in genes associated with atrial fibrillation

BACKGROUND

Wolff-Parkinson-White (WPW) syndrome is a relatively common arrhythmia affecting ~1-3/1,000 individuals. Mutations in PRKAG2 have been described in rare patients in association with cardiomyopathy. However, the genetic basis of WPW in individuals with a structurally normal heart remains poorly understood. Sudden death due to atrial fibrillation (AF) can also occur in these individuals. Several studies have indicated that despite ablation of an accessory pathway, the risk of AF remains high in patients compared to general population.

METHODS

We applied exome sequencing in 305 subjects, including 65 trios, 80 singletons, and 6 multiple affected families. We used de novo analysis, candidate gene approach, and burden testing to explore the genetic contributions to WPW.

RESULTS

A heterozygous deleterious variant in PRKAG2 was identified in one subject, accounting for 0.6% (1/151) of the genetic basis of WPW in this study. Another individual with WPW and left ventricular hypertrophy carried a known pathogenic variant in MYH7. We found rare de novo variants in genes associated with arrhythmia and cardiomyopathy (ANK2, NEBL, PITX2, and PRDM16) in this cohort. There was an increased burden of rare deleterious variants (MAF ≤ 0.005) with CADD score ≥ 25 in genes linked to AF in cases compared to controls (P = .0023).

CONCLUSIONS

Our findings show an increased burden of rare deleterious variants in genes linked to AF in WPW syndrome, suggesting that genetic factors that determine the development of accessory pathways may be linked to an increased susceptibility of atrial muscle to AF in a subset of patients.

Accumulation of sequence variants in genes of Wnt signaling and focal adhesion pathways in human corneas further explains their involvement in keratoconus

Background

Keratoconus (KTCN) is a protrusion and thinning of the cornea, resulting in loss of visual acuity. The etiology of KTCN remains unclear. The purpose of this study was to assess the potential involvement of new genetic variants in KTCN etiology based on both the genomic and transcriptomic findings recognized in the same corneal tissues.

Methods

Corneal tissues derived from five unrelated Polish individuals with KTCN were examined using exome sequencing (ES), followed by enrichment analyses. For comparison purposes, the datasets comprising ES data of five randomly selected Polish individuals without ocular abnormalities and five Polish patients with high myopia were used. Expression levels of selected genes from the overrepresented pathways were obtained from the previous RNA-Seq study.

Results

Exome capture discovered 117 potentially relevant variants that were further narrowed by gene overrepresentation analyses. In each of five patients, the assessment of functional interactions revealed rare (MAF ≤ 0.01) DNA variants in at least one gene from Wnt signaling (VANGL1, WNT1, PPP3CC, LRP6, FZD2) and focal adhesion (BIRC2, PAK6, COL4A4, PPP1R12A, PTK6) pathways. No genes involved in pathways enriched in KTCN corneas were overrepresented in our control sample sets.

Conclusions

The results of this first pilot ES profiling of human KTCN corneas emphasized that accumulation of sequence variants in several genes from Wnt signaling and/or focal adhesion pathways might cause the phenotypic effect and further points to a complex etiology of KTCN.

Clinical, Histopathological, and Molecular Diagnostics in Lethal Lung Developmental Disorders

Lethal lung developmental disorders are a rare but important group of pediatric diffuse lung diseases presenting with neonatal respiratory failure. On the basis of histopathological appearance at lung biopsy or autopsy, they have been termed: alveolar capillary dysplasia with misalignment of the pulmonary veins, acinar dysplasia, congenital alveolar dysplasia, and other unspecified primary pulmonary hypoplasias. However, the histopathological continuum in these lethal developmental disorders has made accurate diagnosis challenging, which has implications for recurrence risk. Over the past decade, genetic studies in infants with alveolar capillary dysplasia with misalignment of the pulmonary veins have revealed the causative role of the dosage-sensitive FOXF1 gene and its noncoding regulatory variants in the distant lung-specific enhancer at chromosome 16q24.1. In contrast, the molecular bases of acinar dysplasia and congenital alveolar dysplasia have remained poorly understood. Most recently, disruption of the TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling pathway has been reported in patients with these lethal pulmonary dysplasias. Application of next-generation sequencing techniques, including exome sequencing and whole-genome sequencing, has demonstrated their complex compound inheritance. These data indicate that noncoding regulatory elements play a critical role in lung development in humans. We propose that for more precise lethal lung developmental disorder diagnosis, a diagnostic pathway including whole-genome sequencing should be implemented.

A de novo 2.2 Mb recurrent 17q23.1q23.2 deletion unmasks novel putative regulatory non-coding SNVs associated with lethal lung hypoplasia and pulmonary hypertension: a case report

BACKGROUND

Application of whole genome sequencing (WGS) enables identification of non-coding variants that play a phenotype-modifying role and are undetectable by exome sequencing. Recently, non-coding regulatory single nucleotide variants (SNVs) have been reported in patients with lethal lung developmental disorders (LLDDs) or congenital scoliosis with recurrent copy-number variant (CNV) deletions at 17q23.1q23.2 or 16p11.2, respectively.

CASE PRESENTATION

Here, we report a deceased newborn with pulmonary hypertension and pulmonary interstitial emphysema with features suggestive of pulmonary hypoplasia, resulting in respiratory failure and neonatal death soon after birth. Using the array comparative genomic hybridization and WGS, two heterozygous recurrent CNV deletions: ~ 2.2 Mb on 17q23.1q23.2, involving TBX4, and ~ 600 kb on 16p11.2, involving TBX6, that both arose de novo on maternal chromosomes were identified. In the predicted lung-specific enhancer upstream to TBX4, we have detected seven novel putative regulatory non-coding SNVs that were absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies.

CONCLUSIONS

Our findings further support a recently reported model of complex compound inheritance of LLDD in which both non-coding and coding heterozygous TBX4 variants contribute to the lung phenotype. In addition, this is the first report of a patient with combined de novo heterozygous recurrent 17q23.1q23.2 and 16p11.2 CNV deletions.

Bi-allelic Pathogenic Variants in TUBGCP2 Cause Microcephaly and Lissencephaly Spectrum Disorders

Lissencephaly comprises a spectrum of malformations of cortical development. This spectrum includes agyria, pachygyria, and subcortical band heterotopia; each represents anatomical malformations of brain cortical development caused by neuronal migration defects. The molecular etiologies of neuronal migration anomalies are highly enriched for genes encoding microtubules and microtubule-associated proteins, and this enrichment highlights the critical role for these genes in cortical growth and gyrification. Using exome sequencing and family based rare variant analyses, we identified a homozygous variant (c.997C>T [p.Arg333Cys]) in TUBGCP2, encoding gamma-tubulin complex protein 2 (GCP2), in two individuals from a consanguineous family; both individuals presented with microcephaly and developmental delay. GCP2 forms the multiprotein γ-tubulin ring complex (γ-TuRC) together with γ-tubulin and other GCPs to regulate the assembly of microtubules. By querying clinical exome sequencing cases and through GeneMatcher-facilitated collaborations, we found three additional families with bi-allelic variation and similarly affected phenotypes including a homozygous variant (c.1843G>C [p.Ala615Pro]) in two families and compound heterozygous variants consisting of one missense variant (c.889C>T [p.Arg297Cys]) and one splice variant (c.2025-2A>G) in another family. Brain imaging from all five affected individuals revealed varying degrees of cortical malformations including pachygyria and subcortical band heterotopia, presumably caused by disruption of neuronal migration. Our data demonstrate that pathogenic variants in TUBGCP2 cause an autosomal recessive neurodevelopmental trait consisting of a neuronal migration disorder, and our data implicate GCP2 as a core component of γ-TuRC in neuronal migrating cells.

A Patient with Berardinelli-Seip Syndrome, Novel AGPAT2 Splicesite Mutation and Concomitant Development of Non-diabetic Polyneuropathy

Primary polyneuropathy in the context of Seip-Berardinelli type 1 seipinopathy, or congenital generalized lipodystrophy type 1 (CGL1) has not been previously reported. We report the case history of a 27 year old female CGL1 patient presenting with an unusual additional development of non-diabetic peripheral neuropathy and learning disabilities in early adolescence. Whole exome sequencing (WES) of the patient genome identified a novel variant, homozygous for a 52 bp intronic deletion in the AGPAT2 locus, coding for 1-acylglycerol-3-phosphate O-acyltransferase 2, which is uniquely associated with CGL1 seipinopathies, with no molecular evidence for dual diagnosis. Functional studies using RNA isolated from patient peripheral blood leucocytes showed abnormal RNA splicing resulting in the loss of 25 amino acids from the patient AGPAT2 protein coding sequence. Stability and transcription levels for the misspliced AGPAT2 mRNA in our patient nonetheless remained normal. Any AGPAT2 protein produced in our patient is therefore likely to be dysfunctional. However, formal linkage of this deletion to the neuropathy observed remains to be shown. The classical clinical presentation of a patient with AGPAT2-associated lipodystrophy shows normal cognition and no development of polyneuropathy. Cognitive disabilities and polyneuropathy are features associated exclusively with clinical CGL type 2 arising from seipin (BSCL2) gene mutations. This case study suggests that in some genetic contexts, AGPAT2 mutations can also produce phenotypes with primary polyneuropathy.

SeQuiLa-cov: A fast and scalable library for depth of coverage calculations

BACKGROUND

Depth of coverage calculation is an important and computationally intensive preprocessing step in a variety of next-generation sequencing pipelines, including the analysis of RNA-sequencing data, detection of copy number variants, or quality control procedures.

RESULTS

Building upon big data technologies, we have developed SeQuiLa-cov, an extension to the recently released SeQuiLa platform, which provides efficient depth of coverage calculations, reaching >100× speedup over the state-of-the-art tools. The performance and scalability of our solution allow for exome and genome-wide calculations running locally or on a cluster while hiding the complexity of the distributed computing with Structured Query Language Application Programming Interface.

CONCLUSIONS

SeQuiLa-cov provides significant performance gain in depth of coverage calculations streamlining the widely used bioinformatic processing pipelines.

A clinical survey of mosaic single nucleotide variants in disease-causing genes detected by exome sequencing

Background

Although mosaic variation has been known to cause disease for decades, high-throughput sequencing technologies with the analytical sensitivity to consistently detect variants at reduced allelic fractions have only recently emerged as routine clinical diagnostic tests. To date, few systematic analyses of mosaic variants detected by diagnostic exome sequencing for diverse clinical indications have been performed.

Methods

To investigate the frequency, type, allelic fraction, and phenotypic consequences of clinically relevant somatic mosaic single nucleotide variants (SNVs) and characteristics of the corresponding genes, we retrospectively queried reported mosaic variants from a cohort of ~ 12,000 samples submitted for clinical exome sequencing (ES) at Baylor Genetics.

Results

We found 120 mosaic variants involving 107 genes, including 80 mosaic SNVs in proband samples and 40 in parental/grandparental samples. Average mosaic alternate allele fraction (AAF) detected in autosomes and in X-linked disease genes in females was 18.2% compared with 34.8% in X-linked disease genes in males. Of these mosaic variants, 74 variants (61.7%) were classified as pathogenic or likely pathogenic and 46 (38.3%) as variants of uncertain significance. Mosaic variants occurred in disease genes associated with autosomal dominant (AD) or AD/autosomal recessive (AR) (67/120, 55.8%), X-linked (33/120, 27.5%), AD/somatic (10/120, 8.3%), and AR (8/120, 6.7%) inheritance. Of note, 1.7% (2/120) of variants were found in genes in which only somatic events have been described. Nine genes had recurrent mosaic events in unrelated individuals which accounted for 18.3% (22/120) of all detected mosaic variants in this study. The proband group was enriched for mosaicism affecting Ras signaling pathway genes.

Conclusions

In sum, an estimated 1.5% of all molecular diagnoses made in this cohort could be attributed to a mosaic variant detected in the proband, while parental mosaicism was identified in 0.3% of families analyzed. As ES design favors breadth over depth of coverage, this estimate of the prevalence of mosaic variants likely represents an underestimate of the total number of clinically relevant mosaic variants in our cohort.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0658-2) contains supplementary material, which is available to authorized users.

Comparison of kNN and k-means optimization methods of reference set selection for improved CNV callers performance

BACKGROUND

There are over 25 tools dedicated for the detection of Copy Number Variants (CNVs) using Whole Exome Sequencing (WES) data based on read depth analysis. The tools reported consist of several steps, including: (i) calculation of read depth for each sequencing target, (ii) normalization, (iii) segmentation and (iv) actual CNV calling. The essential aspect of the entire process is the normalization stage, in which systematic errors and biases are removed and the reference sample set is used to increase the signal-to-noise ratio. Although some CNV calling tools use dedicated algorithms to obtain the optimal reference sample set, most of the advanced CNV callers do not include this feature. To our knowledge, this work is the first attempt to assess the impact of reference sample set selection on CNV detection performance.

METHODS

We used WES data from the 1000 Genomes project to evaluate the impact of various methods of reference sample set selection on CNV calling performance of three chosen state-of-the-art tools: CODEX, CNVkit and exomeCopy. Two naive solutions (all samples as reference set and random selection) as well as two clustering methods (k-means and k nearest neighbours (kNN) with a variable number of clusters or group sizes) have been evaluated to discover the best performing sample selection method.

RESULTS AND CONCLUSIONS

The performed experiments have shown that the appropriate selection of the reference sample set may greatly improve the CNV detection rate. In particular, we found that smart reduction of reference sample size may significantly increase the algorithms' precision while having negligible negative effect on sensitivity. We observed that a complete CNV calling process with the k-means algorithm as the selection method has significantly better time complexity than kNN-based solution.

The role of FREM2 and FRAS1 in the development of congenital diaphragmatic hernia

Congenital diaphragmatic hernia (CDH) has been reported twice in individuals with a clinical diagnosis of Fraser syndrome, a genetic disorder that can be caused by recessive mutations affecting FREM2 and FRAS1. In the extracellular matrix, FREM2 and FRAS1 form a self-stabilizing complex with FREM1, a protein whose deficiency causes sac CDH in humans and mice. By sequencing FREM2 and FRAS1 in a CDH cohort, and searching online databases, we identified five individuals who carried recessive or double heterozygous, putatively deleterious variants in these genes which may represent susceptibility alleles. Three of these alleles were significantly enriched in our CDH cohort compared with ethnically matched controls. We subsequently demonstrated that 8% of Frem2ne/ne and 1% of Fras1Q1263*/Q1263* mice develop the same type of anterior sac CDH seen in FREM1-deficient mice. We went on to show that development of sac hernias in FREM1-deficient mice is preceded by failure of anterior mesothelial fold progression resulting in the persistence of an amuscular, poorly vascularized anterior diaphragm that is abnormally adherent to the underlying liver. Herniation occurs in the perinatal period when the expanding liver protrudes through this amuscular region of the anterior diaphragm that is juxtaposed to areas of muscular diaphragm. Based on these data, we conclude that deficiency of FREM2, and possibly FRAS1, are associated with an increased risk of developing CDH and that loss of the FREM1/FREM2/FRAS1 complex, or its function, leads to anterior sac CDH development through its effects on mesothelial fold progression.

Complex Compound Inheritance of Lethal Lung Developmental Disorders Due to Disruption of the TBX-FGF Pathway

Primary defects in lung branching morphogenesis, resulting in neonatal lethal pulmonary hypoplasias, are incompletely understood. To elucidate the pathogenetics of human lung development, we studied a unique collection of samples obtained from deceased individuals with clinically and histopathologically diagnosed interstitial neonatal lung disorders: acinar dysplasia (n = 14), congenital alveolar dysplasia (n = 2), and other lethal lung hypoplasias (n = 10). We identified rare heterozygous copy-number variant deletions or single-nucleotide variants (SNVs) involving TBX4 (n = 8 and n = 2, respectively) or FGF10 (n = 2 and n = 2, respectively) in 16/26 (61%) individuals. In addition to TBX4, the overlapping ∼2 Mb recurrent and nonrecurrent deletions at 17q23.1q23.2 identified in seven individuals with lung hypoplasia also remove a lung-specific enhancer region. Individuals with coding variants involving either TBX4 or FGF10 also harbored at least one non-coding SNV in the predicted lung-specific enhancer region, which was absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies. The occurrence of rare coding variants involving TBX4 or FGF10 with the putative hypomorphic non-coding SNVs implies a complex compound inheritance of these pulmonary hypoplasias. Moreover, they support the importance of TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling in human lung organogenesis and help to explain the histopathological continuum observed in these rare lethal developmental disorders of the lung.

Identifying Genes Whose Mutant Transcripts Cause Dominant Disease Traits by Potential Gain-of-Function Alleles

Premature termination codon (PTC)-bearing transcripts are often degraded by nonsense-mediated decay (NMD) resulting in loss-of-function (LoF) alleles. However, not all PTCs result in LoF mutations, i.e., some such transcripts escape NMD and are translated to truncated peptide products that result in disease due to gain-of-function (GoF) effects. Since the location of the PTC is a major factor determining transcript fate, we hypothesized that depletion of protein-truncating variants (PTVs) within the gene region predicted to escape NMD in control databases could provide a rank for genic susceptibility for disease through GoF versus LoF. We developed an NMD escape intolerance score to rank genes based on the depletion of PTVs that would render them able to escape NMD using the Atherosclerosis Risk in Communities Study (ARIC) and the Exome Aggregation Consortium (ExAC) control databases, which was further used to screen the Baylor-Center for Mendelian Genomics disease database. This analysis revealed 1,996 genes significantly depleted for PTVs that are predicted to escape from NMD, i.e., PTVesc; further studies provided evidence that revealed a subset as candidate genes underlying Mendelian phenotypes. Importantly, these genes have characteristically low pLI scores, which can cause them to be overlooked as candidates for dominant diseases. Collectively, we demonstrate that this NMD escape intolerance score is an effective and efficient tool for gene discovery in Mendelian diseases due to production of truncated or altered proteins. More importantly, we provide a complementary analytical tool to aid identification of genes associated with dominant traits through a mechanism distinct from LoF.

Comprehensive genomic analysis of patients with disorders of cerebral cortical development

Malformations of cortical development (MCDs) manifest with structural brain anomalies that lead to neurologic sequelae, including epilepsy, cerebral palsy, developmental delay, and intellectual disability. To investigate the underlying genetic architecture of patients with disorders of cerebral cortical development, a cohort of 54 patients demonstrating neuroradiologic signs of MCDs was investigated. Individual genomes were interrogated for single-nucleotide variants (SNV) and copy number variants (CNV) with whole-exome sequencing and chromosomal microarray studies. Variation affecting known MCDs-associated genes was found in 16/54 cases, including 11 patients with SNV, 2 patients with CNV, and 3 patients with both CNV and SNV, at distinct loci. Diagnostic pathogenic SNV and potentially damaging variants of unknown significance (VUS) were identified in two groups of seven individuals each. We demonstrated that de novo variants are important among patients with MCDs as they were identified in 10/16 individuals with a molecular diagnosis. Three patients showed changes in known MCDs genes  and a clinical phenotype beyond the usual characteristics observed, i.e., phenotypic expansion, for a particular known disease gene clinical entity. We also discovered 2 likely candidate genes, CDH4, and ASTN1, with human and animal studies supporting their roles in brain development, and 5 potential candidate genes. Our findings emphasize genetic heterogeneity of MCDs disorders and postulate potential novel candidate genes involved in cerebral cortical development.

Benchmarking distributed data warehouse solutions for storing genomic variant information

Database URL

https://github.com/ZSI-Bio/variantsdwh.

Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor–Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17–50% of pathogenic CNVs in different disease cohorts where 7.1–11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

Identification of novel candidate disease genes from de novo exonic copy number variants

BACKGROUND

Exon-targeted microarrays can detect small (<1000 bp) intragenic copy number variants (CNVs), including those that affect only a single exon. This genome-wide high-sensitivity approach increases the molecular diagnosis for conditions with known disease-associated genes, enables better genotype-phenotype correlations, and facilitates variant allele detection allowing novel disease gene discovery.

METHODS

We retrospectively analyzed data from 63,127 patients referred for clinical chromosomal microarray analysis (CMA) at Baylor Genetics laboratories, including 46,755 individuals tested using exon-targeted arrays, from 2007 to 2017. Small CNVs harboring a single gene or two to five non-disease-associated genes were identified; the genes involved were evaluated for a potential disease association.

RESULTS

In this clinical population, among rare CNVs involving any single gene reported in 7200 patients (11%), we identified 145 de novo autosomal CNVs (117 losses and 28 intragenic gains), 257 X-linked deletion CNVs in males, and 1049 inherited autosomal CNVs (878 losses and 171 intragenic gains); 111 known disease genes were potentially disrupted by de novo autosomal or X-linked (in males) single-gene CNVs. Ninety-one genes, either recently proposed as candidate disease genes or not yet associated with diseases, were disrupted by 147 single-gene CNVs, including 37 de novo deletions and ten de novo intragenic duplications on autosomes and 100 X-linked CNVs in males. Clinical features in individuals with de novo or X-linked CNVs encompassing at most five genes (224 bp to 1.6 Mb in size) were compared to those in individuals with larger-sized deletions (up to 5 Mb in size) in the internal CMA database or loss-of-function single nucleotide variants (SNVs) detected by clinical or research whole-exome sequencing (WES). This enabled the identification of recently published genes (BPTF, NONO, PSMD12, TANGO2, and TRIP12), novel candidate disease genes (ARGLU1 and STK3), and further confirmation of disease association for two recently proposed disease genes (MEIS2 and PTCHD1). Notably, exon-targeted CMA detected several pathogenic single-exon CNVs missed by clinical WES analyses.

CONCLUSIONS

Together, these data document the efficacy of exon-targeted CMA for detection of genic and exonic CNVs, complementing and extending WES in clinical diagnostics, and the potential for discovery of novel disease genes by genome-wide assay.

An Organismal CNV Mutator Phenotype Restricted to Early Human Development

De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the "CNV-mutator state," and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology.

Iterative Sequencing and Variant Screening (ISVS) as a novel pathogenic mutations search strategy - application for TMPRSS3 mutations screen

Autosomal recessive diseases (ARD) are typically caused by a limited number of mutations whose identification is challenged by their low prevalence. Our purpose was to develop a novel approach allowing an efficient search for mutations causing ARD and evaluation of their pathogenicity without a control group. We developed Iterative Sequencing and Variant Screening (ISVS) approach based on iterative cycles of gene sequencing and mutation screening, and ISVS Simulator software ( http://zsibio.ii.pw.edu.pl/shiny/isvs/ ) for assessment of detected variants' significance. As shown by simulations, ISVS efficiently identifies and correctly classifies pathogenic mutations except for cases where the gene of interest has extremely high number of low frequency nonpathogenic variants. By applying ISVS, we found 4 known and 9 novel (p.C73Y, p.S124L, p.C194Mfs*17, c.782 + 2 T > A, c.953-5 A > G, p.L325Q, p.D334Mfs*24, p.R436G, p.M448T) TMPRSS3 variants among deaf patients. For 3 known and 5 novel variants the disease association was supported by ISVS Simulator odds >90:1. Pathogenicity of 6 novel mutations has been supported by in-silico predictions of variants' deleteriousness. By directly comparing variant prevalence in patients and controls, disease association was demonstrated only for two variants and it was relatively weak (P < 0.05). Summarizing, ISVS strategy and ISVS Simulator are useful for detection of genetic variants causing AR diseases.

Lessons learned from additional research analyses of unsolved clinical exome cases

BACKGROUND

Given the rarity of most single-gene Mendelian disorders, concerted efforts of data exchange between clinical and scientific communities are critical to optimize molecular diagnosis and novel disease gene discovery.

METHODS

We designed and implemented protocols for the study of cases for which a plausible molecular diagnosis was not achieved in a clinical genomics diagnostic laboratory (i.e. unsolved clinical exomes). Such cases were recruited to a research laboratory for further analyses, in order to potentially: (1) accelerate novel disease gene discovery; (2) increase the molecular diagnostic yield of whole exome sequencing (WES); and (3) gain insight into the genetic mechanisms of disease. Pilot project data included 74 families, consisting mostly of parent-offspring trios. Analyses performed on a research basis employed both WES from additional family members and complementary bioinformatics approaches and protocols.

RESULTS

Analysis of all possible modes of Mendelian inheritance, focusing on both single nucleotide variants (SNV) and copy number variant (CNV) alleles, yielded a likely contributory variant in 36% (27/74) of cases. If one includes candidate genes with variants identified within a single family, a potential contributory variant was identified in a total of ~51% (38/74) of cases enrolled in this pilot study. The molecular diagnosis was achieved in 30/63 trios (47.6%). Besides this, the analysis workflow yielded evidence for pathogenic variants in disease-associated genes in 4/6 singleton cases (66.6%), 1/1 multiplex family involving three affected siblings, and 3/4 (75%) quartet families. Both the analytical pipeline and the collaborative efforts between the diagnostic and research laboratories provided insights that allowed recent disease gene discoveries (PURA, TANGO2, EMC1, GNB5, ATAD3A, and MIPEP) and increased the number of novel genes, defined in this study as genes identified in more than one family (DHX30 and EBF3).

CONCLUSION

An efficient genomics pipeline in which clinical sequencing in a diagnostic laboratory is followed by the detailed reanalysis of unsolved cases in a research environment, supplemented with WES data from additional family members, and subject to adjuvant bioinformatics analyses including relaxed variant filtering parameters in informatics pipelines, can enhance the molecular diagnostic yield and provide mechanistic insights into Mendelian disorders. Implementing these approaches requires collaborative clinical molecular diagnostic and research efforts.

Haploinsufficiency of the E3 ubiquitin-protein ligase gene TRIP12 causes intellectual disability with or without autism spectrum disorders, speech delay, and dysmorphic features

Impairment of ubiquitin-proteasome system activity involving ubiquitin ligase genes UBE3A, UBE3B, and HUWE1 and deubiquitinating enzyme genes USP7 and USP9X has been reported in patients with neurodevelopmental delays. To date, only a handful of single-nucleotide variants (SNVs) and copy-number variants (CNVs) involving TRIP12, encoding a member of the HECT domain E3 ubiquitin ligases family on chromosome 2q36.3 have been reported. Using chromosomal microarray analysis and whole-exome sequencing (WES), we have identified, respectively, five deletion CNVs and four inactivating SNVs (two frameshifts, one missense, and one splicing) in TRIP12. Seven of these variants were found to be de novo; parental studies could not be completed in two families. Quantitative PCR analyses of the splicing mutation showed a dramatically decreased level of TRIP12 mRNA in the proband compared to the family controls, indicating a loss-of-function mechanism. The shared clinical features include intellectual disability with or without autistic spectrum disorders, speech delay, and facial dysmorphism. Our findings demonstrate that E3 ubiquitin ligase TRIP12 plays an important role in nervous system development and function. The nine presented pathogenic variants further document that TRIP12 haploinsufficiency causes a childhood-onset neurodevelopmental disorder. Finally, our data enable expansion of the phenotypic spectrum of ubiquitin-proteasome dependent disorders.

De Novo Disruption of the Proteasome Regulatory Subunit PSMD12 Causes a Syndromic Neurodevelopmental Disorder

Degradation of proteins by the ubiquitin-proteasome system (UPS) is an essential biological process in the development of eukaryotic organisms. Dysregulation of this mechanism leads to numerous human neurodegenerative or neurodevelopmental disorders. Through a multi-center collaboration, we identified six de novo genomic deletions and four de novo point mutations involving PSMD12, encoding the non-ATPase subunit PSMD12 (aka RPN5) of the 19S regulator of 26S proteasome complex, in unrelated individuals with intellectual disability, congenital malformations, ophthalmologic anomalies, feeding difficulties, deafness, and subtle dysmorphic facial features. We observed reduced PSMD12 levels and an accumulation of ubiquitinated proteins without any impairment of proteasome catalytic activity. Our PSMD12 loss-of-function zebrafish CRISPR/Cas9 model exhibited microcephaly, decreased convolution of the renal tubules, and abnormal craniofacial morphology. Our data support the biological importance of PSMD12 as a scaffolding subunit in proteasome function during development and neurogenesis in particular; they enable the definition of a neurodevelopmental disorder due to PSMD12 variants, expanding the phenotypic spectrum of UPS-dependent disorders.

Primary immunodeficiency diseases: Genomic approaches delineate heterogeneous Mendelian disorders

BACKGROUND

Primary immunodeficiency diseases (PIDDs) are clinically and genetically heterogeneous disorders thus far associated with mutations in more than 300 genes. The clinical phenotypes derived from distinct genotypes can overlap. Genetic etiology can be a prognostic indicator of disease severity and can influence treatment decisions.

OBJECTIVE

We sought to investigate the ability of whole-exome screening methods to detect disease-causing variants in patients with PIDDs.

METHODS

Patients with PIDDs from 278 families from 22 countries were investigated by using whole-exome sequencing. Computational copy number variant (CNV) prediction pipelines and an exome-tiling chromosomal microarray were also applied to identify intragenic CNVs. Analytic approaches initially focused on 475 known or candidate PIDD genes but were nonexclusive and further tailored based on clinical data, family history, and immunophenotyping.

RESULTS

A likely molecular diagnosis was achieved in 110 (40%) unrelated probands. Clinical diagnosis was revised in about half (60/110) and management was directly altered in nearly a quarter (26/110) of families based on molecular findings. Twelve PIDD-causing CNVs were detected, including 7 smaller than 30 Kb that would not have been detected with conventional diagnostic CNV arrays.

CONCLUSION

This high-throughput genomic approach enabled detection of disease-related variants in unexpected genes; permitted detection of low-grade constitutional, somatic, and revertant mosaicism; and provided evidence of a mutational burden in mixed PIDD immunophenotypes.