Evaluation of differentially expressed immune-related genes in intestine of Pelodiscus sinensis after intragastric challenge with lipopolysaccharide based on transcriptome analysis

The impact of aquaculture system on the microbiome and gut metabolome of juvenile Chinese softshell turtle (Pelodiscus sinensis)

The commercial aquatic animal microbiome may markedly affect the successful host's farming in various aquaculture systems. However, very little was known about it. Here, two different aquaculture systems, the rice-fish culture (RFC) and intensive pond culture (IPC) systems, were compared to deconstruct the skin, oral, and gut microbiome, as well as the gut metabolome of juvenile Chinese softshell turtle (Pelodiscus sinensis). Higher alpha-diversity and functional redundancy of P. sinensis microbial community were found in the RFC than those of the IPC. The aquaculture systems have the strongest influence on the gut microbiome, followed by the skin microbiome, and finally the oral microbiome. Source-tracking analysis showed that the RFC's microbial community originated from more unknown sources than that of the IPC across all body regions. Strikingly, the RFC's oral and skin microbiome exhibited a significantly higher proportion of generalists and broader habitat niche breadth than those of the IPC, but not the gut. Null model analysis revealed that the RFC's oral and skin microbial community assembly was governed by a significantly greater proportion of deterministic processes than that of the IPC, but not the gut. We further identified the key gene and microbial contribution to five significantly changed gut metabolites, 2-oxoglutarate, N-acetyl-d-mannosamine, cis-4-hydroxy-d-proline, nicotinamide, and l-alanine, which were significantly correlated with important categories of microbe-mediated processes, including the amino acid metabolism, GABAergic synapse, ABC transporters, biosynthesis of unsaturated fatty acids, as well as citrate cycle. Moreover, different aquaculture systems have a significant impact on the hepatic lipid metabolism and body shape of P. sinensis. Our results provide new insight into the influence of aquaculture systems on the microbial community structure feature and assembly mechanism in an aquatic animal, also highlighting the key microbiome and gene contributions to the metabolite variation in the gut microbiome-metabolome association.

Transcriptome Analysis of Immune Responses and Metabolic Regulations of Chinese Soft-Shelled Turtle (Pelodiscus sinensis) against Edwardsiella tarda Infection

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an important aquatic species in southern China that is threatened by many serious diseases. Edwardsiella tarda is one of the highly pathogenic bacteria that cause the white abdominal shell disease. Yet, little is known about the immune and metabolic responses of the Chinese soft-shelled turtle against E. tarda infection. In the paper, gene expression profiles in the turtle liver were obtained to study the immune responses and metabolic regulations induced by E. tarda infection using RNA sequencing. A total of 3908 differentially expressed unigenes between the experimental group and the control group were obtained by transcriptome analysis, among them, were the significantly upregulated unigenes and downregulated unigenes 2065 and 1922, respectively. Further annotation and analysis revealed that the DEGs were mainly enriched in complement and coagulation cascades, phagosome, and steroid hormone biosynthesis pathways, indicating that they were mainly associated with defense mechanisms in the turtle liver against E. tarda four days post infection. For the first time, we reported on the gene profile of anti-E. tarda response in the soft-shelled turtle, and our research might provide valuable data to support further study on anti-E. tarda defense mechanisms in turtles

De novo assembly and functional annotation of blood transcriptome of loggerhead turtle, and in silico characterization of peroxiredoxins and thioredoxins

The aim of this study was to generate and analyze the atlas of the loggerhead turtle blood transcriptome by RNA-seq, as well as identify and characterize thioredoxin (Tnxs) and peroxiredoxin (Prdxs) antioxidant enzymes of the greatest interest in the control of peroxide levels and other biological functions. The transcriptome of loggerhead turtle was sequenced using the Illumina Hiseq 2000 platform and de novo assembly was performed using the Trinity pipeline. The assembly comprised 515,597 contigs with an N50 of 2,631 bp. Contigs were analyzed with CD-Hit obtaining 374,545 unigenes, of which 165,676 had ORFs encoding putative proteins longer than 100 amino acids. A total of 52,147 (31.5%) of these transcripts had significant homology matches in at least one of the five databases used. From the enrichment of GO terms, 180 proteins with antioxidant activity were identified, among these 28 Prdxs and 50 putative Tnxs. The putative proteins of loggerhead turtles encoded by the genes Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip were predicted and characterized in silico. When comparing Prdxs and Txns of loggerhead turtle with homologous human proteins, they showed 18 (9%), 52 (18%) 94 (43%), 36 (16%), 35 (33%) and 74 (19%) amino acid mutations respectively. However, they showed high conservation in active sites and structural motifs (98%), with few specific modifications. Of these, Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip presented 0, 25, 18, three, six and two deleterious changes. This study provides a high quality blood transcriptome and functional annotation of loggerhead sea turtles.

