Transcriptome Analysis of Paralichthys olivaceus Erythrocytes Reveals Profound Immune Responses Induced by Edwardsiella tarda Infection

Differential tissue immune stimulation through immersion in bacterial and viral agonists in the Antarctic Notothenia rossii

The genome evolution of Antarctic notothenioids has been modulated by their extreme environment over millennia and more recently by human-caused constraints such as overfishing and climate change. Here we investigated the characteristics of the immune system in Notothenia rossii and how it responds to 8 h immersion in viral (Poly I:C, polyinosinic: polycytidylic acid) and bacterial (LPS, lipopolysaccharide) proxies. Blood plasma antiprotease activity and haematocrit were reduced in Poly I:C-treated fish only, while plasma protein, lysozyme activity and cortisol were unchanged with both treatments. The skin and duodenum transcriptomes responded strongly to the treatments, unlike the liver and spleen which had a mild response. Furthermore, the skin transcriptome responded most to the bacterial proxy (cell adhesion, metabolism and immune response processes) and the duodenum (metabolism, response to stress, regulation of intracellular signal transduction, and immune system responses) to the viral proxy. The differential tissue response to the two proxy challenges is indicative of immune specialisation of the duodenum and the skin towards pathogens. NOD-like and C-type lectin receptors may be central in recognising LPS and Poly I:C. Other antimicrobial compounds such as iron and selenium-related genes are essential defence mechanisms to protect the host from sepsis. In conclusion, our study revealed a specific response of two immune barrier tissue, the skin and duodenum, in Notothenia rossii when exposed to pathogen proxies by immersion, and this may represent an adaptation to pathogen infective strategies.

Establishment and identification of a skin cell line from Chinese tongue sole (Cynoglossus semilaevis) and analysis of the changes in its transcriptome upon LPS stimulation

The Chinese tongue sole (Cynoglossus semilaevis) holds significant economic importance within the fishing industry along the eastern coasts of China. In recent years, the frequent outbreaks of bacterial diseases have become a common concern as the aquaculture scale expands. The majority of the diseased fish exhibit symptoms such as skin congestion, damage and skin ulceration. As the skin serves as the first line of defense against bacterial infections, establishing a skin cell line for immunological research on Chinese tongue sole's response to bacterial infection is of utmost importance. In this study, a cell line named CSS (derived from the skin of the Chinese tongue sole) was successfully established. The cells have demonstrated stability during passages and exhibit a multipolar fibroblast-like morphology. They were cultured in L-15 medium with 20% serum and have been successfully passed through 60 passages over a period of 20 months. The identification of the mitochondrial CO1 gene confirmed that the cell originated from Chinese tongue sole. The karyotype detection revealed that the cell had a chromosome number of 2n = 42. After being stored in liquid nitrogen for 15 months, the cells can maintain more than 75% viability upon recovery. After transfecting with cy3-labeled scramble siRNA and pEGFP-N3 plasmid, clear fluorescence was observed in the transfected cells. We observed that lipopolysaccharide (LPS) from Escherichia coli significantly upregulate the gene expression of various immune-related pathways at 2 h in the CSS cell line. Additionally, the differentially expressed genes showed a higher enrichment in immune-related pathways at 2 and 6 h after stimulation compared to the 24 h point. Moreover, we identified 347 genes that exhibited a gradual increase in expression during the 0-24 h stimulation period. These genes were primarily enriched in pathways related to Autophagy, GABAergic synapse, Apelin signaling and Ferroptosis. In general, the CSS cell line established in this study exhibits stable growth and can serve as a valuable tool for in vitro studies of immunology and other basic biologies of Chinese tongue sole.

Immune Activation and Inflammatory Response Mediated by the NOD/Toll-like Receptor Signaling Pathway-The Potential Mechanism of Bullfrog (Lithobates catesbeiana) Meningitis Caused by Elizabethkingia miricola

