Dual RNA-Seq Uncovers Metabolic Amino Acids Dependency of the Intracellular Bacterium Piscirickettsia salmonis Infecting Atlantic Salmon

Transcriptome profiling of lumpfish (Cyclopterus lumpus) head kidney to Renibacterium salmoninarum at early and chronic infection stages

Renibacterium salmoninarum causes Bacterial Kidney Disease (BKD) in several fish species. Atlantic lumpfish, a cleaner fish, is susceptible to R. salmoninarum. To profile the transcriptome response of lumpfish to R. salmoninarum at early and chronic infection stages, fish were intraperitoneally injected with either a high dose of R. salmoninarum (1 × 109 cells dose-1) or PBS (control). Head kidney tissue samples were collected at 28- and 98-days post-infection (dpi) for RNA sequencing. Transcriptomic profiling identified 1971 and 139 differentially expressed genes (DEGs) in infected compared with control samples at 28 and 98 dpi, respectively. At 28 dpi, R. salmoninarum-induced genes (n = 434) mainly involved in innate and adaptive immune response-related pathways, whereas R. salmoninarum-suppressed genes (n = 1537) were largely connected to amino acid metabolism and cellular processes. Cell-mediated immunity-related genes showed dysregulation at 98 dpi. Several immune-signalling pathways were dysregulated in response to R. salmoninarum, including apoptosis, alternative complement, JAK-STAT signalling, and MHC-I dependent pathways. In summary, R. salmoninarum causes immune suppression at early infection, whereas lumpfish induce a cell-mediated immune response at chronic infection. This study provides a complete depiction of diverse immune mechanisms dysregulated by R. salmoninarum in lumpfish and opens new avenues to develop immune prophylactic tools to prevent BKD.

KLF17 is an important regulatory component of the transcriptomic response of Atlantic salmon macrophages to Piscirickettsia salmonis infection

Piscirickettsia salmonis is the most important health problem facing Chilean Aquaculture. Previous reports suggest that P. salmonis can survive in salmonid macrophages by interfering with the host immune response. However, the relevant aspects of the molecular pathogenesis of P. salmonis have been poorly characterized. In this work, we evaluated the transcriptomic changes in macrophage-like cell line SHK-1 infected with P. salmonis at 24- and 48-hours post-infection (hpi) and generated network models of the macrophage response to the infection using co-expression analysis and regulatory transcription factor-target gene information. Transcriptomic analysis showed that 635 genes were differentially expressed after 24- and/or 48-hpi. The pattern of expression of these genes was analyzed by weighted co-expression network analysis (WGCNA), which classified genes into 4 modules of expression, comprising early responses to the bacterium. Induced genes included genes involved in metabolism and cell differentiation, intracellular transportation, and cytoskeleton reorganization, while repressed genes included genes involved in extracellular matrix organization and RNA metabolism. To understand how these expression changes are orchestrated and to pinpoint relevant transcription factors (TFs) controlling the response, we established a curated database of TF-target gene regulatory interactions in Salmo salar, SalSaDB. Using this resource, together with co-expression module data, we generated infection context-specific networks that were analyzed to determine highly connected TF nodes. We found that the most connected TF of the 24- and 48-hpi response networks is KLF17, an ortholog of the KLF4 TF involved in the polarization of macrophages to an M2-phenotype in mammals. Interestingly, while KLF17 is induced by P. salmonis infection, other TFs, such as NOTCH3 and NFATC1, whose orthologs in mammals are related to M1-like macrophages, are repressed. In sum, our results suggest the induction of early regulatory events associated with an M2-like phenotype of macrophages that drives effectors related to the lysosome, RNA metabolism, cytoskeleton organization, and extracellular matrix remodeling. Moreover, the M1-like response seems delayed in generating an effective response, suggesting a polarization towards M2-like macrophages that allows the survival of P. salmonis. This work also contributes to SalSaDB, a curated database of TF-target gene interactions that is freely available for the Atlantic salmon community.

Alternative splicing in Atlantic salmon head kidney and SHK-1 cell line during the Piscirickettsia salmonis infection: A comparative transcriptome survey

Piscirickettsia salmonis, an intracellular bacterium in salmon aquaculture, is a big challenge because it is responsible for 54.2% of Atlantic salmon mortalities. In recent years, the high relevance of Alternative Splicing (AS) as a molecular mechanism associated with infectious conditions and host-pathogen interaction processes, especially in host immune activation, has been observed. Several studies have highlighted the role of AS in the host's immune response during viral, bacterial, and endoparasite infection. In the present study, we evaluated AS transcriptome profiles during P. salmonis infection in the two most used study models, SHK-1 cell line and salmon head kidney tissue. First, the SHK-1 cell line was exposed to P. salmonis infection at 0-, 7-, and 14-days post-infection (dpi). Following, total RNA was extracted for Illumina sequencing. On the other hand, RNA-Seq datasets of Atlantic salmon head kidney infected with the same P. salmonis strayingwase used. For both study models, the highest number of differentially alternative splicing (DAS) events was observed at 7 dpi, 16,830 DAS events derived from 9213 DAS genes in SHK-1 cells, and 13,820 DAS events from 7684 DAS genes in salmon HK. Alternative first exon (AF) was the most abundant AS type in the three infection times analyzed, representing 31% in SHK-1 cells and 228.6 in salmon HK; meanwhile, mutually exclusive exon (MX) was the least abundant. Notably, functional annotation of DAS genes in SHK-1 cells infected with P. salmonis showed a high presence of genes related to nucleotide metabolism. In contrast, the salmon head kidney exhibited many GO terms associated with immune response. Our findings reported the role of AS during P. salmonis infection in Atlantic salmon. These studies would contribute to a better understanding of the molecular bases that support the pathogen-host interaction, evidencing the contribution of AS regulating the transcriptional host response.

