We performed a quantitative trait locus (QTL) analysis of epicuticular hydrocarbon variation in 1650 F 2 males from crosses of Baja California and mainland Mexico populations of Drosophila mojavensis cultured on two major host cacti. Principal component (PC) analysis revealed five PCs that accounted for 82% of the total epicuticular hydrocarbon variation. Courtship trials with mainland females were used to characterize hydrocarbon profiles of mated and unmated F 2 males, and logistic regression analysis showed that cactus substrates, two PCs, and a PC by cactus interaction were associated with mating success. Multiple QTLs were detected for each hydrocarbon PC and seven G × E (cactus) interactions were uncovered for the X, second, and fourth chromosomes. Males from the courtship trials and virgins were used, so "exposure to females" was included as a factor in QTL analyses. "Exposed" males expressed significantly different hydrocarbon profiles than virgins for most QTLs, particularly for the two PCs associated with mating success. Ten QTLs showed G × E (exposure) interactions with most resulting from mainland genotypes expressing altered hydrocarbon amounts when exposed to females compared to Baja genotypes. Many cactus × exposure interaction terms detected across QTL and all PCs confirmed that organ pipe-reared males expressed significantly lower hydrocarbon amounts when exposed to females than when reared on agria cactus. Epicuticular hydrocarbon variation in D. mojavensis is therefore a multigenic trait with some epistasis, multiple QTLs exhibited pleiotropy, correlated groups of hydrocarbons and cactus substrates determined courtship success, and males altered their hydrocarbon profiles in response to females.

Few studies have examined genotype by environment (GxE) effects on premating reproductive isolation and associated behaviors, even though such effects may be common when speciation is driven by adaptation to different environments. In this study, mating success and courtship song differences among diverging populations of Drosophila mojavensis were investigated in a two-environment quantitative trait locus (QTL) analysis. Baja California and mainland Mexico populations of D. mojavensis feed and breed on different host cacti, so these host plants were used to culture F 2 males to examine host-specific QTL effects and GxE interactions influencing mating success and courtship songs. Linear selection gradient analysis showed that mainland females mated with males that produced songs with significantly shorter L(long)-IPIs, burst durations, and interburst intervals. Twenty-one microsatellite loci distributed across all five major chromosomes were used to localize effects of mating success, time to copulation, and courtship song components. Male courtship success was influenced by a single detected QTL, the main effect of cactus, and four GxE interactions, whereas time to copulation was influenced by three different QTLs on the fourth chromosome. Multiplelocus restricted maximum likelihood (REML) analysis of courtship song revealed consistent effects linked with the same fourth chromosome markers that influenced time to copulation, a number of GxE interactions, and few possible cases of epistasis. GxE interactions for mate choice and song can maintain genetic variation in populations, but alter outcomes of sexual selection and isolation, so signal evolution and reproductive isolation may be slowed in diverging populations. Understanding the genetics of incipient speciation in D. mojavensis clearly depends on cactus-specific expression of traits associated with courtship behavior and sexual isolation.

We carried out a three-tiered genetic analysis of egg-to-adult development time and viability in ancestral and derived populations of cactophilic Drosophila mojavensis to test the hypothesis that evolution of these life-history characters has shaped premating reproductive isolation in this species. First, a common garden experiment with 11 populations from Baja California and mainland Mexico and Arizona reared on two host species revealed significant host plant X region and population interactions for viability and development time, evidence for host plant adaptation. Second, replicated line crosses with flies reared on both hosts revealed autosomal, X chromosome, cytoplasmic, and autosome X cactus influences on development time. Viability differences were influenced by host plants, autosomal dominance, and X chromosomal effects. Many of the F 1 , F 2 , and backcross generations showed evidence of heterosis for viability. Third, a QTL analysis of male courtship song and epicuticular hydrocarbon variation based on 1688 Baja × mainland F 2 males also revealed eight QTL influencing development time differences. Mainland alleles at six of these loci were associated with longer development times, consistent with population-level differences. Eight G × E interactions were also detected caused by longer development times of mainland alleles expressed on a mainland host with smaller differences among Baja genotypes reared on the Baja host plant. Four QTL influenced both development time and epicuticular hydrocarbon differences associated with courtship success, and there was a significant QTL-based correlation between development time and cuticular hydrocarbon variation. Thus, the regional shifts in life histories that evolved once D. mojavensis invaded mainland

