Search

Female mate choice plays an important role in sexual isolation. The present study examined sexual isolation using D. simulans and D. mauritiana in conditions where females had no opportunity to compare males versus where females were able to choose males. The sound produced by wing vibration in males (courtship song) affects female receptivity in Drosophila. Females of both species copulated with intact conspecific males more than intact heterospecific males and wingless conspecific and heterospecific males. Drosophila mauritiana females copulated only with intact conspecific males within 30-min observations period without comparing other males, suggesting that absolute criteria are used for decision-making to accept courting males. Females of D. simulans copulated with intact D. mauritiana males as well as wingless D. simulans or D. mauritiana males in no-choice conditions. In a choice situation, D. simulans females copulated with intact D. mauritiana males as well as wingless D. mauritiana males when the females were courted by both types of males, suggesting that D. simulans females accept intact D. mauritiana males as if they are mute. Females of D. simulans copulated with intact D. simulans males as well as intact D. mauritiana males when they were courted by males of either type, whereas they copulated with intact D. simulans males more than intact D. mauritiana males in true choice situations. These results suggest that females make a comparative review of courting males before accepting a male and that conspecific song is a factor in criteria affecting female selectivity.

Behavioral isolation is considered to be the primary mode of species isolation, and the lack of identification of individual genes for behavioral isolation has hindered our ability to address fundamental questions about the process of speciation. One of the major questions that remains about behavioral isolation is whether the genetic basis of isolation between species also varies within a species. Indeed, the extent to which genes for isolation may vary across a population is rarely explored. Here, we bypass the problem of individual gene identification by addressing this question using a quantitative genetic comparison. Using strains from eight different populations of Drosophila simulans, we genetically mapped the genomic regions contributing to behavioral isolation from their closely related sibling species, Drosophila mauritiana. We found extensive variation in the size of contribution of different genomic regions to behavioral isolation among the different strains, in the location of regions contributing to isolation, and in the ability to redetect loci when retesting the same strain.

Significance The morphology of male genitalia evolves rapidly, probably driven by sexual selection. However, little is known about the genes underlying genitalia differences between species. Identifying these genes is key to understanding how sexual selection acts to produce rapid phenotypic change. We have found that the gene tartan underlies differences between male Drosophila mauritiana and Drosophila simulans in the size and bristle number of the claspers—genital projections that grasp the female during copulation. Moreover, since tartan encodes a protein that is involved in cell interactions, this may represent an alternative developmental mechanism for morphological change. Therefore, our study provides insights into the genetic and developmental bases for the rapid evolution of male genitalia and organ size more generally., Male genital structures are among the most rapidly evolving morphological traits and are often the only features that can distinguish closely related species. This process is thought to be driven by sexual selection and may reinforce species separation. However, while the genetic bases of many phenotypic differences have been identified, we still lack knowledge about the genes underlying evolutionary differences in male genital organs and organ size more generally. The claspers (surstyli) are periphallic structures that play an important role in copulation in insects. Here, we show that divergence in clasper size and bristle number between Drosophila mauritiana and Drosophila simulans is caused by evolutionary changes in tartan (trn), which encodes a transmembrane leucine-rich repeat domain protein that mediates cell–cell interactions and affinity. There are no fixed amino acid differences in trn between D. mauritiana and D. simulans, but differences in the expression of this gene in developing genitalia suggest that cis-regulatory changes in trn underlie the evolution of clasper morphology in these species. Finally, analyses of reciprocal hemizygotes that are genetically identical, except for the species from which the functional allele of trn originates, determined that the trn allele of D. mauritiana specifies larger claspers with more bristles than the allele of D. simulans. Therefore, we have identified a gene underlying evolutionary change in the size of a male genital organ, which will help to better understand not only the rapid diversification of these structures, but also the regulation and evolution of organ size more broadly.

Genes involved in immune defense against pathogens provide some of the most well-known examples of both directional and balancing selection. Antimicrobial peptides (AMPs) are innate immune effector genes, playing a key role in pathogen clearance in many species, including Drosophila. Conflicting lines of evidence have suggested that AMPs may be under directional, balancing, or purifying selection. Here, we use both a linear model and control-gene-based approach to show that balancing selection is an important force shaping AMP diversity in Drosophila. In Drosophila melanogaster, this is most clearly observed in ancestral African populations. Furthermore, the signature of balancing selection is even more striking once background selection has been accounted for. Balancing selection also acts on AMPs in Drosophila mauritiana, an isolated island endemic separated from D. melanogaster by about 4 Myr of evolution. This suggests that balancing selection may be broadly acting to maintain adaptive diversity in Drosophila AMPs, as has been found in other taxa.

