Meta-analysis to Unravel Core Transcriptomic Responses in Penaeus vannamei Exposed to Biotic and Abiotic Stresses

Shrimp farming, a dominant economic activity in coastal areas, is affected by different abiotic and biotic stress factors. These stressors, under poor management conditions, could affect growth and health of farmed animals. Understanding the common gene expressions in response to stress, regardless of the specific stress factor, holds significant importance in the field of functional genomics. Scope of this study is to identify the core transcriptomic responses in the shrimp species Penaeus vannamei exposed to various abiotic and biotic stress conditions and to decipher their functional importance. To achieve our objective, we gathered and analyzed multiple RNA-seq datasets related to twelve abiotic and nine biotic stress conditions. Through the in silico meta-analysis, we predicted 961 differentially expressed genes (meta-DEGs) for abiotic stress conditions and 517 meta-DEGs for biotic stress conditions, respectively. These meta-DEGs represent genes that are commonly expressed across different stress factors and are indicative of the organism's general response to stress. The annotation of nineteen core up-regulated meta-DEGs revealed their diverse functions in detoxification, cell adhesion, metal ion binding, and oxidative phosphorylation. These genes play a crucial role in stress response and immune defense. For abiotic stress, significant pathways associated with the stress response include tryptophan metabolism, starch and sucrose metabolism, fatty acid degradation, carbohydrate digestion and absorption, phenylalanine metabolism, drug metabolism-other enzymes, arachidonic acid metabolism, and fatty acid elongation. Similarly, for biotic stress, metabolism of xenobiotics by cytochrome P450, pentose and glucuronate interconversions, steroid hormone biosynthesis, and drug metabolism-cytochrome P450 were found to be significant pathway associations. In addition, the study also predicted 17 stress regulatory motifs present in the identified meta-DEGs. These motifs have significance in identifying the stress responses of the organism. The metabolic pathways and regulatory motifs associated with abiotic and biotic stress factors identified through this study could be a valuable resource for developing stress management approaches in shrimp aquaculture.

Transcriptomic variations associated with salinity stress in Penaeus indicus

BACKGROUND

The Indian white shrimp, Penaeus indicus a native species of India, is important brackishwater aquaculture species. Shrimps are euryhaline in nature and they regulate osmotic and ionic concentrations by osmoregulatory process. However, variations in abiotic factors such as salinity result in stress to the shrimps during culture period affecting their growth and immunity.

METHODS AND RESULTS

To understand the adaptive mechanism to stress in low salinity conditions, RNA-seq was used to compare the transcriptomic response of P. indicus upto 3 weeks. De novo assembly using Trinity assembler generated a total of 173,582 transcripts. The assembly had a mean length of 854 bp, N50 value of 1243 bp and GC content of 42.33%. Differential gene expression analysis, resulted in identification of 2130, 3090, and 5351 DEGs in 7 days, 14 days and 21 days respectively of salinity stress period. The pathway prediction of the assembled trinity transcripts using KEGG database showed total number of 329 pathways linking 12,430 transcripts. KEGG pathway enrichment analyses led to the identification of several enriched pathways related to lipid metabolism, amino acid metabolism, glycolysis, signalling pathways etc. Selected genes involved in osmoregulatory process and immune response in shrimps were validated and analysed for the gene expression levels by quantitative real-time PCR (qPCR).

CONCLUSION

This study on the adaptive transcriptomic response of P. indicus to low salinity, will further help in our understanding of the molecular mechanisms underlying osmoregulation mechanism in shrimps.

Exploration of the candidate beneficial bacteria for Penaeus vannamei culture by core microbiome analysis using amplicon sequencing

