PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(9): e0184358.
Published online 2017 September 12. doi:  10.1371/journal.pone.0184358
PMCID: PMC5595425

Genome-wide analysis of parent-of-origin interaction effects with environmental exposure (PoOxE): An application to European and Asian cleft palate trios

Øystein A. Haaland, Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Software, Visualization, Writing – original draft, Writing – review & editing,1,* Astanand Jugessur, Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing,1,2,3, Miriam Gjerdevik, Conceptualization, Formal analysis, Methodology, Software, Visualization, Writing – review & editing,1,2 Julia Romanowska, Formal analysis, Software, Writing – review & editing,1 Min Shi, Data curation, Writing – review & editing,4 Terri H. Beaty, Conceptualization, Funding acquisition, Resources, Writing – review & editing,5 Mary L. Marazita, Conceptualization, Funding acquisition, Resources, Writing – review & editing,6 Jeffrey C. Murray, Conceptualization, Funding acquisition, Resources, Writing – review & editing,7 Allen J. Wilcox, Funding acquisition, Resources, Writing – review & editing,8 Rolv T. Lie, Conceptualization, Formal analysis, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing,1,9, and Håkon K. Gjessing, Conceptualization, Formal analysis, Funding acquisition, Methodology, Software, Supervision, Visualization, Writing – review & editing1,3,
Andrew T DeWan, Editor

Abstract

Cleft palate only is a common birth defect with high heritability. Only a small fraction of this heritability is explained by the genetic variants identified so far, underscoring the need to investigate other disease mechanisms, such as gene-environment (GxE) interactions and parent-of-origin (PoO) effects. Furthermore, PoO effects may vary across exposure levels (PoOxE effects). Such variation is the focus of this study. We upgraded the R-package Haplin to enable direct tests of PoOxE effects at the genome-wide level. From a previous GWAS, we had genotypes for 550 case-parent trios, of mainly European and Asian ancestry, and data on three maternal exposures (smoking, alcohol, and vitamins). Data were analyzed for Europeans and Asians separately, and also for all ethnicities combined. To account for multiple testing, a false discovery rate method was used, where q-values were generated from the p-values. In the Europeans-only analyses, interactions with maternal smoking yielded the lowest q-values. Two SNPs in the ‘Interactor of little elongation complex ELL subunit 1’ (ICE1) gene had a q-value of 0.14, and five of the 20 most significant SNPs were in the ‘N-acetylated alpha-linked acidic dipeptidase-like 2’ (NAALADL2) gene. No evidence of PoOxE effects was found in the other analyses. The connections to ICE1 and NAALADL2 are novel and warrant further investigation. More generally, the new methodology presented here is easily applicable to other traits and exposures in which a family-based study design has been implemented.

Introduction

With a prevalence of 0.5 per 1000 live births, cleft palate only (CPO) is a common birth defect in humans [1, 2]. It is broadly categorized according to whether it occurs as an isolated defect or together with additional congenital anomalies. In this paper, we focus on isolated CPO.

The particularly high heritability and recurrence risk of orofacial clefts [38] have spurred long-standing efforts to identify genetic variants controlling risk to these common birth defects. However, as with most other complex traits, the genetic variants identified thus far explain only a small fraction of the total heritability and familial recurrence, underscoring the need to examine etiologic mechanisms beyond simple child effects alone. One alternative is to investigate the effect of a risk-allele or haplotype based on whether it is inherited from the mother or the father (i.e., parent-of-origin (PoO) effects). A difference in effect by parent of origin could occur, for example, with genes that are subject to genomic imprinting [9], which occurs when the allele from one parent is silenced but the allele from the other parent is expressed. This possibility is especially relevant for perinatal disorders because the mother defines the prenatal environment of the fetus.

Another popular approach is to explore the role of environmental factors, either independently or in combination with specific genetic variants (GxE effects). Although animal models have long demonstrated that environmental factors are important in clefting (reviewed in [10, 11]), the evidence from human studies is less conclusive. Among a wide array of environmental factors, maternal periconceptional smoking has been consistently associated with increased risk of clefting [1214]. Since most environmental factors are modifiable, identifying GxE effects may help to target genetically susceptible subgroups of the population. A third, yet unexplored approach is to study PoO effects in interaction with environmental exposures (PoOxE); i.e., whether PoO effects vary according to the exposure status of the fetus. With the notable exception of Wang et al. (2011) [15], who assessed differential imprinting across environmental exposures in childhood asthma, the literature on PoOxE effect estimation is sparse. To address this gap, we have developed a comprehensive and user-friendly methodology that is not restricted by assumptions pertaining to imprinting. The theoretical foundation for these new methods has been presented by Skare et al. (2012) [14] and Gjerdevik et al. (2017) [16], and the methods themselves are available in the R-package Haplin [17]. The mathematics behind the PoOxE analyses is outlined in Materials and methods.

