The High-Quality Genome Sequence of the Oceanic Island Endemic Species Drosophila guanche Reveals Signals of Adaptive Evolution in Genes Related to Flight and Genome Stability

Drosophila guanche is a member of the obscura group that originated in the Canary Islands archipelago upon its colonization by D. subobscura. It evolved into a new species in the laurisilva, a laurel forest present in wet regions that in the islands have only minor long-term weather fluctuations. Oceanic island endemic species such as D. guanche can become model species to investigate not only the relative role of drift and adaptation in speciation processes but also how population size affects nucleotide variation. Moreover, the previous identification of two satellite DNAs in D. guanche makes this species attractive for studying how centromeric DNA evolves. As a prerequisite for its establishment as a model species suitable to address all these questions, we generated a high-quality D. guanche genome sequence composed of 42 cytologically mapped scaffolds, which are assembled into six super-scaffolds (one per chromosome). The comparative analysis of the D. guanche proteome with that of twelve other Drosophila species identified 151 genes that were subject to adaptive evolution in the D. guanche lineage, with a subset of them being involved in flight and genome stability. For example, the Centromere Identifier (CID) protein, directly interacting with centromeric satellite DNA, shows signals of adaptation in this species. Both genomic analyses and FISH of the two satellites would support an ongoing replacement of centromeric satellite DNA in D. guanche.

Molecular population genetics of the Polycomb genes in Drosophila subobscura

Polycomb group (PcG) proteins are important regulatory factors that modulate the chromatin state. They form protein complexes that repress gene expression by the introduction of posttranslational histone modifications. The study of PcG proteins divergence in Drosophila revealed signals of coevolution among them and an acceleration of the nonsynonymous evolutionary rate in the lineage ancestral to the obscura group species, mainly in subunits of the Pcl-PRC2 complex. Herein, we have studied the nucleotide polymorphism of PcG genes in a natural population of D. subobscura to detect whether natural selection has also modulated the evolution of these important regulatory genes in a more recent time scale. Results show that most genes are under the action of purifying selection and present a level and pattern of polymorphism consistent with predictions of the neutral model, the exceptions being Su(z)12 and Pho. MK tests indicate an accumulation of adaptive changes in the SU(Z)12 protein during the divergence of D. subobscura and D. guanche. In contrast, the HKA test shows a deficit of polymorphism at Pho. The most likely explanation for this reduced variation is the location of this gene in the dot-like chromosome and would indicate that this chromosome also has null or very low recombination in D. subobscura, as reported in D. melanogaster.

Evidence of neofunctionalization after the duplication of the highly conserved Polycomb group gene Caf1-55 in the obscura group of Drosophila

Drosophila CAF1-55 protein is a subunit of the Polycomb repressive complex PRC2 and other protein complexes. It is a multifunctional and evolutionarily conserved protein that participates in nucleosome assembly and remodelling, as well as in the epigenetic regulation of a large set of target genes. Here, we describe and analyze the duplication of Caf1-55 in the obscura group of Drosophila. Paralogs exhibited a strong asymmetry in evolutionary rates, which suggests that they have evolved according to a neofunctionalization process. During this process, the ancestral copy has been kept under steady purifying selection to retain the ancestral function and the derived copy (Caf1-55dup) that originated via a DNA-mediated duplication event ~18 Mya, has been under clear episodic selection. Different maximum likelihood approaches confirmed the action of positive selection, in contrast to relaxed selection, on Caf1-55dup after the duplication. This adaptive process has also taken place more recently during the divergence of D. subobscura and D. guanche. The possible association of this duplication with a previously detected acceleration in the evolutionary rate of three CAF1-55 partners in PRC2 complexes is discussed. Finally, the timing and functional consequences of the Caf1-55 duplication is compared to other duplications of Polycomb genes.

Dense gene physical maps of the non-model species Drosophila subobscura

The comparative analysis of genetic and physical maps as well as of whole genome sequences had revealed that in the Drosophila genus, most structural rearrangements occurred within chromosomal elements as a result of paracentric inversions. Genome sequence comparison would seem the best method to estimate rates of chromosomal evolution, but the high-quality reference genomes required for this endeavor are still scanty. Here, we have obtained dense physical maps for Muller elements A, C, and E of Drosophila subobscura, a species with an extensively studied rich and adaptive chromosomal polymorphism. These maps are based on 462 markers: 115, 236, and 111 markers for elements A, C, and E, respectively. The availability of these dense maps will facilitate genome assembly and will thus greatly contribute to obtaining a good reference genome, which is a required step for D. subobscura to attain the model species status. The comparative analysis of these physical maps and those obtained from the D. pseudoobscura and D. melanogaster genomes allowed us to infer the number of fixed inversions and chromosomal evolutionary rates for each pairwise comparison. For all three elements, rates inferred from the more closely related species were higher than those inferred from the more distantly related species, which together with results of relative-rate tests point to an acceleration in the D. subobscura lineage at least for elements A and E.

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.

