Cyprinid herpesvirus 2 infection changes microbiota and metabolites in the gibel carp (Carassius auratus gibelio) midgut

Abstract

Cyprinid herpesvirus 2 (CyHV-2) infects gibel carp (Carassius auratus gibelio) and causes severe losses. Microbiota in animal guts involves nutrition intake, development, immunity, and disease resistance. However, the relationship between gibel carp gut microbiota and CyHV-2 infection is not well known. Herein, we analyzed the gut microbiota composition and metabolite profiles in CyHV-2-infected and -uninfected fish using high-throughput sequencing and gas chromatography/mass spectrometry. Results showed that CyHV-2 infection significantly changed gut microbiota and metabolite profiles (p < 0.05). High-throughput sequencing demonstrated that the relative abundance of Aeromonas in the midgut increased dramatically while Cetobacterium decreased. Time-course analysis showed that the number of Aeromonas in the midgut of infected fish increased more than 1,000 times within 5 days post infection. Metabolome analysis illustrated that CyHV-2 infection significantly altered 24 metabolites in the midgut of gibel carp, annotating to the anomaly of digestion and metabolisms of amino acids, carbohydrates, and lipids, such as tryptophan (Trp) metabolism. The Mantel test demonstrated that gut microbiota and metabolite profiles were well related (r = 0.89). Furthermore, Trp metabolism responded to CyHV-2 infection closely was taken as one example to prove the correlation among CyHV-2 infection, metabolites and microbiota in the midgut, and host immunity. Results showed that modulating Trp metabolism could affect the relative abundance of Aeromonas in the midgut of fish, transcription of antiviral cytokines, and CyHV-2 infection. Therefore, we can conclude that CyHV-2 infection significantly perturbed the gut microbiome, disrupted its’ metabolic functions, and caused the proliferation of the opportunistic pathogen Aeromonas. This study also suggests that modulation of the gut microbiome will open a therapeutic opportunity to control CyHV-2 infection in gibel carp.

1 Introduction

Gibel carp (Carassius auratus gibelio), an omnivorous freshwater fish, is one of the leading aquaculture species worldwide. More than 2,700,000 tons of gibel carps were produced in 2020 with an economic value of over 5 billion dollars (Ministry of Agriculture and Rural Affairs of China, 2021), and most of its production farms are located in Yancheng, Province Jiangsu, China (Liu et al., 2018). However, the disease caused by Cyprinid herpesvirus 2 (CyHV-2), also known as the goldfish hematopoietic necrosis virus (GFHNV), limits the sustainable aquaculture of gibel carp (Lu et al., 2022; Zhu et al., 2022). CyHV-2 infection in gibel carp was first reported in 1999 in China (Chang et al., 1999), and the massive mortality reached over 90% in many breeding ponds (Wang et al., 2012). The CyHV-2-infected fish displays lethargic, anorexic, and widespread bleeding in the body, gill, and swim bladder. When the gill cover opens and closes or the fish’s body is stressed or jumping, blood will flow out from the gills (some infected fish may not have gill-bleeding symptoms). Some infected fish have erythema on the gill cover, which is called a beauty spot. Internal organs are congested or bleeding, the fish’s swim bladder has obvious bleeding points compared to the hemorrhages caused by Aeromonas spp., and the tail fin is whitish (Wu et al., 2013). CyHV-2 can be naturally detected in the water and sediments in the culture ponds, and gibel carp is sensitive to this virus (Yao et al., 2016). Gills are the critical organ for the initial invasion of CyHV-2 into fish, while the kidney and spleen are the primary organs for viral replication (Ding et al., 2014; Zhu et al., 2015; Yao et al., 2016). However, there is still no effective strategy to control this infection, which threatens this species’ sustainable culture industry.

Recent studies reveal that modulating the gut microbiota is an effective strategy to control disease development in human beings (Liu et al., 2013), but few studies were reported on the fish (Ran et al., 2021). Gut microbiota is considered as a good indicator of animal health, and plays a vital role in regulating host immune homeostasis and metabolic function (Guarner and Malagelada, 2003; Hu, 2017). The following four roles related to gut microflora have been reported: 1) metabolic activities in the fermentation of non-digestible dietary residue and endogenous mucus: the salvage of energy as short-chain fatty acids, the production of vitamin K, and the absorption of ions; 2) critical trophic effects on intestinal epithelial cell growth and differentiation, and interactions between intestinal bacteria and host immunity at the mucosal interface; 3) protection or mucous integrity-dependent barrier against pathogen invasion by alien microbes; and 4) relation to certain pathological disorders, including multisystem organ failure and inflammatory diseases (Guarner and Malagelada, 2003; Shapira, 2016; Cong et al., 2022). Thus, determining causal links of gut microbiota to diseases can contribute to developing new therapeutic interventions (Guarner and Malagelada, 2003; Li et al., 2017; Agus et al., 2018).

Although the gut microbiota of fish is affected by many factors, including host factors (e.g., genetics, gender, weight, and immunity), environmental factors (e.g., water, salinity, diet, and medicine), and microbial factors (Wang et al., 2018), the composition of the gut core microbiota is relatively stable in fish (Wang et al., 2018). Among them, water, diet, and feeding habits have been studied the most. Aeromonas, Pseudomonas, and Bacteroides type A dominate the gut microbiota of freshwater fish species. Plesiomonas, Enterobacteriaceae, Micrococcus, Acinetobacter, Clostridium, Bacteroides type B, and Fusariumas are less abundant. Additionally, many studies have used the combined metabolomic and metagenomic methods to investigate the correlation between metabolic alterations and gut microbiota changes, providing new insights into the disease pathology and therapeutic intervention (Guarner and Malagelada, 2003; Wang et al., 2018). For gibel carp, the gut microbiota of fish post-CyHV-2 infection was changed, and Plesiomonas was highly abundant in the infected fish Rong et al., 2017. CyHV-2 and CyHV-3 infections ruin the gut integrity of carps. These findings raised the possibility of the direct interaction between CyHV-2 and commensal bacteria (Rong et al., 2017; Ran et al., 2021). However, few studies have been conducted to investigate the effect of CyHV-2 infection on the gut microbial ecology by using combined metabolomic and metagenomic analysis. It will hinder the control and prevention development of CyHV-2 infection in gibel carp culture.

Thus, in this study, the changes of microbiota and metabolites in the midgut of CyHV-2-infected gibel carp and the correlation among them were firstly investigated. Then, the time-dependent number changes of target bacteria associated with CyHV-2 infection were analyzed through the culture method. Furthermore, tryptophan (Trp) metabolism was taken as one of the examples to prove the correlation between CyHV-2 infection and the metabolites since 1) the changes of microbiota and metabolites in the fish midgut post-CyHV-2 infection in this study had been annotated to Trp metabolism; 2) Trp metabolism has been reported to play a crucial role in gut microbiota–host cross-talk in health and disease via the aryl hydrocarbon receptor (AhR) (Hao and Whitelaw, 2013; Agus et al., 2018; Natividad et al., 2018; Sun et al., 2019; Wyant and Moslehi, 2022); 3) AhR, a nuclear protein upon association with certain endogenous and exogenous ligands, can translocate signals into the nucleus, bind DNA, and regulate gene expression. Trp metabolites are one of the most important endogenous AhR ligands (Ghiboub et al., 2020; Willem et al., 2022); and 4) Trp metabolism is known to regulate the intestinal barrier function and host immunity, and Trp metabolites/AhR signaling modulation are a therapeutic perspective in human diseases (Lamas et al., 2016; Lamas et al., 2018; Trikha and Lee, 2019; Platten et al., 2019; Ghiboub et al., 2020; Rosser et al., 2020). Taken together, these results about microbiome-modulated metabolites from this study will provide valuable information to control CyHV-2 infection from a disrupted equilibrium to clinical opportunities.

2 Materials and methods

2.1 Ethics statement

In the present study, all the animal experiments were handled in accordance with the guidelines for the Care and Use of Experimental Animals of China. The Committee for the Welfare and Ethics of Laboratory Animals was approved by the Animal Care and Use Committee of the Center for Applied Aquatic Genomics at the Chinese Academy of Fishery Sciences.

2.2 Fish acclimation

Gibel carp (43.17 ± 2.9 g) was obtained from a fish farm at Yancheng City, Jiangsu province China. After bringing the fish to the lab, we first performed PCR to make sure that the fish did not carry the CyHV-2 and then did further research (Figure S1). The primers for PCR are listed in Table S1. PCR in the reaction volume of 20 μl was performed following the protocols as one cycle of 95°C for 5 min, followed by 30 cycles at 95°C for 30 s, 58°C for 30 s, 72°C for 30 s, and final elongation at 72°C for 5 min (Zhu et al., 2015). The PCR products were checked by 1% agarose gel electrophoresis.

Then, the fish were acclimated in a tank with static water at a water temperature of 23°C–25°C, pH of 7.8–8.2, and dissolved oxygen (DO) higher than 5 mg L⁻¹ and provided with 24-h continuous aeration circulation for 2 weeks prior to the experiments. Gibel carp was fed three times each day with commercial diets (Tongwei, Yancheng, China) to satiation. The diets contained 32.25% crude protein, 5.90% crude lipid, 1.18% calcium, and 1.23% total phosphorus.

2.3 CyHV-2 solution preparation and viral quantification

The CyHV-2 solution was prepared as the previous description by Qiao et al. (2020). In brief, the kidney from naturally CyHV-2-infected gibel carp with typical signs of hemorrhage in the gill and swim bladder was sampled, homogenized by mixing with 10 times volume of phosphate-buffered saline (PBS, 0.01 M, pH7.2), centrifuged at 12,000 × g for 30 min and filtered by a 0.45/0.22-μm biofilter to remove bacteria.