Symbiosis is a major force of evolutionary change, influencing virtually all aspects of biology, from population ecology and evolution to genomics and molecular/biochemical mechanisms of development and reproduction. A remarkable example is Wolbachia endobacteria, present in some parasitic nematodes and many arthropod species. Acquisition of genomic data from diverse Wolbachia clades will aid in the elucidation of the different symbiotic mechanisms(s). However, challenges of de novo assembly of Wolbachia genomes include the presence in the sample of host DNA: nematode/vertebrate or insect. We designed biotinylated probes to capture large fragments of Wolbachia DNA for sequencing using PacBio technology (LEFT-SEQ: Large Enriched Fragment Targeted Sequencing). LEFT-SEQ was used to capture and sequence four Wolbachia genomes: the filarial nematode Brugia malayi, wBm, (21-fold enrichment), Drosophila mauritiana flies (2 isolates), wMau (11-fold enrichment), and Aedes albopictus mosquitoes, wAlbB (200-fold enrichment). LEFT-SEQ resulted in complete genomes for wBm and for wMau. For wBm, 18 single-nucleotide polymorphisms (SNPs), relative to the wBm reference, were identified and confirmed by PCR. A limit of LEFT-SEQ is illustrated by the wAlbB genome, characterized by a very high level of insertion sequences elements (ISs) and DNA repeats, for which only a 20-contig draft assembly was achieved.

Small noncoding piRNAs act as sequence-specific guides to repress complementary targets in Metazoa. Prior studies in Drosophila ovaries have demonstrated the function of the piRNA pathway in transposon silencing and therefore genome defense. However, the ability of the piRNA program to respond to different transposon landscapes and the role of piRNAs in regulating host gene expression remain poorly understood. Here, we comprehensively analyzed piRNA expression and defined the repertoire of their targets in Drosophila melanogaster testes. Comparison of piRNA programs between sexes revealed sexual dimorphism in piRNA programs that parallel sex-specific transposon expression. Using a novel bioinformatic pipeline, we identified new piRNA clusters and established complex satellites as dual-strand piRNA clusters. While sharing most piRNA clusters, the two sexes employ them differentially to combat the sex-specific transposon landscape. We found two piRNA clusters that produce piRNAs antisense to four host genes in testis, including CG12717/pirate, a SUMO protease gene. piRNAs encoded on the Y chromosome silence pirate, but not its paralog, to exert sex- and paralog-specific gene regulation. Interestingly, pirate is targeted by endogenous siRNAs in a sibling species, Drosophila mauritiana, suggesting distinct but related silencing strategies invented in recent evolution to regulate a conserved protein-coding gene.

The compound eyes of insects exhibit extensive variation in ommatidia number and size, which affects how they see and underlies adaptations in their vision to different environments and lifestyles. However, very little is known about the genetic and developmental bases underlying differences in compound eye size. We previously showed that the larger eyes ofDrosophila mauritianacompared toD. simulansis caused by differences in ommatidia size rather than number. Furthermore, we identified an X-linked chromosomal region inD. mauritianathat results in larger eyes when introgressed intoD. simulans. Here, we used a combination of fine-scale mapping and gene expression analysis to further investigate positional candidate genes on the X chromosome. We found thatorthodenticleis expressed earlier inD. mauritianathan inD. simulansduring ommatidial maturation in third instar larvae, and we further show that this gene is required for the correct organisation and size of ommatidia inD. melanogaster. Using ATAC-seq, we have identified several candidate eye enhancers ofotdas well as potential direct targets of this transcription factor that are differentially expressed betweenD. mauritianaandD. simulans. Taken together, our results suggest that differential timing ofotdexpression contributes to natural variation in ommatidia size betweenD. mauritianaandD. simulans, which provides new insights into the mechanisms underlying the regulation and evolution of compound eye size in insects.

Microbes infect a diversity of species, influencing the performance and fitness of their hosts. Maternally transmitted Wolbachia bacteria infect most insects and other arthropods, making these bacteria some of the most common endosymbionts in nature. Despite their global prevalence, it remains mostly unknown how Wolbachia influence host physiology and behavior to proliferate. We demonstrate pervasive effects of Wolbachia on Drosophila temperature preference. Most hosts infected with A-group Wolbachia prefer cooler temperatures, whereas the one host species infected with divergent B-group Wolbachia prefers warmer temperatures, relative to uninfected genotypes. Changes to host temperature preference generally do not alter Wolbachia abundance in host tissues, but for some A-group strains, adult males have increased Wolbachia titer when shifted to a cooler temperature. This suggests that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results help elucidate the impact of endosymbionts on their hosts amid the global Wolbachia pandemic., Heritable symbionts can modify a range of ecologically important host traits, including behavior. About half of all insect species are infected with maternally transmitted Wolbachia, a bacterial endosymbiont known to alter host reproduction, nutrient acquisition, and virus susceptibility. Here, we broadly test the hypothesis that Wolbachia modifies host behavior by assessing the effects of eight different Wolbachia strains on the temperature preference of six Drosophila melanogaster subgroup species. Four of the seven host genotypes infected with A-group Wolbachia strains (wRi in Drosophila simulans, wHa in D. simulans, wSh in Drosophila sechellia, and wTei in Drosophila teissieri) prefer significantly cooler temperatures relative to uninfected genotypes. Contrastingly, when infected with divergent B-group wMau, Drosophila mauritiana prefers a warmer temperature. For most strains, changes to host temperature preference do not alter Wolbachia titer. However, males infected with wSh and wTei tend to experience an increase in titer when shifted to a cooler temperature for 24 h, suggesting that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results indicate that Wolbachia modifications to host temperature preference are likely widespread, which has important implications for insect thermoregulation and physiology. Understanding the fitness consequences of these Wolbachia effects is crucial for predicting evolutionary outcomes of host-symbiont interactions, including how Wolbachia spreads to become common.

