Genetic studies of secondary sexual traits provide insights into whether and how selection drove their divergence among populations, and these studies often focus on the fraction of variation attributable to genes on the X-chromosome. However, such studies may sometimes misinterpret the amount of variation attributable to the X-chromosome if using only simple reciprocal F1 crosses, or they may presume sexual selection has affected the observed phenotypic variation. We examined the genetics of a secondary sexual trait, male sex comb size, in Drosophila subobscura. This species bears unusually large sex combs for its species group, and therefore, this trait may be a good candidate for having been affected by natural or sexual selection. We observed significant heritable variation in number of teeth of the distal sex comb across strains. While reciprocal F1 crosses seemed to implicate a disproportionate X-chromosome effect, further examination in the F2 progeny showed that transgressive autosomal effects inflated the estimate of variation associated with the X-chromosome in the F1. Instead, the X-chromosome appears to confer the smallest contribution of all major chromosomes to the observed phenotypic variation. Further, we failed to detect effects on copulation latency or duration associated with the observed phenotypic variation. Overall, this study presents an examination of the genetics underlying segregating phenotypic variation within species and illustrates two common pitfalls associated with some past studies of the genetic basis of secondary sexual traits.

Speciation, the evolution of reproductive isolation amongst populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting migration locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination- and mutation-rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarise the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a roadmap for future speciation research; a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.

Fine scale meiotic recombination maps have uncovered a large amount of variation in crossover rate across the genomes of many species, and such variation in mammalian and yeast genomes is concentrated to <5 kb regions of highly elevated recombination rates (10–100x the background rate) called “hotspots.” Drosophila exhibit substantial recombination rate heterogeneity across their genome, but evidence for these highly-localized hotspots is lacking. We assayed recombination across a 40 Kb region of Drosophila pseudoobscura chromosome 2, with one 20 kb interval assayed every 5 Kb and the adjacent 20 kb interval bisected into 10 kb pieces. We found that recombination events across the 40 kb stretch were relatively evenly distributed across each of the 5 kb and 10 kb intervals, rather than concentrated in a single 5 kb region. This, in combination with other recent work, indicates that the recombination landscape of Drosophila may differ from the punctate recombination pattern observed in many mammals and yeast. Additionally, we found no correlation of average pairwise nucleotide diversity and divergence with recombination rate across the 20 kb intervals, nor any effect of maternal age in weeks on recombination rate in our sample.

In contrast to the prevailing dogma in the 1990s, recent studies have suggested that an evolutionary history of segregation distortion within species may contribute to sterility in species hybrids. However, this recent work identified segregation distortion exclusively in species hybrids which may never have had an evolutionary history of segregation distortion in either parent species. We expand on previous work by using a strain of Drosophila persimilis exhibiting segregation distortion within species to generate QTL maps for segregation distortion and hybrid sterility in crosses between D. persimilis and D. pseudoobscura. The maps localize regions along the XR contributing to both phenotypes, and they indicate one region of overlap between the two maps. This overlap could provide preliminary evidence for an association between segregation distortion within species and hybrid sterility, but the localizations are currently too broad to have confidence in this conclusion. This work is a first step towards possibly supporting a genetic conflict model of speciation in this system.

Chromosomal inversions disrupt recombination in heterozygotes by both reducing crossing over within inverted regions and increasing it elsewhere in the genome. The reduction of recombination in inverted regions facilitates the maintenance of hybridizing species, as outlined by various models of chromosomal speciation. We present a comprehensive comparison of the effects of inversions on recombination rates and on nucleotide divergence. Within an inversion differentiating Drosophila pseudoobscura and D. persimilis, we detected one double-recombinant among 9739 progeny from F1 hybrids screened, consistent with published double crossover frequencies observed within species. Despite similar rates of exchange within and between species, we found no sequence-based evidence of ongoing gene exchange between species within this inversion, but significant exchange was inferred within species. We also observed greater differentiation at regions near inversion breakpoints between species vs. within species. Moreover, we observed strong ‘interchromosomal effect’ (higher recombination in inversion heterozygotes between species) with up to 9-fold higher recombination rates along collinear segments of chromosome two in hybrids. Further, we observed that regions most susceptible to changes in recombination rates corresponded to regions with lower recombination rates in homokaryotypes. Finally, we showed that interspecies nucleotide divergence is lower in regions with greater increases in recombination rate, potentially resulting from greater interspecies exchange. Overall, we have identified several similarities and differences between inversions segregating within vs. between species in their effects on recombination and divergence. We conclude that these differences are most likely due to lower frequency of heterokaryotypes and to fitness consequences from the accumulation of various incompatibilities between species. Additionally, we have identified possible effects of inversions on interspecies gene exchange that had not been considered previously.

Taxa in the early stages of speciation may bear intraspecific allelic variation at loci conferring barrier traits in hybrids such as hybrid sterility. Additionally, hybridization may spread alleles that confer barrier traits to other taxa. Historically, few studies examine within- and between-species variation at loci conferring reproductive isolation. Here, we test for allelic variation within Drosophila persimilis and within the Bogota subspecies of D. pseudoobscura at regions previously shown to contribute to hybrid male sterility. We also test whether D. persimilis and the USA subspecies of D. pseudoobscura share an allele conferring hybrid sterility in a D. pseudoobscura bogotana genetic background. All loci conferred similar hybrid sterility effects across all strains studied, though we detected some statistically significant quantitative effect variation among D. persimilis alleles of some hybrid incompatibility QTLs. We also detected allelism between D. persimilis and D. pseudoobscura USA at a 2nd chromosome hybrid sterility QTL. We hypothesize that either the QTL is ancestral in D. persimilis and D. pseudoobscura USA and lost in D. pseudoobscura bogotana, or gene flow transferred the QTL from D. persimilis to D. pseudoobscura USA. We discuss our findings in the context of population features that may contribute to variation in hybrid incompatibilities.

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed “recombinational neighborhoods,” where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5–7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.

We carried out a three-tiered genetic analysis of egg-to-adult development time and viability in ancestral and derived populations of cactophilic D. 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 cacti revealed significant host plant X region and population interactions for viability and development time. Second, replicated line crosses with cactus-reared flies revealed autosomal, X chromosome, cytoplasmic, and autosome X cactus influences on development time. Third, a QTL analysis of development time differences on 1688 Baja X mainland F2 males revealed eight QTL. Eight GxE interactions were also detected, caused by longer development times associated with mainland alleles reared a mainland host with smaller differences among Baja genotypes on a Baja host plant. Four QTL influenced both development time and cuticular 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 Mexico from Baja California by shifting host plants were genetically correlated with variation in cuticular hydrocarbon-based mate preferences.