关于人类大脑皮层转录组层级化组织的神经影像学研究报告
作者与发表信息 本研究由Joshua B. Burt（耶鲁大学物理系）、Murat Demirtaş（耶鲁大学医学院精神病学系）、William J. Eckner（耶鲁大学物理系）、Natasha M. Navejar（杜兰大学神经科学项目）、Jie Lisa Ji（耶鲁大学跨部门神经科学项目）、William J. Martin（Blackthorn Therapeutics公司）、Alberto Bernacchia（剑桥大学工程系）、Alan Anticevic（耶鲁大学医学院精神病学系）以及通讯作者John D. Murray（耶鲁大学物理系与精神病学系）共同完成。该研究于2018年发表在《自然-神经科学》（*Nature Neuroscience*）期刊上。
学术背景与研究目标 本研究隶属于系统神经科学与计算神经科学的交叉领域，主要探讨大脑皮层的宏观组织原则。长期以来，对非人灵长类动物的侵入性解剖追踪研究表明，大脑皮层在结构和功能上呈现出层级（hierarchy）组织，即从处理基本感觉信息的初级感觉区（如初级视觉皮层V1）到执行复杂认知功能的高级联合皮层（如前额叶皮层），存在一个有序的等级排列。这种层级结构反映在脑区间连接的模式上：较低层级的区域主要向较高层级发出前馈（feedforward）投射，而较高层级则向较低层级发出反馈（feedback）投射，且这些投射在皮层各层的起源具有特异性。然而，由于无法在人类身上进行侵入性的追踪研究，人类皮层是否以及如何遵循类似的层级组织原则，一直是一个悬而未决的问题。
与此同时，大规模转录组学（transcriptomics）研究，如艾伦人脑图谱（Allen Human Brain Atlas, AHBA），揭示了基因表达在人类大脑皮层不同区域存在显著的异质性。这种转录组架构的差异与皮层功能网络的活动时空结构以及解剖连接模式相关。然而，皮层微观结构（如细胞组成、突触特性）和基因表达的大尺度空间变化，是否也存在一个统一、简洁的组织原则，例如是否也沿着感觉-联合皮层轴呈现层级化的梯度变化，尚不清楚。
因此，本研究旨在解决两个核心科学问题：第一，能否找到一种非侵入性的神经影像学指标，作为人类皮层解剖层级的可靠代理（proxy）？第二，人类皮层的转录组架构（即基因表达的空间模式）是否也表现出层级化的组织？如果是，这种层级化的基因表达梯度是否与皮层微环路（microcircuit）的特化（如特定细胞类型、突触受体的分布）相关，并可能对理解脑疾病的病理生理学有启示？本研究假设，皮层微环路的特化涉及跨皮层区域的基因表达层级梯度，并且这种层级轴是连接皮层解剖、功能和转录组架构的共同组织原则。
详细研究流程 本研究整合了来自人类和猕猴的多种数据类型，包括结构磁共振成像（MRI）数据、解剖学追踪数据以及基因表达数据，流程严谨且多步骤验证。
流程一：验证T1w/T2w图谱作为层级代理的可行性 首先，研究团队需要一种非侵入性指标来量化人类皮层的解剖层级。他们聚焦于MRI衍生的T1加权与T2加权图像比值图（T1w/T2w map）。已有研究表明，该图谱对皮层灰质的髓鞘含量敏感，并能勾勒出已知的细胞构筑和髓鞘构筑边界。研究者观察到，在人类和猕猴的组平均T1w/T2w图谱上，初级感觉皮层（如视觉、体感、听觉皮层）的值较高，而联合皮层（如默认模式网络相关区域）的值较低，且这种模式在个体间稳定、在物种间同源。
为了验证T1w/T2w值能否反映解剖层级，研究者在拥有详细解剖追踪数据的猕猴皮层上进行了关键验证。他们使用了来自公开数据库的猕猴皮层间层状特异性投射数据。对于每一对有连接的脑区，可以计算一个指标：标记神经元中位于浅层（supragranular layers）的比例（SLN）。高SLN值通常对应前馈连接（从低层级到高层级），低SLN值对应反馈连接。研究者发现，相连脑区之间的T1w/T2w值差异与SLN值显著相关，表明T1w/T2w差异能预测连接的层状特异性。
更重要的是，他们采用广义线性模型（Generalized Linear Model, GLM），利用所有成对的SLN数据，拟合出一个全局最优的、自洽的解剖层级排序（即每个脑区被赋予一个层级水平值）。分析结果显示，猕猴皮层T1w/T2w图谱的值与这个模型估计的解剖层级水平之间存在极强的负相关（斯皮尔曼秩相关 rs = -0.76， p < 10^-5）。这意味着T1w/T2w值高的区域（如初级感觉区）位于层级底部，而T1w/T2w值低的区域（如联合皮层）位于层级顶部。
研究还将T1w/T2w图谱与另外两个候选代理指标——皮层厚度图谱和从初级视觉皮层（V1）出发的测地线距离图谱——进行了比较。结果发现，T1w/T2w图谱与解剖层级的关联性显著强于后两者。这些证据共同支持了“T1w/T2w图谱可以作为灵长类（包括人类）皮层解剖层级的非侵入性代理指标”这一结论。
流程二：分析人类皮层基因表达的层级梯度 在确立了T1w/T2w图谱作为层级代理的有效性后，研究者转向分析人类皮层的转录组数据。他们使用了艾伦人脑图谱（AHBA）的微阵列基因表达数据，该数据来自6名正常人类捐赠者左半球大脑的数百个样本。
数据处理流程细致：首先，将每个捐赠者的基因表达样本从其立体定位坐标映射到其个体化的大脑皮层表面网格上。接着，利用人类连接组计划（Human Connectome Project, HCP）的多模态分区（MMP1.0）将左半球皮层划分为180个区域。对于每个基因，在每个捐赠者内部，选择最具代表性的探针，并在每个皮层分区内平均表达值。最后，跨捐赠者平均，得到每个基因在180个分区的组平均表达图谱。经过筛选，最终获得了16,088个基因的可靠表达谱。
核心分析是计算每个基因的表达谱与人类组平均T1w/T2w图谱之间的斯皮尔曼秩相关，即T1w/T2w图谱相关性（T1w/T2w map correlation, TMC）。由于T1w/T2w值与解剖层级负相关，因此，负TMC意味着该基因在高层级（联合皮层）表达更高，呈现正向层级梯度；正TMC意味着该基因在低层级（感觉皮层）表达更高，呈现负向层级梯度。
流程三：将人类转录组梯度与猕猴微观解剖测量关联 为了验证从人类基因表达中推断出的层级梯度的生物学合理性，研究者将其与来自猕猴的、更直接的微观解剖测量结果进行了比较。这构成了跨物种、跨模态的强有力验证。
流程四：探究转录组变异的主导模式 研究者进一步使用主成分分析（Principal Component Analysis, PCA）来识别人类皮层基因表达空间变异的主导模式。他们对所有基因以及几个特定的功能基因集（如脑特异性基因、神经元特异性基因、少突胶质细胞特异性基因、突触组基因）分别进行了分析。
结果发现，对于所有测试的基因集，第一主成分（PC1），即能解释最大比例基因表达方差的空间模式，均与T1w/T2w图谱的地形高度相似（相关系数在0.80-0.81之间）。PC1同样清晰地区分了初级感觉运动皮层和联合皮层。这表明，T1w/T2w图谱所代表的层级轴，捕获了人类皮层转录组变异中最主要、最一致的维度。与皮层厚度、距V1距离等其他候选图谱相比，T1w/T2w图谱与PC1的相关性最强，能解释更多的基因表达方差。
流程五：分析基因表达稳定性与功能、疾病关联 1. 差异稳定性（Differential Stability, DS）：研究者计算了每个基因的皮层差异稳定性（DSc），即个体间基因表达空间图谱相关性的平均值。他们发现，DSc高的基因（即个体间表达模式稳定的基因）其TMC的绝对值也倾向于更大。当过滤掉低DSc基因后，TMC的分布从集中在零附近变为明显的双峰分布（强正TMC和强负TMC），且PC1能解释的方差比例大幅上升。这表明，在个体间稳定表达的基因，更倾向于表现出强烈的层级梯度。 2. 功能与疾病富集分析：研究者检验了具有强TMC（无论是正还是负）的基因是否在某些功能类别或疾病相关基因集中过度出现。结果发现，具有强TMC的基因在多个功能基因本体（如生物过程、细胞组分、分子功能、microRNA结合位点、药物靶点）中均显著富集。更重要的是，对神经精神疾病相关基因集的分析显示，具有强负TMC的基因（即在联合皮层高表达的基因）在精神分裂症、双相情感障碍、自闭症谱系障碍和抑郁症等疾病的基因集中显著富集。例如，阿尔茨海默病的关键风险基因*APOE*和帕金森病的关键基因*SNCA*都表现出强负TMC。
主要研究结果 1. T1w/T2w图谱是非侵入性解剖层级代理：在猕猴皮层上严格验证表明，T1w/T2w值与基于层状投射的模型估计解剖层级高度负相关，且其预测能力优于皮层厚度和距V1距离等其他指标。这为在人类研究中使用该指标奠定了基础。 2. 人类皮层基因表达存在系统性层级梯度：大量与微环路功能相关的基因，其表达沿感觉-联合皮层轴呈现显著的梯度变化。这些基因涉及细胞构筑分层、特定抑制性神经元类型、树突棘密度（通过受体基因间接推断）以及兴奋/抑制性突触受体亚基。 3. 转录组变异的主轴与层级轴对齐：主成分分析揭示，人类皮层基因表达空间变异的最主要模式（PC1）与T1w/T2w图谱高度一致，说明层级是组织皮层转录组架构的一个核心原则。 4. 稳定表达的基因偏好呈现层级模式：在个体间表达模式稳定的基因（高DSc基因），更可能表现出强烈的层级梯度表达，且这些基因构成了转录组低维结构的主体。 5. 层级化基因与脑疾病相关：在联合皮层高表达（强负TMC）的基因，在多种神经精神疾病的遗传风险基因集中显著富集，提示这些疾病可能对皮层层级的不同部分有差异化的影响。
结论与意义 本研究得出结论：人类大脑皮层的转录组和微观解剖结构共享一个共同的主要变异轴，即皮层层级。T1w/T2w图谱作为一种非侵入性神经影像学指标，能够可靠地索引这一层级轴。皮层微环路的特化——包括细胞类型组成、突触特性等——通过基因表达的层级梯度得以实现，这些梯度可能构成了皮层宏观功能特化的微观基础。研究还发现了人类与猕猴在皮层组织原则上的保守性。
其科学价值在于： 1. 提出了统一框架：将皮层解剖、功能和转录组的大尺度组织统一到一个简洁的层级原则之下。 2. 建立了关键桥梁：提供了连接非侵入性神经影像（T1w/T2w）、微观解剖和分子生物学（基因表达）的实证桥梁，使得在活体人类中研究皮层层级及其微观基础成为可能。 3. 揭示了疾病新视角：发现与高阶认知相关的联合皮层高表达基因与多种精神疾病风险基因重叠，为理解这些疾病的皮层特异性病理机制提供了新线索。 4. 启发了新研究方向：为未来研究皮层功能特化的发育轨迹、神经调控药物的靶向性（因受体表达存在层级梯度）以及更精细的空间转录组图谱分析指明了方向。
研究亮点 1. 跨物种、多模态整合：创新性地将人类神经影像、人类转录组数据和猕猴微观解剖数据相结合，进行了严谨的交叉验证，极大地增强了结论的说服力。 2. 方法学创新：明确验证并提出T1w/T2w图谱作为人类皮层解剖层级的非侵入性代理指标，解决了该领域长期缺乏合适工具的问题。 3. 发现转录组主导模式：首次通过PCA系统证明，人类皮层基因表达空间变异的最主要模式与结构层级轴高度一致，揭示了转录组组织的核心原则。 4. 连接微观与宏观：成功地将微观层面的基因表达、细胞类型特征与宏观层面的皮层层级和功能网络联系起来，为系统神经科学提供了分子层面的解释。 5. 临床相关性：将基础科学的发现与神经精神疾病的遗传风险联系起来，提升了研究的转化医学价值。
其他有价值的内容 研究还讨论了T1w/T2w图谱的局限性，指出其不仅对髓鞘敏感，也可能反映其他微结构特性（如细胞大小、密度、树突复杂度、铁含量等）。同时，研究承认存在多种功能层级（如不同感觉模态内），未来的研究可以结合多种神经影像指标（如弥散加权成像）来揭示皮层组织更复杂的多维原则。此外，研究指出AHBA数据空间采样相对稀疏，可能导致对结构-转录关联的估计偏保守。这些讨论体现了研究的严谨性和开放性。