Transcriptome analysis reveals a rich gene set related to innate immunity in the Eastern oyster (Crassostrea virginica)

Transcriptome analysis and microsatellite discovery in the blunt snout bream (Megalobrama amblycephala) after challenge with Aeromonas hydrophila

The blunt snout bream, Megalobrama amblycephala, is a herbivorous freshwater fish species native to China and a major aquaculture species in Chinese freshwater polyculture systems. In recent years, the bacterium Aeromonas hydrophila has been reported to be its pathogen causing great losses of farmed fish. To understand the immune response of the blunt snout bream to A. hydrophila infection, we used the Solexa/Illumina technology to analyze the transcriptomic profile after artificial bacterial infection. Two nonnormalized cDNA libraries were synthesized from tissues collected from control blunt snout bream or those injected with A. hydrophila. After assembly, 155,052 unigenes (average length 692.8 bp) were isolated. All unigenes were annotated using BLASTX relative to several public databases: the National Center for Biotechnology Information nonreduntant (Nr) database, SwissProt, Eukaryotic Orthologous Groups of proteins (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO). The sequence similarity (86%) of the assembled unigenes was to zebrafish based on the Nr database. A number of unigenes (n = 30,482) were assigned to three GO categories: biological processes (25,242 unigenes), molecular functions (26,096 unigenes), and cellular components (22,778 unigenes). 20,909 unigenes were classified into 25 KOG categories and 28,744 unigenes were assigned into 315 specific signaling pathways. In total, 238 significantly differentially expressed unigenes (mapped to 125 genes) were identified: 101 upregulated genes and 24 downregulated genes. Another 303 unigenes were mapped to unknown or novel genes. Among the known expressed genes identified, 53 were immune-related genes and were distributed in 71 signaling pathways. The expression patterns of selected up- and downregulated genes from the control and injected groups were determined with reverse transcription-quantitative PCR (RT-qPCR). Microsatellites (n = 10,877), including di-to pentanucleotide repeat motifs, were also identified in the blunt snout bream transcriptome profiles. This study extends our understanding of the immune defense mechanisms of the blunt snout bream against A. hydrophila and provides useful data for further studies of the immunogenetics of this species.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

Immune responses to infectious diseases in bivalves

Many species of bivalve mollusks (phylum Mollusca, class Bivalvia) are important in fisheries and aquaculture, whilst others are critical to ecosystem structure and function. These crucial roles mean that considerable attention has been paid to the immune responses of bivalves such as oysters, clams and mussels against infectious diseases that can threaten the viability of entire populations. As with many invertebrates, bivalves have a comprehensive repertoire of immune cells, genes and proteins. Hemocytes represent the backbone of the bivalve immune system. However, it is clear that mucosal tissues at the interface with the environment also play a critical role in host defense. Bivalve immune cells express a range of pattern recognition receptors and are highly responsive to the recognition of microbe-associated molecular patterns. Their responses to infection include chemotaxis, phagolysosomal activity, encapsulation, complex intracellular signaling and transcriptional activity, apoptosis, and the induction of anti-viral states. Bivalves also express a range of inducible extracellular recognition and effector proteins, such as lectins, peptidoglycan-recognition proteins, thioester bearing proteins, lipopolysaccharide and β1,3-glucan-binding proteins, fibrinogen-related proteins (FREPs) and antimicrobial proteins. The identification of FREPs and other highly diversified gene families in bivalves leaves open the possibility that some of their responses to infection may involve a high degree of pathogen specificity and immune priming. The current review article provides a comprehensive, but not exhaustive, description of these factors and how they are regulated by infectious agents. It concludes that one of the remaining challenges is to use new "omics" technologies to understand how this diverse array of factors is integrated and controlled during infection.

Domain repeats related to innate immunity in the South African abalone, Haliotis midae

