Evidence for de novo evolution of testis-expressed genes in the Drosophila yakuba/Drosophila erecta clade

The mutational origin and subsequent evolution of de novo genes, which are hypothesized to be genes of recent origin that are not obviously related to ancestral coding sequence, are poorly understood. However, accumulating evidence suggests that such genes may often function in male reproduction. Here we use testis-derived expressed sequence tags (ESTs) from Drosophila yakuba to identify genes that have likely arisen either in D. yakuba or in the D. yakuba/D. erecta ancestor. We found several such genes, which show testis-biased expression and are often X-linked. Comparative data indicate that three of these genes have very short open reading frames, which suggests the possibility that a significant number of testis-biased de novo genes in the D. yakuba/D. erecta clade may be noncoding RNA genes. These data, along with previously published data from D. melanogaster, support the idea that many de novo Drosophila genes function in male reproduction and that a small region of the X chromosome in the melanogaster subgroup may be a hotspot for the evolution of novel testis-biased genes.     

Identical satellite DNA sequences in sibling species of Drosophila

The evolution of simple satellite DNAs was examined by DNA-DNA hybridization of ten Drosophila melanogaster satellite sequences to DNAs of the sibling species, Drosophila simulans and Drosophila erecta. Seven of these repeat types are present in tandem arrays in D. simulans and each of the ten sequences is repeated in D. erecta. In thermal melts, six of the seven satellite sequences in D. simulans and seven of the ten sequences in D. erecta melted within 1 deg.C of the corresponding values in D. melanogaster. The remaining sequences melted within 3 deg.C of the homologous hybrids. Therefore, there is little or no alteration in those satellite sequences held in common, despite a period of about ten million years since the divergence of D. melanogaster and D. simulans from a common ancestor. Simple satellite sequences appear to be more highly conserved than coding regions of the genome, on a per nucleotide basis. Since multiple copies of three satellite sequences could not be detected in D. simulans yet are present in D. erecta, a species more distantly related to D. melanogaster than is D. simulans, these sequences show discontinuities in evolution. There were major quantitative variations between species, showing that satellite DNAs are prone to massive amplification or diminution events over timespans as short as those separating sibling species. In D. melanogaster, these sequences amount to 21% of the genome but only 5% in D. simulans and 0.4% in D. erecta. There was a general trend of lower abundance with evolutionary distance for most satellites, suggesting that the amounts of different satellite sequences do not vary independently during evolution.     

Highly tissue specific expression of Sphinx supports its male courtship related role in Drosophila melanogaster

Sphinx is a lineage-specific non-coding RNA gene involved in regulating courtship behavior in Drosophila melanogaster. The 5' flanking region of the gene is conserved across Drosophila species, with the proximal 300 bp being conserved out to D. virilis and a further 600 bp region being conserved amongst the melanogaster subgroup (D. melanogaster, D. simulans, D. sechellia, D. yakuba, and D. erecta). Using a green fluorescence protein transformation system, we demonstrated that a 253 bp region of the highly conserved segment was sufficient to drive sphinx expression in male accessory gland. GFP signals were also observed in brain, wing hairs and leg bristles. An additional ～800 bp upstream region was able to enhance expression specifically in proboscis, suggesting the existence of enhancer elements. Using anti-GFP staining, we identified putative sphinx expression signal in the brain antennal lobe and inner antennocerebral tract, suggesting that sphinx might be involved in olfactory neuron mediated regulation of male courtship behavior. Whole genome expression profiling of the sphinx knockout mutation identified significant up-regulated gene categories related to accessory gland protein function and odor perception, suggesting sphinx might be a negative regulator of its target genes.     

Drosophila Melanogaster Mitochondrial DNA: Gene Organization and Evolutionary Considerations

