Identification of genes involved in immune response, microsatellite, and SNP markers from expressed sequence tags generated from hemocytes of freshwater pearl mussel (Hyriopsis cumingii)

Single Nucleotide Polymorphism Markers with Applications in Conservation and Exploitation of Aquatic Natural Populations

Simple Summary

In recent decades, societies, states and local authorities have become increasingly aware that for effective long-term management and protection of aquatic ecosystems and populations, it is necessary to take into account the genetic changes occurring in these populations. One type of high-resolution molecular marker suitable for studying the neutral and adaptive genetic diversity of populations is single nucleotide polymorphism (SNP). This review is an attempt to show the benefits of using SNPs to recognize natural populations of aquatic animals and detect the threats to them from accidentally or intentionally released farm animals, fishery and global climate changes. It is postulated that conservation actions should protect not only pristine natural populations that are endangered or overfished, but also protect populations of non-threatened species from unnecessarily released semi-domesticated animals. The enhancement of natural populations with farmed material usually reduces their genetic diversity. Experimental size-selective catches of artificially created populations have caused evolutionary changes in the life cycles of fishes. However, fishery-induced evolution in natural populations is difficult to observe. The negative measurable effects on populations can be expected when the number of breeding individuals is reduced below 100, which occurs very rarely in the sea and more often in fragmented freshwater streams, ponds and seasonal rivers.

Abstract

An increasing number of aquatic species have been studied for genetic polymorphism, which extends the knowledge on their natural populations. One type of high-resolution molecular marker suitable for studying the genetic diversity of large numbers of individuals is single nucleotide polymorphism (SNP). This review is an attempt to show the range of applications of SNPs in studies of natural populations of aquatic animals. In recent years, SNPs have been used in the genetic analysis of wild and enhanced fish and invertebrate populations in natural habitats, exploited migratory species in the oceans, migratory anadromous and freshwater fish and demersal species. SNPs have been used for the identification of species and their hybrids in natural environments, to study the genetic consequences of restocking for conservation purposes and the negative effects on natural populations of fish accidentally escaping from culture. SNPs are very useful for identifying genomic regions correlated with phenotypic variants relevant for wildlife protection, management and aquaculture. Experimental size-selective catches of populations created in tanks have caused evolutionary changes in life cycles of fishes. The research results have been discussed to clarify whether the fish populations in natural conditions can undergo changes due to selective harvesting targeting the fastest-growing fishes.

The immune system of the freshwater zebra mussel, Dreissena polymorpha, decrypted by proteogenomics of hemocytes and plasma compartments

The immune system of bivalves is of great interest since it reflects the health status of these organisms during stressful conditions. While immune molecular responses are well documented for marine bivalves, few information is available for continental bivalves such as the zebra mussel, Dreissena polymorpha. A proteogenomic approach was conducted on both hemocytes and plasma to identified immune proteins of this non-model species. Combining transcriptomic sequences with mass spectrometry data acquired on proteins is a relevant strategy since 3020 proteins were identified, representing the largest protein inventory for this sentinel organism. Functional annotation and gene ontology (GO) analysis performed on the identified proteins described the main molecular players of hemocytes and plasma in immunity. GO analysis highlights the complementary immune functions of these two compartments in the management of micro-organisms. Functional annotation revealed new mechanisms in the immune defence of the zebra mussel. Proteins rarely observed in the hemolymph of bivalves were pinpointed such as natterin-like and thaumatin-like proteins. Furthermore, the high abundance of complement-related proteins observed in plasma suggested a strong implication of the complement system in the immune defence of D. polymorpha. This work brings a better understanding of the molecular mechanisms involved in zebra mussel immunity. SIGNIFICANCE: Although the molecular mechanisms of marine bivalves are widely investigated, little information is known for continental bivalves. Moreover, few proteomic studies described the complementarity of both hemolymphatic compartments (cellular and plasmatic) in the immune defence of invertebrates. The recent proteogenomics concept made it possible to discover proteins in non-model organisms. Here, we propose a proteogenomic strategy with the zebra mussel, a key sentinel species for biomonitoring of freshwater, to identify and describe the molecular actors involved in the immune system in both hemocytes and plasma compartments. More widely, this study provided new insight into bivalve immunity.

Transcriptome Analysis Provides Insights into the Markers of Resting and LPS-Activated Macrophages in Grass Carp (Ctenopharyngodon idella)

