 Proceedings
 Open Access
 Published:
X chromosome effects and their interactions with mitochondrial effects
BMC Genetics volume 6, Article number: S157 (2005)
Abstract
We report a simple and rapid method for detecting additive genetic variance due to Xlinked 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 sexasymmetric genetic covariance). When applied to data from the Collaborative Study on the Genetics of Alcoholism, this method found evidence of Xchromosomal 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 malefemale heterogeneity in environmental variance may be overly conservative in detection of Xchromosomal effects.
Background
Although the X chromosome is one of the largest in the human genome, genetic analysis software packages have been slow to incorporate routines for Xlinked loci. This gap is closing: experimental twopoint Xlinkage analysis routines are now available in MERLIN [1] and SOLAR [2], and Ekstrøm [3] 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 [3], 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 malemale, femalefemale, and malefemale 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 [4] dosage compensation model, which is applicable to the majority of Xchromosomal loci in humans and other eutherian mammals. Direct comparison of our simplified approach with that of Ekstrøm [3] is beyond the scope of this paper. However, based on a simulation study that does make this comparison by Kent et al. [5], we believe that the simplified method has equal or greater power to detect dosage compensated loci, and is conservative for nonlyonizing loci.
We used the method as a tool address the more general question of the interaction of sexrelated effects in complex traits. Both Xlinkage 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 [6] examines possible mitochondrial contributions to a set of alcoholismrelated phenotypes from the Collaborative Study on the Genetics of Alcoholism (COGA); here we extend this investigation to Xchromosomal effects.
Methods
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 [7].
Effect of the X chromosome
We have implemented an extension to the genetic analysis package SOLAR [2] to estimate the additive genetic effect of Xlinked loci. For humans and many other placental mammals in which most Xlinked loci are subject to random inactivation in females, it should be a reasonable first approximation to employ a simple dosagecompensation model [4]. The additive contribution of an Xlinked 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 σ^{2}_{Xf},
Ω = 2Φ σ^{2}_{g} + 2Ψ σ^{2}_{Xf} + I_{f}σ^{2}_{ef} + I_{m}σ^{2}_{em}, (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 [3]) is a new matrix of Xchromosomespecific coefficients of relationship for malemale, femalefemale, and malefemale 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 sexspecific identity matrices and random environmental variances (I_{f}σ^{2}_{ef} and I_{m}σ^{2}_{em}) in Equation 1 is desirable because the different σ^{2}_{Xf} and σ^{2}_{Xm} may appear as different environmental variances in the restricted model where σ^{2}_{Xf} = σ^{2}_{Xm} = 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 variancecomponents 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 σ^{2}_{Xf} to that of the restricted model. Because σ^{2}_{Xf} 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 [9].
Mitochondrial effect
We used the mitochondrial identitybydescent (IBD) calculation implemented in SOLAR [10] (and see discussion in Lease et al. [6]). The additive genetic variance due to the mitochondrion can be added to the basic polygenic model as
Ω = 2Φ σ^{2}_{g} + M σ^{2}_{mt} + I_{f}σ^{2}_{ef} + I_{m}σ^{2}_{em}, (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 (σ^{2}_{mt} = 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 [6] and measured in the same individuals). Note, however, that a modified Bonferroni correction that accommodates these genetic correlations, described by Lease et al. [6], 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 [7]. 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 Xeffect 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 pvalues 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 [6], but no significant effect of the X chromosome. In each case the new parameter accounts for the entire additive genetic effect (i.e., h^{2}_{g} = 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 log_{e}transformed: κ before(after) transformation: MAXDRINKS 9.18 (0.06); CIGPKYR 9.69 (1.46). We also examined 13 electrophysiology traits (ecb21, ntth14, ttdt14, ttth14). 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 pvalues 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 sexrelated effects
The initial screen of traits followed a stepwise 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 (σ^{2}_{mt} = 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 σ^{2}_{Xf} 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 σ^{2}_{Xf} does not yield a falsely significant estimate of σ^{2}_{mt} (model F vs. model B, p = 0.06).
Conclusion
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 Xchromosomal effect that constrain males and females to have equal environmental variances (as, for example, analyses of discrete traits) may be overly conservative.
Abbreviations
 COGA:

Collaborative Study on the Genetics of Alcoholism
 EMF:

Model estimating separate environmental variances by sex
 IBD:

Identitybydescent
 MIT:

Model estimating mitochondrial variance component
 POL:

polygenic model
 XFX:

Model estimating Xlinked variance component
References
 1.
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlinrapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002, 30: 97101. 10.1038/ng786.
 2.
Almasy L, Blangero J: Multipoint quantitativetrait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 11981211. 10.1086/301844.
 3.
Ekstrøm CT: Multipoint linkage analysis of quantitative traits on sexchromosomes. Genet Epidemiol. 2004, 26: 218230. 10.1002/gepi.10310.
 4.
Bulmer MG: The Mathematical Theory of Quantitative Genetics. 1985, New York: Oxford University Press
 5.
Kent JW, Dyer TD, Blangero J: Estimating the additive genetic effect of the X chromosome. Genet Epidemiol.
 6.
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): S15810.1186/147121566S1S158.
 7.
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): S25S30.
 8.
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]
 9.
Self SG, Liang KY: Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J Am Stat Assoc. 1987, 82: 605610. 10.2307/2289471.
 10.
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): S779S782.
Acknowledgements
This work was supported in part by NIH/NIMH grant MH05949006 (JB, principal investigator). We appreciate comments on the first draft by two anonymous reviewers.
Author information
Additional information
Authors' contributions
JWK, TDD, and JB developed the methodology for analyzing the effect of the X chromosome. JWK conceived the study, performed the Xeffect 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.
Rights and permissions
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/147121566S1S157
Published
DOI
Keywords
 Additive Genetic Variance
 Dosage Compensation
 Restricted Model
 Additive Genetic Effect
 Continuous Trait