Genomics of ecological adaptation in cactophilic Drosophila.

Cactophilic Drosophila species provide a valuable model to study gene-environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii-D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii-D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

Diversity of transposable elements and repeats in a 600 kb region of the fly Calliphora vicina.

BACKGROUND: Transposable elements (TEs) are a very dynamic component of eukaryotic genomes with important implications (e.g., in evolution) and applications (e.g., as transgenic tools). They also represent a major challenge for the assembly and annotation of genomic sequences. However, they are still largely unknown in non-model species. RESULTS: Here, we have annotated the repeats and transposable elements present in a 600 kb genomic region of the blowfly Calliphora vicina (Diptera: Calliphoridae) which contains most of the achaete-scute gene complex of this species. This is the largest genomic region to be sequenced and analyzed in higher flies outside the Drosophila genus. We find that the repeat content spans at least 24% of the sequence. It includes 318 insertions classified as 3 LTR retrotransposons, 21 LINEs, 14 cut-and-paste DNA transposons, 4 helitrons and 33 unclassified repeats. CONCLUSIONS: This is the most detailed description of TEs and repeats in the Calliphoridae to date. This contribution not only adds to our knowledge about TE evolution but will also help in the annotation of repeats on Dipteran whole genome sequences.

An arthropod cis-regulatory element functioning in sensory organ precursor development dates back to the Cambrian.

BACKGROUND: An increasing number of publications demonstrate conservation of function of cis-regulatory elements without sequence similarity. In invertebrates such functional conservation has only been shown for closely related species. Here we demonstrate the existence of an ancient arthropod regulatory element that functions during the selection of neural precursors. The activity of genes of the achaete-scute (ac-sc) family endows cells with neural potential. An essential, conserved characteristic of proneural genes is their ability to restrict their own activity to single or a small number of progenitor cells from their initially broad domains of expression. This is achieved through a process called lateral inhibition. A regulatory element, the sensory organ precursor enhancer (SOPE), is required for this process. First identified in Drosophila, the SOPE contains discrete binding sites for four regulatory factors. The SOPE of the Drosophila asense gene is situated in the 5' UTR. RESULTS: Through a manual comparison of consensus binding site sequences we have been able to identify a SOPE in UTR sequences of asense-like genes in species belonging to all four arthropod groups (Crustacea, Myriapoda, Chelicerata and Insecta). The SOPEs of the spider Cupiennius salei and the insect Tribolium castaneum are shown to be functional in transgenic Drosophila. This would place the origin of this regulatory sequence as far back as the last common ancestor of the Arthropoda, that is, in the Cambrian, 550 million years ago. CONCLUSIONS: The SOPE is not detectable by inter-specific sequence comparison, raising the possibility that other ancient regulatory modules in invertebrates might have escaped detection.

The transposon Galileo generates natural chromosomal inversions in Drosophila by ectopic recombination.

BACKGROUND: Transposable elements (TEs) are responsible for the generation of chromosomal inversions in several groups of organisms. However, in Drosophila and other Dipterans, where inversions are abundant both as intraspecific polymorphisms and interspecific fixed differences, the evidence for a role of TEs is scarce. Previous work revealed that the transposon Galileo was involved in the generation of two polymorphic inversions of Drosophila buzzatii. METHODOLOGY/PRINCIPAL FINDINGS: To assess the impact of TEs in Drosophila chromosomal evolution and shed light on the mechanism involved, we isolated and sequenced the two breakpoints of another widespread polymorphic inversion from D. buzzatii, 2z(3). In the non inverted chromosome, the 2z(3) distal breakpoint was located between genes CG2046 and CG10326 whereas the proximal breakpoint lies between two novel genes that we have named Dlh and Mdp. In the inverted chromosome, the analysis of the breakpoint sequences revealed relatively large insertions (2,870-bp and 4,786-bp long) including two copies of the transposon Galileo (subfamily Newton), one at each breakpoint, plus several other TEs. The two Galileo copies: (i) are inserted in opposite orientation; (ii) present exchanged target site duplications; and (iii) are both chimeric. CONCLUSIONS/SIGNIFICANCE: Our observations provide the best evidence gathered so far for the role of TEs in the generation of Drosophila inversions. In addition, they show unequivocally that ectopic recombination is the causative mechanism. The fact that the three polymorphic D. buzzatii inversions investigated so far were generated by the same transposon family is remarkable and is conceivably due to Galileo's unusual structure and current (or recent) transpositional activity.

Evolution of the achaete-scute complex in insects: convergent duplication of proneural genes.

Proneural genes encode transcriptional activators of the basic Helix-loop-helix class that are involved in neuronal specification and differentiation. We have used the recent availability of genome sequences of multiple distant insect species to study the evolution of a family of proneural genes, the achaete-scute genes, and to examine their genomic organization and evolution. We document independent evolution of multiple copies of achaete-scute homologues and argue that this might have contributed to morphological diversity in Diptera and Lepidoptera.

