Family-based genome-wide association study in Patagonia confirms the association of theDMDlocus and cleft lip and palate

MMP16 as NSCL ± P Susceptible Gene in Western Han Chinese

OBJECTIVE

The role of MMP16 in lip development is unclear. This study aimed to identify nonsyndromic cleft lip with or without palate (NSCL ± P) susceptible loci of MMP16 in western Han Chinese.

DESIGN

We performed targeted sequencing around MMP16 combined with a 2-phase association analysis on common variants. Phase 2 association analysis was performed with NSCL ± P specific subphenotypes (NSCL and NSCLP). Then we used rare variants burden analysis and genotyping, accompanied by motif analysis.

SETTING

This study was completed in a tertiary medical center.

PATIENTS, PARTICIPANTS

Phase 1 targeted sequencing included 159 patients with NSCL ± P and 542 normal controls; phase 2 included 1626 patients with NSCL ± P (1047 NSCL and 579 NSCLP) and 2255 normal controls.

INTERVENTIONS

Venous blood samples were collected from patients and used to extract DNA.

MAIN OUTCOME MEASURES

After Bonferroni correction, phase 1 significant threshold of p-value was 4.28 × 10-5 (0.05/1167 single nucleotide polymorphisms [SNPs]), and phase 2 was .00025 (0.05/200 SNPs). Burden analysis significant threshold p-value was .05.

RESULTS

Common variants phase 1 association analysis identified 11 statistically significant SNPs (lowest p = 1.90 × 10-9, odds ratio (OR) = 0.27, 95% CI: 0.17-0.44), phase 2 replication identified 16 SNPs in NSCL ± P (lowest p = 6.26 × 10-6, OR = 0.77, 95% CI: 0.69-0.86) and 9 in NSCL (lowest p = 8.44 × 10-5, OR = 0.76, 95% CI: 0.66-0.87). Rare variants burden analysis showed no significant results, genotyping results showed they were maternally inherited.

CONCLUSIONS

Our study identified MMP16 susceptible SNPs in NSCL ± P and NSCL, emphasizing its potential role in lip development. Our study also highlighted the importance to perform association analysis with subphenotypes divided.

Association Between SNPs in 1q32.2 and NSCL ± P in Han Chinese Population

OBJECTIVES

Non-syndromic cleft lip with or without cleft palate (NSCL ± P) is one of the most common birth malformations. Currently, numerous susceptibility SNPs have been reported by GWA studies, however, the replications of them among NSCL ± P from Han Chinese were very limited.

DESIGN

In this study, we selected 16 SNPs around 1q32.2 based on the published GWA studies and replicated them among 302 trios with NSCL ± P from Han Chinese Population. The genotypic data was analyzed with FBAT, PLINK and R package.

SETTING

The study was conducted in a tertiary medical center.

PATIENTS, PARTICIPANTS

302 patients with CL ± P and their parents.

MAIN OUTCOME MEASURES

To ascertain the genetic variants in 1q32.2 in patients with CL ± P in Han Chinese Population.

INTERVENTIONS

Blood samples were collected.

RESULTS

We found T allele (Z = 4.26, p = 0.00002) and T/T homozygotes (Z = 4.4, p = 0.000011) at rs12063989 was significantly over-transmitted among non-syndromic cleft lip with or without cleft palate (NSCL ± P).

CONCLUSIONS

We found rs12063989 exhibited significant association with the occurrence of NSCL ± P, which would provide new evidence for the future study in the etiology of NSCL ± P.

SNPs at TP63 gene was specifically associated with right-side cleft lip in Han Chinese population

OBJECTIVES

Non-syndromic cleft lip with or without palate is one of the most common birth malformations. TP63 and GREM1 were recently reported to be associated with NSCL/P. However, there were few studies focused on their associations in non-syndromic cleft lip only (NSCLO).

DESIGN

