Identification of reference genes for qRT-PCR analysis in Yesso scallop Patinopecten yessoensis

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora

BACKGROUND

Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation.

RESULTS

In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR.

CONCLUSIONS

This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

Gene delivery available in molluscan cells by strong promoter discovered from bivalve-infectious virus

Understanding gene functions in marine invertebrates has been limited, largely due to the lack of suitable assay systems. Such a system requires investigative methods that are reproducible and can be quantitatively evaluated, such as a cell line, and a strong promoter that can drive high expression of a transgene. In this study, we established primary cell culture from a marine bivalve mollusc, Mizuhopecten yessoensis. Using scallop primary cells, we optimized electroporation conditions for transfection and carried out a luciferase-based promoter activity assay to identify strong promoter sequences that can drive expression of a gene of interest. We evaluated potential promoter sequences from genes of endogenous and exogenous origin and discovered a strong viral promoter derived from a bivalve-infectious virus, ostreid herpesvirus-1 (OsHV-1). This promoter, we termed OsHV-1 promoter, showed 24.7-fold and 16.1-fold higher activity than the cytomegalovirus immediate early (CMV IE) promoter and the endogenous EF1α promoter, the two most commonly used promoters in bivalves so far. Our GFP assays showed that the OsHV-1 promoter is active not only in scallop cells but also in HEK293 cells and zebrafish embryos. The OsHV-1 promoter practically enables functional analysis of marine molluscan genes, which can contribute to unveiling gene-regulatory networks underlying astonishing regeneration, adaptation, reproduction, and aging in marine invertebrates.

Quantitative proteomic analysis reveals the molecular mechanism of the Yesso scallop (Patinopecten yessoensis) in response to Polydora infection

The Yesso scallop is a large and ancient molluscan group with great economic value; however, it has recently suffered severe cases of Polydora infection. Polydora parasitizes the shells of scallops, badly damaging shell structures and affecting growth and mortality. To investigate the molecular mechanism of Yesso scallops’ response to Polydora infection, proteomic profiling changes in the mantle tissues of Polydora-infected (diseased) and healthy scallops were systematically analysed by tandem mass tags (TMT) labelling technology in this study. A total of 519 differentially expressed proteins (DEPs) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed most innated immune-related functions and pathways were significantly downregulated in diseased scallops, except the phagocytosis pathway, indicating an important role of phagocytosis in response to Polydora infection. DEPs involved in the phagocytosis pathway were associated with phagocytic receptor recognition, phagosome biogenesis and pathogen degradation, and they were further verified by quantitative real-time PCR. The results elucidate the molecular components of phagocytosis in molluscs for the first time. Polydora can be encapsulated by melanization with an obvious appearance in shells; indeed, melanization-related DEPs were upregulated in diseased scallops. Inhibition of apoptosis and nervous modulation may be also involved in the response mechanism, with some highly associated proteins significantly differentially expressed. Finally, a protein–protein interaction network was constructed to provide a global view of the interaction relationships of the DEPs. The study predicts the molecular response mechanism of Yesso scallops to Polydora infection, and lays a theoretical foundation for Polydora disease control.

Gonadal somatic cell-specific transforming growth factor-β superfamily member in the Yesso scallop reveals gonadal somatic cell distribution during the reproductive phase

The objective of this study was to identify the gonadal somatic cells in the Yesso scallop using a novel molecular marker. This study is the first to identify the bone morphogenetic protein 2a (Bmp2a) gene as a gonadal somatic cell-specific gene in this bivalve. We performed a transcriptomic survey to identify the transforming growth factor-β (TGFβ) superfamily members that act in Yesso scallop gonad development. BLAST survey, phylogenetic tree, and RT-PCR analyses screened BMP molecules (i.e., bmp2a and bmp10a), which are members of the TGFβ superfamily that show gonad-specific expression. Among the BMPs from the Yesso scallop, in situ hybridization accompanied by RNAscope assay identified that bmp2a mRNA was specifically expressed in the gonadal somatic cells localized in the interspace between germ cells. Real-time quantitative PCR (qPCR) analysis revealed that bmp2a mRNA expression increased during the reproductive phase. The relative expression of bmp2a mRNA was lowest at the beginning of the growing stage and peaked at the mature stage in both sexes. These observations indicate that bmp2a-positive gonadal somatic cells support germ cell growth and differentiation during the reproductive phase for both sexes. This study provides new insights into gonadal somatic cell biology in marine invertebrates and we propose that TGFβ signaling is necessary for gonad development in bivalves.