DNA from the kidney was extracted using the viral DNA extraction kit (TaKaRa BIO Co., Ltd., Dalian, China), and the extracted DNA was assessed and adjusted to 50 ng μl⁻¹ by Biophotometer Plus (Eppendorf). The primers (CyHV-2-RT-F and CyHV-2-RT-R) for quantitative real-time PCR (qPCR) are listed in Table S1. The qPCR reactions were performed using the SYBR^® Premix ExTaq™ kit (Takara Bio Co., Ltd., Dalian, China) in a total volume of 25 μl, containing 12.5 μl SYBR Premix Ex Taq II (TaKaRa, Japan), 2 μl template DNA, 1 μl each primer and 8.5 μl ddH₂O. Distilled water replacing the extracted DNA was set as a negative control for each new run. The thermal cycling condition was as follows: one cycle of 95°C for 3 min, followed by 40 cycles at 95°C for 5 s and 55°C for 30 s and then by dissociation curve analysis (65°C–95°C: increment 0.5°C for 5 s) (Qiao et al., 2020). After amplification, melting curve analysis was performed, and Ct values were obtained. The standard curve of the viral DNA copy number is shown in Figure S2, which was generated as described by Xu et al. (2014). Plasmid DNA containing the sequence of the CyHV-2 helicase gene was selected to serve as the standard for virus quantification. The amplified DNA fragment was gel-purified using a Fermentas Gel Extraction Kit (Thermo). The resulting DNA fragment was inserted into the pMD19-T plasmid to produce pMD-CyHV-2. A 10-fold dilution series of pMD-CyHV-2 was used as the standard template of CyHV-2 in the qPCR. Finally, a viral stock solution at a concentration of 9.7 × 10⁷ copies μl^-1 was prepared.

2.4 Fish infection with CyHV-2 through gills (simulated to natural infection)

The gut microbiota diversity of gibel carp is related to the environment, diets, and feeding habits (Chen et al., 2021), but the composition of core microbiota in the gut is relatively stable under the same cultural circumstance (Wang et al., 2018; Chen et al., 2021). Therefore, allogynogenetic gibel carp, after acclimation for 14 days in the lab, was house-infected by CyHV-2 to avoid the effect of gender, water, and diet on gut microbiota. After acclimation in part 2.2, fish were divided randomly into two groups, namely, a control group with PBS and the CyHV-2-infection group. Each group included six tanks, three of them (volume of 100 L) with 8 individuals were used to calculate cumulative mortality (CM) and another three (volume of 400 L) with 30 individuals was set to sample for gut microbiota and metabolite analysis.

The ‘per-gill infection method’ mocked natural infection, which directly exposes virus to only gills, had been used to infect gibel carp like KHV and CHNV (Miyazaki et al., 2008; Somamoto et al., 2015) since gills are one of the first organs of aquatic animals exposed to pathogens (Murray et al., 2017; Andrews et al., 2010) and key invasion portal for CyHV-2 in gibel carp (Ding et al., 2014). Briefly, fish were anesthetized with tricaine methane sulfonate (MS-222) at 100 mg L⁻¹ and inoculated with a CyHV-2 working solution (9.7×10⁶ copies μl⁻¹) at 10 μl fish^-1 into their gills (both sides) in air. An equal volume of PBS was inoculated into gills as the control group. After gill exposure to the viral solution or PBS, fish were wrapped with wet papers and kept in air for 5–7 min at 25°C to allow the virus to adsorb into the gill tissue. Then, the fish were returned to the tank and maintained at 25°C with 24-h continuous aeration circulation. During the challenge test, no diets were fed (Qiao et al., 2020).

2.5 CyHV-2 infection confirmation

2.5.1 Clinical symptom observation and cumulative mortality

Clinical symptoms were observed daily, focusing on lethargy, anorexia, and hemorrhage in the body, eye, gill, and swim bladder. CM within 14 days post infection (dpi) was counted.

2.5.2 CyHV-2 load measurement

The fish with typical hemorrhagic symptoms at the 5th dpi were sampled. The viral load at 5 dpi in different tissues, including the gill, foregut, midgut, hindgut, liver, kidney, spleen, brain, and muscle, was detected by qPCR as described in part 2.3. Viral load was converted into copies per ng DNA.

2.6 Microbiome analysis of fish gut post-CyHV-2 infection

2.6.1 Midgut microbiota analysis through high-throughput sequencing

The midgut of three individuals with typical hemorrhagic symptoms at 5 dpi from each tank was collected and pooled as one sample. Resultantly, a total of nine fish for each group and three biological replicates were conducted for gut microbiome analysis. Microbial DNA was extracted from the midgut using the HiPure Soil DNA kits (Magen, Guangzhou, China) according to the manufacturer’s protocols. DNA samples were amplified to target the V4 region of bacterial 16S rRNA (Qiao et al., 2020). PCR products were purified with Qiagen Gel Extraction Kit (Qiagen, Suzhou, China). The primers 341F and 806R were used, and sequences are listed in Table S1 (Guo et al., 2017). Purified amplicons were pooled in the equimolar and paired-end sequences on an Illumina platform (MiSeq, Illumina Inc., Guangzhou, China) (2 × 250) following the standard protocols. The effective tags were clustered into operational taxonomic units (OTUs) of ≥97% similarity using the UPARSE (version 9.2.64) pipeline (Edgar, 2013). The tag sequence with the highest abundance was selected as a representative sequence within each cluster. A Venn diagram using the R package (version 1.6.16) to analyze and draw plots was performed in R project UpSetR package (version 1.3.3) to identify unique and common OTUs (Conway et al., 2017). The species comparison between groups was calculated by Welch’s t-test. Principal component analysis (PCA) was performed to distinguish the gut microbiota profiles between control and CyHV-2-infection groups. All α-diversity indexes were calculated in QIIME (version 1.9.1), and the alpha index comparison between groups was calculated by Welch’s t-test (Caporaso et al., 2010). Sequence alignment was performed using muscle (version 3.8.31), and phylogenetic tree was constructed using FastTree (version 2.1) (Edgar, 2004; Price et al., 2010); then, the weighted UniFrac distance matrix was generated by R project GuniFrac package (version 1.0) (Lozupone and Knight, 2006). The Kyoto Encyclopedia of Genes and Genomes KEGG pathway analysis of the OTUs was inferred using PICRUSt (version 2.1.4), and the analysis of functional difference between groups was calculated by Welch’s t-test (Langille et al., 2013). Canonical correspondence analysis (CCA) was executed in the R Project Vegan package (version 2.5.3) to clarify the influence of treatment on community composition (Oksanen et al., 2011). The biomarker features in each group were screened by the CCA effect size (LEfSe) software (version 1.0) (Edgar, 2013). Student’s t-test was used to evaluate the statistical differences of biological parameters between control and CyHV-2-infection groups. The significant difference was set at p < 0.05.

2.6.2 Verification of time-dependent Aeromonas abundance within 5 dpi by the culture method

Based on the results from part 2.6.1 and Aeromonas frequent infection in gibel carp culture practice, the time-course abundance of Aeromonas in the midgut of fish after infected by CyHV-2 was detected. The midgut of fish at 0, 1, 4, and 5 dpi was collected and homogenized in precold PBS. Then, the homogenate was diluted and spread on the RS selective medium, which was cultured at 28°C for 24–48 h. Bacterial colonies were counted and converted into the CFU g^-1 midgut. Meanwhile, the number of Aeromonas in culture water during the whole experiment was counted.

2.7 Metabolomic analysis of fish gut post-CyHV-2 infection

2.7.1 Metabolomic analysis of fish midgut by Jiangsu Liquid Chromatography-Tandem Mass Spectrometry LC-MS/MS

The midgut from each fish with typical hemorrhagic symptoms at 5 dpi was sampled, cut off, and then flushed with PBS to collect metabolites. Six biological replicates were set for each group. LC-MS/MS analysis was performed using a Ultra High Performance Liquid Chromatography UHPLC system (1290, Agilent Technologies). Mass Spectrometer MS raw data (.raw) files were converted to the mzML format using ProteoWizard and processed by R package XCMS (version 3.2), including retention time alignment, peak detection, and peak matching. The OSI-SMMS (version 1.0, Dalian Chem Data Solution Information Technology Co. Ltd.) was used for peak annotation after data processing with an in-house MS/MS database. Partial least squares discriminant analysis (PLS-DA) is a supervised dimensionality reduction method and was applied in comparison groups using R package models (http://www.r-project.org/) (Worley and Powers, 2013). The metabolite resonances were identified according to the Human Metabolome Database (Deng et al., 2014). Significantly changed metabolites were identified based on the following criteria: fold change >2.00 or <0.56, VIP > 1, and p < 0.05.

2.7.2 Detection of time-dependent concentration of aryl hydrocarbon receptor (AhR) within 5 days post-CyHV-2 infection

The results from parts 2.6 and 2.7.1 showed that 1) CyHV-2 infection increased the relative abundance of Aeromonas spp. and decreased Cetobacterium in the midgut significantly; 2) the concentration of indoleacetaldehyde decreased significantly; and 3) a strong positive correlation between them was also found. Additionally, indole derivatives such as indoleacetaldehyde are the products of Trp metabolism and activate the AhR signaling pathway and its downstream signals (Agus et al., 2018). AhR activation has been reported to play important roles in host health and diseases through modulating T-cell differentiation, the expression of cytokines, transcription factors, autoimmunity, and inflammation (Quintana et al., 2008; Kerkvliet et al., 2009; Apetoh et al., 2010; Jeremy et al., 2016). Thus, the time-dependent concentration of AhR in the midgut of infected fish at 0, 6, 24, 72, and 120 hpi was detected using the AhR assay kit (Jiangsu Jianglai Biotechnology Co. Ltd., Suzhou, China) according to the manufacturers’ protocol.

2.7.3 Effect verification of AhR expression on CyHV-2 infection

