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

Opsin expression varies across larval development and taxa in pteriomorphian bivalves

Introduction

Many marine organisms have a biphasic life cycle that transitions between a swimming larva with a more sedentary adult form. At the end of the first phase, larvae must identify suitable sites to settle and undergo a dramatic morphological change. Environmental factors, including photic and chemical cues, appear to influence settlement, but the sensory receptors involved are largely unknown. We targeted the protein receptor, opsin, which belongs to large superfamily of transmembrane receptors that detects environmental stimuli, hormones, and neurotransmitters. While opsins are well-known for light-sensing, including vision, a growing number of studies have demonstrated light-independent functions. We therefore examined opsin expression in the Pteriomorphia, a large, diverse clade of marine bivalves, that includes commercially important species, such as oysters, mussels, and scallops.

Methods

Genomic annotations combined with phylogenetic analysis show great variation of opsin abundance among pteriomorphian bivalves, including surprisingly high genomic abundance in many species that are eyeless as adults, such as mussels. Therefore, we investigated the diversity of opsin expression from the perspective of larval development. We collected opsin gene expression in four families of Pteriomorphia, across three distinct larval stages, i.e., trochophore, veliger, and pediveliger, and compared those to adult tissues.

Results

We found larvae express all opsin types in these bivalves, but opsin expression patterns are largely species-specific across development. Few opsins are expressed in the adult mantle, but many are highly expressed in adult eyes. Intriguingly, opsin genes such as retinochrome, xenopsins, and Go-opsins have higher levels of expression in the later larval stages when substrates for settlement are being tested, such as the pediveliger.

Conclusion

Investigating opsin gene expression during larval development provides crucial insights into their intricate interactions with the surroundings, which may shed light on how opsin receptors of these organisms respond to various environmental cues that play a pivotal role in their settlement process.

A widespread picornavirus affects the hemocytes of the noble pen shell (Pinna nobilis), leading to its immunosuppression

Introduction

The widespread mass mortality of the noble pen shell (Pinna nobilis) has occurred in several Mediterranean countries in the past 7 years. Single-stranded RNA viruses affecting immune cells and leading to immune dysfunction have been widely reported in human and animal species. Here, we present data linking P. nobilis mass mortality events (MMEs) to hemocyte picornavirus (PV) infection. This study was performed on specimens from wild and captive populations.

Methods

We sampled P. nobilis from two regions of Spain [Catalonia (24 animals) and Murcia (four animals)] and one region in Italy [Venice (6 animals)]. Each of them were analyzed using transmission electron microscopy (TEM) to describe the morphology and self-assembly of virions. Illumina sequencing coupled to qPCR was performed to describe the identified virus and part of its genome.

Results and discussion

In 100% of our samples, ultrastructure revealed the presence of a virus (20 nm diameter) capable of replicating within granulocytes and hyalinocytes, leading to the accumulation of complex vesicles of different dimensions within the cytoplasm. As the PV infection progressed, dead hemocytes, infectious exosomes, and budding of extracellular vesicles were visible, along with endocytic vesicles entering other cells. The THC (total hemocyte count) values observed in both captive (eight animals) (3.5 × 104-1.60 × 105 ml-1 cells) and wild animals (14 samples) (1.90-2.42 × 105 ml-1 cells) were lower than those reported before MMEs. Sequencing of P. nobilis (six animals) hemocyte cDNA libraries revealed the presence of two main sequences of Picornavirales, family Marnaviridae. The highest number of reads belonged to animals that exhibited active replication phases and abundant viral particles from transmission electron microscopy (TEM) observations. These sequences correspond to the genus Sogarnavirus-a picornavirus identified in the marine diatom Chaetoceros tenuissimus (named C. tenuissimus RNA virus type II). Real-time PCR performed on the two most abundant RNA viruses previously identified by in silico analysis revealed positive results only for sequences similar to the C. tenuissimus RNA virus. These results may not conclusively identify picornavirus in noble pen shell hemocytes; therefore, further study is required. Our findings suggest that picornavirus infection likely causes immunosuppression, making individuals prone to opportunistic infections, which is a potential cause for the MMEs observed in the Mediterranean.

The MYST Family Histone Acetyltransferase SasC Governs Diverse Biological Processes in Aspergillus fumigatus