Globally, Penaeus vannamei is the vital species in aquaculture production. Beneficial bacterial exploration of gut, sediment, and water were investigated in P. vannamei culture using Illumina Miseq sequencing of 16S RNA V3-V4 hypervariable regions. Predominant phyla identified were Proteobacteria, Tenericutes, Bacteroidetes in gut; Proteobacteria, Bacteroidetes, Planctomycetes in sediment and Cyanobacteria, Proteobacteria, and Planctomycetes in water. In total, 46 phyla, 509 families and 902 genera; 70 phyla, 735 families and 1255 genera; 55 phyla, 580 families and 996 genera were observed in gut, sediment and water, respectively. Diversity of microbial communities in respect of observed Operational Taxonomic Units, diversity indices (Shannon and Simpson), richness index (Chao1) were significantly high P (<0.05) in 60 DoC in gut and 30 DoC in sediment. Beta diversity indicated separate clusters for bacterial communities in gut, sediment and water samples and formation of distinct community profiles. Core microbiome in P. vannamei rearing ponds over a time consisted of 9, 21, and 20 OTUs in gut, rearing water and sediment, respectively. This study helps to intervene with suitable beneficial microbes to establish an aquaculture system thereby contributes to enhance the productivity, improve water quality and pond bottom condition, and control the pathogenic agents at each stage of the culture.

MRF: a tool to overcome the barrier of inconsistent genome annotations and perform comparative genomics studies for the largest animal DNA virus

BACKGROUND

The genome of the largest known animal virus, the white spot syndrome virus (WSSV) responsible for huge economic losses and loss of employment in aquaculture, suffers from inconsistent annotation nomenclature. Novel genome sequence, circular genome and variable genome length led to nomenclature inconsistencies. Since vast knowledge has already accumulated in the past two decades with inconsistent nomenclature, the insights gained on a genome could not be easily extendable to other genomes. Therefore, the present study aims to perform comparative genomics studies in WSSV on uniform nomenclature.

METHODS

We have combined the standard mummer tool with custom scripts to develop missing regions finder (MRF) that documents the missing genome regions and coding sequences in virus genomes in comparison to a reference genome and in its annotation nomenclature. The procedure was implemented as web tool and in command-line interface. Using MRF, we have documented the missing coding sequences in WSSV and explored their role in virulence through application of phylogenomics, machine learning models and homologous genes.

RESULTS

We have tabulated and depicted the missing genome regions, missing coding sequences and deletion hotspots in WSSV on a common annotation nomenclature and attempted to link them to virus virulence. It was observed that the ubiquitination, transcription regulation and nucleotide metabolism might be essentially required for WSSV pathogenesis; and the structural proteins, VP19, VP26 and VP28 are essential for virus assembly. Few minor structural proteins in WSSV would act as envelope glycoproteins. We have also demonstrated the advantage of MRF in providing detailed graphic/tabular output in less time and also in handling of low-complexity, repeat-rich and highly similar regions of the genomes using other virus cases.

CONCLUSIONS

Pathogenic virus research benefits from tools that could directly indicate the missing genomic regions and coding sequences between isolates/strains. In virus research, the analyses performed in this study provides an advancement to find the differences between genomes and to quickly identify the important coding sequences/genomes that require early attention from researchers. To conclude, the approach implemented in MRF complements similarity-based tools in comparative genomics involving large, highly-similar, length-varying and/or inconsistently annotated viral genomes.

In silico prediction of potential indigenous microbial biomarkers in Penaeus vannamei identified through meta-analysis and genome-scale metabolic modelling

BACKGROUND

Understanding the microbiome is crucial as it contributes to the metabolic health of the host and, upon dysbiosis, may influence disease development. With the recent surge in high-throughput sequencing technology, the availability of microbial genomic data has increased dramatically. Amplicon sequence-based analyses majorly profile microbial abundance and determine taxonomic markers. Furthermore, the availability of genome sequences for various microbial organisms has prompted the integration of genome-scale metabolic modelling that provides insights into the metabolic interactions influencing host health. However, the analysis from a single study may not be consistent, necessitating a meta-analysis.

RESULTS

