Genetics: Published Articles Ahead of Print, published on October 18, 2007 as 10.1534/genetics.107.076190

Extensive additivity of expression differentiates

subspecies of the house mouse

Ruth Rottscheidt & Bettina Harr

Institute for Genetics, University of Cologne, Zülpicher Strasse 47, 50674 Cologne,

Germany

Address for correspondence: Dr. Bettina Harr, Institute for Genetics, University of

Cologne, Zülpicher Strasse 47, 50674 Cologne, Germany

Email: [email protected]

Tel: ++49 221 470 6617

Fax: ++49 221 470 5975

Keywords: house mouse, , F1 hybrid, speciation

Running title: Additivity of gene expression Abstract

We have studied different subspecies of the house mouse and their reciprocal F1 hybrids to estimate the within-locus mode of inheritance for subspecies differences in gene expression in three tissues (brain, liver, testis) of male mice. This study investigates mode of inheritance in crosses at a larger taxonomic distance than have been previously systematically investigated. We found the vast majority of transcripts to be additively expressed with only a few transcripts showing dominance or overdominance in expression, except for one direction of one cross, which showed large mis-expression in the testis. We suggest that as time passes, more come to influence expression, and if there is no directional dominance, additivity becomes increasingly more likely, up to a threshold beyond which there is F1 hybrid breakdown. Some previous studies on different organisms have found a large degree of dominance, commonly at shorter taxonomic differences. We surveyed these findings and show that the most consistent association exists between the amount of additivity detected in a study and expression analysis method (in particular microarray platform), suggesting that at least some of the differences among studies might be methodological. Most studies agree with ours in that within-locus additivity seems to be general mode of inheritance for transcript expression.

Differentially expressed transcripts identified in our screen among subspecies of house mice are candidate genes that could be involved in reproductive isolation between these subspecies. Introduction

The methods of quantitative genetics are now being applied to transcript abundance, as measured using microarrays. “Expression Quantitative Trait Locus” mapping (eQTL) identifies the genetic basis of expression differences between organisms (ROCKMAN and

KRUGLYAK 2006). Such mapping studies are being used to ask how many loci contribute to a quantitative trait, what is the size of the effect of the loci, and how individual loci interact to generate a quantitative trait (MACKAY 2001). Treating gene expression levels as measured by microarrays as quantitative traits enables the simultaneous assessment of thousands of traits in parallel. eQTL studies can potentially also be used to ask whether regulatory variants, especially those that are important in evolution, are located in cis (i.e. directly linked to the gene that shows the differential expression) or in trans (i.e. somewhere else in the genome) (HOEKSTRA and COYNE 2007; WRAY 2007).

Two approaches have been used to study the genetics of expression traits. In the first, DNA sequence polymorphisms are correlated with expression differences (BREM et al. 2002; DOSS et al. 2005; SCHADT et al. 2003; STOREY et al. 2005). A key finding from these studies is that expression traits are often affected by multiple underlying loci and interactions among them.

The second approach makes use of F1 hybrids generated from a cross between divergent taxa to ask if a transcript’s expression is intermediate (“additive”) to that of the two parents. These studies have produced varying results. GIBSON et al. (2004) crossed two lines of Drosophila melanogaster and found most transcripts to be “nonadditivly” (i.e. the hybrids showed expression levels that were more closely to either on of the parents) expressed in F1 hybrids. Widespread nonadditivity was also identified in the

Pacific Oyster (HEDGECOCK et al. 2007) and to a lesser extend in Arabidopsis

(VUYLSTEKE et al. 2005). However, later studies on Drosophila (HUGHES et al. 2006), maize (STUPAR and SPRINGER 2006; SWANSON-WAGNER et al. 2006), and laboratory strains of the house mouse (CUI et al. 2006) were consistent with mostly additivity of expression. HUGHES et al. (2006) suggested the discrepancy between the different

Drosophila studies might be explained in two ways. First, inbreeding might affect the results. Much within-locus additivity was observed using natural populations of D. melanogaster (HUGHES et al. 2006) while nonadditivity was observed with strongly inbred lines (GIBSON et al. 2004). It is well known that crosses between inbred strains produce offspring that exhibit greater biomass, speed of development, and fertility than both parents (heterosis) (COMINGS and MACMURRAY 2000). Heterosis is attributed to overdominance (superiority of heterozygotes at genes affecting fitness) or dominance

(masking of recessive deleterious mutations) and plausibly could result in overdominant or dominant expression at the transcript level. However the inbreeding explanation seems to be an unlikely general explanation because additional recent studies have found strong evidence for additivity despite the use of inbred parental lines (CUI et al. 2006; STUPAR and SPRINGER 2006; SWANSON-WAGNER et al. 2006). A second explanation for the

Drosophila differences is that genetic architecture of an expression trait might depend critically on the taxonomic level at which the variation is investigated (i.e., between populations as in the GIBSON et al. (2004) study as opposed to within populations as in the HUGHES et al. (2006) study). These authors assume that the more highly diverged the parents are the greater the nonadditivity in gene expression that will be observed.

However, we might expect additivity of expression to become more likely as taxa diverge and the trait becomes influenced by more and more genes, provided that they do not show directional dominance. This would apply up to some certain threshold, beyond which incompatible genes lead to massive mis-expression in F1 hybrids.

Here we investigate patterns of expression in the house mouse subspecies complex of Mus musculus. We used parental mouse strains that were collected in the wild and inbred in the laboratory for several generations. The strains belong to the three different subspecies of house mice, M. m. musculus, M. m. domesticus and M. m. castaneus, with estimated divergence times of between 300,000 to 1 Million years

(BOURSOT et al. 1993). Previous systematic studies focused on comparisons within a single species or used house mouse laboratory strains. In contrast, our study employs natural-derived lines that are rather highly diverged. We found predominantly additive gene expression differences between these subspecies. Materials and Methods

Mouse strains

Wild-derived strains of domesticus (JPC 2705, from Germany), musculus (JPC 2821, from the Czech Republic) and castaneus (CIM, from India) were used for the microarray experiment. The domesticus and musculus strains were provided by J. Pialek. These strains were collected in the wild and inbred by brother-sister matings for ~13 generations in the laboratory of J. Pialek in the Department of Population Biology in

Studenec, Czech Republic. The castaneus strain was provided by A. Orth and F.

Bonhomme. This strain has been kept in the Laboratory Génome, Populations,

Interactions, Adaptation, Montpellier, France for more than 30 generations in a closed colony. We set up reciprocal crosses between domesticus and musculus and between musculus and castaneus to obtain F1 hybrids. For each parental strain and each reciprocal cross two male individuals were analyzed. Altogether, we performed 42 microarray experiments, 14 in each of three tissues (brain, liver and testis). An overview of the samples used in this study is given in Table 1. All mice were raised under identical standard laboratory conditions and were sacrificed at the age of 8-10 weeks.

RNA extraction

We extracted RNA from three different tissues (brain, liver and testis) using Trizol®

(Invitrogen, Carlsbad, CA) following the manufacturer’s protocol. Quality and integrity of the total RNA was controlled by using the Agilent Technologies 2100 Bioanalyzer and the RNA 6000 Nano LabChip® Kit (Agilent Technologies; Waldbronn, Germany). Microarray

Expression profiles were determined for over 39,000 mouse transcripts using the Mouse

Genome 430 2.0 Affymetrix® GeneChip ®. For biotin-labeled target synthesis starting from 3 µg of total RNA we used standard protocols supplied by the manufacturer

(Affymetrix; Santa Clara, CA). After hybridization the GeneChips ® were washed, stained with SA-PE and read using an Affymetrix ® GeneChip ® fluidic station and scanner.

Data analysis

All data processing and statistical analyses were performed using the statistical language

R. Raw signal intensities were normalized and summarized according to the standard

Affymetrix MA Suite 5.0 algorithm using the program Bioconductor