Microbial signature profiles of Penaeus vannamei larvae in low-survival hatchery tanks affected by vibriosis

Vibriosis is caused by some pathogenic Vibrio and produces significant mortality in Pacific white shrimp Penaeus (Litopenaeus) vannamei larvae in commercial hatcheries. Acute hepatopancreatic necrosis disease (AHPND) is an emerging vibriosis affecting shrimp-producing countries worldwide. Zoea 2 syndrome is another type of vibriosis that affects the early stages of P. vannamei larvae. Although the pathogenesis of AHPND and zoea 2 syndrome is well known, there is scarce information about microbial composition and biomarkers of P.vannamei larvae affected by AHPND, and there is no study of the microbiome of larvae affected by zoea 2 syndrome. In this work, we characterized the microbiome of P. vannamei larvae collected from 12 commercial hatchery tanks by high-throughput sequencing. Seven tanks were affected by AHPND, and five tanks were affected by zoea 2 syndrome. Subsequently, all samples were selected for sequencing of the V3-V4 region of the16S rRNA gene. Similarity analysis using the beta diversity index revealed significant differences in the larval bacterial communities between disease conditions, particularly when Vibrio was analyzed. Linear discriminant analysis with effect size determined specific microbial signatures for AHPND and zoea 2 syndrome. Sneathiella, Cyclobacterium, Haliea, Lewinella, among other genera, were abundant in AHPND-affected larvae. Meanwhile, Vibrio, Spongiimonas, Meridianimaribacter, Tenacibaculum, among other genera, were significantly abundant in larvae affected by zoea 2 syndrome. The bacterial network at the phylum level for larvae collected from tanks affected by AHPND showed greater complexity and connectivity than in samples collected from tanks affected by zoea 2 syndrome. The bacterial connections inter Vibrio genera were higher in larvae from tanks affected by zoea 2 syndrome, also presenting other connections between the genera Vibrio and Catenococcus. The identification of specific biomarkers found in this study could be useful for understanding the microbial dynamics during different types of vibriosis.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Rejection of Lepeophtheirus salmonis driven in part by chitin sensing is not impacted by seawater acclimitization in Coho salmon (Oncorhynchus kisutch)

There is tremendous variation in life-history strategies among anadromous salmonids. Species that enter the ocean environment at small sizes (< 20 g) are likely under more physiological pressure from pathogens; however, little data is available on responses at these early stages. With this in mind, we performed salmon louse challenges with Coho salmon either immediately after seawater entry (SW; ca. 10 g) or after 30 days in SW (ca. 20 g). Irrespective of size or time in SW, parasites were rapidly rejected by the host, with > 90% of all parasites lost by 16 days post-infection (dpi). Rejection was concomitant with host epithelial granulomatous infiltrations that initially targeted the embedded frontal filament (4 dpi) and the entire parasite by 10 dpi. Illumina sequencing, followed by functional enrichment analysis, revealed a concerted defense response in the fin within 1 dpi that included multiple innate and adaptive immunity components. Strikingly, early indications of an allergic-type inflammatory response were associated with chitin sensing pathways orchestrated by early overexpression of the IgE-receptor, fcer1g. Additionally, there was profound overexpression of several classes of c-type lectin receptors, including dectin-2, mincle, and dc-sign at 1 dpi onward. These profiles and upregulation of cellular effector markers were corroborated by histopathological evaluation, revealing the simultaneous presence of mast cell/eosinophilic granular cells, sacciform cells, macrophages/histiocytes, and granulocytes in fin. At 10 dpi and concurrent with parasite expulsion, there was evidence of immunoregulation in addition to tissue remodelling pathways. At 16 dpi, the response was effectively abrogated. Simultaneous profiling of the parasite transcriptome revealed early induction of chitin metabolism and immunomodulation, toxin production and ECM degradation; however, after 7 dpi, these were replaced with overexpression of stress and immune defense genes. These data present the first evidence for Coho salmon demonstrating chitin- and sugar moiety-sensing as key drivers of salmon louse rejection.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Why vaccines fail against Piscirickettsiosis in farmed salmon and trout and how to avoid it: A review

Piscirickettsiosis is the most severe, persistent, and damaging disease that has affected the Chilean salmon industry since its origins in the 1980s. As a preventive strategy for this disease, different vaccines have been developed and used over the last 30 years. However, vaccinated salmon and trout frequently die in the sea cages and the use of antibiotics is still high demonstrating the low efficiency of the available vaccines. The reasons why the vaccines fail so often are still debated, but it could involve different extrinsic and intrinsic factors. Among the extrinsic factors, mainly associated with chronic stress, we can distinguish: 1) biotic including coinfection with sea lice, sealions attacks or harmful algal blooms; 2) abiotic including low oxygen or high temperature; and 3) farm-management factors including overcrowding or chemical delousing treatments. Among the intrinsic factors, we can distinguish: 1) fish-related factors including host's genetic variability (species, population and individual), sex or age; 2) pathogen-related factors including their variability and ability to evade host immune responses; and 3) vaccine-related factors including low immunogenicity and poor matches with the circulating pathogen strain. Based on the available evidence, in order to improve the development and the efficacy of vaccines against P. salmonis we recommend: a) Do not perform efficacy evaluations by intraperitoneal injection of pathogens because they generate an artificial protective immune response, instead cohabitation or immersion challenges must be used; b) Evaluate the diversity of pathogen strains in the field and ensure a good antigenic match with the vaccines; c) Investigate whether host genetic diversity could be improved, e.g. through selection, in favor of better and longer responses to vaccination; d) To reduce the stressful effects at the cage level, controlling the co-infection of pathogens and avoiding fish overcrowding. To date, we do not know the immunological mechanisms by which the vaccines against P. salmonis may or may not generate protection. More studies are required to identify what type of response, cellular or molecular, is required to develop effective vaccines.

