《Horticulturae》:The Development of a Core SSR Marker Set for DNA Fingerprinting and Accession Discrimination in Durian (Durio zibethinus)
编辑推荐:
榴莲(Durio zibethinus)种质的名称可能存在错误、标签误标或遗传误鉴,尤其在引种、交换和苗木流通环节。形态学鉴定还受环境条件和植物发育阶段影响。研究人员开发了用于榴莲种质鉴定的核心简单序列重复(simple sequence repeat,SSR
榴莲(Durio zibethinus)种质的名称可能存在错误、标签误标或遗传误鉴,尤其在引种、交换和苗木流通环节。形态学鉴定还受环境条件和植物发育阶段影响。研究人员开发了用于榴莲种质鉴定的核心简单序列重复(simple sequence repeat,SSR)标记集。基于榴莲参考基因组进行全基因组挖掘,共鉴定出106,776个SSR位点;从160对候选引物中经琼脂糖凝胶筛选和M13尾标荧光毛细管电泳筛选,选出18个稳定且具多态性的标记。在58份榴莲种质中,这18个SSR位点共检测到142个等位基因,每位点3–15个等位基因,多态性信息含量(polymorphism information content,PIC)值为0.422–0.849。研究人员对所有标记组合进行穷举评估,确定最小六位点核心SSR标记集Y6、Y12、Y15、Y75、Y91和Y108可为全部58份种质生成唯一多位点基因型,累积一致概率(probability of identity,PID)为7.99×10?8。STRUCTURE分析中K代表推断遗传簇数,结果支持K=8,并显示若干份种质存在混合祖源。该研究开发的核心SSR标记集可区分58份榴莲种质并用于构建其DNA指纹;在更大且遗传更多样的种质集合中验证后,方可用于常规商业化部署。
研究背景方面,榴莲(Durio zibethinus)为锦葵科热带果树,有“热带水果之王”之称,近年受中国市场拉动贸易快速扩张。榴莲品种在市场价值、果实品质和育种潜力上差异显著,但常规形态学鉴定受环境条件、栽培措施和发育阶段影响;引种、交换与苗木流通中还会出现同物异名、同名异物和标签错误。苗木未到生殖成熟前又难以凭果实用性状可靠鉴别,因此亟需可重复、低成本的分子鉴定体系。已有RAPD、ISSR和SSR用于遗传多样性研究,但系统开发核心SSR标记集、标准化分型流程和DNA指纹仍有限。该论文发表于《Horticulturae》,研究人员旨在利用毛细管电泳开发基因组来源SSR标记,筛选能区分58份种质的紧凑核心SSR标记集,构建DNA指纹与分子身份码,并分析遗传多样性和群体结构。
主要关键技术方法方面,样本为海南三亚某榴莲种质资源圃58份种质,其中52份记起来源于马来西亚、5份泰国、1份越南。研究人员从榴莲参考基因组用MISA平台鉴定SSR并设计引物;先以琼脂糖凝胶初筛扩增特异性和稳定性,再以M13尾标荧光策略结合毛细管电泳(capillary electrophoresis,CE)对全部种质分型;用穷举组合搜索最少位点核心标记集,并按等位基因频率计算单位点与多位点PID;群体分析采用PopGene32、NTSYS-pc、GenAlEx和STRUCTURE等软件进行遗传多样性、UPGMA聚类、分子方差分析(analysis of molecular variance,AMOVA)和贝叶斯聚类。
3.1 全基因组SSR鉴定与引物开发,研究人员从735.19 Mb榴莲参考基因组鉴定出106,776个SSR位点,密度145.24个/Mb,二核苷酸重复最多占54.66%,三核苷酸占29.53%,还检出12,940个复合SSR;据此设计160对候选引物。
3.2 SSR标记筛选与多态性评估,160对引物中56对初筛稳定特异,再选18对做毛细管电泳;18个位点共142个等位基因,观察等位基因数(Na)3–15,有效等位基因数(Ne)2.13–7.26,Nei基因多样性0.5299–0.8622,香农指数0.8179–2.2297,PIC 0.4222–0.8485,15个位点PIC大于0.50;5份种质独立提DNA重测显示核心位点结果100%一致。
3.3 核心SSR标记集判别能力,研究人员穷举18个位点全部262,143种组合,1–5个位点最佳组合最多区分56份,不能全部分离;六位点Y6+Y12+Y15+Y75+Y91+Y108产生58个唯一基因型,判别率100%,累积PID为7.99×10?8;删除任一核心位点后唯一判别数降为51–56份,证明六位点均不可省。
3.4 DNA指纹数据库与分子身份码构建,3.4.1 研究人员以18个位点建立58份种质基因型库,用六核心位点生成黑白方格表示的二进制指纹,每份种质都有唯一六位点基因型。3.4.2 按Y6、Y12、Y15、Y75、Y91、Y108固定顺序将“小等位基因/大等位基因”连字符连接,生成58个互不相同的分子身份码。
3.5 聚类分析与分子方差分析,3.5.1 UPGMA基于18位点等位基因有无矩阵,在Dice相似系数约0.42处描述性分为6簇,大簇分别含19份和23份种质。3.5.2 AMOVA按UPGMA六簇划分时,簇间变异占2%,簇内个体间17%,个体内81%,FST为0.016;按记录来源把越南1份除外后,马来西亚与泰国两组间变异占5%,FST为0.049,说明地理来源解释力有限。
3.6 STRUCTURE推断群体结构,研究人员用18位点进行K=2至10贝叶斯混合模型分析,Evanno法ΔK在K=8最大,多数种质主要归入单簇,若干种质具两簇及以上成员,显示不同程度的混合。
3.7 主坐标分析(principal coordinate analysis,PCoA),基于18位点Dice距离,前两个坐标分别解释19.74%和13.17%,累计32.91%;按K=8主要祖源着色后,部分组间有一定分离,但II、III、VI、VII组左中部重叠明显,说明亚结构弱且无清晰边界。
讨论部分总结,研究人员指出榴莲基因组SSR密度高于部分棉属报道,但受组装质量与判定标准影响,宜描述性看待;18标记平均PIC 0.6639高于既往榴莲SSR研究,这来自候选标记按多态性预筛以及毛细管电泳比琼脂糖凝胶分辨率更高,并非创造新生物学多态性。已有榴莲SSR研究多面向多样性或小样本,本文则用全基因组挖掘加毛细管电泳验证18个多态位点并给出覆盖58份种质的唯一身份码。六核心位点分布在不同染色体且缺一则判别下降,说明信息互补;但58份中52份马来西亚、5份泰国、1份越南,地理不平衡,100%判别仅限本面板。UPGMA六簇与STRUCTURE K=8不矛盾,前者按成对相似聚类,后者允许混合祖源;低组间变异与PCoA重叠说明无强分化。核心SSR集可用于种质核查、指纹库、品种真伪和苗木追溯,但是否适用于更广种质还需大样本、地理更平衡集合验证;研究也未专门评估无效等位基因(null alleles)导致等位基因丢失和假纯合的风险。
结论部分翻译,本研究在榴莲参考基因组鉴定106,776个SSR位点,并从160个候选中筛选出18个稳定且具多态性的SSR引物对。在58份榴莲种质中,这些位点检出142个等位基因,平均PIC为0.6639,显示其对遗传多样性分析具有高信息量。穷举组合分析确定由Y6、Y12、Y15、Y75、Y91和Y108组成的核心SSR标记集,该集为每份种质生成唯一多位点基因型,并可构建DNA指纹与分子身份码。遗传分析表明多数变异存在于个体内和组内个体间,记录地理来源仅解释本面板小部分变异。这些结果为一个紧凑且具前景的标记框架提供支持,用于所分析种质面板内的榴莲种质鉴定、指纹库构建和 accession 管理,在更大型、地理更广且更均衡的种质集合验证后方可更广泛应用。