(http://www.bioconductor.org/). Signal intensities were ln-transformed prior to statistical analyses. MA Suite 5.0 expression values have been submitted to the Gene Expression

Omnibus (accession number XX). We define a “group” as one of domesticus, musculus, castaneus, F1 hybrids from one direction of the cross, or F1 hybrids from the other direction of the cross. We called a transcript “expressed” in this group if the average expression level among the replicate samples within this group was > 500. Applying a threshold of >500 yields about 11,0000 transcripts expressed per group per tissue, a value that corresponds to previous studies (SU et al. 2004).

First, we studied all transcripts that are differentially expressed between the parental subspecies and assessed additivity vs. nonadditivity of expression of these transcripts in the F1 hybrids (see Figure 1A, 1B). Second, we investigated those transcripts which do not differ in expression between the parental lines, but which are expressed differently in the F1 hybrids relative to both their parents (Figure 1C):

1) Additivity/nonadditivity of differentially expressed transcripts

For each parental comparison (i.e., musculus vs. domesticus and musculus vs. castaneus) standard One-Way ANOVA was performed for each tissue separately to identify differentially expressed transcripts. We applied the False Discovery Rate (FDR;

(REINER et al. 2003)) to the ANOVA P-values and only selected transcripts with an FDR

< 10%. Applying this stringent criterion to the musculus-domesticus comparison corresponds to a significance threshold of P < 0.00025 in the brain, P < 0.0000053 in the liver and P < 0.002 in the testis. In the musculus-castaneus comparison this corresponds to P < 0.002 in the testis; no transcript at all is significant in the brain and liver under a

10% False Discovery Rate (see Results). However, preliminary experiments using quantitative Real-Time-PCR of one gene differentially expressed between musculus and domesticus in the liver (Ces1) suggest that transcripts up to a P-value level of 0.06 on the microarray can be confirmed on independent samples (data not shown). Therefore, we repeated the selection of differentially expressed genes for each tissue using arbitrarily significance values set at P < 0.01 and P < 0.05.

Significant transcripts were considered separately according to whether they were expressed (>500) in both subspecies or just in one of the subspecies within a given comparison. Expression levels <500 likely include transcripts that are not expressed at all. In this case, additivity of transcript expression is meaningless. However, if the transcript is absent in one parental subspecies but present in the other, one can minimally ask whether the expression status in the hybrid is consistent with dominance or not. We used three methods to determine the mode of inheritance for the differentially expressed transcripts between the parental subspecies using the F1 hybrids between them. The musculus-domesticus comparison was processed in parallel to the musculus-castaneus comparison using identical methods. We define the average expression level across replicates of the subspecies with lower expression to be x, the average expression level of

the subspecies with higher expression to be y, and the expression in the hybrid to be hij, where the subscript i indicates the cross (i.e., either subspecies a female with subspecies b male or subspecies a male with subspecies b female; i = 1, 2) and the subscript j indicates the replicate within each cross (j = 1, 2). First, following GIBSON et al. (2004) we defined additivity as y / 1.25 > h > x * 1.25, for all h (note that y, x and h are measured on the untransformed scale). In this model, only transcripts with a difference between the parental subspecies of > 2.5 can show additivity, and hence the analysis is restricted to only these transcripts. Dominance is defined as 1.25 * x > h > x / 1.25 or 1.25 * y > h > y

/ 1.25 for all h, and overdominance as h < x / 1.25 or h > y * 1.25 for all h. We call this the “fixed threshold method.”

The problem with the fixed threshold method is that as the fold difference between two subspecies becomes greater, the zone where additivity will be accepted becomes larger. Accordingly, we devised a second measure, which we term the

“fractional threshold method,” based on the relative difference in expression levels between the parental subspecies. We defined additivity to be (x+y)/2 - ε < h < (x+y)/2 + ε, for all h, where ε was set to 5%, 10%, 15% and 20% of the difference between x and y

(e.g., ε = (y-x)*0.15). Dominance was set to be x - ε/2 < h < x + ε/2 or y - ε/2 < h < y +

ε/2 for all h, so that for each transcript the interval over which either dominance or additivity is accepted is equivalent. Overdominance was set to be outside of the range for dominance, i.e., h < x - ε/2 or h > y + ε/2 for all h. To be comparable to the analysis using a fixed threshold of 1.25, we only included genes where the expression difference (fold- change) between the parents was larger than 2.5-fold. For both methods we determined the mode of inheritance by (1) treating F1 hybrids from both directions of the cross together (i.e., all four replicates had to be consistent with a certain mode of inheritance) and (2) treating both directions of the cross separately to account for possible maternal effects, sex-linkage and directional dominance.

In addition to classifying transcripts into arbitrary mode-of-inheritance classes, we analyzed the distribution of dominance effects, which was calculated as d/a, where a is half the difference in expression level between the parental strains and d is the difference in expression level between the F1 hybrid (average across replicates independent of the direction of cross) and the average of the parental strains (FALCONER and MACKAY 1996;

GIBSON et al. 2004). If the expression level of the transcript in the hybrid is exactly intermediate between the parental strains d = 0. A d/a value of 0 corresponds to perfect within-locus additivity, |d/a| = 1 to complete dominance and |d/a| > 1 to overdominance.

We performed this analysis separately for each tissue (brain, liver and testis) and both comparisons (i.e., musculus-domesticus and musculus-castaneus) including all differentially expressed transcripts between parental strains (P<0.01 and P<0.05), with average expression levels >500 in at least one of the parental subspecies. Note that in contrast to the analyses where transcripts are assigned to discrete additive and non- additive classes we do not require that the differentially expressed transcript shows a fold-change > 2.5 between the parental subspecies, thus allowing us to more generally assess dominance effects.

2) Non-additive expression of non-differentially expressed transcripts

To assay parent – F1 hybrid offspring expression differences at the large number of transcripts not differentially expressed between the parental subspecies we used pairwise

Tukey post hoc tests. We selected transcripts in which F1 hybrids have significant (P <

0.05) and 1.25-fold higher or lower expression compared to both parents. We performed this analysis for F1 hybrids from both directions of the cross together as well as from each direction of the cross separately. For those transcripts that had expression levels <

500 in both the parental subspecies we only analyzed those crosses in which the F1 hybrids from at least one direction of the cross had average expression levels > 500.

Determination of subspecies origin of each probeset:

The probesets of the Mouse Genome 430 2.0 Affymetrix® GeneChip were designed based on the sequence of the laboratory inbred strain C57BL/6J. The genome of this strain is a mixture of genetic contributions from all three subspecies, domesticus, musculus and castaneus (WADE and DALY 2005; WADE et al. 2002; YANG et al. 2007).

According to the most recent analyses (YANG et al. 2007) up to 92% of the genome of laboratory inbred strains (and thus also the genome of C57BL/6J) is derived from domesticus, about 7% is from musculus and < 1% is from castaneus. Whenever DNA sequences show differences between these subspecies, we might expect probesets to hybridize more efficiently to that subspecies from which the particular genomic region was derived from. While sequence divergence between the subspecies is generally small

(typically 1% for non-coding DNA at autosomal loci (HARR 2006)) it could potentially enhance or alleviate expression differences. This problem is likely to be more prevalent when other subspecies are compared to domesticus, the subspecies that contributed most of the genetic material to C57BL/6J. Thus, we specifically address this problem for the comparison between domesticus and musculus. We used > 8 Million SNP (Single

Nucleotide Polymorphism) loci distributed all over the genome downloaded from the

