Adaptive selection and coevolution at the proteins of the Polycomb repressive complexes in Drosophila.

Polycomb group (PcG) proteins are important epigenetic regulatory proteins that modulate the chromatin state through posttranslational histone modifications. These interacting proteins form multimeric complexes that repress gene expression. Thus, PcG proteins are expected to evolve coordinately, which might be reflected in their phylogenetic trees by concordant episodes of positive selection and by a correlation in evolutionary rates. In order to detect these signals of coevolution, the molecular evolution of 17 genes encoding the subunits of five Polycomb repressive complexes has been analyzed in the Drosophila genus. The observed distribution of divergence differs substantially among and along proteins. Indeed, CAF1 is uniformly conserved, whereas only the established protein domains are conserved in other proteins, such as PHO, PHOL, PSC, PH-P and ASX. Moreover, regions with a low divergence not yet described as protein domains are present, for instance, in SFMBT and SU(Z)12. Maximum likelihood methods indicate an acceleration in the nonsynonymous substitution rate at the lineage ancestral to the obscura group species in most genes encoding subunits of the Pcl-PRC2 complex and in genes Sfmbt, Psc and Kdm2. These methods also allow inferring the action of positive selection in this lineage at genes E(z) and Sfmbt. Finally, the protein interaction network predicted from the complete proteomes of 12 Drosophila species using a coevolutionary approach shows two tight PcG clusters. These clusters include well-established binary interactions among PcG proteins as well as new putative interactions.

A molecular perspective on a complex polymorphic inversion system with cytological evidence of multiply reused breakpoints.

Genome sequence comparison across the Drosophila genus revealed that some fixed inversion breakpoints had been multiply reused at this long timescale. Cytological studies of Drosophila inversion polymorphism had previously shown that, also at this shorter timescale, some breakpoints had been multiply reused. The paucity of molecularly characterized polymorphic inversion breakpoints has so far precluded contrasting whether cytologically shared breakpoints of these relatively young inversions are actually reused at the molecular level. The E chromosome of Drosophila subobscura stands out because it presents several inversion complexes. This is the case of the E1+2+9+3 arrangement that originated from the ancestral Est arrangement through the sequential accumulation of four inversions (E1, E2, E9 and E3) sharing some breakpoints. We recently identified the breakpoints of inversions E1 and E2, which allowed establishing reuse at the molecular level of the cytologically shared breakpoint of these inversions. Here, we identified and sequenced the breakpoints of inversions E9 and E3, because they share breakpoints at sections 58D and 64C with those of inversions E1 and E2. This has allowed establishing that E9 and E3 originated through the staggered-break mechanism. Most importantly, sequence comparison has revealed the multiple reuse at the molecular level of the proximal breakpoint (section 58D), which would have been used at least by inversions E2, E9 and E3. In contrast, the distal breakpoint (section 64C) might have been only reused once by inversions E1 and E2, because the distal E3 breakpoint is displaced >70 kb from the other breakpoint limits.

Multilocus analysis of nucleotide variation in Drosophila madeirensis, an endemic species of the Laurisilva forest in Madeira.

Drosophila madeirensis is an endemic species of Madeira that inhabits the island Laurisilva forest. Nucleotide variation in D. madeirensis is analysed in six genomic regions and compared to that previously reported for the same regions in Drosophila subobscura, an abundant species in the Palearctic region that is closely related to D. madeirensis. The gene regions analysed are distributed along the O(3) inversion. The O(3) arrangement is monomorphic in D. madeirensis, and it was present in ancestral populations of D. subobscura but went extinct in this species after the origin of the derived O(ST) and O(3+4) arrangements. Levels of nucleotide polymorphism in D. madeirensis are similar to those present in the O(ST) and O(3+4) arrangements of D. subobscura, and the frequency spectrum is skewed towards rare variants. Purifying selection against deleterious nonsynonymous mutations is less effective in D. madeirensis. Although D. madeirensis and D. subobscura coexist at present in Madeira, no clear evidence of introgression was detected in the studied regions.

Multilocus analysis of variation using a large empirical data set: phenylpropanoid pathway genes in Arabidopsis thaliana.