Elizabethkingia miricola is an emerging opportunistic pathogen that is highly pathogenic in both immunocompromised humans and animals. Once the disease occurs, treatment can be very difficult. Therefore, a deep understanding of the pathological mechanism of Elizabethkingia miricola is the key to the prevention and control of the disease. In this study, we isolated the pathogenic bacteria from bullfrogs with dark skin color, weak limbs, wryneck, and cataracts. Via subsequent morphological observations and a 16S rRNA gene sequence analysis, the pathogen was identified as Elizabethkingia miricola. The histopathological and transmission electron microscopy analysis revealed that the brain was the main target organ. Therefore, brain samples from diseased and healthy bullfrogs were used for the RNA-Seq analysis. The comparative transcriptome analysis revealed that the diseased bullfrog brain was characterized by the immune activation and inflammatory response, which were mediated by the "NOD-like receptor signaling pathway" and the "Toll-like receptor signaling pathway". We also performed qRT-PCR to examine the expression profile of inflammation-related genes, which further verified the reliability of our transcriptome data. Based on the above results, it was concluded that the NOD/Toll-like receptor-related networks that dominate the immune activation and inflammatory response were activated in the brain of Elizabethkingia miricola-infected bullfrogs. This study contributes to the search for therapeutic targets for bullfrog meningitis and provides basic information for establishing effective measures to prevent and control bullfrog meningitis.

Red blood cells of Labeo rohita express Toll-like receptors, NOD- like receptors, interleukins, and interferon-I in response to Gram-negative bacterial infections and lipopolysaccharide stimulations

Red blood cells (RBCs) are the most abundant cell types in the circulatory system of vertebrates. In fish, RBCs retain their nuclei throughout their lifetime and remain transcriptionally and translationally active. While their primary function is typically associated with gas exchange, recent reports indicate that nucleated red blood cells can play a significant role in regulating the body's innate immune response. The current article describes the innate immune role of red blood cells in rohu (Labeo rohita), a freshwater fish species that holds significant commercial importance in India and South-East Asian nations. From the whole blood and mucosal surface RBCs have been isolated through density gradient centrifugation with HiSep™LSM 1077 (density 1.007 ± 0.0010) and their purity has been confirmed by the Giemsa staining followed by microscopical observations. Toll-like receptors (TLR2, 3, 4, 5) and nucleotide oligomerization domain (NOD)-like receptors (NOD1 and NOD2) in RBCs of rohu fingerlings were observed to be significantly activated (P < 0.05) on infection with Aeromonas hydrophila and Edwardsiella tarda. This activation resulted in increased expression of interleukins (IL-8, IL-1β) and interferon (IFN)-I genes. The activation of TLR4, NOD1 and NOD2, as well as the expression of interleukins and IFN-I genes have been observed in both in vivo and in vitro stimulation of rohu RBCs with lipopolysaccharides. These findings highlight the importance of fish RBCs in enhancing innate immunity against various pathogenic invasions in rohu.

Red blood cells in proliferative kidney disease-rainbow trout (Oncorhynchus mykiss) infected by Tetracapsuloides bryosalmonae harbor IgM+ red blood cells

The myxozoan parasite Tetracapsuloides bryosalmonae is the causative agent of proliferative kidney disease (PKD)-a disease of salmonid fishes, notably of the commercially farmed rainbow trout Oncorhynchus mykiss. Both wild and farmed salmonids are threatened by this virulent/deadly disease, a chronic immunopathology characterized by massive lymphocyte proliferation and hyperplasia, which manifests as swollen kidneys in susceptible hosts. Studying the immune response towards the parasite helps us understand the causes and consequences of PKD. While examining the B cell population during a seasonal outbreak of PKD, we unexpectedly detected the B cell marker immunoglobulin M (IgM) on red blood cells (RBCs) of infected farmed rainbow trout. Here, we studied the nature of this IgM and this IgM⁺ cell population. We verified the presence of surface IgM via parallel approaches: flow cytometry, microscopy, and mass spectrometry. The levels of surface IgM (allowing complete resolution of IgM^- RBCs from IgM⁺ RBCs) and frequency of IgM⁺ RBCs (with up to 99% of RBCs being positive) have not been described before in healthy fishes nor those suffering from disease. To assess the influence of the disease on these cells, we profiled the transcriptomes of teleost RBCs in health and disease. Compared to RBCs originating from healthy fish, PKD fundamentally altered RBCs in their metabolism, adhesion, and innate immune response to inflammation. In summary, RBCs play a larger role in host immunity than previously appreciated. Specifically, our findings indicate that the nucleated RBCs of rainbow trout interact with host IgM and contribute to the immune response in PKD.

