《Viruses》：Genetic Variability of IFNAR2 Gene and Risk of COVID-19 Critical Illness Among the Tunisian Population
编辑推荐：
专业译文
IFNAR2基因编码I型干扰素受体的一个必需组分，对机体抗病毒防御至关重要。已有证据表明，IFNAR2的遗传变异可能影响COVID-19的严重程度，但在非洲人群中相关结果尚无定论。因此，本研究旨在探讨三个IFNAR2单核苷酸多态性（SNP）与突尼斯人
专业译文
IFNAR2基因编码I型干扰素受体的一个必需组分，对机体抗病毒防御至关重要。已有证据表明，IFNAR2的遗传变异可能影响COVID-19的严重程度，但在非洲人群中相关结果尚无定论。因此，本研究旨在探讨三个IFNAR2单核苷酸多态性（SNP）与突尼斯人群COVID-19危重型的关联。研究共纳入413例COVID-19患者，分为住院危重型（重症监护室（ICU）组，n=46）和非住院组（n=367）。采用TaqMan实时荧光定量PCR（RT-PCR）技术对三个IFNAR2多态性位点（rs1051393、rs2236757、rs9984273）进行基因分型。结果显示，rs2236757（比值比（OR）=0.31，95%置信区间（CI）：0.16–0.60，p=4×10?4）和rs9984273（OR=0.41，95%CI：0.22–0.76，p=0.0055）与COVID-19严重程度存在关联，且具有潜在的保护效应。然而，rs1051393未显示与疾病严重程度显著相关。此外，探索性单倍型分析提示，TAC单倍型可能具有保护性关联（OR=0.11，95%CI：0.02–0.86，p=0.036），但该结果基于极少的携带者数量，且经多重检验校正后不再显著。GAC单倍型在ICU患者中完全缺失，经连续性校正后得到OR值为0.03（95%CI：[0.002–0.53]）；该估计值依赖于零细胞计数，仅为假设生成性质。这些初步发现来自一项单中心队列，其中危重症患者数量有限，尚需在更大规模、独立、多中心的人群中进行重复验证。
论文解读文章
一、研究背景与目的
COVID-19大流行由严重急性呼吸综合征冠状病毒2型（SARS-CoV-2）引起，其临床表现从无症状到危重症差异巨大。尽管高龄、男性、合并症（如糖尿病、肥胖、高血压）等已知风险因素可部分解释病情差异，但许多无上述危险因素的年轻患者仍发展为重症，提示内在遗传因素可能在免疫应答中发挥关键作用。I型干扰素（IFN）通路是抗病毒先天免疫的核心，而IFNAR2基因编码I型干扰素受体亚基，其功能异常可能直接影响干扰素信号效率。此前，多项全基因组关联研究（GWAS）及候选基因研究已在欧洲、巴西、墨西哥、巴勒斯坦等不同人群中报道了IFNAR2单核苷酸多态性（SNP）与COVID-19严重程度的关联，但结果方向不一致，且北非人群特别是突尼斯人群的数据极度缺乏。突尼斯拥有独特的遗传多样性，可能影响免疫应答模式。因此，本研究旨在评估三个IFNAR2 SNP（rs1051393、rs2236757、rs9984273）与突尼斯人群COVID-19危重型的关联，填补该地区的研究空白。
二、研究方法概述
这是一项回顾性遗传关联研究。研究对象来源于突尼斯Ben Hassan Salem医学分析实验室，共纳入413例于2021年3月至8月期间确诊的、未接种疫苗的成年COVID-19患者，分为非住院组（367例）和住院危重型（需入住重症监护室（ICU）的缺氧性肺炎患者，46例）。研究人员从EDTA抗凝全血中提取DNA，采用盐析法纯化，并使用NanoDrop 2000测定浓度和纯度。基因分型采用TaqMan实时荧光定量PCR（RT-PCR）探针法完成。统计分析借助SPSS 23和SNPStats在线平台，包括哈迪-温伯格平衡检验、Fisher精确检验、Logistic回归（校正年龄和性别）以及单倍型分析（基于期望最大化算法）。多重比较采用Benjamini-Hochberg错误发现率（FDR）校正。
三、研究结果
3.1 临床与人口学特征
研究共纳入413例参与者。非住院组平均年龄（51.5±16.65岁）略高于住院组（48.5±17.2岁），但差异无统计学意义（p=0.251）。两组性别分布无显著差异（p=0.871），女性均占多数（非住院组65.9%，住院组67.4%）。生活方式因素（烟草、水烟、酒精）及各合并症（高血压、糖尿病、心血管疾病等）均未显示与危重症显著关联。这一与传统风险因素不符的结果被归因于住院组样本量过小。
3.2 IFNAR2变异与COVID-19严重程度的关联
三个SNP的基因型分布在两组中均符合哈迪-温伯格平衡。Fisher精确检验显示，rs2236757和rs9984273的基因型频率在两组间存在显著差异。经年龄和性别校正的Logistic回归分析中，依据最小赤池信息准则（AIC）和贝叶斯信息准则（BIC），rs2236757的最佳拟合模型为显性模型：以GG基因型为参照，AA+AG基因型在住院组中频率显著降低（37% vs. 62.4%），与COVID-19危重症风险降低相关（OR=0.31，95%CI：0.16–0.60，p=4×10?4），且经FDR校正后仍显著（FDR=0.0036）。rs9984273同样适用显性模型：CC+TC基因型在住院组中频率较低（50% vs. 71.9%），与保护效应相关（OR=0.41，95%CI：0.22–0.76，p=0.0055，FDR=0.0165）。rs1051393未发现显著关联。交互作用分析未发现性别对上述关联有显著修饰效应。
3.3 连锁不平衡与单倍型分析
连锁不平衡（LD）分析显示，rs2236757与rs9984273之间（D′=0.0741，r2=0.0028）、rs9984273与rs1051393之间（D′=0.0524，r2=0.0012）均不存在LD；rs2236757与rs1051393之间存在微弱LD（D′=0.4648，r2=0.18，p<0.001）。单倍型分析识别出8种单倍型，整体分布差异显著（全局检验p=0.00018）。以最常见单倍型TGT为参照，TAC单倍型在住院组中频率较低（1.25% vs. 7.40%），名义上与危重症风险降低相关（OR=0.11，95%CI：0.02–0.86，p=0.036），但该比较经多重检验校正后不再显著，且基于极少数携带者。GAC单倍型在ICU患者中完全缺失（0%），在非住院组中出现频率为10.36%；经Haldane-Anscombe连续性校正后获得OR=0.03（95%CI：0.002–0.53，p=0.00001），但该估计依赖零细胞计数，仅为假设生成性质。
四、讨论与结论
讨论部分指出，本研究发现rs2236757的A等位基因和rs9984273的C等位基因在突尼斯人群中表现出保护效应，这与部分文献（如Safari等人的研究及Jalkanen等人的ARDS/COVID-19队列）一致，但与GenOMICC研究（英国GWAS，认为A等位基因增加风险）、巴勒斯坦队列（A等位基因与呼吸困难加重相关，OR=4.70）及墨西哥队列（A等位基因在死亡者中频率更高）的结果相反。研究人员详细分析了这种“翻转”关联的可能原因：不同人群间LD模式的差异、群体特异性基因-环境互作、表型定义不一致（如ICU入院vs.院内死亡）、研究设计差异（GWAS vs.候选基因、前瞻性vs.回顾性）以及统计效力不足。特别强调，本研究中rs2236757的等位基因编码遵循dbSNP正向链惯例，排除了链方向人为误差。对于rs2236757-A可能的功能机制，引用了相关研究指出该等位基因与更高的IFNAR2表达水平相关，从而可能增强抗病毒信号，但这在本研究中未经实验验证。rs9984273位于3′UTR区域，可能参与转录后调控，但其功能也未在本研究中评估。rs1051393作为错义突变（苯丙氨酸→缬氨酸），未发现关联，可能与统计效力不足有关。
研究人员坦承多项局限性：住院组仅46例，统计效力极低，易产生假阴性和假阳性；单中心招募可能导致选择偏倚，无法推广至全国；样本采集于2021年疫苗接种前，结果不适用于当前流行株或已接种人群；未进行人群分层校正；未系统性筛查继发感染；缺乏功能性验证。因此，所有结论均为假设生成性质，而非确认性结论。
研究结论部分重申：本候选基因关联研究提供了初步证据，表明IFNAR2的rs2236757和rs9984273与突尼斯人群COVID-19危重症风险降低相关，而rs1051393未见独立关联；TAC单倍型可能具有保护趋势，GAC单倍型因零细胞计数仅为假设生成。但由于样本量小、病例对照设计、单中心回顾性、缺乏功能验证及继发感染筛查，这些发现仅可作为未来研究的起点，需在更大规模、独立、多中心队列中重复验证，并结合功能研究阐明机制，方可考虑潜在的临床应用。
