Exploring the Heat-Responsive miRNAs and their Target Gene Regulation in Ruditapes philippinarum Under Acute Heat Stress.

This study aimed to investigate the inherent molecular regulatory mechanisms of Ruditapes philippinarum in response to extremely high-temperature environments and to enhance the sustainable development of the R. philippinarum aquaculture industry. In this study, we established a differential expression profile of miRNA under acute heat stress and identified a total of 46 known miRNAs and 80 novel miRNAs, three of which were detected to be significantly differentially expressed. We analyzed the functions of target genes regulated by differentially expressed miRNAs (DEMs) of R. philippinarum. The findings of the KEGG enrichment analysis revealed that 29 enriched pathways in the group were subjected to acute heat stress. Notably, fatty acid metabolism, FoxO signaling pathway, TGF-ß signaling pathway, and ubiquitin-mediated proteolysis were found to play significant roles in response to acute heat stress. We established a regulatory map of DEMs and their target genes in response to heat stress and constructed the miRNA-mRNA regulation network. This study provides valuable insights into the response of R. philippinarum to high temperature, helping to understand its underlying molecular regulatory mechanisms under high-temperature stress.

Physiological and transcriptomic analysis provides new insights into osmoregulation mechanism of Ruditapes philippinarum under low and high salinity stress.

The Manila clam (Ruditapes philippinarum) is a commercially important marine bivalve, which inhabits the estuarine and mudflat areas. The osmoregulation is of great significance for molluscs adaptation to salinity fluctuations. In this study, we investigated the effects of low salinity (10 psu) and high salinity (40 psu) stress on survival and osmoregulation of the R. philippinarum. The results of physiological parameters showed that the ion (Na+, K+, Cl-) concentrations and Na+/K+-ATPase (NKA) activity of R. philippinarum decreased significantly under low salinity stress, but increased significantly under high salinity stress, indicating that there are differences in physiological adaptation of osmoregulation of R. philippinarum. In addition, we conducted the transcriptome analysis in the gills of R. philippinarum exposed to low (10 psu) and high (40 psu) salinity challenge for 48 h using RNA-seq technology. A total of 153 and 640 differentially expressed genes (DEGs) were identified in the low salinity (LS) group and high salinity (HS) group, respectively. The immune (IAP, TLR6, C1QL4, Ank3), ion transport (Slc34a2, SLC39A14), energy metabolism (PCK1, LDLRA, ACOX1) and DNA damage repair-related genes (Gadd45g, HSP70B2, GATA4) as well as FoxO, protein processing in endoplasmic reticulum and endocytosis pathways were involved in osmoregulation under low salinity stress of R. philippinarum. Conversely, the ion transport (SLC6A7, SLC6A9, SLC6A14, TRPM2), amino acid metabolism (GS, TauD, ABAT, ALDH4A1) and immune-related genes (MAP2K6, BIRC7A, CTSK, GVIN1), and amino acid metabolism pathways (beta-Alanine, Alanine, aspartate and glutamate, Glutathione) were involved in the process of osmoregulation under high salinity stress. The results obtained here revealed the difference of osmoregulation mechanism of R. philippinarum under low and high salinity stress through physiological and molecular levels. This study contributes to the assessment of salinity adaptation of bivalves in the context of climate change and provides useful information for marine resource conservation and aquaculture.

Genome-wide identification and expression analysis of Sox gene family in the Manila clam (Ruditapes philippinarum).

Sox transcription factors are vital in numerous fundamental biological processes. In this study, nine Sox gene family members were discovered in the Ruditapes philippinarum genome, classified into the SoxB1, SoxB2, SoxC, SoxD, SoxE, and SoxF groups, marking the first genome-wide identification of this gene family in R. philippinarum. Analyses of phylogeny, exon-intron structures, and domains bolster the support for their categorization and annotation. Furthermore, transcriptomic analyses across various developmental stages revealed that RpSox4, RpSox5, RpSox9, and RpSox11 were significantly expressed in the D-larval stage. Additionally, investigations into transcriptomes of clams with different shell colors indicated that most sox genes exhibited their highest expression levels in orange clams, followed by zebra, white zebra, and white clams, and the results of transcriptomes analysis in different tissues indicated that 8 Sox genes (except RpSox17) were highly expressed in the mantle tissue. Moreover, qPCR was used to detect the expression of Sox gene in R. philippinarum at different developmental periods, different shell colors and different tissues, and the results showed consistency with those of the transcriptomes. This study's findings lay the groundwork for additional exploration into the role of the Sox gene in melanin production in R. philippinarum shells.

The lgi-miR-2d is Potentially Involved in Shell Melanin Synthesis by Targeting mitf in Manila Clam Ruditapes philippinarum.

