X-linked markers in the Duchenne muscular dystrophy gene associated with oral clefts

Genome-wide Gene-by-Sex Interaction Studies Identify Novel Nonsyndromic Orofacial Clefts Risk Locus

Risk loci identified through genome-wide association studies have explained about 25% of the phenotypic variations in nonsyndromic orofacial clefts (nsOFCs) on the liability scale. Despite the notable sex differences in the incidences of the different cleft types, investigation of loci for sex-specific effects has been understudied. To explore the sex-specific effects in genetic etiology of nsOFCs, we conducted a genome-wide gene × sex (GxSex) interaction study in a sub-Saharan African orofacial cleft cohort. The sample included 1,019 nonsyndromic orofacial cleft cases (814 cleft lip with or without cleft palate and 205 cleft palate only) and 2,159 controls recruited from 3 sites (Ethiopia, Ghana, and Nigeria). An additive logistic model was used to examine the joint effects of the genotype and GxSex interaction. Furthermore, we examined loci with suggestive significance (P < 1E-5) in the additive model for the effect of the GxSex interaction only. We identified a novel risk locus on chromosome 8p22 with genome-wide significant joint and GxSex interaction effects (rs2720555, p_2df = 1.16E-08, p_GxSex = 1.49E-09, odds ratio [OR] = 0.44, 95% CI = 0.34 to 0.57). For males, the risk of cleft lip with or without cleft palate at this locus decreases with additional copies of the minor allele (p < 0.0001, OR = 0.60, 95% CI = 0.48 to 0.74), but the effect is reversed for females (p = 0.0004, OR = 1.36, 95% CI = 1.15 to 1.60). We replicated the female-specific effect of this locus in an independent cohort (p = 0.037, OR = 1.30, 95% CI = 1.02 to 1.65), but no significant effect was found for the males (p = 0.29, OR = 0.86, 95% CI = 0.65 to 1.14). This locus is in topologically associating domain with craniofacially expressed and enriched genes during embryonic development. Rare coding mutations of some of these genes were identified in nsOFC cohorts through whole exome sequencing analysis. Our study is additional proof that genome-wide GxSex interaction analysis provides an opportunity for novel findings of loci and genes that contribute to the risk of nsOFCs.

Genome-wide interaction studies identify sex-specific risk alleles for nonsyndromic orofacial clefts

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect in humans and is notable for its apparent sexual dimorphism where approximately twice as many males are affected as females. The sources of this disparity are largely unknown, but interactions between genetic and sex effects are likely contributors. We examined gene-by-sex (G × S) interactions in a worldwide sample of 2,142 NSCL/P cases and 1,700 controls recruited from 13 countries. First, we performed genome-wide joint tests of the genetic (G) and G × S effects genome-wide using logistic regression assuming an additive genetic model and adjusting for 18 principal components of ancestry. We further interrogated loci with suggestive results from the joint test ( p < 1.00 × 10 -5 ) by examining the G × S effects from the same model. Out of the 133 loci with suggestive results ( p < 1.00 × 10 -5 ) for the joint test, we observed one genome-wide significant G × S effect in the 10q21 locus (rs72804706; p = 6.69 × 10 -9 ; OR = 2.62 CI [1.89, 3.62]) and 16 suggestive G × S effects. At the intergenic 10q21 locus, the risk of NSCL/P is estimated to increase with additional copies of the minor allele for females, but the opposite effect for males. Our observation that the impact of genetic variants on NSCL/P risk differs for males and females may further our understanding of the genetic architecture of NSCL/P and the sex differences underlying clefts and other birth defects.

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.

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.

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.

Family-Based Multi-SNP X Chromosome Analysis Using Parent Information

We propose a method for association analysis of haplotypes on the X chromosome that offers both improved power and robustness to population stratification in studies of affected offspring and their parents if all three have been genotyped. The method makes use of assumed parental haplotype exchangeability (PHE), a weaker assumption than Hardy-Weinberg equilibrium (HWE). PHE requires that in the source population, of the three X chromosome haplotypes carried by the two parents, each is equally likely to be carried by the father. We propose a pseudo-sibling approach that exploits that exchangeability assumption. Our method extends the single-SNP PIX-LRT method to multiple SNPs in a high linkage block. We describe methods for testing the PHE assumption and also for determining how apparent violations can be distinguished from true fetal effects or maternally-mediated effects. We show results of simulations that demonstrate nominal type I error rate and good power. The methods are then applied to dbGaP data on the birth defect oral cleft, using both Asian and Caucasian families with cleft.

Learning about the X from our parents

The X chromosome is generally understudied in association studies, in part because the analyst has had limited methodological options. For nuclear-family-based association studies, most current methods extend the transmission disequilibrium test (TDT) to the X chromosome. We present a new method to study association in case-parent triads: the parent-informed likelihood ratio test for the X chromosome (PIX-LRT). Our method enables estimation of relative risks and takes advantage of parental genotype information and the sex of the affected offspring to increase statistical power to detect an effect. Under a parental exchangeability assumption for the X, if case-parent triads are complete, the parents of affected offspring provide an independent replication sample for estimates based on transmission distortion to their affected offspring. For each offspring sex we combine the parent-level and the offspring-level information to form a likelihood ratio test statistic; we then combine the two to form a combined test statistic. Our method can estimate relative risks under different modes of inheritance or a more general co-dominant model. In triads with missing parental genotypes, the method accounts for missingness with the Expectation-Maximization algorithm. We calculate non-centrality parameters to assess the power gain and robustness of our method compared to alternative methods. We apply PIX-LRT to publically available data from an international consortium of genotyped families affected by the birth defect oral cleft and find a strong, internally-replicated signal for a SNP marker related to cleft lip with or without cleft palate.

Classification, epidemiology, and genetics of orofacial clefts

Orofacial clefts (OFCs) include a broad range of facial conditions that differ in cause and disease burden. In the published literature, there is substantial ambiguity in both terminology and classification of OFCs. This article discusses the terminology and classification of OFCs and the epidemiology of OFCs. Demographic, environmental, and genetic risk factors for OFCs are described, including suggestions for family counseling. This article enables clinicians to counsel families regarding the occurrence and recurrence of OFCs. Although much of the information is detailed, it is intended to be accessible to all health professionals for use in their clinical practices.