Intraspecific variations in oyster (Magallana gigas) ploidy does not affect physiological responses to microplastic pollution.

Recent advances in genetic manipulation such as triploid breeding and artificial selection, have rapidly emerged as valuable hatchery methodologies for enhancing seafood stocks. The Pacific oyster Magallana gigas is a leading aquaculture species worldwide and key ecosystem engineer that has received particular attention in this field of science. In light of the growing recognition of the ecological effects of intraspecific variation, oyster polyploids provide a valuable opportunity to assess whether intraspecific diversity affects physiological responses to environmental stressors. While the responses of diploid and triploid oysters to climate change have been extensively investigated, research on their sensitivity to environmental pollution remains scarce. Here, we assess whether genotypic (i.e., ploidy) variation within Magallana gigas affects physiological responses to microplastic pollution. We show that diploid and triploid M. gigas have similar clearance rates and ingest similar amounts of microplastics under laboratory-controlled condition. In addition, they exhibited similar heart rates after prolonged exposure to microplastic leachates. Our findings suggest that intraspecific variations within M. gigas ploidy does not affect oyster responses to microplastic pollution. However, regardless of ploidy, our work highlights significant adverse effects of microplastic leachates on the heart rate of M. gigas and provides evidence of microplastic ingestion in the laboratory.

Intergenerational toxicity of 17α-ethinylestradiol (EE2): Effects of parental exposure on early larval development and transcriptomic profiles in the Sydney rock oyster, Saccostrea glomerata.

This study exposed adult Sydney rock oysters, of either sex or both, to the synthetic estrogen 17α-ethinylestradiol (EE2) at 50 ng/L for 21 days, followed by an examination of developmental endpoints and transcriptomic responses in unexposed larvae. Reduced survival was observed at 1 day post-fertilisation (dpf) in larvae from bi-parental exposure (FTMT). Motile larvae at 2 dpf were fewer from maternal (FTMC), paternal (FCMT), and FTMT exposures. Additionally, shell length at 7 dpf decreased in larvae from FTMC and FTMT parents. RNA sequencing (RNA-seq) revealed 1064 differentially expressed genes (DEGs) in 1-dpf larvae from FTMT parents, while fewer DEGs were detected in larvae from FTMC and FCMT parents, with 258 and 7, respectively. GO and KEGG analyses showed significant enrichment of DEGs in diverse terms and pathways, with limited overlap among treatment groups. IPA results indicated potential inhibition of pathways regulating energy production, larval development, transcription, and detoxification of reactive oxygen species in FTMT larvae. qRT-PCR validation confirmed significant downregulation of selected DEGs involved in these pathways and relevant biological processes, as identified in the RNA-seq dataset. Overall, our results suggest that the intergenerational toxicity of EE2 is primarily maternally transmitted, with bi-parental exposure amplifying these effects.

An amplicon panel for high-throughput and low-cost genotyping of Pacific oyster.

Maintaining genetic diversity in cultured shellfish can be challenging due to high variance in individual reproductive success, founder effects, and rapid genetic drift, but is important to retain adaptive potential and avoid inbreeding depression. To support broodstock management and selective breeding in cultured Pacific oysters (Crassostrea (Magallana) gigas), we developed an amplicon panel targeting 592 genomic regions and SNP variants with an average of 50 amplicons per chromosome. Target SNPs were selected based on elevated observed heterozygosity or differentiation in Pacific oyster populations in British Columbia, Canada. The use of the panel for parentage applications was evaluated using multiple generations of oysters from a breeding program on Vancouver Island, Canada (n = 181) and families selected for Ostreid herpesvirus-1 resistance from the Molluscan Broodstock Program in Oregon, USA (n = 136). Population characterization was evaluated using wild, naturalized, farmed, or hatchery oysters sampled throughout the Northern Hemisphere (n = 189). Technical replicates showed high genotype concordance (97.5%; n = 68 replicates). Parentage analysis found suspected pedigree and sample handling errors, demonstrating the panel's value for quality control in breeding programs. Suspected null alleles were identified and found to be largely population dependent, suggesting population-specific variation impacting target amplification. Null alleles were identified using existing data without the need for pedigree information, and once they were removed, assignment rates increased to 93.0 and 86.0% of possible assignments in the two breeding program datasets. A pipeline for analyzing the amplicon sequence data from sequencer output, amplitools, is also provided.

