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 log2(TPM) of RNA-Seq data (Ct = - 0.94 log2(TPM) + 29.67, R2 = 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.

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.

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.

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.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation.

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

The Effect of Temperature on Gonadal Sex Differentiation of Yesso Scallop Patinopecten yessoensis.

Many marine organisms are generally poikilotherms, making seawater temperature one of the most important environmental factors affecting gonadal sex differentiation. Mollusca is the second-largest animal phylum with diverse reproductive systems, but studies on the impact of temperature on sex differentiation are limited to a few sequential hermaphrodites. By combining morphological and molecular analyses, we investigated the effect of temperature on gonadal sex differentiation of a commercially important gonochoristic scallop Patinopecten yessoensis in the field and under laboratory conditions. Based on the relative expression of FoxL2 and Dmrt1L in the gonads of 6- to 12 month-old scallops, we found the scallops start to differentiate at 7 months old in September when the seawater temperature was 21°C. To eliminate the effect of factors other than temperature on sex differentiation, we compared the gonadal development of juvenile scallops at different temperatures (21, 16 and 11°C) under laboratory conditions. After 50 days of treatment, the 11°C group contain more germ cell types, and have higher sex differentiation rates than the 21°C group. But no obvious sex bias was observed. These results suggest that high temperature (21°C) inhibits sex differentiation, whereas low temperature (11°C) accelerates sex differentiation by 2 months for this cold-water species. It also supports juvenile P. yessoensis is gonochoristic rather than protandrous hermaphroditic. Our study addresses for the first time an environmental influence associated with genetic controls on scallop sex differentiation. It will facilitate a better understanding of how environmental factors affect gonadal development in poikilotherms, especially in the less studied molluscs.

Dynamics of DNA Methylation and DNMT Expression During Gametogenesis and Early Development of Scallop Patinopecten yessoensis.

DNA methylation reprograms during gametogenesis and embryo development, which is essential for germ cell specification and genomic imprinting in mammals. Corresponding process remains poorly investigated in molluscs. Here, we examined global DNA methylation level in the gonads of scallop Patinopecten yessoensis during gametogenesis and in embryos/larvae at different stages. DNA methylation level fluctuates during gametogenesis and early development, peaking at proliferative stage of ovary, growing stage of testis, and in blastulae. To understand the mechanisms underlying these changes, we conducted genome-wide characterization of DNMT family and investigated their expression profiles based on transcriptomes and in situ hybridization. Three genes were identified, namely PyDNMT1, PyDNMT2, and PyDNMT3. Expression of PyDnmt3 agrees with DNA methylation level during oogenesis and early development, suggesting PyDNMT3 may participate in de novo DNA methylation that occurs mainly at proliferative stage of ovary and testis, and in blastulae and gastrulae. PyDnmt1 expression is positively correlated with DNA methylation level during spermatogenesis, and is higher at maturation stage of ovary and in 2-8 cell embryos than other stages, implying possible involvement of PyDNMT1 in DNA methylation maintenance during meiosis and embryonic development. This study will facilitate better understanding of the developmental epigenetic reprogramming in bivalve molluscs.

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.

FOXL2 and DMRT1L Are Yin and Yang Genes for Determining Timing of Sex Differentiation in the Bivalve Mollusk Patinopecten yessoensis.

Sex determination and differentiation have long been a research hotspot in metazoans. However, little is known about when and how sex differentiation occurs in most mollusks. In this study, we conducted a combined morphological and molecular study on sex differentiation in the Yesso scallop Patinopecten yessoensis. Histological examination on gonads from 5- to 13-month-old juveniles revealed that the morphological sex differentiation occurred at 10 months of age. To determine the onset of molecular sex differentiation, molecular markers were screened for early identification of sex. The gonadal expression profiles of eight candidate genes for sex determination or differentiation showed that only two genes displayed sexually dimorphic expression, with FOXL2 being abundant in ovaries and DMRT1L in testes. In situ hybridization revealed that both of them were detected in germ cells and follicle cells. We therefore developed LOG10(DMRT1L/FOXL2) for scallop sex identification and confirmed its feasibility in differentiated individuals. By tracing its changes in 5- to 13-month-old juveniles, molecular sex differentiation time was determined: some scallops differentiate early in September when they are 7 months old, and some do late in December when they are 10 months old. Two kinds of coexpression patterns were found between FOXL2 and DMRT1L: expected antagonism after differentiation and unexpected coordination before differentiation. Our results revealed that scallop sex differentiation co-occurs with the formation of follicles, and molecular sex differentiation is established prior to morphological sex differentiation. Our study will assist in a better understanding of the molecular mechanism underlying bivalve sex differentiation.

Genome-wide identification and expression profiling of the SOX gene family in a bivalve mollusc Patinopecten yessoensis.

SOX family is composed of transcription factors that play vital roles in various developmental processes. Comprehensive understanding on evolution of the SOX family requires full characterization of SOX genes in different phyla. Mollusca is the second largest metazoan phylum, but till now, systematic investigation on the SOX family is still lacking in this phylum. In this study, we conducted genome-wide identification of the SOX family in Yesso scallop Patinopecten yessoensis and profiled their tissue distribution and temporal expression patterns in the ovaries and testes during gametogenesis. Seven SOX genes were identified, including SOXB1, B2, C, D, E, F and H, representing the first record in protostomes with SOX members identical to that proposed to exist in the last common ancestor of chordates. Genomic structure analysis identified relatively conserved exon-intron structures, accompanied by intron insertion. Quantitative real-time PCR analysis revealed possible involvement of scallop SOX in various functions, including neuro-sensory cell differentiation, hematopoiesis, myogenesis and gametogenesis. This study represents the first systematic characterization of SOX gene family in Mollusca. It will assist in a better understanding of the evolution and function of SOX family in metazoans.

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.

Scallop genome provides insights into evolution of bilaterian karyotype and development.

Reconstructing the genomes of bilaterian ancestors is central to our understanding of animal evolution, where knowledge from ancient and/or slow-evolving bilaterian lineages is critical. Here we report a high-quality, chromosome-anchored reference genome for the scallop Patinopecten yessoensis, a bivalve mollusc that has a slow-evolving genome with many ancestral features. Chromosome-based macrosynteny analysis reveals a striking correspondence between the 19 scallop chromosomes and the 17 presumed ancestral bilaterian linkage groups at a level of conservation previously unseen, suggesting that the scallop may have a karyotype close to that of the bilaterian ancestor. Scallop Hox gene expression follows a new mode of subcluster temporal co-linearity that is possibly ancestral and may provide great potential in supporting diverse bilaterian body plans. Transcriptome analysis of scallop mantle eyes finds unexpected diversity in phototransduction cascades and a potentially ancient Pax2/5/8-dependent pathway for noncephalic eyes. The outstanding preservation of ancestral karyotype and developmental control makes the scallop genome a valuable resource for understanding early bilaterian evolution and biology.

Transcriptome Sequencing and Comparative Analysis of Ovary and Testis Identifies Potential Key Sex-Related Genes and Pathways in Scallop Patinopecten yessoensis.

Bivalve mollusks have fascinatingly diverse modes of reproduction. However, research investigating sex determination and reproductive regulation in this group of animals is still in its infancy. In this study, transcriptomes of three ovaries and three testes of Yesso scallop were sequenced and analyzed. Transcriptome comparison revealed that 4394 genes were significantly different between ovaries and testes, of which 1973 were ovary-biased (upregulated in the ovaries) and 2421 were testis-biased. Crucial sex-determining genes that were previously reported in vertebrates and putatively present in bivalves, namely FOXL2, DMRT, SOXH, and SOXE, were investigated. The genes all possessed conserved functional domains and were detected in the gonads. Except for PySOXE, the other three genes were significantly differentially expressed between the ovaries and testes. PyFOXL2 was ovary-biased, and PyDMRT and PySOXH were testis-biased, suggesting that these three genes are likely to be key candidates for scallop sex determination/differentiation. Furthermore, GO and KEGG enrichment analyses were conducted for both ovary- and testis-biased genes. Interestingly, both neurotransmitter transporters and GABAergic synapse genes were overrepresented in the ovary-biased genes, suggesting that neurotransmitters, such as GABA and glycine, are likely to participate in scallop ovary development. Our study will assist in better understanding of the molecular mechanisms underlying bivalve sex determination and reproductive regulation.

Identification, characterization and expression profiling of the Tollip gene in Yesso scallop (Patinopecten yessoensis).

Toll-interacting protein (Tollip) is a critical regulator of Toll-like receptor (TLR)-mediated innate immune responses. However, the Tollip gene has not been systematically characterized in shellfish. In this study, we identified and characterized a Tollip gene, PyTollip, in Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine its sequence identities and evolutionary relationships. Compared with Tollip genes from other invertebrate and vertebrate species, the PyTollip gene is highly conserved in its sequence and structural features, except that a unique asparagine residue was found at a conserved site in the C2 domain of PyTollip. Quantitative real-time PCR was used to investigate the expression profiles of PyTollip in different developmental stages, healthy adult tissues, and in hemolymph after Micrococcus luteus and Vibrio anguillarum bacterial infection. Real-time PCR analysis demonstrated differential expression of PyTollip at the acute phase (3 h) after infection with Gram-negative (V. anguillarum) and Gram-positive (M. luteus) bacteria. A second strong response of PyTollip expression was observed 24 h after challenge with V. anguillarum. Collectively, these results provide novel insights into the specific role and response of Tollip and TLR signaling pathways in host immune responses against different bacterial pathogens in bivalves.