Geographic Distribution and Inheritance of Three Cytoplasmic Incompatibility Types in Drosophila Simulans

Wolbachia-like microorganisms have been implicated in unidirectional cytoplasmic incompatibility between strains of Drosophila simulans. Reduced egg eclosion occurs when females from uninfected strains (type W) are crossed with males from infected strains (type R). Here we characterize a third incompatibility type (type S) which is also correlated with the presence of Wolbachia-like microorganisms. Despite the fact that the symbionts cannot be morphologically distinguished, we observed complete bidirectional incompatibility between R and S strains. This indicates that the determinants of incompatibility are different in the two infected types. S/W incompatibility is unidirectional and similar to R/W incompatibility. A worldwide survey of D. simulans strains showed that type S incompatibility was found only in insular populations which harbor the mitochondrial type SiI. Both W and R types were found among mainland and island populations harboring the worldwide mitochondrial type SiII. Type S incompatibility could be involved in the reinforcement of the geographical isolation of SiI populations.     

Divergence of mitochondrial dna is not corroborated by nuclear dna, morphology, or behavior in Drosophila simulans

We ask whether the observed mitochondrial DNA (mtDNA) population subdivision of Drosophila simulans is indicative of organismal structure or of specific processes acting on the mitochondrial genome. Factors either intrinsic or extrinsic to the host genome may influence the evolutionary dynamics of mtDNA. Potential intrinsic factors include adaptation of the mitochondrial genome and of nucleomitochondrial gene complexes specific to the local environment. An extrinsic force that has been shown to influence mtDNA evolution in invertebrates is the bacterial endosymbiont Wolbachia. Evidence presented in this study suggests that mtDNA is not a good indicator of organismal subdivision in D. simulans. Furthermore, there is no evidence to suggest that Wolbachia causes any reduction in nuclear gene flow in this species. The observed differentiation in mtDNA is not corroborated by data from NADH: ubiquinone reductase 75kD subunit precursor or the Alcohol dehydrogenase-related loci, from the shape or size of the male genital arch, or from assortative premating behavior. We discuss these results in relation to a mitochondrial genetic species concept and the potential for Wolbachia-induced incompatibility to be a mechanism of speciation in insects. We conclude with an iterated appeal to include phylogenetic and statistical tests of neutrality as a supplement to phylogenetic and population genetic analyses when using mtDNA as an evolutionary marker.     

Variability within the Seychelles Cytoplasmic Incompatibility System in Drosophila Simulans

In Drosophila simulans, we described a cytoplasmic incompatibility (CI) system (Seychelles) restricted to insular populations that harbor the mitochondrial type SiI. Since then, these populations have been shown to be heterogeneous, some being infected by one Wolbachia genetic variant only (wHa), while others are infected simultaneously by wHa and by another variant (wNo) always found in association with wHa. We have experimentally obtained two D. simulans strains only infected by the wNo variant. This variant determines its own cytoplasmic incompatibility type. In particular, the cross between wNo-bearing flies and wHa-bearing ones is bidirectionally incompatible. The Seychelles CI type, stricto sensu, is distinguished by being determined by the simultaneous presence of two Wolbachia variants that we found to be mutually incompatible. In addition, we observed incomplete maternal transmission of the Wolbachia.     

The evolution of reproductive isolation in the ectomycorrhizal Hebeloma crustuliniforme aggregate (Basidiomycetes) in northwestern Europe: a phylogenetic approach