This study is based on the case-parent trio study design, which is applicable to a wide range of etiologic scenarios pertinent to perinatal disorders [18]. We had GWAS data as well as information on periconceptional exposures from the mother (cigarette smoking, alcohol intake, vitamin use) and ethnicity (European, Asian, other) for the largest collection of CPO trios to date [19]. Our aim is to identify PoOxE effects in this data set.

Results

We conducted three sets of analyses: pooled analyses including all participants; analyses restricted to Europeans only; and analyses restricted to Asians only. The remaining ethnic groups in our data set were too small to justify separate analyses (Table 1). Given the phenotypic consistency in clefting across ethnicities, it is reasonable to assume that a proportion of the causal variants for clefting is shared across all ethnicities. Accordingly, we present the results of the pooled analyses first, followed by the Europeans-only and Asians-only analyses. The combination of three environmental exposures and the above subgroup analyses yielded a large amount of results. For simplicity, we chose to focus on the top 20 SNPs (sorted by observed p-value) from each analysis. Details about these SNPs, including relative risk ratios (RRRs), are provided in Table 2 and Fig 1, Table 3 and Fig 2, and Table 4 and Fig 3. The corresponding Manhattan plots are provided as supplementary online material (S1 to S3 Figs). Table 5 contains the full names of all the genes mentioned in Tables Tables22 to to44.

Fig 1
Pooled analyses of all ethnicities combined.
Fig 2
Analyses of the European sample.
Fig 3
Analyses of vitamin use in the Asian sample.
Table 1
Population distribution according to ethnicity and trio completeness.
Table 2
Top 20 SNPs sorted according to p-value in the pooled analysis.
Table 3
Top 20 SNPs sorted according to p-value in the Europeans-only analysis.
Table 4
Top 20 SNPs sorted according to p-value in the Asians-only analysis.
Table 5
Full gene names*.

To adjust for multiple testing, we used a false discovery rate method where q-values are calculated from observed p-values [20]. We used a q-value of 0.1 to assess statistical significance, which means that at least 90% of the significant SNPs are expected to be true positives. Across all analyses, several SNPs had q-values ranging from 0.1 to 0.5 (Tables (Tables22 to to4).4). This corresponds to a false discovery rate between 10% and 50%, implying that many of these SNPs are potentially associated with PoOxE effects. Fig 1 shows QQ-plots for the pooled analyses, comprising all ethnicities. All of the most significant SNPs are within the 95% confidence band at the upper right corner of the distribution. The lowest q-values were 0.8 for rs1116099 for maternal smoking, 0.5 for rs6092934 for maternal alcohol intake, and 0.5 for rs2830634 for maternal vitamin use (Table 2).

QQ-plots for the Europeans-only analyses are shown in Fig 2. The plot for smoking is particularly notable because all the top 12 SNPs had lower p-values than expected, even though most of them were located within the 95% confidence band. Specific p-values and q-values for these SNPs are provided in Table 3. All of these q-values were below 0.5 for the top 12 SNPs, but markedly higher for the remaining SNPs. Among these 12 SNPs, both rs2964447 and rs2964137 had a q-value of 0.14 (RRR = 0.09, 95% CI: 0.04–0.23). For alcohol intake and vitamin use, the top SNPs were rs6092934 (q = 0.8, RRR = 8.0, 95% CI: 3.2–19.8) and rs1400316 (q = 0.4, RRR = 10.1, 95% CI: 4.0–25.6), respectively.

The Asians-only analyses were uninformative due to the low number of trios in which the mother had smoked or consumed alcohol (Table 6). Consequently, tests for interaction had less power than the other analyses. For vitamin use, the QQ-plot did not deviate appreciably from the expected pattern (Fig 3). Table 4 shows the p-values and q-values for the top 20 SNPs. All the SNPs in the Asians-only analyses had q-values equal to one.

Table 6
Characteristics of maternal exposures according to ethnicity.

