Evolution of osmosensing signal transduction in Metazoa: stress-activated protein kinases p38 and JNK

Sponges (Porifera) represent the most basal branch of the Metazoa alive today. We show that two central stress-activated protein kinases involved in the osmosensing pathway, p38 mitogen-activated protein kinase (MAPK) and JNK, can complement for the ancestral MAPK Hog1 in the yeast Saccharomyces cerevisiae. S. cerevisiae mutants lacking Hog1 (hog1-Delta 1) have been complemented with the sponge SDJNK and SDp38 genes. Western blotting has revealed that, after transformation, the hog1-Delta 1+ SDJNK(sense) and hog1-Delta 1+ SDp38(sense) clones express the sponge proteins. Functional studies have demonstrated that the complemented clones grow under hyperosmotic conditions (0.6 M NaCl). Furthermore, the expressed sponge kinases undergo phosphorylation in S. cerevisiae at 0.6 M NaCl. This report documents that the metazoan signal transduction kinases, p38 and JNK, which were originally derived from an common ancestor with yeast HOG1, have retained their function after their specification.     

Nucleotide sequence of the triosephosphate isomerase gene from Aspergillus nidulans: implications for a differential loss of introns

A functional cDNA from Aspergillus nidulans encoding triosephosphate isomerase (TPI) was isolated by its ability to complement a tpi1 mutation in Saccharomyces cerevisiae. This cDNA was used to obtain the corresponding gene, tpiA. Alignment of the cDNA and genomic DNA nucleotide sequences indicated that tpiA contains five introns. The intron positions in the tpiA gene were compared with those in the TPI genes of human, chicken, and maize. One intron is present at an identical position in all four organisms, two other introns are located in similar positions in A. nidulans and maize, and the remaining two introns are unique to A. nidulans. These Aspergillus-specific introns are located in regions of the protein that were predicted to be interrupted by introns based on analysis of a Go plot of chicken TPI. These comparisons are discussed in relation to the evolution of introns within TPI genes.     

Conserved intron positions in FGFR genes reflect the modular structure of FGFR and reveal stepwise addition of domains to an already complex ancestral FGFR

We have analyzed the evolution of fibroblast growth factor receptor (FGFR) tyrosine kinase genes throughout a wide range of animal phyla. No evidence for an FGFR gene was found in Porifera, but we tentatively identified an FGFR gene in the placozoan Trichoplax adhaerens. The gene encodes a protein with three immunoglobulin-like domains, a single-pass transmembrane, and a split tyrosine kinase domain. By superimposing intron positions of 20 FGFR genes from Placozoa, Cnidaria, Protostomia, and Deuterostomia over the respective protein domain structure, we identified ten ancestral introns and three conserved intron groups. Our analysis shows (1) that the position of ancestral introns correlates to the modular structure of FGFRs, (2) that the acidic domain very likely evolved in the last common ancestor of triploblasts, (3) that splicing of IgIII was enabled by a triploblast-specific insertion, and (4) that IgI is subject to substantial loss or duplication particularly in quickly evolving genomes. Moreover, intron positions in the catalytic domain of FGFRs map to the borders of protein subdomains highly conserved in other serine/threonine kinases. Nevertheless, these introns were introduced in metazoan receptor tyrosine kinases exclusively. Our data support the view that protein evolution dating back to the Cambrian explosion took place in such a short time window that only subtle changes in the domain structure are detectable in extant representatives of animal phyla. We propose that the first multidomain FGFR originated in the last common ancestor of Placozoa, Cnidaria, and Bilateria. Additional domains were introduced mainly in the ancestor of triploblasts and in the Ecdysozoa.     

Bruton tyrosine kinase-like protein, BtkSD, is present in the marine sponge Suberites domuncula

Sponges, the simplest and most ancient phylum of Metazoa, encode in their genome complex and highly sophisticated proteins that evolved together with multicellularity and are found only in metazoan animals. We report here the finding of a Bruton tyrosine kinase (BTK)-like protein in the marine sponge Suberites domuncula (Demospongiae). The nucleotide sequence of one sponge cDNA predicts a 700-aa-long protein, which contains all of the characteristic domains for the Tec family of protein tyrosine kinases (PTKs). The highest homology (38% identity, 55% overall similarity) was found with human BTK and TEC PTKs. Sponge PTK was therefore named BtkSD. Human BTK is involved in the maturation of B cells and mutations in the BTK gene cause X-linked agammaglobulinemia. Kinases from the Tec family are not present in Caenorhabditis elegans and, until now, they were found only in insects and higher animal taxa. Our finding implies that the BTK/TEC genes are of a very ancient origin.     

