Genome-wide identification and characterization of five MyD88 duplication genes in Yesso scallop (Patinopecten yessoensis) and expression changes in response to bacterial challenge

A newly identified scallop MyD88 interacts with TLR and functions in innate immunity

Myeloid differentiation factor-88 (MyD88) is a key adaptor of the toll-like receptor (TLR) signaling pathway and plays a crucial role in innate immune signal transduction in animals. However, the MyD88-mediated signal transduction mechanism in shellfish has not been well studied. In this study, a new MyD88 gene, CfMyD88-2, was identified in the Zhikong scallop, Chlamys farreri. The 1779 bp long open reading frame encodes 592 amino acids. The N-terminus of CfMyD88-2 contains a conserved death domain (DD), followed by a TIR (TLR/Interleukin-1 Receptor) domain. The results of the multi-sequence comparison showed that the TIR domain sequences were highly conserved. Phylogenetic analysis revealed that CfMyD88-2 was first associated with Mizuhopecten yessoensis MyD88-4 and Argopecten irradians MyD88-4. CfMyD88-2 mRNA was expressed in all scallop tissues, as detected by qRT-PCR, and the expression level was the highest in the mantle and hepatopancreas. In addition, CfMyD88-2 mRNA expression significantly increased after pathogen-associated molecular patterns (PAMPs, such as lipopolysaccharide, peptidoglycan, or polyinosinic-polycytidylic acid) stimulation. The results of the co-immunoprecipitation experiments in HEK293T cells showed that both CfMyD88-1 and CfMyD88-2 interacted with the TLR protein of scallops, suggesting the existence of more than one functional TLR-MyD88 signaling axis in scallops. Dual luciferase reporter gene assays indicated that the overexpressed CfMyD88-2 in HEK293T cells activated interferon (IFN) α, IFN-β, IFN-γ, and NF-κB reporter genes, indicating that the protein has multiple functions. The results of the subcellular localization experiment uncovered that CfMyD88-2 was mainly localized in the cytoplasm of human cells. In summary, the novel identified CfMyD88-2 can respond to the challenge of PAMPs, participate in TLR immune signaling, and may activate downstream effector genes such as NF-κB gene. These research results will be useful in advancing the theory of innate immunity in invertebrates and provide a reference for the selection of disease-resistant scallops in the future.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Molecular Characterization of MyD88 in Anodonta woodiana and Its Involvement in the Innate Immune Response to Bacterial Infection

Myeloid differentiation factor 88 (MyD88) is a key adapter molecule in Toll-like receptor signal transduction that triggers downstream immune cascades involved in the host defense response to exogenous pathogens. However, the function of MyD88s in mollusks, especially in freshwater shellfish, remains poorly understood. In this study, a novel freshwater shellfish MyD88 (denoted AwMyD88) was characterized from Anodonta woodiana. The present AwMyD88 protein consists of 474 amino acids and contains a conserved a typical death domain (DD) and a conservative Toll/IL-1R (TIR) domain with three typical boxes. Quantitative real-time PCR (qRT-PCR) analysis showed that AwMyD88 was broadly expressed in all the examined tissues, and the highest expression level was observed in hemocytes of A. woodiana. When challenged with Aeromonas hydrophila and lipopolysaccharide (LPS), the mRNA expression levels of AwMyD88 were significantly induced in hemocytes of A. woodiana in vivo and in vitro. In addition, in vivo injection experiments revealed that MyD88 signaling pathway genes showed strong responsiveness to A. hydrophila challenge, and their expression levels were significantly upregulated in hemocytes. Knockdown of AwMyD88 reduced the transcript levels of immune related transcription factors (AwNF-κB and AwAP-1) and effectors (AwTNF, AwLYZ, AwDefense and AwAIF) during A. hydrophila infection. Moreover, subcellular localization analysis indicated that AwMyD88 was mainly localized to the cytoplasm in HEK293T cells. Finally, luciferase reporter assays revealed that AwMyD88 associates with AwTLR to activate the NF-κB and AP-1 signaling pathways in HEK293T cells. These results suggested that AwMyD88 might be involved in the host defense response to bacterial challenge, providing new insight into the immune function of the MyD88 signaling pathway in freshwater shellfish.

Identification, functional characterization and expression pattern of myeloid differentiation factor 88 (MyD88) in Nibea albiflora

Myeloid differentiation factor 88 (MyD88), composed of an N-terminal death domain and a C-terminal Toll/interleukin (IL)-IR homology domain, is a key connector protein in the TLR signal transduction pathway. In this study a novel isoform of MyD88 in Nibea albiflora (named as NaMyD88) was identified and functionally characterized (GenBank accession no. MN384261.1). Its complete cDNA sequence was 1672 bp and contained an open reading frame of 879 bp encoding 292 amino acid residues, which was similar to its teleost fish counterparts in the length. The theoretical molecular mass was 33.63 kDa and the isoelectric point was 5.24. BLASTp analysis suggested that the deduced amino acids sequence of NaMyD88 shared high identity to the known MyD88, for instance, 94.77% identity with Collichthys lucidus. Sequence analysis showed that NaMyD88 protein was consistent with MyD88 protein of other species at three conserved domains, N-terminal DD, short middle domain and C-terminal TIR, and the TIR domain contained three highly conserved motifs: Box1, Box2, and Box3. NaMyD88 and red fluorescent protein (Dsred) were fused and expressed in the cytoplasm of the epithelioma papulosum cyprini (EPC cells). The NaTLR9-TIR-EGFP fusion protein, which was obtained in our previous studies, showed green fluorescence and mainly distributed in the cytoplasm. After co-transfection, NaMyD88-Dsred and NaTLR9-TIR-EGFP obviously overlapped and displayed orange-yellow color. The results showed that the homologous MyD88-Dsred could interact with NaTLR9-TIR-EGFP. Based on this result pcMV-NaMyD88-TIR-Myc plasmids and the pcDNA3.1-NaTLR9-TIR-flag were constructed and co-transfected into 293T cells for the immunoprecipitation test. According to Western blot, the protein eluted by Flag-beads could be detected by anti-Flag-tag antibody and anti-Myc tag antibody respectively, while the protein without NaTLR9-TIR could not be found, which further proved that TLR and MyD88 could interact each other. The prokaryotic plasmid of MyD88-TIR domain was constructed, expressed in BL21 (DE3) and purified by Ni-NAT super flow resin conforming to the expected molecular weight of 27 kDa with the corresponding active sites for its conferring protein-protein interaction functions. Real-time fluorescence quantitative PCR showed that NaMyD88 could be expressed in intestine, stomach, liver, kidney, gill, heart and spleen, with the highest in the kidney, and it was up-regulated after being infected with Polyinosinic:polycytidylic acid - Poly (I:C) and Pseudomonas plecoglossicida, which showed that NaMyD88 was involved in the immune response of N.albiflora. These data afforded a basis for understanding the role of NaMyD88 in the TLR signaling pathway of N.albiflora.

