Reactive Neutrophil Responses Dependent on the Receptor Tyrosine Kinase c-MET Limit Cancer Immunotherapy.

Inhibitors of the receptor tyrosine kinase c-MET are currently used in the clinic to target oncogenic signaling in tumor cells. We found that concomitant c-MET inhibition promoted adoptive T cell transfer and checkpoint immunotherapies in murine cancer models by increasing effector T cell infiltration in tumors. This therapeutic effect was independent of tumor cell-intrinsic c-MET dependence. Mechanistically, c-MET inhibition impaired the reactive mobilization and recruitment of neutrophils into tumors and draining lymph nodes in response to cytotoxic immunotherapies. In the absence of c-MET inhibition, neutrophils recruited to T cell-inflamed microenvironments rapidly acquired immunosuppressive properties, restraining T cell expansion and effector functions. In cancer patients, high serum levels of the c-MET ligand HGF correlated with increasing neutrophil counts and poor responses to checkpoint blockade therapies. Our findings reveal a role for the HGF/c-MET pathway in neutrophil recruitment and function and suggest that c-MET inhibitor co-treatment may improve responses to cancer immunotherapy in settings beyond c-MET-dependent tumors.

TMPRSS2 isoform 1 downregulation by G-quadruplex stabilization induces SARS-CoV-2 replication arrest.

BACKGROUND: SARS-CoV-2 infection depends on the host cell factors angiotensin-converting enzyme 2, ACE2, and the transmembrane serinprotease 2, TMPRSS2. Potential inhibitors of these proteins would be ideal targets against severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) infection. Our data opens the possibility that changes within TMPRSS2 can modulate the outcome during a SARS-CoV-2 infection. RESULTS: We reveal that TMPRSS2 acts not only during viral entry but has also an important role during viral replication. In addition to previous functions for TMPRSS2 during viral entry, we determined by specific downregulation of distinct isoforms that only isoform 1 controls and supports viral replication. G-quadruplex (G4) stabilization by chemical compounds impacts TMPRSS2 gene expression. Here we extend and in-depth characterize these observations and identify that a specific G4 in the first exon of the TMPRSS2 isoform 1 is particular targeted by the G4 ligand and affects viral replication. Analysis of potential single nucleotide polymorphisms (SNPs) reveals that a reported SNP at this G4 in isoform 1 destroys the G4 motif and makes TMPRSS2 ineffective towards G4 treatment. CONCLUSION: These findings uncover a novel mechanism in which G4 stabilization impacts SARS-CoV-2 replication by changing TMPRSS2 isoform 1 gene expression.

Predicting the pathogenicity of missense variants using features derived from AlphaFold2.

MOTIVATION: Missense variants are a frequent class of variation within the coding genome, and some of them cause Mendelian diseases. Despite advances in computational prediction, classifying missense variants into pathogenic or benign remains a major challenge in the context of personalized medicine. Recently, the structure of the human proteome was derived with unprecedented accuracy using the artificial intelligence system AlphaFold2. This raises the question of whether AlphaFold2 wild-type structures can improve the accuracy of computational pathogenicity prediction for missense variants. RESULTS: To address this, we first engineered a set of features for each amino acid from these structures. We then trained a random forest to distinguish between relatively common (proxy-benign) and singleton (proxy-pathogenic) missense variants from gnomAD v3.1. This yielded a novel AlphaFold2-based pathogenicity prediction score, termed AlphScore. Important feature classes used by AlphScore are solvent accessibility, amino acid network related features, features describing the physicochemical environment, and AlphaFold2's quality parameter (predicted local distance difference test). AlphScore alone showed lower performance than existing in silico scores used for missense prediction, such as CADD or REVEL. However, when AlphScore was added to those scores, the performance increased, as measured by the approximation of deep mutational scan data, as well as the prediction of expert-curated missense variants from the ClinVar database. Overall, our data indicate that the integration of AlphaFold2-predicted structures can improve pathogenicity prediction of missense variants. AVAILABILITY AND IMPLEMENTATION: AlphScore, combinations of AlphScore with existing scores, as well as variants used for training and testing are publicly available.

Msx1 deficiency interacts with hypoxia and induces a morphogenetic regulation during mouse lip development.

