The identification of genes with large effects on sexual isolation and speciation is an important link between classic evolutionary genetics and molecular biology. Few genes that affect sexual isolation and speciation have been identified, perhaps because many traits influencing sexual isolation are complex behaviors. Cuticular hydrocarbons (CHs) of species of the Drosophila melanogaster group play a large role in sexual isolation by functioning as contact pheromones influencing mate recognition. Some of the genes that play key roles in determining species-specific CHs have been identified. We have performed separate quantitative trait locus (QTL) analyses of 7-tricosene (7-T) and 7,11-heptacosadiene (7,11-HD), the two major female CHs differing between D. simulans and D. sechellia. We find that 40% of the phenotypic variance in each CH is associated with two to four chromosomal regions. A region on the right arm of chromosome 3 contains QTL that affect both traits, but other QTL are in distinct chromosomal regions. Epistatic interactions were detected between two pairs of QTL for 7,11-HD such that if either were homozygous for the D. simulans allele, the fly was similar to D. simulans in phenotype, with a low level of 7,11-HD. We discuss the location of these regions with regard to candidate genes for CH production, including those for desaturases.

---

STUDIES of the genetics of speciation have, until very recently, been hindered by a paucity of candidate genes for detailed analysis (SWANSON and VACQUIER 2002; NOOR 2003). Reproductive isolation, responsible for speciation, is likely to involve complex, coevolved, polygenic traits (leading to a "type I" genetic architecture in the terminology of TEMPLETON 1981, i.e., numerous genes of small effect). Most empirical studies of reproductive isolation, especially of sexual isolation, have confirmed polygenic effects (HOLLOCHER et al. 1997; RITCHIE and PHILLIPS 1998; TING et al. 2001); however, a few large effect genes have been identified for both postmating (TING et al. 1998; BARBASH et al. 2003; PRESGRAVES et al. 2003) and premating isolation (WHEELER et al. 1991; GREENBERG et al. 2003). In Drosophila, most of these large effect genes have been identified by mutagenesis of D. melanogaster or by association analysis with geographic variation. Quantitative trait loci (QTL) studies, in contrast, can identify the minimum number and location of genes influencing complex traits (MACKAY 2001) and therefore identify the traits that are likely to be influenced by genes of large effect underlying naturally occurring variation. QTL studies may also indicate the potential candidate genes underlying genetic differentiation.

Cuticular hydrocarbons (CHs) affect mating behavior in many Drosophila species and are traits with the potential for a type II (large gene effect) genetic architecture (TEMPLETON 1981). CHs are long-chain hydrocarbons located on the cuticle surface that prevent desiccation and can function as contact pheromones (BLOMQUIST et al. 1987). During courtship, male Drosophila perform several different behaviors, including orienting toward a female, vibrating a wing (producing courtship "song"), and tapping her abdomen with his foreleg. In tapping the female, the male detects her CHs. Drosophila species are either monomorphic, with both males and females having the same CH composition (e.g., D. simulans and D. mauritiana, JALLON and DAVID 1987), or dimorphic with differing profiles (e.g., D. melanogaster and D. sechellia, COBB et al. 1989). CH profiles affect interspecific mating asymmetrically: males of monomorphic species will not court females of dimorphic species, whereas males of dimorphic species will court females of monomorphic species (COBB and JALLON 1990).

The major CH affecting mate recognition in D. simulans is 7-tricosene (7-T, PÉCHINÉ et al. 1985), which is produced by both males and females. Male D. sechellia have several tricosenes, including the same isomer, whereas females have very little 7-T and high levels of 7,11-heptacosadiene (7,11-HD, COBB et al. 1989). In crowding experiments, in which one female of a species was housed with females of the other species, the CH of the abundant females was transferred to the single female. In subsequent mating tests, the behavior of the males toward the females was changed such that D. simulans females housed with D. sechellia were courted less frequently than control D. simulans females by D. simulans males. In contrast, D. sechellia females housed with D. simulans were courted more vigorously than D. sechellia controls (COYNE et al. 1994).

