• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Debaryomyces hansenii supplementation in low fish meal diets promotes growth,modulates microbiota and enhances intestinal condition in juvenile marine fish

    2024-03-14 13:19:52IgnasiSanahujaAlbertoRuizJoanaFirminoFelipeReyespezJuanOrtizDelgadoEvaVallejosVidalLluisTortDarielTovarRamrezIsabelCerezoMiguelMoriigoCarmenSarasqueteandEnricGisbert

    Ignasi Sanahuja,Alberto Ruiz,Joana P.Firmino,Felipe E.Reyes-López,Juan B.Ortiz-Delgado,Eva Vallejos-Vidal,Lluis Tort,Dariel Tovar-Ramírez,Isabel M.Cerezo,Miguel A.Mori?igo,Carmen Sarasquete^ and Enric Gisbert*

    Abstract Background The development of a sustainable business model with social acceptance,makes necessary to develop new strategies to guarantee the growth,health,and well-being of farmed animals.Debaryomyces hansenii is a yeast species that can be used as a probiotic in aquaculture due to its capacity to i) promote cell proliferation and differentiation,ii) have immunostimulatory effects,iii) modulate gut microbiota,and/or iv) enhance the digestive function.To provide inside into the effects of D.hansenii on juveniles of gilthead seabream (Sparus aurata) condition,we integrated the evaluation of the main key performance indicators coupled with the integrative analysis of the intestine condition,through histological and microbiota state,and its transcriptomic profiling.Results After 70 days of a nutritional trial in which a diet with low levels of fishmeal (7%) was supplemented with 1.1% of D.hansenii (17.2 × 105 CFU),an increase of ca. 12% in somatic growth was observed together with an improvement in feed conversion in fish fed a yeast-supplemented diet.In terms of intestinal condition,this probiotic modulated gut microbiota without affecting the intestine cell organization,whereas an increase in the staining intensity of mucins rich in carboxylated and weakly sulphated glycoconjugates coupled with changes in the affinity for certain lectins were noted in goblet cells.Changes in microbiota were characterized by the reduction in abundance of several groups of Proteobacteria,especially those characterized as opportunistic groups.The microarrays-based transcriptomic analysis found 232 differential expressed genes in the anterior-mid intestine of S.aurata,that were mostly related to metabolic,antioxidant,immune,and symbiotic processes.Conclusions Dietary administration of D.hansenii enhanced somatic growth and improved feed efficiency parameters,results that were coupled to an improvement of intestinal condition as histochemical and transcriptomic tools indicated.This probiotic yeast stimulated host-microbiota interactions without altering the intestinal cell organization nor generating dysbiosis,which demonstrated its safety as a feed additive.At the transcriptomic level,D.hansenii promoted metabolic pathways,mainly protein-related,sphingolipid,and thymidylate pathways,in addition to enhance antioxidant-related intestinal mechanisms,and to regulate sentinel immune processes,potentiating the defensive capacity meanwhile maintaining the homeostatic status of the intestine.

    Keywords Debaryomyces hansenii,Intestine condition,Low fish meal diet,Microbiota,Transcriptomics,Yeast probiotic

    Background

    Yeasts are single-celled members of the fungi kingdom,widely distributed in nature,capable of living in extreme environments.Its biological relevance is determined by 1,500 species described so far,with a high impact on both plant and animal life [1],being also a pillar in the human food industry.Used for millennia due to their fermentative activity to produce highly demanded foods (bread,cheese,wine,beer,among others),for some years now yeasts have also been in the crosshair of the pharmaceutical industry[2].Some authors have suggested their use as transport vehicles for bioactive compounds [3],drug delivery [4],or what is more widely considered,as natural probiotics for humans,livestock,and aquaculture species [5,6].Some of them have been isolated from the intestine of several species as part of their commensal microbiota,including fish [7].In this sense,the reasons why some yeast species are highly interesting for their use as probiotics is because they may promote somatic growth,modulate the intestinal microbiota,as well as enhance host health and condition due to their contents in immunomodators like β-glucans among other bioactive compounds [8-10].

    Several studies have recommended the inclusion of yeasts in livestock and aquaculture feeds,as an additive,due to their beneficial effects in key performance indicators related to growth and feed efficiency,and in terms of antioxidant performance,and immunity enhancement,among which is the resistance to infections [11-13].On the other hand,yeasts and their metabolic products from food industry by-products could also be used to produce high quality alternative protein sources within a sustainable and circular economy approach [14-16].It is known that yeasts are a rich source of nutrients and proteinrelated compounds with a wide range of bioactive activities.In this sense,yeasts cells contain different types of wall-related compounds,like glucans and polyamines(biologically active amines),which are essential for the maintenance of life [17].Specifically,amines have been observed to be relevant to the normal cell development and are involved in several cellular processes,along with systemic benefits for the animals [18-20].

    Among all yeast species,Saccharomyces cerevisiaemight be the most studied yeast as a probiotic or alternative protein source.However,another yeast species that has received a lot of attention for its probiotic effect isDebaryomyces hansenii[21].Contrary to other species,D.hanseniibelongs to the normal microbiota of wild and cultivated animals [22],and it has been reported that it also promotes growth and development in various species at different stages,results that have been attributed to its glucans and polyamine content [21].This probiotic has been tested at early life stages in aquatic species with promising results;for instance,larvae ofSeriola rivoliana[23],Danio rerio[24],Dicentrarchus labrax[25]andPenaeus monodon[26] showed improved survival and development when fed a diet supplemented withD.hansenii.Similarly,other studies focused on juveniles and adults showed that a diet supplemented with this yeast had a positive effect on the host in terms of its physiological and immune responses,showing its capacity to counteract possible stress impacts as well as to enhance disease resistance [17,27,28],along with an improvement of the intestinal health and function [29].However,little is known about the effects on fish performance,and especially about the modes of action on the intestinal mucosa,beyond a purely immunological aspect.

    The aim of this study was to evaluate the effect of a feed supplemented withD.hanseniion gilthead seabream (Sparus aurata) on common key performance indicators (KPI)associated to somatic growth and feed efficiency coupled with the analysis of the transcriptomic profile and the histological organization of the anterior-mid intestine,and its autochthonous microbiota in order to provide inside into the effects of this yeast on host condition.The convergence offered in this study by the integration of these analytics in the intestinal mucosa in one of the most highly valued species in the Mediterranean [30],is highly relevant for the introduction and management of sustainable feeding strategies for the growing aquaculture industry.

    Material and methods

    Diet composition

    To evaluate the effects ofD.hanseniias a probiotic in aquafeeds,two diets were formulated with similar levels of protein content (48.4% crude protein),lipid content (17.2% crude fat),and energy levels (21.7 MJ/kg gross energy) (Table 1).Both diets,named Control and Yeast diets,were formulated with low levels of fishmeal (7%) where fishmeal was replaced by a blend of plant-based meals (wheat,corn,pea and soy) (Table 1).Diets only differed in the inclusion ofD.hansenii(CBS8339) at an inclusion rate of 1.1% (17.2 × 105CFU).D.hanseniiwas provided by CIBNOR (La Paz,Mexico)and grown as described in Tovar-Ramírez et al.[25].

    Table 1 Ingredients and proximal composition of the experimental diets

    Both diets (2 mm pellet size) were manufactured by SPAROS Lda (Olh?o,Portugal) in a temperature-controlled low-shear extruder as detailed in Gisbert et al.[31].Samples of each diet were taken for proximate composition,and feeds were stored at 4 °C along the experimental feeding trial to prevent their oxidation.

    Experimental model and design

    A total of 500S.auratajuveniles (body weight,BW=14.5 ± 1.2 g;mean ± SD) were obtained from a commercial farm (PISCIMAR S.L.;AVRAMAR,Burriana,Spain).Fish were transported to IRTA-La Ràpita research facilities (La Ràpita,Spain),and acclimated for two weeks in two tanks of 2,000 L connected to a IRTAmar?recirculating system under natural photoperiod and constant temperature (23.0 ± 1.0 °C).During acclimatation,fish were fed with a commercial diet containing 48.5% protein,18% lipid,and 18.5 MJ/kg digestible energy (Optibream,Skretting).

    Before the trial,200 fish (25 fish per tank;4 tank replicates per diet) were gently anesthetized (50 mg/L tricaine methanesulfonate,MS-222,Sigma-Aldrich,Madrid,Spain)and individually measured in BW and standard length(SL) to the nearest 0.1 g and 1 mm,respectively.The trial lasted for 70 d,during which fish were fed at 3.5% of the stocked biomass by means of automatic feeders (ARVOTEC T Drum 2000;Arvotec,Finland).The daily feed ratio was evenly distributed in 4 meals at 08:00,11:00,13:00 and 16:00 h (at each meal,the corresponding feed ration was distributed during 1 h).One hour after each meal,uneaten pellets were collected,dried overnight (120 °C),and weighted (g) for calculating daily feed ingesta values,while feed ration was adjusted to guarantee 10%-15% of uneaten pellets;thus,confirming that fish were fed at satiation.This trial was run in a recirculating system IRTAmar?under natural photoperiod conditions (14 h light/8 h dark),constant temperature (23.3 ± 1.3 °C),salinity (35-36 ppt),dissolved oxygen (5.7 ± 0.2 mg/L) (OXI330,Crison Instruments,Spain),and pH (8.2 ± 0.1) (pHmeter 507,Crison Instruments).Ammonia (0.13 ± 0.1 mg) and nitrite levels (0.18 ± 0.1 mg) (HACH DR 900 Colorimeter,Hach Company,Spain) were weekly controlled.

    Growth and feed performance indicators were calculated using the following formulae:

    - are the final and initial mean BW of fish.

    - Survival rate(SR,%)=100 × (Final number of fish ∕Initial number of fish).

    - Feed conversion ratio(FCR)=Feed ingesta(g)∕Increase in fish biomass(g).

    At the end of the trial,all fish were individually captured,anaesthetised (100 mg/L MS222) and their BW and SL measured as previously described.In addition,fifteen fish per tank were euthanised (300 mg/L MS222) and their digestive tract removed for further transcriptomic and microbiome studies as described below.In this sense,the anterior-mid intestine was chosen for transcriptomic purposes because of its immunological relevance in comparison with other intestinal sections [32,33].

    Histological and histochemical analysis

    Samples from the anterior-mid intestine (n=12 randomly selected fish per experimental diet;3 fish per tank) were embedded in paraffin and sagittally sectioned (5-6 μm).A total of 576 sections (2 per each sample × 24 samples × 12 histochemical techniques)were used for histological and histochemical purposes.Two sections per each sample were stained with haematoxylin-eosin for descriptive purposes;the rest were used for evaluating the histochemical properties of epithelial and mucous cells.In brief,Schiff,Periodic Acid Schiff (PAS),diastase-PAS and Alcian Blue (AB) pH 2.5,1 and 0.5 (carboxylated and sulphated glycoconjugates/glycoproteins) techniques were used for studying carbohydrate distribution.Furthermore,several horseradish peroxidase (HRP) conjugated lectins (Sigma-Aldrich,Spain) were used for proper characterization of different glucidic residues bound to the glycoconjugates;in particular,Canavalia ensiformes/ConA (mannose and/or glucose),Ulex europeus/UEA-I (L-fucose),Triticum vulgaris/WGA (N-acetyl-D-glucosamine and/or N-acetyl-neuraminic acid,NeuNAc/sialic acid/NANA),Glycine max/SBA (a-N-acetyl-D-galactosamine) andSambucus nigra/SNA (NeuNAc/sialic acid/NANA).Lectin concentrations ranged between 15 μg/L to 30 μg/mL.Regarding negative controls,omission of the respective lectin,substitution of lectin-HPR conjugates by TBS and treatments with different enzymes were performed according to Sarasquete et al.[34].The peroxidase activity was visualized with 3,3-diaminobenzidine tetra hydrochloride/DAB and hydrogen peroxide(0.05%).All the techniques were performed according to Underwood [35] and following proper standardized techniques and protocols [36].All reagents were purchased from Sigma-Aldrich (St Louis,MO,USA).

    Histological images (600 dpi) were taken with a Leitz Wetzlar microscope with a built-in SPOT Insight Color camera (Ernst Leitz Wetzlar GmbH,Germany).Results were manually registered using a semi-quantitative assessment scoring based on color intensity scores(0,negative;1,weak;2,moderate;3,intense;4,very intense) from four independent observers,comparing the sections of the control with the experimental diet.The mucous cell count was determined in four different sites of each histological section,and the number of cells expressed per length unit of the basal lamina of the mucosal epithelium (1 mm) according to Yamamoto et al.[37].

    Intestinal microbiota

    DNA extraction and sequencing by Illumina MiSeq technology

    The complete intestine ofS.aurata(n=8 per dietary condition) was taken for microbiome analyses.Samples were stored at -20 °C,thawed gradually on ice and the mucus intestinal contents (autochthonous microbiota)were collected for further analyses.A 1-mL aliquot per sample of intestinal contents were collected with 1 mL PBS.Total DNA was extracted from each sample according to protocol described in Tapia-Paniagua et al.[38],with minor modifications.DNA concentration and purity were quantified fluorometrically with Qubit?dsDNA HS Assay Kit (Thermo Fisher Scientific,Waltham,MA,USA) and by spectrophotometric and electrophoretic methods to study the degree of purity,quality,and DNA integrity.Isolated DNA was stored at -20 °C until further processing and 30 ng used for subsequent analyses.16S rRNA of samples was sequenced on Illumina MiSeq platform (Illumina,San Diego,CA,USA) with 2 × 300 bp paired-end sequencing in the Ultrasequencing Service of Malaga (Malaga,Spain).Sequencing was carried out using the primers 341F and 805R directed to the variable regions V3-V4 of the 16S rRNA gene [39].Raw read sequences of the 16S rRNA gene fromS.aurataintestinal microbiota are publicly available in the NCBI SRA depository (SRS16609731—SRS16609735) within BioProject PRJNA914904,with BioSample accession numbers SAMN32318396-SAMN32318408.

    Bioinformatic analysis

    Primers were removed and quality control were performed withcutadaptandFastQCsoftware respectively.Barcoding sequences were processed using RDADA2library.Briefly,forward and reverse reads were truncated with decreasing quality metrics while maintaining sequence overlap.Paired reads were assembled after error modelling and correction,creating amplicon sequence variants (ASVs).Chimeric ASVs were removed by reconstruction against more abundant parent ASVs.Taxonomy was assigned to representative sequence variants using SILVA release 138 database,clustered at 99%identity,and trimmed to the amplified region.Abundance of ASV of the intestinal microbiota was processed usingphyloseqandveganlibrary in R statistical package.Shannon’s alpha diversity indices,to evaluate taxonomic diversity and PCoA to beta diversity were performed.ASVs with abundance minor to 10 reads in at least of 10%of samples were filtered to taxonomy analyses.R libraryggplotwas used to represent the abundance of different taxonomy categories and differential abundance of taxa was carried out using the R package DESeq2 with a false discovery rate (FDR < 0.05).

    Gut transcriptional analysis

    RNA isolation and quality control

    Total RNA from the anterior-mid intestine of nine randomly selected fish per dietary treatment (n=4 fish per tank) was extracted using the RNeasy?Mini Kit (Qiagen,Germany).Total RNA was eluted in a final volume of 35 μL nucleasefree water and treated with DNAse (DNA-free?DNA Removal Kit;Invitrogen,Lithuania).Total RNA concentration and purity were measured using Nanodrop-2000?spectrophotometer (Thermo Scientific,USA) and stored at -80 °C until analysis.Prior to hybridization with microarrays,RNA samples were diluted to 133.33 ng/μL concentration,checked for RNA integrity (Agilent 2100 Bioanalyzer;Agilent Technologies,Spain) and selected by the criteria of a RIN value > 8.5.Three different pools of samples per dietary treatment were established (n=4 fish each pool);this choice of pooling individual samples resulted in the loss of information derived from individual sample variability.

    Microarray hybridization and analysis

    Anterior-mid intestine transcriptional analysis from both experimental groups was carried out using the AquagenomicsSparus aurataOligonucleotide Microarray v2.0(4 × 44 K) (SAQ) platform.Detailed information and transcriptomic raw data are available at the Gene Expression Omnibus (GEO) public repository at the US National Center for Biotechnology Information (NCBI),accession numbers GPL13442,and GSE162504,respectively.Sample labelling,hybridization,washes,and scanning was performed as described in Firmino et al.[40].Briefly,a onecolor RNA labelling was used (Agilent One-Color RNA Spike-In kit;Agilent Technologies,USA).RNA from each sample pool (200 ng) was reverse-transcribed with spikein.Then,total RNA was used as template for Cyanine-3(Cy3) labelled cRNA synthesis and amplified with the Quick Amp Labelling kit (Agilent Technologies).cRNA samples were purified using the RNeasy?micro kit (Qiagen).Dye incorporation and cRNA yield were checked(NanoDrop ND-2000?spectrophotometer).Then,Cy3-labeled cRNA (1.5 mg) with specific activity > 6.0 pmol Cy3/mg cRNA was fragmented at 60 °C for 30 min and hybridized with the array in presence of hybridization buffer (Gene expression hybridization kit,Agilent Technologies) at 65 °C for 17 h.For washes,microarrays were incubated with Gene expression wash buffers,and stabilization and drying solution according to manufacturer instructions (Agilent Technologies).Microarray slides were then scanned (Agilent G2505B Microarray Scanner System),and spot intensities and other quality control features extracted (Agilent Feature Extraction software version 10.4.0.0).The extracted raw data were imported and analysed with GeneSpring (version 14.5 GX software,Agilent Technologies).To standardize the arrays intensity,75% percentile normalization was used,and data were filtered by expression levels.The differential expressed genes (DEGs) were obtained from a gene-level differential expression analysis.Expression values with up-value < 0.05 were considered statistically significant.The DEGs were grouped according to their fold-change value (P-value < 0.05) and represented using the GraphPad softwarev7.0 for Windows.The 3D principal component analysis (PCA) was carried out using GeneSpring software(Agilent),four eigenvectors were calculated to describe the aggrupation of diets control and yeast in a 3D plot.The gene expression values (log-expression ratios) were represented by a hierarchical clustering heatmap analysis using GeneSpring software (Agilent).

    The Search Tool for the Retrieval of Interacting Genes(STRING) public repository version 11.0 (https://stringdb.org),and the application ClueGO (v2.5.9) and CluePedia (v1.5.9) using Cytoscape (v3.9.1) software,were used to generate the anterior-mid intestine transcripteractome and cluster aggrupation for the DEGs of fish fed the yeastsupplemented diet.A protein-protein interaction (PPI)Networks Functional Enrichment Analysis for all the differentially expressed genes (DEGs) was conducted with a medium-confidence interaction score (0.4) using Homo sapiens as model organism [40].Gene Ontology (GO)and Kyoto Encyclopedia of Genes and Genomes (KEGG)enrichment analysis of all the DEGs obtained were also assessed through STRING (P<0.05).In order to confirm match of gene acronyms between bothHomo sapiensandS.aurataspecies,human orthology identification based on gene/protein name was accessed through the Genecards (www.genec ards.org) [41] and Uniprot (www.unipr ot.org) databases.Additionally,protein-protein BLAST(BLASTp) was run (E < 10-7;query cover > 95%).

    Statistical analysis

    Data on growth,feed performance indicators and goblet cell density in the intestinal epithelium are expressed as mean ± standard deviation,and a probability value ofP<0.05 was considered as significant by means of a Student’st-test.Regarding microbiota studies,F diversity indexes were compared between both dietary groups by means of a Student’st-test (P<0.05).

    Results

    Growth and feed performance indicators

    At the end of the trial,fish fed low fish meal diet supplemented withD.hansenii(1.1%,17.2 × 105CFU;Yeast diet) had significant differences compared with the fish fed the Control diet in terms of somatic growth.Data showed improved BWf,BWG and the SGR values in fish fed the Yeast diet compared with those fed the Control diet (Table 2;P<0.05).Furthermore,FCR values were significantly lower in fish fed the Yeast diet when compared to the control ones (1.26 ± 0.07 vs.1.06 ± 0.05,respectively).No significant differences in the FI and SR were found between both groups (Table 2;P>0.05).

    Table 2 Growth and feed performance indicators in gilthead sea bream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet)

    Morphological organization of the anterior-mid intestine and histochemical properties of mucins produced by goblet cells

    The mucosal layers from the anterior-mid intestine of fish fed both experimental diets were similar.In brief,both were lined by a simple columnar epithelium with basal nuclei,basophilic cytoplasm,and prominent brush border,with scattered goblet cells.No signals of histological alterations associated with inflammatory processes nor lymphocyte infiltrations in the anterior-mid intestine were observed in fish fed the Yeast diet compared to the control group (Fig.1a and b).No significant differences in goblet cell density were detected when comparing both treatments (Fig.2).

    Fig.1 Histological organization of the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) (a) or devoid of the probiotic (Control diet) (b).Histochemical properties of mucins produced by intestinal goblet cells regarding their content on carboxylated and/or sulphated acidic groups from fish fed the Control diet (c) and the Yeast diet (d).Abbreviations: mc,mucous cells;bb,brush border.Staining: hematoxylin-eosin (a,b),Alcian Blue pH=2.5 (c,d).Scale bar=50 μm

    Fig.2 Box-plot graph of goblet cell density in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).No statistically significant differences were detected among dietary treatments (Student’s t-test;P>0.05).Mean line within the boxplot indicate the average value of the series;O and =outlier points of the series

    Regarding the histochemical properties of mucins produced by goblet cells from the anterior-mid intestine,results showed a variable richness of neutral glycoproteins (PAS and diastase-PAS positive) (Table 3).In addition,mucins showed a mixture of carboxylated (AB pH=2.5) and sulphated acidic groups (weak and strongly ionized;AB pH=1.0 and 0.5,respectively) (Table 3 and Fig.1c).Regarding lectin histochemistry,a specific affinity for WGA,SNA,ConA and SBA lectins was detected in the mucinous content of goblet cells (Table 3 and Fig.3).Moreover,no variations were detected in the distribution of mucosal cell glycoconjugates between the upper and the bottom areas of the intestinal folds.When comparing both dietary groups,the dietary administration ofD.hanseniimodified the composition of glycoproteins of mucins produced by goblet cells,with an increase in the staining intensity of mucins rich in carboxylated and weakly sulphated glycoconjugates (Table 3 and Fig.1c and d).In addition,an increase in affinity for SBA lectin and a decrease in the affinity for the WGA lectin was found in the mucinous content of goblet cells,whereas no changes were detected regarding ConA,SNA and UEA-I lectins (Table 3 and Fig.3).

    Fig.3 Histochemical localization of glycoconjugates containing sugar residues in the anterior-mid intestine of gilthead seabream (Sparus aurata)fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).Presence of glycoconjugates containing α/β-N-acetyl-D-galactosamine in mucous cells of S.aurata fed a Control (a) or a Yeast diet (b).Results denote a moderate increase in affinity for the SBA lectin in the mucous cells and in the intestinal epithelium from Yeast diet.Glycoconjugates containing N-acetyl-D-glucosamine and/or N-acetylneuraminic acid residues in mucous cells of S.aurata fed the Control or the Yeast diet(d).Results denote a moderate decrease in affinity for the WGA lectin in the mucous cells from the Yeast diet.Histochemical detection of glycoconjugates containing N-acetylneuraminic acid/sialic acid residues in the intestine from the Control (e) or Yeast (f) groups.Note the decrease in affinity for the SNA lectin in the anterior-mid intestinal epithelium of S.aurata fed the Yeast diet.Glycoconjugates containing a-mannose/a-glucose residues in anterior-mid intestine from S.aurata fed Control (g) or Yeast diets (h).Observe the increase in affinity for the ConA lectin in the intestinal epithelium of S.aurata fed the Diet Y.Abbreviations: e: epithelium;mc: mucous cells.Scale bar=50 μm

    Table 3 Histochemical results and affinity for different lectins in the anterior-mid intestinal mucosa of gilthead seabream(Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet)

    Intestinal microbiota

    Sequencing data analysis

    After bioinformatic processing,an average of 117,038.2 reads per sample were obtained (minimum reading 71,378 and maximum 162,602 reading) that were classified into 1,594 ASVs.

    The Shannon and Simpson alpha diversity indexes were calculated for each sample of experimental groups(Table 4).The alfa diversity values in the experimental groups did not show significant differences (P>0.05).The PCoA analysis to check the beta diversity between the detected sequences did not show a clearly separated grouping of samples based on the diet (Additional file 1:Fig.S1).

    Table 4 Alpha diversity indices of microbiota intestine population of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet)

    Taxonomic composition of the intestinal microbiota

    At phylum level,the predominant phyla detected in all fish fed both experimental diets were Proteobacteria,followed by Firmicutes (Fig.4A),whereas the most abundant classes in all cases were alpha-and gamma-Proteobacteria(Fig.4B).However,no significant differences were found among diets.A comparative analysis at genus level (abundance > 0.5%) showed a very similar abundance of sequences between specimens from both diets (Fig.5),whereasPseudoalteromonas,Acinetobacter,Pseudomonas,ShewanellaandVibriowere the most abundant genera.

    Fig.4 Values of relative abundance at the phylum (A) and class (B) level,obtained for the anterior-mid intestine microbiota of gilthead seabream(Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic(Control diet)

    Fig.5 Values of relative abundance at genus level,obtained for the anterior-mid intestine microbiota of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet)

    DESeq2 analysis was applied to evaluate the significant differences (FDR < 0.05) between the abundance of the microbial taxa.Significant differences were detected regarding to the genera included in Protobacteria,and with abundance higher than 0.5% such asAcinetobacter,ComamonasandPseudomonas,which showed decreased abundances in fish from the Yeast diet in comparison to those fed the Control diet(Fig.6).Other significant differences were also detected within the Proteobacteria phylum,and particular in bacteria from the genusAnaerococcus,Ascidiaceihabitans,HydrogenophagaandVariovoraxthat displayed lower abundances of ASVs in fish from the Yeast diet,whereas an increase of ASVs related to the genusBacilluswere found in the intestine of fish fed the Yeast diet(Fig.6).

    Fig.6 Logarithmic fold change of differentially abundant ASVs (DESeq2,FDR < 0.05) at the genus level in the anterior-mid intestine microbiota of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet)

    Transcriptomics of the anterior-mid intestine

    Transcriptomic profile

    To evaluate the modulatory effects of the dietary administration ofD.hanseniiat 1.1% upon intestine transcriptome,we performed a microarray-based transcriptomic analysis.As shown in Fig.7 and listed in Additional file 1: Table S1,a total of 232 differentially-expressed genes (DEGs) were found in the anterior-mid intestine when comparing samples from both diets (P<0.05).Among these,177 DEGs were up-regulated,and 55 DEGs were down-regulated showing both different levels of fold-change (FC) (Fig.7A).The percentage expression pattern was similar in both up-and down-regulated genes,although in numerical terms most of the DEGs found were up-regulated with a FC > 2 (12 DEGs upregulated vs.1 DEG down-regulated),even some of these being expressed with a FC of 4,6 or 10 (Fig.7A and Additional file 1: Table S1).The 3D PCA analysis showed a clear grouping for both the Control and Yeast diets(Fig.7B),meanwhile the hierarchical clustering heatmap confirmed the conservative response of the arrays,showing a clear differential expression profile of DEGs between both dietary groups (Fig.7C).

    Fig.7 Transcriptomic response of the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).A The green (up-regulated) and red(down-regulated) colour scheme indicate the gene modulation according to its magnitude interval (fold-change).B 3D principal component analysis (PCA).Control (red spheres) and D.hansenii groups (yellow spheres) are represented (each array included in the study is represented by one sphere).C Hierarchical clustering heatmap representing the 232 DEGs.The results of each microarray analysed for the control and yeast groups are shown.The green (up-regulation) and red (down-regulation) colour scheme indicates the gene normalized intensity values according to its magnitude interval

    A biological network (transcripteractome) considering the complete DEGs annotated list was prepared in order to have an overview of the relationships between them (Fig.8).The gene network resulted in 157 nodes/DEGs and 210 edges/interactions (with an expected number of 145 edges),excluding from the analysis the remaining DEGs (n=75;32% of the total DEGs),classified as unknown genes.

    Fig.8 Biological network of the differentially expressed genes (DEGs) in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).Node fill with a continuous mapping colour from red (down-regulated DEGs) to green (up-regulated DEGs),represents the fold change (FC) intensity.Edge score was represented by doted lines with a continuous mapping from 0.35 (light Blue) to 1.00 (dark blue).The transcripteractome was obtained using Cytoscape (v3.9.1) platform with String (v1.7.1).For details,please refer to Additional file 1: Table S1

    Functional enrichment analysis

    The functional enrichment analysis of the transcripteractome helped us to identify several biological processes and pathways for the DEGs in the anterior-mid intestine ofS.auratafed the Yeast diet in comparison to their congeners from the control group (Control diet).Thus,54 GO were identified (P<0.05;Additional file 1: Table S2),whereas the top-10 representative processes per number of DEGs were as listed: “Nitrogen compound metabolic process” (60.5% of the total DEGs: 83 up-regulated vs.12 down-regulated;GO:0006807),“Organic substance biosynthetic process” (30.6% of the total DEGs:42 up-regulated vs.6 down-regulated;GO:1901576),“Cellular response to chemical stimulus” (28.2% of the total DEGs: 35 up regulated vs.9 down-regulated;GO:0070887),“Nucleic acid metabolic process” (24.8%of the total DEGs: 37 up-regulated vs.2 down-regulated;GO:0090304),“Negative regulation of gene expression”(20.4% of the total DEGs: 27 up-regulated vs.5 downregulated;GO:0010629),“Cellular catabolic process”(18.5% of the total DEGs: 28 up-regulated vs.1 downregulated;GO:0044248),“Heterocycle biosynthetic process” (14.0% of the total DEGs: 21 up-regulated vs.1 down-regulated;GO:0,018130),“Cellular response to oxygen-containing compound” (13.4% of the total DEGs:16 up-regulated vs.5 down-regulated;GO:1901701),“DNA metabolic process” (12.1% of the total DEGs: 17 up-regulated vs.2 down-regulated;GO:0006259) and,“Symbiotic process” (11.5% of the total DEGs: 15 up-regulated vs.3 down-regulated;GO:0044403) (Additional file 1: Table S3).

    According to the results reported by the enrichment analysis,representative clusters of genes were classified as shown in Figs.9,10,11,12 and 13 and Additional file 1:Table S4.A significant number of genes were observed to be shared among the metabolic category: Protein Metabolism (31.2% of known genes;Fig.9),Nucleotide Metabolism (8.3% of known genes;Fig.10),and Lipid Metabolism (5.7% of known genes;Fig.11).The other clustered categories contained genes related to the Regulation of Antioxidant Capacity (6.4% of known genes;Fig.12),and to the Regulation of Immune Response(10.8% of known genes;Fig.13).This clusters suggest a strong relationship between processes that favours metabolic and immune responses to the dietary inclusion ofD.hansenii.

    Fig.9 Protein metabolism cluster of the functional enrichment network for biological processes based on the related functions of the differentially expressed genes (DEGs) in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).The node size determinate the significance(Bonferroni step-down correction;small P<0.05;medium P<0.005;large P<0.0005),the node colour determinate the principal GO process aggrupation,and the node numbers indicate different GO process classified in Additional file 1: Table S4

    Fig.10 Nucleotide metabolism cluster of the functional enrichment network for biological processes based on the related functions of the differentially expressed genes (DEGs) in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).The node size determinate the significance (Bonferroni step-down correction;small P<0.05;medium P<0.005;large P<0.0005),the node colour determinate the principal GO process aggrupation,and the node numbers indicate different GO process classified in Additional file 1: Table S4

    Fig.11 Lipid metabolism cluster of the functional enrichment network for biological processes based on the related functions of the differentially expressed genes (DEGs) in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).The node size determinate the significance(Bonferroni step-down correction;small P<0.05;medium P<0.005;large P<0.0005),the node colour determinate the principal GO process aggrupation,and the node numbers indicate different GO process classified in Additional file 1: Table S4

    Fig.12 Regulation of Antioxidant Capacity cluster of the functional enrichment network for biological processes based on the related functions of the differentially expressed genes (DEGs) in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).The node size determinate the significance (Bonferroni step-down correction;small P<0.05;medium P<0.005;large P<0.0005),the node colour determinate the principal GO process aggrupation,and the node numbers indicate different GO process classified in Additional file 1: Table S4

    Fig.13 Regulation of Immune Response cluster of the functional enrichment network for biological processes based on the related functions of the differentially expressed genes (DEGs) in the anterior-mid intestine of gilthead seabream (Sparus aurata) fed experimental low fishmeal diets containing Debaryomyces hansenii (1.1%,17.2 × 105 CFU;Yeast diet) or devoid of the probiotic (Control diet).The node size determinate the significance (Bonferroni step-down correction;small P<0.05;medium P<0.005;large P<0.0005),the node colour determinate the principal GO process aggrupation,and the node numbers indicate different GO process classified in Additional file 1: Table S4

    Three significative cellular components GO groups(80.2% of the total DEGs in GO:0 043227;73.9% in GO:0 043231;10.2% in GO:0005740) were associated with the single or double lipid bilayer membrane,which indicated that most of DEGs belong to the above-mentioned structures.On the other hand,it is also worth noting that some of these DEGs participate in mitochondrial processes (GO:0005739),proteasome complex (GO:0000502) and DNA replication (replication fork: GO:0 005657;replisome: GO:0030894).For further details see Additional file 1: Table S5.

    Discussion

    Low fishmeal along with high plant-content meal diets are of common practices for the aquaculture feed industry.However,reducing the dependence on marine raw materials is still associated with variable side effects on fish nutrient metabolism,mucosal health,and disease resistance [42],especially due to the presence of antinutritional factors and the lack of essential amino acids and fatty acids in those alternative sources.Thus,it is necessary to develop complementary nutritional tools to reduce the negative-side effects associated with the substitution of traditional ingredients.In this sense,the use of functional feeds seems to be an appropriate strategy [43].When considering new additives,yeasts are among those promising candidates showing a wide range of benefits on animal health and growth [21,44].In this sense,the present study aimed to determine the effects on growth and intestinal condition of a functional feed with a low fishmeal formula (7%) that was supplemented withD.hanseniiinS.aurata,the main marine fish species farmed in the Mediterranean basin [30].

    Effects of D.hansenii on fish performance

    Under present experimental conditions,we found a significant body weight increase in fish fed the Yeast diet coupled to a lower feed conversion rate and similar feed intake values compared with those fish fed the Control diet.Considering that feeding costs represent up to 50%-70% of total production costs in intensive fish farms [45],and fishmeal being amongst the most expensive raw materials in aquafeeds [46];the nutritional strategy tested in our study that consisted of combining a diet containing low fishmeal levels (7%) with the supplementation of a live probiotic seemed very promising when KPIs related to growth and feed utilization were considered.Furthermore,these results are of special relevance since this is the first study in fish juveniles showing the benefits of supplementing low fish meal diets withD.hanseniiin terms of somatic growth and feed efficiency parameters.In this context,several previous studies had reported that the dietary inclusion ofD.hansenii(intact or hydrolysed) in aquafeeds with high levels of fishmeal (39%-41%) did not result in significant differences in terms of fish growth performance in Atlantic salmon (Salmo salar) and tropical gar(Atractosteus tropicus) [27,47].Present findings when using low fishmeal diets may be due to the presence inD.hanseniiof growth factors like biogenic amines,capable of promoting growth and intestinal functionality [21,29,48].Biogenic amines or polyamines are essential metabolites found in almost all living organisms [49].In particular,D.hanseniiis especially rich in putrescine,spermidine,spermine and cadaverine[17,50].These polyamines are essential for many cell functions and play a key role in vertebrate development and cell growth [17,51].In particular,Tovar-Ramírez et al.[25] observed in European sea bass (Dicentrarchus labrax) an increase of larval survival and an enhancement of the digestive function,whereas Teles et al.[23]found that longfin yellowtail (Seriola rivoliana) larvae had higher growth and survival rates,and an earlier digestive tract maturation when fed with a diet supplemented withD.hansenii.These results found in fish larvae were similar to those found in the current study regardless of the developmental stage considered.Furthermore,although it is described that polyamines in high concentrations can generate adverse effects on animal health,their probiotic administration throughD.hanseniiin feed might be of great benefit to fish as reported in several studies [21,25,29,48].In addition,by participating in the protein metabolism,these polyamines could be of special interest in carnivorous species with high protein demand [51],since their dietary provision byD.hanseniimay modulate the amino acid metabolism in fish,which could result in the promotion of fish growth by a better use of dietary proteins [52],as shown in our study.

    Effect of D.hansenii diet on intestinal condition

    The intestine is a complex organ composed of the digestive epithelium with its typical tissue organization in different mucosal layers,its associated lymphoid tissue,and the mucus layer covering the simple intestinal epithelium with its commensal microbiota.Furthermore,enterocytes form a selectively permeable barrier for water,electrolytes,and nutrients,while maintaining an effective defence against pathogens and tolerance towards dietary antigens [53].Thus,studying the effect of the diet on tissue condition and its commensal microbiota is of importance when evaluating new nutritional strategies.In this context,the intestinal microbiota has been extensively studied in a multitude of fish species since microbial gut communities are essential for proper gastrointestinal development and are implicit in aiding digestive and immunological functions and disease resistance [54,55].However,how host-microbiota symbiosis and/or how feed and feed additives can modulate the normal microbiota in the host remains poorly understood,since gut microbial communities are strongly affected by dietary or environmental factors [56,57].Despite that,results in terms of phyla obtained in our study were in agreement within those from other authors [40,58,59].In those studies,Proteobacteria has been reported as the most abundant phylum in the intestine ofS.aurata,followed by Firmicutes,Actinobacteria,and Bacteroidetes.This conservative pattern at phyla level of microbial communities in the intestine was maintained in our study regardless of the diet evaluated.

    At a lower taxonomic level,Pseudomonasspp.andAcinetobacterspp.were the dominant genera in the gut microbiome from both dietary groups.Pseudomonasspp.andAcinetobacterspp.are a part of the typical fish microbiota [60],whereas some strains have been reported as antagonistic against pathogenic microorganisms [61,62].However,under stressful conditions such as malnutrition and overcrowding their abundance and the presence of some opportunistic and pathogenic species may increase [63-65].Under current experimental condition,the reduction in the relative abundance of these genera in fish fed the Yeast diet suggested that the dietary administration ofD.hanseniimay be a strategy for dealing with infections caused by opportunistic and antimicrobial-resistant bacteria,such asPseudomonasspp.andAcinetobacterspp.[66-70].Furthermore,D.hanseniisupplementation also resulted in an increase of the relative abundances of ASVs related to the genusBacillusin comparison to fish fed the Control diet.This Gram-positive endosporeforming genus has been widely tested as a probiotic on several aquatic cultured organisms,enhancing growth,modulating immune response,and promoting disease resistance [71,72].Similar to other species,D.hanseniiare normally found as cohabitants in the gut of wild and farmed fish [21,73],and seems to interact with the host,maintaining intestine microbiota homeostasis,and aiding to its recover in case of dysbiosis [74,75].Although in the current study,the low fishmeal diet did not result in any dietarily-induced inflammatory disorder of the intestinal mucosa or dysbiosis of the gut microbiome,these interactions in terms of DEGs related to membrane-linked organelles and symbiotic biological processes were identified by the transcriptomic functional analysis conducted,and further discussed in "Effect ofD.hanseniidiet on gut metabolism,homeostasis,and immunity".

    The mucus layer that covers enterocytes is considered as the first line of innate host defence as well as the site where commensal microbiota attaches,preventing the colonisation of pathogenic organisms by occupying empty niches.The dynamic relationship between commensal bacteria and mucins secreted by goblet cells plays an important role in gut health and condition [76].Commensal microbiota is adapted to the glycan rich environment of mucus,stablishing host-microbiota interactions,influencing on the quantitatively and qualitatively mucin glycosylation [77].In the current study,the dietary administration ofD.hanseniidid not modify the histochemical properties of mucins produced by goblet cells in terms of neutral glycoproteins,but resulted in a slightly higher richness of carboxylated and weakly ionised sulphated glycoconjugates in fish fed the Yeast diet.These results are of special relevance since the pattern of O-glycosylation in mucin glycoproteins and their polypeptide backbone structure linked to some physical properties of the mucus layer as its viscosity and depth,are key elements in gut health and condition modulation [78].Thus,the increase in carboxylated and weakly ionized acid sulphated grlycoproteins produced by goblet cells in fish fed the Yeast diet may be interpreted as beneficial due to an increase in mucus viscosity that would protect the host form opportunistic pathogenic bacterial [79,80].Furthermore,the inclusion ofD.hanseniiin the diet changed the glycosylation patterns of intestinal mucins,being of special relevance the increase in SBA lectin affinity rather than WGA,changes that may be partially associated to different rates of mucin glycosylation and goblet cell differentiation along the intestinal mucosa [34].However,these results are of special relevance since lectins with affinity for SBA,which bind to GalNAc residues may contain sialic acid residues in the backbone region of the mucin [81],have been postulated to have an antiparasitic effect in turbot (Scophthalmus maximus) [82] and pathogenic organisms [78,83].Furthermore,the abovementioned changes in lectin affinity may be supported by differential expression ofb4galt3andb3galt2genes,which both participates in the elongation and branching biosynthetic steps of glycosylation pathways [84],and the up-expression of a Golgi mannosyl-oligosaccharide 1,2-α-mannosidases gene (man1a1),whose activity is a prerequisite for the formation of complex or hybrid glycans of mucins [85].

    Overall,the results obtained showed the ability ofD.hanseniito modulate the gut microbiota without affecting the intestine cell organization inS.aurata,mainly reducing several groups of Proteobacteria abundances,especially those opportunistic groups,and changing glycosylation patterns of intestinal mucins.It is suggested that these results,which can be attributed to the potential role ofD.hanseniias a probiotic,involve some of the yeast cell wall components,such as β-glucan,polyamines and mannan molecules that have,apart from these attributes,the ability to modulate the response of the immune system [21,86,87].However,the functional benefits of these findings need further attention,including the role ofD.hanseniiand the microbiota or their interactions,in terms of glycosylation patterns variation.

    Effect of D.hansenii diet on gut metabolism,homeostasis,and immunity

    In parallel with the improved performance and the modulation of the intestinal condition,several authors have suggested that the dietary administration ofD.hanseniiresult in a concomitant modification of systemic physiological responses,and especially those related to the digestive function and immunity [17,88].In this sense,the dietary inclusion ofD.hanseniimodulated the expression of different genes in the anterior-mid intestine ofS.auratathat were associated with biological processes related to gut metabolism,antioxidant defence,and immunity.

    Compared with terrestrial animals,fish have higher protein requirements due to their biological lifestyle characteristics [89].For this reason,it is not surprising that most of the DEGs with metabolic functions were mainly related to protein biosynthesis,turnover and degradation to cope with their increased growth rate.In our study,we found that the ubiquitin-proteasome system(UPS) was modulated by the regulation of some ubiquitin enzyme and ligases genes (uchl3,trim 25,trim39,ubr7andube2i) and proteasome subunits (psma4,psma6,psmb10,psmc2andrad23a),which act as protein binding signalling and protein degradation,respectively.Both activities are suggested to be important mechanisms of aging,immune activity,cell quality control and protection,tissue regeneration,or DNA repair and gene transduction [90,91].On the other hand,we found several DEGs related to protein synthesis and modification,especially genes involved in ribosomal subunit components (rpl22,rpl38,mrps2,mrps9,andmrpl54),and amide-related biosynthetic processes.The regulation of several genes involved in amide biosynthesis may be attributed to the amine and polyamine content ofD.hanseniiincluded in the diet [17],which would support the above-mentioned benefits in terms of growth performance as described in several species [21].

    Linked to the increased gut metabolic activity,we found other DEGs not directly related to protein metabolism that were associated with nucleotide metabolism and DNA elongation and repair.According to current results,dietaryD.hanseniienhanced thymidylate nucleotide metabolism as indicated by the up-regulation ofdtymk,tyms,dut,anddckgenes.Thymidylate is necessary to synthetise deoxythymidine triphosphate (dTTP),an essential building block of DNA [92].In particular,studies in larvae showed that TYMK,the protein codded bydtymkgene,and its deoxyribonucleotide triphosphate(dNTP) product,are vital in early embryogenesis and essential for neurodevelopment,and genome integrity in adult fish [93,94].These results would be in agreement with the reported beneficial effect ofD.hanseniilarval development and morphogenesis [23,25].

    The third main metabolic cluster identified in our study was related to lipid metabolism,mainly involving DEGs related to sphingolipid or ceramide metabolic processes(aldh3a2,elovl6,asah1,degs1,b3galt2andb4galt3).Sphingolipids and ceramides are constituents of eucaryotic cell membranes with a wide variety of functions and properties.It is reported that sphingolipids have antimicrobial and immunomodulatory properties [95],along with bioactive signalling functions [96].The DEGasah1encodes for an acid ceramidase N-acylsphingosine amidohydrolase 1 (ASAH1),which catalyses the degradation of ceramide into sphingosine and free fatty acid and suggested to maintain ceramide homeostasis in the lysosome[97].Furthermore,participates in sphingosine recycling pathway along withdegs1gene,which controls the step from dihydroceramides to ceramides,and suggested to be a promising biomarker of redox state [98].

    An increase in the metabolism could lead to increased production of oxidative radicals,that in an unbalanced situation with the antioxidant defence capacity would produce oxidative damage [99].In parallel to the abovementioned enhancement of gut metabolic activity,we also found an increased expression of some genes from the electron chain and others also involved in mitochondrial respiration.Their activity produces chemical energy for the normal cell behaviour,but as a result the superoxide anion (),the precursor of most other reactive oxygen species (ROS),are also obtained [100].To possibly counteract with this situation,we found the up-regulation ofmgst3andgsto1genes in fish fed the Yeast diet,which both together are related with the glutathione cycle,in charge of detoxifying hydrogen peroxide (H2O2) [101],andselenot,which mediates oxidoreductase functions involved in redox homeostasis[102].Magalh?es et al.[103] reported inS.auratadifferent antioxidant levels depending on the organ analysed,being the enzymatic activity of catalase (CAT),superoxide dismutase (SOD) and glutathione reductase (GR) higher in the intestine,and glutathione peroxidase (GPx) in liver.The authors discussed that this pattern is expected due to the elevated cell turnover rate in the intestine that favours oxidative processes in a non-stressful situation,which is also proposed in our study by the supposedly increased gut cell metabolism linked to the increased somatic growth observed.As H2O2plays important cellular functions,e.g.,as a promoter for cell cycle progression,the activity of the glutathione cycle enzymes and cytochrome P450 (cyp1a1),also takes part in metabolic pathways like arachidonic,leukotriene,prostaglandin metabolism,or in estrogen metabolism [51,104],as well as protecting lipids and proteins to be oxidised [105].In that way,the upregulation of malate/oxoglutarate carrier (slc25a11)that is related with the TCA cycle and also plays a key role maintaining GSH levels in the mitochondria [106],would reinforce this hypothesis.Thus,the up-regulation of these genes could be more related to metabolism than to a possible oxidative damage produced by the inclusion of yeast in the diet.This idea is also supported by some down-regulated genes related to oxidative stress like:serpinb1,which is suggested to protect the cell from proteases released into the cytoplasm during stress [107];sin3a,whose expression is promoted under hypoxia stress [108];and the anti-inflammatory cytokine IL-10,that is modulated under redox/oxidant perturbations [109].Regarding antioxidant functions,dietary live yeast administration inD.labraxlarvae,induced lower activity and expression levels of glutathione peroxidase and superoxide dismutase compared to fish fed control diet,assuming a possible involvement of superoxide anion retention in fish larvae,which could represent importance to the host to increase cell or tissue responsiveness to growth-and/or differentiation-enhancing factors [110].

    As reviewed in Angulo et al.[21],D.hanseniidietary supplementation has been shown to promote gut protection in several aquatic and terrestrial species at immunological and gene expression levels.An exhaustive analysis of the DEGs promoted by the Yeast diet in our study showed different regulatory responses of theS.auratagut immune system.Among them,we found regulatory viral genes up-regulated like thelsm14a,a sensor for both viral DNA and RNA which is important to initiate the IFN-β induction [111];thetrim25gene,that encodes the TRIM25 E3 ubiquitin ligase protein,which mediates an essential step to initiate the antiviral responses in the cell through the pattern recognition receptor RIG-I [112];and thephb2,a gene that encodes a prohibitin 2 protein related to the mitochondrial-mediated antiviral innate immunity by activating the RIG-I signalling pathway [113].Yeast dietary inclusion also promoted the up-regulation of regulatory genes required for T and B cell development and proinflammatory mediators (nsd2,rabl3,inhbb,andparp1),a lymphocyte activator (slamf8),and genes involved in proteasomal generation of antigenic peptides (psma4andpsma6).All these genes might participate in the first steps of the immune response,promoting the regulation of sentinel immune responses,maintaining the immune cells quiescent to avoid autoimmunity,but prepared for a rapid activation [114].Contrary to those up-regulated genes,we found down-regulated genes directly related with the activation of immune system processes:serpinb1,which is suggested to be activated upon an infection[115];thetyrobpgene,which participates in the activation of immune cells involved in inflammatory reactions and as a component of NK cells with anti-viral functions[116];and thecfhgene,which regulates the alternative pathway of the complement system [117].We also found two cytokines down-regulated,a pro-inflammatoryil15gene,which promotes Th1 responses and T-cell maturation [118],andil10,an anti-inflammatory gene that is expressed under infection with stimulatory and inhibitory activities on several immune cells [119].Altogether,it may be postulated that yeast supplementation in a low fish meal feed could promote the regulation of sentinel immune processes by increasing the potentiation of the innate immunity rather than activating a response of the immune markers.A good example of this hypothesis was shown by the down-regulation ofpsmb10,which is an IFN-inducible immunoproteasome,mainly activated during an immune response [120].In that sense,it has been postulated thatD.hanseniimight modulate host’s immunity through its wall-related β-glucan and polyamine content,increasing functional and decreasing deleterious immune responses in fish [121].How these components act or in which pathways are involved is still poorly studied.However,as β-glucans are considered pathogenassociated molecular patterns,there seems to be a consensus in its effects potentiating cellular responses and gene expression signalling pathways,activating the communication and activity of the adaptive immune system in several marine species [122-124],which could explain the vast response of the intestine at transcriptomic level.

    Conclusions

    The present multi-integrative study of the intestinal response to the probioticD.hanseniiwhen administered toS.auratajuveniles in diets with low fishmeal levels provided new insights into the use of this yeast and its metabolic and immunomodulatory components on farmed fish.Under current experimental conditions,we demonstrated that dietary administration ofD.hansenii(1.1%,17.2 × 105CFU) enhanced somatic growth and improved feed efficiency parameters,results that were coupled to an improvement of intestinal condition as histochemical,microbiota and transcriptomic tools indicated.In particular,D.hanseniistimulated hostmicrobiota interactions without altering the intestinal cell organization nor generating dysbiosis,which demonstrated the safety use of the yeast as a probiotic.At the transcriptomic level,D.hanseniipromoted metabolic pathways,mainly protein-related,sphingolipid,and thymidylate pathways,in addition to enhance antioxidantrelated intestinal mechanisms,and to regulate sentinel immune processes,potentiating the defensive capacity meanwhile maintaining the homeostatic status of the intestine.These findings confirm the probiotic benefits ofD.hanseniiin the diet,and its use on cultured marine species beyond their larval stages.

    Abbreviations

    AB Alcian blue

    ASAH1 N-acylsphingosine amidohydrolase 1

    ASVs Amplicon sequence variants

    BW Body weight

    CAT Catalase

    CFU Colony-forming unit

    D.hanseniiDebaryomyces hansenii

    DEGs Differential expressed genes

    dNTP Deoxyribonucleotide triphosphate

    dTTP Deoxythymidine triphosphate

    FC Fold-change

    FCR Feed conversion ratio

    FDR False discovery rate

    FI Feed intake

    GEO Gene Expression Omnibus

    GR Glutathione reductase

    GO Gene Ontology

    GPx Glutathione peroxidase

    GSH Reduced glutathione

    HRP Horseradish peroxidase

    IFN-β Interferon beta

    IL-10 Interleukin 10

    KEGG Kyoto Encyclopedia of Genes and Genomes

    KPI Key performance indicators

    NCBI National Center for Biotechnology Information

    NK Natural killer

    PAS Periodic Acid Schiff

    PCA Principal component analysis

    PCoA Principal coordinate analysis

    PPI Protein-protein interaction

    RIG-I Retinoic acid-inducible gene I

    S.aurataSparus aurata

    SBA Soybean agglutinin

    SD Standard deviation

    SGR Specific growth rate

    SOD Superoxide dismutase

    TCA Tricarboxylic acid cycle

    TRIM25 Tripartite motif containing 25

    UPS Ubiquitin-proteasome system

    WGA Wheat Germ Agglutinin (Triticum vulgarislectin)

    Supplementary Information

    The online version contains supplementary material available at https://doi.org/10.1186/s40104-023-00895-4.

    Additional file 1: Fig.S1:Principal coordinate analysisof the anterior-mid intestine microbiota similarity.Table S1:Summary of the differentially expressed genesobtained from the microarray-based transcriptomic analysis for the gut of juvenile gilthead seabream fed theD.hanseniidiet.Table S2:Functional enrichment analysis of the DEGs.Biological processidentified in the intestine of juvenile gilthead seabreamfed theD.hanseniidiet.Table S3:Top-10 biological processof the DEGs identified in the intestine of juvenile gilthead seabreamfed theD.hanseniidiet.Table S4:Representative cluster of the functional enrichment network for biological processes based on the related functions of the differentially expressed genesin the anterior-mid intestine.Table S5:Cellular componentand annotated keywordsof the DEGs identified in the intestine of juvenile gilthead seabreamfed theD.hanseniidiet.

    Acknowledgements

    We thank the technicians at the IRTA facilities for their valuable help during fish rearing.

    Authors’ contributions

    EG designed and carried out the experiments.Biological samplings were performed by EG,JPF,AR,and IS.The transcriptomic data analysis and interpretation were performed by IS,EV-V,and FER-L.IMC and MAM carried out the microbiota analyses.JBO-D and CS carried out the histological and histochemical analysis.DT-R reviewed and validated the methodology used in the study.The study was supervised by EG,IS,and LT.IS wrote the original draft.All the authors provided the critical feedback,read,and agreed to the published version of the manuscript.

    Funding

    This work has been financed through the DIETAplus project of JACUMAR(Junta de Cultivos Marinos,MAPAMA;Spanish government),which is cofunded with FEMP funds (EU).Furthermore,this research was funded by means of grants from the Spanish Government: PID2019-106878RB-I00 and IS was granted with a Postdoctoral fellowship (FJC2020-043933-I).EVV and FERL thank the support of Fondecyt iniciación (project number 11221308)and Fondecyt regular (project number 11221308) grants (Agencia Nacional de Investigacióny Desarrollo de Chile,Government of Chile),respectively.Collaboration between Ibero-American researchers has been done under the framework of the network LARVAplus “Strategies for the development and im-provement of fish larvae production in Ibero-America” (117RT0521) funded by the Ibero-American Program of Science and Technology for Development(CYTED,Spain).

    Availability of data and materials

    Detailed information and transcriptomic raw data are available at the Gene Expression Omnibus (GEO) public repository at the US National Center for Biotechnology Information (NCBI),accession numbers GPL13442,and GSE162504,respectively.Data on intestinal microbiota are publicly available in the NCBI SRA depository (SRS16609731-SRS16609735) within BioProject PRJNA914904,with BioSample accession numbers SAMN32318396 and SAMN32318408.

    Declarations

    Ethics approval and consent to participate

    This study was conducted following the Guiding Principles for Biomedical Research Involving Animals (EU2010/63),the guidelines of the Spanish laws(law 32/2007 and RD 53/2013) and authorized by the Ethical Committee of the Institute for Research and Technology in Food and Agriculture (Spain),and Generalitat de Catalunya government for the use of laboratory animals(FUE-2020-01719167).

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1Aquaculture Program,Institute of Agrifood Research and Technology (IRTA),La Ràpita 43540,Spain.2Centro de Biotecnología Acuícola,Universidad de Santiago de Chile,Santiago,Chile.3Instituto de Ciencias Marinas de Andalucía(ICMAN-CSIC),Avda.República Saharaui n° 2,Campus Universitario Río San Pedro,Puerto Real,Cádiz 11510,Spain.4Núcleo de Investigaciones Aplicadas en Ciencias Veterinarias y Agronómicas,Facultad de Medicina Veterinaria y Agronomía,Universidad de Las Américas,Santiago,Chile.5Department of Cell Biology,Physiology,and Immunology,Universitat Autonoma de Barcelona,Barcelona,Spain.6Centro de Investigaciones Biológicas del Noroeste SC,CIBNOR,La Paz,México.7Department of Microbiology,Instituto de Biotecnología Y Desarrollo Azul (IBYDA),Faculty of Sciences,University of Malaga,29010 Malaga,Spain.8SCBI,Bioinformatic Unit,University of Malaga,29590 Malaga,Spain.

    Received:20 February 2023 Accepted:11 May 2023

    日韩国内少妇激情av| 亚洲熟妇熟女久久| 欧美不卡视频在线免费观看 | 免费高清视频大片| 国产精品影院久久| 亚洲午夜精品一区,二区,三区| xxxhd国产人妻xxx| 国产成人欧美在线观看| 69av精品久久久久久| 国产不卡一卡二| 午夜精品国产一区二区电影| 亚洲专区国产一区二区| 国产亚洲精品一区二区www| 亚洲 欧美一区二区三区| 99久久国产精品久久久| 欧美在线一区亚洲| 亚洲av日韩精品久久久久久密| 亚洲精品一二三| 91成人精品电影| 国产精品亚洲一级av第二区| 韩国av一区二区三区四区| 免费av中文字幕在线| 色婷婷av一区二区三区视频| 最近最新中文字幕大全免费视频| 国产一区二区激情短视频| 日韩欧美一区视频在线观看| 午夜激情av网站| 一级毛片精品| 精品日产1卡2卡| 日本撒尿小便嘘嘘汇集6| 99riav亚洲国产免费| 久久中文字幕一级| 午夜免费观看网址| 涩涩av久久男人的天堂| 黑人巨大精品欧美一区二区mp4| 亚洲 欧美一区二区三区| 欧美中文综合在线视频| 亚洲色图 男人天堂 中文字幕| 久久久久久久午夜电影 | 亚洲国产精品一区二区三区在线| 国产免费现黄频在线看| 久热这里只有精品99| 精品一区二区三区视频在线观看免费 | 免费在线观看完整版高清| 亚洲精品成人av观看孕妇| 久久热在线av| 一区在线观看完整版| 黄色视频,在线免费观看| 午夜免费观看网址| 纯流量卡能插随身wifi吗| 色婷婷久久久亚洲欧美| 法律面前人人平等表现在哪些方面| x7x7x7水蜜桃| 亚洲人成77777在线视频| cao死你这个sao货| 亚洲激情在线av| 99久久久亚洲精品蜜臀av| 欧美精品啪啪一区二区三区| 日日夜夜操网爽| 成人三级做爰电影| 日本一区二区免费在线视频| 亚洲aⅴ乱码一区二区在线播放 | 精品少妇一区二区三区视频日本电影| 黑人操中国人逼视频| 久久国产精品人妻蜜桃| 人人澡人人妻人| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕色久视频| 亚洲免费av在线视频| 国产精品秋霞免费鲁丝片| 久9热在线精品视频| a在线观看视频网站| 久久亚洲真实| 日日干狠狠操夜夜爽| 亚洲狠狠婷婷综合久久图片| 丰满的人妻完整版| 黑人巨大精品欧美一区二区蜜桃| 手机成人av网站| 国产成+人综合+亚洲专区| 久久人妻av系列| 亚洲精品国产色婷婷电影| 国产亚洲欧美98| 97超级碰碰碰精品色视频在线观看| 视频区图区小说| 无人区码免费观看不卡| 亚洲aⅴ乱码一区二区在线播放 | 深夜精品福利| 亚洲熟女毛片儿| 丝袜美足系列| 国产熟女xx| 中文字幕精品免费在线观看视频| av免费在线观看网站| 亚洲人成网站在线播放欧美日韩| 50天的宝宝边吃奶边哭怎么回事| 99久久精品国产亚洲精品| 亚洲色图 男人天堂 中文字幕| 欧美在线一区亚洲| 1024香蕉在线观看| 神马国产精品三级电影在线观看 | 国产亚洲欧美精品永久| 欧美日本中文国产一区发布| 精品国产一区二区三区四区第35| 日本欧美视频一区| 免费av毛片视频| 在线天堂中文资源库| 女性生殖器流出的白浆| 在线天堂中文资源库| 国产精品免费一区二区三区在线| 人人妻,人人澡人人爽秒播| 女同久久另类99精品国产91| xxxhd国产人妻xxx| 免费不卡黄色视频| 国产精品偷伦视频观看了| 久久人妻福利社区极品人妻图片| 欧美黑人精品巨大| 久久精品aⅴ一区二区三区四区| 日韩欧美三级三区| 女人高潮潮喷娇喘18禁视频| 亚洲午夜精品一区,二区,三区| 男男h啪啪无遮挡| 国产精品电影一区二区三区| 久久草成人影院| svipshipincom国产片| 怎么达到女性高潮| 一级毛片高清免费大全| 国产成人啪精品午夜网站| 国产97色在线日韩免费| 国产色视频综合| 一个人观看的视频www高清免费观看 | 女人被狂操c到高潮| 免费高清在线观看日韩| 亚洲视频免费观看视频| 一a级毛片在线观看| 婷婷精品国产亚洲av在线| 夜夜看夜夜爽夜夜摸 | 国产精品98久久久久久宅男小说| 欧美日韩瑟瑟在线播放| 欧美日韩瑟瑟在线播放| 老司机在亚洲福利影院| 在线十欧美十亚洲十日本专区| netflix在线观看网站| 真人做人爱边吃奶动态| 91精品国产国语对白视频| av天堂久久9| 夜夜夜夜夜久久久久| 一进一出抽搐动态| x7x7x7水蜜桃| 老司机午夜十八禁免费视频| 黄色片一级片一级黄色片| 男人的好看免费观看在线视频 | 男女下面进入的视频免费午夜 | 91av网站免费观看| 免费观看人在逋| av超薄肉色丝袜交足视频| 午夜精品久久久久久毛片777| 精品欧美一区二区三区在线| 黄色怎么调成土黄色| 精品一区二区三区av网在线观看| 人妻久久中文字幕网| 黄色片一级片一级黄色片| 男人的好看免费观看在线视频 | 亚洲av日韩精品久久久久久密| 久久精品人人爽人人爽视色| 日韩欧美免费精品| 精品免费久久久久久久清纯| 亚洲国产欧美一区二区综合| 色精品久久人妻99蜜桃| 亚洲aⅴ乱码一区二区在线播放 | 一区二区日韩欧美中文字幕| 免费少妇av软件| 日韩成人在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 欧美乱色亚洲激情| 亚洲色图综合在线观看| 国产成人免费无遮挡视频| 日韩免费高清中文字幕av| 亚洲av成人一区二区三| 久久精品国产综合久久久| 国产成人精品久久二区二区免费| 久久影院123| 99精品在免费线老司机午夜| 亚洲国产欧美网| 国产伦一二天堂av在线观看| av电影中文网址| 91在线观看av| 欧美日韩精品网址| 妹子高潮喷水视频| 亚洲国产毛片av蜜桃av| 亚洲精品国产精品久久久不卡| 久久精品影院6| 女性生殖器流出的白浆| 香蕉丝袜av| 久久久久久久精品吃奶| 免费看十八禁软件| 亚洲人成电影观看| 两性夫妻黄色片| 欧美日韩一级在线毛片| 亚洲人成伊人成综合网2020| 亚洲精品国产一区二区精华液| svipshipincom国产片| 看免费av毛片| 国产99久久九九免费精品| 亚洲一区中文字幕在线| 久久精品亚洲精品国产色婷小说| 午夜91福利影院| 亚洲欧美日韩另类电影网站| 老汉色∧v一级毛片| 午夜视频精品福利| 国产野战对白在线观看| 日韩欧美一区二区三区在线观看| 国产精品自产拍在线观看55亚洲| 亚洲国产欧美一区二区综合| 国产亚洲欧美精品永久| 丁香六月欧美| 午夜91福利影院| 欧美乱色亚洲激情| 亚洲成人免费电影在线观看| 无遮挡黄片免费观看| 91精品三级在线观看| 国产伦人伦偷精品视频| 免费人成视频x8x8入口观看| 国产精品免费一区二区三区在线| 国产精品爽爽va在线观看网站 | 热99re8久久精品国产| 国产亚洲欧美98| 亚洲精品av麻豆狂野| 国产伦一二天堂av在线观看| 在线国产一区二区在线| 岛国在线观看网站| 老司机亚洲免费影院| 久热爱精品视频在线9| 一区二区三区激情视频| 18禁观看日本| 久久国产精品影院| 动漫黄色视频在线观看| 亚洲成人久久性| 日韩欧美国产一区二区入口| 亚洲国产欧美日韩在线播放| 一级a爱视频在线免费观看| 9191精品国产免费久久| 国产一区二区三区综合在线观看| 麻豆久久精品国产亚洲av | 亚洲欧美激情在线| xxxhd国产人妻xxx| 桃色一区二区三区在线观看| 亚洲熟妇熟女久久| 在线观看66精品国产| 变态另类成人亚洲欧美熟女 | 日韩大尺度精品在线看网址 | 在线观看免费视频网站a站| 婷婷精品国产亚洲av在线| ponron亚洲| 涩涩av久久男人的天堂| 色在线成人网| 国产蜜桃级精品一区二区三区| 久久久水蜜桃国产精品网| 亚洲狠狠婷婷综合久久图片| 黄色视频,在线免费观看| 男女之事视频高清在线观看| 最近最新中文字幕大全免费视频| 叶爱在线成人免费视频播放| 免费女性裸体啪啪无遮挡网站| 久久99一区二区三区| 电影成人av| 一进一出抽搐动态| 国产单亲对白刺激| 欧美性长视频在线观看| 91成人精品电影| 日韩欧美一区视频在线观看| 国产一区二区三区综合在线观看| 男女之事视频高清在线观看| 日本vs欧美在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 黄色a级毛片大全视频| 一本综合久久免费| 国产精品一区二区精品视频观看| 精品国产超薄肉色丝袜足j| 午夜免费观看网址| 亚洲av美国av| 亚洲成人免费av在线播放| 国产xxxxx性猛交| 精品久久久久久,| 欧美成人午夜精品| 动漫黄色视频在线观看| 超碰97精品在线观看| 国产视频一区二区在线看| 精品久久久精品久久久| 精品高清国产在线一区| 一级毛片高清免费大全| 水蜜桃什么品种好| 变态另类成人亚洲欧美熟女 | 日韩精品免费视频一区二区三区| 国产三级在线视频| av中文乱码字幕在线| 99久久久亚洲精品蜜臀av| 不卡一级毛片| 91字幕亚洲| 韩国精品一区二区三区| 91成人精品电影| 新久久久久国产一级毛片| 成人永久免费在线观看视频| 国产熟女xx| 国产乱人伦免费视频| 女人高潮潮喷娇喘18禁视频| 午夜影院日韩av| 国产精品香港三级国产av潘金莲| 99re在线观看精品视频| 亚洲av成人av| 欧美激情极品国产一区二区三区| 欧美激情极品国产一区二区三区| 国产精品香港三级国产av潘金莲| 亚洲精品成人av观看孕妇| 国产成人av激情在线播放| 久久久精品国产亚洲av高清涩受| 日本wwww免费看| 日本一区二区免费在线视频| 久久久精品国产亚洲av高清涩受| 三级毛片av免费| 十分钟在线观看高清视频www| 色播在线永久视频| 亚洲国产欧美网| 国产精品乱码一区二三区的特点 | 琪琪午夜伦伦电影理论片6080| 日韩大码丰满熟妇| 国产熟女午夜一区二区三区| 真人一进一出gif抽搐免费| 久久久久久亚洲精品国产蜜桃av| 亚洲成国产人片在线观看| 成人亚洲精品av一区二区 | svipshipincom国产片| 91大片在线观看| 99精品在免费线老司机午夜| 一级片'在线观看视频| 搡老乐熟女国产| 久久久久久久久免费视频了| 国产精品电影一区二区三区| av国产精品久久久久影院| 国产亚洲精品综合一区在线观看 | 欧美在线黄色| 亚洲成人国产一区在线观看| 国产又色又爽无遮挡免费看| 中文字幕人妻丝袜一区二区| 国产三级黄色录像| 亚洲熟妇中文字幕五十中出 | 日韩三级视频一区二区三区| 亚洲欧美一区二区三区黑人| 国产欧美日韩一区二区三| 日韩三级视频一区二区三区| 亚洲七黄色美女视频| 中国美女看黄片| 女警被强在线播放| 成人亚洲精品av一区二区 | √禁漫天堂资源中文www| 两个人看的免费小视频| 久久青草综合色| 日韩免费高清中文字幕av| 黑人巨大精品欧美一区二区mp4| 欧美+亚洲+日韩+国产| 国产高清视频在线播放一区| 久久精品亚洲熟妇少妇任你| 两个人看的免费小视频| 脱女人内裤的视频| 久久精品亚洲精品国产色婷小说| 老熟妇乱子伦视频在线观看| 9热在线视频观看99| 免费在线观看影片大全网站| 99久久人妻综合| 国产蜜桃级精品一区二区三区| 女性生殖器流出的白浆| 五月开心婷婷网| a级毛片在线看网站| 高清av免费在线| 韩国av一区二区三区四区| 久久香蕉精品热| 黄片小视频在线播放| 国内毛片毛片毛片毛片毛片| 国产精华一区二区三区| 欧美一级毛片孕妇| 一级毛片高清免费大全| 亚洲精品在线观看二区| 美女高潮喷水抽搐中文字幕| 精品免费久久久久久久清纯| 日韩欧美免费精品| 久久久久久久午夜电影 | 国产野战对白在线观看| 国产av在哪里看| 国产精品免费一区二区三区在线| 久久中文字幕一级| 大香蕉久久成人网| 热99re8久久精品国产| 老司机午夜福利在线观看视频| 精品日产1卡2卡| 亚洲第一欧美日韩一区二区三区| 亚洲国产中文字幕在线视频| 夜夜躁狠狠躁天天躁| 久久 成人 亚洲| 香蕉久久夜色| 一二三四在线观看免费中文在| 一二三四在线观看免费中文在| 亚洲午夜精品一区,二区,三区| 亚洲人成网站在线播放欧美日韩| 女人被躁到高潮嗷嗷叫费观| 99国产极品粉嫩在线观看| 神马国产精品三级电影在线观看 | 丰满的人妻完整版| 久久狼人影院| 精品国产国语对白av| 高清av免费在线| 怎么达到女性高潮| 亚洲自偷自拍图片 自拍| 亚洲成人免费电影在线观看| 婷婷精品国产亚洲av在线| 久热这里只有精品99| 国产不卡一卡二| 国产午夜精品久久久久久| 狠狠狠狠99中文字幕| 一区在线观看完整版| svipshipincom国产片| 国产亚洲精品综合一区在线观看 | 水蜜桃什么品种好| 无遮挡黄片免费观看| 在线观看一区二区三区| 国产蜜桃级精品一区二区三区| 成人亚洲精品一区在线观看| 成人亚洲精品一区在线观看| 丝袜人妻中文字幕| 国产国语露脸激情在线看| 两人在一起打扑克的视频| 久久久久精品国产欧美久久久| 午夜日韩欧美国产| 国产99久久九九免费精品| 国产乱人伦免费视频| 亚洲专区中文字幕在线| 亚洲激情在线av| 大陆偷拍与自拍| 国产精品成人在线| 老司机福利观看| 一区二区三区国产精品乱码| 欧美日韩av久久| 免费在线观看影片大全网站| 午夜a级毛片| 国产精品自产拍在线观看55亚洲| 麻豆国产av国片精品| 欧美精品一区二区免费开放| 50天的宝宝边吃奶边哭怎么回事| 精品国产乱子伦一区二区三区| 色综合婷婷激情| 精品电影一区二区在线| 国产xxxxx性猛交| 老汉色av国产亚洲站长工具| 亚洲欧美一区二区三区黑人| 国产精华一区二区三区| www国产在线视频色| 久久青草综合色| 欧美日韩av久久| 亚洲自拍偷在线| 国产亚洲欧美在线一区二区| av网站免费在线观看视频| 在线观看日韩欧美| 精品一区二区三区四区五区乱码| 成人影院久久| 91九色精品人成在线观看| 久久久久久人人人人人| 久久亚洲精品不卡| 满18在线观看网站| 久久精品国产99精品国产亚洲性色 | 黑人巨大精品欧美一区二区蜜桃| 宅男免费午夜| 久久亚洲真实| 久久久国产成人精品二区 | 亚洲va日本ⅴa欧美va伊人久久| 精品人妻1区二区| 亚洲中文av在线| 成人手机av| 久久影院123| 激情视频va一区二区三区| 午夜精品国产一区二区电影| 母亲3免费完整高清在线观看| www国产在线视频色| 又黄又粗又硬又大视频| 久久久久亚洲av毛片大全| 午夜福利在线免费观看网站| 成人18禁在线播放| 一级黄色大片毛片| 国产一卡二卡三卡精品| 99riav亚洲国产免费| 久久影院123| 亚洲中文字幕日韩| 91av网站免费观看| 亚洲avbb在线观看| 久久国产精品影院| av欧美777| 日韩欧美免费精品| 国产人伦9x9x在线观看| 亚洲 欧美 日韩 在线 免费| 巨乳人妻的诱惑在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美国产一区二区入口| 久久精品亚洲av国产电影网| 午夜久久久在线观看| 久久狼人影院| 在线看a的网站| 久久久国产一区二区| 国产深夜福利视频在线观看| 久久久久久久午夜电影 | 久久青草综合色| 12—13女人毛片做爰片一| 亚洲熟妇熟女久久| 99久久国产精品久久久| 成人18禁在线播放| 在线观看www视频免费| 精品乱码久久久久久99久播| www.www免费av| 不卡av一区二区三区| 国产又爽黄色视频| 满18在线观看网站| 午夜成年电影在线免费观看| 女人高潮潮喷娇喘18禁视频| 国产成人精品无人区| 久久久国产一区二区| 日韩精品中文字幕看吧| 午夜免费成人在线视频| 在线观看免费日韩欧美大片| 精品日产1卡2卡| 亚洲熟女毛片儿| 欧美日本中文国产一区发布| 日本撒尿小便嘘嘘汇集6| 女同久久另类99精品国产91| 国产精品国产av在线观看| 亚洲av熟女| 久久精品91无色码中文字幕| 国产成人啪精品午夜网站| 中文亚洲av片在线观看爽| 男女做爰动态图高潮gif福利片 | 亚洲熟妇熟女久久| 欧美人与性动交α欧美精品济南到| 真人做人爱边吃奶动态| 一级毛片高清免费大全| 亚洲伊人色综图| 日韩欧美一区二区三区在线观看| 一级,二级,三级黄色视频| 久久影院123| 少妇粗大呻吟视频| 国产精品日韩av在线免费观看 | 国产片内射在线| 在线观看www视频免费| 精品久久久久久久毛片微露脸| av国产精品久久久久影院| 成人三级黄色视频| 亚洲欧美精品综合久久99| 三上悠亚av全集在线观看| 露出奶头的视频| 成人精品一区二区免费| 精品久久久久久久毛片微露脸| 欧美久久黑人一区二区| 99re在线观看精品视频| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 黄片大片在线免费观看| 欧美大码av| 亚洲国产精品sss在线观看 | 久久精品国产清高在天天线| 亚洲黑人精品在线| 超碰成人久久| 欧美黄色片欧美黄色片| x7x7x7水蜜桃| 成人精品一区二区免费| 国产欧美日韩精品亚洲av| 女人精品久久久久毛片| 国产亚洲av高清不卡| 99国产精品免费福利视频| 高清在线国产一区| 亚洲成人免费av在线播放| 精品免费久久久久久久清纯| 欧美大码av| 纯流量卡能插随身wifi吗| 69精品国产乱码久久久| 欧美人与性动交α欧美软件| 久久精品国产清高在天天线| 久9热在线精品视频| 咕卡用的链子| 精品人妻1区二区| 少妇裸体淫交视频免费看高清 | 新久久久久国产一级毛片| av视频免费观看在线观看| 最近最新中文字幕大全免费视频| 中文字幕最新亚洲高清| 亚洲人成网站在线播放欧美日韩| 成人亚洲精品一区在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲精品美女久久av网站| 国产单亲对白刺激| 国产99白浆流出| 久久人妻福利社区极品人妻图片| 多毛熟女@视频| 亚洲国产欧美日韩在线播放| 国产精品爽爽va在线观看网站 | av电影中文网址| 国产91精品成人一区二区三区| 十分钟在线观看高清视频www| 美女福利国产在线| 国产精品九九99| www.熟女人妻精品国产| 涩涩av久久男人的天堂| av网站免费在线观看视频| 窝窝影院91人妻| 精品一区二区三卡| 久久中文字幕一级| 在线十欧美十亚洲十日本专区| www.999成人在线观看| aaaaa片日本免费| 在线观看www视频免费| 国产精华一区二区三区| av电影中文网址|