Characterization of the breakpoints of a polymorphic inversion complex detects strict and broad breakpoint reuse at the molecular level

Inversions are an integral part of structural variation within species, and they play a leading role in genome reorganization across species. Work at both the cytological and genome sequence levels has revealed heterogeneity in the distribution of inversion breakpoints, with some regions being recurrently used. Breakpoint reuse at the molecular level has mostly been assessed for fixed inversions through genome sequence comparison, and therefore rather broadly. Here, we have identified and sequenced the breakpoints of two polymorphic inversions-E1 and E2 that share a breakpoint-in the extant Est and E1 + 2 chromosomal arrangements of Drosophila subobscura. The breakpoints are two medium-sized repeated motifs that mediated the inversions by two different mechanisms: E1 via staggered breaks and subsequent repair and E2 via repeat-mediated ectopic recombination. The fine delimitation of the shared breakpoint revealed its strict reuse at the molecular level regardless of which was the intermediate arrangement. The occurrence of other rearrangements in the most proximal and distal extended breakpoint regions reveals the broad reuse of these regions. This differential degree of fragility might be related to their sharing the presence outside the inverted region of snoRNA-encoding genes.

Structure and population genetics of the breakpoints of a polymorphic inversion in Drosophila subobscura

Drosophila subobscura is a paleartic species of the obscura group with a rich chromosomal polymorphism. To further our understanding on the origin of inversions and on how they regain variation, we have identified and sequenced the two breakpoints of a polymorphic inversion of D. subobscura--inversion 3 of the O chromosome--in a population sample. The breakpoints could be identified as two rather short fragments (∼300 bp and 60 bp long) with no similarity to any known transposable element family or repetitive sequence. The presence of the ∼300-bp fragment at the two breakpoints of inverted chromosomes implies its duplication, an indication of the inversion origin via staggered double-strand breaks. Present results and previous findings support that the mode of origin of inversions is neither related to the inversion age nor species-group specific. The breakpoint regions do not consistently exhibit the lower level of variation within and stronger genetic differentiation between arrangements than more internal regions that would be expected, even in moderately small inversions, if gene conversion were greatly restricted at inversion breakpoints. Comparison of the proximal breakpoint region in species of the obscura group shows that this breakpoint lies in a small high-turnover fragment within a long collinear region (∼300 kb).

Nonrecombining genes in a recombination environment: the Drosophila "dot" chromosome

Rate of recombination is a powerful variable affecting several aspects of molecular variation and evolution. A nonrecombining portion of the genome of most Drosophila species, the "dot" chromosome or F element, exhibits very low levels of variation and unusual codon usage. One lineage of Drosophila, the willistoni/saltans groups, has the F element fused to a normally recombining E element. Here, we present polymorphism data for genes on the F element in two Drosophila willistoni and one D. insularis populations, genes previously studied in D. melanogaster. The D. willistoni populations were known to be very low in inversion polymorphism, thus minimizing the recombination suppression effect of inversions. We first confirmed, by in situ hybridization, that D. insularis has the same E + F fusion as D. willistoni, implying this was a monophyletic event. A clear gradient in codon usage exists along the willistoni F element, from the centromere distally to the fusion with E; estimates of recombination rates parallel this gradient and also indicate D. insularis has greater recombination than D. willistoni. In contrast to D. melanogaster, genes on the F element exhibit moderate levels of nucleotide polymorphism not distinguishable from two genes elsewhere in the genome. Although some linkage disequilibrium (LD) was detected between polymorphic sites within genes (generally <500 bp apart), no long-range LD between F element loci exists in the two willistoni group species. In general, the distribution of allele frequencies of F element genes display the typical pattern of expectations of neutral variation at equilibrium. These results are consistent with the hypothesis that recombination allows the accumulation of nucleotide variation as well as allows selection to act on synonymous codon usage. It is estimated that the fusion occurred ∼20 Mya and while the F element in the willistoni lineage has evolved "normal" levels and patterns of nucleotide variation, equilibrium may not have been reached for codon usage.

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

