Transcriptome-wide identification and characterization of toll-like receptors response to Vibrio anguillarum infection in Manila clam (Ruditapes philippinarum).

The Manila clam (Ruditapes philippinarum), one of the major marine aquaculture species in China, is susceptible to infection with the pathogen Vibrio, which results in massive mortality and economic losses. Toll-like receptors (TLRs) are significant pattern recognition receptors (PRRs) of innate immunity that are involved in immune regulation against pathogenic invasion. Molecular characterization of Manila clam TLRs and investigations of their immune functions are essential to prevent and control Vibrio infection. In the present research, eight cDNA sequences of R. philippinarum TLRs (RpTLRs) were identified from previous transcriptome libraries and then classified into four groups, namely, P-TLR (one sequence), V-TLR (one sequence), Ls-TLR (two sequences) and sP-TLR (four sequences), based on the corresponding LRR domain arrangement of their protein structures within the typical TLR motifs. A selective pressure test firstly suggested that the molluscan P-TLR, V-TLR, Ls-TLR and sP-TLR families underwent positive selection, and different numbers of positive selection sites (PSSs) were identified in different domains of the four types of RpTLRs, as determined by PAML and analysis of website data. These findings indicated that the evolution of RpTLRs may be associated with their immune recognition and function. Furthermore, tissue-specific expression analysis showed that all RpTLRs were ubiquitously expressed in all test tissues and were dominant in hemocytes. Quantitative real-time PCR revealed that the cDNA expression of all eight RpTLRs was upregulated after injection with Vibrio anguillarum (P < 0.01) in R. philippinarum hemocytes, revealing that these RpTLRs play important roles in responding to pathogenic stimulation. In summary, these findings provide a foundation for future investigations of the molecular classification and evolutionary patterns of Toll-like receptors in invertebrates, and the innate immune responses of TLR signaling pathways in Mollusca.

Transcriptome analysis of Ruditapes philippinarum hepatopancreas provides insights into immune signaling pathways under Vibrio anguillarum infection.

The Manila clam, Ruditapes philippinarum, is one of the most economically important aquatic clams that are harvested on a large scale by the mariculture industry in China. However, increasing reports of bacterial pathogenic diseases have had a negative effect on the aquaculture industry of R. philippinarum. In the present study, the two transcriptome libraries of untreated (termed H) and challenged Vibrio anguillarum (termed HV) hepatopancreas were constructed and sequenced from Manila clam using an Illumina-based paired-end sequencing platform. In total, 75,302,886 and 66,578,976 high-quality clean reads were assembled from 101,080,746 and 99,673,538 raw data points from the two transcriptome libraries described above, respectively. Furthermore, 156,116 unigenes were generated from 210,685 transcripts, with an N50 length of 1125 bp, and from the annotated SwissProt, NR, NT, KO, GO, KOG and KEGG databases. Moreover, a total of 4071 differentially expressed unigenes (HV vs H) were detected, including 903 up-regulated and 3168 down-regulated genes. Among these differentially expressed unigenes, 226 unigenes were annotated using KEGG annotation in 16 immune-related signaling pathways, including Toll-like receptor, NF-kappa B, MAPK, NOD-like receptor, RIG-I-like receptor, and the TNF and chemokine signaling pathways. Finally, 20,341 simple sequence repeats (SSRs) and 214,430 potential single nucleotide polymorphisms (SNPs) were detected from the H and HV transcriptome libraries. In conclusion, these studies identified many candidate immune-related genes and signaling pathways and conducted a comparative analysis of the differentially expressed unigenes from Manila clam hepatopancreas in response to V. anguillarum stimulation. These data laid the foundation for studying the innate immune systems and defense mechanisms in R. philippinarum.

The TLR13-MyD88-NF-κB signalling pathway of Cyclina sinensis plays vital roles in innate immune responses.

