Research progress in lotus (Nelumbo) genetics and molecular biology
-
摘要:
莲(Nelumbo)是我国特色的水生经济作物,具有很高的观赏、食用、药用等价值。近年来,在莲基因组学、遗传学和分子生物学等研究领域取得一系列重要进展,绘制了莲属两个种的参考基因组图谱,构建了莲高密度遗传连锁图谱,开展了莲观赏性状(花色、开花时间和花型)、食用品质(莲子和莲藕营养成分)、药用成分(生物碱和类黄酮)和逆境胁迫(水淹、低温、弱光、重金属和盐碱)响应的分子机制研究。本文从莲基因组测序、遗传图谱构建、重要性状基因挖掘和功能分析等方面进行了系统的综述,并对今后研究的重点和发展方向进行了展望,以期为莲的基础研究和遗传改良提供参考。
Abstract:Lotus (Nelumbo), a significant aquatic crop in China, possesses considerable ornamental, edible, and medicinal value. Over the past decade, significant progress has been achieved in the domains of genomics, genetics, and molecular biology pertaining to lotus. The reference genomes of two Nelumbo species have been constructed, along with high-density genetic linkage maps. Furthermore, various studies have explored the molecular mechanisms of ornamental traits (flower color, flowering time, and flower type), food quality (seed and rhizome nutrients), medicinal components (alkaloids and flavonoids), and stress responses (flooding, low temperature, low light, heavy metals, and salinity) of lotus. In this paper, the genome sequencing, genetic mapping, gene mining, and functional analysis of important traits in lotus are systematically reviewed. Potential goals and future directions are also discussed to provide a reference for research and genetic improvement of lotus.
-
Keywords:
- Nelumbo /
- Genome sequence /
- Important traits /
- Gene function
-
生物多样性作为地球上生命体的重要特征之一,是人类社会赖以生存和发展的物质基础[1, 2]。全球范围内的人类活动与快速的气候变化已经造成了生态系统失衡和持续的物种灭绝,导致生物多样性遭受严重的威胁[3, 4]。遗传多样性是生物多样性的重要组成部分,评估遗传多样性已经逐渐成为生态系统稳定性、自然资源与濒危物种保护和利用研究中必不可少的环节之一[5, 6]。因此,保护生物多样性尤其是遗传多样性事关人类福祉,是人类社会可持续发展的必然要求[7]。在此背景下,指导生物多样性保护和维持物种进化潜力的学科—保护生物学(Conservation biology)应运而生[8-10]。随着全球生物多样性面临的危机越来越严峻,该学科已成为生命科学领域中的一个重要分支。
保护生物学旨在研究人类活动或其他因素对包括物种、群落和生态系统在内的全球生物多样性的潜在影响[8]。早期主要从遗传学角度探究物种亲缘关系、近交繁殖、遗传漂变、基因流、种群动态、遗传结构、遗传多样性等方面[11, 12]。经过数十年分子生物学技术的发展,保护生物学从最初利用形态学、细胞学、同工酶等评估多样性水平,发展到如今利用DNA分子标记鉴定遗传多样性[13]。检测遗传多样性常用的分子标记包括限制性内切酶片段长度多态性(Restriction fragment length polymorphism,RFLP)标记、扩增片段长度多态性标记(Amplified fragment length polymorphism,AFLP),以聚合酶链式反应(Polymerase chain reaction, PCR)技术为核心的随机扩增多态性DNA标记(Random amplified polymorphic DNA, RAPD)、单引物扩增反应(Single primer amplification reaction,SPAR)、简单重复序列(Simple sequence repeat,SSR)以及简单序列重复间区的DNA标记(Inter simple sequence repeat,ISSR)技术等[14]。虽然分子标记技术推动了保护生物学的发展,但由于传统的分子标记数目较少,揭示的遗传变异仅代表物种全基因组水平的极少一部分。此外,常用的分子标记多为中性位点,无法解析自然选择和局域适应驱动的适应性遗传变异的时空格局与机理[15-17]。
下一代测序技术(Next-generation sequencing,NGS)的发展和成熟推动了如简化基因组测序(Restriction-site associated DNA sequencing, RAD)、基因分型测序(Genotyping-by-sequencing,GBS)、全基因组重测序(Whole genome resequencing)等技术的出现[18]。研究人员通过数以万计的单核苷酸多态性(SNP)位点,可检测物种不同种群在基因组水平上的遗传变异,对深入解析物种适应性进化具有重要意义[19]。同时,随着分子生物学和基因组学的发展,能够解析种群遗传变异空间分布特征的景观基因组学得以快速发展,推动适应性进化的分子机理研究跨入了一个新的时代[20-22]。景观基因组学依靠基因型-环境互作关系(Genotype-environment associations,GEAs),通过整合遗传变异和生物的空间分布模型,可从基因组水平研究物种的适应性遗传变异与景观特征之间的相互作用,并可与驱动物种自然选择的重要环境因子进行关联,解析环境因子在塑造物种遗传分布格局中的作用[20, 23]。通过鉴定基因组水平上由于环境变化引起的适应性遗传变异[20, 24],可以深入挖掘当前及未来气候变化下功能基因在适应性进化过程中发挥的作用,这为生物多样性保护带来了新的机遇[25, 26]。越来越多的研究基于景观基因组学的手段解析物种遗传多样性分布格局、遗传适应潜力以及濒危机理等科学问题,并进一步提出科学合理的保护策略[27-30]。
本文首先介绍了景观基因组学方法如何基于基因型-环境关联对物种适应性遗传变异的分布格局进行解析,接着阐述了近年来景观基因组学方法在保护生物学研究中应用的案例,最后,进一步针对景观基因组学方法在保护生物学研究中存在的问题及未来发展趋势提出了建议。
1. 景观基因组学研究方法
1.1 异常位点检测
目前,异常位点(Outlier loci)主要基于遗传分化指数(FST)和基因型-环境互作关系两类方法进行[31-34]。前者通过比较种群内或者种群间的FST,筛选超出中性选择预期的较高或较低的遗传分化位点。这些位点即为可能受到适应性进化影响的异常位点[31, 35]。目前,这类方法常用的分析软件包括BAYESCAN[32]、OutFLANK[36]和PCadapt[37, 38]等。
BAYESCAN使用贝叶斯方法将FST系数分解为特定种群和特定遗传位点的成分,并估计遗传位点偏离哈迪-温伯格平衡(Hardy-Weinberg equilibrium)的后验概率。显著偏离的位点即被检测为异常位点[32]。该方法考虑有效种群大小和种群的迁移率,但对错误发现率(False discovery rate,FDR)的设置较为敏感,较低的FDR会导致假阳性。OutFLANK是另一种基于FST的方法[36]。该方法通过计算所有个体的FST分布,去除可能受多样性(Diversity selection)或平衡(Balancing selection)选择导致的异常高以及异常低的位点,得到假定的中性位点的零分布(Null distribution)模型,之后再通过该模型在所有位点中检测异常位点。该方法针对复杂进化历史的种群具有较高的检测效率。PCadapt则是一种基于非约束性排序的主成分分析方法,通过对所有可能的主成分进行碎石图(Scree plot)检测,确定可能的种群结构,并选择合适的主成分(K)值[37, 38],然后根据马氏距离(Mahalanobis distance)筛选与种群结构显著相关的异常位点。在种群中存在杂合个体或种群结构连续的情况下,相比其他两类方法,该方法具有更高的检测效率。除了上述分析方法外,ARLEQUIN[31]、FDIST[39]、FLK[40]、XTX[41]以及最新发展的全基因组分化扫描(Genome wide differentiation scan,GWDS)方法[42]也被用于筛选遗传分化异常位点。
环境因素在不同种群的适应性进化中可能发挥重要作用,上述方法在检测异常位点时并未考虑种群间的环境异质性[43]。基于GEAs的异常位点检测方法可考虑环境因素在塑造不同种群适应性进化中的作用,该方法通过解析种群遗传变异与环境梯度之间的关系,可鉴定与环境梯度存在显著关联的异常位点[19, 24]。环境变量信息(如降水、温度、太阳辐射等)可通过WorldClim气候网站(Https://worldclim.org/)进行下载,并利用R语言[44]或者地理信息系统软件ArcGIS,基于种群经纬度信息提取对应的环境变量。由于环境变量之间的共线性(Collinearity)会对结果造成干扰,在进行分析前,可利用R语言中的usdm函数包剔除共线性较高的环境变量[45]。GEAs常用的分析方法有BAYENV[41, 46]和潜在因素混合模型(Latent factor mixed modelling,LFMM)[33, 47]等。BAYENV是一种贝叶斯方法,通过中性位点建立等位基因频率矩阵作为零分布模型,并基于该模型检验所有遗传位点的等位基因频率与环境变量的相关性。该方法可排除种群历史动态对检测结果的干扰[41]。LFMM基于贝叶斯混合模型,在检测环境变量与遗传变异的相关性时,可将种群结构作为潜在因子引入模型,从而有效地估计由种群进化历史和距离隔离模式引起的随机效应[33]。在更新的2.0版本中,可基于最小二乘法对模型潜在因子进行更准确地评估,在计算速度上也有所提升[47]。除了上述环境关联方法外,BayPass[48]、冗余分析(Redundancy analysis,RDA)[49]、BayeScENV[50]和Samβada[34]等方法也可用于检测与环境关联的异常位点。
1.2 距离/环境隔离检测
距离隔离模式(Isolation by distance,IBD)反映种群扩散受地理梯度影响而导致的种群遗传分化的局部积累情况[51]。而环境隔离模式(Isolation by environment)则反映不同栖息地的环境选择压力对种群遗传分化的影响[19]。因此,距离和环境隔离模式检测可解析地理或环境梯度在塑造物种遗传变异的空间分布和种群遗传结构中发挥重要作用。
蒙特尔检测(Mantel test)和冗余分析是进行距离和环境隔离检测应用较广的两类方法[52, 53]。蒙特尔检测是基于种群间的遗传、地理和环境的分化矩阵,分析其相关性来评估种群遗传分化的驱动因素。遗传分化矩阵可利用R语言中的hierfstat函数包进行计算[54]。地理矩阵通常计算种群间的欧氏(Euclidean)距离,可基于种群经纬度信息利用geosphere函数包实现[55]。计算环境矩阵时,考虑到环境变量间的大小差异,需首先对环境变量进行中心化(Center)和标准化(Scale),并基于R语言ecodist函数包中的布雷柯蒂斯(Bray-curtis)方法计算环境矩阵[56]。由于地理和环境变量之间通常存在一定的相关性,因此,在评估距离(环境)矩阵与遗传分化矩阵的相关性时,会进行偏蒙特尔检测(Partial Mantel test)以控制环境(距离)矩阵的影响。蒙特尔检测由于涉及到模型的统计假设问题,检测效率不高,并可能会导致假阳性的发生[57]。RDA是一种基于约束排序(Constrained ordination)的多变量统计方法,可检验遗传变异与地理、气候的多元关系,与蒙特尔检测相比,假阳性率较低[49, 58]。该方法可通过R语言中的vegan函数包实现[59]。冗余分析假设响应变量(Y)和解释变量(X)之间存在线性关系,在应用时将遗传数据作为响应变量矩阵,环境数据(地理、气候)作为解释变量矩阵,并通过环境数据来解释个体之间的遗传变异。考虑到进行冗余分析时,可能存在对环境变量的过拟合(Overfitting)现象,分析前需对环境变量通过前置选择(Forward selection)进行筛选[60]。与偏蒙特尔检测相似,冗余分析通常也进行偏冗余分析(Partial RDA)来控制协变量对遗传分化的影响[58]。此外,典型相关分析(Canonical correlation analysis,CCA)也可用于估算空间和环境变量的相对贡献[61]。
1.3 遗传变异对环境梯度的响应
解析物种对环境梯度的响应,并进一步确定可用于预测生物多样性模式的重要环境变量是生态学研究和关注的热点之一[62]。在此背景下,研究人员将遗传和环境数据关联,进行空间建模框架的搭建,从而深入解析景观尺度上物种环境适应性的遗传基础[63]。基于机器学习(Machine-learning)的梯度森林(Gradient forest,GF)[64]和基于矩阵回归(Matrix regression)的广义相异建模(Generalized dissimilarity modelling,GDM)[65]是近年来应用较广的生物多样性建模技术,可鉴定个体或种群的等位基因频率沿着环境梯度的变化趋势。这两类方法在基因-环境关系的探究中适用于明显的非线性关系,在应用于大型基因组数据集时,可处理许多低频等位基因[65]。同时,两者均可将物种不同种群的适应性遗传变异在景观尺度上进行映射。最重要的是,GF和GDM模型可评估未来气候情景下物种不同种群的遗传偏移(Genetic offset),即在当前和未来气候情景下遗传组成的不匹配(Mismatch)程度。遗传偏移越高的种群预计更易受到未来气候变化的影响,种群的遗传脆弱性越大。最近的一项研究将迁移和扩散纳入这两种模型,通过预测种群的正向(Forward)、反向(Reverse)和局域(Local)遗传偏移,可评估种群在原地(in situ)留存或迁移到新的合适栖息地的可能性,新模型可更全面地解析物种应对未来气候变化的适应能力[66]。这两类方法基于个体或种群的等位基因频率与环境数据,可利用R语言中的gradientForest函数包(Http://gradientforest.r-forge.r-project.org/)与gdm函数包[67]实现。
1.4 未来气候下的遗传脆弱性
在当今气候快速变化的背景下,解析物种在未来气候下的适应性进化,是生物多样性保护的关键。因此,我们亟需基于基因型-环境关系建立物种当前气候与种群适应性遗传变异的景观模型,然后基于该模型在空间(物种分布)和时间尺度(未来气候)上进行投影,进一步预测在未来气候变化下物种不同种群的遗传脆弱性[68]。未来气候数据可通过WorldClim气候网站(Https://worldclim.org/)下载,依据温室气体排放量的不同,网站提供了4种共享社会经济路径(Shared socioeconomic pathways,SSPs)情景(SSPs 126、 245、370、585),分别代表碳排放量从低到高。此外,网站提供了未来4个时期(2021-2040、2041-2060、2061-2080、2081-2100)的气候预测。GF、GDM和非适应性风险分析(Risk of non-adaptedness,RONA)是预测未来气候情景下遗传脆弱性应用较广的方法[64, 65, 69]。GF和GDM通过评估种群的遗传偏移程度来预测物种在未来气候下的遗传脆弱性[64, 65],这两类方法在上一节已做过介绍。非适应性风险分析基于当前基因型-环境关系建立线性回归,并以此预测未来气候下的等位基因频率,当前与未来等位基因频率的差值即为RONA值,表示物种为适应未来气候所需的等位基因频率的平均改变。该值越高,表明种群不适应未来气候的风险越高,即遗传脆弱性越大。RONA分析可基于Pina-Martins等[70]开发的Python脚本实现。此外, RDA也可基于建立的基因型-环境关系进行未来气候下的景观预测[58]。
2. 景观基因组学方法在保护生物学中的应用
保护生物学是20世纪80年代兴起于美国的一门综合性学科,主要目标是评估人类活动对生物多样性的潜在影响,并提出相关保护策略[8]。基于基因型-环境关联的景观基因组学方法的出现极大推动了保护生物学的发展[20, 22]。
森林是许多陆地生态系统的主体,在保护生物多样性、维持全球碳循环平衡、应对气候变化和全球变暖方面发挥着重要作用[20, 66, 71]。森林树种的边缘种群相对分布于中心的种群来言,可能存在较低的遗传多样性水平[72]。这些种群由于缺少应对未来气候变化的遗传变异,可能具有较高的脆弱性和灭绝风险,应加强对其的保护和管理。如,Yuan等[73]基于全基因组重测序对中国广泛分布的栎属(Quercus)植物麻栎(Q. acutissima Carruthers)27个种群117个个体在未来气候下的适应潜力进行了探究,通过遗传-环境关联分析在麻栎种群中鉴定到与气候适应性相关的多个功能基因,证明了麻栎的多基因适应模式。此外,作者借助GDM分析量化了麻栎应对气候变化的遗传偏移,发现边缘种群面临更高的局部灭绝风险。该研究在基因组水平上为麻栎的环境适应潜力提供了证据,并为物种遗传多样性保护提供了理论依据。Du等[74]对我国西南部生物多样性中心的建群树种栎属川滇高山栎(Q. aquifolioides Rehder & E. H. Wilson)进行了局域适应机理及未来气候变化的适应潜力研究。作者利用65个与干旱胁迫相关的候选基因,检测到川滇高山栎分布范围内不同的遗传变异模式,发现西藏谱系的遗传变异与地理距离相关,而喜马拉雅-横断山脉/川西高原谱系的遗传变异与环境相关,并进一步利用RONA方法检测到种群分布范围最东侧的边缘种群脆弱性较高。该研究为解析树种应对气候选择压力的响应模式提供了参考。珍稀濒危植物作为生物多样性的重要组成部分,一直以来都是保护生物学研究的重点和热点之一。Cao等[75]对东亚第三纪孑遗濒危植物领春木属(Euptelea)现存两物种响应环境变化的基因组变异及未来气候下的脆弱性进行了研究,对中国分布的大果领春木(E. pleiosperma J. D. Hooker & Thomson)和日本分布的多花领春木(E. polyandra Siebold & Zucc.)共171个个体进行简化基因组测序(RAD-seq)。结果表明,中新世晚期的气候震荡和海平面升降促进了领春木属的物种形成,4个与温度相关的气候因子驱使大果领春木种群的局域适应形成,而地理隔离是日本多花领春木种群遗传分化的主要驱动因素。梯度森林进一步检测发现现存的中国中东部种群最易受未来气候变化影响,种群脆弱性较高。该研究促进了人们对东亚孑遗森林树种的物种形成和适应性进化的认识,表明在全球气候变暖背景下我国东部地区孑遗森林树种的保育亟待加强,并为其他孑遗植物类群适应性潜力和灭绝风险区域的评估提供了方法学指导。借助景观基因组学方法,在草本和灌木类群的保护研究中也取得了一定进展。譬如,借助GEAs方法,研究人员在草本植物崖爬藤属(Tetrastigma)[76]、报春花属(Primula)[77]以及灌木植物金菀木属(Ericameria)[78]等物种中鉴定了与环境关联的适应性遗传位点,评估了植物类群在不同气候梯度下的基因组变异以及基因组对未来气候变化的脆弱性,阐明了这些类群的气候适应性遗传机理,为保护管理工作提供了重要基础。
除植物外,基于景观基因组学方法,在动物保护生物学方面也取得了一定的进展。Maier等[79]利用景观基因组学方法对加利福尼亚州内华达山脉的特有蟾蜍属物种进行了气候变化下适应潜力的预测。作者基于双酶切测序(Double-digest,dd-RAD),通过RDA和贝叶斯关联方法筛选到24个响应气候选择的候选位点,候选位点相关基因参与了对环境变化的响应。接着通过预测在未来气候下的选择压力,并利用基因型-环境关联来估计种群的适应性,确定了蟾蜍分布在北部、东部以及西南低海拔区域等3个适应未来气候的进化单元(Evolutionary units)。该研究为利用景观基因组学方法保护濒危物种蟾蜍提供了一种全面的策略,在其他濒危物种中也具有广泛的适用性。Bay等[80]基于RAD测序研究了对北美候鸟黄莺(Setophaga petechia L.)气候适应的基因组变异,发现在黄莺分布范围内存在显著的距离隔离和环境隔离效应,通过LFMM分析,还检测到降水与黄莺遗传变异的关系最为密切。此外,北美洲西部地区的种群表现出较高的基因组脆弱性,气候变化对这些地区的种群影响较大。该研究为解析鸟类及其他野生动物的基因组脆弱性提供了参考。Jaffé等[81]对南美洲分布的重要传粉者无刺蜜蜂(Melipona subnitida Ducke)的遗传结构及局域适应模式进行了评估,检测到4个遗传谱系,通过检测与温度、降水和森林覆盖度相关的基因组变异,发现这些变异在空间上呈纬度和海拔的分布模式。该研究强调了维持局域适应以及改善不同生境之间的连通性等保护行动对未来蜜蜂授粉的重要作用。气候变化、人类活动等因素对生物多样性构成了严重威胁,基于景观基因组学方法,对物种濒危机理、遗传脆弱性及局域适应等进行研究,可为基因组脆弱区解析、脆弱性种群评估、遗传拯救方案(辅助基因流、辅助迁移等)制定、种质资源储备及保护单元划分等提供关键信息,从而推动保护生物学的快速发展[82-84]。
3. 问题与展望
测序技术的发展和基因组学时代的来临,极大地推动了景观基因组学方法在保护生物学研究中的应用。然而,在应用过程中也存在一些问题。首先,在应用景观基因组学方法解决某一类问题时,有多种方法(如上述异常位点筛选方法)可采纳[22, 85]。由于不同方法基于不同的假设,会存在不同的假阳性现象。因此,通常需综合利用多种方法以获得可靠的结果。其次,以往的景观基因组学研究大多关注基因分化,而较少关注不同种群间的表型分化[86]。未来需要加强同质园和交互移植实验表型数据的获得,并通过基因型-表型关联分析鉴定表型变异的遗传基础,为生物多样性保护和濒危物种的种群复壮提供理论基础。第三,目前多从单一物种角度研究气候变化对生物多样性的影响[19],但同一区域不同物种可能采取不同的气候响应策略,为了探究气候变化对群落生物多样性的影响,更好地指导群落恢复和生境重建工作,未来需利用景观群落基因组学(Landscape community genomics)方法,在群落水平上进行多物种气候响应模式的比较研究[87]。最后,解析种群遗传变异的成因(过去)、遗传多样性的分布格局(当前)及响应未来气候变化的分子机理(未来),可全面了解种群的演化动态,是科学合理地进行生物多样性评估、保护和恢复的基础。景观基因组学研究聚焦于解析物种遗传多样性分布格局,以及响应未来气候变化的分子机理,往往忽略了结合种群遗传变异的历史成因。谱系地理学(Phylogeography)利用基因谱系关系及其时空分布来追溯种群的进化历史,可深入理解种群地理分布格局以及遗传变异的历史成因[88]。只有在全面了解物种进化历史的基础上,对遗传变异与重要环境因子进行关联,才能对物种当前遗传多样性分布格局及未来气候响应机理进行科学的探讨。目前,借助谱系地理学在物种进化历史方面已做了大量的研究和综述,系统阐述了相关类群的谱系地理式样及影响其种群分化形成的地质历史与气候因素[89, 90]。因此,后续工作可整合景观基因组学和谱系地理学,从种群进化历史、适应性遗传变异分布格局及未来气候响应机理方面入手,综合全面地开展保护生物学相关研究。
总之,未来在基于景观基因组学方法进行保护生物学研究时,首先需了解物种的进化历史(有效种群大小、谱系分化、基因渐渗等),在此基础上,通过综合比较,选择合适的方法在群落水平上解析同一区域或同一属内多个物种的环境响应机理,挖掘机理背后共同的进化动态。同时,应注意物种不同种群间表型数据的获取,以解析其表型变异的遗传基础。最后,进一步预测物种在未来气候情景下的遗传脆弱性和适应风险,评估不同物种的适应潜力,并以此来开展保护单元划分、遗传拯救方案制定、优先保护种群等保护行动。近十年来 ,高通量测序(High-throughput sequencing)技术得到了蓬勃发展,但简化基因组以及基因分型测序由于本质上仍然是基因型数据,缺少物种全基因组信息,且数据可靠性依然存在一定的问题[91]。而全基因组重测序作为目前分辨率最高的基因组学方法,为保护生物学研究提供了充分保障[18]。未来随着测序技术成本的不断降低,必将有越来越多的物种获得全基因组信息,保护生物学研究也将全面进入基因组时代。目前正处于第六次大规模物种灭绝中,全球生物多样性面临严重的危机[92]。通过全面解析物种过去、当前、未来的种群演化动态,并基于此提出生物多样性保护策略,可为保护生物学研究带来新的机遇。
-
表 1 已绘制的莲基因组图谱
Table 1 Constructed genome of Nelumbo
表 2 已构建的莲遗传连锁图谱
Table 2 Constructed linkage map of Nelumbo
标记类型
Marker type亲本
Parental material图谱长度
Length of map / cM标记个数
No. of markers连锁群数
No. of linkage groups平均图距
Average length / cM参考文献
ReferenceSSR、SRAP 中国莲与美洲黄莲 366.67 47 7 7.78 [12] SSR、SRAP 中国莲与美洲黄莲 524.51 177 11 2.96 [12] SSR、SNP 中国莲与美洲黄莲 656.90 96 8 6.60 [13] SNP 中国莲与美洲黄莲 494.30 4098 9 0.70 [13] SNP 中国莲 581.30 791 bin(8 971 SNP) 8 0.70 [15] SSR 中国莲与美洲黄莲 420.70 88 8 4.80 [14] SNP 中国莲和泰国莲 789.54 2371 bin(217 577 SNP) 8 0.33 [9] SNP 中国莲 772.92 1939 bin (79 596 SNP) 8 0.49 [16] SNP 中国莲 621.40 1475 bin(12 113 SNP) 8 0.42 [17] SNP 温带莲和热带莲 896.11 2935 bin(236 840 SNP) 8 0.27 [18] -
[1] Li HT,Luo Y,Gan L,Ma PF,Gao LM,et al. Plastid phylogenomic insights into relationships of all flowering plant families[J]. BMC Biol,2021,19:232. doi: 10.1186/s12915-021-01166-2
[2] 张行言. 中国荷花新品种图志[M]. 陈龙清, 译. 北京: 中国林业出版社, 2011: 1-100. [3] 中国科学院武汉植物研究所. 中国莲[M]. 北京: 科学出版社, 1987: 1-100. [4] Diao Y,Chen L,Yang GX,Zhou MQ,Song YC,et al. Nuclear DNA C-values in 12 species in Nymphaeales[J]. Caryologia,2006,59 (1):25−30. doi: 10.1080/00087114.2006.10797894
[5] Ming R,VanBuren R,Liu YL,Yang M,Han YP,et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn. )[J]. Genome Biol,2013,14 (5):R41. doi: 10.1186/gb-2013-14-5-r41
[6] Meng Z,Hu XX,Zhang ZL,Li ZJ,Lin QF,et al. Chromosome nomenclature and cytological characterization of sacred lotus[J]. Cytogenet Genome Res,2017,153 (4):223−231. doi: 10.1159/000486777
[7] Wang Y,Fan GY,Liu YM,Sun FM,Shi CC,et al. The sacred lotus genome provides insights into the evolution of flowering plants[J]. Plant J,2013,76 (4):557−567. doi: 10.1111/tpj.12313
[8] Shi T,Rahmani RS,Gugger PF,Wang MH,Li H,et al. Distinct expression and methylation patterns for genes with different fates following a single whole-genome duplication in flowering plants[J]. Mol Biol Evol,2020,37 (8):2394−2413. doi: 10.1093/molbev/msaa105
[9] Gui ST,Peng J,Wang XL,Wu ZH,Cao R,et al. Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements[J]. Plant J,2018,94 (4):721−734. doi: 10.1111/tpj.13894
[10] Zheng XW,Wang T,Cheng T,Zhao LL,Zheng XF,et al. Genomic variation reveals demographic history and biological adaptation of the ancient relictual,lotus(Nelumbo Adans. )[J]. Hortic Res,2022,9:29.
[11] Zheng P,Sun H,Liu J,Lin JS,Zhang XT,et al. Comparative analyses of American and Asian lotus genomes reveal insights into petal color,carpel thermogenesis and domestication[J]. Plant J,2022,110 (5):1498−1515. doi: 10.1111/tpj.15753
[12] Yang M,Han YN,Vanburen R,Ming R,Xu LM,et al. Genetic linkage maps for Asian and American lotus constructed using novel SSR markers derived from the genome of sequenced cultivar[J]. BMC Genomics,2012,13:653. doi: 10.1186/1471-2164-13-653
[13] Zhang Q,Li LT,Vanburen RL,Liu YL,Yang M,et al. Optimization of linkage mapping strategy and construction of a high-density American lotus linkage map[J]. BMC Genomics,2014,15:372. doi: 10.1186/1471-2164-15-372
[14] 陈岳,张微微,莫海波,付乃峰,田代科. EST-SSR标记构建莲(Nelumbo Adans. )遗传连锁图谱[J]. 分子植物育种,2017,15(6):2265−2273. Chen Y,Zhang WW,Mo HB,Fu NF,Tian DK. Construction of genetic linkage map of lotus (Nelumbo Adans. ) using EST-SSR markers[J]. Molecular Plant Breeding,2017,15 (6):2265−2273.
[15] Liu ZW,Zhu HL,Liu YP,Kuang J,Zhou K,et al. Construction of a high-density,high-quality genetic map of cultivated lotus (Nelumbo nucifera) using next-generation sequencing[J]. BMC Genomics,2016,17:466. doi: 10.1186/s12864-016-2781-4
[16] 严寒松. 中国莲高密度遗传图谱的构建及花瓣数控制基因的定位[D]. 福州: 福建农林大学, 2019: 1-10. [17] Huang LY,Li M,Cao DD,Yang PF. Genetic dissection of rhizome yield-related traits in Nelumbo nucifera through genetic linkage map construction and QTL mapping[J]. Plant Physiol Biochem,2021,160:155−165. doi: 10.1016/j.plaphy.2021.01.020
[18] Liu YL,Song HY,Zhang MH,Yang D,Deng XB,et al. Identification of QTLs and a putative candidate gene involved in rhizome enlargement of Asian lotus (Nelumbo nucifera)[J]. Plant Mol Biol,2022,110 (1-2):23−36. doi: 10.1007/s11103-022-01281-w
[19] Wang YJ,Chen YQ,Yuan M,Xue ZY,Jin QJ,Xu YC. Flower color diversity revealed by differential expression of flavonoid biosynthetic genes in sacred lotus[J]. J Am Soc Hortic Sci,2016,141 (6):573−582. doi: 10.21273/JASHS03848-16
[20] Zhu HH,Yang JX,Xiao CH,Mao TY,Zhang J,Zhang HY. Differences in flavonoid pathway metabolites and transcripts affect yellow petal colouration in the aquatic plant Nelumbo nucifera[J]. BMC Plant Biol,2019,19 (1):277. doi: 10.1186/s12870-019-1886-8
[21] Deng J,Li M,Huang LY,Yang M,Yang PF. Genome-Wide analysis of the R2R3 MYB subfamily genes in lotus (Nelumbo nucifera)[J]. Plant Mol Biol Rep,2016,34 (5):1016−1026. doi: 10.1007/s11105-016-0981-3
[22] Sun SS,Gugger PF,Wang QF,Chen JM. Identification of a R2R3-MYB gene regulating anthocyanin biosynthesis and relationships between its variation and flower color difference in lotus (Nelumbo Adans. )[J]. PeerJ,2016,4:e2369. doi: 10.7717/peerj.2369
[23] Deng J,Li JJ,Su MY,Lin ZY,Chen L,Yang PF. A bHLH gene NnTT8 of Nelumbo nucifera regulates anthocyanin biosynthesis[J]. Plant Physiol Biochem,2021,158:518−523. doi: 10.1016/j.plaphy.2020.11.038
[24] 李静,刘正位,朱红莲,匡晶,彭静,等. 莲花色性状QTL定位及分子标记开发[J]. 植物科学学报,2022,40(6):820−828. doi: 10.11913/PSJ.2095-0837.2022.60820 Li J,Liu ZW,Zhu HL,Kuang J,Peng J,et al. Mapping of quantitative trait loci and development of closely linked markers of flower color in lotus (Nelumbo nucifera Gaertner)[J]. Plant Science Journal,2022,40 (6):820−828. doi: 10.11913/PSJ.2095-0837.2022.60820
[25] Liu J, Wang YX, Deng XB, Zhang MH, Sun H, et al. Transcription factor NnMYB5 controls petal color by regulating GLUTATHIONE S-TRANSFERASE2 in Nelumbo nucifera[J]. Plant Physiol, 2023: kiad363. doi: 10.1093/plphys/kiad363.
[26] 杨小凤,李小蒙,廖万金. 植物开花时间的遗传调控通路研究进展[J]. 生物多样性,2021,29(6):825−842. doi: 10.17520/biods.2020370 Yang XF,Li XM,Liao WJ. Advances in the genetic regulating pathways of plant flowering time[J]. Biodiversity Science,2021,29 (6):825−842. doi: 10.17520/biods.2020370
[27] 郭蓓. 荷花成花相关基因NnAP1的克隆与表达[D]. 郑州: 河南农业大学, 2014: 1-10. [28] Yang M,Zhu L,Xu L,Pan C,Liu Y. Comparative transcriptomic analysis of the regulation of flowering in temperate and tropical lotus (Nelumbo nucifera) by RNA-Seq[J]. Ann Appl Biol,2014,165 (1):73−95. doi: 10.1111/aab.12119
[29] 胡裕凤. 莲花期与地下茎膨大相关基因的功能分析[D]. 北京: 中国科学院大学, 2019: 1-10. [30] 宋贺云. 莲FT基因调控开花的功能研究[D]. 北京: 中国科学院大学, 2020: 1-10. [31] Zhang L,Zhang F,Liu FB,Shen J,Wang JX,et al. The lotus NnFTIP1 and NnFT1 regulate flowering time in Arabidopsis[J]. Plant Sci,2021,302:110677. doi: 10.1016/j.plantsci.2020.110677
[32] Yang M,Zhu LP,Xu LM,Liu YL. Population structure and association mapping of flower-related traits in lotus (Nelumbo Adans. ) accessions[J]. Sci Horticult,2014,175:214−222. doi: 10.1016/j.scienta.2014.06.017
[33] 李玲. 莲高密度遗传连锁图谱的构建及花期QTL定位[D]. 北京: 中国科学院大学, 2017: 1-10. [34] 黄秀,田代科,张微微,曾宋君,莫海波. 荷花“重瓣化”的花器官形态发育比较观察[J]. 植物分类与资源学报,2014,36(3):303−309. Huang X,Tian DK,Zhang WW,Zeng SJ,Mo HB. Comparison of floral organ morphological development between single and double flowers in Nelumbo nucifera[J]. Plant Diversity and Resources,2014,36 (3):303−309.
[35] Lin ZY,Damaris RN,Shi T,Li JJ,Yang PF. Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera)[J]. BMC Genomics,2018,19 (1):554. doi: 10.1186/s12864-018-4950-0
[36] Lin ZY,Cao DD,Damaris RN,Yang PF. Comparative transcriptomic analysis provides insight into carpel petaloidy in lotus (Nelumbo nucifera)[J]. PeerJ,2021,9:e12322. doi: 10.7717/peerj.12322
[37] Sun H,Li JJ,Song HY,Yang D,Deng XB,et al. Comprehensive analysis of AGPase genes uncovers their potential roles in starch biosynthesis in lotus seed[J]. BMC Plant Biol,2020,20 (1):457. doi: 10.1186/s12870-020-02666-z
[38] Chen HH,Chu P,Zhou YL,Ding Y,Li Y,et al. Ectopic expression of NnPER1,a Nelumbo nucifera 1-cysteine peroxiredoxin antioxidant,enhances seed longevity and stress tolerance in Arabidopsis[J]. Plant J,2016,88 (4):608−619. doi: 10.1111/tpj.13286
[39] Zhang D,Liu T,Sheng JY,Lv S,Ren L. TMT-Based quantitative proteomic analysis reveals the physiological regulatory networks of embryo dehydration protection in lotus (Nelumbo nucifera)[J]. Front Plant Sci,2021,12:792057. doi: 10.3389/fpls.2021.792057
[40] Chen L,Xin J,Song HY,Xu F,Yang H,et al. Genome-wide study and functional characterization elucidates the potential association of late embryogenesis abundant (LEA) genes with lotus seed development[J]. Int J Biol Macromol,2023,226:1−13. doi: 10.1016/j.ijbiomac.2022.11.301
[41] Li JJ,Shi T,Huang LY,He DL,Nyong'A TM,Yang PF. Systematic transcriptomic analysis provides insights into lotus (Nelumbo nucifera) seed development[J]. Plant Growth Regul,2018,86 (3):339−350.
[42] Wang L,Fu JL,Li M,Fragner L,Weckwerth W,Yang PF. Metabolomic and proteomic profiles reveal the dynamics of primary metabolism during seed development of lotus (Nelumbo nucifera)[J]. Front Plant Sci,2016,7:750.
[43] Gu C,Wang L,Zhang LY,Liu YL,Yang M,et al. Characterization of genes encoding granule-bound starch synthase in sacred lotus reveals phylogenetic affinity of Nelumbo to Proteales[J]. Plant Mol Biol Rep,2013,31 (5):1157−1165. doi: 10.1007/s11105-013-0605-0
[44] Zhu FL,Sun H,Diao Y,Zheng XW,Xie KQ,Hu ZL. Genetic diversity,functional properties and expression analysis of NnSBE genes involved in starch synthesis of lotus (Nelumbo nucifera Gaertn. )[J]. PeerJ,2019,7:e7750. doi: 10.7717/peerj.7750
[45] 刘正位,郭丹丹,彭静,朱红莲,匡晶,等. 莲子产量相关性状的QTL定位[J]. 园艺学报,2020,47(8):1565−1576. doi: 10.16420/j.issn.0513-353x.2019-0665 Liu ZW,Guo DD,Peng J,Zhu HL,Kuang J,et al. QTL mapping of six seed yield related traits in lotus[J]. Acta Horticulturae Sinica,2020,47 (8):1565−1576. doi: 10.16420/j.issn.0513-353x.2019-0665
[46] 韩玉燕. 莲藕NnWOX1-1、NnWOX4-3和NnWOX5-1的克隆、表达及功能初步分析[D]. 扬州: 扬州大学, 2021: 1-10. [47] Yang M,Zhu LP,Pan C,Xu LM,Liu YL,et al. Transcriptomic analysis of the regulation of rhizome formation in temperate and tropical lotus (Nelumbo nucifera)[J]. Sci Rep,2015,5:13059. doi: 10.1038/srep13059
[48] Cao DD,Lin ZY,Huang LY,Damaris RN,Li M,Yang PF. A CONSTANS-LIKE gene of Nelumbo nucifera could promote potato tuberization[J]. Planta,2021,253 (3):65. doi: 10.1007/s00425-021-03581-9
[49] 程立宝,李淑艳,李岩,尹静静,陈学好,李良俊. 莲藕根状茎膨大过程中淀粉合成相关基因的表达[J]. 中国农业科学,2012,45(16):3330−3336. doi: 10.3864/j.issn.0578-1752.2012.16.012 Cheng LB,Li SY,Li Y,Yin JJ,Chen XH,Li LJ. Expression of genes related to starch synthesis during rhizome swelling of lotus root (Nelumbo nucifera Gaertn. )[J]. Scientia Agricultura Sinica,2012,45 (16):3330−3336. doi: 10.3864/j.issn.0578-1752.2012.16.012
[50] 沈王俊. 莲藕颗粒结合淀粉合成酶基因NnGBSS调控直链淀粉合成的初步探究[D]. 扬州: 扬州大学, 2021: 1-10. [51] 张莉,印荔,杨见秋,程立宝,李良俊. 莲藕可溶性淀粉合成酶基因LrSSS的克隆与表达特性分析[J]. 园艺学报,2015,34(3):496−504. Zhang L,Yin L,Yang JQ,Cheng LB,Li LJ. Cloning and expression profiling analysis of soluble starch synthase gene in lotus rhizome[J]. Acta Horticulturae Sinica,2015,34 (3):496−504.
[52] 李婷婷,李玉蕊,鄢敏丽,张海星,王瑞红,等. 莲藕PPO基因家族成员鉴定与分析[J]. 农业生物技术学报,2022,30(1):38−49. doi: 10.3969/j.issn.1674-7968.2022.01.004 Li TT,Li YR,Yan ML,Zhang HX,Wang RH,et al. Identification and analysis of PPO gene family members in lotus (Nelumbo nucifera)[J]. Journal of Agricultural Biotechnology,2022,30 (1):38−49. doi: 10.3969/j.issn.1674-7968.2022.01.004
[53] 郝晓燕. 莲藕多酚氧化酶基因的克隆与表达分析[D]. 杨凌: 西北农林科技大学, 2011: 1-10. [54] Sharma BR,Gautam LNS,Adhikari D,Karki R. A comprehensive review on chemical profiling of Nelumbo nucifera:potential for drug development[J]. Phytother Res,2017,31 (1):3−26. doi: 10.1002/ptr.5732
[55] Menéndez-Perdomo IM,Facchini PJ. Benzylisoquinoline alkaloids biosynthesis in sacred lotus[J]. Molecules,2018,23 (11):2899. doi: 10.3390/molecules23112899
[56] Chen S,Fang LC,Xi HF,Guan L,Fang JB,et al. Simultaneous qualitative assessment and quantitative analysis of flavonoids in various tissues of lotus (Nelumbo nucifera) using high performance liquid chromatography coupled with triple quad mass spectrometry[J]. Anal Chim Acta,2012,724:127−135. doi: 10.1016/j.aca.2012.02.051
[57] Deng XB,Zhu LP,Fang T,Vimolmangkang S,Yang D,et al. Analysis of isoquinoline alkaloid composition and wound-induced variation in Nelumbo using HPLC-MS/MS[J]. J Agric Food Chem,2016,64 (5):1130−1136. doi: 10.1021/acs.jafc.5b06099
[58] Hagel JM,Facchini PJ. Benzylisoquinoline alkaloid metabolism:a century of discovery and a brave new world[J]. Plant Cell Physiol,2013,54 (5):647−672. doi: 10.1093/pcp/pct020
[59] Kashiwada Y,Aoshima A,Ikeshiro Y,Chen YP,Furukawa H,et al. Anti-HIV benzylisoquinoline alkaloids and flavonoids from the leaves of Nelumbo nucifera,and structure-activity correlations with related alkaloids[J]. Bioorg Med Chem,2005,13 (2):443−448. doi: 10.1016/j.bmc.2004.10.020
[60] Menéndez-Perdomo IM,Facchini PJ. Elucidation of the(R)-enantiospecific benzylisoquinoline alkaloid biosynthetic pathways in sacred lotus (Nelumbo nucifera)[J]. Sci Rep,2023,13:2955. doi: 10.1038/s41598-023-29415-0
[61] Deng XB,Zhao L,Fang T,Xiong YQ,Ogutu C,et al. Investigation of benzylisoquinoline alkaloid biosynthetic pathway and its transcriptional regulation in lotus[J]. Hortic Res,2018,5:29. doi: 10.1038/s41438-018-0035-0
[62] Yu YT,Liu Y,Dong GQ,Jiang JZ,Leng L,et al. Functional characterization and key residues engineering of a regiopromiscuity O-methyltransferase involved in benzylisoquinoline alkaloid biosynthesis in Nelumbo nucifera[J]. Hortic Res,2023,10 (2):uhac276. doi: 10.1093/hr/uhac276
[63] Feng CY,Li SS,Taguchi G,Wu Q,Yin DD,et al. Enzymatic basis for stepwise C-glycosylation in the formation of flavonoid di-C-glycosides in sacred lotus(Nelumbo nucifera Gaertn. )[J]. Plant J,2021,106 (2):351−365. doi: 10.1111/tpj.15168
[64] 王海峰. 不同季节长期水淹对几种陆生植物的存活、生长和恢复生长的影响[D]. 重庆: 西南大学, 2008: 1-10. [65] 李祥志,刘兆磊,陈发棣,丁跃生,高迎,王宏辉. 荷花耐深水评价体系及耐深水鉴定[J]. 安徽农业科学,2014,42(3):679−682. doi: 10.3969/j.issn.0517-6611.2014.03.011 Li XZ,Liu ZL,Chen FD,Ding YS,Gao Y,Wang HH. Study on establishment of evaluation system for deepwater tolerance and its identification of Nelumbo nucifera Gaertn.[J]. Journal of Anhui Agricultural Sciences,2014,42 (3):679−682. doi: 10.3969/j.issn.0517-6611.2014.03.011
[66] Deng XB,Yang D,Sun H,Liu J,Song HY,et al. Time-course analysis and transcriptomic identification of key response strategies of Nelumbo nucifera to complete submergence[J]. Hortic Res,2022,9:uhac001. doi: 10.1093/hr/uhac001
[67] Wang B,Jin QJ,Zhang X,Mattson NS,Ren HH,et al. Genome-wide transcriptional analysis of submerged lotus reveals cooperative regulation and gene responses[J]. Sci Rep,2018,8:9187. doi: 10.1038/s41598-018-27530-x
[68] Jin QJ,Xu YC,Mattson NS,Li X,Wang B,et al. Identification of submergence-responsive microRNAs and their targets reveals complex miRNA-mediated regulatory networks in lotus (Nelumbo nucifera Gaertn. )[J]. Front Plant Sci,2019,8:6.
[69] 赵瑞霜. 荷花抗寒品种筛选及生理研究[D]. 南京: 南京农业大学, 2014: 1-10. [70] 黄勇琴. 自然越冬过程中冬荷抗寒及休眠生理生化特性研究[D]. 郑州: 河南农业大学, 2012: 1-10. [71] Li X,Zhang D,Xu J,Jiang J,Jiang HW. The protective effect of cold acclimation on the low temperature stress of the lotus (Nelumbo nucifera)[J]. Hortic Sci,2022,49 (1):29−37. doi: 10.17221/62/2020-HORTSCI
[72] 冷寒冰,秦俊,叶康,奉树成,高凯. 不同光照环境下荷花叶片光合光响应模型比较[J]. 应用生态学报,2014,25(10):2855−2860. doi: 10.13287/j.1001-9332.20140801.010 Leng HB,Qin J,Ye K,Feng SC,Gao K. Comparison of light response models of photosynthesis in Nelumbo nucifera leaves under different light conditions[J]. Chinese Journal of Applied Ecology,2014,25 (10):2855−2860. doi: 10.13287/j.1001-9332.20140801.010
[73] 张晓,任惠惠,曹婧,苗其军,金奇江,等. 弱光导致荷花花芽败育的机制探析[J]. 南京农业大学学报,2019,42(6):1040−1049. doi: 10.7685/jnau.201904027 Zhang X,Ren HH,Cao J,Miao QJ,Jin QJ,et al. Analysis of the mechanism of lotus flower bud abortion caused by low light[J]. Journal of Nanjing Agricultural University,2019,42 (6):1040−1049. doi: 10.7685/jnau.201904027
[74] Sheng JY,Wang GD,Liu T,Xu Z,Zhang D. Comparative transcriptomic and proteomic profiling reveals molecular models of light signal regulation of shade tolerance in bowl lotus (Nelumbo nucifera)[J]. J Proteomics,2022,257:104455. doi: 10.1016/j.jprot.2021.104455
[75] 张正庆,鲍美娥,陈嘉斌,郝海艳. 植物对重金属的耐性机制[J]. 甘肃科技,2013,29(5):69−71. doi: 10.3969/j.issn.1000-0952.2013.05.024 [76] 司卫静. 荷花对铅尾矿渗出液的耐性及修复研究[D]. 南昌: 江西财经大学, 2013: 1-10. [77] 于辉,彭佳师,严明理. 镉在莲各器官中累积规律的研究[J]. 生态科学,2021,40(1):82−85. doi: 10.14108/j.cnki.1008-8873.2021.01.011 Yu H,Peng JS,Yan ML. Research on the accumulation of cadmium in lotus(Nelumbo nucifera Gaertn.) organs[J]. Ecological Science,2021,40 (1):82−85. doi: 10.14108/j.cnki.1008-8873.2021.01.011
[78] 裴康康. 镉和锌胁迫对荷花种子萌发和生理生化的影响[D]. 郑州: 河南农业大学, 2010: 1-10. [79] 王海波. 荷花在硝酸铅环境下的生理变化及对铅吸收效应的研究[D]. 郑州: 河南农业大学, 2009: 1-10. [80] 苏少文,刘莹,黄志远,吴芳芳,刘艺平. 不同荷花对盐碱胁迫的响应[J]. 北方园艺,2020(23):52−59. Su SW,Liu Y,Huang ZY,Wu FF,Liu YP. Response of different lotus flowers to saline-alkali stress[J]. Northern Horticulture,2020 (23):52−59.
[81] 刘艺平,张一琪,苏少文,刘红利,贺丹,孔德政. 混合盐碱胁迫下不同抗性荷花品种比较转录组分析[J]. 江苏农业科学,2022,50(3):1−8. doi: 10.15889/j.issn.1002-1302.2022.03.001 Liu YP,Zhang YQ,Su SW,Liu HL,He D,Kong DZ. Comparative transcriptome analysis of different resistant lotus(Nelumbo nucifera) cultivars under mixed saline-alkali stress[J]. Jiangsu Agricultural Sciences,2022,50 (3):1−8. doi: 10.15889/j.issn.1002-1302.2022.03.001
[82] Liu RJ,Shi HT,Wang YP,Chen S,Deng J,et al. Comparative physiological analysis of lotus(Nelumbo nucifera) cultivars in response to salt stress and cloning of NnCIPK genes[J]. Sci Horticult,2014,173:29−36. doi: 10.1016/j.scienta.2014.04.032
[83] Cheng LB,Yang JJ,Yin L,Hui LC,Qian HM,et al. Transcription factor NnDREB1 from lotus improved drought tolerance in transgenic Arabidopsis thaliana[J]. Biol Plant,2017,61 (4):651−658. doi: 10.1007/s10535-017-0718-7
[84] Liu XJ,Du FF,Li NW,Chang YJ,Yao DR. Gene expression profile in the long-living lotus:insights into the heat stress response mechanism[J]. PLoS One,2016,11 (3):e0152540. doi: 10.1371/journal.pone.0152540
-
期刊类型引用(0)
其他类型引用(2)