Genome-wide association analysis of nutrient traits in the oyster Crassostrea gigas: genetic effect and interaction network

Network analysis combined with genome-wide association study helps identification of genes related to amino acid contents in soybean

BACKGROUND

Additional to total protein content, the amino acid (AA) profile is important to the nutritional value of soybean seed. The AA profile in soybean seed is a complex quantitative trait controlled by multiple interconnected genes and pathways controlling the accumulation of each AA. With a total of 621 soybean germplasm, we used three genome-wide association study (GWAS)-based approaches to investigate the genomic regions controlling the AA content and profile in soybean. Among those approaches, the GWAS network analysis we implemented takes advantage of the relationships between specific AAs to identify the genetic control of AA profile.

RESULTS

For Approach I, GWAS were performed for the content of 24 single AAs under all environments combined. Significant SNPs grouping into 16 linkage disequilibrium (LD) blocks from 18 traits were identified. For Approach II, the individual AAs were grouped by five families according to their metabolic pathways and were examined based on the sum, ratios, and interactions of AAs within the same biochemical family. Significant SNPs grouping into 35 LD blocks were identified, with SNPs associated with traits from the same biochemical family often positioned on the same LD blocks. Approach III, a correlation-based network analysis, was performed to assess the empirical relationships among AAs. Two groups were described by the network topology, Group 1: Ala, Gly, Lys, available Lys (Alys), and Thr and Group 2: Ile and Tyr. Significant SNPs associated with a ratio of connected AAs or a ratio of a single AA to its fully or partially connected metabolic groups were identified within 9 LD blocks for Group 1 and 2 LD blocks for Group 2. Among 40 identified QTL for AA or AA-derived traits, three genomic regions were novel in terms of seed composition traits (oil, protein, and AA content). An additional 24 regions had previously not been specifically associated with the AA content.

CONCLUSIONS

Our results confirmed loci identified from previous studies but also suggested that network approaches for studying AA contents in soybean seed are valuable. Three genomic regions (Chr 5: 41,754,397-41,893,109 bp, Chr 9: 1,537,829-1,806,586 bp, and Chr 20: 31,554,795-33,678,257 bp) were significantly identified by all three approaches. Yet, the majority of associations between a genomic region and an AA trait were approach- and/or environment-specific. Using a combination of approaches provides insights into the genetic control and pleiotropy among AA contents, which can be applied to mechanistic understanding of variation in AA content as well as tailored nutrition in cultivars developed from soybean breeding programs.

Genome-wide association study reveals the underlying regulatory mechanisms of red blood traits in Anadara granosa

BACKGROUND

Anadara granosa, commonly known as the blood clam, exhibits the unusual characteristic of having red blood among invertebrates. There is significant individual variation in blood color intensity among blood clams; individuals with vibrant red blood are deemed healthier and exhibit stronger stress resistance. However, the molecular basis underlying these red blood traits (RBTs) remains poorly understood.

RESULTS

In this study, we performed genome-wide association studies (GWAS) in a population of 300 A. granosa individuals, focusing on RBTs as measured by hemoglobin concentration (HC), total hemocyte count (THC), and heme concentration (HEME). Our analysis identified 18 single nucleotide polymorphisms (SNPs) correlated with RBTs, subsequently selected 117 candidate genes within a 100 kb flanking region of these SNPs, potentially involved in the RBTs of A. granosa. Moreover, we discovered two haplotype blocks specifically associated with THC and HEME. Further analysis revealed eight genes (Septin7, Hox5, Cbfa2t3, Avpr1b, Hhex, Eif2ak3, Glrk, and Rpl35a) that significantly influence RBTs. Notably, a heterozygous A/T mutation in the 3'UTR of Cbfa2t3 was found to promote blood cell proliferation. These genes suggest that the hematopoietic function plays a significant role in the variability of RBTs in A. granosa.

CONCLUSIONS

Our findings reveal a conservation of the regulatory mechanisms of RBTs between blood clams and vertebrates. The results not only provide a scientific basis for selective breeding in blood clams, but also offer deeper insights into the evolutionary mechanisms of RBTs in invertebrates.

