Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

Microsatellite variation suggests a recent fine-scale population structure of <Emphasis Type="Italic">Drosophila sechellia</Emphasis>, a species endemic of the Seychelles archipelago

Drosophila sechellia is closely related to the cosmopolitan and widespread model species, D. simulans. This species, endemic to the Seychelles archipelago, is specialized on the fruits of Morinda citrifolia, and harbours the lowest overall genetic diversity compared to other species of Drosophila. This low diversity is associated with a small population size. In addition, no obvious population structure has been evidenced so far across islands of the Seychelles archipelago. Here, a microsatellite panel of 17 loci in ten populations from nine islands of the Seychelles was used to assess the effect of the D. sechellia’s fragmented distribution on the fine-scale population genetic structure, the migration pattern, as well as on the demography of the species. Contrary to previous results, also based on microsatellites, no evidence for population contraction in D. sechellia was found. The results confirm previous studies based on gene sequence polymorphism that showed a long-term stable population size for this species. Interestingly, a pattern of Isolation By Distance which had not been described yet in D. sechellia was found, with evidence of first-generation migrants between some neighbouring islands. Bayesian structuring algorithm results were consistent with a split of D. sechellia into two main groups of populations: Silhouette/Mahé versus all the other islands. Thus, microsatellites suggest that variability in D. sechellia is most likely explained by local genetic exchanges between neighbouring islands that have recently resulted in slight differentiation of the two largest island populations from all the others.     

Hybrid disadvantage in the larval foraging behaviour of the two neotropical species of <Emphasis Type="Italic">Drosophila pavani</Emphasis> and <Emphasis Type="Italic">Drosophila gaucha</Emphasis>

Flies from two populations of the Chilean endemic neotropical species Drosophila pavani and two populations of its sibling species Drosophila gaucha were crossed reciprocally to obtain intra- and interspecific hybrids. The developmental pathways of locomotor activity and feeding rate were analysed for eleven of twelve possible genotype groups. The hybrids showed reduced fitness indicated by a decrease in the measured traits. Hybrid disadvantage was strongest in interspecific hybrids, especially with respect to feeding behaviour. This evidence supports the contention that D. pavani and D. gaucha have evolved different coadapted gene pools controlling the developmental pathways for behavioural traits expressed during larval foraging; but genetic divergence affecting these behaviours has also taken place between locally adapted populations within each species.     

Association of Misexpression with Sterility in Hybrids of <Emphasis Type="Italic">Drosophila simulans</Emphasis>and <Emphasis Type="Italic">D. mauritiana</Emphasis>

Recent studies have identified genes associated with hybrid sterility and other hybrid dysfunctions, but the consequences of introgressions of these speciation genes are often poorly understood. Previously, we identified a panel of genes that are underexpressed in sterile male hybrids of Drosophila simulans and D. mauritiana relative to pure species. Here, we build on this reverse-genetics approach to demonstrate that the underexpression of at least five of these genes in hybrids is associated with hybrid sterility and that these five genes are coordinately regulated. We map one upstream regulator of these genes to a region previously shown to harbor one or more factors causing hybrid sterility. Finally, we show that the genes underexpressed in hybrids are often highly conserved, as might be predicted for downstream targets of the genetic changes that cause hybrid sterility. This approach integrates forward genetics with reverse genetics to show a proximate consequence of the introgression of particular hybrid sterility-conferring regions between species: underexpression of genes necessary for normal spermatogenesis.[Reviewing Editor: Martin Kreitman]     

Nucleotide variation at the dopa decarboxylase (<Emphasis Type="Italic">Ddc</Emphasis>) gene in natural populations of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

We studied nucleotide sequence variation at the gene coding for dopa decarboxylase (Ddc) in seven populations of Drosophila melanogaster. Strength and pattern of linkage disequilibrium are somewhat distinct in the extensively sampled Spanish and Raleigh populations. In the Spanish population, a few sites are in strong positive association, whereas a large number of sites in the Raleigh population are associated nonrandomly but the association is not strong. Linkage disequilibrium analysis shows presence of two groups of haplotypes in the populations, each of which is fairly diverged, suggesting epistasis or inversion polymorphism. There is evidence of two forms of natural selection acting on Ddc. The McDonald-Kreitman test indicates a deficit of fixed amino acid differences between D. melanogaster and D. simulans, which may be due to negative selection. An excess of derived alleles at high frequency, significant according to the H-test, is consistent with the effect of hitchhiking. The hitchhiking may have been caused by directional selection downstream of the locus studied, as suggested by a gradual decrease of the polymorphism-to-divergence ratio. Altogether, the Ddc locus exhibits a complicated pattern of variation apparently due to several evolutionary forces. Such a complex pattern may be a result of an unusually high density of functionally important genes.     