Identifying conserved polychaete molecular markers of metal exposure: Comparative analyses using the Alitta virens (Annelida, Lophotrochozoa) transcriptome

Polychaetes are vital for evaluating the effects of toxic metals in marine systems, and sensitive molecular biomarkers should be integral to monitoring efforts. However, the few polychaete markers that exist are inconsistent, even within the same species, failing to identify gene expression changes in metal-exposed animals incurring clear metabolic costs. Comparing previously characterised polychaete metal-responsive genes with those of another carefully selected species could identify biomarkers applicable across polychaetes. The ragworm Alitta virens (Sars, 1835) is particularly suited for such comparisons due to its dominance of fully saline coastal areas, widespread distribution, large biomass, and its phylogenetic position relative to other polychaete 'omic' resources. A transcriptome atlas for A. virens was generated and an RNASeq-qPCR screening approach was used to characterise the response to chronic exposures of environmentally relevant concentrations of copper and zinc in controlled mesocosms. Genes presenting dramatic expression changes in A. virens were compared with known metal-responsive genes in other polychaetes to identify new possible biomarkers and assess those currently used. This revealed some current markers should probably be abandoned (e.g. Atox1), while others, such as GST-Omega, should be used with caution, as different polychaete species appear to upregulate distinct GST-Omega orthologues. In addition, the comparisons give some indication of genes that are induced by metal exposure across phylogenetically divergent polychaetes, including a suite of haemoglobin subunits and linker chains that could play conserved roles in metal-stress response. Although such newly identified markers need further characterisation, they offer alternatives to current markers that are plainly insufficient.

Screening of reference genes using real-time quantitative PCR for gene expression studies in Neoseiulus barkeri Hughes (Acari: Phytoseiidae)

A stable reference gene is a key prerequisite for accurate assessment of gene expression. At present, the real-time reverse transcriptase quantitative polymerase chain reaction has been widely used in the analysis of gene expression in a variety of organisms. Neoseiulus barkeri Hughes (Acari: Phytoseiidae) is a major predator of mites on many important economically crops. Until now, however, there are no reports evaluating the stability of reference genes in this species. In view of this, we used GeNorm, NormFinder, BestKeeper, and RefFinder software tools to evaluate the expression stability of 11 candidate reference genes in developmental stages and under various abiotic stresses. According to our results, β-ACT and Hsp40 were the top two stable reference genes in developmental stages. The Hsp60 and Hsp90 were the most stable reference genes in various acaricides stress. For alterations in temperature, Hsp40 and α-TUB were the most suitable reference genes. About UV stress, EF1α and α-TUB were the best choice, and for the different prey stress, β-ACT and α-TUB were best suited. In normal conditions, the β-ACT and α-TUB were the two of the highest stable reference genes to respond to all kinds of stresses. The current study provided a valuable foundation for the further analysis of gene expression in N. barkeri.

Phenotypic Stability of Sex and Expression of Sex Identification Markers in the Adult Yesso Scallop Mizuhopecten yessoensis throughout the Reproductive Cycle

Simple Summary

Bivalve sex is thought to fluctuate depending on environmental conditions. So far, there has been no investigation on the phenotypic stability of sex in the commercially important Yesso scallop Mizuhopecten yessoensis. The present study revealed that the sex of the Yesso scallop is stable after initial sex differentiation and that this species maintains a sex-stable maturation system throughout its life. In addition, gonad differentiation for each sex was precisely characterized by using molecular markers throughout the maturational cycle.

Abstract

The objective of the present study was to analyze the phenotypic stability of sex after sex differentiation in the Yesso scallop, which is a gonochoristic species that has been described as protandrous. So far, no study has investigated in detail the sexual fate of the scallop after completion of sex differentiation, although bivalve species often show annual sex change. In the present study, we performed a tracking experiment to analyze the phenotypic stability of sex in scallops between one and two years of age. We also conducted molecular marker analyses to describe sex differentiation and gonad development. The results of the tracking experiment revealed that all scallops maintained their initial sex phenotype, as identified in the last reproductive period. Using molecular analyses, we characterized my-dmrt2 and my-foxl2 as sex identification markers for the testis and ovary, respectively. We conclude by proposing that the Yesso scallop is a sex-stable bivalve after its initial sex differentiation and that it maintains a sex-stable maturation system throughout its life. The sex-specific molecular markers identified in this study are useful tools to assess the reproductive status of the Yesso scallop.