Macrophages are very versatile immune cells, with the characteristics of a proinflammatory phenotype in response to pathogen-associated molecular patterns. However, the specific activation marker genes of macrophages have not been systematically investigated in teleosts. In this work, leukocytes (WBC) were isolated using the Percoll gradient method. Macrophages were enriched by the adherent culture of WBC, then stimulated with lipopolysaccharide (LPS). Macrophages were identified by morphological features, functional activity and authorized cytokine expression. Subsequently, we collected samples, constructed and sequenced transcriptomic libraries including WBC, resting macrophage (Mø) and activated macrophage (M(LPS)) groups. We gained a total of 20.36 Gb of clean data including 149.24 million reads with an average length of 146 bp. Transcriptome analysis showed 708 differential genes between WBC and Mø, 83 differentially expressed genes between Mø and M(LPS). Combined with RT-qPCR, we proposed that four novel cell surface marker genes (CD22-like, CD63, CD48 and CD276) and two chemokines (CXCL-like and CCL39.3) would be emerging potential marker genes of macrophage in grass carp. Furthermore, CD69, CD180, CD27, XCL32a.2 and CXCL8a genes can be used as marker genes to confirm whether macrophages are activated. Transcriptome profiling reveals novel molecules associated with macrophages in C. Idella, which may represent a potential target for macrophages activation.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.

Proteomic characterization of mucosal secretions in the eastern oyster, Crassostrea virginica

The soft body surface of marine invertebrates is covered by a layer of mucus, a slippery gel secreted by mucocytes lining epithelia. The functions of this gel are diverse including locomotion, cleansing, food particles processing and defense against physicochemical injuries and infectious agents. In oysters, mucus covering pallial organs has been demonstrated to have a major importance in the processing of food particles and in the interactions with waterborne pathogens. Given the limited information available on mucus in bivalves and the apparent wide spectra of activity of bioactive molecules present in this matrix, the characterization of these mucosal secretions has become a research priority. In this study, mucus was separately collected from the mantle, gills and labial palps of the eastern oyster (Crassostrea virginica) and analyzed by liquid chromatography and tandem mass spectrometry. Results showed the presence of a wide variety of molecules involved in host-microbe interactions, including putative adhesion molecules (e.g. c-type lectins) confirming that transcripts previously identified in epithelial cells are translated into proteins secreted in mucus. Mucus composition was different among samples collected from different organs. These results generate a reference map for C. virginica pallial mucus to better characterize the various physiological functions of mucosal secretions.

Microsatellite analysis of genetic diversity and population structure of freshwater mussel (Lamprotula leai)

Lamprotula leai is one of the most commercially important freshwater pearl mussels in China, but there is limited data on its genetic diversity and population structure. In the present study, 119 individuals from four major geographical populations were investigated using 15 microsatellite loci identified via cross-species amplification. A total of 114 alleles were detected, with an average of 7.6 alleles per locus (range: 2 to 21). Among the four stocks, those from Hung-tse Lake and Poyang Lake had the lowest (0.412) and highest (0.455) observed heterozygosity respectively. The polymorphism information content (PIC) ranged from 0.374 to 0.927 (mean: 0.907). AMOVA showed that 12.56% and 44.68% genetic variances were among populations and within individuals, respectively. Pairwise Fst ranged from 0.073 to 0.146, indicating medium genetic differentiation among the populations. In aggregate, our results suggest that inbreeding is a crucial factor accounting for deviations from Hardy-Weinberg equilibrium at 12 loci. Moreover, the genetic distance among four stocks ranged from 0.192 to 0.890. Poyang Lake and Hung-tse Lake were clustered together, joined with Dongting Lake and Anqing Lake. Given that specimens from Hung-tse Lake showed the highest average allele richness, expected heterozygosity and PIC, this location may be the source of the highest quality germplasm resources and the stock from this area may be the best for future breeding efforts.

Characterization of the LECT2 gene and its associations with resistance to the big belly disease in Asian seabass