Quantitative genetics of speciation: additive and non-additive genetic differentiation between <Emphasis Type="Italic">Drosophila madeirensis</Emphasis> and <Emphasis Type="Italic">Drosophila subobscura</Emphasis>

The role of dominance and epistasis in population divergence has been an issue of much debate ever since the neoDarwinian synthesis. One of the best ways to dissect the several genetic components affecting the genetic architecture of populations is line cross analysis. Here we present a study comparing generation means of several life history-traits in two closely related Drosophila species: Drosophila subobscura, D. madeirensis as well as their F ₁ and F ₂ hybrids. This study aims to determine the relative contributions of additive and non-additive genetic parameters to the differentiation of life-history traits between these two species. The results indicate that both negative dominance and epistatic effects are very important in the differentiation of most traits. We end with considerations about the relevance of these findings for the understanding of the role of non-additive effects in speciation.     

Fine-scale genetic mapping of a hybrid sterility factor between <Emphasis Type="Italic">Drosophila simulans</Emphasis> and <Emphasis Type="Italic">D. mauritiana</Emphasis>: the varied and elusive functions of "speciation genes"

Background  

Hybrid male sterility (HMS) is a usual outcome of hybridization between closely related animal species. It arises because interactions between alleles that are functional within one species may be disrupted in hybrids. The identification of genes leading to hybrid sterility is of great interest for understanding the evolutionary process of speciation. In the current work we used marked P-element insertions as dominant markers to efficiently locate one genetic factor causing a severe reduction in fertility in hybrid males of Drosophila simulans and D. mauritiana.     

Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

Molecular Phylogeny of the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> Species Subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba–D. teissieri and D. erecta–D. orena species pairs relative to the D. melanogaster–D. simulans–D. mauritiana–D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed vectorette PCR to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba–D. teissieri and D. erecta–D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica–D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models. Present address (Ryan M. David): University of Texas Health Science Center at San Antonio, 7703 Floyd Curl Drive, San Antonio, TX 78284, USA     

High genetic differentiation between an African and a non-African strain of <Emphasis Type="Italic">Drosophila simulans</Emphasis> revealed by segregation distortion and reduced crossover frequency

Drosophila simulans strains originating from Madagascar and nearby islands in the Indian Ocean often differ from those elsewhere in the number of sex comb teeth and the degree of morphological anomaly in hybrids with D. melanogaster. Here, we report a strong segregation distortion in the F1 intercross between two D. simulans strains originating from Madagascar and the US, possibly at both the gametic and zygotic levels. Strong bias against alleles of the Madagascar strain was observed for all ten marker loci distributed over the entire second chromosome in the F1 intercross, but only a few showed a weak distortion in the isogenic backgrounds of either strains. Significant deviations of genotype frequencies from Hardy–Weinberg proportions were consistently observed for the second chromosome. By contrast, the X and third chromosomes did not show any strong segregation distortion. Crossover frequency on the second chromosome was uniformly reduced in isogenic backgrounds whereas the map lengths in the F1 intercross were comparable to or larger than that of the standard D. melanogaster map. We discuss these findings in relation to previous studies on other traits and interspecific differences between D. mauritiana, which is endemic to Mauritius Island, and D. simulans.     

Population dynamics of a maternally-transmitted <Emphasis Type="Italic">Spiroplasma</Emphasis> infection in <Emphasis Type="Italic">Drosophila hydei</Emphasis>