Molluscs predominately use the cellular defence system as the primary mechanism of defence against pathogenic infection, in which haemocytes play a pivotal role. Haliotis midae is a commercially important South African species that it is susceptible to bacterial pathogens, fungal and yeast infections in the farming environment. The current study aims to enrich the current knowledge regarding H. midae innate immunity by investigating the presence and evolution of domain repeats. The bioinformatics approach used in this study, detected five repeat families in the H. midae transcriptome. These repeats families include mixed alpha and beta (leucine-rich and ankyrin), spectrin repeats, beta-propellers (WD40) and alfa-structure repeat (TPR-like). The expansion of key gene families related to host defence may be important to abalone adaptation to life in a pathogen-rich environment.

Functional marker detection and analysis on a comprehensive transcriptome of large yellow croaker by next generation sequencing

Large yellow croaker (Larimichthys crocea) is an important economic fish in China and Eastern Asia. Because of the exhaustive fishing and overdense aquaculture, the wild population and the mariculture of the species are facing serious challenges on germplasm degeneration and susceptibility to infectious disease agents. However, a comprehensive transcriptome from multi-tissues of the species has not been reported and functional molecular markers have not yet been detected and analyzed. In this work, we applied RNA-seq with the Illumina Hiseq2000 platform for a multi-tissue sample of large yellow croaker and assembled the transcriptome into 88,103 transcripts. Of them, 52,782 transcripts have been successfully annotated by nt/nr, InterPro, GO and KEGG database. Comparing with public fish proteins, we have found that 34,576 protein coding transcripts are shared in large yellow croaker with zebrafish, medaka, pufferfish, and stickleback. For functional markers, we have discovered 1,276 polymorphic SSRs and 261, 000 SNPs. The functional impact analysis of SNPs showed that the majority (~75%) of small variants cause synonymous mutations in proteins, followed by variations in 3' UTR region. The functional enrichment analysis illuminated that transcripts involved in DNA bindings, enzyme activities, and signal pathways prominently exhibit less single-nucleotide variants but genes for the constituent of the muscular tissue, the cytoskeleton, and the immunity system contain more frequent SNP mutations, which may reflect the structural and functional selections of the translated proteins. This is the first work for the high-throughput detection and analysis of functional polymorphic SSR and SNP markers in a comprehensive transcriptome of large yellow croaker. Our study provides valuable transcript sequence and functional marker resources for the quantitative trait locus (QTL) identification and molecular selection of the species in the research community.

Highly diverse fibrinogen-related proteins in the Pacific oyster Crassostrea gigas

Fibrinogen-related proteins (FREPs) are a family of proteins with high sequence diversity, and they play crucial roles in invertebrate immune response. However, few studies have characterized this diversity at the whole-genome level. In the present study, approximately 190 predicted FREPs with more than 200 fibrinogen-like (FBG) domains were identified in the genome of the Pacific oyster (Crassostrea gigas), suggesting a historical expansion of this protein family. A sequence analysis showed high numbers of polymorphisms in C. gigas FREP (CgFREP) genes, which may contribute to the versatile immune function of FREPs. A phylogenetic analysis of molluscan FREP sequences indicated lineage-specific duplication of these genes in C. gigas. Additionally, several CgFREP mRNAs were highly expressed in the gills, digestive glands, and hemocytes. Taken together, these findings will help elucidate FREP immune function and facilitate studies of the functional validation of this gene family.

The genome of the Pacific oyster Crassostrea gigas brings new insights on the massive expansion of the C1q gene family in Bivalvia

C1q domain-containing (C1qDC) proteins are regarded as important players in the innate immunity of bivalve mollusks and other invertebrates and their highly adaptive binding properties indicate them as efficient pathogen recognition molecules. Although experimental studies support this view, the molecular data available at the present time are not sufficient to fully explain the great molecular diversification of this family, present in bivalves with hundreds of C1q coding genes. Taking advantage of the fully sequenced genome of the Pacific oyster Crassostrea gigas and more than 100 transcriptomic datasets, we: (i) re-annotated the oyster C1qDC loci, thus identifying the correct genomic organization of 337 C1qDC genes, (ii) explored the expression pattern of oyster C1qDC genes in diverse developmental stages and adult tissues of unchallenged and experimentally treated animals; (iii) investigated the expansion of the C1qDC gene family in all major bivalve subclasses. Overall, we provide a broad description of the functionally relevant features of oyster C1qDC genes, their comparative expression levels and new evidence confirming that a gene family expansion event has occurred during the course of Bivalve evolution, leading to the diversification of hundreds of different C1qDC genes in both the Pteriomorphia and Heterodonta subclasses.