Molecular allocation of PC4s provides implications for deciphering thermal response in Zhikong scallop (Chlamys farreri)

The increasing sea temperature caused by global warming has led to serious death of Zhikong scallop (Chlamys farreri) and improving its thermal tolerance has become an active research area in scallop aquaculture industry. Gene transcriptional coactivator p15 (PC4) plays pivotally multi-faced roles in most vertebrates and some invertebrates, but the systematic identification and characterization of PC4 genes have less been reported in scallops. In this study, 15 PC4 genes (CfPC4s) were identified in Zhikong scallop through whole-genome scanning, including two pairs of tandem duplicate genes located in the same scaffold (CF-19495.9 and CF-19495.10, CF-6819.1 and CF-6819.2). Protein structural and phylogenetic analyses were performed to verify identities and evolutionary relationships of these genes. Spatiotemporal expression patterns were determined at different development stages and in healthy adult tissues, as well as expression regulations in selected tissues (mantles, gills, hemocytes and hearts) after high temperatures challenge (27 °C) with different durations (3 h, 6 h, 12 h, 24 h, 3 d, 6 d, 15 d and 30 d). Spatiotemporal expressions of CfPC4s were ubiquitous but exhibited different patterns, suggesting the functional roles of CfPC4s in all stages of growth and development of the scallop. Expression regulations of CfPC4s and their functional related factors (TFIIA, TFIID, TFIIH and RNAPII) in pre-initiation complex (PIC) in various tissues displayed up- and/or down-regulated responses at different time points, showing time- and/or tissue-dependent expression patterns with function allocation upon different thermal durations. Collectively, this study demonstrated that gene allocation of CfPC4s provided implications for deciphering thermal response in Zhikong scallop and potentially helped in developing strategies for long-term healthy sustainable Zhikong scallop culture.

Molecular identification and functional analysis of MyD88 in giant freshwater prawn (Macrobrachium rosenbergii) and expression changes in response to bacterial challenge

Myeloid differentiation factor 88 (MyD88) is a crucial adaptor protein for Toll-like receptor (TLR)-mediated signaling pathways and plays an important role in immune response. In this study, the full-length cDNA of MyD88 from Macrobrachium rosenbergii (MRMyD88) was cloned. The MRMyD88 cDNA is 1758 bp long and contains a 1398-bp open reading frame. Multiple sequence alignment and phylogenetic analysis revealed that the amino acid sequence of MRMyD88 shared high identity with the known MyD88 proteins. The MRMyD88 mRNA was widely expressed in all examined tissues, with highest level in intestine, followed by gonad and pleopod. Furthermore, the MRMyD88 promoter region, spanning 1622 bp, contains several transcription factor-binding sites, including nine GATA-1 box motifs. Electrophoretic mobility shift assay showed that Gfi-1, SRF, and Oct-1 bind to the upstream region of MRMyD88. Additionally, the results showed that the expression levels of TLR1, TLR2 and TLR3 were different in response to Vibrio anguillarum, Lactobacillus plantarum and Aeromonas hydrophila infections. However, these bacteria significantly increased the expression levels of MyD88 and prophenoloxidase. These data suggest that the TLR-mediated signaling pathway is MyD88-dependent in response to pathogenic and probiotic bacteria in M. rosenbergii.

Microplastics induce dose-specific transcriptomic disruptions in energy metabolism and immunity of the pearl oyster Pinctada margaritifera

A combined approach integrating bioenergetics and major biological activities is essential to properly understand the impact of microplastics (MP) on marine organisms. Following experimental exposure of polystyrene microbeads (micro-PS of 6 and 10 μm) at 0.25, 2.5, and 25 μg L^-1, which demonstrated a dose-dependent decrease of energy balance in the pearl oyster Pinctada margaritifera, a transcriptomic study was conducted on mantle tissue. Transcriptomic data helped us to decipher the molecular mechanisms involved in P. margaritifera responses to micro-PS and search more broadly for effects on energetically expensive maintenance functions. Genes related to the detoxification process were impacted by long-term micro-PS exposure through a decrease in antioxidant response functioning, most likely leading to oxidative stress and damage, especially at higher micro-PS doses. The immune response was also found to be dose-specific, with a stress-related activity stimulated by the lowest dose present after a 2-month exposure period. This stress response was not observed following exposure to higher doses, reflecting an energy-limited capacity of pearl oysters to cope with prolonged stress and a dramatic shift to adjust to pessimum conditions, mostly limited and hampered by a lowered energetic budget. This preliminary experiment lays the foundation for exploring pathways and gene expression in P. margaritifera, and marine mollusks in general, under MP exposure. We also propose a conceptual framework to properly assess realistic MP effects on organisms and population resilience in future investigations.

Molecular cloning and expression analysis of myd88 from oriental weatherfish (Misgurnus anguillicaudatus) in response to bacterial challenge

Myeloid differentiation factor 88 (Myd88) plays an important role in both innate and adaptive immune response. In this study, the full-length complementary DNA (cDNA) of myd88 from Misgurnus anguillicaudatus was characterized. The myd88 cDNA is 1333 bp in length and contains an 855 bp open reading frame encoding a predicted protein of 284 amino acids. The predicted protein possesses typical Myd88 domain structural features including a death domain in the N-terminus, and box 1, 2, and 3 motifs of the Toll/IL-1 receptor domain in the C-terminus. Quantitative real-time PCR (qRT-PCR) revealed that myd88 messenger RNA (mRNA) was ubiquitously expressed in all examined tissues, especially highly in brain, kidney, blood, intestines and liver. qRT-PCR and western blotting were used to determine the mRNA and protein levels of Myd88 after Aeromonas veronii challenge, respectively. The Myd88 was remarkably upregulated in response to infection of A. veronii. These results suggested that Myd88 may play a vital role during the immune response of M. anguillicaudatus against bacterial infection.