Nonsyndromic clefts of the lip and palate are common birth defects resulting from gene-gene and gene-environment interactions. Mutations in human MSX1 have been linked to orofacial clefting and we show here that Msx1 deficiency causes a growth defect of the medial nasal process (Mnp) in mouse embryos. Although this defect alone does not disrupt lip formation, Msx1-deficient embryos develop a cleft lip when the mother is transiently exposed to reduced oxygen levels or to phenytoin, a drug known to cause embryonic hypoxia. In the absence of interacting environmental factors, the Mnp growth defect caused by Msx1 deficiency is modified by a Pax9-dependent 'morphogenetic regulation', which modulates Mnp shape, rescues lip formation and involves a localized abrogation of Bmp4-mediated repression of Pax9 Analyses of GWAS data revealed a genome-wide significant association of a Gene Ontology morphogenesis term (including assigned roles for MSX1, MSX2, PAX9, BMP4 and GREM1) specifically for nonsyndromic cleft lip with cleft palate. Our data indicate that MSX1 mutations could increase the risk for cleft lip formation by interacting with an impaired morphogenetic regulation that adjusts Mnp shape, or through interactions that inhibit Mnp growth.

Common, low-frequency, rare, and ultra-rare coding variants contribute to COVID-19 severity.

The combined impact of common and rare exonic variants in COVID-19 host genetics is currently insufficiently understood. Here, common and rare variants from whole-exome sequencing data of about 4000 SARS-CoV-2-positive individuals were used to define an interpretable machine-learning model for predicting COVID-19 severity. First, variants were converted into separate sets of Boolean features, depending on the absence or the presence of variants in each gene. An ensemble of LASSO logistic regression models was used to identify the most informative Boolean features with respect to the genetic bases of severity. The Boolean features selected by these logistic models were combined into an Integrated PolyGenic Score that offers a synthetic and interpretable index for describing the contribution of host genetics in COVID-19 severity, as demonstrated through testing in several independent cohorts. Selected features belong to ultra-rare, rare, low-frequency, and common variants, including those in linkage disequilibrium with known GWAS loci. Noteworthily, around one quarter of the selected genes are sex-specific. Pathway analysis of the selected genes associated with COVID-19 severity reflected the multi-organ nature of the disease. The proposed model might provide useful information for developing diagnostics and therapeutics, while also being able to guide bedside disease management.

Genome-wide association study reveals new insights into the heritability and genetic correlates of developmental dyslexia.

Developmental dyslexia (DD) is a learning disorder affecting the ability to read, with a heritability of 40-60%. A notable part of this heritability remains unexplained, and large genetic studies are warranted to identify new susceptibility genes and clarify the genetic bases of dyslexia. We carried out a genome-wide association study (GWAS) on 2274 dyslexia cases and 6272 controls, testing associations at the single variant, gene, and pathway level, and estimating heritability using single-nucleotide polymorphism (SNP) data. We also calculated polygenic scores (PGSs) based on large-scale GWAS data for different neuropsychiatric disorders and cortical brain measures, educational attainment, and fluid intelligence, testing them for association with dyslexia status in our sample. We observed statistically significant (p < 2.8 × 10-6) enrichment of associations at the gene level, for LOC388780 (20p13; uncharacterized gene), and for VEPH1 (3q25), a gene implicated in brain development. We estimated an SNP-based heritability of 20-25% for DD, and observed significant associations of dyslexia risk with PGSs for attention deficit hyperactivity disorder (at pT = 0.05 in the training GWAS: OR = 1.23[1.16; 1.30] per standard deviation increase; p = 8 × 10-13), bipolar disorder (1.53[1.44; 1.63]; p = 1 × 10-43), schizophrenia (1.36[1.28; 1.45]; p = 4 × 10-22), psychiatric cross-disorder susceptibility (1.23[1.16; 1.30]; p = 3 × 10-12), cortical thickness of the transverse temporal gyrus (0.90[0.86; 0.96]; p = 5 × 10-4), educational attainment (0.86[0.82; 0.91]; p = 2 × 10-7), and intelligence (0.72[0.68; 0.76]; p = 9 × 10-29). This study suggests an important contribution of common genetic variants to dyslexia risk, and novel genomic overlaps with psychiatric conditions like bipolar disorder, schizophrenia, and cross-disorder susceptibility. Moreover, it revealed the presence of shared genetic foundations with a neural correlate previously implicated in dyslexia by neuroimaging evidence.

Systematic assays and resources for the functional annotation of non-coding variants.