Splenic protection network revealed by transcriptome analysis in inactivated vaccine-immunized flounder (Paralichthys olivaceus) against Edwardsiella tarda infection

The protective immune response produced by fish after vaccination is crucial for vaccine effectiveness. Our previous studies have shown inactivated vaccine against Edwardsiella tarda can induce immune response in flounder (Paralichthys olivaceus). To elucidate the protective immune response at the genetic level, in this study, flounder was immunized with inactivated E. tarda for 5 weeks, and then they were challenged with E. tarda. The spleen was dissected at 7th day post immunization, 1st and 7th day post challenge, respectively. Transcriptome analysis showed that average of 46 million clean reads were obtained per library, while percentage of clean reads being mapped to reference genome was more than 89% in all cases, which suggested good quality of samples. As for differentially expressed genes (DEGs) identification in inactivated E. tarda groups, at 7th day post immunization, 1422 DEGs were identified and significantly enriched in innate immune-related pathways, such as Phagosome, Cell adhesion molecules and NF-kappa B signaling pathway; At 1st post challenge day, 1210 DEGs were identified and enriched to Antigen processing and presentation and Cell adhesion molecules, indicating that the pathogen was rapidly recognized and delivered; At 7th post challenge day, 1929 DEGs were identified, belonged to Toll-like receptor signaling pathway, Antigen processing and presentation, Th1 and Th2 cell differentiation and Th17 cell differentiation. Compared to 7th post immunization day, 73 immune-associated DEGs were identified at 1st post challenge day. Protein-protein interaction networks analysis revealed 11 hub genes (TLR7, TLR3, CXCR4, IFIH1, TLR8 etc), associated with recognition of pathogens and activation of innate immunity; while for 7th post challenge day, 141 immune-associated DEGs were identified. 30 hub genes (IL6, STAT1, HSP90A.1, TLR7, IL12β etc) were associated with stimulation of lymphocyte differentiation and activation of cellular immunity. Ten immune-related genes were randomly selected for RT-qPCR validation at each time point. In conclusion, data revealed protection of flounder against E. tarda infection by inactivated vaccine is mediated via immediate recognition of pathogen and subsequently activation of cellular immunity. Results give new aspect for vaccine protection cascades, is good references for vaccine evaluation.

Temperature induces metabolic reprogramming in fish during bacterial infection

Water temperature elevation as a consequence of global warming results in increased incidence of bacterial disease, such as edwardsiellosis, in fish farming. Edwardsiellosis is caused by the bacterial pathogen Edwardsiella tarda and affects many farmed fish including flounder (Paralichthys olivaceus). Currently, the effect of temperature on the metabolic response of flounder to E. tarda infection is unclear. In this study, we found that compared to low temperature (15°C), high temperature (23°C) enhanced E. tarda dissemination in flounder tissues. To examine the impact of temperature on the metabolism of flounder induced by E. tarda, comparative metabolomics were performed, which identified a large number of metabolites responsive to E. tarda invasion and temperature alteration. During E. tarda infection, the metabolic profile induced by elevated temperature was mainly featured by extensively decreased amino acids and TCA intermediates such as succinate, a proven immune regulator. Further, 38 potential metabolite markers of temperature effect (MMTE) in association with bacterial infection were identified. When used as exogenous supplements, two of the MMTE, i.e., L-methionine and UDP-glucose, effectively upregulated the expression of pro-inflammatory cytokines and suppressed E. tarda infection in flounder leukocytes. Taken together, the results of this study indicate an important influence of temperature on the metabolism of flounder during bacterial infection, which eventually affects the survivability of the fish.