Several of the top 20 SNPs were the same across the three main analyses (pooled, Europeans-only, and Asians-only). The pooled and Europeans-only analyses had eight of the top SNPs in common for PoOxSmoke, three for PoOxAlcohol, and one for PoOxVitamin (Table 2). Similarly, the pooled and Asians-only analyses had three of the top SNPs in common for PoOxVitamin (Table 2). As several of the top 20 SNPs were located in the gene for ‘N-acetylated alpha-linked acidic dipeptidase-like 2’ (NAALADL2), we generated a regional association plot for rs4243412, which was the SNP in NAALADL2 with the lowest p-value in the Europeans-only analysis (Fig 4). We created a similar plot for rs2964137, which was the SNP with the lowest p-value in the pooled analysis (Fig 5). This SNP is located near the ‘Interactor of little elongation complex ELL subunit 1’ (ICE1) gene, and was also found among the top 20 SNPs in the Europeans-only analysis (Table 2).

Fig 4
Regional association plot for rs4243412 in NAALADL2.
Fig 5
Regional association plot for rs2964137 near ICE1.

Because PoO effects and maternal effects may be mutually confounded [21], we performed sensitivity analyses on the above-mentioned top 20 SNPs, and adjusted for potential maternal effects in each stratum of exposure. In these analyses, the RRRs were similar to those in Tables Tables22 to to4,4, and the Bonferroni corrected p-values for the interaction between maternal and environmental effects were all equal to 1.

Discussion

Our study used data from the largest collection of CPO trios to date [19] to investigate the hitherto untested possibility of interactions between PoO effects and maternal environmental exposures that have previously shown associations with clefts. We introduce new methodology that not only tests for PoOxE effects but also quantifies them as ratios of relative risks. All analyses were implemented in the R-package Haplin, which accommodates a wide range of etiologic scenarios suitable for family-based study designs. An example code for PoOxE analysis is provided in S1 Appendix.

Pooled analyses

For PoOxSmoke, all p-values were higher in the pooled analyses than in the Europeans-only analyses, suggesting a dilution of effects after pooling data. This reduction of the effect estimate in the pooled analyses may reflect heterogeneity of effect among the subgroups. The opposite was true for PoOxAlcohol, which may indicate a more consistent effect of alcohol across ethnicities. Regarding maternal smoking, multiple SNPs in NAALADL2 indicated the presence of PoOxSmoke effects. No genes or SNPs stood out in the PoOxVitamin analysis.

Europeans-only analyses

We found suggestive evidence of a PoOxSmoke effect for rs2964137 and rs2964447. Although neither SNP is located within any known gene, both lie near ICE1 and are only 2–15 kb from three copy-number variant regions (CNVRs). As in the pooled PoOxSmoke analysis, several top SNPs are located in NAALADL2. Previous analyses of data from genome rearrangements, GWAS, and gene-expression have linked this gene to various disorders, including mild retardation [22] and cancer [23]. We have not been able to find a connection between clefting and ‘Glucocorticoid induced transcript 1’ (GLCCI1), ‘Islet cell autoantigen 1’ (ICA1), or ‘Zinc finger and homeobox 2’ (ZHX2).

Regarding PoOxAlcohol effects, ‘Nitric oxide synthase 1’ (NOS1) and ‘Dipeptidyl-peptidase 6’ (DPP6) were among the most interesting genes. NOS1 acts as a physiological modulator of skeletal muscle function and DPP6 is involved in embryonic craniofacial development [24, 25]. Another member of the nitric oxide synthase family, NOS3, is involved in the folate pathway and has previously been linked to a higher risk of isolated CL/P in a non-Hispanic white population [26]. Furthermore, analysis of biopsies of soft palate muscle tissues from children with isolated clefts showed that NOS1 immunoreactivity in the muscle fibers was strongly influenced by the cleft itself [27].

In the PoOxVitamin analysis, three SNPs were located in the ‘Discs, large homolog 2’ (DLG2) gene on chromosome 11q14.1. One of these SNPs in DLG2, rs1400316, had the lowest q-value (0.4). Little has been reported about its role in clefting. Three other genes, ‘Guanylate cyclase activator 1C’ (GUCA1C), ‘TBC1 domain family, member 22A’ (TBC1D22A) and ‘Cytochrome P450, family 4, subfamily F, member 3’ (CYP4F3), each contain two of the top 20 SNPs from this analysis. Based on the literature, however, GUCA1C and TBC1D22A do not appear to have any connections to clefting. In contrast, CYP4F3 belongs to the cytochrome P450 gene family, which is known to be involved in the biotransformation of endobiotics and xenobiotics [28], and may be relevant for clefting. Still, the q-values for SNPs in CYP4F3 were 0.8 or higher.