Toll-like receptors, the best known pattern recognition receptors, play important roles in recognizing non-self molecules and binding pathogen-associated molecular patterns in the innate immune system. In the present research, the cDNA and protein characterization of the TLR signalling pathway genes including IRAK4, TRAK6 and IKKα (named CsIRAK4, CsTRAF6 and CsIKKα, respectively) with the typical motifs from Cyclina sinensis showed significant similarity with their homologues from other shellfish. Furthermore, the mRNA transcripts of these three genes are ubiquitously expressed in all tissues tested and are dominantly expressed in C. sinensis haemocytes (P < 0.05). Moreover, IRAK4, TRAK6 and IKKα cDNA expression levels were all up-regulated after injection with Vibrio anguillarum, Micrococcus luteus and poly I:C (P < 0.01) as shown by quantitative real-time PCR, indicating that they were involved in responding to pathogenic stimulation. We explored the function of the TLR13-MyD88-NF-κB signalling pathway in the innate immune responses of C. sinensis by RNA interference and immune challenges. The results suggested the mRNA expression patterns of CsMyD88, CsIRAK4, CsTRAF6, CsIKKα, CsIκB, CsNF-κB, CsC-LYZ and CsAMP were all down-regulated (P < 0.01) in normal and stimulated C. sinensis haemocytes, revealing the involvement of the TLR13-MyD88-NF-κB signalling pathway in innate immunity by positively adjusting internal signalling factors and immune-related genes. In summary, a TLR13-MyD88-NF-κB signalling pathway exists and plays vital roles in innate immune responses in C. sinensis. These findings collectively lay the foundation for studying the functional characterization of internal signalling factors and establishing a regulatory network for the TLR signalling pathway in molluscs.

Molecular cloning, characterization and functional analysis of a heat shock protein 70 gene in Cyclina sinensis.

Heat shock protein 70 (HSP70) is an important member of the heat shock protein superfamily and is involved in protecting organisms against various stressors. In the present study, we used RACE to clone a full-length Cyclina sinensis HSP70 cDNA termed CsHSP70. The full length of the CsHSP70 cDNA was 2308 bp, with a 5' untranslated region (UTR) of 42 bp, a 3' UTR of 268 bp, and an open reading frame (ORF) of 1998 bp encoding a polypeptide of 655 amino acids with an estimated molecular mass of 72.75 kDa and an estimated isoelectric point of 5.48. Quantitative real-time PCR was employed to analyze the tissue distribution and temporal expression of the CsHSP70 gene after bacterial challenge and cadmium (Cd) exposure. The CsHSP70 mRNA transcript was expressed ubiquitously in five examined tissues, with the highest expression in hemocytes (P < 0.05) and with the lowest expression in the hepatopancreas. Furthermore, the expression level of CsHSP70 in hemocytes at 3 h after Vibrio anguillarum challenge was extremely significantly up-regulated (P < 0.01). Moreover, the CsHSP70 transcript was up-regulated significantly following exposure to a safe Cd concentration (0.1 mg/L). Finally, after the CsHSP70 gene was silenced by RNA interference, the expression of the CsTLR13 and CsMyD88 genes were extremely significantly decreased (P < 0.01). The results indicated that CsHSP70 could play an important role in mediating the environmental stress and immune responses, and regulating TLR signaling pathway in C. sinensis.

Identification and functional characterization of three TLR signaling pathway genes in Cyclina sinensis.

Toll-like receptors (TLRs) are an ancient family of pattern recognition receptors that play a critical role in initiating and activating the innate immune system. In this study, we identified two TLR genes (CsTLR4 and CsTLR13) and the MyD88 (CsMyD88) gene using a transcriptome library from Cyclina sinensis. The sequence features and mRNA expression profiles of the genes were characterized, and their functions in the immune response were investigated to validate the TLR signaling pathway and its potential role in immune defense. The expression patterns of CsTLR4, CsTLR13 and CsMyD88 were detected in all the tissues examined from healthy clams and were primarily expressed in the hemocytes (P < 0.05), as shown by real-time PCR. Upon challenge with Vibrio anguillarum and Micrococcus luteus, they were significantly increased in hemocytes (P < 0.01), whereas only CsTLR13 and CsMyD88 were up-regulated (P < 0.01) by poly (I:C) challenge. In addition, the mRNA expression level of CsC-LYZ and CsAMP was down-regulated at 72 h (P < 0.01) after injection with CsMyD88 RNAi. These findings might be valuable for understanding the innate immune signaling pathways of C. sinensis and enabling future studies on host-pathogen interactions.

Identification and expression analysis of IκB and NF-κB genes from Cyclina sinensis.

With the increasing economic importance of Cyclina sinensis aquaculture, interest in its defense mechanisms against pathogenic infection has grown in recent years. Inhibitor of nuclear factor-kappaB (IκB) and nuclear factor-kappaB (NF-κB) are proteins with central roles in many important physiological and pathological processes, such as innate immune responses. In this study, we identified CsIκB and CsNF-κB genes from a C. sinensis transcriptome library. In healthy adult clams, CsIκB and CsNF-κB genes were widely expressed in various tissues and highly expressed in hemocytes. Further, the expression levels of these genes were significantly increased in hemocytes challenged by Vibrio anguillarum, Micrococcus luteus and poly I:C. Inhibition of CsMyD88 expression by RNAi technology significantly altered the mRNA expression patterns of CsIκB and CsNF-κB as measured using quantitative real-time PCR. These results collectively indicated that the NF-κB signaling pathway, including CsIκB and CsNF-κB genes, might be involved in early innate immune responses and may be regulated by a MyD88-dependent signaling pathway in C. sinensis.

