8个番石榴品种的转录组测序

张丽梅; 陈钟佃; 张朝坤; 魏秀清; 许家辉

doi:10.19303/j.issn.1008-0384.2021.11.008

8个番石榴品种的转录组测序

doi: 10.19303/j.issn.1008-0384.2021.11.008

1.
福建省农业科学院果树研究所, 福建　福州　350013
2.
福建省农业科学院农业生态研究所, 福建　福州　350013
3.
福建省漳州市农业科学研究所, 福建　漳州　363000

基金项目: 福建省农业科学院科技创新团队项目（CXTD2021008-2）；中国热带农业科学院基本科研业务费专项项目（1630032020013）；福建省科技计划区域发展项目（2016N3021)

详细信息

作者简介:
张丽梅（1974−），女，副研究员，研究方向：果树栽培、生理生化与分子技术（E-mail: lm155195@sina.com）

通讯作者:
许家辉（1971−），男，硕士，研究员，研究方向：果树育种、栽培及生理生化（E-mail: xjhui577@163.com）

中图分类号: S 667
计量
- 文章访问数: 539
- HTML全文浏览量: 98
- PDF下载量: 25
- 被引次数: 0
出版历程
- 收稿日期: 2021-06-27
- 修回日期: 2021-12-08
- 网络出版日期: 2021-12-30
- 刊出日期: 2021-11-28

Transcriptome Sequencing on Eight Psidium guajava Cultivars

1.
Fruit Research Institute, Fujian Academy of Agricultural Sciences, Fuzhou, Fujian　350013, China
2.
Agricultural Ecology Institute, Fujian Academy of Agricultural Sciences, Fuzhou, Fujian　350013, China
3.
Zhangzhou Institute of Agricultural Sciences of Fujian Province, Zhangzhou, Fujian　363000, China

摘要

摘要: 目的获得番石榴转录本数据，探讨不同番石榴品种转录本之间差异。方法研究通过Illumina HiSeq 2500平台对红心（HX）、彩虹（CH）、紫红（ZH）、水晶（SJ）、秀珍（XZ）、彩叶变（CYB）、彩叶（CY）和紫星（ZX）等8个番石榴品种的8～9成熟果肉进行转录组测序，通过组装及数据库比对，获得大量番石榴转录本注释信息。结果组装获得了126 979条 Unigene，平均长度766 bp，N₅₀为1654 bp，N₉₀为267 bp；8个番石榴品种间共有27 274个基因为共表达基因。通过主成分分析与差异基因表达聚类分析，发现品种水晶（SJ）和彩叶变（CYB）之间基因功能表达差异较小，品种紫红（ZH）与其他品种之间基因表达差异较大。聚类分析将品种水晶（SJ）和彩叶变（CYB）聚为一类，品种红心（HX）、彩虹（CH）、秀珍（XZ）、彩叶（CY）、紫星（ZX）归为一类，但品种紫红（ZH）与其他品种间距离较远。结论本研究结果为后续番石榴种质资源分类、功能基因表达、功能基因鉴定提供了大量可靠的转录组数据。
- 番石榴 /
- 转录组 /
- 测序 /
- 功能注释
Abstract: Objective Transcriptomes of different guava cultivars were obtained for comparison. Methods Transcriptome sequencing by Illumina Hwaseq 2 500 was conducted on 8-9 mature fruit pulp from 8 varieties of Psidium guajava Linn, namely Red Heart (HX), Rainbow (CH), Purplish Red (ZH), Crystal (SJ), Shuzhen (XZ), Colored-leaf B (CYB), Colored-leaf (CY), and Purple Star (ZX). Results By comparing to the database, 126979 Unigenes of an average length of 766 bp with a N₅₀ of 1654 bp and a N₉₀ of 267 bp were identified. There were 27274 genes co-expressed in the 8 cultivars. The PCA and differential gene expression cluster analyses showed that the differences on the function expressions between the genes of SJ and CYB was less than those between ZH and other cultivars. SJ and CYB were clustered in a category and HX, CH, XZ, CY, and ZX in a different category, while ZH significantly differentiated from the other cultivars. Conclusion The substantial amount of transcription data secured from this study could be useful for the classification as well as the expression and identification of functional genes of guava germplasms.
- Psidium guajava Linn. /
- transcriptome /
- sequencing /
- functional annotation

HTML全文

图 1 Transcript（A）和Unigene（B）长度分布

Figure 1. Distributions of transcript length (A) and Unigene length (B)

下载: 全尺寸图片幻灯片

图 2 番石榴Unigene NR数据库注释分析