Abstract. To reconstruct the evolution of reproductive isolation in the ectomycorrhizal Hebeloma crustuliniforme aggregate (Basidiomycetes), phylogenetic relationships were determined between strains that belong to a clade consisting of nine intercompatibility groups (ICGs, biological species). Four of these nine ICGs are partially compatible and belong to the H. crustuliniforme aggregate. Different levels of partial compatibility have been found between these four ICGs. Between ICGs 3 and 4, 15% of the combinations were compatible. One strain was compatible with all isolates of both ICGs 3 and 4 and also with one isolate of ICG 2. Both a nuclear phylogeny, based on ribosomal IGS sequence data, and a mitochondrial phylogeny, based on a group-I intron located in the large subunit ribosomal RNA gene (LrRNA), were reconstructed. The level of incompatibility was compared with the phylogenetic history of individuals belonging to this clade.
Different relationships were found between the level of compatibility and the relative age of the most recent common ancestor (MRCA) for different ICGs. On the one hand, the evolution of incompatibility between ICGs 2 and 3/4 is most consistent with the class of "divergence- first" models because a positive correlation was found between the relative age of the MRCA and the level of incompatibility for ICG 2 versus 3/4. On the other hand, the lack of such a correlation for ICGs 3 and 4 shows that (partial) incompatibility between these ICGs has arisen without strong divergence. The ecological (and to a lesser extent geographical) differences found between ICGs 3 and 4 suggest that selection for incompatibility, associated with host tree preference, has been important in the evolution of incompatibility between these two ICGs. The incongruence between the nuclear and mitochondrial trees for ICG 1 could be explained by a hybrid origin of this ICG, with different donors of the mitochondrial and nuclear sequences.     

The phylogeny of the social wasp subfamily Polistinae: evidence from microsatellite flanking sequences,mitochondrial COI sequence,and morphological characters

Background  

Social wasps in the subfamily Polistinae (Hymenoptera: Vespidae) have been important in studies of the evolution of sociality, kin selection, and within colony conflicts of interest. These studies have generally been conducted within species, because a resolved phylogeny among species is lacking. We used nuclear DNA microsatellite flanking sequences, mitochondrial COI sequence, and morphological characters to generate a phylogeny for the Polistinae (Hymenoptera) using 69 species.     

Phylogeny, diversity, and classification of the Accipitridae based on DNA sequences of the RAG-1 exon

The avian family Accipitridae has historically been divided into subfamilies or tribes based on features such as general resemblance, feeding ecology, and behavior. Consequently, the monophyly of those groups has been questionable. Recently, three phylogenetic analyses of a majority of the genera have appeared, one based on osteology, one on DNA sequences from a single mitochondrial gene, and the third on mitochondrial plus nuclear DNA sequences, and the resulting phylogenies were in substantial disagreement concerning the composition and basal branching patterns of the clades and hence require further analysis and confirmation. Here we use DNA sequences from the large nuclear RAG-1 exon to investigate the phylogenetic relationships of these birds. Our results largely corroborated the prior study that included nuclear genes. We found strong support for a monophyletic clade comprising the secretarybird Sagittarius serpentarius , the osprey Pandion haliaetus , and the traditional accipitrids. However, every one of the traditionally recognized subfamilies of accipitrids was found to be polyphyletic. The most basal nodes in the phylogeny separate small clades of insectivorous and scavenger species, such as kites and Old World vultures, from the rest of the family. The speciose genera of bird and mammal predators are all relatively derived (terminal) in the phylogeny. Many of the basal clades are cosmopolitan in their distributions, consistent with the great mobility of these raptors. A new classification is proposed that eliminates the problem of polyphyletic intrafamilial taxa.     

Expression of cytoplasmic incompatibility in Drosophila simulans and its impact on infection frequencies and distribution of Wolbachia pipientis

The aim of this study is to examine the expression of cytoplasmic incompatibility and investigate the distribution and population frequencies of Wolbachia pipientis strains in Drosophila simulans. Nucleotide sequence data from 16S rDNA and a Wolbachia surface protein coding sequence and cytoplasmic incompatibility assays identify four distinct Wolbachia strains: wHa, wRi, wMa, and wAu. The levels of cytoplasmic incompatibility between six lines carrying these strains of bacteria and three control lines without bacteria are characterized. Flies infected with wHa and wRi are bidirectionally incompatible, and males that carry either strain can only successfully produce normal numbers of offspring with females carrying the same bacterial strain. Males infected with wAu do not express incompatibility. Males infected with the wMa strain express intermediate incompatibility when mated to females with no bacteria and no incompatibility with females with any other Wolbachia strain. We conduct polymerase chain reaction/restriction fragment length polymorphism assays to distinguish the strain of Wolbachia and the mitochondrial haplotype to survey populations for each type and associations between them. Drosophila simulans is known to have three major mitochondrial haplotypes (siI, sill, and siIII) and two subtypes (siIIA and siIIB). All infected lines of the sil haplotype carry wHa, wNo, or both; wMa and wNo are closely related and it is not clear whether they are distinct strains or variants of the same strain. Infected lines with the silIA haplotype harbor wRi and the siIIB haplotype carries wAu. The wMa infection is found in siIII haplotype lines. The phenotypic expression of cytoplasmic incompatibility and its relation to between-population differences in frequencies of Wolbachia infection are discussed.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.     