Initial screening and replication in large cohorts were used to locate the susceptible SNPs of TP63 and GREM1. Firstly, variations were screened among 192 NSCLO cases by the Sanger sequencing. Then, we selected five associated SNPs in initial screening phase and replicated among 1,006 NSCLO cases and 1,823 normal controls.

RESULTS

Initial chi-square test showed that rs7653848, rs7624324, rs6790167, and rs1345186 in TP63 and rs2280738 in GREM1 achieved statistical significance (p < .05); the subsequent replication analysis showed that rs1345186 was specifically significant in right-side cleft lip (RCL; p = .017, OR = 1.33, and 95% CI: 1.05-1.69).

CONCLUSION

This study firstly used the subphenotype of cleft lip samples to verify the association between TP63 and GREM1, which indicated that TP63 is a promising susceptible gene for RCL in Chinese population. And further confirmed the different etiology in the right-sided cleft lip, left-sided cleft lip, and bilateral cleft lip of NSCLO. This will give new reference for the future research and genetic counseling.

Distinct Expression of miR-378 in Nonsyndromic Cleft Lip and/or Cleft Palate: A Cogitation of Skewed Sex Ratio in Prevalence

OBJECTIVE

Nonsyndromic cleft lip and/or cleft palate (NSCL/P) is an isolated phenotype of orofacial clefts with skewed sex ratio in prevalence. This study aims to identify differentially expressed genes (DEGs) and microRNAs (DEMs) of NSCL/P by integrated bioinformatics analysis, revealing mechanisms for sexual dimorphism in prevalence.

MATERIALS AND METHODS

First, we downloaded the expression profile data from Gene Expression Omnibus database to identify DEGs and DEMs. Second, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses performed DEGs' functions. Then, clustered DEGs were identified through protein-protein interaction networks. Combining clustered DEGs with key genes searched in GeneCards enlarged NSCL/P-related genes. Moreover, the genes were linked by transcription factors (TFs). Subsequently, connected by the above TFs, DEMs and genes were used to establish the miRNA-TF-messenger RNA (mRNA) regulatory networks.

RESULTS

The DEGs in sex-ignored group, female-only group, and male-only group were obtained, respectively. Among the DEMs, miR-378 was downregulated in females but upregulated in males. In female-only group, the miRNA-TF-mRNA regulatory networks showed miR-378-SP1-POLE2/CDK6/EZR regulatory axis was found to be key candidates of NSCL/P.

CONCLUSIONS

Our findings suggest that different expression of miR-378 is consistent with the skewed sex ratio in the prevalence of NSCL/P.

Homozygote C/C at rs12543318 was risk factor for non-syndromic cleft lip only from Western Han Chinese population

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a complex disorder, and it results from both of the genetic modifiers and environmental factors, with genetic modifiers contributes to it more than environmental factors. GWASs made great progress in identifying the candidate genes for NSCL/P, but the findings need to be replicated in other populations. In this study, we selected eleven SNPs from recent GWASs and GWAS meta-analysis to investigate their associations among 308 NSCL/P trios (134 non-syndromic cleft lip only (NSCLO) trios and 174 non-syndromic cleft lip with cleft palate (NSCLP) trios) from Han Chinese population. All SNPs were genotyped using SNPscan method and analyzed the data with FBAT, PLINK, and R package. Allelic TDT analysis showed that allele A at rs12543318 was associated with NSCLO trios (P = .0032, OR = 0.57, 95% CI: 0.39-0.83), and parent-of-origin effect analysis indicated that allele A at rs12543318 was significantly maternally undertransmitted among NSCLO (P = .0046), which implied the potential influence of genomic imprinting; global TDT further confirmed this association. Individual genotypic TDT showed homozygote C/C at rs12543318 was overtransmitted among NSCLO (Z = 3.79, P = .00015) and NSCL/P groups (Z = 3.83, P = .00013), which indicated that it could increase the risk to have cleft babies. Our findings indicated that rs12543318 was associated with NSCLO from Western Han Chinese population, which will give new scientific evidence for later researches in the etiology of NSOCs.