The conserved MYST proteins form the largest family of histone acetyltransferases (HATs) that acetylate lysines within the N-terminal tails of histone, enabling active gene transcription. Here, we have investigated the biological and regulatory functions of the MYST family HAT SasC in the opportunistic human pathogenic fungus Aspergillus fumigatus using a series of genetic, biochemical, pathogenic, and transcriptomic analyses. The deletion (Δ) of sasC results in a drastically reduced colony growth, asexual development, spore germination, response to stresses, and the fungal virulence. Genome-wide expression analyses have revealed that the ΔsasC mutant showed 2402 significant differentially expressed genes: 1147 upregulated and 1255 downregulated. The representative upregulated gene resulting from ΔsasC is hacA, predicted to encode a bZIP transcription factor, whereas the UV-endonuclease UVE-1 was significantly downregulated by ΔsasC. Furthermore, our Western blot analyses suggest that SasC likely catalyzes the acetylation of H3K9, K3K14, and H3K29 in A. fumigatus. In conclusion, SasC is associated with diverse biological processes and can be a potential target for controlling pathogenic fungi.

Genome-wide identification and transcriptome-based expression profiling of E2 gene family: Implication for potential roles in gonad development of Crassostrea gigas

In this study, we investigated the role of E2 ubiquitin conjugating enzymes (E2) in the Pacific oyster Crassostrea gigas, with a focus on their involvement in gonad development. We identified 34 E2 genes clustered into nine subgroups and 24 subfamilies. The gene structure and intron-exon location were conserved within the same subfamily, but motif variation suggested functional diversity. Tissue transcriptome analyses revealed that most E2 genes were broadly expressed, with UBE2CL showing specific expression in the female gonad. Expression profiling of E2 genes during early embryo-larvae development stages suggested that five E2 genes were highly expressed in early embryo development, indicating their involvement in cell division processes. Furthermore, by profiling the expression of E2 genes in different gonadal developmental stages, we observed a gradual increase in expression for four genes during gametogenesis, with significantly higher expression in the female gonad at the maturation stage. Similarly, five E2 genes displayed elevated expression levels in the male gonad at the maturation stage, indicating their crucial roles in gonadal development and gametogenesis. Our study provides valuable insights into the potential functions of the E2 gene family in C. gigas, shedding light on the molecular mechanisms underlying gonad development in oysters.

Investigation of the highly endangered Pinna nobilis' mass mortalities: Seasonal and temperature patterns of health status, antioxidant and heat stress responses

Recently, P. nobilis populations have suffered a tremendous reduction, with pathogens potentially playing a crucial role. Considering its highly endangered status, mechanisms leading to mass mortalities were examined in one or multiple pathogens infected populations. Thus, seasonal antioxidant enzymatic activities, hsp70 and catalase mRNA levels, were investigated in two different Greek populations, during mass mortality events in summer of 2020. Samples were collected from Fthiotis and Lesvos during February (T^oC 14 ± 1.2 and 15 ± 1 respectively), April (T^oC 18 ± 1.2 and 17 ± 1.3 respectively), and June (T^oC 24.5 ± 1.5 and 21.5 ± 1.5 respectively) 2020. In July of the same year (T^oC 26.5 ± 1.7 in Fthiotis and 24.5 ± 1.7 in Lesvos), no live specimens were found. All biochemical parameters and phylogenetic analysis suggest that pathogen infection increases P. nobilis sensitivity to water temperature, subsequently leading to mass mortality. The latter was obvious in Fthiotis individuals, in which Haplosporidium pinnae was also observed with Mycobacterium spp., compared to Lesvos individuals.

Functional Characterization of the GNAT Family Histone Acetyltransferase Elp3 and GcnE in Aspergillus fumigatus

Post-translational modifications of chromatin structure by histone acetyltransferase (HATs) play a pivotal role in the regulation of gene expression and diverse biological processes. However, the function of GNAT family HATs, especially Elp3, in the opportunistic human pathogenic fungus Aspergillus fumigatus is largely unknown. To investigate the roles of the GNAT family HATs Elp3 and GcnE in the A. fumigatus, we have generated and characterized individual null Δelp3 and ΔgcnE mutants. The radial growth of fungal colonies was significantly decreased by the loss of elp3 or gcnE, and the number of asexual spores (conidia) in the ΔgcnE mutant was significantly reduced. Moreover, the mRNA levels of the key asexual development regulators were also significantly low in the ΔgcnE mutant compared to wild type (WT). Whereas both the Δelp3 and ΔgcnE mutants were markedly impaired in the formation of adherent biofilms, the ΔgcnE mutant showed a complete loss of surface structure and of intercellular matrix. The ΔgcnE mutant responded differently to oxidative stressors and showed significant susceptibility to triazole antifungal agents. Furthermore, Elp3 and GcnE function oppositely in the production of secondary metabolites, and the ΔgcnE mutant showed attenuated virulence. In conclusion, Elp3 and GcnE are associated with diverse biological processes and can be potential targets for controlling the pathogenic fungus.