Genetic Control of Mitochondria in Paramecium

Mitochondrial mutations for resistance to various antibiotics (erythromycin, chloramphenicol, spiramycin, mikamycin) have been obtained in Paramecium aurelia and their properties are reviewed. Using these mitochondrial markers, the interactions between nucleus and mitochondria have been studied in two ways: by microinjection of mitochondria from one stock or species into other stocks and species of P. aurelia and by a genetic study of a nuclear mutation affecting mitochondrial multiplication. Both types of experiments show: (1) that there may exist incompatibility between a given type of mitochondria and the cell into which they are introduced and (2) that through multiplication in the host cell, mitochondrial properties can be modified. The possible basis for incompatibility and host-induced modifications is discussed.     

Linking phylogenetics with population genetics to reconstruct the geographic origin of a species

Reconstructing ancestral geographic origins is critical for understanding the long-term evolution of a species. Bayesian methods have been proposed to test biogeographic hypotheses while accommodating uncertainty in phylogenetic reconstruction. However, the problem that certain taxa may have a disproportionate influence on conclusions has not been addressed. Here, we infer the geographic origin of Drosophila simulans using 2,014 bp of the period locus from 63 lines collected from 18 countries. We also analyze two previously published datasets, alcohol dehydrogenase related and NADH:ubiquinone reductase 75 kDa subunit precursor. Phylogenetic inferences of all three loci support Madagascar as the geographic origin of D. simulans. Our phylogenetic conclusions are robust to taxon resampling and to the potentially confounding effects of recombination. To test our phylogenetically derived hypothesis we develop a randomization test of the population genetics prediction that sequences from the geographic origin should contain more genetic polymorphism than those from derived populations. We find that the Madagascar population has elevated genetic polymorphism relative to non-Madagascar sequences. These data are corroborated by mitochondrial DNA sequence data.     

Exploitation of mitochondrial nad6 as a complementary marker for studying population variability in Lepidoptera

The applicability of mitochondrial nad6 sequences to studies of DNA and population variability in Lepidoptera was tested in four species of economically important moths and one of wild butterflies. The genetic information so obtained was compared to that of cox1 sequences for two species of Lepidoptera. nad6 primers appropriately amplified all the tested DNA targets, the generated data proving to be as informative and suitable in recovering population structures as that of cox1. The proposal is that, to obtain more robust results, this mitochondrial region can be complementarily used with other molecular sequences in studies of low level phylogeny and population genetics in Lepidoptera.     

Evolution of the A+T-rich region of mitochondrial DNA in the melanogaster species subgroup of Drosophila

We determined the nucleotide sequences of two regions in the A+T-rich region of mitochondrial DNA (mtDNA) in the siI and siII types of D. simulans, the maII type of D. mauritiana, and D. sechellia. The sequences were aligned with those of the corresponding regions of siIII of D. simulans and maI of D. mauritiana, D. melanogaster, and D. yakuba. The type I and type II elements and the T-stretches were detected in all eight of the mtDNA types compared, indicating that the three elements are essential in the A+T-rich region of this species subgroup. The alignment revealed several short repetitive sequences and relatively large deletions in the central portions of the region. In the highly conserved sequence elements in the type II elements, the substitution rates were not uniform among lineages and acceleration in the substitution rate might have been due to loss of functional constraint in the stem–loop-forming sequences predicted in the type II elements. Patterns of nucleotide substitutions observed in the A+T-rich region were further compared with those in the coding regions and in the intergenic regions of mtDNA. Substitutions between A and T were particularly repressed in the highly conserved sequence elements and in the intergenic regions compared with those in the A+T-rich region excluding the highly conserved sequence elements and in the fourfold degenerate sites in the coding regions. The functional and structural characteristics of the A+T-rich region that might be involved in this substitutional bias are discussed.     

