Exome Survey and Candidate Gene Re-Sequencing Identifies Novel Exstrophy Candidate Genes and Implicates LZTR1 in Disease Formation.

BACKGROUND: The bladder exstrophy-epispadias complex (BEEC) is a spectrum of congenital abnormalities that involves the abdominal wall, the bony pelvis, the urinary tract, the external genitalia, and, in severe cases, the gastrointestinal tract as well. METHODS: Herein, we performed an exome analysis of case-parent trios with cloacal exstrophy (CE), the most severe form of the BEEC. Furthermore, we surveyed the exome of a sib-pair presenting with classic bladder exstrophy (CBE) and epispadias (E) only. Moreover, we performed large-scale re-sequencing of CBE individuals for novel candidate genes that were derived from the current exome analysis, as well as for previously reported candidate genes within the CBE phenocritical region, 22q11.2. RESULTS: The exome survey in the CE case-parent trios identified two candidate genes harboring de novo variants (NR1H2, GKAP1), four candidate genes with autosomal-recessive biallelic variants (AKR1B10, CLSTN3, NDST4, PLEKHB1) and one candidate gene with suggestive uniparental disomy (SVEP1). However, re-sequencing did not identify any additional variant carriers in these candidate genes. Analysis of the affected sib-pair revealed no candidate gene. Re-sequencing of the genes within the 22q11.2 CBE phenocritical region identified two highly conserved frameshift variants that led to early termination in two independent CBE males, in LZTR1 (c.978_985del, p.Ser327fster6) and in SLC7A4 (c.1087delC, p.Arg363fster68). CONCLUSIONS: According to previous studies, our study further implicates LZTR1 in CBE formation. Exome analysis-derived candidate genes from CE individuals may not represent a frequent indicator for other BEEC phenotypes and warrant molecular analysis before their involvement in disease formation can be assumed.

Identification of de novo variants in nonsyndromic cleft lip with/without cleft palate patients with low polygenic risk scores.

BACKGROUND: Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome-wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants. METHODS: To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent-trios with a low polygenic risk for the phenotype (discovery). We prioritized DNV-carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population-matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome-wide association data, expression, protein-protein-interactions), were used for final prioritization. CONCLUSION: In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re-sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.

Allele-specific transcription factor binding in a cellular model of orofacial clefting.

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a frequent congenital malformation with multifactorial etiology. While recent genome-wide association studies (GWAS) have identified several nsCL/P risk loci, the functional effects of the associated non-coding variants are largely unknown. Furthermore, additional risk loci remain undetected due to lack of power. As genetic variants might alter binding of transcription factors (TF), we here hypothesized that the integration of data from TF binding sites, expression analyses and nsCL/P GWAS might help to (i) identify functionally relevant variants at GWAS loci, and (ii) highlight novel risk variants that have been previously undetected. Analysing the craniofacial TF TFAP2A in human embryonic palatal mesenchyme (HEPM) cells, we identified 2845 TFAP2A ChIP-seq peaks, several of which were located near nsCL/P candidate genes (e.g. MSX1 and SPRY2). Comparison with independent data suggest that 802 of them might be specific to craniofacial development, and genes near these peaks are enriched in processes relevant to nsCL/P. Integration with nsCL/P GWAS data, however, did not show robust evidence for co-localization of common nsCL/P risk variants with TFAP2A ChIP-seq peaks. This data set represents a new resource for the analyses of craniofacial processes, and similar approaches with additional cell lines and TFs could be applied to generate further insights into nsCL/P etiology.

Re-sequencing of candidate genes FOXF1, HSPA6, HAAO, and KYNU in 522 individuals with VATER/VACTERL, VACTER/VACTERL-like association, and isolated anorectal malformation.

