Click
here to close Hello! We notice that
you are using Internet Explorer, which is not supported by Echinobase
and may cause the site to display incorrectly. We suggest using a
current version of Chrome,
FireFox,
or Safari.
Front Genet
2021 Jan 01;12:795706. doi: 10.3389/fgene.2021.795706.
Show Gene links
Show Anatomy links
Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions.
Chan J
,
Wang L
,
Li L
,
Mu K
,
Bushek D
,
Xu Y
,
Guo X
,
Zhang G
,
Zhang L
.
Abstract
Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.
FIGURE 1. Schematic diagram of experimental design. FSW, filtered seawater.
FIGURE 2. Transcriptional changes induced by P. marinus in C. virginica.
(A), Volcano plots of transcriptional change under short-term and long-term P. marinus challenge. Blue points represent downregulated genes, red points represent upregulated genes, and black points represent non-differentially expressed genes. (B), Functional annotation of differentially expressed genes under P. marinus challenge. (C), Immune-related genes significantly up-regulated after dermo challenge.
FIGURE 3. Comparative transcriptome analysis between C. gigas and C. virginica under P. marinus challenge. (A), Number of immune-related genes identified (Left) and number of immune-related genes differentially expressed (Right) in C. gigas (Cg) and C. virginica (Cv). (B), Venn diagrams of differentially expressed immune genes under P. marinus challenge between the two oyster species. (C), Unique and differentially expressed immune genes induced by P. marinus in two oyster species. (D), Transcriptional change of differentially expressed immune genes shared by two oyster species, including pattern recognition receptors (e.g., TLRs) and immune effectors (e.g., TRAFs).
FIGURE 4. GO enrichment analysis of the differentially expressed genes after Dermo challenge between C. virginica
(A) and C. gigas
(B). Red and blue columns represent functional enrichments of up-regulated and down-regulated genes, respectively.
FIGURE 5. Expansion of TLRs and divergent expression patterns under P. marinus challenge between C. gigas and C. virginica.
(A), Phylogenetic tree constructed with the maximum likelihood method showing lineage-specific expansion of TLRs in C. gigas (red font) and C. virginica (blue font). Red and blue asterisks represent the differentially expressed TLRs in C. gigas and C. virginica, respectively. (B), Comparison of TLRs with different domain architectures in: human (Homo sapiens-Hs), sea urchin (Strongylocentrotus purpuratus-SP), fruit fly (Drosophila melanogaster-Dm), Pacific oyster (Crassostrea gigas-Cg), Eastern oyster (Crassostrea virginica-Cv), and staghorn coral (Acropora digitifera-Ad). LRRCT in red, LRRNT in blue and LRR in yellow. Toll/interleukin-1 receptor (TIR) domains are shown as gold triangles. (C), Diverse expression pattern of the differentially expressed TLRs, marked with asterisks in (A), during P. marinus challenge in C. gigas
(Top) and C. virginica
(Bottom).
FIGURE 6. Tandem duplication and positive selection of TLRs in C. virginica. (A), Tandem repeats as the major source of TLRs duplication in C. virginica.
(B), The largest tandem repeat cluster (Top) and divergent expression patterns (Bottom) of tandemly linked TLRs in NC_035786.1 under P. marinus challenge. The wide arrows represent TLR genes, and genes with red background are differentially expressed under P. marinus challenge. A pseudogene LOC111105057 was detected in this repeat cluster. (C), Evolutionary relationships of duplicated TLRs. The calculated dN/dS (ω) values (gold) and bootstrap values are shown for each branch. The branches with ω values >1.0 are marked with red rectangle. The differential genes under P. marinus challenge that are positively selected are boxed with red dash lines. (D), Structural modeling of TLR depicting sites (blue) under positive selection. Positive selection sites (Asp27, Glu62, Asn199) are represented by black arrows.
FIGURE 7. Expansion of C1qDCs and divergent expression patterns under P. marinus challenge between C. virginica and C. gigas. (A), Phylogenetic tree constructed with the maximum likelihood method showing lineage-specific expansion of C1qDCs in C. virginica (red font) and C. gigas (gray font). Different color backgrounds represent different types of C1qDC domain composition, and the red and gray asterisks represent differentially expressed C1qDCs under P. marinus challenge in C. virginica and C. gigas, respectively. (B), Domain composition of the “galactose-binding lectin” containing C1qDC gene. (C), A C1qDC gene (LOC111120224) with “galactose-binding lectin” domain was significantly upregulated in C. virginica after P. marinus challenge.
FIGURE 8. Expansion and phylogeny of FBGDC gene family in C. virginica (red font) and C. gigas (gray font). The phylogenetic tree is constructed with the maximum likelihood method showing lineage-specific expansion of FBGDC in C. virginica (red font) and C. gigas (gray font). Different color backgrounds represent different domain composition types, and the red and gray asterisks represent differentially expressed FBGDCs under P. marinus challenge in C. virginica and C. gigas, respectively.
Ahmed,
Knockdown of a galectin-1-like protein in zebrafish (Danio rerio) causes defects in skeletal muscle development.
2009, Pubmed
Ahmed,
Knockdown of a galectin-1-like protein in zebrafish (Danio rerio) causes defects in skeletal muscle development.
2009,
Pubmed
Akira,
Pathogen recognition and innate immunity.
2006,
Pubmed
Aladaileh,
Induction of phenoloxidase and other immunological activities in Sydney rock oysters challenged with microbial pathogen-associate molecular patterns.
2007,
Pubmed
Alcaide,
Molecular evolution of the toll-like receptor multigene family in birds.
2011,
Pubmed
Allam,
Early host-pathogen interactions in marine bivalves: evidence that the alveolate parasite Perkinsus marinus infects through the oyster mantle during rejection of pseudofeces.
2013,
Pubmed
Anders,
HTSeq--a Python framework to work with high-throughput sequencing data.
2015,
Pubmed
Areal,
Signatures of positive selection in Toll-like receptor (TLR) genes in mammals.
2011,
Pubmed
Burge,
Climate change influences on marine infectious diseases: implications for management and society.
2014,
Pubmed
Capella-Gutiérrez,
trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses.
2009,
Pubmed
Carnegie,
A rapid phenotype change in the pathogen Perkinsus marinus was associated with a historically significant marine disease emergence in the eastern oyster.
2021,
Pubmed
Chan,
Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions.
2021,
Pubmed
Chintala,
Comparison of in vitro-cultured and wild-type Perkinsus marinus. II. Dosing methods and host response.
2002,
Pubmed
Chiu,
Molecular Dynamics Simulations on High-Performance Reconfigurable Computing Systems.
2010,
Pubmed
Conant,
Probabilistic cross-species inference of orthologous genomic regions created by whole-genome duplication in yeast.
2008,
Pubmed
Darfour-Oduro,
Adaptive Evolution of Toll-Like Receptors (TLRs) in the Family Suidae.
2015,
Pubmed
Dodds,
Recognition events and host-pathogen co-evolution in gene-for-gene resistance to flax rust.
2009,
Pubmed
Emms,
OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy.
2015,
Pubmed
Figueras,
Genomics and immunity of the Mediterranean mussel Mytilus galloprovincialis in a changing environment.
2019,
Pubmed
Ford,
Comparison of in vitro-cultured and wild-type Perkinsus marinus. I. Pathogen virulence.
2002,
Pubmed
Fornůsková,
Contrasted evolutionary histories of two Toll-like receptors (Tlr4 and Tlr7) in wild rodents (MURINAE).
2013,
Pubmed
Gerdol,
An updated molecular basis for mussel immunity.
2015,
Pubmed
Gertz,
Composition-based statistics and translated nucleotide searches: improving the TBLASTN module of BLAST.
2006,
Pubmed
Groner,
Managing marine disease emergencies in an era of rapid change.
2016,
Pubmed
Grueber,
Episodic positive selection in the evolution of avian toll-like receptor innate immunity genes.
2014,
Pubmed
Guo,
Infectious diseases of marine molluscs and host responses as revealed by genomic tools.
2016,
Pubmed
Guo,
Immune and stress responses in oysters with insights on adaptation.
2015,
Pubmed
Gutiérrez-Rivera,
Differential expression of serine protease inhibitors 1 and 2 in Crassostrea corteziensis and C. virginica infected with Perkinsus marinus.
2015,
Pubmed
Hanington,
The primary role of fibrinogen-related proteins in invertebrates is defense, not coagulation.
2011,
Pubmed
Harvell,
Emerging marine diseases--climate links and anthropogenic factors.
1999,
Pubmed
He,
Mutation in promoter region of a serine protease inhibitor confers Perkinsus marinus resistance in the eastern oyster (Crassostrea virginica).
2012,
Pubmed
Huang,
Genomic analysis of the immune gene repertoire of amphioxus reveals extraordinary innate complexity and diversity.
2008,
Pubmed
,
Echinobase
Hughes,
Apoptosis as a host defense mechanism in Crassostrea virginica and its modulation by Perkinsus marinus.
2010,
Pubmed
Kenjo,
Cloning and characterization of novel ficolins from the solitary ascidian, Halocynthia roretzi.
2001,
Pubmed
Kim,
Noble tandem-repeat galectin of Manila clam Ruditapes philippinarum is induced upon infection with the protozoan parasite Perkinsus olseni.
2008,
Pubmed
Lan,
De novo transcriptome assembly and positive selection analysis of an individual deep-sea fish.
2018,
Pubmed
Li,
Conservation and divergence of mitochondrial apoptosis pathway in the Pacific oyster, Crassostrea gigas.
2017,
Pubmed
Matsushita,
A novel human serum lectin with collagen- and fibrinogen-like domains that functions as an opsonin.
1996,
Pubmed
Mistry,
Pfam: The protein families database in 2021.
2021,
Pubmed
Moreau,
Autophagy plays an important role in protecting Pacific oysters from OsHV-1 and Vibrio aestuarianus infections.
2015,
Pubmed
Nei,
Evolution by the birth-and-death process in multigene families of the vertebrate immune system.
1997,
Pubmed
Phillips,
Immunogenetic novelty confers a selective advantage in host-pathogen coevolution.
2018,
Pubmed
Proestou,
Variation in global transcriptomic response to Perkinsus marinus infection among eastern oyster families highlights potential mechanisms of disease resistance.
2020,
Pubmed
Qu,
Characterization of an inhibitor of apoptosis protein in Crassostrea gigas clarifies its role in apoptosis and immune defense.
2015,
Pubmed
Quevillon,
InterProScan: protein domains identifier.
2005,
Pubmed
Rabinovich,
Role of galectins in inflammatory and immunomodulatory processes.
2002,
Pubmed
Rast,
Genomic insights into the immune system of the sea urchin.
2006,
Pubmed
,
Echinobase
Reece,
Molecular epizootiology of Perkinsus marinus and P. chesapeaki infections among wild oysters and clams in Chesapeake Bay, USA.
2008,
Pubmed
Ren,
Unusual conservation of mitochondrial gene order in Crassostrea oysters: evidence for recent speciation in Asia.
2010,
Pubmed
Royle,
Stimulation of Toll-like receptor 4 by lipopolysaccharide during cellular invasion by live Salmonella typhimurium is a critical but not exclusive event leading to macrophage responses.
2003,
Pubmed
Sackton,
Comparative genomics and transcriptomics of host-pathogen interactions in insects: evolutionary insights and future directions.
2019,
Pubmed
Saha,
Rotavirus vaccines performance: dynamic interdependence of host, pathogen and environment.
2021,
Pubmed
Song,
The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia.
2021,
Pubmed
Sutherland,
Shifting white pox aetiologies affecting Acropora palmata in the Florida Keys, 1994-2014.
2016,
Pubmed
Tanguy,
Discovery of genes expressed in response to Perkinsus marinus challenge in Eastern (Crassostrea virginica) and Pacific (C. gigas) oysters.
2004,
Pubmed
Tasumi,
A galectin of unique domain organization from hemocytes of the Eastern oyster (Crassostrea virginica) is a receptor for the protistan parasite Perkinsus marinus.
2007,
Pubmed
Trapnell,
TopHat: discovering splice junctions with RNA-Seq.
2009,
Pubmed
Wang,
The oyster immunity.
2018,
Pubmed
Wang,
A C1q domain containing protein from scallop Chlamys farreri serving as pattern recognition receptor with heat-aggregated IgG binding activity.
2012,
Pubmed
Wang,
Microarray analysis of gene expression in eastern oyster (Crassostrea virginica) reveals a novel combination of antimicrobial and oxidative stress host responses after dermo (Perkinsus marinus) challenge.
2010,
Pubmed
Wang,
A tandem-repeat galectin involved in innate immune response of the pearl oyster Pinctada fucata.
2011,
Pubmed
Woolhouse,
Biological and biomedical implications of the co-evolution of pathogens and their hosts.
2002,
Pubmed
Xu,
Widespread positive selection on cetacean TLR extracellular domain.
2019,
Pubmed
Yang,
Bayes empirical bayes inference of amino acid sites under positive selection.
2005,
Pubmed
Zhang,
The oyster genome reveals stress adaptation and complexity of shell formation.
2012,
Pubmed
Zhang,
A fibrinogen-related protein from bay scallop Argopecten irradians involved in innate immunity as pattern recognition receptor.
2009,
Pubmed
Zhang,
Massive expansion and functional divergence of innate immune genes in a protostome.
2015,
Pubmed
Zhang,
Gene discovery, comparative analysis and expression profile reveal the complexity of the Crassostrea gigas apoptosis system.
2011,
Pubmed
Zhang,
Transcriptome analysis reveals a rich gene set related to innate immunity in the Eastern oyster (Crassostrea virginica).
2014,
Pubmed
Zhang,
Diversification of Ig superfamily genes in an invertebrate.
2004,
Pubmed
Zhang,
The first evidence of positive selection in peptidoglycan recognition protein (PGRP) genes of Crassostrea gigas.
2013,
Pubmed
Zong,
A novel globular C1q domain containing protein (C1qDC-7) from Crassostrea gigas acts as pattern recognition receptor with broad recognition spectrum.
2019,
Pubmed