欢迎访问作物学报，今天是
作物学报

作物学报 ›› 2024, Vol. 50 ›› Issue (6): 1503-1513.doi: 10.3724/SP.J.1006.2024.34172

• 作物遗传育种·种质资源·分子遗传学 • 上一篇    下一篇

不同马铃薯淀粉含量差异的转录组学解析

赵娜1,2(), 刘宇曦1,2, 张朝澍1,2,*(), 石瑛1,2,*()   

  1. 1东北农业大学农学院, 黑龙江哈尔滨 150030
    2寒地粮食作物种质创新与生理生态教育部重点实验室, 黑龙江哈尔滨 150030
  • 收稿日期:2023-10-20 接受日期:2024-01-31 出版日期:2024-06-12 网络出版日期:2024-02-21
  • 通讯作者: * 石瑛, E-mail: shiying01@163.com; 张朝澍, E-mail: zhangchaoshu@126.com
  • 作者简介:E-mail: zhaona992024@163.com
  • 基金资助:
    黑龙江省“揭榜挂帅”科技攻关项目(2022ZXJ06B02);黑龙江省自然科学基金项目(LH2021C027)

Transcriptomic analysis of differences in the starch content of different potatoes

ZHAO Na1,2(), LIU Yu-Xi1,2, ZHANG Chao-Shu1,2,*(), SHI Ying1,2,*()   

  1. 1College of Agronomy, Northeast Agricultural University, Harbin 150030, Heilongjiang, China
    2Key Laboratory of Germplasm Enhancement, Physiology and Ecology of Food Crops in Cold Region, Ministry of Education, Harbin 150030, Heilongjiang, China
  • Received:2023-10-20 Accepted:2024-01-31 Published:2024-06-12 Published online:2024-02-21
  • Contact: * E-mail: shiying01@163.com;E-mail: zhangchaoshu@126.com
  • Supported by:
    “Open Bidding for Selecting the Best Candidates” Scientific and Technological Research Project of Heilongjiang Province(2022ZXJ06B02);Natural Science Foundation of Heilongjiang Province(LH2021C027)

RichHTML

24

PDF

724

摘要：

淀粉是马铃薯重要的品质性状之一, 广泛应用于食品、医药、石油化工等行业, 市场上对高淀粉马铃薯品种的需求逐年提高。为探究马铃薯块茎淀粉积累特性及调控的关键基因, 本研究通过RNA测序分析, 分别对高、低淀粉含量马铃薯品种大西洋(DXY)、定薯1号(DS1)的匍匐茎(a)、块茎膨大前期(b)、块茎膨大中期(c)、块茎膨大后期(d)及成熟期(e)进行转录组谱分析, 5个时期共获得9494个差异表达基因(DEGs), 结合功能注释发现差异基因主要富集在结合绑定和催化活性等分子功能中, 代谢通路分析发现差异基因主要富集在碳水化合物途径中, 其中有137个DEGs与淀粉和蔗糖代谢相关。筛选获得9个调控淀粉合成关键基因, 其中蔗糖合酶基因PGSC0003DMG400013547在大西洋整个生育期表达量较高, 且在各时期表达量均高于定薯1号; 果糖激酶基因PGSC0003DMG400026916在a、b时期大西洋的表达量显著高于定薯1号; 定薯1号中淀粉酶基因PGSC0003DMG400009891、PGSC0003DMG400001549和葡聚糖内切-1,3-β-葡萄糖苷酶基因PGSC0003DMG400024642、PGSC0003DMG400003181分别在c和e时期表达量显著高于高淀粉品种大西洋。上述基因可能是参与淀粉的合成与积累的关键调控基因。本研究为不同类型马铃薯品种块茎淀粉代谢调控机制研究提供参考。

关键词: 马铃薯, 发育时期, 淀粉合成积累, 转录组测序

Abstract:

Starch is one of the most important quality characteristics of potatoes, which is widely used in the food, medical, petrochemical and other industries. The demand for starchy potato varieties in the market is increasing year by year. To explore the characteristics and key genes of starch accumulation and regulation of potato tuber, in this study, transcriptome profiling of creeping stems (a), pre-tuber expansion (b), mid-tuber expansion (c), late tuber expansion (d), and mature stage (e) of high and low starch potato cultivars DXY and DS1 were conducted by RNA sequencing analysis, and a total of 9494 differentially expressed genes (DEGs) were identified at five stages. Binding function annotation revealed that the differential genes were mainly enriched in molecular functions such as binding and catalytic activity. Metabolic pathway analysis revealed that the differential genes were mainly enriched in carbohydrate-related metabolic pathways, and 137 DEGs were associated with starch and sucrose metabolism. Nine key genes regulating starch synthesis were examined, and the relative expression level of the sucrose synthase gene PGSC0003DMG400013547 was higher during growth period of DXY than that at any stage. The relative expression level of the fructokinase gene PGSC0003DMG400026916 was significantly higher in DXY than in DS1 during the a and b periods. The relative expression level of amylase gene PGSC0003DMG400009891, PGSC0003DMG400001549, and glucan endo-1,3-β-glucosidase gene PGSC0003DMG400024642, PGSC0003DMG400003181 in DS1 was significantly higher than that in DXY with starch-rich cultivars during the c and e periods. The described gene might be a key regulatory gene for starch synthesis and accumulation. This study provides a clue for the investigation of the regulatory mechanism of tuber starch metabolism in different potato cultivars.

Key words: potato, developmental period, accumulation of starch synthesis, transcriptome sequencing

表1

引物序列"

基因编号Gene ID 正向引物Forward primer (5'-3') 反向引物Reverse primer (5'-3')
PGSC0003DMG400026916 TGATCAAGGTCAGCGATGTG GTCTTCACGTGGAATCCTCC
PGSC0003DMG400016481 ACAGTTCTGCGTTTCAGAGC AGCAGCATGGTAAAGCTCCAA
PGSC0003DMG400013547 ACAAGGGGTGCTTTCGTTCA GTAGCCTTGTCGCCCTGATT
PGSC0003DMG400003181 CGCTGAGAGGAATTTCGGGT TTTGGTTGCCCTTTGTTGCC
PGSC0003DMG400024642 TACAAACGGACCTTGCTTTA AGTAGCAGGATAAACACAGC
PGSC0003DMG400009891 GATCATATTGCAGGCATTCG CAAGCAGTTGTGAAACCAG
EF1a GATGTTGTGCCAAAGGATGT AACTTGGTCAATGCGAGA

图1

定薯1号、大西洋块茎不同发育时期图 Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。"

图2

定薯1号vs大西洋不同时期淀粉含量和淀粉含量比率的变化 A: 淀粉含量柱状图; B: 淀粉含量比率折线图。Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。星号表示通过单因素方差分析与定薯1号相比具有显著性(** P < 0.01, *** P < 0.001), ns表示无显著差异(P > 0.05)。"

表2

测序数据质量统计"

