Genetic and evolutionary patterns of innate immune genes in the Pacific oyster Crassostrea gigas

Morphologic, cytometric, quantitative transcriptomic and functional characterisation provide insights into the haemocyte immune responses of Pacific abalone (Haliotis discus hannai)

In recent years, the abalone aquaculture industry has been threatened by the bacterial pathogens. The immune responses mechanisms underlying the phagocytosis of haemocytes remain unclear in Haliotis discus hannai. It is necessary to investigate the immune mechanism in response to these bacterial pathogens challenges. In this study, the phagocytic activities of haemocytes in H. discus hannai were examined by flow cytometry combined with electron microscopy and transcriptomic analyses. The results of Vibrio parahaemolyticus, Vibrio alginolyticus and Staphylococcus aureu challenge using electron microscopy showed a process during phagosome formation in haemocytes. The phagocytic rate (PP) of S. aureus was higher than the other five foreign particles, which was about 63%. The PP of Vibrio harveyi was about 43%, the PP peak of V. alginolyticus in haemocyte was 63.7% at 1.5 h. After V. parahaemolyticus and V. alginolyticus challenge, acid phosphatase, alkaline phosphatase, total superoxide dismutase, lysozyme, total antioxidant capacity, catalase, nitric oxide synthase and glutathione peroxidase activities in haemocytes were measured at different times, differentially expressed genes (DEGs) were identified by quantitative transcriptomic analysis. The identified DEGs after V. parahaemolyticus challenge included haemagglutinin/amebocyte aggregation factor-like, supervillin-like isoform X4, calmodulin-like and kyphoscoliosis peptidase-like; the identified DEGs after V. alginolyticus challenge included interleukin-6 receptor subunit beta-like, protein turtle homolog B-like, rho GTPase-activating protein 6-like isoform X2, leukocyte surface antigen CD53-like, calponin-1-like, calmodulin-like, troponin C, troponin I-like isoform X4, troponin T-like isoform X18, tumor necrosis factor ligand superfamily member 10-like, rho-related protein racA-like and haemagglutinin/amebocyte aggregation factor-like. Some immune-related KEGG pathways were significantly up-regulated or down-regulated after challenge, including thyroid hormone synthesis, Th17 cell differentiation signalling pathway, focal adhesion, melanogenesis, leukocyte transendothelial migration, inflammatory mediator regulation of TRP channels, ras signalling pathway, rap1 signalling pathway. This study is the first step towards understanding the H. discus hannai immune system by adapting several tools to gastropods and providing a first detailed morpho-functional study of their haemocytes.

Identification, diversity, and evolution analysis of thioester-containing protein family in Pacific oyster (Crassostrea gigas) and immune response to biotic and abiotic stresses

Thioester-containing proteins (TEPs) play a vital role in the innate immune response to biotic and abiotic stresses. In this study, the TEPs in C. gigas were identified, and their gene structure, phylogenetic relationships, collinearity relationships, expression profiles, sequence diversity, and alternative splicing were analyzed. Eight Tep genes were identified in C. gigas genome. Functional analysis and evolutionary relationships indicated a high level of homology to other mollusks TEPs. The transcriptome quantitative analysis results showed that the Tep genes in C. gigas respond to heat stress and Vibrio stress. Alternative splicing analysis revealed four Tep genes (designated A2M_1, CD109_3, CD109_5, complement C3) encode multiple alternative splice variants. Analysis of gene structure and multiple alignments revealed that seven CD109_5 variants are produced through the alternative splicing of the 19th exon, which encodes the highly variable central region. Sequence diversity analysis revealed thirteen missense variants within the 19th exon region of these seven CD109_5 alternative splice variants. Furthermore, the differential alternative splicing analysis showed significant induction of CD109_5, A2M_1 and A2M_2 variants after infection with V. parahaemolyticus. This study explores the Tep genes of C. gigas, providing insights into the molecular mechanisms underlying the involvement of C. gigas TEPs in innate immunity.

Discovering the deep evolutionary roots of serum amyloid A protein family

