999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

海南黑山羊與海×努雜交山羊腸道微生物組成的比較

2025-04-07 00:00:00唐嘉陽(yáng)陳韜羽李世元蒙勇陳巧玲陳思高宏巖杜麗王鳳陽(yáng)滿初日嘎
關(guān)鍵詞:分析

DOI:10.3969/j.issn.2095-1191.2025.01.026

摘要:【目的】揭示雜交山羊擁有良好生長(zhǎng)性狀的內(nèi)在機(jī)理,為海南黑山羊的品種選育及海×努雜交山羊的飼養(yǎng)管理提供參考依據(jù)。【方法】選取出生4 d內(nèi)的海南黑山羊(HN組)和海×努雜交山羊(HY組)各6只,在統(tǒng)一的飼養(yǎng)管理標(biāo)準(zhǔn)下飼養(yǎng)至120日齡左右,采集糞便樣品進(jìn)行16S rRNA測(cè)序,分析2組山羊的腸道微生物群落組成差異,并對(duì)腸道微生物與山羊體重及平均日增重進(jìn)行Pearson相關(guān)分析。【結(jié)果】2組山羊的腸道微生物群落Alpha多樣性指數(shù)無(wú)顯著差異(Pgt;0.05),但基于布雷柯蒂斯(Bray-Curtis)距離的主坐標(biāo)分析(PCoA)發(fā)現(xiàn)2組山羊的腸道微生物組成存在一定差異。在門分類水平上,海南黑山羊與海×努雜交山羊腸道微生物均以厚壁菌門、擬桿菌門、疣微菌門、彎曲菌門、螺旋菌門及變形菌門等為優(yōu)勢(shì)菌門;在屬分類水平上,2組山羊腸道微生物組成相似,存在多種能發(fā)酵膳食纖維的菌屬;Metastats分析發(fā)現(xiàn),海×努雜交山羊腸道中的Monoglobus相對(duì)豐度顯著高于海南黑山羊(Plt;0.05,下同)。海南黑山羊腸道微生物的生物標(biāo)志物有腸球菌科、腸球菌屬及Family_XIII_UCG_001等;海×努雜交山羊腸道微生物的生物標(biāo)志物為短促生乳桿菌屬和顫螺旋菌科未知菌屬。海×努雜交山羊腸道微生物中,厚壁菌門相對(duì)豐度與體重及平均日增重呈顯著正相關(guān);海南黑山羊腸道微生物中,另枝菌屬相對(duì)豐度與體重呈極顯著正相關(guān)(Plt;0.01),與平均日增重呈顯著正相關(guān);將2組山羊合并為一個(gè)群體進(jìn)行分析,發(fā)現(xiàn)厚壁菌門相對(duì)豐度與山羊平均日增重呈顯著正相關(guān),Monoglo-bus相對(duì)豐度與山羊體重呈顯著正相關(guān)。【結(jié)論】海南黑山羊與海×努雜交山羊的腸道微生物組成存在一定差異,其中Monoglobus在海×努雜交山羊中的相對(duì)豐度顯著高于海南黑山羊。此外,腸道微生物中的厚壁菌門、另枝菌屬和Monoglobus可能是調(diào)控山羊生長(zhǎng)性能的關(guān)鍵菌群,為育肥山羊靶向型飼料添加劑的選擇提供了新方向。

關(guān)鍵詞:海南黑山羊;海×努雜交山羊;腸道微生物;16S rRNA測(cè)序;體重增長(zhǎng)

中圖分類號(hào):S827文獻(xiàn)標(biāo)志碼:A文章編號(hào):2095-1191(2025)01-0295-10

Comparison of intestinal microbial composition of Hainan black goat and hybrid of Hainan black goat×Nubian goat

TANG Jia-yang,CHEN Tao-yu,LI Shi-yuan,MENG Yong,CHEN Qiao-ling,CHEN Si,GAO Hong-yan,DU Li,WANG Feng-yang,Manchuriga*

(College of Tropical Agriculture and Forestry,Hainan University/Key Laboratory of Tropical Animal Reproductionamp; Breeding and Epidemic Disease Research of Hainan,Haikou,Hainan 570228,China)

Abstract:【Objective】To reveal the internal mechanism of hybrid goats with good growth traits,which could provide reference basis for breedselection of Hainan black goat and feeding management of hybrid of Hainan black goat×Nubian goat.【Method】Six Hainan black goats(HN group)and six hybrids of Hainan black goat×Nubian goat(HY group)born within 4 d were selected and raised under uniform feeding management standards until about 120 days old.Fecal samples were collected for 16S rRNA sequencing to analyze the differences in composition of the intestinal microbial community between the 2 groups of goats,and Pearson correlation analysis was performed to assess relationships between intestinal microbe community and body weight and daily weight gain.【Result】There was no significant difference in the Alpha di-versity index of the intestinal microbe between the 2 groups of goats(Pgt;0.05).However,principal coordinate analysis(PCoA)based on Bray-Curtis distances revealed distinct microbial community structures between the 2 groups of goats.At phylum level,Hainan black goat and hybrid of Hainan black goat×Nubian goat exhibited dominant phyla including Firmicutes,Bacteroidota,Verrucomicrobiota,Campylobacterota,Spirochaetota and Proteobacteria.At the genus level,the intestinal microbial composition was similar in the 2 groups,with a variety of fiber-fermenting genera.Metastats analysis identified significantly higher relative abundance of Monoglobus in hybrid of Hainan black goat×Nubian goat compared to Hainan black goat(Plt;0.05,the same below).The biomarkers of intestinal microbes in Hainan black goat in-cluded Enterococcaceae,Enterococcus,F(xiàn)amily_XIII_UCG_001,etc.The biomarkers of intestinal microbes in hybrid of Hainan black goat×Nubian goat were Levilactobacillus and unclassified_Oscillospiraceae.The relative abundance of Fir-micutes among the intestinal microbes of hybrid of Hainan black goat×Nubian goat was significantly and positively corre-lated with body weight and daily weight gain.In the intestinal microbes of Hainan black goat,there was extremely signifi-cant positive correlation(Plt;0.01)between the relative abundance of Alistipes and body weight,as well as a significant positive correlation with daily weight average gain.Combining 2 groups of goats into a population for analysis,it wasfound that the relative abundance of Firmicutes was significantly and positively correlated with average daily weight gain of goats,and the relative abundance of Monolobus was significantly and positively correlated with body weight of goats.【Conclusion】There are certain differences in the intestinal microbial composition between Hainan black goat and hybrid of Hainan black goat and Nubian goat,among which Monolobus has significantly higher relative abundance in Hainan black goat compared to hybrid of Hainan black goat×Nubian goat.In addition,F(xiàn)irmicutes,Alistipes and Monolobus inintestinal microbes might be key microbial communities that may regulate goat growth performance,which provides new directions for selection of targeted feed additives for goat fattening.