Figure 2. Annotation of NR database on guava Unigenes

下载: 全尺寸图片幻灯片

图 3 番石榴Unigene Go分类

Figure 3. GO functional classification of guava Unigenes

下载: 全尺寸图片幻灯片

图 4 番石榴Unigene的KOG数据库分类

Figure 4. KOG classification of guava Unigenes

下载: 全尺寸图片幻灯片

图 5 番石榴Unigene的KEGG数据库分类

Figure 5. KEGG classification of guava Unigenes

下载: 全尺寸图片幻灯片

图 6 番石榴PCA主成分分析

Figure 6. PCA on guava

下载: 全尺寸图片幻灯片

图 7 样本遗传距离分析

Figure 7. Heat map of distance analysis on guava samples

下载: 全尺寸图片幻灯片

图 8 番石榴基因共表达韦恩图

Figure 8. TPM Venn diagram of guava cultivars

下载: 全尺寸图片幻灯片

表 1 样本质量控制数据统计

Table 1. Statistics on quality control data of samples

品种 Cultivars	总碱基数 Number of bases/bp	Reads/bp	Q₃₀/%	GC/%
红心HX	7123935342	145	95.96	49.98
彩虹CH	8029343436	144	95.47	50.68
紫红ZH	6300972075	145	95.57	50.10
水晶SJ	7681478290	145	95.60	51.23
秀珍XZ	5184950713	144	95.31	50.27
彩叶变CYB	7824848106	146	95.82	50.65
彩叶CY	7607494989	144	95.35	50.44
紫星ZX	6364083214	144	95.50	50.49

下载: 导出CSV

表 2 番石榴8个品种间基因表达差异统计分析

Table 2. Statistical analysis on gene expressions of 8 guava cultivars

对照组 CK	试验组 Experimental cultivars	上调基因数 Up-regulatory genes	下调基因数 Down-regulatory genes	KEGG富集途径 KEGG pathway
对照组 CK	试验组 Experimental cultivars	上调基因数 Up-regulatory genes	下调基因数 Down-regulatory genes	1	2	3	4
SJ	ZX	2601	8289	Y	Y	Y	N
SJ	ZH	5959	11471	Y	Y	Y	N
SJ	CH	1996	7354	Y	Y	Y	N
SJ	XZ	1994	7178	Y	Y	Y	N
CYB	CY	2640	6327	Y	Y	Y	N
SJ	HX	2187	8443	N	Y	Y	Y
ZX	ZH	4876	4232	Y	Y	N	Y
CH	HX	13355	1896	Y	N	N	N
注：KEGG富集途径1：核糖体；2：碳代谢；3：氨基酸生物合成；4：淀粉和蔗糖代谢。Y表示富集到该代谢途径，N表示未富集到该代谢途径。 Note: KEGG pathway 1: ribosome; 2: carbon metabolism; 3: amino acid biosynthesis; 4: starch and sucrose metabolisms; Y: enriched; N: non-enriched.

下载: 导出CSV

参考文献(35)