Molecular evolution of calmodulin-like domain protein kinases (CDPKs) in plants and protists

Many genes for calmodulin-like domain protein kinases (CDPKs) have been identified in plants and Alveolate protists. To study the molecular evolution of the CDPK gene family, we performed a phylogenetic analysis of CDPK genomic sequences. Analysis of introns supports the phylogenetic analysis; CDPK genes with similar intron/exon structure are grouped together on the phylogenetic tree. Conserved introns support a monophyletic origin for plant CDPKs, CDPK-related kinases, and phosphoenolpyruvate carboxylase kinases. Plant CDPKs divide into two major branches. Plant CDPK genes on one branch share common intron positions with protist CDPK genes. The introns shared between protist and plant CDPKs presumably originated before the divergence of plants from Alveolates. Additionally, the calmodulin-like domains of protist CDPKs have intron positions in common with animal and fungal calmodulin genes. These results, together with the presence of a highly conserved phase zero intron located precisely at the beginning of the calmodulin-like domain, suggest that the ancestral CDPK gene could have originated from the fusion of protein kinase and calmodulin genes facilitated by recombination of ancient introns. Received: 11 July 2000 / Accepted: 18 April 2001     

Gene structure of two-domain arginine kinases from Anthopleura japonicus and Pseudocardium sachalinensis

Unusual two-domain arginine kinases (AKs) arose independently at least two times during molecular evolution of phosphagen kinases: AKs from the primitive sea anemone Anthopleura japonicus and from the clam Pseudocardium sachalinensis. To elucidate its unusual evolution, the structures of Anthopleura and Pseudocardium AK genes have been determined. The Anthopleura gene consisted of 4 exons and 3 introns: two domains are linked by a bridge intron, and each domain contains one intron in different positions. On the other hand, the Pseudocardium gene consisted of 10 exons and 9 introns: two domains are also linked by a bridge intron, and domains 1 and 2 contains 3 and 5 introns, respectively, of which 3 introns are located in exactly same positions. Since the two domains of Pseudocardium AK are estimated to have diverged about 290 million years ago, the 3 introns have been conserved at least for this long. Comparison of intron positions in Anthopleura, Pseudocardium and C. elegans AK genes indicates that there is no intron conserved through the three AK lineages, in sharp contrast to relatively conservative intron positions in creatine kinase (CK) gene family.     

Intron loss and gain during evolution of the catalase gene family in angiosperms.

Angiosperms (flowering plants), including both monocots and dicots, contain small catalase gene families. In the dicot, Arabidopsis thaliana, two catalase (CAT) genes, CAT1 and CAT3, are tightly linked on chromosome 1 and a third, CAT2, which is more similar to CAT1 than to CAT3, is unlinked on chromosome 4. Comparison of positions and numbers of introns among 13 angiosperm catalase genomic sequences indicates that intron positions are conserved, and suggests that an ancestral catalase gene common to monocots and dicots contained seven introns. Arabidopsis CAT2 has seven introns; both CAT1 and CAT3 have six introns in positions conserved with CAT2, but each has lost a different intron. We suggest the following sequence of events during the evolution of the Arabidopsis catalase gene family. An initial duplication of an ancestral catalase gene gave rise to CAT3 and CAT1. CAT1 then served as the template for a second duplication, yielding CAT2. Intron losses from CAT1 and CAT3 followed these duplications. One subclade of monocot catalases has lost all but the 5''-most and 3''-most introns, which is consistent with a mechanism of intron loss by replacement of an ancestral intron-containing gene with a reverse-transcribed DNA copy of a fully spliced mRNA. Following this event of concerted intron loss, the Oryza sativa (rice, a monocot) CAT1 lineage acquired an intron in a novel position, consistent with a mechanism of intron gain at proto-splice sites.     

The mitogen-activated protein kinase p38 pathway is conserved in metazoans: cloning and activation of p38 of the SAPK2 subfamily from the sponge Suberites domuncula