We conducted a meta-analysis and integrated with constraint-based metabolic modelling approach, focusing on the microbiome of pacific white shrimp Penaeus vannamei, an extensively cultured marine candidate species. Meta-analysis revealed that Acinetobacter and Alteromonas are significant indicators of "health" and "disease" specific taxonomic biomarkers, respectively. Further, we enumerated metabolic interactions among the taxonomic biomarkers by applying a constraint-based approach to the community metabolic models (4416 pairs). Under different nutrient environments, a constraint-based flux simulation identified five beneficial species: Acinetobacter spWCHA55, Acinetobacter tandoii SE63, Bifidobacterium pseudolongum 49 D6, Brevundimonas pondensis LVF1, and Lutibacter profundi LP1 mediating parasitic interactions majorly under sucrose environment in the pairwise community. The study also reports the healthy biomarkers that can co-exist and have functionally dependent relationships to maintain a healthy state in the host.

CONCLUSIONS

Toward this, we collected and re-analysed the amplicon sequence data of P. vannamei (encompassing 117 healthy and 142 disease datasets). By capturing the taxonomic biomarkers and modelling the metabolic interaction between them, our study provides a valuable resource, a first-of-its-kind analysis in aquaculture scenario toward a sustainable shrimp farming.

In Silico Prediction of Novel Probiotic Species Limiting Pathogenic Vibrio Growth Using Constraint-Based Genome Scale Metabolic Modeling

The prevalence of bacterial diseases and the application of probiotics to prevent them is a common practice in shrimp aquaculture. A wide range of bacterial species/strains is utilized in probiotic formulations, with proven beneficial effects. However, knowledge of their role in inhibiting the growth of a specific pathogen is restricted. In this study, we employed constraint-based genome-scale metabolic modeling approach to screen and identify the beneficial bacteria capable of limiting the growth of V. harveyi, a common pathogen in shrimp culture. Genome-scale models were built for 194 species (including strains from the genera Bacillus, Lactobacillus, and Lactococcus and the pathogenic strain V. harveyi) to explore the metabolic potential of these strains under different nutrient conditions in a consortium. In silico-based phenotypic analysis on 193 paired models predicted six candidate strains with growth enhancement and pathogen suppression. Growth simulations reveal that mannitol and glucoronate environments mediate parasitic interactions in a pairwise community. Furthermore, in a mannitol environment, the shortlisted six strains were purely metabolite consumers without donating metabolites to V. harveyi. The production of acetate by the screened species in a paired community suggests the natural metabolic end product's role in limiting pathogen survival. Our study employing in silico approach successfully predicted three novel candidate strains for probiotic applications, namely, Bacillus sp 1 (identified as B. licheniformis in this study), Bacillus weihaiensis Alg07, and Lactobacillus lindneri TMW 1.1993. The study is the first to apply genomic-scale metabolic models for aquaculture applications to detect bacterial species limiting Vibrio harveyi growth.

Phylogenetic relations and mitogenome-wide similarity metrics reveal monophyly of Penaeus sensu lato

Splitting of the genus Penaeus sensu lato into six new genera based on morphological features alone has been controversial in penaeid shrimp taxonomy. Several studies focused on building phylogenetic relations among the genera of Penaeus sensu lato. However, they lack in utilizing full mitochondrial DNA genome of shrimp representing all the six controversial genera. For the first time, the present study targeted the testing of all the six genera of Penaeus sensu lato for phylogenetic relations utilizing complete mitochondrial genome sequence. In addition, the study reports for the first time about the complete mitochondrial DNA genome sequence of Fenneropenaeus indicus, an important candidate species in aquaculture and fisheries, and utilized it for phylogenomics. The maximum likelihood and Bayesian approaches were deployed to generate and comprehend the phylogenetic relationship among the shrimp in the suborder, Dendrobranchiata. The phylogenetic relations established with limited taxon sampling considered in the study pointed to the monophyly of Penaeus sensu lato and suggested collapsing of the new genera to a single genus. Further, trends in mitogenome-wide estimates of average amino acid identity in the order Decapoda and the genus Penaeus sensu lato supported restoration of the old genus, Penaeus, rather promoting the creation of new genera.