Email updates

Keep up to date with the latest news and content from Molecular Pain and BioMed Central.

Open Access Highly Accessed Methodology

A computational model for sex-specific genetic architecture of complex traits in humans: Implications for mapping pain sensitivity

Chenguang Wang1, Yun Cheng1, Tian Liu1, Qin Li1, Roger B Fillingim2, Margaret R Wallace3, Roland Staud3, Lee Kaplan3 and Rongling Wu1*

Author Affiliations

1 Department of Statistics, University of Florida, Gainesville, FL 32611 USA

2 Department of Community Dentistry and Behavioral Science, University of Florida, Gainesville, FL 32611, USA

3 Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL 32611, USA

For all author emails, please log on.

Molecular Pain 2008, 4:13  doi:10.1186/1744-8069-4-13

The electronic version of this article is the complete one and can be found online at: http://www.molecularpain.com/content/4/1/13


Received:29 November 2007
Accepted:16 April 2008
Published:16 April 2008

© 2008 Wang et al; licensee 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.

Abstract

Understanding differences in the genetic architecture of complex traits between the two sexes has significant implications for evolutionary studies and clinical diagnosis. However, our knowledge about sex-specific genetic architecture is limited largely because of a lack of analytical models that can detect and quantify the effects of sex on the complexity of quantitative genetic variation. Here, we derived a statistical model for mapping DNA sequence variants that contribute to sex-specific differences in allele frequencies, linkage disequilibria, and additive and dominance genetic effects due to haplotype diversity. This model allows a genome-wide search for functional haplotypes and the estimation and test of haplotype by sex interactions and sex-specific heritability. The model, validated by simulation studies, was used to detect sex-specific functional haplotypes that encode a pain sensitivity trait in humans. The model could have important implications for mapping complex trait genes and studying the detailed genetic architecture of sex-specific differences.

Background

Differences in males and females (sexual dimorphism) is ubiquitous in many biological aspects [1-3]. In humans, sexually dimorphic traits include those from morphological shapes and body size to brain development to disease susceptibility [4,5]. Substantial differences are also observed in sensitivities to pain and pain-killing drugs, and susceptibility to developing chronic pain between men and women [6-8]. All these sex-specific differences are due to varying expression of genes on the X/Y chromosome and autosomes, thought to result from differences in cellular and hormonal environments between the two sexes [9]. A growing body of research has been conducted to elucidate the genetic control of sexual dimorphism in various complex phenotypes by gene mapping approaches [4,5,10,11]. Despite these efforts, however, little is known about the genetic architecture underlying sex-related variation in a quantitative trait.

Since sex is easily determined, the effects of sex on morphological, developmental and pathological traits can be directly observed. However, characterizing the impacts of sex on the genetic architecture of these traits has been challenged by a lack of powerful statistical approaches. The motivation of this article is to develop a statistical and computational model that can systematically search for sex-specific genes contributing to quantitative variation and formulate testable hypotheses regarding the interplay between sexes and gene expression. Our model is principally different from those used in many previous studies that are aimed to detect sex-specific quantitative trait loci (QTLs) based on linkage or linkage disequilibrium analysis [2,4,5,12,13]. Our model will be founded on the statistical framework constructed by Liu et al. [14] to detect the effects and diversity of haplotypes constructed by single nucleotide polymorphisms (SNPs) that are genotyped at candidate genes or genome-wide [15]. Our model has been generalized to allow the test of sex differences in haplotype frequencies, allele frequencies and linkage disequilibria between different SNPs as well as additive and dominant effects of haplotypes on complex traits. It has power to identify sex-specific DNA sequence variants that encode complex phenotypes in men and women.

Model

Notation

Suppose there is a diversity of haplotypes constructed by two SNPs each with two alleles designated as 1 and 0. Let p and q be the 1-allele frequencies for the first and second SNP, respectively. Thus, the 0-allele frequencies at these two different SNPs will be 1 - p and 1 - q. The two SNPs that are segregating in a natural human population form four haplotypes, [11], [10], [01], and [00], whose frequencies are constructed by allele frequencies and linkage disequilibrium (D) between the two SNPs, i.e.,

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M1">View MathML</a>

(1)