dRNASb: a systems biology approach to decipher dynamics of host-pathogen interactions using temporal dual RNA-seq data

Infection triggers a dynamic cascade of reciprocal events between host and pathogen wherein the host activates complex mechanisms to recognise and kill pathogens while the pathogen often adjusts its virulence and fitness to avoid eradication by the host. The interaction between the pathogen and the host results in large-scale changes in gene expression in both organisms. Dual RNA-seq, the simultaneous detection of host and pathogen transcripts, has become a leading approach to unravelling complex molecular interactions between the host and the pathogen and is particularly informative for intracellular organisms. The amount of in vitro and in vivo dual RNA-seq data is rapidly growing, which demands computational pipelines to effectively analyse such data. In particular, holistic, systems-level, and temporal analyses of dual RNA-seq data are essential to enable further insights into the host-pathogen transcriptional dynamics and potential interactions. Here, we developed an integrative network-driven bioinformatics pipeline, dRNASb, a systems biology-based computational pipeline to analyse temporal transcriptional clusters, incorporate molecular interaction networks (e.g. protein-protein interactions), identify topologically and functionally key transcripts in host and pathogen, and associate host and pathogen temporal transcriptome to decipher potential between-species interactions. The pipeline is applicable to various dual RNA-seq data from different species and experimental conditions. As a case study, we applied dRNASb to analyse temporal dual RNA-seq data of Salmonella-infected human cells, which enabled us to uncover genes contributing to the infection process and their potential functions and to identify putative associations between host and pathogen genes during infection. Overall, dRNASb has the potential to identify key genes involved in bacterial growth or host defence mechanisms for future uses as therapeutic targets.

Stress response and virulence factors in bacterial pathogens relevant for Chilean aquaculture: current status and outlook of our knowledge

The study of the stress responses in bacteria has given us a wealth of information regarding the mechanisms employed by these bacteria in aggressive or even non-optimal living conditions. This information has been applied by several researchers to identify molecular targets related to pathogeny, virulence, and survival, among others, and to design new prophylactic or therapeutic strategies against them. In this study, our knowledge of these mechanisms has been summarized with emphasis on some aquatic pathogenic bacteria of relevance to the health and productive aspects of Chilean salmon farming (Piscirickettsia salmonis, Tenacibaculum spp., Renibacterium salmoninarum, and Yersinia ruckeri). This study will aid further investigations aimed at shedding more light on possible lines of action for these pathogens in the coming years.

Why Does Piscirickettsia salmonis Break the Immunological Paradigm in Farmed Salmon? Biological Context to Understand the Relative Control of Piscirickettsiosis

Piscirickettsiosis (SRS) has been the most important infectious disease in Chilean salmon farming since the 1980s. It was one of the first to be described, and to date, it continues to be the main infectious cause of mortality. How can we better understand the epidemiological situation of SRS? The catch-all answer is that the Chilean salmon farming industry must fight year after year against a multifactorial disease, and apparently only the environment in Chile seems to favor the presence and persistence of Piscirickettsia salmonis. This is a fastidious, facultative intracellular bacterium that replicates in the host’s own immune cells and antigen-presenting cells and evades the adaptive cell-mediated immune response, which is why the existing vaccines are not effective in controlling it. Therefore, the Chilean salmon farming industry uses a lot of antibiotics—to control SRS—because otherwise, fish health and welfare would be significantly impaired, and a significantly higher volume of biomass would be lost per year. How can the ever-present risk of negative consequences of antibiotic use in salmon farming be balanced with the productive and economic viability of an animal production industry, as well as with the care of the aquatic environment and public health and with the sustainability of the industry? The answer that is easy, but no less true, is that we must know the enemy and how it interacts with its host. Much knowledge has been generated using this line of inquiry, however it remains insufficient. Considering the state-of-the-art summarized in this review, it can be stated that, from the point of view of fish immunology and vaccinology, we are quite far from reaching an effective and long-term solution for the control of SRS. For this reason, the aim of this critical review is to comprehensively discuss the current knowledge on the interaction between the bacteria and the host to promote the generation of more and better measures for the prevention and control of SRS.

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.

Dual transcriptomic analysis reveals metabolic changes associated with differential persistence of human pathogenic bacteria in leaves of Arabidopsis and lettuce