SUMMARYSmall non-coding piRNAs act as sequence-specific guides to repress complementary targets in Metazoa. Prior studies inDrosophilaovaries have demonstrated the function of piRNA pathway in transposon silencing and therefore genome defense. However, the ability of piRNA program to respond to different transposon landscape and the role of piRNAs in regulating host gene expression remain poorly understood. Here, we comprehensively analyzed piRNA expression and defined the repertoire of their targets inDrosophila melanogastertestes. Comparison of piRNA programs between sexes revealed sexual dimorphism in piRNA programs that parallel sex-specific transposon expression. Using a novel bioinformatic pipeline, we identified new piRNA clusters and established complex satellites as dual-strand piRNA clusters. While sharing most piRNA clusters, two sexes employ them differentially to combat sex-specific transposon landscape. We found several host genes targeted by piRNAs in testis, includingCG12717/pita, a SUMO protease gene. piRNAs encoded on Y chromosome silencepita, but not its paralog, to exert sex- and paralog-specific gene regulation. Interestingly,pitais targeted by endogenous siRNAs in a sibling species,Drosophila mauritiana, suggesting distinct but related silencing strategies invented in recent evolution to regulate a conserved protein-encoding gene.

Rapid divergence in genital structures among nascent species has been posited to be an early-evolving cause of reproductive isolation, although evidence supporting this idea as a widespread phenomenon remains mixed. Using a collection of interspecific introgression lines between twoDrosophilaspecies that diverged ∼240,000 years ago, we tested the hypothesis that even modest divergence in genital morphology can result in substantial fitness losses. We studied the reproductive consequences of variation in the male epandrial posterior lobes betweenDrosophila mauritianaandD. sechelliaand found that divergence in posterior lobe morphology has significant fitness costs on several pre-fertilization and post-copulatory reproductive measures. Males with divergent posterior lobe morphology also significantly reduced the life span of their mates. Interestingly, one of the consequences of genital divergence was decreased oviposition and fertilization, which suggests that a sensory bias for posterior lobe morphology could exist in females, and thus posterior lobe morphology may be the target of cryptic female choice in these species. Our results provide evidence that divergence in genitalia can in fact give rise to substantial reproductive isolation early during species divergence, and they also reveal novel reproductive functions of the external male genitalia inDrosophila.

Maternally transmittedWolbachiabacteria infect about half of all insect species. ManyWolbachiacause cytoplasmic incompatibility (CI), reduced egg hatch when uninfected females mate with infected males. Although CI produces a frequency-dependent fitness advantage that leads to high equilibriumWolbachiafrequencies, it does not aidWolbachiaspread from low frequencies. Indeed, the fitness advantages that produce initialWolbachiaspread and maintain non-CIWolbachiaremain elusive.wMauWolbachiainfectingDrosophila mauritianado not cause CI, despite being very similar to CI-causingwNo fromD. simulans(0.068% sequence divergence over 682,494 bp), suggesting recent CI loss. Using draftwMau genomes, we identify a deletion in a CI-associated gene, consistent with theory predicting that selection within host lineages does not act to increase or maintain CI. In the laboratory,wMau shows near-perfect maternal transmission; but we find no significant effect on host fecundity, in contrast to published data. IntermediatewMau frequencies on the island Mauritius are consistent with a balance between unidentified small, positive fitness effects and imperfect maternal transmission. Our phylogenomic analyses suggest that group-BWolbachia, includingwMau andwPip, diverged from group-AWolbachia, such aswMel andwRi, 6–46 million years ago, more recently than previously estimated.