Systematic identification and validation of the reference genes from 60 RNA-Seq libraries in the scallop Mizuhopecten yessoensis

Background

Reverse transcription quantitative PCR (RT-qPCR) is widely used for gene expression analysis in various organisms. Its accuracy largely relies on the stability of reference genes, making reference gene selection a vital step in RT-qPCR experiments. However, previous studies in mollusks only focused on the reference genes widely used in vertebrates.

Results

In this study, we conducted the transcriptome-wide identification of reference genes in the bivalve mollusk Mizuhopecten yessoensis based on 60 transcriptomes covering early development, adult tissues and gonadal development. A total of 964, 1210 and 2097 candidate reference genes were identified, respectively, resulting in a core set of 568 genes. Functional enrichment analysis showed that these genes are significantly overrepresented in Gene Ontology (GO) terms or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to ribosomes, energy production, etc. Six genes (RS23, EF1A, NDUS4, SELR1, EIF3F, and OLA1) were selected from the candidate genes for RT-qPCR validation, together with 6 commonly used reference genes (ACT, CYTC, HEL, EF1B, GAPDH and RPL16). Stability analyses using geNorm, NormFinder and the comparative delta-Ct method revealed that the new candidate reference genes are more stable than the traditionally used genes, and ACT and CYTC are not recommended under either of the three circumstances. There was a significant correlation between the Ct of RT-qPCR and the log₂(TPM) of RNA-Seq data (Ct = − 0.94 log₂(TPM) + 29.67, R² = 0.73), making it easy to estimate the Ct values from transcriptome data prior to RT-qPCR experiments.

Conclusion

Our study represents the first transcriptome-wide identification of reference genes for early development, adult tissues, and gonadal development in the Yesso scallop and will benefit gene expression studies in other bivalve mollusks.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5661-x) contains supplementary material, which is available to authorized users.

A carotenoid oxygenase is responsible for muscle coloration in scallop

As lipid microconstituents mainly of plant origin, carotenoids are essential nutrients for humans and animals, and carotenoid coloration represents an important meat quality parameter for many farmed animals. Currently, the mechanism of carotenoid bioavailability in animals is largely unknown mainly due to the limited approaches applied, the shortage of suitable model systems and the restricted taxonomic focus. The mollusk Yesso scallop (Patinopecten yessoensis) possessing orange adductor muscle with carotenoid deposition, provides a unique opportunity to research the mechanism underlying carotenoid utilization in animals. Herein, through family construction and analysis, we found that carotenoid coloration in scallop muscle is inherited as a recessive Mendelian trait. Using a combination of genomic approaches, we mapped this trait onto chromosome 8, where PyBCO-like 1 encoding carotenoid oxygenase was the only differentially expressed gene between the white and orange muscles (FDR = 2.75E-21), with 11.28-fold downregulation in the orange muscle. Further functional assays showed that PyBCO-like 1 is capable of degrading β-carotene, and inhibiting PyBCO-like 1 expression in the white muscle resulted in muscle coloration and carotenoid deposition. In the hepatopancreas, which is the organ for digestion and absorption, neither the scallop carotenoid concentration nor PyBCO-like 1 expression were significantly different between the two scallops. These results indicate that carotenoids could be taken up in both white- and orange-muscle scallops and then degraded by PyBCO-like 1 in the white muscle. Our data suggest that PyBCO-like 1 is the essential gene for carotenoid metabolism in scallop muscle, and its downregulation leads to carotenoid deposition and muscle coloration.

Molecular identification of steroidogenesis-related genes in scallops and their potential roles in gametogenesis