A Kazal-type serine proteinase inhibitor from Cyclina sinensis is involved in immune response and signal pathway initiation.

Serine protease inhibitors (SPIs) are an important group of protease inhibitors involved in a variety of biological processes. In the present study, a Kazal-type serine protease inhibitor homolog gene (designated as CsKPI) was identified from a Cyclina sinensis cDNA library. The open reading frame consists of 456 bp and encodes a protein of 151 amino acid residues with a theoretical molecular mass of 16.85 kDa and an isoelectric point of 5.74. Furthermore, using quantitative real-time PCR, we focused on the expression patterns of CsKPI found in tissues and on the stimulation of this gene's expression by bacteria. The results show that a higher-level mRNA expression of CsKPI was detected in hemocytes (P < 0.05) and was significantly upregulated at 3 h (P < 0.01) upon receiving bacterial challenges with Vibrio anguillarum. In addition, after the CsKPI gene was silenced by RNA interference, the expression of the CsTLR2 and CsMyD88 genes was extremely significantly decreased (P < 0.01) in C. sinensis. Finally, the recombinant CsKPI (rCsKPI) protein was purified and shown to exhibit less inhibitory activity than C-lyz against V. anguillarum in vitro. Hence, we propose that CsKPI plays an important role in the innate immunity and mediates TLR2 and MyD88-dependent pathway initiation in C. sinensis.

Molecular characterization and expression analysis of Toll-like receptor 21 cDNA from Paralichthys olivaceus.

Toll-like receptor (TLR) is believed to play crucial role in host defense of pathogenic microbes in innate immune system. In the present study, the full-length cDNA of Paralichthys olivaceus Toll-like receptor 21 (Po-TLR21) was cloned by homology cloning and rapid amplification of cDNA ends (RACE) technique. The Po-TLR21 cDNA sequence was 3687 bp, containing an open reading frame of 2922 bp encoding 973 amino acids. TMHMM and SMART program analysis indicated that protein contained one transmembrane domain, eighteen leucine-rich repeats (LRRs), and one Toll/IL-1 receptor homology domain (TIR). Multiple alignment analysis of the Po-TLR21 protein-coding sequence with other known TLR21 from grouper, pufferfish, zebrafish, cod, catfish, carp and chicken showed the homology of 67%, 63%, 54%, 52%, 51%, 49%, and 39%, respectively. The Po-TLR21 mRNA expression patterns were measured by real-time PCR. The results revealed that TLR21 is widely expressed in various tested healthy tissues, and highly expressed in spleen and gill. In vivo immunostimulation experiments revealed that expression of TLR21 is modulated by Vibrio anguillarum (V. anguillarum), CpG oligodeoxynucleotides (CpG ODN) and poly I:C. Moreover, the inhibitor of homodimerization of myeloid differentiation factor 88 (MyD88) could significantly reduce the up-regulation of TLR21, MyD88, and tumor necrosis factor (TNF) expression in CpG ODN or poly I:C-treated head kidney cells in vitro. These results indicate that TLR21 may be involved in the pathogen recognition in the early innate immune.

Molecular characterization of a novel p38 MAPK cDNA from Cyclina sinensis and its potential immune-related function under the threat of Vibrio anguillarum.

The p38 mitogen-activated protein kinase (MAPK) is one important member of MAPK family and reported to serve a predominant function in regulating innate immunity after the occurrence of certain infection. In the present study, one novel p38 MAPK gene was acquired from Cyclina sinensis based on the RNA-seq analysis and designated as Csp38 MAPK. This novel gene contained a full length of 1781 bp, 1104 bp of which was deemed as open reading frames and gave rise to a peptide of 367 amino acids with a predicted molecular weight of 42.31 KDa. A conserved serine/threonine protein kinase (S_Tkc) region along with a Thr-Gly-Tyr motif was discovered in the deduced sequence. According to the phylogenetic analysis, there was a close relationship between this kinase and Meretrix petechialis p38 MAPK. As for the expression pattern, this newly-identified Csp38 MAPK was ubiquitously distributed in several tissues throughout the body but with varied abundance. After the challenge of Vibrio anguillarum, both the transcription and phosphorylation level of Csp38 MAPK in hemolymph were coordinately altered with a time-dependent manner. Besides, with the application of double strand RNA homologous to myeloid differentiation factor 88 (MyD88) of C. sinensis, the activation of Csp38 MAPK was found to obviously decrease in hemolymph after the pathogen stimulation. Hence, our experimental data presented evidence for the potential involvement of p38 MAPK in response to bacterial invaders in C. sinensis, possibly facilitating the understanding for pathogen-induced innate immunity in clams.

