Introduction

Until recent discoveries1,2 it was assumed that the physiological cell death is restricted to multicellular organisms, which have separate germ and somatic cells.3 Originally it was suggested to divide the physiological cell death into two processes: (a) ‘programmed cell death’, describing the developmentally regulated elimination of specific cells during embryogenesis,4 and (b) ‘apoptotis’, describing morphological changes of dying cells.5 While programmed cell death is characterized by the expression of new genes, apoptosis can proceed without new gene expression (reviewed in6). At present, these terms are used interchangeably; therefore, we apply the term apoptosis. In the last two years it has been elucidated that apoptosis is not restricted to Metazoa with separated cell lines but evolved already during the transition from the common ancestor of all metazoan phyla to the phylogenetically oldest metazoan taxon, the sponges [Porifera] (reviewed in7).

Based on the identification and cloning of characteristic metazoan proteins, e.g. cell surface receptors (reviewed in8), extracellular matrix proteins (reviewed in9), immune molecules (reviewed in10), it became evident that all metazoans, including the sponges, are of monophyletic origin.11 This fact suggested also that the key molecules involved in apoptosis are already present in sponges. Subsequent studies elucidated that in sponges, with the two demosponges Geodia cydonium and Suberites domuncula as the most thoroughly investigated examples, polypeptide sequences exist which comprise high sequence similarities to members of the wider Bcl-2 family1,2 as well as to pro-apoptotic molecules containing the death domain.2

Members of the Bcl-2 family are defined by the conservation of two domains, termed Bcl-2 homology domain-1 and -2 (BH1 and BH2).12,13 The Bcl-2 family comprises pro-survival molecules, e.g. Bcl-214 as well as pro-death proteins, e.g. Bax.13 The number of the Bcl-2 family members is rapidly growing.15 Structure analyses revealed that BH1, BH2 and BH3 are involved in binding to other proteins from the Bcl-2 family.16,17 Hence, the members of the Bcl-2 family interact with each other under formation of a complex network of homo- and heterodimers.18 It is assumed that the ratio of pro-survival versus pro-apoptotic dimers is crucial for the resistance of cells to apoptosis.19 While heterodimerization is not required for the pro-survival function20 such a process is thought to be essential for the pro-apoptotic activity, especially via the BH3 region.21 Members of the Bcl-2 family that comprise only domains BH1 and BH2 belong to the pro-survival proteins, e.g. the GRS,22 the A123 and LMW5-HL.24

The first cDNAs encoding putative Bcl-2 homologous proteins from animals evolutionary older than the triploblastic Metazoa have been identified from the sponges G. cydonium and S. domuncula.1 They are characterized by BH1 and BH2 domains and one transmembrane segment. It was the aim of the present study to elucidate the function of the putative Bcl-2 homologous proteins in sponges. For these experiments a Bcl-2 molecule isolated from G. cydonium was chosen. Like the previous sponge Bcl-2's the new molecule, termed GCBHP2 (deduced polypeptide BHP2_GC), comprises the BH1 and the BH2 domains as well as the transmembrane segment.1 Exposure of sponge cells to the organic derivative of the heavy metal tin, tributyltin [TBT] and heat shock causes an increased expression of GCBHP2. TBT also causes an [Ca2+]i overload, cytoskeletal damage and mitochondrial failure (see25); and it was already earlier found that TBT induces apoptosis in sponges26 and also in vertebrate cells.25,27 Here, the molecular marker septin was used to establish the selective expression of GCBHP2. Septin is a molecule involved in the orientation of the cleavage planes during cell division.28 In order to support the finding in the homologous system, HEK-293 cells as well as mouse NIH-3T3 cells were stably transfected with the sponge gene encoding the putative Bcl-2 homologous protein. HEK-293 cells have been selected for the studies since these cells are suitable for the investigation of the Fas-induced apoptotic pathway,29 and hence will also be used for further studies on the role of the sponge Bcl-2 homologous proteins in this metabolism. In addition, transfection studies have been performed with NIH-3T3 cells. It was found that the sponge GCBHP2 gene confers resistance to HEK-293 cells towards the two inducers of apoptosis, TBT and to serum starvation25,30 after transfection. Furthermore, this report describes for the first time that a sponge gene causes a functional expression also in mammalian cells.

Results

Cloning of the Bcl-2 homologous protein GCBHP2 from G. cydonium

The complete cDNA GCBHP2 encoding the putative anti-apoptotic protein BHP2_GC has been isolated after screening the G. cydonium cDNA library with GCBHP1 as a probe. The cDNA for the second Bcl-2 homologous protein (BHP) from G. cydonium is 935 nt long (excluding the poly(A) tail). The potential ORF starting at nt110 to nt112 to the stop codon at nt815–817 encodes a 235 aa long polypeptide (Figure 1A). The deduced protein sequence of this second BHP has a putative size (Mr) of 25 783 and an estimated isoelectric point (pl) of 6.38. The instability index has been computed to be 33.57, indicating that BHP2_GC is a stable protein.31 Northern blot analysis performed with the sponge GCBHP2 clone as a probe yielded one prominent band of approximately 1.2 kb, indicating that the full-length cDNA was isolated (see below). In addition, the exact transcription initiation site for GCBHP2 was determined experimentally by reverse transcription-PCR; it was found to be 109 nts upstream of the putative start-methionine. Like the previously described BHP BHP1_GC also BHP2_GC comprises the two characteristic Bcl-2 protein family domains BH1 (aa125 to aa145) and BH2 (aa174 to aa193) and the COOH-terminal transmembrane segment (aa214 to aa234); Figure 1A. Similarity searches for BH1 and BH2 domains were performed according to Yang and Korsmeyer.32 No similarities to BH3 and BH4 domains could be detected in BHP2_GC.20,33 Therefore, the number of potential helices towards the NH2-terminal end of BHP2_GC was determined. The computer analysis predicted two α-helix regions (Figure 1A) spanning from aa24 to aa46 and from aa94 to aa129.