Sex steroids are crucial for controlling gametogenesis and germ cell maturation in vertebrates. It has been proposed that Yesso scallop (Mizuhopecten yessoensis) has the same sex steroids as those animals, but the scallop biosynthetic pathway is unclear. In this study, we characterized several steroidogenesis-related genes in M. yessoensis and proposed a putative biosynthetic pathway for sex steroids that is similar to that of vertebrates. Specifically, we identified several steroidogenesis-related gene sequences that encode steroid metabolizing enzymes: StAR-related lipid transfer (START) protein, 17α-hydroxylase, 17,20-lyase (cyp17a), 17β-hydroxysteroid dehydrogenase (hsd17b), and 3β-hydroxysteroid dehydrogenase (hsd3b). We sampled adult scallops throughout their reproductive phase to compare their degree of maturation with their intensity of mRNA expression. Semi-quantitative RT-PCR analysis revealed a ubiquitous expression of transcripts for steroid metabolizing enzymes (i.e., star, cyp17a, hsd17b, and hsd3b) in peripheral and gonadal tissues. Real-time PCR analysis revealed a high level of expression of star3 and cyp17a genes in gonadal tissues at the early stage of cell differentiation in scallops. Interestingly, mRNA expression of hsd3b and hsd17b genes showed a synchronous pattern related to degree of gonad maturity. These results indicate that both hsd3b and hsd17b genes are likely involved in steroidogenesis in scallops. We therefore believe that these steroid-metabolizing enzymes allow scallops to endogenously produce sex steroids to regulate reproductive events.

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.

The scallop IGF2 mRNA-binding protein gene PyIMP and association of a synonymous mutation with growth traits

Insulin-like growth factor 2 mRNA-binding proteins (IMPs) function in localization, stability and translational control of their target RNAs. In this study, we identified an IMP gene (PyIMP) from Yesso scallop, Patinopecten yessoensis. The complete DNA sequence of PyIMP was 22,875 bp, consisting of seventeen exons and sixteen introns. The full-length cDNA sequence was 3,293 bp, with an open reading frame of 1,776 bp, encoding 592 amino acids. PyIMP exhibited characters typical of IMPs, namely two RNA recognition motifs and four hnRNP K homology domains. Real-time quantitative reverse transcription PCR analysis indicated that PyIMP was universally expressed, with higher expression levels in the gonad of adult scallops, and in gastrulae and trochophore larvae at developmental stages. A synonymous mutation SNP, c.852A>G, which showed significant associations with growth traits of Yesso scallop, was identified in this gene. Scallops with the AA genotype at this locus had significantly higher trait values than those with the GG genotype for shell length, shell height, body weight, soft tissue weight and striated muscle weight (P < 0.05). Meanwhile, the expression of PyIMP in AA type scallops was significantly higher than that in the GG type, implying a positive effect of PyIMP on scallop growth. PyIMP represents the first mRNA-binding protein gene characterized in mollusks, and SNP c.852A>G will be useful for a better understanding of the role of mRNA-binding proteins in bivalves and for scallop breeding.

Gene expression analysis of Ruditapes philippinarum haemocytes after experimental Perkinsus olseni zoospore challenge and infection in the wild

The production of Manila clam (Ruditapes philippinarum) is seriously threatened by the protistan parasite Perkinsus olseni. We characterized and compared gene expression of Manila clam haemocytes in response to P. olseni in a time-course (10 h, 24 h, 8 d) controlled laboratory challenge (LC), representing the first step of infection, and in a more complex infection in the wild (WI), using a validated oligo-microarray containing 11,232 transcripts, mostly annotated. Several immune-genes involved in NIK/NF-kappaB signalling, Toll-like receptor signalling and apoptosis were activated at LC-10 h. However, down-regulation of genes encoding lysozyme, histones, cathepsins and heat shock proteins indicated signals of immunodepression, which persisted at LC-24 h, when only down-regulated genes were detected. A rebound of haemocyte activity occurred at LC-8 d as shown by up-regulation of genes involved in cytoskeleton organization and cell survival. The WI study showed a more complex picture, and several immune-relevant processes including cytoskeleton organization, cell survival, apoptosis, encapsulation, cell redox- and lipid-homeostasis were activated, illustrating the main mechanism of host response. Our results provide useful information, including potential biomarkers, to develop strategies for controlling Manila clam perkinsosis.

Characterization of the TRAF3IP1 gene in Yesso scallop (Patinopecten yessoensis) and its expression in response to bacterial challenge