Asians-only analyses

Compared with European women, Asian women generally consume little alcohol and tobacco [29, 30], which would be expected to be even less among those who are pregnant or planning to be pregnant. This was also observed in our data (Table 6). Even though a lack of observations was not a problem for the PoOxVitamin analyses, all the q-values were equal to one and there were no convincing associations overall for this ethnic group. Regarding ethnic specificity and generalizability, none of the top SNPs in the Asians-only analyses were among the top SNPs in the Europeans-only analyses (Tables (Tables33 and and4),4), which suggests ethnic-specific effects. Still, the lack of markers in common was somewhat unexpected, as GxE effects have previously been reported across the two ethnicities in the same sample population studied here [31]. However, that study used a different approach; the pooled sample was analyzed first and the top SNPs were verified to see whether the results were consistent across ethnicities. Additionally, the authors did not consider PoOxE.

Methodological considerations

The case-parent trio study design coupled with a large data set provided an excellent opportunity to explore PoOxE effects. The study design protects against false positives due to population substructure, because it aims at detecting asymmetries in allele transmission from parents to the affected child (proband), as opposed to considering only differences in allele frequencies at a population level. Still, if populations of different ethnicities react differently to a given exposure, such that there is a PoOxE effect in one population but not in the other, this effect may be muted or even go undetected in the combined population. It is therefore judicious to stratify analysis by ethnicity.

PoO effects may be seen when a gene associated with a given phenotype is also subjected to genomic imprinting [32, 33]. Through DNA methylation, the expression of a particular gene can be upregulated or downregulated depending on its parental origin [9, 34]. It is thus reasonable to assume that maternal environmental exposures capable of influencing methylation patterns might also influence the phenotype differently for maternally and paternally inherited alleles. Hence, it is conceivable that looking specifically for PoOxE effects rather than standard PoO or GxE effects alone might increase the chance of finding gene effects that are indicative of, for instance, genomic imprinting.

While PoOxE searches combine PoO searches with ordinary GxE searches in a natural way, there is a price to pay in the form of added complexity. Nevertheless, the total PoOxE effect at a locus with two alleles and a dichotomous environmental exposure can be measured as a single ratio of relative risks (RRR). We have

RRR=RRRPoO(1)RRRPoO(0)=RRmat(1)/RRpat(1)RRmat(0)/RRpat(0),
(1)

where RRmat(S) and RRpat(S) are as explained in Materials and Methods, and RRR is the ratio of PoO effects in the two strata. If RRR > 1, the interpretation is that the PoO effect RRmat(1)/RRpat(1) in stratum 1 is larger than the corresponding RRmat(0)/RRpat(0) in stratum 0. Note that this may come about in different ways. For example, consider an allele that increases the risk only when inherited from exposed mothers, so that RRmat(1) > 1. Because the other RRs are equal to 1, RRR would be larger than 1. Similarly, if the allele is protective when inherited from unexposed mothers but has no effect in other situations, RRmat(0) < 1, and again RRR > 1. One might also observe more complex patterns, such as an increased risk when the allele is inherited from the mother, where this effect is larger among the exposed than the unexposed; that is,

RRmat(1)  >  RRmat(0) and RRpat(1)  =  RRpat(0), 

and we would again have RRR > 1. The actual direction of the effect may depend on which allele and exposure group are chosen as reference, which is a general problem when assessing GxE in case-only designs.

While ordinary PoO analyses consider the ratio RRmat/RRpat for both strata combined, and ordinary GxE analyses consider RR(1)/RR(0) without accounting for parental origin, the full PoOxE RRR involves comparing four quantities—the effects of maternally and paternally derived alleles computed in both strata separately. Thus, a certain loss of power would be expected relative to the standard tests for PoO and GxE effects. This is indeed what we observe in the power simulations (Fig 6, right panel). We therefore decided not to include maternal genomic effects in the full GWAS analysis, since this is likely to further reduce power to detect PoO effects [21]. Instead, we performed sensitivity analyses to remove any positive confounding from maternal effects for the 20 most promising SNPs in each set of analyses (shown in Tables Tables22 to to4).4). It is not particularly likely that any of the genes involved in the sensitivity analyses would operate through maternal effects. Complex, but less likely scenarios where maternal effects cancel out PoO effects may be missed by this approach, however.

