Skip to main content

Sex dependent imprinting effects on complex traits in mice



Genomic imprinting is an epigenetic source of variation in quantitative traits that results from monoallelic gene expression, where commonly either only the paternally- or the maternally-derived allele is expressed. Imprinting has been shown to affect a diversity of complex traits in a variety of species. For several such quantitative traits sex-dependent genetic effects have been discovered, but whether imprinting effects also show such sex-dependence has yet to be explored. Moreover, theoretical work on the evolution of sex-dependent genomic imprinting effects makes specific predictions about the phenotypic patterns of such effects, which, however, have not been assessed empirically to date.


Using a genome-scan for loci affecting a set of complex growth and body composition traits from an intercross between two divergent mouse strains, we investigated possible sex-dependent imprinting effects. Our results demonstrate for the first time the existence of genomic imprinting effects that depend on sex and are not related to sex-chromosome effects. We detected a total of 13 loci on 11 chromosomes that showed significant differences between the sexes in imprinting effects. Most loci showed imprinting effects in only one sex, with eight imprinted effects found in males and six in females. One locus showed sex-dependent imprinting effects in both sexes for different traits. The absence of an imprinting effect in one sex was not necessarily indicative of the overall inactivity of the locus in that sex, as for several loci a significant additive or dominance effect was detected. Moreover, three loci exhibited significant additive effects in both sexes but their imprinting effect was restricted to one sex.


Our results clearly show that imprinting effects can be sex-dependent and also suggest that new candidate imprinted loci can be detected when taking account of sex-specific imprinting effects. However, predictions made about the evolution of sex-dependent imprinting effects and associated phenotypic patterns cannot be unequivocally supported at present and further research into the selection pressures applied to the strains of mice used in our study is required.


While most genes are expressed from both the paternally and maternally derived allele, imprinted loci show parent-of-origin-dependent gene expression such that either the paternal or the maternal allele only is expressed [1, 2]. This monoallelic pattern of gene expression is usually caused by differential methylation of the paternal and maternal allele [3] and results in phenotypic differences between the two heterozygotes (where the only difference is the parent-of-origin of the parental alleles) at a given locus [4]. Since the discovery of genomic imprinting over twenty years ago [5], to date over 80 imprinted genes and their associated phenotypic effects have been identified and described in mice (; October 2008). Many imprinted genes have been shown to directly or indirectly affect resource transfer between mother and offspring [68] and thus individual development and growth. The traits affected by genomic imprinting are widespread and are found predominantly in mammalian systems, suggesting that imprinting has evolved for different reasons at different loci, or may be maintained by selection because it serves different functions at different stages in development and in different tissues [2, 9]. Among traits affected by imprinting are weight-related [8] and body composition traits [10] cognitive abilities [11, 12] and associated disorders (e.g. Prader-Willi and Angelman syndrome), and obesity [7, 1316], highlighting the impact of epigenetic sources to variation in complex traits.

Sex-dependent genetic effects

Previous research has demonstrated that patterns of genetic variation such as additive and dominance effects (the former being a measure of the independent effects of alleles, generally defined as half the difference between the average phenotypes of homozygotes, and the latter being a measure of interactions between alleles at a locus, and defined as the deviation of the average phenotype of the heterozygote from the midpoint between that of the two homozygotes) of loci can show considerable differences depending on the sex of the individual [17, 18]. Such sex-dependent genetic effects have been reported for a number of different traits, including blood pressure [19], longevity in humans [20] and Drosophila [21], bodyweight [22, 23], and morphology in mice [24], and various disease-related traits in humans [25, 26]. However, we are not aware of any study that has explored whether autosomal sex-specific genomic imprinting effects exist, although sex-dependent imprinting effects have been suggested for the X chromosome [27, 28] as a potential explanation for disorders such as Turner syndrome and autism [2931]. Thus, we currently know very little about the likelihood and patterns of differences in imprinting effects in males and females.

Evolutionary and genetic significance of sex-dependent imprinting effects

Sex-dependent imprinting effects are of interest both from a genetic and an evolutionary perspective. For example, analyses of sex-dependent imprinting effects might reveal loci that either show imprinting in only one sex but bi-parental expression of alleles in the other sex. Alternatively, loci could be imprinted in both sexes but show different patterns of imprinting (e.g., maternal expression in one sex but paternal expression in the other). Such loci would suggest flexibility in expression patterns and possibly the underlying imprinting mechanism that may not be expected from current models of imprinting. It is also possible that some loci might show sex-dependent imprinting because they exhibit strictly sex-limited expression, where expression is limited to either males or females. Given that the sexes across taxa often differ considerably but are virtually genetically identical except for the sex chromosomes, sex-dependent gene expression has been predicted to be a relatively common phenomenon and its existence has been demonstrated in a number of different systems [18, 32]. Therefore, we may expect a priori that some imprinted loci would show sex-dependent expression.

From an evolutionary perspective, sex-dependent imprinting effects are of particular interest because theory predicts the existence of such differences in imprinting between the two sexes [33] when selection acts differentially on traits expressed in males versus females (e.g. on loci affecting size dimorphism). Specifically, Day and Bonduriansky [33] predict silencing of the allele derived from the sex that does not experience selection. For example if males are under selection for increased body size but females are not, loci affecting body size are predicted to show paternal expression in males. Moreover, sex-dependent imprinting effects might act in same manner as has been suggested for sex-dependent genetic effects [[34] but see [35]] to help maintain genetic variation because selection may favour alternative alleles or imprinting patterns in males and females.