Polymorphisms in the cysteine dioxygenase gene and their association with taurine content in the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is rich in taurine, which is crucial for its adaptation to the fluctuating intertidal environment and presents significant potential in improving taurine nutrition and boosting immunity in humans. Cysteine dioxygenase (CDO) is a key enzyme involved in the initial step of taurine biosynthesis and plays a crucial role in regulating taurine content in the body. In the present study, polymorphisms of CDO gene in C. gigas (CgCDO) and their association with taurine content were evaluated in 198 individuals. A total of 24 single nucleotide polymorphism (SNP) loci were identified in the exonic region of CgCDO gene by direct sequencing. Among these SNPs, c.279G>A and c.287C>A were found to be significantly associated with taurine content, with the GG and AA genotype at the two loci exhibiting enhanced taurine accumulation (p < 0.05). Haplotype analysis revealed that the 279GG/287AA haplotype had the highest taurine content of 29.24 mg/g, while the 279AA/287CC haplotype showed the lowest taurine content of 21.19 mg/g. These results indicated that the SNPs of CgCDO gene could influence the taurine content in C. gigas and have potential applications in the selective breeding of high-taurine varieties.

Preliminary analysis of pathways and their implications during salinity stress in abalone

Transcriptome sequencing has offered immense opportunities to study non-model organisms. Abalone is an important marine mollusk that encounters harsh environmental conditions in its natural habitat and under aquaculture conditions; hence, research that increases molecular information to understand abalone physiology and stress response is noteworthy. Accordingly, the study used transcriptome sequencing of the gill tissues of abalone exposed to low salinity stress. The aim is to explore some enriched pathways during salinity stress and the crosstalk and functions of the genes involved in the candidate biological processes for future further analysis of their expression patterns. The data suggest that abalone genes such as YAP/TAZ, Myc, Nkd, and Axin (involved in the Hippo signaling pathway) and PI3K/Akt, SHC, and RTK (involved in the Ras signaling pathways) might mediate growth and development. Thus, deregulation of the Hippo and Ras pathways by salinity stress could be a possible mechanism by which unfavorable salinities influence growth in abalone. Furthermore, PEPCK, GYS, and PLC genes (mediating the Glucagon signaling pathway) might be necessary for glucose homeostasis, reproduction, and abalone meat sensory qualities; hence, a need to investigate how they might be influenced by environmental stress. Genes such as MYD88, IRAK1/4, JNK, AP-1, and TRAF6 (mediating the MAPK signaling pathway) could be useful in understanding abalone's innate immune response to environmental stresses. Finally, the aminoacyl-tRNA biosynthesis pathway hints at the mechanism by which new raw materials for protein biosynthesis are mobilized for physiological processes and how abalone might respond to this process during salinity stress. Low salinity clearly regulated genes in these pathways in a time-dependent manner, as hinted by the heat maps. In the future, qRT-PCR verification and in-depth study of the various genes and proteins discussed would provide enormous molecular information resources for the abalone biology.

Genome-wide association analysis explores the genetic loci of amino acid content in duck's breast muscle

BACKGROUND

Amino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat.

RESULT

In the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55-56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2-21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways.

CONCLUSION

In this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

Molecular cloning, characterisation and molecular modelling of two novel T-synthases from mollusc origin

The glycoprotein-N-acetylgalactosamine β1,3-galactosyltransferase, known as T-synthase (EC 2.4.1.122), plays a crucial role in the synthesis of the T-antigen, which is the core 1 O-glycan structure. This enzyme transfers galactose from UDP-Gal to GalNAc-Ser/Thr. The T-antigen has significant functions in animal development, immune response, and recognition processes. Molluscs are a successful group of animals that inhabit various environments, such as freshwater, marine, and terrestrial habitats. They serve important roles in ecosystems as filter feeders and decomposers but can also be pests in agriculture and intermediate hosts for human and cattle parasites. The identification and characterization of novel carbohydrate active enzymes, such as T-synthase, can aid in the understanding of molluscan glycosylation abilities and their adaptation and survival abilities. Here, the T-synthase enzymes from the snail Pomacea canaliculata and the oyster Crassostrea gigas are identified, cloned, expressed, and characterized, with a focus on structural elucidation. The synthesized enzymes display core 1 β1,3-galactosyltransferase activity using pNP-α-GalNAc as substrate and exhibit similar biochemical parameters as previously characterised T-synthases from other species. While the enzyme from C. gigas shares the same structural parameters with the other enzymes characterised so far, the T-synthase from P. canaliculata lacks the consensus sequence CCSD, which was previously considered indispensable.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis)

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Differentiation and decreased genetic diversity in field contaminated oysters Crassostrea hongkongensis: Identification of selection signatures