Immune and sex-biased gene expression in the threatened Mojave desert tortoise, Gopherus agassizii

The immune system of ectotherms, particularly non-avian reptiles, remains poorly characterized regarding the genes involved in immune function, and their function in wild populations. We used RNA-Seq to explore the systemic response of Mojave desert tortoise (Gopherus agassizii) gene expression to three levels of Mycoplasma infection to better understand the host response to this bacterial pathogen. We found over an order of magnitude more genes differentially expressed between male and female tortoises (1,037 genes) than differentially expressed among immune groups (40 genes). There were 8 genes differentially expressed among both variables that can be considered sex-biased immune genes in this tortoise. Among experimental immune groups we find enriched GO biological processes for cysteine catabolism, regulation of type 1 interferon production, and regulation of cytokine production involved in immune response. Sex-biased transcription involves iron ion transport, iron ion homeostasis, and regulation of interferon-beta production to be enriched. More detailed work is needed to assess the seasonal response of the candidate genes found here. How seasonal fluctuation of testosterone and corticosterone modulate the immunosuppression of males and their susceptibility to Mycoplasma infection also warrants further investigation, as well as the importance of iron in the immune function and sex-biased differences of this species. Finally, future transcriptional studies should avoid drawing blood from tortoises via subcarapacial venipuncture as the variable aspiration of lymphatic fluid will confound the differential expression of genes.

Oxidative Stress Biomarkers in Erythrocytes of Captive Pre-Juvenile Loggerhead Turtles Following Acute Exposure to Methylmercury

This study describes the use of erythrocytes (RBCs) of loggerhead turtles as in vitro models for evaluating their toxicity to methylmercury. Blood samples of loggerhead turtles that were born in the Colombian Caribbean were used. The LC50 of RBCs to methylmercury was determined at 96 h using methylmercury concentrations of 0.5–100 mg L−1. Next, the viability of the RBCs and the activity of the enzymes superoxide dismutase (SOD), glutathione S-transferase (GST), and lipid peroxidation by malondialdehyde (MDA) at 6 and 12 h of exposure to acute concentrations of 0, 1, and 5 mg L−1 were evaluated. The LC50 for loggerhead turtle RBCs was 8.32 mg L−1. The cell viability bioassay of RBCs exposed for 12 h only showed 100% cell viability. Increasing in vitro MeHg concentrations caused a corresponding increase in MDA concentration as well as decreases in the activities of SOD and GST. The RBCs represent an excellent model for ecotoxicological studies and SOD, GST, and MDA are biomarkers of environmental pollution and oxidative stress in loggerhead turtles. This was the first study conducted on loggerhead turtle where the response of RBCs to MeHg-induced oxidative stress is evaluated.

Antioxidant system of soiny mullet (Liza haematocheila) is responsive to dietary poly-β-hydroxybutyrate (PHB) supplementation based on immune-related enzyme activity and de novo transcriptome analysis

As a dietary supplement, poly-β-hydroxybutyrate (PHB) has been reported to positively influence growth, boost the immune system and enhance disease resistance in fish and shellfish. However, the protective mechanism is little known. Thus, the present study was conducted to evaluate the effect of PHB supplementation on immune-related enzyme activity and transcriptome-based gene expression in soiny mullet (Liza haematocheila). Results showed that dietary PHB supplementation could increase antioxidant enzyme activity, including total antioxidant capacity, catalase and superoxide dismutase. A total of 7,082,094,175 and 7,650,341,357 raw reads with mean length of 757 bp were obtained from control and PHB (dietary PHB supplementation at 2%) groups, respectively. There were 46,106 differentially expressed genes (DEGs) between control and PHB groups, including 21,828 upregulated and 24,278 downregulated DEGs. All the DEGs were classified into three gene ontology categories, and 312 DEGs related with immune system process and 760 with the response to a stimulus. Additionally, all DEGs were allocated to 261 Kyoto Encyclopedia of Gene and Genome pathways, and major immune-related pathways were detected, including MAPK/PI3K-Akt/TNF/NF-κB/TCR/TLR signaling pathways. Moreover, the regulation of several observed immune-related genes was confirmed by qRT-PCR. Altogether, this study suggests that antioxidant system is more effective for dietary PHB supplementation and lays the foundation for further study on the precise immunostimulatory mechanism of PHB. Hopefully, it provides insights into exploring biomarker for assessment of immunostimulants in fish culture.