Perlegene website (http://mouse.perlegen.com/mouse/download.html). These SNPs have been typed in 15 commonly used inbred strains of the laboratory house mouse. The data from these strains can be combined with the known sequence of strain C57BL/6J

(WATERSTON et al. 2002) at the respective SNP position. Two out of the 15 strains are so-called “wild-derived” strains, one of which belongs to musculus (PWD/PhJ) and one to domesticus (WSB/EiJ). We first used the SNPs that distinguished the musculus from the domesticus strain to assign the subspecies origin of each SNP in the C57BL/6J strain.

We then calculated the frequency of musculus-like or domesticus-like SNPs in a 20 kb region surrounding the location of the transcript that is targeted by a specific probe set. A probe set was called musculus-like if >60% of the SNPs in the region matched the musculus strain and domesticus-like if >60% of the SNPs matched the domesticus strain.

We then partitioned the data into two sets: transcripts that are located in “musculus-like” regions but show higher expression in domesticus and transcripts that are located in “domesticus-like” regions but show higher expression in musculus. For these transcripts, sequence divergence is unlikely to explain the expression differences. Due to the absence of a wild-derived castaneus strain in the Perlegen dataset we cannot assign regions in the genome of C57BL/6J that were contributed by this subspecies. Given that only a very small proportion of the genome derives from castaneus, this is probably not a very serious limitation.

Chromosomal location

The chromosomal location (according to NCBI Build 36) of each additively expressed transcript (identified by the fractional (15%) method in F1 hybrids from both directions of the cross (P < 0.05, fold-change between the parental subspecies >2.5) was obtained from the annotation file (“Mouse430_2.na21.annot.csv”) downloaded from the

Affymetrix website in January 2007 (www.affymetrix.com) and transformed into a

MySQL database. We counted the number of additively expressed transcripts on each . If two alternative transcripts of the same gene were additively expressed, we only counted it once. For each chromosome we also obtained the total number of known coding genes from Ensembl (www.ensembl.org). To test for differential representation of additively expressed transcripts on any particular autosome, we first calculated the ratio of additive transcripts to the total number of protein coding genes on this autosome. We then compared this to the ratio of the total number of additively expressed transcripts on all autosomes relative to the total number of protein coding genes on all autosomes, using a Fisher’s exact-test. To test for an overrepresentation of additively expressed transcripts on the X-chromosome we counted the number of additive transcripts on the X-chromosome relative to the total number of protein coding genes on the X-chromosome. This number is compared to the number of additive transcripts on all autosomes relative to the total number of protein coding genes on all autosomes. All analyses were repeated for each tissue separately.

Clustering of additively expressed transcripts along a chromosome

For each chromosome, we determined the chromosomal location of each additively expressed gene and sorted them according to chromosomal position (NCBI Build 36).

Next we calculated the distance between each gene and its nearest downstream neighbor and expressed this distance relative to the total length of that particular chromosome. We then asked whether this distance is smaller than expected by chance. The expectation was derived by randomly drawing two values from a uniform distribution. The absolute distance between these two numbers was calculated and the lower 2.5 percentile was derived from repeating this procedure 100,000 times. We asked whether the observed distance is smaller than the 2.5 percentile from the simulation.

Functional annotation

We used “PANTHER” (http://www.pantherdb.org/tools/genexAnalysis.jsp) to analyze the additively expressed transcripts in each tissue for overrepresentation (relative to the full gene content of the mouse genome) of certain biological functions. We performed this analysis for additively expressed transcripts identified by the fractional (15%) method in F1 hybrids from both directions of the cross (setting P < 0.05, > 2.5-fold difference between the parental subspecies). In addition we analyzed the transcripts differentially expressed between the parental subspecies independent of the mode of inheritance in F1 hybrids (again setting the same criteria P < 0.05, fold-change between parental subspecies > 2.5). Results

About 1/4 of all the transcripts represented on the Affymetrix® GeneChip have expression levels >500 in at least one of the three tissues from at least one of the three subspecies.

The number of expressed transcripts did not differ greatly among tissue (brain: 13,137 expressed transcripts, liver: 11,155 expressed transcripts, testis: 11,608 expressed transcripts). In Table 2 we show the pattern of differential transcript expression between musculus and domesticus and between musculus and castaneus respectively. At a significance level of P < 0.05, we find roughly similar numbers of transcripts in each of the three organs. Given the ~12,000 transcripts expressed in at least one subspecies and tissue, 600 should be assigned significantly different by chance at P < 0.05. Instead we find roughly three times as many significant transcripts in each tissue. The number of transcripts differentially expressed between the subspecies drops considerably when only transcripts with a larger than 2-fold change are considered (i.e., musculus-domesticus:

445 (brain), 683 (liver) and 847 (testis), musculus-castaneus: 435 (brain), 691 (liver) and

827 (testis)). However, even when parental expression differences are constrained to differ by 5-fold, we find a relatively large number of significantly different transcripts for each tissue (between 96 and 218 transcripts, dependent on tissue and subspecies comparison, Table 2).

If the significance is set to a more stringent level, a difference between the tissues becomes apparent (Figure 2). At the extreme, when the P-value level corresponds to a

False Discovery Rate of 10%, we find in the musculus-domesticus comparison 9 transcripts differentially expressed in the liver, one transcript differentially expressed in the brain and 218 transcripts differentially expressed in the testis. In the musculus- castaneus comparison significant transcripts were only found in the testis (213 transcripts). It is clear that there is an excess of differential transcript expression in the testis compared to the other organs (Figure 2).

For both comparisons, we also determined the overlap of differentially expressed transcripts between subspecies among the three different tissues. The vast majority of transcripts are differentially expressed in just one tissue within a comparison (Figure 3A), even though about half of the differentially expressed transcripts in each tissue are detectably expressed in at least one additional tissue (i.e., a large proportion of differentially expressed genes are not tissue specific, data not shown). Between 1/3 and

1/2 (depending on the tissue) of the transcripts that are differentially expressed in the musculus-domesticus comparison are also differentially expressed in the musculus- castaneus comparison (Figure 3B) suggesting that the expression phenotype at these transcripts is specific to musculus.

Proportion of differentially expressed genes showing additivity, dominance and overdominance

As described in the methods we used three different methods to assay the mode of inheritance for transcript expression, two of which assign transcripts into discrete classes

(i.e., additivity, dominance and overdominance): the fixed threshold method and the fractional threshold method. In Tables 3 and 4 we show the number of transcripts identified for each mode of inheritance according to these two methods, using an ε corresponding to 15% of the parental difference in expression for the fractional method.

We scored these numbers separately for the transcripts with expression levels >

500 in both parental strains and those that with an expression level >500 in only one of the two parents. For the latter group, we cannot unequivocally determine whether a transcript is additively expressed in the F1 hybrid because it may have no expression in one of the parental subspecies. For this reason we place numbers for additively expressed transcripts in parentheses in Tables 3 and 4.

When compared to the fixed threshold method (GIBSON et al. 2004) the fractional threshold method generally yields fewer transcripts that can be classified to one or other mode of inheritance. Nevertheless, for both the musculus-domesticus and musculus- castaneus comparisons, the fraction of additive versus nonadditive effects is almost identical for both methods, and changes little with varying values of ε in the fractional method. In all cases transcripts were classified into the same mode-of-inheritance categories by both methods.

The relative frequency of additive effects is generally independent of whether additivity was identified in F1 hybrids from only one direction of the cross or from both directions of the cross combined. The percentage of transcripts showing additivity was also similar across tissues. In all but one comparison > 80% of the transcripts are additively expressed, with about 10% being dominantly and 10% overdominantly expressed. The one exception to this pattern is in the musculus mother - castaneus father cross. In this cross, the testis (but not the other tissues) shows a much higher fraction of nonadditivity (~50% of all the assigned transcripts).

Similar results come from the analysis of transcripts with expression levels >500 in only one of the parental subspecies. Dominance and overdominance among these transcripts are rare, and expression levels of most transcripts are consistent with additivity. The fact that we find such a large fraction of transcripts consistent with additivity suggests that those transcripts with expression levels <500 may not be absent, but rather are expressed at a low level.

The distribution of dominance values is shown in Figure 4. Unlike the previous tests which only included genes showing a >2.5 fold difference in expression between the parental subspecies, this test includes all transcripts that are differentially expressed (at the given P-value threshold). The vast majority of d/a values fall within the -0.5 - +0.5 interval, again suggesting additivity to be predominant (recall that a value of d = 0 indicates exact intermediacy of expression).

Contribution of sequence divergence to expression change

Divergence in the sequences that represent the probes on the microarray may influence the results. The Affymetrix probe sequences are based on predominantly domesticus

DNA (i.e. the lab strain C57BL/6J (WADE and DALY 2005; WADE et al. 2002; YANG et al. 2007)) and we used this microarray to compare domesticus to musculus. In a direct comparison between domesticus and musculus, probesets should hybridize more efficiently in domesticus than in musculus in all regions of the genome where these two subspecies have accumulated sequence changes. To quantify the effect of sequence divergence of probesets for the musculus-domesticus comparison we obtained counts of the number of transcripts differentially expressed, separately for the cases where musculus shows higher expression than domesticus and where domesticus shows higher expression than musculus (Table 5). At low significance there is no difference between these numbers but at high significance about twice as many significant transcripts show higher expression in domesticus than in musculus. This suggests that sequence divergence might contribute to some extend to the detection of differentially expressed transcripts.

We performed two tests to determine whether such an effect could contribute to our estimate of the proportion of additively expressed transcripts. First, as described in the methods section, we assigned each probeset on the microarray as “musculus-like” or as “domesticus-like”. Then, we restricted the analysis to two sets of transcripts: 1) differentially expressed transcripts where the probe set is musculus-like but the expression change is in the direction of higher expression in domesticus, and 2) differentially expressed transcripts where the probe set is domesticus-like but the expression change is in the direction of higher expression in musculus. In both cases, sequence divergence in the probe is not likely to explain the measured expression divergence. As shown in Table 6, the proportions of additively expressed transcripts among the differentially expressed transcripts in this subset of genes do not differ from the analysis including all probe sets. As a second test, we considered the reciprocal crosses between musculus and castaneus. These two subspecies are more closely related to each other than either is to domesticus (PRAGER et al. 1998). Thus, this comparison should suffer much less from an ascertainment bias because both of these subspecies are equally diverged from the subspecies that contributes most sequence to the probe (i.e. domesticus). As shown above

(Table 4), with one exception (one direction of the musculus-castaneus cross, hybrid testis), we can confirm the large surplus of additively expressed transcripts. Thus, we conclude that the excess of additive transcript expression is not substantially affected by divergence in the sequences that correspond to the probes on the Affymetrix microarray.

We also analyzed all differentially expressed transcripts between the parental subspecies with respect to maternal or paternal effects. To identify these transcripts we analyzed both directions of the cross separately and looked for transcripts that were dominantly expressed in both directions of the cross. We classified transcripts as maternal effect transcripts if the expression level in the F1 hybrids would always match that one of the mother, while for paternal effects transcripts F1 hybrids would always match their father. This analysis was restricted to transcripts that are not sex-linked.

Interestingly we found one region on chromosome 7 (around 59 Mb, NCBI Build 36, gene Snrnp) that seems to show a paternal effect in the brain of F1 hybrids from the musculus-domesticus cross as well as in F1 hybrids from the musculus-castaneus cross.

The former cross does show the same paternal effect also in the liver. No maternal effect transcripts were found. Nonadditive expression in F1 hybrids at transcripts that are not differentially expressed between the parents

We included in this analysis all transcripts that did not show a significant difference in expression between the parental subspecies (i.e., P > 0.05) (see Figure 1C). However, if the transcript had an expression level < 500 in both parental subspecies for a given comparison, the transcript had to have an expression level >500 in the respective F1 hybrid to be included. Among the ~10,000 transcripts meeting the criteria for this analysis (i.e. ~12,000 minus the ~2,000 transcripts that are differentially expressed between two subspecies at P < 0.05) we found a few (between 0.01% and 0.5%) that differed significantly in expression between parents and F1 hybrids (Table 7). There is again the familiar exception, i.e., the musculus mother – castaneus father cross in the testis, where ~1/3 of all transcripts show significant nonadditivity.

Thus, apart from the one exception, also in this analysis, nonadditivity of transcript expression is rare.

Chromosomal distribution of additively expressed transcripts

In Supplementary Table 8 we list the chromosomal location (based on NCBI Build 36) of additively expressed transcripts in F1 hybrids that are found in both directions of the cross. For musculus-domesticus a significant overrepresentation of additively expressed transcripts was found for chromosome 13 in the testis (Fisher’s exact test, P=0.01) and for chromosome 12 in the liver (P=0.014). In the brain, no chromosome was significantly overrepresented. The musculus-castaneus comparison yielded a significant overrepresentation of additively expressed genes only for chromosome 10 in the brain

(P=0.03). The significant effect of all disappeared when P-values were

Bonferroni corrected. We found either no transcripts (musculus-domesticus) or only one transcript (musculus-castaneus) additively expressed in F1 hybrids on the X- chromosome; although these are small numbers they are not significantly different from expectation (P>0.05). Due to small sample sizes, we did not test transcripts showing dominance or overdominance.

Clustering of additively expressed transcripts

We analyzed the additively expressed transcripts for chromosomal clustering. Significant clustering might be expected if genes that are involved in the same pathway co-localize along the chromosome (HURST et al. 2004), or alternatively if differential expression is caused by deletion or duplication of a region encompassing more than one gene. We found a few instances of significant clustering of differentially expressed transcripts that show additivity. For the musculus-domesticus comparison we found chromosome 17

(Glo1, Glpr, P=0.008) and chromosome 5 (Camkk2, Tmed2, P=0.024) in the brain, chromosome 13 (Akr1c18, Akr1c12, P=0.0006) and chromosome 7 (Aldoa, Hirip3,

P=0.0006) in the liver and chromosome 2 (Rims2, 4930429F24Rik, P=0.001) and chromosome 3 (Gm440, Tm2d2, P=0.0012) in the testis. For the musculus-castaneus comparison only the brain (chromosome 10 (un-annotated gene and Fuca1, P=0.005)) and liver (chromosome 6 (1700019G17Rik and 2810474O19Rik, P=0.0005)) showed evidence for clustering of additively expressed transcripts. Functional categories of transcripts

For brain and testis no overrepresentation of any biological process was detected among the additively expressed transcripts. In the liver, metabolic functions were significantly over-represented among the additively expressed transcripts in musculus-domesticus hybrids and oxygenase and acetyltransferase functions were over-represented in musculus-castaneus hybrids (both assessed after Bonferroni correction).

We also performed this analysis for all transcripts differentially expressed between the parental subspecies, ignoring the expression status in F1 hybrids. As shown in Supplementary Table 9, we found a significant overrepresentation of certain biological processes in all tissues and in both comparisons of subspecies. Comparing musculus with domesticus, similar biological processes were overrepresented in the brain and testis

(these included transcripts involved in intracellular protein traffic), and in the liver we found a strong and highly significant overrepresentation of transcripts involved in lipid, fatty acid and steroid metabolism. The musculus-castaneus comparison showed similar patterns in the liver, but brain and testis showed different and distinct patterns of overrepresentation (i.e. Carbohydrate metabolism (brain) and detoxification (testis),

Supplementary Table 10).

Finally, we functionally annotated the large number of nonadditively expressed transcripts that did not show an expression difference between the parental subspecies in the testis from the musculus mother - castaneus father direction of the cross (see Supplementary Table 11). Various processes seem to be affected. “Spermatogenesis and motility” and “gametogenesis” rank fifth and sixth in the list with an overrepresentation

P-value of 5 x 10-6. Discussion

We assayed expression levels in subspecies of the house mouse and their reciprocal F1 hybrids in three different tissues: brain, liver and testis. These tissues were chosen as representative of different aspects of the phenotype of an organism. Previously we suggested that expression changes in the brain between the subspecies could reflect behavioral differences, expression changes in the liver could indicate general metabolic differences, while changes in the testis pertain to changes in reproduction (VOOLSTRA et al. 2007). We found a strikingly high preponderance of transcripts additively expressed in

F1 hybrids between the subspecies in all three tissues. Additively expressed transcripts were found more or less evenly distributed among all the chromosomes of the house mouse and functional annotation of the transcripts did not reveal an overrepresentation of specific functional categories. While we did see an overrepresentation of metabolic functions for additively expressed transcripts in the liver, the same overrepresentation is also seen when all transcripts differentially expressed between the parents are functionally annotated, suggesting that this is liver-specific rather than a phenomenon associated with additive expression.

We used three different methods to identify additivity and nonadditivity in transcript expression, two of which classify transcript expression into discrete mode-of inheritance categories, i.e. the fixed and the fractional threshold method. Both these methods have different forms of bias. The fixed threshold method biases the result towards detecting more additivity as the difference in expression between the parental strains rises, whereas the fractional method is biased towards more overdominance as the differences between the parents get smaller. Because both methods affect the result in different directions, yet both of them yield the same result, we find strong support for the notion that additive transcript expression is common. The third method that assesses additivity/nonadditivity in transcript expression level, using a continuous distribution of dominance effects, was also consistent with our interpretation of most expression changes being additively expressed (Figure 4).

Among the seven studies in addition to ours that so far have addressed this question rigorously, two found evidence for substantial nonadditivity in gene expression

(GIBSON et al. 2004; HEDGECOCK et al. 2007) while four studies match ours, showing a strong preponderance of additivity (Drosophila (HUGHES et al. 2006), maize (STUPAR and

SPRINGER 2006; SWANSON-WAGNER et al. 2006), the (CUI et al.

2006)). One study in Arabidopsis (VUYLSTEKE et al. 2005) found relatively similar numbers of additively and nonadditively expressed transcripts. Several factors could explain the discrepancy among studies: 1) Degree of divergence might play a critical role.