HOM-C evolution in Drosophila: is there a need for Hox gene clustering?

The conservation of Homeotic (Hox) gene clustering and colinearity in many metazoans indicates that functional constraints operate on this genome organization. However, several studies have questioned its relevance in Drosophila. Here, we analyse the genomic organization of Hox and Hox-derived genes in 13 fruitfly species and the mosquito Anopheles gambiae. We found that at least seven different Homeotic complex (HOM-C) arrangements exist among Drosophila species, produced by three major splits, five microinversions and six gene transpositions. This dynamism contrasts with the stable organization of the complex in many other taxa. Although there is no evidence of an absolute requirement for Hox gene clustering in Drosophila, we found that strong functional constraints act on the individual genes.

Fast sequence evolution of Hox and Hox-derived genes in the genus Drosophila.

BACKGROUND: It is expected that genes that are expressed early in development and have a complex expression pattern are under strong purifying selection and thus evolve slowly. Hox genes fulfill these criteria and thus, should have a low evolutionary rate. However, some observations point to a completely different scenario. Hox genes are usually highly conserved inside the homeobox, but very variable outside it. RESULTS: We have measured the rates of nucleotide divergence and indel fixation of three Hox genes, labial (lab), proboscipedia (pb) and abdominal-A (abd-A), and compared them with those of three genes derived by duplication from Hox3, bicoid (bcd), zerknüllt (zen) and zerknüllt-related (zen2), and 15 non-Hox genes in sets of orthologous sequences of three species of the genus Drosophila. These rates were compared to test the hypothesis that Hox genes evolve slowly. Our results show that the evolutionary rate of Hox genes is higher than that of non-Hox genes when both amino acid differences and indels are taken into account: 43.39% of the amino acid sequence is altered in Hox genes, versus 30.97% in non-Hox genes and 64.73% in Hox-derived genes. Microsatellites scattered along the coding sequence of Hox genes explain partially, but not fully, their fast sequence evolution. CONCLUSION: These results show that Hox genes have a higher evolutionary dynamics than other developmental genes, and emphasize the need to take into account indels in addition to nucleotide substitutions in order to accurately estimate evolutionary rates.

Conservation of regulatory sequences and gene expression patterns in the disintegrating Drosophila Hox gene complex.

Homeotic (Hox) genes are usually clustered and arranged in the same order as they are expressed along the anteroposterior body axis of metazoans. The mechanistic explanation for this colinearity has been elusive, and it may well be that a single and universal cause does not exist. The Hox-gene complex (HOM-C) has been rearranged differently in several Drosophila species, producing a striking diversity of Hox gene organizations. We investigated the genomic and functional consequences of the two HOM-C splits present in Drosophila buzzatii. Firstly, we sequenced two regions of the D. buzzatii genome, one containing the genes labial and abdominal A, and another one including proboscipedia, and compared their organization with that of D. melanogaster and D. pseudoobscura in order to map precisely the two splits. Then, a plethora of conserved noncoding sequences, which are putative enhancers, were identified around the three Hox genes closer to the splits. The position and order of these enhancers are conserved, with minor exceptions, between the three Drosophila species. Finally, we analyzed the expression patterns of the same three genes in embryos and imaginal discs of four Drosophila species with different Hox-gene organizations. The results show that their expression patterns are conserved despite the HOM-C splits. We conclude that, in Drosophila, Hox-gene clustering is not an absolute requirement for proper function. Rather, the organization of Hox genes is modular, and their clustering seems the result of phylogenetic inertia more than functional necessity.

A new split of the Hox gene complex in Drosophila: relocation and evolution of the gene labial.

Hox genes encode transcription factors involved in the specification of segment identity in the early metazoan embryo. These genes are usually clustered and arranged in the same order as they are expressed along the anteroposterior body axis. This conserved genomic organization has suggested the existence of functional constraints acting on the genome organization. Partial disassembly of the Hox gene complex (HOM-C) in Caenorhabditis elegans and in two different Drosophila lineages, however, calls into question whether this cluster organization is absolutely required for proper function. Here we report a new split of the HOM-C discovered in the species of the Drosophila repleta group, which relocated the most anterior gene of the complex, lab, to a distant chromosomal site near the two most posterior Hox genes, abd-A and Abd-B. To investigate the evolutionary consequences of natural rearrangements of the Hox gene complex, the gene lab has been cloned and sequenced in D. buzzatii, a member of the D. repleta group with the split, and in D. virilis, a member of a different species group without the split. The results show that the structure of lab in D. buzzatii is intact and place the breakpoint at least 8 kb from its transcription start site. The nucleotide sequence evolution of lab in the genus Drosophila has been investigated by means of maximum likelihood methods. No significant variation has been observed among lineages in the rate of nucleotide substitution or in the nonsynonymous/synonymous substitution ratio. Seemingly, the relocation of lab has not induced a change in evolution rate or degree of functional constraint. Nevertheless, further work is needed to ascertain whether the lab-pb split has had any effects on gene expression.