Key words:Hainan black goat;hybrid of Hainan black goat×Nubian goat;intestinal microbes;16S rRNA sequen-cing;body weight gain

Foundation items:China Agriculture Research System(CARS-38);Hainan Academician Innovation Platform Pro-ject(YSPTZX202153)

0引言

【研究意義】海南黑山羊是我國(guó)海南省重要的地方畜牧品種資源,具有耐粗飼、耐高溫高濕、抗病力強(qiáng)、性成熟早及肉質(zhì)鮮美等優(yōu)勢(shì)(葉玉秀等,2016),但也存在體重偏小、生長(zhǎng)速度緩慢等缺點(diǎn)(孫瑞萍等,2015)。努比亞山羊?qū)儆谌槿饧嬗眯蜕窖颍蚱浼∪庳S滿、體格較大、抗病力好、繁殖力強(qiáng),已廣泛應(yīng)用于山羊雜交育種(Stemmer et al.,2009;胡艷等,2021),以海南黑山羊?yàn)槟副尽⑴葋喩窖驗(yàn)楦副镜碾s交后代——海×努雜交山羊具有生長(zhǎng)迅速、體型較大、產(chǎn)肉量高等特點(diǎn)。腸道微生物是指動(dòng)物胃腸道內(nèi)數(shù)量龐大的微生物群(李世杰等,2024;徐蘭夢(mèng)等,2024),雜交引起的腸道微生物區(qū)系改變可通過(guò)調(diào)節(jié)宿主代謝方式,而對(duì)宿主的生長(zhǎng)性狀產(chǎn)生影響(Ber-gamaschietal.,2020)。因此,從腸道微生物層面入手分析雜交后海南黑山羊腸道微生物組成的變化,有利于揭示海×努雜交山羊擁有良好生長(zhǎng)性狀的內(nèi)在機(jī)理,為海南黑山羊的品種選育及飼養(yǎng)管理提供參考依據(jù)。【前人研究進(jìn)展】動(dòng)物腸道微生物主要由細(xì)菌、真菌、古細(xì)菌及病毒等組成(Mabwi et al.,2021)。據(jù)報(bào)道,人體內(nèi)擁有數(shù)萬(wàn)億的腸道微生物(Sender et al.,2016),且這些腸道微生物在機(jī)體代謝及免疫調(diào)節(jié)等方面發(fā)揮著至關(guān)重要的作用(Zhang et al.,2019;Aron-Wisnewsky et al.,2021)。除了飲食和生活環(huán)境是影響腸道微生物群落結(jié)構(gòu)的重要因素外(Singhetal.,2017;Trinhetal.,2018),遺傳因素對(duì)于腸道微生物組成也有影響。Doms等(2022)通過(guò)研究小鼠基因組與腸道黏膜微生物組成,并繪制高分辨率遺傳圖譜,結(jié)果發(fā)現(xiàn)存在一些具有高度遺傳性的腸道微生物;He等(2023)研究證實(shí),與屯昌豬相比,巴克夏豬與屯昌豬雜交子代的腸道微生物組成已發(fā)生變化,其中,乳酸菌屬(Lactobacillus)相對(duì)豐度的降低可能會(huì)導(dǎo)致雜交子代腸道免疫能力變?nèi)酰绽孜质暇鷮伲≒revotella)相對(duì)豐度的增加賦予雜交子代更好的營(yíng)養(yǎng)代謝和吸收功能。但也有研究表明,宿主的基因差異對(duì)腸道微生物組成幾乎沒(méi)有影響。Jiang等(2020)研究發(fā)現(xiàn),在相同的環(huán)境和飼養(yǎng)條件下,海南黑山羊與薩嫩山羊雜交F1代中黑山羊和白山羊的腸道微生物組成差異不明顯;楊蓮等(2021)研究證實(shí),純種柯樂(lè)豬和巴×柯雜交豬的腸道微生物組成相似。此外,Qi等(2021)研究發(fā)現(xiàn),無(wú)菌仔豬的體重和去脂體重較正常仔豬約降低40%,且伴隨著肌肉發(fā)育不良;Cheng等(2022)將高體重綿羊和低體重綿羊的腸道微生物分別移植到小鼠體內(nèi),結(jié)果表明,接受高體重綿羊腸道微生物移植的小鼠擁有更高的體重及日增重,可能與腸道中Prevotellaceae_UCG-001和Lachnospiraceae_NK4A1 36_group的相對(duì)豐度更高有關(guān);同時(shí)有諸多研究證實(shí),豬腸道微生物中的消化球菌屬(Peptococcus)、真桿菌屬(Eubacterium)、月單胞菌屬(Selenomo-nas)、莫拉氏菌屬(Moraxella)、柯林斯菌屬(Collin-sella)、Ruminococcaceae_UCG_014和Ruminococca-ceae_UCG_008等與其體重增長(zhǎng)呈正相關(guān)(Ding et al.,2019;Oh etal.,2020;姜長(zhǎng)津等,2022)。可見(jiàn),腸道微生物組成與宿主的生長(zhǎng)性狀密切相關(guān)。【本研究切入點(diǎn)】至今,針對(duì)山羊雜交后代腸道微生物組成變化的相關(guān)研究鮮見(jiàn)報(bào)道,腸道微生物組成變化對(duì)山羊生長(zhǎng)性狀的影響機(jī)制也尚未明確。【擬解決的關(guān)鍵問(wèn)題】采集海南黑山羊和海×努雜交山羊的糞便樣本,通過(guò)16S rRNA測(cè)序分析雜交后山羊腸道微生物組成的變化,并探究腸道微生物與山羊體重及日增重的相關(guān)性,旨在揭示雜交山羊擁有良好生長(zhǎng)性狀的內(nèi)在機(jī)理,為海南黑山羊的品種選育及海×努雜交山羊的飼養(yǎng)管理提供參考依據(jù)。

1材料與方法

1.1試驗(yàn)動(dòng)物及飼養(yǎng)管理

選取出生4 d內(nèi)的海南黑山羊和海×努雜交山羊各6只(4只公羊和2只母羊),海南黑山羊設(shè)為HN組,海×努雜交山羊設(shè)為HY組(表1)。所有山羊打耳標(biāo)標(biāo)記,采用常規(guī)免疫程序進(jìn)行免疫,并按相同的標(biāo)準(zhǔn)統(tǒng)一飼養(yǎng)管理,每天8:00、14:00和20:00投喂3次,自由飲水。每天觀察山羊健康情況。

1.2體重及日增重測(cè)定

