Assessing Maternal Genetic Associations

Genetic Analysis of the Single-Nucleotide Polymorphisms rs880810, rs545793, rs80094639, and rs13251901 in Nonsyndromic Oral Clefts: A Case–Parent Trio Study

Oral clefts, including cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP), are the most common types of congenital anomalies of the human face. Various genetic and environmental factors play a role in developing oral clefts. Several studies have shown the association of the PAX7 gene and the 8q24 region with these oral clefts in different populations worldwide. However, there are no reported studies on the possible connection between the PAX7 gene and the 8q24 region nucleotide variants and the risk of developing nonsyndromic oral clefts (NSOC) in the Indian population. Hence, this study aimed to test the possible association between PAX7 gene single-nucleotide polymorphisms (SNPs) rs880810, rs545793,rs80094639, and rs13251901 of the 8q24 region using a case-parent trio design. Forty case-parent trios were selected from the CLP center. Genomic DNA was isolated from the cases and their parents. The rs880810, rs545793, rs80094639, and rs13251901 were genotyped by the MassARRAY technique. PLINK software was used for statistical analysis. All the SNPs were tested for Hardy-Weinberg equilibrium. No statistical significance was found with any SNPs, as none of the genotyped SNPs showed a p -value of less than 0.05. Hence, the rs880810, rs545793, and rs80094639 of the PAX7 gene, and rs13251901 of the 8q24 region are not associated with NSOC in the Indian population.

Evidence of gene-environment interaction for two genes on chromosome 4 and environmental tobacco smoke in controlling the risk of nonsyndromic cleft palate

Nonsyndromic cleft palate (CP) is one of the most common human birth defects and both genetic and environmental risk factors contribute to its etiology. We conducted a genome-wide association study (GWAS) using 550 CP case-parent trios ascertained in an international consortium. Stratified analysis among trios with different ancestries was performed to test for GxE interactions with common maternal exposures using conditional logistic regression models. While no single nucleotide polymorphism (SNP) achieved genome-wide significance when considered alone, markers in SLC2A9 and the neighboring WDR1 on chromosome 4p16.1 gave suggestive evidence of gene-environment interaction with environmental tobacco smoke (ETS) among 259 Asian trios when the models included a term for GxE interaction. Multiple SNPs in these two genes were associated with increased risk of nonsyndromic CP if the mother was exposed to ETS during the peri-conceptual period (3 months prior to conception through the first trimester). When maternal ETS was considered, fifteen of 135 SNPs mapping to SLC2A9 and 9 of 59 SNPs in WDR1 gave P values approaching genome-wide significance (10(-6)<P<10(-4)) in a test for GxETS interaction. SNPs rs3733585 and rs12508991 in SLC2A9 yielded P = 2.26×10(-7) in a test for GxETS interaction. SNPs rs6820756 and rs7699512 in WDR1 also yielded P = 1.79×10(-7) and P = 1.98×10(-7) in a 1 df test for GxE interaction. Although further replication studies are critical to confirming these findings, these results illustrate how genetic associations for nonsyndromic CP can be missed if potential GxE interaction is not taken into account, and this study suggest SLC2A9 and WDR1 should be considered as candidate genes for CP.

Fetal and maternal genes' influence on gestational age in a quantitative genetic analysis of 244,000 Swedish births

Although there is increasing evidence that genetic factors influence gestational age, it is unclear to what extent this is due to fetal and/or maternal genes. In this study, we apply a novel analytical model to estimate genetic and environmental contributions to pregnancy history records obtained from 165,952 Swedish families consisting of offspring of twins, full siblings, and half-siblings (1987-2008). Results indicated that fetal genetic factors explained 13.1% (95% confidence interval (CI): 6.8, 19.4) of the variation in gestational age at delivery, while maternal genetic factors accounted for 20.6% (95% CI: 18.1, 23.2). The largest contribution to differences in the timing of birth were environmental factors, of which 10.1% (95% CI: 7.0, 13.2) was due to factors shared by births of the same mother, and 56.2% (95% CI: 53.0, 59.4) was pregnancy specific. Similar models fit to the same data dichotomized at clinically meaningful thresholds (e.g., preterm birth) resulted in less stable parameter estimates, but the collective results supported a model of homogeneous genetic and environmental effects across the range of gestational age. Since environmental factors explained most differences in the timing of birth, genetic studies may benefit from understanding the specific effect of fetal and maternal genes in the context of these yet-unidentified factors.