Understanding the molecular determinants underlying the interaction between the leaf and human pathogenic bacteria is key to provide the foundation to develop science-based strategies to prevent or decrease the pathogen contamination of leafy greens. In this study, we conducted a dual RNA-sequencing analysis to simultaneously define changes in the transcriptomic profiles of the plant and the bacterium when they come in contact. We used an economically relevant vegetable crop, lettuce (Lactuca sativa L. cultivar Salinas), and a model plant, Arabidopsis thaliana Col-0, as well as two pathogenic bacterial strains that cause disease outbreaks associated with fresh produce, Escherichia coli O157:H7 and Salmonella enterica serovar Typhimurium 14028s (STm 14028s). We observed commonalities and specificities in the modulation of biological processes between Arabidopsis and lettuce and between O157:H7 and STm 14028s during early stages of the interaction. We detected a larger alteration of gene expression at the whole transcriptome level in lettuce and Arabidopsis at 24 h post inoculation with STm 14028s compared to that with O157:H7. In addition, bacterial transcriptomic adjustments were substantially larger in Arabidopsis than in lettuce. Bacterial transcriptome was affected at a larger extent in the first 4 h compared to the subsequent 20 h after inoculation. Overall, we gained valuable knowledge about the responses and counter-responses of both bacterial pathogen and plant host when these bacteria are residing in the leaf intercellular space. These findings and the public genomic resources generated in this study are valuable for additional data mining.

Loss of Dihydroxyacid Dehydratase Induces Auxotrophy in Bacillus anthracis

Anthrax disease is caused by infection with the bacteria Bacillus anthracis which, if left untreated, can result in fatal bacteremia and toxemia. Current treatment for infection requires prolonged administration of antibiotics. Despite this, inhalational and gastrointestinal anthrax still result in lethal disease. By identifying key metabolic steps that B. anthracis uses to grow in host-like environments, new targets for antibacterial strategies can be identified. Here, we report that the ilvD gene, which encodes dihydroxyacid dehydratase in the putative pathway for synthesizing branched chain amino acids, is necessary for B. anthracis to synthesize isoleucine de novo in an otherwise limiting microenvironment. We observed that ΔilvD B. anthracis cannot grow in media lacking isoleucine, but growth is restored when exogenous isoleucine is added. In addition, ΔilvD bacilli are unable to utilize human hemoglobin or serum albumin to overcome isoleucine auxotrophy, but can when provided with the murine forms. This species-specific effect is due to the lack of isoleucine in human hemoglobin. Furthermore, even when supplemented with physiological levels of human serum albumin, apotransferrin, fibrinogen, and IgG, the ilvD knockout strain grew poorly relative to nonsupplemented wild type. In addition, comparisons upon infecting humanized mice suggest that murine hemoglobin is a key source of isoleucine for both WT and ΔilvD bacilli. Further growth comparisons in murine and human blood show that the auxotrophy is detrimental for growth in human blood, not murine. This report identifies ilvD as necessary for isoleucine production in B. anthracis, and that it plays a key role in allowing the bacilli to effectively grow in isoleucine poor hosts. IMPORTANCE Anthrax disease, caused by B. anthracis, can cause lethal bacteremia and toxemia, even following treatment with antibiotics. This report identifies the ilvD gene, which encodes a dihydroxyacid dehydratase, as necessary for B. anthracis to synthesize the amino acid isoleucine in a nutrient-limiting environment, such as its mammalian host. The use of this strain further demonstrated a unique species-dependent utilization of hemoglobin as an exogenous source of extracellular isoleucine. By identifying mechanisms that B. anthracis uses to grow in host-like environments, new targets for therapeutic intervention are revealed.

Nutrient Scarcity in a New Defined Medium Reveals Metabolic Resistance to Antibiotics in the Fish Pathogen Piscirickettsia salmonis

Extensive use of antibiotics has been the primary treatment for the Salmonid Rickettsial Septicemia, a salmonid disease caused by the bacterium Piscirickettsia salmonis. Occurrence of antibiotic resistance has been explored in various P. salmonis isolates using different assays; however, P. salmonis is a nutritionally demanding intracellular facultative pathogen; thus, assessing its antibiotic susceptibility with standardized and validated protocols is essential. In this work, we studied the pathogen response to antibiotics using a genomic, a transcriptomic, and a phenotypic approach. A new defined medium (CMMAB) was developed based on a metabolic model of P. salmonis. CMMAB was formulated to increase bacterial growth in nutrient-limited conditions and to be suitable for performing antibiotic susceptibility tests. Antibiotic resistance was evaluated based on a comprehensive search of antibiotic resistance genes (ARGs) from P. salmonis genomes. Minimum inhibitory concentration assays were conducted to test the pathogen susceptibility to antibiotics from drug categories with predicted ARGs. In all tested P. salmonis strains, resistance to erythromycin, ampicillin, penicillin G, streptomycin, spectinomycin, polymyxin B, ceftazidime, and trimethoprim was medium-dependent, showing resistance to higher antibiotic concentrations in the CMMAB medium. The mechanism for antibiotic resistance to ampicillin in the defined medium was further explored and was proven to be associated to a decrease in the bacterial central metabolism, including the TCA cycle, the pentose-phosphate pathway, energy production, and nucleotide metabolism, and it was not associated with decreased growth rate of the bacterium or with the expression of any predicted ARG. Our results suggest that nutrient scarcity plays a role in the bacterial antibiotic resistance, protecting against the detrimental effects of antibiotics, and thus, we propose that P. salmonis exhibits a metabolic resistance to ampicillin when growing in a nutrient-limited medium.

RNA-Seq-Based Analysis of Cortisol-Induced Differential Gene Expression Associated with Piscirickettsia salmonis Infection in Rainbow Trout (Oncorhynchus mykiss) Myotubes