記錄好每只山羊的出生重,至山羊約120日齡后禁食12 h,測(cè)量其空腹體重,然后計(jì)算平均日增重(表1)。

1.3糞便樣本收集

在120日齡左右,收集12只山羊的糞便樣本,裝入做好標(biāo)記的15 mL離心管中,立即投入液氮中保存,運(yùn)回實(shí)驗(yàn)室后轉(zhuǎn)移至-80℃冰箱保存。HN組山羊及其糞便樣本編號(hào)分別為HN1、HN2、HN3、HN4、HN5和HN6,HY組山羊及其糞便樣本編號(hào)分別為HY1、HY2、HY3、HY4、HY5和HY6。

1.4腸道微生物16S rRNA測(cè)序及數(shù)據(jù)處理

采用TGuide S96磁珠法土壤/糞便基因組DNA提取試劑盒提取山羊糞便DNA,根據(jù)全長(zhǎng)引物序列合成攜帶Barcode的特異引物27F(5'-AGRGTTTGA TYNTGGCTCAG-3')和1492R(5'-TASGGHTACCTT GTTASGACTT-3')并進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系30.0μL:DNA模板1.5μL,KOD OneTM PCR Master Mix 15.0μL,上、下游引物各1.5μL,Nuclease-free Water 10.5μL。擴(kuò)增程序:95℃預(yù)變性2 min;98℃10 s,55℃30 s,72℃1 min 30 s,進(jìn)行25個(gè)循環(huán);72℃延伸2 min。PCR擴(kuò)增產(chǎn)物經(jīng)Qubit定量分析及1.8%瓊脂糖凝膠電泳檢測(cè)合格后混樣,以AMPure PB磁珠進(jìn)行純化回收,并通過(guò)Agilent 2100生物分析儀和Illumina文庫(kù)定量試劑盒進(jìn)行評(píng)估,上機(jī)測(cè)序委托北京百邁客生物科技有限公司完成。

下機(jī)數(shù)據(jù)以CCS文件呈現(xiàn),通過(guò)Lima v1.7.0根據(jù)Barcode對(duì)CCS文件進(jìn)行識(shí)別,得到Raw-CCS數(shù)據(jù)并進(jìn)行質(zhì)控過(guò)濾:使用Cutadapt 1.9.1去除數(shù)據(jù)中的引物(Martin,2011),得到Clean-CCS數(shù)據(jù);通過(guò)UCHIME v4.2鑒定并去除嵌合體序列(Edgar et al.,2011),得到Effective-CCS數(shù)據(jù);然后以USEARCH處理Reads(Edgar,2013),在97.0%的相似水平下進(jìn)行聚類分析以獲得OTUs。以SILVA為參考數(shù)據(jù)庫(kù),通過(guò)樸素貝葉斯分類器結(jié)合比對(duì)的方法對(duì)OTUs特征序列進(jìn)行分類學(xué)注釋,得到每個(gè)特征序列對(duì)應(yīng)的物種分類信息,分別在門和屬分類水平上統(tǒng)計(jì)各樣本的微生物群落組成。

1.5統(tǒng)計(jì)分析

通過(guò)QIIME 2 2020.6對(duì)山羊腸道微生物群落進(jìn)行Alpha多樣性分析(Bolyen et al.,2019),使用主坐標(biāo)分析(PCoA)進(jìn)行Beta多樣性分析,并利用LEfSe分析篩選出不同組別間具有統(tǒng)計(jì)學(xué)差異的生物標(biāo)志物(Segataetal.,2011)。此外,以SPSS 22.0對(duì)腸道微生物與山羊體重及日增重進(jìn)行Pearson相關(guān)分析,通過(guò)Metastats在不同分類水平上進(jìn)行組間差異分析。

2結(jié)果與分析

2.1山羊腸道微生物16S rRNA測(cè)序數(shù)據(jù)質(zhì)控結(jié)果

對(duì)12份山羊糞便樣本微生物進(jìn)行16S rRNA測(cè)序,共獲得723879條Raw-CCS,經(jīng)質(zhì)控過(guò)濾,去除嵌合體后得到599689條Effective-CCS,各樣本的Raw-CCS比率均在80.0%以上。其中,HN組得到301589條Effective-CCS,HY組得到298100條Effective-CCS。從圖1可看出,當(dāng)采樣序列數(shù)達(dá)10000時(shí),稀釋曲線逐漸趨于平緩。在97.0%的相似水平下對(duì)所有Effective-CCS進(jìn)行OTU聚類分析,結(jié)果(圖2)顯示,共獲得2385個(gè)OTUs,其中,共有OTUs為1841個(gè),HN組特有OTUs為265個(gè),HY組特有OTUs為279個(gè)。

2.2山羊腸道微生物群落Alpha多樣性分析結(jié)果

山羊腸道微生物群落Alpha多樣性分析結(jié)果(表2)顯示,HY組的ACE指數(shù)、Chao1指數(shù)、Simp-son指數(shù)和Shannon指數(shù)均高于HN組,但差異不顯著(rgt;0.05,下同)。此外,所有山羊糞便樣本的測(cè)序結(jié)果覆蓋率均在99.00%以上,即測(cè)序結(jié)果可代表海南黑山羊和海×努雜交山羊腸道微生物的真實(shí)情況。

2.3山羊腸道微生物群落Beta多樣性分析結(jié)果

基于布雷柯蒂斯(Bray-Curtis)距離矩陣的PCoA分析結(jié)果(圖3)顯示,第一主坐標(biāo)(PCoA1)解釋了32.70%的差異性,第二主坐標(biāo)(PCoA2)解釋了13.62%的差異性。2組山羊的大部分樣本坐標(biāo)呈聚集趨勢(shì),但有小部分樣本偏離主群分布,說(shuō)明HN組和HY組山羊的腸道微生物組成存在潛在差異。

2.4山羊腸道微生物群落組成分析結(jié)果

2.4.1門分類水平上的山羊腸道微生物組成在2組山羊糞便樣本中共檢測(cè)到25種菌門,相對(duì)豐度排前10的菌門如圖4-A所示。其中,厚壁菌門(Fir-micutes)、擬桿菌門(Bacteroidota)、疣微菌門(Verru-comicrobiota)、彎曲菌門(Campylobacterota)、螺旋菌門(Spirochaetota)和變形菌門(Proteobacteria)為六大優(yōu)勢(shì)菌門,其相對(duì)豐度總和均占HN組和HY組山羊腸道微生物總量的99.00%以上。經(jīng)Metastats分析發(fā)現(xiàn),在HN組與HY組間各菌門的相對(duì)豐度差異均不顯著。