In this study, we set out to explore whether sex-dependent imprinting effects on complex traits exist in mice by investigating sex differences in the imprinting effects of quantitative trait loci (QTL) in males and females. We characterize the phenotypic patterns caused by genomic imprinting in males and females separately and also determine whether such loci exhibit sex-specific additive or dominance effects. Furthermore, we estimate the proportion of phenotypic variance explained by sex-dependent imprinting effects. Results demonstrate that genomic imprinting effects can be different between the sexes and reveal the existence of 13 QTL with such patterns located on 11 autosomal chromosomes.


Overall, we detected a total of 13 QTL on 11 chromosomes affecting both body and organ weight measures that showed a significant interaction of sex and genomic imprinting effect (SbiQTL, where 'Sbi' refers to sex by imprinting interaction). For all QTL, we confirmed that these sex interactions were true genomic imprinting effects rather than maternal genetic effects (see Methods).

Nine of the ten measured weekly bodyweights were affected by an interaction between sex and imprinting as was postweaning growth (Additional file 1). In addition, the QTL exhibiting sex-dependent imprinting effects also significantly affected reproductive fatpad, spleen, kidney, heart and liver weights (Additional file 1). Two loci (Sbi15.1 and Sbi17.1) were found to affect all weight traits except for week 1 and week 3 weight, whereas six other SbiQTL affected between one and seven weight traits each. Four SbiQTL influencing organ weights, by contrast, showed an effect on only one specific trait. Focusing on the traits themselves, we found that week 5 body weight was modulated by the highest number of SbiQTL (seven) followed by week 4 weight (six). Organ weights were generally affected by fewer SbiQTL. While liver weight was affected by two SbiQTL, all other organ traits were found to be influenced by three loci.

Sex-dependent genomic imprinting effects

For most loci, the interaction effect was due to the imprinting effect occurring in one sex only. At three loci, the imprinting effect was restricted to females for all affected traits and in five other loci this effect was found in males only. For example, Sbi3.1 showed a bipolar imprinting pattern in males but has no imprinting effect in females (Figure 1) whereas at Sbi6.1 no imprinting effect occurred in males but only in females with a paternal expression pattern (Figure 2). At several loci, no significant imprinting effect in either sex was found for some or all traits despite the fact that the locus shows a significant interaction effect (i.e., a significant difference between the sexes in the imprinting effect; see Additional file 2). In these cases, the estimated imprinting effects are opposite in sign to each other in the two sexes but not strong enough in either sex to be significant. Sbi15.1 showed a somewhat different pattern where the imprinting effect occurred in both sexes for separate traits. While the sex by imprinting interaction effect was very strong for many traits for this locus (exceeding genome-wide significance levels), the imprinting effects in the separate sexes were often only marginally significant.

Figure 1
figure 1

Genomic imprinting effects at Sbi3.1 for week 3 body weight in males (black bars) and females (grey bars) for the four ordered genotypes. There is no difference between female heterozygotes but a significant difference in males showing a bipolar imprinting pattern. Shown are the mean trait values for the four ordered genotypes with their standard errors (SE).

Figure 2
figure 2

Genomic imprinting effects at Sbi6.1 for week 6 body weight in males (black bars) and females (grey bars) for the four ordered genotypes. Here, there is no difference between male heterozygotes but a significant difference in females showing a paternal imprinting pattern. Shown are the mean trait values for the four ordered genotypes with their standard errors (SE).

Interestingly, the fact that an imprinting effect occured in one sex only is not necessarily due to the overall inactivity of the locus in the other sex (i.e. no detectable additive or dominance effect). Rather, at four loci that do not show an imprinting effect in males, we found a significant additive or dominance effect in males (Additional file 2). The reverse was found in females only at two loci. Another indication that overall sex-specific expression of the locus does not sufficiently explain the occurrence of sex-dependent imprinting effects is given by the fact that three loci displayed a significant main additive effect in both sexes but the imprinting effect was restricted to either males or females.

Imprinting patterns and variance explained

The imprinting patterns found at a locus can generally be classified into three canonical categories (though other more complex patterns that deviate from these canonical patterns can occur), parental expression (paternal or maternal expression), bipolar dominance and polar over- or underdominance [10, 36], depending on the pattern of variation among the four ordered genotypes at a locus (i.e., alleles are ordered by the parent-of-origin of alleles, thus allowing the distinction between the alternative heterozygotes). Parental expression occurs when the genotypes sharing the same allele inherited from the same parent have the same phenotype. For example, with paternal expression, individuals that inherited the same paternal allele show the same average phenotype. Bipolar dominance is defined as the pattern when the two heterozygotes are significantly different from each other but the two homozygotes are not, such that the genotypic value of one heterozygote is larger than that of the homozygotes while the other is smaller. Finally, polar over- or underdominance is defined by the pattern where one of the two heterozygotes is significantly different from the three other genotypes (either larger, in the case of polar overdominance or smaller, in the case of polar underdominance), which are not significantly different from each other. Where imprinting appeared in only one sex, we found paternal and maternal expression, polar under- and overdominance and bipolar expression (Additional file 1). Among the loci that showed imprinting in males only, four showed paternal expression, three bipolar dominance and one polar overdominance where the SL heterozygote was larger than all other three genotypes. Imprinting patterns in females were paternal expression for three loci, and polar over- and underdominance as well as bipolar dominance once each in females.

The proportion of phenotypic variance explained by the imprinting and sex interaction ranged from 0.58% to 3.6% depending on the trait (Additional file 1). For example, the interaction at Sbi6.1 accounted for 3.28% of phenotypic variance of week 4 body weight but this value decreased to just 0.75% for week 10 body weight.