Salmonid rickettsial septicemia (SRS) is the major infectious disease of the Chilean salmonid aquaculture industry caused by Piscirickettsia salmonis. Intensive farming conditions generate stress and increased susceptibility to diseases, being skeletal muscle mainly affected. However, the interplay between pathogen infection and stress in muscle is poorly understood. In this study, we perform an RNA-seq analysis on rainbow trout myotubes that are pretreated for 3 h with cortisol (100 ng/mL) and then infected with P. salmonis strain LF-89 for 8 h (MOI 50). Twelve libraries are constructed from RNA samples (n = 3 per group) and sequenced on Illumina HiSeq 4000. A total of 704,979,454 high-quality reads are obtained, with 70.25% mapped against the reference genome. In silico DETs include 175 total genes-124 are upregulated and 51 are downregulated. GO enrichment analysis reveals highly impacted biological processes related to apoptosis, negative regulation of cell proliferation, and innate immune response. These results are validated by RT-qPCR of nine candidate transcripts. Furthermore, cortisol pretreatment significantly stimulated bacterial gene expression of ahpC and 23s compared to infection. In conclusion, for the first time, we describe a transcriptomic response of trout myotubes infected with P. salmonis by inducing apoptosis, downregulating cell proliferation, and intrinsic immune-like response that is differentially regulated by cortisol.

Transcriptomic analysis reveals a Piscirickettsia salmonis-induced early inflammatory response in rainbow trout skeletal muscle

Skeletal muscle is the most abundant tissue in teleosts and is essential for movement and metabolism. Recently, it has been described that skeletal muscle can express and secrete immune-related molecules during pathogen infection. However, the role of this tissue during infection is poorly understood. To determine the immunocompetence of fish skeletal muscle, juvenile rainbow trout (Oncorhynchus mykiss) were challenged with Piscirickettsia salmonis strain LF-89. P. salmonis is the etiological agent of piscirickettsiosis, a severe disease that has caused major economic losses in the aquaculture industry. This gram-negative bacterium produces a chronic systemic infection that involves several organs and tissues in salmonids. Using high-throughput RNA-seq, we found that 60 transcripts were upregulated in skeletal muscle, mostly associated with inflammatory response and positive regulation of interleukin-8 production. Conversely, 141 transcripts were downregulated in association with muscle filament sliding and actin filament-based movement. To validate these results, we performed in vitro experiments using rainbow trout myotubes. In myotubes coincubated with P. salmonis strain LF-89 at an MOI of 50, we found increased expression of the proinflammatory cytokine il1b and the pattern recognition receptor tlr5s 8 and 12 h after infection. These results demonstrated that fish skeletal muscle is an immunologically active organ that can implement an early immunological response against P. salmonis.

Exploring genetic resistance to infectious salmon anaemia virus in Atlantic salmon by genome-wide association and RNA sequencing

Background

Infectious Salmonid Anaemia Virus (ISAV) causes a notifiable disease that poses a large threat for Atlantic salmon (Salmo salar) aquaculture worldwide. There is no fully effective treatment or vaccine, and therefore selective breeding to increase resistance to ISAV is a promising avenue for disease prevention. Genomic selection and potentially genome editing can be applied to enhance host resistance, and these approaches benefit from improved knowledge of the genetic and functional basis of the target trait. The aim of this study was to characterise the genetic architecture of resistance to ISAV in a commercial Atlantic salmon population and study its underlying functional genomic basis using RNA Sequencing.

Results

A total of 2833 Atlantic salmon parr belonging to 194 families were exposed to ISAV in a cohabitation challenge in which cumulative mortality reached 63% over 55 days. A total of 1353 animals were genotyped using a 55 K SNP array, and the estimate of heritability for the trait of binary survival was 0.13–0.33 (pedigree-genomic). A genome-wide association analysis confirmed that resistance to ISAV was a polygenic trait, albeit a genomic region in chromosome Ssa13 was significantly associated with resistance and explained 3% of the genetic variance. RNA sequencing of the heart of 16 infected (7 and 14 days post infection) and 8 control fish highlighted 4927 and 2437 differentially expressed genes at 7 and 14 days post infection respectively. The complement and coagulation pathway was down-regulated in infected fish, while several metabolic pathways were up-regulated. The interferon pathway showed little evidence of up-regulation at 7 days post infection but was mildly activated at 14 days, suggesting a potential crosstalk between host and virus. Comparison of the transcriptomic response of fish with high and low breeding values for resistance highlighted TRIM25 as being up-regulated in resistant fish.

Conclusions

ISAV resistance shows moderate heritability with a polygenic architecture, but a significant QTL was detected on chromosome 13. A mild up-regulation of the interferon pathway characterises the response to the virus in heart samples from this population of Atlantic salmon, and candidate genes showing differential expression between samples with high and low breeding values for resistance were identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07671-6.

Dual RNAseq highlights the kinetics of skin microbiome and fish host responsiveness to bacterial infection

Background

Tenacibaculum maritimum is a fish pathogen known for causing serious damage to a broad range of wild and farmed marine fish populations worldwide. The recently sequenced genome of T. maritimum strain NCIMB 2154^T provided unprecedented information on the possible molecular mechanisms involved in the virulence of this species. However, little is known about the dynamic of infection in vivo, and information is lacking on both the intrinsic host response (gene expression) and its associated microbiota. Here, we applied complementary omic approaches, including dual RNAseq and 16S rRNA gene metabarcoding sequencing using Nanopore and short-read Illumina technologies to unravel the host–pathogen interplay in an experimental infection system using the tropical fish Platax orbicularis as model.

Results

