Mitochondrial DNA evolution in the obscura species subgroup of Drosophila

Introgression in the Drosophila subobscura--D. Madeirensis sister species: evidence of gene flow in nuclear genes despite mitochondrial differentiation

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.

Phylogenetic relationships between Drosophila subobscura, D. guanche and D. madeirensis based on Southern analysis of heat shock genes

A Southern analysis of genomic DNA using Drosophila melanogaster probes for the major heat shock protein genes (Hsp82, Hsp70, Hsps encoding small proteins) was made to study the phylogenetic relationships between three Drosophila species belonging to the obscura group (D. subobscura, D. guanche, and D. madeirensis). The phylogenetic trees showed that D. madeirensis and D. subobscura are the most closely related species in the triad, while D. guanche is the most distantly related one. As in other Drosophila species, Hsp82 is a single copy gene in D. subobscura, D. guanche, and D. madeirensis, while Hsp70 and Hsps, which encode small proteins, are genic families. At least four sequences homologous to D. melanogaster Hsp70 were found in the obscura group species. These species have sequences which showed similarity with the four small Hsps of D. melanogaster.

Distribution of gypsy sequences in Drosophila species of the obscura subgroup

Eight Drosophila species of the obscura subgroup were screened for sequences homologous to the gypsy retrotransposon of D. melanogaster. Molecular characterization of gypsy sequences was first approached through digesting genomic DNAs from these obscura species with appropriate restriction enzymes and subjecting them to Southern blot analysis. The results of this analysis indicate that gypsy-homologous sequences are well conserved among species of the obscura subgroup. With the exception of D. guanche, all other species bear a 7 kb Xho I fragment that represents the complete element in D. melanogaster. Lower molecular weight fragments that could be deleted elements, are shared by different species. Both types of element probably existed before the divergence of this subgroup. Two different species clusters could be established on the basis of hybridization patterns, one represented by D. subobscura and its relative species D. guanche and D. madeirensis, and the other, in which D. obscura, D. tristis and D. subsilvestris are included.

Evolutionary dynamics of the SGM transposon family in the Drosophila obscura species group

SGM (Drosophila subobscura, Drosophila guanche, and Drosophila madeirensis) transposons are a family of transposable elements (TEs) in Drosophila with some functional and structural similarities to miniature inverted-repeat transposable elements (MITEs). These elements were recently active in D. subobscura and D. madeirensis (1-2 MYA), but in D. guanche (3-4 MYA), they gave rise to a species-specifically amplified satellite DNA making up approximately 10% of its genome. SGM elements were already active in the common ancestor of all three species, giving rise to the A-type specific promoter section of the P:-related neogene cluster. SGM sequences are similar to elements found in other obscura group species, such as the ISY elements in D. miranda and the ISamb elements in Drosophila ambigua. SGM elements are composed of different sequence modules, and some of them, i.e., LS and LS-core, are found throughout the Drosophila and Sophophora radiation with similarity to more distantly related TEs. The LS-core module is highly enriched in the noncoding sections of the Drosophila melanogaster genome, suggesting potential regulatory host gene functions. The SGM elements can be considered as a model system elucidating the evolutionary dynamics of mobile elements in their arms race with host-directed silencing mechanisms and their evolutionary impact on the structure and composition of their respective host genomes.

Evolution of gypsy endogenous retrovirus in the Drosophila obscura species group

The Ty3/gypsy family of retroelements is closely related to retroviruses, and some of their members have an open reading frame resembling the retroviral gene env. Sequences homologous to the gypsy element from Drosophila melanogaster are widely distributed among Drosophila species. In this work, we report a phylogenetic study based mainly on the analysis of the 5' region of the env gene from several species of the obscura group, and also from sequences already reported of D. melanogaster, Drosophila virilis, and Drosophila hydei. Our results indicate that the gypsy elements from species of the obscura group constitute a monophyletic group which has strongly diverged from the prototypic D. melanogaster gypsy element. Phylogenetic relationships between gypsy sequences from the obscura group are consistent with those of their hosts, indicating vertical transmission. However, D. hydei and D. virilis gypsy sequences are closely related to those of the affinis subgroup, which could be indicative of horizontal transmission.