Characterization of larval shell formation and CgPOU2F1, CgSox5, and CgPax6 gene expression during shell morphogenesis in Crassostrea gigas

Shell formation is a dynamic process involving organic matrix secretion and calcification. In this study, we characterized shell morphogenesis during larval development in Crassostrea gigas. Using scanning electron microscopy (SEM) and fluorescence staining, we demonstrated that shell field, the first morphologically distinguishable shell-forming tissue, became visible soon after enlargement of the blastopore at the anterior end of the trochophore. Shell organic matrix namely protein polysaccharides and calcified structure appeared as a slit at the dorsal side of the embryo. The early shell field began to extend along the dorsal side of the trochophore larvae, and became a saddle shaped shell field that gave rise to the prodissoconch I embryonic shell in the early D-shaped larvae. Subsequently, prodissoconch II shell was formed in the late D-shaped larvae with a characteristic appearance of growth lines. To identify gene expression markers for studying shell formation, we isolated three potential larval shell formation genes CgPOU2F1, CgSox5, and CgPax6 and analyzed their expression during shell morphogenesis. The three potential shell formation genes possessed a similar pattern of expression. Their expression was detected in the shell gland and shell field regions in early D-shaped larvae, hereafter, their expression was detected at the larval mantle edge in the calcified shell stages. Together, these studies provide knowledge of shell morphogenesis in pacific oyster and molecular markers for studying the molecular regulation of biomineralization and shell formation.

Molecular characterization of Pax7 and its role in melanin synthesis in Crassostrea gigas

The paired-box 7 (Pax7) is a transcription factor crucial for skin color polymorphism. However, the mechanism underlying the pigmentation associated with Pax7 in mollusks have yet to be elucidated. In this study, the cDNA sequence of Pax7 in the Pacific oyster Crassostrea gigas (CgPax7) was characterized. Phylogenetically, the identity of deduced amino acid sequence was similar to that of other mollusks and contained 463 amino acids, with conserved features of paired domain (PRD), homeobox domain (HD) and octapeptide. Gene expression analysis revealed that CgPax7 was markedly increased at D-shaped larvae stage and ubiquitously expressed in six examined tissues in adult oyster. The result of whole-mount in situ hybridization (WMISH) showed a restricted pattern of CgPax7 expression on margins of shell valves at D-shaped and umbo larvae stages. Additionally, although CgPax7 silencing had no significant effect on CgMitf expression, it significantly inhibited the expressions of CgPax7, CgTyr, CgTyrp1, CgTyrp2 and CgCdk2, genes involved in Tyr-mediated melanin synthesis. Furthermore, CgPax7 knockdown obviously decreased the tyrosinase activity. Less brown-granules at mantle edge was detected by micrographic examination and melanosomes defect was observed by transmission electron microscopy. It was demonstrated that CgPax7 play a key role in melanin synthesis by regulating Tyr-pathway in C. gigas. These findings indicated the potential framework by which mollusks pigmentation.

The Lysine Demethylases KdmA and KdmB Differently Regulate Asexual Development, Stress Response, and Virulence in Aspergillus fumigatus

Histone demethylases govern diverse cellular processes, including growth, development, and secondary metabolism. In the present study, we investigated the functions of two lysine demethylases, KdmA and KdmB, in the opportunistic human pathogenic fungus Aspergillus fumigatus. Experiments with mutants harboring deletions of genes encoding KdmA (ΔkdmA) and KdmB (ΔkdmB) showed that KdmA is necessary for normal growth and proper conidiation, whereas KdmB negatively regulates vegetative growth and conidiation. In both mutant strains, tolerance to H2O2 was significantly decreased, and the activities of both conidia-specific catalase (CatA) and mycelia-specific catalase (Cat1) were decreased. Both mutants had significantly increased sensitivity to the guanine nucleotide synthesis inhibitor 6-azauracil (6AU). The ΔkdmA mutant produced more gliotoxin (GT), but the virulence was not changed significantly in immunocompromised mice. In contrast, the production of GT and virulence were markedly reduced by the loss of kdmB. Comparative transcriptomic analyses revealed that the expression levels of developmental process-related genes and antioxidant activity-related genes were downregulated in both mutants. Taken together, we concluded that KdmA and KdmB have opposite roles in vegetative growth, asexual sporulation, and GT production. However, the two proteins were equally important for the development of resistance to 6AU.