The dynamics of the roo transposable element in mutation-accumulation lines and segregating populations of Drosophila melanogaster

We estimated the number of copies for the long terminal repeat (LTR) retrotransposable element roo in a set of long-standing Drosophila melanogaster mutation-accumulation full-sib lines and in two large laboratory populations maintained with effective population size approximately 500, all of them derived from the same isogenic origin. Estimates were based on real-time quantitative PCR and in situ hybridization. Considering previous estimates of roo copy numbers obtained at earlier stages of the experiment, the results imply a strong acceleration of the insertion rate in the accumulation lines. The detected acceleration is consistent with a model where only one (maybe a few) of the approximately 70 roo copies in the ancestral isogenic genome was active and each active copy caused new insertions with a relatively high rate ( approximately 10(-2)), with new inserts being active copies themselves. In the two laboratory populations, however, a stabilized copy number or no accelerated insertion was found. Our estimate of the average deleterious viability effects per accumulated insert [E(s) < 0.003] is too small to account for the latter finding, and we discuss the mechanisms that could contain copy number.

Chromosomal evolution of elements B and C in the Sophophora subgenus of Drosophila: evolutionary rate and polymorphism

The locations of 77 markers along the chromosomal elements B (41 markers) and C (36 markers) of Drosophila subobscura, D. pseudoobscura, and D. melanogaster were obtained by in situ hybridization on polytene chromosomes. In comparisons between D. subobscura and D. pseudoobscura, 10 conserved segments (accounting for 32% of the chromosomal length) were detected on element B and eight (17% of the chromosomal length) on element C. The fixation rate of paracentric inversions inferred by a maximum likelihood approach differs significantly between elements. Muller's element C (0.17 breakpoints/Mb/million years) is evolving two times faster than element B (0.08 breakpoints/Mb/million years). This difference in the evolutionary rate is paralleled by differences in the extent of chromosomal polymorphism in the corresponding lineages. Element C is highly polymorphic in D. subobscura, D. pseudoobscura, and in other obscura group species such as D. obscura and D. athabasca. In contrast, the level of polymorphism in element B is much lower in these species. The fixation rates of paracentric inversions estimated in the present study between species of the Sophophora subgenus are the highest estimates so far reported in the genus for the autosomes. At the subgenus level, there is also a parallelism between the high fixation rate and the classical observation that the species of the Sophophora subgenus tend to be more polymorphic than the species of the Drosophila subgenus. Therefore, the detected relationship between level of polymorphism and evolutionary rate might be a general characteristic of chromosomal evolution in the genus Drosophila.

Sequences upstream of the homologous cis-elements of the Adh adult enhancer of Drosophila are required for maximal levels of Adh gene transcription in adults of Scaptodrosophila lebanonensis

The evolution of cis-regulatory elements is of particular interest for our understanding of the evolution of gene regulation. The Adh gene of Drosophilidae shows interspecific differences in tissue-specific expression and transcript levels during development. In Scaptodrosophila lebanonensis adults, the level of distal transcripts is maximal between the fourth and eighth day after eclosion and is around five times higher than that in D. melanogaster Adh(S). To examine whether these quantitative differences are regulated by sequences lying upstream of the distal promoter, we performed in vitro deletion mutagenesis of the Adh gene of S. lebanonensis, followed by P-element-mediated germ-line transformation. All constructs included, as a cotransgene, a modified Adh gene of D. melanogaster (dAdh) in a fixed position and orientation that acted as a chromosomal position control. Using this approach, we have identified a fragment of 1.5 kb in the 5' region, 830 bp upstream of the distal start site, which is required to achieve maximal levels of distal transcript in S. lebanonensis. The presence of this fragment produces a 3.5-fold higher level of distal mRNA (as determined by real time quantitative PCR) compared with the D. melanogaster dAdh cotransgene. This region contains the degenerated end of a minisatellite sequence expanding farther upstream and does not correspond to the Adh adult enhancer (AAE) of D. melanogaster. Indeed, the cis-regulatory elements of the AAE have been identified by phylogenetic footprinting within the region 830 bp upstream of the distal start site of S. lebanonensis. Furthermore, the deletions Delta-830 and Delta-2358 yield the same pattern of tissue-specific expression, indicating that all tissue-specific elements are contained within the region 830 bp upstream of the distal start site.

A minisatellite with fold-back structure is included in the 5'-flanking region of the Adh gene of Scaptodrosophila lebanonensis