Identification of genetic variation in individual genomes is now a routine procedure in human genetic research and diagnostics. For many variants, however, insufficient evidence is available to establish a pathogenic effect, particularly for variants in non-coding regions. Furthermore, the sheer number of candidate variants renders testing in individual assays virtually impossible. While scalable approaches are being developed, the selection of methods and resources, and the application of a given framework to a particular disease or trait remain major challenges. This limits the translation of results from both genome-wide association studies and genome sequencing. Here, we discuss computational and experimental approaches available for functional annotation of non-coding variation.

Extending the allelic spectrum at noncoding risk loci of orofacial clefting.

Genome-wide association studies (GWAS) have generated unprecedented insights into the genetic etiology of orofacial clefting (OFC). The moderate effect sizes of associated noncoding risk variants and limited access to disease-relevant tissue represent considerable challenges for biological interpretation of genetic findings. As rare variants with stronger effect sizes are likely to also contribute to OFC, an alternative approach to delineate pathogenic mechanisms is to identify private mutations and/or an increased burden of rare variants in associated regions. This report describes a framework for targeted resequencing at selected noncoding risk loci contributing to nonsyndromic cleft lip with/without cleft palate (nsCL/P), the most frequent OFC subtype. Based on GWAS data, we selected three risk loci and identified candidate regulatory regions (CRRs) through the integration of credible SNP information, epigenetic data from relevant cells/tissues, and conservation scores. The CRRs (total 57 kb) were resequenced in a multiethnic study population (1061 patients; 1591 controls), using single-molecule molecular inversion probe technology. Combining evidence from in silico variant annotation, pedigree- and burden analyses, we identified 16 likely deleterious rare variants that represent new candidates for functional studies in nsCL/P. Our framework is scalable and represents a promising approach to the investigation of additional congenital malformations with multifactorial etiology.

Evaluating shared genetic influences on nonsyndromic cleft lip/palate and oropharyngeal neoplasms.

It has been hypothesised that nonsyndromic cleft lip/palate (nsCL/P) and cancer may share aetiological risk factors. Population studies have found inconsistent evidence for increased incidence of cancer in nsCL/P cases, but several genes (e.g., CDH1, AXIN2) have been implicated in the aetiologies of both phenotypes. We aimed to evaluate shared genetic aetiology between nsCL/P and oral cavity/oropharyngeal cancers (OC/OPC), which affect similar anatomical regions. Using a primary sample of 5,048 OC/OPC cases and 5,450 controls of European ancestry and a replication sample of 750 cases and 336,319 controls from UK Biobank, we estimate genetic overlap using nsCL/P polygenic risk scores (PRS) with Mendelian randomization analyses performed to evaluate potential causal mechanisms. In the primary sample, we found strong evidence for an association between a nsCL/P PRS and increased odds of OC/OPC (per standard deviation increase in score, odds ratio [OR]: 1.09; 95% confidence interval [CI]: 1.04, 1.13; p = .000053). Although confidence intervals overlapped with the primary estimate, we did not find confirmatory evidence of an association between the PRS and OC/OPC in UK Biobank (OR 1.02; 95% CI: 0.95, 1.10; p = .55). Mendelian randomization analyses provided evidence that major nsCL/P risk variants are unlikely to influence OC/OPC. Our findings suggest possible shared genetic influences on nsCL/P and OC/OPC.

Investigating the shared genetics of non-syndromic cleft lip/palate and facial morphology.

There is increasing evidence that genetic risk variants for non-syndromic cleft lip/palate (nsCL/P) are also associated with normal-range variation in facial morphology. However, previous analyses are mostly limited to candidate SNPs and findings have not been consistently replicated. Here, we used polygenic risk scores (PRS) to test for genetic overlap between nsCL/P and seven biologically relevant facial phenotypes. Where evidence was found of genetic overlap, we used bidirectional Mendelian randomization (MR) to test the hypothesis that genetic liability to nsCL/P is causally related to implicated facial phenotypes. Across 5,804 individuals of European ancestry from two studies, we found strong evidence, using PRS, of genetic overlap between nsCL/P and philtrum width; a 1 S.D. increase in nsCL/P PRS was associated with a 0.10 mm decrease in philtrum width (95% C.I. 0.054, 0.146; P = 2x10-5). Follow-up MR analyses supported a causal relationship; genetic variants for nsCL/P homogeneously cause decreased philtrum width. In addition to the primary analysis, we also identified two novel risk loci for philtrum width at 5q22.2 and 7p15.2 in our Genome-wide Association Study (GWAS) of 6,136 individuals. Our results support a liability threshold model of inheritance for nsCL/P, related to abnormalities in development of the philtrum.