Genome-wide identification of the HSP70 genes in Pacific oyster Magallana gigas and their response to heat stress.

Heat shock protein 70 (HSP70), the most prominent and well-characterized stress protein in animals, plays an important role in assisting animals in responding to various adverse conditions. In the present study, a total of 113 HSP70 gene family members were identified in the updated genome of Magallana gigas (designated MgHSP70) (previously known as Crassostrea gigas). There were 75, 12, 11, and 8 HSP70s located in the cytoplasm, nucleus, mitochondria, and endoplasmic reticulum, respectively, and 7 HSP70s were located in both the nucleus and cytoplasm. Among 113 MgHSP70 genes, 107 were unevenly distributed in 8 chromosomes of M. gigas with the greatest number in chromosome 07 (61 genes, 57.01%). The MgHSP70 gene family members were mainly assigned into five clusters, among which the HSPa12 subfamily underwent lineage-specific expansion, consisting of 89 members. A total of 68 MgHSP70 genes (60.18%) were tandemly duplicated and formed 30 gene pairs, among which 14 gene pairs were under strong positive selection. In general, the expression of MgHSP70s was tissue-specific, with the highest expression in labial palp and gill and the lowest expression in adductor muscle and hemocytes. There were 35, 31, and 47 significantly upregulated genes at 6, 12, and 24 h after heat shock treatment (28 °C), respectively. The expression patterns of different tandemly duplicated genes exhibited distinct characteristics after shock treatment, indicating that these genes may have different functions. Nevertheless, genes within the same tandemly duplicated group exhibit similar expression patterns. Most of the tandemly duplicated HSP70 gene pairs showed the highest expression levels at 24 h. This study provides a comprehensive description of the MgHSP70 gene family in M. gigas and offers valuable insights into the functions of HSP70 in the mollusc adaptation of oysters to environmental stress.

The development of multiplex PCR assays for the rapid identification of multiple Saccostrea species, and their practical applications in restoration and aquaculture.

BACKGROUND: The ecology and biology of oysters (Ostreidae) across the tropics is poorly understood. Morphological plasticity and shared characteristics among oysters have resulted in the misidentification of species, creating challenges for understanding basic species-specific biological information that is required for restoration and aquaculture. Genetic barcoding has proven essential for accurate species identification and understanding species geographic ranges. To reduce the costs of molecular species identification we developed multiplex assays using the cytochrome c oxidase subunit I (COI or cox1) barcoding gene for the rapid identification of five species of oysters within the genus Saccostrea that are commonly found in Queensland, Australia: Saccostrea glomerata, Saccostrea lineage B, Saccostrea lineage F, Saccostrea lineage G, and Saccostrea spathulata (lineage J). RESULTS: Multiplex assays were successful in species-specific amplification of targeted species. The practical application of these primers was tested on wild spat collected from a pilot restoration project in Moreton Bay, Queensland, with identified species (S. glomerata, lineage B and lineage G) validated by Sanger sequencing. DNA sampling by extraction of oyster pallial fluid was also tested on adult oysters collected from the Noosa estuary in Queensland to assess whether oysters were able to be identified non-destructively. DNA concentrations as low as 1 ng/ µL still amplified in most cases, allowing for identification, and mortality at 6 weeks post pallial fluid collection was low (3 out of 104 sampled oysters). CONCLUSION: These multiplex assays will be essential tools for species identification in future studies, and we successfully demonstrate their practical application in both restoration and aquaculture contexts in Queensland. The multiplex assays developed in this study outline easily replicable methods for the development of additional species-specific primer sets for the rapid identification of other species of Saccostrea found across the Indo-Pacific, which will be instrumental in unravelling the taxonomic ambiguities within this genus in tropical regions.

OysterDB: A Genome Database for Ostreidae.