6-Formylindolo[3,2-b]carbazole (FICZ), a Trp photoproduct postulated as a candidate physiological ligand of AhR (Rannug et al., 1987; Oberg et al., 2005; Jeong et al., 2012), was used to evaluate the effect of AhR expression on the bacterial diversity in the midgut, the expression of antiviral cytokines, and viral infection progress.

2.7.3.1 Optimal FICZ working concentration assessment in gibel carp

FICZ (MedChemExpress, State of New Jersey USA) at the concentration of 1 μg fish^-1 was injected by I.P., and an equal volume of PBS (100 μl fish^-1) was injected as control. At 0, 12, 24, and 48 hpi, the fish midgut from two groups was sampled, immersed in RNA later (Sigma, Shanghai, China), and immediately stored at -80°C. Total RNA extraction, complementary DNA cDNA synthesis, and RT-qPCR were conducted as described in our previous study (Qiao et al., 2020). Briefly, RNA was extracted with the RNeasy mini kit (Qiagen, Valencia, CA, USA), and genomic DNA was erased through DNase I (Qiagen). A Nanodrop ND-1000 spectrophotometer was used to detect the quantity of RNA and concentration. Then, cDNA was synthesized using the oligo^dT primer and Prime-script™ first-strand cDNA synthesis kit (Takara Bio, Dalian, China). Genes AhR1 and AhR2 and downstream gene Cyp1A1 (Cytochrome P450, family 1, member A1) were chosen to investigate the effect of FICZ on the expression of key genes involved in the AhR signaling pathway (Jeong et al., 2012). The gene-specific primers are listed in Table S1. The RT-qPCR was performed through SYBR^® Premix Ex Taq™ (Takara, Dalian, China) in a CFX96 Real-Time PCR Detection System (Bio-Rad, Shanghai, China) (Qiao et al., 2020). The relative fold changes of a specific gene in the FICZ injection group were compared to that in the control PBS-injected group using the 2^-ΔΔCt method (Livak and Schmittgen, 2001).

2.7.3.2 Effect of FICZ on the expression of cytokines

FICZ was injected at 1 μg fish^-1 by I.P., and an equal volume of PBS (100 μl fish^-1) injected as control. At 12 hpi based on the results from part 2.7.3.1, the head kidney was sampled, immersed in RNA later, and immediately stored at -80°C. Total RNA extraction, cDNA synthesis, RT-qPCR, and the relative expression calculation of cytokines, including IFN-γ (interferon), MX1 (myxovirus resistance 1), ISG15 (interferon-stimulated gene 15), JAK (Janus kinase), TNF-α (tumor necrosis factor-α), MAPK3 (mitogen-activated protein kinase 3), IL-4 (interleukin-4), Gata3 (Gata binding protein 3), and T-Bet1 (T-box expressed in T cells 1), were performed as described in part 2.7.3.1. The genes and gene-specific primers are listed in Table S1.

2.7.3.3 Effect of FICZ on the cumulative mortality of fish post-CyHV-2 infection

There were four groups, namely, the PBS single-injection group (10 μl fish^-1), FICZ single-injection group, ‘FICZ+CyHV-2’-infection group, and ‘PBS+CyHV-2’- infection group. FICZ (1 μg fish^-1) was injected into fish by I.P., PBS, and CyHV-2 immersed fish through the ‘per-gill’ method as described in part 2.4. Additionally, two viral concentrations were set, which were 9.7×10⁶ and 2.4×10⁷ copies μl⁻¹. Each group included triplicates. Death within 14 days was recorded, and CM was calculated.

2.7.3.4 Effect of FICZ on the viral replication in vivo post-CyHV-2 infection

The fish from both the FICZ (1 μg fish^-1, I.P.) group and control group (PBS, 100 μl fish^-1) were infected by CyHV-2 through the ‘per-gill’ method as described in part 2.4. At 6, 12, 24, 48, 72, and 120 hpi, the kidney from ‘FICZ +CyHV-2-infected’ fish and ‘PBS +CyHV-2-infected’ fish was sampled, and viral load was quantified by qPCR as described in part 2.3.

2.8 Correlation between the gut microbiome and metabolites