Evidence of gene-environment interaction for the RUNX2 gene and environmental tobacco smoke in controlling the risk of cleft lip with/without cleft palate

This study examined the association between 49 markers in the Runt-related transcription factor 2 (RUNX2) gene and nonsyndromic cleft lip with/without cleft palate (CL/P) among 326 Chinese case-parent trios, while considering gene-environment (GxE) interaction and parent-of-origin effects. Five single-nucleotide polymorphisms (SNPs) showed significant evidence of linkage and association with CL/P and these results were replicated in an independent European sample of 825 case-parent trios. We also report compelling evidence for interaction between markers in RUNX2 and environmental tobacco smoke (ETS). Although most marginal SNP effects (i.e., ignoring maternal exposures) were not statistically significant, eight SNPs were significant when considering possible interaction with ETS when testing for gene (G) and GxE interaction simultaneously or when considering GxE alone. Independent samples from European populations showed consistent evidence of significant GxETS interaction at two SNPs (rs6904353 and rs7748231). Our results suggest genetic variation in RUNX2 may influence susceptibility to CL/P through interacting with ETS.

MI-GWAS: a SAS platform for the analysis of inherited and maternal genetic effects in genome-wide association studies using log-linear models

BACKGROUND

Several platforms for the analysis of genome-wide association data are available. However, these platforms focus on the evaluation of the genotype inherited by affected (i.e. case) individuals, whereas for some conditions (e.g. birth defects) the genotype of the mothers of affected individuals may also contribute to risk. For such conditions, it is critical to evaluate associations with both the maternal and the inherited (i.e. case) genotype. When genotype data are available for case-parent triads, a likelihood-based approach using log-linear modeling can be used to assess both the maternal and inherited genotypes. However, available software packages for log-linear analyses are not well suited to the analysis of typical genome-wide association data (e.g. including missing data).

RESULTS

An integrated platform, Maternal and Inherited Analyses for Genome-wide Association Studies (MI-GWAS) for log-linear analyses of maternal and inherited genetic effects in large, genome-wide datasets, is described. MI-GWAS uses SAS and LEM software in combination to appropriately format data, perform the log-linear analyses and summarize the results. This platform was evaluated using existing genome-wide data and was shown to perform accurately and relatively efficiently.

CONCLUSIONS

The MI-GWAS platform provides a valuable tool for the analysis of association of a phenotype or condition with maternal and inherited genotypes using genome-wide data from case-parent triads. The source code for this platform is freely available at http://www.sph.uth.tmc.edu/sbrr/mi-gwas.htm.

Evidence of gene-environment interaction for the IRF6 gene and maternal multivitamin supplementation in controlling the risk of cleft lip with/without cleft palate

Although multiple genes have been identified as genetic risk factors for isolated, non-syndromic cleft lip with/without cleft palate (CL/P), a complex and heterogeneous birth defect, interferon regulatory factor 6 gene (IRF6) is one of the best documented genetic risk factors. In this study, we tested for association between markers in IRF6 and CL/P in 326 Chinese case-parent trios, considering gene-environment interaction for two common maternal exposures, and parent-of-origin effects. CL/P case-parent trios from three sites in mainland China and Taiwan were genotyped for 22 single nucleotide polymorphisms (SNPs) in IRF6. The transmission disequilibrium test was used to test for marginal effects of individual SNPs. We used PBAT to screen the SNPs and haplotypes for gene-environment (G×E) interaction and conditional logistic regression models to quantify effect sizes for SNP-environment interaction. After Bonferroni correction, 14 SNPs showed statistically significant association with CL/P. Evidence of G×E interaction was found for both maternal exposures, multivitamin supplementation and environmental tobacco smoke (ETS). Two SNPs showed evidence of interaction with multivitamin supplementation in conditional logistic regression models (rs2076153 nominal P=0.019, rs17015218 nominal P=0.012). In addition, rs1044516 yielded evidence for interaction with maternal ETS (nominal P=0.041). Haplotype analysis using PBAT also suggested interaction between SNPs in IRF6 and both multivitamin supplementation and ETS. However, no evidence for maternal genotypic effects or significant parent-of-origin effects was seen in these data. These results suggest IRF6 gene may influence risk of CL/P through interaction with multivitamin supplementation and ETS in the Chinese population.

Maternal immunoglobulin E and childhood leukemia

