The 5S ribosomal genes in the Drosophila melanogaster species subgroup. Nucleotide sequence of a 5S unit from Drosophila simulans and Drosophila teissieri

THE REPRODUCTIVE RELATIONSHIPS OF DROSOPHILA SECHELLIA WITH D. MAURITIANA, D. SIMULANS, AND D. MELANOGASTER FROM THE AFROTROPICAL REGION

Hybridization tests among the four sibling species of the Drosophila melanogaster complex were made to determine the reproductive status of the recently discovered D. sechellia (which is endemic to a few islands and islets of the Seychelles archipelago) with regard to its three close relatives, D. mauritiana (endemic to Mauritius) and Afrotropical strains of the two cosmopolitan species D. melanogaster and D. simulans. Interstrain variation in the ability to hybridize with other species was also analyzed for D. melanogaster and D. simulans. D. mauritiana and D. simulans appear to be more weakly isolated from each other than either species is from D. sechellia. A striking unilateral mating success is observed in the cross of D. sechellia with D. simulans. The most extreme isolation is between D. melanogaster and its three siblings. Variation in the ability of strains to hybridize is observed in heterospecific crosses between D. simulans and either D. melanogaster or D. mauritiana.

RACE FORMATION, SPECIATION, AND INTROGRESSION WITHIN DROSOPHILA SIMULANS, D. MAURITIANA, AND D. SECHELLIA INFERRED FROM MITOCHONDRIAL DNA ANALYSIS

Mitochondrial DNA cleavage maps from three chromosomally homosequential species Drosophila simulans, D. mauritiana, and D. sechellia, were established for 12 restriction enzymes. One isofemale strain was studied in D. sechellia (se), 13 in D. simulans, and 17 in D. mauritiana: in the last two species, respectively, three (siI, II, and III) and two (maI and II) cleavage morphs were found. The evolutionary relationships based on mtDNA cleavage map comparisons show that the maI and se mtDNAs are internal branches of the phylogenetic tree of the D. simulans mtDNA. D. mauritiana and D. sechellia species appear to be derived from a population of D. simulans which carried an ancestral form of the current siI mtDNA type. In addition, two cleavage morphs (siIII [only present in D. simulans from Madagascar] and maI) appeared to be identical, although found in different species. We present a speculative interpretation of data on biogeography and hybridization which is consistent with the hypothesis of a recent introgression of mitochondrial DNA of D. simulans from Madagascar into D. mauritiana.

Abridged 5S rDNA units in sea beet (Beta vulgaris subsp. maritima)

Amplification by polymerase chain reaction of the 5S rDNA repeat units of Beta vulgaris subsp. maritima resulted in a 350-bp product corresponding to the full-length 5S unit, but also revealed 4 abridged unit classes, each with a deletion that removed most of the spacer and 12–76 bp of the coding sequence. Each abridged type lacks at least 1 of the conserved elements involved in transcription of the 5S gene, and so appear to be nonfunctional. Network analysis revealed that the abridged units are evolving in the same manner as the full-length versions.Key words: 5S rDNA, Beta vulgaris subsp. maritima, network analysis, sea beet.

Evolution of 5S rRNA gene families in Drosophila

In Drosophila virilis, the three clusters of 5S rRNA genes on chromosome 5 comprise two different gene families (B and C), which differ profoundly in the organization of their spacer sequences. While C-type genes, which are found in two of the clusters, exhibit a true repetitive character, the B-type genes of the third cluster are each embedded in completely different genomic environments. Southern blots of genomic DNA of different D. virilis subspecies, D. hydei and D. melanogaster probed with 5S rRNA gene spacer and coding sequences demonstrate the specificity of C-type sequences for the D. virilis species group. The comparative analysis of flanking sequences of 5S rRNA genes of D. virilis, members of the D. melanogaster species subgroup and of the blowfly Calliphora erythrocephala reveals the existence of conserved sequence motifs both in the 5' upstream and 3' downstream flanking regions. Their possible roles in the control of expression and processing of the 5S rRNA precursor molecule are discussed.

Polytene chromosomes as indicators of phylogeny in several species groups of Drosophila

BACKGROUND

Polytene chromosome banding patterns have long been used by Drosophila evolutionists to infer degree of relatedness among taxa. Recently, nucleotide sequences have preempted this traditional method. We place the classical Drosophila evolutionary biology tools of polytene chromosome inversion analysis in a phylogenetic context and assess their utility in comparison to nucleotide sequences.