Transcriptome sequencing and characterization of Japanese scallop Patinopecten yessoensis from different shell color lines

Shell color is an important trait that is used in breeding the Japanese scallop Patinopecten yessoensis, the most economically important scallop species in China. We constructed four transcriptome libraries from different shell color lines of P. yessoensis: the left and right shell mantles of ordinary strains of P. yessoensis and the left shell mantles of the ‘Ivory’ and ‘Maple’ strains. These four libraries were paired-end sequenced using the Illumina HiSeq 2000 platform and contained 54,802,692 sequences, 40,798,962 sequences, 74,019,262 sequences, and 44,466,166 sequences, respectively. A total of 214,087,082 expressed sequence tags were assembled into 73,522 unigenes with an average size of 1,163 bp. When the data were compared against the public Nr and Swiss-Prot databases using BlastX, nearly 30.55% (22,458) of the unigenes were significantly matched to known unique proteins. Gene Ontology annotation and pathway mapping analysis using the Kyoto Encyclopedia of Genes and Genomes categorized unigenes according to their diverse biological functions and processes and identified candidate genes that were potentially involved in growth, pigmentation, metal transcription, and immunity. Expression profile analysis was performed on all four libraries and many differentially expressed genes were identified. In addition, 5,772 simple sequence repeats were obtained from the P. yessoensis transcriptomes, and 464,197, 395,646, and 310,649 single nucleotide polymorphisms were revealed in the ordinary strains, the ‘Ivory’ strain, and the ‘Maple’ strain, respectively. These results provide valuable information for future genomic studies on P. yessoensis and improve our understanding of the molecular mechanisms involved in the growth, immunity, shell coloring, and shell biomineralization of this species. These resources also may be used in a variety of applications, such as trait mapping, marker-assisted breeding, studies of population genetics and genomics, and work on functional genomics.

Immune-directed support of rich microbial communities in the gut has ancient roots

The animal gut serves as a primary location for the complex host-microbe interplay that is essential for homeostasis and may also reflect the types of ancient selective pressures that spawned the emergence of immunity in metazoans. In this review, we present a phylogenetic survey of gut host-microbe interactions and suggest that host defense systems arose not only to protect tissue directly from pathogenic attack but also to actively support growth of specific communities of mutualists. This functional dichotomy resulted in the evolution of immune systems much more tuned for harmonious existence with microbes than previously thought, existing as dynamic but primarily cooperative entities in the present day. We further present the protochordate Ciona intestinalis as a promising model for studying gut host-bacterial dialogue. The taxonomic position, gut physiology and experimental tractability of Ciona offer unique advantages in dissecting host-microbe interplay and can complement studies in other model systems.

Identification of potential general markers of disease resistance in American oysters, Crassostrea virginica through gene expression studies