Nine newly identified individuals refine the phenotype associated with MYT1L mutations.

Both point mutations and deletions of the MYT1L gene as well as microdeletions of chromosome band 2p25.3 including MYT1L are associated with intellectual disability, obesity, and behavioral problems. Thus, MYT1L is assumed to be the-at least mainly-causative gene in the 2p25.3 deletion syndrome. Here, we present comprehensive descriptions of nine novel individuals bearing MYT1L mutations; most of them single nucleotide variants (SNVs). This increases the number of known individuals with causative deletions or SNVs of MYT1L to 51. Since eight of the nine novel patients bear mutations affecting MYT1L only, the total number of such individuals now nearly equals the number of individuals with larger microdeletions affecting additional genes, allowing for a comprehensive phenotypic comparison of these two patient groups. For example, 55% of the individuals with mutations affecting MYT1L only were overweight or obese as compared to 86% of the individuals with larger microdeletions. A similar trend was observed regarding short stature with 5 versus 35%, respectively. However, these differences were nominally significant only after correction for multiple testing, further supporting the hypothesis that MYT1L haploinsufficiency is central to the 2p25.3 deletion phenotype. Most importantly, the large number of individuals with MYT1L mutations presented and reviewed here allowed for the delineation of a more comprehensive clinical picture. Seizures, postnatal short stature, macrocephaly, and microcephaly could be shown to be over-represented among individuals with MYT1L mutations.

Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.

Sequencing the GRHL3 Coding Region Reveals Rare Truncating Mutations and a Common Susceptibility Variant for Nonsyndromic Cleft Palate.

Nonsyndromic cleft lip with/without cleft palate (nsCL/P) and nonsyndromic cleft palate only (nsCPO) are the most frequent subphenotypes of orofacial clefts. A common syndromic form of orofacial clefting is Van der Woude syndrome (VWS) where individuals have CL/P or CPO, often but not always associated with lower lip pits. Recently, ∼5% of VWS-affected individuals were identified with mutations in the grainy head-like 3 gene (GRHL3). To investigate GRHL3 in nonsyndromic clefting, we sequenced its coding region in 576 Europeans with nsCL/P and 96 with nsCPO. Most strikingly, nsCPO-affected individuals had a higher minor allele frequency for rs41268753 (0.099) than control subjects (0.049; p = 1.24 × 10(-2)). This association was replicated in nsCPO/control cohorts from Latvia, Yemen, and the UK (pcombined = 2.63 × 10(-5); ORallelic = 2.46 [95% CI 1.6-3.7]) and reached genome-wide significance in combination with imputed data from a GWAS in nsCPO triads (p = 2.73 × 10(-9)). Notably, rs41268753 is not associated with nsCL/P (p = 0.45). rs41268753 encodes the highly conserved p.Thr454Met (c.1361C>T) (GERP = 5.3), which prediction programs denote as deleterious, has a CADD score of 29.6, and increases protein binding capacity in silico. Sequencing also revealed four novel truncating GRHL3 mutations including two that were de novo in four families, where all nine individuals harboring mutations had nsCPO. This is important for genetic counseling: given that VWS is rare compared to nsCPO, our data suggest that dominant GRHL3 mutations are more likely to cause nonsyndromic than syndromic CPO. Thus, with rare dominant mutations and a common risk variant in the coding region, we have identified an important contribution for GRHL3 in nsCPO.

Meta-analysis Reveals Genome-Wide Significance at 15q13 for Nonsyndromic Clefting of Both the Lip and the Palate, and Functional Analyses Implicate GREM1 As a Plausible Causative Gene.