The extent to which chemical contamination affects the population structure and genetic diversity of natural populations remains elusive. Here, we used the whole-genome resequencing and transcriptome to diagnose the effects of long-term exposure to multiple elevated chemical pollutants on the population differentiation and genetic diversity in oysters Crassostrea hongkongensis in a typically polluted Pearl River Estuary (PRE) of Southern China. Population structure revealed an obvious differentiation between the PRE oysters and those collected from a nearby clean Beihai (BH) individuals, while no significant differentiation was observed among individuals collected from the three pollution sites within PRE due to the high gene flow. The decreased genetic diversity in the PRE oysters reflected the long-term effects of chemical pollutants. Selective sweeps between BH and PRE oysters revealed that chemical defensome genes, including glutathione S-transferase, zinc transporter, were responsible for their differentiation, sharing common metabolic process of other pollutants. Combined with the genome-wide association analysis, 25 regions containing 77 genes were identified to be responsible for the direct selection regions of metals. Linkage disequilibrium blocks and haplotypes within these regions provided the biomarkers of permanent effects. Our results provide important insights to the genetic mechanisms underlying the rapid evolution under chemical contamination in marine bivalves.

A core of functional complementary bacteria infects oysters in Pacific Oyster Mortality Syndrome

BACKGROUND

The Pacific oyster Crassostrea gigas is one of the main cultivated invertebrate species worldwide. Since 2008, oyster juveniles have been confronted with a lethal syndrome known as the Pacific Oyster Mortality Syndrome (POMS). POMS is a polymicrobial disease initiated by a primary infection with the herpesvirus OsHV-1 µVar that creates an oyster immunocompromised state and evolves towards a secondary fatal bacteremia.

RESULTS

In the present article, we describe the implementation of an unprecedented combination of metabarcoding and metatranscriptomic approaches to show that the sequence of events in POMS pathogenesis is conserved across infectious environments. We also identified a core bacterial consortium which, together with OsHV-1 µVar, forms the POMS pathobiota. This bacterial consortium is characterized by high transcriptional activities and complementary metabolic functions to exploit host's resources. A significant metabolic specificity was highlighted at the bacterial genus level, suggesting low competition for nutrients between members of the core bacteria.

CONCLUSIONS

Lack of metabolic competition between the core bacteria might favor complementary colonization of host tissues and contribute to the conservation of the POMS pathobiota across distinct infectious environments.

Analyzing the mechanism by which oyster peptides target IL-2 in melanoma cell apoptosis based on RNA-seq and m6A-seq

Melanoma is a kind of skin cancer with high malignancy and strong proliferation and invasion abilities. Chemotherapy drugs in the clinic have the disadvantages of high price and high toxicity. Peptides are natural active ingredients that have many functions and are safe and effective. Previous studies have shown that oysters are rich in protein and have antitumor effects. In this study, a high-throughput strategy combined with MALDI TOF/TOF-MS and molecular docking was developed to screen peptides with antitumor functions from oyster hydrolysate. Three dominant peptides were predicted to have similar functions to IL-2 via molecular docking. Then, the activity of the peptides was confirmed in B16 cells, and we found that the three peptides increased the apoptosis of B16 cells. Furthermore, via RNA-seq and m6A-seq of B16 cells treated with the peptides, we found that ILADSAPR downregulates the expression of Pcna, Tlr4, and Ncbp2 and upregulates the expression of Bax, Bad, Pak4, Rasa2, Cct6, and Gbp2. ILADSAPR inhibited B16 cell proliferation and promoted cell apoptosis by regulating the expression of these genes. In addition, the result of metabolic pathway analysis also proved this point. This study provides a preliminary reference for antitumor research on oyster peptides.

Genome-Wide Association and Genomic Prediction of Growth Traits in the European Flat Oyster (Ostrea edulis)

