- 1Key Open Laboratory of Cold Water Fish Germplasm Resources and Breeding of Heilongjiang Province, Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin, China
- 2College of Animal Science, Northeast Agricultural University, Harbin, China
- 3College of Animal Science, Guizhou University, Guiyang, China
- 4Fishery Technical Extension Station of Jilin Province, Changchun, China
This study aimed to demonstrate the effects of dietary glutathione (GSH) on growth, intestinal antioxidant capacity, histology, gene expression, and microbiota in juvenile triploid rainbow trout (Oncorhynchus mykiss). Different diets (G0-control, G100, G200, G400, and G800) containing graded levels of GSH (0, 100, 200, 400, and 800mgkg−1) were fed to triplicate groups of 30 fish (initial mean weight 4.12±0.04g) for 56days. G400 had significantly improved weight gain and feed conversion rate. Based on the broken-line regression analysis, the optimum dietary GSH level was 447.06mgkg−1. Catalase and superoxide dismutase activities were significantly higher in G200–G800. G200 had significantly lower malondialdehyde content. The height of the intestinal muscular layer in G400 was significantly higher than that of the control group. Intestinal PepT1 and SLC1A5 gene expression was significantly increased, and the highest was observed in G400. TNF-α, IL-1β, IL-2, and IL-8 expression were significantly decreased than that of G0. Next-generation sequencing of the 16S rDNA showed a significant difference in alpha diversity whereas no differences in beta diversity. On the genus level, LefSe analysis of indicator OTUs showed Ilumatobacter, Peptoniphilus, Limnobacter, Mizugakiibacter, Chelatococcus, Stella, Filimonas, and Streptosporangium were associated with the treatment diet, whereas Arcobacter, Ferrovibrio, Buchnera, Chitinophaga, Stenotrophobacter, Solimonadaceae, Polycyclovorans, Rhodococcus, Ramlibacter, and Azohydromonas were associated with the control diet. In summary, feeding juvenile triploid O. mykiss 200–800mgkg−1 GSH improved growth and intestinal health.
Introduction
Fish in intensive aquaculture are frequently subjected to a range of negative environmental stresses, including high temperatures, overcrowding, deteriorating water quality, pathogen invasion, and disinfectant and antibiotic (Ming et al., 2015). This causes bacterial resistance and residues in aquatic products, transfer of antimicrobial resistance genes into the environment and food, immunosuppression in fish and increased susceptibility to different diseases (Liu et al., 2010), which can lead to significant economic losses. As a consequence, one of the methods is to improve fish immune and stress resistance through dietary approaches.
Glutathione (GSH), a tripeptide consisting of glutamate, cysteine, and glycine residues (Kosower and Kosower, 1978) can eliminate unnecessary free oxygen radicals from cells. It also has many physiological roles, such as increasing antioxidant activity (Doyotte et al., 1997), protecting liver cells (Ponsoda et al., 1999), maintaining DNA synthesis, enhancing immunity (Will, 1999), and alleviating neuron intoxication (Raghunathan et al., 2007). Several GSH studies in freshwater fish and shrimp found growth, antioxidant ability, and immunity could be considerably improved in European bass (Dicentrarchus labrax; Zambonino-Infante et al., 1997), Nile tilapia (Oreochromis niloticus; Zhou et al., 2013), grass carp (Ctenopharyngodon idella; Ming et al., 2015), Atlantic salmon (Salmo salar; Ma et al., 2019), and whiteleg shrimp (Litopenaeus vannamei; Xia and Wu, 2018).
Rainbow trout (Oncorhynchus mykiss) is a typical coldwater fish. As one of China’s major farmed species, triploid O. mykiss is grown mainly in cold-water regions. Induced triploidy through chromosome set manipulation leads to an additional chromosome package in each somatic cell that makes fish sterile (Thorgaard, 1983). The negative effects of gonad growth with female triploid salmonids can often benefit recirculating aquaculture producers (Good and Davidson, 2016). Due to physiological changes, there may be variations in nutritional needs between ploidies (Fjelldal et al., 2016). This may be either because of increased growth potential or the odd number of chromosomes inherent in genetic variations (Ren et al., 2017). Currently, there is little information on the nutritional needs of triploid O. mykiss.
The intestine is a delicate tissue that plays an important role in fish health and nutrition. Maintaining intestinal homeostasis is important to enhance the growth performance and fish health status (Rombout et al. 2011). Intestinal health is determined by host (immunity, mucosal barrier), nutritional, microbial, and environmental factors (Kelly and Salinas, 2017). Thus, intestinal health could be improved through dietary approaches. GSH may benefit intestinal health in aquatic animals. In shrimp, dietary glutathione (150–250mgkg−1) improved the jejunum wall thickness and villus height of intestine (Wang et al., 2018). To date, there is limited published research on the intestinal health of dietary GSH in fish. Furthermore, the published data were mainly based on the purified diets. Further studies based on the practical diet are needed for the application of GSH or related ingredients in fish feeds, especially for the commercial fish feeds. Therefore, the aim of this research was to study the effects of GSH supplementation in practical diets on growth, intestinal antioxidant capacity, histology, gene expression, and microbiota in triploid O. mykiss diets.
Materials and Methods
Diets
Experimental diets (G0, G100, G200, G400, and G800) with five GSH levels (0, 100, 200, 400, and 800mgkg−1) were prepared (Table 1). The spectrophotometric approach (Ming et al., 2015) revealed that the dietary GSH levels were 8.52, 106.36, 210.32, 409.51, and 796.55mgkg−1, respectively. Fish meal, soybean protein concentrate, chicken meal, and extruded soybean were the main protein sources in the experimental diets, with fish oil and soybean oil as lipid sources and wheat middling as a carbohydrate source. Ingredients were finely ground before mixing (<250μm) and then blended with minerals and vitamins. After adding the lipid source, all the ingredients were thoroughly mixed for 15min before being mixed again for 10min. The dough was shaped into feed pellets with a diameter of 1.2mm using a small-scale extruder (G250; Machine Factory of Muyang, China). After pelleting, feed pellets were dried for approximately 12h in a ventilated oven at 45°C, then sieved. The pellets were then frozen at −20°C until use.
Feeding Management
Before feeding the basal diet, the fish were acclimatized in the laboratory (Chinese Academy of Fishery Sciences Coldwater Fish Experimental Station, Mudanjiang, China) for 14days. At the start of the experiment, the fish were starved for 24h before being pooled. In 15 tanks, 450 fish (initial weight 4.12±0.04g) were distributed (size: 300L). Each group had 3 tanks, and each tank was used as a replicate. The fish in each tank were weighed in batches. During the 56-day feeding trial, the fish were hand-fed four times a day (08:00, 11:00, 14:00, and 17:00) until they were satiated. Fish were raised in a water flow-through system (flow rate: 0.2Ls−1). Water quality was measured (YSI 6600 V2-2, Ohio State, United States) daily during the experimental period, water temperature (11.3–15.8°C), pH (7.2–7.5), dissolved oxygen (7.8–9.2mgL−1), and ammonia nitrogen (<0.2mgL−1). After a 24-h starvation period, the fish in each tank were weighed again as a batch at the end of the feeding trial.
Sample Collection
At the end of the experiment, all fish had been starved for 24h. Ten fish were randomly selected from each tank and anesthetized (tricaine methane sulfonate MS-222, 75mgL−1) before weighing. The weight (precision 0.01g) and fork length (precision 0.01cm) of the fish were then determined to calculate the condition factor (CF). Further, the total weight of the fish in each tank was measured to calculate the weight gain rate (WGR), specific growth rate (SGR), feed conversion rate (FCR), daily feed intake (DFI). For body composition analysis, four fish were sampled from each tank and stored at −80°C. Three fish were collected from each tank and aseptically killed in an ice bath. The body surfaces of the sampled fish were washed with 70% ethanol, and the fish were dissected using sterile surgical scissors. Then their mid-intestines (located right after the pyloric ceca, 1 inch from the stomach) were collected until being tested for antioxidant enzymes, and faecal samples from the mid-intestine were collected in sterile tubes from each dietary replicated tank and stored at −80°C for microbial analysis. The research protocol was handled following the Chinese Animal Health Protection Law and the Scientific Laboratory Animal Permit Approval (Ethical Approval No. SCXK(YU)2005-0001).
Chemical Analysis
Normal procedures were used to determine the crude protein, crude lipid, ash, moisture, and gross energy of feed and whole body (AOAC, 2012). After acid digestion using a Kjeltec system, crude protein (N×6.25) was measured using the Kjeldahl method (KDN-102C Autoanalyzer, Xianjian, China). The ether-extraction process was used to measure crude lipid with the Soxtec System (SXT-06-analyzer, Hongji, China). Moisture was measured by oven drying for 6h at 105°C. Ash was placed in a muffle furnace at 550°C for 12h. The energy content of the diet was determined using bomb calorimetry (XRY-1A, Jingmi, China).
Intestinal Antioxidant Capacity
Intestinal samples were weighed and mixed with an ice-cold buffer in a 1:9 ratio (0.86% normal saline). The extract was then centrifuged for 15min at 7,700g and 4°C, and the supernatant was used to determine superoxide dismutase (SOD), and catalase (CAT), glutathione reductase (GR) activities, GSH, and malondialdehyde (MDA) content. The standard was bovine serum albumin, and the spectrophotometer absorbance was determined at 750nm. Spectrophotometric kits were bought from the Chinese Nanjing Jiancheng Institute of Bioengineering and were used to analyze the MDA content (cat. no. A003-1), GSH content (cat. no. A006-1-1), and the activity of GR (cat. no. A062-1-1), SOD (cat. no. A001-1), and CAT (cat. no. A007-1). The concentrations of MDA, SOD, CAT, GR and GSH were then calculated according to the instructions provided with the respective kits as described by Deng et al. (2015).
Intestinal Gene Expression
Total RNA was extracted from mid-intestines according to the manufacturer’s instructions using RNAiso Plus Reagent (TaKaRa, Dalian, China). A spectrophotometer was used to examine the absorbance at 260nm to determine the RNA concentration. Using agarose gel electrophoresis, the RNA integrity was determined, and the absorbance ratio at A260 nm/A280 nm ranged from 1.8 to 2.0. Tumor necrosis factor (TNF-α), interleukin 1 (IL-1β), interleukin 2 (IL-2), interleukin 8 (IL-8), solute carrier family 1 member 5 (SLC1A5), peptide transporter 1 (Pep T1), and β-actin expression levels were determined using quantitative real-time PCR (ABI 7500, USA) with a reaction length of 20μl, including 10μl of 2 SYBR® Premix Ex Taq (TaKaRa, Dalian, China), 0.8μl for quantitative real-time PCR, specific primers were constructed based on sequences cloned and published in the O. mykiss gene bank (Table 2). The cycling conditions were 95°C for 30s followed by 35cycles of 95°C for 5s, 59°C for 10s, and 72°C for 30s. The housekeeping gene (β-actin) was chosen as a reference gene to normalize the results. 2−ΔΔCt was used to measure the expression values.
Intestine Histology
The intestinal samples were fixed in Bouin’s solution for 24h, rinsed multiple times with water to extract the fixative, dehydrated in ethanol, placed in xylene, embedded in wax for 2h at 60°C, sectioned to 6μm thick, with a microtome (Leica-RM2235), stained with hematoxylin–eosin, and sealed with neutral resin. Intestinal morphology was measured using Motic Images Plus 2.0 software after being photographed (Nikon, DS-Ri2). The muscular layer and villus height were determined by randomly selecting 10 villi per slide.
Intestine Bacterial DNA Extraction and Sequencing
Complete DNA was extracted from 0.2g of intestinal faeces using a DNA Extraction Kit (Beijing Tiangen Biochemical Technology Co. Ltd., China) according to the guidelines for total DNA extraction from intestinal microorganisms. DNA was isolated from the collected intestinal microorganisms using specific primers of the forward sequence 341F (5′-CCTAYGGGRBGCASCAG-3′) and reverse sequence 806R (5′-GGACTACNNGGGTATCTAAT-3′; Sun et al., 2020). The optimized conditions for amplification were as follows: one pre-denaturation cycle at 95°C for 5min, 27cycles of denaturation at 95°C for 30s, annealing at 55°C for 30s, elongation at 72°C for 45s, and a final extension at 72°C for 10min. The resulting amplicons were purified from 0.7% agarose gel, measured concentration using Qubit dsDNA broad-range assay kit (Life Technologies, United States), and equal concentration (20ngμl−1) of amplicons were pooled together, and sequenced on the Illumina MiSeq platform (300bp paired-end reads; Beijing Baimaike Technology Co. Ltd., China).
Intestinal Microbiota Bioinformatics Analysis
Using QIIME (Version 1.8.0, http://qiime.org/), low-quality reads with quality scores <20e were filtered out and raw reads were sorted according to their Barcode sequences (Caporaso et al., 2010). The reads were first cut in poor quality and then separated based on Barcode from the subsequent reads. Through binding the raw reads to the Barcode, the raw reads were then separated from the subsequent reads. To detect chimera, the Barcode and primer sequences were removed from the quality control, and the read sequences were removed using VSearch software (https://github.com/torognes/vsearch/; Li et al., 2015). The chimeric sequences were compared to the species annotation database, and the remaining chimeric sequences were deleted. The chimeric sequences were deleted from the database, leaving only the clean reads. The final correct reads were collected using the Uparse program (Version 7.0.1001, http://www.drive5.com/uparse/) to group all clean reads from all samples (Haas et al., 2011). The sequences were grouped into operative taxonomic units (OTUs) with 97% sequence similarity. The most frequent sequences were chosen as OTUs based on the algorithm (Edgar, 2010). According to the algorithm, the most commonly occurring sequences were chosen as symbolic OTUs. The Mothur system was used to species-annotate the OTUs sequences. Taxon resembling chloroplasts, mitochondria, unknowns, archaea, and eukaryotes were removed. The RDP Classifier Bayesian algorithm (http://sourceforge.net/projects/rdpclassifier/; Wang et al., 2007) was used to cluster the samples into OTUs, as was the GreenGene database1 for species annotation analysis, and OTUs for abundance. Raw reads were saved to the National Center for Biotechnology Information’s (NCBI) Sequence Read Archive (SRA) database (Accession Number: PRJNA714809). At each taxonomic rank, the species abundance, diversity index (Chao1 index, Shannon index and ACE index), and population structure were studied (Lozupone and Knight, 2005).
Calculations and Statistical Analysis
Weight gain rate (WGR, %)=100×(weight gain, g)/(initial weight, g); Survival rate (%)=100×(final number of fish)/(initial number of fish); Feed conversion ratio (FCR)=(dry dietary intake, g)/(weight gain, g). Specific growth rate (SGR, % per day)=100×[Ln (final weight, g)−Ln (initial weight, g)]/duration (days); Condition factor (CF)=100×(body weight, g)/(body length, cm)3; survival rate (SR, %)=100×(final fish number/initial fish number); Daily feed intake (DFI)=(feed consumed, g)/[(initial weight+final weight)/2, g]×(days, d).
All data in tables and figures are expressed as mean±standard deviation (SD). After normality and homogeneity checking, one-way variance analysis (ANOVA) and Duncan’s multiple range tests were used to examine the data. p values <0.05 were considered significantly different. The SPSS statistical package 23.0 (SPSS Inc., Chicago, IL, USA) was used for statistical analysis. The sigma plot software version 14.0 is used to draw column graphs and curves.
Results
Growth Performance
The effects of dietary GSH level on weight gain, feed conversion rate, and survival rate are shown in Table 3. With increasing dietary GSH level, weight gain rate and specific growth rate gradually increased, reaching the highest level in G400, showing a significant difference with G0 and G100 (p<0.05). The daily feed intake of G200 and G400 was significantly higher than that of the control group (p<0.05). The feed conversion rate decreased with increasing GSH, and G400 was the lowest, which was significantly different from G0 and G100 (p<0.05). The survival rate of G400 was the highest, and significantly higher than that of G0, G100, and G200 (p<0.05). There were no significant differences in condition factor among the groups (p>0.05). According to the relationship between weight gain and dietary GSH levels, the broken line model was selected as a good fit for the model, and the optimal dietary GSH level was 447.06mgkg−1 (Figure 1).
Figure 1. Relationship between weight gain and dietary GSH for O. mykiss as described by a broken line regression (n=3). The breakpoint in the broken line is 447.06mgkg−1.
Body Composition
Crude protein, crude lipid, ash, and moisture levels did not change significantly among the groups (p>0.05; Table 4).
Antioxidant Capacity
GSH content and GR activity did not vary considerably with increased dietary GSH supplementation (Table 5). Compared to the control group, SOD and CAT activities were significantly higher in G200, G400, and G800 (p<0.05). MDA content was significantly lower in G200 than in the other groups (p<0.05).
Table 5. Effects of dietary GSH on the antioxidant capacity of mid-intestine in O. mykiss (mean±SD, n=9).
Intestine Histology
Table 6 and Figure 2 show the mid-intestine morphology. No pathological differences in the intestines were found between the different groups. When the dietary GSH level was 0–200mgkg−1, the height of intestinal muscular layer was not significantly different (p>0.05); however, it was significantly higher in G400 (p<0.05). The intestinal villus height and width were also not significantly different between groups (p>0.05).
Figure 2. The light micrograph in the triploid O. mykiss fed GSH diets. VH, VW and MT represent villus height, villus width and muscular thickness.
Gene Expression
Figure 3 shows the mid-intestine gene expression. Compared to that of the control group, PepT1 and SLC1A5 gene expression in each GSH feeding group were significantly increased (p<0.05); G400 had the highest expression, and there was no significant difference between G100, G200, and G800 (p>0.05). TNF-α, IL-1β, IL-2, and IL-8 gene expressions were significantly reduced than that of the control group, respectively.
Figure 3. The gene expression of the intestine in each group (mean±SD, n=3). (A) SLC1A5, (B) pepT1, (C) TNF-α, (D) IL-1β, (E) IL-2, and (F) IL-8 [β-actin was chosen as a reference gene to normalize the results, different letters indicate significant differences between the different groups (p<0.05)].
Microbial Diversity
The 16S rDNA sequencing produced 1.18 million reads with a mean of 68,310 clean tags per sample. Chloroplast, mitochondria and eukaryotes were removed and reduced the amount of reads by 66.7% to a total of 1.0 million reads. Table 7 shows the Chao1, ACE, and Shannon indexes. Of all experimental treatments, fish fed a G400 diet had the highest values of Chao1, ACE, and Shannon index. On the phylum level, dietary GSH resulted in reduced Proteobacteria abundance and increased Bacteroidetes, Firmicutes, Acideobacteria, and Actinobacteria abundance (Figure 4). On the genus level, faeces were mainly composed of Cetobacterium, Nicotiana_otophora, Sphingomonas, Bacteroides, Pseudomonas, Bradyrhizobium, and Candidatus_Branchiomona (Figure 5). Figure 6 depicts a heatmap study of species abundance clustering at the phylum stage. Clustering results showed that the intestinal microbiota composition was classified into three classes (G0, G100; G400, G800; G200).
Table 7. Effects of dietary GSH on alpha diversity index of intestinal microbiota in O. mykiss (mean±SD, n=3).
Figure 4. Composition and relative abundance of bacterial community-based 16S rDNA sequences in phylum level (n=3).
Figure 5. Composition and relative abundance of bacterial community-based 16S rDNA sequences in genus level (n=3).
Figure 6. Heatmap analysis of the species abundance clustering in the top 20 on the phylum level (mean±SD, n=3).
Principle coordinate analysis (PCoA) illustrated beta diversity. Samples in PCoA using weighted and unweighted UniFrac distances were evenly distributed among the groups (Figure 7). The findings revealed that the intestinal microbiota composition of G0 and G100 are similar. However, there was no significant difference in beta diversity using ANOSIM (R=0.049, p=0.328).
Figure 7. Principle coordinate analysis (PCoA) based on weighted-unifrac (A) and unweighted-unifrac (B) analysis of bacterial profiles from intestines of O. mykiss (n=3).
Analysis with LefSe noted several indicator bacteria species associated with each diet (Table 8). In the faeces, dietary GSH 200–800mgkg−1 increased the abundance of Ilumatobacter, Peptoniphilus, Limnobacter, Mizugakiibacter, Chelatococcus, Stella, Filimonas, and Streptosporangium, while the control diet had increased abundance of Ferrovibrio, Buchnera, Chitinophaga, Stenotrophobacter, Solimonadaceae, Polycyclovorans, Rhodococcus, Ramlibacter, Azohydromonas, and Arcobacter.
Table 8. Linear discriminant analysis effect size (Lefse) of indicator bacteria species that were significantly (p<0.05) associated with each group.
Microbial Function
The microbial functions of the intestine were predicted with PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States; Figure 8). The intestine microbiota was enriched with functions related to carbohydrate metabolism, amino acid metabolism, energy metabolism, the metabolism of cofactors and vitamins, nucleotide metabolism, membrane transport, translation, lipid metabolism, and xenobiotics biodegradation and metabolism. However, the abundant among the groups had little change, and it is speculated that the microbial functions of the intestine were similar.
Discussion
Growth Performance
The results showed that dietary GSH can significantly improve the weight gain rate and survival rate of triploid O. mykiss, and reduce the feed conversion rate, which is similar to results obtained for grass carp (Ming et al., 2015) and Atlantic salmon (Ma et al., 2019). GSH can enhance the growth performance of European bass larvae by improving amylase and protease activity (Zambonino-Infante et al., 1997). In a study of Japanese flounder (Paralichthys olivaceus), the deamination product of cysteine and hydrolysate of GSH was found to be a component of coenzyme A (Wang et al., 2011). This can also enhance grass carp growth by destroying the growth inhibitory molecules, and promoting the secretion of growth hormone (Xiao and Lin, 2003). It was found that GSH can promote the secretion of pituitary growth hormone and liver insulin-like growth factor-1 level, to promote the synthesis of protein and improve the nutrient utilization rate (Ming et al., 2015). Furthermore, GSH can protect intestinal mucosa from the damage of toxins and peroxides, protect the digestive system, and help the intestinal absorption of nutrients in animals (Aw et al., 1992), which is in line with reports of improving growth and feed conversion rate.
In the present study, the optimal addition of GSH in the feed of triploid O. mykiss was 447.06mgkg−1, which is close to that of grass carp (381mgkg−1; Ming et al., 2015), but much higher than that of Japanese flounder (368.92mg·kg−1; Wang et al., 2011) and Nile tilapia (355.13mgkg−1; Zhou et al., 2013). The differences in these results may be related to the culture environment, species and size, feed formula, feeding strategy, and raw material processing technology. Adding an appropriate amount of GSH to the feed can improve the growth performance, but excessive GSH will result in negative effects. The weight and protein efficiency of yellow catfish (Pelteobagrus fulvidraco) first increased and then decreased, and reached a maximum value at 300mgkg−1, and the high GSH levels possibly had a toxic effect on the fish (Zhou et al., 2017). A low GSH concentration can inhibit lipid peroxidation, whereas a high GSH concentration can promote mitochondrial lipid peroxidation, which is related to the mutual transformation of reduced glutathione and oxidized glutathione in the body (Gao et al., 1998). Moreover, GSH can be synthesized with many compounds, such as aromatic epoxides, halogenated hydrocarbons, to produce toxic metabolites, which cause DNA damage and have adverse effects on the body (Kleider et al., 2016).
Body Composition
The effects of GSH on fish body composition varied with different species. There was no significantly different effect of GSH on the dry matter, crude lipid, and ash content of grass carp (Ming et al., 2015). Similarly, there were no significant differences among the groups for body composition in this study. However, crude lipid and protein in yellow catfish were higher than that of the control group (Zhou et al., 2017). This may be related to GSH participating in the transport of amino acids and small peptides, which promotes amino acid absorption (Meister and Anderson, 1983).
Histology
Increased intestinal villus height can enhance digestion and absorption by improving the interaction between the intestine and nutrients (Long et al., 2018). The results in this study were similar to previous reports of Atlantic salmon (Refstie et al., 2000) and meager (Argyrosomus regius; Ribeiro et al., 2015). The small peptide can improve the intestinal development of fish, and increase the height of small intestinal villi (Murashita et al., 2015). It has been reported that adding GSH to the feed of Pacific white shrimp improved intestinal development, promoted intestinal villi growth, and thickened the muscle layer (Wang et al., 2018). It was found that when fish were fed dietary GSH at 400mgkg−1, the myometrium thickness and intestinal villus height increased significantly. The results showed that GSH could enhance intestinal tract development by affecting the tissue structure of the small intestine, improve the utilization rate of feed nutrition, and enhance the growth and development of triploid O. mykiss.
Antioxidant Capacity
In the present study, the SOD, CAT, GSH, and glutathione reductase (GR) activities in the mid-intestine when fish fed GSH were higher than those of the control group. Similarly, the SOD activity also increased as dietary GSH levels increased (Zhou et al., 2012). It has been speculated that exogenous GSH can reduce the oxygen-free radical damage to the body and antioxidant stress of O. mykiss (Pena-Llopis et al., 2001). However, the SOD activity of Japanese flounder was not affected by the amount of GSH added to the feed (Wang et al., 2011). It is possible that the effect of GSH on SOD activity is different for different species. With the increase of GSH, CAT activity in tilapia increased first and then decreased; therefore, CAT could possibly eliminate hydrogen peroxide in the body and protect cells from its toxicity (Zhou et al., 2012). The change in CAT activity was due to the activation of the antioxidant enzyme system in the body. In this study, the mid-intestine CAT activity in triploid O. mykiss was higher than that of the control group. Therefore, adding appropriate GSH could regulate body antioxidant status, decrease the concentration of hydrogen peroxide in the cells, and reduce the degree of cell damage. In this experiment, the mid-intestine GSH and GR activities in triploid O. mykiss are higher than those in the control group. It is speculated that exogenous GSH could activate a reaction and metabolism related to GSH and increase the demand for GSH. Therefore, the GR activity that can convert oxidized GSH into GSH (Ming et al., 2015) also increased accordingly. Furthermore, intestinal MDA content when fish were fed GSH was lower than that of the control group, which is associated with previous study findings (Wang et al., 2011). It is possible that an appropriate amount of GSH could reduce the toxic effect of lipid peroxidation in cells, reducing the degree of cell damage (Birnie-Gauvin et al., 2017).
Gene Expression
When protein ingestion increases, the PepT1 expression in the brush border of the small intestine increases, transport efficiency of small peptides increases, and the use of protein, tripeptides, and free amino acids will be improved through enhanced amino acid uptake efficiency (Ostaszewska et al., 2010). SLC1A5, an amino acid transporter, transports amino acids using the concentration gradient of Na+ inside and outside the cell membrane and is the most important transporter for cells to absorb exogenous glutamine. Compared to that of the control group, the PepT1 and SLC1A5 gene expression in the intestinal oligopeptide transporter by dietary GSH supplementation increased significantly, suggesting that dietary glutamine utilization was improved, which enhances intestinal development.
Proinflammatory factors are necessary to signal cytokines to initiate and regulate this reaction, which mainly includes tumor necrosis factor-α, interleukin-1β, interleukin-2, and interleukin-8. Interleukin-1β can promote the proliferation and activation of immune cells, such as thymocytes and T cells, and also promote the synthesis and secretion of immune proteins by B cells, and mediate the inflammatory response (Sims and Smith, 2010). Interleukin-2 is primarily generated by activated T cells, which can promote the growth, proliferation, and differentiation of lymphocytes, enhance NK cell function and play an important role in immune response and antiviral infection (Koreth et al., 2011). Interleukin-8 regulates immunity, promoting cell mitosis, and stimulating capillary formation (Hoffmann et al., 2002). Tumor necrosis factor-α is a member of the cytokine family, which can induce apoptosis of tumor cells and coordinate non-specific immune response (Paul, 1985), which is produced in the early stage of inflammatory response. When the level of intestinal pro-inflammatory factors increases, it promotes the production of the inflammatory response, causing damage to intestinal mucosa tissue cells, destroying the barrier function of the intestinal epithelium, increasing intestinal permeability, making pathogens and endotoxins enter the blood circulation, causing functional damage to more organs, resulting in a systemic inflammatory response of the body (Costantini et al., 2010). It has been found that GSH can reduce intestinal injury, the expression of intestinal inflammatory factors, incidence of pathogen and endotoxin translocation in the intestinal tract of rats with acute necrotizing pancreatitis, and play a role in protecting the intestinal mucosa (Aw et al., 1992). In addition to interleukin-8, the expression of proinflammatory factors in the intestinal tract in the treatment group was significantly lower than that of the other groups. The mechanism of action may be related to the NF-κB/MLCK pathway inhibition by GSH’s nitroso derivatives, to protect the structure and function of closely connected intestinal epithelial cells (Koeberle et al., 2020).
Intestinal Microbiota
Feed composition is a key aspect influencing fish intestinal microbiota (Ringø et al., 2006). In this study, when fish were fed GSH, the Chao1 index, ACE index, and Shannon index increased, which indicated that the abundance and evenness of the intestinal microbiota of triploid O. mykiss increased. The abundance of the intestinal microbiota reached a significant level when fish were fed 200–800mgkg−1 GSH. Results have shown that several species with low microbiotal abundance in the intestines could benefit from dietary GSH supplementation. The change in alpha diversity might be related to the antibacterial properties of GSH, and the mechanism needs further study. In the present study, the predominant intestinal microbiota were Proteobacteria, Firmicutes, and Bacteroidetes, which were similar to those of O. mykiss (Huyben et al., 2017; Lyons et al., 2017; Huyben et al., 2021). Some studies have shown that Firmicutes and Bacteroides in the intestine are related to fat deposition, and can ferment more short-chain fatty acids, and promote fat deposition when the proportion of Bacteroides increases (Turnbaugh et al., 2006). This study showed that the abundance of Firmicutes/Bacteroides in each group increased, and the proportion of Actinobacteria with the addition of GSH increased compared with that of the control group. Proteobacteria includes many pathogenic bacteria, such as Escherichia coli, Vibrio cholerae, and Salmonella enterica (Holben et al., 2002), which are usually present in the intestines. In this study, although there were no significant differences in beta diversity among the groups, dietary GSH 200–800mgkg−1 increased the abundance of Ilumatobacter, Peptoniphilus, Limnobacter, etc. Moreover, the control diet had an increased abundance of Arcobacter, which is responsible for causing diseases in fish. The results in the present study indicated that 200–800mgkg−1 dietary GSH may decrease the probability of triploid O. mykiss being infected by pathogenic bacteria. This may be related to the role of GSH in scavenging intracellular peroxides and free radicals, protecting cells from oxidative damage. This would maintain the intestinal mucosal barrier, promoting immune cell proliferation, and inhibiting the expression of pro-inflammatory factors, resulting in increased abundance and homogeneity of microbiota in the intestines of triploid O. mykiss. However, the microbial functions of the intestine were similar through the functional analysis. Therefore, the dynamics of the microbiota structure and mechanism of action between the altered microbiota and intestine function are still not well understood and need more research.
Conclusion
In conclusion, dietary GSH is beneficial for triploid O. mykiss and may be related to the ability of GSH to scavenge peroxides and free radicals in cells, protect cells from oxidative damage, maintain the intestinal mucosal barrier, promote immune cell proliferation, inhibit the expression of the pro-inflammatory factors, and increase the abundance and evenness of intestinal microbiota. The mechanism of interaction is still unclear, and further research is needed. In addition, more studies are required to elucidate the exact regulatory mechanisms involved by which GSH causes these coordinated responses and explore products (e.g., yeast, wheat germ, and animal liver which contains approximately 1–100mgg−1 GSH) or use feed ingredients with high glutathione content appropriately to generate substantial economic returns for the aquaculture industry.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/genbank/, PRJNA714809.
Ethics Statement
The animal study was reviewed and approved by Chinese Animal Health Protection Law and the Scientific Laboratory Animal Permit Approval [Ethical Approval No. SCXK(YU)2005-0001]. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
CW, HL and YY designed the study. BS, YL, and SH carried out the rearing work. CW, BS, HJ, and ZL tested the samples and analyzed the results. CW wrote the manuscript with contributions from the other authors. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Natural Science Funds of Heilongjiang (YQ2019C036), China Agriculture Research System of MOF and MARA (CARS-46), China Scholarship Council (Grant No.202003260012), the Science and Technology Project of Guizhou Province (20162502 and 20162511), the Guizhou Science and Technology Plan Project (QKHZC20172532), the Guizhou Technology Innovation Team Project (QKHRCTD20154016), and the Beijing Sturgeon & Trout Innovation Team (BAIC08-2018).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors thank the participants who gave their time to the trial.
Footnotes
References
AOAC, 2012. Official Methods for Analysis, 19th edn. Association of Official Analytical Chemists, Washington, DC.
Aw, T. Y., Williams, M. W., and Gray, L. (1992). Absorption and lymphatic transport of peroxidized lipids by rat small intestine in vivo: role of mucosal GSH. Am. J. Physiol. Gastrointest. Liver Physiol. 262, G99–G106. doi: 10.1152/ajpgi.1992.262.1.G99
Birnie-Gauvin, K., Costantini, D., Cooke, S. J., and Willmore, W. G. (2017). A comparative and evolutionary approach to oxidative stress in fish: a review. Fish Fish. 18, 928–942. doi: 10.1111/faf.12215
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Costantini, T. W., Bansal, V., Krzyzaniak, M., Putnam, J. G., Peterson, C. Y., Loomis, W. H., et al. (2010). Vagal nerve stimulation protects against burn-induced intestinal injury through activation of enteric glia cells. Am. J. Physiol. Gastrointest. Liver Physiol. 299, G1308–G1318. doi: 10.1152/ajpgi.00156.2010
Deng, J., Mai, K., Chen, L., Mi, H., and Zhang, L. (2015). Effects of replacing soybean meal with rubber seed meal on growth, antioxidant capacity, non-specific immune response, and resistance to Aeromonas hydrophila in tilapia (Oreochromis niloticus × O. aureus). Fish Shellfish Immunol. 44, 436–444. doi: 10.1016/j.fsi.2015.03.018
Doyotte, A., Cossu, C., Jacquin, M., Babut, M., and Vasseur, P. (1997). Antioxidant enzymes, glutathione and lipid peroxidation as relevant biomarkers of experimental or field exposure in the gills and the digestive gland of the freshwater bivalve Unio tumidus. Aquat. Toxicol. 39, 93–110. doi: 10.1016/S0166-445X(97)00024-6
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Fjelldal, P. A., Hansen, T. J., Lock, E. J., Wargelius, A., Fraser, T., Sambraus, F., et al. (2016). Increased dietary phosphorous prevents vertebral deformities in triploid Atlantic salmon (Salmo salar L.). Aquac. Nutr. 22, 72–90. doi: 10.1111/anu.12238
Gao, S., Liu, D., and Luo, G. (1998). ESR study on the mechanism of anti-lipid peroxidation of antioxidants. Chin. J. Magn. Reson. 15, 139–143.
Good, C., and Davidson, J. (2016). A review of factors influencing maturation of Atlantic salmon, Salmo salar, with focus on water recirculation aquaculture system environments. J. World Aquacult. Soc. 47, 605–632. doi: 10.1111/jwas.12342
Haas, B. J., Gevers, D., Earl, A. M., Feldgarden, M., Ward, D. V., Giannoukos, G., et al. (2011). Chimeric 16S rRNA sequence formation and detection in sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. doi: 10.1101/gr.112730.110
Hoffmann, E., Dittrich Breiholz, O., Holtmann, H., and Kracht, M. (2002). Multiple control of interleukin-8 gene expression. J. Leukoc. Biol. 47, 605–632. doi: 10.1111/jwas.12342
Holben, W. E., Williams, P., Saarinen, M., Särkilahti, L. K., and Apajalahti, J. H. (2002). Phylogenetic analysis of intestinal microflora indicates a novel mycoplasma phylotype in farmed and wild salmon. Microb. Ecol. 44, 175–185. doi: 10.1007/s00248-002-1011-6
Huyben, D., Chiasson, M., Lumsden, J. S., Pham, P. H., and Chowdhury, M. A. K. (2021). Dietary microencapsulated blend of organic acids and plant essential oils affects intestinal morphology and microbiome of rainbow trout (Oncorhynchus mykiss). Microorganisms 9:2063. doi: 10.3390/microorganisms9102063
Huyben, D., Nyman, A., Vidaković, A., Passoth, V., Moccia, R., Kiessling, A., et al. (2017). Effects of dietary inclusion of the yeasts Saccharomyces cerevisiae and Wickerhamomyces anomalus on gut microbiota of rainbow trout. Aquaculture 473, 528–537. doi: 10.1016/j.aquaculture.2017.03.024
Kelly, C., and Salinas, I. (2017). Under pressure: interactions between commensal microbiota and the teleost immune system. Front. Immunol. 8:559. doi: 10.3389/fimmu.2017.00559
Kleider, C., Esch, H. L., and Lehmann, L. (2016). Reactivity against DNA bases despite previous conjugation to glutathione–A new mechanism of toxicity for α, β-unsaturated carbonyls? Toxicol. Lett. 258:S254. doi: 10.1016/j.toxlet.2016.06.1899
Koeberle, S. C., Gollowitzer, A., Laoukili, J., Kranenburg, O., Werz, O., Koeberle, A., et al. (2020). Distinct and overlapping functions of glutathione peroxidases 1 and 2 in limiting NF-κB-driven inflammation through redox-active mechanisms. Redox Biol. 28:101388. doi: 10.1016/j.redox.2019.101388
Koreth, J., Matsuoka, K., Kim, H. T., McDonough, S. M., Bindra, B., Alyea, E. P. III, et al. (2011). Interleukin-2 and regulatory T cells in graft-versus-host disease. New Engl. J. Med. 365, 2055–2066. doi: 10.1056/NEJMoa1108188
Kosower, N. S., and Kosower, E. M. (1978). The glutathione status of cells. Int. Rev. Cytol. 54, 109–160. doi: 10.1016/S0074-7696(08)60166-7
Li, S., Zhang, C., Gu, Y., Chen, L., Ou, S., Wang, Y., et al. (2015). Lean rats gained more body weight than obese ones from a high-fibre diet. Brit. J. Nutr. 114, 1188–1194. doi: 10.1017/S0007114515002858
Liu, B., Ge, X., He, Y., Xie, J., Xu, P., He, Y., et al. (2010). Effects of anthraquinones extracted from Rheum officinale Bail on the growth, non-specific immune response of Macrobrachium rosenbergii. Aquaculture 310, 13–19. doi: 10.1016/j.aquaculture.2010.09.020
Long, M., Yang, S., Li, P., Song, X., Pan, J., He, J., et al. (2018). Combined use of C. butyricum Sx-01 and L. salivarius C-1-3 improves intestinal health and reduces the amount of lipids in serum via modulation of gut microbiota in mice. Nutrients 10:810. doi: 10.3390/nu10070810
Lozupone, C., and Knight, R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005
Lyons, P. P., Turnbull, J. F., Dawson, K. A., and Crumlish, M. (2017). Exploring the microbial diversity of the distal intestinal lumen and mucosa of farmed rainbow trout Oncorhynchus mykiss (Walbaum) using next generation sequencing (NGS). Aquac. Res. 48, 77–91. doi: 10.1111/are.12863
Ma, J., Zhang, J., Sun, G., Lou, Y., and Li, Y. (2019). Effects of dietary reduced glutathione on the growth and antioxidant capacity of juvenile Atlantic salmon (Salmo salar). Aquac. Nutr. 25, 1028–1035. doi: 10.1111/anu.12920
Meister, A., and Anderson, M. E. (1983). Glutathione. Annu. Rev. Biochem. 52, 711–760. doi: 10.1146/annurev.bi.52.070183.003431
Ming, J., Ye, J., Zhang, Y., Xu, P., and Xie, J. (2015). Effects of dietary reduced glutathione on growth performance, non-specific immunity, antioxidant capacity and expression levels of IGF-I and HSP70 mRNA of grass carp (Ctenopharyngodon idella). Aquaculture 438, 39–46. doi: 10.1016/j.aquaculture.2014.12.038
Murashita, K., Fukada, H., Takahashi, N., Hosomi, N., Matsunari, H., Furuita, H., et al. (2015). Effect of feed ingredients on digestive enzyme secretion in fish. Bull. Fish. Res. Agency 40, 69–74.
Ostaszewska, T., Dabrowski, K., Kamaszewski, M., Grochowski, P., Verri, T., Rzepkowska, M., et al. (2010). The effect of plant protein-based diet supplemented with dipeptide or free amino acids on digestive tract morphology and PepT1 and PepT2 expressions in common carp (Cyprinus carpio L.). Comp. Biochem. Physiol. A 157, 158–169. doi: 10.1016/j.cbpa.2010.06.162
Paul, W. E. (1985). Fundamental immunology. Fundam. Immun. 14, 549–550. doi: 10.3109/08820138509022683
Pena-Llopis, S., Pena, J. B., Sancho, E., Fernandez-Vega, C., and Ferrando, M. D. (2001). Glutathione-dependent resistance of the European eel Anguilla anguilla to the herbicide molinate. Chemosphere 45, 671–681. doi: 10.1016/S0045-6535(00)00500-2
Ponsoda, X., Bort, R., Jover, R., Gómez-Lechón, M. J., and Castell, J. V. (1999). Increased toxicity of cocaine on human hepatocytes induced by ethanol: role of GSH. Biochem. Pharmacol. 58, 1579–1585. doi: 10.1016/S0006-2952(99)00249-X
Raghunathan, V. K., Ellis, E. M., Tettey, J. N., and Grant, M. H. (2007). Involvement of reduced glutathione and glutathione reductase in the chronic toxicity of hexavalent chromium to monocytes in vitro. Toxicology 231, 105–106. doi: 10.1016/j.tox.2006.11.030
Refstie, S., Korsøen, Ø. J., Storebakken, T., Baeverfjord, G., Lein, I., and Roem, A. J. (2000). Differing nutritional responses to dietary soybean meal in rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar). Aquaculture 190, 49–63. doi: 10.1016/S0044-8486(00)00382-3
Ren, J., Wu, P., Trampe, B., Tian, X., Lübberstedt, T., and Chen, S. (2017). Novel technologies in doubled haploid line development. Plant Biotechnol. J. 15, 1361–1370. doi: 10.1111/pbi.12805
Ribeiro, L., Moura, J., Santos, M., Colen, R., Rodrigues, V., Bandarra, N., et al. (2015). Effect of vegetable based diets on growth, intestinal morphology, activity of intestinal enzymes and haematological stress indicators in meagre (Argyrosomus regius). Aquaculture 447, 116–128. doi: 10.1016/j.aquaculture.2014.12.017
Ringø, E., Sperstad, S., Myklebust, R., Refstie, S., and Krogdahl, Å. (2006). Characterisation of the microbiota associated with intestine of Atlantic cod (Gadus morhua L.): the effect of fish meal, standard soybean meal and a bioprocessed soybean meal. Aquaculture 261, 829–841. doi: 10.1016/j.aquaculture.2006.06.030
Rombout, J. H., Abelli, L., Picchietti, S., Scapigliati, G., and Kiron, V. (2011). Teleost intestinal immunology. Fish Shellfish Immunol. 31, 616–626. doi: 10.1016/j.fsi.2010.09.001
Sims, J. E., and Smith, D. E. (2010). The IL-1 family: regulators of immunity. Nat. Rev. Immunol. 10, 89–102. doi: 10.1038/nri2691
Sun, X., Li, M., Xia, L., Fang, Z., Yu, S., Gao, J., et al. (2020). Alteration of salivary microbiome in periodontitis with or without type-2 diabetes mellitus and metformin treatment. Sci. Rep. 10, 15363–15314. doi: 10.1038/s41598-020-72035-1
Thorgaard, G. H. (1983). 8 chromosome set manipulation and sex control in fish. Fish Physiol. 9, 405–434. doi: 10.1016/S1546-5098(08)60308-8
Turnbaugh, P. J., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., and Gordon, J. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031. doi: 10.1038/nature05414
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, X., Xu, W., Zhou, H., Zhang, Y., Gao, W., Zhang, W., et al. (2018). Reduced glutathione supplementation in practical diet improves the growth, anti-oxidative capacity, disease resistance and gut morphology of shrimp Litopenaeus vannamei. Fish Shellfish Immunol. 73, 152–157. doi: 10.1016/j.fsi.2017.11.043
Wang, F., Zhang, W., Mai, K., Xu, W., Wang, X., and Liufu, Z. (2011). Effects of dietary reduced glutathione on growth and antioxidative capability of Japanese flounder, Paralichthys olivaceus. Period. Ocean Univ. China 4, 55–60. doi: 10.3969/j.issn.1672-5174.2011.04.008
Will, Y. (1999). Overview of glutathione function and metabolism. Curr. Protoc. Toxicol. 1, 1–6. doi: 10.1002/0471140856.tx0601s00
Xia, Z., and Wu, S. (2018). Effects of glutathione on the survival, growth performance and non-specific immunity of white shrimps (Litopenaeus vannamei). Fish Shellfish Immunol. 73, 141–144. doi: 10.1016/j.fsi.2017.12.015
Xiao, D., and Lin, H. (2003). Effects of cysteamine—a somatostatin-inhibiting agent—on serum growth hormone levels and growth in juvenile grass carp (Ctenopharyngodon idellus). Comp. Biochem. Physiol. A 134, 93–99. doi: 10.1016/S1095-6433(02)00191-5
Zambonino-Infante, J. L., Cahu, C. L., and Peres, A. (1997). Partial substitution of di-and tripeptides for native proteins in sea bass diet improves Dicentrarchus labrax larval development. J. Nutr. 127, 608–614. doi: 10.1093/jn/127.4.608
Zhou, T., Cao, J., Huang, Y., Wang, G., Zhao, H., Sun, Z., et al. (2012). Effects of dietary glutathione on the mRNA expression of growth and antioxidant related genes in liver of GIFT Oreochromis niloticus. Guang. Agric. Sci. 14, 152–155. doi: 10.3969/j.issn.1004-874X.2012.14.045
Zhou, W., Liang, Q., Xu, Y., Gurr, G. M., Bao, Y., Zhou, X., et al. (2013). Genomic insights into the glutathione S-transferase gene family of two rice planthoppers, Nilaparvata lugens (Stål) and Sogatella furcifera (Horváth)(Hemiptera: Delphacidae). PLoS One 8:e56604. doi: 10.1371/journal.pone.0056604
Zhou, Y., Sun, Y., Cao, J., Huang, Y., Wang, G., Chen, B., et al. (2017). Effects of dietary glutathione on growth performance, body composition, serum biochemical indices and anti-ammonia-nitrogen stress ability of juvenile yellow catfish (Pelteobagrus fulvidraco). Chin. J. Ani. Nutr. 29, 4189–4197. doi: 10.3969/j.issn.1006-267x.2017.11.043
Keywords: triploid O. mykiss, glutathione, intestinal health, growth, microbiota
Citation: Wang C, Su B, Lu S, Han S, Jiang H, Li Z, Liu Y, Liu H and Yang Y (2021) Effects of Glutathione on Growth, Intestinal Antioxidant Capacity, Histology, Gene Expression, and Microbiota of Juvenile Triploid Oncorhynchus mykiss. Front. Physiol. 12:784852. doi: 10.3389/fphys.2021.784852
Edited by:
Serhat Turkmen, University of Alabama at Birmingham, United StatesReviewed by:
Beth Cleveland, Cool and Cold Water Aquaculture Research (USDA-ARS), United StatesDavid Huyben, University of Guelph, Canada
Copyright © 2021 Wang, Su, Lu, Han, Jiang, Li, Liu, Liu and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chang’an Wang, gordoncase@126.com; Hongbai Liu, liuhongbai@hrfri.ac.cn; Yuhong Yang, yuhongyang68@163.com