Genetic dissection of hybrid incompatibilities between Drosophila simulans and D. mauritiana. II. Mapping hybrid male sterility loci on the third chromosome

Hybrid male sterility (HMS) is a rapidly evolving mechanism of reproductive isolation in Drosophila. Here we report a genetic analysis of HMS in third-chromosome segments of Drosophila mauritiana that were introgressed into a D. simulans background. Qualitative genetic mapping was used to localize 10 loci on 3R and a quantitative trait locus (QTL) procedure (multiple-interval mapping) was used to identify 19 loci on the entire chromosome. These genetic incompatibilities often show dominance and complex patterns of epistasis. Most of the HMS loci have relatively small effects and generally at least two or three of them are required to produce complete sterility. Only one small region of the third chromosome of D. mauritiana by itself causes a high level of infertility when introgressed into D. simulans. By comparison with previous studies of the X chromosome, we infer that HMS loci are only approximately 40% as dense on this autosome as they are on the X chromosome. These results are consistent with the gradual evolution of hybrid incompatibilities as a by-product of genetic divergence in allopatric populations.     

Quantitative trait loci for sexual isolation between Drosophila simulans and D. mauritiana

Sexual isolating mechanisms that act before fertilization are often considered the most important genetic barriers leading to speciation in animals. While recent progress has been made toward understanding the genetic basis of the postzygotic isolating mechanisms of hybrid sterility and inviability, little is known about the genetic basis of prezygotic sexual isolation. Here, we map quantitative trait loci (QTL) contributing to prezygotic reproductive isolation between the sibling species Drosophila simulans and D. mauritiana. We mapped at least seven QTL affecting discrimination of D. mauritiana females against D. simulans males, three QTL affecting D. simulans male traits against which D. mauritiana females discriminate, and six QTL affecting D. mauritiana male traits against which D. simulans females discriminate. QTL affecting sexual isolation act additively, are largely different in males and females, and are not disproportionately concentrated on the X chromosome: The QTL of greatest effect are located on chromosome 3. Unlike the genetic components of postzygotic isolation, the loci for prezygotic isolation do not interact epistatically. The observation of a few QTL with moderate to large effects will facilitate positional cloning of genes underlying sexual isolation.     

A quantitative genetic analysis of male sexual traits distinguishing the sibling species Drosophila simulans and D. sechellia

A quantitative trait locus (QTL) genetic analysis of morphological and reproductive traits distinguishing the sibling species Drosophila simulans and D. sechellia was carried out in a backcross design, using 38 markers with an average spacing of 8.4 cM. The direction of QTL effects for the size of the posterior lobe was consistent across the identified QTL, indicating directional selection for this trait. Directional selection also appears to have acted on testis length, indicating that sexual selection may have influenced many reproductive traits, although other forms of directional selection cannot be ruled out. Sex comb tooth number exhibited high levels of variation both within and among isofemale lines and showed no evidence for directional selection and, therefore, may not have been involved in the early speciation process. A database search for genes associated with significant QTL revealed a set of candidate loci for posterior lobe shape and size, sex comb tooth number, testis length, tibia length, and hybrid male fertility. In particular, decapentaplegic (dpp), a gene known to influence the genital arch, was found to be associated with the largest LOD peak for posterior lobe shape and size.     

Genetic architecture of variation in sex-comb tooth number in Drosophila simulans

The sex comb on the forelegs of Drosophila males is a secondary sexual trait, and the number of teeth on these combs varies greatly within and between species. To understand the relationship between the intra- and interspecific variation, we performed quantitative trait locus (QTL) analyses of the intraspecific variation in sex-comb tooth number. We used five mapping populations derived from two inbred Drosophila simulans strains that were divergent in the number of sex-comb teeth. Although no QTLs were detected on the X chromosome, we identified four QTLs on the second chromosome and three QTLs on the third chromosome. While identification and estimated effects of the second-chromosome QTLs depend on genetic backgrounds, significant and consistent effects of the two third-chromosome QTLs were found in two genetic backgrounds. There were significant epistatic interactions between a second-chromosome QTL and a third-chromosome QTL, as well as between two second-chromosome QTLs. The third-chromosome QTLs are concordant with the locations of the QTLs responsible for the previously observed differences in sex-comb tooth number between D. simulans and D. mauritiana.     

The genetics of mating recognition between Drosophila simulans and D. sechellia

