全世界范围内近视的患病率急剧升高，预计2050年近视患者数量将约50亿，达到全球人口的50%，其中高度近视患者可能增加至约10亿^[1]。近视程度的加深伴随着各种高度近视并发症的发生，如后巩膜葡萄肿、脉络膜视网膜萎缩、脉络膜新生血管、近视牵引性黄斑病变等，这些变化通常与视网膜和脉络膜厚度或血管改变相关，最终可导致视力丧失^[2]。

随着近视的加重和眼轴长度（axial length，AL）的增长，眼球形状变为长椭圆体^[3]。最近的一项研究^[4]报道，后极部眼球的形状对近视牵引性黄斑病变的影响甚至比后巩膜葡萄肿更大，视网膜越陡峭，影响视力的近视牵引性黄斑病变越容易发生。这表明异常的视网膜形态不仅可以引起单纯的近视，还可以通过生物力学拉伸引起其他眼部并发症。评估近视患者视网膜形态与眼底结构改变的关系，对理解近视的病理生理学机制有重要意义。以往的研究主要关注黄斑区视网膜厚度与AL的关系，黄斑区视网膜厚度与视网膜形态间的关系尚不明确。因此，本研究应用光学相干断层扫描（optical coherence tomography，OCT）获得较大范围的眼球后极部视网膜图像和视网膜参数，探讨近视眼视网膜形态与黄斑区视网膜厚度的关系。

1 材料与方法 1.1 研究对象和分组

本研究为横断面研究，从2021年中国医科大学附属第一医院眼科招募120例20~30岁的志愿者中纳入符合标准的101例（101只眼）。由于同一个体双眼的对称性，选取右眼进行分析。若由于固定不佳或信号强度 < 8，导致OCT图像质量差，则选择左眼进行分析。最终纳入90只右眼，11只左眼。

纳入标准：（1）经裂隙灯、前置镜和OCT等完整的眼科学检查，为健康眼；（2）最佳矫正视力≥0.8；（3）3次测量的眼压值均 < 21 mmHg。排除标准：（1）存在病理性近视，如严重的后巩膜葡萄肿、脉络膜新生血管等；（2）存在非生理性散光（柱镜超过0.75 D）；（3）任何眼科疾病，如斜视、糖尿病性视网膜病变、青光眼、神经眼科疾病等；（4）有角膜塑形镜佩戴史、屈光手术史、内眼手术史。

根据REZAPOUR等^[5]提出的轴性近视分类，根据AL，将所有研究对象分为3组：无轴性近视组（AL≤24 mm，27例27只眼）、轻度轴性近视组（AL > 24~26 mm，42例42只眼）和高度轴性近视组（AL > 26 mm，32人32只眼）。本研究获得中国医科大学附属第一医院伦理委员会批准，所有研究对象签署知情同意书。

1.2 方法

1.2.1 一般临床资料的收集

收集研究对象的性别、年龄，测量身高和体重，取3次测量值的平均值，计算体重指数（body mass index，BMI），BMI=体重/身高²。

1.2.2 常规眼科检查

使用KR-800电脑验光仪（日本拓普康公司）测量屈光度，取3次测量值的平均值，等效球镜度数（spherical equivalent，SE）=球镜度数+1/2柱镜度数。使用IOLMaster500（德国蔡司公司）测量AL，取5次合格测量值的平均值。使用CY-800电脑非接触式眼压计（日本佳能公司）测量眼压。

1.2.3 视网膜厚度的测量

使用扫频源OCT（VG200S，中国视微影像公司）测量视网膜厚度，黄斑中心凹的位置通过患者固视后机器自动识别来确定，位置不佳者，通过手动调整来确定。基于早期治疗糖尿病视网膜病变研究的环形网格，收集以黄斑中心凹反光点为中心、直径0~1 mm的黄斑中心凹以及直径1~6 mm的中心凹周围环形区域内上方、颞侧、下方、鼻侧共5个区域的视网膜厚度。见图 1。

1.2.4 视网膜形态定量方法

使用扫频源OCT，选取同时经过中心凹和视盘中心的16 mm扫描线用于分析。将经过AL校正后的视网膜图像导入ImageJ程序，获得坐标数据集，导入Anaconda科学计算环境，使用Python编程语言，实现利用最小二乘法拟合椭圆曲线的算法，并在Jupyter Notebook中执行。拟合过程采用一种多维无约束非线性最小化算法，并结合迭代Nelder-Mead算法。每次迭代都构造一个椭圆，计算拟合椭圆与实际曲线间的均方根误差，一直迭代到均方根误差最小值时终止程序，此时可以得到与实际视网膜形态最相符的拟合曲线及其相关参数。该拟合椭圆的离心率（e）越大，后极部视网膜形态越陡峭。

1.3 统计学分析

采用SPSS 26.0软件进行统计分析。采用Shapiro-Wilk检验进行正态性检验，符合正态分布且方差齐性的计量资料用x±s表示，多组间比较采用单因素方差分析，组间两两比较采用LSD-t检验；不符合正态分布的计量资料用M（P₂₅~P₇₅）表示，多组间比较采用非参数Kruskall-Wallis H检验，组间两两比较采用Bonferroni校正。计数资料用率（%）表示，采用χ²检验进行比较。采用Pearson相关分析对AL、SE与视网膜厚度的相关性进行分析；采用Spearman相关分析对e与视网膜厚度的相关性进行分析。将相关性分析中P < 0.1的变量进一步进行多元线性回归分析（逐步法），分析变量与视网膜厚度的关系。P < 0.05为差异有统计学意义。