We used whole transcriptome microarrays to assess changes in gene expression and monitored mortality rates and epicuticular hydrocarbons (CHCs) in response to desiccation stress in four natural populations of Drosophila mojavensis from Baja California and mainland Mexico. Desiccation had the greatest effect on gene expression, followed by biogeographical variation at regional and population levels. Genes involved in environmental sensing and cuticular structure were up-regulated in dry conditions, while genes involved in transcription itself were down-regulated. Flies from Baja California had higher expression of reproductive and mitochondrial genes, suggesting that these populations have greater fecundity and higher metabolic rates. Host plant differences had a surprisingly minor effect on the transcriptome. In most cases, desiccation-caused mortality was greater in flies reared on fermenting cactus tissues than laboratory media. Water content of adult females and males was significantly different, and was lower in Baja California males. Different groups of CHCs simultaneously increased and decreased in amounts due to desiccation exposure of 9 and 18 hr and were population-specific and dependent on larval rearing substrates. Overall, we observed that changes in gene expression involved a coordinated response of behavioral, cuticular and metabolic genes. Together with differential expression of cuticular hydrocarbons, this study revealed some of the mechanisms that have allowed D. mojavensis to exploit its harsh desert conditions. Certainly, for D. mojavensis that uses different host plants, population-level understanding of responses to stressors associated with future climate change in desert regions must be evaluated across geographical and local ecological scales.

BackgroundWe characterized variation and chemical composition of epicuticular hydrocarbons (CHCs) in the seven species of the Drosophila buzzatii cluster with gas chromatography/mass spectrometry. Despite the critical role of CHCs in providing resistance to desiccation and involvement in communication, such as courtship behavior, mating, and aggregation, few studies have investigated how CHC profiles evolve within and between species in a phylogenetic context. We analyzed quantitative differences in CHC profiles in populations of the D. buzzatii species cluster in order to assess the concordance of CHC differentiation with species divergence.ResultsThirty-six CHC components were scored in single fly extracts with carbon chain lengths ranging from C29 to C39, including methyl-branched alkanes, n-alkenes, and alkadienes. Multivariate analysis of variance revealed that CHC amounts were significantly different among all species and canonical discriminant function (CDF) analysis resolved all species into distinct, non-overlapping groups. Significant intraspecific variation was found in different populations of D. serido suggesting that this taxon is comprised of at least two species. We summarized CHC variation using CDF analysis and mapped the first five CHC canonical variates (CVs) onto an independently derived period (per) gene + chromosome inversion + mtDNA COI gene for each sex. We found that the COI sequences were not phylogenetically informative due to introgression between some species, so only per + inversion data were used. Positive phylogenetic signal was observed mainly for CV1 when parsimony methods and the test for serial independence (TFSI) were used. These results changed when no outgroup species were included in the analysis and phylogenetic signal was then observed for female CV3 and/or CV4 and male CV4 and CV5. Finally, removal of divergent populations of D. serido significantly increased the amount of phylogenetic signal as up to four out of five CVs then displayed positive phylogenetic signal.ConclusionsCHCs were conserved among species while quantitative differences in CHC profiles between populations and species were statistically significant. Most CHCs were species-, population-, and sex-specific. Mapping CHCs onto an independently derived phylogeny revealed that a significant portion of CHC variation was explained by species' systematic affinities indicating phylogenetic conservatism in the evolution of these hydrocarbon arrays, presumptive waterproofing compounds and courtship signals as in many other drosophilid species.