Tumor necrosis factor receptor-associated factor 3 (TRAF3) is an important adaptor that transmits upstream activation signals to induce innate immune responses. TRAF3 interacting protein 1 (TRAF3IP1) interacts specifically with TRAF3, but its function in innate immunity remains unclear, especially in marine invertebrates. In this study, to better understand the functions of TRAFs in innate immune responses, we identified and characterized the first bivalve TRAF3IP1 gene, PyTRAF3IP1, from Yesso scallop (Patinopecten yessoensis), one of the most important mollusk species for aquaculture. The PyTRAF3IP1 cDNA is 2,367 bp, with an open reading frame of 1,629 bp encoding 542 amino acids. Phylogenetic and protein structural analysis confirmed the gene's identity and revealed that PyTRAF3IP1 was more similar to vertebrate TRAF3IP1s than to those of invertebrates. PyTRAF3IP1 was expressed in all the adult tissues and developmental stages sampled, implying that it plays versatile roles in many biological processes. Furthermore, PyTRAF3IP1 expression was dramatically induced in the acute phase (3-6 h) after infection with both Gram-positive (Micrococcus luteus) and Gram-negative (Vibrio anguillarum) bacteria, even stronger induction being observed after V. anguillarum challenge. This is the first report of the characterization and immune response involvement of TRAF3IP1 in marine invertebrates, and suggests that TRAF3IP1 contributes to innate immunity in bivalves.

Validation of reference genes for RT-qPCR in marine bivalve ecotoxicology: Systematic review and case study using copper treated primary Ruditapes philippinarum hemocytes

The appropriate selection of reference genes for the normalization of non-biological variance in reverse transcription real-time quantitative PCR (RT-qPCR) is essential for the accurate interpretation of the collected data. The use of multiple validated reference genes has been shown to substantially increase the robustness of the normalization. It is therefore considered good practice to validate putative genes under specific conditions, determine the optimal number of genes to be employed, and report the method or methods used. Under this premise, we assessed the current state of reference gene based normalization in RT-qPCR bivalve ecotoxicology studies (post 2011), employing a systematic quantitative literature review. A total of 52 papers met our criteria and were analysed for genes used, the use of multiple reference genes, as well as the validation method employed. We further critically discuss methods for reference gene validation based on a case study using copper exposed primary hemocytes from the marine bivalve Ruditapes philippinarum; including the established algorithms geNorm, NormFinder and BestKeeper, as well as the popular online tool RefFinder. We identified that RT-qPCR normalization is largely performed using single reference genes, while less than 40% of the studies attempted to experimentally validate the expression stability of the genes used. 18s rRNA and β-Actin were the most popular genes, yet their un-validated use did introduce artefactual variance that altered the interpretation of the resulting data. Our findings further suggest that combining the results from multiple individual algorithms and calculating the overall best-ranked gene, as computed by the RefFinder tool, does not by default lead to the identification of the most suitable reference genes.

Assessment of housekeeping genes as internal references in quantitative expression analysis during early development of oyster

The early development of mollusks exhibits important characteristics from the developmental and evolutionary perspective. With the increasing number of genome-wide studies, accurate analyses of quantitative gene expression during development are impeded by the lack of validated reference genes. To improve the situation, in this study, we analyzed the expression stability of seven candidate housekeeping genes during early development of the Pacific oyster Crassostrea gigas: actin, glyceraldehyde-3-phosphate dehydrogenase (gapdh), α subunit of elongation factor 1 (elf1α), adp-ribosylation factor 1 (arf1), heterogeneous nuclear ribonucleoprotein q, ubiquitin-conjugating enzyme e2d2 and ribosomal protein s18. We focused on 11 stages from oocyte to D-veliger, which include crucial developmental processes such as axis determination, gastrulation and shell formation. Gene expression stabilities were assessed with the three commonly used programs geNorm, NormFinder and BestKeeper. Although the results obtained with the three programs varied to some extent, in general, arf1, elf1α and gapdh were highly ranked and actin was poorly ranked. This analysis also indicated that multiple genes should be used for normalization, and we concluded that arf1-elf1α-gapdh should be used as internal references. The findings of this study will help researchers to obtain accurate results in future quantitative gene expression analysis of development in bivalve mollusks.

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.

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.

Evaluation of housekeeping genes as references for quantitative real-time PCR analysis of gene expression in the murrel Channa striatus under high-temperature stress

