1.?larINDEL_FZ与棉花种子短绒和叶茸毛性状显著相关

本研究中，作者利用215份亚洲棉自然群体，在前人基于SNP标记的GWAS基础上（Du et al.?2018），基于SV（结构变异）标记对种子短绒和叶茸毛性状进行GWAS，鉴定得到位于A08染色体末端的一个非编码插入片段larINDEL_FZ与短绒和叶茸毛表型显著相关（图1a, b）；连锁不平衡分析发现，在关联区域有4个连锁不平衡(LD)位点，SNP_FZ不位于这些位点内（图1c, d）；对SNP_FZ和larINDEL_FZ基因型与种子短绒的相关性进行研究，发现只有larINDEL_FZ基因型与无茸毛表型高度相关（图1e）。以上结果表明，插入片段larINDEL_FZ与种子短绒和叶茸毛性状显著相关。

图1 larINDEL_FZ与棉花种子短绒和叶茸毛性状显著相关

2.?larINDEL_FZ的遗传结构及其与种子短绒和叶茸毛的相关性

为了验证插入片段larINDEL_FZ，作者利用19个重叠引物在具有代表性的短绒和无绒材料（GA0146和GA0149）中进行序列分析，发现短绒材料中有一个大小为6226 bp的插入片段larINDEL_FZ（图2a）；进一步分析发现，该插入片段由一个约5?kb的重复序列和一个约1.2?kb的易位片段组成（图2b）；作者在GWAS群体以及GA0146和GA0149构建的包含1212单株的F₂群体中分析larINDEL_FZ和表型的相关性，结果表明，位于亚洲棉第8染色体上的6226?bp的larINDEL_FZ与种子短绒和叶茸毛显著相关（图2c）。

?图2 larINDEL_FZ的遗传机制

3.larINDEL_FZ是激活GaFZ表达的关键调控增强子

作者对关联区域内8个基因的表达模式分析，发现在短绒发育起始阶段，无绒材料中GaFZ基因表达显著高于短绒材料（图3a）；从GWAS群体中随机挑选各6个短绒和无绒材料进行基因和启动子序列分析，发现larINDEL_FZ存在于所有的无绒材料中（图3b, c）；荧光素酶活性分析实验结果表明，larINDEL_FZ可能扮演远端增强子角色特异激活其下游短绒发育关键基因的GaFZ?表达（图3d, e）

?图3 larINDEL_FZ激活GaFZ基因表达

4.GaFZ基因转基因植株的表现

为了进一步验证GaFZ基因的功能，作者分别在拟南芥和有短绒的陆地棉中过表达GaFZ基因，发现转基因拟南芥叶茸毛的发育收到了明显的抑制，同时转基因棉花株系的茎毛、叶茸毛和短绒的发育也均受到了不同程度的抑制（图a-h）；但是，转基因棉花和拟南芥中根毛的数量均没有发生变化；qRT-PCR分析显示，过表达株系中GaFZ基因显著高于野生型（图4i-j）；结果表明，GaFZ基因可以抑制叶茸毛和短绒纤维的发育，而对根毛的发育没有抑制作用。

?图4 GaFZ基因转基因棉花的表现

5.GaFZ基因调控亚洲棉短绒发育的遗传机制

作者进一步在短绒和无绒棉花材料和过表达的拟南芥材料中分析全基因共表达谱，在拟南芥和棉花的短绒关键发育阶段，共鉴定出1343个和367个差异表达基因，这些基因主要参与次生代谢，如脂肪酸伸长和角质、栓质和蜡质的生物合成；结合互作分析发现，GaFZ基因可能通过抑制超长链脂肪酸伸长途径中关键基因的表达，进一步抑制脂肪酸延伸的下游途径，从而影响短绒/叶茸毛的发育。作者基于基因功能分析,提出了一个亚洲棉短绒发育的可能机制（图5）。

?图5 larINDEL_FZ调控亚洲棉短绒发育的可能机制