1. 领域背景与文献
文献英文标题:Structure and evolution of homomorphic UV sex chromosomes in the liverwort Blasia pusilla;
发表期刊:Genome Biology;影响因子:17.9(2024年);研究领域:植物基因组学、性染色体进化
性染色体起源于一对常染色体,初始性别决定区域(SDR)的出现启动重组抑制,随时间推移重组抑制扩展可导致性染色体异型化。在二倍体XY/ZW系统中,非重组的Y/W染色体因重组抑制发生退化,形成异型性染色体;而单倍体UV系统中,雌性(U)和雄性(V)性染色体的重组抑制是对称的,二者退化程度相似,可能长期保持细胞学同型。在光合真核生物类群中,绿藻、红藻、褐藻及苔藓植物均存在UV性染色体,其中地钱模式物种Marchantia polymorpha的异型UV性染色体已被深入解析,其基因组结构、基因组成及性别调控机制研究较为系统。然而,同型UV性染色体仅通过细胞学观察被报道,其基因组序列、结构及基因组成尚未明确,这一空白限制了对苔藓植物性染色体进化关系(同型与异型性染色体的祖先/衍生关系)的理解。本研究通过测序组装地钱Blasia pusilla的雌雄基因组,首次解析同型UV性染色体的基因组结构,为苔藓植物性染色体进化提供关键证据。

2. 文献综述解析
作者对领域内现有研究的分类维度主要包括性染色体的进化模式(二倍体XY/ZW系统vs单倍体UV系统)、研究类群(藻类vs苔藓植物)及性染色体形态(同型vs异型)。
现有研究在二倍体性染色体进化方面已明确重组抑制扩展与染色体退化的关联,而单倍体UV性染色体的研究主要集中在异型类群,如M. polymorpha的微型性染色体,其研究成果包括性染色体的基因组成、重组抑制区域的进化及性别决定调控模块的鉴定,技术方法上采用了长读长测序、Hi-C scaffolding等基因组学技术,优势在于能获得染色体级别的组装并解析精细结构,但局限性在于对同型UV性染色体的研究仅停留在细胞学层面,缺乏基因组层面的解析,无法明确同型与异型UV性染色体的进化起源关系。
本研究的创新价值在于填补了同型UV性染色体基因组研究的空白,首次发现苔藓植物中同型UV性染色体具有明确的假常染色体区域(PAR)和性别决定区域(SDR)的结构,这一结构此前仅在藻类的UV性染色体中被报道,而苔藓植物中尚未发现。通过与M. polymorpha的异型性染色体比较,本研究为Marchantiidae亚纲性染色体的进化提供了两种可能的模型,推动了对苔藓植物性染色体进化轨迹的理解。
3. 研究思路总结与详细解析
本研究的整体目标是解析Blasia pusilla同型UV性染色体的基因组结构,明确其与异型UV性染色体的进化关系;核心科学问题是同型UV性染色体的结构特征、基因组成及进化起源;技术路线遵循“基因组测序组装→性染色体鉴定→结构与进化分析→性别调控模块验证”的闭环逻辑,通过多组学技术结合生物信息学分析,系统解析同型UV性染色体的特征。
3.1 雌雄基因组测序与组装
实验目的是获得Blasia pusilla雌雄株的高质量染色体级基因组组装,为后续性染色体分析提供基础。方法:选取芬兰雄株和中国雌株,雄株采用PacBio长读长测序(约30 Gbp)结合Illumina短读长测序(约307 Gbp),并通过Dovetail Omni-C技术进行scaffolding;雌株采用ONT长读长测序(约115 Gbp)结合Illumina短读长测序,通过Hi-C技术进行scaffolding。组装过程中进行污染去除(利用BlobTools、MEGAN等工具去除非植物序列)、组装polishing(POLCA、Racon、Pilon)及scaffolding优化(3D-DNA、Juicer),最终对雌株组装进行参考引导scaffolding(基于雄株的9条染色体)。结果:雄株组装为9条染色体,基因组大小377 Mbp,符合k-mer估计的基因组大小;雌株组装为14条scaffold,基因组大小390 Mbp,端粒组装在雌株中更完整。BUSCO分析显示基因空间覆盖度良好,基因组重复序列占比约61%,显著高于M. polymorpha的32%,其中Helitrons转座子异常富集(约占注释转座子的16.81%),未检测到全基因组复制信号。产品关联:实验所用关键产品:Omega Bio-Tek的植物EZ DNA提取试剂盒,PacBio SMRT测序平台,Oxford Nanopore Technologies的PromethION测序平台,Illumina NovaSeq测序平台,Dovetail Omni-C试剂盒,Hi-C文库构建试剂盒。
3.2 性染色体鉴定与结构分析
实验目的是鉴定Blasia pusilla的性染色体并解析其精细结构。方法:将Blasia pusilla基因组与M. polymorpha基因组进行共线性分析,比较各染色体的基因密度、重复序列密度,分析雌雄株同源基因对的同义替换率(Ks)以判断重组抑制程度。结果:共线性分析显示Blasia pusilla与M. polymorpha染色体存在多次断裂、融合和易位,未检测到与M. polymorpha UV性染色体的共线性块,但9号染色体的末端区域与其他染色体特征不同。进一步分析发现9号染色体可分为两个区域:约2/3的区域(0-21.41 Mbp)为假常染色体区域(PAR),其基因密度、重复序列密度与常染色体类似,雌雄同源基因对的同义替换率低,显示高重组率;约1/3的末端区域(21.41 Mbp至染色体末端)为性别决定区域(SDR),基因密度显著降低至4.90基因/Mbp(n=雄株9号染色体,P<0.05,与常染色体比较),Copia和Gypsy类LTR反转录转座子富集,雌雄同源基因对的同义替换率高,显示重组抑制。雌株组装中性别决定区域与假常染色体区域位于不同scaffold,但二者的基因密度和重复序列组成与雄株相似。产品关联:文献未提及具体生物信息学软件的品牌,领域常规使用Minimap2、MCScanX、Circos、ggplot2等工具进行基因组比较与可视化。