Reevaluation of phylogeny in the Drosophila obscura species group based on combined analysis of nucleotide sequences

The Drosophila obscura species group has served as an important model system in many evolutionary and population genetic studies. Despite the amount of study this group has received, some phylogenetic relationships remain unclear. While individual analysis of different nuclear, mitochondrial, allozyme, restriction fragment, and morphological data partitions are able to discern relationships among closely related species, they are unable to resolve relationships among the five obscura species subgroups. A combined analysis of several nucleotide data sets is able to provide resolution and support for some nodes not seen or well supported in analyses of individual loci. A phylogeny of the obscura species group based on combined analysis of nucleotide sequences from six mitochondrial and five nuclear loci is presented here. The results of several different combined analyses indicate that the Old World obscura and subobscura subgroups form a monophyletic clade, although they are unable to resolve the relationships among the major lineages within the obscura species group.

GEM, a cluster of repetitive sequences in the Drosophila subobscura genome

GEM is a new family of repetitive sequences detected in the D. subobscura genome. Two of the four described GEM elements encompass a heterogeneous central module, with no detectable ORF, flanked by two long inverted repeats. These elements are composed of a set of repetitive modules, which are inverted repeat (IR), direct repeat (DR), palindromic sequence (PS), long sequence (LS) and short sequence (SS). These five modules can be found either clustered or dispersed as single modules in the D. subobscura genome, in euchromatic and heterochromatic regions. In addition to the 3' region of Adh retrosequences, single IR and LS blocks were found associated with the promoter region of different genes, in particular, LS-like blocks have also been found associated with functional genes in D. melanogaster and D. virilis. Conversely, the DR block is highly similar to satellite DNAs from some other species of the obscura group. In addition, GEM elements share some structural features with IS elements described in different Drosophila species. It is likely that both GEM and IS sequences would be vestiges of an ancestral transposable element.

The phylogenetic position of Drosophila eskoi deduced from P element and Adh sequence data

PCR screening with primers specific for the T-, M-, and O-type P element subfamilies was performed to investigate the interspecific distribution in 18 species and to reconstruct the phylogenetic history of the various types within the obscura species group. T-type elements occur in D. ambigua, D. tristis, D. obscura, D. subsilvestris, and D. eskoi. In the genomes of D. subobscura, D. madeirensis, and D. guanche they are present in the form of terminally truncated T-type derivatives. The wide distribution suggests that the T-type subfamily had a long evolutionary history in the obscura lineage. In contrast, the patchy occurrence of M- and O-type elements can be ascribed to four independent events of horizontal invasion of different lineages. The cladogenesis of the obscura group was investigated using a partial sequence of the Adh gene as a marker. In contrast to earlier findings, the position of D. eskoi had to be revised. D. eskoi appears as the closest relative of the D. ambigua clade, whereas D. tsukubaensis is the sister taxon of the species pair D. bifasciata/D. imaii. This result is in good accordance with the P element data, where high sequence similarity (95%) was found among the T-type elements of D. eskoi and those of D. ambigua and D. tristis.

Molecular and chromosomal phylogeny in the obscura group of Drosophila inferred from sequences of the rp49 gene region

A region of approximately 1.6 kb encompassing the ribosomal protein 49 gene (rp49) has been sequenced and compared in nine species of the obscura group of Drosophila: four species belonging to the obscura subgroup, three to the pseudoobscura subgroup, and two to the affinis subgroup. Our data provide strong support that the nearctic species (pseudoobscura and affinis subgroups) are monophyletic and place D. bifasciata with the other species of the obscura subgroup. Nucleotide sequence information at the rp49 gene region (located very close to one of the breakpoints of inversion O3) has also been used to infer the phylogeny of the O chromosome in the subobscura species cluster. Analysis based both on parsimony-informative sites and on genetic distances confirms that the O3 gene arrangement, present in D. guanche (together with inversion g) and in D. madeirensis, is ancestral to gene arrangements O3 + 4 and Ost present in extant populations of D. subobscura.

Identifying and distinguishing sibling species in the Tetrahymena pyriformis complex (Ciliophora, Oligohymenophorea) using PCR/RFLP analysis of nuclear ribosomal DNA