Interactive effects of salinity variation and exposure to ZnO nanoparticles on the innate immune system of a sentinel marine bivalve, Mytilus edulis

ZnO nanoparticles (nZnO) are released into the coastal environment from multiple sources, yet their toxicity to marine organisms is not well understood. We investigated the interactive effects of salinity (normal 15, low 5, and fluctuating 5-15) and nZnO (100 μg l-1) on innate immunity of the blue mussels Mytilus edulis from a brackish area of the Baltic Sea. Exposure to ionic Zn (100 μg l-1) was used to test whether the toxic effects of nZnO can be attributed to the potential release of Zn2+. Functional parameters and the expression of key immune-related genes were investigated in the mussels exposed to nZnO or ionic Zn under different salinity regimes for 21 days. nZnO exposures elevated hemocyte mortality, suppressed adhesion, stimulated phagocytosis, and led to an apparent increase in lysosomal volume. At salinity 15, nZnO suppressed the mRNA expression of the Toll-like receptors TLRb and c, C-lectin, and the complement system component C3q indicating impaired ability for pathogen recognition. In contrast, the mRNA levels of an antimicrobial peptide defensin increased during nZnO exposure at salinity 15. At fluctuating salinity (5-15), nZnO exposure increased expression of multiple immune-related genes in hemocytes including the complement system components C1 and C3q, and the Toll-like receptors TLRa, b and c. Low salinity (5) had strong immunosuppressive effects on the functional and molecular immune traits of M. edulis that overshadowed the effects of nZnO. The salinity-dependent modulation of immune response to nZnO cannot be attributed to the differences in the aggregation or solubility of nZnO, and likely reflects the interaction of the toxic effects of nanoparticles and physiological effects of the osmotic stress. These findings have implications for the environmental risk assessment of nanomaterials and the development of the context-specific biomarker baselines for coastal pollution monitoring.

Evolutionary and functional analysis of MyD88 genes in pearl oyster Pinctada fucata martensii

Myeloid differentiation factor 88 (MyD88) is an adapter protein that links toll-like receptor and interleukin 1 receptor-mediated signal transduction. In this study, we identified 20 MyD88 genes from eight mollusk genomes and found that MyD88 was expanded in bivalves. This expansion tends to be tandem duplication. Phylogenetic analysis suggested that the tandem duplication of MyD88 was formed before bivalve differentiation. All of the identified MyD88 contained both of death domain (DD) and toll/interleukin-1 receptor (TIR) domain, and 13 mollusks MyD88 have low complexity regions (LCRs), which were not found in the MyD88 from humans and zebrafish. The genomic structure showed that most of the mollusk MyD88 (14 of 19) contained five conserved introns, four of which were found in humans and zebrafish. Furthermore, the cDNA full length of PfmMyD88-2 (one of the two identified MyD88 in Pincatada fucata martensii) was obtained with 1591 bp, including 260 bp of 5'UTR, 257 bp of 3'UTR, and 1077 bp of open reading frame encoding 358 amino acids. Quantitative real-time PCR analysis demonstrated that PfmMyD88-2 mRNA was widely expressed in all detected tissues. The highest expression level was in the gills and followed by hepatopancreas and feet. After lipopolysaccharide stimulation, PfmMyD88-2 expression level increased and reached the highest level at 12 h and then gradually declined to the normal level. Over-expression of PfmMyD88-2 in HEK293T increased the luciferase activity of the pNF-κB-Luc reporter. We also identified that PfmmiR-4047 could regulate the expression of PfmMyD88-2. These results help us elucidate the mechanism underlying mollusk immune response.

Genome-wide identification, characterization of RLR genes in Yesso scallop (Patinopecten yessoensis) and functional regulations in responses to ocean acidification

Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), are crucial sensors with a conserved structure in cytoplasm, inducing the production of cytokines, chemokines and host restriction factors which mediate a variety of intracellular activities to interfere with distinct PAMPs (pathogen-associated molecular patterns) for eliminating pathogens in innate immune system. Although RLR genes have been investigated in most vertebrates and some invertebrates, the systematic identification and characterization of RLR genes have not been reported in scallops. In this study, four RLR genes (PY-10413.4, PY-10413.5, PY-443.7 and PY-443.8, designated PyRLRs) were identified in Yesso scallop (Patinopecten yessoensis) through whole-genome scanning through in silico analysis, including two pairs of tandem duplicate genes located on the same scaffold (PY-10413.4 and PY-10413.5, PY-443.7 and PY-443.8, respectively). Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of these genes. The expression profiles of PyRLRs were determined in all developmental stages, in healthy adult tissues, and in mantles that simulated ocean acidification (OA) exposure (pH = 6.5 and 7.5) at different time points (3, 6, 12 and 24 h). Spatiotemporal expression patterns suggested the functional roles of PyRLRs in all stages of development and growth of the scallop. Regulation expressions revealed PY-10413.4 and PY-10413.5 with one or two CARD(s) (caspase activation and recruitment domain) were up-regulated expressed at most time points, whereas PY-443.8 and PY-10413.4 without CARD were significantly down-regulated at each time points, suggesting functional differentiations in the two pairs of PyRLRs based on the structural differences in response to OA. Collectively, this study demonstrated gene duplication of RLR family genes and provide primary analysis for versatile roles in the response of the bivalve innate immune system to OA challenge.

Alternative splicing, spatiotemporal expression of TEP family genes in Yesso scallop (Patinopecten yessoensis) and their disparity in responses to ocean acidification

