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.      

Nucleotide Polymorphism at the rp49 region of Drosophila subobscura: Lack of geographic subdivision within chromosomal arrangements in Europe

Restriction-site polymorphism at the rp49 gene region has been studied in 234 lines of Drosophila subobscura representing different gene arrangements for the O chromosome. The population surveyed (El Pedroso, Spain) was sampled four times in each of two consecutive years. The data indicate that the two chromosomal classes studied, O_[ST] and O_{[3 + 4]}, are genetically differentiated in El Pedroso. Comparison of the present results with those previously obtained for two other populations further supports that, for a given chromosomal class, European populations are not genetically differentiated. This lack of differentiation at the rp49 region within O_[ST] and within O_{[3 + 4]} stands in contrast to the clear latitudinal clines found in Europe for these arrangements.     

Tracing the colonization of Madeira and the Canary Islands by Drosophila subobscura through the study of the rp49 gene region

Nucleotide variation at the nuclear ribosomal protein 49 (rp49) gene region has been analysed by fine restriction mapping in a sample of 47 lines from a population from Madeira. Five restriction-site (out of 37 sites scored) and 3 length polymorphisms have been detected, resulting in 14 different haplotypes. This population shows less variation than both continental and Canary Island populations. The population from Madeira shows some differentiation from mainland populations, which does not favor the idea of extensive migration between the continent and Madeira. Chromosomal and restriction-map variation of the rp49 region in D. subobscura populations, together with data on sequence comparison of this nuclear region in D. guanche and D. madeirensis clearly indicate that the Canary Islands underwent at least two colonization events from the nearby continent. Although the data for Madeira are compatible with a single colonization event by a continental sample polymorphic for gene arrangements O₃ and O_{3 + 4}, an alternative scenario with at least two colonization events seems more likely.     

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.     

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.     

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.     

Inferring the demographic history of Drosophila subobscura from nucleotide variation at regions not affected by chromosomal inversions

Drosophila subobscura presents a rich and complex chromosomal inversion polymorphism. It can thus be considered a model system (i) to study the mechanisms originating inversions and how inversions affect the levels and patterns of variation in the inverted regions and (ii) to study adaptation at both the single‐gene and chromosomal inversion levels. It is therefore important to infer its demographic history as previous information indicated that its nucleotide variation is not at mutation–drift equilibrium. For that purpose, we sequenced 16 noncoding regions distributed across those parts of the J chromosome not affected by inversions in the studied population and possibly either by other selective events. The pattern of variation detected in these 16 regions is similar to that previously reported within different chromosomal arrangements, suggesting that the latter results would, thus, mainly reflect recent demographic events rather than the partial selective sweep imposed by the origin and frequency increase of inversions. Among the simple demographic models considered in our Approximate Bayesian Computation analysis of variation at the 16 regions, the model best supported by the data implies a population size expansion soon after the penultimate glacial period. This model constitutes a better null model, and it is therefore an important resource for subsequent studies aiming among others to uncover selective events across the species genome. Our results also highlight the importance of introducing the possibility of multiple hits in the coalescent simulations with an outgroup.     

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.     

Molecular population genetics of the OBP83 genomic region in Drosophila subobscura and D. guanche: contrasting the effects of natural selection and gene arrangement expansion in the patterns of nucleotide variation

Chromosomal inversion polymorphism play a major role in the evolutionary dynamics of populations and species because of their effects on the patterns of genetic variability in the genomic regions within inversions. Though there is compelling evidence for the adaptive character of chromosomal polymorphisms, the mechanisms responsible for their maintenance in natural populations is not fully understood. For this type of analysis, Drosophila subobscura is a good model species as it has a rich and extensively studied chromosomal inversion polymorphism system. Here, we examine the patterns of DNA variation in two natural populations segregating for chromosomal arrangements that differentially affect the surveyed genomic region; in particular, we analyse both nucleotide substitutions and insertion/deletion variations in the genomic region encompassing the odorant-binding protein genes Obp83a and Obp83b (Obp83 region). We show that the two main gene arrangements are genetically differentiated, but are consistent with a monophyletic origin of inversions. Nevertheless, these arrangements interchange some genetic information, likely by gene conversion. We also find that the frequency spectrum-based tests indicate that the pattern of nucleotide variation is not at equilibrium; this feature probably reflects the rapid increase in the frequency of the new gene arrangement promoted by positive selection (that is an adaptive change). Furthermore, a comparative analysis of polymorphism and divergence patterns reveals a relaxation of the functional constraints at the Obp83b gene, which might be associated with particular ecological or demographic features of the Canary island endemic species D. guanche     

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.     

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).     

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.