Figure 1
figure 1

Sponge Bcl-2 homologue proteins and their phylogenetic relationships. (A) Alignments of the two deduced aa sequences of the G. cydonium GCBHP1 (protein BHP1_GC; accession number Y19156), GCBHP2 (BHP2_GC). The alignment was performed using CLUSTAL W program. Similar aa residues are shown in white letters on black background. The Bcl-2 domains, BH1 and BH2, as well as the putative transmembrane segment (TM) are indicated. (B) Phylogenetic relationships of the sponge BH1 domains, present in BHP1_GC and BHP2_GC. (B-1) Alignment of the domains in BHP1_GC and BHP2_GC (BHP1_GC, aa segment 144–165; and BHP2_GC, aa 125–145) with the corresponding domains from the pro-survival proteins (marked: S in B-2 and C-2) from higher metazoan phyla, Ced-9 from C. elegans (CED9_CAEEL; P41958, aa 158–180), NR-13 from quail (NR13_COTJA; Q90343, aa 73–95), frog R11 (AR11_XENLA; Q91828, aa 100–121), human MCL-1 (MCL1_HUMAN; Q07820, aa 251–273), human A1-BFL-1 (BFL1_HUMAN; Q16548, aa 76–98), mouse Bcl-xL (BCLX_MOUSE; AAC53459, aa 128–149), human Bcl-2 (BCL2_HUMAN; P10415, aa 135–156) and Bcl-W (BCLW_HUMAN; Q92843, aa 84–105) and the following apoptotic proteins (marked: A in B-2 and C-2) the mouse Mtd (MTD_MOUSE; AAC53582, aa 112–133), mouse Bax (BAXA_MOUSE; Q07813, aa 97–119), human Bak (BAK_HUMAN; 1096082, aa 116–137), mouse Bad (BAD_MOUSE; A55671; aa 137–159) and the human activator of apoptosis Hrk (HRK_HUMAN; AAC34931, aa 1–20). (B-2) Unrooted tree constructed from sequences listed in B-1 by neighbour-joining as described under Materials and Methods. The numbers at the nodes are an indication of the level of confidence (given in percentage) for the branches as determined by bootstrap analysis [1000 bootstrap replicates]. The scale bar indicates an evolutionary distance of 0.1 aa substitutions per position in the sequence. (C) Phylogenetic relationships of the BH2 domains in the sponge BHP1_GC and BHP2_GC sequences (BHP1_GC, aa 195–216; BHP2_GC, aa 173–193) with the BH2 domains found in the pro-survival proteins Ced-9 from C. elegans (CED9_CAEEL, aa 210–231), NR-13 from quail (NR13_COTJA, aa 123–144), frog R11 (AR11_XENLA, aa 149–170), human Mcl-1 (MCL1_HUMAN, aa 301–322), the mouse Bcl-2 sequences Bcl-xL (BCLX_MOUSE, aa 177–198), human Bcl-W (BCLW_HUMAN, aa 133–154) and human Bcl-2 (BCL2_HUMAN, aa 184–205) as well as human A1-Bfl-1 (BFL1_HUMAN, aa 129–150), and the following apoptotic proteins mouse Bax (BAXA_MOUSE, aa 127–149), mouse Mtd (MTD_MOUSE, aa 161–182), mouse Bad (BAD_MOUSE, aa 172–193), human Bak (BAK_HUMAN, aa 166–186), human apoptosis inducer Bik (BIK_HUMAN; Q13323, aa 118–140), human apoptosis inducer Nbk (NBK_HUMAN; AAC79124, aa 118–140), and rat Bcl-2-related ovarian killer protein Bok (BOK_RAT; AAC61928, aa 118–139). C-1 Alignment and C-2 unrooted tree

The new sponge BHP protein identified has been aligned with the most similar polypeptides listed in the data banks. All belong to the Bcl-2 family comprising the BH1 and BH2 domains (data not shown). Highest similarity was found to the G. cydonium BHP1_GC with 10% of identical aa and 31% of similar aa. Somewhat lower are similarities to the vertebrate members of the Bcl-2 family, both promoting cell survival (pro-survival) and facilitating apoptosis (pro-apoptotic), as well as the Caenorhabditis elegans Ced-9 sequence (≈8% identical aa ≈25% similar aa).

Phylogenetic analysis of the sponge BHP domains

The BH1- (Figure 1B) as well as the BH2 domains (Figure 1C) from BHP2_GC were compared with the corresponding domains found in the most similar sequences of the Bcl-2 family. The analysis revealed that the sponge molecules share (almost) all aa defined to be characteristic for BH1 domains32,33 with other members of the Bcl-2 family, with ≈30–45% of identical aa and ≈50–72% of similar aa. In contrast, the similarity to the C. elegans domain is only 17% (identical aa) and ≈47% (similar aa). The phylogenetic tree (unrooted) constructed from the alignment (Figure 1 B-1) from these BH1 segments, revealed that neither the sponge domains from BHP1_GC and BHP2_GC nor the BH1 domains from the pro-survival and pro-apoptotic sequences can be grouped to distinct clusters (Figure 1 B-2).

In the same way the similarities of the sponge BH2 domains have been analyzed (Figure 1C). Based on the alignment (Figure 1 C-1) the sponge domains share ≈20–30% of identical aa and ≈25–55% of similar aa with the corresponding pro-survival and pro-apoptotic sequences from vertebrates. Again the similarity to the C. elegans domain is low (9% identical aa and ≈18% similar aa). The unrooted tree shows again that no distinction between pro-survival and pro-apoptotic sequences can be made which is based on the aa similarity of the BH2 domains (Figure 1 C-2).

