Nucleotide sequence comparison of the rp49 gene region between Drosophila subobscura and D. melanogaster

A 1.6-kb fragment encompassing the rp49 gene, which codes for a ribosomal protein, has been cloned and sequenced in Drosophila subobscura. The rp49 coding region has accumulated 46 nucleotide differences out of 402 bp since D. subobscura diverged from D. melanogaster. Forty-three percent of the effectively silent sites have changed since both species diverged. Both silent and replacement differences are distributed at random between the two exons of the gene. The frequency of silent differences in exons does not differ from that observed in the 5' leader sequence and in the intron. The frequency of silent differences in exon and intron sites is much greater than the number of amino acid replacement differences. This observation indicates strong purifying selection against amino acid replacements.      

Adaptation of Drosophila subobscura chromosomal inversions to climatic variables: the Balkan natural population of Avala

Genetica - The adaptive value of chromosomal inversions continues raising relevant questions in evolutionary biology. In many species of the Drosophila genus, different inversions have been...     

Genetic load and coadaptation of chromosomal inversions. II. O-chromosomes in Drosophila subobscura populations

We have analysed the inversion polymorphism and genetic load of O-chromosomes in three populations of D. subobscura from southeastern Europe. As expected for a central populations the inversion polymorphism was extensive. In a like fashion, the genetic load, in particular the frequency of lethals, was heavy in all three populations. There were significant differences in the frequency of moderately deleterious genes. These differences in viability can be attributed to balancing selection. A comparison of these two kinds of genetic polymorphism indicates that there are differences in mean viability among different gene arrangements of O-chromosomes in the three populations. The differences observed are due to an unequal distribution of various viability classes among O-chromosome gene arrangements. We here show for the first time a specific distribution of lethal genes among these arrangements within the Palearctic distribution area of D. subobscura. The lethal allelism test showed lethals are non-randomly associated with the Ost gene arrangement. The amount of genetic load is heavy in gene arrangements with a high frequency, in comparison with the ones with a low frequency. Lethal genes may be protected in combinations of low and moderate frequency gene arrangements that harbor more lethal genes, as the Ost in the one population. Some arrangements that are less protected against recombination have a higher load than ones that are more protected against recombination. This can be taken as evidence for coadaptation.     

DNA variation at the rp49 gene region in Drosophila madeirensis and D. subobscura from Madeira: inferences about the origin of an insular endemic species

An ~1.6 ‐ kb fragment spanning the rp49 gene was sequenced in 16 lines of Drosophila subobscura from Madeira and in 22 lines of the endemic species D. madeirensis. Nucleotide diversity in D. subobscura from Madeira (π=0.0081) was similar to that in lines from Spain carrying the O₃₊₄ chromosomal arrangement (π=0.0080). No significant genetic differentiation was detected between insular and continental O₃₊₄ lines of D. subobscura. These results are compatible both with a rather recent and massive colonization, and with multiple colonization events from the continent. Nucleotide diversity in D. madeirensis (π=0.0076) was similar to that in D. subobscura, which deviates from the expectation, under strict neutrality, of a lower level of variation in an insular species with a small population size. The observed numbers of shared polymorphisms and of fixed differences between D. madeirensis and D. subobscura are compatible with the isolation model of speciation, where shared polymorphisms are due to common ancestry.     

Nucleotide divergence of the rp49 gene region between Drosophila melanogaster and two species of the Obscura group of Drosophila

Summary A 2.1-kb SStI fragment including the rp49 gene and the 3 end of the -serendipity gene has been cloned and sequenced in Drosophila pseudoobscura. rp49 maps at region 62 on the tip of chromosome II of this species. Both the coding and flanking regions have been aligned and compared with those of D. subobscura. There is no evidence for heterogeneity in the rate of silent substitution between the rp49 coding region and the rate of substitutions in flanking regions, the overall silent divergence per site being 0.19. Noncoding regions also differ between both species by different insertions/deletions, some of which are related to repeated sequences. The rp49 region of D. pseudoobscura shows a strong codon bias similar to those of D. subobscura and D. melanogaster. Comparison of the rates of silent (K _S ) and nonsilent (K _a ) substitutions of the rp49 gene and other genes completely sequenced in D. pseudoobscura and D. melanogaster confirms previous results indicating that rp49 is evolving slowly both at silent and nonsilent sites. According to the data for the rp49 region, D. pseudoobscura and D. subobscura lineages would have diverged some 9 Myr ago, if one assumes a divergence time of 30 Myr for the melanogaster and obscura groups.Offprint requests to: C. Segarra     

DNA variation at the rp49 gene region of Drosophila simulans: evolutionary inferences from an unusual haplotype structure.