Several diseases have a significant impact on American oyster populations in the Atlantic coasts of North America. Knowledge about the responses of oysters to pathogenic challenge could help in identifying potential markers of disease resistance and biomarkers of the health status of an oyster population. A previous analysis of the transcriptome of resistant and susceptible American oysters in response to challenge with the bacterial pathogen Roseovarius crassostreae, as well as sequencing of suppression subtractive hybridization libraries from oysters challenged with the protozoan parasite Perkinsus marinus, provided a list of genes potentially involved in disease resistance or susceptibility. We investigated the patterns of inducible gene expression of several of these genes in response to experimental challenge with the oyster pathogens R. crassostreae, Vibrio tubiashii, and P. marinus. Oysters showing differential susceptibility to R. crassostreae demonstrated differential patterns of expression of genes coding for immune (serine protease inhibitor-1, SPI1) and stress-related (heat shock protein 70, HSP70; arginine kinase) proteins 30 days after challenge with this bacterial pathogen. Differential patterns of expression of immune (spi1, galectin and a matrix metalloproteinase) and stress-related (hsp70, histone H4, and arginine kinase) genes was observed in hemocytes from adult oysters challenged with P. marinus, but not with V. tubiashii. While levels of spi1 expression in hemocytes collected 8 and 21 days after P. marinus challenge were negatively correlated with parasite load in oysters tissues at the end of the challenge (62 days), levels of expression of hsp70 in hemocytes collected 1-day after challenge were positively correlated with oyster parasite load at 62 days. Our results confirm previous research on the role of serine protease inhibitor-1 in immunity and disease resistance in oysters. They also suggest that HSP70 and histone H4 could be used as a markers of health status or disease susceptibility in oysters.

Repurposed transcriptomic data facilitate discovery of innate immunity toll-like receptor (TLR) Genes across Lophotrochozoa

The growing volume of genomic data from across life represents opportunities for deriving valuable biological information from data that were initially collected for another purpose. Here, we use transcriptomes collected for phylogenomic studies to search for toll-like receptor (TLR) genes in poorly sampled lophotrochozoan clades (Annelida, Mollusca, Brachiopoda, Phoronida, and Entoprocta) and one ecdysozoan clade (Priapulida). TLR genes are involved in innate immunity across animals by recognizing potential microbial infection. They have an extracellular leucine-rich repeat (LRR) domain connected to a transmembrane domain and an intracellular toll/interleukin-1 receptor (TIR) domain. Consequently, these genes are important in initiating a signaling pathway to trigger defense. We found at least one TLR ortholog in all but two taxa examined, suggesting that a broad array of lophotrochozoans may have innate immune systems similar to those observed in vertebrates and arthropods. Comparison to the SMART database confirmed the presence of both the LRR and the TIR protein motifs characteristic of TLR genes. Because we looked at only one transcriptome per species, discovery of TLR genes was limited for most taxa. However, several TRL-like genes that vary in the number and placement of LRR domains were found in phoronids. Additionally, several contigs contained LRR domains but lacked TIR domains, suggesting they were not TLRs. Many of these LRR-containing contigs had other domains (e.g., immunoglobin) and are likely involved in innate immunity.

Transcriptome of American oysters, Crassostrea virginica, in response to bacterial challenge: insights into potential mechanisms of disease resistance

The American oyster Crassostrea virginica, an ecologically and economically important estuarine organism, can suffer high mortalities in areas in the Northeast United States due to Roseovarius Oyster Disease (ROD), caused by the gram-negative bacterial pathogen Roseovarius crassostreae. The goals of this research were to provide insights into: 1) the responses of American oysters to R. crassostreae, and 2) potential mechanisms of resistance or susceptibility to ROD. The responses of oysters to bacterial challenge were characterized by exposing oysters from ROD-resistant and susceptible families to R. crassostreae, followed by high-throughput sequencing of cDNA samples from various timepoints after disease challenge. Sequence data was assembled into a reference transcriptome and analyzed through differential gene expression and functional enrichment to uncover genes and processes potentially involved in responses to ROD in the American oyster. While susceptible oysters experienced constant levels of mortality when challenged with R. crassostreae, resistant oysters showed levels of mortality similar to non-challenged oysters. Oysters exposed to R. crassostreae showed differential expression of transcripts involved in immune recognition, signaling, protease inhibition, detoxification, and apoptosis. Transcripts involved in metabolism were enriched in susceptible oysters, suggesting that bacterial infection places a large metabolic demand on these oysters. Transcripts differentially expressed in resistant oysters in response to infection included the immune modulators IL-17 and arginase, as well as several genes involved in extracellular matrix remodeling. The identification of potential genes and processes responsible for defense against R. crassostreae in the American oyster provides insights into potential mechanisms of disease resistance.