Detecting the signature of adaptation on nucleotide variation is often difficult in species that like Arabidopsis thaliana might have a complex demographic history. Recent re-sequencing surveys in this species provided genome-wide information that would mainly reflect its demographic history. We have used a large empirical data set (LED) as well as multilocus coalescent simulations to analyse sequence variation at loci involved in the phenylpropanoid pathway of this species. We surveyed and examined DNA sequence variation at nine of these loci (about 19.7 kb) in 23 accessions of A. thaliana and one accession of its closely related species Arabidopsis lyrata. Nucleotide variation was lower at nonsynonymous sites than at silent sites in all loci, indicating generalized functional constraint at the protein level. No association between variation and position in the metabolic pathway was detected. When the data were contrasted against the standard neutral model, significant deviations for silent variation were detected with Tajima's D, Fu's F(S) and Fay and Wu's H multilocus test statistics. These deviations were in the same direction than in previous large-scale multilocus analyses, suggesting a genome-wide effect. When the nine-locus data set was contrasted against the large empirical data set, the level (Watterson's theta) and pattern of variation (Tajima's D) detected in these loci did not deviate either at the single-locus or multilocus level from the corresponding empirical distributions. These results would support an important role of the demographic history of A. thaliana in shaping nucleotide variation at the nine studied phenylpropanoid loci. The potential and limitations of the empirical distribution approach are discussed.

Restriction map variation at the adh locus of Drosophila melanogaster in inverted and noninverted chromosomes.

Restriction map variation among 39 Standard and 40 In(2L)t chromosomes extracted from a Spanish natural population of Drosophila melanogaster was investigated for a 2.7-kb region encompassing the Adh locus with ten four-cutter restriction enzymes. A total of 20 polymorphisms were detected, representing 15 restriction site polymorphisms, 4 length polymorphisms and the allozyme polymorphism. Variation at the DNA level was compared among St-Adh(F), St-Adh(S) and t-Adh(S) chromosomes. t-Adh(S) chromosomes show a higher level of variation than St-Adh(F) chromosomes. This suggests that In(2L)t arose before the fast/slow allozyme divergence in the evolutionary history of D. melanogaster.

Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

Nucleotide sequence variation at the Acp29AB gene region has been surveyed in Drosophila melanogaster from Spain (12 lines), Ivory Coast (14 lines), and Malawi (13 lines) and in one line of D. simulans. The approximately 1.7-kb region studied encompasses the Acp29AB gene that codes for a male accessory gland protein and its flanking regions. Seventy-seven nucleotide and 8 length polymorphisms were detected. Nonsynonymous polymorphism was an order of magnitude lower than synonymous polymorphism, but still high relative to other non-sex-related genes. In D. melanogaster variation at this region revealed no major genetic differentiation between East and West African populations, while differentiation was highly significant between the European and the two African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites showed an excess of fixed nonsynonymous changes, which indicates that the evolution of the Acp29AB protein has been driven by directional selection at least after the split of the D. melanogaster and D. simulans lineages. The pattern of variation in extant populations of D. melanogaster favors a scenario where the fixation of advantageous replacement substitutions occurred in the early stages of speciation and balancing selection is maintaining variation in this species.

Different forces drive the evolution of the Acp26Aa and Acp26Ab accessory gland genes in the Drosophila melanogaster species complex.

The Acp26Aa and Acp26Ab genes that code for male accessory gland proteins are tandemly arranged in the species of the Drosophila melanogaster complex. An approximately 1.6-kb region encompassing both genes has been sequenced in 10, 24, and 18 lines from Spain, Ivory Coast, and Malawi, respectively; the previously studied 10 lines from North Carolina have also been included in the analyses. A total of 110 nucleotide and 4 length polymorphisms were detected. Silent variation for the whole Acp26A region was slightly higher in African than in non-African populations, while for both genes nonsynonymous variation was similar in all populations studied. Based on Fst estimates no major genetic differentiation was detected between East and West Africa, while in general non-African populations were strongly differentiated from both African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites revealed that directional selection acting on amino acid replacement changes has driven the evolution of the Acp26Aa protein in the last 2.5 myr.

Molecular organization of the X chromosome in different species of the obscura group of Drosophila.

Nine single copy regions located on the X chromosome have been mapped by in situ hybridization in six species of the obscura group of Drosophila. Three Palearctic species, D. subobscura, D. madeirensis and D. guanche, and three Nearctic species, D. pseudoobscura, D. persimilis and D. miranda, have been studied. Eight of the regions include known genes from D. melanogaster (Pgd, zeste, white, cut, vermilion, RNA polymerase II 215, forked and suppressor of forked) and the ninth region (lambda DsubF6) has not yet been characterized. In all six species, as in D. melanogaster, all probes hybridize to a single site. Established chromosomal arm homologies of Muller's element A are only partly supported by present results since two of the probes (Pgd and zeste) hybridize at the proximal end of the XR chromosomal arm in the three Nearctic species. In addition to the centric fusion of Muller's A (= XL) and D (= XR) elements, the metacentric X chromosome of the Nearctic species requires a pericentric inversion to account for this result. Previously proposed homologies of particular chromosomal regions of the A (= X) chromosome in the three species of the D. subobscura cluster and of the XL chromosomal arm in the three species of the D. pseudoobscura cluster are discussed in light of the present results. Location of the studied markers has changed drastically not only since the divergence between the melanogaster and obscura groups but also since the Palearctic and Nearctic species of the obscura group diverged.

Evolutionary history of the sex-peptide (Acp70A) gene region in Drosophila melanogaster.

In Drosophila the products of the seminal fluid stimulate oviposition and suppress remating in the female. Of all the accessory gland peptides (Acp's) involved in these two responses, the sex-peptide (coded by the Acp70A gene) is among the best characterized at the functional level. A 1.2-kb fragment encompassing the Acp70A gene of nine lines from a natural population of D. melanogaster and one allele of D. sechellia was sequenced to study the forces shaping nucleotide variation within and between species. The coding region of D. simulans and D. mauritiana was also sequenced. A Ser to Ala replacement polymorphism at the last position of the signal peptide was detected in D. melanogaster. The Ser and Ala alleles are at intermediate frequencies. The level of nucleotide variation is lower for the derived Ala allele, which is compatible with a recent origin and an increase in frequency due to positive selection. Variation at the 5' flanking region is structured in two major highly differentiated haplotypes, whose distribution does not conform to neutral expectations. Selective and/or historical factors could contribute to the observed overall patterning of nucleotide variation at the Acp70A region.

Nucleotide variation at the CHALCONE ISOMERASE locus in Arabidopsis thaliana.

An approximately 1.9-kb region encompassing the CHI gene, which encodes chalcone isomerase, was sequenced in 24 worldwide ecotypes of Arabidopsis thaliana (L.) Heynh. and in 1 ecotype of A. lyrata ssp. petraea. There was no evidence for dimorphism at the CHI region. A minimum of three recombination events was inferred in the history of the sampled ecotypes of the highly selfing A. thaliana. The estimated nucleotide diversity theta(TOTAL) = 0.004, theta(SIL) = 0. 005 was on the lower part of the range of the corresponding estimates for other gene regions. The skewness of the frequency spectrum toward an excess of low-frequency polymorphisms, together with the bell-shaped distribution of pairwise nucleotide differences at CHI, suggests that A. thaliana has recently experienced a rapid population growth. Although this pattern could also be explained by a recent selective sweep at the studied region, results from the other studied loci and from an AFLP survey seem to support the expansion hypothesis. Comparison of silent polymorphism and divergence at the CHI region and at the Adh1 and ChiA revealed in some cases a significant deviation of the direct relationship predicted by the neutral theory, which would be compatible with balancing selection acting at the latter regions.

Nucleotide polymorphism at the RpII215 gene in Drosophila subobscura. Weak selection on synonymous mutations.

Nucleotide variation in an 8.1-kb fragment encompassing the RpII215 gene, which encodes the largest subunit of the RNA polymerase II complex, is analyzed in a sample of 11 chromosomes from a natural population of Drosophila subobscura. No amino acid polymorphism was detected among the 157 segregating sites. The observed numbers of preferred and unpreferred derived synonymous mutations can be explained by neutral mutational processes. In contrast, preferred mutations segregate at significantly higher frequency than unpreferred mutations, suggesting the action of natural selection. The polymorphism to divergence ratio is different for preferred and unpreferred changes, in agreement with their beneficial and deleterious effects on fitness, respectively. Preferred and unpreferred codons are nonrandomly distributed in the RpII215 gene, leading to a heterogeneous distribution of polymorphic to fixed synonymous differences across this coding region. This intragenic variation of the polymorphism/divergence ratio cannot be explained by different patterns of gene expression, mutation, or recombination rates, and therefore it indicates that selection coefficients for synonymous mutations can vary extensively across a coding region. The application of nucleotide composition stationarity tests in coding and flanking noncoding regions, assumed to behave neutrally, allows the detection of the action of natural selection when stationarity holds in the noncoding region.

Evidence of extensive genetic exchange in the rp49 region among polymorphic chromosome inversions in Drosophila subobscura.

Restriction map variation in 107 lines extracted from two natural populations of Drosophila subobscura was investigated with seven four-nucleotide-recognizing enzymes in a 1.6-kb region including the rp49 gene, that is located very close to the proximal breakpoint of inversion O3. Fourteen restriction site and 8 length polymorphisms, resulting in 73 haplotypes, have been identified. Estimated heterozygosity per nucleotide, pi = 0.0045, is comparable to the average nucleotide variation observed in Drosophila melanogaster. Because of the location of the rp49 region in D. subobscura, variation in three different gene arrangements-Ost, O3 + 4 and O3 + 4 + 8-has been compared. Out of 14 restriction site polymorphisms, 3 are shared by Ost, O3 + 4 and O3 + 4 + 8, and 3 additional ones are shared by Ost and O3 + 4, evidencing extensive genetic exchange among these polymorphic inversions. In agreement with previous data, the higher level of variation of O3 + 4 (as measured by haplotype diversity and nucleotide variation) suggests that O3 + 4 may be ancestral in relationship to extant gene arrangements.

Synonymous rates at the RpII215 gene of Drosophila: variation among species and across the coding region.

The region encompassing the RpII215 gene that encodes the largest component of the RNA polymerase II complex (1889 amino acids) has been sequenced in Drosophila subobscura, D. madeirensis, D. guanche, and D. pseudoobscura. Nonsynonymous divergence estimates (Ka) indicate that this gene has a very low rate of amino acid replacements. Given its low Ka and constitutive expression, synonymous substitution rates are, however, unexpectedly high. Sequence comparisons have allowed the molecular clock hypothesis to be tested. D. guanche is an insular species and it is therefore expected to have a reduced effective size relative to D. subobscura. The significantly higher rate of synonymous substitutions detected in the D. guanche lineage could be explained if synonymous mutations behave as nearly neutral. Significant departure from the molecular clock hypothesis for synonymous and nonsynonymous substitutions was detected when comparing the D. subobscura, D. pseudoobscura, and D. melanogaster lineages. Codon bias and synonymous divergence between D. subobscura and D. melanogaster were negatively correlated across the RpII215 coding region, which indicates that selection coefficients for synonymous mutations vary across the gene. The C-terminal domain (CTD) of the RpII215 protein is structurally and functionally differentiated from the rest of the protein. Synonymous substitution rates were significantly different in both regions, which strongly indicates that synonymous mutations in the CTD and in the non-CTD regions are under detectably different selection coefficients.

Excess polymorphism at the Adh locus in Drosophila melanogaster.

The evolutionary history of a region of DNA encompassing the Adh locus is studied by comparing patterns of variation in Drosophila melanogaster and its sibling species, D. simulans. An unexpectedly high level of silent polymorphism in the Adh coding region relative to the 5' and 3' flanking regions in D. melanogaster is revealed by a populational survey of restriction polymorphism using a four-cutter filter hybridization technique as well as by direct sequence comparisons. In both of these studies, a region of the Adh gene encompassing the three coding exons exhibits a frequency of polymorphism equal to that of a 4-kb 5' flanking region. In contrast, an interspecific sequence comparison shows a two-fold higher level of divergence in the 5' flanking sequence compared to the structural locus. Analysis of the patterns of variation suggest an excess of polymorphism within the D. melanogaster Adh locus, rather than lack of polymorphism in the 5' flanking region. An approach is outlined for testing neutral theory predictions about patterns of variation within and between species. This approach indicates that the observed patterns of variation are incompatible with an infinite site neutral model.

Divergence of the yellow gene between Drosophila melanogaster and D. subobscura: recombination rate, codon bias and synonymous substitutions.

The yellow (y) gene maps near the telomere of the X chromosome in Drosophila melanogaster but not in D. subobscura. Thus the strong reduction in the recombination rate associated with telomeric regions is not expected in D. subobscura. To study the divergence of a gene whose recombination rate differs between two species, the y gene of D. subobscura was sequenced. Sequence comparison between D. melanogaster and D. subobscura revealed several elements conserved in noncoding regions that may correspond to putative cis-acting regulatory sequences. Divergence in the y gene coding region between D. subobscura and D. melanogaster was compared with that found in other genes sequenced in both species. Both, yellow and scute exhibit an unusually high number of synonymous substitutions per site (ps). Also for these genes, the extent of codon bias differs between both species, being much higher in D. subobscura than in D. melanogaster. This pattern of divergence is consistent with the hitchhiking and background selection models that predict an increase in the fixation rate of slightly deleterious mutations and a decrease in the rate of fixation of slightly advantageous mutations in regions with low recombination rates such as in the y-sc gene region of D. melanogaster.

Molecular evolution of the Cecropin multigene family in Drosophila. functional genes vs. pseudogenes.

Approximately 4 kb of the Cecropin cluster region have been sequenced in nine lines of Drosophila melanogaster and one line of the sibling species D. simulans, D. mauritiana, and D. sechellia. This region includes three functional genes (CecA1, CecA2, and CecB), which are involved in the insect immune response, and two pseudogenes (CecPsi1 and CecPsi2). The level of silent polymorphism in the three Cec genes is rather high (0.028), and there is no excess of nonsynonymous polymorphism. There is no evidence of gene conversion in the history of these genes. The interspecific comparison has revealed that in the three species of the simulans cluster the CecA2 gene is partially deleted and has therefore lost its function and become a pseudogene; in each of the species, subsequent deletions have accumulated. Divergence estimates indicate that the CecPsi1 and CecPsi2 pseudogenes are highly diverged, both between themselves and relative to the other three Cec genes. However, both CecPsi1 and CecPsi2 have conserved transcriptional signals and splice sites, and they present an open reading frame; also, correctly spliced transcripts have been detected for both CecPsi1 and CecPsi2. The data support that these genes are either active genes with some null alleles or young pseudogenes.

Differentiation of Muller's chromosomal elements D and E in the obscura group of Drosophila.

Twenty-two markers located on Muller's elements D or E have been mapped by in situ hybridization in six species of the obscura group of Drosophila and in D. melanogaster. The obscura species can be grouped into a Palearctic cluster (D. subobscura, D. madeirensis and D. guanche) and a Nearctic one (D. pseudoobscura, D. persimilis and D. miranda). Eleven of the probes contain known genes: E74, Acp70A, Est5, hsp28/23, hsp83, emc, hsp70, Xdh, Acph-1, Cec and rp49. The remaining probes are recombinant phages isolated from a D. subobscura genomic library. All these markers hybridize to the putative homologous chromosome or chromosomal arm of elements D and E. Thus, these elements have conserved their genic content during species divergence. Chromosomal homologies proposed previously for each element among the species of the same cluster have been compared with the present results. The distribution of markers within each element has changed considerably as inferred from pairwise comparisons of obscura species included in the two different clusters. Only chromosomal segments defined by closely linked markers have been conserved: one such segment has been detected in element D and three in element E between D subobscura and D. pseudoobscura.

Synonymous substitutions in the Xdh gene of Drosophila: heterogeneous distribution along the coding region.

The Xdh (rosy) region of Drosophila subobscura has been sequenced and compared to the homologous region of D. pseudoobscura and D. melanogaster. Estimates of the numbers of synonymous substitutions per site (Ks) confirm that Xdh has a high synonymous substitution rate. The distributions of both nonsynonymous and synonymous substitutions along the coding region were found to be heterogeneous. Also, no relationship has been detected between Ks estimates and codon usage bias along the gene, in contrast with the generally observed relationship among genes. This heterogeneous distribution of synonymous substitutions along the Xdh gene, which is expression-level independent, could be explained by a differential selection pressure on synonymous sites along the coding region acting on mRNA secondary structure. The synonymous rate in the Xdh coding region is lower in the D. subobscura than in the D. pseudoobscura lineage, whereas the reverse is true for the Adh gene.

A test of neutral molecular evolution based on nucleotide data.

The neutral theory of molecular evolution predicts that regions of the genome that evolve at high rates, as revealed by interspecific DNA sequence comparisons, will also exhibit high levels of polymorphism within species. We present here a conservative statistical test of this prediction based on a constant-rate neutral model. The test requires data from an interspecific comparison of at least two regions of the genome and data on levels of intraspecific polymorphism in the same regions from at least one species. The model is rejected for data from the region encompassing the Adh locus and the 5' flanking sequence of Drosophila melanogaster and Drosophila sechellia. The data depart from the model in a direction that is consistent with the presence of balanced polymorphism in the coding region.

Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster.

Restriction map variation in 64 X chromosome lines extracted from three different populations of Drosophila melanogaster was investigated with seven six-nucleotide-recognizing restriction enzymes for a 106-kb region encompassing the yellow gene and the achaete-scute complex that is located in the region of reduced crossing over close to the telomere. Nine restriction site polymorphisms (out of 176 sites scored) and 19 length polymorphisms (15 insertions and 4 deletions) were detected. The estimated level of heterozygosity per nucleotide, H = 0.0003, is much lower than that reported for autosomal and sex-linked loci located in regions with normal levels of crossing over. The overall frequency of polymorphic restriction sites is reduced. Six out of nine restriction site polymorphisms are unique and the other three have frequencies less than 0.17. Some large insertions have reached relatively high frequencies, 0.08 to 0.17. Consistent with the theoretically predicted negative relationship between crossing over and the magnitude of linkage disequilibrium, an increase in the relative number of nonrandom associations was observed in the y-ac-sc region.