An approximately 1.3-kb region including the rp49 gene plus its 5' and 3' flanking regions was sequenced in 24 lines of Drosophila simulans (10 from Spain and 14 from Mozambique). Fifty-four nucleotide and 8 length polymorphisms were detected. All nucleotide polymorphisms were silent: 52 in noncoding regions and 2 at synonymous sites in the coding region. Estimated silent nucleotide diversity was similar in both populations (pi = 0.016, for the total sample). Nucleotide variation revealed an unusual haplotype structure showing a subset of 11 sequences with a single polymorphism. This haplotype was present at intermediate frequencies in both the European and the African samples. The presence of such a major haplotype in a highly recombining region is incompatible with the neutral equilibrium model. This haplotype structure in both a derived and a putatively ancestral population can be most parsimoniously explained by positive selection. As the rate of recombination in the rp49 region is high, the target of selection should be close to or within the region studied.     

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long‐standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion‐size and inversion‐shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.     

Genetic loads and coadaptivity of chromosomal inversions I. Inversion polymorphism and genetic load on chromosome O in a Drosophila subobscura population from Petnica

Inversion polymorphism on chromosome O and polymorphism for the viability of determining genes have been studied in a natural population of Drosophila subobscura from Petnica (Serbia). The range of inversion polymorphism and the abundance of particular gene arrangements in the study population agree with a general pattern of inversion polymorphism of D. subobscura in Europe. The data obtained on the amount of genetic loads show that the D. subobscura population from Petnica displays a moderate degree of that polymorphism, compared to the other studied populations of these species. Therefore, the D. subobscura population from Petnica could be tentatively classified as an ecologically central population. Examination association of chromosomal, thus, inversion polymorphism with gene polymorphism, in the form of genetic loads show that differences exist in the mean viability among certain gene arrangements. The distribution of deleterious genes among chromosome O gene arrangements were non-random.     

The genetics of Drosophila subobscura populations XVIII. Multiple insemination and sperm displacement in Drosophila subobscura

The multiple insemination of Drosophila subobscura females in nature was studied by the analysis of single female offspring using two highly polymorphic enzyme systems. On qualitative grounds it was found that 23% of these females were inseminated at least twice. The power of detection of one enzyme system was estimated using the already known multi-inseminated females detected by the other system. This permitted the true frequency estimation of multi-inseminated females, in nature, which is 0.276 (for 95% confidence limits it ranges from 0.186 to 0.417). In a laboratory experiment it was found that sperm from the second inseminating male displaces the sperm from the first male.     

Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajima's D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.     

Allozyme divergence between different gene arrangements in Drosophila subobscura

Allozyme and chromosomal polymorphism were analysed simultaneously in three population samples of Drosophila subobscura from Scotland, Germany and Italy. The study revealed significant genic divergence between the various chromosomal gene arrangements for the loci Acph, Lap and Pep, which are located inside inversions. No difference was found with the loci Me and Odh, which are both on the same chromosome but proximal to the structurally polymorphic region. Samples of a given gene arrangement taken from the different localities are genetically uniform. The data indicate that the observed non-random associations between enzyme genes and chromosomal inversions are the result of inhibited recombination in heterokaryotypes and gradual differentiation of the separated gene pools, rather than selection and coadaptation.     

Molecular population genetics of Drosophila subtelomeric DNA

DNA sequence surveys in yeast and humans suggest that the forces shaping telomeric polymorphism and divergence are distinctly more dynamic than those in the euchromatic, gene-rich regions of the chromosomes. However, the generality of this pattern across outbreeding, multicellular eukaryotes has not been determined. To characterize the structure and evolution of Drosophila telomeres, we collected and analyzed molecular population genetics data from the X chromosome subtelomere in 58 lines of North American Drosophila melanogaster and 29 lines of African D. melanogaster. We found that Drosophila subtelomeres exhibit high levels of both structural and substitutional polymorphism relative to linked euchromatic regions. We also observed strikingly different patterns of variation in the North American and African samples. Moreover, our analyses of the polymorphism data identify a localized hotspot of recombination in the most-distal portion of the X subtelomere. While the levels of polymorphism decline sharply and in parallel with rates of crossing over per physical length over the distal first euchromatic megabase pairs of the X chromosome, our data suggest that they rise again sharply in the subtelomeric region (approximately 80 kbp). These patterns of historical recombination and geographic differentiation indicate that, similar to yeast and humans, Drosophila subtelomeric DNA is evolving very differently from euchromatic DNA.     

Relationship between chromosomal polymorphism and wing size in a natural population of Drosophila subobscura

Chromosomal polymorphism and wing size (as a measure of body size) were analysed simultaneously in two samples of Drosophila subobscura from Barcelona, Spain. The very rich chromosomal polymorphism of this species makes it difficult to detect the relationship of this polymorphism with any phenotypical character. However, a positive significant regression of wing size on the percentage of the autosome length with standard arrangement was found. Furthermore, for each polymorphic chromosome, except for the J chromosome, an association between the most frequent arrangements and wing size was observed. This trend, which was the same in the two samples, was that expected according to the latitudinal clines of both characters.     

The relationship between allozyme and chromosomal polymorphism inferred from nucleotide variation at the Acph-1 gene region of Drosophila subobscura.