Heme-Peroxidase 2 Modulated by POU2F1 and SOX5 is Involved in Pigmentation in Pacific Oyster (Crassostrea gigas)

Color polymorphism is frequently observed in molluscan shellfish, while the molecular regulation of shell pigmentation is not well understood. Peroxidase is a key enzyme involved in melanogenesis. Here, we identified a heme-peroxidase 2 gene (CgHPX2), and characterized the expression patterns and transcriptional regulation of CgHPX2 in the Pacific oyster Crassostrea gigas. Tissues expression analysis showed that CgHPX2 was a mantle-specific gene and primarily expressed in the edge mantle in black shell color oyster compared with white shell oyster. In situ hybridization showed that strong signals for CgHPX2 were detected in the both inner and outer surface of the outer fold of mantle in the black shell color oyster, whereas positive signals in white shell oyster were mainly localized in the outer surface of the outer fold of mantle. In the embryos and larvae, a high expression level of CgHPX2 was detected in the trochophore stage in both black and white shell color oysters. The temporal localization of CgHPX2 was mainly detected in the shell gland and edge mantle of trochophore and calcified shell larvae, respectively. In addition, a 2227 bp of 5' flanking region sequence of CgHPX2 was cloned, which contained a presumed core promoter region and many potential transcription factor binding sites. Further luciferase assay experiment confirmed that POU domain, class 2, transcription factor 1 (POU2F1), and SRY-box transcription factor 5 (SOX5) were involved in transcriptional regulation of CgHPX2 gene through binding to its specific promoter region. After CgPOU2F1 and CgSOX5 RNA interference, the CgHPX2 gene expression was significantly decreased. These results suggested that CgPOU2F1 and CgSOX5 might be two important transcription factors that positively regulated the expression of CgHPX2 gene, improving our understanding of the transcriptional regulation of molluscan shell pigmentation.

Acute hypoxic exposure: Effect on hemocyte functional parameters and antioxidant potential in gills of the pacific oyster, Crassostrea gigas

Bivalve mollusks are frequently subjected to fluctuations of dissolved oxygen concentration in the environment which can represent a significant threat to bivalve antioxidant status. In this work the effects of hypoxia on hemocyte reactive oxygen species (ROS) production and level of mitochondrial potential as well as the activity and expression level of catalase (CAT) and superoxide dismutase (SOD) in gills of Crassostrea gigas were investigated after 24 h and 72 h exposure. 24 h hypoxia promoted an increase of mitochondrial membrane potential in agranulocytes and induced ROS accumulation in granulocytes. 72 h exposure substantially decreased hemocyte mitochondrial potential and intracellular ROS level in all hemocyte types. No significant changes in the activity of CAT in gills were observed following both 24 h and 72 h exposure periods compared to control. SOD activity in gills decreased after 72 h exposure to hypoxia but did not change under 24 h hypoxia. Significant up-regulation of SOD gene and no changes in expression level of CAT were observed in all experimental groups. The results indicate an overall shift in antioxidant status in gills and hemocytes of the Pacific oyster that may act as compensatory mechanisms to maintain redox homeostasis after a short-term (24 h) exposure and represent the occurrence of oxidative stress conditions at the end of 72 h hypoxia.

Evaluation of candidate reference genes for quantitative RTqPCR analysis in goldfish (Carassius auratus L.) in healthy and CyHV-2 infected fish