The Evolution of a Specialized, Highly Virulent Fish Pathogen through Gene Loss and Acquisition of Host-Specific Survival Mechanisms

Photobacterium damselae comprises two subspecies, P. damselae subsp. damselae and P. damselae subsp. piscicida, that contrast remarkably despite their taxonomic relationship. The former is opportunistic and free-living but can cause disease in compromised individuals from a broad diversity of taxa, while the latter is a highly specialized, primary fish pathogen. Here, we employ new closed curated genome assemblies from Australia to estimate the global phylogenetic structure of the species P. damselae. We identify genes responsible for the shift from an opportunist to a host-adapted fish pathogen, potentially via an arthropod vector as fish-to-fish transmission was not achieved in repeated cohabitation challenges despite high virulence for Seriola lalandi. Acquisition of ShdA adhesin and of thiol peroxidase may have allowed the environmental, generalist ancestor to colonize zooplankton and to occasionally enter in fish host sentinel cells. As dependence on the host has increased, P. damselae has lost nonessential genes, such as those related to nitrite and sulfite reduction, urea degradation, a type 6 secretion system (T6SS) and several toxin-antitoxin (TA) systems. Similar to the evolution of Yersinia pestis, the loss of urease may be the crucial event that allowed the pathogen to stably colonize zooplankton vectors. Acquisition of host-specific genes, such as those required to form a sialic acid capsule, was likely necessary for the emergent P. damselae subsp. piscicida to become a highly specialized, facultative intracellular fish pathogen. Processes that have shaped P. damselae subsp. piscicida from subsp. damselae are similar to those underlying evolution of Yersinia pestis from Y. pseudotuberculosis.

IMPORTANCEPhotobacterium damselae subsp. damselae is a ubiquitous marine bacterium and opportunistic pathogen of compromised hosts of diverse taxa. In contrast, its sister subspecies P. damselae subsp. piscicida (Pdp) is highly virulent in fish. Pdp has evolved from a single subclade of Pdd through gene loss and acquisition. We show that fish-to-fish transmission does not occur in repeated infection models in the primary host, Seriola lalandi, and present genomic evidence for vector-borne transmission, potentially via zooplankton. The broad genomic changes from generalist Pdd to specialist Pdp parallel those of the environmental opportunist Yersinia pseudotuberculosis to vector-borne plague bacterium Y. pestis and demonstrate that evolutionary processes in bacterial pathogens are universal between the terrestrial and marine biosphere.

Molecular Identification of Nocardia seriolae and Comparative Analysis of Spleen Transcriptomes of Hybrid Snakehead (Channa maculata Female × Channa argus Male) With Nocardiosis Disease

Hybrid snakehead (Channa maculata female × Channa argus male) is a new freshwater aquaculture fish species in southern China. During intensive aquaculture, hybrid snakeheads are often infected by Nocardia seriolae. In this study, hybrid snakehead infected suspiciously by N. seriolae in an artificial breeding pond were examined. Diseased hybrid snakeheads swam slowly without food intake, and the clinical symptoms included skin wound, anal swelling and ascites, and white granulomatous in liver, spleen, and kidney of fish. Through bacterial isolation, 16S rDNA sequencing, fluorescence in situ hybridization (FISH) and artificial infection experiment, the pathogen was identified as N. seriolae. Furthermore, the spleen samples from diseased and healthy male hybrid snakeheads in the same pond were used for RNA-Seq analysis. A total of 3,512 unique transcripts (unigenes) were identified as differentially expressed genes (DEGs), and 1,886 of them were up-regulated in diseased fish. The expression patterns of 20 DEGs were verified by quantitative polymerase chain reaction (qPCR). Several immune-related pathways and many immune-related genes were identified. qPCR results showed that the expression patterns of immune-related genes in the liver and kidney of diseased fish were comparable to that in the spleen. This study provides deep-sequencing data of hybrid snakehead spleen and will help understand the immune response of hybrid snakehead to N. seriolae. It is also helpful for the biomarker screening of fish-borne Nocardia spp. and the breeding of nocardiosis-resistant fish species.