2 结果 2.1 一般临床资料的比较

本研究共纳入101例研究对象的101只眼。3组比较，性别、年龄、体重、BMI均无统计学差异（均P > 0.05），具有可比性，见表 1。

3组比较，SE、AL和e有统计学差异（均P < 0.001）。无轴性近视组的AL、e明显小于轻度轴性近视组和高度轴性近视组，轻度轴性近视组的AL、e明显小于高度轴性近视组（均P < 0.05）。无轴性近视组的SE明显高于轻度轴性近视组和高度轴性近视组，轻度轴性近视组的SE明显高于高度轴性近视组（均P < 0.05），表明无轴性近视组的近视程度低于轻度轴性近视组和高度轴性近视组，轻度轴性近视组的近视程度低于高度轴性近视组。见表 1。

2.2 视网膜厚度的比较

3组比较，5个区域的视网膜厚度均无统计学差异（P > 0.05）。见表 2。

2.3 视网膜厚度与AL、e、SE的相关性分析

在黄斑中心凹区域内，视网膜厚度与AL、e、SE均无相关性（均P > 0.1）。在中心凹周围视网膜上方、下方、鼻侧3个区域，视网膜厚度均与AL呈负相关（均P < 0.1），在上方、下方、鼻侧、颞侧4个区域，视网膜厚度均与e呈负相关，与SE呈正相关（均P < 0.1）。见表 3。

2.4 视网膜厚度与AL、e、SE的多元线性回归分析

基于相关性分析结果，为矫正混杂关系，将相关性分析中P < 0.1的变量纳入多元线性回归模型中。逐步回归结果显示，e是影响中心凹周围上方（β =-0.345，P < 0.001）、颞侧（β =-0.339，P = 0.001）、下方（β =-0.337，P = 0.001）、鼻侧（β =-0.298，P = 0.003）视网膜厚度的独立相关因素。见表 4。

3 讨论

多项研究表明，近视主要由眼轴延长导致，眼轴过度伸长会导致视网膜和脉络膜发生改变^[6]，这也是导致近视性器质改变的主要原因。本研究对黄斑区视网膜厚度进行分析，结果发现，在根据AL分成的3组间，5个区域的视网膜厚度均无统计学差异。

本研究显示，黄斑中心凹（0~1 mm）的视网膜厚度与AL、e、SE均无显著相关性。目前，对于黄斑中心凹视网膜厚度与AL关系的研究仍存在争议。一项使用扫频源OCT的研究^[7]测量了146名健康志愿者的中心凹视网膜厚度，结果发现，其与AL无关。还有研究发现，中度以上近视眼的中心凹视网膜厚度比正视眼或低度近视眼厚^[8]，认为黄斑中心凹视网膜厚度与AL呈正相关^[9]。本研究结果支持中心凹视网膜厚度与AL无关，且结果进一步表明中心凹视网膜厚度与视网膜形态也无相关性。提示可能存在一种机制来维持中心凹视网膜厚度的稳定性，如黄斑Bruch膜的长度不随近视的增加而增加^[10]，也印证了之前研究^[11]报道的无近视相关视网膜病变的眼睛最佳矫正视力与近视程度无关。在近视眼球轴向发展的过程中，近视相关的视网膜变薄可能始于中心凹周边区域^[12]。

本研究中，对于中心凹周围（1~6 mm）视网膜，上方、下方、鼻侧视网膜厚度与AL、e均呈负相关。这与普遍认为的周边视网膜厚度随AL的增加而减小^[9]结果一致。此外，还有研究^[13]表明，随着AL的增加，黄斑区视网膜神经节细胞复合体的厚度也明显变薄。以往的研究认为，AL增加是导致视网膜结构和功能改变的主要原因。因此，为了得到准确的视网膜厚度与视网膜形态间的关系，本研究在单因素分析的基础上进行多元逐步回归分析，结果显示，e是视网膜厚度的独立相关因素。与AL相比，视网膜形态对周边视网膜厚度的影响更大。这可能是由于视网膜形态能够比眼轴更好地描述眼球形态。还有研究发现，随着近视的加深，视网膜形态在后极附近越陡峭，赤道部越平坦。这种趋势是由于陡峭视网膜的周边受到相对远视的刺激可能会导致补偿性眼球生长，进而使正视眼产生轴向增长的倾向^[14]。这一结果也充分强调了视网膜形态在近视发展中的重要作用。由于本研究为横断面回顾性研究，未来需要进行进一步的纵向研究来完善和证实研究结果。

综上所述，本研究在分析黄斑区视网膜厚度与AL关系的基础上，进一步探索了黄斑区视网膜厚度与视网膜形态间的关系。黄斑中心凹视网膜厚度与视网膜形态无关，中心凹周围视网膜厚度变薄与视网膜形态陡峭更相关。这提示临床上应当重点关注部分近视程度较轻且AL在正常范围但视网膜形态陡峭的人群，进行早期筛查和病情监测，预防病理性近视等并发症的发生。