Differential mantle transcriptomics and characterization of growth-related genes in the diploid and triploid pearl oyster Pinctada fucata

Comparative transcriptomics identifies genes underlying growth performance of the Pacific black-lipped pearl oyster Pinctada margaritifera

BACKGROUND

In bivalves, the rate at which organisms grow is a major functional trait underlying many aspects of their commercial production. Growth is a highly polygenic trait, which is typically regulated by many genes with small to moderate effects. Due to its complexity, growth variability in such shellfish remains poorly understood. In this study, we aimed to investigate differential gene expression among spat of the pearl oyster Pinctada margaritifera with distinct growth phenotypes.

RESULTS

We selected two groups of P. margaritifera spat belonging to the same F2 cohort based on their growth performance at 5.5 months old. Transcriptome profile analysis identified a total of 394 differentially expressed genes between these Fast-growing (F) and Slow-growing (S) phenotypes. According to functional enrichment analysis, S oysters overexpressed genes associated with stress-pathways and regulation of innate immune responses. In contrast, F oysters up-regulated genes associated with cytoskeleton activity, cell proliferation, and apoptosis. Analysis of genome polymorphism identified 16 single nucleotide polymorphisms (SNPs) significantly associated with the growth phenotypes. SNP effect categorization revealed one SNP identified for high effect and annotated for a stop codon gained mutation. Interestingly, this SNP is located within a gene annotated for scavenger receptor class F member 1 (SRF1), which is known to modulate apoptosis. Our analyses also revealed that all F oysters showed up-regulation for this gene and were homozygous for the stop-codon mutation. Conversely, S oysters had a heterozygous genotype and a reduced expression of this gene.

CONCLUSIONS

Altogether, our findings suggest that differences in growth among the same oyster cohort may be explained by contrasted metabolic allocation between regulatory pathways for growth and the immune system. This study provides a valuable contribution towards our understanding of the molecular components associated with growth performance in the pearl oyster P. margaritifera and bivalves in general.

Exploration of Molecular Mechanisms of Immunity in the Pacific Oyster (Crassostrea gigas) in Response to Vibrio alginolyticus Invasion

Over the years, oysters have faced recurring mass mortality issues during the summer breeding season, with Vibrio infection emerging as a significant contributing factor. Tubules of gill filaments were confirmed to be in the hematopoietic position in Crassostrea gigas, which produce hemocytes with immune defense capabilities. Additionally, the epithelial cells of oyster gills produce immune effectors to defend against pathogens. In light of this, we performed a transcriptome analysis of gill tissues obtained from C. gigas infected with Vibrio alginolyticus for 12 h and 48 h. Through this analysis, we identified 1024 differentially expressed genes (DEGs) at 12 h post-injection and 1079 DEGs at 48 h post-injection. Enrichment analysis of these DEGs revealed a significant association with immune-related Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. To further investigate the immune response, we constructed a protein-protein interaction (PPI) network using the DEGs enriched in immune-associated KEGG pathways. This network provided insights into the interactions and relationships among these genes, shedding light on the underlying mechanisms of the innate immune defense mechanism in oyster gills. To ensure the accuracy of our findings, we validated 16 key genes using quantitative RT-PCR. Overall, this study represents the first exploration of the innate immune defense mechanism in oyster gills using a PPI network approach. The findings provide valuable insights for future research on oyster pathogen control and the development of oysters with enhanced antimicrobial resistance.

Comparative Transcriptome Profiling of CMS-D2 and CMS-D8 Systems Characterizes Fertility Restoration Genes Network in Upland Cotton

