Accumulation of rare coding variants in genes implicated in risk of human cleft lip with or without cleft palate

The heterogeneous genetic architectures of orofacial clefts

Orofacial clefts (OFCs) are common, affecting 1:1000 live births. OFCs occur across a phenotypic spectrum - including cleft lip (CL), cleft lip and palate (CLP), or cleft palate (CP) - and can be further subdivided based on laterality, severity, or specific structures affected. Herein we review what is known about the genetic architecture underlying each of these subtypes, considering both shared and subtype-specific risks. While there are more known genetic similarities between CL and CLP than CP, recent research supports both shared and subtype-specific genetic risk factors within and between phenotypic classifications of OFCs. Larger sample sizes and deeper phenotyping data will be of increasing importance for the discovery of novel genetic risk factors for OFCs and various subtypes going forward.

Identification of Novel Risk Variants of Non-Syndromic Cleft Palate by Targeted Gene Panel Sequencing

Non-syndromic cleft palate (ns-CP) has a genetically heterogeneous aetiology. Numerous studies have suggested a crucial role of rare coding variants in characterizing the unrevealed component of genetic variation in ns-CP called the "missing heritability". Therefore, this study aimed to detect low-frequency variants that are implicated in ns-CP aetiology in the Polish population. For this purpose, coding regions of 423 genes associated with orofacial cleft anomalies and/or involved with facial development were screened in 38 ns-CP patients using the next-generation sequencing technology. After multistage selection and prioritisation, eight novel and four known rare variants that may influence an individual's risk of ns-CP were identified. Among detected alternations, seven were located in novel candidate genes for ns-CP, including COL17A1 (c.2435-1G>A), DLG1 (c.1586G>C, p.Glu562Asp), NHS (c.568G>C, p.Val190Leu-de novo variant), NOTCH2 (c.1997A>G, p.Tyr666Cys), TBX18 (c.647A>T, p.His225Leu), VAX1 (c.400G>A, p.Ala134Thr) and WNT5B (c.716G>T, p.Arg239Leu). The remaining risk variants were identified within genes previously linked to ns-CP, confirming their contribution to this anomaly. This list included ARHGAP29 (c.1706G>A, p.Arg569Gln), FLNB (c.3605A>G, Tyr1202Cys), IRF6 (224A>G, p.Asp75Gly-de novo variant), LRP6 (c.481C>A, p.Pro161Thr) and TP63 (c.353A>T, p.Asn118Ile). In summary, this study provides further insights into the genetic components contributing to ns-CP aetiology and identifies novel susceptibility genes for this craniofacial anomaly.

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.

Identification of novel susceptibility genes for non-syndromic cleft lip with or without cleft palate using NGS-based multigene panel testing

For non-syndromic cleft lip with or without cleft palate (ns-CL/P), the proportion of heritability explained by the known risk loci is estimated to be about 30% and is captured mainly by common variants identified in genome-wide association studies. To contribute to the explanation of the "missing heritability" problem for orofacial clefts, a candidate gene approach was taken to investigate the potential role of rare and private variants in the ns-CL/P risk. Using the next-generation sequencing technology, the coding sequence of a set of 423 candidate genes was analysed in 135 patients from the Polish population. After stringent multistage filtering, 37 rare coding and splicing variants of 28 genes were identified. 35% of these genetic alternations that may play a role of genetic modifiers influencing an individual's risk were detected in genes not previously associated with the ns-CL/P susceptibility, including COL11A1, COL17A1, DLX1, EFTUD2, FGF4, FGF8, FLNB, JAG1, NOTCH2, SHH, WNT5A and WNT9A. Significant enrichment of rare alleles in ns-CL/P patients compared with controls was also demonstrated for ARHGAP29, CHD7, COL17A1, FGF12, GAD1 and SATB2. In addition, analysis of panoramic radiographs of patients with identified predisposing variants may support the hypothesis of a common genetic link between orofacial clefts and dental abnormalities. In conclusion, our study has confirmed that rare coding variants might contribute to the genetic architecture of ns-CL/P. Since only single predisposing variants were identified in novel cleft susceptibility genes, future research will be required to confirm and fully understand their role in the aetiology of ns-CL/P.

The genetic factors contributing to the risk of cleft lip-cleft palate and their clinical utility

Cleft lip and cleft palate (CL/P) are among the most common congenital malformations in neonates and have syndromic or nonsyndromic forms. Nonsyndromic forms of malformation are being reported to be associated with chromosomal DNA modification by teratogenic exposure and to complex genetic contributions of multiple genes. Syndromic forms are shown to be related to chromosomal aberrations or monogenic diseases. There is a growing body of data illustrating the association of several genes with risk of developing this malformation, including genetic defects in T-box transcription factor-22 (TBX22), interferon regulatory factor-6 (IRF6), and poliovirus receptor-like-1 (PVRL1), responsible for X-linked cleft palate, cleft lip/palate-ectodermal dysplasia syndrome, and Van der Woude and popliteal pterygium syndromes, respectively. Genetic variants in MTR, PCYT1A, ASS1, SLC 25A13, GSTM1, GSTT1, SUMO1 BHMT1, and BHMT2 are being reported to be linked with CL/P risk. The etiology of nonsyndromic CLP is still remained to be unknown, although mutations in candidate genes have been found. Here, we provide an overview about the potential variants to be associated with CL/P for identification of the relative risk of CLP with respect to the basis of genetic background and environmental factors (e.g., dietary factors, alcohol use).

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

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.