The complement system constitutes a highly sophisticated and powerful body defense machinery acting in the innate immunity of both vertebrates and invertebrates. As central components of the complement system, significant effects of thioester-containing protein (TEP) family members on immunity have been reported in most vertebrates and in some invertebrates, but the spatiotemporal expression and regulatory patterns of TEP family genes under environmental stress have been less widely investigated in scallops. In this study, expression profiling of TEP family members in the Yesso scallop Patinopecten yessoensis (designated PyTEPs) was performed at all developmental stages, in different healthy adult tissues, and in mantles during exposure to different levels of acidification (pH = 6.5 and 7.5) for different time points (3, 6, 12 and 24 h); this profiling was accomplished through in silico analysis of transcriptome and genome databases. Spatiotemporal expression patterns revealed that PyTEPs had specific functional differentiation in all stages of growth and development of the scallop. Expression analysis confirmed the inducible expression patterns of PyTEPs during exposure to acidification. Gene duplication and alternative splicing events simultaneously occurred in PyTEP1. Seven different cDNA variants of PyTEP1 (designated PyTEP1-A-PyTEP1-G) were identified in the scallop mantle transcriptome during acidic stress. These variants were produced by the alternative splicing of seven differentially transcribed exons (exons 18-24), which encode the highly variable central region. The responses to immune stress may have arisen through the gene duplication and alternative splicing of PyTEP1. The sequence diversity of PyTEP1 isoforms and their different expression profiles in response to ocean acidification (OA) suggested a mechanism used by scallops to differentiate and regulate PyTEP1 gene expression. Collectively, these results demonstrate the gene duplication and alternative splicing of TEP family genes and provide valuable resources for elucidating their versatile roles in bivalve innate immune responses to OA challenge.

Northeast Chinese lamprey (Lethenteron morii) MyD88: Identification, expression, and functional characterization

Myeloid differentiation factor 88 (MyD88) is a key adaptor of Toll-like receptors (TLR), an important pattern recognition receptor of the innate immune system. To study the origin and evolution of the vertebrate TLR signaling pathway in innate immune systems, we analyzed the biological characteristics and functions of the MyD88 gene in Northeast Chinese lamprey (Lethenteron morii) using PCR amplification, real-time PCR analysis, dual luciferase reporter gene assay, immunofluorescence assay, and other methods. Bioinformatics analysis showed that LmMyD88 has a modular structure consisting of Toll/IL-1R domain (TIR) and death domain (DD), which is typical of the MyD88 family. A phylogenetic tree showed that the homology of LmMyD88 was consistent with the phylogenetic status of lampreys. Tissue expression analysis indicated that the mRNA expression was expressed in some normal tissues of larval and adult L. morii. Real-time PCR analysis showed that the expression of LmMyD88 in tissues, such as gill and kidney, of the adult increased significantly after infection by Pseudomonas aeruginosa. Subcellular localization results showed that LmMyD88 was expressed in the nucleus, cytoplasm, and other parts. A dual luciferase reporter assay indicated that LmMyD88 activated nuclear factor kappa B downstream of the TLR signaling pathway. This study suggested that LmMyD88 might play an important role in the innate immune signal transduction process of L. morii.

Identification and functional characterization of three myeloid differentiation factor 88 (MyD88) isoforms from thick shell mussel Mytilus coruscus

Myeloid differentiation factor 88 (MyD88) is a pivotal adapter protein that involved in interleukin-1 receptor/toll-like receptor (IL-1R/TLR) signal transduction, which could spur downstream cascades and eventually drawn into innate immune response. MyD88 has been extensively studied in vertebrates, however, the information ascribe to MyD88 in invertebrates is still very scarce especially its function annotation remains extremely obscure. At here, three novel MyD88 isoforms termed McMyD88a, McMyD88b and McMyD88c were firstly cloned from thick shell mussel Mytilus coruscus. McMyD88a, McMyD88b and McMyD88c shared domain topology containing the Death domain (DD) and TIR domain (TIR) with its counterparts in mammals. All three McMyD88s were ubiquitously expressed in examined tissues in thick shell mussel, with the higher expression levels in immune-related tissues such as haemocytes, gills and digestive glands. Upon Vibrio alginolyticus, polyinosine-polycytidylic acid (poly I:C) and lipopolysaccharide (LPS) challenge, McMyD88a, McMyD88b and McMyD88c transcripts were significantly induced in haemocytes despite of differential expression levels and responsive time points. Overexpression of McMyD88a, McMyD88b and McMyD88c showed a dose-dependent induction to NF-κB or ISRE in mammalian cell lines. Taken together, these results suggested that McMyD88a, McMyD88b and McMyD88c are members of MyD88 family and play potential roles in innate immune response to pathogenic invasions in thick shell mussel. Moreover, these results suggested indirectly the existence of a MyD88-dependent signaling pathway in thick shell mussel, and provide insight into the immunoregulatory effect in molluscs.

Expression of Cathepsin F in response to bacterial challenges in Yesso scallop Patinopecten yessoensis

Cathepsin F is a unique papain cysteine proteinase with highly conserved structures: catalytic triad and a cystatin domain contained in the elongated N-terminal pro-region. It has been reported that cathepsin F is associated with the establishment of innate immune in several vertebrate including fish in aquaculture, but not known in bivalves. In this study, we firstly identified and characterized cathepsin F in the Yesso scallop (Patinopecten yessoensis). The protein structural and phylogenetic analyses were then conducted to determine its identity and evolutionary position. We've also investigated the expression levels of cathepsin F gene at different embryonic developmental stages, in healthy adult tissues and especially in the hemocytes and hepatopancreas after Gram-positive (Micrococcus luteus) and negative (Vibrio anguillarum) challenges using quantitative real-time PCR (qPCR). Cathepsin F was significantly up-regulated 3 h after infection of V. anguillarum in hemocytes, suggesting its participation in immune response. Our findings have provided strong evidence that cathepsin F may be a good target for enhancing the immune activity in Yesso scallop.

The GRP94 gene of Yesso scallop (Patinopecten yessoensis): Characterization and expression regulation in response to thermal and bacterial stresses