RESULTS

A simultaneous analysis framework was used to examine the congruence of the chromosomal inversion data with more recent DNA sequence data in four Drosophila species groups - the melanogaster, virilis, repleta, and picture wing. Inversions and nucleotides were highly congruent with one another based on incongruence length difference and partitioned Bremer support values. Inversion phylogenies were less resolved because of fewer numbers of characters. Partitioned Bremer supports, corrected for the number of characters in each matrix, were higher for inversion matrices.

CONCLUSIONS

Polytene chromosome data are highly congruent with DNA sequence data and, when placed in a simultaneous analysis framework, are shown to be more information rich than nucleotide data.

Structural evolution of the Drosophila 5S ribosomal genes

We compare the 5S gene structure from nine Drosophila species. New sequence data (5S genes of D. melanogaster, D. mauritiana, D. sechellia, D. yakuba, D. erecta, D. orena, and D. takahashii) and already-published data (5S genes of D. melanogaster, D. simulans, and D. teissieri) are used in these comparisons. We show that four regions within the Drosophila 5S genes display distinct rates of evolution: the coding region (120 bp), the 5'-flanking region (54-55 bp), the 3'-flanking region (21-22 bp), and the internal spacer (149-206 bp). Intra- and interspecific heterogeneity is due mainly to insertions and deletions of 6-17-bp oligomers. These small rearrangements could be generated by fork slippages during replication and could produce rapid sequence divergence in a limited number of steps.

Deletions and amplifications of tandemly arranged ribosomal 5S genes internal to a P element occur at a high rate in a dysgenic context

We observed unusual kinds of rearrangements within tandemly clustered 5S genes internal to a P element in dysgenic context. Rearranged P transposons, initially containing eight 5S genes, were found to display discrete numbers of 5S genes, from 4 up to 17 units. Precise deletions and amplifications occurred at a high rate (40%), at both original and new insertion sites. These events can be explained by a "cut and paste" transposition model. Possible links between rearrangements due to dysgenic-like processes and concerted evolution are discussed.

I elements and the Drosophila genome

LINEs are a large class of transposable elements in eukaryotes. They transpose by reverse transcription of an RNA intermediate. I elements of Drosophila melanogaster belong to this class and are responsible for the I-R system of hybrid dysgenesis. Many results indicate that at the beginning of the century natural populations of this species were devoid of active I elements and that they were invaded by functional I elements in the last decades. Many Drosophila species contain both defective and active I elements. It seems that the latter were lost in Drosophila melanogaster before its spread throughout the world, and that the recent invasion results from the spread of functional elements originating either from another species by horizontal transfer or from an isolated population of the same species. These data are discussed, as well as their significance in evolutionary processes.

Two highly conserved transcribed regions in the 5S DNA repeats of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae

The 5S RNA genes of Caenorhabditis briggsae consist of approximately 65 copies of a 1 kb repeat unit and 20 copies of a related 0.7 kb repeat unit, organized in separate tandem clusters. DNA sequence comparisons with the 1kb 5S DNA repeat from the closely related nematode C. elegans show that the 5S RNA coding region is perfectly conserved. Both C. briggsae 1 kb and 0.7 kb repeats are also efficiently transcribed in vitro,suggesting that both represent functional 5S RNA genes. Surprisingly, a second block of 118 bp is also perfectly conserved between the 1 kb repeats,and is less well conserved in the 0.7 kb repeat. In C. elegans, this DNA is transcribed to produce and abundant 100 nt transcript (SL RNA) which participates in a trans-splicing process (Krause and Hirsh, Cell 49:753, 1987). This SL RNA region of the C. briggsae 1 kb 5S DNA repeat also appears to be transcribed in vivo, while the corresponding region of the 0.7 kb repeat is not.

A long interspersed repetitive element--the I factor of Drosophila teissieri--is able to transpose in different Drosophila species

