Evolution of the Larval Cuticle Proteins Coded by the Secondary Sex Chromosome Pair: X2 and Neo-Y of Drosophila miranda: II. Comparison at the Amino Acid Sequence Level

The larval cuticle proteins (LCPs) are encoded by a multigene family, Lcp1–4, located at the right arm of the metacentric autosome 2 (2R) in Drosophila melanogaster. Due to a chromosome fusion the Lcp locus of Drosophila miranda is situated on a pair of secondary sex chromosomes, the X2 and neo-Y chromosomes. Comparing the deduced amino acid sequences of the autosomal D. melanogaster loci with the sex-chromosomal loci of D. miranda, we were able to trace the evolution of the Lcp loci with respect to their different chromosomal inheritance. The length of the signal peptide is conserved in all four LCPs, while the size of the mature LCPs varies. Conserved protein motifs became obvious from the alignment, indicating regions of structural and functional importance. Analyzing intra- and interspecific sequence similarities of the Lcp gene families allowed us to reconstruct the phylogeny of the gene cluster. Alignment with cuticle amino acid sequences originating from divergent insect species reveals motifs already present in the primordial insect LCPs. These motifs indicate different levels of constraint acting during the evolution of the LCPs. Received: 27 December 1995 / Accepted: 30 April 1996     

Evolutionary changes in the organization of the major LCP gene cluster during sex chromosomal differentiation in the sibling speciesDrosophila persimilis,D. pseudoobscura andD. miranda

The major larval cuticle protein (LCP) genes I–IV ofDrosophila melanogaster are clustered on the right arm of the second chromosome. By cross-hybridization we cloned the corresponding genes from three different members of theobscura group:D. persimilis, D. pseudoobscura andD. miranda. InD. pseudoobscura andD. persimilis the gene cluster maps to autosome3. In contrast, inD. miranda it was found on theX2 andY sex chromosome. Hence, this exceptional karyotypic situation offers a unique opportunity to analyse the molecular processes underlying the phenomenon of chromosome degeneration. Comparison of LCP genes I–IV in theX2 andY chromosomal region inD. miranda revealed extensive DNA rearrangements at the latter. TheY chromosomal LCP cluster is characterized by DNA insertions which are absent in the correspondingX2 chromosomal DNA, suggesting that these DNA sequences must have invaded this area. In addition, part of the analysedY chromosomal region is duplicated.     

Absence of Protein Polymorphism in the Ras Genes of Drosophila melanogaster

Sequence analysis of 27 alleles of each of the three Ras-related genes in Drosophila melanogaster indicates that they all have low levels of polymorphism but may experience slightly different evolutionary pressures. No amino acid replacement substitutions were indicated in any of the sequences, or in the sibling species D. simulans and D. mauritiana. The Dras1 gene, which is the major ras homologue in Drosophila, has less within-species variation in D. melanogaster relative to the amount of divergence from the sibling species than does Dras2, although the contrast was not significant by the HKA test. Dras2 appears to be maintaining two classes of haplotype in D. melanogaster, one of which is closer to the alleles observed in the sibling species, suggesting that this is not likely to be a pseudogene despite the absence of a mutant phenotype. Although differences in level of expression may affect the function of the genes, it is concluded that genetic variation in the Ras signal transduction pathways cannot be attributed to catalytic variation in the Ras proteins. Received: 5 November 1998 / Accepted: 26 March 1999     

Histone H1 Genes and Histone Gene Clusters in the Genus Drosophila