2.4.2屬分類水平上的山羊腸道微生物組成在2組山羊糞便樣本中共檢測(cè)到339種菌屬,相對(duì)豐度排前25的菌屬如圖4-B所示。山羊腸道微生物群落中相對(duì)豐度排名前10的優(yōu)勢(shì)菌屬分別是毛螺菌科未知菌屬(unclassified_Lachnospiraceae)、毛螺菌科NK4A136菌群(Lachnospiraceae_NK4A136_group)、UCG_005、uncultured_rumen_bacterium、理研菌科RC9腸道菌群(Rikenellaceae_RC9_gut_group)、擬桿菌屬(Bacteroides)、毛螺菌科AC2044菌群(Lachno-spiraceae_AC2044_group)、克里斯滕森菌科R_7菌群(Christensenellaceae R 7 group)、瘤胃球菌屬(Ru-minococcus)及unclassified_[Eubacterium]_coprosta-noligenes_group,占HN組山羊腸道微生物總量的58.67%以上,占HY組山羊腸道微生物總量的54.80%以上。經(jīng)Metastats分析發(fā)現(xiàn),Monoglobus在HY組的相對(duì)豐度顯著高于HN組(Plt;0.05,下同)。

2.5山羊腸道微生物組間差異分析結(jié)果

LEfSe分析結(jié)果(圖5和圖6)表明,當(dāng)LDA分值≥2時(shí),HN組山羊腸道微生物的生物標(biāo)志物有腸球菌科(Enterococcaceae)、腸球菌屬(Enterococcus)、Family_XIII_UCG_001、Clostridium_methylpentosum_group和unclassified_Clostridium_methylpentosum_group;HY組山羊腸道微生物的生物標(biāo)志物為短促生乳桿菌屬(Levilactobacillus)和顫螺旋菌科未知菌屬(unclassified_Oscillospiraceae)。

2.6腸道微生物與山羊生長(zhǎng)速度的相關(guān)分析結(jié)果

選取山羊腸道微生物群落中相對(duì)豐度排前10的菌門和排名前30的菌屬,分別與山羊體重及平均日增重進(jìn)行Pearson相關(guān)分析,結(jié)果(圖7)表明:在HN組,另枝菌屬(Alistipes)相對(duì)豐度與山羊的體重呈極顯著正相關(guān)(Plt;0.01),與平均日增重呈顯著正相關(guān)(圖7-B)。在HY組,厚壁菌門相對(duì)豐度與山羊的體重及平均日增重呈顯著正相關(guān)(圖7-C)。將2組山羊群體合并為一個(gè)群體進(jìn)行分析,結(jié)果顯示,厚壁菌門相對(duì)豐度與山羊的平均日增重呈顯著正相關(guān)(圖7-E),Monoglobus相對(duì)豐度與山羊的體重呈顯著正相關(guān)(圖7-F)。

3討論

腸道微生物在保障宿主健康及其免疫力方面發(fā)揮著至關(guān)重要的作用。腸道微生物通過(guò)利用蛋白質(zhì)和碳水化合物進(jìn)行代謝,得到短鏈脂肪酸,而短鏈脂肪酸不僅能供給機(jī)體能量,還可通過(guò)結(jié)合G蛋白偶聯(lián)受體及促進(jìn)組蛋白乙酰化等方式,調(diào)控機(jī)體健康水平,具有抗炎、維持腸道穩(wěn)態(tài)和預(yù)防疾病的功能作用(Simpson and Campbell,2015;Fu et al.,2019;Zhao et al.,2019;van der Hee and Wells,2021)。此外,腸道微生物區(qū)系受飲食、遺傳及外界環(huán)境等因素的影響,而腸道微生物區(qū)系的改變又能進(jìn)一步影響宿主的表型性狀。山羊作為草食動(dòng)物,其飲食結(jié)構(gòu)中膳食纖維占比較高,導(dǎo)致山羊腸道內(nèi)存在大量消化膳食纖維的菌群,以便于將膳食纖維代謝為短鏈脂肪酸,而有利于機(jī)體吸收。本研究中,2組山羊雖然采用統(tǒng)一的飼養(yǎng)管理標(biāo)準(zhǔn),較好地限制了飲食及環(huán)境等因素對(duì)機(jī)體腸道微生物組成的影響,但其腸道微生物組成仍存在一定差異,由此推測(cè)是努比亞山羊基因的引入致使雜交山羊腸道微生物組成較海南黑山羊有所改變,且這些腸道微生物組成變化或許對(duì)雜交山羊優(yōu)良性狀的形成有促進(jìn)作用。

為了從腸道微生物層面揭示雜交山羊擁有良好生長(zhǎng)性狀的內(nèi)在機(jī)理,本研究分別采集海南黑山羊和海×努雜交山羊的糞便樣本進(jìn)行16S rRNA測(cè)序分析,結(jié)果顯示,在門分類水平上,2組山羊腸道微生物均以厚壁菌門和擬桿菌門為絕對(duì)優(yōu)勢(shì)菌門,與Cao等(2023)的研究結(jié)果一致。厚壁菌門和擬桿菌門作為哺乳動(dòng)物腸道中最具優(yōu)勢(shì)的菌門(Ley et al.,2008),均可分解利用膳食纖維,而增強(qiáng)宿主對(duì)營(yíng)養(yǎng)物質(zhì)的吸收利用(McKee et al.,2021;Sun et al.,2023)。在屬分類水平上,2組山羊腸道的優(yōu)勢(shì)菌屬包括毛螺菌科未知菌屬、毛螺菌科NK4A136菌群、UCG_005及理研菌科RC9腸道菌群等,與Zhi等(2022)的研究結(jié)果相似。毛螺菌科各菌屬和理研菌科RC9腸道菌群均可發(fā)酵膳食纖維及其他糖類,而產(chǎn)生大量短鏈脂肪酸,進(jìn)而調(diào)節(jié)宿主生理功能并提供能量,供宿主生長(zhǎng)所需及日常活動(dòng)消耗(Ahmad etal.,2022;Zaplana et al.,2023)。Monoglo-bus于2017年從人類糞便中首次分離獲得,是一種具有果膠及膳食纖維代謝功能的腸道微生物,可能與機(jī)體的健康相關(guān)(Kimetal.,2017;Li etal.,2024)。本研究發(fā)現(xiàn),Monoglobus在海×努雜交山羊腸道中的相對(duì)豐度顯著高于海南黑山羊。在海南黑山羊腸道生物標(biāo)志物中,腸球菌科及其腸球菌屬包含多種致病性細(xì)菌(Fiore et al.,2019);Clostridium_methyl-pentosum_group能以戊糖或甲基戊糖為底物進(jìn)行代謝,而產(chǎn)生短鏈脂肪酸(Himelbloom and Canale-Parola,1989)。在海×努雜交山羊腸道生物標(biāo)志物中,短促生乳桿菌屬是一種常見(jiàn)的益生菌,通常存在于母乳中,可產(chǎn)生大量γ-氨基丁酸,具有增強(qiáng)免疫的作用(Wu and Shah,2017;Duraisamyet al.,2022),故推測(cè)海×努雜交山羊腸道免疫水平高于海南黑山羊。γ-氨基丁酸是一種抑制性神經(jīng)遞質(zhì),當(dāng)其與受體結(jié)合后可引起鉀離子外流及氯離子內(nèi)流,而導(dǎo)致突觸后膜超極化,降低神經(jīng)元的興奮性。γ-氨基丁酸還具有擴(kuò)張血管、抗氧化、抗抑郁、促進(jìn)生長(zhǎng)激素分泌及抗腫瘤細(xì)胞增殖等功能(Sarasa etal.,2020)。顫螺旋菌科未知菌屬則在膳食纖維的降解過(guò)程中發(fā)揮關(guān)鍵作用,常被認(rèn)為是丁酸鹽的生產(chǎn)者(Rang et al.,2023),可為海×努雜交山羊提供了更好的膳食纖維消化能力。