A tandem repetitive sequence with a repeat unit of 12 bp has been found 1.3 kb upstream of the Adh gene of Scaptodrosophila lebanonensis. This repetitive sequence extends over 4.3 kb and consists of two inverted arrays (a fold-back segment). The repeated unit with a consensus sequence GAATACAGAATA is highly conserved and the nucleotide substitutions are not distributed randomly among the 12 bp. In situ hybridization in S. lebanonensis polytene chromosomes revealed two signals, one at the 60A section, the Adh locus, and a second site in the same chromosome at the 60C section close to the telomere. This same pattern of hybridization is obtained in all the analyzed strains including the subspecies S. lebanonensis casteeli. The minisatellite sequence accounts for about 0.03-0.04% of the S. lebanonensis genome and showed intraspecific variability in tandem repeat numbers. Possible functions of this sequence are discussed.

Evolutionary change in the structure of the regulatory region that drives tissue and temporally regulated expression of alcohol dehydrogenase gene in Drosophila funebris

The Adh locus of Drosophilidae is organized as a single gene transcribed from two spatially and temporally regulated promoters except in species of the repleta group, which have two single promoter genes. Here we show that in Drosophila funebris the Adh gene is transcribed from a single promoter, in both larva and adult, with qualitative and quantitative species specific-differences in tissue distribution. The gene is expressed in larval fat body but in other tissues such as gastric caeca, midgut and Malpighian tubules its expression is reduced compared to most Drosophilidae species, and in adults it is almost limited to the fat body. The comparative analysis of gene expression of two strains, which differ by a duplication, indicates that the cis elements necessary for this pattern of expression in larvae are included in the region of 1.55 kb upstream of the transcription initiation site. This new organization reveals the evolution of a different regulatory strategy to express the Adh gene in the subgenus Drosophila.

Chromosomal homologies between Drosophila melanogaster and D. funebris determined by in-situ hybridization

Seventeen biotin-labeled DNA sequences were hybridized to polytene chromosomes of Drosophila melanogaster and D. funebris in order to establish chromosomal homologies between these species. Ten probes correspond to cloned DNA sequences from D. melanogaster (RpII 215, MHC, H3-H4, Tor, hsp 68, hsp 28/23, hsp 83, PP1alpha, RpII 140, and ey), four are clones isolated from a D. subobscura genomic library (Xdh, lambdaDsubS3, lambdaDsubG3, and lambdaDsubG4), two are clones from D. funebris (F2 and Adh) and one from D. virilis (ci). The probes were chosen in order to cover all the autosomes, since X-chromosome homologies have already been studied by linkage analysis of morphological mutants. Most probes gave a unique hybridization signal; consequently, our results allow unambiguous inferences about chromosomal homologies. The results show extensive gene rearrangement within all chromosomal elements, probably due to paracentric inversions, but are consistent with Muller's proposal that chromosomal elements have conserved their genetic content during the evolution of Drosophila.

Fate of dot chromosome genes in Drosophila willistoni and Scaptodrosophila lebanonensis determined by in situ hybridization

One modification of the primitive karyotype of the Drosophilidae is the absence of dot chromosomes. The origin of this modification is diverse. In some cases, the fate of the dot chromosomes can be directly inferred from cytogenetic analysis but in other cases a genetic or a combined molecular and cytogenetic analysis is needed, as occurs in Drosophila willistoni and Scaptodrosophila lebanonensis. We determined the location of four dot chromosome sequences in D. willistoni and S. lebanonensis using in situ hybridization. Drosophila melanogaster and D. virilis, which possess dot chromosomes, were used as a control. The in situ hybridization results show that dot chromosome genes of D. melanogaster and D. virilis are closely linked in chromosome 3 of D. willistoni and in chromosome X of S. lebanonensis. These results suggest an autosome-dot fusion in D. willistoni and an X-dot fusion in S. lebanonensis, two different paths in the evolution of dot chromosomes.

Nucleotide sequence of the genomic region encompassing Adh and Adh-dup genes of D. lebanonensis (Scaptodrosophila): gene expression and evolutionary relationships