Resolving the genetic basis of fertility restoration for cytoplasmic male sterility (CMS) can improve the efficiency of three-line hybrid breeding. However, the genetic determinants of male fertility restoration in cotton are still largely unknown. This study comprehensively compared the full-length transcripts of CMS-D2 and CMS-D8 systems to identify potential genes linked with fertility restorer genes Rf₁ or Rf₂. Target comparative analysis revealed a higher percentage of differential genes in each restorer line as compared to their corresponding sterile and maintainer lines. An array of genes with specific expression in the restorer line of CMS-D2 had functional annotations related to floral development and pathway enrichments in various secondary metabolites, while specifically expressed genes in the CMS-D8 restorer line showed functional annotations related to anther development and pathway enrichment in the biosynthesis of secondary metabolites. Further analysis identified potentially key genes located in the target region of fertility restorer genes Rf₁ or Rf₂. In particular, Ghir_D05G032450 can be the candidate gene related to restorer gene Rf₁, and Ghir_D05G035690 can be the candidate gene associated with restorer gene Rf₂. Further gene expression validation with qRT-PCR confirmed the accuracy of our results. Our findings provide useful insights into decoding the potential regulatory network that retrieves pollen fertility in cotton and will help to further reveal the differences in the genetic basis of fertility restoration for two CMS systems.

Transcriptome profiling explores the immune defence mechanism of triploid Pacific oyster (Crassostrea gigas) blood against Vibrio alginolyticus based on protein interaction networks

Triploid oysters have provided the oyster industry with many benefits, such as fast growth rates, meat quality improvement, and increased oyster production and economic benefits, since the first report on triploid oysters was published. The development of polyploid technology has remarkably increased the output of triploid oysters to meet the increasing demand of consumers for Crassostrea gigas in the past decades. At present, research on triploid oyster has mainly focused on breeding and growth, but studies on the immunity of triploid oysters are limited. According to recent reports, Vibrio alginolyticus is a highly virulent strain that can cause disease and death in shellfish, shrimp, as well as serious economic losses. V. alginolyticus may be a reason why oysters die during summer. Therefore, using V. alginolyticus to explore the resistance and immune defense mechanisms of triploid oysters against pathogens presents practical significance. Transcriptome analysis of gene expression was performed in triploid C. gigas at 12 and 48 h after infection with V. alginolyticus, and the respective 2257 and 191 differentially expressed genes (DEGs) were identified. The results of GO and KEGG enrichment analyses showed that multiple significantly enriched GO terms and KEGG signaling pathways are associated with immunity. A protein-protein interaction network was constructed to investigate the interaction relationship of immune-related genes. Finally, we verified the expression situation of 16 key genes using quantitative RT-PCR. This study is the first to use the PPI network in exploring the immune defense mechanism of triploid C. gigas blood to fill the gap in the immune mechanism of triploid oysters and other mollusks, and provide valuable reference for future triploid farming and pathogen prevention and control.

Mantle Transcriptome Provides Insights into Biomineralization and Growth Regulation in the Eastern Oyster (Crassostrea virginica)

Growth of the eastern oyster Crassostrea virginica, a major aquaculture species in the USA, is highly variable and not well understood at molecular levels. As growth of mollusks is confined in shells constructed by the mantle, mantle transcriptomes of large (fast-growing) and small (slow-growing) eastern oysters were sequenced and compared in this study. Transcription was observed for 31,186 genes, among which 104 genes were differentially expressed between the large and small oysters, including 48 upregulated and 56 downregulated in large oysters. Differentially expressed genes (DEGs) included genes from diverse pathways highlighting the complexity of shell formation and growth regulations. Seventeen of the 48 upregulated DEGs were related to shell matrix formation, most of which were upregulated in large oysters, indicating that large oysters are more active in biomineralization and shell formation. Genomic and transcriptomic analyses identified 22 genes encoding novel polyalanine containing proteins (Pacps) with characteristic motifs for matrix function that are tandemly duplicated on one chromosome, all specifically expressed in mantle and at higher levels in large oysters, suggesting that these expanded Pacps play important roles in shell formation and growth. Analysis of sequence variation identified 244,964 SNPs with 328 associated with growth. This study provides novel candidate genes and markers for shell formation and growth, and suggests that genes related to shell formation are important for the complex regulation of growth in the eastern oyster and possibly other bivalve mollusks. Results of this study show that both transcriptional modulation and functional polymorphism are important in determining growth.

Transcriptomic analysis provides insights into candidate genes and molecular pathways involved in growth of Manila clam Ruditapes philippinarum