The 94-kDa glucose-regulated protein (GRP94) belonging to the HSP90 family is an endoplasmic reticulum (ER) chaperone. It plays critical roles in ER quality control, and has been implicated as a specialized immune chaperone to regulate both innate and adaptive immunity. In this study, we identified and characterized a GRP94 gene (PyGRP94) from Yesso scallop (Patinopecten yessoensis). The protein sequence of PyGRP94 is highly conserved with its homologs in vertebrates, with a signal sequence in N-terminal, an ER retrieval signal sequence in C-terminal and a HATPase_c domain. Expression analysis suggests that PyGRP94 transcripts in early embryos are maternally derived and the zygotic expression is started from D-shaped larvae. This gene is also expressed in almost all the adult tissues examined except smooth muscle, with the highest expression level in hemocytes. Besides, PyGRP94 was demonstrated to be induced by heat shock and both Gram-positive (Micrococcus luteus) and Gram-negative (Vibrio anguillarum) bacterial infection, with much more dramatic changes being observed after V. anguillarum challenge. Our results suggest the involvement of PyGRP94 in response to thermal stress, and that it might play an important role in the innate immune defense of scallop.

Identification and characterization of TEP family genes in Yesso scallop (Patinopecten yessoensis) and their diverse expression patterns in response to bacterial infection

Thioester-containing protein (TEP) family members are characterized by their unique intrachain β-cysteinyl-γ-glutamyl thioesters, and they play important roles in innate immune responses. Although significant effects of TEP members on immunity have been reported in most vertebrates, as well as certain invertebrates, the complete TEP family has not been systematically characterized in scallops. In this study, five TEP family genes (PyC3, PyA2M, PyTEP1, PyTEP2 and PyCD109) were identified from Yesso scallop (Patinopecten yessoensis) through whole-genome scanning, including one pair of tandem duplications located on the same scaffold. Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of the five genes (PyTEPs). The vast distribution of PyTEPs in TEP subfamilies confirmed that the Yesso scallop contains relatively comprehensive types of TEP members in evolution. The expression profiles of PyTEPs were determined in hemocytes after bacterial infection with gram-positive (Micrococcus luteus) and gram-negative (Vibrio anguillarum) using quantitative real-time PCR (qRT-PCR). Expression analysis revealed that the PyTEP genes exhibited disparate expression patterns in response to the infection by gram bacteria. A majority of PyTEP genes were overexpressed after bacterial stimulation at most time points, especially the notable elevation displayed by duplicated genes after V. anguillarum challenge. Interestingly, at different infection times, PyTEP1 and PyTEP2 shared analogous expression patterns, as did PyC3 and PyCD109. Taken together, these results help to characterize gene duplication and the evolutionary origin of PyTEPs and supplied valuable resources for elucidating their versatile roles in bivalve innate immune responses to bacterial pathogen challenges.

Identification and Characterization of Neuropeptides by Transcriptome and Proteome Analyses in a Bivalve Mollusc Patinopecten yessoensis

Neuropeptides play essential roles in regulation of reproduction and growth in marine molluscs. But their function in marine bivalves – a group of animals of commercial importance – is largely unexplored due to the lack of systematic identification of these molecules. In this study, we sequenced and analyzed the transcriptome of nerve ganglia of Yesso scallop Patinopecten yessoensis, from which 63 neuropeptide genes were identified based on BLAST and de novo prediction approaches, and 31 were confirmed by proteomic analysis using the liquid chromatography-tandem mass spectrometry (LC-MS/MS). Fifty genes encode known neuropeptide precursors, of which 20 commonly exist in bilaterians and 30 are protostome specific. Three neuropeptides that have not yet been reported in bivalves were identified, including calcitonin/DH31, lymnokinin and pleurin. Characterization of glycoprotein hormones, insulin-like peptides, allatostatins, RFamides, and some reproduction, cardioactivity or feeding related neuropeptides reveals scallop neuropeptides have conserved molluscan neuropeptide domains, but some (e.g., GPB5, APGWamide and ELH) are characterized with bivalve-specific features. Thirteen potentially novel neuropeptides were identified, including 10 that may also exist in other protostomes, and 3 (GNamide, LRYamide, and Vamide) that may be scallop specific. In addition, we found neuropeptides potentially related to scallop shell growth and eye functioning. This study represents the first comprehensive identification of neuropeptides in scallop, and would contribute to a complete understanding on the roles of various neuropeptides in endocrine regulation in bivalve molluscs.

Genetic and molecular basis of the immune system in the brachiopod Lingula anatina

The extension of comparative immunology to non-model systems, such as mollusks and annelids, has revealed an unexpected diversity in the complement of immune receptors and effectors among evolutionary lineages. However, several lophotrochozoan phyla remain unexplored mainly due to the lack of genomic resources. The increasing accessibility of high-throughput sequencing technologies offers unique opportunities for extending genome-wide studies to non-model systems. As a result, the genome-based study of the immune system in brachiopods allows a better understanding of the alternative survival strategies developed by these immunologically neglected phyla. Here we present a detailed overview of the molecular components of the immune system identified in the genome of the brachiopod Lingula anatina. Our findings reveal conserved intracellular signaling pathways as well as unique strategies for pathogen detection and killing in brachiopods.

The various components implied the diversified Toll-like receptor (TLR) signaling pathway in mollusk Chlamys farreri

Toll-like receptor (TLR) signaling pathway, composed of various components, plays pivotal roles in host innate immune defense mechanism. In the present study, twenty-nine TLR signaling pathway components, including receptors, adaptors, transduction molecules and immune effectors, were identified in Zhikong scallop Chlamys farreri via assembling and screening public available transcriptomic data and expression sequence tags (ESTs). These identified TLR signaling pathway components were constitutively expressed and detectable in various tissues, and almost all of them were highly expressed in gill and hepatopancreas. These results indicated the presence of TLR signaling pathways in both MyD88-dependent and MyD88-independent forms in scallop, and implied the diversified TLR signaling pathway in mollusk C. farreri.

Molecular identification and functional analysis of two variants of myeloid differentiation factor 88 (MyD88) from disk abalone (Haliotis discus discus)