為進(jìn)一步探究山羊腸道微生物與生長(zhǎng)性狀之間的相關(guān)性,本研究還分別對(duì)2組山羊腸道微生物與體重及平均日增重進(jìn)行Pearson相關(guān)分析,結(jié)果發(fā)現(xiàn),海×努雜交山羊腸道微生物中的厚壁菌門與體重及平均日增重呈顯著正相關(guān),可能是厚壁菌門給海×努雜交山羊帶來(lái)更好的膳食纖維代謝功能所致,膳食纖維是山羊飼料最重要的組成部分,較高的厚壁菌門相對(duì)豐度有助于山羊?qū)︼暳系南铡4送猓癖诰T與擬桿菌門比值較高的腸道微生物結(jié)構(gòu)能更好地幫助宿主利用能量,與宿主的肥胖息息相關(guān)(DiBaise etal.,2008)。在海南黑山羊腸道微生物中,另枝菌屬與體重呈極顯著正相關(guān),與平均日增重呈顯著正相關(guān)。另枝菌屬是一種革蘭氏陰性厭氧菌,其主要代謝產(chǎn)物為乙酸鹽和丙酸鹽,通常與健康的腸道微生物代謝狀態(tài)相關(guān)(Parker et al.,2020;Xu et al.,2022)。另枝菌屬相對(duì)豐度與體重增長(zhǎng)呈正相關(guān),可能是由于海南黑山羊正處于生長(zhǎng)發(fā)育期,擁有健康穩(wěn)定的代謝狀態(tài)是確保發(fā)育期山羊體重增長(zhǎng)的關(guān)鍵。本研究將2組山羊合并為一個(gè)群體進(jìn)行分析,結(jié)果發(fā)現(xiàn),厚壁菌門相對(duì)豐度與山羊的平均日增重呈顯著正相關(guān),Monoglobus相對(duì)豐度與山羊的體重呈顯著正相關(guān),故推測(cè)Monoglobus的膳食纖維及果膠代謝能力對(duì)山羊生長(zhǎng)發(fā)育有促進(jìn)作用,但其代謝功能及對(duì)宿主機(jī)體的影響仍需進(jìn)一步驗(yàn)證。

4結(jié)論

海南黑山羊與海×努雜交山羊的腸道微生物組成存在一定差異,其中Monoglobus在海×努雜交山羊中的相對(duì)豐度顯著高于海南黑山羊。此外,腸道微生物中的厚壁菌門、另枝菌屬和Monoglobus可能是調(diào)控山羊生長(zhǎng)性能的關(guān)鍵菌群,為育肥山羊靶向型飼料添加劑的選擇提供了新方向。

參考文獻(xiàn)(References):

胡艷,張叁保,張瑜,高小童,江雨航,申玉建,韋英明,蔣欽楊.2021.努比亞山羊KISS1基因多態(tài)性與繁殖性狀的關(guān)聯(lián)分析[J].南方農(nóng)業(yè)學(xué)報(bào),52(10):2880-2886.[Hu Y,Zhang S B,Zhang Y,Gao X T,Jiang Y H,Shen Y J,Wei Y M,Jiang Q Y.2021.Association analysis on KISS1 gene polymorphism with reproductive traits in Nubian goats[J].Journal of Southern Agriculture,52(10):2880-2886.]doi:10.3969/j.issn.2095-1191.2021.10.029.

姜長(zhǎng)津,陳云,潘鵬丞,陳寶劍,關(guān)志惠,盧慧林,謝炳坤,覃兆鮮.2022.育肥豬糞便菌群結(jié)構(gòu)及其與生長(zhǎng)速度的相關(guān)性分析[J].中國(guó)畜牧獸醫(yī),49(3):924-931.[Jiang C J,Chen Y,Pan P C,Chen B J,Guan Z H,Lu H L,Xie B K,Qin Z X.2022.Fecal microflora structure and its correla-tion with growth rate in finishing pigs[J].China Animal Husbandryamp;Veterinary Medicine,49(3):924-931.]doi:10.16431/j.cnki.1671-7236.2022.03.014.

李世杰,張金輝,王英民,李嫚嫚,趙航,王榮軍,李俊強(qiáng),李曉迎,許會(huì)艷,菅復(fù)春.2024.不同性別湖羊早期定植腸道菌群分析[J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào),58(2):249-258.[Li S J,Zhang J H,Wang Y M,Li M M,Zhao H,Wang R J,Li J Q,Li X Y,Xu H Y,Jian F C.2024.Analysis of intestinal flora early colonizing in Hu-sheep of different genders[J].Journal of Henan Agricultural University,58(2):249-258.]doi:10.16445/j.cnki.1000-2340.20231120.001.

孫瑞萍,魏立民,劉圈煒,鄭心力,宋奇,王峰.2015.復(fù)合酶制劑對(duì)育肥前期海南黑山羊生產(chǎn)性能和血清生化指標(biāo)的影響[J].安徽農(nóng)業(yè)科學(xué),43(19):102-103.[Sun R P,Wei L M,Liu Q W,Zheng X L,Song Q,Wang F.2015.Effect of complex enzyme preparation on the growth perfor-mance and serum biochemical indicators in Hainan black goat in early fattening stage[J].Journal of Anhui Agricul-tural Sciences,43(19):102-103.]doi:10.13989/j.cnki.0517-6611.2015.19.040.