Pearson correlation coefficient between gut microbial composition and experimental treatment was calculated with the R-project psych package (version 1.8.4) (Hornik, 2012). A network of correlation coefficients was generated using Omicsmart, a dynamic real-time interactive online platform for data analysis (http://www.omicsmart.com). In order to further prove the change of metabolites on the relative abundance change of main bacteria species, the effect of AhR expression on the relative abundance of Aeromonas in the midgut of fish post-CyHV-2 infection was investigated. Briefly, based on the results from part 2.7.3.1, the optimal working model of FICZ in vivo was injected by I.P. at 1 μg fish^-1 for 12 h. Four groups were set, including ‘PBS+ PBS’ single injection group (100 μl fish^-1), ‘FICZ+ PBS’ injection group, ‘FICZ+CyHV-2’ injection group, and ‘PBS+CyHV-2’ injection group. FICZ and PBS were injected into fish by I.P. and CyHV-2-infected fish through the ‘per-gill’ method as described in part 2.4. They were conducted at 12-h intervals. At 0, 1, 4, 5, and 7 dpi, the midgut of fish in four groups was sampled and homogenized in precold PBS. Then, the homogenate was diluted and spread on TSA and RS-selective medium, which was cultured at 28°C for 24–48 h. The TSA medium was used to count the total culturable bacteria, and the RS medium was used for Aeromonas spp. Bacterial colonies were counted and converted into the CFU g^-1 midgut.

2.9 Statistical analysis

Results are expressed as mean ± standard deviation and subjected to one-way analysis of variance using the statistical software program Statistical Product and Service Solutions SPSS version 17.0 (SPSS Inc., IL, USA). The significant difference was determined using Tukey’s multiple comparison test for more than three experimental groups and Student’s t-test for two experimental groups. Statistical significance was considered at P-values less than 0.05, and the results were expressed as mean ± SE (standard error).

3 Results

3.1 Viral load in vivo post-CyHV-2 infection through ‘per-gill’ method

No CyHV-2 was detected by PCR in the fish before the experiment (Figure S1A). Gibel carp could be infected by CyHV-2 through the ‘per-gill’ method in the lab. Like natural infection, hemorrhage in the fish skin and bulging eyes was observed at 5 dpi (Figure 1A). Severe hemorrhage was observed in the gills and swim bladder wall (Figures 1B, C). In the control group, no CyHV-2 was detected and no clinical symptoms were observed (Figures S1A, B). CM in the CyHV-2-infection (9.7 × 10⁷ copies per fish) group was 50%, and no mortality was observed in the PBS control group (Table S2). The standard curve of viral load using the qPCR method is shown in Figure S2. CyHV-2 loads among various tissues were significantly different (Figure 2). The viral loads in the kidney, spleen, midgut, hindgut, and gill were 4.4 × 10⁵, 3.5 × 10⁵, 3.0 × 10⁵, 1.6 × 10⁵, and 1.3 × 10⁵ copies ng^-1 DNA, respectively. Less than 1.0 × 10² copies ng^-1 DNA were detected in the liver, foregut, muscle, and brain. Thus, different organs had different viral loads.

Figure 1

Figure 2

3.2 CyHV-2 infection changed gut microbiome

3.2.1 OTU statistics and α-diversity analysis

We examined the microbiome in the midgut of infected and uninfected fish. The raw data about the microbiome and metabolites have been deposited in the NCBI (https://www.ncbi.nlm.nih.gov/sra/PRJNA655390), and the accession number is PRJNA655390. The abundance statistics of OTUs are shown in Table S3. Chao1 richness estimators in control and CyHV-2-infection groups were 926.60 and 985.54, respectively. The Simpson index in the CyHV-2-infection group was higher than that in the control group, and the Shannon index was lower. The less Simpson index, the higher the bacterial community diversity, while the Simpson index is inverse. The present results suggested that there was no significant difference in bacterial diversity between CyHV-2-infection and control groups (Table S4). Following the microbiota diversity analysis, we classified bacteria in fish midguts. A total of 383 OTUs were shared by the control and CyHV-2-infection groups. A total of 264 OTUs in control and 273 OTUs in the CyHV-2-infection group were unique, respectively (Figure S3).

3.2.2 Microbiota composition

At the phylum level, the microbiota in both control and CyHV-2-infection groups were dominated by Proteobacteria, Fusobacteria, Firmicutes, and Bacteroidetes (Figure 3A). Notably, compared to the control group, the relative abundance of Proteobacteria and Firmicutes in the CyHV-2-infection group was significantly increased, while the abundance of Fusobacteria decreased (p < 0.05). At the genus level, the gut microbiota in both control and CyHV-2-infection groups were dominated by Aeromonas, Cetobacterium, ZOR0006, Flavobacterium, Shewanella, Vibrio, Gemmobacter, Bacteroides, Acinetobacter, and Pseudomonas. Compared to the control group, the relative abundance of Aeromonas in the CyHV-2-infection group increased from 41.81% to 67.00%, while Cetobacterium decreased from 38.36% to 4.41% (Figure 3B). Multivariate statistical analysis–Principal Co-ordinates Analysis PCoA was used to reveal the difference in midgut microbiome patterns that responded to CyHV-2 infection. The gut microbiota in the CyHV-2-infection group was clearly separated from the control group, with 69.12% and 22.65% variation explained by PC1 and PC2 at the genus level, respectively (Figure S4A). The result from hierarchical clustering analysis with unweighted pair-group method with arithmetic means UPGMA was consistent with the PCoA plot, i.e., all OTUs in both control and CyHV-2-infection groups were clustered in their own cluster except sample 2 by CyHV-2 infection (Figure S4B). Additionally, the analysis of similarities ANOSIM analysis showed that the R-value between the control and CyHV-2-infection groups was 0.63, illustrating that the middle difference between the two groups was found (Figure S4C). These differences of midgut microbiome patterns in different groups can be directly attributed to the relative abundance changes of dominant species. Therefore, CyHV-2 infection changed the abundance of different classes of bacteria in the midgut of fish.

Figure 3

LEfSe and CCA analysis showed that the specific main species in each group could be found. The most relevant species in the midgut-related CyHV-2 infection was Aeromonas. In the control group, Fusobacteria and Cetobacterium among different individuals were donated (Figure 4). Aeromonas presented a closely positive association with the viral infection (Figure 5). Aeromonas accounted for 41.81% in the control group and 67.00% in the CyHV-2-infection group (Figure 3B). A time-dependent abundance of Aeromonas in the midgut of gibel carp post-CyHV-2 infection by the culture method increased from 3.8 × 10⁴ CFU g^-1 midgut at 0 dpi to 6.4 × 10⁵ CFU g^-1 midgut at 1 dpi, 2.8 × 10⁶ CFU g^-1 midgut at 4 dpi, and 7.6 × 10⁷ CFU g^-1 midgut at 5 dpi, respectively. The abundance of Aeromonas at 5 dpi increased almost 1,000 times compared to 0 hpi, which was in accordance with the microbiome analysis by high-throughput sequencing. There was no significant difference in Aeromonas (less than 8 CFU ml^-1 water) in culture water between the two groups during the whole challenge test. These results illustrated that CyHV-2 infection increased the relative abundance of Aeromonas and decreased Cetobacterium.

Figure 4

Figure 5

3.2.3 Functional prediction of gut microbiota through PICRUSt

A total of 35 KEGG functions of all the midgut microbiota from two groups were annotated through PICRUSt. The major functions were donated into three categories, namely, metabolism, genetic information processing, and environmental information processing. The functions regulated by the top 10 gut microbiota were membrane transport, carbohydrate metabolism, amino acid metabolism, replication and repair, energy metabolism, translation, metabolism of cofactors and vitamins, nucleotide metabolism, cell motility, and lipid metabolism (Figure S5A). The different microbiota between the control and CyHV-2-infection groups were annotated into 16 functions (Figure S5B).

3.3 CyHV-2 infection altered metabolite profiles in fish gut

3.3.1 Basic analysis of midgut metabolites from two groups

The PLS-DA model was conducted on the NMR data sets of all individuals, and the first two primary components PC1 and PC2 were selected to identify the discrimination between the control and CyHV-2-infection groups. Results showed that the metabolite profiles between the two groups were separated with 24% and 46% variation explained by PC1 and PC2, respectively (Figure 6).

Figure 6

3.3.2 Different metabolites between two groups

Compared to the control, the concentrations of the main 24 metabolites were significantly changed in the CyHV-2-infection group. Among them, the concentrations of 8 metabolites (norizalpinin, D-maltose, (dibromomethylene)-17beta-hydroxy-androst-4-en-3-one propiote, deoxy-D-altro-heptulose 7-phosphate, ellagic acid, alanyl-valine, N-a-acetyl-L-arginine, and phosphophosphite) increased and 16 metabolites decreased significantly (Table 1). These metabolites are related to various metabolism pathways, such as the digestive system (p-cresol, D-glucose, and D-maltose), amino acid metabolism (2-oxo-4-methylthiobutanoic acid, phenylpyruvic acid, hydroxyphenyllactic acid, ketoleucine, O-phospho-L-serine, picolinic acid, 2-oxoarginine, and 2-oxoadipate), metabolism of other amino acids (L-selenocysteine), carbohydrate metabolism (D-maltose), lipid metabolism (phytosphingosine), xenobiotics biodegradation, and metabolism (cyclophosphamide).

3.3.3 Function annotation of metabolites

All of the metabolites from both control and CyHV-2-infection groups were enriched, and the top 20 KEGG pathways are shown in Figure 7A. The first seven pathways were metabolic pathways, 2-oxocarboxilic acid metabolism, carbohydrate digestion and absorption, Trp metabolism, starch and sucrose metabolism, taste transduction, and biosynthesis of amino acids (Figure 7A). Different metabolites between the two groups were enriched in 11 signaling pathways (p < 0.05) (Figure 7B). Among them, the most obviously different signaling pathways were Trp metabolism, starch and sucrose metabolism, carbohydrate digestion and absorption, metabolic pathways, and taste transduction (Figure 7). Accordingly, different metabolites were mainly annotated to the digestive system, amino acid metabolism, carbohydrate metabolism, and lipid metabolism (Table 2). Taken together with metabolites and their related signaling pathways, we found that five metabolites were involved in these significantly different signaling pathways, including picolinic acid (PLA), 2-oxoadipate, and indoleacetaldehyde participated in Trp metabolism, phytosphingosine entangled in non-alcoholic fatty liver disease, D-maltose and D-glucose annotated in carbohydrate digestion and absorption, starch and sucrose metabolism, taste transduction, and ATP-binding cassette ABC transporters. Therefore, CyHV-2 infection changed midgut metabolites significantly, further affecting metabolisms, such as Trp metabolism.

Figure 7

3.4 Verification of the correlation between AhR and CyHV-2 infection

3.4.1 CyHV-2 infection decreased the concentration of AhR in the midgut

In the control group, there was no significant difference on the concentration of AhR among 0, 6, 24, 72, and 120 hpi. However, in the CyHV-2-infection group, the concentration decreased from 5.38 pM (0 hpi) to 3.77 pM (120 hpi) (Figure 8A).

Figure 8

3.4.2 AhR expression affected the transcription of cytokines

The optimal FICZ injection (I.P.) was 1 μg fish^-1 for 12 hpi (Figure 8B). At 12 hpi, the highest expression of AhR1 and Cyp1A1 was observed, and AhR2 upregulated significantly (Figure 8B). The transcriptional levels of cytokines JAK, IFN-γ, ISG15, MX1, IL-4, Gata3, and T-Bet1 were upregulated significantly by 9.75-fold, 9.92-fold, 7.05-fold, 5.56-fold, 2.03-fold, 4.14-fold, and 2.75-fold, respectively. However, the mRNA expression of genes TNF-α and MAPK3 was not regulated significantly (Figure 8C).

3.4.3 Higher AhR expression reduced the cumulative mortality of fish and inhibited CyHV-2 replication in vivo

Compared to the ‘PBS+CyHV-2’ infection group, the CM of fish post-CyHV-2 infection in the ‘FICZ+CyHV-2’ group was much lower. CM in the ‘FICZ+CyHV-2’ group at the viral concentrations of 2.4 × 10⁷ and 9.7×10⁶ copies μl^-1 was 41.7% and 33.3%, respectively. They were 83.8% and 58.3% in the ‘PBS+CyHV-2’ group. The fish in the ‘FICZ+CyHV-2’ group died later than that in the ‘PBS+CyHV-2’ group. No death was observed in PBS single and FICZ single injection groups (Table 3).

Compared to the ‘PBS+CyHV-2’-infection group, the viral replication in the ‘FICZ+CyHV-2’ group was slower, suggesting that higher AhR expression inhibited CyHV-2 replication in vivo. The viral load in the ‘FICZ+CyHV-2’ group at 6 hpi (6.9 × 10³ copies ng^-1 DNA) and 120 hpi (2.3 × 10⁴ copies ng^-1 DNA) was lower than that in the ‘PBS+CyHV-2’ group, which were 1.4 × 10⁴ copies ng^-1 DNA at 6 hpi and 5.5 × 10⁴ copies ng^-1 DNA at 120 hpi, respectively (Figure 8D). These findings and in-depth verification suggested that Trp metabolism is one of the cross-talk ways between a viral infection and a host.

3.5 Correlation between the altered microbiota and metabolites

In the CyHV-2-infection group, the relative abundance of Aeromonas significantly increased, while Cetobacterium decreased. We generated a scatter plot using Pearson correlation coefficient analysis to investigate the potential association between metabolites and bacterial composition. The heatmap of the correlation between gut metabolites and microbiota is shown in Figure 9. The apparent correlations were identified between the perturbed bacteria genus and altered metabolites (Mantel test, r = 0.89) (Figure 10 and Table 2). For example, the OTUs assigned to Aeromonas were positively correlated with glycyrrhetinic acid, sedoheptulose 1-phosphate, donepezil, and 2-oxoarginine and negatively correlated with indoleacetaldehyde, 1-kestose, and 8-iso-15-keto-PGE2. OTUs belonging to Cetobacterium were positively correlated with indoleacetaldehyde, adenosine monophosphate, 4-ketocyclophosphamide, LysoPC (24:0), D-glucono, cis-aconitatecis-4-carboxymethylenebut-2-en-4-olide, 4-pyridoxate, 4-hydroxybenzaldehyde, pyrrolidonecarboxylic acid, 13-oxoODE, L-cysteine, ketoleucine, stearoyl sphingomyelin, deoxyribose 5-phosphate, phytosphingosine, 2-oxoarginine and 4-hydroxyphenylpyruvic acid (Figure 10 and Table 2). In this study, the most noteworthy functional genes were involved in ABC transporters and glycerophospholipid metabolism related to membrane synthesis and gut integrity, followed by primary immunodeficiency (Figure S5B). Our data demonstrated that CyHV-2 infection changed both the gut microbiome and metabolome, and these two changes were correlated.

Figure 9

Figure 10

3.6 Alteration of AhR expression affected microbiota in the midgut of fish post CyHV-2 infection

Compared to the PBS injection group, CyHV-2-infection group, and ‘FICZ+CyHV-2’ group, FICZ stimulation at 4 and 5 dpi increased the number of total culturable bacteria by approximately two folds significantly. There was no significant difference among the PBS injection group, CyHV-2-infection group, and ‘FICZ+CyHV-2’-infection group (Figure 11A). However, compared to PBS injection group and FICZ injection group, the number of Aeromonas spp. in the CyHV-2-infection group and ‘FICZ+CyHV-2’- infection group from 1 to 7 dpi increased significantly. There was still a significant difference between the CyHV-2-infection group and the ‘FICZ+CyHV-2’ infection group from 1 to 7 dpi. The number of Aeromonas spp. in the CyHV-2-infection group was significantly higher than that in the ‘FICZ+CyHV-2’ infection group from 1 to 7 dpi. FICZ injection could inhibit the replication of Aeromonas spp. in the midgut of fish post-CyHV-2 infection (Figure 11B).

Figure 11

4 Discussion

Gibel carp can be artificially infected by CyHV-2 through the ‘per-gill’ method. Like natural infection, the typical clinical syndromes of hemorrhagic gills and the swim bladder were observed post-CyHV-2 infection through the ‘per-gill’ method, and CyHV-2 could be quantified in different inner organs. The hemorrhage of gills is mainly observed during the transition from spring to summer and in autumn, similar to the epidemiological characteristics of the hematopoietic necrosis (HVHN) of goldfish infected with CyHV-2 (Jung and Miyazaki, 1995; Zhu et al., 2015). Until now, there is no effective method to control this disease, although some studies on a vaccine are undergoing in the lab (Takafumi and Yukio, 2015; Yan et al., 2020). Early prevention before disease outbreaks becomes more important to control this epidemic disease. Growing pieces of evidence have demonstrated that the gut microbiota plays an important role in health and diseases. The gut microbiota in fish can regulate the expression of 212 genes, of which some are related to the promotion of nutrient metabolism, stimulation of epithelial proliferation, and immunity. The absence of gut microbiota in fish may lead to epithelial cell dysfunction and weaker immune responses (Wang et al., 2018). Colonization by commensal in newly hatched zebrafish (Denio rerio) primes neutrophils and induces several genes encoding proinflammatory and antiviral mediators (Galindo-Villegas et al., 2012).

To our knowledge, this study was the first to address the relationship among CyHV-2 infection, gut microbiota, and metabolite profiles in the gibel carp. The present results showed (1) gut microbiota dysbiosis, especially the relative abundance increment of Aeromonas and decrement of Cetobacterium; (2) metabolome alternation regarding amino acid metabolism, carbohydrate metabolism, and lipid metabolism, especially Trp metabolism; (3) the strong correlation between gut microbiota and metabolite profiles; and (4) the interference of viral infection through modulating gut ecology. The predominant genera in the gibel carp were Aeromonas, Cetobacterium, Flavobacterium, Shewanella, Vibrio, Gemmobacteria, Bacteroides, Acinetobacter, and Pseudomonas. Accordingly, Aeromonas, Cetobacterium, Pseudomonas, Acinetobacter, and Bacteroides were reported as the predominant genera in freshwater fish (Li et al., 2017; Wang et al., 2018). However, there are some differences in genus composition; Flavobacterium, Shewanella, and Vibiro were found in our experimental fish, while these genera are mainly reported in marine fish (Onarheim et al., 1994; Blanch et al., 1997). It is supposed that this difference is mainly attributed to the culture environment, especially salinity, since our experimental fish was obtained from Yancheng and the culture sediment in the pond was closer to the mudflat (Chen et al., 2021). The culture farm in Yancheng is mainly located along the coastline with the longest and biggest mudflat in China. Consistent to the previous report (Zhang et al., 2016), the abundance of dominant species such as Cetobacterium increases along with higher salinity.

The microbiota depends on a functional equilibrium related to the host, environment, and dietary factors. The gut microbiome is responsible for maintaining gut metabolic homeostasis and coordinating immune responses (Kau et al., 2011; Shapira, 2016). A disturbance in the functional equilibrium can lead to changes in the microbial diversity and abundance of certain bacteria, which are beneficial or harmful to fish (Liu et al., 2013; Ling et al., 2014; Lukens et al., 2014). Thus, we used a laboratory infection model instead of naturally infected fish to study the effect of CyHV-2 infection on the gut microbiota and microbiome. After acclimation, the gut microbiota is more stable (Chen et al., 2021). In this study, the midgut of fish was chosen as the target organ to study the microbiome and metabolite changes post-CyHV-2 infection since the viral load in the midgut was higher than that in the foregut and hindgut (Figure 2). Accordingly, the microbial abundance and gene gata5 expression in the midgut of zebrafish are higher than that in the foregut and hindgut (Hu, 2017; Rong et al., 2017). The gata5 gene regulates gut mucus secretion and gut epithelial differentiation (Hu, 2017).

CyHV-2 infection changed the gut microbiota of fish. The abundance of genera Aeromonas, Pseudomonas, and Flavobacterium were significantly increased in the CyHV-2-infected fish. Particularly, the relative abundance of Aeromonas in the CyHV-2-infected fish increased significantly from 48.74% to 67.00%, consistent with reports in crucian carp (C. carassius) (Caruso et al., 2017). Accordingly, a time-dependent abundance of Aeromonas in the midgut of gibel carp increased along with the CyHV-2 infection progress. However, it is different from the previous report (Rong et al., 2017), which showed that Plesiomonas was highly abundant in CyHV-2-infected gibel carp and could be used as a microbial biomarker for CyHV-2 infection. In this study, no dominant Plesiomonas was detected in both uninfected and infected gibel carp, as in the previous report by Chen et al. (2021). This difference might be related to the environment salinity, water, and developmental stage, which have been widely reported to affect the fish gut microbiota (Wang et al., 2018). This study suggests that the abundance of Aeromonas could be attractively monitored during the disease development progress. Aeromonas are ubiquitous inhabitants of freshwater and estuary environments as potential etiological agents in diseases of fish, terrestrial animals, and humans. It can cause motile aeromonad septicemia with high mortality and severe economic losses (Ran et al., 2018).

The gut metabolic disorders were observed in the CyHV-2-infected fish. For instance, three digestive pathways related to metabolites of p-cresol, D-glucose, and D-maltose and eight amino acid metabolism–related pathways to the metabolites of 2-oxo-4-methylthiobutanoic acid, phenylpyruvic acid, hydroxyphenyllactic acid, ketoleucine, O-phospho-L-serine, PLA, 2-oxoarginine, and 2-oxoadipate were changed in the gut of CyHV-2-infected fish. PLA has broad and effective antimicrobial activity against both bacteria and fungi (Dieuleveux et al., 1998; Lavermicocca et al., 2003). PLA is a major metabolic biomarker for oxidative damage to the cerebral cortex, phenylketonuria, and alcohol-induced liver disease (Manna et al., 2011). More recently, PLA has been well documented to play a positive role in the progression of different cancers (Lavermicocca et al., 2003; Fong et al., 2011) as shown in this study. Glucose as a good carbon source for cofactor regeneration in Trp metabolism is consistent with the phenomenon that the glucose content increased, while the corresponding maltose content decreased significantly in this study. D-maltose and D-glucose are involved in the digestive system. Sugars are the primary carbon source for bacterial cells, and they play important roles in bacterial pathogenesis (Pacheco et al., 2012; Ng et al., 2013). Both pathogenic and commensal bacteria compete for carbon sources such as glucuronate, mannose, fucose, and ribose to colonize and proliferate in the gut (Lustri et al., 2017). Sugar availability influences the microbiota composition and expansion of bacterial pathogens (Lustri et al., 2017). The change of glucose concentration in the environment regulates the expression of the locus of enterocyte effacement and catabolic, and osmotic stress (Njoroge et al., 2012; Njoroge et al., 2013).

In this study, an apparent correlation between specific microbes and metabolites was observed. For instance, Cetobacterium OTUs are negatively correlated with phytosphingosine, which is involved in membrane synthesis. The concentration increment of phytosphingosine raises the permeability of endothelial cells and leads to the barrier function decline (Luther et al., 2013). Another interesting metabolite is 1-kestose, the smallest component of Fructooligosaccharide FOS in the probiotics from the human gut. These findings suggest that 1-kestose is a potential new prebiotic targeting drug (Tochio et al., 2018). It has been shown that the hepatocyte membrane has abundant glycyrrhetinic acid (GA) receptors (Shan et al., 2017). Some documents have authenticated that GA is combined with the GA receptors of hepatocyte membrane to enhance the targeted therapy of liver (Shan et al., 2017; Wang et al., 2018). A more interesting finding is that Cetobacterium OTUs are correlated with indoleacetaldehyde positively (Natividad et al., 2018). Both indoleacetaldehyde and PLA from this study were annotated to Trp metabolism. Actually, most indole derivatives are formed from Trp metabolism (Lamas et al., 2018). Indole substances are the endogenous ligands of AhR and thus participate in many important biological processes of the body, such as cell differentiation, apoptosis, and inflammation (Agus et al., 2018). Indole and Trp metabolites are a major source of endogenous AhR ligand precursors (Bjeldanes et al., 1991; Perdew and Babbs, 1991). It has been reported that the gut microbiota affects health and disease through the regulation of Trp metabolism and the activation of AhR (Lamas et al., 2016; Lamas et al., 2018; Agus et al., 2018; Natividad et al., 2018; Trikha and Lee, 2019). CARD9 was reported to impact colitis by altering the gut microbiota metabolism of Trp into AhR ligands (Lamas et al., 2016).

AhR is a member of the basic helix–loop–helix–(bHLH) superfamily of transcription factors, which are associated with cellular responses to environmental stimuli. Early AhR studies focused on understanding the role of AhR in mediating the toxicity and carcinogenesis properties of the prototypic ligand 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Recently, AhR has been highly receptive to a wide array of endogenous and exogenous ligands. Its activation leads to a myriad of crucial host physiological functions, such as intestinal barrier function and immune cells, as well as intestinal homeostasis (Lamas et al., 2018). The byproducts of Trp photooxidation possess a high AhR-binding capacity, and are able to induce the expression of Cyp1A1 and other AhR target genes (Paine, 1976; Paine and Francis, 1980; Goerz et al., 1996). In fish, two kinds of AhR (AhR1 and AhR2) have been reported mostly (Prasch et al., 2003; Aluru et al., 2011). Thus, the AhR pathway related to Trp metabolism in this study was chosen to prove the effect of gut micro-ecology changes on the viral infection process through FICZ stimulation, which is a prime example of a Trp photoproduct and can agonistically stimulate AhR activity in as low as picomolar ranges (Jeong et al., 2012). The present results showed that 1) CyHV-2 infection reduced AhR concentration in the midgut (Figure 8A) and serum (data not shown), 2) higher AhR expression upregulated the mRNA expression of antiviral genes in the head kidney; 3) higher AhR expression reduced the relative abundance of Aeromonas spp. in the midgut of fish post-CyHV-2 infection; 4) higher AhR expression reduced the 14-day CM of fish post-CyHV-2 infection and inhibited viral replication in vivo. In mammals, AhR is known to express by different intestinal immune cells, such as intestinal epithelial cells (IECs), Th17 cells, innate lymphoid cells (ILCs), macrophages, and neutrophils. The AhR signaling pathway is pivotal to the differentiation and proliferation of T cells and the regulation of mucosal intestinal immune responses (Qiu et al., 2012; Goudot et al., 2017). Although the correlation between metabolite changes and CyHV-2 infection focusing on Trp metabolism was proven, further in-depth studies on the correlation among gut microbes, gut barrier disturbance, Aeromonas translocation, Trp metabolism, immunity, and antiviral infection need to be conducted. Especially, regulating Trp metabolism in the gut by diets including probiotics will attract more attention as reported in mammals (Whitfield-Cargile et al., 2016; Cervantes-Barragan et al., 2017). These actionable studies will probably provide an excellent immunotherapeutic intervention to control carp diseases and thus advance the sustainable culture.