Analysis of Parent-of-Origin Effects on the X Chromosome in Asian and European Orofacial Cleft Triads Identifies Associations with DMD, FGF13, EGFL6, and Additional Loci at Xp22.2

Background: Although both the mother's and father's alleles are present in the offspring, they may not operate at the same level. These parent-of-origin (PoO) effects have not yet been explored on the X chromosome, which motivated us to develop new methods for detecting such effects. Orofacial clefts (OFCs) exhibit sex-specific differences in prevalence and are examples of traits where a search for various types of effects on the X chromosome might be relevant. Materials and Methods: We upgraded our R-package Haplin to enable genome-wide analyses of PoO effects, as well as power simulations for different statistical models. 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European case-parent triads of isolated OFCs were available from a previous GWAS. For each ethnicity, cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO) were analyzed separately using two X-inactivation models and a sliding-window approach to haplotype analysis. In addition, we performed analyses restricted to female offspring. Results: Associations were identified in "Dystrophin" (DMD, Xp21.2-p21.1), "Fibroblast growth factor 13" (FGF13, Xq26.3-q27.1) and "EGF-like domain multiple 6" (EGFL6, Xp22.2), with biologically plausible links to OFCs. Unlike EGFL6, the other associations on chromosomal region Xp22.2 had no apparent connections to OFCs. However, the Xp22.2 region itself is of potential interest because it contains genes for clefting syndromes [for example, "Oral-facial-digital syndrome 1" (OFD1) and "Midline 1" (MID1)]. Overall, the identified associations were highly specific for ethnicity, cleft subtype and X-inactivation model, except for DMD in which associations were identified in both CPO and CL/P, in the model with X-inactivation and in Europeans only. Discussion/Conclusion: The specificity of the associations for ethnicity, cleft subtype and X-inactivation model underscores the utility of conducting subanalyses, despite the ensuing need to adjust for additional multiple testing. Further investigations are needed to confirm the associations with DMD, EGF16, and FGF13. Furthermore, chromosomal region Xp22.2 appears to be a hotspot for genes implicated in clefting syndromes and thus constitutes an exciting direction to pursue in future OFCs research. More generally, the new methods presented here are readily adaptable to the study of X-linked PoO effects in other outcomes that use a family-based design.

Third molar agenesis as a potential marker for craniofacial deformities

The identification of clinical patterns of tooth agenesis in individuals born with craniofacial deformities may be a useful tool for risk determination of these defects. We hypothesize that specific craniofacial deformities are associated with third molar agenesis.

OBJECTIVE

The aim of this study was to identify if third molar agenesis could have a relation with other craniofacial structure alterations, such as cleft lip and palate, skeletal malocclusion, or specific growth patterns in humans.

DESIGN

Data were obtained from 550 individuals ascertained as part of studies aiming to identify genetic contributions to oral clefts. 831 dental records of patients aged over eight years seeking orthodontic treatment were also included. SN-GoGn angle were used to classify the growth pattern (hypo-divergent, normal and hyper-divergent), and the ANB angle was used to verify the skeletal malocclusion pattern (Class I, II and III). Panoramic radiographs were used to determine third molar agenesis.

RESULTS

A high frequency of third molar agenesis among individuals born with cleft lip with or without cleft palate (55%), as well as among their relatives (93.5%) was found. Third molar agenesis was not associated to skeletal malocclusion or growth pattern.

CONCLUSION

It appears that third molar agenesis is associated with the disturbances that lead to cleft lip and palate.

A new approach to chromosome-wide analysis of X-linked markers identifies new associations in Asian and European case-parent triads of orofacial clefts

Background

GWAS discoveries on the X-chromosome are underrepresented in the literature primarily because the analytical tools that have been applied were originally designed for autosomal markers. Our objective here is to employ a new robust and flexible tool for chromosome-wide analysis of X-linked markers in complex traits. Orofacial clefts are good candidates for such analysis because of the consistently observed excess of females with cleft palate only (CPO) and excess of males with cleft lip with or without cleft palate (CL/P).