Leukocyte cell-derived chemotaxin-2 (LECT2) is an important protein of the innate immune system for the defense against bacterial infection. We cloned and characterized the LECT2 gene from Asian seabass (Lates calcarifer). Its complete cDNA consisted of an open reading frame of 459 bp encoding a protein of 152 amino acids. The genomic DNA sequence of this gene consists of four exons and three introns. Quantitative real-time PCR revealed that the LECT2 gene was expressed predominantly in liver while its expression was moderate in spleen and heart, and weak in other tissues. The LECT2 transcript was up-regulated in the kidney, spleen and liver in response to a challenge with a pathogenic bacterium Vibrio harveyi. In addition, we identified three single nucleotide polymorphisms (SNPs) in the LECT2 gene, and found significant associations between these polymorphisms and resistance to the big belly disease. These results suggest that the LECT2 gene play an important role in resistance to bacterial pathogens in fish. The SNP markers in the gene associated with the resistance to bacterial pathogens may facilitate selecting Asian seabass resistant to bacterial diseases.

Effects of benzo(a)pyrene on differentially expressed genes and haemocyte parameters of the clam Venerupis philippinarum

In this study a suppression subtractive hybridisation method was employed to identify differentially expressed genes of the clam Venerupis philippinarum exposed to benzo(a)pyrene (BaP). Nineteen known transcripts and seven predicted proteins were found from the subtractive cDNA library of the clam, which could provide more sequence information for further study. Seven of the differentially expressed genes were selected for mRNA expression analysis. Real-time PCR analysis revealed that the expression level of the selected cDNAs of clams was up-regulated to varying degrees by different concentration of BaP. They are suggested as potential molecular biomarkers for polycyclic aromatic hydrocarbons (PAHs) pollution monitoring in aquatic ecosystems. In addition, haemocyte parameters were also measured, and a decrease of total haemocyte counts and suppression of antibacterial and bacteriolytic activities were detected in BaP-stressed clams. We suggest that the modulation of the expression of the selected genes caused by PAHs probably leads to the disturbance of the immune defense of the clam. Meanwhile, the adverse effects of PAHs on haemocyte parameters caused the suppression of the immune defense and susceptibility to infectious diseases. Therefore, it is inferred that PAHs pollutants could interact with components of the immune system and interferes with defense functions of the clam V. philippinarum.

Transcriptome-wide single nucleotide polymorphisms (SNPs) for abalone (Haliotis midae): validation and application using GoldenGate medium-throughput genotyping assays

Haliotis midae is one of the most valuable commercial abalone species in the world, but is highly vulnerable, due to exploitation, habitat destruction and predation. In order to preserve wild and cultured stocks, genetic management and improvement of the species has become crucial. Fundamental to this is the availability and employment of molecular markers, such as microsatellites and single nucleotide (SNPs). Transcriptome sequences generated through sequencing-by-synthesis technology were utilized for the in vitro and in silico identification of 505 putative SNPs from a total of 316 selected contigs. A subset of 234 SNPs were further validated and characterized in wild and cultured abalone using two Illumina GoldenGate genotyping assays. Combined with VeraCode technology, this genotyping platform yielded a 65%–69% conversion rate (percentage polymorphic markers) with a global genotyping success rate of 76%–85% and provided a viable means for validating SNP markers in a non-model species. The utility of 31 of the validated SNPs in population structure analysis was confirmed, while a large number of SNPs (174) were shown to be informative and are, thus, good candidates for linkage map construction. The non-synonymous SNPs (50) located in coding regions of genes that showed similarities with known proteins will also be useful for genetic applications, such as the marker-assisted selection of genes of relevance to abalone aquaculture.

Significant association between SNPs in the superoxide dismutase 3, extracellular (SOD3) gene and resistance to Aeromonas hydrophila in the freshwater mussel Hyriopsis cumingii

Extracellular superoxide dismutase (SOD3) is a major antioxidant enzyme that protects organs from damage by reactive oxygen species (ROS). In this study, the SOD3 gene was identified and characterized from the freshwater mussel Hyriopsis cumingii (Hc-SOD3). The cDNA sequence consists of 763 bp, encoding a protein of 208 amino acids. The amino acid sequence possesses two CuZnSOD signature sequences, and amino acids required for binding of Cu (His-93, -95, -110 and -169) and Zn (His-110, -118, -129 and Asp-132) were conserved in Hc-SOD3. The Hc-SOD3 genomic sequence was 9165 bp in length, containing four exons and three introns. Eighteen single nucleotide polymorphisms were detected in the Hc-SOD3 gene from resistant stock (RS) and susceptible stock (SS) of H. cumingii to Aeromonas hydrophila. The genotype and allele distribution were examined in resistant and susceptible stocks. Among them, a C/G substitution at the g.7994C>G locus and G/C substitution at the g.8087G>C locus were significantly associated with resistance/susceptibility of H. cumingii to A. hydrophila, both in genotype (P = 0.017, P = 0.004 respectively) and allele frequency (P = 0.021, P = 0.006 respectively). Linkage disequilibrium analysis revealed that g.7994C>G, g.8001A>G, g.8035G>A, g.8087G>C and g.8191T>A were in linkage disequilibrium. The results suggest that the two polymorphic loci, g.7994C>G and g.8087G>C, could be potential genetic markers for future molecular selection of strains that are resistant to diseases.