徐蘭夢(mèng),袁巖聰,彭昕,張城,張宏亮,方成偉,萬(wàn)堃,何航,張傳師,劉安芳,章杰.2024.復(fù)合菌制劑對(duì)麻羽肉雞性能和腸道微生物的影響[J].西南大學(xué)學(xué)報(bào)(自然科學(xué)版),46(5):87-98.[Xu L M,Yuan Y C,Peng X,Zhang C,Zhang H L,F(xiàn)ang C W,Wan K,He H,Zhang C S,Liu AF,Zhang J.2024.Effects of compound bacterial preparation on performance and gut microbiota of Mayu broilers[J].Journal of Southwest University(Natural Science Edition),46(5):87-98.]doi:10.13718/j.cnki.xdzk.2024.05.008.

楊蓮,燕志宏,黃維江,楊仕鈺,吳光松,顧麗菊,林鵬飛,楊蓉,李平,任麗群,張蕓.2021.純種柯樂(lè)豬與巴×柯雜交豬腸道菌群結(jié)構(gòu)的研究[J].動(dòng)物營(yíng)養(yǎng)學(xué)報(bào),33(3):1359-1371.[Yang L,Yan Z H,Huang W J,Yang S Y,Wu G S,Gu L J,Lin P F,Yang R,Li P,Ren L Q,Zhang Y.2021.Study on intestinal microflora structure of pure-bred Kele pigs and Berkshire×Kele hybrid pigs[J].Chinese Journal of Animal Nutrition,33(3):1359-1371.]doi:10.3969/j.issn.1006-267x.2021.03.019.

葉玉秀,雷湘蘭,程文科,晁哲.2016.海南黑山羊研究進(jìn)展與產(chǎn)業(yè)發(fā)展思路[J].熱帶農(nóng)業(yè)科學(xué),36(10):114-118.[Ye Y X,Lei X L,Cheng W K,Chao Z.2016.Research prog-ress and some thoughts of industrial development of Hainan black goats[J].Chinese Journal of Tropical Agri-culture,36(10):114-118.]doi:10.12008/j.issn.1009-2196.2016.10.024.

Ahmad A A,Zhang J B,Liang Z Y,Du M,Yang Y Y,Zheng J S,Yan P,Long R J,Tong B,Han J L,Ding X Z.2022.Age-dependent variations in rumen bacterial community of Mongolian cattle from weaning to adulthood[J].BMC Mi-crobiology,22:213.doi:10.1186/s 12866-022-02627-6.

Aron-Wisnewsky J,Warmbrunn M V,Nieuwdorp M,Clément K.2021.Metabolism and metabolic disorders and the microbiome:The intestinal microbiota associated with obesity,lipid metabolism,and metabolic health-patho-physiology and therapeutic strategies[J].Gastroenterology,160(2):573-599.doi:10.1053/j.gastro.2020.10.057.

Bergamaschi M,Maltecca C,Schillebeeckx C,McNulty N P,Schwab C,Shull C,F(xiàn)ix J,Tiezzi F.2020.Heritability and genome-wide association of swine gut microbiome fea-tureswith growth and fatness parameters[J].Scientific Re-ports,10:10134.doi:10.1038/s41598-020-66791-3.

Bolyen E,Rideout J R,Dillon M R,Bokulich NA,Abnet C C,Al-Ghalith G A,Alexander H,Alm E J,Arumugam M,Asnicar F,Bai Y,Bisanz J E,Bittinger K,BrejnrodA,Bri-slawn C J,Brown C T,Callahan B J,Caraballo-Rodríguez A M,Chase J,Cope E K,Da Silva R,Diener C,Dorrestein P C,Douglas G M,Durall D M,Duvallet C,Edwardson C F,Ernst M,Estaki M,F(xiàn)ouquier J,Gauglitz J M,Gibbons S M,Gibson D L,Gonzalez A,Gorlick K,Guo J R,Hill-mann B,Holmes S,Holste H,Huttenhower C,Huttley G A,Janssen S,Jarmusch A K,Jiang L J,Kaehler B D,Kang K B,Keefe C R,Keim P,Kelley S T,Knights D,Koester I,Kosciolek T,Kreps J,Langille M G I,Lee J,Ley R,Liu Y X,Loftfield E,Lozupone C,Maher M,Marotz C,Mar-tin B D,McDonald D,McIver L J,Melnik A V,Metcalf J L,Morgan S C,Morton J T,Naimey A T,Navas-Molina J A,Nothias L F,Orchanian S B,Pearson T,Peoples S L,Petras D,Preuss M L,Pruesse E,Rasmussen L B,Rivers A,Robeson M S,Rosenthal P,Segata N,Shaffer M,Shiffer A,Sinha R,Song S J,Spear J R,Swafford A D,Thompson L R,Torres P J,Trinh P,Tripathi A,Turnbaugh P J,Ul-Hasan S,van der Hooft J J J,Vargas F,Vázquez-Baeza Y,Vogtmann E,von Hippel M,Walters W,Wan Y H,Wang M X,Warren J,Weber K C,Williamson C H D,Willis A D,Xu Z Z,Zaneveld J R,Zhang Y L,Zhu Q Y,Knight R,Caporaso J G.2019.Reproducible,interactive,scalable and extensible microbiome data science using QIIME 2[J].Nature Biotechnology,37:852-857.doi:10.1038/s41587-019-0209-9.

Cao Y H,F(xiàn)eng T,Wu Y J,Xu Y X,Du L,Wang T,Luo Y H,Wang Y,Li Z P,Xuan Z Y,Chen S M,Yao N,Gao N L,Xiao Q,Huang K W,Wang X B,Cui K Q,Rehman S U,Tang X F,Liu D W,Han H B,Li Y,Chen W H,Liu Q Y.2023.The multi-kingdom microbiome of the goat gastroin-testinal tract[J].Microbiome,11(1):219.doi:10.1186/s40168-023-01651-6.

Cheng J B,Zhang X X,Zhang D Y,Zhang Y K,Li X L,Zhao Y,Xu D,Zhao L M,Li W X,Wang J H,Zhou B B,Lin C C,Yang X B,Zhai R,Cui P P,Zeng X W,Huang Y L,Ma Z W,Liu J,Wang W M.2022.Sheep fecal transplantation affects growth performance in mouse models by altering gut microbiota[J].Journal of Animal Science,100(11):skac303.doi:10.1093/jas/skac303.

DiBaise J K,Zhang H S,Crowell M D,Krajmalnik-Brown R,Anton Decker G,Rittmann B E.2008.Gut microbiota and its possible relationship with obesity[J].Mayo Clinic Pro-ceedings,83(4):460-469.doi:10.4065/83.4.460.

Ding X H,Lan W S,Liu G,Ni H J,Gu J D.2019.Exploring possible associations of the intestine bacterial microbiome with the pre-weaned weight gaining performance of piglets in intensive pig production[J].Scientific Reports,9:15534.doi:10.1038/s41598-019-52045-4.