HUGHES et al. (2006) suggested that more recently diverged lines could show a higher fraction of additivity. However, one might expect that as more and more genes come to influence a trait, provided those genes have no net directional dominance, additivity of expression will increase with time, at least up to a certain threshold. 2) The fraction of additivity might depend on whether or not the parental lines are inbred. Given the fact that many crosses from highly inbred lines exhibit heterosis (i.e., F1 hybrid fitness- related traits lie outside the range of the two parents) one might expect nonadditively expressed transcripts causally relate to the heterosis phenotype in studies using inbred lines. 3) There could be methodological reasons for the discrepancy among these studies.

We pursue this possibility in more detail in the next paragraphs.

Several methodological differences between studies are relevant. First, the definition of additivity depends on the measurement scale, i.e., on some scales an expression level might be called additive and on others not. In our study of house mouse subspecies we have tried to minimize such an effect by using three different analysis methods, all of which gave very similar results and supported the view of pervasive additivity in gene expression in F1 hybrid mice. In addition, while we used only a small number of replicates in our study, our result is very robust with respect to changes in P- value and fold-change levels (data not shown) and thus is unlikely explained by the measurement scale.

Second, the outcome of a study might depend on the biological material that is used for RNA extractions. Some studies use single organs (such as the study of laboratory strains of the house mouse and our study) while other studies pool different tissues (i.e. a whole plant, or a whole fly is processed). Third, studies might differ due to technical reasons, most notable microarray platform. Recently there has been an interest in comparing the performance of different methods to measure gene expression level on the genome-wide scale. While some studies did not find difference in performance among platforms (PATTERSON et al. 2006), other studies conclude that generally one- color microarrays (such as Affmetrix or ABI) outperform two-color microarrays (such as homemade cDNA-arrays or Agilent arrays) (DE REYNIES et al. 2006; KUO et al. 2006). In Table 12 we have summarized all studies that have systematically assayed patterns of inheritance in gene expression. We have also included four additional studies that have less rigorously treated the issue, but these studies allow a qualitative evaluation.