Nonsyndromic orofacial clefts are common birth defects with multifactorial etiology. The most common type is cleft lip, which occurs with or without cleft palate (nsCLP and nsCLO, respectively). Although genetic components play an important role in nsCLP, the genetic factors that predispose to palate involvement are largely unknown. In this study, we carried out a meta-analysis on genetic and clinical data from three large cohorts and identified strong association between a region on chromosome 15q13 and nsCLP (P = 8.13 × 10(-14) for rs1258763; relative risk (RR): 1.46, 95% confidence interval (CI): 1.32-1.61)) but not nsCLO (P = 0.27; RR: 1.09 (0.94-1.27)). The 5 kb region of strongest association maps downstream of Gremlin-1 (GREM1), which encodes a secreted antagonist of the BMP4 pathway. We show during mouse embryogenesis, Grem1 is expressed in the developing lip and soft palate but not in the hard palate. This is consistent with genotype-phenotype correlations between rs1258763 and a specific nsCLP subphenotype, since a more than two-fold increase in risk was observed in patients displaying clefts of both the lip and soft palate but who had an intact hard palate (RR: 3.76, CI: 1.47-9.61, Pdiff<0.05). While we did not find lip or palate defects in Grem1-deficient mice, wild type embryonic palatal shelves developed divergent shapes when cultured in the presence of ectopic Grem1 protein (P = 0.0014). The present study identified a non-coding region at 15q13 as the second, genome-wide significant locus specific for nsCLP, after 13q31. Moreover, our data suggest that the closely located GREM1 gene contributes to a rare clinical nsCLP entity. This entity specifically involves abnormalities of the lip and soft palate, which develop at different time-points and in separate anatomical regions.

Genome-wide association study and meta-analysis identify ISL1 as genome-wide significant susceptibility gene for bladder exstrophy.

The bladder exstrophy-epispadias complex (BEEC) represents the severe end of the uro-rectal malformation spectrum, and is thought to result from aberrant embryonic morphogenesis of the cloacal membrane and the urorectal septum. The most common form of BEEC is isolated classic bladder exstrophy (CBE). To identify susceptibility loci for CBE, we performed a genome-wide association study (GWAS) of 110 CBE patients and 1,177 controls of European origin. Here, an association was found with a region of approximately 220kb on chromosome 5q11.1. This region harbors the ISL1 (ISL LIM homeobox 1) gene. Multiple markers in this region showed evidence for association with CBE, including 84 markers with genome-wide significance. We then performed a meta-analysis using data from a previous GWAS by our group of 98 CBE patients and 526 controls of European origin. This meta-analysis also implicated the 5q11.1 locus in CBE risk. A total of 138 markers at this locus reached genome-wide significance in the meta-analysis, and the most significant marker (rs9291768) achieved a P value of 2.13 × 10-12. No other locus in the meta-analysis achieved genome-wide significance. We then performed murine expression analyses to follow up this finding. Here, Isl1 expression was detected in the genital region within the critical time frame for human CBE development. Genital regions with Isl1 expression included the peri-cloacal mesenchyme and the urorectal septum. The present study identified the first genome-wide significant locus for CBE at chromosomal region 5q11.1, and provides strong evidence for the hypothesis that ISL1 is the responsible candidate gene in this region.

Compound heterozygosity of low-frequency promoter deletions and rare loss-of-function mutations in TXNL4A causes Burn-McKeown syndrome.

Mutations in components of the major spliceosome have been described in disorders with craniofacial anomalies, e.g., Nager syndrome and mandibulofacial dysostosis type Guion-Almeida. The U5 spliceosomal complex of eight highly conserved proteins is critical for pre-mRNA splicing. We identified biallelic mutations in TXNL4A, a member of this complex, in individuals with Burn-McKeown syndrome (BMKS). This rare condition is characterized by bilateral choanal atresia, hearing loss, cleft lip and/or palate, and other craniofacial dysmorphisms. Mutations were found in 9 of 11 affected families. In 8 families, affected individuals carried a rare loss-of-function mutation (nonsense, frameshift, or microdeletion) on one allele and a low-frequency 34 bp deletion (allele frequency 0.76%) in the core promoter region on the other allele. In a single highly consanguineous family, formerly diagnosed as oculo-oto-facial dysplasia, the four affected individuals were homozygous for a 34 bp promoter deletion, which differed from the promoter deletion in the other families. Reporter gene and in vivo assays showed that the promoter deletions led to reduced expression of TXNL4A. Depletion of TXNL4A (Dib1) in yeast demonstrated reduced assembly of the tri-snRNP complex. Our results indicate that BMKS is an autosomal-recessive condition, which is frequently caused by compound heterozygosity of low-frequency promoter deletions in combination with very rare loss-of-function mutations.