Molecular and expression characterization of Toll-like receptor family genes from the Anadara sativa (Bivalvia, Arcidae) transcriptome.

Innate immunity plays an important role in invertebrates because it provides the first line of protection by recognizing invading microbial pathogens and then activating downstream signaling pathways. However, until now, increasing reports of clam diseases did not include those of Anadara sativa, which are widely distributed and economically important maritime clams. In the present study, transcriptome libraries of untreated (termed H) and Vibrio anguillarum-challenged (termed HV) A. sativa hepatopancreases were constructed and sequenced using the Illumina HiSeq4000 platform. In total, we obtained 78,012,510 and 84,937,516 clean reads from 80,006,030 to 86,871,742 raw data reads, respectively, assembled by different software programs. Furthermore, 150,274 unigenes were generated from 196,003 transcripts, with an N50 length of 1088 bp, and then annotated with the SwissProt, NR, NT, PFAM, KO, GO, KOG and KEGG databases. Moreover, 3982 differentially expressed unigenes (H vs HV) were determined, with 3583 upregulated and 399 downregulated genes. Among these differentially expressed unigenes, 207 unigenes were found using KEGG annotation in 16 immune-related signaling pathways, such as Toll-like receptor (TLR), NOD-like receptor (NLR), and RIG-I-like receptor (RLR) signaling pathways. Finally, we selected 11 full-length TLRs and classified them into 3 groups, namely, one V-TLR, four Ls-TLR and six sP-TLR; furthermore, we validated the increased expression patterns of the 11 TLRs in response to LPS injection. In summary, these results revealed multiple findings on potential immune-related genes, such as the differential expression analysis and annotation based on the A. sativa transcriptome in response to V. anguillarum stimulation, and explored the molecular and expression characterization of A. sativa TLRs, which provide new insights into the innate immune responses and defense mechanisms in shellfish.

Evolutionary and functional analysis of Cyclina sinensis c-Jun AP-1 gene in response to LPS stimulation.

The transcription factor activator protein-1 (AP-1) plays an essential and critical role in the regulation of numerous downstream genes involved in various physiological and chemical responses. In this study, we identified a full-length cDNA of the c-Jun AP-1 gene (termed Csc-Jun) from the transcriptome library in Cyclina sinensis. The cDNA contains an 825-bp open reading frame that encodes a 274-amino acid protein sequence, including a characteristic Jun transcription factor domain and a highly conserved basic leucine zipper (bZIP) signature that shares 90% identity to that of Ruditapes philippinarum. Furthermore, a phylogenetic analysis using MrBayes and PhyML software (with Bayesian and maximum likelihood approaches, respectively) revealed that the c-Jun AP-1 family genes might be involved in adapting to various environments in different invertebrates. We implemented the PAML software with the maximum likelihood method to further select and verify the positive selection sites (PSSs) in the Mollusca c-Jun AP-1 genes, and we detected four PSSs located in the Jun transcription factor domain. In addition, a spatial expression analysis showed that the Csc-Jun cDNA transcript was ubiquitously expressed in all of the tested tissues and was strongly expressed in the hepatopancreas and weakly expressed in the tissues of the hemocytes, gill filaments, mantle and adductor muscle. Quantitative real-time PCR showed that the expression profiles of Csc-Jun were significantly upregulated at different times in all of the tested tissues when challenged with lipopolysaccharide (LPS). Furthermore, knockdown of Csc-Jun by RNA interference resulted in a higher mortality of C. sinensis following LPS exposure. Finally, we explored the function of the TLR13-MyD88 signaling pathway in the innate immunity of C. sinensis by RNA interference and immune challenges. The results revealed that the mRNA expression levels of Csc-Jun were all decreased (P < 0.01) in normal and stimulated C. sinensis hemocytes. These data collectively indicated that the c-Jun AP-1 gene might play vital roles in innate immunity and provide new evidence for the evolutionary patterns of innate immune genes in Mollusca.

Comparative and evolutionary analysis of an adapter molecule MyD88 in invertebrate metazoans.