As shown in the table, neither the level of inbreeding nor differences in divergence time are associated with the frequency of additivity. This is for example illustrated by the study by RANZ et al. (2004) who compared Drosophila melanogaster to Drosophila simulans and their female F1 hybrids. Using heads as source for the RNA, they found a high proportion of additivity, but the whole body of these flies showed a surplus of nonadditivity. In this example, complexity of the tissue used as RNA source might explain the difference, although this cannot apply in other cases, such as the study by

HUGHES et al. (2006) who also used whole flies as RNA source and found a strong preponderance of additivity (Table 12). The most striking correlate of additivity is the platform used. Among the eight studies that support additivity as the major mode of inheritance, five used Affymetrix arrays, two were based on cDNA arrays and one study used differential display. None of the four studies that found a preponderance of nonadditive expression used the Affymetrix platform but instead used either 2-color microarray or sequencing-based expression measures. Thus, the combined evidence seem to suggest that methodological issues are a main cause of different results. This implies that additivity may be a general feature of most (but not all) genes differentially expressed between taxa, independent of divergence level or inbreeding, and most likely also independent of complexity of the tissue. Expression changes and reproductive isolation

House mice are increasingly becoming a model system to study the genetics of incipient speciation (GALTIER et al. 2004; GUENET and BONHOMME 2003; HARR 2006).

Under the biological species concept (MAYR 1963) the different taxa of house mice have been considered subspecies (BOURSOT et al. 1993) because they are only partially reproductively isolated from each other (i. e., male and female hybrids are fully viable and only some crosses between subspecies show hybrid male sterility, while females are fully fertile). VOOLSTRA et al. (2007) have previously tested whether the evolution of reproductive isolation between house mouse subspecies could be explained by an exceptionally fast rate of divergence in gene expression between the subspecies in reproductive organs (i.e. testis), as predicted by some models of sexual selection and sexual conflict (RICE 1998). While VOOLSTRA et al. (2007) found a large number of genes differing in expression in the testis between two species of Mus (Mus spretus and

Mus musculus), comparisons between house mouse subspecies showed most differences in expression in the combined liver and kidney tissue and only very few genes differentially expressed in the testis. In contrast, we find here that most gene expression changes between subspecies of M. musculus are in the testis, especially at high significance threshold. Some reasons that could explain the discrepancy between the current study and VOOLSTRA et al. (2007) are methodological. VOOLSTRA et al. (2007) used two-color microarrays while the current study uses Affymetrix microarrays.

VOOLSTRA et al. (2007) also used individuals that were collected directly in the wild and standardized to laboratory conditions for a few days only, while all mice used in this study have been born in the laboratory. Finally, VOOLSTRA et al. (2007) used multiple unrelated individuals within each subspecies as replicates in the statistical tests, whereas we used only a single (and to a large extent inbred) strain within each subspecies. It is possible that genes affecting hybrid sterility are segregating within populations. If different individuals in Voolstra et al.’s sample were of different genotypes, the consequent within-subspecies variability would obscure between subspecies differences.

In support of this explanation are experiments by VYSKOCILOVA et al. (2005), who crossed laboratory strain C57BL/6J to several different musculus strains from various locations in the Czech Republic and found several hybrid sterility loci segregating within natural populations. Interestingly, the only other study we are aware of that has mapped hybrid sterility loci in recently diverged taxa (in two species of Mimulus) with some precision and determined the frequency of these alleles in natural populations also found several hybrid sterility loci segregating within natural populations but no fixed differences between the species (SWEIGART et al. 2007).

We have not measured the fertility status of the F1 hybrid males we studied.

Some may be sterile. In one cross, between a musculus mother and a castaneus father, an extraordinary large number of transcripts were “mis-expressed” in F1 hybrids relative to both parents in the testis. The fact that this high number of mis-expressed transcripts is restricted to the testis and is not observed in brain or liver rules out methodological/technical problems to have caused this mis-expression. If this is linked to sterility, it indicates that a very large number of genes are affected, perhaps a result of a few master control genes influencing many downstream targets. The observed case of pervasive mis-expression in the testis could represent the crossing of a threshold whereby genetic incompatibilities between taxa result in low hybrid fitness. Such incompatibilities are often thought to affect fertility before viability (COYNE and ORR 2004).

It is less clear how the relatively large number of differentially expressed transcripts between the parental strains that are additively expressed in F1 hybrids would relate to F1 hybrid fitness. While intermediate expression level might already be deleterious in F1 hybrids these genes might also be important isolating genes in F2 or later generation backcrosses, a phenomenon that has been observed in some between subspecies crossing schemes (OKA et al. 2007). While the contribution of the identified genes to reproductive isolation at this stage is highly speculative, our study highlights valuable candidate genes (both, the additive and the non-additively expressed transcripts) on which future expression, phenotypic and population genetic studies using larger samples of multiple unrelated individuals within each subspecies could be targeted.

Ultimately, these efforts will help to determine the genetic basis of reproductive isolation in the house mouse in nature.

Acknowledgements

We are especially grateful to J. Pialek for providing the mouse strains of M. m. musculus and M. m. domesticus and A. Orth and F. Bonhomme for providing M. m. castaneus. B.

Schmitz contributed expert technical help and C. Voolstra helpful experimental advice. P.

Nuernberg and C. Becker from the Center for Functional Genomics in Cologne performed the microarray experiments. We would like to thank T. Price for valuable discussion and J. Baines and R. Scavetta and two anonymous reviewers for helpful comments on the manuscript. The work has been supported by Emmy-Noether and

SFB680 funds by the DFG to B. Harr. Literature

BOURSOT, P., J. C. AUFFRAY, J. BRITTONDAVIDIAN and F. BONHOMME, 1993 The

Evolution of House Mice. Annual Review of Ecology and Systematics 24: 119-

152.

BREM, R. B., G. YVERT, R. CLINTON and L. KRUGLYAK, 2002 Genetic dissection of

transcriptional regulation in budding yeast. Science 296: 752-755.

COMINGS, D. E., and J. P. MACMURRAY, 2000 Molecular heterosis: a review. Mol Genet