BACKGROUND: The acronym VATER/VACTERL association describes the combination of at least three component features (CFs): vertebral defects (V), anorectal malformations (ARM) (A), cardiac defects (C), tracheoesophageal fistula with or without esophageal atresia (TE), renal malformations (R), and limb defects (L). Individuals presenting two CFs have been termed VATER/VACTERL-like. Recently, FOXF1, HSPA6, HAAO, KYNU, TRAP1, and ZIC3 have been proposed as candidate genes for VATER/VACTERL, VATER/VACTERL-like, and ARM. Re-sequencing studies identified disease-causing variants in TRAP1 and ZIC3, the contribution of other genes was not independently investigated. One affected variant carrier in FOXF1 was previously identified. Here we re-sequenced FOXF1, HSPA6, HAAO, and KYNU in 522 affected individuals. METHODS: Using molecular inversion probe (MIP) technology, re-sequencing was performed in 63 individuals with VATER/VACTERL association, 313 with VATER/VACTERL-like association, and 146 with ARM. All individuals were of European ethnicity. Variant filtering considered variants with a minor allele frequency (MAF) ≤0.01 for putative recessive disease-genes HSPA6, HAAO, and KYNU. For the putative dominant disease-gene FOXF1 we considered variants with a MAF ≤0.0001. In silico prediction tools were used for further prioritization. RESULTS: Only two variants in FOXF1 in two independently affected individuals [c.443G>T, p.(Cys148Phe); c.850T>C, p.(Tyr284His)] passed our filter criteria. One individual presented with ARM, the second presented with TE and C comprising atrial and ventricular septal defects. Sanger sequencing confirmed both variants but also their inheritance from the healthy mother. CONCLUSION: Our analysis suggests that FOXF1, HSPA6, HAAO and KYNU do not play a major role in the formation of VACTER/VACTERL phenotypes or ARM.

Molecular Genetic Screening in Patients With ACE Inhibitor/Angiotensin Receptor Blocker-Induced Angioedema to Explore the Role of Hereditary Angioedema Genes.

Angioedema is a relatively rare but potentially life-threatening adverse reaction to angiotensin-converting enzyme inhibitors (ACEi) and angiotensin receptor blockers (ARBs). As with hereditary forms of angioedema (HAE), this adverse reaction is mediated by bradykinin. Research suggests that ACEi/ARB-induced angioedema has a multifactorial etiology. In addition, recent case reports suggest that some ACEi/ARB-induced angioedema patients may carry pathogenic HAE variants. The aim of the present study was to investigate the possible association between ACEi/ARB-induced angioedema and HAE genes via systematic molecular genetic screening in a large cohort of ACEi/ARB-induced angioedema cases. Targeted re-sequencing of five HAE-associated genes (SERPING1, F12, PLG, ANGPT1, and KNG1) was performed in 212 ACEi/ARB-induced angioedema patients recruited in Germany/Austria, Sweden, and Denmark, and in 352 controls from a German cohort. Among patients, none of the identified variants represented a known pathogenic variant for HAE. Moreover, no significant association with ACEi/ARB-induced angioedema was found for any of the identified common [minor allele frequency (MAF) >5%] or rare (MAF < 5%) variants. However, several non-significant trends suggestive of possible protective effects were observed. The lowest p-value for an individual variant was found in PLG (rs4252129, p.R523W, p = 0.057, p.adjust > 0.999, Fisher's exact test). Variant p.R523W was found exclusively in controls and has previously been associated with decreased levels of plasminogen, a precursor of plasmin which is part of a pathway directly involved in bradykinin production. In addition, rare, potentially functional variants (MAF < 5%, Phred-scaled combined annotation dependent depletion score >10) showed a nominally significant enrichment in controls both: 1) across all five genes; and 2) in the F12 gene alone. However, these results did not withstand correction for multiple testing. In conclusion, our results suggest that HAE-associated mutations are, at best, a rare cause of ACEi/ARB-induced angioedema. Furthermore, we were unable to identify a significant association between ACEi/ARB-induced angioedema and other variants in the investigated genes. Further studies with larger sample sizes are warranted to draw more definite conclusions concerning variants with limited effect sizes, including protective variants.

TBK1 and TNFRSF13B mutations and an autoinflammatory disease in a child with lethal COVID-19.

Among children, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections are typically mild. Here, we describe the case of a 3.5-year-old girl with an unusually severe presentation of coronavirus disease (COVID-19). The child had an autoinflammatory disorder of unknown etiology, which had been treated using prednisolone and methotrexate, and her parents were half cousins of Turkish descent. After 5 days of nonspecific viral infection symptoms, tonic-clonic seizures occurred followed by acute cardiac insufficiency, multi-organ insufficiency, and ultimate death. Trio exome sequencing identified a homozygous splice-variant in the gene TBK1, and a homozygous missense variant in the gene TNFRSF13B. Heterozygous deleterious variants in the TBK1 gene have been associated with severe COVID-19, and the variant in the TNFRSF13B gene has been associated with common variable immunodeficiency (CVID). We suggest that the identified variants, the autoinflammatory disorder and its treatment, or a combination of these factors probably predisposed to lethal COVID-19 in the present case.