Evolution of Chromatin-Remodeling Complexes: Comparative Genomics Reveals the Ancient Origin of “Novel” Compensasome Genes

Dosage compensation in Drosophila is mediated by a complex, called compensasome, composed of at least five proteins and two noncoding RNAs. Genes encoding compensasome proteins have been collectively named male-specific lethals or msls. Recent work showed that three of the Drosophila msls (msl-3, mof, and mle) have an ancient origin. In this study, I describe likely orthologues of the two remaining msls, msl-1 and msl-2, in several invertebrates and vertebrates. The MSL-2 protein is the only one found in Drosophila and vertebrate genomes that contains both a RING finger and a peculiar type of CXC domain, related to the one present in Enhancer of Zeste proteins. MSL-1 also contains two evolutionarily conserved domains: a leucine zipper and a second characteristic region, described here for the first time, which I have called the PEHE domain. These two domains are present in the likely orthologues of MSL-1 as well as in other genes in several invertebrate and vertebrate species. Although it cannot be excluded that the compensasome complex is a recent evolutionary novelty, these results shows that all msls are found in mammals, suggesting that protein complexes related to the compensasome may be present in mammalian species. Metazoans that lack several of the msls, such as Caenorhabditis elegans, cannot contain compensasomes. The evolutionary relationships of the compensasome and the NuA4 complex, another chromatin-remodeling complex that contains related subunits, are discussed.     

Drosophila sex-lethal inhibits the stable association of the 40S ribosomal subunit with msl-2 mRNA

The inhibition of male-specific lethal-2 (msl-2) mRNA translation in female flies is essential for X chromosome dosage compensation in Drosophila melanogaster. Translational repression of msl-2 requires sex-lethal (SXL) binding to uridine-rich sequences in both the 5' and 3' untranslated regions (UTRs) of the message. We delineate the msl-2 mRNA sequence elements that are important for regulation by SXL and identify functionally critical sequences adjacent to regulatory SXL binding sites. We demonstrate that SXL inhibits translation initiation and prevents the stable association of the 40S ribosomal subunit with the mRNA in a manner that does not require the presence of a cap structure at the 5' end of the mRNA. These results elucidate a critical regulatory step for dosage compensation in Drosophila melanogaster.     

Multiple tandem gene duplications in a neutral lipase gene cluster in Drosophila

We have examined a highly dynamic section of the Drosophila melanogaster genome which contains neutral lipase family genes that have undergone multiple tandem duplication events. We have identified the orthologous clusters, encoding between five and eight apparently functional lipases, in other Drosophila genomes: yakuba, ananassae, pseudoobscura, virilis, mojavensis, persimilis, grimshawi and willistoni. We examined their gene structure, duplication and pseudogene formation, and the presence of transposable elements. Based on phylogenetic comparisons, the lipase genes contained in each of the clusters fall into four distinct clades. Clades I and II have distinct evolutionary constraints to clades III and IV. Multiple gene duplications have occurred in different lineages of clades I and II while clades III and IV contain a single lipase gene from each species. Compared with lipases from other clades, clade IV genes contain an additional 3' domain of tandemly repeated sequence of varying length and composition, and a substitution in the residue adjacent to the key catalytic serine in the encoded proteins. A comparison of non-synonymous to synonymous nucleotide substitution (dN/dS) rates within each clade showed the highest rate of divergence was between paralogous lipase gene pairs suggesting selection pressure on duplicated genes. Analysis of the encoded lipase protein sequences within each species using PAML identified positively selected sites; structure homology modeling based on human pancreatic lipase indicated many of these residues formed part of the active site of the enzyme. As some of the cluster lipase genes are known to be expressed in the insect midgut and respond to changes in dietary components, we propose that the lipase cluster has undergone dynamic evolutionary changes to maximize absorption of lipid nutrients from the diet.     

A 9.6 kb intervening sequence in D. virilis rDNA, and sequence homology in rDNA interruptions of diverse species of Drosophila and other diptera.

A large proportion of the 28S ribosomal RNA genes in Drosophila virilis are interrupted by a DNA sequence 9.6 kilobase pairs long. As regards both its presence and its position in the 28S gene (about two thirds of the way in), the D. virilis rDNA intervening sequence is similar to that found in D. melanogaster rDNA, but lengths differ markedly between the two species. Degrees of nucleotide sequence homology have been detected bewteen rDNA interruptions of the two species. This homology extends to putative rDNA intervening sequences in diverse higher diptera (other Drosophila species, the house fly and the flesh fly), but hybridization of cloned D. melanogaster and D. virilis rDNA interruption segments to DNA of several lower diptera has been negative. As is the case with melanogaster rDNA interruptions, segments of the virilis rDNA intervening sequence hybridize with non-rDNA components of the virilis genome, and interspecific homology may involve these non-rDNA sequences as well as rDNA interruptions. There is, however, evidence from buoyant density fractionation of DNA that the distributions of interruption-related sequences are distinct in D. melanogaster and D. virilis genomes. Moreover, thermal denaturation studies have indicated differing extents of homology between hybridizable sequences in D. virilis DNA and different segments of the D. melanogaster rDNA intervening sequence. We infer from our studies that rDNA intervening sequences are prevalent among higher diptera; that in the course of the evolution of these organisms, elements of the intervening sequences have been moderately to highly conserved; and that this conservation extends in at least two distantly related species of Drosophila to similar sequences found elsewhere in the genomes.     