We showed that the infection of the host is characterised by an enhancement of functions associated with antibiotic and glucans catabolism functions but a reduction of sulfate assimilation process in T. maritimum. The fish host concurrently displays a large panel of immune effectors, notably involving innate response and triggering acute inflammatory response. In addition, our results suggest that fish activate an adaptive immune response visible through the stimulation of T-helper cells, Th17, with congruent reduction of Th2 and T-regulatory cells. Fish were, however, largely sensitive to infection, and less than 25% survived after 96 hpi. These surviving fish showed no evidence of stress (cortisol levels) or significant difference in microbiome diversity compared with controls at the same sampling time. The presence of T. maritimum in resistant fish skin and the total absence of any skin lesions suggest that these fish did not escape contact with the pathogen, but rather that some mechanisms prevented pathogens entry. In resistant individuals, we detected up-regulation of specific immune-related genes differentiating resistant individuals from controls at 96 hpi, which suggests a possible genomic basis of resistance, although no genetic variation in coding regions was found.

Conclusion

Here we focus in detail on the interplay between common fish pathogens and host immune response during experimental infection. We further highlight key actors of defence response, pathogenicity and possible genomic bases of fish resistance to T. maritimum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00097-1.

Chimeric Protein IPath® with Chelating Activity Improves Atlantic Salmon's Immunity against Infectious Diseases

Infection processes displayed by pathogens require the acquisition of essential inorganic nutrients and trace elements from the host to survive and proliferate. Without a doubt, iron is a crucial trace metal for all living organisms and also a pivotal component in the host–parasite interactions. In particular, the host reduces the iron available to face the infectious disease, increasing iron transport proteins’ expression and activating the heme synthesis and degradation pathways. Moreover, recent findings have suggested that iron metabolism modulation in fish promotes the immune response by reducing cellular iron toxicity. We hypothesized that recombinant proteins related to iron metabolism could modulate the fish’s immune system through iron metabolism and iron-responsive genes. Here a chimeric iron transport protein (IPath^®) was bioinformatically designed and then expressed in a recombinant bacterial system. The IPath^® protein showed a significant chelating activity under in vitro conditions and biological activity. Taking this evidence, a vaccine candidate based on IPath^® was evaluated in Atlantic salmon challenged with three different fish pathogens. Experimental trials were conducted using two fish groups: one immunized with IPath^® and another injected with adjutant as the control group. After 400 accumulated thermal units (ATUs), two different infection trials were performed. In the first one, fish were infected with the bacterium Aeromonas salmonicida, and in a second trial, fish were exposed to the ectoparasite Caligus rogercresseyi and subsequently infected with the intracellular bacterium Piscirickettsia salmonis. Fish immunized with IPath^® showed a significant delay in the mortality curve in response to A. salmonicida and P. salmonis infections. However, no significant differences between infected and control fish groups were observed at the end of the experiment. Notably, sea lice burden reduction was observed in vaccinated Atlantic salmon. Transcriptional analysis evidenced a high modulation of iron-homeostasis-related genes in fish vaccinated with IPath^® compared to the control group during the infection. Moreover, increasing expression of Atlantic salmon IgT was associated with IPath^® immunization. This study provides evidence that the IPath^® protein could be used as an antigen or booster in commercial fish vaccines, improving the immune response against relevant pathogens for salmon aquaculture.

In vivo transcriptome analysis provides insights into host-dependent expression of virulence factors by Yersinia entomophaga MH96, during infection of Galleria mellonella

The function of microbes can be inferred from knowledge of genes specifically expressed in natural environments. Here, we report the in vivo transcriptome of the entomopathogenic bacterium Yersinia entomophaga MH96, captured during initial, septicemic, and pre-cadaveric stages of intrahemocoelic infection in Galleria mellonella. A total of 1285 genes were significantly upregulated by MH96 during infection; 829 genes responded to in vivo conditions during at least one stage of infection, 289 responded during two stages of infection, and 167 transcripts responded throughout all three stages of infection compared to in vitro conditions at equivalent cell densities. Genes upregulated during the earliest infection stage included components of the insecticidal toxin complex Yen-TC (chi1, chi2, and yenC1), genes for rearrangement hotspot element containing protein yenC3, cytolethal distending toxin cdtAB, and vegetative insecticidal toxin vip2. Genes more highly expressed throughout the infection cycle included the putative heat-stable enterotoxin yenT and three adhesins (usher-chaperone fimbria, filamentous hemagglutinin, and an AidA-like secreted adhesin). Clustering and functional enrichment of gene expression data also revealed expression of genes encoding type III and VI secretion system-associated effectors. Together these data provide insight into the pathobiology of MH96 and serve as an important resource supporting efforts to identify novel insecticidal agents.

Global Proteomic Profiling of Piscirickettsia salmonis and Salmon Macrophage-Like Cells during Intracellular Infection

Piscirickettsiasalmonis is an intracellular bacterial fish pathogen that causes piscirickettsiosis, a disease with numerous negative impacts in the Chilean salmon farming industry. Although transcriptomic studies of P. salmonis and its host have been performed, dual host–pathogen proteomic approaches during infection are still missing. Considering that gene expression does not always correspond with observed phenotype, and bacteriological culture studies inadequately reflect infection conditions, to improve the existing knowledge for the pathogenicity of P. salmonis, we present here a global proteomic profiling of Salmon salar macrophage-like cell cultures infected with P. salmonis LF-89. The proteomic analyses identified several P. salmonis proteins from two temporally different stages of macrophages infection, some of them related to key functions for bacterial survival in other intracellular pathogens. Metabolic differences were observed in early-stage infection bacteria, compared to late-stage infections. Virulence factors related to membrane, lipopolysaccharide (LPS) and surface component modifications, cell motility, toxins, and secretion systems also varied between the infection stages. Pilus proteins, beta-hemolysin, and the type VI secretion system (T6SS) were characteristic of the early-infection stage, while fimbria, upregulation of 10 toxins or effector proteins, and the Dot/Icm type IV secretion system (T4SS) were representative of the late-infection stage bacteria. Previously described virulence-related genes in P. salmonis plasmids were identified by proteomic assays during infection in SHK-1 cells, accompanied by an increase of mobile-related elements. By comparing the infected and un-infected proteome of SHK-1 cells, we observed changes in cellular and redox homeostasis; innate immune response; microtubules and actin cytoskeleton organization and dynamics; alteration in phagosome components, iron transport, and metabolism; and amino acids, nucleoside, and nucleotide metabolism, together with an overall energy and ATP production alteration. Our global proteomic profiling and the current knowledge of the P. salmonis infection process allowed us to propose a model of the macrophage–P. salmonis interaction.