Genome-wide characterization of gap junction (connexins and pannexins) genes in turbot (Scophthalmus maximus L.): evolution and immune response following Vibrio anguillarum infection

Gap junction (GJ), a special intercellular junction between different cell types, directly connects the cytoplasm of adjacent cells, allows various molecules, ions and electrical impulses to pass through the intercellular regulatory gate, and plays vital roles in response to bacterial infection. Up to date, the information about the GJ in turbot (Scophthalmus maximus L.) is still limited. In current study, 43 gap junction genes were identified in turbot, phylogeny analysis suggested that gap junctions from turbot and other species were clustered into six groups, GJA, GJB, GJC, GJD, GJE and PANX, and turbot GJs together with respective GJs from Japanese flounder, half-smooth tongue sole and large yellow croaker, sharing same ancestors. In addition, these 43 GJ genes distributed in different chromosomes unevenly. According to gene structure and domain analysis, these genes (in GJA-GJE group) were highly conserved in that most of them contain the transmembrane area, connexin domain (CNX) and cysteine-rich domain (connexin CCC), while PANXs contain Pfam Innexin. Although only one tandem duplication was identified in turbot gap junction gene, 235 pairs of segmental duplications were identified in the turbot genome. To further investigate their evolutionary relationships, Ka/Ks was calculated, and results showed that most ratios were lower than 1, indicating they had undergone negative selection. Finally, expression analysis showed that gap junction genes were widely distributed in turbot tissues and significantly regulated after Vibrio anguillarum infection. Taken together, our research could provide valuable information for further exploration of the function of gap junction genes in teleost.

Transcriptomic analysis using dual RNA sequencing revealed a Pathogen-Host interaction after Edwardsiella anguillarum infection in European eel (Anguilla anguilla)

Many studies have explored differentially expressed genes (DEGs) between some pathogens and hosts, but no study has focused on the interaction of DEGs between Edwardsiella anguillarum (Ea) and Anguilla anguilla (Aa). In this study, we examined the interactions of DEGs during Ea infection and Aa anti-infection processes by dual RNA sequencing. Total RNA from in vitro and in vivo (Aa liver) Ea culture was extracted. Using high-throughput transcriptomics, significant DEGs that were expressed between Ea cultured in vitro versus in vivo and those in the liver of the infected group versus control group were identified. Protein-protein interactions between the pathogen and host were explored using Cytoscape according to the HPIDB 3.0 interaction transcription database. The results showed that the liver in the infection group presented with severe bleeding and a large number of thrombi in the hepatic vessels. We found 490 upregulated and 398 downregulated DEGs of Ea in vivo versus Ea cultured in vitro, and 2177 upregulated and 970 downregulated genes in the liver of the infected eels. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the pathogen DEGs revealed that the upregulated genes were mainly enriched in migration, colonization, biofilm formation, and significantly enriched in ABC transport and quorum sensing; the downregulated genes were mainly involved in metabolism, information transduction, organelle formation, enzyme catalysis, molecular transport, and binding. GO of the host DEGs showed that metabolic process, catalytic activity, single organism metabolic process, small molecule binding, nucleotide binding, nucleotide phosphate binding, and anion binding were markedly enriched. Finally, we found that 79 Ea and 148 Aa proteins encoded by these DEGs were involved in an interaction network, and some pathogen (DegP, gcvP, infC, carB, rpoC, trpD, sthA, and FhuB) and host proteins (MANBA, STAT1, ETS2, ZEP1, TKT1, NMI and RBPMS) appear to play crucial roles in infection. Thus, determining the interaction networks revealed crucial molecular mechanisms underlying the process of pathogenic infection and host anti-infection.