Myeloid differentiation factor 88 (MyD88) is a crucial adaptor protein of the Toll-like receptor (TLR)- and interleukin 1 receptor-mediated signaling pathways and is involved in a diverse array of inflammatory responses via NF-κB activation. In the present study, two MyD88 variants were identified from disk abalone (Haliotis discus discus) and designated AbMyD88-2 and AbMyD88-X. The deduced AbMyD88-2 and AbMyD88-X comprised 433 and 354 amino acids with predicted molecular masses of 48.85 kDa and 40.17 kDa, respectively. AbMyD88-2 and AbMyD88-X possessed typical MyD88 domain structural features including an N-terminal death domain (DD) and C-terminal toll interleukin 1 receptor (TIR) domain similar to those in mammals. Expression analysis of AbMyD88-2 and AbMyD88-X mRNA at different early embryonic developmental stages of abalone by qPCR revealed that their constitutive expression at all developmental stages analyzed with the considerably higher values at the 16-cell (AbMyD88-2) and morula stages (AbMyD88-X). In unchallenged disk abalones, AbMyD88-2 was highly expressed in muscles, while AbMyD88-X mRNA was predominantly transcribed in hemocytes. Moreover, AbMyD88-2 and AbMyD88-X mRNA were differentially modulated in abalone hemocytes after a challenge with live bacteria (Vibrio parahaemolyticus, Listeria monocytogenes), virus (viral hemorrhagic septicemia virus), and pathogen-associated molecular patterns (lipopolysaccharides and Poly I:C). Overexpression of AbMyD88-2 and AbMyD88-X in HEK293T cells induced the activation of the NF-κB promoter. AbMyD88-2 and AbMyD88-X involvement in inflammatory responses was characterized by their overexpression in RAW264.7 murine macrophage cells. These results revealed comparatively higher NO (Nitric oxide) production, induction of inflammatory mediator genes (iNOS and COX2), and proinflammatory genes (IL1β, IL6 and TNFα) expression in abalone MyD88s-overexpressing cells than in mock control in the presence or absence of LPS stimulation. Altogether, these results suggest that existence of a MyD88-dependent like signaling pathway in disk abalone and that both AbMyD88-2 and AbMyD88-X might be involved in innate immune and inflammatory responses.

Transcriptome analysis and discovery of genes involved in immune pathways from coelomocytes of Onchidium struma after bacterial challenge

Onchidium struma widely distributes in subtidal and low-tidal zones, which is considered to be an economical species with rich nutrition, a valuable biomonitor for heavy metal pollution and a representative species for evolution from ocean to land. However, there is limited genetic information available for O. struma development. This study compared transcriptomic profiles of coelomocytes from normal and bacteria infected O. struma by Illumina-based paired-end sequencing to explore the molecular immune mechanism of O. struma against bacterial infection. After assembly, a total of 92,450 unigenes with an average length of 1019 bp were obtained. Approximately 34,964 (37.82%) unigenes were annotated in the Nr NCBI database and 40.1% of unigenes were similar with that of Aplysia californica. Among them, 7609 unigenes were classified into three Gene Ontology (GO) categories: biological process (3250 unigenes, 42.7%), cellular component (2,281, 30.0%) and molecular function (2078 unigenes, 27.3%). A total of 22,776 unigenes were aligned to the Clusters of Orthologous Groups (COG) of proteins and classified into 25 functional categories. Following bacterial infection, 10,623 differently expressed unigenes (DEGs) were identified, including 7644 up-regulated and 2979 down-regulated unigenes. Further KEGG analysis annotated 11,681 DEGs to 42 pathways, and 11 pathways were identified to be related with diseases and immune system. To our knowledge, it was first time to analyze transcriptome profiles of O. struma. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on O. struma. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

Genome-wide identification, characterization and expression analyses of TLRs in Yesso scallop (Patinopecten yessoensis) provide insight into the disparity of responses to acidifying exposure in bivalves

Toll-like receptors (TLRs) play a crucial role in innate immunity by recognizing specific pathogen-associated molecular patterns, including lipoproteins, lipopeptides, lipopolysaccharide, flagellin, dsRNA, ssRNA and CpG DNA motifs. Although significant effects of TLRs on immunity have been reported in most vertebrates and some invertebrates, the complete TLR superfamily has not been systematically characterized in scallops. In this study, 18 TLR genes were identified from Yesso scallop (Patinopecten yessoensis) using whole-genome scanning. Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of the 18 genes. Extensive expansion of TLR genes from the Yesso scallop genome indicated gene duplication events. In addition, expression profiling of PyTLRs was performed at different acidifying exposure levels (pH = 6.50, 7.50) with different challenge durations (3, 6, 12 and 24 h) via in silico analysis using transcriptome and genome databases. Our results confirmed the inducible expression patterns of PyTLRs under acidifying exposure, and the responses to immune stress may have arisen through adaptive recruitment of tandem duplications of TLR genes. Collectively, this study provides novel insight into PyTLRs as well as the specific role and response of TLR signaling pathways in host immune responses against acidifying exposure in bivalves.

Diversity and evolution of TIR-domain-containing proteins in bivalves and Metazoa: New insights from comparative genomics

The Toll/interleukin-1 receptor (TIR) domain has a fundamental role in the innate defence response of plants, vertebrate and invertebrate animals. Mostly found in the cytosolic side of membrane-bound receptor proteins, it mediates the intracellular signalling upon pathogen recognition via heterotypic interactions. Although a number of TIR-domain-containing (TIR-DC) proteins have been characterized in vertebrates, their evolutionary relationships and functional role in protostomes are still largely unknown. Due to the high abundance and diversity of TIR-DC proteins in bivalve molluscs, we investigated this class of marine invertebrates as a case study. The analysis of the available genomic and transcriptomic data allowed the identification of over 400 full-length sequences and their classification in protein families based on sequence homology and domain organization. In addition to TLRs and MyD88 adaptors, bivalves possess a surprisingly large repertoire of intracellular TIR-DC proteins, which are conserved across a broad range of metazoan taxa. Overall, we report the expansion and diversification of TIR-DC proteins in several invertebrate lineages and the identification of many novel protein families possibly involved in both immune-related signalling and embryonic development.

Identification of a new calreticulin homolog from Yesso scallop (Patinopecten yessoensis) and its role in innate immunity