Shell color as an important economic trait is also the crucial target trait for breeding and production. MicroRNA (miRNA) is an endogenous small non-coding RNA that can post-transcriptionally regulate the expression of target genes, it plays important roles in many life activities and physiological processes, such as shell color, stress response, and disease traits. In this study, we investigated the function of lgi-miR-2d in shell melanin formation and the expression patterns of lgi-miR-2d and target gene Rpmitf in Manila clam Ruditapes philippinarum. We further explored and verified the relationship between Rpmitf and lgi-miR-2d and identified the expression level of shell color-related gene changes by RNAi and injecting the antagomir of lgi-miR-2d, respectively. Our results indicated that lgi-miR-2d antagomir affected the expression of its target gene Rpmitf. In addition, the dual-luciferase reporter assay was conducted to confirm the direct interaction between lgi-miR-2d and Rpmitf. The results showed that the expression levels of melanin-related genes such as Rpmitf and tyr were significantly decreased in the positive treatment group compared with the blank control group after the Rpmitf dsRNA injection, indicating Rpmitf plays a crucial role in the melanin synthesis pathway. Taken together, we speculated that lgi-miR-2d might be negatively modulating Rpmitf, which might regulate other shell color-related genes, thereby affecting melanin synthesis in R. philippinarum.

Transcriptome analysis provides new insights into the immune response of Ruditapes philippinarum infected with Vibrio alginolyticus.

Manila clam (Ruditapes philippinarum) is a bivalve species with commercial value, but it is easily infected by pathogenic microorganisms in aquaculture, which restricts the shellfish industry. Notably, the impact of Vibrio alginolyticus on clam culture is obvious. In this study, RNA-seq was performed to analyze clam hepatopancreas tissue in 48 h (challenge group, G48h) and 96 h (challenge group, G96h) after infection with V. alginolyticus and 0 h after injection of PBS (control group, C). The results showed that a total of 1670 differentially expressed genes were detected in the G48h vs C group, and 1427 differentially expressed genes were detected in the G96h vs C group. In addition, KEGG analysis showed that DEGs were significantly enriched in pathways such as Lysosome and Mitophagy. Moreover, 15 immune related DEGs were selected for qRT-PCR analysis to verify the accuracy of RNA-seq, and the results showed that the expression level of DEGs was consistent with that of RNA-seq. Therefore, the results obtained in this study provides a preliminary understanding of the immune defense of R. philippinarum and molecular insights for genetic breeding of V. alginolyticus resistance in Manila clam.

Genetic diversity and differentiation of nine populations of the hard clam (Meretrix petechialis) assessed by EST-derived microsatellites

BACKGROUND: Meretrix petechialis is one of the commercially important marine bivalves. In this study, we selected six highly polymorphic EST-derived microsatellite markers to assess the genetic diversity and population differentiation on nine wild populations of Meretrix petechialis. RESULTS: The number of alleles detected per-locus ranged from 4 to 30 (mean NA = 27.5) with a total of 165 alleles. The mean value of observed and expected heterozygosities varied from 0.717 to 0.861 and from 0.797 to 0.856, respectively. Meanwhile, the result of Neighbor-joining and overall FST = 0.214 (P < 0.01) reveled that M. petechialis populations from GX are the farthest populations, a certain degree of genetic variation among individuals in each population and the genetic differentiation is significant. CONCLUSIONS: GX population has high genetic diversity among individual, and there are certain differences in genetic characteristics among different populations. This study will provide a basis for the domestication and cultivation of genetic diversity of M. petechialis population and the protection of clam germplasm resources.

The genetic diversity of wild and cultivated Manila clam (Ruditapes philippinarum) revealed by 29 novel microsatellite markers

Background: Microsatellite loci often used as a genetic tool for estimating genetic diversity population variation in a wide variety of different species. The application of microsatellite markers in genetics and breeding includes investigating the genetic differentiation of wild and cultured populations, assessing and determining the genetic relationship of different populations. The aim of this work is to develop several microsatellite markers via highthroughput sequencing and characterize these markers in commercially important bivalve Ruditapes philippinarum. Results: Among the two populations of R. philippinarum studied, 110 alleles were detected. The number of alleles at the cultured population ranged from 3 to 17 (mean NA = 6.897) and wild population ranged from 2 to 15 (mean NA = 6.793). The observed and expected heterozygosities of cultured population ranged from 0.182 to 0.964, and from 0.286 to 0.900, with an average of 0.647 and 0.692, respectively. The observed and expected heterozygosities of wild population ranged from 0.138 to 1.000, and from 0.439 to 0.906, with an average of 0.674 and 0.693, respectively. The polymorphism information content ranged from 0.341 to 0.910 with an average of 0.687. Sixteen and thirteen microsatellite loci deviated significantly from Hardy­Weinberg equilibrium after correction for multiple tests in cultured and wild population, respectively. Conclusions: Twenty-nine novel microsatellite loci were developed using Illumina paired-end shotgun sequencing and characterized in two population of R. philippinarum.