Growth is one of the most important traits of aquaculture breeding programs. Understanding the mechanisms underlying growth differences between individuals can contribute to improving growth rates through more efficient breeding schemes. Ruditapes philippinarum is an economically important marine bivalve. In order to gain insights into the molecular mechanisms to growth variability in marine shellfish, we conducted the transcriptome sequencing and examined the expression differences in growth-related gene and molecular pathways involved in growth trait of R. philippinarum. In this study, we investigated the molecular and gene expression differences in fast-growing and slow-growing Manila clam and focused on the analysis of the differential expression patterns of specific genes associated with growth by RNA-seq and qPCR analysis. A total of 61 differentially expressed genes (DEGs) were captured significantly differentially expressed, and were categorized into Ras signaling pathway, hedgehog signaling pathway, AMPK signaling pathway, p53 signaling pathway, regulation of actin cytoskeleton, focal adhesion, mTOR signaling pathway, VEGF signaling pathway, and TGF-beta signaling pathway. A total of 34 growth-related genes were validated significantly and up/downregulated at fast growing and slow growing of R. philippinarum. Functional enrichment analysis revealed the insulin signaling pathway, PI3K-Akt signaling pathway, and mTOR signaling pathway play pivotal roles in molecular function and regulation of growth trait in R. philippinarum. The growth-related genes and pathways obtained here provide important insights into the molecular basis of physiological acclimation, metabolic activity, and growth variability in marine bivalves.

PacBio single molecule long-read sequencing provides insight into the complexity and diversity of the Pinctada fucata martensii transcriptome

BACKGROUND

The pearl oyster Pinctada fucata martensii is an economically valuable shellfish for seawater pearl production, and production of pearls depends on its growth. To date, the molecular mechanisms of the growth of this species remain poorly understood. The transcriptome sequencing has been considered to understanding of the complexity of mechanisms of the growth of P. f. martensii. The recently released genome sequences of P. f. martensii, as well as emerging Pacific Bioscience (PacBio) single-molecular sequencing technologies, provide an opportunity to thoroughly investigate these molecular mechanisms.

RESULTS

Herein, the full-length transcriptome was analysed by combining PacBio single-molecule long-read sequencing (PacBio sequencing) and Illumina sequencing. A total of 20.65 Gb of clean data were generated, including 574,561 circular consensus reads, among which 443,944 full-length non-chimeric (FLNC) sequences were identified. Through transcript clustering analysis of FLNC reads, 32,755 consensus isoforms were identified, including 32,095 high-quality consensus sequences. After removing redundant reads, 16,388 transcripts were obtained, and 641 fusion transcripts were derived by performing fusion transcript prediction of consensus sequences. Alternative splicing analysis of the 16,388 transcripts was performed after accounting for redundancy, and 9097 gene loci were detected, including 1607 new gene loci and 14,946 newly discovered transcripts. The original boundary of 11,235 genes on the chromosomes was corrected, 12,025 complete open reading frame sequences and 635 long non-coding RNAs (LncRNAs) were predicted, and functional annotation of 13,482 new transcripts was achieved. Two thousand three hundred eighteen alternative splicing events were detected. A total of 228 differentially expressed transcripts (DETs) were identified between the largest (L) and smallest (S) pearl oysters. Compared with the S, the L showed 99 and 129 significantly up-and down-regulated DETs, respectively. Six of these DETs were further confirmed by quantitative real-time RT-PCR (RT-qPCR) in independent experiment.

CONCLUSIONS

Our results significantly improve existing gene models and genome annotations, optimise the genome structure, and in-depth understanding of the complexity and diversity of the differential growth patterns of P. f. martensii.

Characterization of Chinese hamster ovary cells with disparate chromosome numbers: Reduction of the amount of mRNA relative to total protein