The myeloid differentiation factor 88 (MyD88) is an essential adapter in Toll-like receptor (TLR) signalling pathways, with TLR the first pattern-recognition receptor (PRR) that was discovered in Drosophila. In the present study, a MyD88 gene was identified and characterized from a commercially important shellfish, Scapharca subcrenata, including a DEATH domain and TIR domain conserved within other molluscs. Furthermore, comparative genomic evidence revealed that MyD88 was of different lengths and contained quantitative exon and intron regions, which might be involved in specific mechanisms. To further explore the phylogenetic relationships of invertebrate metazoan MyD88, we applied MrBayes and PhyML software to construct phylogenetic trees using Bayesian and maximum likelihood approaches, respectively, which suggested that the MyD88 of Arthropoda was closely related to lower invertebrates, in contrast to morphological taxonomy. Finally, we investigated the evolutionary patterns and location of positive selection sites (PSSs) in the MyD88 gene from Arthropoda, Mollusca and Insecta using PAML software with the maximum likelihood method. The data showed that positive selection sites were detected in these groups, and partial sites were located in the TIR domain but were not found in the DEATH domain. To summarize, in this study, we report on the diversification of MyD88 in invertebrate metazoans, the specific evolutionary position of Arthropoda MyD88, and the positive selection pressures on MyD88 of Arthropoda, Mollusca and Insecta. These results are a valuable contribution to understand and clarify the evolutionary pattern of TLR/MyD88 signalling pathways in invertebrate and vertebrate taxa.

A new record and two new species of Cryptochironomus Kieffer, 1918 from China (Diptera, Chironomidae).

Two new species of the genus Cryptochironomus Kieffer, C. maculus Yan & Wang sp. n. and C. protuberans Yan & Wang sp. n. are described and illustrated as adult males. C. albofasciatus (Staeger) is recorded from China for the first time. A key to the males of Cryptochironomus in China is presented and generic diagnosis is emended.

Review of the genus Harnischia Kieffer from China (Diptera, Chironomidae), with description of one new species.

The genus Harnischia Kieffer, 1921 from China is reviewed and one new species, Harnischia parallela Yan & Wang, sp. n., is described and illustrated as adult male. Harnischia okilurida Sasa is recognized as a new synonym of Harnischia longispuria Wang & Zheng. The pupae of Harnischia fuscimana Kieffer and Harnischia curtilamellata Malloch are redescribed from material collected in China, and an unplaced pupa is described. Key to male adults and pupae of known species of Harnischia from China is given.

De novo RNA-Seq analysis of the venus clam, Cyclina sinensis, and the identification of immune-related genes.

The Venus clam, Cyclina sinensis, is one of the most important bivalves in China. In recent years, increasing expansive morbidity has occurred in breeding areas, imposing significant losses on the national economy. To understand the molecular mechanisms of immune-related genes, we analyzed and sequenced hemolymph samples that were injected with two pathogenic microorganisms using the Illumina Miseq system. After trimming, more than 12 M PE reads with an average length greater than 410 bp were assembled into 70,079 transcripts with a mean length of 980 bp. Using a homology analysis, 102 (135 transcripts) potentially immune-related genes were identified, and most of them exhibited a similar pattern in both samples. These data indicated that the response of the clam to both types of bacterial infection might follow a similar molecular mechanism. Using the TreeFam method, 9,904 gene families and 1,031 unique families of the clam were preliminarily classified in comparison to five related species. A significant number of SSRs were identified, which could facilitate the identification of polymorphisms in Venus clam populations. These datasets will improve our knowledge of the molecular mechanisms driving the immune response to bacterial infection in clam populations and will provide basic data about clam breeding and disease control.

Parachironomus lenz from china and Japan (Diptera, chironomidae).

Members of the genus Parachironomus Lenz known from China and Japan are revised, and a key to their male adults is given. Parachironomuspoyangensis sp. n. is described in this life stage. Parachironomusfrequens (Johannsen) and Parachironomusmonochromus (van der Wulp) are recorded from China for the first time, thus are redescribed from Chinese specimens. Parachironomuskamaabeus Sasa & Tanaka and Parachironomustoneabeus Sasa & Tanaka are new junior synonyms of Parachironomusfrequens. Three Chinese or Japanese species formerly placed in Parachironomus are transferred to other genera, resulting in the new combinations Cryptochironomusinafegeus (Sasa, Kitami & Suzuki), Demicryptochironomus (Irmakia) lobus (Yan, Sæther, Jin & Wang), and Microchironomuslacteipennis (Kieffer). Chironomussauteri Kieffer, Parachironomuskisobilobalis Sasa & Kondo and Parachironomuskuramaexpandus Sasa are removed from Parachironomus; the last of these three denotes a valid species of uncertain generic placement, the first two are nomina dubia.