The Acph-1 gene region was sequenced in 51 lines of Drosophila subobscura. Lines differ in their chromosomal arrangement for segment I of the O chromosome (O(st) and O(3+4)) and in the Acph-1 electrophoretic allele (Acph-1(100), Acph-1(054), and Acph-1(>100)). The ACPH-1 protein exhibits much more variation than previously detected by electrophoresis. The amino acid replacements responsible for the Acph-1(054) and Acph-1(>100) electrophoretic variants are different within O(st) and within O(3+4), which invalidates all previous studies on linkage disequilibrium between chromosomal and allozyme polymorphisms at this locus. The Acph-1(>100) allele within O(3+4) has a recent origin, while both Acph-1(054) alleles are rather old. Levels of nucleotide variation are higher within the O(3+4) than within the O(st) arrangement except for nonsynonymous sites. The McDonald and Kreitman test shows a significant excess of nonsynonymous polymorphisms within O(st) when D. guanche is used as the outgroup. According to the nearly neutral model of molecular evolution, this excess is consistent with a smaller effective size of O(st) relative to O(3+4) arrangements. A smaller population size, a lower recombination, and a more recent bottleneck might be contributing to the smaller effective size of O(st).     

Linkage disequilibria between mtDNA haplotypes and chromosomal arrangements in a natural population of Drosophila subobscura

The association between mtDNA haplotypes and chromosomal arrangements in a natural population of Drosophila subobscura from Calvia (Balearic Islands, Spain) was studied in order to search for linkage disequilibria, in an attempt to explain the populational dynamics of the mtDNA haplotypes of this species in nature. The presence of Wolbachia was not detected. Two main haplotypes (I and II) were found, as well as a series of less common ones. The Tajima D-test seemed to indicate some kind of seasonal population bottleneck. An analysis of linkage disequilibrium and factorial analysis of correspondences detected an association between haplotype I and the J(ST) inversion and haplotype II and the J(1) inversion.     

Molecular evolution of a duplication: the sex-peptide (Acp70A) gene region of Drosophila subobscura and Drosophila madeirensis

In Drosophila melanogaster, the Acp70A gene, which is involved in thepostmating reactions of the female, is a single-copy gene. However, inDrosophila subobscura, the gene is duplicated and both copies aretranscribed. To study the molecular evolution of the duplication, a 2.1- kbfragment encompassing both copies of the duplication was sequenced for 10lines of D. subobscura and one line of Drosophila madeirensis. Estimates ofthe divergence between the two copies of the duplicated region and betweenthe two species studied, D. subobscura and D. madeirensis, revealed thatboth copies of the Acp70a gene had evolved independently since theirduplication. The ratio of nonsynonymous to silent divergence between copieswas generally higher than one. The McDonald and Kreitman test revealed anexcess of nonsynonymous changes fixed since the duplication and before thesplit of the D. subobscura and D. madeirensis lineages. These results pointto natural selection driving protein evolution after the duplication.Specifically, adaptive evolution appears to have caused the initialdifferentiation between copies of the N-terminal parts of the proteins,while purifying selection could be responsible for the high conservation ofthe C- terminal parts.     

Molecular population genetics of X-linked genes in Drosophila pseudoobscura

This article presents a nucleotide sequence analysis of 500 bp determined in each of five X-linked genes, runt, sisterlessA, period, esterase 5, and Heat-shock protein 83, in 40 Drosophila pseudoobscura strains collected from two populations. Estimates of the neutral migration parameter for the five loci show that gene flow among D. pseudoobscura populations is sufficient to homogenize inversion frequencies across the range of the species. Nucleotide diversity at each locus fails to reject a neutral model of molecular evolution. The sample of 40 chromosomes included six Sex-ratio inversions, a series of three nonoverlapping inversions that are associated with a strong meiotic drive phenotype. The selection driven by the Sex-ratio meiotic drive element has not fixed variation across the X chromosome of D. pseudoobscura because, while significant linkage disequilibrium was observed within the sisterlessA, period, and esterase 5 genes, we did not find evidence for nonrandom association among loci. The Sex-ratio chromosome was estimated to be 25,000 years old based on the decomposition of linkage disequilibrium between esterase 5 and Heat-shock protein 83 or 1 million years old based on the net divergence of esterase 5 between Standard and Sex-ratio chromosomes. Genetic diversity was depressed within esterase 5 within Sex-ratio chromosomes, while the four other genes failed to show a reduction in heterozygosity in the Sex-ratio background. The reduced heterogeneity in esterase 5 is due either to its location near one of the Sex-ratio inversion breakpoints or that it is closely linked to a gene or genes responsible for the Sex-ratio meiotic drive system.     

Molecular and chromosomal polymorphism in continental and insular populations from the southwestern range of Drosophila subobscura

Eight insular and continental populations from the south-western range of Drosophila subobscura have been studied with regard to molecular and inversion polymorphisms. Heterogeneity between populations with respect to allele frequencies of 4 gene loci (Amy, Est-8, Est-9 and Pep-1) is the highest known for natural populations of the species. Moreover, the most common allele for these loci is not the same in all populations. Cladograms based on UPGMA clustering of the genetic distances based on allele frequency do not coincide with those constructed with inversion data. The allele-frequency differences between the populations may be due to non-random associations between enzyme alleles and gene arrangements and to founder effects appearing in the insular populations.