The European flat oyster (Ostrea edulis) is a bivalve mollusc that was once widely distributed across Europe and represented an important food resource for humans for centuries. Populations of O. edulis experienced a severe decline across their biogeographic range mainly due to overexploitation and disease outbreaks. To restore the economic and ecological benefits of European flat oyster populations, extensive protection and restoration efforts are in place within Europe. In line with the increasing interest in supporting restoration and oyster farming through the breeding of stocks with enhanced performance, the present study aimed to evaluate the potential of genomic selection for improving growth traits in a European flat oyster population obtained from successive mass-spawning events. Four growth-related traits were evaluated: total weight (TW), shell height (SH), shell width (SW) and shell length (SL). The heritability of the growth traits was in the low-moderate range, with estimates of 0.45, 0.37, 0.22, and 0.32 for TW, SH, SW and SL, respectively. A genome-wide association analysis revealed a largely polygenic architecture for the four growth traits, with two distinct QTLs detected on chromosome 4. To investigate whether genomic selection can be implemented in flat oyster breeding at a reduced cost, the utility of low-density SNP panels was assessed. Genomic prediction accuracies using the full density panel were high (> 0.83 for all traits). The evaluation of the effect of reducing the number of markers used to predict genomic breeding values revealed that similar selection accuracies could be achieved for all traits with 2K SNPs as for a full panel containing 4,577 SNPs. Only slight reductions in accuracies were observed at the lowest SNP density tested (i.e., 100 SNPs), likely due to a high relatedness between individuals being included in the training and validation sets during cross-validation. Overall, our results suggest that the genetic improvement of growth traits in oysters is feasible. Nevertheless, and although low-density SNP panels appear as a promising strategy for applying GS at a reduced cost, additional populations with different degrees of genetic relatedness should be assessed to derive estimates of prediction accuracies to be expected in practical breeding programmes.

Population Genomics, Transcriptional Response to Heat Shock, and Gut Microbiota of the Hong Kong Oyster Magallana hongkongensis

The Hong Kong oyster Magallana hongkongensis, previously known as Crassostrea hongkongensis, is a true oyster species native to the estuarine-coast of the Pearl River Delta in southern China. The species—with scientific, ecological, cultural, and nutritional importance—has been farmed for hundreds of years. However, there is only limited information on its genetics, stress adaptation mechanisms, and gut microbiota, restricting the sustainable production and use of oyster resources. Here, we present population structure analysis on M. hongkongensis oysters collected from Deep Bay and Lantau Island in Hong Kong, as well as transcriptome analysis on heat shock responses and the gut microbiota profile of M. hongkongensis oysters collected from Deep Bay. Single nucleotide polymorphisms (SNPs), including those on the homeobox genes and heat shock protein genes, were revealed by the whole genome resequencing. Transcriptomes of oysters incubated at 25 °C and 32 °C for 24 h were sequenced which revealed the heat-induced regulation of heat shock protein pathway genes. Furthermore, the gut microbe community was detected by 16S rRNA sequencing which identified Cyanobacteria, Proteobacteria and Spirochaetes as the most abundant phyla. This study reveals the molecular basis for the adaptation of the oyster M. hongkongensis to environmental conditions.

MultiGWAS: An integrative tool for Genome Wide Association Studies in tetraploid organisms

The genome-wide association studies (GWASs) are essential to determine the genetic bases of either ecological or economic phenotypic variation across individuals within populations of the model and nonmodel organisms. For this research question, the GWAS replication testing different parameters and models to validate the results' reproducibility is common. However, straightforward methodologies that manage both replication and tetraploid data are still missing. To solve this problem, we designed the MultiGWAS, a tool that does GWAS for diploid and tetraploid organisms by executing in parallel four software packages, two designed for polyploid data (GWASpoly and SHEsis) and two designed for diploid data (GAPIT and TASSEL). MultiGWAS has several advantages. It runs either in the command line or in a graphical interface; it manages different genotype formats, including VCF. Moreover, it allows control for population structure, relatedness, and several quality control checks on genotype data. Besides, MultiGWAS can test for additive and dominant gene action models, and, through a proprietary scoring function, select the best model to report its associations. Finally, it generates several reports that facilitate identifying false associations from both the significant and the best-ranked association Single Nucleotide Polymorphisms (SNPs) among the four software packages. We tested MultiGWAS with public tetraploid potato data for tuber shape and several simulated data under both additive and dominant models. These tests demonstrated that MultiGWAS is better at detecting reliable associations than using each of the four software packages individually. Moreover, the parallel analysis of polyploid and diploid software that only offers MultiGWAS demonstrates its utility in understanding the best genetic model behind the SNP association in tetraploid organisms. Therefore, MultiGWAS probed to be an excellent alternative for wrapping GWAS replication in diploid and tetraploid organisms in a single analysis environment.