We describe a riboprinting strategy for identifying and distinguishing among sibling species in the Tetrahymena pyriformis complex. It involves use of the polymerase chain reaction to amplify a large segment of the nuclear ribosomal DNA and internal transcribed spacers, and digestion of this DNA with restriction enzymes. Unique restriction fragment length patterns or haplotypes were then used to distinguish species into: (1) six taxa that were identifiable to the species level, (2) eight taxa that were separated into four pairs, and (3) a group of eight taxa that were identical to each other. The latter result indicates that a more variable molecule is needed to distinguish the most closely related species in the complex. There was no intraspecific variation between two strains from one species (Tetrahymena thermophila) nor among multiple isolates from another species (Tetrahymena empidokyrea). This approach provides an alternative to traditional techniques for identifying T. pyriformis species that require living reference specimens and/or that reveal high levels of intraspecific variation.

Structure and organization of the P element related sequences in Drosophila madeirensis

The P element homologous sequences of the two closely related species Drosophila guanche and Drosophila subobscura represent a very special case of transposable-element derivatives. Although they have lost the regions known to be essential for P transposition by random mutations, all of them have selectively conserved the coding capacity for "P-repressor-like" proteins during the past few millions years. In both species, they are tandemly amplified in a single euchromatic gene cluster at equivalent chromosomal positions. In contrast, Drosophila madeirensis, an endemic species that is very closely related to both D. subobscura and D. guanche, harbours an additional P homologous site. Several mechanisms can be invoked to explain the generation of the new site in this species. In this work we present several molecular and cytological data in order to elucidate the possible evolutionary origin of the P derivatives of D. madeirensis.

The Drosophila montium subgroup species. Phylogenetic relationships based on mitochondrial DNA analysis

Mitochondrial DNA (mtDNA) restriction site maps for three Drosophila montium subgroup species of the melanogaster species group, inhabiting Indian and Afrotropical montium subgroup territories, were established. Taking into account previous mtDNA data concerning six oriental montium species, a phylogeny was established using distance-matrix and parsimony methods. Both genetic diversity and mtDNA size variations were found to be very narrow, suggesting close phylogenetic relationships among all montium species studied. The phylogenetic trees that were constructed revealed three main lineages for the montium subgroup species studied: one consisting of the Afrotropical species Drosophila seguyi, which is placed distantly from the other species, one comprising the north-oriental (Palearctic) species, and one comprising the southwestern (south-oriental, Australasian, Indian, and Afrotropical) species. The combination of the mtDNA data presented here with data from other species belonging to the melanogaster and obscura subgroups revealed two major clusters: melanogaster and obscura. The melanogaster cluster is further divided into two compact lineages, comprising the montium subgroup species and the melanogaster complex species; the species of the other complex of the melanogaster subgroup, yakuba, disperse among the obscura species. The above grouping is in agreement with the mtDNA size variations of the species. Overall, among all subgroups studied, the species of the montium subgroup seem to be the most closely related.

Structure and expression of clustered P element homologues in Drosophila subobscura and Drosophila guanche

Sequence relationships and functional aspects were analysed in the P element homologues of Drosophila subobscura (Ds) and D. guanche (Dg). In both species, the P homologues are clustered at a single genomic position. They lack the characteristic terminal structures of actively transposing P elements, but they have the coding capacity for a 66-kDa 'repressor-like' protein. Two different types of cluster units (G-type and A-type) can be distinguished. The A-type unit, which is present in multiple copies, is transcribed in adult flies. In contrast, the G-type unit has a much lower copy number and is apparently not expressed. In Dg, the isolated G-type sequence carries a 420-bp insertion in the promoter region, which is probably responsible for inactivation. Sequence comparisons of different cluster units show that differentiation of the two types precedes the lineage split of these species. Substitution rates of the deduced proteins reveal two distinct subregions: high variability at the N terminus and strong sequence conservation in the rest of the protein. The variable region contains motifs characteristic of DNA-binding proteins. Adaptive diversification of the cluster units towards specific binding properties might be a plausible explanation for variability in the N-termini. Both unit types have lost the weak promoter region characteristic of P transposons. In the A-type unit, a new promoter has been formed which is apparently composed of parts of insertion sequences derived from two different mobile elements.(ABSTRACT TRUNCATED AT 250 WORDS)