Long interspersed repetitive elements (LINEs) are transposable elements present in many species. In mammals they are difficult to study because most of them are defective and their transposition frequency is low. The I factor of Drosophila melanogaster is a LINE element that is particularly interesting because its transposition occurs at high frequency during I-R hybrid dysgenesis. This phenomenon occurs when males from the class of inducer strains are crossed with females from the class of reactive strains. Inducer strains contain several complete 5.4-kilobase I factors at various sites on the chromosomal arms. Reactive strains are devoid of complete I factors. Many results indicate that active I factors have invaded the D. melanogaster genome recently. To study the evolutionary history of I elements, we have cloned and sequenced a potentially active I factor from Drosophila teissieri. It is flanked by a target-site duplication and terminates at the 3' end by tandem repeats of the sequence TAA. When introduced into the germ line of a reactive strain of D. melanogaster by P element-mediated transformation, it is able to transpose and induces hybrid dysgenesis. This strengthens the hypothesis of a recent reinvasion of the D. melanogaster genome by active I factors giving rise to the inducer strains. They could have originated by horizontal transfer from another species. Such events also could occur for other LINE elements and might explain the spread of new variants in mammalian genomes. Moreover, the results give a further insight into I factor functional organization.

An approach to study the evolution of the Drosophila 5S ribosomal genes using P-element transformation

The P-element-mediated gene transfer system was used to introduce Drosophila teissieri 5S genes into the Drosophila melanogaster genome. Eight transformed D. melanogaster strains that carry D. teissieri 5S mini-clusters consisting of 9-21 adjacent 5S units were characterized. No genetic exchanges between D. melanogaster and D. teissieri 5S clusters were detected over a 2-year survey of the eight strains. The occurrence of small rearrangements within the D. melanogaster 5S cluster was demonstrated in one of the transformed strains.

Bipartite structure of the 5S ribosomal gene family in a Drosophila melanogaster strain, and its evolutionary implications

Knowledge of multigenic family organization should provide insight into their mode of evolution. Accordingly, we characterized the 5S ribosomal gene family in the Drosophila melanogaster strain ry506. The 5S genes in this strain display a striking HindIII restriction difference compared to the "standard" D. melanogaster 5S genes. The sequence of three ry506 5S genes was determined. We show that the HindIII restriction site heterogeneity within the ry506 5S family most probably results from the same point mutation, suggesting that a single 5S variant was propagated into the 5S cluster of this strain. Furthermore, we demonstrate that the structural organization of the 5S genes in ry506 is a bipartite structure, i.e., that about 40% of the 5S genes constitute a HindIII+/HindIII- mixed cluster, while those remaining constitute an homogeneous HindIII- cluster. The events which might lead to such an heterogeneous pattern are discussed from an evolutionary point of view.

5S-rRNA genes in rice embryos

The 5S-rRNA from the ungerminated and 48-h-germinated rice embryos differs from the wheat, rye and maize by two nucleotides. The 48-h-germinated embryos contain another species of 5S-rRNA which differs by 3 nucleotides from the ungerminated embryos, thereby showing the expression of two 5S-rRNA genes during germination. The 5S-rRNA genes are present in tandem repeats of a 0.3-kb sequence with some length heterogeneity in the rice genome. The 5S-rRNA gene that was sequenced is identical to that of wheat and maize, except for two nucleotides, C and T, which are interchanged at positions 107 and 117. The insert of continuous 5S-rRNA gene in pBR322 was transcribed in vitro much more efficiently than the discontinuous gene. There was no homology between the 184-bp spacer sequence of 5S-rRNA genes in rice and other systems except the presence of the oligo(T) transcription terminator sequence.

Conserved 5' flank homologies in dipteran 5S RNA genes that would function on 'A' form DNA

We have sequenced the 480 base pair (bp) repeating unit of the 5S RNA genes of the Dipteran fly Calliphora erythrocephala and compared this sequence to the three known 5S RNA gene sequences from the Dipteran Genus Drosophila (1,2). A striking series of five perfectly conserved homologies identically positioned within the 5' flanks of all four Dipteran 5S RNA coding regions has thus been identified. The spacing (12-13 bp) between all of these homologies is typical of A form rather than B form DNA. Given that the eukaryotic 5S RNA gene specific initiation factor TFIIIA (3) is a DNA unwinding protein (4), a role for these Dipteran 5' flank homologies in initiation site selection on 5S RNA genes transiently unwound for transcription is suggested. One of the Dipteran homology blocks is highly conserved in sequence and position in all but one of the eukaryotic 5S RNA gene sequences known to date (17/18 genes). Its sequence (consensus: TATAAG) and position (average center: -26 bp) are highly reminiscent of the polymerase II gene 'TATA' box (5).