The sequence of a 8351-nucleotide mitochondrial DNA (mtDNA) fragment has been obtained extending the knowledge of the Drosophila melanogaster mitochondrial genome to 90% of its coding region. The sequence encodes seven polypeptides, 12 tRNAs and the 3' end of the 16S rRNA and CO III genes. The gene organization is strictly conserved with respect to the Drosophila yakuba mitochondrial genome, and different from that found in mammals and Xenopus. The high A + T content of D. melanogaster mitochondrial DNA is reflected in a reiterative codon usage, with more than 90% of the codons ending in T or A, G + C rich codons being practically absent. The average level of homology between the D. melanogaster and D. yakuba sequences is very high (roughly 94%), although insertion and deletions have been detected in protein, tRNA and large ribosomal genes. The analysis of nucleotide changes reveals a similar frequency for transitions and transversions, and reflects a strong bias against G + C on both strands. The predominant type of transition is strand specific.     

Sequences sufficient for correct regulation of Sgs-3 lie close to or within the gene.

The Drosophila melanogaster 68C chromosomal locus is the site of a prominent polytene chromosome puff that harbors the genes Sgs-3, Sgs-7 and Sgs-8. These genes code for proteins that are part of the salivary glue that Drosophila larvae secrete as a means of fixing themselves to an external substrate for the duration of the pre-pupal and pupal period. The 68C glue genes are regulated by the steroid hormone ecdysterone, with the hormone required for both initiation and cessation of gene expression during the third larval instar. Previous work has defined sequences sufficient for expression of abundant levels of Sgs-3 mRNA at the correct time and in the correct tissue. We show here that sequences sufficient for normal tissue- and stage-specific accumulation of Sgs-3 RNA, but adequate only for low levels of expression, lie within 130 bp of the 5' end of the gene, or within the gene.     

Divergence and heterogeneity of the histone gene repeating units in the Drosophila melanogaster species subgroup

The repeating units of the histone gene cluster containing the H1, H2A, H2B and H4 genes were amplified by PCR from the Drosophila melanogaster species subgroup, i.e., D. yakuba, D. erecta, D. sechellia, D. mauritiana, D. teissieri and D. orena. The PCR products were cloned and their nucleotide sequences of about 4.6-4.8kbp were determined to elucidate the mechanism of molecular evolution of the histone gene family. The heterogeneity among the histone gene repeating units was 0.6% and 0.7% for D. yakuba and D. sechellia, respectively, indicating the same level of heterogeneity as in the H3 gene region of D. melanogaster. Divergence of the genes among species even in the most closely related ones was much greater than the heterogeneity among family members, indicating a concerted mode of evolution for the histone gene repeating units. Among the species in the D. melanogaster species subgroup, the histone gene regions as well as 3rd codon position of the coding region showed nearly the same GC contents. These results suggested that the previous conclusion on analysis of the H3 gene regions, the gene family evolution in a concerted fashion, holds true for the whole histone gene repeating unit.     

Cooperative enhancement at the Drosophila Sgs-3 locus

The Drosophila glue gene Sgs-3 is specifically expressed in the secretory cells of the salivary glands of third instar larvae. We have assayed the expression of gene fusions to determine the role of cis-acting Sgs-3 sequences in conferring this pattern of expression. These experiments define two regulatory regions required for expression of reporter genes from the Sgs-3 promoter. One region, between 106 and 56 bp upstream of the Sgs-3 mRNA 5' end is sufficient for low but correct tissue- and stage-specific expression. A second region, lying between 629 and 130 bp 5' of the RNA start site is functionally equivalent; that is, it alone will also direct low level, specific expression. These two regions act synergistically to give high level expression. More distant upstream regions function to further increase levels of expression. These two regulatory elements can confer a salivary gland-specific pattern of expression on a heterologous promoter and are also sufficient to drive gene expression in other Drosophila species, implying conservation of regulators.     

Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements

Horizontal transmission has been well documented as a major mechanism for the dissemination of mariner-like elements (MLEs) among species. Less well understood are mechanisms that limit vertical transmission of MLEs resulting in the "spotty" or discontinuous distribution observed in closely related species. In this article we present evidence that the genome of the common ancestor of the melanogaster species subgroup of Drosophila contained an MLE related to the mellifera (honey bee) subfamily. Horizontal transmission, approximately 3-10 MYA, is strongly suggested by the observation that the sequence of the MLE in Drosophila erecta is 97% identical in nucleotide sequence with that of an MLE in the cat flea, Ctenocephalides felis. The D. erecta MLE has a spotty distribution among species in the melanogaster subgroup. The element has a high copy number in D. erecta and D. orena, a moderate copy number in D. teissieri and D. yakuba, and was apparently lost ("stochastic loss") in the lineage leading to D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In D. erecta, most copies are concentrated in the heterochromatin. Two copies from D. erecta, denoted De12 and De19, were cloned and sequenced, and they appear to be nonfunctional ("vertical inactivation"). It therefore appears that the predominant mode of MLE evolution is vertical inactivation and stochastic loss balanced against occasional reinvasion of lineages by horizontal transmission.      

The hemoglobin genes of Drosophila

We recently reported the unprecedented occurrence of a hemoglobin gene (glob1) in the fruitfly Drosophila melanogaster. Here we investigate the structure and evolution of the glob1 gene in other Drosophila species. We cloned and sequenced glob1 genes and cDNA from D. pseudoobscura and D. virilis, and identified the glob1 gene sequences of D. simulans, D. yakuba, D. erecta, D. ananassae, D. mojavensis and D. grimshawi in the databases. Gene structure (introns in helix positions D7.0 and G7.0), gene synteny and sequence of glob1 are highly conserved, with high ds/dn ratios indicating strong purifying selection. The data suggest an important role of the glob1 protein in Drosophila, which may be the control of oxygen flow from the tracheal system. Furthermore, we identified two additional globin genes (glob2 and glob3) in the Drosophilidae. Although the sequences are highly derived, the amino acids required for heme- and oxygen-binding are conserved. In contrast to other known insect globin, the glob2 and glob3 genes harbour both globin-typical introns at positions B12.2 and G7.0. Both genes are conserved in various drosophilid species, but only expression of glob2 could be demonstrated by western blotting and RT-PCR. Phylogenetic analyses show that the clade leading to glob2 and glob3, which are sistergroups, diverged first in the evolution of dipteran globins. glob1 is closely related to the intracellular hemoglobin of the botfly Gasterophilus intestinalis, and the extracellular hemoglobins from the chironomid midges derive from this clade.     

cis-acting sequences required for expression of the divergently transcribed Drosophila melanogaster Sgs-7 and Sgs-8 glue protein genes.

The Sgs-7 and Sgs-8 glue genes at 68C are divergently transcribed and are separated by 475 bp. Fusion genes with Adh or lacZ coding sequences were constructed, and the expression of these genes, with different amounts of upstream sequences present, was tested by a transient expression procedure and by germ line transformation. A cis-acting element for both genes is located asymmetrically in the intergenic region between -211 and -43 bp relative to Sgs-7. It is required for correct expression of both genes. This element can confer the stage- and tissue-specific expression pattern of glue genes on a heterologous promoter. An 86-bp portion of the element, from -133 to -48 bp relative to Sgs-7, is shown to be capable of enhancing the expression of a truncated and therefore weakly expressed Sgs-3 fusion gene. Recently described common sequence motifs of glue gene regulatory elements (T. Todo, M. Roark, K. Vijay Raghavan, C. A. Mayeda, and E.M. Meyerowitz, Mol. Cell. Biol. 10:5991-6002, 1990) are located within this 86-bp region.     

A screen for immunity genes evolving under positive selection in Drosophila

Genes involved in the immune system tend to have higher rates of adaptive evolution than other genes in the genome, probably because they are coevolving with pathogens. We have screened a sample of Drosophila genes to identify those evolving under positive selection. First, we identified rapidly evolving immunity genes by comparing 140 loci in Drosophila erecta and D. yakuba. Secondly, we resequenced 23 of the fastest evolving genes from the independent species pair D. melanogaster and D. simulans, and identified those under positive selection using a McDonald-Kreitman test. There was strong evidence of adaptive evolution in two serine proteases (persephone and spirit) and a homolog of the Anopheles serpin SRPN6, and weaker evidence in another serine protease and the death domain protein dFADD. These results add to mounting evidence that immune signalling pathway molecules often evolve rapidly, possibly because they are sites of host-parasite coevolution.