《Human Molecular Genetics》：Genomic profiling of microphthalmia trios reveals novel pathogenic drivers and repetitive element instability
编辑推荐：
先天性小眼球（CMIC）是一种具有显著遗传异质性的罕见眼部畸形。为探究其基因组景观，研究人员对10个患病先证者-双亲三联家系进行了全基因组测序（WGS），并整合分析了遗传变异、候选胚系新生变异（gDNM）及串联重复序列变异。经质量、群体频率及新生突变概率过
先天性小眼球（CMIC）是一种具有显著遗传异质性的罕见眼部畸形。为探究其基因组景观，研究人员对10个患病先证者-双亲三联家系进行了全基因组测序（WGS），并整合分析了遗传变异、候选胚系新生变异（gDNM）及串联重复序列变异。经质量、群体频率及新生突变概率过滤后，共优先筛选出涉及37个基因的45个候选胚系新生位点，包括POMT1、MYH11和NR4A1中的变异。ANXA9和SCARA5中的罕见纯合变异在各自家系中呈现隐性共分离。整合基因表达谱与已发表的功能证据，凸显了其与眼部发育、钙稳态、铁转运及细胞应激反应之间潜在的生物学联系。重复序列分析进一步在先证者5、9和10中识别出显著的重复长度离群信号，其中包括先证者10中一个与SOX9/SOX9-AS1位点重叠的候选AACCCT重复长度改变。综上，这些发现拓宽了CMIC候选基因组变异的谱系，并提名了可能参与其分子发病机制的基因与调控位点。其具体致病作用与功能效应尚需在独立队列和实验模型中进一步验证。
# 四、论文解读
先天性小眼球（congenital microphthalmos, CMIC）是一种严重的眼部发育畸形，以眼球体积显著缩小为主要特征，发病率约为1.5/10000活产儿，与无眼畸形、眼组织缺损统称为MAC，是儿童致盲的重要原因之一，约占11%。该病可孤立发生，也可作为综合征表现之一，约20%–45%的病例伴随系统性发育异常。其病因复杂，遗传易感性与致畸环境暴露均为主要 contributor。目前已知与MAC直接或间接相关的基因超过100个，包括ALDH1A3、BCOR、GJA8、PAX6、OTX2和SOX2等，但这些已知基因的致病变异仅能解释部分病例，临床遗传学诊断率仅为20%–30%，重症或双侧病例的诊断率相对较高。既往研究提示，家系研究有助于深入理解遗传风险并更精准定位致病因素，例如通过对11个正视眼（nanophthalmos）三联家系的研究在PRSS56和MFRP两个基因中鉴定出3个新生突变（DNM）。然而，CMIC的基因组变异图谱，尤其是非编码区及串联重复序列层面的变异，仍未被充分刻画。基于此，研究人员对2011至2017年间招募的10个CMIC先证者-双亲三联家系开展全基因组测序（WGS），整合遗传变异、候选胚系新生变异与全基因组串联重复变异分析，并结合基因表达谱、功能注释及已发表生物学证据，对可能与眼部发育相关的基因和调控位点进行优先级排序，以拓宽CMIC候选变异谱并为分子发病机制研究提供新的遗传与调控候选靶点。
在技术方法层面，研究人员采集10个先证者及其双亲共30人的外周血样本（来自北京同仁医院，2011—2017年），提取基因组DNA后构建文库，在Illumina HiSeq 2500平台进行双端测序。原始reads经trimming后以BWA v0.7.12比对至hg19参考基因组，Picard标记PCR重复，GATK v3.4完成变异检测与初步过滤，ANNOVAR结合RefSeq基因模型与dbNSFP v4.2a进行注释。候选变异要求在各群体数据库（1KGP、PGG Han、gnomAD）中缺失或AF<0.01，并采用五种计算机预测工具（MVP、MetaLR、MetaSVM、REVEL、VEST4）进行排序。遗传变异与孟德尔不一致基因型通过自编Perl脚本及VCFtools鉴定，候选新生变异经DeNovoGear（pp_dnm>0.8）筛选；家系关系用Somalier评估。重复序列分析采用ExpansionHunter Denovo离群模式（以70例北方汉族1KGP样本为对照），并以ExpansionHunter进行三联家系靶向基因分型，MSIsensor2计算探索性微卫星不稳定（MSI）评分，TelSeq估计端粒长度。眼部表达谱来自人类蛋白质图谱（HPA）单细胞数据、眼组织数据库（OTDB）及GEO数据集GSE119320（人视网膜类器官发育表达谱）。
结果部分：（1）样本采集与全基因组测序——完成10个三联家系共30人的WGS，测序与比对质量指标汇总于补充材料。（2）全基因组变异检测、过滤与注释——样本层面注释计数约446万（先证者9）至933万（先证者2），约90%注释位于基因间区或内含子区，外显子与ncRNA外显子合计约占0.55%；约90%变异可注释于dbSNP，未注释变异平均约80万/样本（先证者9为36万）。（3）遗传变异与孟德尔违背识别——每个三联家系纳入遗传变异分类的位点为805095至1019595个（均值923601），其中先证者纯合变异平均49491个/家系；初始三联比较识别出每家系14087至335190个孟德尔基因型不一致，家系5、9、10比例最高（13.33%、10.79%、28.65%），其余家系为1.36%–3.89%；经群体频率过滤与功能效应评估后保留6248个变异用于下游优先级排序，最终聚焦ANXA9、SCARA5、POMT1、MYH11和NR4A1。（4）CMIC三联家系中的罕见纯合遗传候选变异——家系9中ANXA9的I64T变异（rs143089777）为先证者纯合、双亲杂合，一般人群频率0.000416，1KGP与PGG汉族人群中未见纯合个体，该替换位于第一内部重复结构域且受影响异亮氨酸残基保守；家系10中SCARA5的R43W变异（rs138480005）为先证者纯合，人群频率0.000596，SCARA5介导铁蛋白摄取与视网膜铁转运，为其优先级排序提供生物学背景。（5）孟德尔基因型不一致与候选胚系新生变异——过滤后保留1484个不一致位点用于新生突变评估，除家系5、9、10外每家系为11–74个，不一致数量与父母年龄无关，仅1.67%位点（对应5.58%受累基因）为两个及以上家系共享；经DeNovoGear及人工审校优先筛选出37个基因中的45个候选gDNM位点（48个先证者-变异观察），家系5中包括POMT1 c.396A>T（p.Leu132Phe）与MYH11 c.1301G>A（p.Arg434His），均被ClinVar注释为意义未明变异；POMT1介导α-抗肌萎缩相关蛋白聚糖糖基化，其致病变异可导致伴眼部畸形的肌营养不良糖基化病（如Walker-Warburg综合征），MYH11表达于血管壁细胞与角膜内皮细胞，STRING分析提示其与ACTB、ACTG1存在功能关联（综合评分0.94与0.91）；NR4A1中两个错义变异p.Gln5Arg与p.Gly258Arg见于先证者5，NR4A参与DNA双链断裂修复，NR4A1与R环形成、染色质可及性及复制应激转录反应相关。此外，研究人员汇总了43个基因的文献眼部证据，其中10个具有眼部发育或表达证据的基因被重点总结，CHD7为已知CMIC相关基因，其余基因证据主要为实验性或表达层面。（6）CMIC三联家系中的重复长度变异——ExpansionHunter Denovo在先证者5、9、10中识别出较其余先证者及对照数据集更丰富的重复离群信号，端粒样基序占较大比例（三家系平均比例0.794），全基因组重复离群信号与细胞遗传学脆性位点存在空间重叠，与GC丰富区及串联重复区的相似系数分别为0.614与0.612；经ExpansionHunter靶向分型与过滤，三家系共保留26个候选重复长度改变，其中先证者10中chr17:70120823–70120829（hg19）的候选AACCCT重复长度改变估计基因型为1/38（双亲均为1/1），定位于SOX9/SOX9-AS1区域并重叠外显子、3′UTR及增强子注释；SOX9表达于视网膜祖细胞并为正常Müller胶质细胞发育所必需，还调控视网膜色素上皮的视觉循环基因表达，Odd Sex小鼠模型中视网膜色素上皮的不适当Sox9表达可导致伴白内障的小眼球，为重复长度变异影响眼部发育提供了具体假设。
讨论部分总结：本研究通过整合10个三联家系的遗传变异、候选胚系新生变异与重复长度谱分析，将候选发现拓展至单个编码变异之外，主要贡献在于识别出可在其家系与发育背景中进一步研究的特定基因与基因组位点。ANXA9与SCARA5的罕见纯合基因型（先证者纯合、未受累双亲杂合）支持隐性模型下的优先级排序，SCARA5在铁蛋白介导的铁转运中的作用已在实验中得到证实。候选新生变异中，POMT1将分析与α-抗肌萎缩相关蛋白聚糖糖基化在伴眼部畸形疾病中的已知作用相联系；MYH11的表达与网络分析为细胞骨架及血管过程研究提供依据；NR4A1在DNA修复与转录应激调控中的功能提示先证者5候选变异的独立研究方向。重复长度发现提示常规编码变异分析之外的变异亦可能与CMIC相关，SOX9/SOX9-AS1重叠候选因提供具体发育假设而尤具信息价值，DNA修复通路参与重复长度维持，为结合细胞应激相关候选研究重复变异提供背景，但先证者5的NR4A1发现并不能确立重复离群信号家系存在共同的DNA修复缺陷。主要局限性包括队列小、单中心、临床表型与亚型分层不完整、缺乏独立重复验证及正交与功能验证，队列对正式基因负荷检验效力有限；横断面血液WGS无法判定重复长度信号是否参与疾病、是否为继发或偶然/技术性变异，也无法确定其在胚胎眼组织中的效应。解读文章至此，以下翻译研究结论部分：综上，本研究拓宽了CMIC的候选基因组景观，并为其后续遗传学与功能研究提供了聚焦的基础；优先筛选的变异及其所提出的机制有待独立队列、正交基因分型及发育模型加以检验与区分。