The molluscan family Ostreidae, commonly known as oysters, is an important molluscan group due to its economic and ecological importance. In recent years, an abundance of genomic data of Ostreidae species has been generated and available in public domain. However, there is still a lack of a high-efficiency database platform to store and distribute these data with comprehensive tools. In this study, we developed an oyster genome database (OysterDB) to consolidate oyster genomic data. This database includes eight oyster genomes and 208,923 protein-coding gene annotations. Bioinformatic tools, such as BLAST and JBrowse, are integrated into the database to provide a user-friendly platform for homologous sequence searching, visualization of genomes, and screen for candidate gene information. Moreover, OysterDB will be continuously updated with ever-growing oyster genomic resources and facilitate future studies for comparative and functional genomic analysis of oysters ( http://oysterdb.com.cn/ ).

Exploring the Feasibility of 16S rRNA Short Amplicon Sequencing-Based Microbiota Analysis for Microbiological Safety Assessment of Raw Oyster.

16S rRNA short amplicon sequencing-based microbiota profiling has been thought of and suggested as a feasible method to assess food safety. However, even if a comprehensive microbial information can be obtained by microbiota profiling, it would not be necessarily sufficient for all circumstances. To prove this, the feasibility of the most widely used V3-V4 amplicon sequencing method for food safety assessment was examined here. We designed a pathogen (Vibrio parahaemolyticus) contamination and/or V. parahaemolyticus-specific phage treatment model of raw oysters under improper storage temperature and monitored their microbial structure changes. The samples stored at refrigerator temperature (negative control, NC) and those that were stored at room temperature without any treatment (no treatment, NT) were included as control groups. The profiling results revealed that no statistical difference exists between the NT group and the pathogen spiked- and/or phage treated-groups even when the bacterial composition was compared at the possible lowest-rank taxa, family/genus level. In the beta-diversity analysis, all the samples except the NC group formed one distinct cluster. Notably, the samples with pathogen and/or phage addition did not form each cluster even though the enumerated number of V. parahaemolyticus in those samples were extremely different. These discrepant results indicate that the feasibility of 16S rRNA short amplicon sequencing should not be overgeneralized in microbiological safety assessment of food samples, such as raw oyster.

The first reporting of prevalence Vibrio species and expression of HSP genes in rayed pearl oyster (Pinctada radiata) under thermal conditions.

The main goal of the present study was to evaluate the influence of thermal exposure on Vibrio population and HSP genes expression (HSP 90, HSP70, and HSP20) in rayed pearl oyster (P. radiata). To this end, the oysters were reared for 30 days at temperatures of 22 °C (control), 25 °C, 27 °C, and 29 °C. The results showed that five dominate Vibrio strains including Vibrio hepatarius, V. harveyi, V. alginolyticus, V. parahaemolyticus, and V. rotiferianus were identified. The highest population of V. parahaemolyticus, V. alginolyticus, and V. harveyi, was found in 29οC group. According to real-time PCR, mantle exhibited the highest expression levels of HSP20, HSP70, and HSP90 genes. A higher level of HSP20 expression was observed at high temperatures (25 °C, 27 °C, and 29 °C) in the gonad and mantle compared to the control group (22 °C) while decrease in HSP90 expression level was recorded in 25 °C, 27 °C, and 29 °C groups. HSP20 expression level in adductor muscle was remarkably down-regulated in 27 °C and 29 °C groups. In this tissue, HSP70 was detected at highest levels in the 29οC group. In mantle, HSP90 gene expression was lowest at 22 °C water temperature. Several Vibrio strains have been identified from pearl Gulf oyster that haven't been previously reported. The identification of dominant Vibrio species is essential for epidemiological management strategies to control and prevent Vibrio outbreaks in pearl oyster farms. The expression pattern of HSP genes differs in rayed pearl oyster tissues due to differences in their thermal tolerance capability and physiological and biological characteristics. The present study provides useful molecular information for the ecological adaptation of rayed pearl oysters after exposure to different temperature levels.

Early gonadal differentiation is associated with the antagonistic action of Foxl2 and Dmrt1l in the Pacific oyster.

As the second largest phylum in the zoological kingdom next to arthropods, the mechanism of gonadal differentiation in mollusca is quite complex. Currently, although much has been carried out on gonadal differentiation in the Pacific oyster, there is still unknown information that needs to be further explored. Here, analysis of the Foxl2 and Dmrt1l expression in samples at different development periods of male and female gonads as well as in annual gonad samples revealed that Log10 (Foxl2/Dmrt1l) values were an effective method for sex identification in oysters. In differentiated gonadal tissue, Log10 (Foxl2/Dmrt1l) values greater than 2 were females and less than 1 for males. Subsequent sequential sampling of the same individuals verified that Log10 (Foxl2/Dmrt1l) values greater than 2 for resting gonads would develop as females and less than 1 would develop as males in the future. Relative expression analysis of Foxl2 and Dmrt1l in the annual samples revealed a negative correlation between Log10 (Foxl2) and Log10 (Dmrt1l). Double fluorescence reporter validation results showed that DMRT1L protein was able to bind the Foxl2 promoter and repress its activity with a weak dosage effect. Antagonism between Dmrt1l and Foxl2 is therefore not restricted to vertebrates, and the competing regulatory networks are of great significance in the maintenance of gonadal sex in oysters after sexual differentiation. This study provides novel ideas and insights into the study of early gonadal differentiation in the adult oyster.

Teneurin and TCAP Phylogeny and Physiology: Molecular Analysis, Immune Activity, and Transcriptomic Analysis of the Stress Response in the Sydney Rock Oyster (Saccostrea glomerata) Hemocytes.

Teneurin C-terminal associated peptide (TCAP) is an ancient bioactive peptide that is highly conserved in metazoans. TCAP administration reduces cellular and behavioral stress in vertebrate and urochordate models. There is little information for invertebrates regarding the existence or function of a TCAP. This study used the Sydney rock oyster (SRO) as a molluscan model to characterize an invertebrate TCAP, from molecular gene analysis to its physiological effects associated with hemocyte phagocytosis. We report a single teneurin gene (and 4 teneurin splice variants), which encodes a precursor with TCAP that shares a vertebrate-like motif, and is similar to that of other molluscan classes (gastropod, cephalopod), arthropods and echinoderms. TCAP was identified in all SRO tissues using western blotting at 1-2 different molecular weights (~22 kDa and ~37kDa), supporting precursor cleavage variation. In SRO hemolymph, TCAP was spatially localized to the cytosol of hemocytes, and with particularly high density immunoreactivity in granules. Based on 'pull-down' assays, the SRO TCAP binds to GAPDH, suggesting that TCAP may protect cells from apoptosis under oxidative stress. Compared to sham injection, the intramuscular administration of TCAP (5 pmol) into oysters modulated their immune system by significantly reducing hemocyte phagocytosis under stress conditions (low salinity and high temperature). TCAP administration also significantly reduced hemocyte reactive oxygen species production at ambient conditions and after 48 h stress, compared to sham injection. Transcriptomic hemocyte analysis of stressed oysters administered with TCAP demonstrated significant changes in expression of genes associated with key metabolic, protective and immune functions. In summary, this study established a role for TCAP in oysters through modulation of physiological and molecular functions associated with energy conservation, stress and cellular defense.

Migration of repetitive DNAs during evolution of the permanent translocation heterozygosity in the oyster plant (Tradescantia section Rhoeo).

Due to translocation heterozygosity for all chromosomes in the cell complement, the oyster plant (Tradescantia spathacea) forms a complete meiotic ring. It also shows Rabl-arrangement at interphase, featured by polar centromere clustering. We demonstrate that the pericentromeric regions of the oyster plant are homogenized in concert by three subtelomeric sequences: 45S rDNA, (TTTAGGG)n motif, and TSrepI repeat. The Rabl-based clustering of pericentromeric regions may have been an excellent device to combine the subtelomere-pericentromere sequence migration (via inversions) with the pericentromere-pericentromere DNA movement (via whole arm translocations) that altogether led to the concerted homogenization of all the pericentromeric domains by the subtelomeric sequences. We also show that the repetitive sequence landscape of interstitial chromosome regions contains many loci consisting of Arabidopsis-type telomeric sequence or of TSrepI repeat, and it is extensively heterozygous. However, the sequence arrangement on some chromosomal arms suggest segmental inversions that are fully or partially homozygous, a fact that could be explained if the inversions started to create linkages already in a bivalent-forming ancestor. Remarkably, the subterminal TSrepI loci reside exclusively on the longer arms that could be due to sharing sequences between similarly-sized chromosomal arms in the interphase nucleus. Altogether, our study spotlights the supergene system of the oyster plant as an excellent model to link complex chromosome rearrangements, evolution of repetitive sequences, and nuclear architecture.

Detection and Genetic Correlation Analysis of Diarrhea Cases and Norovirus in Oysters in Yantai, China.

Background: This study aimed to assess the correlation between Norovirus (NoV), diarrhea, and raw oysters from the eastern coastal areas of Yantai, Shandong, China. Methods: Marine oysters were selected from the three aquatic markets in Laishan district, Yantai City, in March 2019. Meanwhile, 100 fecal samples were collected from patients with diarrhea from the same areas during the same period. Nucleic acids were extracted from these samples and detected by employing reverse transcription polymerase chain reaction (RT-PCR) for NoV GI/GII. The VP1 gene of the coat protein of NoV was amplified by semi-nested RT-PCR and sequenced. Sequence comparison of VP1 was performed with BioEdit software, and the evolutionary tree was constructed with Mega7.0 software. Results: Of the 151 oysters, 42 (27.8%) were positive for NoV. Among them, 32 (21.2%) were GII-positive, 10 (6.6%) were GI-positive, and one GI VP1 sequence was obtained in the oyster samples. Of 100 fecal samples from patients with diarrhea, 38 were GII-positive and 17 were GI-positive. Totally, 19 GII VP 1 sequences and eight GI VP 1 sequences were obtained. Two G1 VP 1 sequences in two fecal samples showed 98.7% nucleotide sequence identity and 99.1% amino acid sequence identity G1 VP 1 acquired in the oyster sample. Conclusions: The results suggest that oysters may be responsible for the spread of NoV in Yantai, Shandong province, China.

Variation in Survival and Gut Microbiome Composition of Hatchery-Grown Native Oysters at Various Locations within the Puget Sound.

The Olympia oyster (Ostrea lurida) of the Puget Sound suffered a dramatic population crash, but restoration efforts hope to revive this native species. One overlooked variable in the process of assessing ecosystem health is association of bacteria with marine organisms and the environments they occupy. Oyster microbiomes are known to differ significantly between species, tissue type, and the habitat in which they are found. The goals of this study were to determine the impact of field site and habitat on the oyster microbiome and to identify core oyster-associated bacteria in the Puget Sound. Olympia oysters from one parental family were deployed at four sites in the Puget Sound both inside and outside of eelgrass (Zostera marina) beds. Using 16S rRNA gene amplicon sequencing of the oyster gut, shell, and surrounding seawater and sediment, we demonstrate that gut-associated bacteria are distinct from the surrounding environment and vary by field site. Furthermore, regional differences in the gut microbiota are associated with the survival rates of oysters at each site after 2 months of field exposure. However, habitat type had no influence on microbiome diversity. Further work is needed to identify the specific bacterial dynamics that are associated with oyster physiology and survival rates. IMPORTANCE This is the first exploration of the microbial colonizers of the Olympia oyster, a native oyster species to the West Coast, which is a focus of restoration efforts. The patterns of differential microbial colonization by location reveal microscale characteristics of potential restoration sites which are not typically considered. These microbial dynamics can provide a more holistic perspective on the factors that may influence oyster performance.

A microcosm study of microbial community profiles during sediment remediation using pyrolyzed oyster shells.

The accumulation of organic and inorganic components in sediments leads to a deterioration in the environment and an imbalance in the coastal ecosystem. Currently, capping is the most effective technology for remediating polluted sediment and restoring ecosystems. A microcosm experiment was designed using pyrolyzed oyster shell (POS). These were mixed in with coastal sediment or added as a capping layer. The results showed that POS effectively decreased pollutants, including PO4-P and NH4-N. Metagenomics analysis was performed using 16S rRNA gene sequencing and the most abundant phyla identified in the POS treated and untreated sediments were Proteobacteria, followed by Firmicutes, Bacteroidetes, Chloroflexi, Fusobacteria, Nitrospirae, and Spirochaetes. The relative abundance of Proteobacteria members of the Class Gammaproteobacteria significantly increased, but Deltaproteobacteria gradually decreased throughout the experiment in POS-covered sediment. This suggests that the POS effectively promoted a shift from anaerobic to facultative anaerobic or aerobic microbial communities in the sediment. Dominant species of facultative anaerobic or microaerophilic bacteria from the order Chromatiales and phylum Nitrospirae were observed in the POS-covered sediment. Based on these study results, it can be concluded that POS is an effective covering material for sediment remediation and restores the microbial communities in sediments.

Diversification of the aquaporin family in geographical isolated oyster species promote the adaptability to dynamic environments.

BACKGROUND: The diversified aquaporin (AQP) family that was derived from gene duplication and subsequent functional differentiation play critical roles in multiple physiological processes and in adaptation to the dynamic environments during the evolutionary process. Oysters are a group of bivalve fauna in Mollusca that were widely distributed around the world and show extraordinary adaptation to harsh environments. However, knowledge is lacking with the diversity and evolution of the AQP family in oysters, even in molluscs. RESULTS: Here, we performed a comprehensive analysis of the AQP family in three geographical isolated oyster species that are native to different environments. Genome distribution and phylogenetic analysis revealed that the expansion of the AQP family in oysters were attributed to tandem duplication. Synteny analysis indicated that large-scale inversions lead to the independent duplication or deletion of the AQPs after speciation. As a consequence, these independent duplication events contributed to the diversification of the AQP family in different oysters. Pore pattern analysis suggested that the duplicated AQPs in oysters were highly diversified in inner surface profiles, implying the subsequent functional differentiation. The comparison conducted based on the transcriptome data demonstrated that the functional differentiated AQP family members in oysters may play critical roles in maintaining the balance between the stationary homeostasis and dynamic environments. CONCLUSIONS: Our observation provides evidence for the correlation between the duplicated and functional differentiated AQP family and the adaptation to stationary life under dynamic environments in oysters. Additionally, it also broadens our knowledge of the evolution of AQP family in molluscs.

Revealing the composition of the eukaryotic microbiome of oyster spat by CRISPR-Cas Selective Amplicon Sequencing (CCSAS).

BACKGROUND: The microbiome affects the health of plants and animals, including humans, and has many biological, ecological, and evolutionary consequences. Microbiome studies typically rely on sequencing ribosomal 16S RNA gene fragments, which serve as taxonomic markers for prokaryotic communities; however, for eukaryotic microbes this approach is compromised, because 18S rRNA gene sequences from microbial eukaryotes are swamped by contaminating host rRNA gene sequences. RESULTS: To overcome this problem, we developed CRISPR-Cas Selective Amplicon Sequencing (CCSAS), a high-resolution and efficient approach for characterizing eukaryotic microbiomes. CCSAS uses taxon-specific single-guide RNA (sgRNA) to direct Cas9 to cut 18S rRNA gene sequences of the host, while leaving protistan and fungal sequences intact. We validated the specificity of the sgRNA on ten model organisms and an artificially constructed (mock) community of nine protistan and fungal pathogens. The results showed that > 96.5% of host rRNA gene amplicons were cleaved, while 18S rRNA gene sequences from protists and fungi were unaffected. When used to assess the eukaryotic microbiome of oyster spat from a hatchery, CCSAS revealed a diverse community of eukaryotic microbes, typically with much less contamination from oyster 18S rRNA gene sequences than other methods using non-metazoan or blocking primers. However, each method revealed taxonomic groups that were not detected using the other methods, showing that a single approach is unlikely to uncover the entire eukaryotic microbiome in complex communities. To facilitate the application of CCSAS, we designed taxon-specific sgRNA for ~16,000 metazoan and plant taxa, making CCSAS widely available for characterizing eukaryotic microbiomes that have largely been neglected. CONCLUSION: CCSAS provides a high-through-put and cost-effective approach for resolving the eukaryotic microbiome of metazoa and plants with minimal contamination from host 18S rRNA gene sequences. Video Abstract.

Elevated seasonal temperature disrupts prooxidant-antioxidant homeostasis and promotes cellular apoptosis in the American oyster, Crassostrea virginica, in the Gulf of Mexico: a field study.

One of the major impacts of climate change has been the marked rise in global temperature. Recently, we demonstrated that high temperatures (1-week exposure) disrupt prooxidant-antioxidant homeostasis and promote cellular apoptosis in the American oyster. In this study, we evaluated the effects of seasonal sea surface temperature (SST) on tissue morphology, extrapallial fluid (EPF) conditions, heat shock protein-70 (HSP70), dinitrophenyl protein (DNP, an indicator of reactive oxygen species, ROS), 3-nitrotyrosine protein (NTP, an indicator of RNS), catalase (CAT), superoxide dismutase (SOD) protein expressions, and cellular apoptosis in gills and digestive glands of oysters collected on the southern Texas coast during the winter (15 °C), spring (24 °C), summer (30 °C), and fall (27 °C). Histological observations of both tissues showed a notable increase in mucus production and an enlargement of the digestive gland lumen with seasonal temperature rise, whereas biochemical analyses exhibited a significant decrease in EPF pH and protein concentration. Immunohistochemical analyses showed higher expression of HSP70 along with the expression of DNP and NTP in oyster tissues during summer. Intriguingly, CAT and SOD protein expressions exhibited significant upregulation with rising seasonal temperatures (15 to 27 °C), which decreased significantly in summer (30 °C), leaving oysters vulnerable to oxidative and nitrative damage. qRT-PCR analysis revealed a significant increase in HSP70 mRNA levels in oyster tissues during the warmer seasons. In situ TUNNEL assay showed a significant increase in apoptotic cells in seasons with high temperature. These results suggest that elevated SST induces oxidative/nitrative stress through the overproduction of ROS/RNS and disrupts the antioxidant system which promotes cellular apoptosis in oysters.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters.

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6-26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.

Reconstruction of the evolutionary biogeography reveal the origins and diversification of oysters (Bivalvia: Ostreidae).

Oysters (Bivalvia: Ostreidae Rafinesque, 1815) live in the intertidal and shallow subtidal areas worldwide. Despite their long evolutionary histories, abundant fossil records, global distribution, and ecological significance, a systematic time-dependent biogeographical analysis of this family is still lacking. Using combined mitochondrial (COI and 16S rRNA) and nuclear (18S rRNA, 28S rRNA, H3 and ITS2) gene makers for 80% (70/88) of the recognized extant Ostreidae, we reconstructed the global phylogenetic and biogeographical relationships throughout the evolutionary history of oysters. The result provided a holistic view of the origin, migration and dispersal patterns of Ostreidae. The phylogenetic results and fossil evidence indicated that Ostreidae originated from the circum-Arctic region in the Early Jurassic. The widening of the Atlantic Ocean and changes in the Tethys Ocean further facilitated their subsequent diversification during the Cretaceous and the Palaeogene periods. In particular, Crassostrea and Saccostrea exhibited relatively low dispersal abilities and their major diversifications were consistent with the tectonic events. Environmental adaptations and reproductive patterns, therefore, should play key roles in the formation of oyster distribution patterners, rather than the dispersal ability of their planktonic larvae. The diversity dynamics inferred by standard phylogenetic are consistent with the fossil record, however, further systematic classification, especially for fossil genus Ostrea, would enhance our understanding on extant and fossil oysters. The present study of the historical biogeography of oysters provides new insights into the evolution and speciation of oysters. Our findings also provide a foundation for the assessment of evolutionary patterns and ecological processes in intertidal and inshore life.

Evidence from oyster suggests an ancient role for Pdx in regulating insulin gene expression in animals.

Hox and ParaHox genes encode transcription factors with similar expression patterns in divergent animals. The Pdx (Xlox) homeobox gene, for example, is expressed in a sharp spatial domain in the endodermal cell layer of the gut in chordates, echinoderms, annelids and molluscs. The significance of comparable gene expression patterns is unclear because it is not known if downstream transcriptional targets are also conserved. Here, we report evidence indicating that a classic transcriptional target of Pdx1 in vertebrates, the insulin gene, is a likely direct target of Pdx in Pacific oyster adults. We show that one insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue. Transcriptomic comparison suggests that this tissue plays a similar role to the vertebrate pancreas. Using ATAC-seq and ChIP, we identify an upstream regulatory element of the cgILP gene which shows binding interaction with cgPdx protein in oyster hepatopancreas and demonstrate, using a cell culture assay, that the oyster Pdx can act as a transcriptional activator through this site, possibly in synergy with NeuroD. These data argue that a classic homeodomain-target gene interaction dates back to the origin of Bilateria.