Molecular identification and phylogeny of parasitic wasp species (Hymenoptera: Trichogrammatidae) by mitochondrial DNA RFLP and RAPD markers

Mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP) and random amplified polymorphic DNA (RAPD) were tested for their ability to discriminate between seven species of minute parasitic wasps belonging to the genus Trichogramma. They proved to be reliable species-diagnostic molecular markers. Pairwise comparisons of the mtDNA restriction maps revealed considerable differentiation among the seven species. Percentage of common restriction sites ranged from 30% to 83%. Phylogenetic analyses performed either on the mtDNA nucleotide distance matrix or on the matrix of the restriction site-state generated a tree congruent with those based on allozymes and morphology. RAPD procedures also revealed species-specific banding patterns and seem promising for a rapid and easy identification of Trichogramma species. Moreover, for some Trichogramma species, RAPD banding patterns might be informative of the phylogenetic relatedness.

Phylogeny of the Drosophila obscura species group deduced from mitochondrial DNA sequences

Approximately 2 kb corresponding to different regions of the mtDNA of 14 different species of the obscura group of Drosophila have been sequenced. In spite of the uncertainties arising in the phylogenetic reconstruction due to a restrictive selection toward a high mtDNA A+T content, all the phylogenetic analysis carried out clearly indicate that the obscura group is formed by, at least, four well-defined lineages that would have appeared as the consequence of a rapid phyletic radiation. Two of the lineages correspond to monophyletic subgroups (i.e., affinis and pseudoobscura), whereas the obscura subgroup remains heterogeneous assemblage that could be reasonably subdivided into at least two complexes (i.e., subobscura and obscura).

Two distinct P element subfamilies in the genome of Drosophila bifasciata

The genome of Drosophila bifasciata harbours two distinct subfamilies of P-homologous sequences, designated M-type and O-type elements based on similarities to P element sequences from other species. Both subfamilies have some general features in common: they are of similar length (M-type: 2935 bp, O-type: 2986 bp), are flanked by direct repeats of 8 bp (the presumptive target sequence), contain terminal inverted repeats, and have a coding region consisting of four exons. The splice sites are at homologous positions and the exons have the coding capacity for proteins of 753 amino acids (M-type) and 757 amino acids (O-type). It seems likely that both types of element represent functional transposons. The nucleotide divergence of the two P element subfamilies is high (31%). The main structural difference is observed in the terminal inverted repeats. Whereas the termini of M-type elements consists of 31 bp inverted repeats, the inverted repeats of the O-type elements are interrupted by non-complementary stretches of DNA, 12 bp at the 5' end and 14 bp at the 3' end. This peculiarity is shared by all members of the O-type subfamily. Comparison with other P element sequences indicates incongruities between the phylogenies of the species and the P transposons. M-type and O-type elements apparently have no common origin in the D. bifasciata lineage. The M-type sequence seems to be most closely related to the P element from Scaptomyza pallida and thus could be considered as a more recent invader of the D. bifasciata gene pool. The origin of the O-type elements cannot be unequivocally deduced from the present data.(ABSTRACT TRUNCATED AT 250 WORDS)

Speciation in the Artemia genus: mitochondrial DNA analysis of bisexual and parthenogenetic brine shrimps

From the cloned mitochondrial DNAs (mtDNAs) isolated from two bisexual species, one Mediterranean, Artemia salina, and one American, Artemia franciscana, and two parthenogenetic (diploid and tetraploid) strains of Artemia parthenogenetica collected in Spain, physical maps have been constructed and compared. They are extremely different among themselves, much more than the differences between Drosophila melanogaster and D. yakuba and in the same range of different mammalian species such as mouse/rat or man/cow. The nucleotide sequences of two regions of mtDNA encoding parts of the cytochrome c oxidase subunit I (COI) and cytochrome b (Cytb) genes have been determined in the two bisexual species and the two parthenogenetic strains. Comparisons of these sequences have revealed a high degree of divergence at the nucleotide level, averaging more than 15%, in agreement with the differences found in the physical maps. The majority of the nucleotide changes are silent and there is a strong bias toward transitions, with the C<==>T substitutions being highly predominant. The evolutionary distance between the two Artemia parthenogenetica is high and there is no clear relationship with any of the bisexual species, including the one present nowadays in Spain. Using a combination of molecular (mtDNA) and morphological markers it is possible to conclude that all of these Artemia isolates should be actually considered as belonging to different species, even the two Artemia parthenogenetica diploidica and tetraploidica.