Deep evolutionary origin of the conserved animal serum amyloid A (SAA) apolipoprotein family leading to yet unknown highly similar SAA-like sequences occurring in certain bacterial genomes is demonstrated in this contribution. Horizontal gene transfer event of corresponding genes between gut bacteria and non-vertebrate animals was discovered in the reconstructed phylogenetic tree obtained with maximum likelihood and neighbor-joining methods, respectively. This detailed phylogeny based on totally 128 complete sequences comprised diverse serum amyloid A isoforms from various animal vertebrate and non-vertebrate phyla and also corresponding genes coding for highly similar proteins from animal gut bacteria. Typical largely conserved sequence motifs and a peculiar structural fold consisting mainly of four α-helices in a bundle within all reconstructed clades of the SAA protein family are discussed with respect to their supposed biological functions in various organisms that contain corresponding genes.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

Virulence of Vibrio alginolyticus Accentuates Apoptosis and Immune Rigor in the Oyster Crassostrea hongkongensis

Vibrio species are ubiquitously distributed in marine environments, with important implications for emerging infectious diseases. However, relatively little is known about defensive strategies deployed by hosts against Vibrio pathogens of distinct virulence traits. Being an ecologically relevant host, the oyster Crassostrea hongkongensis can serve as an excellent model for elucidating mechanisms underlying host-Vibrio interactions. We generated a Vibrio alginolyticus mutant strain (V. alginolyticus^△vscC ) with attenuated virulence by knocking out the vscC encoding gene, a core component of type III secretion system (T3SS), which led to starkly reduced apoptotic rates in hemocyte hosts compared to the V. alginolyticus^WT control. In comparative proteomics, it was revealed that distinct immune responses arose upon encounter with V. alginolyticus strains of different virulence. Quite strikingly, the peroxisomal and apoptotic pathways are activated by V. alginolyticus^WT infection, whereas phagocytosis and cell adhesion were enhanced in V. alginolyticus^△vscC infection. Results for functional studies further show that V. alginolyticus^WT strain stimulated respiratory bursts to produce excess superoxide (O2^•-) and hydrogen peroxide (H₂O₂) in oysters, which induced apoptosis regulated by p53 target protein (p53tp). Simultaneously, a drop in sGC content balanced off cGMP accumulation in hemocytes and repressed the occurrence of apoptosis to a certain extent during V. alginolyticus^△vscC infection. We have thus provided the first direct evidence for a mechanistic link between virulence of Vibrio spp. and its immunomodulation effects on apoptosis in the oyster. Collectively, we conclude that adaptive responses in host defenses are partially determined by pathogen virulence, in order to safeguard efficiency and timeliness in bacterial clearance.

Genomic Landscape of Mutational Biases in the Pacific Oyster Crassostrea gigas

Mutation is a driving force of evolution that has been shaped by natural selection and is universally biased. Previous studies determined genome-wide mutational patterns for several species and investigated the heterogeneity of mutational patterns at fine-scale levels. However, little evidence of the heterogeneity of mutation rates over large genomic regions was shown. Hence, the mutational patterns of different large-scale genomic regions and their association with selective pressures still need to be explored. As the second most species-rich animal phylum, little is known about the mutational patterns in Mollusca, especially oysters. In this study, the mutational bias patterns are characterized by using whole-genome resequencing data in the Crassostrea gigas genome. I studied the genome-wide relative rates of the pair mutations and found that the predominant mutation is GC -> AT, irrespective of the genomic regions. This analysis reveals that mutational biases were associated with gene expression levels across the C. gigas genome. Genes with higher expression levels and breadth expression patterns, longer coding length, and more exon numbers had relatively higher GC -> AT rates. I also found that genes with larger dN/dS values had relatively higher GC -> AT rates. This work represents the first comprehensive research on the mutational biases in Mollusca species. Here, I comprehensively investigated the relationships between mutational biases with some intrinsic genetic factors and evolutionary indicators and proposed that selective pressures are important forces shaping the mutational biases across the C. gigas genome.

First de novo transcriptome assembly of Iwagaki oyster, Crassostrea nippona, and comparative evolutionary analysis of salinity-stress response genes in Crassostrea oysters

