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

    檳榔不同發(fā)育時期果實轉(zhuǎn)錄組特征分析

    2020-09-02 06:54:35押輝遠陳葉張巖松許啟泰
    熱帶作物學報 2020年7期
    關(guān)鍵詞:檳榔

    押輝遠 陳葉 張巖松 許啟泰

    摘 ?要:檳榔(Areca catechuL.)果實是四大南藥之一。檳榔果的研究主要集中在生理生化、生防菌、有效成分及藥理、加工和利用等方面,對檳榔果的發(fā)育及其次生物質(zhì)形成的分子機制尚不清楚。本研究對不同發(fā)育時期的檳榔果皮和果核進行轉(zhuǎn)錄組測序,鑒定檳榔果不同發(fā)育時期的關(guān)鍵基因,以探討果實發(fā)育相關(guān)基因的表達特征及次生物質(zhì)形成有關(guān)的基因調(diào)控。結(jié)果顯示,檳榔果皮中檢測到4491個差異基因,其中617個差異基因共參與了111條KEGG代謝通路,生物過程代謝類有82個通路,共257個差異基因被注釋,參與次生代謝途徑共有5個,共27個差異基因。檳榔果核中檢測到5443個差異基因,其中898個差異基因共參與了118條通路,466個差異基因被注釋在生物代謝類通路上,共涉及89條通路,參與次生代謝相關(guān)的基因有53個,參與次生代謝途徑共7條。進一步分析表明:隨著果實的發(fā)育,果皮中80%次級代謝通路差異相關(guān)基因呈下調(diào)表達趨勢;而果核中71.4%次級代謝通路差異相關(guān)基因呈上調(diào)表達趨勢。本研究結(jié)果在轉(zhuǎn)錄組水平揭示了檳榔果發(fā)育的生物學過程,發(fā)現(xiàn)了不同時期檳榔果皮和果核中次級代謝相關(guān)調(diào)控基因表達的變化規(guī)律,也為檳榔的遺傳育種研究奠定了基礎。

    關(guān)鍵詞:檳榔;果實發(fā)育時期;轉(zhuǎn)錄組中圖分類號:S31??????文獻標識碼:A

    Analysis of Transcriptome Characteristics of ArecaatDifferent Developmental Stages

    YA Huiyuan1, CHEN Ye1, ZHANG Yansong1, XU Qitai1,2

    1. School of Food and Medicine, Luoyang Normal University, Luoyang, Henan 471934,?China; 2. Hainan Green Areca Science &?Technology Development Co., Ltd., Dingan, Hainan 571200, China

    Abstract: ArecaAreca catechu L.)?is?one of the four primary medicinal plants in south China.?In the study, the high-throughput sequencing technology was used to sequence the transcriptome of the peel and kernel in different periods to find differentially expressed genes in different developmental stages. Among the peels, 4491 differential genes were divided, of which 617 differential genes were involved in 111 KEGG metabolic pathways,257 differential genes with 82 pathways in the biological process metabolism class, and a total of 27 genes with 5 genes involved in the secondary metabolic pathway. There were 5443 differential genes in the betel nut kernel, according to the KEGG pathway annotation results, 898 differential genes were involved in 118 pathways, 466 differential genes were annotated on biological metabolic pathways for 89 pathways, 53 genes involved in secondary metabolism with 7 secondary metabolic pathways involved. Further analysis showed that with the development of fruit, 80% of the secondary metabolic pathways in the pericarp showed a down-regulated expression of the genes, while 71.4% of the secondary metabolic pathways in the kernel showed up-regulated expression. The results of the study preliminarily revealed the overall characteristics of the transcriptome of different tissues and different developmental stages of areca nut, and found that the expression of secondary metabolism-related regulatory genes in betel nut and pit were observed in different stages, which was the medicinal development and secondary of betel nut.

    Keywords: Areca catechuL.;?fruit development period; transcriptome

    DOI: 10.3969/j.issn.1000-2561.2020.07.001

    檳榔(Areca catechuL.)是棕櫚科檳榔屬常綠喬木,廣泛分布于南亞和東南亞等國家[1]。在我國,檳榔主要種植于云南、海南和臺灣等地,種植面積居世界第五,產(chǎn)量居世界第三,是我國熱帶、亞熱帶地區(qū)僅次于橡膠的第二大產(chǎn)業(yè)[2-4]。檳榔果營養(yǎng)價值豐富,是一種天然的藥食兩用食品,被稱為四大南藥之首[5]

    目前已從檳榔中分離鑒定出生物堿、黃酮類、單寧、三萜類、類固醇、脂肪酸等多種化學成分。生物堿具有顯著的生物活性,有驅(qū)蟲、消積、利氣利水、利濕除疸等功效[6]。植物生物堿分為萜類吲哚生物堿、哌啶生物堿、芐基異喹啉生物堿、茛菪堿和煙堿、嘌呤生物堿等。已有研究表明,檳榔堿是檳榔中的主要生物堿,屬于哌啶類生物堿,其含量為0.3%~0.7%[7-8]。目前,國內(nèi)外對檳榔生物堿的研究主要集中在檳榔堿及其衍生物的分離鑒定與藥理、生理作用、免疫等方面[5]?,F(xiàn)代醫(yī)學研究表明,檳榔堿在促進胃腸道運動、抗腫瘤、鎮(zhèn)痛、消炎、抗氧化、保護心血管系統(tǒng)等方面均有作用[1]。

    迄今為止,檳榔堿及其多種衍生物已被分離和鑒定[9]。檳榔的生物堿含量與其果實成熟度相關(guān)[10],成熟果胚乳中的檳榔堿含量高于青果胚乳,而成熟果果皮的含量低于青果果皮[11]。不同發(fā)育時期檳榔花序的檳榔堿含量也存在顯著差異,其含量隨著檳榔花的成熟而不斷降低[12]。然而,不同的檳榔品種由于成分也有較大差異,所以其藥用成分含量差異也較大,這種差異性極大地影響了檳榔用于中藥的藥效穩(wěn)定性[13-15]。關(guān)于檳榔生物堿等次級代謝產(chǎn)物生物合成的遺傳控制和分子調(diào)控至今還不清楚,這種狀況極大地阻礙了檳榔新品種的選育和藥用成分的有效利用。

    轉(zhuǎn)錄組測序能夠從整體水平研究基因表達量以及基因結(jié)構(gòu),揭示特定生物學過程中的分子機理,目前已廣泛應用于基礎研究、臨床診斷、藥物研發(fā)和分子育種等領(lǐng)域。檳榔果是藥食兩用的組織部位,研究檳榔果果實發(fā)育的遺傳基礎、果實內(nèi)次生物質(zhì)合成的分子調(diào)控,將十分有助于檳榔分子育種的快速發(fā)展。本研究對海南檳榔的檳榔果皮和檳榔果核進行轉(zhuǎn)錄組測序,分析篩選不同時期果皮和果核的差異基因,探究檳榔果生長發(fā)育過程中生物堿等次級代謝產(chǎn)物的分子調(diào)控機理。

    1??材料與方法

    1.1材料

    本研究于海南省定安縣采集坐果30 d和坐果180 d的檳榔(圖1)作為實驗材料用于轉(zhuǎn)錄組測序,每組樣品設置3個實驗重復(表1)。

    1.2方法

    1.2.1 ?檳榔轉(zhuǎn)錄組RNA提取及轉(zhuǎn)錄組測序??材料總RNA提取方法參照天根(TIANGEN)公司的TRNzol Universal總RNA提取試劑使用說明書,分別提取檳榔果的果皮和果核的總RNA。用安捷倫2100系統(tǒng)檢測總RNA完整性,Qubit?RNA測定試劑盒純度質(zhì)量,將樣品送至百邁客生物公司在Illumina Hiseq平臺進行轉(zhuǎn)錄組的測序。

    1.2.2??檳榔轉(zhuǎn)錄組數(shù)據(jù)的功能注釋和代謝通路分析??使用BLAST[16]軟件將Unigene序列與NR[17]、Swiss-Prot[18]、Gene Ontology(GO)[19]、Clusters of Orthologous Groups(COG)[20]、euKaryotic Orthologous Groups(KOG)[21]、eggNOG?4.5[22]、Kyoto Encyclopedia of Genes and Genomes(KEGG)[23]數(shù)據(jù)庫比對,使用KOBAS?2.0軟件[24]得到Unigene在KEGG中的KEGG Orthology結(jié)果,預測Unigene的氨基酸序列后,使用HMMER[25]軟件與Protein family(Pfam)[26]數(shù)據(jù)庫比對,獲得Unigene的注釋信息。

    1.2.3??檳榔轉(zhuǎn)錄組表達差異分析??將皮爾遜相關(guān)系數(shù)r(pearson correlation coefficient)作為樣品間相關(guān)性的評估指標[27]。r2越接近1,說明2個樣品的相關(guān)性越強。采用DESeq[28]進行樣品組間的差異表達分析,獲得2個條件之間的差異表達基因集;在差異表達分析過程中采用了公認有效的Benjamini-Hochberg方法對原有假設檢驗得到的顯著性P值(P-value)進行校正,并最終采用校正后的P值,即False Discovery Rate(FDR)作為差異表達基因篩選的關(guān)鍵指標,以降低對大量基因的表達值進行獨立的統(tǒng)計假設檢驗帶來的假陽性。在篩選過程中,將FDR小于0.01且差異倍數(shù)Fold Change(FC)大于或等于2作為篩選標準。其中,F(xiàn)C表示兩樣品(組)間表達量的比值。

    1.2.4??差異基因富集分析與關(guān)鍵功能基因篩選??差異表達基因GO功能富集,差異表達基因COG分類,差異表達基因KEGG注釋,差異表達基因KEGG通路富集分析,差異表達基因蛋白互作網(wǎng)絡。

    2??結(jié)果與分析

    2.1轉(zhuǎn)錄組測序數(shù)據(jù)統(tǒng)計

    選擇坐果30?d和180?d的檳榔果果皮和果核各3次重復,共12個樣品,獲得104.78?Gb Clean Data,各樣品Clean Data均達到8.73?Gb,Q30堿基百分比在88.78%及以上(表2)。

    轉(zhuǎn)錄組數(shù)據(jù)集已保藏在NCBI SRA數(shù)據(jù)庫中,編號:PRJNA590547。

    組裝后共獲得259?401條轉(zhuǎn)錄本序列,平均長1779.75?bp,N50長2827?bp。獲得94?562條轉(zhuǎn)錄本序列組裝單基因序列,其中長度在1?kb以上的Unigene有31?381條,對轉(zhuǎn)錄本序列組裝單基因序列進行功能注釋,包括與NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam數(shù)據(jù)庫的比對,共獲得35?806條單基因序列的注釋結(jié)果(表3)。

    2.2單基因功能注釋

    通過BLAST軟件將Unigene序列與NR、Swiss-Prot、GO、COG、KOG、eggNOG?4.5、KEGG數(shù)據(jù)庫比對,共得到35 806個基因獲得注釋(表4)。有8664(24.20%)條基因在COG數(shù)據(jù)庫中獲得了注釋,有17?155(47.91%)條基因在GO數(shù)據(jù)庫中獲得了注釋,有10 510(29.35%)條基因在KEGG數(shù)據(jù)庫中獲得了注釋,有17 917(50.04%)條基因在KOG數(shù)據(jù)庫中獲得了注釋,有18?072(50.47%)條基因在Pfam數(shù)據(jù)庫中獲得了注釋,有17 224(48.10%)條基因在Swiss-Prot數(shù)據(jù)庫中獲得了注釋,有30 231(84.43%)條基因在eggNOG?4.5數(shù)據(jù)庫中獲得了注釋,有34 334(95.89%)條基因在NR數(shù)據(jù)庫中獲得了注釋。

    2.3檳榔果皮不同發(fā)育時期差異表達基因篩選和富集分析

    以AcF3為實驗組,AcF1為對照組,共篩選到4491個差異基因,上調(diào)差異基因1991個,下調(diào)差異基因2500個。1563個差異基因在GO功能得到注釋,共分成3大類,45個小類。生物學過程(biological process)主要涉及細胞過程、代謝過程、單一機體過程、應激反應和生物調(diào)節(jié)5種生物學過程;細胞組分(cellular components)主要涉及細胞核、細胞器和質(zhì)膜3種細胞組分;分子功能(molecular function)主要涉及蛋白結(jié)合和催化活性2種分子功能(表5)。

    一般認為,經(jīng)過校正的P值(P-FDR)≤0.05 時,此GO功能存在顯著富集情況。在檳榔果皮差異基因的GO注釋結(jié)果顯示,生物學過程中差異基因顯著有單一生物體過程、生物調(diào)節(jié)、細胞成分的組織或者生物合成、信號、多生物過程、免疫生物進程和生物粘附。細胞組分中細胞器、細胞器部分、大分子復合物、細胞外區(qū)和細胞連接的差異基因有顯著差異。分子功能的催化活性、結(jié)構(gòu)分子活性和轉(zhuǎn)運活性的差異基因數(shù)有顯著差異。

    將果皮差異基因與KEGG數(shù)據(jù)庫進行比對,共有617個差異基因注釋在111個通路上,包括138個上調(diào)表達基因,479個下調(diào)表達基因。KEGG生化代謝通路共分成5個生物過程,包括代謝(metabolism)、遺傳信息過程(genetic informa tion processing)、細胞學過程(cellular processes),環(huán)境信息過程(environmental information processing)和有機體系統(tǒng)(organismal systems)。

    代謝類通路共257個差異基因被注釋,參與了82個通路,主要涉及氨基酸代謝、碳水化合物代謝、能量代謝、脂質(zhì)代謝和其他次級代謝產(chǎn)物的生物合成等。其中,參與次生代謝相關(guān)的基因有27個,參與次生代謝途徑共有5個,苯丙烷生物合成途徑(phenylpropanoid biosynthesis),類黃酮生物合成途徑(flavonoid biosynthesis),芪類、二芳基庚烷和姜醇生物合成途徑(stilbenoid, diarylheptanoid and gingerol biosynthesis),異喹啉類生物堿生物合成途徑(isoquinoline alkaloid biosynthesis),萜類、哌啶和吡啶生物堿合成途徑(tropane, piperidine and pyridine alkaloid biosynthesis)。其中,上調(diào)差異表達只有苯丙烷生物合成途徑,共3個基因包含的酶有4-香豆酸-CoA連接酶、肉桂酰輔酶A還原酶。下調(diào)差異表達基因分布在5個次生代謝途徑中:苯丙烷生物合成途徑,類黃酮生物合成途徑,芪類、二芳基庚烷和姜醇生物合成途徑,異喹啉類生物堿生物合成途徑,萜類、哌啶和吡啶生物堿合成途徑。下調(diào)基因共10個。調(diào)控酶:反式肉桂酸4-單加氧酶、咖啡酰-CoA O-甲基轉(zhuǎn)移酶、雙功能天冬氨酸氨基轉(zhuǎn)移酶和谷氨酸/天冬氨酸-預苯酸氨基轉(zhuǎn)移酶、伯胺氧化酶;協(xié)同調(diào)控苯丙烷生物合成途徑和萜類、哌啶和吡啶生物堿合成途徑共17個基因,包含β-葡萄糖苷酶、肉桂醇脫氫酶、過氧化物酶、托品酮還原酶I(表6)。

    2.4檳榔果核不同發(fā)育時期差異表達基因篩選和富集分析

    以AcN3為實驗組,AcN1為對照組,共篩選到5443個差異基因,上調(diào)差異基因2958個,下調(diào)差異基因2485個。這些差異基因有2452個在GO數(shù)據(jù)庫得到注釋,共分成3大類,48個小類。生物學過程(biological process)主要涉及細胞過程、代謝過程、單一機體過程、應激反應和生物調(diào)節(jié)基因5種,其中,呈顯著差異的有單一生物體過程、對刺激的響應、生物調(diào)節(jié)、復制、生物學階段。細胞組分(cellular components)主要涉及細胞核、細胞器和質(zhì)膜3種,細胞器和大分子復合物的差異基因數(shù)有顯著差異。分子功能(molecular function)主要涉及蛋白結(jié)合和催化活性2種,有顯著差異的包含催化活性和轉(zhuǎn)運活性(表7)

    將果核差異基因與KEGG數(shù)據(jù)庫進行比對,共有898個差異基因注釋在118個通路上,包括595個上調(diào)表達基因,303個下調(diào)表達基因,其中466個差異基因被注釋在代謝類通路上,共涉及89條通路。

    參與次生代謝相關(guān)的基因53個,參與次生代謝途徑共7條,包括咖啡因生物合成途徑(caffeine metabolism),黃酮和黃酮醇生物合成途徑(flavone and flavonol biosynthesis),類黃酮生物合成途徑,異喹啉類生物堿生物合成途徑,苯丙類生物合成途徑,芪類、二芳基庚烷和吡啶生物合成途徑,萜類、哌啶和吡啶生物堿合成途徑。其中,參與上調(diào)表達的基因41個,1個咖啡因生物合成途徑基因,12個苯丙類生物合成途徑基因,15個類黃酮生物合成途徑基因,5個黃酮和黃酮醇生物合成途徑基因,4個異喹啉類生物堿生物合成途徑基因,3個芪類、二芳基庚烷和吡啶生物合成途徑基因,3個萜類、哌啶和吡啶生物堿合成途徑基因,上調(diào)表達涉及的酶有尿酸氧化酶、伯胺氧化酶、4-羥基苯丙酮酸雙加氧酶、3-羥基丁?;?CoA脫氫酶、苯丙氨酸氨裂解酶、反式肉桂酸4-單加氧酶、4-香豆酸-輔酶A連接酶、反式肉桂酸4-單加氧酶、反式肉桂酸4-單加氧酶、查爾酮合成酶、查爾酮異構(gòu)酶類、黃酮3-單加氧酶、類黃酮3,5-羥化酶、白細胞花青素還原酶、花青素還原酶、類黃酮3,5-羥化酶、類黃酮3-單加氧酶、酪氨酸氨基轉(zhuǎn)移酶、酪氨酸脫羧酶(表8)。

    參與下調(diào)表達的基因2個,一個是苯丙類生物合成途徑基因,另一個是萜類、哌啶和吡啶生物堿合成途徑基因。下調(diào)表達涉及的酶為組氨醇-磷酸氨基轉(zhuǎn)移酶。協(xié)同調(diào)控的基因共19個,涉及的酶有雙功能天冬氨酸氨基轉(zhuǎn)移酶和谷氨酸/天冬氨酸-預苯酸氨基轉(zhuǎn)移酶、過氧化物酶、天冬氨酸氨基轉(zhuǎn)移酶。

    3討論

    檳榔果是檳榔的藥用部位。本研究利用高通量測序技術(shù)對2個不同時期的檳榔果皮和果核的轉(zhuǎn)錄組進行測序,初步揭示了檳榔轉(zhuǎn)錄組的整體表達特征,獲得了不同時期檳榔果皮和果核中次級代謝相關(guān)調(diào)控基因表達的變化規(guī)律。本研究結(jié)果為檳榔的遺傳多樣性、分子遺傳育種、遺傳進化等研究提供大量的參考信息。

    目前,已有研究表明在檳榔果實或種子中分離和鑒定出59種化合物,涉及生物堿、黃酮類化合物、單寧、三萜類和類固醇、脂肪酸等化合物,并證實其生物堿為檳榔藥理活性主要成分[29-32]。本研究通過轉(zhuǎn)錄組測序技術(shù)對其次級代謝物相關(guān)基因進行了檢索分析,從果皮中檢索到次級代謝相關(guān)差異基因共33個基因,主要涉及苯丙烷生物合成途徑,類黃酮生物合成途徑,芪類、二芳基庚烷和姜醇生物合成途徑,異喹啉類生物堿生物合成途徑和萜類、哌啶和吡啶生物堿合成途徑。從果核中檢索到次級代謝相關(guān)差異基因共78個基因,主要涉及咖啡因生物合成途徑,苯丙烷生物合成途徑,類黃酮生物合成途徑,黃酮和黃酮醇生物合成途徑,芪類、二芳基庚烷和姜醇生物合成途徑,異喹啉類生物堿生物合成途徑和萜類、哌啶和吡啶生物堿合成途徑。在檳榔果皮和果核中分別找到了參與生物堿代謝相關(guān)的差異基因,果皮中參與萜類、哌啶和吡啶生物堿合成途徑的基因有4個(2個基因協(xié)同表達,2個基因呈下調(diào)表達);果核中參與萜類、哌啶和吡啶生物堿合成途徑的基因為5個(2個基因協(xié)同表達,2個基因呈上調(diào)表達,1個下調(diào)表達)。這些研究對檳榔生物堿的生成調(diào)控機理研究具有重要意義。

    有關(guān)檳榔轉(zhuǎn)錄組測序分析的研究目前僅有一則相關(guān)報道。Manimekalai等[33]對檳榔葉片進行了轉(zhuǎn)錄組分析,在葉片中發(fā)現(xiàn)了7種高度表達的類黃酮和萜類化合物生物合成的酶。本研究在果核中僅檢測到查爾酮異構(gòu)酶、類黃酮3-單加氧酶、白細胞花青素還原酶隨時間的變化呈上調(diào)趨勢。此外,Manimekalai等[33]在葉片中檢測到泛醌和其他萜類代謝途徑、黃酮類代謝途徑和苯丙氨酸途徑的部分基因呈現(xiàn)高表達量,而本研究中果皮和果核的差異基因所涉及的酶類并沒有找到相應的酶,這可能與所用研究材料不同相關(guān)。

    Wang等[10]在1997年就曾報道過檳榔的生物堿含量與其果實成熟度相關(guān),劉蕊[11]的研究結(jié)果顯示,成熟果胚乳中的檳榔堿含量高于青果胚乳,而成熟果果皮的含量低于青果果皮,這與本研究中萜類、哌啶和吡啶生物堿合成途徑差異基因表達的結(jié)果基本一致。根據(jù)表6和表8的統(tǒng)計結(jié)果顯示:隨著果實的發(fā)育,果皮中次級代謝調(diào)控相關(guān)的下調(diào)表達基因數(shù)明顯多于上調(diào)表達的基因,80%次級代謝通路差異相關(guān)基因呈下調(diào)表達趨勢;而果核中的次級代謝調(diào)控相關(guān)基因上調(diào)表達基因數(shù)多于下調(diào)表達的基因,71.4%次級代謝通路差異相關(guān)基因呈上調(diào)表達趨勢。上述比對結(jié)果表明,實驗從轉(zhuǎn)錄組水平很好地驗證了前人的研究結(jié)果,為檳榔生物堿合成的分子調(diào)控機理的下一步研究提供了數(shù)據(jù)支持。

    [20]?Tatusov R L, Galperin M Y, Natale D A,et al. The COG database: a tool for genome scale analysis of protein functions and evolution[J]. Nucleic Acids Research,?2000,?28(1):?33-36.

    [21]?Koonin E?V, Fedorova N?D, Jackson J?D,et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes[J]. Genome Biology, 2004, 5(2): R7.

    [22]?Jaime H C, Damian S, Kristoffer F,et al. eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequen ces[J]. Nucleic Acids Research, 2016, 44(Database issue): D286-D293.

    [23]?Kanehisa M, Goto S, Kawashima S,et al. The KEGG resource for deciphering the genome[J]. Nucleic Acids Research,?2004,?32(Database issue): D277-D280.

    [24]?Xie C, Mao X, Huang J,et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases[J]. Nucleic Acids Research, 2011, 39: W316- W322.

    [25]?Eddy S?R. Profile hidden Markov models[J]. Bioinformatics, 1998, 14(9): 755-763.

    [26]?Finn R?D, Bateman A, Clements J,et al. Pfam: the protein families database[J]. Nucleic Acids Research,?2014, 42(Database issue): D222-230.

    [27]?Anders S, Huber W. Differential expression analysis for sequence count data[J]. Genome Biology, 2010, 11(10): R106.

    [28]?Leng N, Dawson J A, Thomson J A,et al. EBSeq: An empirical bayes hierarchical model for inference in RNA-seq experiments[J]. Bioinformatics, 2013, 29(16): 2073.

    [29]?Holdsworth D?K, Jones R?A, Self R.?Volatile alkaloids fromAreca catechu[J].?Phytochemistry,?1998, 48(3):?581-582.

    [30]?張??興,?梅文莉,?曾艷波,?等. 檳榔果實的酚類化學成分與抗菌活性的初步研究[J].?熱帶亞熱帶植物學報,?2009,?17(1):?74-76.

    [31]?楊文強, 王紅程, 王文婧, 等. 檳榔化學成分研究[J]. 中藥材, 2012, 35(3): 400-403.

    [32]?何細新, 李亞軍, 胡小鵬,?等. 檳榔殼中三萜及蒽醌類成分的分離與結(jié)構(gòu)鑒定[J]. 中藥新藥與臨床藥理, 2010, 21(6): 634-636.

    [33]?Manimekalai R, Nair S, Naganeeswaran A,et al.?Transcriptome sequencing and de novo assembly in arecanut,Areca catechuL.?elucidates the secondary metabolite pathway genes[J].?Biotechnology Reports,?2018, 17:?63-69.

    猜你喜歡
    檳榔
    叫停檳榔廣告很有必要
    人民周刊(2021年18期)2021-10-22 07:15:26
    有利有害的檳榔
    5月檳榔市場監(jiān)測分析
    8月海南氣候動蕩 農(nóng)作物減產(chǎn)劇增
    6月檳榔市場監(jiān)測分析
    燙漂和烘干過程對檳榔中檳榔堿含量的影響
    檳榔之味《粵中見聞》:“以檳榔肉兼食,味厚而芳。”
    讀書文摘(2017年10期)2017-10-16 14:26:57
    檳榔醉,愛心頭
    椰城(2016年2期)2016-09-21 09:18:04
    夢回檳榔園
    民族音樂(2016年1期)2016-08-28 20:02:52
    蒙藥檳榔十三味丸的抗抑郁作用
    在线观看午夜福利视频| 国产白丝娇喘喷水9色精品| av在线亚洲专区| 18+在线观看网站| 国产久久久一区二区三区| 国产在视频线在精品| 国产精品av视频在线免费观看| 精品久久久久久成人av| 99在线视频只有这里精品首页| av在线天堂中文字幕| 男女做爰动态图高潮gif福利片| 午夜福利在线观看免费完整高清在 | 熟女人妻精品中文字幕| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 欧美最新免费一区二区三区| 国产高潮美女av| 精华霜和精华液先用哪个| 色哟哟·www| 六月丁香七月| 久99久视频精品免费| 久久精品国产亚洲网站| 人妻少妇偷人精品九色| 精品无人区乱码1区二区| 女人被狂操c到高潮| 91麻豆精品激情在线观看国产| 亚洲自拍偷在线| 国内揄拍国产精品人妻在线| 国产精品一及| 国产成人福利小说| 听说在线观看完整版免费高清| 女人十人毛片免费观看3o分钟| 性色avwww在线观看| 非洲黑人性xxxx精品又粗又长| 校园春色视频在线观看| 性欧美人与动物交配| 99热网站在线观看| 九九久久精品国产亚洲av麻豆| 精品国产三级普通话版| 亚洲av男天堂| 国产成人福利小说| 99国产精品一区二区蜜桃av| 亚洲欧美日韩卡通动漫| 夜夜爽天天搞| 日韩亚洲欧美综合| 成年免费大片在线观看| 日本在线视频免费播放| 国产av麻豆久久久久久久| 国内精品一区二区在线观看| 国产女主播在线喷水免费视频网站 | 久久国内精品自在自线图片| 一本精品99久久精品77| 一级黄色大片毛片| 男女做爰动态图高潮gif福利片| 天天躁日日操中文字幕| 成年av动漫网址| 老熟妇乱子伦视频在线观看| 国产精品av视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 如何舔出高潮| 伦理电影大哥的女人| 免费观看a级毛片全部| 免费观看人在逋| 国产精品一二三区在线看| 亚洲性久久影院| 男女啪啪激烈高潮av片| 久久久精品94久久精品| 人妻系列 视频| kizo精华| 村上凉子中文字幕在线| 亚洲av一区综合| 欧美3d第一页| 久久鲁丝午夜福利片| av国产免费在线观看| 国产高清有码在线观看视频| 人妻久久中文字幕网| www.av在线官网国产| 少妇被粗大猛烈的视频| 六月丁香七月| 国产亚洲av片在线观看秒播厂 | 99久久成人亚洲精品观看| 欧美不卡视频在线免费观看| 国产蜜桃级精品一区二区三区| 午夜福利在线观看吧| 久久精品国产鲁丝片午夜精品| 亚洲成人av在线免费| 色吧在线观看| 久久久久九九精品影院| 亚洲在线自拍视频| 天天躁夜夜躁狠狠久久av| 久久久久久国产a免费观看| 一级黄片播放器| 日韩在线高清观看一区二区三区| 少妇人妻一区二区三区视频| 91在线精品国自产拍蜜月| 欧美极品一区二区三区四区| 免费无遮挡裸体视频| 丰满的人妻完整版| 午夜福利视频1000在线观看| 精品国内亚洲2022精品成人| 日日啪夜夜撸| 亚洲av第一区精品v没综合| 中文字幕制服av| 成年版毛片免费区| 免费电影在线观看免费观看| 在线免费十八禁| 午夜免费男女啪啪视频观看| 村上凉子中文字幕在线| 免费搜索国产男女视频| 久久人妻av系列| 麻豆久久精品国产亚洲av| 国产一区二区在线av高清观看| 51国产日韩欧美| 亚洲国产日韩欧美精品在线观看| 婷婷色综合大香蕉| 亚洲精品久久国产高清桃花| 啦啦啦韩国在线观看视频| 91午夜精品亚洲一区二区三区| 中文字幕熟女人妻在线| 亚洲精品456在线播放app| www.av在线官网国产| 日日撸夜夜添| 成人性生交大片免费视频hd| 亚洲无线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 一本久久精品| 99久久精品一区二区三区| 2022亚洲国产成人精品| 99国产精品一区二区蜜桃av| 精品久久久久久久久久久久久| 成年女人看的毛片在线观看| 九九热线精品视视频播放| 色视频www国产| 九草在线视频观看| 久久久久久久久大av| 久久久色成人| 亚洲成av人片在线播放无| 哪个播放器可以免费观看大片| 一个人观看的视频www高清免费观看| 精品人妻一区二区三区麻豆| 国产高潮美女av| 深爱激情五月婷婷| 校园春色视频在线观看| 三级国产精品欧美在线观看| 狂野欧美激情性xxxx在线观看| 久久久久久大精品| 干丝袜人妻中文字幕| 大香蕉久久网| 色综合亚洲欧美另类图片| 亚洲av免费高清在线观看| 18+在线观看网站| 国产三级在线视频| 久久精品综合一区二区三区| 最近手机中文字幕大全| 桃色一区二区三区在线观看| 亚洲在线自拍视频| h日本视频在线播放| 欧美日韩在线观看h| or卡值多少钱| 国产片特级美女逼逼视频| 国产在视频线在精品| 婷婷亚洲欧美| 2021天堂中文幕一二区在线观| 哪里可以看免费的av片| 午夜老司机福利剧场| 亚洲精品乱码久久久v下载方式| 欧美+日韩+精品| 亚洲精品日韩在线中文字幕 | 天堂影院成人在线观看| 亚洲经典国产精华液单| 久久精品久久久久久噜噜老黄 | 欧美一级a爱片免费观看看| 久久久午夜欧美精品| 99久久无色码亚洲精品果冻| 精品久久久噜噜| 久久99蜜桃精品久久| 亚洲欧美日韩卡通动漫| 亚洲人成网站高清观看| 又粗又爽又猛毛片免费看| 69人妻影院| 老女人水多毛片| 久久久久久九九精品二区国产| 一级毛片久久久久久久久女| 中文字幕人妻熟人妻熟丝袜美| 身体一侧抽搐| 69人妻影院| 国产 一区 欧美 日韩| av在线蜜桃| 乱人视频在线观看| 国产乱人偷精品视频| 大型黄色视频在线免费观看| 欧美日韩乱码在线| 久久久国产成人免费| 成人永久免费在线观看视频| 亚洲av男天堂| 欧美色欧美亚洲另类二区| 国产白丝娇喘喷水9色精品| 天堂√8在线中文| 国产精品蜜桃在线观看 | 欧美激情在线99| 成人综合一区亚洲| 哪个播放器可以免费观看大片| 国产亚洲av嫩草精品影院| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 亚洲av成人av| 久久久国产成人精品二区| 欧美丝袜亚洲另类| 免费av不卡在线播放| 丰满乱子伦码专区| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 夜夜爽天天搞| 日韩在线高清观看一区二区三区| 久久精品国产清高在天天线| 美女cb高潮喷水在线观看| 精品久久久噜噜| 22中文网久久字幕| 亚洲精品国产成人久久av| 性欧美人与动物交配| 黄色配什么色好看| 久久久久久伊人网av| 亚洲国产精品成人久久小说 | 久久久久久大精品| 99久国产av精品| 精品不卡国产一区二区三区| 我的女老师完整版在线观看| 精品欧美国产一区二区三| 啦啦啦观看免费观看视频高清| 欧美一区二区精品小视频在线| 中文字幕熟女人妻在线| 又爽又黄无遮挡网站| 欧美+日韩+精品| 欧美bdsm另类| av卡一久久| 99久久人妻综合| 九九在线视频观看精品| 国产精品人妻久久久影院| 麻豆乱淫一区二区| 白带黄色成豆腐渣| 丰满的人妻完整版| 日本一二三区视频观看| 又黄又爽又刺激的免费视频.| 三级国产精品欧美在线观看| 男女边吃奶边做爰视频| 久久久久性生活片| 91狼人影院| 五月伊人婷婷丁香| 六月丁香七月| 精品久久久噜噜| 波多野结衣高清无吗| 69av精品久久久久久| 黄色欧美视频在线观看| 国内少妇人妻偷人精品xxx网站| 午夜亚洲福利在线播放| 看片在线看免费视频| 久久婷婷人人爽人人干人人爱| 非洲黑人性xxxx精品又粗又长| 中文亚洲av片在线观看爽| 超碰av人人做人人爽久久| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 美女 人体艺术 gogo| 国产老妇女一区| 最近中文字幕高清免费大全6| 国产精品野战在线观看| 久久久久久久久大av| 国产精品久久久久久久久免| 此物有八面人人有两片| 午夜精品在线福利| 国产黄片美女视频| 我要搜黄色片| 国产伦精品一区二区三区四那| 菩萨蛮人人尽说江南好唐韦庄 | 人妻久久中文字幕网| 欧美又色又爽又黄视频| 麻豆成人午夜福利视频| 插逼视频在线观看| 人人妻人人澡欧美一区二区| 欧美色视频一区免费| 国产黄片美女视频| 日韩成人av中文字幕在线观看| 联通29元200g的流量卡| 久久精品国产清高在天天线| 中文字幕免费在线视频6| 亚洲精品日韩在线中文字幕 | 中文欧美无线码| 成人性生交大片免费视频hd| 天天一区二区日本电影三级| 久久中文看片网| 国产亚洲精品久久久com| 国产在视频线在精品| 美女脱内裤让男人舔精品视频 | av女优亚洲男人天堂| 国产伦理片在线播放av一区 | 亚洲精品久久国产高清桃花| 啦啦啦韩国在线观看视频| 国产三级在线视频| 亚洲高清免费不卡视频| 欧美最黄视频在线播放免费| 日本五十路高清| 欧美激情在线99| 亚洲欧美日韩高清在线视频| 99国产极品粉嫩在线观看| 女人十人毛片免费观看3o分钟| 国产色爽女视频免费观看| 国产精品久久久久久亚洲av鲁大| 免费人成视频x8x8入口观看| 成人综合一区亚洲| 亚洲av.av天堂| 亚洲人成网站在线播| 久久99精品国语久久久| 在线免费观看不下载黄p国产| 精品一区二区三区人妻视频| 亚洲精品久久国产高清桃花| 亚洲经典国产精华液单| 深爱激情五月婷婷| 18+在线观看网站| 国产黄a三级三级三级人| 国产又黄又爽又无遮挡在线| 丰满的人妻完整版| 直男gayav资源| 精品久久久久久久久亚洲| 国产日本99.免费观看| 12—13女人毛片做爰片一| 美女cb高潮喷水在线观看| 夫妻性生交免费视频一级片| 国产精品久久久久久av不卡| 婷婷亚洲欧美| 女人十人毛片免费观看3o分钟| 熟妇人妻久久中文字幕3abv| 在线观看美女被高潮喷水网站| 亚洲av中文字字幕乱码综合| 久久久久久国产a免费观看| 大又大粗又爽又黄少妇毛片口| 美女黄网站色视频| 国产 一区 欧美 日韩| 日韩制服骚丝袜av| 精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 久久精品久久久久久久性| 此物有八面人人有两片| 非洲黑人性xxxx精品又粗又长| 日韩人妻高清精品专区| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| 禁无遮挡网站| 人人妻人人澡人人爽人人夜夜 | 欧美激情国产日韩精品一区| 91久久精品国产一区二区三区| 久久久a久久爽久久v久久| 哪个播放器可以免费观看大片| 美女 人体艺术 gogo| 亚洲一区二区三区色噜噜| 嫩草影院入口| 婷婷亚洲欧美| 狠狠狠狠99中文字幕| 国产精品不卡视频一区二区| АⅤ资源中文在线天堂| 亚洲人成网站高清观看| 免费不卡的大黄色大毛片视频在线观看 | 国产免费一级a男人的天堂| 午夜福利在线观看吧| 九草在线视频观看| 一夜夜www| 91精品一卡2卡3卡4卡| 少妇的逼好多水| 真实男女啪啪啪动态图| 亚洲综合色惰| 国产精品久久电影中文字幕| 国产精品福利在线免费观看| 三级毛片av免费| 12—13女人毛片做爰片一| 亚洲精品国产av成人精品| 大香蕉久久网| 一区二区三区四区激情视频 | 深爱激情五月婷婷| 亚洲精华国产精华液的使用体验 | 久久久欧美国产精品| 亚洲av成人av| 欧美一级a爱片免费观看看| av免费在线看不卡| 国产av一区在线观看免费| avwww免费| 干丝袜人妻中文字幕| 日日啪夜夜撸| 人妻少妇偷人精品九色| 成人毛片a级毛片在线播放| 久久久久九九精品影院| 熟女电影av网| 最近手机中文字幕大全| 插逼视频在线观看| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 亚洲七黄色美女视频| 欧美+日韩+精品| 18禁在线播放成人免费| 日韩在线高清观看一区二区三区| 日韩强制内射视频| 麻豆一二三区av精品| 麻豆国产97在线/欧美| 联通29元200g的流量卡| 少妇的逼好多水| 国产色婷婷99| 国产成人福利小说| 免费一级毛片在线播放高清视频| 国产国拍精品亚洲av在线观看| 丝袜喷水一区| 久久婷婷人人爽人人干人人爱| 久久精品夜色国产| 99视频精品全部免费 在线| 久久韩国三级中文字幕| 波多野结衣巨乳人妻| 亚洲美女视频黄频| 国产91av在线免费观看| 男人舔女人下体高潮全视频| 成年av动漫网址| 欧美bdsm另类| 亚洲不卡免费看| 国内揄拍国产精品人妻在线| 别揉我奶头 嗯啊视频| 亚洲欧美成人综合另类久久久 | 99热精品在线国产| 亚洲欧洲日产国产| 亚洲国产日韩欧美精品在线观看| 国产精品女同一区二区软件| 欧美日韩国产亚洲二区| 校园人妻丝袜中文字幕| 看免费成人av毛片| av视频在线观看入口| 亚洲国产精品合色在线| 国产精品久久视频播放| 亚洲中文字幕日韩| 免费观看的影片在线观看| 国产熟女欧美一区二区| 男人和女人高潮做爰伦理| 色哟哟哟哟哟哟| 日本欧美国产在线视频| 亚洲中文字幕一区二区三区有码在线看| 毛片女人毛片| 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 亚洲激情五月婷婷啪啪| 久久99热6这里只有精品| 国产精品福利在线免费观看| 晚上一个人看的免费电影| av免费观看日本| 晚上一个人看的免费电影| 国产高清激情床上av| 国产精品麻豆人妻色哟哟久久 | 天堂√8在线中文| 欧美高清成人免费视频www| 嫩草影院新地址| 国产精品一区二区在线观看99 | 中文精品一卡2卡3卡4更新| av专区在线播放| 一级毛片久久久久久久久女| av福利片在线观看| 国产老妇伦熟女老妇高清| 国产av在哪里看| 欧美另类亚洲清纯唯美| 国产精品一区二区三区四区久久| 精品久久国产蜜桃| 国产91av在线免费观看| 国产av在哪里看| 日韩 亚洲 欧美在线| 精品一区二区三区人妻视频| 欧美成人一区二区免费高清观看| 久久鲁丝午夜福利片| 久久综合国产亚洲精品| 国产精品人妻久久久久久| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 深夜精品福利| 99久久成人亚洲精品观看| 亚洲美女视频黄频| 成熟少妇高潮喷水视频| 亚洲真实伦在线观看| 99久国产av精品| 国产私拍福利视频在线观看| 在线国产一区二区在线| 99国产精品一区二区蜜桃av| 国产精品一区二区性色av| 成年免费大片在线观看| 亚洲av男天堂| 亚洲精华国产精华液的使用体验 | 午夜久久久久精精品| 国产一区二区在线av高清观看| 内射极品少妇av片p| 精品久久久久久久末码| 日韩一本色道免费dvd| 久久精品国产99精品国产亚洲性色| 男女下面进入的视频免费午夜| 最近最新中文字幕大全电影3| 亚洲精品粉嫩美女一区| www.av在线官网国产| 给我免费播放毛片高清在线观看| 神马国产精品三级电影在线观看| 欧美高清性xxxxhd video| 天堂√8在线中文| 国产精品99久久久久久久久| 国产精华一区二区三区| 日本与韩国留学比较| 亚洲图色成人| 成人特级av手机在线观看| 大香蕉久久网| 国产精品久久久久久精品电影小说 | 热99在线观看视频| 日韩一区二区视频免费看| av在线播放精品| 成人特级黄色片久久久久久久| 久久久久性生活片| 久久久久久伊人网av| 在线国产一区二区在线| 欧美xxxx黑人xx丫x性爽| 91午夜精品亚洲一区二区三区| 午夜精品国产一区二区电影 | 一本精品99久久精品77| 乱人视频在线观看| 插逼视频在线观看| 久久99热6这里只有精品| 99热网站在线观看| 国产精品不卡视频一区二区| 国产成人影院久久av| 国产av在哪里看| 天堂影院成人在线观看| 欧美性感艳星| 亚洲人成网站在线播放欧美日韩| 久久99热6这里只有精品| 色吧在线观看| 亚洲,欧美,日韩| 99riav亚洲国产免费| 国产av在哪里看| 乱码一卡2卡4卡精品| 99九九线精品视频在线观看视频| 在线天堂最新版资源| 夜夜夜夜夜久久久久| 久久久久久伊人网av| 国产亚洲精品av在线| 一本久久精品| 国产精品无大码| 国产成人精品一,二区 | 国产成人精品久久久久久| 久久精品国产99精品国产亚洲性色| 国产成年人精品一区二区| 一个人看视频在线观看www免费| 日本欧美国产在线视频| www日本黄色视频网| 欧美最黄视频在线播放免费| 免费人成视频x8x8入口观看| 又粗又爽又猛毛片免费看| 亚洲国产精品国产精品| 欧美性猛交╳xxx乱大交人| 成人特级黄色片久久久久久久| 亚洲一区高清亚洲精品| 亚洲熟妇中文字幕五十中出| 国产精品一区二区三区四区免费观看| 国产美女午夜福利| 国产精品久久久久久精品电影| 91精品一卡2卡3卡4卡| 最近手机中文字幕大全| 日本色播在线视频| 可以在线观看毛片的网站| 69av精品久久久久久| 久久精品人妻少妇| 亚洲天堂国产精品一区在线| 亚洲av一区综合| 我的老师免费观看完整版| 99热这里只有精品一区| 欧美一区二区国产精品久久精品| 成人亚洲欧美一区二区av| 日本黄色片子视频| 亚洲av成人av| 一个人看的www免费观看视频| 51国产日韩欧美| 我的女老师完整版在线观看| 国产免费男女视频| 亚洲欧美中文字幕日韩二区| 99久久精品热视频| 91精品一卡2卡3卡4卡| 给我免费播放毛片高清在线观看| 3wmmmm亚洲av在线观看| 国产午夜精品一二区理论片| 午夜免费激情av| 免费观看a级毛片全部| 国产精品野战在线观看| 日韩欧美精品免费久久| 亚洲精品日韩av片在线观看| 免费一级毛片在线播放高清视频| 在线观看午夜福利视频| 91av网一区二区| 天堂av国产一区二区熟女人妻| 性插视频无遮挡在线免费观看| 只有这里有精品99| 国产在线精品亚洲第一网站| 欧美在线一区亚洲| 国产三级中文精品| 久久人妻av系列| 成人午夜高清在线视频| 日韩人妻高清精品专区| 国产又黄又爽又无遮挡在线| 成人午夜高清在线视频| 亚洲精品久久久久久婷婷小说 | 麻豆成人av视频| 国产精品久久久久久精品电影小说 | 天堂√8在线中文| 免费av不卡在线播放| 国产单亲对白刺激| 九色成人免费人妻av| 在线观看免费视频日本深夜| 91精品国产九色| a级毛片a级免费在线| av.在线天堂| 女的被弄到高潮叫床怎么办| 噜噜噜噜噜久久久久久91|