Novel mutations in LRP6 highlight the role of WNT signaling in tooth agenesis.

PURPOSE: We aimed to identify a novel genetic cause of tooth agenesis (TA) and/or orofacial clefting (OFC) by combining whole-exome sequencing (WES) and targeted resequencing in a large cohort of TA and OFC patients. METHODS: WES was performed in two unrelated patients: one with severe TA and OFC and another with severe TA only. After deleterious mutations were identified in a gene encoding low-density lipoprotein receptor-related protein 6 (LRP6), all its exons were resequenced with molecular inversion probes in 67 patients with TA, 1,072 patients with OFC, and 706 controls. RESULTS: We identified a frameshift (c.4594delG, p.Cys1532fs) and a canonical splice-site mutation (c.3398-2A>C, p.?) in LRP6, respectively, in the patient with TA and OFC and in the patient with severe TA only. The targeted resequencing showed significant enrichment of unique LRP6 variants in TA patients but not in nonsyndromic OFC patients. Of the five variants in patients with TA, two affected the canonical splice site and three were missense variants; all variants segregated with the dominant phenotype, and in one case the missense mutation occurred de novo. CONCLUSION: Mutations in LRP6 cause TA in humans.Genet Med 18 11, 1158-1162.

Replication analysis of 15 susceptibility loci for nonsyndromic cleft lip with or without cleft palate in an italian population.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is one of the most common congenital malformations in humans. Its average global incidence is 1.7 per 1000 live births, with wide variation according to geographical location and ethnicity. Its etiology involves both genetic and environmental factors. The aim of the present study was to confirm genetic association of a selection of 15 candidate nsCL/P loci using an independent sample of the Italian population. METHODS: At least one single-nucleotide polymorphism (SNP) for each locus was genotyped in 380 nuclear trios. RESULTS: Transmission disequilibrium analysis revealed significant associations for three variants at two loci (8q24 and 1p22). Two SNPs at 8q24 showed the strongest level of association, the rs987525 (PTDT = 6.81 × 10(-6) ; homozygous relative risk = 3.60 [95% confidence interval, 2.12-6.13]), and the rs17241253 (PTDT = 1.03 × 10(-5) ; homozygous relative risk = 3.75 [95% confidence interval, 2.10-6.67]). Four additional loci (at 1q32, 3q12, 8q21, and 10q25) achieved nominally significant p-values. Two SNPs at 1p36 achieved p-values of < 0.1. The present data suggest that 6 of the 15 analyzed nsCL/P risk loci contribute significantly to nsCL/P risk in the Italian population. These include the 1p22 locus, which previous research has implicated predominantly in Asian populations. CONCLUSION: Different loci, including 8q24 and 1p22 have been found associated with nsCL/P in multiple populations. Further efforts are needed to identify causative variants and transfer knowledge to clinical application, such as personal genetic risk assessment.

Further evidence for deletions in 7p14.1 contributing to nonsyndromic cleft lip with or without cleft palate.

BACKGROUND: Nonsyndromic cleft with or without cleft palate (nsCL/P) is a common birth defect. Although genome-wide association studies (GWAS) have identified numerous risk variants, a considerable fraction of the genetic heritability remains unknown. The aim of the present study was to replicate a previous finding that de novo deletions in a 62 kb region of chromosome 7p14 are a risk factor for nsCL/P, using an independent cohort. METHODS: Data from a published case-control GWAS cohort of 399 patients and 1318 controls were used. Copy number variant (CNV) detection in the 62 kb candidate region of 7p14 was performed using QuantiSNP. Putative CNVs in probands were verified and validated by quantitative polymerase chain reaction. Segregation analyses were performed in family members for whom DNA was available. RESULTS: Within the 62 kb candidate region, a deletion of 7.4 kb showed association with nsCL/P (13/387 cases, 20/1300 controls, plowest = 0.024, odds ratio = 2.22). In all families with a sporadic case (n = 3), the deletion occurred de novo. In multiplex families, both incomplete segregation and incomplete penetrance were observed. CONCLUSION: The present data support the hypothesis that deletions at 7p14 are a common risk factor for nsCL/P. Genome-wide CNV analyses in nsCL/P cohorts are warranted to explore the functional relevance of these deletions and their contribution to nsCL/P, and to determine exact breakpoints. Birth Defects Research (Part A) 106:767-772, 2016. © 2016 Wiley Periodicals, Inc.