5 Conclusion

This study found a significant difference in microbiota-regulated metabolites in the fish midgut post-CyHV-2 infection. Notably, the relative abundance of Aeromonas increased while Cetobacterium decreased. The different metabolites were annotated to the metabolism of the digestive system, amino acid, carbohydrate, lipid, and xenobiotic biodegradation, especially Trp metabolism. An in-depth analysis showed that some of these altered metabolites were highly correlated with the genera Aeromonas and Cetobacterium. Our results supported that CyHV-2 infection–induced changes in the midgut, bacterial community structure disturbed the metabolic functions of the gut microbiome, and the viral infection process could be interrupted through the modulation of gut microecology. These findings will provide a new therapeutic strategy to prevent disease and thus advance the sustainable culture.

References

• 1

  AgusA.PlanchaisJ.SokolH. (2018). Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe23, 716–724. doi: 10.1016/j.chom.2018.05.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AluruN.KarchnerS. I.HahnM. E. (2011). Role of DNA methylation of AHR1 and AHR2 promoters in differential sensitivity to PCBs in Atlantic killifish (Fundulus heteroclitus). Aquat. Toxicology.101, 288–294. doi: 10.1016/j.aquatox.2010.10.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AndrewsM.BattagleneS.CobcroftJ.AdamsM.NogaE.NowakB. (2010). Host response to the chondracanthid copepod (Chondracanthus goldsmidi), a gill parasite of the striped trumpeter Latris lineata (Forster), in Tasmania. J. Fish Diseases.33, 211–220. doi: 10.1111/j.1365-2761.2009.01107.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  ApetohL.QuintanaF. J.PotC.JollerN.XiaoS.KumarD.et al. (2010). The aryl hydrocarbon receptor interacts with c-maf to promote the differentiation of type 1 regulatory T cells induced by IL-27. Nat. Immunol.11, 854–861. doi: 10.1038/ni.1912

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BjeldanesL. F.KimJ. Y.GroseK. R.BartholomewJ. C.BradfieldC. A. (1991). Aromatic hydrocarbon responsiveness-receptor agonists generated from indole-3-carbinol in vitro and in vivo: comparisons with 2,3,7,8-tetrachlorodibenzo-p-dioxin. PNAS Proc. Natl. Acad. Sci. U.S.A.88, 9543–9547. doi: 10.1073/pnas.88.21.9543

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BlanchA. R.AlsinaM.SimonM.JofreJ. (1997). Determination of bacteria associated with reared turbot (Scophthalmus maximus) larvae. J. Appl. Microbiol.82, 729–734. doi: 10.1046/j.1365-2672.1997.00190.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods7, 335–336. doi: 10.1038/nmeth.f.303

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CarusoC.PastorinoP.BurioliE. A. V.PelettoS.RighettiM.BonaC.et al. (2017). CyHV-2 outbreak associated with aeromonas spp. in crucian carp (Carassius carassius) in piedmont (Italy)Bamberg, Banagiotis Berillis. 307–314.

  • Google Scholar
  • View reference in article