Crassostrea nippona is a commercially important oyster species in East Asia for it is edible during the summer when the other oyster species are unavailable. Salinity is one of the important limiting factors to the survival and distribution of this stenohaline species. In this study, 535 million reads (74G data) from C. nippona were produced and assembled into 66,742 transcripts. The number of 19,253 differentially expressed genes (DEGs) under salinity stress were identified as salinity stress-response genes. Through comparative evolutionary analysis in five Crassostrea species from East Asia, salinity stress-response genes were noticed to have higher adaptive evolution rates than other genes. This study presents the first de novo transcriptome of C. nippona. Furthermore, comparative evolutionary analysis implies that salinity plays an important role in speciation of Crassostrea species.

Adaptive Evolution Patterns in the Pacific Oyster Crassostrea gigas

Estimation of adaptive evolution rates at the molecular level is important in evolutionary genomics. However, knowledge of adaptive evolutionary patterns in Mollusca is very scarce, especially for oysters. Such information would help clarify how oysters adapt to pathogen-rich and dynamically changing intertidal environments. In this study, we characterized the patterns of adaptive evolution in the Crassostrea gigas genome, using population diversity analysis and congeneric comparison. Our analysis revealed that gene expression patterns were positively associated with adaptive evolution rates, which suggested that positive selection played an important role in gene evolution. The genes with more exons and alternative splicing events had higher adaptive evolution rates. The rates of adaptive evolution in immune-related and stress-response genes were higher than those in other genes, suggesting that these groups of genes experienced strong positive selection. This study represents the first analysis of adaptive evolution rates in oysters and the first comprehensive study of a Mollusca species. These results provide a system-level investigation of association between adaptive evolution rates with some intrinsic genetic factors. They also suggest that adaptation to pathogens and environmental stressors are important forces driving the adaptive evolution of genes.

Targeted Gene Disruption in Pacific Oyster Based on CRISPR/Cas9 Ribonucleoprotein Complexes

The Pacific oyster (Crassostrea gigas) is a representative bivalve mollusc that is widely cultured in the world. In recent years, it has become an important model species for ecological, evolutionary, and developmental studies because of its ability to survive in extreme environmental conditions as a sessile filter feeder and its classical mosaic pattern of development. Although the complete genome sequence of C. gigas is available and omics data have been rapidly generated for the past few years, the genetic tools for gene functional studies have thus far been limited to RNA interference technology. In this study, we developed a gene editing system for C. gigas based on CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 ribonucleoprotein complexes. Two candidate genes, myostatin (MSTN) and Twist, were selected as targets. After microinjecting CRISPR/Cas9 ribonucleoprotein complexes into fertilized eggs, CRISPR-induced indel mutations were detected in the target genes. The CRISPR/Cas9-induced mutations were predominantly small indel mutations ranging in size from 1 to 24 bp in these two target genes. These results demonstrate that CRISPR/Cas9 can be successfully used as an effective targeted gene editing system in C. gigas. The method reported here provides a powerful tool for gene functional studies in oysters and other marine bivalves, and potentially as a new technology for genetic engineering to improve oyster traits for aquaculture.

Relationship Among Intron Length, Gene Expression, and Nucleotide Diversity in the Pacific Oyster Crassostrea gigas

Crassostrea gigas is a model mollusk, but its genetic features have not been studied comprehensively. In this study, we used whole-genome resequencing data to identify and characterize nucleotide diversity and population recombination rate in a diverse collection of 21 C. gigas samples. Our analyses revealed that C. gigas harbors both extremely high genetic diversity and recombination rates across the whole genome as compared with those of the other taxa. The noncoding regions, introns, intergenic spacers, and untranslated regions (UTRs) showed a lower level diversity than the synonymous sites. The larger introns tended to have lower diversity. Moreover, we found a negative association of the non-synonymous diversity with gene expression, which suggested that purifying selection played an important role in shaping genetic diversity. The nucleotide diversity at the 100- and 50-kb levels was positively correlated with population recombination rates, which was expected if the diversity was shaped by purifying selection or hitchhiking of advantageous mutants. Our work gives a general picture of the oyster's polymorphism pattern and its association with recombination rates.