Calreticulin (CRT) is a multifunctional and highly conserved Ca²⁺-binding protein shared among vertebrates and invertebrates. In this study, we cloned and characterized a CRT gene, PyCRT, from Yesso scallop, Patinopecten yessoensis. The full-length cDNA of PyCRT was 1830 bp, including a 1242 bp open reading frame (ORF), a 29 bp 5'-untranslated region and a 559 bp 3'-untranslated region. PyCRT was consisted of three distinct structural and functional domains (N-, P- and C-domains), a signal peptide and an endoplasmic reticulum (ER) retrieval signal sequence (HDEL). Tissue specific expression analysis showed that PyCRT was distributed widely in Yesso scallop, and was highly expressed in the mantle and hemocytes. After Vibrio anguillarum challenge, the expression of PyCRT in hemocytes had a significant increase and reached the maximum level at 12 h post-infection. We also demonstrated for the first time in mollusc that the recombinant PyCRT (rPyCRT) could bind to the Gram-negative bacterium V. anguillarum, Escherichia coli and the Gram-positive bacterium Staphylococcus aureus. Our results suggested that the CRT gene from Yesso scallop possessed immune-related regulatory functions in the innate immune system in scallops.

Physical mapping of immune-related genes in Yesso scallop (Patinopecten yessoensis) using fluorescent in situ hybridization

The innate immune system plays a pivotal role in defending invasion of microorganisms for scallops. Previous studies on immune-related genes in the Yesso scallop, Patinopecten yessoensis (Jay, 1857) have mainly focused on characterization and expression pattern in response to bacterial challenge, no research has been carried out on the cytogenetic level yet. In the present study, eight fosmid clones containing the sequences of key immune-related genes (PyNFkB, PyTRAF2, PyTRAF4, PyTRAF7, PyMyd88-1, PyMyd88-3, PyMKK-7 and PyTNFR) were isolated and seven of them were successfully mapped on chromosomes of Patinopecten yessoensis utilizing fluorescence in situ hybridization. Wherein, PyMyd88-1, PyMyd88-3 and PyMKK-7 located on the same chromosome pair with adjacent positions and the other genes were mapped on four non-homologous chromosome pairs, showing a similar distribution to another five model species. The isolation and mapping of such genes of the Yesso scallop will lay a foundation for studies such as assignment of interested genes to chromosomes, construction cytogenetic maps and so on.

Effects of dietary oxidized konjac glucomannan sulfates (OKGMS) and acidolysis-oxidized konjac glucomannan (A-OKGM) on the immunity and expression of immune-related genes of Schizothorax prenanti

In the present study, konjac glucomannan (KGM) was degraded by H2O2, and then used trisulfonated sodium amine and HCl, individually, to obtain two kinds of derivatives: oxidized konjac glucomannan sulfates (OKGMS) and acidolysis-oxidized konjac glucomannan (A-OKGM). The effects of two OKGM modified products on the immune parameters and expressions of toll-like receptor 22 (TLR22), myeloid differentiation factor 88 (MyD88) and interferon regulatory factors 7 (IRF7) genes in Schizothorax prenanti were determined. The alternative haemolytic complement (ACH50) activity was found to be significantly increased by the OKGMS diets. The immunoglobulin M (IgM) level was significantly enhanced by the OKGMS diets. The lysozyme activity was significantly increased by both OKGMS and A-OKGM diets. The superoxide dismutase (T-SOD) activity in fish fed with all doses of OKGMS diets was significantly higher than that in fish fed with basal diet. The glutathione peroxidase (GSH-PX) activity in fish fed with 0.8% and 1.6% A-OKGM diets was significantly higher than control group. The malondialdehyde (MDA) level was significantly decreased by both OKGMS and A-OKGM diets. The 0.8% A-OKGM diet significantly up-regulated TLR22 gene expression in the head kidney and spleen. TLR22 gene expression was significantly promoted by all OKGMS diets in the mesonephros and liver. The MyD88 mRNA level in 1.6% A-OKGM group significantly increased in the head kidney. The low dose of OKGMS significantly induced the MyD88 gene expression in the mesonephros, gut and liver, while 0.8% A-OKGM group also showed a significantly enhanced MyD88 mRNA expression in the gut. High dose of OKGMS significantly increased the IRF7 mRNA expression in the mesonephros and spleen. Fish fed with low dose of A-OKGM showed significantly higher expression of IRF7 in the gut and liver. Present study suggested that OKGMS and A-OKGM can act as immunostimulant to improve the immune indexes and up-regulate the immune-related gene expressions.

Characterization of three mitogen-activated protein kinases (MAPK) genes reveals involvement of ERK and JNK, not p38 in defense against bacterial infection in Yesso scallop Patinopecten yessoensis

Mitogen-activated protein kinases (MAPKs) are protein Ser/Thr kinases that play a vital role in innate immune responses by converting extracellular stimuli into a wide range of cellular responses. Although MAPKs have been extensively studied in various vertebrates and invertebrates, our current understanding of MAPK signaling cascade in scallop is in its infancy. In this study, three MAPK genes (PyERK, PyJNK, and Pyp38) were identified from Yesso scallop Patinopecten yessoensis. The open reading frame of PyERK, PyJNK, and Pyp38 was 1104, 1227, and 1104 bp, encoding 367, 408, and 367 amino acids, respectively. Conservation in some splicing sites was revealed across the three PyMAPKs, suggesting the common descent of MAPKs genes. The expression profiles of PyMAPKs over the course of ten different developmental stages showed that they had different expression patterns. In adult scallops, PyMAPKs were primarily expressed in muscles, hemocytes, gill, and mantle. To gain insights into their role in innate immunity, we investigated their expression profiles after infection with Gram-positive bacteria (Micrococcus luteus) and Gram-negative bacteria (Vibrio anguillarum). Significant difference in gene expression was only found in PyERK and PyJNK, but not Pyp38, suggesting Pyp38 may not participate in immune response to bacterial infection. Besides, PyERK and PyJNK exhibited more drastic change against the invasion of V. anguillarum than M. luteus, suggesting they could be more sensitive to Gram-negative bacteria than Gram-positive bacteria. This study provides valuable resource for elucidating the role of MAPK signal pathway in bivalve innate immune response.

The categorization and mutual modulation of expanded MyD88s in Crassostrea gigas