Quantitative real-time polymerase chain reaction is the most advanced method of quantifying gene expression studies; however, the significance of the obtained results strongly depends on the normalization of the data to compensate for differences between the samples. In the present study, expression analysis of six different constitutively expressed genes viz. 18S ribosomal RNA, glyceraldehyde-3-phosphate dehydrogenase (gapdh), beta actin (βactin), ribosomal binding protein L13, tubulin and TATA-box-binding protein (tbp) were carried out to test their efficacy as reference genes in three different tissues, namely liver, gill and muscle of murrel Channa striatus exposed to high temperature for variable time periods. The stability and suitability of the genes were determined by using bioinformatic tools: GeNorm, NormFinder and BestKeeper. Based on the results, tub/βactin could be used as the reference genes for liver and gill tissues and βactin/gapdh could be the reference genes for muscle tissues in Channa striatus under both short- and long-term thermal stress.

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.

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

Myeloid differentiation factor 88 (MyD88) is a pivotal adaptor in the TLR/IL-1R signaling pathway, which plays an important role in activating the innate immune system. Although MyD88 genes have been identified in a variety of species, they have not been systematically characterized in scallops. In this study, five MyD88 genes were identified in Yesso scallop (Patinopecten yessoensis), PyMyD88-1, PyMyD88-2a, PyMyD88-2b, PyMyD88-3 and PyMyD88-4, which consisted of two pairs of tandem duplications located on the same chromosome. To our knowledge, this is the largest number of MyD88 genes found in an invertebrate. Phylogenetic and protein structural analyses were carried out to determine the identities and evolutionary relationships of these genes. PyMyD88s have highly conserved structures compared to MyD88 genes from other invertebrate species, except for PyMyD88-4, which contains only a DD domain, suggesting the evolutionarily conserved form of this particular gene member. We investigated the expression profiles of PyMyD88 genes at different developmental stages and in healthy adult tissues and hemocytes after Micrococcus luteus and Vibrio anguillarum infection using quantitative real-time PCR (qRT-PCR). The expression of most PyMyD88s was significantly induced in the acute phase (3-6 h) after infection with both gram-positive (M. luteus) and gram-negative (V. anguillarum) bacteria, with much more dramatic changes in PyMyD88 expression being observed after V. anguillarum challenge. Collectively, the abundance of MyD88s and their specific expression patterns provide insight into their versatile roles in the response of the bivalve innate immune system to gram-negative bacterial pathogens.

Identification and Characterization of Reference Genes for Normalizing Expression Data from Red Swamp Crawfish Procambarus clarkii

qRT-PCR is a widely used technique for rapid and accurate quantification of gene expression data. The use of reference genes for normalization of the expression levels is crucial for accuracy. Several studies have shown that there is no perfect reference gene that is appropriate for use in all experimental conditions, and research on suitable reference genes in red swamp crawfish (Procambarus clarkii) is particularly scarce. In this study, eight commonly used crustacean reference genes were chosen from P. clarkii transcriptome data and investigated as potential candidates for normalization of qRT-PCR data. Expression of these genes under different experimental conditions was examined by qRT-PCR, and the stability of their expression was evaluated using three commonly used statistical algorithms, geNorm, NormFinder and BestKeeper. A final comprehensive ranking determined that EIF and 18S were the optimal reference genes for expression data from different tissues, while TBP and EIF were optimal for expression data from different ovarian developmental stages. To our knowledge, this is the first systematic analysis of reference genes for normalization of qRT-PCR data in P. clarkii. These results will facilitate more accurate and reliable expression studies of this and other crustacean species.

The genome-wide identification of mitogen-activated protein kinase kinase (MKK) genes in Yesso scallop Patinopecten yessoensis and their expression responses to bacteria challenges