A maternally-inherited spiroplasma endosymbiont of Drosophila hydei does not exert apparent phenotypes on both sexes of its host and is prevalent in natural populations of D. hydei. Our previous experiments using a laboratory stock of D. hydei revealed that low temperatures (such as 15°C and 18°C) dramatically lower the vertical transmission rates of this spiroplasma. Therefore, we hypothesized that, in temperate regions, the infection frequencies may decrease in cool seasons but increase in the summer season. To clarify the temporal population dynamics of the spiroplasma infection, D. hydei were collected from two Japanese populations in 2006–2008 from May to early August, representing the only period when a number of D. hydei are collectable in Japan, and examined for spiroplasma infection. Within each year, the frequency of spiroplasma infection fluctuated considerably in both populations. Consistent with our hypothesis, the infection frequency showed an increasing trend in both populations in 2007. However, the data in 2006 and 2008 did not show consistent patterns of increase. The population dynamics of spiroplasma infection may be affected but not critically determined by temperature. Moreover, despite the fluctuation within each year, the infection frequencies seemed to be stable across the years. The frequencies of spiroplasma infection in D. hydei populations may be stabilized by multiple factors. One of these factors may involve a context-dependent positive effect of spiroplasma on the fitness of D. hydei, as was recently observed in laboratory experiments.     

Similar rates of protein adaptation in <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. melanogaster</Emphasis>, two species with different current effective population sizes

Background  

Adaptive protein evolution is common in several Drosophila species investigated. Some studies point to very weak selection operating on amino-acid mutations, with average selection intensities on the order of N _e s ~ 5 in D. melanogaster and D. simulans. Species with lower effective population sizes should undergo less adaptation since they generate fewer mutations and selection is ineffective on a greater proportion of beneficial mutations.     

Natural hybridization of two mussel species <Emphasis Type="Italic">Dreissena polymorpha</Emphasis> (Pallas, 1771) and <Emphasis Type="Italic">Dreissena bugensis</Emphasis> (Andrusov, 1897)

Dreissenids display a high diversity of shell morphology, and it is frequently difficult to ascribe some individuals from mixed populations to one of the two species, Dreissena polymorpha (Pallas, 1771) or D. bugensis (Andrusov, 1897). Presumably, such individuals may be interspecific hybrids. We have analyzed species-specific allozyme loci of the typical representatives of these two mussel species and putative interspecific hybrids. A natural interspecific hybrid between D. polymorpha and D. bugensis was discovered for the first time by genetic methods. It has been demonstrated that D. bugensis was a maternal parent.     

Male accessory gland secretory protein polymorphism in natural populations of <Emphasis Type="Italic">Drosophila nasuta nasuta</Emphasis> and <Emphasis Type="Italic">Drosophila sulfurigaster neonasuta</Emphasis>

Male accessory gland secretory protein polymorphism was analysed in natural populations of Drosophila nasuta nasuta and D. sulfurigaster neonasuta for the first time, using SDS-PAGE to score polymorphism of these proteins in 2788 individuals of D. n. nasuta and 2232 individuals of D. s. neonasuta from 12 different populations from southern India. A total of 25 and 18 variant protein phenotypes were identified in D. n. nasuta and D. s. neonasuta, respectively. Protein fractions of group III were more polymorphic than those from groups I and II. The results show that accessory gland secretory proteins show high levels of polymorphism, irrespective of species or habitat. Moreover, we have used the variation in the accessory gland proteins to assess the extent of divergence between the species and to infer their population structure. The study suggests that though both D. n. nasuta and D. s. neonasuta belong to the same subgroup, they differ in population structure, as far as accessory gland protein polymorphism is concerned.     

Distinct electrophoretic polymorphism pattern at alcohol dehydrogenase (<Emphasis Type="Italic">Adh</Emphasis>) locus of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> natural populations from Turkey

Small number of Drosophila melanogaster populations from two distinct geographical regions of Turkey, Central Anatolia and Black Sea, were studied. Populations sampled were electrophoresed for a single locus, alcohol dehydrogenase (Adh) to assess population differentiation. Both the magnitude of genetic differentiation levels and the population structure based on hierarchical F-statistics allow populations to be grouped on two genetically divergent area, Central Anatolian and Black Sea. One ecological correlate, average yearly maximum rainfall, R _year, seems to track this Adh genetic variation pattern. The study also shows that a typical pattern of geographical Adh polymorphism can emerge with a handfull of populations sampled across a relatively small distance. The text was submitted by the author in English.     

Genetic differentiation in sympatric wood ants, <Emphasis Type="Italic">Formica rufa</Emphasis> and <Emphasis Type="Italic">F. polyctena</Emphasis>