样品名称
Sample name		 有效数据
Clean reads		 基因组比对
Mapped reads		 基因组比对率
Mapped reads (%)		 GC含量
GC content (%)		 ≥Q20 (%) ≥Q30 (%)
Da-1 41,710,042 35,971,385 86.24 42.61 97.77 93.98
Da-2 43,243,492 37,736,387 87.26 42.56 97.67 93.48
Da-3 38,555,946 34,064,046 88.35 42.86 97.96 94.21
Db-1 40,147,012 34,518,026 85.98 42.47 97.95 94.06
Db-2 41,452,450 35,287,769 85.13 43.29 97.64 93.87
Db-3 41,485,436 35,915,281 86.57 42.84 98.07 94.60
Dc-1 38,488,084 30,296,884 78.72 41.65 97.88 93.99
Dc-2 38,988,180 30,707,745 78.76 41.46 97.99 94.24
Dc-3 40,364,158 34,737,460 86.06 42.59 98.10 94.55
Dd-1 39,833,928 31,152,181 78.21 41.68 97.96 94.32
Dd-2 41,331,412 32,215,855 77.95 41.58 97.97 94.34
Dd-3 41,865,794 33,566,589 80.18 41.71 97.85 93.99
De-1 41,161,928 35,155,400 85.41 44.47 97.31 93.57
De-2 38,363,322 29,762,706 77.58 42.16 97.46 93.60
De-3 42,150,428 32,190,017 76.37 42.13 97.71 93.98
Xa-1 40,619,164 35,357,453 87.05 42.54 98.28 94.87
Xa-2 40,601,792 35,382,889 87.15 42.68 97.95 94.07
Xa-3 41,200,382 35,910,567 87.16 42.65 97.67 93.37
Xb-1 44,599,490 38,997,447 87.44 42.85 98.26 94.79
Xb-2 41,602,324 35,537,658 85.42 42.59 98.12 94.52
Xb-3 42,202,702 35,250,119 83.53 43.60 97.99 94.49
Xc-1 40,274,444 31,599,634 78.46 42.08 98.18 94.72
Xc-2 39,511,076 32,460,237 82.15 42.20 97.88 94.03
Xc-3 40,493,952 31,485,795 77.75 41.53 98.09 94.57
Xd-1 40,558,872 32,781,195 80.82 41.94 98.16 94.59
Xd-2 40,557,808 32,991,680 81.34 41.88 98.31 94.91
Xd-3 40,031,264 30,777,525 76.88 41.10 98.02 94.18
Xe-1 40,503,764 33,015,497 81.51 42.09 98.24 94.84
Xe-2 41,759,668 34,419,851 82.42 42.56 98.13 94.62
Xe-3 39,888,416 31,042,204 77.82 41.33 98.10 94.51

图3

定薯1号vs大西洋不同时期差异基因比较 A: 差异基因火山图; B: 差异基因数目统计图; C: 差异基因韦恩图。Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。"

图4

定薯1号vs大西洋差异基因富集分析 A: 差异表达基因的GO分类; B: 前20条最显著富集的KEGG途径。"

图5

定薯1号vs大西洋不同时期差异基因KEGG富集分析 Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。"

图6

不同时期淀粉和蔗糖代谢通路相关差异基因相对表达量热图 用DESeq2 log2 (定薯1号 vs 大西洋)表示各基因的相对表达量, 浅蓝色表示上调, 深蓝色表示下调。Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。"

图7

两品种不同时期差异基因表达量热图分析 用log2 (FPKM)表示基因表达量, 红色表示高表达, 蓝色表示低表达。Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。"

图8

差异基因qRT-PCR验证 A: qRT-PCR与RNA-seq结果比较, B: qRT-PCR与RNA-seq相关性分析。Da~De: 定薯1号定植后25 d、40 d、61 d、82 d和103 d; Xa~Xe: 大西洋定植后25 d、40 d、61 d、82 d和103 d。"