Cloning of the sponge septin

As a first marker for proliferation in sponges, the cDNA GCSEP1, encoding the putative septin SEP1_GC, was cloned; the sequence has been deposited (AJ293509).

Preparation of recombinant Bcl-2 homologous protein and antibodies

GCBHP2 was expressed in E. coli and the purified recombinant oligohistidine-rBHP2GC fusion protein was prepared. The bacteria remained either uninduced or were induced by IPTG. The fusion protein was purified by affinity chromatography using Ni-NTA-agarose resin; the recombinant protein preparation is almost completely pure. A size of 23 kDa was determined for the recombinant protein. The calculated Mr for the deduced aa sequence of the cDNA BHP2_GC (inclusive the histidine tag) is 23.5 kDa (not shown). Polyclonal antibodies have been prepared against rBHP2GC. Western blot analysis revealed only one band corresponding to rBHP2GC. Control studies were performed with antibodies which had been adsorbed with rBHP2GC; under those conditions no staining was seen (not shown).

Exposure of cells from G. cydonium to selected stressors

Two stressors have been chosen to study their effect on the expression of GCBHP2 in the homologous system: tributyltin (TBT) and heat shock.

Tissue samples were exposed to 1 μM of TBT for up to 8 h. Subsequently, RNA was prepared from the samples and subjected to Northern blotting to estimate the levels of gene expression. It was found that the steady-state level of septin (size: 1.4 kb), which was chosen as a reference gene to monitor the viability and proliferation of the cells, did not alter, while the expression levels for the genes encoding the putative anti-apoptotic protein BHP2_GC GCBHP2 (1.2 kb) and the heat shock protein-70 GCHSP70 (accession number X94985; 2.1 kb) increased drastically from 0% at time zero (not detectable under the conditions used) to 100% after incubation for 8 h in TBT (Figure 2A). At concentrations above 10 μM of TBT the induction of the two ‘stress-responsive’ genes was abolished (not shown).

Figure 2
figure 2

Effect of two stressors on tissue samples from G. cydonium. (A) Effect of TBT on the expression of GCBHP2 (encoding for GCBHP2), GCHSP70 (HSP-70) and GCSEP1 (septin). A concentration of 1 μM was chosen. Tissue samples were taken after an incubation period of 0–8 h (as indicated). Then RNA was prepared and subjected to Northern blot analysis. 5 μg of RNA per slot were loaded onto the gel. After separation and blot transfer, hybridization was performed as described under Materials and Methods. The signals were visualized applying the chemiluminescence procedure. The intensities of the transcripts for GCBHP2, GCHSP70 and GCSEP1 are correlated with the highest level of expression for the respective gene (this value was set to 100%). To assure that the same amount of RNA was loaded, the gel was stained with toluidine blue to visualize the rRNA. M: size markers used. (B) Effect of heat shock on the expression of the sponge genes GCBHP2, GCHSP70 and GCSEP1 in tissue samples. The tissue samples were taken 0–5 h after heat shock, as described under Materials and Methods. RNA was isolated, blot transferred and probed with GCBHP2, GCHSP70 and GCSEP1. (C) Effect of TBT on BHP2_GC protein expression in G. cydonium tissue. Samples were treated with TBT (1 μM) for 0–8 h, protein was extracted and subjected to Western blot analysis. Equal amounts of protein were loaded onto the gel. After separation and blot transfer, the blots were incubated with PoAb-BHP2GC and processed as described under Materials and Methods

In a parallel line of experiments the tissue samples were treated with heat at 26°C (9°C above the ambient temperature) as described under Materials and Methods. It was found again that the expression of the 1.2 kb GCBHP2 transcript and the 2.1 kb GCHSP70 transcript are significantly altered. While the expression of septin (1.4 kb) remained unchanged, the levels of GCBHP2 and GCHSP70 expression strongly increased from 0% (time zero) to 100% after an incubation for 5 h under heat stress (Figure 2B). If the tissue samples were treated at higher temperature stress, e.g. 14°C above the ambient temperature for 5 h, no induction of gene expression was measured (not shown). These data show that the two stressors TBT and heat shock cause a selective expression of the genes encoding the stress proteins BHP2_GC and heat shock protein-70.

The degree of expression of GCBHP2 was also monitored on protein level by the Western blot technique. Tissue samples were treated both with TBT or heat as the protocols describe above. Both stressors caused an increase of BHP2_GC as checked. The data of the TBT induced expression are documented in Figure 2C. It is evident that the amount of BHP2_GC strongly increased after 1 h of incubation. A prolonged incubation for 8 h even enhanced the degree of expression as can be deduced from the strength of the 26 kDa band [calculated size of the sponge protein: Mr 25 783; see above] reflecting the BHP2_GC protein.

Transfection of sponge GCBHP2 into mammalian cells

The cDNA construct with the G. cydonium GCBHP2 was introduced into human HEK-293 cells. Stable transfectants were selected by geneticine, G418, and the expression of GCBHP2 was confirmed by the following techniques.

Identification of the transcripts

Using the primers described under Materials and Methods, the 590 bp long PCR product, reflecting the expression of GCBHP2, could be identified in transfected cells; in mock-transfected cells this band was missing. As a control, the expression of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene was determined; the expected 471-bp long PCR fragment was obtained both in transfected and non-transfected cells (data not shown).

Western blotting

Antibodies, raised against the recombinant G. cydonium BHP2_GC, rBHP2GC, were applied to demonstrate the presence of the expressed sponge gene in transfected HEK-293 cells. The Western blot analysis of the extracts from transfected cells revealed one band corresponding to a size of 26 kDa (data not shown). Control studies have been performed with an extract from mock-transfected cells. In this sample no protein was recognized by the PoAb-BHP2GC antibody.