Summary Direct observations have suggested that the closely related wood ants Formica polyctena and F. rufa represent different social organizations, with high queen number in F. polyctena and a high frequency of monogynous nests in F. rufa. We examined social organization and genetic population structure in a setup where populations of the two species are sympatric and gene flow between the species is possible. Our aim was to compare social organization in the species, and study evolutionary relationships between them. The observed relatedness among colony workers suggested that the difference in the level of polygyny is quantitative rather than qualitative, with a higher queen number in F. polyctena. The observed difference in polygyny was not accompanied by a difference in spatial genetic differentiation which was weak in both species. The genetic distance between the species is consistent with limited interspecific gene flow. Identification of a few possible F. rufa migrants in F. polyctena populations suggests potential interspecific gene flow. Thus, reproductive isolation of the species may not be complete when they are sympatric.Received 14 March 2003; revised 10 October 2003; accepted 20 October 2003.     

Localization of <Emphasis Type="Italic">Dras1</Emphasis> gene in polytene chromosomes of sibling species of <Emphasis Type="Italic">Drosophila virilis</Emphasis> group

The Dras1 gene was mapped by in situ hybridization to polytene chromosomes of several sibling species of the Drosophila virilis group and their hybrids. A 1037-bp fragment of Dras1 gene from the D. virilis genome was used as the probe. The gene sequence was localized in the region of a 25 A-B disk in chromosome 2 (in accordance with the D. virilis polytene chromosome map (Gubenko and Evgen’ev, 1984).     

A New Retroposed Gene in <Emphasis Type="Italic">Drosophila</Emphasis> Heterochromatin Detected by Microarray-Based Comparative Genomic Hybridization

A genomic pattern of new gene origination is often dependent on a genomic method that can efficiently identify a statistically adequate number of recently originated genes. The heterochromatic regions have often been viewed as genomic deserts with low coding potential and thus a low flux of new genes. However, increasing reports revealed unexpected roles of heterochromatic regions in the evolution of genes and genomes. We identified recently retroposed genes that originated in heterochromatic regions in Drosophila, by developing microarray-based comparative genomic hybridization (CGH) with multiple species. This new gene family, named Ifc-2h, originated in the common ancestor of the clade of D. simulans, D. mauritiana, and D. sechellia. The sequence features and phylogenetic distribution indicated that Ifc-2h resulted from the retroposition from its parental gene, Infertile crescent (Ifc), and integrated into heterochromatic region of common ancester of the three sibling species 2 million years ago. Expression analysis revealed that Ifc-2h had developed a new expression pattern by recruiting a putative regulatory element from its target sequence. The distribution of indel variation in Ifc-2h of D. simulans and D. mauritiana revealed a significant sequence constraint, suggesting that the Ifc-2h gene may be functional. These analyses cast fresh insight into the evolution of heterochromatin and the origin of its coding regions. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

Molecular characterization of the <Emphasis Type="Italic">singed wings</Emphasis> locus of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Background  

Hormones frequently guide animal development via the induction of cascades of gene activities, whose products further amplify an initial hormonal stimulus. In Drosophila the transformation of the larva into the pupa and the subsequent metamorphosis to the adult stage is triggered by changes in the titer of the steroid hormone 20-hydroxyecdysone. singed wings (swi) is the only gene known in Drosophila melanogaster for which mutations specifically interrupt the transmission of the regulatory signal from early to late ecdysone inducible genes.     

Post-pollination hybridization barriers in <Emphasis Type="Italic">Nicotiana</Emphasis> section <Emphasis Type="Italic">Alatae</Emphasis>

Nicotiana section Alatae contains eight species with variable flower sizes and morphologies. Section members readily hybridize in the glasshouse, but no hybrids have been observed in natural sympatric and parapatric populations. To investigate interspecific crossing relationships with respect to mechanisms preventing hybridization, all members of section Alatae were intercrossed in a complete diallel. We found positive correlation between the pistil length of the pollen donor and interspecific seed set relative to the conspecific cross. Pollen tube growth rate and pollen donor pistil length were positively correlated as well. Furthermore, pollen from short-pistil members of section Alatae could only grow a maximum distance proportional to, but greater than, their own pistil lengths. Our results show that pollen tube growth capacity (i.e., rate and distance), provides a hybridization barrier in long-pistil species × short-pistil species crosses. We also found another hybridization barrier not specifically related to pollen tube growth capacity in short-pistil species × long-pistil species. Taken together, these barriers can generally be described by a ‘pistil-length mismatch’ rule; in section Alatae, pollen has the most success fertilizing ovules from species with pistil lengths similar to their own. This rule could contribute to hybridization barriers in Section Alatae because the species display dramatically different pistil lengths. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.