[1]	GUTIÉRREZ R M P, MITCHELL S, SOLIS R V. Psidium guajava: A review of its traditional uses, phytochemistry and pharmacology [J]. Journal of Ethnopharmacology, 2008, 117(1): 1−27. doi: 10.1016/j.jep.2008.01.025
[2]	HUANG C S, YIN M C, CHIU L C. Antihyperglycemic and antioxidative potential of Psidium guajava fruit in streptozotocin-induced diabetic rats [J]. Food and Chemical Toxicology, 2011, 49(9): 2189−2195. doi: 10.1016/j.fct.2011.05.032
[3]	CHENG F C, SHEN S C, WU J S B. Effect of guava (Psidium guajavaL. ) leaf extract on glucose uptake in rat hepatocytes [J]. Journal of Food Science, 2009, 74(5): H132−H138. doi: 10.1111/j.1750-3841.2009.01149.x
[4]	UBOH F E, OKON I E, EKONG M B. Effect of aqueous extract of Psidium guajava leaves on liver enzymes, histological integrity and hematological indices in rats [J]. Gastroenterology Research, 2010, 3(1): 32−38.
[5]	HAN E H, HWANG Y P, KIM H G, et al. Ethyl acetate extract of Psidium guajava inhibits IgE-mediated allergic responses by blocking FcεRI signaling [J]. Food and Chemical Toxicology, 2011, 49(1): 100−108. doi: 10.1016/j.fct.2010.10.003
[6]	RATTANACHAIKUNSOPON P, PHUMKHACHOM P. Contents and antibacterial activity of flavonoids extracted from leaves of Psidium guajava [J]. Journal of Medicinal Plants Research, 2010, 4(5): 393−396.
[7]	易富, 蔡敏琪, 李娅, 等. 番石榴叶水提取物对11种细菌的体外抑菌实验 [J]. 西南国防医药, 2013, 23(11):1168−1169. doi: 10.3969/j.issn.1004-0188.2013.11.003 YI F, CAI M Q, LI Y, et al. Study on antimicrobial effects of guava leaf water extract on 11 kinds of bacteria in vitro [J]. Medical Journal of National Defending Forces in Southwest China, 2013, 23(11): 1168−1169.(in Chinese) doi: 10.3969/j.issn.1004-0188.2013.11.003
[8]	ADEYEMI O S, SYKES M L, AKANJI MA, et al. Anti-trypanosomal and cytotoxic activity of ethanolic extracts of Psidium guajava leaves in Alamar Blue based assays [J]. Veterinarski Arhiv, 2011, 81(5): 623−633.
[9]	LEE W C, MAHMUD R, NOORDIN R, et al. Free radicals scavenging activity, cytotoxicity and Anti-parasitic activity of essential oil of Psidium guajava L. leaves against Toxoplasma gondii [J]. Journal of Essential Oil Bearing Plants, 2013, 16(1): 32−38. doi: 10.1080/0972060X.2013.764196
[10]	BHAGAVATHY S, MAHENDIRAN C, KANCHANA R. Identification of glucosyl transferase inhibitors from Psidium guajava against Streptococcus mutans in dental caries [J]. Journal of Traditional and Complementary Medicine, 2019, 9(2): 124−137. doi: 10.1016/j.jtcme.2017.09.003
[11]	MEHTA V V, RAJESH G, RAO A, et al. Antimicrobial efficacy of Punica granatum mesocarp, Nelumbo nucifera leaf, Psidium guajava leaf and Coffea canephora extract on common oral pathogens: An In-vitro study [J]. Journal of Clinical and Diagnostic Research, 2014, 8(7): ZC65−ZC68.
[12]	郭胜男, 刘洪斌, 李东华, 等. 番石榴叶总黄酮对糖尿病小鼠肝脏葡萄糖代谢及胰岛素信号通路的影响 [J]. 中国实验方剂学杂志, 2015, 21(4):166−170. GUO S N, LIU H B, LI D H, et al. Effects of flavonoids from Psidium guajava leaves on hepatic glucose metabolism and insulin signaling pathway in diabetic mice [J]. Chinese Journal of Experimental Traditional Medical Formulae, 2015, 21(4): 166−170.(in Chinese)
[13]	EL-AHMADY S H, ASHOUR M L, WINK M. Chemical composition and anti-inflammatory activity of the essential oils of Psidium guajava fruits and leaves [J]. Journal of Essential Oil Research, 2013, 25(6): 475−481. doi: 10.1080/10412905.2013.796498
[14]	BATISTA SILVA W, COSME SILVA G M, SANTANA D B, et al. Chitosan delays ripening and ROS production in guava (Psidium guajava L. ) fruit [J]. Food Chemistry, 2018, 242: 232−238. doi: 10.1016/j.foodchem.2017.09.052
[15]	SANTOS T M, SOUZA FILHO M D S M, SILVA E D O, et al. Enhancing storage stability of guava with tannic acid-crosslinked zein coatings [J]. Food Chemistry, 2018, 257: 252−258. doi: 10.1016/j.foodchem.2018.03.021
[16]	NAIR M S, SAXENA A, KAUR C. Effect of chitosan and alginate based coatings enriched with pomegranate peel extract to extend the postharvest quality of guava (Psidium guajava L. ) [J]. Food Chemistry, 2018, 240: 245−252. doi: 10.1016/j.foodchem.2017.07.122
[17]	DE AQUINO A B, BLANK A F, DE AQUINO SANTANA L C L. Impact of edible chitosan-cassava starch coatings enriched with Lippia gracilis Schauer genotype mixtures on the shelf life of guavas (Psidium guajava L. ) during storage at room temperature [J]. Food Chemistry, 2015, 171: 108−116. doi: 10.1016/j.foodchem.2014.08.077
[18]	WANG H, DU Y J, SONG H C. Α-Glucosidase and α-amylase inhibitory activities of guava leaves [J]. Food Chemistry, 2010, 123(1): 6−13. doi: 10.1016/j.foodchem.2010.03.088
[19]	NUNES J C, LAGO M G, CASTELO-BRANCO V N, et al. Effect of drying method on volatile compounds, phenolic profile and antioxidant capacity of guava powders [J]. Food Chemistry, 2016, 197: 881−890. doi: 10.1016/j.foodchem.2015.11.050
[20]	SACCHETTI G, MAIETTI S, MUZZOLI M, et al. Comparative evaluation of 11 essential oils of different origin as functional antioxidants, antiradicals and antimicrobials in foods [J]. Food Chemistry, 2005, 91(4): 621−632. doi: 10.1016/j.foodchem.2004.06.031
[21]	FERNANDES M R V, DIAS A L T, CARVALHO R R, et al. Antioxidant and antimicrobial activities of Psidium guajava L. spray dried extracts [J]. Industrial Crops and Products, 2014, 60: 39−44. doi: 10.1016/j.indcrop.2014.05.049
[22]	张宏康, 林奕楠, 李德荣, 等. 番石榴泡腾片固体饮料研制 [J]. 食品工业, 2015, 36(1):4−7. ZHANG H K, LIN Y N, LI D R, et al. Preparation of solid beverage of effervescent tablet from the extract of Psidium guajava [J]. The Food Industry, 2015, 36(1): 4−7.(in Chinese)
[23]	刘蒙佳, 周强, 陈淑娣. 番石榴汁茶酒的发酵工艺研究 [J]. 茶叶科学, 2014, 34(1):21−28. doi: 10.3969/j.issn.1000-369X.2014.01.004 LIU M J, ZHOU Q, CHEN S D. Study on fermentation technology of tea wine with guava juice [J]. Journal of Tea Science, 2014, 34(1): 21−28.(in Chinese) doi: 10.3969/j.issn.1000-369X.2014.01.004
[24]	BOLGER A M, LOHSE M, USADEL B. Trimmomatic: a flexible trimmer for Illumina sequence data [J]. Bioinformatics (Oxford, England), 2014, 30(15): 2114−2120. doi: 10.1093/bioinformatics/btu170
[25]	HAAS B J, PAPANICOLAOU A, YASSOUR M, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis [J]. Nature Protocols, 2013, 8(8): 1494−1512. doi: 10.1038/nprot.2013.084
[26]	ALTSCHUL S F, MADDEN T L, SCHÄFFER A A, et al. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs [J]. Nucleic Acids Research, 1997, 25(17): 3389−3402. doi: 10.1093/nar/25.17.3389
[27]	LANGMEAD B, SALZBERG S L. Fast gapped-read alignment with Bowtie 2 [J]. Nature Methods, 2012, 9(4): 357−359. doi: 10.1038/nmeth.1923
[28]	WANG L G, WANG S Q, LI W. RSeQC: quality control of RNA-seq experiments [J]. Bioinformatics, 2012, 28(16): 2184−2185. doi: 10.1093/bioinformatics/bts356
[29]	QUINLAN A R, HALL I M. BEDTools: a flexible suite of utilities for comparing genomic features [J]. Bioinformatics, 2010, 26(6): 841−842. doi: 10.1093/bioinformatics/btq033
[30]	PATRO R, DUGGAL G, LOVE M I, et al. Salmon provides fast and bias-aware quantification of transcript expression [J]. Nature Methods, 2017, 14(4): 417−419. doi: 10.1038/nmeth.4197
[31]	LANGFELDER P, HORVATH S. WGCNA: an R package for weighted correlation network analysis [J]. BMC Bioinformatics, 2008, 9: 559. doi: 10.1186/1471-2105-9-559
[32]	LOVE M I, HUBER W, ANDERS S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 [J]. Genome Biology, 2014, 15(12): 550. doi: 10.1186/s13059-014-0550-8
[33]	YU G C, WANG L G, HAN Y Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters [J]. Omics, 2012, 16(5): 284−287. doi: 10.1089/omi.2011.0118
[34]	GHOSH R, TIWARY B K, KUMAR A, et al. Guava leaf extract inhibits quorum-sensing and Chromobacterium violaceum induced Lysis of human hepatoma cells: Whole transcriptome analysis reveals differential gene expression [J]. PLoS One, 2014, 9(9): e107703. doi: 10.1371/journal.pone.0107703
[35]	张丽梅, 张朝坤, 陈洪彬, 等. 番石榴种质资源果实性状的聚类分析 [J]. 中国南方果树, 2019, 48(6):58−63. ZHANG L M, ZHANG C K, Chen H B, et al. Cluster analysis of fruit characters of guava germplasm resources [J]. South China Fruit, 2019, 48(6): 58−63.(in Chinese)