During courtship, visual and chemical signals are often exchanged between the sexes. The proper exchange of such signals ensures intraspecific recognition. We have examined the genetic basis of interspecific differences in male mating behaviour and pheromone concentration between Drosophila simulans and D. sechellia by using Drosophila simulans/D. sechellia introgression lines. Our results show a majority of quantitative trait loci (QTLs) explaining variation in both male mating behaviour and pheromone concentration to be located on the third chromosome. One QTL found on the third chromosome explains variation in time needed to start courtship and copulation as well as time spent courting. The position of such QTL (approximately 84A-88B) with effects on courtship and copulation aspects of mating includes the candidate sex determination gene doublesex (84E5-6) and Voila (86E1-2), a gene that affects male courtship in D. melanogaster. One additional third chromosome QTL explained variation in 7-tricosene pheromone concentrations among males. The interval mapping position of this QTL (approximately 68E-76E) did not overlap with the position detected for differences in mating behaviour and the intervals did not include candidate genes previously identified as having an effect on D. melanogaster cuticular hydrocarbon production. We did not detect any directionality of the effect of Drosophila sechellia allele introgressions in male mating recognition.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

Genetics of sexual isolation in females of the Drosophila simulans species complex.

Genetic analysis of hybrids between Drosophila simulans and D. sechellia shows that sexual isolation in females is caused by at least two genes, one on each major autosome, while the X chromosome has no effect. These results are similar to those of a previous study of hybrids between D. simulans and another sibling species, D. mauritiana. In this latter hybridization, each arm of the second chromosome carries genes causing sexual isolation in females, implying a total divergence of at least three loci. The genetic similarity between the D. simulans/D. mauritiana and D. simulans/D. sechellia hybridizations probably results from independent evolution and not phylogenetic artifacts, because the dominance relationships and behavioural interactions differ between the two hybridizations. The lack of an X-chromosome effect on sexual isolation contrasts with genetic studies of post-zygotic reproductive isolation, which invariably show strong effects of this chromosome.     

The genetic basis of prezygotic reproductive isolation between Drosophila santomea and D. yakuba due to mating preference

Sexual isolating mechanisms that act before fertilization are often considered the most important genetic barriers leading to speciation in animals. While progress has been made toward understanding the genetic basis of the postzygotic isolating mechanisms of hybrid sterility and inviability, little is known about the genetic basis of prezygotic sexual isolation. Here, we map quantitative trait loci (QTL) contributing to prezygotic reproductive isolation between the sibling species Drosophila santomea and D. yakuba. We mapped at least three QTL affecting discrimination of D. santomea females against D. yakuba males: one X-linked and one autosomal QTL affected the likelihood of copulation, and a second X chromosome QTL affected copulation latency. Three autosomal QTL also affected mating success of D. yakuba males with D. santomea. No epistasis was detected between QTL affecting sexual isolation. The QTL do not overlap between males and females and are not disproportionately concentrated on the X chromosome. There was some overlap in map locations of QTL affecting sexual isolation between D. santomea and D. yakuba with QTL affecting sexual isolation between D. simulans and D. mauritiana and with QTL affecting differences in pigmentation between D. santomea and D. yakuba. Future high-resolution mapping and, ultimately, positional cloning, will reveal whether these traits do indeed have a common genetic basis.     

Mapping resistance to cereal aphids in barley

 A set of 150 doubled-haploid (DH) barley (Hordeum vulgare L.) lines derived from the cross of Harrington/TR306 was used to determine the number and chromosomal location of quantitative trait loci (QTLs) controlling resistance to cereal aphids. The experiments were conducted under natural infestation in the field during two growing seasons: 1994 and 1995. Aphid resistance was measured by counting the number of aphids per plot. Counts were made on a weekly basis. Each year at the time of maximum aphid density there was an obvious difference in reaction between the parental genotypes. The DH lines showed continuous variation for aphid density. Simple interval mapping and simplified composite interval mapping revealed that the principal QTL determining cereal aphid resistance is on the distal region of the short arm of chromosome 1. This aphid-resistance QTL is linked with a heading-date QTL. At the time of highest aphid infestation, this QTL accounted for 31% and 22% of the total variance of aphid density in 1994 and 1995, respectively. A QTL on chromosome 5 was also detected but only by simplified composite interval mapping. However, the largest consistent effect was due to the QTL on the short arm of chromosome 1. This QTL may be a useful target for marker-assisted selection for adult plant cereal aphid resistance in barley. Received: 10 September 1996/Accepted: 11 October 1996     

Do quantitative trait loci (QTL) for a courtship song difference between Drosophila simulans and D. sechellia coincide with candidate genes and intraspecific QTL?