• 9

  Cervantes-BarraganL.ChaiJ. N.TianeroM. D.DiLucciaB.AhernP. P.MerrimanJ.et al. (2017). Lactobacillus reuteri induces gut intraepithelial CD4⁺ CD8αα⁺ T cells. Science357, 806–810. doi: 10.1126/science.aah5825

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  ChangP. H.Shu-HwaeL.ChiangH. C.JongM. H. (1999). Epizootic of herpes-like virus infection in goldfish (Carassius auratus) in Taiwan. Fish Pathology.34, 209–210.

  • Google Scholar
  • View reference in article

• 11

  ChenP.SunQ. R.ZhangH. H.WangS.WangJ. W.XuY.et al. (2021). Intestinal microbiota analysis in gibel carp (Carassius auratus gibelio) based on 16S rRNA gene sequence. Fisheries Science.8, 1–10. doi: 10.16378/j.cnki.1003-1111.20208

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  CongJ.ZhouP.ZhangR. Y. (2022). Intestinal microbiota-derived short chain fatty acids in host health and disease. Nutrients14:1977. doi: 10.3390/nu14091977

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  ConwayJ. R.LexA.GehlenborgN. (2017). UpSetR: an r package for the visualization of intersecting sets and their properties. Bioinformatics18, 2938–2940. doi: 10.1093/bioinformatics/btx364

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  DengY. F.ZhangY.ZhangR.WuB.DingL.XuK.et al. (2014). Mice in vivo toxicity studies for monohaloacetamides emerging disinfection byproducts based on metabolomic methods. Environ. Sci. Technology.48, 8212. doi: 10.1021/es502010v

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  DieuleveuxV.LemarinierS.GuéguenM. (1998). Antimicrobial spectrum and target site of d-3-phenyllactic acid. Int. J. Food Microbiol.40, 177–183. doi: 10.1016/S0168-1605(98)00031-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  DingZ. F.XiaS. Y.ZhaoZ. M.XiaA. J.ShenM. F.TangJ. Q.et al. (2014). Histopathological characterization and fluorescence in situ hybridization of cyprinid herpesvirus 2 in cultured Prussian carp (Carassius auratus gibelio) in China. J. Virological Methods206, 76–83. doi: 10.1016/j.jviromet.2014.05.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  EdgarR. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res.32, 1792–1797. doi: 10.2460/ajvr.69.1.82

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  EdgarR. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods10, 996–998. doi: 10.1038/nmeth.2604

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  FongM. Y.McDunnJ.KakarS. S. (2011). Identification of metabolites in the normal ovary and their transformation in primary and metastatic ovarian cancer. PloS One6, e19963. doi: 10.1371/journal.pone.0019963

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  Galindo-VillegasJ.Garcia-MorenoD.OliveiraS. D.MeseguerJ.MuleroV. (2012). Regulation of immunity and disease resistance by commensal microbes and chromatin modifications during zebrafish development. PNAS Proc. Natl. Acad. Sci. U.S.A.109 (39), E2605–E2614.

  • Google Scholar
  • View reference in article

• 21

  GhiboubM.VerburgtC. M.SovranB.BenningaM. A.LimbergenJ. (2020). Nutritional therapy to modulate tryptophan metabolism and aryl hydrocarbon-receptor signaling activation in human diseases. Nutrients12 (9), 2846. doi: 10.3390/nu12092846

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  GoerzG.BarnstorfW.WinnekendonkG.BolsenK.FritschC.KalkaK.et al. (1996). Influence of UVA and UVB irradiation on hepatic and cutaneous P450 isoenzymes. Arch. Dermatol. Res.289, 46–51. doi: 10.1007/s004030050151

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  GoudotC.CoillardA.VillaniA. C.GueguenP.CrosA.SarkizovaS.et al. (2017). Aryl hydrocarbon receptor controls monocyte differentiation into dendritic cells versus macrophages. Immunity47, 582–596. doi: 10.1016/j.immuni.2017.08.016

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  GuarnerF.MalageladaJ. R. (2003). Gut flora in health and disease. Lancet361, 512–519. doi: 10.1016/S0140-6736(03)12489-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  GuoM. J.WuF. H.HaoG. E.QiQ.LiR.LiN.et al. (2017). Bacillus subtilis improves immunity and disease resistance in rabbits. Front. Immunol.8. doi: 10.3389/fimmu.2017.00354

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  HaoN.WhitelawM. L. (2013). The emerging roles of AhR in physiology and immunity. Biochem. Pharmacol.86, 561–570. doi: 10.1016/j.bcp.2013.07.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  HornikK. (2012). The comprehensive r archive network Vol. 4 (Wiley Interdisciplinary Reviews-Computational Statistics), 394–398. doi: 10.1002/wics.1212

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  HuJ. (2017). The homeostasis of gut adhesive microbiota as a marker for evaluation of fish health. [dissertation’s thesis] (WuHan: Huazhong Agricultural University).

  • Google Scholar
  • View reference in article