Whereas the genomes of many organisms contain several nonallelic types of linker histone genes, one single histone H1 type is known in Drosophila melanogaster that occurs in about 100 copies per genome. Amplification of H1 gene sequences from genomic DNA of wild type strains of D. melanogaster from Oregon, Australia, and central Africa yielded numerous clones that all exhibited restriction patterns identical to each other and to those of the known H1 gene sequence. Nucleotide sequences encoding the evolutionarily variable domains of H1 were determined in two gene copies of strain Niamey from central Africa and were found to be identical to the known H1 sequence. Most likely therefore, the translated sequences of D. melanogaster H1 genes do not exhibit intragenomic or intergenomic variations. In contrast, three different histone H1 genes were isolated from D. virilis and found to encode proteins that differ remarkably from each other and from the H1 of D. melanogaster and D. hydei. About 40 copies of H1 genes are organized in the D. virilis genome with copies of core histone genes in gene quintets that were found to be located in band 25F of chromosome 2. Another type of histone gene cluster is present in about 15 copies per genome and contains a variable intergenic sequence instead of an H1 gene. The H1 heterogeneity in D. virilis may have arisen from higher recombination rates than occur near the H1 locus in D. melanogaster and might provide a basis for formation of different chromatin subtypes. Received: 2 March 2000 / Accepted: 1 June 2000     

Rapid Concerted Evolution via Gene Conversion at the <Emphasis Type="BoldItalic">Drosophila hsp70</Emphasis> Genes

We analyzed nucleotide variation in the hsp70 genes of Drosophila melanogaster (five genes) and D. simulans (four genes) to characterize the homogenizing and diversifying roles of gene conversion in their evolution. Gene conversion within and between the 87A7 and 87C1 gene clusters homogenize the hsp70 coding regions; in both D. melanogaster and D. simulans, same-cluster paralogues are virtually identical, and large intercluster conversion tracts diminish 87A7/87C1 divergence. Same-cluster paralogues share many polymorphisms, consistent with frequent intracluster conversion. Shared polymorphism is highly biased toward silent variation; homogenizing conversion interacts with purifying selection. In contrast to the coding regions, some hsp70 flanking regions show conversion-mediated diversification. Strong reductions of nucleotide variability and linkage disequilibria among conversion-mediated sites in hsp70Ab and hsp70Bb alleles sampled from a single natural population are consistent with a selective sweep. Comparison of the D. melanogaster and D. simulans hsp70 genes reveals whole-family fixed differences, consistent with rapid propagation of novel mutations among duplicate genes. These results suggest that the homogenizing and diversifying roles of conversion interact to drive dynamic concerted evolution of the hsp70 genes. Received: 25 June 2001 / Accepted: 10 October 2001     

Identification and Characterization of the Cecropin Antibacterial Protein Gene Locus in Drosophila virilis

Cecropin is a type of antibacterial peptide that is synthesized in response to infection and has been characterized in many insect species and one mammal. The Cecropin locus of Drosophila melanogaster also contains the gene Andropin, which has been identified only in this species and encodes a male-specific antibacterial peptide. As a first step in studying the molecular evolution of the cecropin and andropin genes among Drosophila species, we have isolated genomic clones that cover the Cecropin locus in Drosophila virilis. The cloned region totals approximately 25 kb, within which a 9-kb fragment contains four cecropin genes and one pseudogene. All four genes have a high level of sequence homology to D. melanogaster Cecropin, about 80% identity in the coding regions, and the intron positions are conserved. As in D. melanogaster and other insects, κB-related cis-regulatory elements are found upstream of these cecropin genes. An Andropin-related sequence was not identified in D. virilis; however, genome Southern hybridizations suggest that Andropin-related sequences are present in at least the melanogaster species subgroup. Analysis of 19 insect cecropin genes identifies a common ancestral Cecropin before the divergence of Diptera and Lepidoptera. In addition, D. melanogaster and D. virilis can be identified by monophyletic clades for Cecropin. In contrast, the Lepidopteran species show polyphyletic relationships for duplicated cecropin genes. Received: 12 August 1996 / Accepted: 18 October 1996     

Dosage compensation of a retina-specific gene in Drosophila miranda