Immune response of apoptosis-related cysteine peptidases from the red abalone Haliotis rufescens (HrCas8 and HrCas3): molecular characterization and transcription expression

Caspases play an important role in the different stages of programmed cell death, or apoptosis, which has been related to the immune response in multicellular organisms. The present study characterized an initiator caspase (HrCas8) and an effector caspase (HrCas3) from the red abalone Haliotis rufescens using the RACE method and qPCR analysis. HrCas8 showed a complete sequence of 2529 base pairs (bp) with an open-reading frame (ORF) of 1911 bp, a 5'UTR of 201 bp, and a 3'UTR of 417 bp. The estimated molecular mass for the 636 amino acids from HrCas8 was 71.5 kDa with an isoelectric point of 6.2. The HrCas8 sequence had two death-effector domains (DEDs) and the subunits p20 and p10, in addition to an active site characteristic of cysteine proteins. Meanwhile, the effector caspase HrCas3 showed a complete sequence of 1404 bp, a 5'UTR of 82 bp, and a 3'UTR of 574 bp. The ORF of this caspase had 747 bp that coded for 248 residues. Moreover, the predicted molecular mass of HrCas3 was 29.4 kDa; the theoretical isoelectric point was 5.70, and the sequence evidenced a conserved caspase recruitment domain (CARD). The distribution of the caspases in distinct tissues revealed that HrCas8 was principally expressed in the hemolymph, while HrCas3 had a higher expression in the gills. A basal level of expression was found for both caspases in muscle tissue. The immune response of caspases in H. rufescens was evaluated through an injection of Vibrio anguillarum. The results showed an increase in the transcription of HrCas8 post-challenge, as well as an activation of HrCas3, which together suggest the initiation of apoptosis as a response to bacterial infection in H. rufescens.

Transcriptomic analysis of candidate osmoregulatory genes in the eastern oyster Crassostrea virginica

BACKGROUND

The eastern oyster, Crassostrea virginica, is a euryhaline species that can thrive across a wide range of salinities (5-35). As with all estuarine species, individual oysters must be able to regulate their osmotic balance in response to constant temporal variation in salinity. At the population level, recurrent viability selection may be an additional mechanism shaping adaptive osmoregulatory phenotypes at the margins of oyster salinity tolerance. To identify candidate genes for osmoregulation, we sequenced, assembled, and annotated the transcriptome of wild juvenile eastern oysters from 'high' and 'low' salinity regimes. Annotations and candidates were mostly based on the Pacific oyster (Crassostrea gigas) genome sequence so osmoregulatory relevance in C. virginica was explored by testing functional enrichment of genes showing spatially discrete patterns of expression and by quantifying coding sequence divergence.

RESULTS

The assembly of sequence reads and permissive clustering of potentially oversplit alleles resulted in 98,729 reftigs (contigs and singletons). Of these, 50,736 were annotated with 9,307 belonging to a set of candidate osmoregulatory genes identified from the C. gigas genome. A total of 218,777 SNPs (0.0185 SNPs/bp) were identified in annotated reftigs of C. virginica. Amino acid divergence between translations of C. virginica annotated reftigs and C. gigas coding sequence averaged 23.2 % with an average dN/dS ratio of 0.074, suggesting purifying selection on protein sequences. The high and low salinity source oysters each expressed a subset of genes unique to that group, and the functions for these annotated genes were consistent with known molecular mechanisms for osmotic regulation in molluscs.

CONCLUSIONS

Most of the osmoregulatory gene candidates experimentally identified in C. gigas are present in this C. virginica transcriptome. In general these congeners show coding sequence divergence too high to make the C. gigas genome a useful reference for C. virginica bioinformatics. However, strong purifying selection is characteristic of the osmoregulatory candidates so functional annotations are likely to correspond. An initial examination of C. virginica presence/absence expression patterns across the salinity gradient in a single estuary suggests that many of these candidates have expression patterns that co-vary with salinity, consistent with osmoregulatory function in C. virginica.