• <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
    蒙藥檳榔十三味丸的抗抑郁作用
    黑人猛操日本美女一级片| 国产亚洲av高清不卡| 色婷婷av一区二区三区视频| 99国产综合亚洲精品| 在线视频色国产色| 国产精品久久久久成人av| 免费不卡黄色视频| 久久人妻熟女aⅴ| 他把我摸到了高潮在线观看| 亚洲欧美日韩另类电影网站| 50天的宝宝边吃奶边哭怎么回事| 69精品国产乱码久久久| 久久香蕉国产精品| 亚洲欧美日韩无卡精品| 一区二区三区激情视频| 欧美最黄视频在线播放免费 | 欧美日本中文国产一区发布| 身体一侧抽搐| 亚洲午夜理论影院| www.熟女人妻精品国产| 桃红色精品国产亚洲av| 一二三四社区在线视频社区8| 老司机靠b影院| 亚洲av片天天在线观看| 99久久综合精品五月天人人| 人成视频在线观看免费观看| 岛国视频午夜一区免费看| 伊人久久大香线蕉亚洲五| 超碰成人久久| 看免费av毛片| 国产区一区二久久| 一边摸一边抽搐一进一小说| 啦啦啦 在线观看视频| 巨乳人妻的诱惑在线观看| 成人永久免费在线观看视频| 99久久综合精品五月天人人| 午夜精品在线福利| 久久人妻熟女aⅴ| 亚洲精品美女久久av网站| 国产在线观看jvid| 在线观看免费视频日本深夜| 国产精品电影一区二区三区| 亚洲精品国产一区二区精华液| 亚洲色图 男人天堂 中文字幕| 日本免费a在线| 国产一区二区激情短视频| 国产又爽黄色视频| 新久久久久国产一级毛片| 亚洲色图av天堂| 国产伦一二天堂av在线观看| 国产精品亚洲av一区麻豆| 日韩高清综合在线| 啪啪无遮挡十八禁网站| 色哟哟哟哟哟哟| 国产aⅴ精品一区二区三区波| 又大又爽又粗| 日韩av在线大香蕉| av在线播放免费不卡| 精品福利观看| 国产1区2区3区精品| 亚洲欧美精品综合久久99| 日韩人妻精品一区2区三区| 黑丝袜美女国产一区| 黄色女人牲交| 久久欧美精品欧美久久欧美| 久热爱精品视频在线9| av电影中文网址| 国产精品 欧美亚洲| 在线观看免费午夜福利视频| 久久久久久久精品吃奶| 久久人妻福利社区极品人妻图片| 国产在线观看jvid| 波多野结衣av一区二区av| 又黄又爽又免费观看的视频| videosex国产| 亚洲 国产 在线| 欧美激情久久久久久爽电影 | 久久精品成人免费网站| av在线天堂中文字幕 | 一夜夜www| 免费在线观看完整版高清| 欧美一区二区精品小视频在线| 757午夜福利合集在线观看| 美女福利国产在线| 大码成人一级视频| 99香蕉大伊视频| 亚洲av熟女| 99热只有精品国产| 91成人精品电影| 精品一区二区三区视频在线观看免费 | 99在线人妻在线中文字幕| 国产不卡一卡二| 超碰97精品在线观看| 韩国av一区二区三区四区| 麻豆国产av国片精品| 亚洲av日韩精品久久久久久密| 国产成年人精品一区二区 | 如日韩欧美国产精品一区二区三区| 无人区码免费观看不卡| 久久午夜亚洲精品久久| 国产色视频综合| 国产精品成人在线| 1024视频免费在线观看| 亚洲五月天丁香| 啪啪无遮挡十八禁网站| 国产黄色免费在线视频| 国产精品一区二区三区四区久久 | 女性生殖器流出的白浆| 精品乱码久久久久久99久播| 精品电影一区二区在线| 国产亚洲精品久久久久久毛片| 最近最新中文字幕大全电影3 | 色婷婷av一区二区三区视频| 国产区一区二久久| 亚洲人成伊人成综合网2020| 首页视频小说图片口味搜索| 亚洲欧美精品综合久久99| 麻豆国产av国片精品| av有码第一页| 久久国产精品男人的天堂亚洲| 国产精品综合久久久久久久免费 | 中文字幕人妻熟女乱码| 高清毛片免费观看视频网站 | 久久精品影院6| 老熟妇仑乱视频hdxx| 三上悠亚av全集在线观看| 老汉色av国产亚洲站长工具| 村上凉子中文字幕在线| 在线观看www视频免费| 亚洲精品在线美女| 激情视频va一区二区三区| 动漫黄色视频在线观看| 一边摸一边做爽爽视频免费| 色哟哟哟哟哟哟| 超碰97精品在线观看| 午夜视频精品福利| 欧美日韩亚洲高清精品| 色哟哟哟哟哟哟| 亚洲精品国产区一区二| 国产精品影院久久| 精品人妻1区二区| 久久国产精品男人的天堂亚洲| 亚洲专区中文字幕在线| 少妇粗大呻吟视频| 日韩欧美国产一区二区入口| 天堂俺去俺来也www色官网| 老司机亚洲免费影院| 精品国产乱子伦一区二区三区| 女人精品久久久久毛片| 老司机午夜福利在线观看视频| 国产亚洲欧美精品永久| 国产精品永久免费网站| 亚洲精品久久午夜乱码| 长腿黑丝高跟| 在线观看一区二区三区激情| 男女床上黄色一级片免费看| www.www免费av| 色综合欧美亚洲国产小说| 91精品国产国语对白视频| 欧美 亚洲 国产 日韩一| 欧美日韩av久久| 亚洲专区字幕在线| 国产精品美女特级片免费视频播放器 | 国产精华一区二区三区| 成人黄色视频免费在线看| 亚洲一卡2卡3卡4卡5卡精品中文| av福利片在线| 亚洲人成伊人成综合网2020| 欧美中文日本在线观看视频| 性欧美人与动物交配| www国产在线视频色| 免费久久久久久久精品成人欧美视频| 熟女少妇亚洲综合色aaa.| 校园春色视频在线观看| 成人三级做爰电影| av网站在线播放免费| 成熟少妇高潮喷水视频| 黄色丝袜av网址大全| 久久人人爽av亚洲精品天堂| 久久热在线av| 亚洲色图综合在线观看| 人妻丰满熟妇av一区二区三区| 熟女少妇亚洲综合色aaa.| 人妻久久中文字幕网| 国产精品永久免费网站| 在线播放国产精品三级| 制服人妻中文乱码| 人人妻人人添人人爽欧美一区卜| 亚洲在线自拍视频| 美女福利国产在线| 日韩国内少妇激情av| 侵犯人妻中文字幕一二三四区| 九色亚洲精品在线播放| 侵犯人妻中文字幕一二三四区| 色哟哟哟哟哟哟| 日本免费一区二区三区高清不卡 | 变态另类成人亚洲欧美熟女 | 无限看片的www在线观看| 欧美性长视频在线观看| 一二三四在线观看免费中文在| 亚洲成人免费电影在线观看| 久久香蕉精品热| 国产欧美日韩一区二区三区在线| 99久久人妻综合| 精品国产一区二区三区四区第35| 日韩欧美一区二区三区在线观看| 宅男免费午夜| 欧美激情极品国产一区二区三区| 亚洲av电影在线进入| 男女之事视频高清在线观看| av网站免费在线观看视频| 91国产中文字幕| 在线av久久热| svipshipincom国产片| 午夜免费鲁丝| 日韩免费高清中文字幕av| 午夜激情av网站| 午夜精品国产一区二区电影| 国产成人系列免费观看| av天堂在线播放| 亚洲国产精品999在线| 国产真人三级小视频在线观看| 国产精品国产av在线观看| 国产一区二区三区综合在线观看| 黄色视频,在线免费观看| 色精品久久人妻99蜜桃| 国产av又大| 可以免费在线观看a视频的电影网站| √禁漫天堂资源中文www| 正在播放国产对白刺激| 日韩大尺度精品在线看网址 | 国产高清国产精品国产三级| 人妻丰满熟妇av一区二区三区| 欧美日韩国产mv在线观看视频| 免费日韩欧美在线观看| 午夜两性在线视频| 天天躁夜夜躁狠狠躁躁| 精品一区二区三卡| 亚洲熟妇中文字幕五十中出 | 最近最新中文字幕大全电影3 | 日韩中文字幕欧美一区二区| www.自偷自拍.com| 精品一品国产午夜福利视频| 亚洲久久久国产精品| 欧美成人性av电影在线观看| 午夜精品在线福利| 婷婷丁香在线五月| 亚洲人成电影免费在线| 国产区一区二久久| 91老司机精品| 美女福利国产在线| 侵犯人妻中文字幕一二三四区| 亚洲成人久久性| 亚洲第一av免费看| 欧美老熟妇乱子伦牲交| 大陆偷拍与自拍| 一区二区三区精品91| 久久这里只有精品19| 色尼玛亚洲综合影院| 亚洲国产欧美网| 精品少妇一区二区三区视频日本电影| 悠悠久久av| 在线观看一区二区三区激情| 深夜精品福利| 国产伦人伦偷精品视频| 国产精品久久视频播放| 色婷婷av一区二区三区视频| 久久久久久久久免费视频了| 国产精品久久久久成人av| 男人舔女人的私密视频| 18禁国产床啪视频网站| 久久久精品欧美日韩精品| 桃红色精品国产亚洲av| 国产精品久久电影中文字幕| av电影中文网址| 韩国精品一区二区三区| 国产蜜桃级精品一区二区三区| 99国产综合亚洲精品| 高清毛片免费观看视频网站 | 看免费av毛片| 国产精品一区二区精品视频观看| 91成人精品电影| 免费高清视频大片| 国产精品久久久久久人妻精品电影| 免费在线观看黄色视频的| 国产aⅴ精品一区二区三区波| 日韩国内少妇激情av| 亚洲精品一区av在线观看| 一二三四社区在线视频社区8| 亚洲aⅴ乱码一区二区在线播放 | 校园春色视频在线观看| 免费在线观看视频国产中文字幕亚洲| 国产精品综合久久久久久久免费 | 天堂俺去俺来也www色官网| 男女做爰动态图高潮gif福利片 | 搡老岳熟女国产| 免费高清在线观看日韩| 99热国产这里只有精品6| 精品第一国产精品| 少妇 在线观看| 精品久久久久久,| a在线观看视频网站| 欧美中文综合在线视频| 精品卡一卡二卡四卡免费| 黄色 视频免费看| 欧美午夜高清在线| 69av精品久久久久久| 欧美丝袜亚洲另类 | 真人做人爱边吃奶动态| 动漫黄色视频在线观看| 日韩欧美一区二区三区在线观看| 欧美日韩乱码在线| 亚洲成a人片在线一区二区| 日日摸夜夜添夜夜添小说| 亚洲中文字幕日韩| 国产一区在线观看成人免费| 极品教师在线免费播放| 三级毛片av免费| 99久久99久久久精品蜜桃| 亚洲精品国产区一区二| 国产真人三级小视频在线观看| av有码第一页| 后天国语完整版免费观看| 丰满人妻熟妇乱又伦精品不卡| 国产免费男女视频| 国产男靠女视频免费网站| 琪琪午夜伦伦电影理论片6080| 欧美成狂野欧美在线观看| 午夜久久久在线观看| 国产伦一二天堂av在线观看| 男人舔女人的私密视频| 亚洲精品一二三| 麻豆av在线久日| 成人亚洲精品av一区二区 | 日本黄色日本黄色录像| 亚洲精华国产精华精| 久久久国产成人精品二区 | 99热只有精品国产| 最近最新免费中文字幕在线| 亚洲一区高清亚洲精品| 精品高清国产在线一区| 少妇粗大呻吟视频| 亚洲精品美女久久av网站| 这个男人来自地球电影免费观看| 亚洲男人的天堂狠狠| 午夜福利影视在线免费观看| 亚洲av成人不卡在线观看播放网| 亚洲男人天堂网一区| 国产成人影院久久av| www.精华液| 热99国产精品久久久久久7| 村上凉子中文字幕在线| 俄罗斯特黄特色一大片| 日韩一卡2卡3卡4卡2021年| 美女福利国产在线| 99久久久亚洲精品蜜臀av| 亚洲成人国产一区在线观看| 久久天躁狠狠躁夜夜2o2o| 人妻久久中文字幕网| 手机成人av网站| 国产蜜桃级精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 操美女的视频在线观看| 亚洲中文日韩欧美视频| 日韩精品青青久久久久久| 在线天堂中文资源库| 国产熟女午夜一区二区三区| 人妻丰满熟妇av一区二区三区| 天堂俺去俺来也www色官网| 18禁国产床啪视频网站| www.www免费av| 国产精品一区二区在线不卡| 女同久久另类99精品国产91| 欧美午夜高清在线| 日韩高清综合在线| 黄片播放在线免费| 动漫黄色视频在线观看| 亚洲av美国av| 久久九九热精品免费| 一级毛片女人18水好多| 黄色 视频免费看| 国产午夜精品久久久久久| 久久久国产成人免费| xxxhd国产人妻xxx| 欧美日韩亚洲高清精品| 亚洲av成人av| 日韩有码中文字幕| 成人精品一区二区免费| 亚洲精品在线观看二区| 90打野战视频偷拍视频| 性少妇av在线| 又黄又粗又硬又大视频| 国产真人三级小视频在线观看| 满18在线观看网站| 久久久久国产一级毛片高清牌| 99国产精品一区二区三区| 久久狼人影院| 麻豆av在线久日| 国产日韩一区二区三区精品不卡| 夜夜看夜夜爽夜夜摸 | cao死你这个sao货| 伦理电影免费视频| 国产亚洲精品久久久久5区| 99热国产这里只有精品6| 亚洲精品av麻豆狂野| 亚洲少妇的诱惑av| 大陆偷拍与自拍| 久久久久国产一级毛片高清牌| 国产精品99久久99久久久不卡| 久久精品国产清高在天天线| 欧美色视频一区免费| 成人精品一区二区免费| 精品人妻1区二区| 亚洲av熟女| 国产av一区在线观看免费| 国产无遮挡羞羞视频在线观看| www.999成人在线观看| 精品国产一区二区久久| 日韩免费高清中文字幕av| 美女午夜性视频免费| 这个男人来自地球电影免费观看| 精品一区二区三区av网在线观看| 天天躁夜夜躁狠狠躁躁| 男女之事视频高清在线观看| 成人亚洲精品av一区二区 | 日韩欧美免费精品| 精品久久久久久成人av| 亚洲精品成人av观看孕妇| 天堂√8在线中文| 久久久久久亚洲精品国产蜜桃av| 老司机在亚洲福利影院| 岛国视频午夜一区免费看| 99热国产这里只有精品6| 两性午夜刺激爽爽歪歪视频在线观看 | 色综合婷婷激情| 欧美av亚洲av综合av国产av| 一级片免费观看大全| 久久青草综合色| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 两个人看的免费小视频| 男人舔女人的私密视频| 精品免费久久久久久久清纯| www.熟女人妻精品国产| 亚洲 国产 在线| 性色av乱码一区二区三区2| 国产人伦9x9x在线观看| 老司机靠b影院| 俄罗斯特黄特色一大片| 国内久久婷婷六月综合欲色啪| 身体一侧抽搐| 中文字幕另类日韩欧美亚洲嫩草| 免费久久久久久久精品成人欧美视频| 女同久久另类99精品国产91| 亚洲av成人av| 国产精品亚洲av一区麻豆| 午夜两性在线视频| 久久久精品欧美日韩精品| 欧美中文综合在线视频| 每晚都被弄得嗷嗷叫到高潮| 欧美在线一区亚洲| 夜夜躁狠狠躁天天躁| 久久久久亚洲av毛片大全| a级毛片在线看网站| 窝窝影院91人妻| 国产精品久久久av美女十八| av超薄肉色丝袜交足视频| 黄色怎么调成土黄色| 久久久国产欧美日韩av| 久久久久九九精品影院| 一区二区三区国产精品乱码| 国产精品久久久av美女十八| 欧美激情久久久久久爽电影 | 亚洲国产精品999在线| 在线观看午夜福利视频| 热re99久久国产66热| av欧美777| 99香蕉大伊视频| 色婷婷久久久亚洲欧美| 男女床上黄色一级片免费看| 免费在线观看亚洲国产| 9热在线视频观看99| 亚洲男人天堂网一区| 久久久久国产一级毛片高清牌| 欧美亚洲日本最大视频资源| 日本a在线网址| 怎么达到女性高潮| 亚洲国产毛片av蜜桃av| 青草久久国产| 无限看片的www在线观看| 88av欧美| 亚洲avbb在线观看| 自线自在国产av| 国产成人影院久久av| 久久精品亚洲精品国产色婷小说| 亚洲第一青青草原| 欧美丝袜亚洲另类 | 久久精品91蜜桃| 亚洲欧美一区二区三区黑人| 夜夜躁狠狠躁天天躁| 亚洲成人免费av在线播放| 麻豆国产av国片精品| 国产精品98久久久久久宅男小说| 淫秽高清视频在线观看| 亚洲一区中文字幕在线| 亚洲全国av大片| 亚洲人成电影观看| 97人妻天天添夜夜摸| 黑丝袜美女国产一区| 天堂影院成人在线观看| 在线观看一区二区三区激情| 欧美人与性动交α欧美精品济南到| a级毛片黄视频| 精品久久久久久久毛片微露脸| 国产aⅴ精品一区二区三区波| 精品电影一区二区在线| 这个男人来自地球电影免费观看| 亚洲中文日韩欧美视频| 亚洲视频免费观看视频| 夜夜夜夜夜久久久久| 老鸭窝网址在线观看| 91av网站免费观看| 三级毛片av免费| 老司机深夜福利视频在线观看| 免费在线观看视频国产中文字幕亚洲| 日本免费一区二区三区高清不卡 | 国产精品久久久av美女十八| 久久久久国内视频| 91九色精品人成在线观看| 日本 av在线| 欧美人与性动交α欧美精品济南到| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久av美女十八| 精品熟女少妇八av免费久了| 国产成+人综合+亚洲专区| 国产欧美日韩一区二区三| 性少妇av在线| 久久婷婷成人综合色麻豆| 久久精品人人爽人人爽视色| av网站在线播放免费| 日本黄色视频三级网站网址| 香蕉丝袜av| 50天的宝宝边吃奶边哭怎么回事| 国产主播在线观看一区二区| 极品人妻少妇av视频| 国产高清视频在线播放一区| 亚洲人成电影观看| 看黄色毛片网站| 免费高清视频大片| 男女高潮啪啪啪动态图| 男人舔女人下体高潮全视频| 国产97色在线日韩免费| 纯流量卡能插随身wifi吗| 欧美精品一区二区免费开放| 岛国在线观看网站| 91av网站免费观看| 精品久久久精品久久久| 国产精品九九99| 精品一区二区三区av网在线观看| 啦啦啦免费观看视频1| 高清毛片免费观看视频网站 | 热99国产精品久久久久久7| 亚洲,欧美精品.| 法律面前人人平等表现在哪些方面| 亚洲自拍偷在线| 亚洲一卡2卡3卡4卡5卡精品中文| av免费在线观看网站| 亚洲七黄色美女视频| 亚洲av成人一区二区三| 国产亚洲欧美在线一区二区| 怎么达到女性高潮| 88av欧美| 9色porny在线观看| 久久中文看片网| 级片在线观看| 国产欧美日韩一区二区三| 搡老乐熟女国产| 亚洲欧美日韩无卡精品| 国产伦一二天堂av在线观看| 国内毛片毛片毛片毛片毛片| 久久人人爽av亚洲精品天堂| 久久亚洲精品不卡| 久99久视频精品免费| 午夜免费成人在线视频| 男女高潮啪啪啪动态图| 久久久久久大精品| 久久香蕉精品热| 精品无人区乱码1区二区| 日韩 欧美 亚洲 中文字幕| 免费观看精品视频网站| 免费人成视频x8x8入口观看| 久久久久久人人人人人| svipshipincom国产片| 国产免费现黄频在线看| 一二三四在线观看免费中文在| 99riav亚洲国产免费| 最好的美女福利视频网| 制服人妻中文乱码| 丝袜美腿诱惑在线| 国产亚洲精品久久久久久毛片| 97超级碰碰碰精品色视频在线观看| 精品人妻在线不人妻| 黄色毛片三级朝国网站| 久久精品成人免费网站| 日韩免费av在线播放| 成人亚洲精品av一区二区 | 国产99白浆流出| 中文字幕色久视频| 国产一卡二卡三卡精品| 久久精品国产综合久久久| 中文欧美无线码| 国产不卡一卡二| 99久久99久久久精品蜜桃| 久久午夜亚洲精品久久| 精品一区二区三卡|