The accuracy of quantitative real time PCR (RTqPCR) can be attained only when a suitable reference gene is used. The gene expression for a particular gene may vary within different cells at different conditions. Hence, the suitability and stability of various potential reference genes have to be determined for expression studies. In this study, we have examined the potential of four different reference genes including β-Actin (ACTB), 18S ribosomal RNA (18S), glyceraldehyde-3P-dehydrogenase (GAPDH), and elongation factor 1 alpha (EF1A_A) in seven different tissues including gill, liver, kidney, spleen, heart, muscle and intestine of goldfish (Carassius auratus). The housekeeping genes were analyzed from healthy fish and in CyHV-2 challenged fish. Based upon the real time PCR results the gene expression varied among the genes and in tissues. The expression levels of the housekeeping genes were then compared and evaluated with the RefFinder web tool which analyses results using four different algorithms - BestKeeper, delta Ct, geNorm and NormFinder. EF1A_A was ranked to be the best gene in healthy fish by BestKeeper and geNorm analysis. The delta Ct and NormFinder algorithm have found 18S to be a stable gene in healthy fish but 18S was given to be least expressed in challenged fish. ACTB was also given as a stable gene by geNorm analysis in both healthy and challenged fish. Also, in CyHV-2 challenged fish, EF1A_A was identified as the best gene by all the three analysis except by BestKeeper analysis, where it has ranked GADPH as the best housekeeping gene. Expression of the four candidate reference genes differed across all tissue types tested, inferring that a thorough study of the reference genes is necessary for cross tissue comparison. These results can be further used in the immune gene response study of goldfish infected with any viral pathogen to develop better health strategies in the disease management of goldfish aquaculture.

Characterization of the mbsA Gene Encoding a Putative APSES Transcription Factor in Aspergillus fumigatus

The APSES family proteins are transcription factors (TFs) with a basic helix-loop-helix domain, known to regulate growth, development, secondary metabolism, and other biological processes in Aspergillus species. In the genome of the human opportunistic pathogenic fungus Aspergillus fumigatus, five genes predicted to encode APSES TFs are present. Here, we report the characterization of one of these genes, called mbsA (Afu7g05620). The deletion (Δ) of mbsA resulted in significantly decreased hyphal growth and asexual sporulation (conidiation), and lowered mRNA levels of the key conidiation genes abaA, brlA, and wetA. Moreover, ΔmbsA resulted in reduced spore germination rates, elevated sensitivity toward Nikkomycin Z, and significantly lowered transcripts levels of genes associated with chitin synthesis. The mbsA deletion also resulted in significantly reduced levels of proteins and transcripts of genes associated with the SakA MAP kinase pathway. Importantly, the cell wall hydrophobicity and architecture of the ΔmbsA asexual spores (conidia) were altered, notably lacking the rodlet layer on the surface of the ΔmbsA conidium. Comparative transcriptomic analyses revealed that the ΔmbsA mutant showed higher mRNA levels of gliotoxin (GT) biosynthetic genes, which was corroborated by elevated levels of GT production in the mutant. While the ΔmbsA mutant produced higher amount of GT, ΔmbsA strains showed reduced virulence in the murine model, likely due to the defective spore integrity. In summary, the putative APSES TF MbsA plays a multiple role in governing growth, development, spore wall architecture, GT production, and virulence, which may be associated with the attenuated SakA signaling pathway.

Tissue-specific evaluation of suitable reference genes for RT-qPCR in the pond snail, Lymnaea stagnalis

Reverse transcription quantitative PCR (RT-qPCR) is a robust technique for the quantification and comparison of gene expression. To obtain reliable results with this method, one or more reference genes must be employed to normalize expression measurements among treatments or tissue samples. Candidate reference genes must be validated to ensure that they are stable prior to use in qPCR experiments. The pond snail (Lymnaea stagnalis) is a common research organism, particularly in the areas of learning and memory, and is an emerging model for the study of biological asymmetry, biomineralization, and evolution and development. However, no systematic assessment of qPCR reference genes has been performed in this animal. Therefore, the aim of our research was to identify stable reference genes to normalize gene expression data from several commonly studied tissues in L. stagnalis as well as across the entire body. We evaluated a panel of seven reference genes across six different tissues in L. stagnalis with RT-qPCR. The genes included: elongation factor 1-alpha, glyceraldehyde-3-phosphate dehydrogenase, beta-actin, beta-tubulin, ubiquitin, prenylated rab acceptor protein 1, and a voltage gated potassium channel. These genes exhibited a wide range of expression levels among tissues. The tissue-specific stability of each of the genes was consistent when measured by the standard stability assessment algorithms: geNorm, NormFinder, BestKeeper, and RefFinder. Our data indicate that the most stable reference genes vary among the tissues that we examined (central nervous system, tentacles, lips, penis, foot, mantle). Our results were generally congruent with those obtained from similar studies in other molluscs. Given that a minimum of two reference genes are recommended for data normalization, we provide suggestions for strong pairs of reference genes for single- and multi-tissue analyses of RT-qPCR data in L. stagnalis.