The phylogenetic relationships of introduced Aphelinus (Hymenoptera: Aphelinidae), biological control agents of the Russian wheat aphid (Homoptera: Aphididae)

Abstract  Several species of Aphelinus have been introduced to the US from the Old World for biological control of the Russian wheat aphid, Diuraphis noxia (Modvilko). Reproductive incompatibility has been observed among populations collected from different geographic areas. We examined whether or not the reproductive incompatibility between strains of A. asychis was caused by distant phylogenetic relationships. Ribosomal DNA sequences of internal transcribed spacers 2 (ITS2) were analyzed in several species of Aphelinus collected from multiple sites of Europe and Asia. The phylogenetic analysis showed that strains within the species A. albipodus and A. asychis are not monophyletic, and two clearly divergent clades were revealed among sequenced samples. Our results suggest that the reproductive incompatibility between three exotic strains of A. asychis was more likely caused by divergence of phylogeny than by symbiotic bacteria.     

Evolution of Wolbachia-induced cytoplasmic incompatibility in Drosophila simulans and D. sechellia

Abstract.— The intracellular bacterium Wolbachia invades arthropod host populations through various mechanisms, the most common of which being cytoplasmic incompatibility (CI). CI involves elevated embryo mortality when infected males mate with uninfected females or females infected with different, incompatible Wolbachia strains. The present study focuses on this phenomenon in two Drosophila species: D. simulans and D. sechellia . Drosophila simulans populations are infected by several Wolbachia strains, including w Ha and w No. Drosophila sechellia is infected by only two Wolbachia : w Sh and w Sn. In both Drosophila species, double infections with Wolbachia are found. As indicated by several molecular markers, w Ha is closely related to w Sh, and w No to w Sn. Furthermore, the double infections in the two host species are associated with closely related mitochondrial haplotypes, namely si I (associated with w Ha and w No in D. simulans ) and se (associated with w Sh and w Sn in D. sechellia ). To test the theoretical prediction that Wolbachia compatibility types can diverge rapidly, we injected w Sh and w Sn into D. simulans , to compare their CI properties to those of their sister strains w Ha and w No, respectively, in the same host genetic background. We found that within each pair of sister strains CI levels were similar and that sister strains were fully compatible. We conclude that the short period for which the Wolbachia sister strains have been evolving separated from each other was not sufficient for their CI properties to diverge significantly.     

Distribution of sequences homologous to the impCAB operon of TP110 among bacterial plasmids of different incompatibility groups

Summary Mutagenic DNA repair is a function of many naturally occurring plasmids belonging to several different incompatibility groups. A DNA probe corresponding to the impCAB operon of the IncIl plasmid TP110, which encodes such functions, was used to investigate the distribution of homologous sequences in both related and unrelated plasmids. Southern blotting was used to demonstrate considerable sequence conservation amongst a number of plasmid types, with imp-related sequences being found on plasmids belonging to the I1, I1/B, B and FIV incompatibility groups. However, no homology was detected amongst plasmids of the N and L/M incompatibility groups, many of which carry functionally similar gene clusters. It appears that sequences determining mutagenic repair functions have been largely conserved within any one incompatibility group, but that significant divergent evolution has occurred between groups.     

Rates and Patterns of Scndna and Mtdna Divergence within the Drosophila Melanogaster Subgroup