Metab 71: 19-31.

COYNE, J. A., and H. A. ORR, 2004 Speciation. Sinauer Associates, Sunderland, MA.

CUI, X., J. AFFOURTIT, K. R. SHOCKLEY, Y. WOO and G. A. CHURCHILL, 2006

Inheritance patterns of transcript levels in F1 hybrid mice. Genetics 174: 627-637.

DE REYNIES, A., D. GEROMIN, J. M. CAYUELA, F. PETEL, P. DESSEN et al., 2006

Comparison of the latest commercial short and long oligonucleotide microarray

technologies. BMC Genomics 7: 51.

DOSS, S., E. E. SCHADT, T. A. DRAKE and A. J. LUSIS, 2005 Cis-acting expression

quantitative trait loci in mice. Genome Res 15: 681-691.

FALCONER, D. S., and T. F. C. MACKAY, 1996 Introduction to Quantitative Genetics.

Longman, Essex, UK. GALTIER, N., F. BONHOMME, C. MOULIA, K. BELKHIR, P. CAMINADE et al., 2004 Mouse

biodiversity in the genomic era. Cytogenet Genome Res 105: 385-394.

GIBSON, G., R. RILEY-BERGER, L. HARSHMAN, A. KOPP, S. VACHA et al., 2004 Extensive

sex-specific nonadditivity of gene expression in Drosophila melanogaster.

Genetics 167: 1791-1799.

GUENET, J. L., and F. BONHOMME, 2003 Wild mice: an ever-increasing contribution to a

popular mammalian model. Trends Genet 19: 24-31.

HARR, B., 2006 Genomic islands of differentiation between house mouse subspecies.

Genome Res 16: 730-737.

HEDGECOCK, D., J. Z. LIN, S. DECOLA, C. D. HAUDENSCHILD, E. MEYER et al., 2007

Transcriptomic analysis of growth heterosis in larval Pacific oysters (Crassostrea

gigas). Proc Natl Acad Sci U S A 104: 2313-2318.

HEY, J., and R. M. KLIMAN, 1993 Population genetics and phylogenetics of DNA

sequence variation at multiple loci within the Drosophila melanogaster species

complex. Molecular Biology and Evolution 10: 804-822.

HOEKSTRA, H. E., and J. A. COYNE, 2007 The locus of evolution: evo devo and the

genetics of adaptation. Evolution Int J Org Evolution 61: 995-1016.

HUGHES, K. A., J. F. AYROLES, M. M. REEDY, J. M. DRNEVICH, K. C. ROWE et al., 2006

Segregating variation in the transcriptome: cis regulation and additivity of effects.

Genetics 173: 1347-1355. HURST, L. D., C. PAL and M. J. LERCHER, 2004 The evolutionary dynamics of eukaryotic

gene order. Nat Rev Genet 5: 299-310.

KARN, R. C., A. ORTH, F. BONHOMME and P. BOURSOT, 2002 The complex history of a

gene proposed to participate in a sexual isolation mechanism in house mice. Mol

Biol Evol 19: 462-471.

KUO, W. P., F. LIU, J. TRIMARCHI, C. PUNZO, M. LOMBARDI et al., 2006 A sequence-

oriented comparison of gene expression measurements across different

hybridization-based technologies. Nat Biotechnol 24: 832-840.

MACKAY, T. F., 2001 The genetic architecture of quantitative traits. Annu Rev Genet 35:

303-339.

MAYR, E., 1963 Animal Species and Evolution. The Belknap press, Cambridge, MA.

MICHALAK, P., and M. A. NOOR, 2003 Genome-wide patterns of expression in

Drosophila pure species and hybrid males. Mol Biol Evol 20: 1070-1076.

OKA, A., T. AOTO, Y. TOTSUKA, R. TAKAHASHI, M. UEDA et al., 2007 Disruption of

genetic interaction between two autosomal regions and the x chromosome causes

reproductive isolation between mouse strains derived from different subspecies.

Genetics 175: 185-197.

PATTERSON, T. A., E. K. LOBENHOFER, S. B. FULMER-SMENTEK, P. J. COLLINS, T. M.

CHU et al., 2006 Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project. Nat Biotechnol 24:

1140-1150.

PRAGER, E. M., C. ORREGO and R. D. SAGE, 1998 Genetic variation and phylogeography

of central Asian and other house mice, including a major new mitochondrial

lineage in Yemen. Genetics 150: 835-861

RANZ, J. M., K. NAMGYAL, G. GIBSON and D. L. HARTL, 2004 Anomalies in the

expression profile of interspecific hybrids of Drosophila melanogaster and

Drosophila simulans. Genome Res 14: 373-379.

REILAND, J., and M. A. NOOR, 2002 Little qualitative RNA misexpression in sterile male

F1 hybrids of Drosophila pseudoobscura and D. persimilis. BMC Evol Biol 2: 16.

REINER, A., D. YEKUTIELI and Y. BENJAMINI, 2003 Identifying differentially expressed

genes using false discovery rate controlling procedures. Bioinformatics 19: 368-

375.

RICE, W. R., 1998 Intergenomic conflict, interlocus antagonistic coevolution, and the

evolution of reproductive isolation, pp. 261-270 in Endless Forms: Species and

Speciation, edited by S. BERLOCHER. Oxford University Press, New York,

Oxford, Oxford.

ROCKMAN, M. V., and L. KRUGLYAK, 2006 Genetics of global gene expression. Nat Rev

Genet 7: 862-872. SCHADT, E. E., S. A. MONKS, T. A. DRAKE, A. J. LUSIS, N. CHE et al., 2003 Genetics of

gene expression surveyed in maize, mouse and man. Nature 422: 297-302.

STOREY, J. D., J. M. AKEY and L. KRUGLYAK, 2005 Multiple locus linkage analysis of

genomewide expression in yeast. PLoS Biol 3: e267.

STUPAR, R. M., and N. M. SPRINGER, 2006 Cis-transcriptional variation in maize inbred

lines B73 and Mo17 leads to additive expression patterns in the F1 hybrid.

Genetics 173: 2199-2210.

SU, A. I., T. WILTSHIRE, S. BATALOV, H. LAPP, K. A. CHING et al., 2004 A gene atlas of

the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A

101: 6062-6067.

SWANSON-WAGNER, R. A., Y. JIA, R. DECOOK, L. A. BORSUK, D. NETTLETON et al.,

2006 All possible modes of gene action are observed in a global comparison of

gene expression in a maize F1 hybrid and its inbred parents. Proc Natl Acad Sci U

S A 103: 6805-6810.

SWEIGART, A. L., A. R. MASON and J. H. WILLIS, 2007 Natural variation for a hybrid

incompatibility between two species of Mimulus. Evolution Int J Org Evolution

61: 141-151.

VOOLSTRA, C., D. TAUTZ, P. FARBROTHER, L. EICHINGER and B. HARR, 2007

Contrasting evolution of expression differences in the testis between species and

subspecies of the house mouse. Genome Res 17: 42-49. VUYLSTEKE, M., F. VAN EEUWIJK, P. VAN HUMMELEN, M. KUIPER and M. ZABEAU, 2005

Genetic analysis of variation in gene expression in Arabidopsis thaliana. Genetics

171: 1267-1275.

VYSKOCILOVA, M., Z. TRACHTULEC, J. FOREJT and J. PIALEK, 2005 Does geography

matter in hybrid sterility in house mice? Biol J Linn Soc 84: 663-674.

WADE, C. M., and M. J. DALY, 2005 Genetic variation in laboratory mice. Nat Genet 37:

1175-1180.

WADE, C. M., E. J. KULBOKAS, 3RD, A. W. KIRBY, M. C. ZODY, J. C. MULLIKIN et al.,

2002 The mosaic structure of variation in the laboratory mouse genome. Nature

420: 574-578.

WATERSTON, R. H., K. LINDBLAD-TOH, E. BIRNEY and E. AL., 2002 Initial sequencing

and comparative analysis of the mouse genome. Nature 420: 520-562.

WRAY, G. A., 2007 The evolutionary significance of cis-regulatory mutations. Nat Rev

Genet 8: 206-216.

YANG, H., T. A. BELL, G. A. CHURCHILL and F. PARDO-MANUEL DE VILLENA, 2007 On