The most important results emerging from our study are firstly that genomic imprinting effects on complex traits can depend on sex and, secondly, that new loci can be found in sex-specific analyses that may escape detection in searches for imprinting effects in the population as a whole. Our analysis of sex-dependent imprinting effects on complex traits uncovered a total of 13 SbiQTL where often the imprinting effect was restricted to one sex only. Comparing the locations of these SbiQTL with previously-detected imprinted QTL, i.e. those that exert a significant main effect on a given trait [10, 36], many QTL found here are novel loci with the exception of Sbi7.1 and those located on the proximal region on chromosome 3, which are within the confidence regions of previously described loci. This is a particularly important result since it suggests the existence of more imprinted loci that would remain undetected if males and females were not differentiated.

Imprinting patterns

Our study also showed interesting imprinting patterns in the separate sexes. At five loci, imprinting effects were detected in males only while at three loci such effects were found in females only. In the vast majority of cases, the sign of the imprinting effect was opposite in the two sexes at loci in which the imprinting effect occurred in one sex only (e.g. Sbi1.1: negative in males and positive in females, see Additional file 2). Overall, the imprinting effect patterns for a given sex were diverse, ranging from parental (paternal or maternal) expression to bipolar dominance and polar over- or underdominance. We note that when an imprinting effect was limited to one sex only, other genetic effects (additive and dominance) may still have an effect in the other sex or in both sexes. This suggests that the overall expression of the locus is not sex-specific but rather that the imprinting effect alone is sex-dependent.

Comparison with genetic effects

Compared to the occurrence of sex-dependent additive and dominance effects reported in the F2 SM/J × LG/J cross [23], sex-dependent imprinting effects were more common for weight and growth traits but each individual locus was found to affect fewer traits. Vaughn et al. [23] report five QTL on chromosomes 9, 12, 15 and 16 that showed significant sex-specific additive and dominance interactions for 11 different weight and growth traits with a total of 27 traits affected by all loci (one trait can be affected by multiple loci). This compares with eight SbiQTL affecting a total of nine weight and growth traits where seven loci affected two or more traits each. Of the five loci with sex-dependent additive or dominance effects found by Vaughn et al. [23], four had an additive and dominance effect in males and were not expressed in females at all (i.e. they showed neither a significant additive nor a dominance effect) whereas on chromosome 16 there was a male-specific QTL at the proximal part of the chromosome and a female-specific QTL on the distal portion of the chromosome. Vaughn et al.'s [23] findings differ from ours in that a locus that does not show an imprinting effect in one sex can still have additive effects in the same sex or even both sexes. This latter result may be somewhat surprising, because it suggests that such loci do not show simple sex-specific expression in the sense that overall expression is restricted to one sex. Rather, it appears that the imprinting mechanism is not acting in one sex. Perhaps differential methylation at the locus in males and females accounts for this observation, or there may be some trans control of expression that depends on sex.

Genome locations of SbiQTL

With the exception of Sbi6.1 the locations of our detected SbiQTL do not map to regions where currently known imprinted genes are located, although we note that the confidence intervals of Sbi15.1 and Sbi7.1 are just outside of a large imprinting cluster on chromosome 7 and the region containing Peg13 on chromosome 15 (Additional file 2). The confidence region for Sbi6.1 on chromosome 6 is close to a cluster with several imprinted genes, amongst them paternally expressed gene 10 (Peg10) and Sgce [37]. In the light of our finding that imprinting effects can be sex-dependent it might be worthwhile to investigate whether known imprinted genes such as the aforementioned show similar effects, perhaps by comparing male and female phenotypes of knock-out mutants separately or by looking for differences in expression levels.

It should be noted that the primary objective of this study was not to find imprinted loci per se (see [10, 36] for details on these analyses) but rather to search specifically for sex by imprinting interaction effects. In a simulation study, Luedi et al. [38] predicted a total of 600 imprinted genes across the genome. While some congruence of our confidence intervals and the locations of predicted genes might be expected by chance, the loci identified in this study may yield novel candidate gene locations that could be identified in fine mapping or gene expression studies.

The evolution of sex-dependent imprinting effects