Accumulation of different metals in oyster Crassostrea gigas: Significance and specificity of SLC39A (ZIP) and SLC30A (ZnT) gene families and polymorphism variation

The Zrt/Irt-like proteins (ZIP, SLC39A) and zinc transporters (ZnT, SLC30A) are the two major gene families responsible for the import/export of Zn and other metals. In this study, the mRNA expression levels and genetic variations of eight ZnTs and 14 ZIPs were identified in Crassostrea gigas after exposure to Zn, Cd, Cu, Hg, and Pb. Metal exposure induced reactive oxygen species (ROS) and malondialdehyde (MDA) accumulation and antioxidant enzyme expression. The expanded gene numbers of superoxide dismutase (SOD) in the oysters exhibited diverse expression under exposure to the five metals, and the contrasting expressions of both ZnTs and ZIPs under different metal exposures were observed, revealing their ion-specific responses. Zn and Cu have similar transporters and induce high expression levels of ZnT1, 2, 7, and 9 and ZIP1, 3, 6, 9, 10, 11, and 14. Pb induced high expression levels of ZIP7, and 13 and ZnT5, 6, and 7, which are mainly expressed in the endoplasmic reticulum (ER). Cd induced high expression levels of ZnT1, 2, and 7 and ZIP1, 6, 9, 10, 11, and 13. Hg exposure was found to have little effect on the ZIP and ZnT expression levels. Based on 3784 single nucleotide polymorphisms (SNPs) within the ZnTs and ZIPs, genetic association analysis for Zn accumulation was conducted on 427 oyster samples. The 38 SNPs, which were located within 12 genes, were identified to be associated with Zn content (p < 0.01), explaining the phenotypic variation from 1.61% to 3.37%. One nonsynonymous mutation and related haplotypes were identified within ZIP1, explaining 1.69% of the variation in Zn. Its high expression under Zn exposure revealed its important role in Zn transportation. To the best of our knowledge, this study is the first comprehensive investigation of the transportation mechanisms of ZIPs and ZnTs under different metal exposures and the genetic effect of Zn accumulation in oysters, and provides valuable biomarkers and genetic resources to evaluate environmental metal pollution.

Identification of SNPs involved in Zn and Cu accumulation in the Pacific oyster (Crassostrea gigas) by genome-wide association analysis

Oysters accumulate high concentrations of zinc (Zn) and copper (Cu), which can be transferred to human due to sea food consumption. Breeding new oyster varieties with low Zn and Cu accumulations is one important way to improve food safety. However, the genetic basis for metal accumulation in mollusks is not well understood. To address this issue, oysters collected in the field were used for genome-wide association study (GWAS) and then the identified genes were used for mRNA expressions analysis in laboratory. First, GWAS were conducted for Zn and Cu accumulation in 288 wild Pacific oysters (Crassostrea gigas) farmed in the same ocean environment. The oysters did not show obvious population structure or kinship but exhibited 8.43- and 10.0- fold changes of Zn and Cu contents respectively. GWAS have identified 11 and 12 single nucleotide polymorphisms (SNPs) associated with Zn and Cu, respectively, as well as 16 genes, which were Zn-containing proteins or participated in caveolae-dependent endocytosis. Second, the mRNA expressions of these 16 genes were observed under Zn and Cu exposure. After 9 days of Zn exposure, Zn contents increased 3.1-fold, while the mRNA expression of cell number regulator 3 increased 1.65-fold. Under 9 days of Cu exposure, Cu contents increased 1.97-fold, while the mRNA expression of caveolin-1 decreased 0.61-fold. These provide the evidence for their roles in regulating physiological levels of these two metals. The findings advance our understanding of the genetic basis of Zn and Cu accumulation in mollusks, which can be useful for breeding new, less toxic varieties of oysters.