MyD88 serves as a critical cytosolic adaptor mediating activation of NF-κB in innate immunity. It has been found that there is a considerable expansion of MyD88 in Crassostrea gigas. In the present study, four typical MyD88 genes in Crassostrea gigas (CgMyD88-A to CgMyD88-D) were successfully cloned and their potential functions were investigated together with another two known ones (CgMyD88-T1 and CgMyD88-T2). Multiple alignments revealed that CgMyD88-B and CgMyD88-C remained the conserved DD and TIR domains, while there was a significant variation of E51Q in the DD of CgMyD88-A, and some variations in both DD and TIR domains of CgMyD88-D, respectively. Both truncated CgMyD88-T1 and CgMyD88-T2 lacked Box II in their only TIR domains. Expression pattern analysis showed that CgMyD88-B and CgMyD88-C genes possessed higher expression in normal tissues, compared with the other four. When oysters were under bacteria challenge, CgMyD88-B, CgMyD88-C, CgMyD88-T1 and CgMyD88-T2 were firstly induced, while CgMyD88-A and CgMyD88-D were suppressed. Dual luciferase reporter assays showed that CgMyD88-B and CgMyD88-C could promote the activation of NF-κB signaling pathway, while the other four CgMyD88 genes failed or even suppressed the activities of CgMyD88-B and CgMyD88-C on the activation of NF-κB signaling. It was deduced that after oysters were challenged by bacteria, CgMyD88-B and CgMyD88-C could rapidly and efficiently activate NF-κB signaling pathway to elicit anti-pathogen responses before suppressor CgMyD88 genes (CgMyD88-T1 and CgMyD88-T2) exceeding their expression level. These results suggested that there was mutual modulation of expanded CgMyD88 genes on activating NF-κB signaling pathway in oyster C. gigas.

Genome-wide identification, characterization and expression analyses of two TNFRs in Yesso scallop (Patinopecten yessoensis) provide insight into the disparity of responses to bacterial infections and heat stress in bivalves

Tumor necrosis factors receptors (TNFRs) comprise a superfamily of proteins characterized by a unique cysteine-rich domain (CRD) and play important roles in diverse physiological and pathological processes in the innate immune system, including inflammation, apoptosis, autoimmunity and organogenesis. Although significant effects of TNFRs on immunity have been reported in most vertebrates as well as some invertebrates, the complete TNFR superfamily has not been systematically characterized in scallops. In this study, two different types of TNFR-like genes, including PyTNFR1 and PyTNFR2 genes were identified from Yesso scallop (Patinopecten yessoensis, Jay, 1857) through whole-genome scanning. Phylogenetic and protein structural analyses were carried out to determine the identities and evolutionary relationships of the two genes. The expression profiling of PyTNFRs was performed at different development stages, in healthy adult tissues and in hemocytes after bacterial infection and heat stress. Expression analysis revealed that both PyTNFRs were significantly induced during the acute phase (3 h) after infection with Gram-positive (Micrococcus luteus) and Gram-negative (Vibrio anguillarum) bacteria, though much more dramatic chronic-phase (24 h) changes were observed after V. anguillarum challenge. For heat stress, only PyTNFR2 displayed significant elevation at 12 h and 24 h, which suggests a functional difference in the two PyTNFRs. Collectively, this study provides novel insight into the PyTNFRs and the specific role and response of TNFR-involved pathways in host immune responses against different bacterial pathogens and heat stress in bivalves.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

Evaluation of the 2b-RAD method for genomic selection in scallop breeding

The recently developed 2b-restriction site-associated DNA (2b-RAD) sequencing method provides a cost-effective and flexible genotyping platform for aquaculture species lacking sufficient genomic resources. Here, we evaluated the performance of this method in the genomic selection (GS) of Yesso scallop (Patinopecten yessoensis) through simulation and real data analyses using six statistical models. Our simulation analysis revealed that the prediction accuracies obtained using the 2b-RAD markers were slightly lower than those obtained using all polymorphic loci in the genome. Furthermore, a small subset of markers obtained from a reduced tag representation (RTR) library presented comparable performance to that obtained using all markers, making RTR be an attractive approach for GS purpose. Six GS models exhibited variable performance in prediction accuracy depending on the scenarios (e.g., heritability, sample size, population structure), but Bayes-alphabet and BLUP-based models generally outperformed other models. Finally, we performed the evaluation using an empirical dataset composed of 349 Yesso scallops that were derived from five families. The prediction accuracy for this empirical dataset could reach 0.4 based on optimal GS models. In summary, the genotyping flexibility and cost-effectiveness make 2b-RAD be an ideal genotyping platform for genomic selection in aquaculture breeding programs.

Genome-wide identification and characterization of TRAF genes in the Yesso scallop (Patinopecten yessoensis) and their distinct expression patterns in response to bacterial challenge

The tumor necrosis factor (TNF) receptor associated factors (TRAFs) are the major signal transducers for the TNF receptor superfamily and the interleukin-1 receptor/Toll-like receptor (IL-1R/TLR) superfamily, which regulate a variety of cellular activities and innate immune responses. TRAF genes have been extensively studied in various species, including vertebrates and invertebrates. However, as one of the key component of NF-κB pathway, TRAF genes have not been systematically characterized in marine invertebrates. In this study, we identified and characterized five TRAF genes, PyTRAF2, PyTRAF3, PyTRAF4, PyTRAF6 and PyTRAF7, in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine their identities and evolutionary relationships. In comparison with the TRAF genes from vertebrate species, the structural features were all relatively conserved in the PyTRAF genes. To gain insights into the roles of TRAF genes during scallop innate immune responses, quantitative real-time PCR was used to investigate the expression profiles in the different stages of scallop development, in the healthy adult tissues, and in the hemocytes after bacterial infection with Micrococcus luteus and Vibrio anguillarum. Based on the qRT-PCR analysis, the expression of most of the PyTRAFs was significantly induced in the acute phases (3-6 h) after infection with Gram-positive (M. luteus) and Gram-negative (V. anguillarum) bacteria, and many more dramatic changes in PyTRAFs expression were observed after V. anguillarum challenge. Notably, the strong response in the up-regulation of PyTRAF6 post-bacterial challenge was distinct from that previously reported in scallops and crabs but was similar to that of other shellfish, Echinodermata and even teleost fish. The high level expressions of PyTRAFs in the hemocytes and the gill, and their specific expression patterns after challenges provide insights into the versatile roles and responses of TRAFs in the innate immune system against Gram-negative bacterial pathogens in bivalves.