Our recent data suggest that during auto- and allograft recognition in sponges (Porifera), cytokines are differentially expressed. Since the mitogen-activated protein kinase (MAPK) signal transduction modulates the synthesis and release of cytokines, we intended to identify one key molecule of this pathway. Therefore, a cDNA from the marine sponge Suberites domuncula encoding the MAPK was isolated and analyzed. Its encoded protein is 366 amino acids long (calculated Mr 42 209), has a TGY dual phosphorylation motif in protein kinase subdomain VIII and displays highest overall similarity to the mammalian p38 stress activated protein kinase (SAPK2), one subfamily of MAPKs. The sponge protein was therefore termed p38_SD. The overall homology (identity and similarity) between p38_SD and human p38alpha (CSBP2) kinase is 82%. One feature of the sponge kinase is the absence of threonine at position 106. In human p38alpha MAPK this residue is involved in the interaction with the specific pyridinyl-imidazole inhibitor; T106 is replaced in p38_SD by methionine. Inhibition studies with the respective inhibitor SB 203580 showed that it had no effect on the phosphorylation of the p38 substrate myelin basic protein. A stress responsive kinase Krs_SD similar to mammalian Ste20 kinases, upstream regulators of p38, had already previously been found in S. domuncula. The S. domuncula p38 MAPK is phosphorylated after treatment of the animal in hypertonic medium. In contrast, exposure of cells to hydrogen peroxide, heat shock and ultraviolet light does not cause any phosphorylation of p38. It is concluded that sponges, the oldest and most simple multicellular animals, utilize the conserved p38 MAPK signaling pathway, known to be involved in stress and immune (inflammatory) responses in higher animals.     

Nucleotide sequence and intron structure of the apocytochrome b gene of Neurospora crassa mitochondria

The sequence of the apocytochrome b (cob) gene of Neurospora crassa has been determined. The structural gene is interrupted by two intervening sequences of approximately 1260 bp each. The polypeptide encoded by the exons shows extensive homology with the cob proteins of Aspergillus nidulans and Saccharomyces cerevisiae (79% and 60%, respectively). The two introns are, however, located at sites different from those of introns in the cob genes of A. nidulans and S. cerevisiae (which contain highly homologous introns at the same site within the gene). The introns share several short regions of sequence homology (10-12 bp long) with each other and with other fungal mitochondrial introns. Moreover, the second intron contains a 50 nucleotide long sequence that is highly homologous with sequences within every ribosomal intron of fungal mitochondria sequenced to date. The conserved sequences may allow the formation of a core secondary structure, which is nearly identical in many mitochondrial introns. The conserved secondary structure may be required for intron splicing. The second intron contains an open reading frame, continuous with the preceding exon, of approximately 290 codons. Two stretches of 10 amino acid residues, conserved in many introns, are present in the open reading frame.     

Genomic organization of the human multidrug resistance (MDR1) gene and origin of P-glycoproteins

The MDR1 gene, responsible for multidrug resistance in human cells, encodes a broad specificity efflux pump (P-glycoprotein). P-glycoprotein consists of two similar halves, each half including a hydrophobic transmembrane region and a nucleotide-binding domain. On the basis of sequence homology between the N-terminal and C-terminal halves of P-glycoprotein, we have previously suggested that this gene arose by duplication of a primordial gene. We have now determined the complete intron/exon structure of the MDR1 gene by direct sequencing of cosmid clones and enzymatic amplification of genomic DNA segments. The MDR1 gene includes 28 introns, 26 of which interrupt the protein-coding sequence. Although both halves of the protein-coding sequence are composed of approximately the same number of exons, only two intron pairs, both within the nucleotide-binding domains, are located at conserved positions in the two halves of the protein. The other introns occur at different locations in the two halves of the protein and in most cases interrupt the coding sequence at different positions relative to the open reading frame. These results suggest that the P-glycoprotein arose by fusion of genes for two related but independently evolved proteins rather than by internal duplication.     

A glimpse of a putative pre-intron phase of eukaryotic evolution

Comparison of the exon-intron structures of ancient eukaryotic paralogs reveals the absence of conserved intron positions in these genes. This is in contrast to the conservation of intron positions in orthologous genes from even the most evolutionarily distant eukaryotes and in more recent paralogs. The lack of conserved intron positions in ancient paralogs probably reflects the origination of these genes during the earliest phase of eukaryotic evolution, which was characterized by concomitant invasion of genes by group II self-splicing elements (which were to become introns in the future) and extensive duplication of genes.     

