Both IRAK3 and IRAK1 Activate the MyD88-TRAF6 Pathway in Zebrafish.

IL-1R-associated kinases (IRAKs) are signal transducers of the TLR/IL-1R-MyD88-TRAF6 pathways. Vertebrates possess two IRAK lineages, IRAK1/2/3 and IRAK4. In mammals, IRAK4/IRAK1 and IRAK4/IRAK2 are pathway enhancers, whereas IRAK3 is a repressor. However, in bony fish, IRAK2 is absent, and it remains elusive how fish IRAK1/3/4 functionally differ from their mammalian counterparts. In this study, we explored this using the zebrafish model. First, we showed that in human 293T cells, zebrafish IRAK1 and IRAK4 were components of the Myddosome (MyD88-IRAK4-IRAK1) complex, with IRAK1 serving as a potent pathway enhancer. Then, we discovered two zebrafish IRAK3 variants: one (IRAK3a) contains an N-terminal Death domain, a middle pseudokinase domain, and a C-terminal TRAF6-binding domain, whereas the other (IRAK3b) lost both the kinase and TRAF6-binding domains. This truncation of IRAK3 variants could be a conserved phenomenon in fish, because it is also observed in trout and grass carp. We proceeded to show that zebrafish IRAK3a acts as a pathway enhancer by binding with MyD88 and TRAF6, but its activity is milder than IRAK1, possibly because it has no kinase activity. Zebrafish IRAK3b, however, plays a sheer negative role, apparently because of its lack of kinase and TRAF6-binding domains. Moreover, zebrafish IRAK3a/3b inhibit the activity of IRAK1/4, not by interacting with IRAK1/4 but possibly by competing for MyD88 and TRAF6. Finally, we have verified the essential activities of zebrafish IRAK1/3a/3b/4 in zebrafish cells and embryos. In summary, to our knowledge, our findings provide new insights into the molecular functions of fish IRAKs and the evolution of the IRAK functional modes in vertebrates.

The IL-17 pathway intertwines with neurotrophin and TLR/IL-1R pathways since its domain shuffling origin.

The IL-17 pathway displays remarkably diverse functional modes between different subphyla, classes, and even orders, yet its driving factors remains elusive. Here, we demonstrate that the IL-17 pathway originated through domain shuffling between a Toll-like receptor (TLR)/IL-1R pathway and a neurotrophin-RTK (receptor-tyrosine-kinase) pathway (a Trunk-Torso pathway). Unlike other new pathways that evolve independently, the IL-17 pathway remains intertwined with its donor pathways throughout later evolution. This intertwining not only influenced the gains and losses of domains and components in the pathway but also drove the diversification of the pathway's functional modes among animal lineages. For instance, we reveal that the crustacean female sex hormone, a neurotrophin inducing sex differentiation, could interact with IL-17Rs and thus be classified as true IL-17s. Additionally, the insect prothoracicotropic hormone, a neurotrophin initiating ecdysis in Drosophila by binding to Torso, could bind to IL-17Rs in other insects. Furthermore, IL-17R and TLR/IL-1R pathways maintain crosstalk in amphioxus and zebrafish. Moreover, the loss of the Death domain in the pathway adaptor connection to IκB kinase and stress-activated protein kinase (CIKSs) dramatically reduced their abilities to activate nuclear factor-kappaB (NF-κB) and activator protein 1 (AP-1) in amphioxus and zebrafish. Reinstating this Death domain not only enhanced NF-κB/AP-1 activation but also strengthened anti-bacterial immunity in zebrafish larvae. This could explain why the mammalian IL-17 pathway, whose CIKS also lacks Death, is considered a weak signaling activator, relying on synergies with other pathways. Our findings provide insights into the functional diversity of the IL-17 pathway and unveil evolutionary principles that could govern the pathway and be used to redesign and manipulate it.

SUMO and PIAS repress NF-κB activation in a basal chordate.

Small ubiquitin-like modifier (SUMO) regulates various biological processes, including the MyD88/TICAMs-IRAKs-TRAF6-NF-κB pathway, one of the core immune pathways. However, its functions are inconsistent between invertebrates and vertebrates and have rarely been investigated in lower chordates, including amphioxus and fishes. Here, we investigated the SUMOylation gene system in the amphioxus, a living basal chordate. We found that amphioxus has a SUMOylation system that has a complete set of genes and preserves several ancestral traits. We proceeded to study their molecular functions using the mammal cell lines. Both amphioxus SUMO1 and SUMO2 were shown to be able to attach to NF-κB Rel and to inhibit NF-κB activation by 50-75% in a dose-dependent fashion. The inhibition by SUMO2 could be further enhanced by the addition of the SUMO E2 ligase UBC9. In comparison, while human SUMO2 inhibited RelA, human SUMO1 slightly activated RelA. We also showed that, similar to human PIAS1-4, amphioxus PIAS could serve as a SUMO E3 ligase and promote its self-SUMOylation. This suggests that amphioxus PIAS is functionally compatible in human cells. Moreover, we showed that amphioxus PIAS is not only able to inhibit NF-κB activation induced by MyD88, TICAM-like, TRAF6 and IRAK4 but also able to suppress NF-κB Rel completely in the presence of SUMO1/2 in a dose-insensitive manner. This suggests that PIAS could effectively block Rel by promoting Rel SUMOylation. In comparison, in humans, only PIAS3, but not PIAS1/2/4, has been reported to promote NF-κB SUMOylation. Taken together, the findings from amphioxus, together with those from mammals and other species, not only offer insights into the functional volatility of the animal SUMO system, but also shed light on its evolutionary transitions from amphioxus to fish, and ultimately to humans.