The piRNA pathway is an adaptive mechanism that maintains genome stability by repression of selfish genomic elements. In the male germline of Drosophila melanogaster repression of Stellate genes by piRNAs generated from Supressor of Stellate (Su(Ste)) locus is required for male fertility, but both Su(Ste) piRNAs and their targets are absent in other Drosophila species. We found that D. melanogaster genome contains multiple X-linked non-coding genomic repeats that have sequence similarity to the protein-coding host gene vasa. In the male germline, these vasa-related AT-chX repeats produce abundant piRNAs that are antisense to vasa; however, vasa mRNA escapes silencing due to imperfect complementarity to AT-chX piRNAs. Unexpectedly, we discovered AT-chX piRNAs target vasa of Drosophila mauritiana in the testes of interspecies hybrids. In the majority of hybrid flies, the testes were strongly reduced in size and germline content. A minority of hybrids maintained wild-type array of premeiotic germ cells in the testes, but in them harmful Stellate genes were derepressed due to the absence of Su(Ste) piRNAs, and meiotic failures were observed. Thus, the piRNA pathway contributes to reproductive isolation between D. melanogaster and closely related species, causing hybrid male sterility via misregulation of two different host protein factors.

Natural variation in the courtship song of Drosophila is mapped to the intronic insertion of a retroelement at the slowpoke locus, which encodes an ion channel. An important aspect of courtship behaviour in Drosophila is the male courtship song, generated when the males vibrate their wings. The features of this courtship song and how this varies between Drosophila species have been well-characterized. David Stern and colleagues now map the genetic variation causal for natural variation in courtship song between two wild isolates of D. simulans and D. mauritiana to the insertion of a retroelement at the slowpoke (slo) locus, which encodes an ion channel. Animal species display enormous variation for innate behaviours, but little is known about how this diversity arose. Here, using an unbiased genetic approach, we map a courtship song difference between wild isolates of Drosophila simulans and Drosophila mauritiana to a 966 base pair region within the slowpoke (slo) locus, which encodes a calcium-activated potassium channel1. Using the reciprocal hemizygosity test2, we confirm that slo is the causal locus and resolve the causal mutation to the evolutionarily recent insertion of a retroelement in a slo intron within D. simulans. Targeted deletion of this retroelement reverts the song phenotype and alters slo splicing. Like many ion channel genes, slo is expressed widely in the nervous system and influences a variety of behaviours3,4; slo-null males sing little song with severely disrupted features. By contrast, the natural variant of slo alters a specific component of courtship song, illustrating that regulatory evolution of a highly pleiotropic ion channel gene can cause modular changes in behaviour.

Y chromosomes display population variation within and between species. Co-evolution within populations is expected to produce adaptive interactions between Y chromosomes and the rest of the genome. One consequence is that Y chromosomes from disparate populations could disrupt harmonious interactions between co-evolved genetic elements and result in reduced male fertility, sterility or inviability. Here we address the contribution of 'heterospecific Y chromosomes' to fertility in hybrid males carrying a homozygous region of Drosophila mauritiana introgressed in the Drosophila simulans background. In order to detect Y chromosome-autosome interactions, which may go unnoticed in a single-species background of autosomes, we constructed hybrid genotypes involving three sister species: Drosophila simulans, D. mauritiana, and D. sechellia. These engineered strains varied due to: (i) species origin of the Y chromosome (D. simulans or D. sechellia); (ii) location of the introgressed D. mauritiana segment on the D. simulans third chromosome, and (iii) grandparental genomic background (three genotypes of D. simulans). We find complex interactions between the species origin of the Y chromosome, the identity of the D. mauritiana segment and the grandparental genetic background donating the chromosomes. Unexpectedly, the interaction of the Y chromosome and one segment of D. mauritiana drastically reduced fertility in the presence of Ysim, whereas the fertility is partially rescued by the Y chromosome of D. sechellia when it descends from a specific grandparental genotype. The restoration of fertility occurs in spite of an autosomal and X-linked genome that is mostly of D. simulans origin. These results illustrate the multifactorial basis of genetic interactions involving the Y chromosome. Our study supports the hypothesis that the Y chromosome can contribute significantly to the evolution of reproductive isolation and highlights the conditional manifestation of infertility in specific genotypic combinations.