Large-scale sequencing of flatfish genomes provides insights into the polyphyletic origin of their specialized body plan

The evolutionary and genetic origins of the specialized body plan of flatfish are largely unclear. We analyzed the genomes of 11 flatfish species representing 9 of the 14 Pleuronectiforme families and conclude that Pleuronectoidei and Psettodoidei do not form a monophyletic group, suggesting independent origins from different percoid ancestors. Genomic and transcriptomic data indicate that genes related to WNT and retinoic acid pathways, hampered musculature and reduced lipids might have functioned in the evolution of the specialized body plan of Pleuronectoidei. Evolution of Psettodoidei involved similar but not identical genes. Our work provides valuable resources and insights for understanding the genetic origins of the unusual body plan of flatfishes.

Genetic variants in S-adenosyl-methionine synthesis pathway and nonsyndromic cleft lip with or without cleft palate in Chile

BACKGROUND

The S-adenosyl-methionine (SAM) availability is crucial for DNA methylation, an epigenetic mechanism involved in nonsyndromic cleft lip with or without cleft palate (NSCL/P) expression. The aim of this study was to assess the association between single-nucleotide polymorphisms (SNPs) of genes involved in SAM synthesis and NSCL/P in a Chilean population.

METHODS

In 234 cases and 309 controls, 18 SNPs in AHCY, MTR, MTRR, and MAT2A were genotyped, and the association between them and the phenotype was evaluated based on additive (allele), dominant, recessive and haplotype models, by odds ratio (OR) computing.

RESULTS

Three deep intronic SNPs of MTR showed a protective effect on NSCL/P expression: rs10925239 (OR 0.68; p = 0.0032; q = 0.0192), rs10925254 (OR 0.66; p = 0.0018; q = 0.0162), and rs3768142 (OR 0.66; p = 0.0015; q = 0.0162). Annotations in expression database demonstrate that the protective allele of the three SNPs is associated with a reduction of MTR expression summed to the prediction by bioinformatic tools of its potentiality to modify splicing sites.

CONCLUSIONS

The protective effect against NSCL/P of these intronic MTR SNPs seems to be related to a decrease in MTR enzyme expression, modulating the SAM availability for proper substrate methylation. However, functional analyses are necessary to confirm our findings.

IMPACT

SAM synthesis pathway genetic variants are factors associated to NSCL/P. This article adds new evidence for folate related genes in NSCL/P in Chile. Its impact is to contribute with potential new markers for genetic counseling.

A LCMS-based untargeted lipidomics analysis of cleft palate in mouse

BACKGROUND

Recent studies have shown that lipid metabolism was abnormal during the formation of cleft palate. However, the composition of these lipid species remains unclear.

OBJECTIVE

Aims of this study were to identify the lipid species components and reveal the key lipid metabolic disorders in cleft palate formation.

METHODS

The pregnant mice were divided into experimental group exposed to all-trans retinoic acid (RA-treated group) (n = 12) and control group (n = 12) at embryonic gestation day 10.5 (E0.5). The component of the palatal tissue metabolome was analyzed using a LCMS-based nontargeted lipidomics approach. Multivariate statistical analysis was then carried out to assess the differences between the RA-treated group and the control group.

RESULTS

Twenty-nine lipid species were found to discriminate between RA-treated and control embryos. Among them, 28 lipid species increased and 1 lipid species decreased in the RA-treated group. Among these lipids, 13 were triglycerides, 9 were PEs, 3 were PCs, 2 were PSs, 2 were DGs. Further analysis of the number of carbons and unsaturated bond of triglycerides showed that TGs with high unsaturated bonds constituted a higher fraction in the RA-treated group. A higher amount of triglycerides containing 52, 54, 56, 58, 60 carbons, and 1 to 8 unsaturated bonds. Of note, under RA treatment, TG 50:1, 52:2, 56:6and 60:8 became the most prominent.

CONCLUSION

Lipid metabolism is significantly different in the formation of cleft palate induced by RA, and the unsaturated triglycerides increased in the RA-treated group may play an important role in the formation of cleft palate.

Accumulation of rare coding variants in genes implicated in risk of human cleft lip with or without cleft palate

Cleft lip with/without cleft palate (CLP) is a common craniofacial malformation with complex etiologies, reflecting both genetic and environmental factors. Most of the suspected genetic risk for CLP has yet to be identified. To further classify risk loci and estimate the contribution of rare variants, we sequenced the exons in 49 candidate genes in 323 CLP cases and 211 nonmalformed controls. Our findings indicated that rare, protein-altering variants displayed markedly higher burdens in CLP cases at relevant loci. First, putative loss-of-function mutations (nonsense, frameshift) were significantly enriched among cases: 13 of 323 cases (~4%) harbored such alleles within these 49 genes, versus one such change in controls (p = 0.01). Second, in gene-level analyses, the burden of rare alleles showed greater case-association for several genes previously implicated in cleft risk. For example, BHMT displayed a 10-fold increase in protein-altering variants in CLP cases (p = .03), including multiple case occurrences of a rare frameshift mutation (K400 fs). Other loci with greater rare, coding allele burdens in cases were in signaling pathways relevant to craniofacial development (WNT9B, BMP4, BMPR1B) as well as the methionine cycle (MTRR). We conclude that rare coding variants may confer risk for isolated CLP.