• 29

  JeongK.-T.HwangS.-J.OhG.-S.ParkJ.-H. (2012). FICZ, a tryptophan photoproduct, suppresses pulmonary eosinophilia and Th2-type cytokine production in a mouse model of ovalbumin-induced allergic asthma. Int. Immunopharmacology.13, 377–385. doi: 10.1016/j.intimp.2012.04.014

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  JeremyA. G.RoopaliG.JessicaE. K.AdaY.GopalM.SharmilaS.et al. (2016). AHR activation is protective against colitis driven by t cells in humanized mice. Cell Rep.17, 1318–1329. doi: 10.1016/j.celrep.2016.09.082

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  JungS. J.MiyazakiT. (1995). Herpesviral haematopoietic necrosis of goldfish (Carassius auratus l.). J. Fish Diseases.18, 211–220. doi: 10.1111/j.1365-2761.1995.tb00296.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  KauA. L.AhernP. P.GriffinN. W.GoodmanA. L.GordonJ. I. (2011). Human nutrition, the gut microbiome and the immune system. Nature474, 327–336. doi: 10.1038/nature10213

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  KerkvlietN. I.SteppanL. B.VorachekW.OdaS.FarrerD.WongC. P.et al. (2009). Activation of aryl hydrocarbon receptor by TCDD prevents diabetes in NOD mice and increases Foxp3+ T cells in pancreatic lymph nodes. Immunotherapy1, 539–547. doi: 10.2217/IMT.09.24

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  LamasB.NatividadJ. M.SokolH. (2018). Aryl hydrocarbon receptor and intestinal immunity. Mucosal Immunol.11, 1024–1038. doi: 10.1038/s41385-018-0019-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  LamasB.RichardM. L.LeducqV.PhamH.-P.MichelM.-L.CostaG. D.et al. (2016). CARD9 impacts colitis by altering gut microbiota metabolism of tryptophan into aryl hydrocarbon receptor ligands. Nat. Med.22, 598–605. doi: 10.1038/nm.4102

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  LangilleM. G. I.ZaneveldJ.CaporasoJ. G.McdonaldD.KnightsD.ReyesJ. A.et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol.31, 814–821. doi: 10.1038/nbt.2676

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  LavermicoccaP.ValerioF.ViscontiA. (2003). Antifungal activity of phenyllactic acid against molds isolated from bakery products. Appl. Environ. Microbiol.69, 634–640. doi: 10.1128/AEM.69.1.634-640.2003

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  LiT. T.LongM.LiH.GatesoupeF. J.ZhangX. J.ZhangQ. Q.et al. (2017). Multi-omics analysis reveals a correlation between the host phylogeny, gut microbiota and metabolite profiles in cyprinid fishes. Front. Microbiol.8. doi: 10.3389/fmicb.2017.00454

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LingZ. X.LiuX.JiaX. Y.ChengY. W.LuoY. Q.YuanL.et al. (20147485). Impacts of infection with different toxigenic Clostridium difficile strains on faecal microbiota in children. Sci. Rep.4, 7485. doi: 10.1038/srep07485

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LiuZ. J.CaoA. T.CongY. Z. (2013). Microbiota regulation of inflammatory bowel disease and colorectal cancer. Semin. Cancer Biol.23, 543–552. doi: 10.1016/j.semcancer.2013.09.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  LiuW.WuJ. P.LiZ.DuanZ. Y.WenH. (2018). Effects of dietary coated protease on growth performance, feed utilization, nutrient apparent digestibility, intestinal and hepatopancreas structure in juvenile gibel carp (Carassius auratus gibelio). Aquaculture Nutr.24, 47–55. doi: 10.1111/anu.12531

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2^-△△CT method. Methods25, 402–408.

  • Google Scholar
  • View reference in article

• 43

  LozuponeC.KnightR. (2006). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol.71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  LukensJ. R.GurungP.VogelP.JohnsonG. R.CarterR. A.McGoldrickD. J.et al. (2014). Dietary modulation of the microbiome affects autoinflammatory disease. Nature516, 246–249. doi: 10.1038/nature13788

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  LustriB. C.SperandioV.MoreiraC. G. (2017). Bacterial chat: Intestinal metabolites and signals in host-microbiota-pathogen interactions. Infection Immunity.85, e00476–e00417. doi: 10.1128/IAI.00476-17

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  LuC. Y.TangR. Z.SuM. Z.ZouJ. X.LvL. Q. (2022). Induction of reactive oxygen species is necessary for efficient onset of cyprinid herpesvirus 2 replication: implications for novel antiviral strategy with antioxidants. Front. Microbiol.12. doi: 10.3389/fmicb.2021.792655

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  LutherJ. A.EnesJ.BirrenS. J. (2013). Neurotrophins regulate cholinergic synaptic transmission in cultured rat sympathetic neurons through a p75-dependent mechanism. J. Neurophysiology.109, 485–496. doi: 10.1152/jn.00076.2011

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  MannaS. K.PattersonA. D.YangQ.KrauszW. K.IdleJ. R.FornaceA. J.et al. (2011). UPLC-MS-based urine metabolomics reveals indole-3-lactic acid and phenyllactic acid as conserved biomarkers for alcohol-induced liver disease in the para-null mouse model. J. Proteome Res.10, 4120–4133. doi: 10.1021/pr200310s

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  Ministry of Agriculture and Rural Affairs of China (2021). China Fisheries statistical yearbook, volume 14: Freshwater aquaculture output in each province (Beijing: China Agriculture Press).

  • Google Scholar
  • View reference in article

• 50

  MiyazakiT.KuzuyaY.YasumotoS.YasudaM.KobayashiT. (2008). Histopathological and ultrastructural features of koi herpesvirus (KHV)-infected carp (Cyprinus carpio), and the morphology and morphogenesis of KHV. Dis. Aquat. Organisms.80, 1–11. doi: 10.3354/dao01929

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  MurrayH. M.LeggiadroC. T.DouglasS. E. (2017). Immunocytochemical localization of pleurocidin to the cytoplasmic granules of eosinophilic granular cells from the winter flounder gill. J. Fish Biol.70, 336–345. doi: 10.1111/j.1095-8649.2007.01452.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  NatividadJ. M.AgusA.PlanchaisJ.LamasB.JarryA. C.MartinR.et al. (2018). Impaired aryl hydrocarbon receptor ligand production by the gut microbiota is a key factor in metabolic syndrome. Cell Metab.28, 737–749. doi: 10.1016/j.cmet.2018.07.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  NgK. M.FerreyraJ. A.HigginbottomS. K.LynchJ. B.KashyapP. C.GopinathS.et al. (2013). Microbiota-liberated host sugars facilitate post-antibiotic expansion of enteric pathogens. Nature502, 96–99. doi: 10.1038/nature12503

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  NjorogeJ. W.GruberC.SperandioV. (2013). The interacting cra and kdpe regulators are involved in the expression of multiple virulence factors in enterohemorrhagic Escherichia coli. J. Bacteriology195, 2499–2508. doi: 10.1128/JB.02252-12

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  NjorogeJ. W.NguyenY.CurtisM. M.MoreiraC. G.SperandioV. (2012). Virulence meets metabolism: Cra and kdpe gene regulation in enterohemorrhagic Escherichia coli. Mbio3, e00280–e00212. doi: 10.1128/mBio.00280-12

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  ObergM.BerganderL.HåkanssonH.RannugU.RannugA. (2005). Identification of the tryptophan photoproduct 6-formylindolo[3,2-b]carbazole, in cell culture medium, as a factor that controls the background aryl hydrocarbon receptor activity. Toxicological Sci.85, 935–943. doi: 10.1016/j.inoche.2007.11.024

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  OksanenJ.BlanchetF. G.KindtR.LegendreP.O'haraR. B. (2011). R package version 1.17-6. the comprehensive r archive network website (Vegan: Community Ecology Package).

  • Google Scholar
  • View reference in article

• 58

  OnarheimA. M.WiikR.BurghardtJ.StackebrandtE. (1994). Characterization and identification of two Vibrio species indigenous to the intestine of fish in cold sea water; description of Vibrio iliopiscarius sp. nov. Systematic Appl. Microbiol.17, 370–379. doi: 10.1016/S0723-2020(11)80053-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  PachecoA. R.CurtisM. M.RitchieJ. M.MuneraD.WaldorM. K.MoreiraC. G.et al. (2012). Fucose sensing regulates bacterial intestinal colonization. Nature492, 113–117. doi: 10.1038/nature11623

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  PaineA. J. (1976). Induction of benzo[a]pyrene mono-oxygenase in liver cell culture by the photochemical generation of active oxygen species. evidence for the involvement of singlet oxygen and the formation of a stable inducing intermediate. Biochem. J.158, 109–117. doi: 10.1042/bj1580109

  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  PaineA. J.FrancisJ. E. (1980). The induction of benzo[a]pyrene-3-mono-oxygenase by singlet oxygen in liver cell culture is mediated by oxidation products of histidine. Chemico-Biological Interactions.30, 343–353. doi: 10.1016/0009-2797(80)90057-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  PerdewG. H.BabbsC. F. (1991). Production of ah receptor ligands in rat fecal suspensions containing tryptophan or indole-3-carbinol. Nutr. Cancer.16, 209–218. doi: 10.1080/01635589109514159

  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  PlattenM.NollenE. A. A.RöhrigU. F.FallarinoF.OpitzC. A. (2019). Tryptophan metabolism as a common therapeutic target in cancer, neurodegeneration and beyond. Nat. Rev. Drug Discovery.18, 379–401. doi: 10.1038/s41573-019-0016-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  PraschA. L.HirokiT.CarneyS. A.WuD.TakeoH.StegemanJ. J.et al. (2003). Aryl hydrocarbon receptor 2 mediates 2,3,7,8-tetrachlorodibenzo-p-dioxin developmental toxicity in zebrafish. Toxicological Sci.1, 138–150. doi: 10.1093/toxsci/kfg202

  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  PriceM. N.DehalP. S.ArkinA. P. (2010). Fast tree 2-approximately maximum-likelihood trees for large alignments. PloS One5, e9490.

  • Google Scholar
  • View reference in article