Fig 6
Simulation-based power curves.

As shown in Fig 6, PoOxE analyses will generally have lower power, given similar effect sizes, compared with PoO and GxE analyses. However, because PoOxE effects are measured as ratios of RRRs (see Eq (1)), it is hypothetically possible that PoOxE effects are larger than PoO effects or GxE effects, in particular in the presence of ‘qualitative interactions’, where effects are in opposite directions across strata. This is illustrated in S1 Appendix, and may partly explain some of the large effects in Tables Tables22 to to4.4. Under such scenarios, some of the lost power may be regained. Nevertheless, none of the q-values were lower than 0.14, which suggests that low power may have been an issue in this study. Still, several SNPs had q-values below 0.5, meaning that we expect fewer than half of them to be false positives. SNPs presented in Tables Tables22 to to44 should be interpreted as candidates to be further investigated in other studies. The next steps would be to replicate these candidates in other data sets/populations, followed by targeted functional analyses to help elucidate the importance of these SNPs in the interplay between environmental factors and risk of CPO.

To summarize, this study presents new methodology, implemented in the R-package Haplin, to investigate PoOxE effects in the context of family trios or duos. Our analyses pointed to several SNPs with PoOxSmoke effects in the European sample. We were unable to assess the generalizability of this finding across ethnicities, because few of the Asian mothers smoked cigarettes or consumed alcohol. We did not find any evidence for PoOxAlcohol effects in the European sample, and there were no PoOxVitamin effects in either ethnicity. Still, these analyses highlight the versatility of Haplin in studying complex disease models.

Materials and methods

Study participants

The majority of the participants belonged to one of two major ethnicities (European or Asian). Table 1 outlines the population distribution by ethnicity and trio completeness, and Table 6 summarizes characteristics of the maternal exposures by ethnicity.

Quality control

Genotypes for 569 244 SNPs were available for the current analyses. The PLINK software [35] was used for quality control, with the following criteria applied for excluding SNPs: (i) >5% missing genotype for a given SNP, (ii) minor allele frequency (MAF) <5%, (iii) Hardy-Weinberg equilibrium (HWE) p-value <0.001 for parental alleles, (iv) >10% Mendelian error rate, and finally (v) linkage disequilibrium (LD) of r2 = 1 with other SNPs (to exclude SNPs with redundant information due to complete LD). Overall, genotypes for 550 families with isolated CPO were available for the current analyses. Criteria for excluding individuals were: (vi) >10% missing genotype within an individual, and (vii) >5% Mendelian errors within a family. Table 7 provides the total number of individuals after the above pruning. Because none of the families had Mendelian error rates >5%, they were all retained in the analyses. The total number of SNPs remaining after quality control is shown in Table 8, along with the different criteria used for pruning.

Table 7
Individuals before and after pruning.
Table 8
SNPs before and after pruning.

Statistical analysis

All analyses were conducted using the statistical software package Haplin, http://people.uib.no/gjessing/genetics/software/haplin. Haplin is particularly tailored to the analysis of offspring-parent trios and duos, but is also applicable to case-control data [17]. It is implemented as a package in the statistical programming language R [36]. We applied the function haplinSlide to analyze all SNPs sequentially. For each SNP, a log-linear maximum likelihood model is applied to the trio genotype frequencies, allowing different risk of disease (penetrance) depending on the parent of origin of the allele. The effect of each SNP was assumed to be multiplicative in allele dose, with the most common (major) allele used as reference. Missing alleles were imputed using the EM-algorithm; standard errors and p-values were corrected for this imputation [17].

The following section outlines how the PoOxE effects are computed in Haplin. First, a PoO analysis is performed for each stratum of an exposure, where S = 0 represents the unexposed and S = 1 the exposed. The PoO analysis in stratum S computes two relative risks

RRmat(S)=P(CPO|pat=a,mat=a1,S)P(CPO|pat=a,mat=a0,S)

for a maternally inherited allele, and

RRpat(S)=P(CPO|pat=a1,mat=a,S)P(CPO|pat=a0,mat=a,S)

for a paternally inherited allele, where a0 is the reference allele, a1 is the alternative allele, and “a” denotes any one of the two alleles. The PoO relative risk ratio (RRRPoO) then compares the two separate relative risks, so that