Chromosomes in Chinese hamster ovary (CHO) cells are labile. We have shown that high-chromosome-number CHO cells have greater potential to become robust producers of recombinant proteins. One explanation being the increase in transgene integration sites. However, high-chromosome-number cell clones produce more IgG3 following culture of single-cell clones, even under conditions that yield the same number of integrations as cells with normal chromosome numbers. Here, we characterized high-chromosome-number cells by transcriptome analysis. RNA standards were used to normalize transcriptomes of cells that had different chromosome numbers. Our results demonstrate that the mRNA ratio of β-actin and many other genes in high-chromosome-number cells to that in normal-chromosome-number cells per cell (normalized to RNA standards) was smaller than the equivalent genomic size and cell volume ratios. Many genes encoding membrane proteins are more highly expressed in high-chromosome-number cells, probably due to differences in cell size caused by the increase in chromosomes. In addition, genes related to histone modification and lipid metabolism are differentially expressed. The reduced transcript level required per protein produced in total and the different intracellular signal transductions might be key factors for antibody production.

Comparative Transcriptome Analysis Reveals Molecular Basis Underlying Fast Growth of the Selectively Bred Pacific Oyster, Crassostrea gigas

Fast growth is one of the most desired traits for all food animals, which affects the profitability of animal production. The Pacific oyster, Crassostrea gigas, is an important aquaculture shellfish around the world with the largest annual production. Growth of the Pacific oyster has been greatly improved by artificial selection breeding, but molecular mechanisms underlying growth remains poorly understood, which limited the molecular integrative breeding of fast growth with other superior traits. In this study, comparative transcriptome analyses between the fast-growing selectively bred Pacific oyster and unselected wild Pacific oysters were conducted by RNA-Seq. A total of 1,303 protein-coding genes differentially expressed between fast-growing oysters and wild controls were identified, of which 888 genes were expressed at higher levels in the fast-growing oysters. Functional analysis of the differentially expressed genes (DEGs) indicated that genes involved in microtubule motor activity and biosynthesis of nucleotides and proteins are potentially important for growth in the oyster. Positive selection analysis of genes at the transcriptome level showed that a significant number of ribosomal protein genes had undergone positive selection during the artificial selection breeding process. These results also indicated the importance of protein biosynthesis and metabolism for the growth of oysters. The alternative splicing (AS) of genes was also compared between the two groups of oysters. A total of 3,230 differential alternative splicing events (DAS) were identified, involved in 1,818 genes. These DAS genes were associated with specific functional pathways related to growth, such as “long-term potentiation,” “salivary secretion,” and “phosphatidylinositol signaling system.” The findings of this study will be valuable resources for future investigation to unravel molecular mechanisms underlying growth regulation in the oyster and other marine invertebrates and to provide solid support for breeding application to integrate fast growth with other superior traits in the Pacific oyster.

Towards a better understanding of allograft-induced stress response in the pearl oyster Pinctada fucata martensii: Insights from iTRAQ-based comparative proteomic analysis

Implantation of a spherical nucleus into a recipient oyster is a critical step in artificial pearl production. The implanted nucleus is known to trigger cellular stress responses at several levels, yet the molecular mechanism underpinning physiological adaptation of the pearl oysters to nucleus implantation is still poorly understood. In this study, we took advantage of the iTRAQ-based proteomics and LC-MS/MS approach to look into allograft induced gene regulation at the protein expression level in the pearl oyster Pinctada fucata martensii, across a period of 30 days following nucleus implantation. A wide variety of proteins, including a group of immune-related proteins such as E3 ubiquitin-ligase and heat shock proteins, exhibited differential expression in response to the surgical operation. Further comparisons between different sampling points revealed that GO terms including "translation" and "oxidation-reduction process" and KEGG pathways including "glycolysis/gluconeogenesis" and "pyruvate metabolism" were significantly enriched at several time points, indicating the important roles of these molecular events in the stress response of pearl oysters to nucleus implantation. In addition, considerable discrepancy between protein expression level and gene transcript abundancy was identified, as only a few genes showed at least 2-fold expression changes at both proteomic and transcriptomic levels. The result implies that post-transcriptional gene regulation for the key proteins may represent an important aspect of allograft-induced stress response in the pearl oysters. Taken together, the data obtained would contribute to a deeper understanding of the molecular mechanisms enabling stress adaptation of the pearl oysters in response to nucleus implantation.