Differential roles of the fish chitinous membrane in gut barrier immunity and digestive compartments.

The chitin-based peritrophic matrix (PM) is a structure critical for both gut immunity and digestion in invertebrates. PM was traditionally considered lost in all vertebrates, but a PM-like chitinous membrane (CM) has recently been discovered in fishes, which may increase the knowledge on vertebrate gut physiology and structural evolution. Here, we show that in zebrafish, the CM affects ingestion behavior, microbial homeostasis, epithelial renewal, digestion, growth, and longevity. Young mutant fish without CM appear healthy and are able to complete their life cycle normally, but with increasing age they develop gut inflammation, resulting in gut atrophy. Unlike mammals, zebrafish have no visible gel-forming mucin layers to protect their gut epithelia, but at least in young fish, the CM is not a prerequisite for the antibacterial gut immunity. These findings provide new insights into the role of the CM in fish prosperity and its eventual loss in tetrapods. These findings may also help to improve fish health and conservation, as well as to advance the understanding of vertebrate gut physiology and human intestinal diseases.

Two Amphioxus ApeC-Containing Proteins Bind to Microbes and Inhibit the TRAF6 Pathway.

The apextrin C-terminal (ApeC) domain is a class of newly discovered protein domains with an origin dating back to prokaryotes. ApeC-containing proteins (ACPs) have been found in various marine and aquatic invertebrates, but their functions and the underlying mechanisms are largely unknown. Early studies suggested that amphioxus ACP1 and ACP2 bind to bacterial cell walls and have a role in immunity. Here we identified another two amphioxus ACPs (ACP3 and ACP5), which belong to the same phylogenetic clade with ACP1/2, but show distinct expression patterns and sequence divergence (40-50% sequence identities). Both ACP3 and ACP5 were mainly expressed in the intestine and hepatic cecum, and could be up-regulated after bacterial challenge. Both prokaryotic-expressed recombinant ACP3 and ACP5 could bind with several species of bacteria and yeasts, showing agglutinating activity but no microbicidal activity. ELISA assays suggested that their ApeC domains could interact with peptidoglycan (PGN), but not with lipoteichoic acid (LTA), lipopolysaccharides (LPS) and zymosan A. Furthermore, they can only bind to Lys-type PGN from Staphylococcus aureus, but not to DAP-type PGN from Bacillus subtilis and not to moieties of PGN such as MDPs, NAMs and NAGs. This recognition spectrum is different from that of ACP1/2. We also found that when expressed in mammalian cells, ACP3 could interact with TRAF6 via a conserved non-ApeC region, which inhibited the ubiquitination of TRAF6 and hence suppressed downstream NF-κB activation. This work helped define a novel subfamily of ACPs, which have conserved structures, and have related yet diversified molecular functions. Its members have dual roles, with ApeC as a lectin and a conserved unknown region as a signal transduction regulator. These findings expand our understanding of the ACP functions and may guide future research on the role of ACPs in different animal clades.

Chordate PIAS proteins act as conserved repressors of the TRAF6 self-polyubiquitination.

In mammals, PIAS proteins are important SUMO E3 ligases and act as versatile regulators of over sixty different proteins, including components from the NF-κB pathways. But the PIAS functions are not well-understood due to complicated molecular mechanisms and multiple gene paralogs with overlapping roles, which is especially true in lower vertebrates where dedicated studies are scarce. As a basal chordate with a single PIAS gene, amphioxus is a convenient model to study PIAS from the evolutionary perspective. TRAF6 is a critical adaptor of the NF-κB pathways but it is not known whether TRAF6 is regulated by PIAS. Here we discover that in mammalian cells, amphioxus PIAS inhibited NF-κB activation by co-localizing and binding with TRAF6. The interaction relied on the N-terminal SAP and PINIT domains of PIAS. TRAF6 is an E3 ubiquitin ligase, which initiates downstream NF-κB signaling by promoting its self-ubiquitination. Both amphioxus SUMO1 and Ubc9 (SUMO E2 ligase) could suppress TRAF6 self-ubiquitination and NF-κB activation, suggesting that the SUMOylation activity competed away the ubiquitination activity of TRAF6. However, we show that the wild-type PIAS and the mutant PIAS without SUMO E3 ligase activity both could inhibit TRAF6-mediated NF-κB activation by reducing TRAF6 self-ubiquitination. This implies that SUMO ligase activity is not the only mechanism for PIAS to negatively regulate TRAF6. Finally, we tested the interactions between human PIAS1-4 and TRAF6. It reveals that human PIAS1, 3 and 4, but not 2, were able to repress NF-κB activation by reducing TRAF6 self-ubiquitination. Taken together, our study discovers a conserved regulatory interaction between chordate PIAS and TRAF6. It therefore sheds light on the complicated role of PIAS in immune regulation, and may help to understand the PIAS functions in other lower chordate taxa, such as jawless and jawed fishes.