• 66

  QiaoG.ChenP.SunQ. R.ZhaLivakngM. M.ZhangJ. L.LiZ.et al. (2020). Poly-beta-hydroxybutyrate (PHB) in bioflocs alters intestinal microbial community structure, immune-related gene expression and early cyprinid herpesvirus 2 replication in gibel carp (Carassius auratus gibelio). Fish Shellfish Immunol.97, 72–82. doi: 10.1016/j.fsi.2019.12.045

  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  QiuJ.HellerJ.GuoX.ChenZ. M.FishK.FuY. X.et al. (2012). The aryl hydrocarbon receptor regulates gut immunity through modulation of innate lymphoid cells. Immunity36, 92–104. doi: 10.1016/j.immuni.2011.11.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 68

  QuintanaF. J.BassoA. S.IglesiasA. H.KornT.FarezM. F.BettelliE.et al. (2008). Control of t(reg) and T(H)17 cell differentiation by the aryl hydrocarbon receptor. Nature453, 65–71. doi: 10.1038/nature06880

  • CrossRef
  • Google Scholar
  • View reference in article

• 69

  RanC.LiY.MaX. F.XieY. D.XieM. X.ZhangY. T.et al. (2021). Interactions between commensal bacterial and viral infection: insights for viral disease control in farmed animals. Sci. China-Life Sci.64, 1–12. doi: 10.1007/S11427-020-1721-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 70

  RannugA.RannugU.RosenkranzH. S.WingvistL.WesterholmR.AgurellE.et al. (1987). Certain photooxidized derivatives of tryptophan bind with very high affinity to the ah receptor and are likely to be endogenous signal substances. J. Biol. Chem.262, 15422–15427. doi: 10.1016/S0021-9258(18)47743-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 71

  RanC.QinC. B.XieM. X.ZhangJ. X.LiJ.XieY. D.et al. (2018). Aeromonas veronii and aerolysin are important for the pathogenesis of motile aeromonad septicemia in cyprinid fish. Environ. Microbiol.20, 3442–3456. doi: 10.1111/1462-2920.14390

  • CrossRef
  • Google Scholar
  • View reference in article

• 72

  RongS.LiT. T.LuoD.LiJ. B.YinL. Y.LiH.et al. (2017). Changes in the intestinal microbiota of gibel carp (Carassius gibelio) associated with cyprinid herpesvirus 2 (CyHV-2) infection. Curr. Microbiol.74, 1130–1136. doi: 10.1007/s00284-017-1294-y

  • CrossRef
  • Google Scholar
  • View reference in article

• 73

  RosserE. C.PiperC.MateiD. E.BlairP. A.MauriC. (2020). Microbiota-derived metabolites suppress arthritis by amplifying aryl-hydrocarbon receptor activation in regulatory b cells. Cell Metab.31 (4), 837–851.e10. doi: 10.1016/j.cmet.2020.03.003

  • CrossRef
  • Google Scholar

• 74

  ShanF.LiuY. J.JiangH. Y.TongF. (2017). In vitro and in vivo protein release and anti-ischemic/reperfusion injury properties of bone morphogenetic protein-2-loaded glycyrrhetinic acid-poly (ethylene glycol)-b-poly(l-lysine) nanoparticles. Int. J. Nanomedicine.12, 7613–7625. doi: 10.2147/IJN.S146546

  • CrossRef
  • Google Scholar

• 75

  ShapiraM. (2016). Gut microbiotas and host evolution: scaling up symbiosis. Trends Ecol. Evolution.31, 539–549. doi: 10.1016/j.tree.2016.03.006

  • CrossRef
  • Google Scholar

• 76

  SomamotoT.MiuraY. H.NakanishiT.NakaoM. (2015). Local and systemic adaptive immune responses toward viral infection via gills in ginbuna crucian carp. Dev. Comp. Immunol.52, 81–87. doi: 10.1016/j.dci.2015.04.016

  • CrossRef
  • Google Scholar

• 77

  SunM.MaN.HeT.JohnstonL. J.MaX. (2019). Tryptophan (trp) modulates gut homeostasis via aryl hydrocarbon receptor (AhR). Crit. Rev. Food Sci. Nutr.1-9. doi: 10.1080/10408398.2019.1598334

  • CrossRef
  • Google Scholar

• 78

  TakafumiL.YukioM. (2015). Effect of booster shot and investigation of vaccination efficacy period against herpesviral haematopoietic necrosis (HVHN) in goldfish (Carassius auratus). Veterinary Microbiol.175, 139–144. doi: 10.1016/j.vetmic.2014.11.005

  • CrossRef
  • Google Scholar

• 79

  TochioT.KadotaY.TanakaT.KogaY. (2018). 1-kestose, the smallest fructooligosaccharide component, which efficiently stimulates Faecalibacterium prausnitzii as well as bifidobacteria in humans. Foods7, 140. doi: 10.3390/foods7090140

  • CrossRef
  • Google Scholar

• 80

  TrikhaP.LeeD. A. (2019). The role of ahr in transcriptional regulation of immune cell development and function. Biochim. Biophys. Acta (BBA) - Rev. Cancer.1873, 188335. doi: 10.1016/j.bbcan.2019.188335

  • CrossRef
  • Google Scholar

• 81

  WangL.HeJ. G.LiangL. G.ZhengX.JiaP.ShiX. J.et al. (2012). Mass mortality caused by cyprinid herpesvirus 2 (CyHV-2) in Prussian carp (Carassius gibelio) in China. Bull. Eur. Assoc. Fish Pathologists.32, 164–173.

  • Google Scholar

• 82

  WangA. R.RanC.RingoeE.ZhouZ. G. (2018). Progress in fish gastrointestinal microbiota research. Rev. Aquaculture.10, 626–640. doi: 10.1111/raq.12191

  • CrossRef
  • Google Scholar

• 83

  Whitfield-CargileC. M.CohenN. D.ChapkinR. S.WeeksB. R.DavidsonL. A.GoldsbyJ. S.et al. (2016). The microbiota-derived metabolite indole decreases mucosal inflammation and injury in a murine model of NSAID enteropathy. Gut Microbes7, 246–261. doi: 10.1080/19490976.2016.1156827

  • CrossRef
  • Google Scholar

• 84

  WillemM. D. V.HerbertT.MatthiasV. H.PatriceD. C. (2022). Gut microbiome and health: mechanistic insights. Gut71 (5), 1020–1032. doi: 10.1136/gutjnl-2021-326789

  • CrossRef
  • Google Scholar

• 85

  WorleyB.PowersR. (2013). Multivariate analysis in metabolomics. Curr. Microbiol.1, 92–107. doi: 10.2174/2213235X11301010092

  • CrossRef
  • Google Scholar

• 86

  WuT.DingZ. F.RenM.AnL.XiaoZ. Y.LiuP.et al. (2013). The histo- and ultra-pathological studies on a fatal disease of Prussian carp (Carassius gibelio) in mainland China associated with cyprinid herpesvirus 2 (CyHV-2). Aquaculture412-413, 8–13.

  • Google Scholar

• 87

  WyantG. A.MoslehiJ. (2022). Expanding the therapeutic world of tryptophan metabolism. Circulation24, 1799–1820. doi: 10.1161/CIRCULATIONAHA.122.059812

  • CrossRef
  • Google Scholar

• 88

  XuL. J.PodokP.XieJ.LiL. Q. (2014). Comparative analysis of differential gene expression in kidney tissues of moribund and surviving crucian carp (Carassius auratus gibelio) in response to cyprinid herpesvirus 2 infection. Arch. Virology.159 (8), 1961–1974. doi: 10.1007/s00705-014-2011-9

  • CrossRef
  • Google Scholar

• 89

  YanY. Y.HuoX. C.AiT. S.SuJ. G. (2020). β-glucan and anisodamine can enhance the immersion immune efficacy of inactivated cyprinid herpesvirus 2 vaccine in Carassius auratus gibelio. Fish Shellfish Immunol.98, 285–295. doi: 10.1016/j.fsi.2020.01.025

  • CrossRef
  • Google Scholar

• 90

  YaoZ. F.HuangJ.ZhangJ. M.NieH. H.ZhengX. H.LiL. J.et al. (2016). Investigation of viral load in healthy fish naturally infected with cyprinid herpesvirus II (CyHV-2). Hubei Agric. Sci.4, 980–983. doi: 10.14088/j.cnki.issn0439-8114.2016.04.041

  • CrossRef
  • Google Scholar

• 91

  ZhangM. L.SunY. H.LiuY. K.QiaoF.ChenL. Q.LiuW. T.et al. (2016). Response of gut microbiota to salinity change in two euryhaline aquatic animals with reverse salinity preference. Aquaculture454, 72–80. doi: 10.1016/j.aquaculture.2015.12.014

  • CrossRef
  • Google Scholar

• 92

  ZhuM.DaiY. P.TongX. Y.ZhangY. X.ZhouY.ChengJ. L.et al. (2022). Circ-udg derived from cyprinid herpesvirus 2 promotes viral replication. Microbiol. Spectrum.10, 2165–0497. doi: 10.1128/spectrum.00943-22

  • CrossRef
  • Google Scholar

• 93

  ZhuM.LiuB.CaoG. L.HuX. L.WeiY. H.YiJ. T.et al. (2015). Identification and rapid diagnosis of the pathogen responsible for haemorrhagic disease of the gill of allogynogenetic crucian carp. J. Virological Methods219, 67–74. doi: 10.1016/j.jviromet.2015.03.019

  • CrossRef
  • Google Scholar