GWAS Analysis Indicated Importance of NF-κB Signaling Pathway in Host Resistance Against Motile Aeromonas Septicemia Disease in Catfish

Motile Aeromonas septicemia (MAS) disease caused by a bacterial pathogen, Aeromonas hydrophila, is an emerging but severe disease of catfish. Genetic enhancement of disease resistance is considered to be effective to control the disease. To provide an insight into the genomic basis of MAS disease resistance, in this study, we conducted a genome-wide association study (GWAS) to identify quantitative trait loci (QTL). A total of 1820 interspecific backcross catfish of 7 families were challenged with A. hydrophila, and 382 phenotypic extremes were selected for genotyping with the catfish 690 K SNP arrays. Three QTL on linkage group (LG) 2, 26 and 29 were identified to be significantly associated with MAS resistance. Within these regions, a total of 24 genes had known functions in immunity, 10 of which were involved in NF-κB signaling pathway, suggesting the importance of NF-κB signaling pathway in MAS resistance. In addition, three suggestively significant QTL were identified on LG 11, 17, and 20. The limited numbers of QTL involved in MAS resistance suggests that marker-assisted selection may be a viable approach for catfish breeding.

A transcriptome study on Macrobrachium nipponense hepatopancreas experimentally challenged with white spot syndrome virus (WSSV)

White spot syndrome virus (WSSV) is one of the most devastating pathogens of cultured shrimp, responsible for massive loss of its commercial products worldwide. The oriental river prawn Macrobrachium nipponense is an economically important species that is widely farmed in China and adult prawns can be infected by WSSV. However, the molecular mechanisms of the host pathogen interaction remain unknown. There is an urgent need to learn the host pathogen interaction between M. nipponense and WSSV which will be able to offer a solution in controlling the spread of WSSV. Next Generation Sequencing (NGS) was used in this study to determin the transcriptome differences by the comparison of control and WSSV-challenged moribund samples, control and WSSV-challenged survived samples of hepatopancreas in M. nipponense. A total of 64,049 predicted unigenes were obtained and classified into 63 functional groups. Approximately, 4,311 differential expression genes were identified with 3,308 genes were up-regulated when comparing the survived samples with the control. In the comparison of moribund samples with control, 1,960 differential expression genes were identified with 764 genes were up-regulated. In the contrast of two comparison libraries, 300 mutual DEGs with 95 up-regulated genes and 205 down-regulated genes. All the DEGs were performed GO and KEGG analysis, overall a total of 85 immune-related genes were obtained and these gene were groups into 13 functions and 4 KEGG pathways, such as protease inhibitors, heat shock proteins, oxidative stress, pathogen recognition immune receptors, PI3K/AKT/mTOR pathway, MAPK signaling pathway and Ubiquitin Proteasome Pathway. Ten genes that valuable in immune responses against WSSV were selected from those DEGs to furture discuss the response of host to WSSV. Results from this study contribute to a better understanding of the immune response of M. nipponense to WSSV, provide information for identifying novel genes in the absence of genome of M. nipponense. Furthermore, large number of transcripts obtained from this study could provide a strong basis for future genomic research on M. nipponense.

Identification, characterization and expression analysis of MAVS in Pelodiscus sinensis after challenge with Poly I:C

Pelodiscus sinensis, which is one of the important reptile species in the aquaculture industry in China, frequently suffers from serious infectious diseases caused by viruses. However, there is a lack of biological knowledge about its antiviral innate immunity. In this study, we identified and characterized the open reading frame (ORF) of PsMAVS cDNA in P. sinensis. It consisted of 2691 nucleotides encoding a protein of 896 amino acid residues, which were composed of an N-terminal CARD, a central proline-rich domain and a C-terminal TM domain. Based on the amino acid sequence, phylogenetic analyses revealed a closer relationship of PsMAVS with those of Chelonia. qRT-PCR analysis indicated that PsMAVS was ubiquitously expressed in all of the examined healthy tissues with different expression levels; it was expressed at high levels in spleen, muscle and heart and at moderate levels in kidney, liver, intestine, intestinum crissum and oesophagus. PsMAVS was detected in embryos at 10 days post hatching, and it gradually upregulated with the embryonic development stage. Its expression levels in the examined tissues were all upregulated significantly after challenge with Poly I:C. The PsMAVS protein was detected in the intestinal tissues from both the challenge and the control groups, and it was distributed widely in the cytoplasm of the intestinal cells, suggesting PsMAVS plays multiple roles in the complicated mechanisms of immune defence against virus invasion in P. sinensis.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.