Genome-wide discovery of maternal effect variants
© Kent et al; licensee BioMed Central Ltd. 2009
Published: 15 December 2009
Skip to main content
© Kent et al; licensee BioMed Central Ltd. 2009
Published: 15 December 2009
Many phenotypes may be influenced by the prenatal environment of the mother and/or maternal care, and these maternal effects may have a heritable component. We have implemented in the computer program SOLAR a variance components-based method for detecting indirect effects of maternal genotype on offspring phenotype. Of six phenotypes measured in three generations of the Framingham Heart Study, height showed the strongest evidence (P = 0.02) of maternal effect. We conducted a genome-wide association analysis for height, testing both the direct effect of the focal individual's genotype and the indirect effect of the maternal genotype. Offspring height showed suggestive evidence of association with maternal genotype for two single-nucleotide polymorphisms in the trafficking protein particle complex 9 gene TRAPPC9 (NIBP), which plays a role in neuronal NF-κB signalling. This work establishes a methodological framework for identifying genetic variants that may influence the contribution of the maternal environment to offspring phenotypes.
Many phenotypes may be influenced by the prenatal environment of the mother and/or maternal care, and these maternal effects may have a heritable component. Much research has focused on the impact of measurable properties of the mother (e.g., adiposity, diabetes, alcohol, or tobacco use) on subsequent phenotypes in their children (e.g., birthweight , insulin resistance , cognitive function ). A more general question is: does the mother's measured genotype influence offspring phenotypes, whether or not the intermediate maternal phenotypes are known or measurable? This 'agnostic' (with respect to maternal phenotype) approach has the potential both to identify novel genetic variants of maternal effect and, via 'reverse epidemiology,' to identify novel maternal phenotypes for such effects.
For the purposes of this study, we accept the strict definition of a genetic maternal effect as the indirect effect of maternal genotype on offspring phenotype , as distinct from asymmetric transmission of parental alleles (e.g., mitochondrial inheritance ) or asymmetric expression of alleles in the offspring depending on parent of origin (e.g., imprinting ). Here we develop mixed variance-components models in the computer program SOLAR  to estimate maternal random effects on quantitative phenotypes, and use the best such models as null hypotheses for measured-genotype genome-wide association tests of single-nucleotide polymorphism (SNP) genotypes of individuals and their mothers.
Data include adult quantitative phenotypes and Affymetrix SNP genotypes provided in the Genetic Analysis Workshop (GAW) 16 Framingham Heart Study (FHS) data release (Problem 2). All authors of this study are 'approved users' of these data per the NHLBI Data Use Certification of April 2008. Analysis of these data was approved by the Institutional Review Board of the University of Texas Health Science Center, San Antonio.
Outlying phenotype measurements (more than four standard deviations from the mean) were removed from the lipid measures (10 for total cholesterol, 9 for high-density lipoprotein (HDL), 43 for triglyceride (TG)) on the assumption that these represented assay errors. The data were normal-quantile-transformed before analysis using the SOLAR "inormal" option to meet the distributional assumptions of the variance components and regression methods. The normal quantile ("inverse normal") transformation is robust to a range of departures from normality and also removed scale effects by standardizing the data. Transformations of this type are convenient for batch processing of multiple phenotypes (e.g., Peng et al. ).
Individuals with incomplete genotype data were given imputed genotype scores for the missing markers using the -- infer option in the computer program Merlin [9, 10] Merlin imputes an expected genotype score based on the probability of each possible genotype at a locus given information on marker allele frequency, adjacent markers, and pedigree relationships. We chose not to exclude any SNPs or individuals on the basis of number of incomplete genotypes (unless no genotypes were available at all), given the robustness of imputation from family data . Genotypes were similarly imputed for all SNPs for the non-genotyped implicit mothers of genotyped and phenotyped founders. These maternal genotypes entered the association analysis as properties of their offspring (see "Measured genotype analysis," below); the 'virtual' mothers did not enter the analysis otherwise.
with R = 1 if i = j (equivalent to the identity matrix in Eq. 1), ρ sib ∈ [0,1] if i, j are siblings or half-siblings, ρmo ∈ [-1,1] if i, j are mother and offspring, or 0 otherwise. Our modification from Bijma  was that twins were not treated differently than other siblings because dizygotic twins could not be distinguished in the de-identified FHS data. Our full mixed model also included the fixed effects of relevant covariates and the random effect of mitochondrial inheritance, σ2 mito ; the mitochondrial variance component is structured by a matrix whose elements are 1 if i, j belong to the same matriline or 0 otherwise, as described by Czerwinski et al. .
Measured genotype analysis was conducted for each polymorphic SNP by including its genotype score (the number of copies of the minor allele, range [0,2] with non-integral values for imputed genotypes) as a covariate in the mixed model . Unlike standard association analysis, we included the indirect effect of the mother's genotype in addition to the direct effect of the focal subject's genotype. These effects were tested separately, with an additional test of the mother's genotype conditional on that of the focal subject. The latter test was intended to account for the non-independence of maternal and offspring genotypes: reduction of evidence of maternal association in the conditional test would suggest that the unconditioned maternal effect represented a 'bleed-through' of the direct effect, while an increase in evidence would suggest that the locus affects the trait both directly and indirectly.
FHS cohorts/examination periods used in this study
FHS Cohort/Exam Period
Age in years [mean (SE)]
Original/Exam 4 (1954-1958)
Offspring/Exam 3 (1983-1987)
Generation 3/Exam 1 (2002-2005)
Parameter estimates [mean (SE)] for saturated random-effects model
Comparative evidence (log likelihoods), variance-components models
Polygenic + mitochondrial
Suggestive associations with height
Sorted by OMG
Sorted by MMG
Sorted by CMMG
Several recent studies have undertaken genome-wide association analysis of human height [20–22]. These studies have typically examined very large numbers of individuals (~10,000-25,000, with multi-stage designs), larger than the 6,775 individuals in the FHS cohort available for this study. These studies agree in finding numerous loci associated with height, as may be expected for a trait long assumed to be polygenic. Under these circumstances, it is not surprising that we did not replicate specific SNPs or locations identified in these larger studies. It should be noted, however, that we did find suggestive evidence of association (OMG) with a broad genomic region identified by Gudbjartsson et al. : 1q24-25 (Table 4). Our candidate genes in this region is MPZL1 (OMIM #604376), a protein tyrosine phosphatase involved in cell proliferation and differentiation. Our next four highest 'hits' were in the mucolipin2 gene MCOLN2 on 1p22.
Because the published genome-wide association study on height used unrelated individuals, none reported maternal effects. Interestingly, among our strongest maternal associations are repeated hits in two regions: the TRAPPC9 gene on chromosome 8 and an intergenic region on chromosome 18. The trafficking protein particle complex 9 gene TRAPPC9 (NIBP) plays a role in neuronal NF-κB signalling  but has not, to our knowledge, been associated with stature in any published study. The existence of repeated (albeit suggestive) associations in this gene makes it a candidate for further investigation of the effects of maternal genotype on height.
We have implemented combined random-effects, measured-genotype fixed effects approach for discovery of genetic variants contributing to the indirect effect of maternal genotype on offspring phenotype. We have identified two regions on chromosomes 2 and 8 - with suggestive association at two SNPs in each region - that may contribute to maternal effects on human height. The tools developed here should be of use for a variety of phenotypes and diseases for which an effect of maternal environment is known or suspected, including height, hypertension, birthweight, and the metabolic syndrome.
Body mass index
Conditional maternal measured genotype
Diastolic blood pressure
Framingham Heart Study
Genetic Analysis Workshop 16
Maternal measured genotype
Own measured genotype
Polygenic SBP: Systolic blood pressure
The Genetic Analysis Workshops are supported by NIH grant R01 GM031575 from the National Institute of General Medical Sciences. We thank the NHLBI and the Framingham Heart Study participants and investigators for making these data available. Our work in this study was supported by NIH grant R01 MH59490 (John Blangero, PI) and was performed in computing facilities funded in part by the AT&T Foundation.
This article has been published as part of BMC Proceedings Volume 3 Supplement 7, 2009: Genetic Analysis Workshop 16. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/3?issue=S7.
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.