The parameters contained in equation (1) can be used to describe some important aspects of the genetic structure and diversity of a natural population. Thus, differences in genetic architecture between the two different sexes can be characterized by these sex-specific parameters. Because it is easy to derive the closed forms for estimating haplotype frequencies [14], we will estimate the linkage disequilibrium from the estimated haplotype frequencies.

Let <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M2">View MathML</a> and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M3">View MathML</a> be the vectors of haplotype frequencies among males and females, respectively. All the genotypes for the two SNPs are consistent with diplotypes, except for the double heterozygote, 10/10, that belongs to a diplotype of either [11] [00] or [10] [01] (Table 1). Assuming that the population is at Hardy-Weinberg equilibrium, the frequency of a diplotype is expressed as the product of the frequencies of the two haplotypes that construct it. Table 1 characterizes the differences in diplotype frequencies between the males and females.

Table 1. Diplotypes and their frequencies for each of nine genotypes at two SNPs, and composite diplotypes for one assumed sex-specific risk haplotype chosen from four possible haplotypes.

If haplotypes triggers an effect on a quantitative trait, this means that at least one haplotype performs differently from the rest of the haplotypes. Without loss of generality, let haplotype [11] be such a distinct haplotype, called risk haplotype, designated as A. All the other non-risk haplotypes, [10], [01] and [00], are collectively expressed as <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M31">View MathML</a>. The risk and non-risk haplotypes form three composite diplotypes AA (symbolized as 2), <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M10">View MathML</a> (symbolized as 1) and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M6">View MathML</a> (symbolized as 0). The genotypic values of the three composite diplotypes may be different between the two sexes, arrayed in <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M32">View MathML</a> for the males and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M33">View MathML</a> for the females, respectively. Let (aM, dM) and (aF, dF) be the additive and dominance genetic effects due to the risk and non-risk haplotypes in males and females, respectively.

Likelihoods