Mitogen-activated protein kinase kinases (MKK) are the essential components of the evolutionarily conserved MAPK signaling cascade, which regulates a variety of cellular activities and innate immune responses. Although MKK genes have been extensively studied in various vertebrate and invertebrate species, they have not been systematically characterized in bivalves. In this study, we identified and characterized five MKK genes (PyMKK1/2, PyMKK4, PyMKK5, PyMKK3/6 and PyMKK7) in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine their identities and evolutionary relationships. To gain insights into the possible roles of MKK genes during scallop innate immune responses, quantitative real-time PCR (qRT-PCR) was used to investigate their expression profiles during different developmental stages in samples taken from healthy adult tissues and hemocytes after Micrococcus luteus and Vibrio anguillarum bacterial infections. The Yesso scallop MKKs (PyMKKs) were found to have highly conserved structural features compared to the MKK genes from other invertebrate species. Using qRT-PCR analysis, three distinct expression patterns were detected among the PyMKKs over the course of ten different developmental stages. In adult scallops, the majority of the PyMKKs were highly expressed in mantle, gill, muscle and hemocytes. The differential expression patterns of the five PyMKKs after M. luteus (Gram-positive) and V. anguillarum (Gram-negative) bacterial infections suggested their possible involvement in the innate immune response and provide the foundation and resource for the further study on innate immune response of MAPK signal pathway in mollusk.

Characterizations and expression analyses of NF-κB and Rel genes in the Yesso scallop (Patinopecten yessoensis) suggest specific response patterns against Gram-negative infection in bivalves

Rel/NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) genes are evolutionarily conserved and play a pivotal role in several physiological events. They have been extensively studied from various species, including both vertebrates and invertebrates. However, the Rel/NF-κB genes have not been systematically characterized in bivalves. In this study, we identified and characterized PyNF-κB and PyRel in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine the identities and evolutionary relationships of Rel/NF-κB genes in Yesso scallop. Compared with the Rel/NF-κB genes from vertebrate species, the PyNF-κB and PyRel are relatively conserved in their structural features, but there were no paralogs found in P. yessoensis or other invertebrates. To gain insights into the roles of Rel/NF-κB genes during the innate immune response in scallop, quantitative real-time PCR was used to investigate the expression profiles of these genes at different developmental stages, in healthy adult tissues and in the hemolymph after bacterial infection with Micrococcus luteus and Vibrio anguillarum. The real-time PCR results indicated the abundance of PyNF-κB in the first four embryonic stages, including oocytes, fertilized eggs, morulae and blastulae. By contrast, PyRel was abundantly expressed in blastulae, trochophores and D-shaped larvae. In adult scallops, PyNF-κB and PyRel were ubiquitously expressed in most healthy tissues and highly expressed in most of the immune related tissues. Both genes were significantly up-regulated during the acute phase (3 h) after infection with Gram-positive (M. luteus) and negative (V. anguillarum) bacteria, while the much higher expression level of PyNF-κB suggested the involvement of the extra immune deficiency (IMD)-like pathway against the Gram-negative bacterial infection. The complex pattern of Rel/NF-κB induced expression suggested that PyNF-κB and PyRel both have specific and cooperative roles in the acute immune responses to bacterial infection.

Characterization of reference genes for RT-qPCR in the desert moss Syntrichia caninervis in response to abiotic stress and desiccation/rehydration

Syntrichia caninervis is the dominant bryophyte of the biological soil crusts found in the Gurbantunggut desert. The extreme desert environment is characterized by prolonged drought, temperature extremes, high radiation and frequent cycles of hydration and dehydration. S. caninervis is an ideal organism for the identification and characterization of genes related to abiotic stress tolerance. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) expression analysis is a powerful analytical technique that requires the use of stable reference genes. Using available S. caninervis transcriptome data, we selected 15 candidate reference genes and analyzed their relative expression stabilities in S. caninervis gametophores exposed to a range of abiotic stresses or a hydration-desiccation-rehydration cycle. The programs geNorm, NormFinder, and RefFinder were used to assess and rank the expression stability of the 15 candidate genes. The stability ranking results of reference genes under each specific experimental condition showed high consistency using different algorithms. For abiotic stress treatments, the combination of two genes (α-TUB2 and CDPK) were sufficient for accurate normalization. For the hydration-desiccation-rehydration process, the combination of two genes (α-TUB1 and CDPK) were sufficient for accurate normalization. 18S was among the least stable genes in all of the experimental sets and was unsuitable as reference gene in S. caninervis. This is the first systematic investigation and comparison of reference gene selection for RT-qPCR work in S. caninervis. This research will facilitate gene expression studies in S. caninervis, related moss species from the Syntrichia complex and other mosses.

Identification and validation of suitable reference genes for RT-qPCR analysis in mouse testis development