Male genital structures are among the most rapidly evolving morphological traits and are often the only features that can distinguish closely related species. This process is thought to be driven by sexual selection and may reinforce species separation. However, while the genetic basis of many phenotypic differences have been identified, we still lack knowledge about the genes underlying evolutionary differences in male genital organs and organ size more generally. The claspers (surstyli) are periphallic structures that play an important role in copulation in insects. Here we show that natural variation in clasper size and bristle number betweenDrosophila mauritianaandD. simulansis caused by evolutionary changes intartan (trn), which encodes a transmembrane leucine-rich repeat domain protein that mediates cell-cell interactions and affinity differences. There are no fixed amino acid differences intrnbetweenD. mauritianaandD. simulansbut differences in the expression of this gene in developing genitalia suggest cis-regulatory changes intrnunderlie the evolution of clasper morphology in these species. Finally, analysis of reciprocal hemizyotes that are genetically identical, except for which species the functional allele oftrnis from, determined that thetrnallele ofD. mauritianaspecifies larger claspers with more bristles than the allele ofD. simulans. Therefore we have identified the first gene underlying evolutionary change in the size of a male genital organ, which will help to better understand the rapid diversification of these structures and the regulation and evolution of organ size more broadly.Significance StatementThe morphology of male genital organs evolves rapidly driven by sexual selection. However, little is known about the genes underlying genitalia differences between species. Identifying these genes is key to understanding how sexual selection acts on development to produce rapid phenotypic change. We have found that the genetartanunderlies differences between maleDrosophila mauritianaandD. simulansin the size and bristle number of the claspers - genital projections that grasp the female during copulation. Moreover, sincetartanencodes a protein that is involved in cell affinity, this may represent a new developmental mechanism for morphological change. Therefore, our study provides new insights into genetic and developmental bases for the rapid evolution of male genitalia and organ size more generally.

Haldane’s rule has been the basis of speciation research during the last 30 years. Most studies have focused on the nature of incompatibilities in the hybrid male, but not much attention has been given to the genetic basis of fertility and inviability in hybrid females. Hybridizations between Drosophila simulans and Drosophila mauritiana produce fertile females and sterile males. Here, we re-examined the level of fertility in reciprocal F1 females of these two species and looked for the presence of maternal effects. Our results show that the reciprocal F1 females of D. simulans and D. mauritiana hybridizations are fully fertile and in fact show a significant level of heterosis in the rate of oviposition but display reduced egg hatching in one direction. Reduced egg hatching was observed in the progenies of F1 hybrid females with D. mauritiana as mother, the same cross that showed a stronger negative effect on F1 male fertility. A review of the literature on the hybridizations in Lepidoptera also showed a maternal effect on inviability when reciprocal crosses produced asymmetric results. Our findings point to the importance of maternal effects in the evolution of embryo inviability and thus enhancing the process of speciation through the evolution of hybrid inviability.

Crosses between closely related animal species often result in male hybrids that are sterile, and the molecular and functional basis of genetic factors for hybrid male sterility is of great interest. Here, we report a molecular and functional analysis of HMS1, a region of 9.2 kb in chromosome 3 of Drosophila mauritiana, which results in virtually complete hybrid male sterility when homozygous in the genetic background of sibling species Drosophila simulans. The HMS1 region contains two strong candidate genes for the genetic incompatibility, agt and Taf1. Both encode unrelated DNA-binding proteins, agt for an alkyl-cysteine-S-alkyltransferase and Taf1 for a subunit of transcription factor TFIID that serves as a multifunctional transcriptional regulator. The contribution of each gene to hybrid male sterility was assessed by means of germ-line transformation, with constructs containing complete agt and Taf1 genomic sequences as well as various chimeric constructs. Both agt and Taf1 contribute about equally to HMS1 hybrid male sterility. Transgenes containing either locus rescue sterility in about one-half of the males, and among fertile males the number of offspring is in the normal range. This finding suggests compensatory proliferation of the rescued, nondysfunctional germ cells. Results with chimeric transgenes imply that the hybrid incompatibilities result from interactions among nucleotide differences residing along both agt and Taf1. Our results challenge a number of preliminary generalizations about the molecular and functional basis of hybrid male sterility, and strongly reinforce the role of DNA-binding proteins as a class of genes contributing to the maintenance of postzygotic reproductive isolation.

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.

Odysseus (OdsH) has been identified as a hybrid male sterility gene between Drosophila mauritiana and D. simulans with accelerated evolutionary rate in both expression and DNA sequence. Loss of a testis-specific expression of Od sH causes male fertility defect in D. melanogaster. Yet, the underlying mechanisms at t he cellular level are unknown. In an attempt to identify the possible mechanisms and functional roles of OdsH i n spermatogenesis, the cell numbers at different developmental stages during spermatogenesis between the OdsH nul l mutant and wild-type flies were compared. The results showed that the early developing g erm cells, including spermatogonia and spermatocytes, were reduced in the OdsH mutant males. In addition, the number of germline stem cells in aged males was a lso reduced, presumably due to the disruption of germline stem cell maintenance, which resulted in more severe fertility defect. These results suggest that the f unction of the enhancement of sperm production by OdsH acted across males of all ages.