Intron evolution in a phylogenetic perspective: Divergent trends in the two copies of the duplicated def gene in Impatiens L. (Balsaminaceae)

The history of MADS box genes is well-known in angiosperms. While duplication events and gene losses occur frequently, gene structure and intron positions are very conserved. We investigated all six introns in a duplicated MADS box gene (deficiens, def) in selected Impatiens taxa, thereby assessing intron features. For the first time, our study provides a comparison of molecular changes in all introns of a gene from a phylogenetic perspective. Interestingly, a uniform pattern of molecular evolution in the introns of each copy was not observed, but intron length increases, decreases, and size retention can be found in each copy. A tendency to accumulate long autapomorphic indels is also present, thus, a longer intron length does not reflect a higher number of parsimony-informative characters. Substitution rates vary between introns of each gene copy. While four of the six introns of def1 exhibit a change in their substitution rate, five of the six def2 introns maintain their rates throughout the genus albeit at different levels. In MADS box genes several regulatory sequences are found residing in introns. Thus, presence of putative regulatory motifs was investigated. Most of them are not conserved in position and usually present in only one of the gene copies. In addition, the potential for phylogenetic reconstruction of introns in both def copies is shortly discussed.     

In contrast to the nematode and fruit fly all 9 intron positions of the sea anemone lamin gene are conserved in human lamin genes

We identified the single gene for nuclear lamin in the genome draft of the sea anemone Nematostella vectensis, a member of the cnidaria, a very old metazoan phylum. The gene consists of 10 exons and 9 introns. Strikingly all 9 intron positions are conserved in the human lamin B genes, which have only 1 (lamin B1) or 2 (lamin B2) additional introns. Using the information on neighboring genes we propose that the human lamin B1 gene on chromosome 5 is the true homolog of the Nematostella lamin gene, while the lamin B2 gene on chromosome 19 arose during vertebrate evolution. In marked contrast to this conservation of gene structure are the results in the rapidly evolving genomes of Drosophila and Caenorhabditis elegans. Here the lamin genes have much fewer introns and these occur often at novel positions. In the single nematode lamin gene and the Drosophila lamin Dmo gene no intron position coincides with an intron in the sea anemone lamin gene.     

The genomic locus of the human hemopoietic-specific cell protein tyrosine kinase (PTK)-encoding gene (HCK) confirms conservation of exon-intron structure among human PTKs of the src family.

Protein tyrosine kinases (PTKs) are implicated in the control of cell growth by virtue of their frequent appearance as products of retroviral oncogenes, as intracellular signal transducers, and as growth factor receptors or components thereof. The knowledge of the structure and sequence of family genes encoding PTKs is still limited. To date, the complete genomic structure of human src family members is only available for the C-FGR gene (encoding p55 Fgr, PTK). Sequence analysis and characterization of the intron/exon organization of the human HCK gene, encoding a hemopoietic-specific cell PTK of the src-related family, revealed a length of over 16 kb for the seven 3'-exons. All intron/exon splice junctions agree with the GT/AG rule. In each case where a boundary occurs at a Gly codon, GGG or GGA, the triplet is split between the first and second nucleotide (nt). A total of eight complete and one partial Alu repeats were identified within the introns. The nt sequence of the genomic clones resolves existing discrepancies among two published sequences of HCK cDNAs. Human HCK, C-SRC (encoding p60 Src PTK), C-FGR and LCK (encoding p56 Lck, PTK) genes thus share very similar exon/intron structures for the conserved exons. These results provide additional evidence that the different PTKs of the src-like family most likely arose by duplication of an ancestral src-like gene.     

Structural organization of lower marine nonvertebrate calmodulin genes.

The troponin C (TnC) superfamily genes generally possess five introns, and the positions where they are inserted are well conserved except for the fourth intron. Based on a structural comparison of TnC genes, we proposed that the common ancestor of TnC or TnC superfamily genes had no intron corresponding to the modern fourth intron, and therefore members of the superfamily have gained the fourth intron independently within each lineage. Here, we cloned calmodulin (CaM, one of the members of the TnC superfamily) cDNAs from two lower marine nonvertebrates, the sea anemone, Metridium senile, belonging to the Cnidaria, and the sponge, Halichondria okadai, belonging to the Porifera, and also determined their genomic organization. Chordate CaM genes generally possess five introns, but neither sea anemone nor sponge CaM has anything corresponding to the fourth intron of chordate CaMs, suggesting that the early metazoan CaM must have had only four introns. The modern fourth intron of chordate CaMs was acquired within the chordate lineage after nonvertebrate/chordate divergence. This notion concurs with our proposal explaining the evolution of the TnC superfamily genes.     