[1] 李扬, 王靖, 唐建昭, 张君, 胡琦, 潘志华, 潘学标. 中国马铃薯主产区生产特点、限制因子和对策分析. 中国马铃薯, 2020, 34: 374-382.
Li Y, Wang J, Tang J Z, Zhang J, Hu Q, Pan Z H, Pan X B. Analysis of production characteristics, restrictive factors, and strategies for main potato production areas in China. Chin Potato J, 2020, 34: 374-382. (in Chinese with English abstract)
[2] 田甲春, 胡新元, 田世龙, 葛霞, 李梅. 19个品种马铃薯营养成分分析. 营养学报, 2017, 39(1): 102-104.
Tian J C, Hu X Y, Tian S L, Ge X, Li M. Analysis of nutrient composition of 19 variety potatoes. Acta Nutr Sin, 2017, 39(1): 102-104. (in Chinese)
[3] 周平, 王海玲, 陆燚, 陈军, 吴显, 马杰, 杨姣, 马维. 马铃薯块茎营养品质分析鉴评. 农业科技通讯, 2021, (8): 132-136.
Zhou P, Wang H L, Lu Y, Chen J, Wu X, Ma J, Yang J, Ma W. Analysis and evaluation of nutritional quality of potato tubers. Bull Agric Sci Technol, 2021, (8): 132-136. (in Chinese)
[4] 李志新. 高淀粉马铃薯新品种克新22的选育及配套栽培技术. 黑龙江农业科学, 2011, (1): 142-143.
Li Z X. Breeding and supporting cultivation technology of new high-starch potato variety Kexin 22. Heilongjiang Agric Sci, 2011, (1): 142-143. (in Chinese)
[5] 石瑛, 张丽莉, 魏峭嵘, 秦昕. 淀粉加工型马铃薯新品种东农308的选育. 中国蔬菜, 2014, (2): 54-56.
Shi Y, Zhang L L, Wei Q R, Qin X. A new starch processing type potato variety: ‘Dongnong 308’. China Veget, 2014, (2): 54-56. (in Chinese with English abstract)
[6] 王金明, 石瑛, 梁晓丽, 赵媛媛, 黄越. 钾肥对高淀粉马铃薯块茎淀粉合成相关酶活性的影响. 作物杂志, 2016, (2): 118-123.
Wang J C, Shi Y, Liang X L, Zhao Y Y, Huang Y. Effects of potassium fertilizer on the related enzymes activity of starch synthesis in tubers of high starch potato varieties. Crops, 2016, (2): 118-123. (in Chinese with English abstract)
[7] 金光辉, 孙秀梅, 台莲梅, 姜丽丽, 杨庆东, 马力, 张志军, 李云成. 淀粉加工型马铃薯新品种‘垦薯1号’的选育. 中国马铃薯, 2014, 28: 125-126.
Jin G H, Sun X M, Tai L M, Jiang L L, Yang Q D, Ma L, Zhang Z J, Li Y C. Breeding and selection of starch processing potato variety ‘Kenshu 1’. Chin Potato J, 2014, 28: 125-126. (in Chinese with English abstract)
[8] 王腾, 马爽, 金光辉. 中国马铃薯全粉加工型品种研究进展. 中国马铃薯, 2022, 36: 266-270.
Wang T, Ma S, Jin G H. Research progress in development of potato granule and flake processing varieties in China. Chin Potato J, 2022, 36: 266-270. (in Chinese with English abstract)
[9] 唐珂, 严彩虹, 朱博, 曾子贤. 马铃薯淀粉代谢相关基因及其顺式调控元件的分子研究进展. 农业生物技术学报, 2023, 31: 833-843.
Tang K, Yan C H, Zhu B, Zeng Z X. Advances in molecular research of starch-related genes and their cis regulatory elements in potato (Solanum tuberosum). J Agric Biotechnol, 2023, 31: 833-843. (in Chinese with English abstract)
[10] Tiwari J K, Buckseth T, Challam C, Rasna Z, Nisha B, Dalamu D, Sharmistha N, Anuj K P, Rajesh K S, Satish K L, Vinod K, Manoj K. CRISPR/Cas genome editing in potato: current status and future perspectives. Front Genet, 2022, 13: 827808.
[11] Andersson M, Turesson H, Olsson N, Fält A, Ohlsson P, Gonzalez M N, Samuelsson M, Hofvander P. Genome editing in potato via CRISPR-Cas9 ribonucleoprotein delivery. Physiol Plant, 2018, 164: 378-384.
doi: 10.1111/ppl.12731 pmid: 29572864
[12] 宋波涛, 谢从华, 柳俊. 马铃薯sAGP基因表达对块茎淀粉和还原糖含量的影响. 中国农业科学, 2005, 38: 1439-1446.
Song B T, Xie C H, Liu J. Expression of Potato sAGP gene and its effects on contents of starch and reducing sugar of transgenic potato tubers. Sci Agric Sin, 2005, 38: 1439-1446. (in Chinese with English abstract)
[13] Fernie A R, Swiedrych A, Tauberger E, Lytovchenko A, Trethewey R N, Willmitzer L. Potato plants exhibiting combined antisense repression of cytosolic and plastidial isoforms of phosphoglucomutase surprisingly approximate wild type with respectto the rate of starch synthesis. Plant Physiol Biochem, 2002, 40: 921-927.
[14] McGettigan P A. Transcriptomics in the RNA-seq era. Curr Opin Chem Biol, 2013, 17: 4-11.
doi: 10.1016/j.cbpa.2012.12.008 pmid: 23290152
[15] 刘宇曦. 马铃薯块茎淀粉含量分子标记开发与应用. 东北农业大学硕士学位论文,黑龙江哈尔滨, 2023.
Liu Y X. Development and Application of Molecular Markers for Starch Content in Potato Tubers. MS Thesis of Northeast Agricultural University, Harbin, Heilongjiang, China, 2023. (in Chinese with English abstract)
[16] Li J, Baroja-Fernández E, Bahaji A, Muñoz F J, Ovecka M, Montero M, Sesma M T, Alonso-Casajús N, Almagro G, Sánchez-López A M, Hidalgo M, Zamarbide M, Pozueta-Romero J. Enhancing sucrose synthase activity results in increased levels of starch and ADP-glucose in maize (Zea mays L.) Plant Cell Physiol, 2013, 54: 282-294.
[17] 邓英毅, 郑虚. 作物叶中蔗糖磷酸合成酶的生物学功能与调控的研究进展. 广西科学院学报, 2009, 25(1): 65-71.
Deng Y Y, Zheng X. Advance on biological function and control of sucrose phosphate synthase in crop leaf. J Guangxi Acad Sci, 2009, 25(1): 65-71. (in Chinese with English abstract)
[18] Huber S C, Huber J L. Role and regulation of sucrose-phosphate synthase in higher plants. Annu Rev Plant Biol, 1996, 47: 431-444.
[19] Sun J, Zhang J, Larue C T, Huber S C. Decrease in leaf sucrose synthesis leads to increased leaf starch turnover and decreased RuBP regeneration-limited photosynthesis but not Rubisco- limited photosynthesis in Arabidopsis null mutants of SPSA1. Plant Cell Environ, 2011 34: 592-604.
[20] 叶香媛, 周文彬. 植物果糖激酶研究进展. 科学通报, 2021, 66: 2820-2831.
Ye X Y, Zhou W B. Research advances in plant fructokinases. Chin Sci Bull, 2021, 66: 2820-2831 (in Chinese with English abstract).
[21] Geng L, He X Y, Ye L Z, Zhang G P. Identification of the genes associated with β-glucan synthesis and accumulation during grain development in barley. Food Chem-mol Sci, 2022, 5: 100136.
[22] Burton R A, Collins H M, Kibble N A J, Smith J A, Shirley N J, Jobling S A, Henderson M, Singh R R, Pettolino P, Wilson S M, Bird A R, Topping D L, Bacic A, Fincher G B. Over-expression of specific HvCslF cellulose synthase-likegenes in transgenic barley increases the levels of cell wall (1,3;1,4)-beta-D-glucans and alters their fine structure. Plant Biotechnol J, 2011, 9: 117-135.
[23] Tetlow I J, Morell M K, Emes M J. Recent developments in understanding the regulation of starch metabolism in higher plants. J Exp Bot, 2004, 55: 2131-2145.
doi: 10.1093/jxb/erh248 pmid: 15361536
[24] Kim Y S, Sohn H, Jin U H, Suh S J, Lee S C, Jeon J H, Lee D S, Kim C H, Ko J H. Molecular cloning and analysis of the Thermus caldophilus ADP-glucose pyrophosphorylase. Enzyme Microb Technol, 2007, 41: 423-431.
[25] 贾小霞, 李建武, 齐恩芳, 文国宏, 李高峰, 吕和平, 马胜, 刘石, 黄伟, 张荣. ‘陇薯8号’马铃薯块茎淀粉积累特性及淀粉-蔗糖代谢途径转录组分析. 中国农业大学学报, 2023, 28(2): 23-34.
Jia X X, Li J B, Qi E F, Wen G H, Li G F, Lyu H P, Ma S, Liu S, Huang W, Zhang R. Starch accumulation pattern and transcriptome analysis of starch-sucrose metabolic pathway in potato ‘Longshu 8’ tuber. J China Agric Univ, 2023, 28(2): 23-34. (in Chinese with English abstract)
[26] Yokobayashi K, Misaki A, Harada T. Purification and properties of Pseudomonas isoamylase. Biochim Biophys Acta, 1970, 212: 458-469.
pmid: 5456995
[1] 陈国欢, 张锐, 李艳迪, 赵佳琪, 任湧涛, 张天赐, 郭华春, 李俊, 杨芳. 外源硒叶面喷施对浅紫色马铃薯块茎花青素合成的影响[J]. 作物学报, 2026, 52(6): 1876-1890.
[2] 王文辕, 燕雪嘉, 刘玉霖, 孙晓彤, 李亚楠, 唐鑫华, 石瑛. 耐弱光马铃薯品种筛选及转录因子编码基因StPIF3的克隆与功能分析[J]. 作物学报, 2026, 52(6): 1631-1645.
[3] 邹仪妹, 徐敏, 汪海洋, 姚辉, 王加峰, 刘浩, 任代胜. 两系不育系水稻幼苗根系响应盐胁迫的转录因子调控网络鉴定[J]. 作物学报, 2026, 52(6): 1728-1742.
[4] 张曦, 王广恩, 李邵琦, 刘祎, 李俊兰, 钱玉源. 基于转录组测序解析陆海杂交姊妹系马克隆值差异的形成机制[J]. 作物学报, 2026, 52(5): 1442-1458.
[5] 田立涛, 丁宁, 王树林, 齐恩芳, 张荣, 王蕊蕊, 马利雯, 李建武, 杨江伟. 马铃薯Argonaute基因家族鉴定及晚疫病诱导表达分析[J]. 作物学报, 2026, 52(4): 1116-1126.
[6] 张宇, 刘芳, 蔡诚诚, 杨小华, 吉阿么石扎, 杨元军, 王西瑶. 溴乙烷与赤霉素协同处理破除马铃薯块茎休眠的机理初探[J]. 作物学报, 2026, 52(3): 825-838.
[7] 汪玲, 胡好, 宋家凤, 程洁蓝, 陈颖, 郑婷婷, 吕钊彦, 朱晓彪, 侯华兰. 马铃薯UDP-糖基转移酶基因StUGT52的克隆及功能验证[J]. 作物学报, 2026, 52(3): 665-676.
[8] 马亮, 马璐, 张舒钰, 章慧敏, 王仁明, 宋旭东, 张振良, 冒宇翔, 陆虎华, 陈国清, 郝德荣, 周广飞. 玉米苞叶数目转录组分析及候选基因鉴定[J]. 作物学报, 2026, 52(3): 790-801.
[9] 徐强, 谢奎忠, 胡新元, 岳云, 董博, 罗爱花. 连作对马铃薯根际土壤线虫群落结构与功能的影响[J]. 作物学报, 2026, 52(2): 527-538.
[10] 杨飚, 杜帅康, 张继旺, 石瑛, 张丽莉. 马铃薯III类POD基因家族的全基因组鉴定及其表达谱分析[J]. 作物学报, 2026, 52(2): 405-420.
[11] 姬炫彤, 卞春松, 金黎平, 李森, 秦军红, 李广存. 不同耐旱型马铃薯根际微生物对干旱的响应[J]. 作物学报, 2026, 52(1): 165-177.
[12] 田甲春, 葛霞, 李守强, 李梅, 田世龙, 张亚倩, 程建新, 李玉梅. 低O2高CO2贮藏环境延缓马铃薯块茎衰老的作用机制[J]. 作物学报, 2026, 52(1): 262-278.
[13] 王雅致, 杨飚, 季香林, 石瑛, 张丽莉. 二倍体马铃薯抗旱资源鉴定及抗旱基因初步筛选[J]. 作物学报, 2026, 52(1): 72-84.
[14] 卓峰琦, 唐振三, 雷雨俊, 程李香, 赵甜甜, 吕汰, 杨晨, 张峰. 基于烹饪方式及回生温度筛选低升糖马铃薯品种(系)[J]. 作物学报, 2025, 51(9): 2538-2546.
[15] 朱锦程, 杨秋华, 程李香, 李文丽, 石明明, 李惠霞, 张峰. 马铃薯抗南方根结线虫种质资源筛选及相关生理反应分析[J]. 作物学报, 2025, 51(9): 2307-2317.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
京ICP备15008789号-2