3.3 SDR基因组成与进化分析
实验目的是解析性别决定区域的基因组成、起源及进化特征,验证重组抑制的存在。方法:搜索Blasia pusilla基因组中M. polymorpha UV gametologs的同源基因,构建系统发育树分析基因分化时间,比较性别决定区域、假常染色体区域与常染色体的重复序列景观及基因特征。结果:Blasia pusilla的性别决定区域中包含11个M. polymorpha gametologs的同源基因,其中7个基因的分化早于两个物种的分歧,显示重组抑制的古老起源;性别决定区域的雌雄同源基因对同义替换率分为两个类群,混合分布在区域中,无明显的分层趋势;重复序列景观显示性别决定区域中古老转座子的比例高于假常染色体区域,假常染色体区域则富集年轻基因(孤儿基因比例高于常染色体)。此外,性别决定区域中还包含部分物种特异性基因,可能为近期捕获的基因。产品关联:领域常规使用OrthoFinder、MAFFT、IQ-TREE、RepeatModeler、RepeatMasker等工具进行基因家族分析、系统发育分析及重复序列注释。

3.4 性别决定调控模块验证
实验目的是验证Blasia pusilla中是否存在保守的性别决定调控模块。方法:搜索M. polymorpha中核心性别决定基因(BPCU/V、SUF、FGMYB)的同源基因,分析其基因组定位,并通过RNA-seq检测表达模式。结果:Blasia pusilla的BPCU和BPCV同源基因分别位于U和V性染色体的性别决定区域中,FGMYB同源基因位于1号常染色体,其侧翼基因组序列与M. polymorpha具有共线性;RNA-seq显示雌株中FGMYB表达量低,雄株中存在SUF的反义转录本,与M. polymorpha的性别调控模块表达模式一致,表明该调控模块在Blasia pusilla中保守存在。产品关联:实验所用关键产品:Sigma-Aldrich的Spectrum Plant Total RNA Kit,Illumina Truseq stranded mRNA Library Prep kit,Illumina Novaseq 6000测序平台。

4. Biomarker研究及发现成果
本研究中的Biomarker主要包括性染色体性别决定区域中的核心性别决定基因(BPCU、BPCV)以及同型UV性染色体特有的假常染色体区域-性别决定区域结构,其筛选与验证逻辑为“基因组比较鉴定→系统发育分析验证→表达模式确认”的完整链条。
Biomarker的来源为Blasia pusilla的雌雄基因组及转录组数据,验证方法包括基因组共线性分析、系统发育树构建、RNA-seq表达分析。其中,BPCU和BPCV作为性别决定的关键基因,其特异性定位于性别决定区域,雌雄株中的同源基因显示高同义替换率(文献未明确提供具体数值,基于图表趋势推测),表明重组抑制;假常染色体区域-性别决定区域结构作为同型性染色体的特征,其假常染色体区域的基因密度与常染色体相近(文献未明确提供具体数值,基于图表趋势推测),性别决定区域基因密度为4.90基因/Mbp(n=雄株9号染色体,P<0.05,与常染色体比较),重复序列占比显著高于假常染色体区域和常染色体。
核心成果方面,本研究首次发现苔藓植物同型UV性染色体的假常染色体区域-性别决定区域结构,这一结构为Marchantiidae亚纲性染色体进化提供了两种可能模型:一是同型性染色体的假常染色体区域-性别决定区域结构为祖先状态,异型性染色体通过假常染色体区域的丢失或易位进化而来;二是同型性染色体起源于常染色体与微型性染色体的融合。此外,保守的BPCU/V-SUF-FGMYB性别调控模块在Blasia pusilla中的存在,表明该模块在苔藓植物中具有古老的起源。这些成果填补了同型UV性染色体基因组研究的空白,推动了对苔藓植物性染色体进化轨迹的理解。