The region of the genome of D. lebanonensis that contains the Adh gene and the downstream Adh-dup gene was sequenced. The structure of the two genes is the same as has been described for D. melanogaster. Adh has two promoters and Adh-dup has only one putative promoter. The levels of expression of the two genes in this species are dramatically different. Hybridizing the same Northern blots with a specific probe for Adh-dup, we did not find transcripts for this gene in D. lebanonensis. The level of Adh distal transcript in adults of D. lebanonensis is five times greater than that of D. melanogaster adults. The maximum levels of proximal transcript are attained at different larval stages in the two species, being three times higher in D. melanogaster late-second-instar larvae than in D. lebanonensis first-instar larvae. The level of Adh transcripts allowed us to determine distal and proximal initiation transcription sites, the position of the first intron, the use of two polyadenylation signals, and the heterogeneity of polyadenylation sites. Temporal and spatial expression profiles of the Adh gene of D. lebanonensis show qualitative differences compared with D. melanogaster. Adh and Adh-dup evolve differently as shown by the synonymous and nonsynonymous substitution rates for the coding region of both genes when compared across two species of the melanogaster group, two of the obscura group of the subgenus Sophophora and D. lebanonensis of the victoria group of the subgenus Scaptodrsophila. Synonymous rates for Adh are approximately half those for Adh-dup, while nonsynonymous rates for Adh are generally higher than those for Adh-dup. Adh shows 76.8% identities at the protein level and 70.2% identities at the nucleotide level while Adh-dup shows 83.7% identities at the protein level and 67.5% identities at the nucleotide level. Codon usage for Adh-dup is shown to be less biased than for Adh, which could explain the higher synonymous rates and the generally lower nonsynonymous substitution rates in Adh-dup compared with Adh. Phylogenetic trees reconstructed by distance matrix and parsimony methods show that Sophophora and Scaptodrosophila subgenera diverged shortly after the separation from the Drosophila subgenus.

Chromosomal homologies between Drosophila lebanonensis and D. melanogaster determined by in situ hybridization

Twelve biotin-labelled recombinant DNA probes were hybridized to polytene chromosomes of Drosophila melanogaster and D. lebanonensis. Probes were chosen in order to cover the whole chromosomal complement. Six probes correspond to known genes from D. melanogaster (RpII215, H3-H4, MHC, hsp28/23, hsp83, hsp70), four probes are clones isolated from a D. subobscura library (Xdh, lambda DsubS3, lambdaDsubG3, lambdaDsubG4) and the remaining two probes correspond to the Adh gene of D. lebanonensis and to one sequence (262), not yet characterized, from the same species. The chromosomal homologies obtained from the in situ hybridization results allow us to determine that Muller's C and D chromosomal elements are fused in the karyotype of D. lebanonensis and constitute the large metacentric chromosome. Single pericentric inversions in the E and B elements have generated the medium and small metacentric chromosomes, respectively. No great changes are detected in Muller's A element, which remains acrocentric. The changes detected in the karyotypic evolution of D. lebanonensis are frequently observed in Drosophila evolution, as deduced from chromosomal homologies of several Drosophila species. The results are also consistent with Muller's proposal that chromosomal elements have been conserved during the evolution of Drosophila.

A photographic map of Drosophila madeirensis polytene chromosomes

A photographic map of salivary gland polytene chromosomes of Drosophila madeirensis has been constructed showing homologies and differences with respect to the standard gene arrangement of D. subobscura. Only two paracentric inversions in the X chromosome and some slight minor dissimilarities of one or two bands in the autosomes differentiate the chromosomes of these species.

Genetic analysis of extra sex combs in the hybrids between Drosophila subobscura and D. madeirensis

Drosophila subobscura and D. madeirensis are closely related species, the first distributed over a large area and the latter restricted to the island of Madeira. These species can hybridize in laboratory conditions, yielding fertile females and sterile males. Hybrid offspring show several phenotypic anomalies, including sex combs on the second and third pairs of legs in males. The extra sex comb trait is a homeotic transformation of the mesothoracic and metathoracic legs into prothoracic legs. This anomaly is observed almost exclusively in F1 males with D.madeirensis mothers. Analysis of backcross males shows that D.subobscura and D.madeirensis have diverged at a minimum of four autosomal loci affecting the extra sex comb anomaly. In addition, some incompatibility involving the X chromosome and/or a maternal effect is also implicated.

Differences in chromosome A arrangement between Drosophila madeirensis and Drosophila subobscura

The proximal half of the A (= X) chromosome of D. madeirensis has a gene arrangement very similar to the A1 or A6 inversions found in D. subobscura. Polytene chromosome analysis of hybrids between D. madeirensis and strains of D. subobscura homozygous for such inversions shows, however, that D. madeirensis has a gene arrangement different from any known for D. subobscura. These results provide evidence for a greater differentiation of the X chromosome in these species than has previously been described; it seems that the X chromosome is the only one that has undergone structural variation during the speciation process.