Association between SNPs in interferon regulatory factor 2 (IRF-2) gene and resistance to Aeromonas hydrophila in freshwater mussel Hyriopsis cumingii

Interferon regulatory factor 2 (IRF-2) is a multi-functional transcription factor in the IRF family exhibiting both transcriptional activating and repressing activities. In this study, an IRF-2 gene (HcIRF-2) from Hyriopsis cumingii was identified and characterized. The cDNA sequence consisted of 2688 bp, encoding a 329 amino acid-protein. The amino acid sequence had a highly conserved N-terminal DBD structure, containing characteristic repeats of six tryptophan residues. The 5'-flanking region contained several transcription regulation elements such as AP1, CdxA, HSF, NIT2 and HNF-3b. Nine SNPs were obtained through direct sequencing of HcIRF-2 from resistant and susceptible stock. Only +2365T/C SNP was significantly associated with resistance/susceptibility of H. cumingii to Aeromonas hydrophila both in genotype (P = 0.021) and allele (P = 0.006) analysis. The SNPs +2248T/C and +2365T/C were in high linkage disequilibrium, and haplotype analysis revealed that haplotype TT frequency in the resistant group was significantly higher than in the susceptible group. The mortality in +2248CC genotype individuals was significantly higher than in CT and TT genotype individuals. These results indicated that haplotype TT and genotype +2248CT and +2248GT individuals were resistant to A. hydrophila, which could make them potential markers in selective breeding of H. cumingii.

Comparative analysis of the transcriptome in tissues secreting purple and white nacre in the pearl mussel Hyriopsis cumingii

The triangle sail mussel Hyriopsis cumingii (Lea) is the most important mussel species used for commercial freshwater pearl production in China. Mussel color is an important indicator of pearl quality. To identify genes involved in the nacre coloring, we conducted RNA-seq and obtained 541,268 sequences (298 bp average size) and 440,034 sequences (293 bp average size) in secreting purple and white nacre libraries (P- and W-libraries), respectively. The 981,302 Expressed Sequence Tags (ESTs) were assembled into 47,812 contigs and 289,386 singletons. In BLASTP searches of the deduced protein, 22,495 were proteins with functional annotations. Thirty-three genes involved in pearl or shell formation were identified. Digital expression analysis identified a total of 358 differentially expressed genes, and 137 genes in the P-library and 221 genes in the W-library showed significantly higher expression. Furthermore, a set of SSR motifs and SNPs between the two samples was identified from the ESTs, which provided the markers for genetic linkage, QTL analysis and future breeding. These EST sequences provided valuable information to further understand the molecular mechanisms involved in the formation, color determination and evolution of the pearl or shell.

Six defensins from the triangle-shell pearl mussel Hyriopsis cumingii

Antimicrobial peptides (AMPs) are the first line of defense of invertebrates against invading pathogens. Defensins, unique AMPs, have a cysteine-stabilized α-helix and β-sheet (CSαβ) motif. In invertebrates, defensins have been reported in arthropods and mussels. Recently, six defensins were identified from Hyriopsis cumingii for the first time, and were designated as HcDef1, HcDef2, HcDef3, HcDef4, HcDef5, and HcDef6. HcDef1 and HcDef2 encode a protein containing 61 and 60 amino acids, respectively. HcDef3, HcDef4, and HcDef6 have 65 amino acids each. HcDef5 is longer than the other five defensins, comprising 83 amino acids. HcDef3 and HcDef4 have three pairs of disulfide bonds. HcDef1, HcDef5, and HcDef6 are exceptions; each has four pairs of disulfide bonds. Evolutionary analysis revealed that only purifying selection and no positive selection could be detected in defensin genes; purifying selection might be the major evolutionary driving force in the evolution of defensin genes. The present study reveals for the first time that the defensins from H. cumingii are diverse and phylogenetic analysis showed that these 6 defensins from H. cumingii were clustered into one group. Reverse transcriptase polymerase chain reaction (RT-PCR) analysis showed that HcDef1-HcDef4 could be detected in the hepatopancreas and gills whereas HcDef5-HcDef6 could only be detected in gills. In addition, the expression levels of HcDef2, HcDef3, and HcDef5 in H. cumingii with pearls were higher than that in H. cumingii without pearls. Quantitative RT-PCR analysis showed that HcDef1, HcDef2, HcDef3, and HcDef5 were downregulated by Vibrio anguillarum challenge whereas HcDef4 and HcDef6 were upregulated under Vibrio challenge. Our results suggest the roles of defensins in the innate immunity of H. cumingii.