The X1R chromosome of Drosophila miranda and the 3L autosome of Drosophila melanogaster are thought to have originated from the ancestral D chromosomal element and therefore may contain the same set of genes. It is expected that these genes will be dosage compensated in D. miranda because of their X linkage. To test these possibilities and to study evolution of the dosage compensation mechanism, we used the 3L-linked autosomal head-specific gene 507ml of D. melanogaster to isolate the homologous gene (507 mr) from a D. miranda genomic library. In situ hybridization showed that gene 507 is located at the 12A region of the X1R chromosome of D. miranda, indicating that the chromosomal homology deduced by cytogenetic means is correct. Restriction analysis and cross-specific DNA and RNA blot hybridization revealed the presence of extensive restriction pattern polymorphism and lack of sequence similarity in some areas of the 507 mr and 507 ml DNA, including the 3 portion of the transcribed region. However, the 5 portion of the transcribed region and the DNA sequences, located approximately 0.8 kb upstream and 3 kb downstream from the 507 ml gene showed a high degreee of similarity with the DNA sequences of comparable regions of the 507 mr gene. In both species gene 507 codes for a highly abundant 1.8 kb RNA which is expressed in the retina of the compound eye. Although in D. miranda the males have one and the females have two copies of the 507 gene, the steady-state levels of the 507 mRNA in both sexes were found to be similar, indicating that gene 507 is dosage compensated in D. miranda. Thus, along with the disparate rates of evolution in different areas of the DNA associated with gene 507, in D. miranda this gene has come under the regulation of the X chromosomal dosage compensation mechanism.by M.L. Pardue     

Reconstructing the Diversification of α-Esterases: Comparing the Gene Clusters of Drosophila buzzatii and D. melanogaster

A cluster composed of 10 active α-esterase genes and a pseudogene is distributed over 60 kb in the Drosophila melanogaster genome. This paper describes the corresponding cluster in Drosophila buzzatii, whose lineage diverged from that of D. melanogaster when the subgenera Drosophila and Sophophora diverged about 50 Mya. With three exceptions we find that the composition of the cluster is conserved in the two lineages. The location of αE1 in D. melanogaster differs from that of its nearest relative in D. buzzatii, and αE4 has duplicated independently in the two lineages. The nature of these differences indicates that a mechanism exists whereby copies of genes can be placed in opposite orientation and nonadjacent positions within a gene cluster, although this does not seem to be a feature of earlier events in the cluster's evolution. The rates of amino acid change are not significantly different between orthologs, but the rates differ sevenfold among paralogs, indicating that very different selective forces are acting on the genes of the cluster. Mapping of sequence differences onto a model of the tertiary structure of the enzymes indicates that motifs contributing to substrate binding and catalysis have changed radically in the αE4s and suggest that this subgroup of α-esterases may be evolving into a substantially different functional niche. Received: 4 January 2000 / Accepted: 18 April 2000     

Rapid Evolution of the Male-Specific Antibacterial Protein Andropin Gene in Drosophila

Andropin, which encodes an antibacterial protein, is closely linked to the Cecropin gene cluster of D. melanogaster. Andropin and Cecropins are considered to have originated from one common ancestor. However, the expression pattern of Andropin is distinct from that of Cecropins, being restricted to the adult male ejaculatory duct. To elucidate the evolutionary process of Andropin, we have sequenced Andropin genes from D. melanogaster and its closely related species. In D. melanogaster, the nucleotide diversity of Andropin is remarkably low compared to that of Cecropin. In contrast, nonsynonymous substitutions of Andropin are conspicuously frequent between species. From genomic Southern analysis, Andropin-like genes are present in at least the melanogaster species subgroup. The series of present results suggests that Andropin was born in the course of constructing the Drosophila Cecropin gene family and then started to evolve rapidly, in contrast to Cecropins. Received: 10 August 2001 / Accepted: 29 October 2001     

The Amylase gene cluster on the evolving sex chromosomes of Drosophila miranda.