The genetic architecture of traits influencing sexual isolation can give insight into the evolution of reproductive isolation and hence speciation. Here we report a quantitative trait loci (QTL) analysis of the difference in mean interpulse interval (IPI), an important component of the male courtship song, between Drosophila simulans and D. sechellia. Using a backcross analysis, we find six QTL that explain a total of 40.7% of the phenotypic variance. Three candidate genes are located in the intervals bounded by two of the QTL and there are no significant QTL on the X chromosome. The values of mean IPI for hybrid individuals imply the presence of dominant alleles or epistasis. Because unisexual hybrid sterility prevents an F(2) analysis, we cannot distinguish dominant from additive genetic effects at the scale of QTL. A comparison with a study of QTL for intraspecific variation in D. melanogaster shows that, for these strains, the QTL we have identified for interspecific variation cannot be those that contribute to intraspecific variation. We find that the QTL have bidirectional effects, which indicates that the genetic architecture is compatible with divergence due to genetic drift, although other possibilities are discussed.     

High-resolution mapping of quantitative trait loci for sternopleural bristle number in Drosophila melanogaster.

We have mapped quantitative trait loci (QTL) harboring naturally occurring allelic variation for Drosophila bristle number. Lines with high (H) and low (L) sternopleural bristle number were derived by artificial selection from a large base population. Isogenic H and L sublines were extracted from the selection lines, and populations of X and third chromosome H/L recombinant isogenic lines were constructed in the homozygous low line background. The polymorphic cytological locations of roo transposable elements provided a dense molecular marker map with an average intermarker distance of 4.5 cM. Two X chromosome and six chromosome 3 QTL affecting response to selection for sternopleural bristle number and three X chromosome and three chromosome 3 QTL affecting correlated response in abdominal bristle number were detected using a composite interval mapping method. The average effects of bristle number QTL were moderately large, and some had sex-specific effects. Epistasis between QTL affecting sternopleural bristle number was common, and interaction effects were large. Many of the intervals containing bristle number QTL coincided with those mapped in previous studies. However, resolution of bristle number QTL to the level of genetic loci is not trivial, because the genomic regions containing bristle number QTL often did not contain obvious candidate loci, and results of quantitative complementation tests to mutations at candidate loci affecting adult bristle number were ambiguous.     

The genetic architecture of Drosophila sensory bristle number

We have mapped quantitative trait loci (QTL) for Drosophila mechanosensory bristle number in six recombinant isogenic line (RIL) mapping populations, each of which was derived from an isogenic chromosome extracted from a line selected for high or low, sternopleural or abdominal bristle number and an isogenic wild-type chromosome. All RILs were evaluated as male and female F(1) progeny of crosses to both the selected and the wild-type parental chromosomes at three developmental temperatures (18 degrees, 25 degrees, and 28 degrees ). QTL for bristle number were mapped separately for each chromosome, trait, and environment by linkage to roo transposable element marker loci, using composite interval mapping. A total of 53 QTL were detected, of which 33 affected sternopleural bristle number, 31 affected abdominal bristle number, and 11 affected both traits. The effects of most QTL were conditional on sex (27%), temperature (14%), or both sex and temperature (30%). Epistatic interactions between QTL were also common. While many QTL mapped to the same location as candidate bristle development loci, several QTL regions did not encompass obvious candidate genes. These features are germane to evolutionary models for the maintenance of genetic variation for quantitative traits, but complicate efforts to understand the molecular genetic basis of variation for complex traits.     

Genetic mapping of adaptive wing size variation in Drosophila simulans

Many ecologically important traits exhibit latitudinal variation. Body size clines have been described repeatedly in insects across multiple continents, suggesting that similar selective forces are shaping these geographical gradients. It is unknown whether these parallel clinal patterns are controlled by the same or different genetic mechanism(s). We present here, quantitative trait loci (QTL) analysis of wing size variation in Drosophila simulans. Our results show that much of the wing size variation is controlled by a QTL on Chr 3L with relatively minor contribution from other chromosome arms. Comparative analysis of the genomic positions of the QTL indicates that the major QTL on Chr 3 are distinct in D. simulans and D. melanogaster, whereas the QTL on Chr 2R might overlap between species. Our results suggest that parallel evolution of wing size clines could be driven by non-identical genetic mechanisms but in both cases involve a major QTL as well as smaller effects of other genomic regions.     

Variation in fiber number of a male-specific muscle between Drosophila species: a genetic and developmental analysis