Structural organization of the 5' region of the thyroglobulin gene. Evidence for intron loss and "exonization" during evolution

More than one third of thyroglobulin (1190 residues out of 2750) is made of one peptide motif repeated ten times in tandem. Segments unrelated to the motif interrupt this structure at various places. The corresponding gene region, which extends over 40 x 10(3) bases, was studied in detail. All exon borders and exon/intron junctions were localized precisely and sequenced, and their positions were correlated with the repetitive organization of the protein. When intron positions were compiled on a consensus sequence of all repeats, three categories of introns were observed. Except between repeats numbers 5 and 6, an intron was invariably found within the Cys codon making the limit of each motif. This category of intron most probably reflects the serial duplication events responsible for the evolution of this region of the gene. All other introns, except no. 2, are found at positions were the repetitive structure is disrupted by "inserted" peptides. We present the hypothesis that this second category of introns was already present in the original unit before the first duplication. Thereafter, they would have experienced either complete loss (some units do not contain any intron) or partial or total exonization, resulting in the slipping of intronic material into coding sequence. Intron no. 2, finally, separates motif no. 1 at a position on the boundary between two segments presenting sequence homology. This last type of intron probably reflects an initial duplication event at the origin of a primordial thyroglobulin gene motif. With all these characteristics, the thyroglobulin gene is presented as a paradigm for the analysis of the fate of introns in gene evolution.     

Promoter and exon-intron structure of the protein kinase C gene from the marine sponge Geodia cydonium: evolutionary considerations and promoter activity

We report the gene structure of a key signaling molecule from a marine sponge, Geodia cydonium. The selected gene, which codes for a classical protein kinase C (cPKC), comprises 13 exons and 12 introns; the introns are, in contrast to those found in cPKC from higher Metazoa, small in size ranging from 93 nt to 359 nt. The complete gene has a length of 4229 nt and contains exons which encode the characteristic putative regulatory and catalytic domains of metazoan cPKCs. While in the regulatory domain only one intron is in phase 0, in the catalytic domain most introns are phase 0 introns, suggesting that the latter only rarely undergo module duplication. The 5'-flanking sequence of the sponge cPKC gene contains a TATA-box like motif which is located 35-26 nt upstream from the start of the longest sequenced cDNA. This 5'-flanking sequence was analyzed for promoter activity. The longest fragment (538 nt) was able to drive the expression of luciferase in transient transfections of NIH 3T3 fibroblasts; the strong activity of the sponge promoter was found to be half the one displayed by the SV40 reference promoter. Deletion analysis demonstrates that the AP4 site and the GC box which is most adjacent to the TATA box are the crucial elements for maximal promoter activity. The activity of the promoter is not changed in 3T3 cells which are kept serum starved or in the presence of a phorbol ester. In conclusion, these data present the phylogenetically oldest cPKC gene which contains in the 5'-flanking region a promoter functional in the heterologous mammalian cell system.     

Drosophila melanogaster Methoprene-tolerant (Met) gene homologs from three mosquito species: Members of PAS transcriptional factor family

The Methoprene-tolerant (Met) gene in Drosophila melanogaster has been shown to function in juvenile hormone (JH) action. Met homologs were isolated from three mosquito species, Culex pipiens, Aedes aegypti and Anopheles gambiae. Sequence similarity was found to be high in bHLH and PAS conserved domains, and the majority of the 7-9 introns in AaMet and AgMet are located in either identical or similar positions, indicating evolutionary relatedness. Sequence comparison with Met and the similar germ-cell expressed (gce) gene in D. melanogaster showed that the mosquito genes are more similar to gce than to Met. Moreover, the multiple introns in AgMet and AaMet are more similar in number with the 7 introns in Dmgce than to the single intron in DmMet; in fact, six intron positions in AaMet and AgMet are similar to those in Dmgce. Efforts to identify a second homologous gene in mosquitoes were unsuccessful, suggesting a single gene in lower Diptera, consistent with the single gene uncovered in genomic sequencing of Ae. aegypti and An. gambiae. These results suggest that a gene duplication occurred during the evolution of higher Diptera, resulting in Met and gce.