Page 29 - 《广西植物》2025年第4期
P. 29

4 期                      徐梦阳等: 南酸枣转录组特征分析及 SSR 标记开发                                         6 4 3

            价、品种鉴定、遗传图谱构建、QTL 定位、濒危物种                          物ꎬ以期为其遗传多样性评价、雌雄性别标记开发
            保 护 等 领 域 ( Kuroda et al.ꎬ 2006ꎻ 阎 毛 毛 等ꎬ          及指纹图谱构建等研究奠定基础ꎮ
            2011)ꎮ 通过不同分子标记检测物种多态性位点
            时ꎬSSR 位点所提供的多态性通常远高于其他分                            1  材料与方法
            子标记ꎬ可以更好地反映种群间和种群内的遗传
            分化和亚分化ꎬ提供更准确的物种进化信息ꎮ 与                             1.1 试验材料
            其他分子标记相比ꎬSSR 位点通常提供更高的多                                转录组测序材料采自江西农业大学南酸枣种

            态性ꎬ更准确地反映种群间和种群内的遗传分化ꎮ                             质资源 圃ꎬ 选 取 10 年 生 ( 雌 雄 无 性 系 各 3 株:
            例如ꎬ郭晋宏等(2019) 通过 SSR 分析山西晋南漆                       Female ̄1、Female ̄2、Female ̄3、Male ̄1、Male ̄2、Male ̄
            树(Toxicodendron vernicifluum)群体ꎬ发现其群体内             3)南酸枣新鲜叶片经液氮速冻后ꎬ干冰保存送至

            遗传多样性丰富ꎬ遗传变异主要集中在群体内部ꎻ                             检测 公 司 进 行 转 录 组 测 序ꎮ 采 集 江 西 崇 义
            任重等(2022)基于 SSR 分子标记研究了中国黄连                        (JXCY)、贵州江口( GZJK)、湖南桑植( HNSZ)、云
            木( Pistacia chinensis) 的遗传多样性ꎬ发现其具有                南芒 市 ( YNMS)、 广 西 平 南 ( GXPN)、 湖 北 赤 壁
            中等偏下的遗传分化ꎬ遗传分化主要集中在群体                              (HBCB )、 安 徽 祁 门 ( AHQM ) 和 重 庆 沙 坪 坝
            内部ꎬ划分为两个主要群体ꎮ 南酸枣的分子生物                             (CQSPB)8 个群体各 1 个样品作为引物筛选样品ꎮ
            研究因起步较晚而缺乏基因组和相关序列信息ꎮ                              引物多态性验证所用试验材料取自 3 个省(区) 的
            本研究利用南酸枣雌雄叶片的转录组数据ꎬ分析                              85 份样品(表 1)ꎮ 所有南酸枣样本采其健康无病
            南酸枣雌雄表达差异ꎬ以及其转录组 SSR 位点的                           虫害嫩叶ꎬ经过硅胶迅速干燥之后保存用于 DNA
            分布及序列特征ꎬ并在此基础上开发 SSR 标记引                           提取ꎬ最后进行引物筛选及多态性验证ꎮ


                                               表 1  9 个南酸枣群体采样信息
                              Table 1  Sampling information of nine Choerospondias axillaris populations
                  群体编号                      采样地                    样品数量             经度              纬度
                Population code          Sampling location      Number of samples  Longitude       Latitude
                   HNTY              湖南桃源 Taoyuanꎬ Hunan              10         111°30′00″ E    28°54′36″ N
                   HNYS              湖南永顺 Yongshunꎬ Hunan             10         109°51′36″ E    28°59′24″ N

                   HNXT              湖南湘潭 Xiangtanꎬ Hunan             7          112°57′36″ E    27°46′48″ N
                   HBCY              湖北崇阳 Chongyangꎬ Hubei            10         114°03′00″ E    29°33′36″ N
                   HBXA              湖北咸安 Xian’anꎬ Hubei              10         114°16′48″ E    29°52′12″ N

                   HBCB                湖北赤壁 Chibiꎬ Hubei              9          113°54′36″ E    29°43′48″ N
                   GXZS             广西钟山 Zhongshanꎬ Guangxi           9          111°18′36″ E    24°31′48″ N
                   GXCW              广西苍梧 Cangwuꎬ Guangxi             8          111°33′36″ E    23°51′00″ N

                   GXPN              广西平南 Pingnanꎬ Guangxi            12         110°24′00″ E    23°33′00″ N


            1.2 试验方法                                           的有效序列 ( clean reads) 与参考序列比对ꎬ 得到
            1.2.1 转录组组装及分析  对高通量测序获得的原                         read countꎮ 将 read count 转换为 FPKM (Fragments
            始数据(raw data) 进行过滤ꎬ获得高质量的有效数                       Per Kilobase Million)ꎬ得到基因的表达水平ꎮ

            据(clean data)ꎮ 利用 Trinity 软件(Grabherr et al.ꎬ      1.2.2 SSR 位点挖掘与引物设计  使用 MISA 软件
            2011)对过滤后的高质量数据进行无参基因组组                            (http: / / pgrc.ipk ̄gatersleben. de / misa / misa. html) 对
            装ꎬ获得转录本信息(transcripts)及 Unigenes 序列信               组装后的 Unigenes 进行搜索ꎬ寻找 Unigenes 中的
            息ꎮ 以拼接得到的 Unigenes 作为参考序列ꎬ将样品                      SSR 位点ꎬ设置单核苷酸至六核苷酸重复参数分
   24   25   26   27   28   29   30   31   32   33   34