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.

A drosophila genetic resource of mutants to study mechanisms underlying human genetic diseases.

Invertebrate model systems are powerful tools for studying human disease owing to their genetic tractability and ease of screening. We conducted a mosaic genetic screen of lethal mutations on the Drosophila X chromosome to identify genes required for the development, function, and maintenance of the nervous system. We identified 165 genes, most of whose function has not been studied in vivo. In parallel, we investigated rare variant alleles in 1,929 human exomes from families with unsolved Mendelian disease. Genes that are essential in flies and have multiple human homologs were found to be likely to be associated with human diseases. Merging the human data sets with the fly genes allowed us to identify disease-associated mutations in six families and to provide insights into microcephaly associated with brain dysgenesis. This bidirectional synergism between fly genetics and human genomics facilitates the functional annotation of evolutionarily conserved genes involved in human health.

Human CLP1 mutations alter tRNA biogenesis, affecting both peripheral and central nervous system function.

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis.

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.

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/.

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.

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.

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.

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.

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.

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.

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.

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.

The Thousand Polish Genomes-A Database of Polish Variant Allele Frequencies.

Although Slavic populations account for over 4.5% of world inhabitants, no centralised, open-source reference database of genetic variation of any Slavic population exists to date. Such data are crucial for clinical genetics, biomedical research, as well as archeological and historical studies. The Polish population, which is homogenous and sedentary in its nature but influenced by many migrations of the past, is unique and could serve as a genetic reference for the Slavic nations. In this study, we analysed whole genomes of 1222 Poles to identify and genotype a wide spectrum of genomic variation, such as small and structural variants, runs of homozygosity, mitochondrial haplogroups, and de novo variants. Common variant analyses showed that the Polish cohort is highly homogenous and shares ancestry with other European populations. In rare variant analyses, we identified 32 autosomal-recessive genes with significantly different frequencies of pathogenic alleles in the Polish population as compared to the non-Finish Europeans, including C2, TGM5, NUP93, C19orf12, and PROP1. The allele frequencies for small and structural variants, calculated for 1076 unrelated individuals, are released publicly as The Thousand Polish Genomes database, and will contribute to the worldwide genomic resources available to researchers and clinicians.

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.

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.

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.

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.

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 Ca2+ 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 Ca2+ 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.

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.