A heterochromatic P sequence in the D. subobscura genome

The study of a heterochromatic P sequence of D. subobscura reveals that it is a degraded element, located at the centromeric region of the A chromosome (X chromosome in this species), and that it is strongly diverged from the euchromatic P sequences previously described in this species. This heterochromatic sequence is composed of some P element fragments embedded in undefined beta-heterochromatic sequences. These mosaic P sequences do not show any transcriptional activity and seem to be ancient parasites of the D. subobscura genome. Phylogenetic analyses indicate that both the euchromatic and heterochromatic P sequences of D. subobscura could come from an ancestral element which was present before the divergence of the subobscura species cluster.

Analysis of mitochondrial DNA in two species of the bipectinata species complex of Drosophila

Drosophila bipectinata and D. malerkotliana are members of the bipectinata species complex of the ananassae subgroup in the melanogaster species group of Drosophila. The mtDNA from 18 isofemale strains of these species was analyzed, using 13 restriction endonucleases. Altogether, eight mtDNA haplotypes were detected, of which the haplotype 1B was shared by the two species. Restriction cleavage map of mtDNA of these species was constructed. The net nucleotide substitutions per site calculated between these species was found to be 0.0002. This value appeared to be relatively much lower than expected at the interspecific level, even lower than that reported between two subspecies of Drosophila despite both being good species. This extreme closeness of their mtDNAs is discussed in the light of recent findings.

Differences in gene activity in a Drosophila species cluster belonging to the Obscura group

The polytene chromosome puffing pattern of Drosophila madeirensis was established and compared with those of the related species D. subobscura and D. guanche. A total of 145 loci, active in some of the 12 developmental stages analysed, were described, 38 of which were found to form the puffing pattern characteristic to this species. Taking into account the number of puffs as well as the mean puff expression, D. madeirensis shows a similar activity level to D. guanche, both species being less active than D. subobscura. The low gene activity of D. madeirensis and D. guanche was explained as a consequence of their ecological characteristics.

Characterization of the length polymorphism in the A + T-rich region of the Drosophila obscura group species

In the twelve Drosophila obscura group species studied, belonging to the affinis, obscura, and pseudoobscura subgroups, the mitochondrial DNA length ranges from 15.8 to 17.2 kb. This length polymorphism is mainly due to insertions/deletions in the variable region of the A + T-rich region. In addition, one species (D. tristis) possess a tandem duplication of a 470-bp fragment that contains the replication origin. The same duplication has occurred at least twice in the Drosophila evolutionary history due to the fact that the repetition is analogous to repetitions found in the four species of the D. melanogaster complex. By comparing the nucleotide sequence of the conserved region in D. ambigua, D. obscura, D. yakuba, D. teissieri, and D. virilis, we show the presence of a secondary structure, likely implied in the replication origin, which could favor the generation of this kind of duplications. Finally, we propose that the high A and T content in the variable region of the A + T-rich region favors the formation of less-stable secondary structures, which could explain the generation of minor insertion/deletions found in this region.

Nuclear and mitochondrial ribosomal RNA variability in the obscura group of Drosophila

Parts of 28S (nuclear) and 12S (mitochondrial) ribosomal RNA of Palearctic, Nearctic and African species of the obscura group have been sequenced by the direct method of sequencing. Rates of nucleotide substitutions in both molecules were compared. The nucleotide divergence is higher in the mitochondrial rRNA. Average distances of species taken in pairwise were compared to results obtained with the melanogaster subgroup: the divergence of nuclear rRNA appears lower, that of the mtDNA higher whereas genetic distances (allozymes) and sncDNA distances are similar. Noticeable variability of evolutionary rates can be observed even in low taxonomical levels. Phylogenetic trees for the obscura group are in general agreement with those obtained with other characters.