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.

Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions.

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.

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.

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila.

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.

Nucleotide sequence variation at two genes of the phenylpropanoid pathway, the FAH1 and F3H genes, in Arabidopsis thaliana.

The FAH1 and F3H genes encode ferulate-5-hydroxylase and flavanone-3-hydroxylase, which are enzymes in the pathways leading to the synthesis of sinapic acid esters and flavonoids, respectively. Nucleotide variation at these genes was surveyed by sequencing a sample of 20 worldwide Arabidopsis thaliana ecotypes and one Arabidopsis lyrata spp. petraea stock. In contrast with most previously studied genes, the percentage of singletons was rather low in both the FAH1 and the F3H gene regions. There was, therefore, no footprint of a recent species expansion in the pattern of nucleotide variation in these regions. In both FAH1 and F3H, nucleotide variation was structured into two major highly differentiated haplotypes. In both genes, there was a peak of silent polymorphism in the 5' part of the coding region without a parallel increase in silent divergence. In FAH1, the peak was centered at the beginning of the second exon. In F3H, nucleotide diversity was highest at the beginning of the gene. The observed pattern of variation in both FAH1 and F3H, although suggestive of balancing selection, was compatible with a neutral model with no recombination.

Nucleotide variation at the yellow gene region is not reduced in Drosophila subobscura: a study in relation to chromosomal polymorphism.

In contrast to Drosophila melanogaster and Drosophila simulans, the yellow (y) gene region of Drosophila subobscura is not located in a region with a strong reduction in recombination. In addition, this gene maps very close to the breakpoints of different inversions that segregate as polymorphic in natural populations of D. subobscura. Therefore, levels of variation at the y gene region in this species relative to those found in D. melanogaster and D. simulans may be affected not only by the change in the recombinational environment, but also by the presence of inversion polymorphism. To further investigate these aspects, an approximately 5.4-kb region of the A (=X) chromosome including the y gene was sequenced in 25 lines of D. subobscura and in the closely related species Drosophila madeirensis and Drosophila guanche. The D. subobscura lines studied differed in their A-chromosomal arrangements, A(st), A(2), and A(1). Unlike in D. melanogaster and D. simulans, levels of variation at the y gene region of D. subobscura are not reduced relative to those found at other genomic regions in the same species (rp49, Acp70A, and Acph-1). This result supports the effect of the change in the recombinational environment of a particular gene on the level of neutral variation. In addition, nucleotide variation is affected by chromosomal polymorphism. A strong genetic differentiation is detected between the A(1) arrangement and either A(st) or A(2), but not between A(st) and A(2). This result is consistent with the location of the y gene relative to the breakpoints of inversions A(1) and A(2). In addition, the pattern of nucleotide polymorphism in A(st)+A(2) and A(1) seems to point out that variation at the y gene region within these chromosomal classes is in the phase transient to equilibrium. The estimated ages of these arrangements assuming a star genealogy indicate that their origin cannot predate the D. madeirensis split. Therefore, the present results are consistent with a chromosomal phylogeny where Am(1), which is an arrangement present in D. madeirensis but absent in current populations of D. subobscura, would be the ancestral arrangement.

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.

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.

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

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.

The acid phosphatase-1 gene region in the Drosophila species of the subobscura cluster.

The acid phosphatase-1 (Acph-1) gene region was sequenced in three species of Drosophila: D. subobscura, D. madeirensis and D. guanche. These three closely related species, which are included in the obscura group, form the subobscura cluster. The different functional regions of the gene were identified by similarity with the sequence of D. melanogaster. The structure of Acph-1 is conserved in the four species. Average divergence at synonymous and nonsynonymous sites between D. melanogaster and the species of the subobscura cluster is Ks = 1.1354 and Ka = 0.1743, respectively. The rather high Ka value confirms that ACPH-1 is a rapidly evolving enzyme in Drosophila, as previously suggested by immunological studies. Amino acid replacements are not randomly distributed along the gene. In fact, an excess of replacements is detected in exon I, indicating that the signal peptide encoded by this exon evolves even faster than the rest of the protein. Divergence at the Acph-1 gene region is further evidence that D. madeirensis and D. subobscura are more closely related than D. guanche is to any of them. In addition, both silent divergence in noncoding regions and synonymous divergence in the coding region indicate that the split of the D. guanche lineage is about twice as old as the split of the lineages leading to D. madeirensis and D. subobscura. These phylogenetic relationships are, however, not supported by divergence at nonsynonymous sites since the lowest Ka estimate is between D. guanche and D. subobscura.

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.

Molecular population genetics of the rp49 gene region in different chromosomal inversions of Drosophila subobscura.

Nucleotide variation at the ribosomal protein 49 (rp49) gene region has been studied in 75 lines of Drosophila subobscura belonging to four chromosomal arrangements (Ost, O3+4, O3+4+8, and O3+4+23). The location of the rp49 gene region within the inversion loop differs among heterokaryotypes: it is very close to one of the breakpoints in heterozygotes involving Ost chromosomes, while it is in a more central position in all other heterokaryotypes. The distribution of nucleotide polymorphism in the different arrangements is consistent with a monophyletic origin of the inversions. The data also provide evidence that gene conversion and possibly double crossover are involved in shuffling nucleotide variation among gene arrangements. The analyses reveal that the level of genetic exchange is higher when the region is located in a more central position of the inverted fragment than when it is close to the breakpoints. The pairwise difference distributions as well as the negative values of Tajima's and Fu and Li's statistics further support the hypothesis that nucleotide variation within chromosomal arrangements still reflects expansion after the origin of the inversions. Under the expansion model, we have estimated the time of origin of the studied inversions.

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila.

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.

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.

An evaluation of measures of synonymous codon usage bias.

Synonymous codons are not generally used at equal frequencies, and this trend is observed for most genes and organisms. Several methods have been proposed and used to estimate the degree of the nonrandom use of the different synonymous codons. The estimates obtained by these methods, however, show different levels of both precision and dispersion when coding regions of a finite number of codons are under analysis. Here, we present a study, based on computer simulation, of how the different methods proposed to evaluate the nonrandom use of synonymous codons are affected by the length of the coding region analyzed. The results show that some of these methods are heavily influenced by the number of codons and that the comparison of codon usage bias between coding regions of different lengths shows a methodological bias under different conditions of nonrandom use of synonymous codons. The study of the dispersion of the estimates obtained by the different methods gives, on the other hand, an indication of the methods to be applied to compare values of codon usage bias among coding regions of equivalent length.

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.