Doms S,F(xiàn)okt H,Rühlemann M C,Chung C J,Kuenstner A,Ibrahim S M,F(xiàn)rankeA,Turner L M,Baines J F.2022.Keyfeatures of the genetic architecture and evolution of host-microbe interactions revealed by high-resolution genetic mapping of the mucosa-associated gut microbiome in hybrid mice[J].eLife,11:e75419.doi:10.7554/eLife.75419.

Duraisamy S,Husain F,Balakrishnan S,Sathyan A,Subramani P,Chidambaram P,Arokiyaraj S,Al-Qahtani W H,Raja-bathar J,Kumarasamy A.2022.Phenotypic assessment of probiotic and bacteriocinogenic efficacy of indigenous LAB strains from human breast milk[J].Current Issues in Molecular Biology,44(2):731-749.doi:10.3390/cimb440 20051.

Edgar R C,Haas B J,Clemente J C,Quince C,Knight R.2011.UCHIME improves sensitivity and speed of chimera detec-tion[J].Bioinformatics,27(16):2194-2200.doi:10.1093/bioinformatics/btr381.

Edgar R C.2013.UPARSE:Highly accurate OTU sequences from microbial amplicon reads[J].Nature Methods,10:996-998.doi:10.1038/nmeth.2604.

Fiore E,van Tyne D,Gilmore M S.2019.Pathogenicity of enterococci[J].Microbiology Spectrum,7(4):23.doi:10.1128/microbiolspec.gpp3-0053-2018.

Fu X D,Liu Z M,Zhu C L,Mou H J,Kong Q.2019.Nondi-gestible carbohydrates,butyrate,and butyrate-producing bacteria[J].Critical Reviews in Food Science and Nutri-tion,59(S1):S130-S152.doi:10.1080/10408398.2018.15 42587.

He J Y,Zhang Y C,Li H,Xie Y S,Huang G Q,Peng C,Zhao P J,Wang Z G.2023.Hybridization alters the gut microbial and metabolic profile concurrent with modifying intestinal functions in Tunchang pigs[J].Frontiers in Microbiology,14:1159653.doi:10.3389/fmicb.2023.1159653.

Himelbloom B H,Canale-Parola E.1989.Clostridium methyl-pentosum sp.nov.:A ring-shaped intestinal bacterium that ferments only methylpentoses and pentoses[J].Archives of Microbiology,151(4):287-293.doi:10.1007/bf00406 553.

Jiang S M,Huo D X,You Z K,Peng Q N,Ma C C,Chang H B,Lin X,Wang L,Zhang J H.2020.The distal intestinal microbiome of hybrids of Hainan black goats and Saanen goats[J].PLoS One,15(1):e0228496.doi:10.1371/jour-nal.pone.0228496.

Kim C C,Kelly W J,Patchett M L,Tannock G W,Jordens Z,Stoklosinski H M,Taylor J W,Sims I M,Bell T J,Rosen-dale D I.2017.Monoglobus pectinilyticus gen.nov.,sp.nov.,a pectinolytic bacterium isolated from human faeces[J].International Journal of Systematic and Evolutionary Microbiology,67(12):4992-4998.doi:10.1099/ijsem.0.0 02395.

Ley R E,Hamady M,Lozupone C,Turnbaugh P J,Ramey R R,Bircher J S,Schlegel M L,Tucker T A,Schrenzel M D,Knight R,Gordon J I.2008.Evolution of mammals and their gut microbes[J].Science,320(5883):1647-1651.doi:10.1126/science.1155725.

Li W,Guan S H,Hu X L,Zhao H,Cai J H,Li X H,Zhang X Y,Zhu W,Pan X,Li S X,Tian J K.2024.Lysimachiacapillipes Hemsl.saponins ameliorate colorectal cancer in mice via regulating gut microbiota and restoring metabolic profiles[J].Fitoterapia,175:105959.doi:10.1016/j.fitote.2024.105959.

Mabwi H A,Kim E,Song D G,Yoon H S,Pan C H,Komba E V G,Ko G P,Cha K H.2021.Synthetic gut microbiome:Advances and challenges[J].Computational and Structu-ral Biotechnology Journal,19:363-371.doi:10.1016/j.csbj.2020.12.029.

Martin M.2011.Cutadapt removes adapter sequences from high-throughput sequencing reads[J].EMBnet.Journal,17:10-12.doi:10.14806/EJ.17.1.200.

McKee L S,La Rosa S L,Westereng B,Eijsink V G,Pope P B,Larsbrink J.2021.Polysaccharide degradation by the Bac-teroidetes:Mechanisms and nomenclature[J].Environmen-tal Microbiology Reports,13(5):559-581.doi:10.1111/1758-2229.12980.

Oh J K,Chae J P,Pajarillo E A B,Kim S H,Kwak M J,Eun J S,Chee S W,Whang K Y,Kim S H,Kang D K.2020.Association between the body weight of growing pigs and the functional capacity of their gut microbiota[J].Animal Science Journal,91(1):e13418.doi:10.1111/asj.13418.

Parker B J,Wearsch P A,Veloo A C M,Rodriguez-Palacios A.2020.The genus Alistipes:Gut bacteria with emergingimplications to inflammation,cancer,and mental health[J].Frontiers in Immunology,11:906.doi:10.3389/fim-mu.2020.00906.

Qi R L,Sun J,Qiu X Y,Zhang Y,Wang J,Wang Q,Huang J X,Ge L P,Liu Z H.2021.The intestinal microbiota con-tributes to the growth and physiological state of muscle tis-sue in piglets[J].Scientific Reports,11:11237.doi:10.1038/s41598-021-90881-5.

Rang Y F,Liu H,Cheng X B,Li WY,Shi J,Ou G H,Huang H Y,Chen C Y,Xiao X M,Liu C H.2023.Structural charac-terization of pectic polysaccharides from Amaranth cauda-tus leaves and the promotion effect on hippocampal glucagon-like peptide-1 level[J].International Journal of Biological Macromolecules,242(4):124967.doi:10.1016/j.ijbiomac.2023.124967.

Sarasa S B,Mahendran R,Muthusamy G,Thankappan B,Selta D R F,Angayarkanni J.2020.A brief review on the non-protein amino acid,gamma-amino butyric acid(GABA):Its production and role in microbes[J].Current Microbio-logy,77(4):534-544.doi:10.1007/s00284-019-01839-w.

Segata N,Izard J,Waldron L,Gevers D,Miropolsky L,Garrett W S,Huttenhower C.2011.Metagenomic biomarker dis-covery and explanation[J].Genome Biology,12(6):R60.doi:10.1186/gb-2011-12-6-r60.