Childhood leukemia, particularly acute lymphoblastic leukemia (ALL), has long been hypothesized to be affected by abnormal immune responses to microbial challenges stemming from a lack of immune modulation in early childhood. Studies of allergies suggest that a child's immune development may be modulated by maternal immune status. We conducted a study to explore the relationship between maternal immunoglobulin E (IgE) and childhood leukemia and to investigate whether maternal immune status can influence childhood leukemia risk. Serum total and specific IgE (respiratory and food) were measured in biological mothers of 352 children (193 healthy controls and 159 leukemia cases, including 139 ALL cases) ages <8 years who were enrolled in the Northern California Childhood Leukemia Study. Odds ratios associated with maternal IgE were calculated using unconditional logistic regression adjusted for child's age, sex, race/ethnicity, and annual household income. A positive association between childhood leukemia or ALL and elevated levels of maternal serum total IgE was observed, especially among Hispanics. In addition, a positive association was observed between childhood leukemia or ALL and maternal respiratory or food IgE status. These results suggest that maternal immune function may play a crucial role in the etiology of childhood leukemia, although additional studies need to be conducted to confirm the results of this study and provide a perspective on mechanisms.

Evidence that TGFA influences risk to cleft lip with/without cleft palate through unconventional genetic mechanisms

This study examined the association between markers in transforming growth factor alpha (TGFA) and isolated, non-syndromic cleft lip with/without palate (CL/P) using a case-parent trio design, considering parent-of-origin effects. We also tested for gene-environmental interaction with common maternal exposures, and for gene-gene interaction using markers in TGFA and another recognized causal gene, IRF6. CL/P case-parent trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 17 single nucleotide polymorphisms (SNPs) in TGFA. The transmission disequilibrium test was used to test individual SNPs, and the parent-of-origin likelihood ratio test (PO-LRT) was used to assess parent-of-origin effects. We also screened for possible gene-environment interaction using PBAT, and tested for gene-gene interaction using conditional logistic regression models. When all trios were combined, four SNPs showed significant excess maternal transmission, two of which gave significant PO-LRT values [rs3821261: P = 0.004 and OR(imprinting) = 4.17; and rs3771475: P = 0.027 and OR(imprinting) = 2.44]. Haplotype analysis of these two SNPS also supported excess maternal transmission. We saw intriguing but suggestive evidence of G x E interaction for several SNPs in TGFA when either individual SNPs or haplotypes of adjacent SNPs were considered. Thus, TGFA appears to influence risk of CL/P through unconventional means with an apparent parent-of-origin effect (excess maternal transmission) and possible interaction with maternal exposures.

Differential parental transmission of markers in BCL3 among Korean cleft case-parent trios

OBJECTIVES

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of approximately 1 in 700 live births. The B-Cell Leukemia/lymphoma 3 (BCL3) gene has been suggested as a candidate gene for CL/P based on association and linkage studies in some populations. This study tests for an association between markers in BCL3 and isolated, non-syndromic CL/P using a case-parent trio design, while considering parent-of-origin effects.

METHODS

Forty case-parent trios were genotyped for two single nucleotide polymorphisms (SNPs) in the BCL3 gene. We performed a transmission disequilibrium test (TDT) on individual SNPs, and the FAMHAP package was used to estimate haplotype frequencies and to test for excess transmission of multi-SNP haplotypes.

RESULTS

The odds ratio for transmission of the minor allele, OR (transmission), was significant for SNP rs8100239 (OR=3.50, p=0.004) and rs2965169 (OR=2.08, p=0.027) when parent-of-origin was not considered. Parent-specific TDT revealed that SNP rs8100239 showed excess maternal transmission. Analysis of haplotypes of rs2965169 and rs8100239 also suggested excess maternal transmission.

CONCLUSIONS

BCL3 appears to influence risk of CL/P through a parent-of-origin effect with excess maternal transmission.