We characterize a newly discovered morphological difference between species of the Drosophila melanogaster subgroup. The muscle of Lawrence (MOL) contains about four to five fibers in D. melanogaster and Drosophila simulans and six to seven fibers in Drosophila mauritiana and Drosophila sechellia. The same number of nuclei per fiber is present in these species but their total number of MOL nuclei differs. This suggests that the number of muscle precursor cells has changed during evolution. Our comparison of MOL development indicates that the species difference appears during metamorphosis. We mapped the quantitative trait loci responsible for the change in muscle fiber number between D. sechellia and D. simulans to two genomic regions on chromosome 2. Our data eliminate the possibility of evolving mutations in the fruitless gene and suggest that a change in the twist might be partly responsible for this evolutionary change.     

The genetic basis of behavioral isolation between Drosophila mauritiana and D. sechellia

Understanding how species form is a fundamental question in evolutionary biology. Identifying the genetic bases of barriers that prevent gene flow between species provides insight into how speciation occurs. Here, I analyze a poorly understood reproductive isolating barrier, prezygotic reproductive isolation. I perform a genetic analysis of prezygotic isolation between two closely related species of Drosophila, D. mauritiana and D. sechellia. I first confirm the existence of strong behavioral isolation between D. mauritiana females and D. sechellia males. Next, I examine the genetic basis of behavioral isolation by (1) scanning an existing set of introgression lines for chromosomal regions that have a large effect on isolation; and (2) mapping quantitative trait loci (QTL) that underlie behavioral isolation via backcross analysis. In particular, I map QTL that determine whether a hybrid backcross female and a D. sechellia male will mate. I identify a single significant QTL, on the X chromosome, suggesting that few major-effect loci contribute to behavioral isolation between these species. In further work, I refine the map position of the QTL to a small region of the X chromosome.     

Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species

The Dobzhansky-Muller model posits that postzygotic reproductive isolation results from the evolution of incompatible epistatic interactions between species: alleles that function in the genetic background of one species can cause sterility or lethality in the genetic background of another species. Progress in identifying and characterizing factors involved in postzygotic isolation in Drosophila has remained slow, mainly because Drosophila melanogaster, with all of its genetic tools, forms dead or sterile hybrids when crossed to its sister species, D. simulans, D. sechellia, and D. mauritiana. To circumvent this problem, we used chromosome deletions and duplications from D. melanogaster to map two hybrid incompatibility loci in F(1) hybrids with its sister species. We mapped a recessive factor to the pericentromeric heterochromatin of the X chromosome in D. simulans and D. mauritiana, which we call heterochromatin hybrid lethal (hhl), which causes lethality in F(1) hybrid females with D. melanogaster. As F(1) hybrid males hemizygous for a D. mauritiana (or D. simulans) X chromosome are viable, the lethality of deficiency hybrid females implies that a dominant incompatible partner locus exists on the D. melanogaster X. Using small segments of the D. melanogaster X chromosome duplicated onto the Y chromosome, we mapped a dominant factor that causes hybrid lethality to a small 24-gene region of the D. melanogaster X. We provide evidence suggesting that it interacts with hhl(mau). The location of hhl is consistent with the emerging theme that hybrid incompatibilities in Drosophila involve heterochromatic regions and factors that interact with the heterochromatin.     

The genetic architecture of selection response. Inferences from fine-scale mapping of bristle number quantitative trait loci in Drosophila melanogaster.

Quantitative trait loci (QTL) affecting responses and correlated responses to selection for abdominal and sternopleural bristle number have been mapped with high resolution to the X and third chromosomes. Advanced intercross recombinant isogenic chromosomes were constructed from high and low selection lines in an unselected inbred background, and QTL were detected using composite interval mapping and high density transposable element marker maps. We mapped a total of 26 bristle number QTL with large effects, which were in or immediately adjacent to intervals previously inferred to contain bristle number QTL on these chromosomes. The QTL contributing to response to selection for high bristle number were not the same as those contributing to response to selection for low bristle number, suggesting that distributions of allelic effects per locus may be asymmetrical. Correlated responses were more often attributable to loose linkage than pleiotropy or close linkage. Bristle number QTL mapping to the same locations have been inferred in studies with different parental strains. Of the 26 QTL, 20 mapped to locations consistent with candidate genes affecting peripheral nervous system development and/or bristle number. This facilitates determining the molecular basis of quantitative variation and allele frequencies by associating molecular variation at the candidate genes with phenotypic variation in bristle number in samples of alleles from nature.