Levels of DNA divergence among the eight species of the Drosophila melanogaster subgroup and D. takahashii have been determined using the technique of DNA-DNA hybridization. Two types of DNA were used: single-copy nuclear DNA (scnDNA) and mitochondrial DNA (mtDNA). The major findings are: (1) A phylogeny has been derived for the group based on scnDNA which is congruent with chromosomal data, morphology, and behavior. The three homosequential species, simulans, sechellia, and mauritiana, are very closely related; the scnDNA divergence indicate the two island species are a monophyletic group. (2) The rates of change of scnDNA and mtDNA are not greatly different; if anything scnDNA evolves faster than mtDNA. (3) The rates of scnDNA evolution are not closely correlated to chromosomal (inversion) evolution. (4) The Drosophila genome appears to consist of two distinct classes of scnDNA with respect to rate of evolutionary change, a very rapidly evolving fraction and a relatively conservative fraction. (5) The absolute rate of change was estimated to be at least 1.7% nucleotide substitution per one million years. (6) DNA distance estimates based on restriction site variation are correlated with distances based on DNA-DNA hybridization, although the correlation is not very strong.     

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

Nuclear gene sequences from a late pleistocene sloth coprolite

The determination of nuclear DNA sequences from ancient remains would open many novel opportunities such as the resolution of phylogenies, the sexing of hominid and animal remains, and the characterization of genes involved in phenotypic traits. However, to date, single-copy nuclear DNA sequences from fossils have been determined only from bones and teeth of woolly mammoths preserved in the permafrost. Since the best preserved ancient nucleic acids tend to stem from cold environments, this has led to the assumption that nuclear DNA would be retrievable only from frozen remains. We have previously shown that Pleistocene coprolites stemming from the extinct Shasta sloth (Nothrotheriops shastensis, Megatheriidae) contain mitochondrial (mt) DNA from the animal that produced them as well as chloroplast (cp) DNA from the ingested plants. Recent attempts to resolve the phylogeny of two families of extinct sloths by using strictly mitochondrial DNA has been inconclusive. We have prepared DNA extracts from a ground sloth coprolite from Gypsum Cave, Nevada, and quantitated the number of mtDNA copies for three different fragment lengths by using real-time PCR. We amplified one multicopy and three single-copy nuclear gene fragments and used the concatenated sequence to resolve the phylogeny. These results show that ancient single-copy nuclear DNA can be recovered from warm, arid climates. Thus, nuclear DNA preservation is not restricted to cold climates.     

Parallel evolution in eight-barbel loaches of the genus Lefua (Balitoridae, Cypriniformes) revealed by mitochondrial and nuclear DNA phylogenies

The evolutionary history of eight-barbel loaches of the genus Lefua contains important phylogenetic information that will aid in resolution of the faunal formations and evolutionary histories of Japanese and East Asian freshwater fishes. Our sequencing of the mitochondrial D-loop region in a large number of samples allowed construction of the most comprehensive phylogeny of these loaches to date; we demonstrated monophyly of five Lefua species and identified populations of Lufua. sp. and Lefua echigonia. Loaches inhabiting the Tokai region in Japan were morphologically and ecologically indistinguishable from Lefua sp. However, they were included in the L. echigonia lineage. We determined a novel phylogeny by sequencing the nuclear ribosomal S7 subunit and showed that nuclear DNA phylogeny essentially matched the mitochondrial DNA phylogeny. Loaches from the Tokai region were part of the L. echigonia lineage, indicating parallel evolution between Tokai loaches and Lefua sp. in western Japan. We presented the most robust phylogeny to date using concatenated mitochondrial and nuclear sequences. The wealth of molecular information allowed us to speculate on evolutionary processes in the genus Lefua.     

Unidirectional incompatibility in Drosophila simulans: inheritance, geographic variation and fitness effects

In California, Drosophila simulans females from some populations (type W) produce relatively few adult progeny when crossed to males from some other populations (type R), but the productivity of the reciprocal cross is comparable to within-population controls. These two incompatibility types are widespread in North America and are also present elsewhere. Both types sometimes occur in the same population. Type R females always produce type R progeny irrespective of the father's type. However, matings between R males and females from stocks classified as type W produce type R progeny at low frequency. This suggests rare paternal transmission of the R incompatibility type, as we have found no evidence for segregation of incompatibility types in the W stocks. There is quantitative variation among type R lines for compatibility with W females, but not vice versa. Population cage studies and productivity tests suggest that deleterious side effects are associated with the type R cytoplasm.