麻栎(Quercus acutissima Carr.),隶属于壳斗科(Fagaceae)、栎属(Quercus L.),是东亚温带和亚热带地区的优势森林树种,分布于中国、朝鲜、韩国、日本、越南等多个国家。在我国的东经91°-123°,北纬18°-41°地区均有分布,跨越多种气候环境地带,和属间其他物种又存在着杂交现象,因此是研究环境气候变化对物种进化影响以及物种遗传结构和遗传分化的优良材料。本研究在野外实地采样调查的基础上,利用单亲遗传的cpDNA标记和双亲遗传的nSSR标记,探讨了分布区内麻栎的谱系地理结构和遗传多样性,并结合了该物种的当今分布和气候地质事件,模拟了其在末次盛冰期间可能的生境以及冰后期的扩张。主要的研究结果如下:(1)叶绿体片段的联合分析利用4个叶绿体DNA片段(编码区 ycfl和非编码区trnD-trnE、psbA-trnH、rps16-trnQ)对麻栎55个居群502个个体进行了分析。共检测出26种单倍型,单倍型多样性Hd=0.7628,核苷酸多样性π= 0.61×10-3,平均核苷酸差异K=1.079。麻栎的遗传多样性略高(HT=0.770),同时具有较小的基因流(Nm =0.12)。麻栎55个居群不存在谱系地理结构(NST<GST;p<0.05),IBD结果显示其遗传距离和地理距离没有显著的相关性(r = 0.0353;p>0.05)。AMOVA分析显示遗传变异主要发生在居群间(63.25%),居群内的变异较小。失配分析显示麻栎近期存在着扩张,扩张时间约为11.2 kyr。NETWORK和贝叶斯分析将麻栎分为两个枝系,BEAST计算得出麻栎近期共同祖先的分化时间约为6.55百万年前,两个枝系的分化时间约为4.81百万年前和3.42百万年前。单倍型在地理上的分布模式及各居群的遗传多样性暗示秦巴山区、四川盆地、云贵高原和武夷山、天目山地区为麻栎潜在的冰期避难所。(2)核微卫星位点的联合分析利用10对nSSR引物对麻栎33个居群510个个体进行了分析。结果显示麻栎的遗传多样性水平并不高(Na = 5.264,Ne = 3.504,I= 1.218,Ho = 0.283,He = 0.595,PPL=96.06%)。STRUCTURE结果将麻栎分为两组,TESS结果和UPGMA树也支持这一分组;而Mantel检验显示遗传距离和地理距离没有显著相关性(r = 0.138,p>0.05)。Genetix分析显示麻栎33个居群发生了显著的分化(p<0.05),通过AMOVA分析得出遗传变异主要来自于居群内(84.26%),居群间的变异较少;依据STRUCTURE结果分成两组时,遗传变异依然主要来自于居群内(82.59%),组间的变异仅为3.95%。F统计也显示群体间遗传分化较小(FST = 0.203),并且有较大的基因流(Nm= 1.353)。BOTTLENECK分析显示一些居群经历过瓶颈效应。通过标准化的等位基因丰富度和私有等位基因丰富度分析,我们发现秦巴山区、云贵高原、武夷山和天目山地区的这些指数较高,可能为麻栎的冰期避难所;AIS分析也支持这一推断。另外通过EDENETWORK,我们也发现避难所之间存在较多的基因流,暗示麻栎居群间微弱的隔离效应。(3)生态位模拟和冰期避难所叶绿体数据和核微卫星数据都表明麻栎具有多个避难所,生态位模拟显示秦巴山区和云贵高原地区可能为麻栎的冰期避难所。结合这些分析我们认为麻栎的潜在避难所位于秦巴山区和云贵高原地区。冰期后麻栎从避难所向四周适合的生境扩张,在我国华中地区,各避难所的居群得到交流,可能存在着二次接触现象。(4)保护建议虽然麻栎在我国广泛分布,但我们的研究发现其多样性水平不是很高。人类活动影响的加剧,以及种质资源的不合理使用,都对活动对麻栎的分布产生了巨大的影响。麻栎的现今分布已呈现一种破碎化的趋势,且多地为次生林。因此我们建议应禁止滥砍滥伐,恢复其成林规模。