Convergent intron loss of MRP1 in Drosophila and mosquito species

Previous study revealed that the MRP1 gene ortholog DMRP1/CG6214 of Drosophila melanogaster contains 12 exons in the coding region. In the current study, the genes of DMRP1/CG6214 from D. melanogaster and Drosophila virilis were compared, and the result indicated that D. virilis had an extra intron located in exon 2, implying that intron loss or gain might have occurred at this locus. To track the evolution of the extra intron (Intron Z), orthologous nucleotide sequences of 37 arthropod species were cloned or annotated. Based on phylogenetic analysis, we found that Intron Z should present in the common ancestor of arthropod species, more than 420 Ma. In addition, we found that Sophophora subgenus species and mosquito (Culex pipiens) lost Intron Z independently, showing evolutionary convergence.     

Promoter analysis of the Drosophila melanogaster gene encoding the pterin 4alpha-carbinolamine dehydratase.

Pterin-4alpha-carbinolamine dehydratase (PCD) is a key enzyme in the regeneration pathway of tetrahydrobiopterin. Previously, we isolated and reported the Drosophila melanogaster gene encoding PCD. In the present study, we isolated and characterized the Drosophila virilis gene encoding PCD. The Drosophila virilis PCD gene has two introns and an open reading frame to encode a protein of 101 amino acids. The amino acid sequence of Drosophila virilis PCD shows a 83% homology to that of the Drosophila melanogaster PCD protein. From the alignment of the nucleotide sequence in the 5'-flanking region of the Drosophila melanogaster and Drosophila virilis PCD genes, we found four conserved sequences. Using a transient transfection assay, we showed that one of the conserved sequences (-127 to approximately -115) is critical for expression, also the minimal promoter region between -127 and +51 is necessary for the efficient expression of Drosophila melanogaster PCD.     

Nucleotide variation at the no-on-transient A gene in Drosophila littoralis

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.     

Molecular evolution of odorant-binding protein genes OS-E and OS-F in Drosophila

The Drosophila olfactory genes OS-E and OS-F are members of a family of genes that encode insect odorant-binding proteins (OBPs). OBPs are believed to transport hydrophobic odorants through the aqueous fluid within olfactory sensilla to the underlying receptor proteins. The recent discovery of a large family of olfactory receptor genes in Drosophila raises new questions about the function, diversity, regulation, and evolution of the OBP family. We have investigated the OS-E and OS-F genes in a variety of Drosophila species. These studies highlight potential regions of functional significance in the OS-E and OS-F proteins, which may include a region required for interaction with receptor proteins. Our results suggest that the two genes arose by an ancient gene duplication, and that in some lineages, one or the other gene has been lost. In D. virilis, the OS-F gene shows a different spatial pattern of expression than in D. melanogaster. One of the OS-F introns shows a striking degree of conservation between the two species, and we identify a putative regulatory sequence within this intron. Finally, a phylogenetic analysis places both OS-E and OS-F within a large family of insect OBPs and OBP-like proteins.     

Dosage Compensation Regulatory Proteins and the Evolution of Sex Chromosomes in Drosophila

In the fruitfly Drosophila melanogaster, the four male-specific lethal (msl) genes are required to achieve dosage compensation of the male X chromosome. The MSL proteins are thought to interact with cis-acting sites that confer dosage compensation to nearby genes, as they are detected at hundreds of discrete sites along the length of the polytene X chromosome in males but not in females. The histone H4 acetylated isoform, H4Ac16, colocalizes with the MSL proteins at a majority of sites on the D. melanogaster X chromosome. Using polytene chromosome immunostaining of other species from the genus Drosophila, we found that X chromosome association of MSL proteins and H4Ac16 is conserved despite differences in the sex chromosome karyotype between species. Our results support a model in which cis-acting regulatory sites for dosage compensation evolve on a neo-X chromosome arm in response to the degeneration of its former homologue.     

Expression and translocation of Drosophila nuclear encoded cytochrome b(5) proteins to mitochondria

DNA sequence studies of cytochrome b(5) (Cyt-b) genes from Drosophila melanogaster and Drosophila virilis predict that the Drosophila Cyt-b proteins are extremely hydrophobic and have at least eight potential transmembrane spanning domains. Primary protein sequence analysis also predicts that the Cyt-b proteins have mitochondrial targeting sequences and they contain sites for potential post-translational modification similar to other cytochrome proteins. We report the characterization of the cytochrome b(5) proteins from Drosophila melanogaster and Drosophila virilis. We have used a Drosophila cytochrome b(5) specific antibody to demonstrate that cytochrome b(5) proteins are expressed in muscle-containing tissues in the fly. We also provide evidence that the nuclear encoded cytochrome b(5) protein that contains a mitochondrial targeting sequence is translocated to mitochondria.