The family of amphioxus chitin synthases offers insight into the evolution of chitin formation in chordates.

Chitin is a very important and widely-used biopolymer in fungi and lower metazoans, but mysteriously disappears in mammals. Recent studies reveal that at least lower vertebrates have chitin synthases (CS) and use them to synthesize endogenous chitin. Amphioxus, a basal chordate, therefore becomes critical to understand the evolution of CS, as it occupies the transitional position from invertebrates to vertebrates, and is considered as a good proxy to the chordate ancestor. Here, by exploiting multiple genome assemblies, high-depth RNA-seq data and synteny relations, we identify 11-12 CS genes for each amphioxus species. It represents the largest CS gene pool ever found in eukaryotes so far. As comparison, most metazoans have one or two CSs. Amphioxus is the only chordate that has both the very ancient type-I CS family and the more broadly distributed type-II CS family. Specifically, amphioxus has only one type-II CS but 10-11 type-I CSs, which means that amphioxus is the only metazoan with a greatly expanded type-I CS family. Further analysis suggests that the chordate ancestor have at least one type-II CS and an expanded of type-I CS family. We hypothesize that: these ancient CSs are mostly retained in amphioxus; but the whole type-I CS family was lost in urochordates and vertebrates; the type-II CS was later duplicated into two lineages in vertebrates and followed by stochastic losses, till all type-II CSs were eventually lost in birds and mammals. Finally, our expression profiling and preliminary gene knockout analysis suggest that amphioxus CSs could have highly diverse but mildly overlapping functions in various tissues and organs. Taken together, these findings not only provide insights into the evolution of chordate CSs, lay a foundation for further functional study of the chordate CSs. After all, it is mysterious that our chordate ancestor needed so many isoenzymes for chitin formation.

LanceletDB: an integrated genome database for lancelet, comparing domain types and combination in orthologues among lancelet and other species.

Lancelet (amphioxus) represents the most basally divergent extant chordate (cephalochordates) that diverged from the other two chordate lineages (urochordates and vertebrates) more than half a billion years ago. As it occupies a key position in evolution, it is considered as one of the best proxies for understanding the chordate ancestral state. Thus, the construction of a database with multiple lancelet genomes and gene annotation data, including protein domains, is urgently needed to investigate the loss and gain of domains in orthologues among species, especially ancient domain types (non-vertebrate-specific domains) and novel domain combination, which is helpful for providing new insight into the chordate ancestral state and vertebrate evolution. Here, we present an integrated genome database for lancelet, LanceletDB, which provides reference haploid genome sequence and annotation data for lancelet (Branchiostoma belcheri), including gene models and annotation, protein domain types, gene expression pattern in embryogenesis, different expression sequence tag sets and alternative polyadenylation (APA) sites profiled by the sequencing APA sites method. Especially, LanceletDB allows comparison of domain types and combination in orthologues among type species so as to decode the ancient domain types and novel domain combination during evolution. We also integrated the released diploid lancelet genome annotation data (Branchiostoma floridae) to expand LanceletDB and extend its usefulness. These data are available through the search and analysis page, basic local alignment search tool page and genome browser to provide an integrated display.

The conserved ancient role of chordate PIAS as a multilevel repressor of the NF-κB pathway.

In vertebrates, PIAS genes encode versatile cellular regulators, with functions extremely complex and redundant. Here we try to understand their functions from an evolutionary perspective. we evaluate the sequences, expression and molecular functions of amphioxus PIAS genes and compare them with their vertebrate counterparts. Phylogenetic analysis suggests a single PIAS gene in ancestral chordates, which has been duplicated into four families (PIAS1-4) in vertebrates by 2R-WGD but remained single in a basal chordate (amphioxus). Amphioxus PIAS encodes two variants with and without a Serine/Threonine-rich tail, which are retained in human PIAS1-3 but lost in PIAS4. We show that amphioxus PIAS binds C-terminus of NF-κB Rel and blocks the DNA binding activity. In humans, such function is retained in PIAS1, altered in PIAS4, and lost in PIAS2-3. Instead, PIAS3 has evolved new ability to inhibit Rel by binding RHD and promoting SUMOylation. We show that amphioxus PIAS also inhibits NF-κB by binding with upstream signalling adaptor TICAM-like and MyD88. Finally, we verify that human PIAS1, 3 and 4, but not 2, were capable of these newly-discovered functions. Our study offers insight into the sub- and neo-functionalization of PIAS genes and suggests a conserved ancient role for chordate PIAS in NF-κB signalling.