【原】网页工具的差异分析和功能分析就发了4+分还小2区

Identification of age- and gender-associated long noncoding RNAs in the human brain with Alzheimer’s disease阿尔茨海默病患者大脑年龄及性别相关lncRNA的鉴定

一、研究背景

  阿尔茨海默病(AD)是一种与年龄及性别相关的脑部疾病。研究表明，长链非编码RNA(lncRNAs)在脑部发育、稳态和脑部病理中起到重要的调节作用。因此，鉴定年龄及性别相关lncRNA在AD患者脑中的差异表达不仅可以为进一步的机制研究提供靶点，还可以为AD患者的个体化治疗提供理论基础。

二、分析流程

1. limma包：矫正表达数据

2. GEO2R网页工具：差异分析

3. KEGG富集分析：网页工具(R2 KEGG pathway: http://r2.)

4. 其余：GraphPad Prism

下面给大家推荐2个网页工具，也可以代替上面的R包或GraphPad Prism进行分析：

IDE.91

(http://bioinformatics./idep/)

矫正表达数据和差异分析

易汉博

(http://www./ImageGP/index.php/Home/Index/index.html)

可视化的一大利器

三、结果解读

1.衰老大脑和AD患者大脑的lncRNA表达模式

• 分析正常年轻大脑（≤40岁）和衰老大脑（≥80岁）前额叶皮质的lncRNA表达差异：衰老大脑中22个lncRNA表达上调，22个lncRNA表达下调。

• 分析正常大脑和AD患者大脑的lncRNA表达差异：AD患者大脑中18个lncRNA表达上调，32个lncRNA表达下调。

• 正常衰老大脑和AD患者大脑中有15个表达发生变化的lncRNA重合。

• 重合的15个lncRNA在衰老大脑和AD患者大脑中表达改变模式相同，且LINC01094，NEAT1和LINC00844表达显著上调(Log FC>1)，LINC01007，LINC00507和OIP5-AS1表达显著下调(Log FC<-1)。（如图1）

图1.衰老大脑和正常大脑中重合lncRNA的表达改变模式

2.AD患者大脑中异常表达的年龄相关lncRNA

• 对年龄和lncRNA log2表达量进行Pearson相关分析：23个lncRNA表达量与年龄呈负相关，29个lncRNA表达量与年龄呈正相关（|r|>0.5,p<0.05）。

• 对正常大脑和(nd)和AD患者大脑lncRNA的平均log2表达值与年龄进行Pearson相关分析，同时提供95%的预测区间：在nd中与年龄显著相关的16个lncRNA(|r|>0.75,p<0.05)在AD脑部表达发生异常改变。

  如年龄正相关lncRNA LINC00844，PRR34-AS1在AD中表达异常上调，年龄负相关lncRNAMIR7-3HG，LINC00643在AD中表达异常下调。（如图2）

图2.年龄相关lncRNA在nd和AD额叶皮质的表达量改变比较

3.AD患者大脑中异常调节的性别相关lncRNA

• 分析GSE48350数据集的nd和AD样本（如图3A）：正常性别差异lncRNA中，有3个lncRNA在女性AD中差异表达，11个lncRNA在男性AD中差异表达；患者性别差异lncRNA中，有3个lncRNA在女性nd中差异表达，14个lncRNA在男性nd中差异表达。

• 与nd相比，13个性别相关lncRNA在男性和女性AD中表达均发生了改变（黑点为男性，灰点为女性）。

  双向ANOVA分析（性别和AD两个变量）发现性别因素在SLC25A25-AS1和RP11-45P15.4的表达中起重要作用，表明AD对性别相关lncRNA表达的影响依性别不同而不同。（图3B、E）

• 用CPC和CPAT分析lncRNA的编码潜能，发现AD中表达异常的这16个年龄相关lncRNA和13个性别相关lncRNA均不编码RNA。

图3.性别相关lncRNA在AD患者中表达异常情况

4.年龄及性别相关lncRNA与AD Braak分期及宿主基因的相关性

• 对GSE48350数据集的lncRNA log2平均表达值和Braak分期进行Pearson相关性分析，研究lncRNA表达与疾病分期的关系：年龄相关lncRNA LINC00672与Braak分期负相关，年龄相关lncRNA SNHG19与Braak分期正相关；性别相关lncRNA RNF144A-AS1，LY86-AS1和LINC00639）与Braak分期负相关。

• 研究lncRNA表达是否与AD的重要组织病理缺损之一——神经炎性斑块（NP）密度有关，计算年龄及性别相关lncRNA log2表达量和平均NP密度的相关性，发现有4个lncRNA（RP11-223I10.1，LINC01094，OIP5-AS1和LINC01007）与平均NP密度呈中等程度线性相关(|r|>0.3,p<0.05,n=102)。

5.AD患者神经元中的年龄及性别相关lncRNA的调节异常

• 用含AD患者和正常人额上回神经元基因表达谱的GSE5281数据集，检测lncRNA是否是神经元特异性：发现大多数lncRNA在nd和AD患者中的表达改变模式相同。

• 年龄负相关lncRNA LY86-AS1，RP11-219A15.5，MIR7-3HG，LINC01007和LINC00643在AD患者神经元中显著下调；年龄正相关lncRNA NEAT1，PRR34-AS1和lnc-BLID-2在AD患者神经元中显著上调。

• 年龄负相关lncRNA OIP5-AS1在AD患者神经元中显著上调，表明OIP5-AS1在不同的中枢神经系统细胞中改变模式不同，如在AD神经元中上调，在AD神经胶质细胞中下调。

• AD相关lncRNA反映可以反映AD早期的氧化损伤和DNA损伤修复(DDR)水平。用GSE66333数据集评估AD相关lncRNA表达与DDR强度的关系发现，AD早期的高水平氧化应激和神经元DDR可以诱发AD相关lncRNA的下调。

6.年龄及性别相关lncRNA的功能分析

• 用R2 KEGG pathway finder对lncRNA相关基因进行通路富集分析，研究lncRNA在AD大脑中的潜在功能并列举出前30个显著富集的通路。（图4）

• 年龄及性别相关lncRNA相关基因在AD，HD，PD等神经退行性疾病中均显著富集。

• AD的性别及年龄相关基因均在胞吞和前列腺癌通路中显著富集。

• 正常大脑的lncRNA相关基因在蛋白酶体和长时程增强KEGG通路中均富集。

  

  

图4.年龄及性别相关lncRNA相关基因富集通路（示前30个通路）

小结