the subspecific origin of the laboratory mouse. Nat Genet 39: 1100-1107. Figure legends

Figure 1

Schematic representation of within locus modes of gene action in the comparison between house mouse subspecies and their F1 hybrids. Sp = Subspecies.

Figure 2

Number of differentially expressed transcripts between M. m. musculus and M. m. domesticus and between M. m. musculus and M. m. castaneus in different tissues and at different P-value thresholds.

Figure 3

A: Venn diagram showing the overlap of transcripts differentially expressed between parental subspecies (P < 0.01, fold-change >2.5) among the different tissues.

B: Venn diagram showing the overlap of differentially expressed transcripts between the musculus-domesticus and the musculus-castaneus comparison within a given tissue

(ANOVA P < 0.01, fold-change >2.5 within each comparison). Figure 4

Distribution of dominance values in F1 hybrids for differentially expressed transcripts between M. m. musculus and M. m. domesticus and M. m. musculus and M. m. domesticus for two different P-value thresholds. Table 1

Overview of the samples hybridized separately on the microarray

Cross name Subspecies involved tissue No. of individuals

domesticus Brain, Liver, 2 males Testis

musculus Brain, Liver, 2 males Testis

castaneus Brain, Liver, 2 males Testis dom-mus F1 hybrid Brain, Liver, 2 males (domesticus female x musculus male) Testis mus-dom F1 hybrid Brain, Liver, 2 males (musculus female x domesticus male) Testis mus-cas F1 hybrid Brain, Liver, 2 males (musculus female x castaneus male) Testis cas-mus F1 hybrid Brain, Liver, 2 males (castaneus female x musculus male) Testis Table 2

Number of transcripts differentially expressed in an ANOVA analysis between the parental subspecies, M. m. musculus and M. m. domesticus and between M. m. musculus and M. m. castaneus at different ANOVA P-value thresholds and at different magnitudes of change.

ANOVA P<0.05 ANOVA P<0.01 ANOVA P<10% FDR

musculus musculus musculus musculus musculus musculus vs. vs. vs. vs. vs. vs. domesticus castaneus domesticus castaneus domesticus castaneus

Brain all 1943 1819 573 516 1 0

2-fold 445 435 208 192 0 0

5-fold 96 98 63 54 0 0

Liver all 1777 1808 490 463 9 0

2-fold 683 691 249 228 7 0

5-fold 179 185 90 85 3 0

Testis all 2563 2514 882 822 218 213

2-fold 847 827 431 401 128 130

5-fold 218 197 151 132 55 61 Table 3

Number of additively, dominantly and overdominantly expressed transcripts among the differentially (fold-change >2.5, ANOVA P

<0.01) expressed transcripts between M. m. musculus and M. m. domesticus in different tissues.

mus-dom mus-dom mus-dom AND AND AND dom-mus dom-mus dom-mus mus-dom mus-dom mus-dom dom-mus dom-mus dom-mus

Additive Dominant Overdominant Additive Dominant Overdominant Additive Dominant Overdominant Fixed threshold, expressed transcripts in domesticus AND musculus Brain 39 5 0 41 5 0 31 1 0 Liver 104 17 2 113 10 1 84 1 0 Testis 97 11 1 104 10 0 80 4 0 Fixed threshold, expressed transcripts in domesticus OR musculus Brain (158) 6 0 (156) 13 0 (132) 0 0 Liver (219) 24 2 (200) 21 2 (164) 5 0 Testis (342) 12 0 (334) 22 2 (293) 4 0 Fractional threshold (15%), expressed transcripts in domesticus AND musculus Brain 27 0 0 20 0 0 13 0 0 Liver 57 6 2 56 2 1 24 0 0 Testis 58 0 1 66 1 0 39 0 0 Fractional threshold (15%), expressed transcripts in domesticus OR musculus Brain (90) 4 1 (87) 3 0 (40) 0 0 Liver (89) 6 7 (74) 9 5 (41) 2 2 Testis (172) 2 1 (163) 1 3 (87) 0 0 Numbers in brackets indicate that additive expression of a transcript cannot unequivocally be identified as the transcript is not expressed in one of the subspecies. Table 4

Number of additively, dominantly and overdominantly expressed transcripts among the differentially (fold-change >2.5, ANOVA P

<0.01) expressed transcripts between M. m. musculus and M. m. castaneus in different tissues.

mus-cas mus-cas mus-cas AND AND AND cas-mus cas-mus cas-mus mus-cas mus-cas mus-cas dom-mus dom-mus dom-mus

Additive Dominant Overdominant Additive Dominant Overdominant Additive Dominant Overdominant Fixed threshold, expressed transcripts in castaneus AND musculus Brain 21 1 0 22 0 1 17 0 0 Liver 49 5 0 35 3 0 32 2 0 Testis 43 7 0 31 18 4 21 3 0 Fixed threshold, expressed transcripts in castaneus OR musculus Brain (110) 5 0 (93) 6 1 (83) 1 0 Liver (129) 14 0 (116) 8 0 (105) 3 0 Testis (223) 16 1 (153) 56 10 (126) 5 0 Fractional threshold (15%), expressed transcripts in castaneus AND musculus Brain 10 0 0 13 0 1 4 0 0 Liver 11 1 0 11 1 0 4 0 0 Testis 19 2 0 9 8 10 6 0 0 Fractional threshold (15%), expressed transcripts in castaneus OR musculus Brain (55) 1 1 (46) 1 2 (26) 0 0 Liver (43) 2 0 (42) 4 0 (21) 0 0 Testis (123) 5 1 (38) 22 20 (21) 0 0 Numbers in brackets indicate that additive expression of a transcript cannot unequivocally be identified as the transcript is not expressed in one of the subspecies. Table 5

Number of transcripts differentially expressed between the parental subspecies, M. m. musculus and M. m. domesticus at different P-value thresholds and at different magnitudes of change, separately for the cases where domesticus shows higher expression than musculus (domesticus > musculus) and vice versa (musculus > domesticus).

musculus domesticus musculus domesticus musculus domesticus > > > > > > domesticus musculus domesticus musculus domesticus musculus ANOVA P< 10% ANOVA P<0.05 ANOVA P<0.01 FDR

Brain all 956 987 267 306 1 0 2-fold 159 286 69 139 0 0 5-fold 29 67 21 42 0 0

Liver all 913 864 253 237 3 6 2-fold 315 368 112 137 2 5 5-fold 81 98 42 48 2 1

Testis all 1269 1294 417 465 87 131 2-fold 331 516 164 267 41 87 5-fold 76 142 56 95 17 38 Table 6

Number of additively, dominantly and overdominantly expressed transcripts among the differentially (fold-change >2.5, ANOVA P

<0.01) expressed transcripts between M. m. musculus and M. m. domesticus in different tissues, considering the

“musculus”/”domesticus”-likeness of the probe set, excluding sequence divergence as a contributor to expression differences.

mus-dom mus-dom mus-dom AND AND AND dom-mus dom-mus dom-mus mus-dom mus-dom mus-dom dom-mus dom-mus dom-mus

Additive Dominant Overdominant Additive Dominant Overdominant Additive Dominant Overdominant Fixed threshold, expressed transcripts in domesticus AND musculus Brain 14 1 0 14 2 0 12 1 0 Liver 42 7 0 45 1 0 36 0 0 Testis 33 1 0 33 2 0 28 1 0 Fixed threshold, expressed transcripts in domesticus OR musculus Brain (48) 1 0 (50) 4 0 (43) 0 0 Liver (90) 11 0 (85) 9 2 (70) 3 0 Testis (128) 2 0 (128) 6 0 (114) 0 0 Fractional threshold (15%), expressed transcripts in domesticus AND musculus Brain 9 0 0 9 0 0 6 0 0 Liver 20 3 0 18 1 0 7 0 0 Testis 16 0 0 19 0 0 11 0 0 Fractional threshold (15%), expressed transcripts in domesticus OR musculus Brain (29) 1 0 (28) 1 0 (13) 0 0 Liver (29) 3 3 (25) 3 4 (14) 1 2 Testis (52) 0 0 (55) 0 0 (25) 0 0 Numbers in brackets indicate that additive expression of a transcript cannot unequivocally be identified as the transcript is not expressed in one of the subspecies. Table 7

Number of transcripts nonadditively expressed in F1 hybrids compared to the parental subspecies