Frost (Fst) is a gene associated with cold exposure in Drosophila melanogaster. We used real-time PCR to assess whether cold exposure induces expression of Fst in 10 different life stages of D. melanogaster, and adults of seven other Drosophila species. We exposed groups of individuals to 0 °C (2 h), followed by 1 h recovery (22 °C). Frost was significantly upregulated in response to cold in eggs, third instar larvae, and 2- and 5-day-old male and female adults in D. melanogaster. Life stages in which cold did not upregulate Fst had high constitutive expression. Frost is located on the opposite strand of an intron of Diuretic hormone (DH), but cold exposure did not upregulate DH. Frost orthologues were identified in six other species within the Melanogaster group (Drosophila sechellia, Drosophila simulans, Drosophila yakuba, Drosophila erecta, Drosophila ananassae and Drosophila mauritiana). Frost orthologues were upregulated in response to cold exposure in both sexes in adults of all of these species. The predicted structure of a putative Frost consensus protein shows highly conserved tandem repeats of motifs involved in cell signalling (PEST and TRAF2), suggesting that Fst might encode an adaptor protein involved in acute stress or apoptosis signalling in vivo.

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100 000 years) compared to the origin of the species, around 0.25 MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.

The study of speciation has advanced considerably in the last decades because of the increased application of molecular tools. In particular, the quantification of gene flow between recently diverged species could be addressed. Drosophila simulans and Drosophila mauritiana diverged, probably allopatrically, from a common ancestor approximately 250 000 years ago. However, these species share one mitochondrial DNA (mtDNA) haplotype indicative of a recent episode of introgression. To study the extent of gene flow between these species, we took advantage of a large sample of D. mauritiana and employed a range of different markers, i.e. nuclear and mitochondrial sequences, and microsatellites. This allowed us to detect two new mtDNA haplotypes (MAU3 and MAU4). These haplotypes diverged quite recently from haplotypes of the siII group present in cosmopolitan populations of D. simulans. The mean divergence time of the most diverged haplotype (MAU4) is approximately 127 000 years, which is more than 100 000 years before the assumed speciation time. Interestingly, we also found some evidence for gene flow at the nuclear level because an excess of putatively neutral loci shows significantly reduced differentiation between D. simulans and D. mauritiana. Our results suggest that these species are exchanging genes more frequently than previously thought.

The Dobzhansky–Muller model posits that intrinsic postzygotic reproductive isolation—the sterility or lethality of species hybrids—results from the evolution of incompatible epistatic interactions between species: favorable or neutral alleles that become fixed in the genetic background of one species can cause sterility or lethality in the genetic background of another species. The kind of hybrid incompatibility that evolves between two species, however, depends on the particular evolutionary history of the causative substitutions. An allele that is functionally derived in one species can be incompatible with an allele that is functionally derived in the other species (a derived-derived hybrid incompatibility). But an allele that is functionally derived in one species can also be incompatible with an allele that has retained the ancestral state in the other species (a derived-ancestral hybrid incompatibility). The relative abundance of such derived-derived vs. derived-ancestral hybrid incompatibilities is unknown. Here, we characterize the genetics and evolutionary history of a lethal hybrid incompatibility between Drosophila mauritiana and its two sibling species, D. sechellia and D. simulans. We show that a hybrid lethality factor(s) in the pericentric heterochromatin of the D. mauritiana X chromosome, hybrid lethal on the X (hlx), is incompatible with a factor(s) in the same small autosomal region from both D. sechellia and D. simulans, Suppressor of hlx [Su(hlx)]. By combining genetic and phylogenetic information, we infer that hlx-Su(hlx) hybrid lethality is likely caused by a derived-ancestral incompatibility, a hypothesis that can be tested directly when the genes are identified.

Drosophila focal adhesion kinase (Dfak) gene is a single-copy nuclear gene. Previous study revealed that Drosophila melanogaster and Drosophila simulans had lost an intron precisely within the tyrosine kinase (TyK) domain of this gene. However, this did not happen in several other Drosophila species, including Drosophila elegans, Drosophila ficusphila, Drosophila biarmipes, Drosophila jambulina, Drosophila prostipennis, Drosophila takahashii, and Drosophila pseudoobscura. In the current study, homologous sequences of Drosophila sechellia, Drosophila mauritiana, Drosophila yakuba, Drosophila teissieri, Drosophila santomea, and Drosophila erecta were amplified by polymerase chain reaction, and further sequencing analysis indicated that these species were missing a TyK domain intron, indicating they were closely related. The relationship of the D. melanogaster species group was reconstructed using TyK domain nucleotide sequences. The resulting phylogenetic tree revealed that these 8 species were the most related species in the melanogaster group. These results strongly support previously proposed classifications based on morphological and molecular data.

SUMMARY 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.