RRRPoO(S)=RRmat(S)RRpat(S).

RRRPoO = 1 means a1 increases (or decreases) the risk by the same amount regardless of whether the allele is maternally or paternally inherited. Next, Haplin compares the RRRPoO for all strata. In the case of two strata, S = 0 represents the unexposed and S = 1 the exposed, and Haplin tests whether RRRPoO(0) = RRRPoO(1). The test is performed as a Wald test by exploiting the fact that the estimated log(RRRPoO(0)) and log(RRRPoO(1)) are independent and asymptotically normally distributed, as outlined in Skare et al. (2012) [14] and Gjerdevik et al. (2017) [16].

P-values from the PoOxE analyses were displayed in a QQ-plot, with expected p-values plotted against the observed. Under the null hypothesis of no PoOxE effect, all SNPs should lie along the diagonal line representing a uniform distribution, whereas significant SNPs are expected to appear markedly above the diagonal line and outside the confidence bands.

To visualize the strength of the association signal and regional information flanking the most significant SNPs, we used a modified version of the R-script for regional plots available at http://www.broadinstitute.org/files/shared/diabetes/scandinavs/assocplot.R. The plot also displays the degree of LD between top SNPs and neighboring SNPs, recombination patterns, and positional information about genes in the region [37].

To assess the a priori power to detect PoOxE effects with our model, we performed power simulations based on 1000 replications and a significance level of 0.05 (Fig 6). The black line shows the power for a PoOxE analysis based on 500 case-parent trios (consistent with the sample size in this study), a MAF of 0.20, and equally-sized exposed and unexposed groups. The left panel of Fig 6 depicts different sample sizes and the middle panel depicts different MAFs. The right panel shows the power for different etiologic scenarios (child, PoO, GxE, and PoOxE). The child effect is the direct risk associated with the allele when it is carried by the child, regardless of parental origin or environmental exposures. The PoO effect is the risk associated with maternally-inherited alleles relative to paternally-inherited alleles. The GxE effect is the ratio of RRs in the two exposure groups. Finally, the PoOxE effect is the maternal to paternal risk ratio for the exposed divided by the same ratio for the unexposed.

Ethics approvals

This specific study did not need approval from an ethics committee because ethics approvals for the consortium were obtained from the respective ethics committees at each institution contributing data to the consortium. Details have been provided in our original publication [19].

Supporting information

S1 Fig

Manhattan plots for the different exposures in the analyses of the pooled sample.

SNPs with p-values below 10−5 are in blue.

(TIFF)

S2 Fig

Manhattan plots for the different exposures in the analyses of the European sample.

SNPs with p-values below 10−5 are in blue.

(TIFF)

S3 Fig

Manhattan plots for the different exposures in the analyses of the Asian sample.

(TIFF)

S1 Appendix

Example code for PoOxE analysis.

(DOCX)

Acknowledgments

We are indebted to the families who contributed to this study, and the orofacial cleft consortium as a whole. We also sincerely thank everyone involved in the recruitment process and the genotyping of DNA from the families.

This research was supported by the Bergen Medical Research Foundation, grant 807191 (AJ, HKG, RTL), in part by the Intramural Program of the National Institute of Environmental Health Sciences, NIH/NIEHS (AJW), by NIH grant DE08559 (JCM), and by the Research Council of Norway through its Centres of Excellence funding scheme, project number 262700 (HKG, AJ).

Funding Statement

This research was supported by the Bergen Medical Research Foundation, grant 807191, in part by the Intramural Program of the National Institute of Environmental Health Sciences, National Institutes of Health (NIH/NIEHS), by NIH grant DE08559, and partly by the Research Council of Norway through its Centres of Excellence funding scheme, project number 262700. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