RT-qPCR is a commonly used method for evaluating gene expression; however, its accuracy and reliability are dependent upon the choice of appropriate reference gene(s), and there is limited information available on suitable reference gene(s) that can be used in mouse testis at different stages. In this study, using the RT-qPCR method, we investigated the expression variations of six reference genes representing different functional classes (Actb, Gapdh, Ppia, Tbp, Rps29, Hprt1) in mice testis during embryonic and postnatal development. The expression stabilities of putative reference genes were evaluated using five algorithms: geNorm, NormFinder, Bestkeeper, the comparative delta C(t) method and integrated tool RefFinder. Analysis of the results showed that Ppia, Gapdh and Actb were identified as the most stable genes and the geometric mean of Ppia, Gapdh and Actb constitutes an appropriate normalization factor for gene expression studies. The mRNA expression of AT1 as a test gene of interest varied depending upon which of the reference gene(s) was used as an internal control(s). This study suggested that Ppia, Gapdh and Actb are suitable reference genes among the six genes used for RT-qPCR normalization and provide crucial information for transcriptional analyses in future studies of gene expression in the developing mouse testis.

Identification of two secreted ferritin subunits involved in immune defense of Yesso scallop Patinopecten yessoensis

As an important iron storage protein, ferritin plays a crucial role in the iron-withholding defense system. In this study, two secreted ferritin subunits (PyFerS1 and PyFerS2) were identified from the Yesso scallop, Patinopecten yessoensis. The complete DNA sequences of the two ferritins are 7101 and 5359 bp, consisting of seven and five exons, respectively. The full-length cDNAs of PyFerS1 and PyFerS2 are 960 and 956 bp in length, encoding 228 and 220 amino acids, respectively. They have typical ferritin structures, with four long α-helices, one short α-helix and an L-loop. Signal peptides were found at the N-terminus of both ferritins, and phylogenetic analysis showed that they both clustered with secreted mollusc ferritins. PyFerS1 possesses all seven conserved residues of the ferroxidase center, whereas PyFerS2 only has two. Real-time PCR analysis indicated high expression level of PyFerS2 in the D-shaped larvae, and PyFerS1 in both D-shaped larvae and fertilized eggs. In adult scallops, PyFerS1 was only detected in the hepatopancreas, whereas PyFerS2 was detected in both hepatopancreas and mantle. After the scallops were challenged by iron ion or bacteria Vibrio anguillarum, the expression of both PyFerS1 and PyFerS2 was significantly elevated, suggesting they may play a role in scallop innate immune defense. For the first time, secreted ferritins were cloned and comprehensively characterized in bivalve molluscs. It will assist in better understanding of the role of secreted ferritins in bivalve innate immunity.

A scallop IGF binding protein gene: molecular characterization and association of variants with growth traits

Background

Scallops represent economically important aquaculture shellfish. The identification of genes and genetic variants related to scallop growth could benefit high-yielding scallop breeding. The insulin-like growth factor (IGF) system is essential for growth and development, with IGF binding proteins (IGFBPs) serving as the major regulators of IGF actions. Although an effect of IGF on growth was detected in bivalve, IGFBP has not been reported, and members of the IGF system have not been characterized in scallop.

Results

We cloned and characterized an IGFBP (PyIGFBP) gene from the aquaculture bivalve species, Yesso scallop (Patinopecten yessoensis, Jay, 1857). Its full-length cDNA sequence was 1,445 bp, with an open reading frame of 378 bp, encoding 125 amino acids, and its genomic sequence was 10,193 bp, consisting of three exons and two introns. The amino acid sequence exhibited the characteristics of IGFBPs, including multiple cysteine residues and relatively conserved motifs in the N-terminal and C-terminal domains. Expression analysis indicated that PyIGFBP was expressed in all the tissues and developmental stages examined, with a significantly higher level in the mantle than in other tissues and a significantly higher level in gastrulae and trochophore larvae than in other stages. Furthermore, three single nucleotide polymorphisms (SNPs) were identified in this gene. SNP c.1054A>G was significantly associated with both shell and soft body traits in two populations, with the highest trait values in GG type scallops and lowest in AG type ones.

Conclusion

We cloned and characterized an IGFBP gene in a bivalve, and this report also represents the first characterizing an IGF system gene in scallops. A SNP associated with scallop growth for both the shell and soft body was identified in this gene. In addition to providing a candidate marker for scallop breeding, our results also suggest the role of PyIGFBP in scallop growth.