Assume that a total of n subjects (including nM males and nF females) sampled from the population are phenotyped for a quantitative trait. In each sex, there are nine possible genotypes for the two SNPs, each genotype with an observed number generally expressed as <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M34">View MathML</a> for the males and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M35">View MathML</a> for the females (r1 r'1, r2 r'2, r3 r'3 = 1,0).

Many physiological traits scale with body weight (W) according to a power function with a certain allometric exponent. Thus, we implement this allometric scaling law to describe the phenotypic value of a trait for subject i within male or female subpopulations in terms of the haplotypes considered as

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M36">View MathML</a>

(2)

where (αM, βM) or (αF, βF) are body weight-related allometric coefficients, <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M37">View MathML</a> or <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M38">View MathML</a> are the indicator variables associated with the additive and dominance effects, respectively, and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M39">View MathML</a> or <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M40">View MathML</a> is the residual error, normally distributed as <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M41">View MathML</a> or <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M42">View MathML</a>. The genotypic values of composite diplotypes and variance are arrayed by a quantitative genetic parameter vector <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M43">View MathML</a> for the males and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M44">View MathML</a> for the females, respectively.

The log-likelihood of haplotype frequencies, genotypic values of composite diplotypes and residual variances given sex-specific phenotypic (yM, yF) and SNP data (SM, SF) is factorized into two parts, expressed as

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M45">View MathML</a>

(3)

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M46">View MathML</a>

(4)

where equation (4) is derived from equation (3) because the males and females are assumed to be independent, and

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M47">View MathML</a>

(5)

where <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M48">View MathML</a> is a normal distribution density function of composite diplotype j (j = 2, 1, 0) for sex k, and

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M49">View MathML</a>

(6)

is the relative proportion of diplotype [11] [00] within the double heterozygote for sex k.

It can be seen from equation (3) or (4) that maximizing <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M50">View MathML</a> is equivalent to maximizing log <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M51">View MathML</a> and log <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M52">View MathML</a> individually in equation (5).

The EM algorithm

A closed-form solution for the EM algorithm [14] has been derived to estimate the unknown parameters that maximize the two sex-specific likelihoods of (5). The estimates of sex-specific haplotype frequencies are based on the log-likelihood function <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M51">View MathML</a>, whereas the estimates of sex-specific genotypic values of composite diplotypes and the residual variance are based on the log-likelihood function <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M53">View MathML</a>. These two different types of parameters can be estimated using a two-stage hierarchical EM algorithm (see ref. [14] for a detailed implementation).

Model selection

According to equation (5), the summed likelihood across the sexes, <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M54">View MathML</a>, is formulated by assuming that haplotype [11] is a risk haplotype. However, a real risk haplotype is unknown from raw data (yk, Sk). An additional step for choosing the most likely risk haplotype should be implemented. The simplest way to do so is to calculate the likelihood values by assuming that any one of the four haplotypes can be a risk haplotype (Table 1). Thus, we obtain four possible likelihood values as follows:

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M55">View MathML</a>

The largest likelihood value calculated is thought to correspond to the most likely risk haplotype. Under an optimal risk haplotype, we estimate sex-specific quantitative genetic parameters <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M56">View MathML</a> and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M57">View MathML</a>.

Hypothesis tests

The genetic architecture of a quantitative trait is characterized by population (including haplotype frequencies, allele frequencies, and linkage disequilibria) and quantitative genetic parameters (including haplotype effects and mode of inheritance for haplotypes). The model proposed provides a meaningful way for estimating the genetic architecture of a trait and further testing sex-specific differences in genetic control.

After haplotype frequencies are estimated, allele frequencies and linkage disequilibrium between the two SNPs with each sex can be calculated as

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M58">View MathML</a>

(7)

The influence of haplotypes on a quantitative trait is quantified in terms of the additive (a) and dominant genetic effects (d), and the mode of inheritance (ρ), which are estimated for each sex. Each of these population and quantitative genetic parameters can be tested when appropriate hypotheses are formulated.

Overall genetic control

Haplotype effects on the trait, i.e., the existence of functional haplotypes, in both male and female populations can be tested using the following hypotheses expressed as

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M59">View MathML</a>

(8)

The log-likelihood ratio test statistic (LR) under these two hypotheses can be similarly calculated,

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M60">View MathML</a>

(9)

where the L0 and L1 are the plug-in likelihood values under the null and alternative hypotheses of (8), respectively. Although the critical threshold for determining the existence of a functional haplotype can be based on empirical permutation tests, the LR may asymptotically follow a χ2 distribution with four degrees of freedom, so that the threshold can be obtained from the χ2distribution table.

Sex-specific population genetic architecture

The male and female populations may be different in terms of population genetic parameters. Such sex-specific differences can be tested by formulating the following hypotheses

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M61">View MathML</a>

(10)

for allele frequency at SNP 1,

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M62">View MathML</a>

(11)

for allele frequency at SNP 2, and

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M63">View MathML</a>

(12)

for the linkage disequilibrium between the two SNPs.

For each of the hypotheses (10)–(12), the LR values are calculated, which are each thought to asymptotically follow a χ2-distribution with one degree of freedom. Sex-specific differences in overall population genetic architecture can be tested with the null hypothesis H0: pM = pF, qM = qF, and DM = DF, with the corresponding LR value to be χ 2-distributed with three degrees of freedom.

Sex-specific quantitative genetic architecture

Sex-specific differences in overall quantitative genetic architecture can be tested by formulating the hypotheses

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M64">View MathML</a>

(13)

The LR value calculated under the null and alternative hypotheses is suggested to follow a χ2-distribution with two degrees of freedom. The rejection of the null hypothesis implies that the effects of the same haplotype are different between the two sexes. If there exists a sex-specific difference, the next step is to test whether this difference is due to the additive or dominant genetic effects, or both.

Sex-specific risk haplotypes

In the preceding sections, the same risk haplotype was assumed between the male and female populations. It is possible that the two sexes have different risk haplotypes. Let <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M65">View MathML</a> (jm = 2, 1, 0) and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M66">View MathML</a> (jf = 2, 1, 0) be the genotypic values of composite diplotypes for the males and females constructed by a sex-specific rick haplotype. By reformulating the likelihood log <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M52">View MathML</a> of equation (5) based on sex-specific composite diplotypes, these genotypic values can be estimated with the EM algorithm. A best combination of risk haplotypes between the two sexes can be determined from the AIC values.

Multi-locus haplotyping

Three-SNP model

Consider three associated SNPs, S1, S2, and S3, each with two alleles denoted by 1 and 0. Let p, q and r, and D12, D13, D23 and D123 be the 1-allele frequencies for the three SNPs, and the linkage disequilibria between SNPs 1 and 2, SNPs 1 and 3, SNPs 2 and 3 and among the three SNPs, respectively. Eight haplotypes, [111], [110], [101], [100], [011], [010], [001] and [000], formed by these three SNPs, have sex-specific frequencies arrayed in <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M67">View MathML</a> for sex k. Each of these haplotype frequencies is constructed by allele frequencies at different SNPs and their linkage disequilibria of different orders, expressed as

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M68">View MathML</a>

(14)

for sex k.

Sex-specific population genetic architecture can be tested by comparing the differences in allele frequencies (pk, qk, rk) and linkage disequilibria of different orders <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M69">View MathML</a> between the males and females.

In a natural population, there are 27 genotypes for the three SNPs. The frequency of each genotype is expressed in terms of haplotype frequencies. Some genotypes are consistent with diplotypes, whereas the others that are heterozygous at two or more SNPs are not. Each double heterozygote contains two different diplotypes. One triple heterozygote, i.e., 10/10/10, contains four different diplotypes, [111] [000] (in a probability of <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M70">View MathML</a>), [110] [001] (in a probability of <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M71">View MathML</a>), [101] [010] (in a probability of <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M72">View MathML</a>) and [100] [011] (in a probability of 2p100p011). The relative frequencies of different diplotypes for this double or triple heterozygote are a function of haplotype frequencies (Supporting Information Table 1). The integrative EM algorithm can be employed to estimate the MLEs of haplotype frequencies. A general formula for estimating haplotype frequencies can be derived.

By assuming [111] as a risk haplotype (labeled by A) and all the others as non-risk haplotypes (labelled by <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M31">View MathML</a>), the formulation of genotypic values for three composite diplotypes, μ2 for AA, μ1 for <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M10">View MathML</a> and μ0 for <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M6">View MathML</a> can be derived. Similar procedures described for the two-SNP model can be obtained to estimate and test sex-specific additive and dominance genetic effects when a haplotype contains three SNPs.

L-SNP model

It is possible that the two- and three-SNP models are too simple to characterize genetic variants for quantitative variation. We can develop a model that includes an arbitrary number of SNPs whose sequences are associated with the phenotypic variation. A key issue for the multi-SNP sequencing model is how to distinguish among 2ℓ-1 different diplotypes for the same genotype heterozygous at ℓ loci. The relative frequencies of these diplotypes can be expressed in terms of haplotype frequencies.

Consider a a functional haplotype that contains L SNPs among which there exist linkage disequilibria of different orders. The two alleles, 1 and 0, at each of these SNPs are symbolized by r1, ..., rL, respectively. Let <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M73">View MathML</a> be the allele frequencies for these different SNPs within sex k. A haplotype frequency, denoted as <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M74">View MathML</a>, is decomposed into the following components:

<a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M75">View MathML</a>

where Dk's are the linkage disequilibria of different orders among particular SNPs for sex k.

Sex-specific difference in terms of allele frequencies and linkage disequilibria between different SNPs as well as haplotype additive and dominance effects can be tested by formulating the corresponding hypotheses.

Results

Pain genetics study

The model proposed was used to detect differences in the genetic architecture of pain sensitivity between men and women. Genetic and phenotypic data were from a pain genetics project in which 237 subjects (including 143 men and 94 women) from five different races were sampled for six SNPs at three candidate genes. As a demonstration of the utilization of the model, we will focus on two SNPs, OPRDT80G (with two alleles T and G) and OPRDT921C (with two alleles T and C), at the delta opioid receptor. Pain testing procedures followed Fillingim et al. [16]. The phenotypic values of traits were subtracted by the means for each race to remove the effect due to races.

These two SNPs construct four haplotypes, [TC], [TT], [GC], and [GT], which yield 10 diplotypes, [TC] [TC], [TC] [TT], [TT] [TT], [TC] [GC], [TC] [GT], [TT] [GC], [TT] [GT], [GC] [GC], [GC] [GT], and [GT] [GT] and nine genotypes, TT/CC, TT/CT, TT/TT, TG/CC, TG/CT, TG/TT, GG/CC, GG/CT and GG/TT. Based on the observed numbers of each genotype in the male and female populations, we estimated sex-specific haplotype frequencies (Table 1). The pattern of haplotype distribution is consistent between the two sexes, with haplotypes [TC] and [TT] jointly occupying a majority proportion in the populations. Haplotype [GT] is very rare, with the frequency close to zero. SNP OPRDT80G has a low heterozygosity because the frequency of its commoner allele (T) is closer to 0.90, whereas there is a high heterozygosity for SNP OPRDT921C in terms of its averaged allele frequencies. The two SNPs are highly significantly associated at p = 3.41 × 10-5 for males and p = 9.63 × 10-5 for females, with a normalized linkage disequilibrium of D' = 1.00, because alleles T from OPRDT80G and T from OPRDT921C as well alleles G from OPRDT80G and C from OPRDT921C tend to form the same haplotypes more frequently than at random. There is no sex-specific difference in allele frequencies at the two SNPs and their linkage disequilibrium.

By assuming that one of the haplotypes is a risk haplotype, we estimated the effects of each haplotype on a pain sensitivity trait, assessed with a baseline pressure pain threshold measured at the ulna, in the pooled male and female population. The likelihoods of haplotype [TC], [TT] and [GC] as a risk haplotype are -594.8, -594.5, and -596.1, and thus the most likely risk haplotype is [TT]. The genotypic values of composite diplotypes constructed by this risk haplotype and its non-risk haplotype counterpart were estimated and compared between different sexes. In both males and female, the three composite diplotypes do not display significant genetic differences in the pain trait studied, but the directions of the additive and dominance effects are different between the two sexes (Table 2). In males, the non-risk haplotype tends to increase pressure pain thresholds, and it is overdominant to the risk haplotype, leading to increased pressure pain thresholds at a marginal significance level (p = 0.058) (Fig. 1). By contrast, in females, the non-risk haplotype tends to reduce pressure pain thresholds, and it also tends to be overdominant to the risk haplotype by reducing pressure pain thresholds. These discrepancies in both effect size and direction (Fig. 2) make the overall quantitative genetic architecture of the pain sensitivity trait significantly different between the two sexes (p = 1.49 × 10-7) (Table 2). Although the additive genetic effect displays a gene by sex interaction at the p = 0.03 significance level, a gene by sex interaction for the dominance effect is highly significant at p = 6.16 × 10-7. No significant difference was observed in inheritance mode between males and females.

Table 2. The estimates and tests of population genetic structure for two SNPs, OPRDT80G (with alleles T and G) and OPRDT921C (with two alleles C and T), and quantitative genetic effects of haplotypes constructed by these two SNPs on baseline pressure pain thresholds measured at the ulna in males and females.

thumbnailFigure 1. Different genotypic values of baseline pressure pain thresholds measured at the ulna for composite diplotypes, AA, <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M10">View MathML</a>, and <a onClick="popup('http://www.molecularpain.com/content/4/1/13/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.molecularpain.com/content/4/1/13/mathml/M6">View MathML</a>, constructed by risk haplotype [TT] and non-risk haplotype in males and females. The origin indicates the mean of the genotypic values between the two homozygotes.

thumbnailFigure 2. Different additive (a) and dominant genetic effects (d) of haplotypes on baseline pressure pain thresholds measured at the ulna in males and females.

Monte carlo simulation

Simulation studies were performed to test the statistical properties of the model proposed. Given a certain sample size (n), we simulated two SNPs by assuming different allele frequencies and linkage disequilibria between two sexes. The hypothesized allele frequencies at the two SNPs are pM = 0.5, qM = 0.6 and DM = 0.1 for males and pF = 0.8, qF = 0.9 and DF = 0.06 for females. By postulating one of the four haplotypes constructed by the two SNPs as a risk haplotype, we calculated the genetic variance among three composite diplotypes using the additive effect a = 0.6 and dominance effect d = 0.8, from which the residual variance was calculated when the heritability (H2) of a trait is given. The phenotypic values of the trait were simulated by assuming that they follow a normal distribution under four different simulation designs, (1) n = 100 and H2 = 0.1, (2) n = 400 and H2 = 0.1, (3) n = 100 and H2 = 0.4, and (4) n = 400 and H2 = 0.4.

The new model was used to analyze the simulated SNP and phenotypic data with the results tabulated in Table 3. Population genetic parameters including allele frequencies (p and q) and linkage disequilibrium (D) can well be estimated, with increasing precision when the sample size increases from 100 to 400. The power to detect the given sex-specific differences in allele frequencies and linkage disequilibrium is as high as 0.95 even with a modest sample size (100). Although the traditional model that does not implement sex-specific differences can provide precise estimates of these parameters, the estimates are generally between the true values of males and females.

Table 3. The MLEs of population and quantitative genetic parameters and the standard errors of the estimates obtained by the new model and the power to detect sex-specific differences under different simulation designs. Parameter estimates by a conventional model are also given.

The estimation of quantitative genetic parameters including the additive (a) and dominance effects (d) needs the determination of an optimal risk haplotype. When all possible risk haplotypes were assumed for the simulated data, we found that the true risk haplotype gave the largest likelihood among the four possible cases. In general, quantitative genetic parameters can well be estimated, but the estimation precision increases dramatically with sample size and heritability. Although the additive effect can be obtained with reasonable precision at a modest heritability (0.1) with a modest sample size (100), the precise estimation of the dominance effect relies upon a larger heritability and sample size. Also, the given difference in the additive effect between two sexes can be detected with great power, even when both the sample size and heritability are small. But the same size of sex-specific difference in the dominance effect can be detected with the same power only when the sample size is 400 and heritability is 0.4 (Table 3). The traditional model gave biased estimates of the sex-specific additive and dominance effects regardless of increasing sample size and heritability.

We conducted an additional simulation study, in which the data simulated under the assumption of no sex-specific differences in all genetic parameters were analyzed by the new and traditional model. As expected, both the models provide reasonable estimates of population and quantitative genetic parameters, with estimation precision increasing with increasing sample size and heritability (data not shown). This, in conjunction with the results in Table 3, suggests that the new model provides a general tool for study the genetic architecture of a complex trait, regardless of whether the genetic control of the trait is sex-specific.

Discussion

The genetic architecture of a quantitative trait is complex in terms of interactions between its underlying genetic factors and various environments including sex [1-3]. However, in many current studies, gene by sex interactions are often ignored simply because existing analytical models are not incorporated by environmental factors. While different phenotypic expressions of a trait between the two sexes can be easily measured [4,5], sex-specific discrepancy in the genetic control of the trait can be discerned only when a sophisticated model is used. There is strong evidence for sex-specific genetic influence [17,18] even for traits that display no sexual dimorphism [19]. Although quantitative genetic models are available to estimate sex-specific heritabilities due to aggregative effects of many genes [5,19-22] or map sex-specific QTLs for phenotypic variation [4,5,11-13,23], the model proposed in this article can dissect sex-specific genetic control at the DNA sequence levels.

Our model is founded on the conceptual framework for haplotyping a trait with single nucleotide polymorphisms (SNPs) formulated by Liu et al. [14]. Since haplotypes constructed by physically associated SNPs are thought to affect the expressivity of a complex trait [23,24], it is more precise to characterize such haplotype effects by incorporating gene by gene and gene by sex interactions. Lin and Wu [25] extended Liu et al.'s [14] model to estimate haplotype-haplotype interactions. The new model reported here can not only estimate sex-specific genetic parameters, but also provide a series of statistical procedures for testing sex-specific differences in the genetic architecture of quantitative variation. In a natural population, the structure and pattern of genetic variation can be studied by population genetic parameters, such as haplotype frequencies, allele frequencies and linkage disequilibria. Thus, the understanding of differences in these parameters between the two sexes help to infer the sex-specific genetic structure of a natural population and its evolutionary processes. As shown through simulation studies, our model is alert to discern sex-specific differences in basic population genetic parameters.

Liu et al. [14] assumed a so-called reference or risk haplotype that triggers an effect on complex traits in a different way from the other haplotypes. Thus, the combinations between the risk and non-risk haplotypes (composite diplotypes) will perform differently, depending on the type of combination, i.e., risk by risk, risk by non-risk and non-risk by non-risk. Liu et al. [14] proposed the concepts of the additive effect due to the substitution of the non-risk (or risk) haplotype by the risk (or non-risk) haplotype and the dominance effect due to the interaction between the risk and non-risk haplotypes. These concepts have been integrated into the current model that allows the test of additive by sex and dominance by sex interaction effects. Simulation studies suggest that the new model displays adequate power to detect differences in these quantitative genetic parameters between the sexes, although the detection of sex-specific dominance effects needs a much larger sample size and/or heritability level.

Our model was used to analyze a real data set for pain genetics. Our analyses of a pain sensitivity trait-baseline pressure pain threshold measured at the ulna-by estimating genetic parameters and testing their sex-specific differences in a combination of male and female samples revealed that males and females have different population structure at two SNPs genotyped from a candidate gene (delta opioid receptor) for human pain and that haplotypes exert different genetic effects on the trait between the two sexes. A further test indicates that the risk haplotype [TT] detected by the model exemplifies sex-specific modes of inheritance in affecting the pain trait. While there are no differences among composite diplotypes in males, a significant additive effect was detected in females. Both additive and dominance effects due to the risk haplotype identified are different between the two sexes. Anholt and Mackay [2] described three major mechanisms that explain sex-specific difference in trait control, i.e., sex-specific effects (a gene affects only one sex), sex-biased effects (a gene affects both sexes but to different degrees), and sex-antagonistic effects (a gene affects both sexes but in opposite directions). In our example, the functional haplotype detected affects the pain trait in a sex-antagonistic effect manner, a mechanism thought to help the maintenance of genetic variation in natural populations [26].

In practice, failure to model sex-specific architecture may significantly hamper the ability to detect signals of functional genetic variants in genomewide screens. Although combining male and female data to increase sample size are tempting approaches to increase power, the estimates in this way will be biased from true sex-specific differences. Possible mechanisms that cause sex differences include parent-of-origin effects [27], linkage to or interaction with sex chromosomes, or differences arising from sex-specific hormonal environments. Our gene by sex interaction model that is incorporated by these mechanisms can be modified to consider interactions between genes and any other environments such as life style. Our interaction models should provide a more powerful tool to draw a detailed and precise picture of the genetic architecture of any complex traits that are important to human health.

Authors' contributions

CW, YC, TL and QL wrote the programs and performed data analyses. RF, MW, and RL designed the experiment. LK performed the experiment. RW conceived the idea and wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

The preparation of this manuscript is partially supported by a joint NSF/NIH grant (0540745) to RW and an NIH grant (R01 NS041670-05A1) to RBF.

References

  1. Mackay TF: The genetic architecture of quantitative traits: lessons from Drosophila.

    Curr Opin Genet Dev 2004, 14(3):253-257. PubMed Abstract | Publisher Full Text OpenURL

  2. Anholt RR, Mackay TFC: Quantitative genetic analyses of complex behaviours in Drosophila.

    Nature Reviews Genetics 2004, 5:838-849. PubMed Abstract | Publisher Full Text OpenURL

  3. Foerster K, Coulson T, Sheldon BC, Pemberton JM, Clutton-Brock TH, Kruuk LEB: Sexually antagonistic genetic variation for fitness in red deer.

    Nature 2007, 447:1107-1110. PubMed Abstract | Publisher Full Text OpenURL

  4. Weiss LA, Pan L, Abney M, Ober C: The sex-specific genetic architecture of quantitative traits in humans.

    Nature Genetics 2006, 38:218-222. PubMed Abstract | Publisher Full Text OpenURL

  5. Weiss LA, Abney M, Cook EH, Ober C: Sex-specific genetic architecture of whole blood serotonin levels.

    American Journal of Human Genetics 2005, 76:33-41. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Fillingim RB: Sex-related influences on pain: a review of mechanisms and clinical implications.

    Rehabilitation Psychology 2003, 48:165-174. Publisher Full Text OpenURL

  7. Fillingim RB: Sex, gender and pain: The biopsychosocial model in action XX vs. XY.

    The International Journal of Sex Differences in the Study of Health, Disease and Aging 2003, 1:98-101. OpenURL

  8. Craft RM, Mogil JS, Aloisi AM: Sex differences in pain and analgesia: the role of gonadal hormones.

    European Journal of Pain 2004, 8:397-411. PubMed Abstract | Publisher Full Text OpenURL

  9. Rinn JL, Snyder M: Sexual dimorphisms in mammalian gene expression.

    Trends in Genetics 2005, 21:298-305. PubMed Abstract | Publisher Full Text OpenURL

  10. Mogil JS, Richards SP, O'Toole LA, Helms ML, Mitchell SR, Kest B, Belknap JK: Identification of a sex-specific quantitative trait locus mediating nonopioid stress-induced analgesia in female mice.

    Journal of Neuroscience 1997, 17:7995-8002. PubMed Abstract | Publisher Full Text OpenURL

  11. Zhao W, Ma CX, Cheverud JM, Wu RL: A unifying statistical model for QTL mapping of genotype-sex interaction for developmental trajectories.

    Physiological Genomics 2004, 19:218-227. PubMed Abstract | Publisher Full Text OpenURL

  12. Butterfield RJ, Roper RR, Rhein DM, Melvold RW, Haynes L, Ma RZ, Doerge RW, Teuscher C: Sex-specific QTL govern susceptibility to Theilers murine en-cephalomyelitis virus-induced demyelination (TMEVD).

    Genetics 2003, 163:1041-1046. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Femandez JR, Vogler GP, Tarantino LM, Vignetti S, Plomin R, McClearn GE: Sex-exclusive quantitative trait loci influences in alcohol-related phenotypes.

    American Journal of Medical Genetics 1999, 88:647-652. PubMed Abstract | Publisher Full Text OpenURL

  14. Liu T, Johnson JA, Casella G, Wu RL: Sequencing complex diseases with HapMap.

    Genetics 2004, 168:503-511. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. The International HapMap Consortium: The International HapMap Project.

    Nature 2003, 426:789-794. PubMed Abstract | Publisher Full Text OpenURL

  16. Fillingim RB, Kaplan L, Staud R, Ness TJ, Glover TL, Campbell CM, Mogil JS, Wallace MR: The A118G single nucleotide polymorphism of the mu-opioid receptor gene (OPRM1) is associated with pressure pain sensitivity in humans.

    Journal of Pain 2005, 6:159-167. PubMed Abstract | Publisher Full Text OpenURL

  17. Jensen H, Sather B-E, Ringsby TH, Tufto J, Griffth SC, Ellegren H: Sexual variation in heritability and genetic correlations of morphological traits in house sparrow (Passer domesticus).

    Journal of Evolutionary Biology 2003, 16:1296-1307. PubMed Abstract | Publisher Full Text OpenURL

  18. Leips J, Mackay TFC: Quantitative trait loci for life span in Drosophila melanogaster: interactions with genetic background and larval density.

    Genetics 2002, 155:1773-1788. OpenURL

  19. Schousboe K, Visscher PM, Henriksen JE, Hopper JL, Sorensen TIA, Kyvik KO: Twin study of genetic and environmental influences on glucose tolerance and indices of insulin sensitivity and secretion.

    Diabetologia 2003, 46:1276-1283. PubMed Abstract | Publisher Full Text OpenURL

  20. van Beijsterveldt CE, Bartels M, Hudziak JJ, Boomsma DI: Causes of stability of aggression from early childhood to adolescence: a longitudinal genetic analysis in Dutch twins.

    Behavioral Genetics 2003, 33:591-605. Publisher Full Text OpenURL

  21. Keski-Rahkonen A, Neale BM, Bulik CM, Pietil KH, Pietilainen KH, Rose RJ, Kaprio J, Rissanen A: -bf Intentional weight loss in young adults: Sex-specific genetic and environmental effects.

    Obesity Research 2005, 13:745-753. PubMed Abstract | Publisher Full Text OpenURL

  22. Ober C, Pan L, Phillips N, Parry R, Kurina LM: Sex-specific genetic architecture of asthma-associated quantitative trait loci in a founder population.

    Curr Allergy Asthma Rep 2006, 6(3):241-246. PubMed Abstract | Publisher Full Text OpenURL

  23. Bader JS: The relative power of SNPs and haplotype as genetic markers for association tests.

    Pharmacogenomics 2001, 2(1):11-24. PubMed Abstract | Publisher Full Text OpenURL

  24. Judson R, Stephens JC, Windemuth A: The predictive power of haplotypes in clinical response.

    Pharmacogenomics 2000, 1:5-26. PubMed Abstract | Publisher Full Text OpenURL

  25. Lin M, Wu RL: Detecting sequence-sequence interactions for complex diseases.

    Current Genomics 2006, 7:59-72. Publisher Full Text OpenURL

  26. Rice WR: Sexually antagonistic genes: experimental evidence.

    Science 1992, 256:1436-1439. PubMed Abstract | Publisher Full Text OpenURL

  27. Reik W, Walter J: Genomic imprinting: parental influence on the genome.

    Nature Reviews Genetics 2001, 2:21-32. PubMed Abstract | Publisher Full Text OpenURL