On the basis of chromosomal homology, the Amylase gene cluster in Drosophila miranda must be located on the secondary sex chromosome pair, neo-X (X2) and neo-Y, but is autosomally inherited in all other Drosophila species. Genetic evidence indicates no active amylase on the neo-Y chromosome and the X2-chromosomal locus already shows dosage compensation. Several lines of evidence strongly suggest that the Amy gene cluster has been lost already from the evolving neo-Y chromosome. This finding shows that a relatively new neo-Y chromosome can start to lose genes and hence gradually lose homology with the neo-X. The X2-chromosomal Amy1 is intact and Amy2 contains a complete coding sequence, but has a deletion in the 3''-flanking region. Amy3 is structurally eroded and hampered by missing regulatory motifs. Functional analysis of the X2-chromosomal Amy1 and Amy2 regions from D. miranda in transgenic D. melanogaster flies reveals ectopic AMY1 expression. AMY1 shows the same electrophoretic mobility as the single amylase band in D. miranda, while ectopic AMY2 expression is characterized by a different mobility. Therefore, only the Amy1 gene of the resident Amy cluster remains functional and hence Amy1 is the dosage compensated gene.     

Anopheles stephensi Dox-A2 Shares Common Ancestry with Genes from Distant Groups of Eukaryotes Encoding a 26S Proteasome Subunit and Is in a Conserved Gene Cluster

The sequence of a cloned Anopheles stephensi gene showed 72% inferred amino acid identity with Drosophila melanogaster Dox-A2 and 93% with its putative ortholog in Anopheles gambiae. Dox-A2 is the reported but herein disputed structural locus for diphenol oxidase A2. Database searches identified Dox-A2 related gene sequences from 15 non-insect species from diverse groups. Phylogenetic trees based on alignments of inferred protein sequences, DNA, and protein motif searches and protein secondary structure predictions produced results consistent with expectations for genes that are orthologous. The only inconsistency was that the C-terminus appears to be more primitive in the yeasts than in plants. In mammals, plants, and yeast these genes have been shown to code for a non-ATPase subunit of the PA700 (19S) regulatory complex of 26S proteasome. The analyses indicated that the insect genes contain no divergent structural features, which taken within an appraisal of all available data, makes the reported alternative function highly improbable. A plausible additional role, in which the 26S proteasome is implicated in regulation of phenol oxidase, would also apply to at least the mammalian genes. No function has yet been reported for the other included sequences. These were from genome projects and included Caenorhabiditus elegans, Arabidopsis thaliana, Fugu rubripes, and Toxoplasma gondii. A consensus of the results predicts a protein containing exceptionally long stretches of helix with a hydrophilic C-terminus. Phosphorylation site motifs were identified at two conserved positions. Possible SRY and GATA-1 binding motifs were found at conserved positions upstream of the mosquito genes. The location of A. stephensi Dox-A2 was determined by in situ hybridization at 34D on chromosome arm 3R. It is in a conserved gene cluster with respect to the other insects. However, the A. stephensi cluster contains a gene showing significant sequence identity to human and pigeon carnitine acetyltransferase genes, therefore showing divergence with the distal end of the D. melanogaster cluster. Received: 3 July 1998 / Accepted: 22 December 1999     

The Heat-Shock Gene hsp83 of Drosophila auraria: Genomic Organization, Nucleotide Sequence, and Long Antiparallel Coupled ORFs (LAC ORFs)

The genomic organization of the hsp83 gene of Drosophila auraria, a far-eastern endemic species belonging to the montium subgroup of the melanogaster species group, is presented here. Based on in situ hybridization on polytene chromosomes, cDNA and genomic clone mapping, nucleotide sequencing, and genomic Southern analysis, hsp83 is shown to be present as a single-copy gene at locus 64B on the 3L chromosome arm in D. auraria. This gene is organized into two exons separated by a 929-bp intron. The first exon represents the mRNA leader sequence and is not translated, while the coding region, having a length of 2,151 bp, is solely included in the second exon. Nucleotide sequence comparisons of D. auraria hsp83 with homologous sequences from other organisms show high conservation of the coding region (88–92% identity) in the genus Drosophila, in addition to the conserved genomic organization of two-exons–one-intron, of comparable size and arrangement. A phylogenetic tree based on the protein sequences of homologous genes from representative organisms is in accord with the accredited phylogenetic position of D. auraria. In the hsp83 gene region, a second case of long antiparallel coupled open reading frames (LAC ORFs) for this species was found. The antiparallel to the hsp83 gene ORF is 1,554 bases long, while the two ORFs overlap has a size of 1,548 bp. The anti-hsp83 ORF does not show significant homology to any known gene sequences. In addition, no similar LAC ORF structures were found in homologous gene regions of other organisms. Received: 18 April 1997 / Accepted: 1 August 1997     