Transcriptomic Profiling of the Adaptive and Innate Immune Responses of Atlantic Salmon to Renibacterium salmoninarum Infection

Bacterial Kidney Disease (BKD), which is caused by a Gram-positive, intracellular bacterial pathogen (Renibacterium salmoninarum), affects salmonids including Atlantic salmon (Salmo salar). However, the transcriptome response of Atlantic salmon to BKD remained unknown before the current study. We used a 44K salmonid microarray platform to characterise the global gene expression response of Atlantic salmon to BKD. Fish (~54 g) were injected with a dose of R. salmoninarum (H-2 strain, 2 × 10⁸ CFU per fish) or sterile medium (control), and then head kidney samples were collected at 13 days post-infection/injection (dpi). Firstly, infection levels of individuals were determined through quantifying the R. salmoninarum level by RNA-based TaqMan qPCR assays. Thereafter, based on the qPCR results for infection level, fish (n = 5) that showed no (control), higher (H-BKD), or lower (L-BKD) infection level at 13 dpi were subjected to microarray analyses. We identified 6,766 and 7,729 differentially expressed probes in the H-BKD and L-BKD groups, respectively. There were 357 probes responsive to the infection level (H-BKD vs. L-BKD). Several adaptive and innate immune processes were dysregulated in R. salmoninarum-infected Atlantic salmon. Adaptive immune pathways associated with lymphocyte differentiation and activation (e.g., lymphocyte chemotaxis, T-cell activation, and immunoglobulin secretion), as well as antigen-presenting cell functions, were shown to be differentially regulated in response to BKD. The infection level-responsive transcripts were related to several mechanisms such as the JAK-STAT signalling pathway, B-cell differentiation and interleukin-1 responses. Sixty-five microarray-identified transcripts were subjected to qPCR validation, and they showed the same fold-change direction as microarray results. The qPCR-validated transcripts studied herein play putative roles in various immune processes including pathogen recognition (e.g., tlr5), antibacterial activity (e.g., hamp and camp), regulation of immune responses (e.g., tnfrsf11b and socs1), T-/B-cell differentiation (e.g., ccl4, irf1 and ccr5), T-cell functions (e.g., rnf144a, il13ra1b and tnfrsf6b), and antigen-presenting cell functions (e.g., fcgr1). The present study revealed diverse immune mechanisms dysregulated by R. salmoninarum in Atlantic salmon, and enhanced the current understanding of Atlantic salmon response to BKD. The identified biomarker genes can be used for future studies on improving the resistance of Atlantic salmon to BKD.

Pharmacological iron-chelation as an assisted nutritional immunity strategy against Piscirickettsia salmonis infection

Salmonid Rickettsial Septicaemia (SRS), caused by Piscirickettsia salmonis, is a severe bacterial disease in the Chilean salmon farming industry. Vaccines and antibiotics are the current strategies to fight SRS; however, the high frequency of new epizootic events confirms the need to develop new strategies to combat this disease. An innovative opportunity is perturbing the host pathways used by the microorganisms to replicate inside host cells through host-directed antimicrobial drugs (HDAD). Iron is a critical nutrient for P. salmonis infection; hence, the use of iron-chelators becomes an excellent alternative to be used as HDAD. The aim of this work was to use the iron chelator Deferiprone (DFP) as HDAD to treat SRS. Here, we describe the protective effect of the iron chelator DFP over P. salmonis infections at non-antibiotic concentrations, in bacterial challenges both in vitro and in vivo. At the cellular level, our results indicate that DFP reduced the intracellular iron content by 33.1% and P. salmonis relative load during bacterial infections by 78%. These findings were recapitulated in fish, where DFP reduced the mortality of rainbow trout challenged with P. salmonis in 34.9% compared to the non-treated group. This is the first report of the protective capacity of an iron chelator against infection in fish, becoming a potential effective host-directed therapy for SRS and other animals against ferrophilic pathogens.

Dual RNA-seq provides novel insight into the roles of dksA from Pseudomonas plecoglossicida in pathogen-host interactions with large yellow croakers ( Larimichthys crocea)