Maternal transmission effects of the PAX genes among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of 1 in 700 live births. The paired box (PAX) genes have been suggested as candidate genes for CL/P based largely on mouse models; however, few human studies have focused on this gene family. This study tests for association between markers in four PAX genes and CL/P using a case-parent trio design considering parent-of-origin effects. Trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 34 single nucleotide polymorphisms (SNPs) in the PAX3, PAX6, PAX7, and PAX9 genes. We performed the transmission disequilibrium test (TDT) on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test (TAT) and the parent-of-origin likelihood ratio test (PO-LRT). TDT analysis showed one SNP (rs766325) in PAX7 yielding evidence of linkage and association when parent-of-origin was not considered, with an OR(transmission)=1.62 (P=0.003), and five SNPs in PAX6 (including two pairs in near perfect linkage disequilibrium). TAT analysis of all trios revealed two SNPs in PAX7 and four SNPs in PAX3 showing significant excess maternal transmission. For these six SNPs, the maternal OR(transmission) ranged between 1.74 and 2.40, and PO-LRT was also significant (P-values=0.035-0.012). When this analysis was limited to trios with male cases, SNPs in PAX7 showed higher maternal OR(transmission) and greater significance. PAX genes may influence the risk of CL/P through maternal effects, possibly imprinting, which seems to be stronger among male cases.

Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence around 1 in 700 live births. The Runt-related transcription factor 2 (RUNX2) gene has been suggested as a candidate gene for CL/P based largely on mouse models; however, no human studies have focused on RUNX2 as a risk factor for CL/P. This study examines the association between markers in RUNX2 and isolated, nonsyndromic CL/P using a case-parent trio design, while considering parent-of-origin effects. Case-parent trios from four populations (77 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 24 single nucleotide polymorphisms (SNPs) in the RUNX2 gene. We performed the transmission disequilibrium test on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test and the parent-of-origin likelihood ratio test (PO-LRT). When all trios were combined, the transmission asymmetry test revealed a block of 11 SNPs showing excess maternal transmission significant at the P<0.01 level, plus one SNP (rs1934328) showing excess paternal transmission (P=0.002). For the 11 SNPs showing excess maternal transmission, odds ratios of being transmitted to the case from the mother ranged between 3.00 and 4.00. The parent-of-origin likelihood ratio tests for equality of maternal and paternal transmission were significant for three individual SNPs (rs910586, rs2819861, and rs1934328). Thus, RUNX2 appears to influence risk of CL/P through a parent-of-origin effect with excess maternal transmission.

Excess maternal transmission of markers in TCOF1 among cleft palate case-parent trios from three populations

Isolated cleft palate is among the most common human birth defects. The TCOF1 gene has been suggested as a candidate gene for cleft palate based on animal models. This study tests for association between markers in TCOF1 and isolated, nonsyndromic cleft palate using a case-parent trio design considering parent-of-origin effects. Case-parent trios from three populations (comprising a total of 81 case-parent trios) were genotyped for single nucleotide polymorphisms (SNPs) in the TCOF1 gene. We used the transmission disequilibrium test and the transmission asymmetry test on individual SNPs. When all trios were combined, the odds ratio for transmission of the minor allele, OR(transmission), was significant for SNP rs15251 (OR = 2.88, P = 0.007), as well as rs2255796 and rs2569062 (OR = 2.08, P = 0.03; OR = 2.43, P = 0.041; respectively) when parent of origin was not considered. The transmission asymmetry test also revealed one SNP (rs15251) showing excess maternal transmission significant at the P = 0.005 level (OR = 6.50). Parent-of-origin effects were assessed using the parent-of-origin likelihood ratio test on both SNPs and haplotypes. While the parent-of-origin likelihood ratio test was only marginally significant for this SNP (P = 0.136), analysis of haplotypes of rs2255796 and rs15251 suggested excess maternal transmission. Therefore, these data suggest TCOF1 may influence risk of cleft palate through a parent-of-origin effect.

Analytical approaches to detect maternal/fetal genotype incompatibilities that increase risk of pre-eclampsia

BACKGROUND

In utero interactions between incompatible maternal and fetal genotypes are a potential mechanism for the onset or progression of pregnancy related diseases such as pre-eclampsia (PE). However, the optimal analytical approach and study design for evaluating incompatible maternal/offspring genotype combinations is unclear.

METHODS

Using simulation, we estimated the type I error and power of incompatible maternal/offspring genotype models for two analytical approaches: logistic regression used with case-control mother/offspring pairs and the log-linear regression used with case-parent triads. We evaluated a real dataset consisting of maternal/offspring pairs with and without PE for incompatibility effects using the optimal analysis based on the results of the simulation study.

RESULTS

We identified a single coding scheme for the incompatibility effect that was equally or more powerful than all of the alternative analysis models evaluated, regardless of the true underlying model for the incompatibility effect. In addition, the log-linear regression was more powerful than the logistic regression when the heritability was low, and more robust to adjustment for maternal or fetal effects. For the PE data, this analysis revealed three genes, lymphotoxin alpha (LTA), von Willebrand factor (VWF), and alpha 2 chain of type IV collagen (COL4A2) with possible incompatibility effects.

CONCLUSION

The incompatibility model should be evaluated for complications of pregnancy, such as PE, where the genotypes of two individuals may contribute to the presence of disease.

Using risk-based sampling to enrich cohorts for endpoints, genes, and exposures

Targeting the first-degree relatives of people with a particular complex disease can offer a powerful approach to building a risk-based cohort for prospective studies of etiologic factors. Such a cohort provides both a sizable increase in the rate of accrual of newly incident cases, enriching for risk factors that are known or even unknown, and a high level of motivation among participants. A nationwide study of breast cancer in the United States and Puerto Rico, the Sister Study, made up of women who are each the sister of a woman with breast cancer, exemplifies this approach. In this paper, the authors provide power calculations to aid in the design of such studies and quantify their benefits for detecting both genetic variants related to risk and interactive effects of genetic and environmental factors. While the risk-based cohort can have markedly increased prevalences of rare causative alleles, most of the power advantages for this design is due to the increased rate of accrual of newly incident cases rather than the increase in any one individual allele.

Studying genetic predisposition among small-for-gestational-age newborns

The objective of this report is to provide a summary overview of genetic association studies for the small-for-gestational-age (SGA) pregnancy outcome. Only the thrombophilia and xenobiotic-metabolizing genetic pathways were studied with any frequency. Most studies used case-control designs and analyzed only the maternal genotype. A brief critique of some features of the published studies is presented: it addresses mainly the selection of controls, study power, the need to evaluate gene-environment interaction, and the potential for population stratification bias, believed likely to affect such studies. Alternative designs, not vulnerable to the population structure bias, are also discussed; they include case-parental trios and a mixture of both case trios and case-control data. Aspects that have almost never been considered in the published studies, but that are particularly relevant for adverse pregnancy outcomes, are maternally mediated and parent-of-origin effects. These are defined and methods to evaluate them are briefly presented.

Power calculations for likelihood ratio tests for offspring genotype risks, maternal effects, and parent-of-origin (POO) effects in the presence of missing parental genotypes when unaffected siblings are available

Genotype-based likelihood-ratio tests (LRT) of association that examine maternal and parent-of-origin effects have been previously developed in the framework of log-linear and conditional logistic regression models. In the situation where parental genotypes are missing, the expectation-maximization (EM) algorithm has been incorporated in the log-linear approach to allow incomplete triads to contribute to the LRT. We present an extension to this model which we call the Combined_LRT that incorporates additional information from the genotypes of unaffected siblings to improve assignment of incompletely typed families to mating type categories, thereby improving inference of missing parental data. Using simulations involving a realistic array of family structures, we demonstrate the validity of the Combined_LRT under the null hypothesis of no association and provide power comparisons under varying levels of missing data and using sibling genotype data. We demonstrate the improved power of the Combined_LRT compared with the family-based association test (FBAT), another widely used association test. Lastly, we apply the Combined_LRT to a candidate gene analysis in Autism families, some of which have missing parental genotypes. We conclude that the proposed log-linear model will be an important tool for future candidate gene studies, for many complex diseases where unaffected siblings can often be ascertained and where epigenetic factors such as imprinting may play a role in disease etiology.

The C677T polymorphism in the methylenetetrahydrofolate reductase gene (MTHFR), maternal use of folic acid supplements, and risk of isolated clubfoot: A case-parent-triad analysis

Worldwide, 1-4 per 1,000 births are affected by clubfoot. Clubfoot etiology is unclear, but both genetic and environmental factors are thought to be involved. Low folate status in pregnant women has been implicated in several congenital malformations, and folate metabolism may be affected by polymorphisms in the methylenetetrahydrofolate reductase gene (MTHFR). Using a case-parent-triad design, the authors investigated whether the MTHFR C677T polymorphism, and maternal periconceptional folic acid supplement use, influenced risk of isolated clubfoot. Three hundred seventy-five United Kingdom case-parent triads were recruited in 1998-1999. Among the children, there was a significant trend of decreasing clubfoot risk with increasing number of T alleles: relative risk for CT vs. CC = 0.75, 95% confidence interval: 0.57, 0.97; relative risk for TT vs. CC = 0.57, 95% confidence interval: 0.35, 0.91; p trend = 0.006. This association was not modified by maternal folic acid use. Maternal MTHFR genotype did not influence clubfoot risk for the offspring overall, although a possible interaction with folic acid use was found. This is the first known report of a specific genetic polymorphism associated with clubfoot. The direction of the association is intriguing and suggests that DNA synthesis may be relevant in clubfoot development. However, clubfoot mechanisms are poorly understood, and the folate metabolism pathway is complex. Further research is needed to elucidate these relations.