Sender R,F(xiàn)uchs S,Milo R.2016.Revised estimates for the number of human and bacteria cells in the body[J].PLoS Biology,14(8):e1002533.doi:10.1371/journal.pbio.1002533.

Simpson H L,Campbell B J.2015.Review article:Dietary fibre-microbiota interactions[J].Alimentary Pharmacologyamp;Therapeutics,42(2):158-179.doi:10.1111/apt.13248.

Singh R K,Chang H W,Yan D,Lee K M,Ucmak D,Wong K,Abrouk M,F(xiàn)arahnik B,Nakamura M,Zhu T H,Bhutani T,Liao W.2017.Influence of diet on the gut microbiome and implications for human health[J].Journal of Translational Medicine,15:73.doi:10.1186/s 12967-017-1175-y.

Stemmer A,Siegmund-Schultze M,Gall C,Z?rate A V.2009.Development and worldwide distribution of the Anglo Nubian goat[J].Tropical and Subtropical Agroecosys-tems,11:185-188.

Sun Y G,Zhang S S,Nie Q X,He H J,Tan H Z,Geng F,Ji H H,Hu J L,Nie S P.2023.Gut firmicutes:Relationship with dietary fiber and role in host homeostasis[J].Critical Reviews in Food Science and Nutrition,63(33):12073-12088.doi:10.1080/10408398.2022.2098249.

Trinh P,Zaneveld J R,Safranek S,Rabinowitz P M.2018.One health relationships between human,animal,and environ-mental microbiomes:A mini-review[J].Frontiers in Pub-lic Health,6:235.doi:10.3389/fpubh.2018.00235.

van der Hee B,Wells J M.2021.Microbial regulation of host physiology by short-chain fatty acids[J].Trends in Micro-biology,29(8):700-712.doi:10.1016/j.tim.2021.02.001.

Wu Q L,Shah N P.2017.Highγ-aminobutyric acid production from lactic acid bacteria:Emphasis on Lactobacillus brevis as a functional dairy starter[J].Critical Reviews in Food Science and Nutrition,57(17):3661-3672.doi:10.1080/10408398.2016.1147418.

Xu Z L,Jiang W,Huang W L,Lin Y,Chan F K L,Ng S C.2022.Gut microbiota in patients with obesity and meta-bolic disorders—A systematic review[J].Genesamp;Nutri-tion,17:2.doi:10.1186/s 12263-021-00703-6.

Zaplana T,Miele S,Tolonen A C.2023.Lachnospiraceae are emerging industrial biocatalysts and biotherapeutics[J].Frontiers in Bioengineering and Biotechnology,11:1324 396.doi:10.3389/fbioe.2023.1324396.

Zhang C X,Wang H Y,Chen T X.2019.Interactions betweenintestinal microflora/probiotics and the immune system[J].BioMed Research International,2019:6764919.doi:10.1155/2019/6764919.

Zhao J F,Zhang X Y,Liu H B,Brown M A,Qiao S Y.2019.Dietary protein and gut microbiota composition and func-tion[J].Current Proteinamp;Peptide Science,20(2):145-154.doi:10.2174/1389203719666180514145437.

Zhi W B,Tang K,Yang J S,Yang T S,Chen R,Huang J M,Tan H S,Zhao J G,Sheng Z W.2022.Research on the gut microbiota of Hainan black goat[J].Animals,12(22):3129.doi:10.3390/ani 12223129.

(責(zé)任編輯:蘭宗寶)

猜你喜歡
分析
禽大腸桿菌病的分析、診斷和防治
隱蔽失效適航要求符合性驗(yàn)證分析
電力系統(tǒng)不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
電力系統(tǒng)及其自動(dòng)化發(fā)展趨勢(shì)分析
經(jīng)濟(jì)危機(jī)下的均衡與非均衡分析
對(duì)計(jì)劃生育必要性以及其貫徹實(shí)施的分析
GB/T 7714-2015 與GB/T 7714-2005對(duì)比分析
出版與印刷(2016年3期)2016-02-02 01:20:11
中西醫(yī)結(jié)合治療抑郁癥100例分析
偽造有價(jià)證券罪立法比較分析
在線教育與MOOC的比較分析
主站蜘蛛池模板: 欧美高清三区| 欧美成人免费| 再看日本中文字幕在线观看| 色综合中文| 91久久国产热精品免费| 无码福利日韩神码福利片| аⅴ资源中文在线天堂| 尤物特级无码毛片免费| 日韩 欧美 国产 精品 综合| 国产精品第一区| 亚洲香蕉伊综合在人在线| 国产成人免费视频精品一区二区| 91www在线观看| 国产网站免费看| 国产女人18水真多毛片18精品| 99免费视频观看| 国产一在线| 熟妇丰满人妻| 精品福利一区二区免费视频| 欧美精品影院| 在线观看av永久| 国产在线精彩视频二区| 在线观看网站国产| 伊人成人在线视频| 国产精品入口麻豆| 欧美中文字幕一区二区三区| 青青草综合网| 丁香五月激情图片| 免费一级毛片| 久久综合九色综合97网| 国产高潮流白浆视频| 国产第二十一页| 日韩精品专区免费无码aⅴ| 东京热av无码电影一区二区| 中文字幕日韩欧美| 亚国产欧美在线人成| 99久久这里只精品麻豆| 3D动漫精品啪啪一区二区下载| 中国成人在线视频| 一级在线毛片| 亚洲天堂视频在线播放| a欧美在线| 欧美日韩国产在线人| 在线播放91| 亚洲中文久久精品无玛| 国产在线观看一区精品| 4虎影视国产在线观看精品| 亚洲 日韩 激情 无码 中出| 伊人福利视频| 蜜桃视频一区| 天堂网亚洲综合在线| 久草视频一区| 另类欧美日韩| 一区二区自拍| 国产女人18毛片水真多1| 久久久受www免费人成| a毛片基地免费大全| 久久99国产乱子伦精品免| 在线日韩一区二区| 日本黄网在线观看| 四虎AV麻豆| 国产91精选在线观看| 国产成人精品男人的天堂下载| 毛片视频网址| 国产精品xxx| 久久综合伊人77777| 欧美日韩北条麻妃一区二区| 日韩精品久久无码中文字幕色欲| 欧美日本在线一区二区三区| 国产呦精品一区二区三区网站| www.91在线播放| 亚州AV秘 一区二区三区| 亚洲综合狠狠| 在线播放真实国产乱子伦| 青青青伊人色综合久久| 韩国v欧美v亚洲v日本v| 五月天久久综合国产一区二区| 亚洲精品无码在线播放网站| 91精品在线视频观看| 国产特级毛片aaaaaa| 综合色在线| 国产二级毛片|