An intact mariner-like element was isolated from the genome of tobacco budworm, Heliothis virescens. This is the first report of an intact mariner element after the initial identification of Mos1 from Drosophila mauritiana. The full-length Hvmar1 has 30 bp inverted terminal repeats and a complete 1065 bp open reading frame (ORF) encoding the transposase with a 'D,D(34)D' motif in the catalytic domain. Polymerase chain reaction results show that at least one insertion of the Hvmar1 element is conserved in this Heliothis strain. Phylogenetic analysis demonstrated that Hvmar1 belongs to the irritans subfamily.

Species often produce sterile hybrids early in their evolutionary divergence, and some evidence suggests that hybrid sterility may be associated with deviations or disruptions in gene expression. In support of this idea, many studies have shown that a high proportion of male-biased genes are underexpressed, compared with non-sex-biased genes, in sterile F1 male hybrids of Drosophila species. In this study, we examined and compared patterns of misexpression in sterile F1 male hybrids of Drosophila simulans and 2 of its sibling species, Drosophila mauritiana and Drosophila sechellia, at both the larval and adult life stages. We analyzed hybrids using both commercial Drosophila melanogaster microarrays and arrays we developed from reverse transcriptase-polymerase chain reactions of spermatogenesis and reproduction-related transcripts from these species (sperm array). Although the majority of misexpressed transcripts were underexpressed, a disproportionate number of the overexpressed transcripts were located on the X chromosome. We detected a high overlap in the genes misexpressed between the 2 species pairs, and our sperm array was better at detecting such misexpression than the D. melanogaster array, suggesting possible weaknesses in the use of an array designed from another species. We found only minimal misexpression in the larval samples with the sperm array, suggesting that disruptions in spermatogenesis occur after this life stage. Further study of these misexpressed loci may allow us to identify precisely where disruptions in the spermatogenesis pathway occur.

We determined the complete nucleotide sequence of the A+T-rich region of the maII type of mtDNA in D. mauritiana. The nucleotide sequence was found to contain 3,206 bp. Three types of conserved element, i.e., type I element, type II element, and T-stretch, were included in this sequence, as reported for D. melanogaster. Comparison between the two species revealed that the type I elements were less conserved than the type II elements. However, each of these type I elements contained a G-stretch within a loop of a putative stem-loop-forming sequence, which has also been observed in D. melanogaster. Moreover, in both type I and type II repeat arrays, the elements closest to the T-stretch diverged the most, due to nucleotide substitution and/or the insertion of short repeats. Sequence comparison of the two complete sequences of the A+T-rich region of D. melanogaster and the maII type of D. mauritiana, as well as comparison of partial sequences in other types of mtDNA within the melanogaster complex, suggested that the A+T-rich region in this complex has been maintained by concerted evolution after the duplication of two types of element, i.e., type I and type II.

Crosses of Drosophila melanogaster females to males of its sibling species Drosophila simulans, Drosophila mauritiana and Drosophila sechellia produce no sons and daughters that are viable only at low temperatures. We describe here a novel rescue allele Df(1)EP307-1-2 isolated on the basis of its suppression of high temperature hybrid female lethality. Df(1)EP307-1-2 also rescues hybrid males to the pharate adult stage, the same stage at which it is lethal to D. melanogaster pure species males. Molecular analysis indicates that Df(1)EP307-1-2 is associated with a deletion of about 61 kb in the 9D region of the X chromosome. The structure of Df(1)EP307-1-2 suggests that it was formed by a process similar to P-element induced male recombination.

The genetic basis of hybrid incompatibility in crosses between Drosophila mauritiana and D. simulans was investigated to gain insight into the evolutionary mechanisms of speciation. In this study, segments of the D. mauritiana third chromosome were introgressed into a D. simulans genetic background and tested as homozygotes for viability, male fertility, and female fertility. The entire third chromosome was covered with partially overlapping segments. Many segments were male sterile, while none were female sterile or lethal, confirming previous reports of the rapid evolution of hybrid male sterility (HMS). A statistical model was developed to quantify the HMS accumulation. In comparison with previous work on the X chromosome, we estimate that the X has ∼2.5 times the density of HMS factors as the autosomes. We also estimate that the whole genome contains ∼15 HMS “equivalents”—i.e., 15 times the minimum number of incompatibility factors necessary to cause complete sterility. Although some caveats for the quantitative estimate of a 2.5-fold density difference are described, this study supports the notion that the X chromosome plays a special role in the evolution of reproductive isolation. Possible mechanisms of a “large X” effect include selective fixation of new mutations that are recessive or partially recessive and the evolution of sex-ratio distortion systems.

