Øivind Skare1†
Rolv T. Lie2,3†Øystein A. Haaland2Miriam Gjerdevik2,4Julia Romanowska2,5
Håkon K. Gjessing2,3‡
Astanand Jugessur2,3,4*‡- 1Department of Occupational Medicine and Epidemiology, National Institute of Occupational Health, Oslo, Norway
- 2Department of Global Public Health and Primary Care, University of Bergen, Bergen, Norway
- 3Centre for Fertility and Health (CeFH), Norwegian Institute of Public Health, Oslo, Norway
- 4Department of Genetics and Bioinformatics, Norwegian Institute of Public Health, Oslo, Norway
- 5Computational Biology Unit, University of Bergen, Bergen, Norway
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.
Introduction
Orofacial clefts (OFCs) are the most common craniofacial birth defects in humans, affecting approximately 1–2/1000 live births worldwide (Mossey et al., 2009). They are broadly categorized according to whether the affected region includes the primary palate, the secondary palate, or both, and whether they occur with or without additional congenital anomalies (Dixon et al., 2011; Marazita, 2012; Beaty et al., 2016). OFCs pose a substantial public health burden in terms of the medical costs and sequelae associated with their treatment, which may persist from infancy to childhood and throughout life (Wehby and Cassell, 2010; Berg et al., 2016a,b). OFCs have also been linked to higher risk of specific types of cancer in later life (Zhu et al., 2002; Bille et al., 2005), increased overall mortality well into adulthood (Christensen et al., 2004), and lower academic achievement (Wehby et al., 2014). Despite major progress in surgery and other medical interventions aimed at repairing and managing the cleft itself, current understanding of the biological underpinnings of these relatively common birth defects is still incomplete.
Both genetic and environmental factors have been reported to influence the risk of OFCs, either individually or through their complex interactions in relevant biological pathways (Jugessur et al., 2009a; Dixon et al., 2011; Marazita, 2012; Rahimov et al., 2012; Beaty et al., 2016; Kousa and Schutte, 2016). Genome-wide association studies (GWAS) have successfully identified several genes and loci for OFCs, contributing to an improved understanding of the biological processes underlying these relatively common birth defects (Birnbaum et al., 2009; Grant et al., 2009; Beaty et al., 2010; Mangold et al., 2010; Camargo et al., 2012; Ludwig et al., 2012; Wolf et al., 2015; Leslie et al., 2016a,c). However, most of the genetic variation in OFCs remains unexplained. Given the more than 30-fold increased risk of recurrence in first-degree relatives of patients with OFCs (Sivertsen et al., 2008; Grosen et al., 2010), exploring alternative genetic mechanisms beyond simple fetal or maternal gene-effects alone may be important. This entails investigating higher-order interactions, such as epistasis (Cordell, 2009; Wei et al., 2014) and gene-environment interaction (Thomas, 2010), and studying parent-of-origin (PoO) effects (Ferguson-Smith, 2011; Guilmatre and Sharp, 2012; Lawson et al., 2013; Peters, 2014; Connolly and Heron, 2015; Gjerdevik et al., 2017; Haaland et al., 2017).
A PoO effect describes the situation where the effect of an allele in the offspring differs according to the parental origin of the allele (Guilmatre and Sharp, 2012; Gjerdevik et al., 2017). PoO effects are particularly relevant for birth defects because the mother influences the development of the fetus through the action of her own genes and through providing the prenatal environment for the fetus. To estimate PoO effects, one contrasts the frequency of alleles transmitted to an affected offspring from the mother versus the father, and if transmission distortion to the affected offspring is stronger for mothers than fathers (or vice versa), there is evidence of a PoO effect (Weinberg, 1999; Jugessur et al., 2012b; Connolly and Heron, 2015). PoO effects might also occur on the X chromosome. Given 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), OFCs are good examples of traits that might be caused by various types of effects on the X chromosome (including PoO effects). Moreover, several genes on the X chromosome are known to cause syndromic forms of clefts, and there is growing evidence that X-linked genes might also contribute to isolated clefts (Jugessur et al., 2012a; Patel et al., 2013; Fonseca et al., 2015; Wise et al., 2016; Skare et al., 2017).
A major shortcoming of most previous GWAS has been the systematic exclusion of SNPs on the X chromosome prior to analysis, even though this chromosome comprises approximately 5% of the human genome and many genotyping platforms do include X-linked SNPs. This has led to only a few analyses of X-linked markers for complex traits in general (Wise et al., 2013), and for OFCs in particular, only two studies have explored PoO effects at the genome-wide level (Shi et al., 2012; Garg et al., 2014), and neither investigated PoO effects on the X chromosome.
With these gaps in mind, we upgraded our R-package Haplin (Gjessing and Lie, 2006) with new functionalities to enable PoO analyses at the genome-wide level, as well as an assessment of statistical power for different statistical models. The current analyses are based on the largest collection of case-parent triads of OFCs to date (Beaty et al., 2010). We implemented a sliding-window approach to haplotype analysis and used two X-inactivation models, one with and one without the assumption of X-inactivation in females (Yang et al., 2011), to explore PoO effects on the X chromosome. In addition, we performed separate analyses on female offspring alone to allow for possible sex-specific differences.
Materials and Methods
Study Populations
Characteristics of the study populations, the genotyping platform, and the quality control criteria used for data cleaning have been detailed elsewhere (Beaty et al., 2010; Shi et al., 2012; Patel et al., 2013; Skare et al., 2017). Briefly, genotyping was performed on an Illumina Human610-Quad® platform and genotypes for 589,945 SNPs (99.56% of the attempted SNPs) were released and later deposited in the Database of Genotypes and Phenotypes (dbGaP; http://www.ncbi.nlm.nih.gov/gap) under study accession ID phs000094.v1.p1. Genotypes for 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European case-parent triads of isolated clefts were available for the current analyses. For additional data cleaning, we used PLINK (Purcell et al., 2007) to exclude individuals with more than 10% missing genotypes, SNPs with more than 1% missing genotypes, and SNPs with a minor allele frequency (MAF) less than 0.01. Mendelian errors were not assessed during data cleaning; however, SNPs with more than 30 Mendelian errors were excluded at the analysis stage. After data cleaning, 13,180 X-chromosome SNPs were available for the current PoO analyses. Table 1 shows the number of triads according to ethnicity, cleft subtype and child's sex. For each ethnicity, isolated CL/P and isolated CPO triads were analyzed separately.
Table 1. Number of case-parent triads according to ethnicity, cleft subtype, and child's sex.
Statistical Methods
We recently published a new approach to chromosome-wide analysis of X-linked SNPs using the same dataset as here (Skare et al., 2017). In the current paper, we extend the approach to also cover PoO effects on the X chromosome, as implemented in our statistical software Haplin (Gjessing and Lie, 2006). Two types of analyses were performed on the Asian and European samples: (i) Single-marker analyses, where SNPs were analyzed individually, and (ii) haplotype analyses, where up to four SNPs per sliding window were analyzed together. Haplin fits a log-linear model to genotype data from case-parent triads. It reconstructs haplotypes from multiple-SNP data and estimates the relative risk and confidence interval for one or two copies of a target allele or haplotype. Since the expectation-maximization algorithm is implemented in Haplin, incomplete triads can also be used in the analyses after accounting for missing parental genotypes in the maximum likelihood estimation.
For haplotype analysis, we used the HaplinSlide function in Haplin. For additional information on data formats and the HaplinSlide function, see our recent chromosome-wide analysis of X-linked SNPs (Skare et al., 2017), our website at https://people.uib.no/gjessing/genetics/software/haplin/, or the R help pages for Haplin at https://CRAN.R-project.org/package=Haplin.
Estimating PoO Effects on the X Chromosome Using Haplin
Within the maximum-likelihood framework of Haplin, autosomal PoO effects are derived essentially by contrasting the relative frequency of the variant allele when transmitted from the mother to the affected child versus the relative frequency of the variant allele when transmitted from the father; there is evidence of a PoO effect if the frequencies differ (Weinberg, 1999; Gjerdevik et al., 2017). Haplin provides estimates of RRm and RRf, which are the relative risk increase (or decrease) associated with inheriting the variant allele from the mother and from the father, respectively. The measure of PoO effect is then the ratio of relative risks, RRR = RRm/RRf, i.e., the PoO effect is a comparison of the relative risks derived from the maternally and paternally inherited alleles. A value of RRR = 1 would be obtained from an allele with no PoO effect, i.e., RRm = RRf. Note that even when RRR = 1, RRm and RRf could themselves still be different from 1, meaning that there could be an effect of the allele carried by the fetus even in the absence of a PoO effect. Under a multiplicative model, the risk is assumed to be RRm*RRf when the allele is inherited both from the mother and the father (i.e., a double dose of the variant allele).
Similarly, on the X chromosome, the estimates of RRm and RRf can be obtained by restricting the analyses to girls only. However, since the X chromosome in boys is maternally derived, Haplin allows an increase in power by combining the relative risk estimate RRB from boys with the estimates from girls. The two options for doing so assume either X-inactivation among girls, where RRm*RRf = RRB, or no X-inactivation, where RRm = RRB. Note that in the model assuming X-inactivation, equating the estimate RRB from the boys with the product of RRm and RRf may influence the RRm and RRf estimates themselves, but has relatively little influence on the ratio RRR. When studying the significance of PoO specifically, the hits obtained from investigating girls only will thus be similar to those obtained from the model assuming X-inactivation. The combined model assuming no X-inactivation, however, may well produce other hits. Both combined models allow different baselines risks for boys and girls to be fitted.
Figure 1 provides a detailed explanation of the Haplin model for PoO effects on the X chromosome.
Figure 1. An illustration of the Haplin model for parent-of-origin (PoO) effects on the X chromosome. The red arrows show the relative risks associated with girls inheriting the risk allele “a” from the mother (RRm) or from the father (RRf). Under the multiplicative risk model illustrated here, the relative risk increase for “aa” girls, i.e., girls inheriting allele “a” from both the mother and the father, is RRaa = RRm*RRf. The ratio RRm/RRf = 1.4 is a measure of the PoO effect (blue arrow). The risk increase for boys when inheriting the “a” allele is RRB = 2.0. Under the assumption of X-inactivation (A), the risk increase for girls inheriting “aa” is the same as that for a single “a” in boys, i.e., RRaa = RRB. Under the assumption of no X-inactivation (B), the risk increase for girls inheriting “a” from the mother is the same as that for boys, i.e., RRm = RRB. In the model without X-inactivation, when a girl inherits “a” from the father as well, this may lead to a higher total risk increase for girls than for boys inheriting the one “a” from the mother. In this illustration, RRaa = RRm*RRf = 2.86 > RRB. The model allows different baseline risks for girls and for boys, here 0.4 and 0.5%, respectively.
Post-processing of Results
For each cleft subtype, –log10(p) were plotted against chromosomal position in a Manhattan plot, before all the Manhattan plots were collated and displayed in a single figure. To control for the proportion of falsely rejected hypotheses, we applied a false discovery rate (FDR) method where the original p-values were replaced by “q-values” (Storey and Tibshirani, 2003). For example, among SNPs with a q-value ≤ 0.2, one would expect an FDR of less than 20%.
Electronic Database Information
Haplin version 6.2.1 is implemented as a package in the R statistical software (R Development Core Team, 2014) and can be installed from the CRAN library. More information can be found at our web site (https://people.uib.no/gjessing/genetics/software/haplin/).
Ethics Approval
Ethics approvals for the International Cleft Consortium were obtained from the respective institutional review boards of the participating sites. The consortium was formed in 2007 and each participating institution approved research protocols for the recruitment of case-parent triads from 13 individual sites. All participants have granted their written informed consents. The participating sites included institutions in the US (Johns Hopkins University; University of Iowa; Utah State University; National Institute of Environmental Health Sciences (NIEHS); University of Pittsburgh), Denmark (University of Southern Denmark), Norway (University of Bergen), China (Peking University Health Science Center; Wuhan University; Peking Union Medical College; West China School of Stomatology, Sichuan University; School of Stomatology, Beijing University), Korea (Yonsei University), Taiwan (Chang Gung Memorial Hospital), and Singapore (KK Women's & Children's Hospital; National University of Singapore). For more details on the recruitment sites, the research approvals and protocols, see the online “Supplementary Note” of the original publication (Beaty et al., 2010), as well as the study outline at dbGAP (https://www.ncbi.nlm.nih.gov/gap), under study accession number phs000094.v1.p1.
Results
The results are organized in two main parts: (Part A) presents the results of the analyses of the boys and girls together (Figure 2), and (Part B) presents the results of the analyses of the girls only (Figure 3). For each of the X-inactivation models depicted in Figure 2, the results of the single-marker analyses are presented first, followed by the results of the haplotype analyses. Table 2 provides additional information on the SNPs and haplotypes, along with their relative risks (RRs) and 95% confidence intervals (CIs). Note that the Manhattan plots for the single-marker and haplotype analyses (Figures 2, 3) only show the lead SNPs and haplotypes. A more complete list of SNPs and haplotypes lying above the p-value cutoff of 10−4 in Figures 2, 3 is provided in Table 2.
Figure 2. Single-marker and haplotype analyses in the Asian and European samples without stratification by child's sex. The Manhattan plots show the single-marker and haplotype analyses based on the model without and with X-inactivation in females, respectively. For convenience, we have added a vertical line corresponding to a Bonferroni-corrected p-value cutoff of 10−4.
Figure 3. Single-marker and haplotype analyses in the Asian and European girls only. The Manhattan plots show the single-marker and haplotype analyses of the girls only. For convenience, we have added a vertical line corresponding to a Bonferroni-corrected p-value cutoff of 10−4.
Table 2. Q-values and RR estimates for associations identified with SNPs and haplotypes.
Figure 4 shows the results of power simulations for different statistical models (girls only, without X-inactivation, with X-inactivation), based on sample sizes reflecting those available in the current GWAS dataset. Table 3 provides a synopsis of all the genes in which associations were identified in the current analyses. To determine whether a given SNP or haplotype was located in or near a gene within 20 kb, we used the 1,000 Genomes browser (https://www.ncbi.nlm.nih.gov/variation/tools/1000genomes).
Figure 4. Statistical power for a single SNP. (A) Shows the simulated power for the model restricted to girls only, (B) shows the simulated power for the model without X-inactivation, and (C) shows the simulated power for the model with X-inactivation. The figures on the left depict the power with increasing relative risk ratios (RRRs) (increasing values of RRm; RRf = 1) and differing sample sizes of case-parent triads, using a fixed minor allele frequency (MAF) of 0.2. The figures on the right-hand side show the power with varying RRRs and MAFs, assuming a total of 300 case-parent triads for the girls-only analyses in (A) and 600 case-parent triads in (B,C) (assuming an equal distribution of boys and girls). We used a nominal significance level of 0.05 throughout.
Table 3. Synopsis of the genes in which associations were identified in this study.
(A) Single-Marker and Haplotype Analyses Combining Estimates from Boys and Girls
Figure 2 displays the results of the single-marker and haplotype analyses in the Asian and European samples for each X-inactivation model.
(i) Asian Sample
There were no associations in the single-marker analyses of isolated CPO or isolated CL/P in the model without X-inactivation. By contrast, three SNPs stood out in the model with X-inactivation: rs5935623 in CPO; and both rs12690207 and rs2405829 in CL/P. rs5935623 is located in the gene for “EGF-like domain multiple 6” (EGFL6 at Xp22.2), rs12690207 is located in “Fibroblast growth factor 13” (FGF13 at Xq26.3-q27.1), and rs2405829 is located in “Dystrophin” (DMD at Xp21.2-p21.1).
In haplotype analyses, haplotype rs6527993-rs4969754-rs5990881-rs4969757 was associated with CL/P in the model without X-inactivation. Although the q-value for this haplotype was 0.008 (Table 2), none of the SNPs is located in or near a gene within 20 kb. In the model with X-inactivation, we found associations with two haplotypes: (i) rs4366220-rs5928243-rs988530-rs5928245 in CPO and (ii) rs5979748-rs2162-rs5935409-rs5935412 in CL/P. The first SNP in haplotype (i) lies ~9.7 kb from DMD, and the first SNP in haplotype (ii) lies ~5.3 kb from “Phosphoribosyl pyrophosphate synthetase 2” (PRPS2 at Xp22.2).
(ii) European Sample
In the model without X-inactivation, we found associations with rs4559361, rs5933262, and rs4128285 in CPO (all three SNPs had the same q-value of 0.155; Figure 2, Table 2). rs4559361 lies ~4.4 kb from the gene for “SPANX family, member N4” (SPANXN4 at Xq27.3), rs5933262 is located in “Ubiquitin specific peptidase 26” (USP26 at Xq26.2), and rs4128285 is not located near any gene within 20 kb. In the model with X-inactivation, rs6529005 was associated with CPO. This SNP is not located near any gene within 20 kb.
In haplotype analyses, rs2024759-rs6609469-rs743046-rs2498183 was associated with CPO in the model without X-inactivation. The last SNP in this haplotype is located ~24.2 kb from “Zinc finger protein 157” (ZNF157 at Xp11.3). In the model with X-inactivation, four haplotypes were associated with CPO, but only rs1139916-rs2266856-rs5925077-rs5925079 (q = 0.169) has SNPs located in a specific gene—in “Gamma-aminobutyric acid type A receptor epsilon subunit” (GABRE at Xq28).
(B) Single-Marker and Haplotype Analyses Restricted to Girls Only
Figure 3 displays the results of the single-marker and haplotype analyses restricted to the Asian and European girls only.
(i) Asian Sample
In single-marker analyses, we identified associations with the same three SNPs (rs5935623 in CPO, and both rs12690207 and rs2405829 in CL/P) as in the single-marker analyses based on the model with X-inactivation (Figure 2). In the haplotype analyses, two of the haplotypes that were associated in the analyses based on the model with X-inactivation (Figure 2) were also identified here and in the same cleft categories; notably, rs4366220-rs5928243-rs988530-rs5928245 in CPO and rs5979748-rs2162-rs5935409-rs5935412 in CL/P. There were two additional haplotypes associated with CPO: (i) rs982598-rs12843271-rs4826609-rs5961136 and (ii) rs6639176-rs5978529-rs1874039-rs5978530. The last SNP in haplotype (i) is located in “Inter-alpha-trypsin inhibitor heavy chain family member 6” (ITIH6 at Xp11.22), while the SNPs in haplotype (ii) are located in “FERM and PDZ domain containing 4” (FRMPD4 at Xp22.2).
(ii) European Sample
In single-marker analyses, rs6529005 (q = 0.161) was associated with CPO, which was also identified in the single-marker analyses in the same cleft category based on the model with X-inactivation (Figure 2). There were no associations with CL/P. In the haplotype analyses, the same SNP identified in the single-marker analysis above, rs6529005, appeared to be driving the association here between rs2869525-rs6529005-rs2180212-rs5954651 and CPO. One more haplotype was associated with CPO: rs10521639-rs12115935-rs16979405-rs12009996. Both of these haplotypes were already identified in the haplotype analyses based on the model with X-inactivation (Figure 2), but none of the SNPs in these haplotypes is located near any known gene within 20 kb.
Power Simulations
Haplin includes a complete framework for power simulations (Gjerdevik et al., 2017). Figure 4 displays the a priori power calculations for a single SNP to detect PoO effects on the X chromosome. All calculations were based on 1,000 simulated datasets with a 0.05 nominal significance level (note that a baseline risk of one was used throughout). The power of a single SNP to detect PoO effects on the X chromosome depends on several factors, including the minor allele frequency (MAF), effect size, sample size and family design. The power of the model focusing on girls only is similar to that of the model with X-inactivation in females (Figures 4A,C), which is consistent with the fact that those models provide similar estimates of RRR (as also mentioned in the Materials and Methods section). Overall, the power is sufficient for RRRs > 2, with the range of MAFs and sample sizes presented in our analyses (RRRs ≥ 2.5 for the smallest sample sizes).
Discussion
Our analyses detected possible PoO effects with several SNPs and haplotypes on the X chromosome, some of which were located in or close to genes (summarized in Table 3). Perhaps the most prominent gene on the list is “Dystrophin” (DMD), the second largest gene in humans according to the size of the transcript and protein product (Richards and Hawley, 2011). Dystrophin forms part of the dystrophin-glycoprotein complex whose function is to provide mechanical stability to the plasma membrane of striated muscle cells (Michele and Campbell, 2003; Gao and McNally, 2015). Several lines of evidence point to a link between muscular dystrophy and OFCs. For example, both congenital muscular dystrophy and CL/P have been observed in the rare autosomal-recessive Walker-Warburg syndrome (Moore et al., 1988; Dobyns et al., 1989; Pratap et al., 2007; Vajsar et al., 2008). More recently, we and others have reported associations between OFCs and variants in DMD in several chromosome-wide studies of X-linked SNPs (Patel et al., 2013; Fonseca et al., 2015; Wise et al., 2016; Skare et al., 2017). Even though only the Fonseca et al. (2015) study was based on a different GWAS dataset and study population compared to the other three studies, the repeated identification of DMD nevertheless constitutes an exciting new direction to pursue in future OFCs research.
Besides DMD, we identified PoO effects with several other genes, two of which—EGFL6 and FGF13—offer biologically plausible links to OFCs. EGFL6 belongs to the EGF-repeat containing superfamily of genes known to be involved in cell cycle regulation and proliferation (Buchner et al., 2000a). Expression analyses of murine Egf16 show strong expression in the mesenchymal components of both the hard and soft palate from embryonic stage E12.5 (Buchner et al., 2000b). This spatio-temporal pattern of gene expression supports a role for EGFL6 in the development of the orofacial complex, but as with DMD, the associations with EGFL6 will need to be validated in other OFCs samples.
FGF13 belongs to the FGF family of signaling molecules known to play key roles in embryonic development (Nie et al., 2006; Rahimov et al., 2012; Du et al., 2016). Several FGF ligands and receptors are expressed in the developing facial primordia, and mutations in specific FGF receptors have been reported to cause clefting syndromes (Rahimov et al., 2012). For example, Kallmann syndrome is caused by mutations in FGFR1 (Dodé et al., 2003) whereas Apert syndrome is caused by mutations in FGFR2 (Wilkie et al., 1995). FGF13 is thus a promising gene to pursue in other independent cleft samples.
As opposed to the above genes, there are no obvious connections between OFCs and the remaining genes in Table 3. It is nonetheless noteworthy that both PRPS2 and FRMPD4 are located in chromosomal region Xp22.2, which is also the location of EGFL6 whose link to OFCs was discussed above. Figure 5 is a collage of screenshots showing all the genes in chromosomal region Xp22.2 as they appear in the genome browser of The Ensembl Project (Kersey et al., 2016; Yates et al., 2016). This chromosomal region contains two genes that have been linked to clefting syndromes. The first is “Oral-facial-digital syndrome 1” (OFD1), a gene in which we had previously identified associations with OFCs in a candidate-gene analysis of X-linked markers (Jugessur et al., 2012a). Mutations in OFD1 underlie the X-linked dominant oral-facial digital syndrome type 1, which is characterized by malformations of the face, oral cavity and digits, and lethality in most affected males (Ferrante et al., 2001). The orofacial abnormalities in Oral-facial-digital syndrome 1 include median cleft lip, clefts of the alveolar ridge, and cleft palate. The second gene in chromosomal region Xp22.2 region is “Midline 1” (MID1). Mutations in this gene cause the X-linked Opitz GBBB syndrome, a congenital midline malformation syndrome featuring clefting of the lip and palate as part of the overall clinical picture (Quaderi et al., 1997). Furthermore, associations between isolated CL/P and specific haplotypes in MID1 have been reported in an Italian population (Scapoli et al., 2008).
Figure 5. Genes located in chromosomal region Xp22.2. This collage of screenshots of chromosomal region Xp22.2 was generated using the Ensembl genome browser (http://www.ensembl.org/) (Kersey et al., 2016; Yates et al., 2016). The three genes identified in this study are shown in boxes, in addition to the “Oral-facial-digital syndrome 1” (OFD1) gene, in which we had identified associations in our previous study on candidate genes on the X chromosome (Jugessur et al., 2012a), and the “Midline 1” (MID1) gene.
As seen in Figure 5, MID1, FRMPD4, PRPS2, EGFL6, and OFD1 all lie relatively close to one another on chromosome Xp22.2. Specifically, MID1 is ~1.3 Mb from FRMPD4, FRMPD4 is ~66.8 kb from PRPS2, PRPS2 is ~745.3 kb from EGFL6, and EGFL6 is ~81.8 kb from OFD1. Aside from the potential contributions of the genes themselves, their shared location on Xp22.2 suggests that there might be a specific locus (or loci) in that region driving the associations detected by our PoO analyses—a hypothesis worth pursuing in other OFCs samples. In addition, the apparent excess of genes on the X chromosome that are associated with sexually dimorphic traits (Saifi and Chandra, 1999; Vawter et al., 2004; Berletch et al., 2011) makes genes that escape X-inactivation (Yang et al., 2011; Deng et al., 2014; Balaton and Brown, 2016) particularly attractive candidates in explaining some of the observed male-female differences in development and physiology (Ober et al., 2008; Berletch et al., 2011). Such a scenario would be captured by our model assuming no X-inactivation (Figure 1B). Since both CL/P and CPO consistently exhibit a skewed sex ratio in prevalence, it may be worth exploring this hypothesis in future studies of PoO effects in OFCs. In addition, scenarios involving variable escape from X-chromosome inactivation (Peeters et al., 2014) are also amenable to testing using the X-inactivation models in Haplin.
Identifying associations between isolated OFCs and genes for clefting syndromes is not surprising, as it has long been recognized that Mendelian forms of clefting that phenocopy isolated clefts may provide important clues in identifying genes involved in isolated forms of clefts (Jugessur and Murray, 2005). One of the best examples to date is Van der Woude syndrome (VWS), in which hypodontia and lip pits are the only additional features distinguishing patients with VWS from those with isolated clefts. Mutations in interferon regulatory factor 6 (IRF6) were reported to cause VWS (Kondo et al., 2002), and, subsequently, we and others reported associations between IRF6 variants and isolated clefts (Zucchero et al., 2004; Jugessur et al., 2008, 2009b; Beaty et al., 2013; Leslie et al., 2016b; Moreno Uribe et al., 2017). Our current identification of several genes associated with clefting syndromes, while studying isolated OFCs per se, supports the notion that there may be a wider spectrum of subclinical features beside the overt cleft itself, blurring the distinction between isolated and non-isolated clefts.
Indeed, a handful of subclinical features have been characterized in patients who were initially classified as “isolated” cleft cases, as well as in their immediate family members (Weinberg et al., 2006; Jugessur et al., 2009a; Marazita, 2012). The inclusion of one such subclinical feature—breaks in the orbicularis oris muscle in the mouth—revealed a pattern of familial segregation more akin to that of a Mendelian trait (Marazita, 2007; Neiswanger et al., 2007). This further supports the existence of a wider array of subphenotypes that are associated with the overt cleft, and the identification of genes underlying syndromic clefts while analyzing isolated OFCs samples, like we do in this paper, is consistent with previous descriptions of an “extended cleft phenotype” (Weinberg et al., 2006; Marazita, 2007, 2012; Jugessur et al., 2009a; Rahimov et al., 2012).
Another important factor to consider in genetic association analyses is whether associations identified in one ethnicity are generalizable to other ethnicities. DMD was the only gene in which associations were identified in both CPO and CL/P, in the model with X-inactivation, in both single-marker and haplotype analyses, and in Europeans only. The remaining associations were highly specific for ethnicity, cleft subtype and X-inactivation model, with none of the associations in the Asian sample overlapping with those in the European sample. For example, the associations with EGFL6, FGF13, DMD, and PRPS2 were only detected in the Asian sample and in the model with X-inactivation, whereas the associations with SPANXN4, ZNF157, and GABRE were only observed in the European sample. Furthermore, associations with SPANXN4 and ZNF157 were only detected in the model assuming no X-inactivation, whereas associations with GABRE were only detected in the model with X-inactivation. These findings challenge the common assumption that the remarkable phenotypic consistency of OFCs across different ethnicities is because at least a subset of the causal variants is shared across different ethnicities. The apparent lack of overlap between Asians and Europeans was also observed in our recent chromosome-wide analysis of X-linked variants in the same GWAS dataset as here (Skare et al., 2017). Collectively, these findings underscore the importance of examining different ethnicities, X-inactivation models and cleft subtypes, even though this means adding more multiplicity to the analyses.
We observed a substantial overlap in the associations identified in the analyses based on the model with X-inactivation (Figure 2) and those identified in the analyses restricted to girls only (Figure 3). This was as expected, for the reasons provided in the Materials and Methods section. Regarding statistical power, our simulations indicated that there would be sufficient power for single-SNP analyses based on relative risk ratios (RRR) ≥ 2 (RRR ≥ 2.5 for the smaller sample of CPO triads) (Figure 4). It could be argued that the assumption of a RRR ≥ 2 or RRR ≥ 2.5 is optimistic in the context of OFCs, or most other complex traits for that matter. Our current analyses of X-linked PoO effects are novel and should be regarded as exploratory at this stage. With continuously accruing sample sizes of OFCs through multiple international collaborations, we envisage our new methods to gain wider currency in future analyses of X-linked PoO-effects in the larger pooled datasets.
This study benefited from being based on the largest collection of case-parent triads of OFCs to date, which enabled two major isolated cleft subtypes—CPO and CL/P—to be analyzed separately and without biases due to population stratification. In addition, having access to genotype data on two major ethnicities (Asian and European) enabled a more in-depth exploration of ethnicity- and sex-specific differences for each of the identified associations. The new methods presented here enabled an investigation of the hitherto untested possibility of PoO effects on the X chromosome, through the use of models with or without the assumption of X-inactivation in females. Limitations of the study include a lack of control triads to verify potential population transmission ratio distortions that could affect the PoO estimates. This could be the case if a genetic variant can cause early abortions that are not even registered as births. However, since the PoO RRR estimates are a ratio of two transmission distortions (from the mother and from the father), distortions unrelated to clefting, such as early abortions, would at least need to exhibit a parent-of-origin dependent pattern to substantially bias the RRR value. Furthermore, the lack of a comparable replication cohort precluded a formal validation of the current findings, which currently makes it difficult to interpret their relevance to clinical research.
To summarize, we developed new methods for a robust investigation of PoO-effects on the entire X chromosome and present the first results of such a screening in the largest available collection of OFCs triads to date. Our analyses identified associations with several genes, in particular with DMD, FGF13, and EFG16 that offer biologically plausible links to OFCs. The additional identification of associations with other genes on chromosomal region Xp22.2 deserves further scrutiny because several genes for clefting syndromes appear to cluster in this region. Except for DMD, in which associations were identified in both CPO and CL/P, the rest of the associations were highly specific for ethnicity, cleft subtype and X-inactivation model, highlighting the importance of performing such subanalyses, despite the need to adjust for additional multiple testing. More generally, the new methodology presented here is easily adaptable to the study of X-linked PoO effects in other outcomes that use a family-based study design.
Author Contributions
Conception of the work, study design, data analysis and interpretation, manuscript preparation and final manuscript approval: AJ, ØS, HG, RL, ØH, MG, and JR. Statistical modeling and software design: HG, ØS, MG, JR, and ØH. Funding acquisition: AJ, HG, and RL.
Funding
This research was supported by the Bergen Medical Research Foundation (BMFS) grant 807191, by the Research Council of Norway (RCN) through its Centers of Excellence funding scheme, grant 262700, and by the Biobank Norway II grant 245464/F50 from the RCN.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are grateful to all the participating families, medical staff and field workers at all the 13 recruitment sites for helping to make this study possible.
References
Al-Mubarak, B., Abouelhoda, M., Omar, A., AlDhalaan, H., Aldosari, M., Nester, M., et al. (2017). Whole exome sequencing reveals inherited and de novo variants in autism spectrum disorder: a trio study from Saudi families. Sci. Rep. 7:5679. doi: 10.1038/s41598-017-06033-1
Balaton, B. P., and Brown, C. J. (2016). Escape artists of the X chromosome. Trends Genet. 32, 348–359. doi: 10.1016/j.tig.2016.03.007
Beaty, T. H., Marazita, M. L., and Leslie, E. J. (2016). Genetic factors influencing risk to orofacial clefts: today's challenges and tomorrow's opportunities. F1000Res 5:2800. doi: 10.12688/f1000research.9503.1
Beaty, T. H., Murray, J. C., Marazita, M. L., Munger, R. G., Ruczinski, I., Hetmanski, J. B., et al. (2010). A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4. Nat. Genet. 42, 525–529. doi: 10.1038/ng.580
Beaty, T. H., Taub, M. A., Scott, A. F., Murray, J. C., Marazita, M. L., Schwender, H., et al. (2013). Confirming genes influencing risk to cleft lip with/without cleft palate in a case-parent trio study. Hum. Genet. 132, 771–781. doi: 10.1007/s00439-013-1283-6
Berg, E., Haaland, Ø. A., Feragen, K. B., Filip, C., Vindenes, H. A., Moster, D., et al. (2016a). Health status among adults born with an oral cleft in Norway. JAMA Pediatr. 170, 1063–1070. doi: 10.1001/jamapediatrics.2016.1925
Berg, E., Sivertsen, Å., Ariansen, A. M., Filip, C., Vindenes, H. A., Feragen, K. B., et al. (2016b). Socio-economic status and reproduction among adults born with an oral cleft: a population-based cohort study in Norway. PLoS ONE 11:e0162196. doi: 10.1371/journal.pone.0162196
Berletch, J. B., Yang, F., Xu, J., Carrel, L., and Disteche, C. M. (2011). Genes that escape from X inactivation. Hum. Genet. 130, 237–245. doi: 10.1007/s00439-011-1011-z
Bille, C., Winther, J. F., Bautz, A., Murray, J. C., Olsen, J., and Christensen, K. (2005). Cancer risk in persons with oral cleft–a population-based study of 8,093 cases. Am. J. Epidemiol. 161, 1047–1055. doi: 10.1093/aje/kwi132
Birnbaum, S., Ludwig, K. U., Reutter, H., Herms, S., Steffens, M., Rubini, M., et al. (2009). Key susceptibility locus for nonsyndromic cleft lip with or without cleft palate on chromosome 8q24. Nat. Genet. 41, 473–477. doi: 10.1038/ng.333
Buchner, G., Broccoli, V., Bulfone, A., Orfanelli, U., Gattuso, C., Ballabio, A., et al. (2000a). MAEG, an EGF-repeat containing gene, is a new marker associated with dermatome specification and morphogenesis of its derivatives. Mech. Dev. 98, 179–182. doi: 10.1016/S0925-4773(00)00462-7
Buchner, G., Orfanelli, U., Quaderi, N., Bassi, M. T., Andolfi, G., Ballabio, A., et al. (2000b). Identification of a new EGF-repeat-containing gene from human Xp22: a candidate for developmental disorders. Genomics 65, 16–23. doi: 10.1006/geno.2000.6146
Camargo, M., Rivera, D., Moreno, L., Lidral, A. C., Harper, U., Jones, M., et al. (2012). GWAS reveals new recessive loci associated with non-syndromic facial clefting. Eur. J. Med. Genet. 55, 510–514. doi: 10.1016/j.ejmg.2012.06.005
Christensen, K., Juel, K., Herskind, A. M., and Murray, J. C. (2004). Long term follow up study of survival associated with cleft lip and palate at birth. BMJ 328:1405. doi: 10.1136/bmj.38106.559120.7C
Connolly, S., and Heron, E. A. (2015). Review of statistical methodologies for the detection of parent-of-origin effects in family trio genome-wide association data with binary disease traits. Brief. Bioinformatics 16, 429–448. doi: 10.1093/bib/bbu017
Cordell, H. J. (2009). Detecting gene-gene interactions that underlie human diseases. Nat. Rev. Genet. 10, 392–404. doi: 10.1038/nrg2579
Cunningham, J. T., Moreno, M. V., Lodi, A., Ronen, S. M., and Ruggero, D. (2014). Protein and nucleotide biosynthesis are coupled by a single rate-limiting enzyme, PRPS2, to drive cancer. Cell 157, 1088–1103. doi: 10.1016/j.cell.2014.03.052
Deng, X., Berletch, J. B., Nguyen, D. K., and Disteche, C. M. (2014). X chromosome regulation: diverse patterns in development, tissues and disease. Nat. Rev. Genet. 15, 367–378. doi: 10.1038/nrg3687
Derry, J. M., Jess, U., and Francke, U. (1995). Cloning and characterization of a novel zinc finger gene in Xp11.2. Genomics 30, 361–365. doi: 10.1006/geno.1995.0031
Dixon, M. J., Marazita, M. L., Beaty, T. H., and Murray, J. C. (2011). Cleft lip and palate: understanding genetic and environmental influences. Nat. Rev. Genet. 12, 167–178. doi: 10.1038/nrg2933
Dobyns, W. B., Pagon, R. A., Armstrong, D., Curry, C. J., Greenberg, F., Grix, A., et al. (1989). Diagnostic criteria for Walker-Warburg syndrome. Am. J. Med. Genet. 32, 195–210. doi: 10.1002/ajmg.1320320213
Dodé, C., Levilliers, J., Dupont, J. M., De Paepe, A., Le Dû, N., Soussi-Yanicostas, N., et al. (2003). Loss-of-function mutations in FGFR1 cause autosomal dominant Kallmann syndrome. Nat. Genet. 33, 463–465. doi: 10.1038/ng1122
Du, W., Prochazka, J., Prochazkova, M., and Klein, O. D. (2016). Expression of FGFs during early mouse tongue development. Gene Expr. Patterns 20, 81–87. doi: 10.1016/j.gep.2015.12.003
Fatemi, S. H., Folsom, T. D., Rooney, R. J., and Thuras, P. D. (2013). Expression of GABAA alpha2-, beta1- and epsilon-receptors are altered significantly in the lateral cerebellum of subjects with schizophrenia, major depression and bipolar disorder. Transl. Psychiatry 3:e303. doi: 10.1038/tp.2013.64
Ferguson-Smith, A. C. (2011). Genomic imprinting: the emergence of an epigenetic paradigm. Nat. Rev. Genet. 12, 565–575. doi: 10.1038/nrg3032
Ferrante, M. I., Giorgio, G., Feather, S. A., Bulfone, A., Wright, V., Ghiani, M., et al. (2001). Identification of the gene for oral-facial-digital type I syndrome. Am. J. Hum. Genet. 68, 569–576. doi: 10.1086/318802
Fonseca, R. F., de Carvalho, F. M., Poletta, F. A., Montaner, D., Dopazo, J., Mereb, J. C., et al. (2015). Family-based genome-wide association study in Patagonia confirms the association of the DMD locus and cleft lip and palate. Eur. J. Oral Sci. 123, 381–384. doi: 10.1111/eos.12212
Gao, Q. Q., and McNally, E. M. (2015). The dystrophin complex: structure, function, and implications for therapy. Compr. Physiol. 5, 1223–1239. doi: 10.1002/cphy.c140048
Garg, P., Ludwig, K. U., Böhmer, A. C., Rubini, M., Steegers-Theunissen, R., Mossey, P. A., et al. (2014). Genome-wide analysis of parent-of-origin effects in non-syndromic orofacial clefts. Eur. J. Hum. Genet. 22, 822–830. doi: 10.1038/ejhg.2013.235
Gjerdevik, M., Haaland, O. A., Romanowska, J., Lie, R. T., Jugessur, A., and Gjessing, H. K. (2017). Parent-of-origin-environment interactions in case-parent triads with or without independent controls. Ann. Hum. Genet. doi: 10.1111/ahg.12224. [Epub ahead of print].
Gjessing, H. K., and Lie, R. T. (2006). Case-parent triads: estimating single- and double-dose effects of fetal and maternal disease gene haplotypes. Ann. Hum. Genet.70(Pt 3), 382–396. doi: 10.1111/j.1529-8817.2005.00218.x
Grant, S. F., Wang, K., Zhang, H., Glaberson, W., Annaiah, K., Kim, C. E., et al. (2009). A genome-wide association study identifies a locus for nonsyndromic cleft lip with or without cleft palate on 8q24. J. Pediatr. 155, 909–913. doi: 10.1016/j.jpeds.2009.06.020
Grosen, D., Chevrier, C., Skytthe, A., Bille, C., Mølsted, K., Sivertsen, A., et al. (2010). A cohort study of recurrence patterns among more than 54,000 relatives of oral cleft cases in Denmark: support for the multifactorial threshold model of inheritance. J. Med. Genet. 47, 162–168. doi: 10.1136/jmg.2009.069385
Guilmatre, A., and Sharp, A. J. (2012). Parent of origin effects. Clin. Genet. 81, 201–209. doi: 10.1111/j.1399-0004.2011.01790.x
Haaland, Ø. A., Jugessur, A., Gjerdevik, M., Romanowska, J., Shi, M., Beaty, T. H., et al. (2017). Genome-wide analysis of parent-of-origin interaction effects with environmental exposure (PoOxE): an application to European and Asian cleft palate trios. PLoS ONE 12:e0184358. doi: 10.1371/journal.pone.0184358
Hu, H., Haas, S. A., Chelly, J., Van Esch, H., Raynaud, M., de Brouwer, A. P., et al. (2016). X-exome sequencing of 405 unresolved families identifies seven novel intellectual disability genes. Mol. Psychiatry 21, 133–148. doi: 10.1038/mp.2014.193
Jugessur, A., Farlie, P. G., and Kilpatrick, N. (2009a). The genetics of isolated orofacial clefts: from genotypes to subphenotypes. Oral Dis. 15, 437–453. doi: 10.1111/j.1601-0825.2009.01577.x
Jugessur, A., and Murray, J. C. (2005). Orofacial clefting: recent insights into a complex trait. Curr. Opin. Genet. Dev. 15, 270–278. doi: 10.1016/j.gde.2005.03.003
Jugessur, A., Rahimov, F., Lie, R. T., Wilcox, A. J., Gjessing, H. K., Nilsen, R. M., et al. (2008). Genetic variants in IRF6 and the risk of facial clefts: single-marker and haplotype-based analyses in a population-based case-control study of facial clefts in Norway. Genet. Epidemiol. 32, 413–424. doi: 10.1002/gepi.20314
Jugessur, A., Shi, M., Gjessing, H. K., Lie, R. T., Wilcox, A. J., Weinberg, C. R., et al. (2009b). Genetic determinants of facial clefting: analysis of 357 candidate genes using two national cleft studies from Scandinavia. PLoS ONE 4:e5385. doi: 10.1371/journal.pone.0005385
Jugessur, A., Skare, Ø., Harris, J. R., Lie, R. T., and Gjessing, H. K. (2012b). Using offspring-parent triads to study complex traits: a tutorial based on orofacial clefts. Norweg. J. Epidemiol. 21, 251–267. doi: 10.5324/nje.v21i2.1503
Jugessur, A., Skare, Ø., Lie, R. T., Wilcox, A. J., Christensen, K., Christiansen, L., et al. (2012a). X-linked genes and risk of orofacial clefts: evidence from two population-based studies in Scandinavia. PLoS ONE 7:e39240. doi: 10.1371/journal.pone.0039240
Kersey, P. J., Allen, J. E., Armean, I., Boddu, S., Bolt, B. J., Carvalho-Silva, D., et al. (2016). Ensembl Genomes 2016: more genomes, more complexity. Nucleic Acids Res. 44, D574–D580. doi: 10.1093/nar/gkv1209
Kondo, S., Schutte, B. C., Richardson, R. J., Bjork, B. C., Knight, A. S., Watanabe, Y., et al. (2002). Mutations in IRF6 cause Van der Woude and popliteal pterygium syndromes. Nat. Genet. 32, 285–289. doi: 10.1038/ng985
Kouprina, N., Mullokandov, M., Rogozin, I. B., Collins, N. K., Solomon, G., Otstot, J., et al. (2004). The SPANX gene family of cancer/testis-specific antigens: rapid evolution and amplification in African great apes and hominids. Proc. Natl. Acad. Sci. U.S.A. 101, 3077–3082. doi: 10.1073/pnas.0308532100
Kousa, Y. A., and Schutte, B. C. (2016). Toward an orofacial gene regulatory network. Dev. Dyn. 245, 220–232. doi: 10.1002/dvdy.24341
Lawson, H. A., Cheverud, J. M., and Wolf, J. B. (2013). Genomic imprinting and parent-of-origin effects on complex traits. Nat. Rev. Genet. 14, 609–617. doi: 10.1038/nrg3543
Leslie, E. J., Carlson, J. C., Shaffer, J. R., Feingold, E., Wehby, G., Laurie, C. A., et al. (2016a). A multi-ethnic genome-wide association study identifies novel loci for non-syndromic cleft lip with or without cleft palate on 2p24.2, 17q23 and 19q13. Hum. Mol. Genet. 25, 2862–2872. doi: 10.1093/hmg/ddw104
Leslie, E. J., Koboldt, D. C., Kang, C. J., Ma, L., Hecht, J. T., Wehby, G. L., et al. (2016b). IRF6 mutation screening in non-syndromic orofacial clefting: analysis of 1521 families. Clin. Genet. 90, 28–34. doi: 10.1111/cge.12675
Leslie, E. J., Liu, H., Carlson, J. C., Shaffer, J. R., Feingold, E., Wehby, G., et al. (2016c). A Genome-wide association study of nonsyndromic cleft palate identifies an etiologic missense variant in GRHL3. Am. J. Hum. Genet. 98, 744–754. doi: 10.1016/j.ajhg.2016.02.014
Ludwig, K. U., Mangold, E., Herms, S., Nowak, S., Reutter, H., Paul, A., et al. (2012). Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci. Nat. Genet. 44, 968–971. doi: 10.1038/ng.2360
Mangold, E., Ludwig, K. U., Birnbaum, S., Baluardo, C., Ferrian, M., Herms, S., et al. (2010). Genome-wide association study identifies two susceptibility loci for nonsyndromic cleft lip with or without cleft palate. Nat. Genet. 42, 24–26. doi: 10.1038/ng.506
Marazita, M. L. (2007). Subclinical features in non-syndromic cleft lip with or without cleft palate (CL/P): review of the evidence that subepithelial orbicularis oris muscle defects are part of an expanded phenotype for CL/P. Orthod. Craniofac. Res. 10, 82–87. doi: 10.1111/j.1601-6343.2007.00386.x
Marazita, M. L. (2012). The evolution of human genetic studies of cleft lip and cleft palate. Annu. Rev. Genomics Hum. Genet. 13, 263–283. doi: 10.1146/annurev-genom-090711-163729
Matosin, N., Green, M. J., Andrews, J. L., Newell, K. A., and Fernandez-Enright, F. (2016). Possibility of a sex-specific role for a genetic variant in FRMPD4 in schizophrenia, but not cognitive function. Neuroreport 27, 33–38. doi: 10.1097/WNR.0000000000000491
Michele, D. E., and Campbell, K. P. (2003). Dystrophin-glycoprotein complex: post-translational processing and dystroglycan function. J. Biol. Chem. 278, 15457–15460. doi: 10.1074/jbc.R200031200
Moore, G., Ivens, A., Chambers, J., Bjornsson, A., Arnason, A., Jensson, O., et al. (1988). The application of molecular genetics to detection of craniofacial abnormality. Development 103(Suppl.), 233–239.
Moreno Uribe, L. M., Fomina, T., Munger, R. G., Romitti, P. A., Jenkins, M. M., Gjessing, H. K., et al. (2017). A population-based study of effects of genetic Loci on orofacial clefts. J. Dent. Res. 96, 1322–1329. doi: 10.1177/0022034517716914
Mossey, P. A., Little, J., Munger, R. G., Dixon, M. J., and Shaw, W. C. (2009). Cleft lip and palate. Lancet 374, 1773–1785. doi: 10.1016/S0140-6736(09)60695-4
Neiswanger, K., Weinberg, S. M., Rogers, C. R., Brandon, C. A., Cooper, M. E., Bardi, K. M., et al. (2007). Orbicularis oris muscle defects as an expanded phenotypic feature in nonsyndromic cleft lip with or without cleft palate. Am. J. Med. Genet. A 143A, 1143–1149. doi: 10.1002/ajmg.a.31760
Nie, X., Luukko, K., and Kettunen, P. (2006). FGF signalling in craniofacial development and developmental disorders. Oral Dis. 12, 102–111. doi: 10.1111/j.1601-0825.2005.01176.x
Ober, C., Loisel, D. A., and Gilad, Y. (2008). Sex-specific genetic architecture of human disease. Nat. Rev. Genet. 9, 911–922. doi: 10.1038/nrg2415
Patel, P. J., Beaty, T. H., Ruczinski, I., Murray, J. C., Marazita, M. L., Munger, R. G., et al. (2013). X-linked markers in the Duchenne muscular dystrophy gene associated with oral clefts. Eur. J. Oral Sci. 121, 63–68. doi: 10.1111/eos.12025
Peeters, S. B., Cotton, A. M., and Brown, C. J. (2014). Variable escape from X-chromosome inactivation: identifying factors that tip the scales towards expression. Bioessays 36, 746–756. doi: 10.1002/bies.201400032
Peters, J. (2014). The role of genomic imprinting in biology and disease: an expanding view. Nat. Rev. Genet. 15, 517–530. doi: 10.1038/nrg3766
Pratap, A., Agrawal, A., Tiwari, A., Lakshmi, R., and Rajbanshi, S. (2007). The walker-warburg syndrome with cleft lip and palate. Singapore Med. J. 48, e66–e67.
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Quaderi, N. A., Schweiger, S., Gaudenz, K., Franco, B., Rugarli, E. I., Berger, W., et al. (1997). Opitz G/BBB syndrome, a defect of midline development, is due to mutations in a new RING finger gene on Xp22. Nat. Genet. 17, 285–291. doi: 10.1038/ng1197-285
Rahimov, F., Jugessur, A., and Murray, J. C. (2012). Genetics of nonsyndromic orofacial clefts. Cleft Palate Craniofac. J. 49, 73–91. doi: 10.1597/10-178
R Development Core Team (2014). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Richards, J. E., and Hawley, R. S. (2011). The Human Genome, A User's Guide, 3rd Edn. London: Academic Press.
Saifi, G. M., and Chandra, H. S. (1999). An apparent excess of sex- and reproduction-related genes on the human X chromosome. Proc. Biol. Sci. 266, 203–209. doi: 10.1098/rspb.1999.0623
Scapoli, L., Martinelli, M., Arlotti, M., Palmieri, A., Masiero, E., Pezzetti, F., et al. (2008). Genes causing clefting syndromes as candidates for non-syndromic cleft lip with or without cleft palate: a family-based association study. Eur. J. Oral. Sci. 116, 507–511. doi: 10.1111/j.1600-0722.2008.00574.x
Shi, M., Murray, J. C., Marazita, M. L., Munger, R. G., Ruczinski, I., Hetmanski, J. B., et al. (2012). Genome wide study of maternal and parent-of-origin effects on the etiology of orofacial clefts. Am. J. Med. Genet. A 158A, 784–794. doi: 10.1002/ajmg.a.35257
Sivertsen, A., Wilcox, A. J., Skjaerven, R., Vindenes, H. A., Abyholm, F., Harville, E., et al. (2008). Familial risk of oral clefts by morphological type and severity: population based cohort study of first degree relatives. BMJ 336, 432–434. doi: 10.1136/bmj.39458.563611.AE
Skare, Ø., Gjessing, H. K., Gjerdevik, M., Haaland, Ø. A., Romanowska, J., Lie, R. T., et al. (2017). A new approach to chromosome-wide analysis of X-linked markers identifies new associations in Asian and European case-parent triads of orofacial clefts. PLoS ONE 12:e0183772. doi: 10.1371/journal.pone.0183772
Storey, J. D., and Tibshirani, R. (2003). Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. U.S.A. 100, 9440–9445. doi: 10.1073/pnas.1530509100
Stouffs, K., Lissens, W., Tournaye, H., Van Steirteghem, A., and Liebaers, I. (2005). Possible role of USP26 in patients with severely impaired spermatogenesis. Eur. J. Hum. Genet. 13, 336–340. doi: 10.1038/sj.ejhg.5201335
Thomas, D. (2010). Gene–environment-wide association studies: emerging approaches. Nat. Rev. Genet. 11, 259–272. doi: 10.1038/nrg2764
Trujillano, D., Bertoli-Avella, A. M., Kumar Kandaswamy, K., Weiss, M. E., Köster, J., Marais, A., et al. (2017). Clinical exome sequencing: results from 2819 samples reflecting 1000 families. Eur. J. Hum. Genet. 25, 176–182. doi: 10.1038/ejhg.2016.146
Vajsar, J., Baskin, B., Swoboda, K., Biggar, D. W., Schachter, H., and Ray, P. N. (2008). Walker-Warburg Syndrome with POMT1 mutations can be associated with cleft lip and cleft palate. Neuromuscul. Disord. 18, 675–677. doi: 10.1016/j.nmd.2008.05.014
Vawter, M. P., Evans, S., Choudary, P., Tomita, H., Meador-Woodruff, J., Molnar, M., et al. (2004). Gender-specific gene expression in post-mortem human brain: localization to sex chromosomes. Neuropsychopharmacology 29, 373–384. doi: 10.1038/sj.npp.1300337
Wehby, G. L., and Cassell, C. H. (2010). The impact of orofacial clefts on quality of life and healthcare use and costs. Oral Dis. 16, 3–10. doi: 10.1111/j.1601-0825.2009.01588.x
Wehby, G. L., Collet, B., Barron, S., Romitti, P. A., Ansley, T. N., and Speltz, M. (2014). Academic achievement of children and adolescents with oral clefts. Pediatrics 133, 785–792. doi: 10.1542/peds.2013-3072
Wei, W. H., Hemani, G., and Haley, C. S. (2014). Detecting epistasis in human complex traits. Nat. Rev. Genet. 15, 722–733. doi: 10.1038/nrg3747
Weinberg, C. R. (1999). Methods for detection of parent-of-origin effects in genetic studies of case-parents triads. Am. J. Hum. Genet. 65, 229–235. doi: 10.1086/302466
Weinberg, S. M., Neiswanger, K., Martin, R. A., Mooney, M. P., Kane, A. A., Wenger, S. L., et al. (2006). The Pittsburgh oral-facial cleft study: expanding the cleft phenotype. Background and justification. Cleft Palate Craniofac. J. 43, 7–20. doi: 10.1597/04-122R1.1
Wilkie, A. O., Slaney, S. F., Oldridge, M., Poole, M. D., Ashworth, G. J., Hockley, A. D., et al. (1995). Apert syndrome results from localized mutations of FGFR2 and is allelic with Crouzon syndrome. Nat. Genet. 9, 165–172. doi: 10.1038/ng0295-165
Wise, A. L., Gyi, L., and Manolio, T. A. (2013). eXclusion: toward integrating the X chromosome in genome-wide association analyses. Am. J. Hum. Genet. 92, 643–647. doi: 10.1016/j.ajhg.2013.03.017
Wise, A. S., Shi, M., and Weinberg, C. R. (2016). Family-based multi-SNP X chromosome analysis using parent information. Front. Genet. 7:20. doi: 10.3389/fgene.2016.00020
Wolf, Z. T., Brand, H. A., Shaffer, J. R., Leslie, E. J., Arzi, B., Willet, C. E., et al. (2015). Genome-wide association studies in dogs and humans identify ADAMTS20 as a risk variant for cleft lip and palate. PLoS Genet. 11:e1005059. doi: 10.1371/journal.pgen.1005059
Wosnitzer, M. S., Mielnik, A., Dabaja, A., Robinson, B., Schlegel, P. N., and Paduch, D. A. (2014). Ubiquitin specific protease 26 (USP26) expression analysis in human testicular and extragonadal tissues indicates diverse action of USP26 in cell differentiation and tumorigenesis. PLoS ONE 9:e98638. doi: 10.1371/journal.pone.0098638
Yang, C., Chapman, A. G., Kelsey, A. D., Minks, J., Cotton, A. M., and Brown, C. J. (2011). X-chromosome inactivation: molecular mechanisms from the human perspective. Hum. Genet. 130, 175–185. doi: 10.1007/s00439-011-0994-9
Yates, A., Akanni, W., Amode, M. R., Barrell, D., Billis, K., Carvalho-Silva, D., et al. (2016). Ensembl 2016. Nucleic Acids Res. 44, D710–D716. doi: 10.1093/nar/gkv1157
Zhang, W., Liu, T., Mi, Y. J., Yue, L. D., Wang, J. M., Liu, D. W., et al. (2015). Evidence from enzymatic and meta-analyses does not support a direct association between USP26 gene variants and male infertility. Andrology 3, 271–279. doi: 10.1111/andr.295
Zhu, J. L., Basso, O., Hasle, H., Winther, J. F., Olsen, J. H., and Olsen, J. (2002). Do parents of children with congenital malformations have a higher cancer risk? A nationwide study in Denmark. Br. J. Cancer 87, 524–528. doi: 10.1038/sj.bjc.6600488
Keywords: parent-of-origin, X chromosome, GWAS, case-parent triads, orofacial clefts, birth defects, genetic epidemiology, Haplin
Citation: Skare Ø, Lie RT, Haaland ØA, Gjerdevik M, Romanowska J, Gjessing HK and Jugessur A (2018) 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. Front. Genet. 9:25. doi: 10.3389/fgene.2018.00025
Received: 03 October 2017; Accepted: 22 January 2018;
Published: 22 February 2018.
Edited by:
Dana C. Crawford, Case Western Reserve University, United StatesReviewed by:
Marie-Hélène Roy-Gagnon, University of Ottawa, CanadaJing Dong, Baylor College of Medicine, United States
Copyright © 2018 Skare, Lie, Haaland, Gjerdevik, Romanowska, Gjessing and Jugessur. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Astanand Jugessur, YXN0YW5hbmQuanVnZXNzdXJAZmhpLm5v
†Joint first authors.
‡Joint senior authors.