Identification of candidate reference genes for qRT-PCR normalization studies of salinity stress and injury in Onchidium reevesii

Real-time quantitative reverse transcription-PCR (qRT-PCR) is an undeniably effective tool for measuring levels of gene expression, but the accuracy and reliability of the statistical data obtained depend mainly on the basal expression of selected housekeeping genes in many samples. To date, there have been few analyses of stable housekeeping genes in Onchidium reevesii under salinity stress and injury. In this study, the gene expression stabilities of seven commonly used housekeeping genes, CYC, RPL28S, ACTB, TUBB, EF1a, Ubiq and 18S RNA, were investigated using BestKeeper, geNorm, NormFinder and RefFinfer. Although the results of the four programs varied to some extent, in general, RPL28S, TUBB, ACTB and EF1a were ranked highly. ACTB and TUBB were found to be the most stable housekeeping genes under salinity stress, and EF1a plus TUBB was the most stable combination under injury stress. When analysing target gene expression in different tissues, RPL28S or EF1a should be selected as the reference gene according to the level of target gene expression. Under extreme environmental stress (salinity) conditions, ACTB (0 ppt, 5 ppt, 15 ppt, 25 ppt) and TUBB (35 ppt) are reasonable reference gene choices when expression stability and abundance are considered. Under conditions of 15 ppt salinity and injury stress, our results showed that the best two-gene combination was TUBB plus EF1a. Therefore, we suggest that RPL28S, ACTB and TUBB are suitable reference genes for evaluating mRNA transcript levels. Based on candidate gene expression analysis, the tolerance of O. reevesii to low salinity (low osmotic pressure) is reduced compared to its tolerance to high salinity (high osmotic pressure). These findings will help researchers obtain accurate results in future quantitative gene expression analyses of O. reevesii under other stress conditions.

Reference genes for accurate evaluation of expression levels in Trichophyton interdigitale grown under different carbon sources, pH levels and phosphate levels

Tinea pedis is a type of dermatophytosis caused by anthropophilic keratinolytic fungi such as Trichophyton interdigitale. Quantitative reverse transcription PCR (RT-qPCR) is a reliable and reproducible technique for measuring changes in target gene expression across various biological conditions. A crucial aspect of accurate normalization is the choice of appropriate internal controls. To identify reference genes for accurate evaluation of expression levels in T. interdigitale, the transcription levels of eight candidate reference genes (adp-rf, β-act, ef1-α, gapdh, psm1, sdha, rpl2 and ubc) and one target gene (Tri m4) were analysed by RT-qPCR after growing the dermatophyte under different environmental conditions. The results obtained from expression stability evaluations with NormFinder, geNorm, BestKeeper, and RefFinder software demonstrated that adp-rf and psm1 were the most stable internal control genes across all experimental conditions. The present study constitutes the first report of the identification and validation of reference genes for RT-qPCR normalization for T. interdigitale grown under different environmental conditions resembling the conditions encountered by fungi during invasion of skin.

A SoxC gene related to larval shell development and co-expression analysis of different shell formation genes in early larvae of oyster

Among the potential larval shell formation genes in mollusks, most are expressed in cells surrounding the shell field during the early phase of shell formation. The only exception (cgi-tyr1) is expressed in the whole larval mantle and thus represents a novel type of expression pattern. This study reports another gene with such an expression pattern. The gene encoded a SoxC homolog of the Pacific oyster Crassostrea gigas and was named cgi-soxc. Whole-mount in situ hybridization revealed that the gene was highly expressed in the whole larval mantle of early larvae. Based on its spatiotemporal expression, cgi-soxc is hypothesized to be involved in periostracum biogenesis, biomineralization, and regulation of cell proliferation. Furthermore, we investigated the interrelationship between cgi-soxc expression and two additional potential shell formation genes, cgi-tyr1 and cgi-gata2/3. The results confirmed co-expression of the three genes in the larval mantle of early D-veliger. Nevertheless, cgi-gata2/3 was only expressed in the mantle edge, and the other two genes were expressed in all mantle cells. Based on the spatial expression patterns of the three genes, two cell groups were identified from the larval mantle (tyr1 ⁺/soxc ⁺/gata2/3 ⁺ cells and tyr1 ⁺/soxc ⁺/gata2/3 ^- cells) and are important to study the differentiation and function of this tissue. The results of this study enrich our knowledge on the structure and function of larval mantle and provide important information to understand the molecular mechanisms of larval shell formation.