Transcriptome of Edwardsiella anguillarum in vivo and in vitro revealed two-component system, ABC transporter and flagellar assembly are three pathways pathogenic to European eel (Anguilla anguilla)

Edwardsiella anguillarum is one of the common bacterial pathogens for the cultivated eels in China. The aim of this study was to reveal the cause of E. anguillarum pathogenic to European eel (Anguilla anguilla) from the perspective of the transcriptome. In this study, we first prepared E. anguillarum cultured in vitro and analysed the whole transcriptome after extracting the total RNA. Then, eels were i.p injected with E. anguillarum, and total RNA were extracted from the liver of European eels 48 h after the infection. After sequencing the transcriptome, we obtained average 1.97 × 10⁸ clean reads cultured in vitro and 1.36 × 10⁵ clean reads located in vivo after annotating all reads into the genome of E. anguillarum. The whole transcriptome showed, compared to the E. anguillarum cultured in vitro, 503 significantly up and 657 significantly down-regulated different expressed genes (DEGs) were observed. KEGG analysis showed that 38 DEGs of Two-Component System, 41 DEGs of ABC transporter, and 10 DEGs flagellar assembly pathways were highly upregulated in E. anguillarum located in vivo. Then, we designed primers to analyse the up-regulated DEGs through qRT-PCR and confirmed some up-regulated DEGs. The results of this study provide important reference for the further study of pathogen-host interaction between E. anguillarum and European eel.

Identification and Characterization of lncRNAs Related to the Muscle Growth and Development of Japanese Flounder (Paralichthys olivaceus)

Long noncoding RNAs (lncRNAs) play an important role in many life activities, but the expression pattern and function of lncRNAs in Japanese flounder skeletal muscle are unclear. In this study, 751 lncRNAs were selected from skeletal muscle in different development stages of the Japanese flounder [stage A: larval 7 days post hatching (dph); stage B: juvenile about 90 dph; stage C (female) and stage D (male): adult about 24 months] using coding potential analysis methods. In total, 232, 211, 194, 28, 29, and 14 differentially expressed lncRNAs and 9549, 8673, 9181, 1821, 1080, and 557 differentially expressed mRNAs were identified in comparisons of A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. We identified the cis- and trans-regulatory target genes of differentially expressed lncRNAs, and lncRNA-gene interaction networks were constructed using the Cytoscape program. In total, there were 200, 200, 200, 93, 47, and 11 cis-regulation relationships, and 29, 19, 24, 38, 8, and 47 trans-regulation relationships in the comparisons between A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. These results indicate that lncRNA may participate in the development of Japanese flounder skeletal muscle through cis- or trans-acting mechanisms, thus providing a scientific basis for further study of the biological function of lncRNA in Japanese flounder skeletal muscle. Based on these relationships, functional annotation of the related lncRNAs was performed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Differentially expressed genes associated with muscle development were enriched in multiple pairs of comparisons (e.g., differentially expressed genes LOC109640370, LOC109634180, LOC109643555, rusc1, and LOC109626999 were enriched in the actin-binding term, and differentially expressed genes LOC109640370, was, LOC109644970, LOC109643555, and itga9 were enriched in the regulation of the actin cytoskeleton pathway in the KEGG pathway analysis in the comparison of stages C and D). We predicted lncRNA-mRNA, miRNA-mRNA, and lncRNA-miRNA regulatory relationships and constructed interactive networks using Cytoscape software. Co-expression networks show that most lncRNAs interact with one or two predicted miRNAs involved in muscle growth and development. These results provide a basis for further study of the function of lncRNAs on skeletal muscle in different developmental stages of Japanese flounder.