Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2005, p. 6934-6946, Vol. 71, No. 11
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.11.6934-6946.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
B Pathway in the Beneficial Euprymna scolopes-Vibrio fischeri Light Organ Symbiosis
Joshua V. Troll,1
Todd E. Scheetz,2
Thomas L. Casavant,2,3
M. Bento Soares,4 and
Margaret J. McFall-Ngai1*
Department of Medical Microbiology and Immunology, University of Wisconsin, Madison, Wisconsin 53706,1 Center for Bioinformatics and Computational Biology and Departments of Biomedical Engineering and Opthalmology,2 Department of Electrical and Computer Engineering,3 Department of Pediatrics, University of Iowa, Iowa City, Iowa 522424
Received 1 April 2005/ Accepted 1 July 2005
|
|
|---|
B pathway is a common, evolutionarily conserved innate immune pathway that modulates the responses of animal cells to microbe-associated molecular patterns (MAMPs). Because MAMPs have been implicated as critical elements in the signaling of symbiont-induced development, an expressed sequence tag library from the juvenile light organ of Euprymna scolopes was used to identify members of the Toll/NF-
B pathway. Full-length transcripts were identified by using 5' and 3' RACE PCR. Seven transcripts critical for MAMP-induced triggering of the Toll/NF-
B phosphorylation cascade have been identified, including receptors, signal transducers, and a transcription factor. Further investigations should elucidate the role of the Toll/NF-
B pathway in the initiation of the beneficial symbiosis between E. scolopes and Vibrio fischeri. |
|
|---|
B pathway, that have striking similarities among these phylogenetically diverse organisms (Fig. 1) (5, 30, 51). The highly conserved Toll/NF-
B signaling pathway has been studied extensively for its role in microbial pathogenesis in two of the three major branches of the animal kingdom, the vertebrates (Deuterostomia) and the insects (Ecdysozoa) (1, 21, 64). The NF-
B family proteins exist in unstimulated cells of these animals as homo- or heterodimers that are sequestered in the cytoplasm by a member of the Inhibitor of NF-
B (I
B) proteins. Through interaction with a receptor, MAMPs elicit an intracellular phosphorylation cascade. This phosphorylation cascade converges on, and activates, a complex of serine kinases termed I
B kinases (IKKs) (35). The IKK complex phosphorylates serine residues of I
B, leading to its ubiquitination and degradation by the proteosome (25). Degradation of I
B allows the NF-
B dimer to translocate into the nucleus and activate transcription of target genes possessing the
B promoter site. Although evidence is available for the presence of the NF-
B pathway in a mollusk, the Pacific oyster (20, 28, 48), the signaling pathways activating the immune responses in mollusks and other Lophotrochozoa, the members of which form the third branch of the animal kingdom (2, 29, 59), are largely unknown (Fig. 1).
![]() View larger version (40K): [in a new window] |
FIG. 1. The signaling pathways activating the innate immune responses of plants and animals are similar. NF- B p105 and Relish have a REL homology domain located at their N terminus (N) and an I B-like domain at their C terminus. The C terminus domain can be proteolytically cleaved and degraded. Brackets enclose pathway members identified in the present study. The molecules listed in the figure represent full-length sequences. Ai, Argopecten irradians; At, Arabidopsis thaliana; C, C terminus; Cg, Crassostrea gigas; CpG DNA, unmethylated CG dinucleotides prevalent in bacterial DNA; d, drosophila; ECSIT, evolutionarily conserved signaling intermediate in Toll pathways; Es, Euprymna scolopes; I B, inhibitor of NF- B; IKK, I B kinase; Imd, immune deficiency; IRAK, interleukin-1-receptor-associated kinase; ird, immune response deficient; LPS, lipopolysaccharide; MEKK, mitogen-activated protein kinase/extracellular signal-regulated kinase kinase kinase; MyD88, myeloid differentiation primary response gene 88; N, N terminus; ND*, not determined that MAMPs activate this pathway in this superphylum; NEMO, NF- B essential modulator; NF- B, nuclear factor B; NIK, NF- B inducing kinase; NOD, nucleotide-binding oligomerization domain protein; o, oyster; PGN, peptidoglycan; PGRP, PGN recognition proteins; R gene, plant resistance genes; REL, proteins containing regions of homology to the rel oncogene, the definition of the NF- B family of genes; ssRNA, single-stranded RNA; TAB, TAK-binding protein; TAK, transforming growth factor ß-activated kinase 1; TLR, Toll-like receptor; Tollip, Toll interacting protein; TRAF, tumor necrosis factor receptor-associated factor.
|
Both biochemical and cellular evidence have suggested that, during the complex program of early host light organ development, V. fischeri MAMPs signal and activate host cells through pathways either analogous or homologous to the Toll/NF-
B signaling pathway. For example, mucus secretion, as a response to interactions with bacteria, is controlled by the Toll/NF-
B pathway in other systems (34, 40, 63). In addition, the activities of halide peroxidase (58, 65, 71) and nitric oxide synthase (16) have been implicated in the dynamics of the squid-vibrio system, and the transcription of the genes that encode these proteins is either directly or indirectly activated by NF-
B in other hosts (4, 60). At the cellular level, the recruitment of macrophage-like hemocytes into the light organ (39) and apoptosis of the cells of the superficial ciliated epithelium (24) are types of cellular behaviors modulated by activity of the NF-
B pathway in many host-pathogen interactions (13, 33, 36, 41, 55).
In the present study, we used juvenile light organ expressed sequence tag (EST) libraries to search for transcripts encoding proteins of the Toll/NF-
B pathway and rapid amplification of cDNA ends (RACE) PCR to obtain full-length sequences. Our data provide evidence that the pathway is present in E. scolopes, a finding that paves the way for the study of its regulation in beneficial animal-bacterial interactions using the squid-vibrio model.
|
|
|---|
Molecular analysis.
Candidate sequences were identified by analyses of a light organ EST database that had been constructed by using juvenile E. scolopes light organs (EST sequences available through The Wellcome Trust Sanger Institute and European Bioinformatics Institute [http://trace.ensembl.org/]). The sequences from the EST library were translated and compared to protein sequences already in GenBank by using the BlastX algorithm of the Blastcl3 netclient software (NCBI [ftp://ftp.ncbi.nih.gov/BLAST/blastcl3]) (3). For use in cloning and RACE (rapid amplification of cDNA ends) PCR techniques, total RNA was extracted from 100 juvenile light organs by using the MasterPure complete DNA and RNA purification kit (Epicentre, Madison, WI) according to the manufacturer's protocol. 5'- and 3'-RACE PCR was performed by using the SMART RACE cDNA amplification kit including the Advantage II PCR kit (BD Biosciences, Palo Alto, CA) and the FirstChoice RLM-RACE kit (Ambion, Austin, TX) according to the manufacturers' instructions. 5'- and 3'-RACE-ready cDNA was produced by using 1 µg of symbiotic juvenile light organ mRNA. RACE primers were constructed from the EST sequences identified to have similarity to a transcript of interest. The resulting amplification products were gel purified by using a Geneclean kit (Bio 101, Carlsbad, CA) or the QIAquick gel extraction kit (QIAGEN, Valencia, CA). Purified RACE products were ligated into pGEM-T Easy Vector (Promega, Madison, WI). Calcium-competent cells of Escherichia coli DH5
were transformed with the resultant plasmids. Recombinant colonies were screened for the insert of the predicted length, and plasmids were prepared for sequencing by using QIAprep spin miniprep kit columns (QIAGEN). PCR products were sequenced from the recombinant plasmid by using T7 and SP6 primers (Promega, Madison, WI).
Sequence analysis.
EST library clones that contained transcripts of interest were fully sequenced. In cases where the EST library clone did not contain the full-length transcript, 5' or 3' RACE PCR was used to obtain the complete transcript. For each RACE amplification product, a consensus contiguous sequence was generated by aligning all sequenced clones using MacVector and AssemblyLIGN or the DSgene software package (Accelrys, San Diego, CA), ensuring a three times independent PCR coverage per base. The full sequence of the transcript was determined by (i) the inability of 5' RACE PCR to amplify any further sequence, coupled with the presence of in-frame stop codons upstream of the prospective methionine start site at the 5' end, and (ii) the presence of a poly(A) tail following the prospective stop codon at the 3' end of the transcript. The resulting sequence was analyzed by BLAST searches of GenBank using the default settings (3). Potential open reading frames (ORFs) were predicted by alignment with sequences in GenBank and codon preference analysis (GCG Sequence Analysis; Genetics Computer Group, Inc.). Predicted protein sequences were generated by MacVector or DSgene. The predicted proteins were compared to protein sequences in GenBank by using conserved domain searches within BLAST (44) and SMART (61), secondary structure prediction (Jpred) (14), and alignment. Bootstrapped (1,000 iterations) neighbor-joining trees were constructed based on the Kimura two-parameter correction by using CLUSTAL X (69) with gapped positions included. Maximum-likelihood analyses were also performed on these data using the default parameters of the program proML with global rearrangement in the PHYLIP (version 3.63) suite of programs (22). Trees constructed by the neighbor-joining and maximum-likelihood methods were compared using the Kishino-Hasegawa-Templeton test in proML (38).
|
|
|---|
B pathway (Table 1 and Fig. 1), including molecules of the receptor complex, signaling molecules, and proteins with high similarity to members of the NF-
B family. The transcripts of genes that are regulated by the Toll/NF-
B pathway in other systems are also present (Table 1 and Fig. 1). With the exception of peptidoglycan recognition protein 3 (PGRP3), the 3' ESTs are similar to the 3' region of the most similar molecule in the database, adding to the confidence of their identity. The PGRP3 3' EST was not found to have sequence similarity with sequences in GenBank because it included only a small portion of the ORF. This EST was identified as a candidate PGRP by a BLASTn comparison of E. scolopes PGRP1, -2, and -4 against the entire squid 3' EST database. |
View this table: [in a new window] |
TABLE 1. E. scolopes 3' ESTs most similar to the molecules of the Toll/NF- B pathwaye
|
B pathway were investigated further, namely, a Toll-like receptor (TLR), the four PGRPs, interleukin-1 receptor-associated kinase 4 (IRAK4), TNF-
receptor associated factor 6 (TRAF6), inhibitor of
B kinase
(IKK
), inhibitor of
B
(I
B
), and REL (Table 2). In all cases, sequencing of the entire coding region improved the similarity to a sequence in the database. However, in the case of Toll, the PGRPs, IRAK4, TRAF6, I
B
, and REL, the best hit to the entire transcript was to the same protein but from a different organism than that of the 3' EST (Table 2). Trees constructed using the neighbor-joining and maximum likelihood with global rearrangement methods were not found to be statistically different. |
View this table: [in a new window] |
TABLE 2. Full-length transcripts of E. scolopes most similar to key molecules of the Toll/NF- B pathwayc
|
![]() View larger version (25K): [in a new window] |
FIG. 2. (A) Comparison of Toll-like receptor proteins from E. scolopes, fruit fly, and human. Gray regions indicate the signal peptide, black diagonal shading indicates the LRR region, black regions indicate a transmembrane domain, and black horizontal shading indicates the TIR domain. The number indicates the size of the derived amino acid sequence. Accession numbers are AY956809, NP_524757, and NP_003256, respectively. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between Toll-like protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. Symbols on the branches represent percentage bootstrap values: +, 90 to 100%; , 80 to 89%; , 70 to 79%; , 60 to 69%. Accession numbers clockwise from E. scolopes are as follows: AY956809, XP_393712, AAF86226, AAL37902, NP_524757, AAL37904, XP_393713, BAB85498, AAF86225, NP_649719, NP_523519, NP_733166, AAF86227, NP_003255, NP_112218, NP_003254, NP_006059, NP_003259, NP_942086, NP_001008664, NP_003256, NP_059138, AAF78035, NP_057694, NP_612564, and AAK37544. Abbreviations: Ag, Anopheles gambiae (mosquito); Am, Apis mellifera (honeybee); Bm, Bombyx mori (silk moth); Bt, Bos taurus (cow); Ce, Caenorhabditis elegans (nematode); Dm, Drosophila melanogaster (fruit fly); Es, Euprymna scolopes (Hawaiian bobtail squid); Hs, Homo sapiens (human); Rn, Rattus norvegicus (rat).
|
(i) EsPGRP1.
5' and 3' RACE revealed a 1,362-bp transcript that is predicted to encode a 207-aa protein (Table 2 and Fig. 3A). The putative protein sequence for EsPGRP1 contains a conserved PGRP domain (CDS E value of 1e43). Because no other protein domains, such as a signal peptide, were strongly suggested by in silico analysis of the protein structure, EsPGRP1 is predicted to have a cytosolic cellular localization.
![]() View larger version (73K): [in a new window] |
FIG. 3. (A) Comparison of PGRP proteins from E. scolopes, bay scallop, fruit fly, and human sources. Gray regions indicate the signal peptide, black regions indicate a transmembrane domain, and black diagonal shading indicate the PGRP domain. The number indicates the size of the derived amino acid sequence. Reading from left to right and down from the top, accession numbers are AY956811, AY956812, AY956813, AY956814, AAR92030, NP_572727, NP_731576, NP_005082, and AAL05629, respectively. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between PGRP protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. Symbols on the branches represent percentage bootstrap values: +, 90 to 100%; , 80 to 89%; , 70 to 79%. Values below 70% are not shown. Accession numbers clockwise from E. scolopes PGRP1 are as follows: AY956811, AY956812, AAR92030, AAL05629, XP_395941, NP_731576, AAG32064, NP_005082, NP_729468, AAG23736, NP_788467, AAK00295, NP_572727, BAA77210, CAD89136, CAD89193, NP_443123, EAA10910, AY956814, and AY956813. (C) Alignment of the PGRP domains of the EsPGRPs and those of other organisms. Residues critical to PGRP function are highlighted. Accession numbers from top to bottom are as follows: AY956811, AY956812, AY956813, AY956814, NP_729468, NP_005082, NP_731576, NP_572727, and AAL05629. Abbreviations: Ag, Anopheles gambiae (mosquito); Ai, Argopecten irradians (bay scallop); Am, Apis mellifera (honeybee); Bm, Bombyx mori (silk moth); Dm, Drosophila melanogaster (fruit fly); Es, Euprymna scolopes (Hawaiian bobtail squid); Hs, Homo sapiens (human).
|
(iii) EsPGRP3.
5' RACE PCR completed a 2,156-bp EsPGRP3 transcript encoding a 243-aa protein (Table 2 and Fig. 3A). At the nucleotide level EsPGRP3 does not share significant similarity to any of the other EsPGRPs, but the putative EsPGRP3 protein is strongly similar to that of EsPGRP1 and EsPGRP2 (60% identity) with a somewhat weaker similarity to EsPGRP4 (48% identity). EsPGRP3 protein is predicted to have a signal peptide and a PGRP domain (CDS E value of 6.4e68). It is also strongly predicted to have a glycosylphosphatidylinositol (GPI)-anchor, which would bind the PGRP to the outer surface of the cell membrane, with a possible GPI-anchor cleavage site at position 216. Therefore, we predict that EsPGRP3 is an extracellular extrinsic membrane protein that may be released into the extracellular environment through cleavage of the GPI-anchor.
(iv) EsPGRP4.
The 5' and 3' RACE PCR product sequence indicated that the EsPGRP4 transcript is 1,527 bp encoding a 270-aa protein (Table 2 and Fig. 3A). At the nucleotide level, the coding sequence of EsPGRP4 does not share significant similarity to either EsPGRP1 or EsPGRP2 (using the BLASTn algorithm). The EsPGRP4 protein sequence is about 45 and 52% identical to EsPGRP1 and EsPGRP2, respectively. EsPGRP4 contains a putative PGRP domain (CDS E value of 1e42). The prediction of two N-terminal transmembrane domains suggests that EsPGRP4 is an integral membrane protein.
When the E. scolopes PGRPs (EsPGRPs) are compared to PGRP molecules from other organisms (Fig. 3B), the E. scolopes molecules group together with strong bootstrap support (>70%). Apart from this cluster and the cluster containing the PGRP-LD and PGRP-LAa from D. melanogaster, there is no other strong bootstrap support for the PGRPs of a single organism grouping together. EsPGRPs also did not group with a PGRP from another mollusk, Argopecten irradians (the bay scallop). These data suggest that the PGRPs identified from the E. scolopes light organ were derived from a single ancestral EsPGRP. Such findings do not preclude the possibility of other EsPGRPs that would cluster outside of those already found, and are derived from a different PGRP.
Some PGRPs have T7 lysozyme-like catalytic activity, including PGRP-LB from D. melanogaster (37, 45), which coordinates a zinc ion required for N-acetylmuramic acid L->alanine amidase (NAMLAA) activity with a conserved tetrad of two histidines, a critical cysteine, and tyrosine residue (H42, Y78, H152, and C160, PGRP-LB numbering). EsPGRP1, EsPGRP2, and EsPGRP3 contain all four zinc coordinating residues in the PGRP-binding domain, suggesting zinc-dependent NAMLAA activity in these three proteins (Fig. 3C). At the homologous position of the critical cysteine 160, EsPGRP4 contains a serine substitution similar to D. melanogaster PGRP-SA, which does not coordinate a zinc ion in the binding pocket (56). It is likely that EsPGRP4 does not act as a catalytic PGRP and is limited to PGN recognition and/or signal transduction functions. All four EsPGRPs contain a pair of cysteines C50 and C56 (using PGRP-LB numbering) at synonymous positions with those essential for the structure of the PGN binding pocket and a highly conserved threonine, T158 (using PGRP-LB numbering) that is required in the binding of PGN by PGRP-LB from D. melanogaster (37).
IRAK4.
5' and 3' RACE PCR of juvenile light organ RNA resulted in a 1,791-bp transcript that contains an ORF encoding a 537-aa protein (Table 2 and Fig. 4). The translated transcript exhibited a C-terminal serine/threonine protein kinase domain (CDS E value of 5e34) and at least four N-terminal
-helices of a possible death domain (31), characteristic of IRAK proteins (Fig. 4A and Fig. S2 in the supplemental material). When the E. scolopes IRAK4 is compared to IRAK4-like molecules from other organisms, the E. scolopes molecule clusters with the vertebrate IRAK4 protein sequences with strong bootstrap support (Fig. 4B). However, the branch lengths within this cluster indicate that the E. scolopes sequence is separate from the vertebrate IRAK4 cluster.
![]() View larger version (14K): [in a new window] |
FIG. 4. (A) Comparison of IRAK4-like proteins from E. scolopes, rat, fruit fly, and tomato. Black regions indicate the serine/threonine protein kinase domain, and gray regions indicate the alpha helices of a possible death domain. Number indicates size of the derived amino acid sequence. The accession numbers are AY956810, XP_217026, NP_476971, and T07416, respectively. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between IRAK4 protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. The numbers on the branches are percent bootstrap values. Accession numbers clockwise from E. scolopes are as follows: AY956810, NP_057207, NP_084202, XP_217026, NP_502587, NP_476971, XP_311931, AAF76311, T07416, and AAK11568. Abbreviations: Ag, Anopheles gambiae (mosquito); Ce, Caenorhabditis elegans (nematode); Dm, Drosophila melanogaster (fuit fly); Es, Euprymna scolopes (Hawaiian bobtail squid); Hs, Homo sapiens (human); Le, Lycopersicon esculentum (tomato); Lh, Lycopersicon hirsutum (tomato); Mm, Mus musculus (mouse); Rn, Rattus norvegicus (rat).
|
B activation (7, 12, 15). Comparison of the E. scolopes sequence with TRAF sequences from other organisms supports its classification as TRAF6 (Fig. 5B). E. scolopes TRAF6 clusters with the vertebrate TRAF6 sequences with strong bootstrap support. E. scolopes TRAF6 is 33% similar to human TRAF6, whereas only sharing similarities of 26, 28, 26, and 29% to TRAF2, TRAF3, TRAF4, and TRAF5 of humans, respectively. TRAF-like proteins from D. melanogaster and A. gambiae branch one node away from TRAF6 of E. scolopes and vertebrates, and this grouping is also well supported. The branch lengths indicate that E. scolopes, vertebrates, and insects are separate from each other within the TRAF6-like sequence cluster.
![]() View larger version (21K): [in a new window] |
FIG. 5. (A) Comparison of TRAF6 proteins from E. scolopes, zebrafish, human, and fruit fly with characteristic domains indicated. Black hatched shading indicates the RING finger domain, gray regions indicate the TRAF-type zinc finger domains, and black regions indicate the TRAF homology domain. The number indicates size of the derived amino acid sequence. The accession numbers are AY956816, NP_956115, NP_004611, and AAD47895, respectively. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between TRAF protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. Numbers on the branches are percentage bootstrap values. Accession numbers clockwise from E. scolopes are as follows: AY956816, NP_956115, NP_004611, P70196, XP_230377, XP_311879, AAD34345, AAD47895, AAD34346, NP_991325, NP_004286, NP_033449, XP_220640, CAE88928, CAC82653, CAD69021, NP_066961, XP_231032, NP_989550, NP_004610, NP_035763, CAD57164, NP_035762, and Q13114. Abbreviations: Ag, Anopheles gambiae (mosquito); Ca, Carassius auratus (goldfish); Dm, Drosophila melanogaster (fruit fly); Dr, Danio rerio (zebrafish); Es, Euprymna scolopes (Hawaiian bobtail squid); Gg, Gallus gallus (chicken); He, Hydractinia echinata (hydrozoa); Hs, Homo sapiens (human); Mm, Mus musculus (mouse); Om, Oncorhynchus mykiss (rainbow trout); Rn, Rattus norvegicus (rat).
|
.
proteins requires secondary structure analysis, as opposed to conserved domain searches, because IKK
proteins do not contain unique IKK
domains. IKK
proteins are predominantly helical with large stretches of coiled-coil structure, including a leucine zipper motif and a zinc finger motif at the C terminus (35, 57). These motifs are likely to be involved in protein-protein interactions. Secondary structure analysis of the derived amino acid sequence (Fig. 6A and Fig. S4 in the supplemental material) indicates that the E. scolopes protein is a canonical IKK
(35, 57). Similarly, the E. scolopes sequence clusters with IKK
proteins of vertebrates with strong bootstrap support, although the branch lengths indicate the derived amino acid sequence of E. scolopes is separate from the vertebrate sequences (Fig. 6B). IKK
from D. melanogaster branches one node away and is separated by a long branch from IKK
of E. scolopes and vertebrates. The insect, vertebrate, and mollusk proteins form a well-supported IKK
grouping that is separate from IKK
and IKKß sequences.
![]() View larger version (42K): [in a new window] |
FIG. 6. (A) The characteristic domains of IKK exhibited by the derived amino acid sequence from E. scolopes. The underlined region represents a helical structure containing coiled coil regions (gray shading) and a leucine zipper motif with specific leucines (and hydrophobic isoleucine) marked with an asterisk. The C-terminal zinc finger motif is shaded in black with characteristic amino acids labeled with a "#" symbol. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between IKK protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. Numbers on the branches are percentage bootstrap values. *, The original description does not indicate assignment to , ß, or classes of IKK. Accession numbers clockwise from E. scolopes are as follows: AY956817, NP_499831, AAG02485, AAC05683, NP_001269, AAH18243, NP_001547, NP_034676, AAG02486, NP_989567, AAH63986, NP_776779, NP_003630, and AAM54725. Abbreviations: Bt, Bos taurus (cow); Ce, Caernorhabditis elegans (nematode); Cg, Crassostrea gigas (Pacific oyster); Dm, Drosophila melanogaster (fruit fly); Dr, Danio rerio (zebrafish); Es, Euprymna scolopes (Hawaiian bobtail squid); Gg, Gallus gallus (chicken); Hs, Homo sapiens (human); Mm, Mus musculus (mouse).
|
B.
B family members, the derived amino acid sequence exhibited six ankyrin repeat domains (CDS E value of 9e23), an N-terminal regulatory region, and a C-terminal region rich in Pro, Glu/Asp, Ser, and Thr residues (PEST-rich region) (Fig. 7A and Fig. S5 in the supplemental material). Proteins of the I
B family contain either six or seven ankyrin repeats that mediate binding to the REL homology domain and mask the nuclear localization signal of REL family proteins (35). The N-terminal regulatory region is required for stimulus-induced degradation of the protein, the key step in NF-
B activation. Only the "professional" I
Bs (74), namely, I
B
, I
Bß, I
B
, and Cactus, contain this short peptide stretch centered around two serine residues that can be inducibly phosphorylated (35). Mutations at either serine residue of this region results in the stabilization of the protein (6, 50, 74). The presence of this N-terminal regulatory region, coupled with the inability of 5' RACE PCR to amplify any further sequence and in-frame stop codons upstream of the prospective start site, provides evidence that the sequence is complete and is not the C terminus of a p105-like or RELISH-like transcript that also has I
B-like characteristics (27, 41, 64). However, the sequence of the regulatory region in the E. scolopes protein has a conserved substitution in the third position, i.e., instead of the DSGXXS characteristic of the vertebrate I
Bs and the drosophila homolog Cactus, the sequence of this region of the E. scolopes I
B is DSAXXS. The ubiquitin conjugation site, a lysine residue located 9 to 12 aa N-terminal to the recognition site of professional I
B is, however, shared by the E. scolopes sequence (74). The C-terminal contains a PEST-rich region. PEST sequences are often found in proteins with a rapid turnover rate and may be necessary for the degradation of I
B. For example, in mouse cell lines that were transfected with human I
B
missing the final 41 aa, including the PEST sequence, the mutant I
B
was phosphorylated but not degraded (11). The E. scolopes sequence clustered with I
B
-like proteins (Fig. 7B). Within this I
B
cluster, the putative protein sequence from E. scolopes did not group with either D. melanogaster Cactus or vertebrate I
B
proteins with strong bootstrap support, and branch lengths indicated that the vertebrates, D. melanogaster, and E. scolopes sequences were distinct from each other (Fig. 7B).
![]() View larger version (19K): [in a new window] |
FIG. 7. (A) Comparison of I B family molecules from E. scolopes, rat, and fruit fly. Black diagonal shading represents the N-terminal regulatory region, with the serines of the phosphorylation site and the lysines of the possible ubiquitin conjugation sites underlined and highlighted in boldface. The ankyrin repeat domains are shaded in gray, and the PEST-rich region is shaded in black. The number indicates the size of the derived amino-acid sequence. The accession numbers are AY956818, XP_343066, and NP_476942, respectively. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between I B protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. The numbers on the branches are percent bootstrap values. The accession numbers clockwise from E. scolopes are as follows: AY956818, CAC85086, Q91974, NP_035037, XP_343066, Q08353, CAD12786, O00221, NP_032716, AAB58262, NP_002494, NP_035038, and NP_476942. Abbreviations: At, Arabidopsis thaliana (thale cress); Bt, Bos taurus (cow); Dm, Drosophila melanogaster (fruit fly); Es, Euprymna scolopes (Hawaiian bobtail squid); Gg, Gallus gallus (chicken); Hs, Homo sapiens (human); Mm, Mus musculus (mouse); Om, Oncorhynchus mykiss (rainbow trout); Rn, Rattus norvegicus (rat); Ss, Sus scrofa (pig).
|
B (35). The transactivation domain, C-terminal of the RHD, is less well-conserved and is thought to be involved in the activation of transcription (41). No regions that are critical for transcriptional activation have been identified within the transactivation domain, apart from the presence of the transactivation domain itself (26, 27, 41, 47). The REL family of proteins can be separated into two distinct classes: class I contains the precursor REL proteins with C-terminal domains that have I
B-like characteristics, such as vertebrate and insect p100, p105, and Relish; class II contains the NF-
B subunit REL proteins, including p65/RelA, RelB, c-REL, v-Rel, and Dorsal (27). Compared to the other class II REL proteins, the E. scolopes REL sequence is relatively short (Fig. 8A). The difference appears to be in the length of the transactivation domain, since the size of the RHD of both proteins is similar. Since there is a 3' stop codon and 3' RACE PCR failed to amplify any sequence further 3' of the poly(A) tail and, since the transactivation domain is known to be less conserved, we conclude that the sequence is complete. Compared to class II REL proteins from other organisms, the E. scolopes sequence clusters with the Crassostrea gigas REL sequences with strong bootstrap support (Fig. 8B). Since the insects and the vertebrate class II REL proteins, respectively, also form groupings that are strongly supported and the branch lengths separating the sequences are long, it is concluded that the mollusk sequences are distinct from the vertebrate and insect class II REL proteins.
![]() View larger version (26K): [in a new window] |
FIG. 8. (A) Comparison of REL family proteins from E. scolopes, oyster, human, and fruit fly. The RHD is indicated by gray shading and contains the following: the conserved residues of the DNA-binding motif (black horizontal shading) underlined and highlighted in boldface; the immunoglobulin-like plexin domain, found in all REL family molecules, indicated by black diagonal shading; and the nuclear localization signal indicated by solid black shading, with conserved residues underlined and highlighted in boldface. The unshaded region that is C-terminal to the RHD represents the transactivation domain. The number indicates size of the derived amino acid sequence. Accession numbers are AY956819, AAK72691, NP_068810, and NP_724052, respectively. (B) Unrooted phylogenetic tree produced by the neighbor-joining method, using Kimura two-parameter distances between REL protein sequences of the organisms named on the tree. Gapped positions were included in the analysis. Numbers on the branches are percentage bootstrap values. Accession numbers clockwise from E. scolopes are as follows: AY956819, AAK72690, AAK72691, Q04865, NP_068810, NP_033071, XP_238994, XP_223688, P16236, NP_033072, Q01201, BAB47172, BAB47173, XP_310177, AAP23055, AAG22858, NP_724052, and NP_724054. Abbreviations: Ag, Anopheles gambiae (mosquito); Am, Apis mellifera (honeybee); Cg, Crassostrea gigas (Pacific oyster); Dm, Drosophila melanogaster (fruit fly); Es, Euprymna scolopes (Hawaiian bobtail squid); Hr, Halocynthia roretzi (ascidian); Hs, Homo sapiens (human); Mm, Mus musculus (mouse); Rn, Rattus norvegicus (rat); Tc, Tribolium castaneum (beetle); XI, Xenopus laevis (African clawed frog).
|
|
|
|---|
B pathway in E. scolopes, including receptors, signal transducers, and targets of the pathway (Fig. 9). Although this pathway appears to be present, its characteristics suggest certain unique features. This finding is not surprising, since E. scolopes is phylogenetically distant from the other animals for which this pathway has been described in detail.
![]() View larger version (33K): [in a new window] |
FIG. 9. Putative signaling pathway of the evolutionarily conserved Toll/NF- B pathway utilizing the E. scolopes transcripts identified in the present study. Black shading indicates identification of full-length transcripts, and gray shading indicates identification of partial length transcripts (see Tables 1 and 2). Dashed lines indicate a pathway in which not all of the signaling intermediates have been found in E. scolopes.
|
Within the pathway, several elements critical for MAMP-induced triggering of the NF-
B pathway were identified, including IRAK4, TRAF6, IKK
, I
B-like, and REL-like transcripts (17, 25, 35, 43, 57, 64, 67, 73). In another mollusk, Crassostrea gigas (the Pacific oyster), full-length transcripts of an IKK (20), which appears to be an IKKß, and a REL-like protein (48) were identified, lending further support for the presence of the pathway in the mollusks. An element not identified in the E. scolopes EST database was a MyD88-like transcript, which in other systems encodes a critical protein that links the Toll-like receptor to the IRAK protein. However, a partial sequence of a molluscan MyD88-like transcript has been identified in C. gigas (28), suggesting that E. scolopes may have this gene, although it may not be expressed at high levels in the light organ. Attempts to identify this transcript in other E. scolopes tissues by reverse transcription-PCR amplification were not successful.
The transcripts described in the present study can be used to propose a putative Toll/NF-
B signaling pathway in E. scolopes (Fig. 9). Taken together, the data for E. scolopes and the distantly related C. gigas provide strong evidence that this pathway is likely to occur throughout the Phylum Mollusca, as well as across the entire superphylum containing the mollusks, i.e., the Lophotrochozoa. The phylogenetic trees for the E. scolopes proteins do not clarify whether the E. scolopes sequences group closer to the vertebrates or the insects, although most of the closest BLAST hits are to vertebrate sequences. The relationship among the superphyla of the animals is controversial, with most phylogenies aligning the Lophotrochozoa with the Ecdysozoa, rather than the Deuterostomia (59). Further analysis of the E. scolopes EST library and those of other Lophotrochozoa will be necessary to address this issue rigorously.
The complex patterns of bacteria-induced development in the squid light organ, along with the different functions of MAMPs in the process, suggest that the NF-
B pathway may be modulated differently over time and in location through the early development of the light organ. The high levels of nitric oxide synthase (NOS) (16) and the mucus shedding (52, 54) in response to environmental PGN would be consistent with activity of the pathway in the superficial field of ciliated cells during the initial phases of colonization. However, nitric oxide/NOS levels are attenuated (16), mucus shedding ceases (52), and apoptosis (24) is induced in this field of cells by the subsequent interactions of the host and symbiont within the crypt space, which is several cell layers away. These observations suggest that remote signaling by the symbionts may turn down the activity of the pathway in the superficial field of ciliated cells. In contrast, the MAMPs responsible for development, i.e., derivatives of PGN and LPS, are presented directly to the crypt cells. These receptor-ligand interactions may increase the activity of the pathway in the cells with which the bacteria directly interact. In support of this possibility is the observation that mucus shedding by the crypt cells results from colonization of the crypts with symbionts (52).
The data presented here demonstrate the presence of the NF-
B pathway in the E. scolopes light organ, and the patterns of development and the involvement of MAMPs provide circumstantial evidence for its activity during the onset of the symbiosis. A full understanding of the connection among these findings awaits an in-depth analysis of the activity of the pathway in the various tissues of the light organ.
This research was supported by a grant from the W. M. Keck Foundation (to E. P. Greenberg and M.M.-N.), by NSF grant IBN 0451102 (to M.M.-N. and E. G. Ruby), by NIH grant AI 50611 (to M.M.-N.), by NIH grant RR12294 (to E. G. Ruby and M.M.-N.), and by a Career Development Award from Research to Prevent Blindness (to T.E.S).
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
Present address: CEA Cadarache DEVM/LBC-Bat 156, 13108 Saint Paul-Lez-Durance Cedex, France. ![]()
|
|
|---|
B and I
B proteins: new discoveries and insights. Annu. Rev. Immunol. 14:649-683.[CrossRef][Medline]
B
proteolysis by site-specific, signal-induced phosphorylation. Science 267:1485-1488.
B-dependent inhibition of apoptosis is essential for host cell survival during Rickettsia rickettsii infection. Proc. Natl. Acad. Sci. USA 95:4646-4651.
B activation but is dispensable for c-Jun N-terminal kinase signaling. J. Biol. Chem. 273:24775-24780.
and ß subunits of I
B kinase (IKK) mediate TRAF2-dependent IKK recruitment to tumor necrosis factor (TNF) receptor 1 in response to TNF. Mol. Cell. Biol. 21:3986-3994.
B and Rel proteins: evolutionarily conserved mediators of immune responses. Annu. Rev. Immunol. 16:225-260.[CrossRef][Medline]
B, KBF1, dorsal, and related matters. Cell 62:841-843.[CrossRef][Medline]
B signal transduction pathway: introduction. Oncogene 18:6842-6844.[CrossRef][Medline]
B/NF-
B system: a key determinant of mucosal inflammation and protection. Am. J. Physiol. Cell Physiol. 278:C451-C462.
-Smad signaling pathway cooperates with NF-
B to mediate nontypeable Haemophilus influenzae-induced MUC2 mucin transcription. J. Biol. Chem. 277:45547-45557.
B activity. Annu. Rev. Immunol. 18:621-663.[CrossRef][Medline]
B at the crossroads of life and death. Nat. Immunol. 3:221-227.[CrossRef][Medline]
B via a Src-dependent Ras-MAPK-pp90rsk pathway is required for Pseudomonas aeruginosa-induced mucin overproduction in epithelial cells. Proc. Natl. Acad. Sci. USA 95:5718-5723.
B regulation in the immune system. Nat. Rev. Immunol. 2:725-734.[CrossRef][Medline]
B/I
B story. Adv. Cancer Res. 66:255-292.[Medline]
B homolog characterized in a mollusk, the Pacific oyster Crassostrea gigas. FEBS Lett. 561:75-82.[CrossRef][Medline]
B functions as a tumour promoter in inflammation-associated cancer. Nature 431:461-466.[CrossRef][Medline]
is an essential regulatory subunit of the I
B kinase complex. Nature 395:297-300.[CrossRef][Medline]
B signaling pathways in mammalian and insect innate immunity. Genes Dev. 15:2321-2342.
B kinase complex essential for NF-
B activation. Cell 93:1231-1240.[CrossRef][Medline]
B cellular function via specific targeting of the I
B-ubiquitin ligase. EMBO J. 16:6486-6494.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»