Methods

Genotypes for 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European isolated cleft triads were available from a previously published GWAS. The R-package HAPLIN enables genome-wide–level analyses as well as statistical power simulations for a range of biologic scenarios. We analyzed isolated CL/P and isolated CPO for each ethnicity in HAPLIN, using a sliding-window approach to haplotype analysis and two different statistical models, with and without X-inactivation in females.

Results

There was a larger number of associations in the Asian versus the European sample, and similar to previous reports that have analyzed the same GWAS dataset using different methods, we identified associations with EFNB1/PJA1 and DMD. In addition, new associations were detected with several other genes, among which KLHL4, TBX22, CPXCR1 and BCOR were noteworthy because of their roles in clefting syndromes. A few of the associations were only detected by one particular X-inactivation model, whereas a few others were only detected in one sex.

Discussion/Conclusion

We found new support for the involvement of X-linked variants in isolated clefts. The associations were specific for ethnicity, sex and model parameterization, highlighting the need for flexible tools that are capable of detecting and estimating such effects. Further efforts are needed to verify and elucidate the potential roles of EFNB1/PJA1, KLHL4, TBX22, CPXCR1 and BCOR in isolated clefts.

Craniofacial anomalies associated with hypospadias. Description of a hospital based population in South America

Introduction:

Hypospadias is a congenital abnormality of the penis, in which there is incomplete development of the distal urethra. There are numerous reports showing an increase of prevalence of hypospadias. Association of craniofacial malformations in patients diagnosed with hypospadias is rare. The aim of this study is to describe the association between hypospadias and craniofacial congenital anomalies.

Materials and Methods:

A retrospective review of the Latin-American collaborative study of congenital malformations (ECLAMC) data was performed between January 1982 and December 2011. We included children diagnosed with associated hypospadias and among them we selected those that were associated with any craniofacial congenital anomaly.

Results:

Global prevalence was 11.3 per 10.000 newborns. In this population a total of 809 patients with 1117 associated anomalies were identified. On average there were 1.7 anomalies per patient. Facial anomalies were present in 13.2%. The most commonly major facial anomaly associated to hypospadias was cleft lip/palate with 52 cases. We identified that 18% have an association with other anomalies, and found an association between craniofacial anomalies and hypospadias in 0.59 cases/10.000 newborns.

Discussion:

Hypospadias is the most common congenital anomaly affecting the genitals. Its association with other anomalies is rare. It has been reported that other malformations occur in 29.3% of the cases with hypospadias. The more proximal the meatus, the higher the risk for having another associated anomaly.

Conclusion:

Associated hypospadias are rare, and it is important to identify the concurrent occurrence of craniofacial anomalies to better treat patients that might need a multidisciplinary approach.

Analysis of the genetic ancestry of patients with oral clefts from South American admixed populations

Increased susceptibility to cleft lip, with or without cleft palate (CL±P) has been observed in South America, as related to Amerindian ancestry, using epidemiological data, uniparental markers, and blood groups. In this study, it was evaluated whether this increased risk remains when Amerindian ancestry is estimated using autosomal markers and considered in the predictive model. Ancestry was estimated through genotyping 62 insertion and deletion (INDEL) markers in sample sets of patients with CL±P, patients with cleft palate (CP), and controls, from Patagonia in southern Argentina and Belém in northern Brazil. The Amerindian ancestry in patients from Patagonia with CL±P was greater than in controls although it did not reach statistical significance. The European ancestry in patients with CL±P from Belém and in patients with CP from Belém and Patagonia was higher than in controls and statistically significant for patients with CP who were from Belém. This high contribution of European genetic ancestry among patients with CP who were from Belém has not been previously observed in American populations. Our results do not corroborate the currently accepted risks for CL±P and CP estimated by epidemiological studies in the North American populations and probably reflect the higher admixture found in South American ethnic groups when compared with the same ethnic groups from the North American populations.