BHP2_GC expression in HEK-293 cells and NIH-3T3 cells

GCBHP2 was expressed in HEK-293 cells. Stable transfectants changed their pattern of attachment to the solid phase, irrespectively if the cells were plated on plastics or glass. Control cells, transfected with the pcDNA3 plasmid without insert were shown to spread homogeneously over the surface of the plates. In contrast, the HEK-293 cells transfected with the plasmid containing GCBHP2 clump together and form ‘foci’-like clusters (data not shown).

The expression of GCBHP2 was demonstrated by immunofluorescence. HEK-293 cells which had been stably transfected with the expression plasmid lacking GCBHP2 did not show any staining with the antibody, PoAb-BHP2GC, raised against the recombinant rBHP2GC (Figure 3a,b and e,f). However cells which had been stably transfected with the expression plasmid containing GCBHP2 are brightly stained using PoAb-BHP2GC (Figure 3c,d and g,h).

Figure 3
figure 3

Identification of the expressed G. cydonium GCBHP2 in HEK-293 cells by immunofluorescence studies. Cells were stably transfected using the expression plasmid lacking GCBHP2 (a and b; e and f) or with the plasmid containing GCBHP2 as an insert (c and d; g and h). After fixation the cells were permeabilized and stained with PoAb-BHP2GC followed by incubation with Cy3-labeled anti-rabbit Ig. Inspection with Nomarsky: a, c, e and g; analysis by fluorescence microscopy: b, d, f and h. Magnification: a–d, ×25; and e–f ×150

In a parallel series of experiments NIH-3T3 cells were stably transfected with GCBHP2. Subsequently, the cells were stained both with the monoclonal anti-Hexa-His antibody and with the affinity purified polyclonal antibody against p33/gC1qR (a mitochondrial marker protein), to obtain a more detailed intracellular location in cells larger than HEK-293. The results show that the distribution of the staining with the PoAb-BHP2GC and with anti-p33/gC1qR is almost identical (Figure 4); however a distinct localization of BHP2_GC on the surface of the mitochondria can not be deduced as reported earlier (see34).

Figure 4
figure 4

Immunofluorescence analysis of expressed sponge GCBHP2 in mouse NIH-3T3 cells. These cells were stably transfected with GCBHP2 and subsequently stained with the anti-Hexa-His antibody [polyclonal] (a) and with antibody against a mitochondrial protein, p33/gC1qR [monoclonal] (b) Hexa-His antigen/antibody complexes were identified with Cy3-conjugated secondary antibody and the p33/gC1qR antigen/antibody complexes with an FITC-labeled secondary antibody. Magnification: ×600

Resistance of HEK-293 cells against stressors after transfection with GCBHP2

HEK-293 cells stably transfected with the sponge GCBHP2 have been subjected to two distinct stresses: serum starvation and TBT treatment. The experiments have been performed under conditions described earlier for HEK-293 cells transfected with human Bcl-2 and induced to apoptosis using serum deprivation.34

Serum starvation

In medium, containing no FCS, the viability of the mock-transfected cells was strongly reduced; after 96 h only ≈10% of the cells were still viable. However, the number of viable HEK-293 cells transfected with GCBHP2 was significantly enhanced; after an incubation period of 96 h the per cent of surviving cells was 46% with respect to the number at the beginning of the experiment. After a period of 72 h, 65% of the GCBHP2-transfected cells survived compared to 15% in the assays with mock-transfected cells (Figure 5). These values are almost similar to those reported for HEK-293 cells transfected with human Bcl-2;34 in this report it was found that 74% of the cells transfected with the human gene survived, compared with 17% of mock-transfected cells.

Figure 5
figure 5

Increase of survival of HEK-293 cells, transfected with GCBHP2. The cells were kept under serum-starved conditions and incubated for 0–96 h. After the indicated periods the concentration of viable cells per square mm was determined in GCBHP2-transfected cultures (closed bars) and in mock-transfected ones (open bars). Further details are given under Materials and Methods

TBT treatment

Likewise an increased viability could be determined when GCBHP2-transfected cells were treated with 0.25 μM of TBT. As under serum starvation, the concentration of viable cells in transfected cells was markedly increased. After an incubation for 96 h the per cent of surviving cells was 32% for GCBHP2-transfected cells, compared to 5% of the mock-transfected cells. During a period of 72 h, 52% of the GCBHP2-transfected cells were alive compared to 12% in the assays with mock-transfected cells (data not shown).

In order to demonstrate that the level of GCBHP2-expression in HEK-293 cells was not altered after treatment with the stressors, the cell extracts were analyzed for the expression of the sponge gene. The extracts were subjected to Western blot analysis and reacted with the PoAb–BHP2GC antibody. The 26 kDa signal, corresponding to the sponge (recombinant) BHP2_GC protein, was found not to be significantly changed during the 92 h incubation period in the assays under serum starvation or in cells after TBT treatment (data not shown).

Expression of human Bcl-2 in transfected cells under stress conditions

In parallel, the expression of the cellular Bcl-2 gene was monitored on protein level by Western blotting. The results revealed that both in non-treated cells (time zero) as well as in BHP2_GC-transfected cells the Bcl-2 protein (29 kDa) is present at almost the same level, irrespectively of the presence of serum in the culture medium (Figure 6).

Figure 6
figure 6