Molecular Characterization and Evolution of the Amylase Multigene Family of Drosophila ananassae

Drosophila ananassae is known to produce numerous alpha-amylase variants. We have cloned seven different Amy genes in an African strain homozygous for the AMY1,2,3,4 electrophoretic pattern. These genes are organized as two main clusters: the first one contains three intronless copies on the 2L chromosome arm, two of which are tandemly arranged. The other cluster, on the 3L arm, contains two intron-bearing copies. The amylase variants AMY1 and AMY2 have been assigned to the intronless cluster, and AMY3 and AMY4 to the second one. The divergence of coding sequences between clusters is moderate (6.1% in amino acids), but the flanking regions are very different, which could explain their differential regulation. Within each cluster, coding and noncoding regions are conserved. Two very divergent genes were also cloned, both on chromosome 3L, but very distant from each other and from the other genes. One is the Amyrel homologous (41% divergent), the second one, Amyc1 (21.6% divergent) is unknown outside the D. ananassae subgroup. These two genes have unknown functions. Received: 30 May 2000 / Accepted: 17 July 2000     

Selection and Methionine Accumulation in the Fat Body Protein 2 Gene (FBP2), a Duplicate of the Drosophila Alcohol Dehydrogenase (ADH) Gene

The Drosophila fat body protein 2 gene (Fbp2) is an ancient duplication of the alcohol dehydrogenase gene (Adh) which encodes a protein that differs substantially from ADH in its methionine content. In D. melanogaster, there is one methionine in ADH, while there are 51 (20% of all amino acids) in FBP2. Methionine is involved in 46% of amino acid replacements when Fbp2 DNA sequences are compared between D. melanogaster and D. pseudoobscura. Methionine accumulation does not affect conserved residues of the ADH-ADH^r-FBP2 multigene family. The multigene family has evolved by replacement of mildly hydrophobic amino acids by methionine with no apparent reversion. Its short-term evolution was compared between two Drosophila species, while its long-term evolution was compared between two genera belonging respectively to acalyptrate and calyptrate Diptera, Drosophila and Sarcophaga. The pattern of nucleotide substitution was consistent with an independent accumulation of methionines at the Fbp2 locus in each lineage. Under a steady-state model, the rate of methionine accumulation was constant in the lineage leading to Drosophila, and was twice as fast as that in the calyptrate lineage. Substitution rates were consistent with a slight positive selective advantage for each methionine change in about one-half of amino acid sites in Drosophila. This shows that selection can potentially account for a large proportion of amino acid replacements in the molecular evolution of proteins. Received: 12 December 1994 / Accepted: 15 April 1996     

Frequency of telomere repeat units in the Drosophila miranda genome

Cloned DNA fragments of Drosophila miranda which label all chromosome ends show a basic tandem repeat unit of 4.4 kb. The D. miranda telomere specific tandem repeats do not cross-hybridize with genomic D. melanogaster DNA which itself contains telomere repeat units of 3 kb. For a more detailed analysis of the functional criteria of telomere specific sequences we determined the repetition frequency of the tandem repeat units. As a low estimate we found a repetition frequency of 20 for female D. miranda DNA. This is on average equivalent to 2 telomere repeat units per chromosome end in the female D. miranda karyotype. However, a variable number of tandem repeat units per chromosome end would describe more closely the obtained differences in the labeling intensity between the individual chromosomes (X¹L-5). For the D. miranda male DNA we determined a repetition frequency of 90. The frequency difference of 70 copies between male and female DNA must be due to the Y-chromosome.     

Duplication and divergence of the genes of the α-esterase cluster ofDrosophila melanogaster

The α-esterase cluster of D. melanogaster contains 11 esterase genes dispersed over 60 kb. Embedded in the cluster are two unrelated open reading frames that have sequence similarity with genes encoding ubiquitin-conjugating enzyme and tropomyosin. The esterase amino acid sequences show 37–66% identity with one another and all but one have all the motifs characteristic of functional members of the carboxyl/cholinesterase multigene family. The exception has several frameshift mutations and appears to be a pseudogene. Patterns of amino acid differences among cluster members in relation to generic models of carboxyl/cholinesterase protein structure are broadly similar to those among other carboxyl/cholinesterases sequenced to date. However the α-esterases differ from most other members of the family in: their lack of a signal peptide; the lack of conservation in cysteines involved in disulfide bridges; and in four indels, two of which occur in or adjacent to regions that align with proposed substrate-binding sites of other carboxyl/cholinesterases. Phylogenetic analyses clearly identify three simple gene duplication events within the cluster. The most recent event involved the pseudogene which is located in an intron of another esterase gene. However, relative rate tests suggest that the pseudogene remained functional after the duplication event and has become inactive relatively recently. The distribution of indels also suggests a deeper node in the gene phylogeny that separates six genes at the two ends of the cluster from a block of five in the middle. Received: 18 January 1996 / Accepted: 12 March 1996     

Gene Duplication and Gene Conversion in the Caenorhabditis elegans Genome

A comprehensive analysis of duplication and gene conversion for 7394 Caenorhabditis elegans genes (about half the expected total for the genome) is presented. Of the genes examined, 40% are involved in duplicated gene pairs. Intrachromosomal or cis gene duplications occur approximately two times more often than expected. In general the closer the members of duplicated gene pairs are, the more likely it is that gene orientation is conserved. Gene conversion events are detectable between only 2% of the duplicated pairs. Even given the excesses of cis duplications, there is an excess of gene conversion events between cis duplicated pairs on every chromosome except the X chromosome. The relative rates of cis and trans gene conversion and the negative correlation between conversion frequency and DNA sequence divergence for unconverted regions of converted pairs are consistent with previous experimental studies in yeast. Three recent, regional duplications, each spanning three genes are described. All three have already undergone substantial deletions spanning hundreds of base pairs. The relative rates of duplication and deletion may contribute to the compactness of the C. elegans genome. Received: 30 July 1998 / Accepted: 12 October 1998     

Characterization of the Hydra Lamin and Its Gene: A Molecular Phylogeny of Metazoan Lamins

We report sequences for nuclear lamins from the teleost fish Danio and six invertebrates. These include two cnidarians (Hydra and Tealia), one priapulid, two echinoderms, and the cephalochordate Branchiostoma. Combining these results with earlier data on Drosophila, Caenorhabditis elegans, and various vertebrates, the following conclusions on lamin evolution can be drawn. First, all invertebrate lamins resemble in size the vertebrate B-type lamin. Second, all lamins described previously for amphibia, birds and mammals as well as the first lamin of a fish, characterized here, show a cluster of 7 to 12 acidic residues in the tail domain. Since this acidic cluster is absent from all invertebrate lamins including that of the cephalochordate Branchiostoma, it was acquired with the vertebrate lineage. The larger A-type lamin of differentiated cells must have arisen subsequently by gene duplication and insertion of an extra exon. This extra exon of the vertebrate A-lamins is the only major change in domain organization in metazoan lamin evolution. Third, the three introns of the Hydra and Priapulus genes correspond in position to the last three introns of vertebrate B-type lamin genes. Thus the entirely different gene organization of the C. elegans and Drosophila Dmo genes seems to reflect evolutionary drift, which probably also accounts for the fact that C. elegans has the most diverse lamin sequence. Finally we discuss the possibility that two lamin types, a constitutively expressed one and a developmentally regulated one, arose independently on the arthropod and vertebrate lineages. Received: 4 February 1999 / Accepted: 1 April 1999