Pseudomonas plecoglossicida is a rod-shaped, gram-negative bacterium with flagella. It causes visceral white spot disease and high mortality in Larimichthys crocea during culture, resulting in serious economic loss. Analysis of transcriptome and quantitative real-time polymerase chain reaction (PCR) data showed that dksA gene expression was significantly up-regulated after 48 h of infection with Epinephelus coioides (log 2FC=3.12, P<0.001). RNAi of five shRNAs significantly reduced the expression of dksA in P. plecoglossicida, and the optimal silencing efficiency was 96.23%. Compared with wild-type strains, the symptoms of visceral white spot disease in L. crocea infected with RNAi strains were reduced, with time of death delayed by 48 h and mortality reduced by 25%. The dksA silencing led to a substantial down-regulation in cellular component-, flagellum-, and ribosome assembly-related genes in P. plecoglossicida, and the significant up-regulation of fliC may be a way in which virulence is maintained in P. plecoglossicida. The GO and KEGG results showed that RNAi strain infection in L. crocea led to the down-regulation of inflammatory factor genes in immune-related pathways, which were associated with multiple immune response processes. Results also showed that dksA was a virulence gene in P. plecoglossicida. Compared with the wild-type strains, RNAi strain infection induced a weaker immune response in L. crocea.

Next-generation technologies for studying host-pathogen interactions: a focus on dual transcriptomics, CRISPR/Cas9 screening and organs-on-chips

Pathogens constantly interact with their hosts and the environment, and therefore have evolved unique virulence mechanisms to target and breach host defense barriers and manipulate host immune response to establish an infection. Advances in technologies that allow genome mining, gene editing such as CRISPR/Cas9, genomic, epigenomic and transcriptomic studies such as dual RNA-seq, coupled with bioinformatics, have accelerated the field of host-pathogen interactions within a broad range of infection models. Underpinning of the molecular changes that accompany invasion of eukaryotic cells with pathogenic microorganisms at the intersection of host, pathogen and their local environment has provided a better understanding of infectious disease mechanisms and antimicrobial strategies. The recent evolution of physiologically relevant three-dimensional (3-D) tissue/organ models and microfluidic organ-on-chip devices also provided a window to a more predictive framework of infectious disease processes. These approaches combined hold the potential to highly impact discovery of novel drug targets and vaccine candidates of the future. Here, we review three of the available and emerging technologies-dual RNA-seq, CRISPR/Cas9 screening and organs-on-chips, applicable to the high throughput study and deciphering of interaction networks between pathogens and their hosts that are critical for the development of novel therapeutics.

Transcriptomic Changes of Piscirickettsia salmonis During Intracellular Growth in a Salmon Macrophage-Like Cell Line

Piscirickettsia salmonis is the causative agent of Piscirickettsiosis, a systemic infection of salmonid fish species. P. salmonis infects and survives in its host cell, a process that correlates with the expression of virulence factors including components of the type IVB secretion system. To gain further insights into the cellular and molecular mechanism behind the adaptive response of P. salmonis during host infection, we established an in vitro model of infection using the SHK-1 cell line from Atlantic salmon head kidney. The results indicated that in comparison to uninfected SHK-1 cells, infection significantly decreased cell viability after 10 days along with a significant increment of P. salmonis genome equivalents. At that time, the intracellular bacteria were localized within a spacious cytoplasmic vacuole. By using a whole-genome microarray of P. salmonis LF-89, the transcriptome of this bacterium was examined during intracellular growth in the SHK-1 cell line and exponential growth in broth. Transcriptome analysis revealed a global shutdown of translation during P. salmonis intracellular growth and suggested an induction of the stringent response. Accordingly, key genes of the stringent response pathway were up-regulated during intracellular growth as well as at stationary phase bacteria, suggesting a role of the stringent response on bacterial virulence. Our results also reinforce the participation of the Dot/Icm type IVB secretion system during P. salmonis infection and reveals many unexplored genes with potential roles in the adaptation to intracellular growth. Finally, we proposed that intracellular P. salmonis alternates between a replicative phase and a stationary phase in which the stringent response is activated.

Determination of an Interaction Network between an Extracellular Bacterial Pathogen and the Human Host

Dual RNA sequencing (RNA-seq) offers the promise of determining an interactome at a transcriptional level between a bacterium and the host but has yet to be done on any bacterial infection in human tissue. We performed dual RNA-seq and metabolomics analyses on wounded and infected sites following experimental infection of the arm with H. ducreyi. Our results suggest that H. ducreyi survives in an abscess by utilizing l-ascorbate as an alternative carbon source, possibly taking advantage of host ascorbic acid recycling, and that H. ducreyi also adapts by upregulating genes involved in anaerobic metabolism and inorganic ion and nutrient transport. To our knowledge, this is the first description of an interaction network between a bacterium and the human host at a site of infection.

A major gap in understanding infectious diseases is the lack of information about molecular interaction networks between pathogens and the human host. Haemophilus ducreyi causes the genital ulcer disease chancroid in adults and is a leading cause of cutaneous ulcers in children in the tropics. We developed a model in which human volunteers are infected on the upper arm with H. ducreyi until they develop pustules. To define the H. ducreyi and human interactome, we determined bacterial and host transcriptomic and host metabolomic changes in pustules. We found that in vivoH. ducreyi transcripts were distinct from those in the inocula, as were host transcripts in pustule and wounded control sites. Many of the upregulated H. ducreyi genes were found to be involved in ascorbic acid and anaerobic metabolism and inorganic ion/nutrient transport. The top 20 significantly expressed human pathways showed that all were involved in immune responses. We generated a bipartite network for interactions between host and bacterial gene transcription; multiple positively correlated networks contained H. ducreyi genes involved in anaerobic metabolism and host genes involved with the immune response. Metabolomic studies showed that pustule and wounded samples had different metabolite compositions; the top ion pathway involved ascorbate and aldarate metabolism, which correlated with the H. ducreyi transcriptional response and upregulation of host genes involved in ascorbic acid recycling. These data show that an interactome exists between H. ducreyi and the human host and suggest that H. ducreyi exploits the metabolic niche created by the host immune response.