We compared deleted copies of the seven mauritiana subfamilies of mariner transposable elements in species of the Drosophilidae. All elements were detected by PCR using the inverted terminal repeats of the Mos1 element of Drosophila mauritiana as primers. A higher frequency of breakpoints in the 5prime prime or minute part of the element compared to the 3prime prime or minute part was observed. Of the 27 deletions, 9 (33%) occurred between short direct repeats (SDR) of 5 to 8 bp. The SDRs can be at or close to the breakpoints of the deletion. A deleted copy of D. simulans (St. Martin population) had three repeats of a motif present only once in the complete consensus sequence. The high frequency of SDRs at or near the breakpoints of the deletions strongly suggests that some of them do not occur at random. Mechanisms that might explain these deletions, such as unequal crossing-over, ectopic recombination, and abortive gap repair, are discussed.

Male sexual characters are often among the first traits to diverge between closely related species and identifying the genetic basis of such changes can contribute to our understanding of their evolutionary history. However, little is known about the genetic architecture or the specific genes underlying the evolution of male genitalia. The morphology of the claspers, posterior lobes, and anal plates exhibit striking differences between Drosophila mauritiana and D. simulans. Using QTL and introgression-based high-resolution mapping, we identified several small regions on chromosome arms 3L and 3R that contribute to differences in these traits. However, we found that the loci underlying the evolution of clasper differences between these two species are independent from those that contribute to posterior lobe and anal plate divergence. Furthermore, while most of the loci affect each trait in the same direction and act additively, we also found evidence for epistasis between loci for clasper bristle number. In addition, we conducted an RNAi screen in D. melanogaster to investigate if positional and expression candidate genes located on chromosome 3L, are also involved in genital development. We found that six of these genes, including components of Wnt signaling and male-specific lethal 3 (msl3), regulate the development of genital traits consistent with the effects of the introgressed regions where they are located and that thus represent promising candidate genes for the evolution these traits.

The ability of transposable elements to mobilize across genomes and affect the expression of genes makes them exceptional tools for genetic manipulation methodologies. Several transposon-based systems have been modified and incorporated into shuttle mutagenesis approaches in a variety of organisms. We have found that the Mos1 element, a DNA transposon from Drosophila mauritiana, is suitable and readily adaptable to a variety of strategies to the study of trypanosomatid parasitic protozoa. Trypanosomatids are the causative agents of a wide range of neglected diseases in underdeveloped regions of the globe. In this chapter we describe the basic elements and the available protocols for the in vitro use of Mos1 derivatives in the protozoan parasite Leishmania.

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is ∼50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.

Molecular differentiation between races or closely related species is often incongruent with the reproductive divergence of the taxa of interest. Shared ancient polymorphism and/or introgression during secondary contact may be responsible for the incongruence. At loci contributing to speciation, these two complications should be minimized (1, 2); hence, their variation may more faithfully reflect the history of the species' reproductive differentiation. In this study, we analyzed DNA polymorphism at the Odysseus ( OdsH ) locus of hybrid sterility between Drosophila mauritiana and Drosophila simulans and were able to verify such a prediction. Interestingly, DNA variation only a short distance away (1.8 kb) appears not to be influenced by the forces that shape the recent evolution of the OdsH coding region. This locus thus may represent a test case of inferring phylogeny of very closely related species.

The sexual isolation among the related species Drosophila melanogaster, D. simulans and D. mauritiana is asymmetrical. While D. mauritiana males mate well with both D. melanogaster and D. simulans females, females of D. mauritiana discriminate strongly against males of these two species. Similarly, D. simulans males mate with D. melanogaster females but the reciprocal cross is difficult. Interspecific crosses between several populations of the three species were performed to determine if (i) males and females of the same species share a common sexual isolation genetic system, and (ii) males (or females) use the same genetic system to discriminate against females (or males) of the other two species. Results indicate that although differences in male and female isolation depend on the populations tested, the isolation behaviour between a pair of species is highly correlated despite the variations. However, the rank order of the isolation level along the populations was not correlated in both sexes, which suggests that different genes act in male and female sexual isolation. Neither for males nor for females, the isolation behaviour of one species was paralleled in the other two species, which indicates that the genetic systems involved in this trait are species-pair specific. The implications of these results are discussed.

Fully degenerate polymersae chain reaction (PCR) primers, designed to represent amino-acid conserved regions between transposase genes of the Drosophila mauritiana Tsacas & David and Hyalophora cecropia (L.) mariner transposable elements, were used to screen genomic DNA of 17 insect cell lines for a possible mariner -like element. All cell lines except that of Aedes albopictus (Skuse) were positive, yielding a single appropriately sized band of 490 bp. This PCR product from 2 of the cell lines, Heliothis subflexa (Guenee) and Heliothis virescens (F.), was further characterized by cloning and direct sequencing of PCR-amplified DNA products to confirm whether the product was a mariner -like transposon element. Alignment of the direct sequenced PCR-amplified DNA mariner -like element showed an amino acid identity of 43% to the mariner transposase gene of the Surinam cockroach, Pycnocelus surinamensis (L.).