One possible mechanism underlying our findings involves the role of modifier loci that regulate the expression of imprinted loci affecting the measured traits. In this scenario, the expression of the modifier locus is sex-specific (e.g. comparable to sex-specific effects of additive loci as reported e.g. in Vaughn et al. [23] and occurs during gametogenesis in offspring. Day & Bonduriansky [33] invoke such 'sexually dimorphic imprinting' as part of an intralocus sexual conflict model to explain the evolution of genomic imprinting. Assuming that some loci have positive fitness effects when expressed in one sex but negative when expressed in the other, imprinting is predicted to evolve because selection favours silencing of loci in the sex in which these loci are not under directional selection. Among the traits with potential sexually antagonistic effects proposed by Day & Bonduriansky [33] is growth because size dimorphism observed in many species suggests that loci modulating growth may be under sex-specific selection. Applied to loci that affect growth and assuming antagonistic selection working on males and females, one would predict from this hypothesis that such loci should show paternal expression in males and maternal expression in females. While mice in general show sexual size dimorphism and are, under natural conditions, likely to be under selection for larger males to some degree, the specific selective regime applied to the strains used in our study does not suggest a sex difference. Our results do not show that maternal expression is more common in females than in males and paternal expression is found equally often in males and females and the pattern seen at loci affecting growth (Sbi3.2 and Sbi4.1) is bipolar dominance in males. Thus, at present it seems rather more difficult to assess the predictions made using our data, in particular in the absence of clear evidence for differential selection in males and females in the founding lines.


Results of our study showing sex-dependent imprinting have not been previously reported to our knowledge. Nonetheless, existing data could be analysed to investigate whether such effects occur in other populations. For example, at the phenotypic level one could compare the growth of male and female individuals carrying a mutation at an imprinted locus (e.g. [8]) to the growth of the wildtype and explore whether the growth difference is proportionally greater in one sex than the other. It may also be possible to directly compare differences in expression levels of imprinted genes in both sexes in a given tissue.

We furthermore demonstrate that candidate imprinted loci can be found by analysing sex-specific imprinting effects as opposed to pooled data because the sign of effect can be opposite in the two sexes and phenotypic patterns can be different. The latter also has important implications for research on disease-related phenotypes associated with dysfunctional parent-of-origin-dependent effects and future research in this area should consider sex-dependent imprinting effects.

Further theoretical studies are required to investigate whether sex-specific imprinting effects can be predicted as a result of asymmetrical fitness effects on patrilineal versus matrilineal relatives as outlined by the currently favoured hypothesis for the evolution of imprinting, the conflict or kinship theory of genomic imprinting [39, 40]. In addition, our results should encourage future theoretical work to address the unresolved issue of whether sex-dependent imprinting effects might play a role in the maintenance of genetic variation.


Study population and phenotypes

In this study we analysed variation in the F3 generation of an intercross between the Large (LG/J) and Small (SM/J) inbred mouse strains [41] that were selected for either large or small body weight at 60 days of age [42]. These strains differ by 6–8 standard deviations in size and growth related traits [41] and therefore represent an excellent model system to study imprinting effects arising from genes regulating growth and development. The study population was generated by first mating ten males of the SM/J strain to ten females of the LG/J strain resulting in the F1 population of 52 individuals. These F1 individuals were randomly mated to produce 510 F2 animals, representing our parental generation. Random mating among F2 animals yielded 200 full-sib families of the F3 generation with a total of 1632 individuals. Further details of the husbandry are given in Vaughn et al. [23].

The traits analysed in this study are weekly body weights taken in mice 1 to 10 weeks of age. These measures were used to calculate two growth periods: pre-weaning (weeks 1–3) and post-weaning growth (weeks 3–10). Growth was calculated as the difference between weekly weights (i.e., gain in body weight) such that, for example, the growth from week 1 to week 3 is the difference between week 3 weight and week 1 weight. Pups were weighed weekly using a digital scale with an accuracy of 0.1 g. In addition, we analysed heart, kidney, liver, reproductive fatpad and spleen weights. To obtain these organ weights, mice were sacrificed after 70 days of age or after having produced and reared their offspring to weaning (at three weeks of age). At necropsy, all mice were first weighed to obtain an overall measure of body size and then dissected and their organs weighed to the nearest 0.01 gram using digital scales. The effects of age at necropsy, sex, and litter size at birth were removed from the data and the residuals used in the analysis prior to gene mapping [41].


DNA was extracted from livers of both the F2 and F3 individuals using Qiagen DNeasy tissue kits and samples were scored for 384 SNP markers using the Golden Gate Assay by Illumina, San Diego, USA. These markers were previously found to be polymorphic between the two strains Fifteen loci had to be excluded from the analysis because they were not reliably scored. In addition, 16 loci were scored on the X chromosome, but were not included in this analysis because the statistical model for the X chromosome is currently unresolved. Thus, in total we analysed 353 loci across the 19 autosomes in this study. A genetic map of these markers in cM was produced using R/QTL and validated against the genome coordinate locations in the Ensembl database The average map distance between markers in the F2 generation was 4 cM. Markers were reasonably evenly located throughout the genome except for several regions in which LG/J and SM/J strains have been found to be monomorphic [43].

Haplotype recontruction

Both parental and offspring genotypes were used to reconstruct haplotypes for all mice using PedPhase [44], which produced a set of unordered haplotypes for the F2 parents and a set of ordered haplotypes (i.e. ordered by parent-of-origin of alleles) for their F3 offspring. Thus, we were able to distinguish the four possible genotypes at a given locus, LL, SL, LS or SS where the first allele refers to the paternally-derived allele and the second to the maternally-derived copy.

Analysis of parent-of-origin-dependent effects

The four ordered genotypes at the marker loci (LL, LS, SL and SS) were assigned additive (a), dominance (d) and parent-of-origin (i) genotypic index scores following Mantey et al. [45]. In matrix form, these index scores are given by:

[ L L ¯ L S ¯ S L ¯ S S ¯ ] = [ 1 1 0 0 1 0 1 1 1 0 1 1 1 1 0 0 ] [ r a d i ] yielding estimates of the parameters : [ r a d i ] = [ L L ¯ 2 + S S ¯ 2 L L ¯ 2 S S ¯ 2 L S ¯ 2 + S L ¯ 2 L L ¯ 2 S S ¯ 2 L S ¯ 2 S L ¯ 2 ] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaamWaaeaafaqabeabbaaaaeaadaqdaaqaaiabdYeamjabdYeambaaaeaadaqdaaqaaiabdYeamjabdofatbaaaeaadaqdaaqaaiabdofatjabdYeambaaaeaadaqdaaqaaiabdofatjabdofatbaaaaaacaGLBbGaayzxaaGaeyypa0ZaamWaaeaafaqabeabeaaaaaqaaiabigdaXaqaaiabigdaXaqaaiabicdaWaqaaiabicdaWaqaaiabigdaXaqaaiabicdaWaqaaiabigdaXaqaaiabigdaXaqaaiabigdaXaqaaiabicdaWaqaaiabigdaXaqaaiabgkHiTiabigdaXaqaaiabigdaXaqaaiabgkHiTiabigdaXaqaaiabicdaWaqaaiabicdaWaaaaiaawUfacaGLDbaadaWadaqaauaabeqaeeaaaaqaaiabdkhaYbqaaiabdggaHbqaaiabdsgaKbqaaiabdMgaPbaaaiaawUfacaGLDbaacqqG5bqEcqqGPbqAcqqGLbqzcqqGSbaBcqqGKbazcqqGPbqAcqqGUbGBcqqGNbWzcqqGGaaicqqGLbqzcqqGZbWCcqqG0baDcqqGPbqAcqqGTbqBcqqGHbqycqqG0baDcqqGLbqzcqqGZbWCcqqGGaaicqqGVbWBcqqGMbGzcqqGGaaicqqG0baDcqqGObaAcqqGLbqzcqqGGaaicqqGWbaCcqqGHbqycqqGYbGCcqqGHbqycqqGTbqBcqqGLbqzcqqG0baDcqqGLbqzcqqGYbGCcqqGZbWCcqGG6aGodaWadaqaauaabeqaeeaaaaqaaiabdkhaYbqaaiabdggaHbqaaiabdsgaKbqaaiabdMgaPbaaaiaawUfacaGLDbaacqGH9aqpdaWadaqaauaabeqaeeaaaaqaamaaleaaleaadaqdaaqaaiabdYeamjabdYeambaaaeaacqaIYaGmaaGccqGHRaWkdaWcbaWcbaWaa0aaaeaacqWGtbWucqWGtbWuaaaabaGaeGOmaidaaaGcbaWaaSqaaSqaamaanaaabaGaemitaWKaemitaWeaaaqaaiabikdaYaaakiabgkHiTmaaleaaleaadaqdaaqaaiabdofatjabdofatbaaaeaacqaIYaGmaaaakeaadaWcbaWcbaWaa0aaaeaacqWGmbatcqWGtbWuaaaabaGaeGOmaidaaOGaey4kaSYaaSqaaSqaamaanaaabaGaem4uamLaemitaWeaaaqaaiabikdaYaaakiabgkHiTmaaleaaleaadaqdaaqaaiabdYeamjabdYeambaaaeaacqaIYaGmaaGccqGHsisldaWcbaWcbaWaa0aaaeaacqWGtbWucqWGtbWuaaaabaGaeGOmaidaaaGcbaWaaSqaaSqaamaanaaabaGaemitaWKaem4uamfaaaqaaiabikdaYaaakiabgkHiTmaaleaaleaadaqdaaqaaiabdofatjabdYeambaaaeaacqaIYaGmaaaaaaGccaGLBbGaayzxaaaaaa@B56D@

The vectors of genotypic means are L L ¯ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaa0aaaeaacqWGmbatcqWGmbataaaaaa@2E28@ , L S ¯ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaa0aaaeaacqWGmbatcqWGtbWuaaaaaa@2E36@ , S L ¯ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaa0aaaeaacqWGtbWucqWGmbataaaaaa@2E36@ , S S ¯ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaa0aaaeaacqWGtbWucqWGtbWuaaaaaa@2E44@ ,, r is the reference point for the model (the mid-point between homozygotes), a is the additive genotypic value (half the difference between homozygotes), d is the dominance genotypic value (the difference between the mean of the heterozygotes and the mid-point of the homozygote means), and i is the parent-of-origin or imprinting genotypic value (half the difference between heterozygotes) (cf. [45]).

These index scores were used to build a model for a genome scan for loci showing significant sex by parent-of-origin-dependent interaction effects (i.e. sex by i effects) using the Proc mixed Procedure in SAS 9.1; SAS Institute, Cary, NC, USA using maximum likelihood (see [46] for details). In this model, growth traits or organ weights are the dependent variables and sex, the additive, dominance and imprinting index scores as well as their interactions with sex were the fixed effects and family was a random effect. The mixed model with the fixed genetic effects and random family effect was used to scan the genome to produce a probability distribution for the overall effect of the locus as well as the sex by imprinting interaction effect. The probability was generated by comparing the -2 res log likelihood computed by SAS for the model with the fixed genetic effects and their interaction with sex with a reduced model that did not include these six genetic effects. The difference in the -2 res log likelihoods of the two models is chi-square distributed with 6 degrees of freedom (representing the fact that the two models differ by 6 fixed effects). These probability values were then transformed to a log probability ratio (LPR) in order to make them comparable to the LOD scores commonly seen in QTL analyses (LPR = -log10 [probability]).

Parent-of-origin effects can also appear as a result of maternal genetic effects rather than genomic imprinting [47]. Maternal genetic effects occur if the mother's genotype affects the offspring's phenotype beyond the effects of her genetic contribution to offspring phenotypes. Maternal effects can result in the appearance of parent-of-origin dependent effects because homozygous mothers can produce only one type of heterozygous offspring. By analogy, sex-dependent maternal effects can result in the appearance of sex by i effects. Therefore, we tested all loci with a significant sex by i interaction effect to confirm that the effect was not caused by a sex-dependent maternal genetic effect [cf. [47]]. Using our mixed model framework, we have advanced on our previous approach by testing whether the sex by i effect is dependent on whether individuals are reared by homozygous mothers (where a maternal effect could result in the appearance of a parent-of-origin dependent effect) or heterozygous mothers (where all four ordered genotypes of offspring come from the same genotype of mother and hence maternal effects cannot result in the appearance of a parent-of-origin dependent effect). A parent-of-origin dependent effect was attributed to a maternal effect if the sex by i effect depended significantly on the type of mother (i.e., there was a significant three-way interaction between sex, i and the type of mother, classed as heterozygote or homozygote).

Significance testing of QTL

Our significance testing approach needed to take account of the autocorrelation of siblings within the F3 families. Therefore, we first calculated chromosome-wise and genome-wise threshold LOD scores for the sex-by-imprinting interaction in separate permutation procedures [48] for each of the traits ensuring that the specific family structure in this generation was maintained. For each trait, we achieved this by first calculating deviations of each individual from its family mean and randomly permuting these deviations within each family. We then randomly permuted all F3 family means and reconstructed new values for each individual by adding its permuted deviation to its new mean. Using these new values for each individual, we ran a canonical correlation analysis and computed the highest LPR score on each chromosome. This procedure was repeated 1000 times, and 5% chromosome-wise threshold values were obtained from the 50th highest values generated for each chromosome. The 5% genome-wise threshold value for each trait was obtained from the 50th highest value among the 1000 highest LPR values across all 19 chromosomes in each permutation run [48]. First simulations revealed that the genome- and chromosome-wise thresholds were very similar to those obtained using the effective number of markers method based on the eigenvalues of the marker correlation matrix [49]. Since the latter method allows for direct computation of the thresholds for all traits, whereas the simulation required significant computing time for each trait, we used the thresholds obtained using the effective number of markers. This method calculates the number of independent tests in a genome or chromosome scan and uses the "effective number of markers" in a Bonferonni correction. We used both the conservative genome-wise significance threshold as well as the chromosome-wise thresholds because this approach has been shown to give overall the best results by increasing the discovery of true positives while at the same time avoiding problems using the false discovery rate in gene mapping experiments [50].

After loci were identified using the genome- or chromosome-wide significance thresholds for the interaction model we used post-hoc tests to characterize the phenotypic patterns caused by genomic imprinting. We included all significant effects of QTL using a protected test where pleiotropic effects of QTL are included whenever the effect of a locus on other traits is significant at the pointwise (p < 0.05; LPR > 1.3) level. Thus, while we apply the stringent genome- or chromosome-wide threshold for QTL detection to minimize type I errors, we characterize the distribution of pleiotropic effects across the entire set of traits to minimize type II errors. The imprinting patterns were determined using the mixed model approach and are given by the relationship between the additive (a), dominance (d) and parent-of-origin genotype value (i) [36]. We distinguish a total of three different imprinting patterns: parental (paternal or maternal) expression, bipolar dominance and polar dominance, see above [36].


  1. Bartolomei MS, Tilghman SM: Genomic imprinting in mammals. Annu Rev Genet. 1997, 31: 493-525. 10.1146/annurev.genet.31.1.493.

    Article  CAS  PubMed  Google Scholar 

  2. Wood AJ, Oakey RJ: Genomic imprinting in mammals: Emerging themes and established theories. PLoS Genet. 2006, 2: 1677-1685. 10.1371/journal.pgen.0020147.

    Article  CAS  Google Scholar 

  3. Reik W, Walter J: Genomic imprinting: parental influenceon the genome. Nat Rev Genet. 2001, 2: 21-32. 10.1038/35047554.

    Article  CAS  PubMed  Google Scholar 

  4. Santure AW, Spencer HG: Influence of mom and dad: quantitative genetic models for maternal effects and genomic imprinting. Genetics. 2006, 173: 2297-2361. 10.1534/genetics.105.049494.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Surani MA, Barton SC, Norris ML: Development of reconstituted mouse egges suggests imprinting of the genome during gametogenesis. Nature. 1984, 308: 548-550. 10.1038/308548a0.

    Article  CAS  PubMed  Google Scholar 

  6. Burt A, Trivers R: Genes in Conflict. 2006, Cambridge MA: Harvard University Press

    Book  Google Scholar 

  7. Morison IM, Ramsay JP, Spencer HG: A census of mammalian imprinting. Trends Genet. 2005, 21: 457-465. 10.1016/j.tig.2005.06.008.

    Article  CAS  PubMed  Google Scholar 

  8. Li LL, Keverne EB, Aparicio SA, Ishino F, Barton SC, Surani MA: Regulation of maternal behaviour and offspring growth by paternally expressed Peg3. Science. 1999, 284: 330-333. 10.1126/science.284.5412.330.

    Article  CAS  PubMed  Google Scholar 

  9. Wolf JB, Hager R: A maternal-offspring coadaptation theory for the evolution of genomic imprinting. PLoS Biol. 2006, 4: 2238-2243. 10.1371/journal.pbio.0040380.

    Article  CAS  Google Scholar 

  10. Cheverud JM, Hager R, Roseman C, Fawcett G, Wang B, Wolf JB: Genomic Imprinting Effects on Adult Body Composition in Mice. Proc Natl Acad Sci USA. 2008, 105: 4253-4258. 10.1073/pnas.0706562105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Isles AR, Wilkinson LS: Imprinted genes, cognition and behaviour. Trends Cogn Sci. 2000, 4: 309-318. 10.1016/S1364-6613(00)01504-7.

    Article  PubMed  Google Scholar 

  12. Wilkinson LS, Davies W, Isles AR: Genomic imprinting effects on brain development and function. Nat Rev Neurosci. 2007, 8 (11): 832-843. 10.1038/nrn2235.

    Article  CAS  PubMed  Google Scholar 

  13. Nicholls RD: The impact of genomic imprinting for neurobehavioral and developmental disorders. J Clin Invest. 2000, 105: 413-418. 10.1172/JCI9460.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Constancia M, Kelsey G, Reik W: Resourceful imprinting. Nature. 2004, 432: 53-57. 10.1038/432053a.

    Article  CAS  PubMed  Google Scholar 

  15. Dong C, Li WD, Geller F, Lei L, Li D, et al: Possible genomic imprinting of three human obesity-related genetic loci. Am J Hum Genet. 2005, 76: 427-437. 10.1086/428438.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Heijmans BT, Kremer D, Tobi EW, Boomsma DI, Slagboom PE: Heritable rather than age-related environmental stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus. Hum Mol Genet. 2007, 16: 547-554. 10.1093/hmg/ddm010.

    Article  CAS  PubMed  Google Scholar 

  17. Zhao W, Changxing M, Cheverud JM, Wu R: A unifying statistical model for QTL mapping of genotype x sex interaction for developmental trajectories. Physiol Genomics. 2004, 19: 218-227. 10.1152/physiolgenomics.00129.2004.

    Article  CAS  PubMed  Google Scholar 

  18. Ellengren H, Parsch J: The evolution of sex-biased genes and sex-biased gene expression. Nat Rev Genet. 2007, 8: 689-698. 10.1038/nrg2167.

    Article  Google Scholar 

  19. Franceschini N, MacCluer JW, Göring HHH, Shelley AC, Rose KM, et al: A Quantitative Trait Loci-Specific Gene-by-Sex Interaction on Systolic Blood Pressure Among American Indians. Hypertension. 2006, 46: 266-270. 10.1161/01.HYP.0000231651.91523.7e.

    Article  Google Scholar 

  20. Tan Q, Yashin AI, Bladbjerg EM, de Maat MPM, Andersen-Ranberg K, et al: A Case-Only Approach for Assessing Gene by Sex Interaction in Human Longevity. J Gerontol A Biol Sci Med Sci. 2002, 57 (4): B129-B133.

    Article  PubMed  Google Scholar 

  21. Mackay TCF: Quantitative trait loci in Drospophila. Nat Rev Genet. 2001, 2: 11-20. 10.1038/35047544.

    Article  CAS  PubMed  Google Scholar 

  22. Cheverud JM, Vaughn T, Pletscher LS, Peripato A, Adams E, et al: Genetic architecture of adiposity in the cross of Large (LG/J) and Small (SM/J) inbred mice. Mamm Genome. 2001, 12: 3-12. 10.1007/s003350010218.

    Article  CAS  PubMed  Google Scholar 

  23. Vaughn TT, Pletscher LS, Peripato A, King-Ellison K, Adams E, et al: Mapping quantitative trait loci for murine growth – A closer look at genetic architecture. Genet Res. 1999, 74: 313-322. 10.1017/S0016672399004103.

    Article  CAS  PubMed  Google Scholar 

  24. Parker LT, Leamy LJ: Fluctuating asymmetry of morphometric characters in a randombred population of house mice: The effects of age, sex, and phenotypical extremeness. J Hered. 1991, 82: 145-150.

    CAS  PubMed  Google Scholar 

  25. Butterfield RJ, Roper RJ, Rhein DM, Melvold RW, Haynes L, et al: Sex-specific QTL govern susceptibility to Theiler's murine encephalomyelitis virus-induced demethylation (TMVED). Genetics. 2003, 163: 1041-1046.

    PubMed Central  CAS  PubMed  Google Scholar 

  26. Mogil JS, Richards SP, O'Toole P, Helms ML, Mitchell SR, et al: Identification of a sex-specific quantitative trait locus mediating nonopoid stress-induced analgesia in female mice. J Neurosci. 1997, 17: 7995-8002.

    CAS  PubMed  Google Scholar 

  27. Latham KE: X chromosome imprinting and inactivation in p reimplantation mammalian embryos. Trends Genet. 2005, 21: 120-127. 10.1016/j.tig.2004.12.003.

    Article  CAS  PubMed  Google Scholar 

  28. Davies W, Isles AR, Burgoyne PS, Wilkinson LS: X-linked imprinting: effects on brain and behaviour. Bioessays. 2005, 28: 35-44. 10.1002/bies.20341.

    Article  Google Scholar 

  29. Sagi L, Zuckerman-Levin N, Gawlik A, Ghizzoni L, Buyukgebiz A, et al: Clinical significance of the parental origin of the X chromosome in turner syndrome. J Clin Endocrinol Metab. 2007, 92: 846-52. 10.1210/jc.2006-0158.

    Article  CAS  PubMed  Google Scholar 

  30. Burgoyne PS, Thornhill AR, Boudrean SK, Darling SM, Bishop CE, et al: The genetic basis of XX-XY difference present before gonadal sex differentiation in the mouse. Phil Trans R Soc Lond Ser B. 1995, 350: 253-261. 10.1098/rstb.1995.0159.

    Article  CAS  Google Scholar 

  31. Skuse DH: Imprinting, the X-chromosome, and the male brain: explaining sex differences in the liability to autism. Pediatr Res. 2000, 47: 9-16. 10.1203/00006450-200001000-00006.

    Article  CAS  PubMed  Google Scholar 

  32. Ranz JM, Castillo-Davis CI, Meiklejohn CD, Hartl DL: Sex-dependent gene expression and evolution of the Drosophila transcriptome. Science. 300: 1742-1745. 10.1126/science.1085881.

  33. Day T, Bonduriansky R: Intralocus sexual conflict can drive the evolution of genomic imprinting. Genetics. 167: 1537-1546. 10.1534/genetics.103.026211.

  34. Nuzhdin SV, Pasyukova EG, Dilda CL, Zeng Z, Mackay TFC: Sex-specific quantitative trait loci affecting longevity in Drosophila melanogaster. Proc Natl Acad Sci USA. 1997, 94: 9734-9739. 10.1073/pnas.94.18.9734.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Turelli M, Barton NH: Polygneic variation maintained by balancing selection: pleiotropy, sex-dependent allelic effects and GxE interactions. Genetics. 2004, 166: 1053-1079. 10.1534/genetics.166.2.1053.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Wolf JB, Cheverud JM, Roseman C, Hager R: Genome-wide analysis reveals a complex pattern genomic imprinting. PloS Genet. 2008, 4 (6): e1000091-10.1371/journal.pgen.1000091.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Ono R, Shiura H, Aburatani H, Kohda T, Kaneko-Ishino T, et al: Identification of a large novel imprinted gene cluster on mouse proximal chromosome 6. Genome Res. 2003, 13: 1696-705. 10.1101/gr.906803.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Luedi PP, Hartemink AJ, Jirtle RL: Genome-wide prediction of imprinted murine genes. Genome Res. 2005, 15: 875-884. 10.1101/gr.3303505.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Haig D, Graham C: Genomic imprinting and the strange case of the insulin-like growth factor II receptor. Cell. 1991, 64: 1045-1046. 10.1016/0092-8674(91)90256-X.

    Article  CAS  PubMed  Google Scholar 

  40. Wilkins JF, Haig D: What good is genomic imprinting: The function of parent-specific gene expression. Nat Rev Genet. 2003, 4 (5): 359-368. 10.1038/nrg1062.

    Article  CAS  PubMed  Google Scholar 

  41. Kramer MG, Vaughn TT, Pletscher LS, King-Ellison K, Adams E, et al: Genetic variation for body weight gain and composition in the intercross of Large (LG/J) and Small (SM/J) inbred strains of mice. Genet Mol Biol. 1998, 21: 211-218. 10.1590/S1415-47571998000200007.

    Article  Google Scholar 

  42. Chai CK: Analysis of quantitative inheritance of body size in mice. I. Hybridization and maternal influence. Genetics. 1956, 41 (2): 157-164.

    PubMed Central  CAS  PubMed  Google Scholar 

  43. Hrbek T, de Brito RA, Wang B, Pletscher LS, Cheverud JM: Genetic characterization of a new set of recombinant inbred lines (LGXSM) formed from the intercross of SM/J and LG/J inbred mouse strains. Mamm Genome. 2006, 17: 417-429. 10.1007/s00335-005-0038-7.

    Article  CAS  PubMed  Google Scholar 

  44. Li J, Jiang T: Efficient rule-based haplotyping algorithms for pedigree. J Bioinfo Comp Biol. 2003, 1: 41-69. 10.1142/S0219720003000204.

    Article  CAS  Google Scholar 

  45. Mantey C, Brockmann GA, Kalm E, Reinsch N: Mapping and exclusion mapping of genomic imprinting effects in mouse F2 families. J Hered. 2005, 96: 329-338. 10.1093/jhered/esi044.

    Article  CAS  PubMed  Google Scholar 

  46. Leamy LJ, Routman EJ, Cheverud JM: Quantitative trait loci for early and late developing skull characters in mice: a test of the genetic independence model of morphological integration. Am Nat. 1999, 153: 201-214. 10.1086/303165.

    Article  Google Scholar 

  47. Hager R, Cheverud JM, Wolf JB: Maternal effects as the cause of parent-of-origin-dependent effects that mimic genomic imprinting. Genetics. 2008, 178: 1755-1762. 10.1534/genetics.107.080697.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.

    PubMed Central  CAS  PubMed  Google Scholar 

  49. Li J, Ji L: Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity. 2005, 95: 221-227. 10.1038/sj.hdy.6800717.

    Article  CAS  PubMed  Google Scholar 

  50. Chen L, Storey LD: Relaxed significance criteria for linkage analysis. Genetics. 2006, 173: 2371-2381. 10.1534/genetics.105.052506.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


We thank Will Pitchers for help with haplotype reconstruction and three anonymous referees for helpful comments. This research was supported by grants from the Biotechnology and Biological Sciences Research Council, UK (BBSRC), an Underwood Fellowship from the BBSRC and NIH grant DK055736. RH is supported by a NERC Research Fellowship.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Reinmar Hager.

Additional information

Authors' contributions

JBW and JMC conceived experiment and data analysis. RH analysed data. JBW and LJL conducted significance testing and provided programs used for analyses. RH, JMC, LJM and JBW wrote the paper. All authors contributed to previous drafts of the paper and have approved the final manuscript.

Electronic supplementary material


Additional file 1: Sex-dependent imprinting effect QTL ( Sbi QTL) and their patterns for body weights, growth and organ weights. Coordinates (Mb) are based on mouse genome build 36 and location refers to F2 cM map distances. Patterns are listed as Paternal = paternal expression, Maternal = maternal expression, Bipolar = bipolar imprinting, Under = polar underdominance, Over = polar overdominance, NS = non-significant. Loci were identified when the LPR for the overall effect of a locus exceeded the chromosome-wide or genome-wide threshold and they showed a significant sex by i interaction effect. The genome-wide significance threshold is given at an LPR of 3.41 and the chromosome-wide thresholds are given for each locus (in italics) as is the highest LPR for the locus. Traits shown in bold indicate a QTL effect. (DOC 36 KB)


Additional file 2: The table lists the identified sex by imprinting QTL ( Sbi QTL) and effects for males and females separately. First, the SbiQTL is identified, followed by the trait affected and by marker name, the F2 map location in cM and the genome coordinate based on mouse build 36 along with the confidence intervals. This is followed by the locus LPR and the chromosome-wide threshold (the genome-wide is 3.41), the estimates for the sex by i interaction, the corresponding standard error, the interaction LPR and the r2 for the interaction effect. The next section details for males the additive a, dominance d, and imprinting i effects with their estimates, the estimates divided by the standard deviation and standard error as well as the corresponding LPR. Next, the imprinting patterns are given where 'Paternal' refers to paternal expression, 'Maternal' to maternal expression, 'Bipolar' to bipolar dominance, 'Over' to polar overdominance, 'Under' to polar underdominance and 'NS' to non-significant. This is followed by the genotypic values of the four ordered genotypes with their standard errors SE. We next present in the same order the equivalent information for females. (XLS 84 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Hager, R., Cheverud, J.M., Leamy, L.J. et al. Sex dependent imprinting effects on complex traits in mice. BMC Evol Biol 8, 303 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: