Multiple calmodulin mRNA species are derived from two distinct genes.

We have observed three calmodulin mRNA species in rat tissues. In order to know from how many expressed genes they are derived, we have investigated the genomic organization of calmodulin genes in the rat genome. From a rat brain cDNA library, we obtained two kinds of cDNAs (pRCM1 and pRCM3) encoding authentic calmodulin. DNA sequence analysis of these cDNA clones revealed substitutions of nucleotides at 73 positions of 450 nucleotides in the coding region, although the amino acid sequences of these calmodulins are exactly the same. DNA sequences in the 5' and 3' noncoding regions are quite different between these two cDNAs. From these results, we conclude that they are derived from two distinct bona fide calmodulin genes, CaMI (pRCM1) and CaMII (pRCM3). Total genomic Southern hybridization suggested four distinct calmodulin-related genes in the rat genome. By cloning and sequencing the calmodulin-related genes from rat genomic libraries, we demonstrated that the other two genes are processed pseudogenes generated from the CaMI (lambda SC9) and CaMII (lambda SC8) genes, respectively, through an mRNA-mediated process of insertions. Northern blotting showed that the CaMI gene is transcribed in liver, muscle, and brain in similar amounts, whereas the CaMII gene is transcribed mainly in brain. S1 nuclease mapping indicated that the CaMI gene produced two mRNA species (1.7 and 4 kilobases), whereas the CaMII gene expressed a single mRNA species (1.4 kilobases).     

Concerted evolution within a trypsin gene cluster in Drosophila.

A cluster of four trypsin genes has previously been localized to cytological position 47D-F of the Drosophila melanogaster genome. One of these genes had been sequenced, and the presence of the other three genes was identified by cross-hybridization. Here, we present the DNA sequence of the entire genomic region encoding these four trypsin genes. In addition to the four previously inferred genes, we have identified a fifth trypsin-coding sequence located within this gene cluster. This new gene shows a high degree of sequence divergence (more than 30%) from the other four genes, although it retains all of the functional motifs that are characteristic of trypsin-coding sequences. In order to trace the molecular evolution of this gene cluster, we isolated and sequenced the homologous 7-kb region from the closely related species Drosophila erecta. A comparison of the DNA sequences between the two species provides strong evidence for the concerted evolution of some members of this gene family. Two genes within the cluster are evolving in concert, while a third gene appears to be evolving independently. The remaining two genes show an intermediate pattern of evolution. We propose a simple model, involving chromosome looping and gene conversion, to explain the relatively complex patterns of molecular evolution within this gene cluster.     

Patterns of nucleotide polymorphism and divergence in the odorant-binding protein genes OS-E and OS-F: analysis in the melanogaster species subgroup of Drosophila

The Olfactory Specific-E and -F genes (OS-E and OS-F) belong to the odorant-binding protein gene family, which includes the general odorant-binding proteins and the pheromone-binding proteins. In Drosophila melanogaster, these genes are arranged in tandem in a genomic region near the centromere of chromosome arm 3R. We examined the pattern of DNA sequence variation in an approximately 7-kb genomic region encompassing the two OS genes in four species of the melanogaster subgroup of Drosophila and in a population sample of D. melanogaster. We found that both the OS-E and the OS-F gene are present in all surveyed species. Nucleotide divergence estimates would support that the two genes are functional, although they diverge in their functional constraint. The pattern of nucleotide variation in D. melanogaster also differed between genes. Variation in the OS-E gene region exhibited an unusual and distinctive pattern: (i) a relatively high number of fixed amino acid replacements in the encoded protein and (ii) a peak of nucleotide polymorphism around the OS-E gene. These results are unlikely under the neutral model and suggest the action of natural selection in the evolution of the two odorant-binding protein genes.     

Cloning and characterization of avfA and omtB genes involved in aflatoxin biosynthesis in three Aspergillus species

A collection of lethal and semi-lethal P-element insertions in the 70CD region of chromosome 3 of Drosophila melanogaster was used to investigate genes and gene arrangements by a combination of genetic, cytological, functional and molecular methods. The 12 lethal insertions studied fall into seven complementation groups of six genes. Lethal phases, expression patterns and other phenotypic aspects of these genes were determined. The genes and additional available sequences were placed on cloned genomic DNA fragments and arranged in an EcoRI map of 150kb that covers approximately the bands 70C7-8 to 70D1. Determination of deficiency breakpoints links the genetic, physical and molecular data. The sequences adjacent to seven independent P-element insertions were established after plasmid rescue or polymerase chain reaction. Similarity searches allowed the assignment of the P-element insertions to known mutations, expressed sequence tags, sequence tagged sites, or homologous genes of other species. Among these were identified a putative transacylase, a putative cell cycle gene, and the gene responsible for the dominant Polycomb-suppressor phenotype of devenir. The genomic sequence of the l(3)70Ca/b gene reveals a novel heat shock protein (hsc70Cb). l(3)70Da was identified as a member of the CDC48/PEX1 ATPase family and its coding sequence was determined.     

Rapid evolution of two odorant-binding protein genes, Obp57d and Obp57e, in the Drosophila melanogaster species group

Genes encoding odorant-binding protein (OBP) form a large family in an insect genome. Two OBP genes, Obp57d and Obp57e, were previously identified to be involved in host-plant recognition in Drosophila sechellia. Here, by comparing the genomic sequences at the Obp57d/e locus from 27 Drosophila species, we found large differences in gene number between species. Phylogenetic analysis revealed that Obp57d and Obp57e in the D. melanogaster species group arose by gene duplication of an ancestral OBP gene that remains single in the obscura species group. Further gain and loss of OBP genes were observed in several lineages in the melanogaster group. Site-specific analysis of evolutionary rate suggests that Obp57d and Obp57e have functionally diverged from each other. Thus, there are two classes of gene number differences in the Obp57d/e region: the difference of the genes that have functionally diverged from each other and the difference of the genes that appear to be functionally identical. Our analyses demonstrate that these two classes of differences can be distinguished by comparisons of many genomic sequences from closely related species.     

Isolation of three novel Drosophila melanogaster genes containing repetitive sequences rich in GCN triplets.

Several cDNA and genomic clones were isolated from Drosophila melanogaster gene libraries by hybridization with a region of a mammalian gene that contains a simple repetitive sequence of six GCN repeats. One of the cDNA clones, E6, was completely sequenced and it was shown that it contains a region of 16 GCN repeats; these repeats encode a polyalanine stretch within a long open reading frame. The sequencing of three different genomic clones (A, B, and D) revealed that all the isolated Drosophila clones are similar to one another in a short region containing variable numbers of the GCN repeat. The genomic clone B was found to be the genomic counterpart of the cDNA clone E6. The other genomic clones, A and D, also hybridize with Drosophila cDNA clones at high stringency. These results indicate that the short GCN repetitive sequences, which we have named ala, are found within transcribed regions of the Drosophila genome. These Drosophila genes containing the ala repeat do not show significant sequence similarity to any presently known gene; we have named these novel genes ala-A, ala-B, and ala-D. The cDNA clone from gene ala-B was named ala-E6.     

Structure of a gene for rat calmodulin

The structural organization of the entire rat calmodulin gene was determined by cloning and sequencing overlapping genomic and cDNA clones from rat genomic and brain cDNA libraries. The intron/exon organization was determined by direct comparison of these sequences. Rat calmodulin gene is 9000 bases long and consisted of six exons interrupted by introns of variable sizes. The first intron separates the initiation codon (ATG) from the coding region of the protein. Three out of four intron/exon junctions in the coding region reside in the middle of calcium binding subdomains and do not correlate with the quarterly divided intramolecular homology of the protein. Their positions exactly coincide with those of the corrected version of chicken calmodulin gene. The rat calmodulin gene harbors a stretch of sequences homologous to a rat middle repetitive "identifier sequence" in the middle of the third intron. Analysis of the immediate 5' upstream region detected a TATA box (TATATATAT) and three C-G boxes (CCGCCC) but not a CAT box (CCAAT). A conserved sequence (GCGCCGCGYCYYGGGGGC) was found at -125 for rat and at -204 for chicken calmodulin genes.     

Origin of new genes and source for N-terminal domain of the chimerical gene, jingwei, in Drosophila

This paper deals with a general question posed by the origin of new processed chimerical genes: when a new retrosequence inserts into a new genome position, how does it become activated and acquire novel protein function by recruiting new functional domains and regulatory elements? Jingwei (jgw), a newly evolved functional gene with a chimerical structure in Drosophila, provides an opportunity to examine such questions. The source of its exon encoding C-terminal peptide has been identified as an Adh retrosequence, which extends the concept of exon shuffling from recombination to retroposition as a general molecular mechanism for the origin of a new gene. However, the origin of 5' exons remains unclear. We examined two hypotheses concerning the origin of these non-Adh-derived jgw exons: (i) these exons might originate from a unique genomic sequence that fortuitously evolved a standard intron-exon structure and regulatory sequence for jgw; (ii) these exons might be a duplicate of an unrelated previously existing gene. Genomic Southern analysis, in conjunction with construction and screening of a genomic bookshelf (sub-library), was conducted in a group of Drosophila species. The results demonstrated that there are duplicate genes containing the same structure as the recruited portion of jgw. We name this duplicate gene in Drosophila teissieri and Drosophila yakuba and its orthologous gene in Drosophila melanogaster as yellow-emperor (ymp). Thus, the 5' exons/introns originated from a previously existing gene that provided new modules with specific sub-function to create jgw.     

Low codon bias and high rates of synonymous substitution in Drosophila hydei and D. melanogaster histone genes

We have evaluated codon usage bias in Drosophila histone genes and have obtained the nucleotide sequence of a 5,161-bp D. hydei histone gene repeat unit. This repeat contains genes for all five histone proteins (H1, H2a, H2b, H3, and H4) and differs from the previously reported one by a second EcoRI site. These D. hydei repeats have been aligned to each other and to the 5.0-kb (i.e., long) and 4.8-kb (i.e., short) histone repeat types from D. melanogaster. In each species, base composition at synonymous sites is similar to the average genomic composition and approaches that in the small intergenic spacers of the histone gene repeats. Accumulation of synonymous changes at synonymous sites after the species diverged is quite high. Both of these features are consistent with the relatively low codon usage bias observed in these genes when compared with other Drosophila genes. Thus, the generalization that abundantly expressed genes in Drosophila have high codon bias and low rates of silent substitution does not hold for the histone genes.      

The Evolution of Duplicate Glyceraldehyde-3-Phosphate Dehydrogenase Genes in Drosophila

In Drosophila melanogaster there are two genes which encode the enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Gapdh-43E and Gapdh-13F. We have shown that Gapdh-43E codes for the GAPDH subunit with an apparently larger molecular weight while Gapdh-13F encodes the GAPDH subunit having an apparently smaller molecular weight. Immunoblots of sodium dodecyl sulfate gels were used to survey species from throughout the genus and results indicated that two classes of GAPDH subunits are present only in Drosophila species of the melanogaster and takahashi subgroups of the melanogaster group. Only the smaller subunit is found in species of the obscura group while all other species have only a large subunit. Drosophila hydei was analyzed at the DNA level as a representative species of the subgenus Drosophila. The genome of this species has a single Gapdh gene which is localized at a cytogenetic position likely to be homologous to Gapdh-43 E of D. melanogaster. Comparison of its sequence with the sequence of the D. melanogaster Gapdh genes indicates that the two genes of D. melanogaster are more similar to one another than either is to the gene from D. hydei. The Gapdh gene from D. hydei contains an intron following codon 29. Neither Gapdh gene of D. melanogaster has an intron within the coding region. Southern blots of genomic DNA were used to determine which species have duplicate Gapdh genomic sequences. Gene amplification was used to determine which species have a Gapdh gene that is interrupted by an intron. Species of the subgenus Drosophila have a single Gapdh gene with an intron. Species of the willistoni and saltans groups have a single Gapdh gene that does not contain an intron.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure of rat calmodulin processed genes with implications for a mRNA-mediated process of insertion

Two distinct processed calmodulin genes of rat (lambda SC8 and lambda SC9) were identified, cloned and their DNA sequences determined. The existence of direct repeats of 19 base-pairs for lambda SC8 or 9 base-pairs for lambda SC9 at both ends of the coding plus non-coding regions suggested a possible involvement of a mRNA-mediated process of insertion. Total genomic Southern hybridization suggested the existence of at least three different calmodulin-related genes in the rat genome. The other gene was the bona fide calmodulin gene (lambda SC4) which was split into at least five exons. lambda SC9 contained insertions of one nucleotide and two 17 base-pair direct repeats in the coding region. These insertions cause frameshift mutations probably preventing it from encoding a functional calmodulin. It also carried an insertion of a rat middle repetitive sequence, identifier sequence (IDS: Sutcliffe et al., 1982) in the 3'-non-coding region. Otherwise, it consisted of an almost identical DNA sequence to that of the bona fide calmodulin gene (lambda SC4), including the 3'-non-coding region down to the poly(A) recognition signal, A-A-T-A-A-A. On the other hand, lambda SC8 did not possess frameshift mutations in the coding region, and hence was capable of encoding a functional protein. In fact, a probe specific to the lambda SC8 sequence identified a band in Northern blotting whose size was 300 nucleotides smaller than that of authentic calmodulin mRNA. Comparison of the nucleotide sequences showed that only the coding regions of these two processed genes were homologous, indicating that the divergence of these two processed genes from the common ancestor calmodulin was an ancient event.     

Novel eye-specific calmodulin methylation characterized by protein mapping in Drosophila melanogaster

Post-translational methylation of the epsilon-amino group of lysine residues regulates a number of protein functions. Calmodulin, a key modulator of intracellular calcium signaling, is methylated on lysine 115 in many species. Although the amino acid sequence of calmodulin is highly conserved in eukaryotes, it has been shown that lysine 115 is not methylated in Drosophila calmodulin and no other methylation site has been reported. In this study, we characterized in vivo modification states of Drosophila calmodulin using proteomic methodology involving the protein mapping of microdissected Drosophila tissues on 2-D gels. We found that Drosophila calmodulin was highly expressed in methylated forms in the compound eye, whereas its methylation was hardly detected in other tissues. We identified that lysine 94 located in an EF-hand III is the methylation site in Drosophila calmodulin. The predominance of methylated calmodulin in the compound eye may imply the involvement of calmodulin in photoreceptor-specific functions through methylation.     

Hsp70 duplication in the Drosophila melanogaster species group: how and when did two become five?

To determine how the modern copy number (5) of hsp70 genes in Drosophila melanogaster evolved, we localized the duplication events that created the genes in the phylogeny of the melanogaster group, examined D. melanogaster genomic sequence to investigate the mechanisms of duplication, and analyzed the hsp70 gene sequences of Drosophila orena and Drosophila mauritiana. The initial two-to-four hsp70 duplication occurred 10--15 MYA, according to fixed in situ hybridization to polytene chromosomes, before the origin and divergence of the melanogaster and five other species subgroups of the melanogaster group. Analysis of more than 30 kb of flanking sequence surrounding the hsp70 gene clusters suggested that this duplication was likely a retrotransposition. For the melanogaster subgroup, Southern hybridization and an hsp70 restriction map confirmed the conserved number (4) and arrangement of hsp70 genes in the seven species other than D. melanogaster. Drosophila melanogaster is unique; tandem duplication and gene conversion at the derived cluster yielded a fifth hsp70 gene. The four D. orena hsp70 genes are highly similar and concertedly evolving. In contrast, the D. mauritiana hsp70 genes are divergent, and many alleles are nonfunctional. The proliferation, concerted evolution, and maintenance of functionality in the D. melanogaster hsp70 genes is consistent with the action of natural selection in this species.     

Schistosoma haematobium: analysis of eggshell protein genes and their expression

Eggshell protein genes of Schistosoma mansoni that encode a 14 kDa protein have been shown to be highly conserved and expressed in a sex-, tissue-, and temporal-specific manner. To initiate studies on the eggshell protein genes of S. haematobium, a cDNA probe, pSMf 61-46, representing a S. mansoni eggshell protein mRNA was used to screen a S. haematobium genomic library. Of the seven independent recombinant clones isolated, two (lambda SH 2-1 and lambda SH 6-1) were analyzed and compared to those of S. mansoni. lambda SH 2-1 and lambda SH 6-1 each contain a different genomic copy of the gene encoding a 19.8 and 17.6 kDa protein, respectively. This is due to an additional 78 bp present in the coding region of lambda SH 2-1 relative to lambda SH 6-1. The rest of the coding sequences are identical, and the 5' and 3' untranslated regions are nearly identical. The deduced amino acid sequences of S. haematobium eggshell proteins are very rich in glycine (47 and 50%) when compared to 43.5% glycine in the protein encoded by S. mansoni. Long stretches of glycines, as many as 15 in a row, occur in the S. haematobium sequence. DNA comparison of the eggshell protein genes of the two schistosome species yielded an overall homology of 83.1%. The homology is much higher in the 5' and 3' untranslated regions than in the protein-coding regions. Genomic clones of both species contained second open reading frames, which appeared to be kept open as a consequence of the amino acid composition of the other. There are no introns in S. haematobium or S. mansoni eggshell protein genes, and the genomic Southern data indicated a similar arrangement of these genes in the genome of both species. Primer extension experiments and dideoxynucleotide sequencing of the RNA determined the mRNA cap site sequence as ATCAT and ATCAC in lambda SH 2-1 and lambda SH 6-1, respectively. Northern blot analysis determined the size of the mRNA to be about 1.0 kp. Expression of the RNA from these genes appears to be regulated in a manner similar to the corresponding genes in S. mansoni. mRNA is found only in mature females and first appears at 70 days after infection of hamsters. DNA sequence comparisons of the 5' flanking regions of S. haematobium and S. mansoni eggshell protein genes to each other and to those of silkmoth and Drosophila revealed several short sequence elements that are shared.(ABSTRACT TRUNCATED AT 400 WORDS)