- Open Access
X chromosome effects and their interactions with mitochondrial effects
BMC Genetics volume 6, Article number: S157 (2005)
We report a simple and rapid method for detecting additive genetic variance due to X-linked loci in the absence of marker data for this chromosome. We examined the interaction of this method with an established method for detecting mitochondrial linkage (another source of sex-asymmetric genetic covariance). When applied to data from the Collaborative Study on the Genetics of Alcoholism, this method found evidence of X-chromosomal linkage for one continuous trait (ntth1) and one discrete trait (SPENT). Evidence of mitochondrial contribution was found for one discrete trait (CRAVING) and three continuous traits (ln(CIGPKYR), ecb21, and tth1). Results for ntth1 suggest that methods that do not also allow for male-female heterogeneity in environmental variance may be overly conservative in detection of X-chromosomal effects.
Although the X chromosome is one of the largest in the human genome, genetic analysis software packages have been slow to incorporate routines for X-linked loci. This gap is closing: experimental two-point X-linkage analysis routines are now available in MERLIN  and SOLAR , and Ekstrøm  recently proposed an algorithm for multipoint linkage analysis of the X. Even prior to localizing quantitative trait loci within the X, however, it may be useful to determine the additive genetic effect of the X as a whole, both to allow the investigator to decide whether to pursue a more detailed analysis of the chromosome, and to improve the specificity of variance components models for autosomal linkage.
In Ekstrøm , in part to avoid assumptions about presence or absence of dosage compensation, both the linkage model and the model for the residual additive genetic effect of the X are estimated in terms of separate variance components for male-male, female-female, and male-female relative pairs. Here we offer a simplified mechanism for analyzing the 'X effect' as a single parameter (for ease of both computation and interpretation), even when marker data for this chromosome are unavailable. Our formulation derives from Bulmer's  dosage compensation model, which is applicable to the majority of X-chromosomal loci in humans and other eutherian mammals. Direct comparison of our simplified approach with that of Ekstrøm  is beyond the scope of this paper. However, based on a simulation study that does make this comparison by Kent et al. , we believe that the simplified method has equal or greater power to detect dosage compensated loci, and is conservative for non-lyonizing loci.
We used the method as a tool address the more general question of the interaction of sex-related effects in complex traits. Both X-linkage and mitochondrial linkage display asymmetric patterns of allele transmission from males and females, and either may be confounded with genotype × sex or environment × sex interaction effects on traits influenced by autosomal genes. A fully accurate specification of the variance components model for inheritance of a particular trait should distinguish among these possibilities. A companion study  examines possible mitochondrial contributions to a set of alcoholism-related phenotypes from the Collaborative Study on the Genetics of Alcoholism (COGA); here we extend this investigation to X-chromosomal effects.
The COGA data
We used behavioral and electrophysiological phenotypes in the subset of COGA data generously provided for the Genetic Analysis Workshop 14. The data collection process, including acquisition of written informed consent from all participants and approvals by Institutional Review Boards, are described more fully elsewhere .
Effect of the X chromosome
We have implemented an extension to the genetic analysis package SOLAR  to estimate the additive genetic effect of X-linked loci. For humans and many other placental mammals in which most X-linked loci are subject to random inactivation in females, it should be a reasonable first approximation to employ a simple dosage-compensation model . The additive contribution of an X-linked locus is simply the effect of the single allele in the male, or the single active allele in a homozygous female. As an approximation, the effect of the locus in heterozygous females is just the average of the allele effects, because each female is a mosaic of cells in which one or the other X is inactivated at random. Under these assumptions, the mean X effect is the same in both sexes, while the variance due to the X in males is twice that in females. This linear relationship between the variances allows us to express the additive effect of the X in one sex in terms of the other. For example, in terms of the female variance σ2Xf,
Ω = 2Φ σ2g + 2Ψ σ2Xf + Ifσ2ef + Imσ2em, (1)
where Ω is the matrix of phenotypic correlations between pairs of individuals, Φ is the matrix of autosomal kinship coefficients, and Ψ (using the notation of Ekstrøm ) is a new matrix of X-chromosome-specific coefficients of relationship for male-male, female-female, and male-female pairs (Table 1). Like the elements of Φ, the elements of Ψ can be calculated recursively from pedigree data, as by the algorithm in MINX (MERLIN in X [1, 8]).
Inclusion of sex-specific identity matrices and random environmental variances (Ifσ2ef and Imσ2em) in Equation 1 is desirable because the different σ2Xf and σ2Xm may appear as different environmental variances in the restricted model where σ2Xf = σ2Xm = 0, and also to account for genotype × sex or environment × sex interactions, if any. An important exception is that discrete traits (e.g., CRAVING) are analyzed using heritability parameterization with the phenotypic variance in both sexes fixed at 1 (because the variance-components model for discrete traits assumes an underlying liability with distribution N(0, 1)).
The significance test is straightforward, testing the ratio of the likelihood of the model with the single added parameter σ2Xf to that of the restricted model. Because σ2Xf is tested on its lower bound, the likelihood ratio statistic is distributed as a 1/2:1/2 mixture of χ2 with one degree of freedom, and a point mass at zero .
We used the mitochondrial identity-by-descent (IBD) calculation implemented in SOLAR  (and see discussion in Lease et al. ). The additive genetic variance due to the mitochondrion can be added to the basic polygenic model as
Ω = 2Φ σ2g + M σ2mt + Ifσ2ef + Imσ2em, (2)
where M is the mitochondrial IBD matrix. For comparison to the X effects, we retain separate environmental variance terms for males and females (except for the discrete traits). The significance of the model that includes the mitochondrial variance is tested against the restricted model (σ2mt = 0) by a likelihood ratio test.
Results and Discussion
In the following, evidence of mitochondrial or X effect is listed as significant at p < 0.05, without correction for multiple testing (the traits are genetically and phenotypically correlated  and measured in the same individuals). Note, however, that a modified Bonferroni correction that accommodates these genetic correlations, described by Lease et al. , yields a more stringent experiment wise type I error rate = 0.01.
Mitochondrial and X effects in discrete traits
Phenotypes include 12 discrete behavioral traits related to alcohol consumption . Two traits (BINGE and MORNING) did not show significant heritability in the polygenic model and were not considered further. Of the remaining traits (SMOKER, DESIRE, CRAVING, SPENT, NARROW, GAVEUP, BLACKOUT, WITHDRAWAL, HEALTH, and PSYCHO), data are shown only for those whose model likelihoods were increased by the addition of mitochondrial or X-effect terms (DESIRE, which showed a marginally significant X effect, is also shown). The sequence of analysis was: estimation of a polygenic model (POL); addition of a mitochondrial heritability term (MIT); addition of a heritability term for the additive genetic effect of the X (XFX). At each step the more restricted model was retained unless the model with the added parameter had a likelihood that was significantly greater at p < 0.05. Table 2 reports the p-values for the successive models and point estimates of the parameters retained in the final model for each trait. Trait SPENT showed a significant effect (p = 0.034) of the X chromosome, but the size of this effect (6% of total variance) is small. CRAVING showed a significant mitochondrial effect (p = 0.005), as reported , but no significant effect of the X chromosome. In each case the new parameter accounts for the entire additive genetic effect (i.e., h2g = 0) in the final model.
Mitochondrial and X effects in continuous traits
We examined two behavioral traits, MAXDRINKS and CIGPKYR. Both measures were leptokurtic and were loge-transformed: κ before(after) transformation: MAXDRINKS 9.18 (0.06); CIGPKYR 9.69 (-1.46). We also examined 13 electrophysiology traits (ecb21, ntth1-4, ttdt1-4, ttth1-4). All showed a significant additive genetic component; data are shown only for those traits whose model likelihoods were increased by addition of more parameters. The sequence of analysis was: estimation of a (model POL); estimation of a polygenic model with separate environmental variance terms for the two sexes (EMF); addition of a mitochondrial variance term (MIT); addition of a variance term for the additive genetic effect of the X (XFX). At each step the more restricted model was retained unless the model with the added parameter had a significantly greater likelihood at p < 0.05. Table 3 reports the p-values for the successive models and the estimates of the parameters retained in the final model for each trait. Two traits (LNMXDK and ttdt2) had significantly different environmental variances in males and females, although neither showed significant mitochondrial or X effects. Three traits (LNCIGP, ecb21, and ttth1) showed evidence for mitochondrial effect. One trait (ntth1) showed a significant effect (p = 0.021) of the X.
Interaction of sex-related effects
The initial screen of traits followed a step-wise procedure in which each additional parameter was retained in subsequent analyses only if it significantly increased the likelihood of the model. We were curious if estimating other combinations of parameters might alter the conclusions from this analysis. Table 4 presents results from a series of nested models for continuous traits LNCIGP and ntth1. (In the following, models A, B, etc. for each trait refer to their listing in Table 4.)
LNCIGP showed a significant mitochondrial effect in the initial screen (Table 3). This result is robust; the effect is detectable whether male and female environmental variances are constrained to be the same (model D vs. A, p = 0.037) or allowed to diverge (F vs. B, p = 0.030). In model E (σ2mt = 0) most of the variance actually due to the mitochondrion appears as residual additive genetic variance, with only minor misallocation of variance to the X.
A different pattern appears for ntth1, which showed a significant X effect in the initial screen. Interestingly, a significant X effect is detectable only if male and female environmental variances are allowed to differ (model E vs. B, p = 0.02): the addition of σ2Xf to the polygenic model barely increases the model likelihood (model H vs. A, p = 0.44). Different estimates of male and female environmental variance are retained in models that include the X effect, which may suggest genotype × sex interactions. Exclusion of σ2Xf does not yield a falsely significant estimate of σ2mt (model F vs. model B, p = 0.06).
In this study we have applied a new method for detecting additive genetic variance due to the X chromosome to a representative set of real data. (Simulation studies to validate the method will be reported elsewhere.) We find that real mitochondrial effects (estimated by an established method) are not confounded with X effects. However, it appears that a true X effect can appear as heterogeneous environmental variance in males and females; indeed, in one case (ntth1), detection of an X effect depended on modelling this heterogeneity. This suggests that tests of X-chromosomal effect that constrain males and females to have equal environmental variances (as, for example, analyses of discrete traits) may be overly conservative.
Collaborative Study on the Genetics of Alcoholism
Model estimating separate environmental variances by sex
Model estimating mitochondrial variance component
Model estimating X-linked variance component
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin-rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002, 30: 97-101. 10.1038/ng786.
Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.
Ekstrøm CT: Multipoint linkage analysis of quantitative traits on sex-chromosomes. Genet Epidemiol. 2004, 26: 218-230. 10.1002/gepi.10310.
Bulmer MG: The Mathematical Theory of Quantitative Genetics. 1985, New York: Oxford University Press
Kent JW, Dyer TD, Blangero J: Estimating the additive genetic effect of the X chromosome. Genet Epidemiol.
Lease LR, Winnier DA, Williams JT, Dyer TD, Almasy L, Mahaney MC: Mitochondrial genetic effects on latent class variables associated with susceptibility to alcoholism. BMC Genet. 2005, 6 (Suppl 1): S158-10.1186/1471-2156-6-S1-S158.
Begleiter H, Reich T, Nurnberger J, Li TK, Conneally PM, Edenberg H, Crowe R, Kuperman S, Schuckit M, Bloom F, Hesselbrock V, Porjesz B, Cloninger CR, Rice J, Goate A: Description of the Genetics Analysis Workshop 11 Collaborative Study on the Genetics of Alcoholism. Genet Epidemiol. 1999, 17 (Suppl 1): S25-S30.
MINX: Chromosome X Analyses. The University of Michigan Center for Statistical Genetics MERLIN Reference Sheet, [http://www.sph.umich.edu/csg/abecasis/Merlin/reference.html]
Self SG, Liang K-Y: Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J Am Stat Assoc. 1987, 82: 605-610. 10.2307/2289471.
Czerwinski SA, Williams JT, Demerath EW, Towne B, Siervogel RM, Blangero J: Does accounting for mitochondrial genetic variation improve the fit of genetic models?. Genet Epidemiol. 2001, 21 (Suppl 1): S779-S782.
This work was supported in part by NIH/NIMH grant MH059490-06 (JB, principal investigator). We appreciate comments on the first draft by two anonymous reviewers.
JWK, TDD, and JB developed the methodology for analyzing the effect of the X chromosome. JWK conceived the study, performed the X-effect analyses, and drafted the manuscript. LRL and MCM performed the mitochondrial analyses. MCM and LA provided critical reviews of the methodology. TDD and LA coordinated preparation of the COGA data for investigators at SFBR. All authors have read and approved this manuscript.
About this article
Cite this article
Kent, J.W., Lease, L.R., Mahaney, M.C. et al. X chromosome effects and their interactions with mitochondrial effects. BMC Genet 6, S157 (2005) doi:10.1186/1471-2156-6-S1-S157
- Additive Genetic Variance
- Dosage Compensation
- Restricted Model
- Additive Genetic Effect
- Continuous Trait