Development of type I genetic markers from expressed sequence tags in highly polymorphic species

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: 1) priming site polymorphism in the template leads to inferior or erratic amplification; 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions; and 3) at least occasionally, a PCR primer straddles an exon-intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3' end of coding sequence and 3' UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.

Identification of immune genes of the Agamaki clam (Sinonovacula constricta) by sequencing and bioinformatic analysis of ESTs

The Agamaki clam (Sinonovacula constricta) is an economically important shellfish in Asia. However, genomic research on this species is still in its infancy, and genomic resources are largely unavailable. The objective of this study was to generate expressed sequence tags (ESTs) from a normalized liver complementary DNA library and to identify genes that function in immune defense. A total of 5,296 ESTs were sequenced, from which 540 contigs and 3,473 singletons were identified. BLAST homology analysis indicated that only 20.7% of these ESTs were homologues of known genes while the remaining 79.3% appeared to be novel sequences. Based on sequence similarities, 43 putative immune genes were identified that encode proteases and protease inhibitors, adhesive proteins, stress proteins, lysosomal enzymes, and signal transduction regulators. Our study thus provides both a large collection of novel transcripts and a detailed annotation of immune genes for an important bivalve species.

Generation and analysis of 10 000 ESTs from the half-smooth tongue sole Cynoglossus semilaevis and identification of microsatellite and SNP markers

Three normalized cDNA libraries were constructed, two of which were constructed from reproductive tissues ovary and testis, and the other one from pooled immune tissues including head kidney, intestine, liver and spleen. A total of 10 542 clones were sequenced generating 10 128 expressed sequence tags (ESTs). Cluster analysis indicated a total of 5808 unique sequences including 1712 contigs and 4096 singletons. A total of 4249 (73%) of the unique ESTs had significant hits to the non-redundant protein database, 2253 of which were annotated using Gene Ontology (GO) terms. A total of 311 microsatellites (with 246 having sufficient flanking sequences for primer design) and 6294 putative SNPs were identified. These genome resources provide the material basis for future microarray development, marker validation and genetic linkage and QTL analysis.

Chasing relationships between nutrition and reproduction: A comparative transcriptome analysis of hepatopancreas and testis from Eriocheir sinensis

There is a delicate relationship between nutrition and reproduction of mitten crab (Eriocheir sinensis). The crabs store significant amounts of energy in hepatopancreas, which is prepared for significant energy output and expenditure during reproduction, but the internal molecular mechanism has never been known. Here we present the first relationship between hepatopancreas and testis of E. sinensis. We acquired 6287 high quality expressed sequence tags (EST), representing 3829 unigenes totally, from healthy male mitten crabs of first grade. We investigated the Gene Ontology and the main metabolism processes of hepatopancreas and testis from E. sinensis. Genes most likely expressed more frequently and localized in hepatopancreas, and abundant genes from testis for multiple functions. Many genes important for the nutrition regulation are in the EST resource, including arginine kinase, leptin receptor-like protein, seminal plasma glycoprotein 120, and many kinds of zinc finger proteins. The EST data also revealed genes such as heat shock protein 70, testis enhanced gene transcript (TEGT), Cyclin K, etc. predicted to play important roles in regulation of reproduction mechanisms. Among these genes, alignment of leptin receptor-like protein and vasa-like protein from E. sinensis and other species showed even more genomic information on E. sinensis. We identified seventeen genes relevant to control of nutrition mechanisms and eleven genes involved in regulation of reproduction. And this study provides insights into the genetic and molecular mechanisms of nutrition and reproduction in the crab. Such information would facilitate the optimization of breeding in the aquaculture of mitten crabs.