Through study of the biochemistry of CH production in Drosophila, several genes that affect CH production have been identified. In insects, especially Dipterans, long-chain hydrocarbons with odd numbers of carbons are derived from even-numbered fatty acids by reduction and decarboxylation (BLOMQUIST et al. 1987), a pattern confirmed for the biosynthesis of the major CHs of D. melanogaster and D. simulans (JALLON 1984; CHAN YONG and JALLON 1986; PENNANEC'H et al. 1997). Early pathway steps may be common to both sexes and both species. The main substrates are myristate and palmitate, which are produced by a fatty acid synthetase (DE RENOBALES and BLOMQUIST 1984). From labeling experiments, two types of desaturases are implicated and molecular studies have discovered two closely linked genes that code 9 desaturases, but use different substrates: palmitate for desat1 (WICKER-THOMAS et al. 1997; DALLERAC et al. 2000) and myristate for desat2 (DALLERAC et al. 2000). The desat1 gene is involved in the first desaturation of these saturated fatty acids leading to hydrocarbons bearing a 7 double bond common to males and females of both species (LABEUR et al. 2002). The desat1 gene affects the level of 7,11-HD and the desat2 gene affects that of 5,9-HD in D. melanogaster females (WICKER-THOMAS et al. 1997; COYNE et al. 1999; DALLERAC et al. 2000; TAKAHASHI et al. 2001; FANG et al. 2002; GREENBERG et al. 2003; MARCILLAC et al. 2005). The desat2 gene, functional in African D. melanogaster populations, is not functional in Cosmopolitan strains because of a small deletion in the promoter region and results in females that are rich in 5,9-HD (TAKAHASHI et al. 2001; WICKER-THOMAS and JALLON 2001). The desat1 gene differs between D. melanogaster types in only three amino acid replacements and its function is not greatly altered, leading to the major production of 7-monenes in both sexes and to 7,11-dienes in non-African females (DALLERAC et al. 2000). Environmental selection upon the CH components may have driven the change at desat2 (GREENBERG et al. 2003), but whether this has also influenced assortative mating within D. melanogaster (including the premating isolation seen between Zimbabwe and Cosmopolitan strains) is uncertain (COYNE et al. 1999; GREENBERG et al. 2003; RITCHIE and NOOR 2004).

The diene second double bond of D. melanogaster females requires a second desaturation step, which is probably catalyzed by a different desaturase (WICKER-THOMAS and JALLON 2001). Flies carrying various deficiencies in the third chromosome region 67E–69B have variation in unsaturated hydrocarbon levels in both sexes, although the difference is much larger in females than in males (WICKER-THOMAS and JALLON 2000). One gene in this region, Enhancer of zeste [E(z)], decreases the level of female 7,11-HD and has a correlated increase in 7-heptacosene, depending on the strength of the allele (WICKER-THOMAS and JALLON 2000). The same type of phenotype has been observed in mutations in two temperature-sensitive endocrine genes, ecdysoneless (ecd, WICKER and JALLON 1995) and Dopa decarboxylase (Ddc, MARICAN et al. 2004) when females were bred at restrictive temperature.

Other genes have also been implicated in CH variation and possibly mate choice. Two second chromosome genes, seven pentacosene (sept) and small monoene quantities (smoq), affect the balance between 7-T and 7-pentacosene (7-P) in D. melanogaster males (FERVEUR and JALLON 1996). The nerd locus affects the level of 7-T production in D. melanogaster (FERVEUR and JALLON 1993b) and has been mapped to the third chromosome. Although sept and smoq affect the relative production of 7-T and 7-P in D. melanogaster, in D. simulans, a single locus on the second chromosome, Ngbo, has a dose-dependent effect on production of 7-P. Ngbo has been proposed to be a structural gene that changes enzyme specificity (FERVEUR 1991). Also in D. simulans, the sex-linked kété locus affects overall levels of 7-T, 7-P, and all linear hydrocarbons (FERVEUR and JALLON 1993a).

Genes involved in somatic sexual differentiation may also affect CH composition. The gene doublesex (dsx) affects the production of female-specific dienes in D. melanogaster (JALLON et al. 1988). The female-specific transcript DSX F acts as a dominant regulator of female dienes (WATERBURY et al. 1999). The intersex (ix) locus affects the ratio of production of the female-specific 7,11-HD and 7,11-nonacosadiene (7,11-ND) to the production of 7-T and 5-T, which are more prominent in males (WATERBURY et al. 1999). Differential expression of transformer (tra) affects the production of sex-specific pheromones (SAVARIT et al. 1999) as does transformer-2 (tra2, JALLON et al. 1986) and Sex-lethal (Sxl, TOMPKINS and MCROBERT 1989; TOMPKINS and MCROBERT 1995). The homeotic genes Antennapedia (Antp) and Ultrabithorax (Ubx) also influence the relative levels of monoenes vs. dienes (WICKER-THOMAS and JALLON 2001).

With the completion of the D. melanogaster genome sequence, additional genes are identifiable through domain similarity with genes of known function. Two such genes are Fad2 and infertile crescent (ifc), both of which have stearoyl CoA-desaturase activity (FLYBASE CONSORTIUM 2003). Additional unknown genes with this function include CG8630, CG9743, and CG15531. Other activities include a SUR2-type hydroxylase/desaturase catalytic domain (CG1998, CG11162, and CG30502), acyl-CoA 11-desaturase activity (CG9747), and a fatty acid elongase (CG2781).

In this study, we use quantitative trait loci (QTL) mapping to refine the number and locations of loci affecting 7-T and 7,11-HD, the most important CHs differing between D. simulans and D. sechellia, and compare the position of QTL with that of the candidate genes. Previously, the location of genes involved in the difference in CH composition between D. simulans and D. sechellia has been studied in crossing experiments with limited markers. Using a D. sechellia strain with only one marker per chromosome, COYNE et al. (1994) studied the ratio of 7-T to 7,11-HD in the female progeny of hybrid females backcrossed to D. sechellia. The third chromosome had the largest effect on the trait. The X chromosome had an effect only on the amount of 7-P, a polymorphic CH found only on D. sechellia, which also induced male courtship by D. sechellia males (COBB and FERVEUR 1996; COYNE et al. 1994).

Strains and crosses:

One strain each of the two species was used in this study. A D. sechellia strain was kindly provided by Jean R. David. The strain was inbred for 18 generations of brother-sister mating to produce the strain DavidA4 used in the subsequent crosses. The D. simulans line (f²;nt, pm;st,e), kindly provided by Jerry Coyne, had five morphological markers, one per chromosome arm (Table 1).These are the same strains used in a QTL analysis of courtship song differences between these species (GLEASON and RITCHIE 2004).

View this table: In this window In a new window   TABLE 1 Markers used and their resulting map order

 

View larger version (23K): In this window In a new window Download PPT slide   FIGURE 1.— QTL mapping of cuticular hydrocarbons. The amount of both 7-tricosene (7-T, shaded) and 7,11-heptacosadiene (7,11-HD, solid) in D. simulans-D. sechellia backcross hybrids was mapped using composite interval mapping. The positions of the markers, listed in Table 1, are denoted by circles along the x-axis. Triangles at the top of the graph indicate the positions of each QTL and are numbered for reference to the text. The significance level at P = 0.05 was determined by 1000 permutations of the data set and is shown for 7-T at a likelihood ratio of 14.34 (the level for 7,11-HD is 13.95). The locations of candidate genes are indicated by vertical lines along the x-axis, with the exception of the one candidate gene, tra, that was also a marker and is represented by a shaded circle. Abbreviations for candidate genes are given in Table 6.

 
All fly culturing was at 25° using standard techniques and a 12 hr light:12 hr dark cycle. Female D. simulans were crossed to male D. sechellia and the female progeny were backcrossed to D. simulans males. For each cross, one female was paired with a single male in a vial (95 x 16.5 mm) for 7 days. A total of 44 crosses and 78 backcrosses were used to generate 487 females that were subsequently analyzed for cuticular hydrocarbons. Using the five morphological markers, we attempted to sample equally among the 32 backcross phenotypes to obtain representatives of all recombinants.

Cuticular hydrocarbon extraction and analysis:

Flies were collected within 8 hr of eclosion and phenotyped for the five morphological markers. Females were housed individually in small vials (95 x 16.5 mm). Five days post-eclosion, the female fly was placed for 5 min in 70 µl hexane containing 800 ng hexacosane as a control for recovery. The fly was then removed and DNA was extracted using the single fly extraction method and a total of 100 µl squishing buffer (GLOOR and ENGELS 1992). The hexane containing the cuticular hydrocarbons was dried by evaporation and stored at –20° until gas chromatography.

Before gas chromatography, the sample was resuspended in 40 µl hexane. Five microliters was injected into a Perkin-Elmer (Norwalk, CT) gas chromatograph autosystem with split injection, a flame ionization detector, and a 25QC2/BP1 0.1 column (Scientific Glass Engineering, 25 mm x 0.22 mm x 0.1 µm) programmed to run from 180° to 270° with a 3°/min gradient.

The quantities of two cuticular hydrocarbons, 7-T and 7,11 HD, were calculated from the resulting chromatographs as the area under each peak. These were standardized by dividing by the quantity of hexacosane recovered, to account for differences in extraction and gas chromatography analysis. Eight data outliers were removed and the resulting values were natural log transformed to remove a right skew. Remaining residuals were normally distributed. All subsequent analyses were done with the transformed variables. Final sample sizes were 484 females for 7-T and 482 females for 7,11-HD.

Marker scoring:

Thirty-nine molecular markers were scored for each individual (N = 487, Table 1). These markers were all polymerase chain reaction (PCR) amplified and showed different size fragments for D. sechellia and D. simulans on 2% agarose, 4% Metaphor agarose (Cambex), or 8% acrylamide gels. Size differences were caused by natural variation in sequence length (indels or microsatellites) or differences in restriction enzyme sites (Table 1). Hybrids were easily distinguished from homozygotes.

Genetic mapping and QTL analysis:

With the morphological markers (see Strains and crosses), a total of 44 markers were scored. These markers were mapped using MAPMAKER (LANDER et al. 1987). The map obtained was subsequently used in QTL analyses using QTL Cartographer v1.16c (BASTEN et al. 1997) to map QTL. Epistatic interactions between QTL were detected by a generalized linear model in SPSS 13.0 (SPSS) with Bonferroni correction for multiple tests, using allelic data at the nearest marker to the QTL.

Calculation of effects in terms of parental difference:

All analyzed data were natural log tranformed. The effect, as calculated, is the difference between the mean of the logs for the genotype (MLG) and the mean of the logs for the traits of all backcross individuals (MLT). This difference is equivalent to the log of ratios of the geometric mean for the genotypes (GMG) and the grand geometric mean (GMT); that is, MLG – MLT = log(GMG/GMT). Therefore, if exp(effect)–1 is multiplied by GMG, the result is GMT. Subtracting GMT from GMG yields the effect in the original units. The percentage of the parental difference was then calculated by dividing this effect by the difference in the mean value between the parental strains.

Marker mapping:

The markers were chosen to have an average spacing of approximately 7 cM on the basis of the D. melanogaster map. The average spacing realized between markers was 23.2 cM, because of segregation distortion and experimenter selection for recombinants. Most markers mapped in the same linear order as in D. melanogaster with the exception of two marker pairs on each end of the second chromosome (ex and nt; pm and twi). On the right arm of the third chromosome there is an inversion in D. simulans and D. sechellia relative to D. melanogaster and five markers (Mtn, pros, gl, nos, e) show this inversion with respect to D. melanogaster.

Trait values:

As expected, D. simulans females had more 7-T than D. sechellia and less 7,11-HD (Table 2). Values for hybrid and backcross females were intermediate, but hybrid females were more like D. sechellia for 7-T and more like D. simulans for 7,11-HD [similar to other crosses between D. mauritiana and D. sechellia (COYNE and CHARLESWORTH 1997)], but different from crosses between D. simulans and D. melanogaster (COYNE 1996). Backcross females were more like D. simulans than D. sechellia for both CHs. CH components are potentially interrelated in that they are products of a common production pathway. The correlation between 7-T and 7,11-HD across recombinant females here was only –0.2, indicating the potential for different genes to independently influence the production of the compounds. We therefore analyzed QTL for the two traits separately (previous studies have analyzed the ratio of the two compounds).

View this table: In this window In a new window   TABLE 2 Mean ± standard deviation (in nanograms) and sample sizes (n) for both cuticular hydrocarbons

 

QTL analysis:

The marker on the fourth chromosome (ey) was not significantly associated with either trait and thus results for this chromosome, which composes about 1% of the genome, are not shown. Composite interval mapping (CIM) was performed for the rest of the genome. CIM (JANSEN and STAM 1994; ZENG 1994) combines interval mapping (LANDER et al. 1987) with multiple regression. Each interval flanked by adjacent markers is tested for the presence of a QTL affecting the trait while statistically accounting for the effects of additional segregating QTL outside the interval. Significance levels of P = 0.05 for both CH were calculated from 1000 permutations of the trait data among marker classes (CHURCHILL and DOERGE 1994) and corresponds to a likelihood ratio of 14.33 for 7-T and 13.95 for 7,11-HD.

Parameters potentially affecting the detection of QTL using CIM are the size of the window around the tested interval, within which linked markers are excluded from multiple regression, and the number of background markers used. We tested a range of window sizes from 2.5 to 20 cM and this parameter did not influence the result. Figure 1 depicts the results using a backcross design, the Kosambi map function, all background markers, a walking speed of 2 cM, and a window size of 5 cM. Forward/backward stepwise regression resulted in nine significant markers for 7-T and eight significant markers for 7,11-HD that could be used in CIM. Results varied slightly with the addition of markers. For 7-T, QTL 6 narrowed with the addition of the fifth marker, Mtn, which bounds one edge of the QTL. Otherwise, all QTL were stable after the addition of the second marker for this trait. For 7,11-HD, adding the third marker (which is on the right boundary of QTL2) reduced QTL1 from approximately equal to QTL2 to a value of 26. Adding the sixth marker (which is on the right boundary of QTL 3) reduced it further to the presented value. For QTL3, adding the second marker (Mtn, in the middle of the QTL) narrowed this QTL to where it became stable. For these reasons, using all the markers for the final results gives a conservative estimate of the height of the QTL.

Seven QTL were found for 7-T and three for 7,11-HD (Figure 1, Table 3). For 7-T, QTL 1 and 2 are located on the X chromosome (Figure 1), whereas the other five are on the third chromosome. Six of the QTL are in pairs with one contributing more to the phenotypic variance (V_p) than the other (QTL1 and 2, QTL6 and 7) or both having a small contribution (QTL3 and 4), and QTL5 is independent. Therefore, at least four distinct genomic regions are associated with the trait, with the greatest effects being on the right-hand sides of the first and third chromosomes. For 7-T, the proportion of V_p explained by each QTL varies from 2.29 to 17.68% (Table 3). Together, the QTL explain 44.10% of V_p. The effects for the QTL are all in the same direction, consistent with the trait being expressed in one species and not in the other. The effects, expressed as a percentage of the difference in means between the two parental species, range from 53 to 68% (Table 3).

View this table: In this window In a new window   TABLE 3 QTL locations and effects

 

View this table: In this window In a new window   TABLE 4 Significance of epistatic interaction of 7,11-HD

 
For 7,11-HD, all of the QTL are on the third chromosome (Figure 1). Two QTL are paired, with one contributing more than the other (QTL1 and 2), and the other is independent, so that CIM identifies two regions of large effect. The QTL vary from 2.77 to 18.20% of V_p and together account for a total of 37.97% of V_p (Table 3). Again, all of the effects are in the same direction, which is opposite of that for the 7-T QTL. The effects are approximately the same when expressed as a proportion of the parental difference, ranging from 6.4 to 6.7%. This is a small proportion, heavily influenced by our crossing design (backcrosses to D. simulans). Hence, there was relatively little V_p, centered around the D. simulans parental value, so only a fraction of the parental difference could be explained here.

Epistatic interactions were tested using a generalized linear model in SPSS. No interactions between pairs of QTL for 7-T were significant after Bonferroni correction. For 7,11-HD, there were significant interactions between QTL3 and the other two QTL (Table 4). The interaction effect is such that both QTL regions need to be in a hybrid state to produce a significant amount of 7,11-HD (Table 5).

View this table: In this window In a new window   TABLE 5 Trait values for epistatic interactions of 7,11-HDa

 
Twenty-seven candidate genes for cuticular hydrocarbon biosynthesis were identified from the literature and from the available sequence data. The approximate locations of the candidate genes were estimated by their genetic location in D. melanogaster. Candidate genes do not fall directly with the QTL (Table 6, Figure 1), although several are on the edges of QTL including CG2781, dsx, desat1, and desat2 for 7-T. E(z) and Fad2 are on the edge of a 7,11-HD QTL.

View this table: In this window In a new window   TABLE 6 Cuticular hydrocarbon biosynthesis candidate genes coincidence with cuticular hydrocarbon QTL

 

QTL studies of interspecific differences in adaptive quantitative traits have found QTL with both large (LAURIE et al. 1991; MACDONALD and GOLDSTEIN 1999) and minor effects (e.g., FISHMAN et al. 2002). Many interspecific trait differences have been shown to be polygenic (e.g., KIM and RIESEBERG 1999; ZENG et al. 2000). The time since species divergence has been hypothesized to be positively correlated with the number of QTL found (KIM and RIESEBERG 1999), although another possibility is that some traits are more prone to major gene effects than others. If a threshold of >25% of the phenotypic variance explained is used to designate a major QTL (BRADSHAW et al. 1995, 1998), then neither trait studied here is influenced by major QTL. The majority of the phenotypic variation is not explained, indicating that, despite the large sample size, the study lacked sufficient resolution to detect additional small-effect QTL that also influence this trait difference. However, a high proportion of the parental difference is explained for 7-T. Given the values for the hybrids (Table 2), there are clearly dominance effects of alleles for both traits. For 7-T, D. sechellia alleles appear to be dominant, whereas for 7,11-HD, D. simulans alleles appear to be dominant. Using a backcross design, we cannot distinguish dominance from other genetic effects; because of the sterility of F₁ males, backcross analysis is the only crossing scheme possible with these species. Furthermore, we cannot detect all of the possible effects of D. sechellia alleles because our backcross was to D. simulans.

The effects of the QTL are consistent with the difference between the species: all effects are positive for 7-T and negative for 7,11-HD. This pattern is consistent with directional selection operating on these two species. The question remains open whether this might be from sexual selection, by male preferences for the female trait, or from natural selection, such as that found for different dienes in D. melanogaster females (GREENBERG et al. 2003). The epistatic effect means that both members of interacting loci need to have an allele from D. sechellia to produce increased levels of 7,11-HD. This implies that each locus codes for an essential step in the production of 7,11-HD, which is absent in D. simulans. The combination of D. sechellia alleles at both loci produces a substantial change away from the D. simulans-like phenotype.

The large contribution of QTL on the third chromosome to both traits is probably what was detected by COYNE et al. (1994) in their study with one marker per chromosome. In that study, the ratio of 7-T to 7,11-HD was analyzed rather than each compound separately; thus the effect of the X chromosome, contributing only to 7-T production, was not detected. In a study using introgression lines (CIVETTA and CANTOR 2003), a QTL for 7-T quantity in males of these species was found at approximately the same postion as our 7-T QTL 3 and 4. In contrast to our results, they did not find any QTL for this trait in females, probably because the study had low resolution for D. sechellia (CIVETTA and CANTOR 2003).

A QTL study for a courtship song difference between D. simulans and D. sechellia (GLEASON and RITCHIE 2004) found that most of the QTL do not coincide with previously identified candidate genes. The results here are similar. In Figure 1, the 7-T QTL 6 appears to be adjacent to two fatty acid desaturases, desat1 and desat2. Because alleles of desat2 affect the ratio of 7,11-HD, in D. melanogaster females (GREENBERG et al. 2003), this gene was thought to be a strong candidate gene for the species difference. Using a 2 LOD (9.22 likelihood ratio) confidence interval for the width of 7,11-HD QTL 3, these genes are not even adjacent to that QTL.

The desaturase family of enzymes do not colocalize. The complete D. melanogaster genome has revealed five more fatty acid desaturase genes, in addition to desat1 and desat2, on chromosome 3 and one on chromosome 2 (JALLON and WICKER-THOMAS 2003). Markedly decreased levels of dienes were found in D. melanogaster females carrying deletions in the region 67E–69B [corresponding here to the region surrounding E(z) and Fad2, Figure 1]. Molecular work in progress in that region has characterized Fad2 as an additional desaturase gene, which, expressed only in females of the dimorphic species, acts on unsaturated fatty acids and leads to dienes (T. CHERTEMPS and C. WICKER-THOMAS, unpublished results).

Another group of enzymes (elongases) is involved in the elongation of the hydrocarbon chains and one of these (CG2781) is found on the third chromosome and on the edge of 7-T QTL 5. The same is true for doublesex. Although these two genes are unlikely to contribute to this QTL on the basis of these data, to distinguish the contribution of linked candidate genes to wide QTL peaks is notoriously difficult, and further resolution will require high-level recombination mapping and assessing allelic variation at the candidate loci themselves.

While there is partial overlap in some of the QTL identified for the two CH components studied here, at least three that affect 7-T are not implicated in the production of 7,11-HD. Thus, some steps in the biosynthesis of these compounds are affected by different genes. 7-T and 7,11-HD differ in chain length (23 vs. 27 carbons, respectively) and saturation (one vs. two double bonds, respectively). The genetic effects could therefore be in chain elongation, desaturation, or in overall levels of production. Only by further refined mapping to the genes underlying these traits will we be able to determine the contribution of each gene to hydrocarbon production. Despite the small effects relative to the phenotypic variance, the likelihood ratios of the QTL presented here are very large and indicate that it should be possible to map these traits more finely and identify the genes contributing to this important trait difference.

For technical help with scoring markers, we thank Tanya Hamill and Melanie Edgar. Carrie Adamson assisted with the crosses, cuticular hydrocarbon extractions, and marker scoring. Rosemary Bevan assisted with Drosophila culturing and Terry Gleason assisted with statistical advice. The manuscript was improved thanks to the comments of two anonymous reviewers. This work was supported by a grant (GR3/10786) from the Natural Environment Research Council (UK) to M.G.R.

BARBASH, D. A., D. F. SIINO, A. M. TARONE and J. ROOTE, 2003 A rapidly evolving MYB-related protein causes species isolation in Drosophila. Proc. Natl. Acad. Sci. USA 100: 5302–5307.[Abstract/Free Full Text]

BASTEN, C. J., B. S. WEIR and Z-B. ZENG, 1997 QTL Cartographer: A Reference Manual and Tutorial for QTL Mapping. Department of Statistics, North Carolina State University, Raleigh, NC.

BLOMQUIST, G. J., J. W. DILLWORTH and T. S. ADAMS, 1987 Biosynthesis and endochrine regulation of sex pheomone production in Diptera, pp. 217–250 in Pheromone Biochemistry, edited by G. D. PRESTWITCH and G. J. BLOMQUIST. Academic Press, London.

BRADSHAW, H. D. J., K. G. OTTO, B. E. FREWEN, J. K. MCKAY and D. W. SCHEMSKE, 1998 Quantitative trait loci affecting differences in floral morphology between two species of monkeyflower (Mimulus). Genetics 149: 367–382.[Abstract/Free Full Text]

BRADSHAW, H. D. J., S. M. WILBERT, K. G. OTTO and D. W. SCHEMSKE, 1995 Genetic mapping of floral traits associated with reproductive isolation in monkeyflowers (Mimulus). Nature 376: 762–765.[CrossRef]

CHAN YONG, T. P., and J.-M. JALLON, 1986 Synthèse de novo d'hydrocarbures potentiellement aphrodisiaques ches les Drosophile. C. r: hebd. Séanc. Acad. Sci. 303: 197–202.

CHURCHILL, G. A., and R. W. DOERGE, 1994 Empirical threshold values for quantitative trait mapping. Genetics 138: 963–971.[Abstract]

CIVETTA, A., and E. J. F. CANTOR, 2003 The genetics of mating recognition between Drosophila simulans and D. sechellia. Genet. Res. 82: 117–126.[CrossRef][Medline]

COBB, M., B. BURNET, R. BLIZARD and J.-M. JALLON, 1989 Courtship in Drosophila sechellia: its structure, functional aspects, and relationship to those of other members of the Drosophila melanogaster species subgroup. J. Insect Behav. 2: 63–89.

COBB, M., and J.-F. FERVEUR, 1996 Evolution and genetic control of mate recognition and stimulation in Drosophila. Behav. Proc. 35: 35–54.[CrossRef]

COBB, M., and J.-M. JALLON, 1990 Pheromones, mate recognition and courtship stimulation in the Drosophila melanogaster species subgroup. Anim. Behav. 39: 1058–1067.[CrossRef]

COLSON, I., and D. B. GOLDSTEIN, 1999 Evidence for complex mutations at microsatellite loci in Drosophila. Genetics 152: 617–627.[Abstract/Free Full Text]

COLSON, I., S. J. MACDONALD and D. B. GOLDSTEIN, 1999 Microsatellite markers for interspecific mapping of Drosophila simulans and D. sechellia. Mol. Ecol. 8: 1951–1955.[CrossRef][Medline]

COYNE, J. A., 1996 Genetics of a difference in male cuticular hydrocarbons between two sibling species, Drosophila simulans and D. sechellia. Genetics 143: 1689–1698.[Abstract]

COYNE, J. A., and B. CHARLESWORTH, 1997 Genetics of a pheromonal difference affecting sexual isolation between Drosophila mauritiana and D. sechellia. Genetics 145: 1015–1030.[Abstract]

COYNE, J. A., A. P. CRITTENDEN and K. MAH, 1994 Genetics of a pheromonal difference contributing to reproductive isolation in Drosophila. Science 265: 1461–1464.[Abstract/Free Full Text]

COYNE, J. A., C. WICKER-THOMAS and J.-M. JALLON, 1999 A gene responsible for a cuticular hydrocarbon polymorphism in Drosophila melanogaster. Genet. Res. 73: 189–203.[CrossRef][Medline]

DALLERAC, R., C. LABEUR, J.-M. JALLON, D. C. KNIPPLE and W. L. ROELOFS, 2000 A 9 desaturase gene with a different substrate specificity is responsible for the cuticular diene hydrocarbon polymorphism in Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 97: 9449–9454.[Abstract/Free Full Text]

DE RENOBALES, M., and G. J. BLOMQUIST, 1984 Biosynthesis of medium chain fatty acids in Drosophila melanogaster. Arch. Biochem. Biophys. 228: 407–414.[CrossRef][Medline]

FANG, S., A. TAKAHASHI and C.-I WU, 2002 A mutation in the promoter of desaturase 2 is correlated with sexual isolation between Drosophila behavioral races. Genetics 162: 781–784.[Abstract/Free Full Text]

FERVEUR, J.-F., 1991 Genetic control of pheromones in Drosophila simulans. I. Ngbo, a locus on the second chromosome. Genetics 128: 293–301.[Abstract]

FERVEUR, J.-F., and J.-M. JALLON, 1993a Genetic control of pheromones in Drosophila simulans. II. kété, a locus on the X chromosome. Genetics 133: 561–567.[Abstract]

FERVEUR, J.-F., and J.-M. JALLON, 1993b Nerd, a locus on Chromosome III, affects male reproductive behavior in Drosophila melanogaster. Naturwissenshaften 80: 474–475.[CrossRef][Medline]

FERVEUR, J.-F., and J.-M. JALLON, 1996 Genetic control of male cuticular hydrocarbons in Drosophila melanogaster. Genet. Res. 67: 211–218.[Medline]

FISHMAN, L., A. J. KELLY and J. H. WILLIS, 2002 Minor quantitative trait loci underlie floral taits associated with mating system divergence in Mimulus. Evolution 56: 2138–2155.[CrossRef][Medline]

FLYBASE CONSORTIUM, 2003 The FlyBase database of the Drosophila genome projects and community literature. Nucleic Acids. Res. 31: 172–175.[Abstract/Free Full Text]

GLEASON, J. M., and M. G. RITCHIE, 2004 Do quantitative trait loci for a courtship song difference between Drosophila simulans and D. sechellia overlap with candidate genes and intraspecific QTLs? Genetics 166: 1303–1311.[Abstract/Free Full Text]

GLOOR, G., and W. ENGELS, 1992 Single-fly DNA preps for PCR. Dros. Inf. Ser. 71: 148–149.

GOLDSTEIN, D. B., and A. G. CLARK, 1995 Microsatellite variation in North American populations of Drosophila melanogaster. Nucleic Acids. Res. 23: 3882–3886.[Abstract/Free Full Text]

GREENBERG, A. J., J. R. MORAN, J. A. COYNE and C.-I WU, 2003 Ecological adaptation during incipient speciation revealed by precise gene replacement. Science 302: 1754–1757.[Abstract/Free Full Text]

HOLLOCHER, H., C. T. TING, M. L. WU and C.-I WU, 1997 Incipient speciation by sexual isolation in Drosophila melanogaster: variation in mating preference and correlation between sexes. Evolution 51: 1175–1181.[CrossRef]

JALLON, J.-M., 1984 A few chemical words exchanged by Drosophila during courtship and mating. Behav. Genet. 14: 441–478.[CrossRef][Medline]

JALLON, J.-M., C. ANTONY, T. P. CHAN YONG and S. MANIAR, 1986 Genetic factors controlling the production of aphrodisiac substance in Drosophila. Adv. Invert. Reprod. 4: 445–452.

JALLON, J.-M., and J. R. DAVID, 1987 Variations in cuticular hydrocarbons among the eight species of the Drosophila melanogaster subgroup. Evolution 41: 294–302.[CrossRef]

JALLON, J.-M., G. LAUGE, L. ORSSAUD and C. ANTONY, 1988 Female pheromones in Drosophila melanogaster Genetical Researchare controlled by the doublesex locus. Genet. Res. 51: 17–22.

JALLON, J.-M., and C. WICKER-THOMAS, 2003 Genetic studies on pheromone production in Drosophila, pp. 253–281 in Insect Pheromone Biochemistry and Molecular Biology: The Biosynthesis and Detection of Pheromones and Plant Volatiles, edited by G. BLOMQUIST and R. VOGT. Elsevier/Academic Press, Boston.

JANSEN, R. C., and P. STAM, 1994 High resolution of quantitative traits into multiple loci via interval mapping. Genetics 136: 1447–1455.[Abstract]

KIM, S.-C., and L. H. RIESEBERG, 1999 Genetic architecture of species differences in annual sunflowers: implications for adaptive trait introgression. Genetics 153: 965–977.[Abstract/Free Full Text]

LABEUR, C., R. DALLERAC and C. WICKER-THOMAS, 2002 Involvement of desat1 gene in the control of Drosophila melanogaster pheromone biosynthesis. Genetica 114: 269–274.[CrossRef][Medline]

LANDER, E. S., P. GREEN, J. ABRAHAMSON, A. BARLOW, M. J. DALY et al., 1987 MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics 1: 174–181.[CrossRef][Medline]

LAURIE, C. C., J. T. BRIDGHAM and M. CHOUDHARY, 1991 Associations between DNA sequence variation and variation in expression of the Adh gene in natural population of Drosophila melanogaster. Genetics 129: 489–499.[Abstract]

LIU, J., J. M. MERCER, L. F. STAM, G. C. GIBSON, Z-B. ZENG et al., 1996 Genetic analysis of a morphological shape difference in the male genitalia of Drosophila simulans and D. mauritiana. Genetics 142: 1129–1145.[Abstract]

MACDONALD, S. J., and D. B. GOLDSTEIN, 1999 A quantitative genetic analysis of male sexual traits distinguishing the sibling species Drosophila simulans and D. sechellia. Genetics 153: 1683–1699.[Abstract/Free Full Text]

MACKAY, T. F. C., 2001 Quantitative trait loci in Drosophila. Nat. Rev. Genet. 2: 11–20.[Medline]

MARCILLAC, F., Y. GROSJEAN and J. F. FERVEUR, 2005 A single gene mutation alters production and discrimination of Drosophila sex pheromones. Proc. R. Soc. Lond. Ser. B. 272: 303–309.[Medline]

MARICAN, C., L. DUPORTETS, S. BIRMAN and J.-M. JALLON, 2004 Female-specific regulation of cuticular hydrocarbon biosynthesis by dopamine in Drosophila melanogaster. Insect Biochem. Mol. Biol. 34: 823–830.[Medline]

MICHALAKIS, Y., and M. VEUILLE, 1996 Length variation of CAG/CAA trinucleotide repeats in natural populations of Drosophila melanogaster and its relation to the recombination rate. Genetics 143: 1713–1725.[Abstract]

NOOR, M. A. F., 2003 Evolutionary biology: Genes to make new species. Nature 423: 699–700.[CrossRef][Medline]

PENNANEC'H, M., L. BRICARD, G. KUNESCH and J.-M. JALLON, 1997 Incorporation of fatty acids into cuticular hydrocarbons of male and female Drosophila melanogaster. J. Insect Phys. 43: 1111–1116.[CrossRef][Medline]

PRESGRAVES, D. C., L. BALAGOPALAN, S. M. ABMAYR and H. A. ORR, 2003 Adaptive evolution drives divergence of a hybrid inviability gene between two species of Drosophila. Nature 423: 715–719.[CrossRef][Medline]

PÉCHINÉ, J. M., F. PEREZ, C. ANTONY and J.-M. JALLON, 1985 A further characterization of Drosophila cuticular monoenes using a mass spectrometry method to localize double bonds in complex mixtures. Anal. Biochem. 145: 177–182.[CrossRef][Medline]

RITCHIE, M. G., and M. A. F. NOOR, 2004 Gene replacement and the genetics of speciation. Heredity 93: 1–2.[CrossRef][Medline]

RITCHIE, M. G., and S. D. F. PHILLIPS, 1998 The genetics of sexual isolation, pp. 291–308 in Endless Forms: Species and Speciation, edited by D. J. HOWARD and S. H. BERLOCHER. Oxford University Press, New York.

SAVARIT, F., G. SUREAU, M. COBB and J.-F. FERVEUR, 1999 Genetic elimination of known pheromones reveals the fundamental chemical bases of mating and isolation in Drosophila. Proc. Natl. Acad. Sci. USA 96: 9015–9020.[Abstract/Free Full Text]

SCHUG, M. D., T. F. C. MACKAY and C. F. AQUADRO, 1997 Low mutation rates of microsatellite loci in Drosophila melanogaster. Nat. Genet. 15: 99–102.[CrossRef][Medline]

SCHUG, M. D., K. A. WETTERSTRAND, M. S. GAUDETTE, R. H. LIM, C. M. HUTTER et al., 1998 The distribution and frequency of microsatellite loci in Drosophila melanogaster. Mol. Ecol. 7: 57–70.[CrossRef][Medline]

SWANSON, W. J., and V. D. VACQUIER, 2002 The rapid evolution of reproductive proteins. Nat. Rev. Genet. 3: 137–144.[Medline]

TAKAHASHI, A., S. TSAUR, J. A. COYNE and C.-I. WU, 2001 The nucleotide changes governing cuticular hydrocarbon variation and their evolution in Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 98: 3920–3925.[Abstract/Free Full Text]

TEMPLETON, A. R., 1981 Mechanisms of speciation—a population genetic approach. Ann. Rev. Ecol. and Syst. 12: 23–48.

TING, C., A. TAKAHASHI and C. WU, 2001 Incipient speciation by sexual isolation in Drosophila: concurrent evolution at multiple loci. Proc. Natl. Acad. Sci. USA 98: 6709–6713.[Abstract/Free Full Text]

TING, C. T., S. C. TSAUR, M. L. WU and C.-I WU, 1998 A rapidly evolving homeobox at the site of a hybrid sterility gene. Science 282: 1501–1504.[Abstract/Free Full Text]

TOMPKINS, L., and S. P. MCROBERT, 1989 Regulation of behavioral and pheromonal aspects of sex determination in Drosophila melanogaster by the Sex-lethal gene. Genetics 123: 535–541.[Abstract/Free Full Text]

TOMPKINS, L., and S. P. MCROBERT, 1995 Behavioral and pheromonal phenotypes associated with expression of loss-of-function mutations in the Sex-lethal gene of Drosophila melanogaster. J. Neurogenet. 9: 219–226.[Medline]

WATERBURY, J. A., L. L. JACKSON and P. SCHEDL, 1999 Analysis of the doublesex female protein in Drosophila melanogaster: role in sexual differentiation and behavior and dependence on intersex. Genetics 152: 1653–1667.[Abstract/Free Full Text]

WHEELER, D. A., C. P. KYRIACOU, M. L. GREENACRE, Q. YU, J. E. RUTILA et al., 1991 Molecular transfer of a species-specific behavior from Drosophila simulans to Drosophila melanogaster. Science 251: 1082–1085.[Abstract/Free Full Text]

WICKER, C., and J.-M. JALLON, 1995 Influence of ovary and ecdysteroids on pheromone biosynthesis in Drosophila melanogaster (Diptera: Drosphilidae). Europ. J. Ent. 92: 197–202.

WICKER-THOMAS, C., and J.-M. JALLON, 2000 Role of Enhancer of zeste on the production of Drosophila melanogaster pheromonal hydrocarbons. Naturwissenshaften 87: 76–79.[CrossRef][Medline]

WICKER-THOMAS, C., and J.-M. JALLON, 2001 Control of female pheromones in Drosophila melanogaster by homeotic gene. Genet. Res. 78: 235–242.[Medline]

WICKER-THOMAS, C., C. HENRIET and R. DALLERAC, 1997 Partial characterization of a fatty acid desaturase gene in Drosophila melanogaster. Insect Biochem. Mol. Biol. 27: 963–972.[CrossRef][Medline]

ZENG, Z-B., 1994 Precision mapping of quantitative trait loci. Genetics 136: 1457–1468.[Abstract]

ZENG, Z-B., J. LIU, L. F. STAM, C.-H. KAO, J. M. MERCER et al., 2000 Genetic architecture of a morphological shape difference between two Drosophila species. Genetics 154: 299–310.[Abstract/Free Full Text]

Communicating editor: D. RAND

This article has been cited by other articles:

				B. Foley, S. F. Chenoweth, S. V. Nuzhdin, and M. W. Blows Natural Genetic Variation in Cuticular Hydrocarbon Expression in Male and Female Drosophila melanogaster Genetics, March 1, 2007; 175(3): 1465 - 1477. [Abstract] [Full Text] [PDF]

International Access Link

Online ISSN: 1943-2361  Print ISSN: 0016-6731
Copyright 2009 by the Genetics Society of America
phone: 412-268-1812   fax: 412-268-1813   email: genetics-gsa{at}andrew.cmu.edu