Data are available from the dbGAP database (http://www.ncbi.nlm.nih.gov/gap) under accession ID phs000094.v1.p1.

References

1. Mossey P, Castilla EE. Global registry and database on craniofacial anomalies. World Health Organization, Geneva, 2003.
2. Mossey PA, Little J, Munger RG, Dixon MJ, Shaw WC (2009) Cleft lip and palate. Lancet 374: 1773–85. doi: 10.1016/S0140-6736(09)60695-4 [PubMed]
3. Grosen D, Chevrier C, Skytthe A, Bille C, Molsted 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. Journal of medical genetics 47: 162–8. doi: 10.1136/jmg.2009.069385 [PMC free article] [PubMed]
4. Sivertsen A, Wilcox AJ, Skjaerven R, Vindenes HA, 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 (Clinical research ed 336: 432–4. [PMC free article] [PubMed]
5. Lie RT, Wilcox AJ, Skjaerven R (1994) A population-based study of the risk of recurrence of birth defects. N Engl J Med 331: 1–4. doi: 10.1056/NEJM199407073310101 [PubMed]
6. Christensen K, Mitchell LE (1996) Familial recurrence-pattern analysis of nonsyndromic isolated cleft palate—a Danish Registry study. American journal of human genetics 58: 182–90. [PubMed]
7. Grosen D, Bille C, Pedersen JK, Skytthe A, Murray JC, Christensen K (2010) Recurrence risk for offspring of twins discordant for oral cleft: a population-based cohort study of the Danish 1936–2004 cleft twin cohort. American journal of medical genetics Part A 152A: 2468–74. doi: 10.1002/ajmg.a.33608 [PMC free article] [PubMed]
8. Grosen D, Bille C, Petersen I, Skytthe A, Hjelmborg JB, Pedersen JK, et al. (2011) Risk of oral clefts in twins. Epidemiology 22: 313–9. doi: 10.1097/EDE.0b013e3182125f9c [PMC free article] [PubMed]
9. Lawson HA, Cheverud JM, Wolf JB (2013) Genomic imprinting and parent-of-origin effects on complex traits. Nature reviews 14: 609–17. [PMC free article] [PubMed]
10. Jugessur A, Farlie PG, Kilpatrick N (2009) The genetics of isolated orofacial clefts: from genotypes to subphenotypes. Oral diseases 15: 437–53. doi: 10.1111/j.1601-0825.2009.01577.x [PubMed]
11. Rahimov F, Jugessur A, Murray JC (2012) Genetics of nonsyndromic orofacial clefts. Cleft Palate Craniofac J 49: 73–91. doi: 10.1597/10-178 [PMC free article] [PubMed]
12. Dixon MJ, Marazita ML, Beaty TH, Murray JC (2011) Cleft lip and palate: understanding genetic and environmental influences. Nature reviews Genetics 12: 167–78. doi: 10.1038/nrg2933 [PMC free article] [PubMed]
13. Marazita ML (2012) The Evolution of Human Genetic Studies of Cleft Lip and Cleft Palate. Annual review of genomics and human genetics. [PMC free article] [PubMed]
14. Skare O, Jugessur A, Lie RT, Wilcox AJ, Murray JC, Lunde A, et al. (2012) Application of a Novel Hybrid Study Design to Explore Gene-Environment Interactions in Orofacial Clefts. Ann Hum Genet 76: 221–36. doi: 10.1111/j.1469-1809.2012.00707.x [PMC free article] [PubMed]
15. Wang S, Yu Z, Miller RL, Tang D, Perera FP (2011) Methods for detecting interactions between imprinted genes and environmental exposures using birth cohort designs with mother-offspring pairs. Hum Hered 71: 196–208. doi: 10.1159/000328006 [PMC free article] [PubMed]
16. Gjerdevik M, Haaland ØA, Romanovska J, Lie RT, Jugessur A, Gjessing HK (2017) Parent-of-Origin-Environment Interactions in Case-Parent Triads With or Without Independent Controls. Ann Hum Genet. doi: 10.1111/ahg.12224
17. Gjessing HK, Lie RT (2006) Case-parent triads: estimating single- and double-dose effects of fetal and maternal disease gene haplotypes. Ann Hum Genet 70: 382–96. doi: 10.1111/j.1529-8817.2005.00218.x [PubMed]
18. Jugessur A, Skare Ø, Harris JR, Lie RT, Gjessing HK (2012) Using offspring-parent triads to study complex traits: A tutorial based on orofacial clefts. Norwegian Journal of Epidemiology 21: 251–67.
19. Beaty TH, Murray JC, Marazita ML, Munger RG, Ruczinski I, Hetmanski JB, et al. (2010) A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4. Nature genetics 42: 525–9. doi: 10.1038/ng.580 [PMC free article] [PubMed]
20. Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proceedings of the National Academy of Sciences of the United States of America 100: 9440–5. doi: 10.1073/pnas.1530509100 [PubMed]
21. Hager R, Cheverud JM, Wolf JB (2008) Maternal effects as the cause of parent-of-origin effects that mimic genomic imprinting. Genetics 178: 1755–62. doi: 10.1534/genetics.107.080697 [PubMed]
22. Borg K, Stankiewicz P, Bocian E, Kruczek A, Obersztyn E, Lupski JR, et al. (2005) Molecular analysis of a constitutional complex genome rearrangement with 11 breakpoints involving chromosomes 3, 11, 12, and 21 and a approximately 0.5-Mb submicroscopic deletion in a patient with mild mental retardation. Human genetics 118: 267–75. doi: 10.1007/s00439-005-0021-0 [PubMed]
23. Whitaker HC, Shiong LL, Kay JD, Gronberg H, Warren AY, Seipel A, et al. (2014) N-acetyl-L-aspartyl-L-glutamate peptidase-like 2 is overexpressed in cancer and promotes a pro-migratory and pro-metastatic phenotype. Oncogene 33: 5274–87. doi: 10.1038/onc.2013.464 [PubMed]
24. Du J, Fan Z, Ma X, Gao Y, Wu Y, Liu S, et al. (2011) Expression of Dpp6 in mouse embryonic craniofacial development. Acta Histochem 113: 636–9. doi: 10.1016/j.acthis.2010.08.001 [PubMed]
25. Du J, Fan Z, Ma X, Wu Y, Liu S, Gao Y, et al. (2014) Expression of DPP6 in Meckel's cartilage and tooth germs during mouse facial development. Biotech Histochem 89: 14–8. doi: 10.3109/10520295.2013.795661 [PubMed]
26. Blanton SH, Henry RR, Yuan Q, Mulliken JB, Stal S, Finnell RH, et al. (2011) Folate pathway and nonsyndromic cleft lip and palate. Birth Defects Res A Clin Mol Teratol 91: 50–60. doi: 10.1002/bdra.20740 [PMC free article] [PubMed]
27. Krey KF, Dannhauer KH, Hemprich A, Zaitsev S, Bankfalvi A, Buchwalow IB, et al. (2002) Cytophotometrical and immunohistochemical analysis of soft palate muscles of children with isolated cleft palate and combined cleft lip and palate. Exp Toxicol Pathol 54: 69–75. doi: 10.1078/0940-2993-00235 [PubMed]
28. Corcos L, Lucas D, Le Jossic-Corcos C, Dreano Y, Simon B, Plee-Gautier E, et al. (2012) Human cytochrome P450 4F3: structure, functions, and prospects. Drug Metabol Drug Interact 27: 63–71. doi: 10.1515/dmdi-2011-0037 [PubMed]
29. Ng M, Freeman MK, Fleming TD, Robinson M, Dwyer-Lindgren L, Thomson B, et al. (2014) Smoking prevalence and cigarette consumption in 187 countries, 1980–2012. JAMA 311: 183–92. doi: 10.1001/jama.2013.284692 [PubMed]
30. Yang W, Lu J, Weng J, Jia W, Ji L, Xiao J, et al. (2010) Prevalence of diabetes among men and women in China. N Engl J Med 362: 1090–101. doi: 10.1056/NEJMoa0908292 [PubMed]
31. Beaty TH, Ruczinski I, Murray JC, Marazita ML, Munger RG, Hetmanski JB, et al. (2011) Evidence for gene-environment interaction in a genome wide study of nonsyndromic cleft palate. Genetic epidemiology 35: 469–78. doi: 10.1002/gepi.20595 [PMC free article] [PubMed]
32. Sharp GC, Ho K, Davies A, Stergiakouli E, Humphries K, McArdle W, et al. (2017) Distinct DNA methylation profiles in subtypes of orofacial cleft. Clin Epigenetics 9: 63 doi: 10.1186/s13148-017-0362-2 [PMC free article] [PubMed]
33. Sharp GC, Stergiakouli E, Sandy J, Relton C (2017) Epigenetics and Orofacial Clefts: A Brief Introduction. Cleft Palate Craniofac J. [PubMed]
34. Smith ZD, Meissner A (2013) DNA methylation: roles in mammalian development. Nature reviews 14: 204–20. doi: 10.1038/nrg3354 [PubMed]
35. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics 81: 559–75. doi: 10.1086/519795 [PubMed]
36. R Core Team (2015) R: A Language and Environment for Statistical Computing.
37. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. (2010) LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26: 2336–7. doi: 10.1093/bioinformatics/btq419 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science