Caspase activation and Bcl-2 expression in HEK-293 cells in the absence or presence of serum. Above Time course of caspase activation in mock-transfected (open bars) and GCBHP2-transfected HEK-293 cells (closed bars). HEK-293 cells were incubated in the absence of FCS for 0–92 h. Then the activity of caspase-3 was determined as described under Materials and Methods; the control values were set to 100%. The asterisks indicate the statistical significance (P<0.005) of the difference between the activities measured in the extracts from mock-transfected and those from GCBHP2-transfected HEK-293 cells; n=20. Below Expression of Bcl-2 in GCBHP2-transfected HEK-293. As indicated, the cells were incubated for 0–92 h in the presence or absence of serum as outlined under Materials and Methods. Then total protein was extracted and analyzed for the level of Bcl-2 protein by Western blotting. The positions of marker proteins (M) were given

Effect of the stressors on the activation of caspase-3

Quantitation of caspase-3 activity in HEK-293 cells was performed by measuring cleavage of the fluorescent peptide AC-DEVD-AMC [Ac-Asp-Glu-Val-Asp-(7-amino-4-methylcoumarin). This peptide mimics the target sequence of the group II effector caspases.35 Treatment of mock-transfected HEK-293 under serum starvation resulted in a significant increase in caspase activity above control levels within the first 20 h of incubation. As shown in Figure 6, the increase of caspase activity in BHP2_GC-transfected cells is significantly lower compared to mock-transfected HEK-293 cells. Similar results were observed with BHP2_GC-transfected cells and mock-transfected HEK-293 cells which had been treated with TBT (data not shown). Therefore, we conclude that BHP2_GC-transfected HEK-293 cells are significantly more resistant towards the apoptotic stressors serum starvation and TBT than mock-transfected cells.

Discussion

One key pathway in the development and homeostasis of Metazoa is the active form of cell death, apoptosis. This form of cell death allows the organism to remove unwanted and/or damaged cells that could cause inflammation. In Metazoa, the apoptotic pathway(s) have been described, until recently, only in phyla evolved from the hypothetical ancestor of the Protostomia and Deuterostomia, the Urbilateria (reviewed in36,37). As model systems for Protostomia, the nematode Caenorhabditis elegans and for Deuterostomia the mammalian systems, e.g. Homo sapiens (reviewed in15,36), have been well described. With the recent description of the presence of pro-apoptotic as well as pro-survival molecules,1,2 in sponges, the closest relative to the common ancestral phylum of all Metazoa, the Urmetazoa,38 a new avenue to the understanding of basic apoptotic pathways in all Metazoa phyla became possible.

One suitable approach to prove that the sponge molecules are true homologues to the nucleotide sequence-related molecules identified in other Metazoan phyla, is the performance of transfection studies in heterologous cell systems. In the present work the gene from G. cydonium encoding the putative pro-survival protein BHP2_GC, GCBHP2 was stably transfected into human embryonic kidney cells HEK-293. The deduced polypeptide of BHP2_GC comprises the two characteristic Bcl-2 protein family domains BH1 and BH2. These domains are known to be important for Bcl-2 to form heterodimeric complexes with the family members and to carry out the anti-apoptotic function (reviewed in16). Mutations in BH1 and BH2 prevent Bcl-2 from the formation of heterodimers with Bax, a homologue of Bcl-2, resulting in an abrogation of the Bcl-2-mediated cell survival.16

In the first series of experiments the expression of the sponge Bcl-2-related molecule BHP2_GC in the sponge tissue was studied in response to defined stressors. The capacity of inducers to cause apoptosis depends on the respective species, tissues or cells (reviewed in39). Therefore, in the present study the established stimuli for apoptosis in sponges, TBT26 and heat shock40 have been chosen. With respect to HEK-293 cells serum starvation and TBT, both triggering the cells to apoptosis dose-dependently via the caspase pathway,25,30 have been selected as inducers.

As a measure of apoptosis the TUNEL assay had previously been used in the sponge system.41 Considering the limitations of this assay, functional studies have been applied here. Until now the searches for caspases in sponges have been without success. Therefore, as a molecular marker for viability and proliferation the expression of septin, which is involved in the orientation of the cleavage planes during cell division,28 has been chosen.42,43 The experiments with tissue from G. cydonium revealed that in response to 1 μM of TBT and heat shock (1–5 h at 9°C above ambient temperature), the expression of GCBHP2 is upregulated with almost the same kinetics like the gene GCHSP70, encoding the sponge heat shock protein-70. Under the conditions chosen, the level of expression for the septin gene remains unchanged. These data show that the sponge Bcl-2 homologue GCBHP2 is, like the mammalian Bcl-2 gene (e.g.44), inducible; it is selectively upregulated in response to non-toxic apoptotic stimuli. At higher TBT concentrations (>10 μM) or more severe heat shock (14°C above ambient temperature) the expressions of genes coding for BHP2_GC, HSP-70 and septin decline26,45 (and unpublished).

Transfection of HEK-293 cells and NIH-3T3 cells with the sponge GCBHP2 was successfully performed as shown both on the transcriptional as well as on the protein level. The GCBHP2 transcripts have been identified by PCR analysis, while the protein expression in the transfected cells was demonstrated by immunofluorescence using antibodies raised against the recombinant BHP2_GC protein. The transfected cells showed a different phenotype, a phenomenon which is known already from previous studies in which HEK-293 cells were transfected with the Kin17 gene.46 While the non-transfected cells show the well established evenly distributed growth pattern, the transfected cells clumped together. The staining in transfected cells using an antibody against the sponge BHP2_GC protein and a mitochondrial marker protein, p33/gC1qR, gave almost identical patterns. This result suggests that the sponge protein in mammalian cells is co-localized with mitochondria. The association of Bcl-2 with mitochondria in mammalian cells is established (reviewed in37).

The expression of the sponge GCBHP2 in HEK-293 cells remained at a high level even after a 92 h exposure period in medium completely deprived of fetal calf serum. The activation of a caspase, we have selected caspase-3 since this enzyme is involved in the cleavage of Bcl-XL,47 has been measured in serum starved HEK-293 cells. The activity of the enzyme was found to be strongly increased in mock-transfected cells; a considerably lower increase in caspase activity was measured in GCBHP2-transfected cells. In the final series of experiments it was established that cells exposed to medium in the absence of FCS displayed a significantly increased survival rate; as an example, after 72 h 65% of the GCBHP2-transfected cells were found to be still alive compared to 15% in the assays with mock-transfected cells. This degree of protection against the apoptotic stimulus is almost as high as that seen in HEK-293 cells transfected with the human Bcl-2 gene.34 To verify that HEK-293 cells transfected with the sponge gene are also significantly more resistant to other apoptotic stimuli, the cells were treated with TBT. Again, the GCBHP2-transfected cells displayed a much higher survival rate during the 96 h incubation period compared to the mock-transfected cells.

The data reported here demonstrate that the tools of sponge cultivation48 and molecular biology are powerful means to identify components of complex physiological pathways in the phylogenetically oldest metazoan phylum, in sponges. It can now be concluded that these animals comprise also the characteristic metazoan pathway of apoptosis which had been considered until recently to have emerged only in animals evolved from the common ancestor of Protostomia and Deuterostomia. This conclusion is of considerable importance for the understanding of evolution, since yeast, the kingdom which has a common ancestor with Metazoa,49 shows the phenomenon of apoptosis without having the metazoan apoptotic/anti-apoptotic genes.39 Hence it must be concluded (as shown here for BHP2_GC, a member of the Bcl-2 protein family), that key molecules of apoptosis emerged during evolution of the Fungi to the Metazoa. In a subsequent study the function of the sponge pro-apoptotic molecules comprising two death domains2 in mammalian cells will be investigated. With respect to the Bcl-2 protein family it can now be stated that this apoptosis regulator family is an autapomorphic character of Metazoa. In addition, this study is also the first to demonstrate that sponge gene(s) can display their function in human cells, again supporting the fact that sponges, as the simplest metazoan animals, share common characteristics with other metazoan taxa. This monophyletic view of animals has been previously proposed based on cDNA sequence data.11

Materials and Methods

Materials

Restriction endonucleases and other enzymes for recombinant DNA techniques and vectors were obtained as described earlier.50,51 In addition, FuGENE 6 reagent, PCR-DIG (digoxigenin)-Probe-Synthesis Kit, anti-DIG AP Fab fragments and CDP [disodium 2-chloro-5-(4-methoxyspiro{1,2-dioxetane-3,2′-(5′-chloro)-tricyclo[3.3.1.1.3,7]decan}-4-yl)phenyl phosphate] from Roche Diagnostics (Mannheim; Germany); RPMI-1640 and OptiMEM from Life Technologies (Karlsruhe; Germany); geneticine G418 from Invitrogen (Groningen; The Netherlands); tributyltin and labelled anti-rabbit Ig from Sigma (Deisenhofen; Germany); AC-DEVD-AMC from Peptide Institute (Osaka; Japan); monoclonal anti-Bcl-2 [human] antibodies (mouse) from Santa Cruz (Santa Cruz, Cal; USA); anti-mouse (Cy3 labeled) and anti-rabbit (FITC-labeled) from Dianova (Hamburg; Germany); the mouse anti-Hexa-histidine monoclonal antibody (McAb-anti-His) from Invitrogen (Groningen; The Netherlands). The mammalian cell line was purchased from the American Type Culture Collection (Manassas, VA; USA). The antibody raised against the mitochondrial marker protein p33/gC1qR was a gift of Drs Müller-Esterl and Dedio.52

Sponges

Specimens of the marine sponge Geodia cydonium (Porifera, Demospongiae, Geodiidae) were collected in the Northern Adriatic near Rovinj (Croatia), and then kept in aquaria in Mainz (Germany) at a temperature of 17°C.

Exposure of tissue cubes from G. cydonium to selected stressors

Tissue samples, cubes of ≈0.05 g of G. cydonium remained either untreated in 10 l aquaria at 17°C (ambient temperature) under continuous aeration or were exposed to defined stress conditions. The water was changed once daily. For the heat shock experiments the cubes were treated for 1–5 h at 26°C. Exposure to tributyltin (TBT) was performed at a concentration of 1 μM for 1–8 h. TBT was dissolved in ethanol (final concentration in the cultures 0.01%); the control samples were supplemented without the compound but in the presence of the solvent. For the experiments aliquots of the respective sponge samples were taken and immediately frozen in liquid nitrogen and stored at −80°C.

Cloning of the putative sponge pro-survival cDNA encoding BHP2_GC

The complete sponge cDNA GCBHP2, encoding the putative Bcl-2 homologous proteins (BHP2_GC), was cloned by screening the cDNA library from G. cydonium53 using the GCBHP1 (accession number Y191561) as a probe. Screening of the library was performed under low stringency hybridization as described before.50 Positive clones were detected with an alkaline phosphatase conjugated anti-digoxygenin antibody using BCIP/NBT as substrate.54 All cDNAs have been obtained at least from two different cDNA libraries resulting in at least three independent clones each. DNA sequencing was performed with an automatic DNA sequenator [Li-Cor 4000S].

Cloning of the putative sponge septin

The sponge cDNA, encoding the septin-like molecule, termed GCSEP1, was isolated from the G. cydonium cDNA library by polymerase chain reaction (PCR). The degenerate reverse primer, directed against the conserved aa segments found in the P loop motif of the Drosophila melanogaster pnut sequence (accession number U08103;55 aa148 to aa157) 5′-TTT/CICCIAGICCICAT/CTCICCIACIACCAT-3′, was used in conjunction with the 5′-end vector-specific primer. Further data will be given elsewhere. The clone was termed GCSEP1.

Sequence comparisons

The sequences were analyzed using computer programs BLAST56 and FASTA.57 Multiple alignments were performed with CLUSTAL W Ver. 1.6.58 Phylogenetic trees were constructed on the basis of aa sequence alignments by neighbour-joining, as implemented in the ‘Neighbor’ program from the PHYLIP package.59 The distance matrices were calculated using the Dayhoff PAM matrix model as described.60 The degree of support for internal branches was further assessed by bootstrapping.59 The graphic presentations were prepared with GeneDoc.61 The hydropathy values were calculated according to Kyte and Doolittle.62

Cell transfection

HEK-293

The sponge GCBHP2 cDNA was transfected into the human embryonal kidney cells HEK-293 (ATCC CRL 1573). These cells were grown as monolayer in a humidified 5% CO2 atmosphere at 37°C in RPMI-1640 medium supplemented with 10% fetal calf serum (v/v; FCS).

A fragment of 771 bp (inclusive the His-tag) from GCBHP2, comprising the complete deduced open reading frame (ORF) was obtained by PCR from nt−23 to nt−42 (forward primer: 5′-GGTACCTCACGACGGTCGGCGTATTA-3′; the KpnI site is underlined), located upstream from the start ATG site, to nt676 to nt696 (reverse primer: 5′-TCTAGATCAGTGATGGTGATGGTGATGTTGTCTCATTCCAAGAGCCAG-3′; the XbaI site is underlined and the six histidine codons are in bold type, TCA is the stop codon). After digestion with KpnI and XbaI the fragment was subcloned into the KpnI and XbaI sites of the plasmid pcDNA3 (Invitrogen). Cells were transfected using the FuGENE 6 reagent following the procedure described by Mansour et al.63 using the pcDNA3 containing the GCBHP2 in transcriptional orientation. As a control, mock-transfection was performed with the pcDNA3 without the GCBHP2 insert.

Cells having reached 50–60% confluency were transfected with 2.8 μg/ml of pcDNA3-GCBHP2 or pcDNA3 without insert. Immediately prior transfection cells were transferred into OptiMEM, supplemented with 10% FCS. After 24 h of incubation of the cells with DNA/FuGENE 6 the transfection medium was changed back to the medium routinely used. To select for stably transfected cells, which had acquired neomycin resistance, G418 (500 μg/ml) was added to the culture medium after an additional 24–48 h incubation period. After selection with G418 in 96-well plates, cloned tranfectants were maintained in RPMI 1640, containing 10% FCS. For microscopical inspection the cells were plated onto cover slides. Expression of GCBHP2 was confirmed by Northern blotting, Western blotting and by immunocytochemistry.

For the analysis of the expression of GCBHP2 in HEK-293 cells the method of reverse transcription-PCR64 was applied. RNA was extracted64 and the sponge GCBHP2 transcripts were identified using the primers (forward primer: nt−10 to nt11; reverse primer: nt560 to nt580). As a control, the level of expression of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was determined. Primers (forward primer: 5′-AATGCCTCCTGCACCACCAA-3′ directed against nt445 to nt464; the reverse primer: 5′-GTCGTTGAGGGCAATGCCAGC-3′ against nt895 to nt915 of the human GAPDH (accession number NM 002046) were used to identify the transcripts of GAPDH in mock-transfected and GCBHP2-transfected HEK-293 cells. The PCR was carried out as follows: denaturation at 95°C (3 min), 35 amplification cycles (95°C, 30 s; 58°C 45 s, 72°C 1.5 min) and final extension (72°C, 10 min). The PCR products were size separated on an agarose gel and stained with ethidium bromide.

NIH-3T3

NIH-3T3 (ATCC CRL 1658) cells were transfected with pcDNA3–GCBHP2 as described above, with the exception that 250 μg/ml of G418 were used for selection. These cells were used for additional immunofluorescence studies to determine the intracellular distribution of the sponge protein.

Treatment of HEK-293 cells with two stressors

Stably transfected or mock-transfected HEK-293 cells were treated in two ways to induce apoptosis. (1) Serum starvation: The cells were cultured in RPMI 1640 medium in the absence of FCS. After treatment for 0–96 h cell survival was measured. For the determination of viability the trypan blue-exclusion assay has been applied. Under these conditions the viability in the mock-transfected cells was determined to be ≈10% after a 96 h-incubation period. (2) TBT treatment: The cells were exposed to 0.25 μM of TBT for 0 to 96 h. Subsequently, cell viability was determined. After 96 h the viability in the TBT-treated cultures was ≈5%.

Northern blot analysis

RNA was extracted from liquid-nitrogen pulverized sponge tissue as described.65,66 RNA (total) samples (5 μg) were electrophoresed in a formaldehyde/1% agarose gel and blotted onto Hybond N+ membrane following the manufacturer's instructions (Amersham; Little Chalfont, Buckinghamshire; UK). Hybridization was performed with the 586 nt long GCBHP2 fragment (spanning nt83 to nt668 of the ORF) in High-NaDodSO4 hybridization Buffer (7% NaDodSO4, 50% formamide, 5×SSC, 0.1% Na-laurylsulfate, 50 mM Na3PO4 [pH 7.0] and 2% of blocking reagent [Roche]) at 42°C. The probe had been labeled using the PCR-DIG-Probe Synthesis kit (Roche) according to the ‘Instruction Manual’ (Roche), applying the primers corresponding to nt83 to nt103 and nt648 to nt668. In parallel, the blots were hybridized with the 762 bp long heat shock protein-70 probe from G. cydonium (nt303 to nt1064 GCHSP70; accession number X9498540), or the 799 bp long G. cydonium septin probe (nt−3 to nt796 GCSEP1 cloned in the present study). The signals of the Northern blot were visualized by the chemiluminescence procedure;67 CDP-Star was used as substrate. To quantitate the signals of the Northern blots the screen was scanned with the GS-525 Molecular Imager (Bio-Rad). The relative values for the expressions are given. To document that the same amount of RNA was loaded onto the gels the total RNA was stained with toluidine blue; using this technique both the eukaryotic [28S and 18S rRNA] as well as the prokaryotic rRNA species [23S and 16S] became visible. It is known that the majority of demosponges, including G. cydonium, harbors bacteria.68

GCBHP2 cDNA expression and antibody preparation

A fragment of GCBHP2 ranging from nt83 to nt659 has been used for the expression of the protein; a forward primer, 5′-GGATCCTGATGGAGATGGAGGAGCTCTACAGAA-3′, comprising nt83 to nt103 of the coding region (the BamHI site is underlined); and a reverse primer, 5′-GTCGACATGACACCAATACCGGCTACT-3′; nt639 to nt659 (the Sa/I site is underlined) was chosen for the PCR. This cDNA was used for expression in Escherichia coli as described.53,69 The cDNA was inserted into the bacterial oligohistidine expression vector pQE-32 (Quiagen). E. coli (XL1-blue) were transformed with this plasmid and expression of fusion protein was induced with isopropyl 1-thio-β-D-galactopyranoside (IPTG). Bacteria were extracted with phosphate buffered saline (PBS)/urea and the suspension was centrifuged; the supernatant was collected ‘bacterial crude extract’. The purification of the recombinant oligohistidine-BHP2GC fusion proteins, termed rBHP2GC, was performed by metal-chelate affinity chromatography using Ni-NTA-agarose resin (Qiagen) according to Hochuli et al.70 as well as by the ‘Instruction of Manufacturer’. The purity of the material was checked by 15% polyacrylamide gels containing 0.1% NaDodSO4 [PAGE] according to Laemmli.71

Polyclonal antibodies against rBHP2GC were raised in rabbits as described.72 Purified rBHP2GC (3×10 μg of protein) was injected; after 15-weeks serum was collected and the antibodies were prepared;73 they are termed PoAb-BHP2GC. In control experiments 100 μl of the PoAb-BHP2GC were adsorbed to 50 μg of rBHP2GC (30 min; 4°C) prior to its use.

Western blot analysis

HEK-293 cells or sponge tissue were lysed in 50 mM Tris-HCl (pH 7.5; 150 mM NaCl, 1% Triton, aprotinin [0.10 U/ml), 20 mM leupeptin and 1 mM phenylmethylsulfonyl fluoride). Lysates were separated by NaDodSO4–PAGE (15% gel) and transferred to a nitrocellulose membrane.74 After blocking with Blocking Reagent (Roche) for 1 h, the membrane was incubated with the primary (PoAb-BHP2GC [1 : 10 000 dilution]), then with the secondary antibody (anti-rabbit Ig [1 : 10 000 dilution]) for 1 h each, and finally developed with enhanced chemiluminesence (Amersham).

Immunofluorescence analysis

Transfected cells were washed twice in PBS, fixed and permeabilized with methanol/0.02% (w/v) EGTA for 10 min. After a further washing procedure (three times) the cells were incubated with PoAb–BHP2GC (1 : 1000 dilution in Blocking Reagent) for 60 min at 37°C.74 The cells were again washed three times with PBS and incubated with Cy3-conjugated goat anti-rabbit Ig (1 : 100 dilution in Blocking Reagent) for 30 min 37°C. The cells were embedded into PBS/glycerol (1 : 1; v/v) and inspected with an Olympus AHBT3 microscope.

In a parallel series of experiments, NIH-3T3 cells transfected with pcDNA3–GCBHP2, were stained with an anti-Hexa-His antibody (1 : 200 dilution) and, in addition, with an antibody raised against the mitochondrial marker protein anti-p33/gC1qR (polyclonal antibodies [rabbits]; 1 : 20 dilution).52 The antigen/antibody immunocomplexes were detected either with anti-mouse (Cy3 labeled [identification of the His-labeled BHP2GC]) or with anti-rabbit (FITC labeled [anti-p33/gC1qR]) antibodies.

Caspase-3 assay

For the testing of caspase-3 activity HEK-293 cells were lysed in 10 mM HEPES–KOH [pH 7.4], 2 mM ethylenediaminetetra-acetic acid, 5 mM dithiothreitol (DTT), 1 mM phenylmethylsulfonyl fluoride (PMSF) and 0.1% CHAPS. Subsequently, 50 μl aliquots were added to the 50 μl caspase assay which contained the following components (final concentrations): 20 mM HEPES–KOH [pH 7.4], 10% glycerol, 0.5 mM PMSF, 2 mM DTT and 20 μM of the caspase-3 substrate Ac-DEVD-AMC. Substrate cleavage to release free AMC (excitation at 355 nm; emission at 460 nm) was monitored at 37°C (after 3 h). Fluorescence units measured were normalized with respect to protein concentration. The values measured for the non-treated cells, either mock-transfected or GCBHP2 transfected, at the indicated time were set to 100%. In parallel, Western blot experiments were performed with the same extracts; total tissue extracts (10 μg/lane) were subjected to electrophoresis in 15% polyacrylamide gels containing 0.1% NaDodSO4 as described above. After electrophoresis the proteins were transferred to a nitrocellulose membrane and incubated with anti-Bcl-2 antibodies [1 : 2500 dilution]) and then with the secondary antibody (anti-mouse Ig [1 : 10 000 dilution]).

Further analytical procedure

For protein determination the Fluoram method was used the standard was bovine serum albumin.75

GenBank information

The sequence reported here is deposited in the EMBL/GenBank data base (Accession no. AJ293508 for the Geodia cydonium Bcl-2 homologous protein BHP2 and AJ293509 for the G. cydonium septin).