dom-mus AND dom-mus mus-dom mus-dom Brain 7 19 1 Liver 23 44 5 Testis 15 17 3

mus-cas AND mus-cas cas-mus cas-mus Brain 23 3 2 Liver 15 33 7 Testis 2798 12 16 Supplementary Table 8

Number of additively expressed transcripts (identified by the fractional threshold method) on different chromosomes. Counted are only transcripts that differ significantly in expression (ANOVA P<0.05) between the parental strains and where the parental strains show a 2.5-fold difference in expression level. MD: musculus-domesticus comparison:

MC: musculus-castaneus comparison.

No. of No. of No. of No. of known transcripts$ in transcripts$ in transcripts$ in protein coding brain liver testis chromosome genes MD / MC MD / MC MD / MC 1 1275 3 / 4 3 / 4 4 / 0 2 1903 3 / 2 5 / 2 5 / 3 3 1076 2 / 2 3 / 3 5 / 4 4 1350 4 / 2 3 / 1 2 / 1 5 1322 4 / 2 2 / 2 4 / 2 6 1189 1 / 4 4 / 4 0 / 2 7 1982 4 / 1 5 / 0 3 / 0 8 1110 2 / 3 2 / 2 5 / 2 9 1277 3 / 2 3 / 1 2 / 3 10 1044 4 / 6* 3 / 2 6 / 3 11 1713 5 / 4 2 / 4 4 / 1 12 715 3 / 4 6* / 1 3 / 1 13 881 2 / 2 2 / 2 8* / 1 14 845 1 / 1 0 / 0 4 / 1 15 839 0 / 0 3 / 2 3 / 1 16 699 2 / 1 1 / 0 0 / 1 17 1070 4 / 0 6 / 4 4 / 1 18 551 1 / 1 1 / 2 4 / 0 19 758 3 / 4 1 / 3 2 / 2 Sum 21599 51 / 45 55 / 39 68 / 29

X 906 0 / 0 0 / 0 0 / 1 * significant surplus of transcripts (Fisher’s exact test)

$ a single transcript per gene is used Supplementary Table 9

Overrepresented biological processes among the differentially expressed transcripts between M. m. musculus and M. m. domesticus (P<0.05, fold-change>2.5)

Observed Expected number of number of genes in genes in category category P-value* Brain Intracellular protein traffic 18 7.01 0.00849 Protein modification 22 8.73 0.01070

Liver Lipid, fatty acid and steroid metabolism 39 10.55 <0.00001 Steroid metabolism 23 3.37 <0.00001 Cholesterol metabolism 11 0.97 <0.00001 Electron transport 16 4.3 0.00030 Other metabolism 22 7.52 0.00031 Detoxification 8 1.13 0.00329 Proteolysis 31 13.81 0.00446 Protein metabolism and modification 69 45.83 0.01060 Coenzyme and prosthetic group metabolism 9 2.14 0.01110 Immunity and defense 42 24.16 0.01190 Carbohydrate metabolism 18 7.3 0.01540

Testis Protein targeting and localization 11 3.12 0.01170 Intracellular protein traffic 29 15.01 0.02120 * after Bonferroni correction Supplementary Table 10

Overrepresented biological processes among the differentially expressed transcripts between M. m. musculus and M. m. castaneus (P<0.05, fold-change>2.5)

Observed Expected number of number of genes in genes in category category P-value * Brain Carbohydrate metabolism 16 4.29 0.000257 Other polysaccharide metabolism 7 1.05 0.015200

Liver Lipid, fatty acid and steroid metabolism 39 11.11 <0.000001 Steroid metabolism 18 3.55 0.000005 Other transport 8 0.76 0.000189 Other metabolism 23 7.92 0.000225 Immunity and defense 49 25.43 0.000302 Cholesterol metabolism 8 1.02 0.002260 Coenzyme and prosthetic group metabolism 10 2.25 0.003470 Transport 36 18.38 0.003600 Carbohydrate metabolism 20 7.68 0.003840 Steroid hormone-mediated signaling 6 0.56 0.004890 Amino acid metabolism 11 3.13 0.012000 Protein metabolism and modification 71 48.25 0.017400

Testis Detoxification 10 1.33 0.000195 * after Bonferroni correction Supplementary Table 11

Overrepresented biological processes among the nonadditively expressed transcripts between F1 hybrids (from a cross between a Mus m. musculus mother with a Mus m. castaneus father) and both parents in the testis.

Observed Expected number of number of genes in genes in category category P-value* Intracellular protein traffic 164 83.05 2.26x1014 Cell structure and motility 153 89.74 1.03 x108 Cell cycle 140 79.95 1.12 x108 Protein modification 173 103.42 1.33 x108 Spermatogenesis and motility 37 12.81 4.96 x106 Gametogenesis 51 21.48 5.29 x106 Nucleoside, nucleotide and nucleic acid metabolism 393 306.36 5.50 x106 Other intracellular protein traffic 21 5.25 2.47 x105 Protein metabolism and modification 385 303.81 2.52 x105 Pre-mRNA processing 56 26.41 4.41 x105 Lipid, fatty acid and steroid metabolism 112 69.93 4.66 x105 Protein targeting and localization 40 17.26 5.66 x105 Protein phosphorylation 100 59.67 1.63 x104 Cell structure 90 53.94 4.93 x104 mRNA splicing 43 19.57 5.68 x104 Lipid metabolism 31 12.33 7.71 x104 General vesicle transport 45 21.48 8.17 x104 Miscellaneous 30 13.44 2.00 x103 Mitosis 55 29.75 2.76 x103 Exocytosis 32 14.24 4.78 x103 Other intracellular signaling cascade 38 18.54 8.92 x103 Protein folding 32 14.8 9.65 x103 Intracellular signaling cascade 106 72.31 1.38 x102 Carbohydrate metabolism 73 48.37 1.56 x102 * after Bonferroni correction Table 12

Literature review on studies differentiating additivity and nonadditivity in gene expression. Shaded in grey indicate the studies that have less rigorously investigated the topic of additivity/nonadditivity.

Divergence 1-dye time vs. Inbreeding Tissue Add. vs non- Organism (generations) Reference Method 2-dye status complexity add. Studies consistent with additivity predominant Affymetrix Drosophila within species (HUGHES et al. 2006) microarray 1 no whole fly add>>non-add (STUPAR and Affymetrix Maize within species SPRINGER 2006) microarray 1 high whole plant add>>non-add all above (SWANSON-WAGNER ground Maize within species et al. 2006) cDNA microarray 2 high tissues add>>non-add M. musculus lab Affymetrix strains within species (CUI et al. 2006) microarray 1 high single tissue add>>non-add M. musculus wild Affymetrix strains 2 - 3 Mio* this study microarray 1 partially single tissue add>>non-add D. melanogaster vs. D. simulans 12.5 - 34 Mio# (RANZ et al. 2004) cDNA microarray 2 high head add>non-add D. simulans vs. (MICHALAK and NOOR Affymetrix D. mauritiana 2.5 - 8 Mio# 2003) microarray 1 high whole fly non-add small D. persimilis differential vs. (REILAND and NOOR display PCR (DD- D. pseudoobscura 2.5 -5 Mio# 2002) PCR) 0 high whole fly non-add small Studies consistent with nonadditivity predominant Agilent Drosophila within species (GIBSON et al. 2004) microarray 2 high whole fly non-add>add (VUYLSTEKE et al. Arabidopsis within species 2005) cDNA microarray 2 high first leaf pair add = non-add massive parallel (HEDGECOCK et al. signature partially 1-2 Mio Oyster within species 2007) sequencing 0 (f=0.375) larvae non-add>>add D. melanogaster vs D. simulans 12.5 - 34 Mio# (RANZ et al. 2004) cDNA microarray 2 high body non-add>add * assuming 2-3 generations/year (KARN et al. 2002)

# assuming 5-10 generations/year (HEY and KLIMAN 1993) Figure 1 Figure 2 Figure 3 A

musculus-domesticus musculus-castaneus

Testis (335) Liver (182) Testis (295) Liver (164)

295 15 145 263 9 130 8 8 17 14 15 17

107 97

Brain (146) Brain (137) Figure 3 B

musculus-domesticus musculus-castaneus

Brain 110 36 101

Liver 139 43 121

Testis 138 97 198 Figure 4