CN112218958A - 估计皮肤内在年龄的表观遗传学方法 - Google Patents
估计皮肤内在年龄的表观遗传学方法 Download PDFInfo
- Publication number
- CN112218958A CN112218958A CN201980035763.0A CN201980035763A CN112218958A CN 112218958 A CN112218958 A CN 112218958A CN 201980035763 A CN201980035763 A CN 201980035763A CN 112218958 A CN112218958 A CN 112218958A
- Authority
- CN
- China
- Prior art keywords
- skin
- genomic dna
- age
- sites
- intrinsic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000001973 epigenetic effect Effects 0.000 title description 10
- 210000003491 skin Anatomy 0.000 claims abstract description 59
- 230000030933 DNA methylation on cytosine Effects 0.000 claims abstract description 7
- 210000004927 skin cell Anatomy 0.000 claims abstract description 4
- 108020004414 DNA Proteins 0.000 claims description 25
- 239000000523 sample Substances 0.000 claims description 25
- 108091028043 Nucleic acid sequence Proteins 0.000 claims description 10
- 239000012472 biological sample Substances 0.000 claims description 4
- LSNNMFCWUKXFEE-UHFFFAOYSA-M Bisulfite Chemical compound OS([O-])=O LSNNMFCWUKXFEE-UHFFFAOYSA-M 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 239000003153 chemical reaction reagent Substances 0.000 claims description 2
- 238000009396 hybridization Methods 0.000 claims description 2
- 238000006116 polymerization reaction Methods 0.000 claims description 2
- 238000012175 pyrosequencing Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000032683 aging Effects 0.000 description 17
- 230000011987 methylation Effects 0.000 description 15
- 238000007069 methylation reaction Methods 0.000 description 15
- 238000012360 testing method Methods 0.000 description 15
- 108091029430 CpG site Proteins 0.000 description 13
- 238000012549 training Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 9
- 230000014509 gene expression Effects 0.000 description 8
- 230000007067 DNA methylation Effects 0.000 description 5
- 230000001681 protective effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 230000004224 protection Effects 0.000 description 4
- 108090000623 proteins and genes Proteins 0.000 description 4
- 210000000434 stratum corneum Anatomy 0.000 description 4
- 206010051246 Photodermatosis Diseases 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 210000002615 epidermis Anatomy 0.000 description 3
- 210000000245 forearm Anatomy 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000003908 quality control method Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000009759 skin aging Effects 0.000 description 3
- 108010077544 Chromatin Proteins 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 210000003483 chromatin Anatomy 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 210000004207 dermis Anatomy 0.000 description 2
- 230000037213 diet Effects 0.000 description 2
- 235000005911 diet Nutrition 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000019975 dosage compensation by inactivation of X chromosome Effects 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 230000030279 gene silencing Effects 0.000 description 2
- 238000012226 gene silencing method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000004060 metabolic process Effects 0.000 description 2
- 238000002493 microarray Methods 0.000 description 2
- 230000008845 photoaging Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000007388 punch biopsy Methods 0.000 description 2
- 238000007390 skin biopsy Methods 0.000 description 2
- 230000036555 skin type Effects 0.000 description 2
- 230000035882 stress Effects 0.000 description 2
- 230000036561 sun exposure Effects 0.000 description 2
- 230000037316 sun-exposed skin Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000004543 DNA replication Effects 0.000 description 1
- 238000001712 DNA sequencing Methods 0.000 description 1
- 240000008168 Ficus benjamina Species 0.000 description 1
- 206010063493 Premature ageing Diseases 0.000 description 1
- 208000032038 Premature aging Diseases 0.000 description 1
- 108700025716 Tumor Suppressor Genes Proteins 0.000 description 1
- 102000044209 Tumor Suppressor Genes Human genes 0.000 description 1
- 108700029631 X-Linked Genes Proteins 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000003679 aging effect Effects 0.000 description 1
- 230000004071 biological effect Effects 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000033077 cellular process Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000013020 embryo development Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011985 exploratory data analysis Methods 0.000 description 1
- 238000013100 final test Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 229940029329 intrinsic factor Drugs 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000002690 local anesthesia Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011325 microbead Substances 0.000 description 1
- 238000013188 needle biopsy Methods 0.000 description 1
- 238000000059 patterning Methods 0.000 description 1
- 230000003711 photoprotective effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000009290 primary effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013442 quality metrics Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000037380 skin damage Effects 0.000 description 1
- 230000000392 somatic effect Effects 0.000 description 1
- 238000013518 transcription Methods 0.000 description 1
- 230000035897 transcription Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6813—Hybridisation assays
- C12Q1/6834—Enzymatic or biochemical coupling of nucleic acids to a solid phase
- C12Q1/6837—Enzymatic or biochemical coupling of nucleic acids to a solid phase using probe arrays or probe chips
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/154—Methylation markers
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Zoology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Microbiology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
Description
技术领域
本发明涉及检测和分析基因组DNA中胞嘧啶甲基化模式的方法。更具体地,本发明涉及检测和分析基因组DNA中特定位点的胞嘧啶甲基化模式以确定皮肤内在年龄和健康。
背景技术
众所周知,老化是主要由个体年龄驱动的多因素过程。皮肤老化是由内在和外在因素两者驱动的特别的多因素现象。就内在因素而言,个体的实足年龄是最为熟知的,但其他内在因素如个体的代谢、饮食、应激和潜在健康也对皮肤年龄有影响。除了这些内在因素之外,皮肤暴露于外部挑战如UV辐射、污染、干燥条件和极端温度。这些外在因素因此也影响个体皮肤年龄。
因此,显然存在两种不同形式的皮肤年龄:外在年龄,其由外在因素(即,源于角质层外表面的外部,然后其穿透角质层进入皮肤)引起的老化累积主导,尤其是由日光暴露主导(光老化);以及内在年龄,其为由于内源产生的因素所致的皮肤老化程度;换言之,不是由于外在因素所致的老化。为了理解起见,考虑个体的两种不同类型的皮肤是有帮助的。一种来自通常由衣服保护的部位(如臀部区域或上臂内侧区域)。另一种来自日光暴露部位(如面部或手背)。保护部位会极少暴露于外在老化因素,因此任何老化会归因于内在因素。暴露部位会完全暴露于外在老化因素,因此该区域老化的年龄会归因于内在因素导致的固有内在年龄和外在因素导致的老化两者的组合。
本发明涉及开发估计个体皮肤内在年龄的表观遗传学方法。
DNA甲基化是基因表达的表观遗传学决定因素。CpG甲基化模式是可遗传的、组织特异性的,并且与基因表达相关。甲基化(特别是如果位于基因启动子中)的结果通常是基因沉默。DNA甲基化还与其他细胞过程相关,所述细胞过程包括胚胎发育、染色质结构、基因组印记、雌性中的体细胞X染色体失活、外源DNA的转录和转座的抑制以及DNA复制时序(timing of DNA replication)。当基因被高度甲基化时,其不太可能被表达。因此,鉴别基因组中含有5-meC的位点对于了解基因表达的细胞类型特异性程序以及了解在正常发育、老化和疾病如癌症期间基因表达谱如何改变是重要的。DNA甲基化模式的作图(mapping)对于了解多样生物学过程如印记基因的调节、X染色体失活和人类癌症中肿瘤抑制基因沉默是重要的。
Horvath S.等人的“DNA methylation age of human tissues and cell types”(Genome Biology 14(2103)R115)报道了使用实足年龄的转化版本,其使用惩罚回归模型(弹性网)在CpGs上回归。弹性网回归模型选择被称为表观遗传学时钟CpGs的353个CpGs,因为它们的加权平均值(由回归系数形成)据称相当于表观遗传学时钟。该研究在本专利中称为“Horvath研究”。
然而,我们现在已经发现,对于日光暴露皮肤部位,基于这些353个基因座的预测年龄比它们的实际(“实足”)年龄大约年轻9岁,表明它们没有检测日光诱导的皮肤损伤。另外,发现防晒皮肤样品的年龄比实足年龄年轻4岁,这是对防晒皮肤年龄的低估,预期防晒皮肤年龄与从其采集样品的受试者的实足年龄大致相同。因此,这些353个基因座不能识别光损伤和光保护皮肤类型之间的差异,低估了防晒皮肤的年龄,并且预测光损伤皮肤比光保护皮肤更年轻。因此,可以理解,该模型不能评估不同形式的老化-外在老化和内在老化。
因此,本发明旨在解决这种现有技术老化模型的不良性能,并提供用于评估皮肤内在年龄的改进方法。
发明内容
我们已惊奇地发现,不同的、特定的甲基化位点集为内在皮肤年龄的预测提供提高的准确度。特别地,该位点能够预测保护皮肤的年龄,并且还能够给出令人惊奇地不受外在因素影响的暴露皮肤内在年龄。
因此,在第一方面,本发明提供用于获得可用于确定个体的皮肤内在年龄的信息的方法,所述方法包括以下步骤:
(a)从来源于所述个体的皮肤细胞获得基因组DNA;以及
(b)观察所述基因组DNA中选自下组的>30个CpG基因座的胞嘧啶甲基化:
从而获得可用于确定所述个体的所述皮肤内在年龄的信息。
基因组DNA获自来源于个体的皮肤细胞。皮肤样品优选包含表皮,所述表皮是单独的或与真皮组合。
优选使用来自该组的>40个位点,更优选使用该组的>45、>50、>55、>60、>65、>70、>75、>80、>85个位点,最优选所有89个位点。
优选所观察的基因座是:
更优选所观察的基因座是:
在另一实施方案中,评估基因组DNA中胞嘧啶甲基化,其中所述基因组DNA在以上列出的CpG基因座名称的20kBp内,优选在15kBp内,更优选在10kBp内,还更优选在5kBp内,甚至更优选在1kBp内,最优选在0.5kBp内。
在第二方面,本发明提供用于获得可用于确定个体的皮肤内在年龄的信息的试剂盒,所述试剂盒包含:
-特异于生物样品中的>30个基因组DNA序列的引物或探针,其中所述基因组DNA序列在所述基因组DNA中包含选自下组的CpG基因座,所述组仅由以下CpG基因座名称组成:
用于以下方法的试剂:
基因组DNA聚合方法;
基因组DNA杂交方法;
基因组DNA直接测序方法;
基因组DNA亚硫酸氢盐转化方法;或
基因组DNA焦磷酸测序方法。
优选地,引物或探针特异于生物样品中的>40个基因组DNA序列,更优选>45、>50、>55、>60、>65、>70、>75、>80、>85个,最优选引物或探针特异于该组中的所有89个位点。
优选地,引物或探针特异于皮肤样品中的基因组DNA序列,最优选地,皮肤样品包含表皮,所述表皮是单独的或与真皮组合。
优选地,引物或探针特异于以下CpG基因座名称:
更优选地,引物或探针特异于以下CpG基因座名称:
在另一实施方案中,评估基因组DNA中的胞嘧啶甲基化,其中所述基因组DNA在以上列出的CpG基因座名称的20kBp内,优选在15kBp内,更优选在10kBp内,还更优选在5kBp内,甚至更优选在1kBp内,最优选在0.5kBp内。
优选地,试剂盒包含甲基化微阵列。
优选地,试剂盒包含DNA测序方法。
具体实施方式
如所讨论的,皮肤中的老化过程是也在身体上变化的高度多因素现象。例如,保护皮肤比暴露皮肤受到少得多的损伤,并且因此显然来自同一个体的不同皮肤区域会具有不同程度的损伤,因此具有不同的“年龄”。
在本发明中,我们考虑两种形式的皮肤年龄:内在年龄;和外在年龄。
就内在年龄而言,个体的实足年龄是主导性的,但其他内在因素如个体的代谢、饮食、应激和潜在健康也有影响。因此,在本发明的上下文中,内在年龄意指由内源因素导致的皮肤年龄。
就外在年龄而言,固有年龄仍会是基本组成部分,但此外,外源因素如UV辐射、污染、干燥条件和极端温度也会有影响。因此,在本发明的上下文中,外在年龄意指主要由外源因素导致的皮肤年龄。
为了清楚起见:外在年龄由外在因素(即,源于角质层外表面的外部,然后其穿透角质层进入皮肤)引起的老化累积主导,尤其是由日光暴露主导(光老化);而内在年龄是由于内源产生的因素引起的皮肤老化程度;换言之,不是由于外在因素所致的老化。
本发明涉及开发估计个体皮肤内在年龄的表观遗传学方法。
数据集
本申请使用三个表观遗传学数据集。
-鉴别:第一数据集用于鉴别与皮肤中保护部位和暴露部位相关的甲基化位点。
-训练:第二数据集用于训练数学模型,其中评估从鉴别数据集鉴别的甲基化位点,确定最能够预测皮肤年龄的甲基化位点,并且构建预测模型。
-测试:最后,第三测试数据集用于评估这些甲基化位点在确定皮肤样品年龄中的准确度,以及这些甲基化位点的使用是否比在Horvath研究中鉴别的那些更准确。
第一数据集(鉴别)是涉及24名中国女性和24名高加索女性参与者的单中心、横断面活检研究,其中已经登记了24名年轻女性和24名老年女性。从每个受试者的两个不同区域收集皮肤样品:来自皮肤暴露区域的样品;和来自皮肤保护区域的样品。指定为暴露的部位位于下臂外侧。保护部位位于上臂内侧,通常位于肘和腋窝区域之间的一半处。
第二训练数据集(训练)是公众可获得的数据集(Bormann F.等人:Reduced DNAmethylation patterning and transcriptional connectivity define human skinaging.Aging Cell(2016)1-9.Array express id:EMTAB-4385)。数据集包括总共108个表皮样品,48个样品已经从穿刺活检分离,所述穿刺活检是从24名年轻人(18-27岁)和24名老年人(61-78岁)的前臂外侧获得的。从60名20-79岁的志愿者的前臂外侧获得泡罩顶(suction blister roof)形式的60个样品。所有志愿者均为女性、高加索人,且无病。
最终测试数据集(测试)是公众可获得的数据集(Vandiver A.R.等人:Age andsun exposure-related widespread genomic blocks of hypomethylation innonmalignant skin.Genome Biology(2015)16:80)Gene Expression Omnibus accessionnumber:GSE51954)。数据集包括来自20名高加索受试者的表皮样品(N=38)。在局部麻醉下,从前臂外侧或侧向内眦赘皮(lateral epicanthus)(暴露区域)和上臂内侧(保护区域)收集成对的直径4mm的钻孔活检(punch biopsy)样品。
训练数据集和测试数据集的选择
由以下标准指导数据集的选择。首先,训练和测试数据需要来自表皮皮肤(皮肤活检或者仅表皮)。所选的训练数据(Bormann等人)来自皮肤活检和前臂外侧的泡罩,并且表皮样品可用于测试(Vandiver等人)数据集。其次,训练数据需要在连续年龄上,而测试数据需要在年轻和年老的年龄组中都具有暴露样品和保护样品。第三,训练数据集中的平均年龄(47岁,标准偏差=21)需要并且与测试数据集的平均年龄(51岁,标准偏差=25岁)相当。
甲基化数据质量检查
所有三个数据集均使用杂交至Infinium 450k人甲基化微珠芯片的亚硫酸氢盐转化DNA。
鉴别数据集中所有DNA样品的甲基化数据均通过基于三个阵列质量度量(MAplot,Boxplot,Heatmap)的质量检查。β-值计算为B=R/R+G,并且M-值计算为M=log2(R/G),其中R表示甲基化信号,而G表示未甲基化信号。将60的偏移量加至分母。M-值用于创建表达矩阵。使用分位数归一化对原始数据进行归一化。将β-值用于随后的建模和筛选统计结果。
对‘minfi’R软件包中的raw.idat文件进行训练数据集的质量控制和预处理。使用阵列归一化内的子集分位数(Subset-quantile Within Array Normalization,SWAN)对原始数据进行归一化。
对于测试数据集,执行SWAN所需的raw.idat文件不可用。因此,将提供的经Illumina预处理的β值用于后续分析。还使用‘minfi’R软件包进行应用于数据的质量控制和预处理。
对数据的技术影响
对鉴别数据集使用主成分分析(PCA)进行探索性分析。发现阵列间重复(between-array replicate)不聚集在一起,可能是由于与阵列编号相关的分批效应(batcheffect)。测试数据集的聚集分析揭示了类似的阵列分批效应。在培训数据集上未看到技术分批效应。
分批效应校正数据
使用ComBat方法(Johnson W.E.等人:Adjusting batch effects in microarrayexpression data using empirical Bayes methods.Biostatistics 8(1)(2007)118-127),按照质量控制、归一化和阵列内重复的平均,对在鉴别数据集和测试数据集中观察到的阵列分批效应进行调整。在分批校正之后得到的数据集显示在阵列上没有聚集。其余生物效应仍然存在,并倾向于成为数据中的主要效应。
CpG基因座鉴别
如本文所用,CpG基因座是指在Illumina CpG基因座数据库中找到的独特标识符(如Technical Note:Epigenetics,CpG Loci Identification ILLUMINAInc.2010,https://www.illumina.com/documents/products/technotes/technote_cpg_loci_identification.pdf中所述)。这些CpG位点标识符因此提供一致性和确定性的CpG基因座数据库以确保甲基化数据报告的一致性。
Horvath表观遗传学时钟在预测防晒皮肤年龄方面的性能
来自Horvath研究(其使用以上讨论的353个CpG位点)的年龄预测值针对测试数据集的暴露(Se)样品和保护(Sp)样品运行。使用线性回归评估Horvath模型的性能,从线性回归获得R2(“pho”或“ρ”)。还计算中值误差(预测年龄相对于实际年龄)。结果在表1中提供。
*Se和Sp样品未配对。暴露受试者的年龄为34岁,保护受试者的年龄为30岁。
表1–使用来自Horvath研究的预测值预测的暴露皮肤样品和保护皮肤样品年龄
可以看到,对于19个受试者中的15个受试者,Horvath模型计算的暴露样品比保护样品更年轻,这是不正确的,因为预期经受暴露如UV辐射的样品比免于UV损害的样品更老。
对预测年龄的平均年龄加速显示日光暴露皮肤样品比实足年龄年轻9岁,这与暴露(尤其是日光暴露)引起皮肤过早老化的已知生理学相反。此外,对于仅与内在老化相关的模型,预期这会为保护样品和暴露样品给出大致相似的年龄。
另外,发现保护皮肤样品的年龄比实足年龄年轻4岁,这是对保护皮肤年龄的低估,预期保护皮肤年龄与从其采集样品的人的实足年龄大致相同。
因此可以得出结论,来自Horvath研究的353个CpG位点不能识别暴露皮肤类型和保护皮肤类型之间的差异,也不能识别暴露皮肤中的内在老化效应,不正确地预测晒伤皮肤比防晒皮肤更年轻,并且低估了保护样品的年龄。
还发现通过Horvath研究鉴别的353个CpG位点就保护样品的准确度评分而言性能较差。
保护样品的准确度评分为:
ρ=0.93(误差=16.6岁)。
因此可以理解,需要用于确定皮肤内在年龄的改进的表观遗传学方法。
与保护部位相关的甲基化位点的鉴别(来自鉴别数据集)
为了产生外在年龄列表和内在年龄列表,使用不同线性模型对归一化的分批校正数据进行总共5个比较(表2)。应用在多个测试校正列表(调整P值-adjP,benjaminiHochberg)<0.05的统计截断设置以及δ-β>=0.05。
对于暴露部位年轻相对于年老的比较(比较1:n=10,649),检测大量差异甲基化CpG位点。对于年龄组相对于部位相互作用的比较(比较5:n=233),鉴别相对较少的差异甲基化CpG位点。
表2–统计结果。5个比较中的每一个的差异甲基化位点的数目,其中调整的截断值p为0.05
内在部位列表
为了鉴别仅捕获内在老化的CpG位点,筛选比较2(年轻相对于年老保护部位)结果以去除在年轻或年老中因部位改变的探针(比较3和4),以去除保护皮肤中可能另外受外在因素影响的任何老化变化。
得到的列表是1,575个CpG位点。对这些1,575个位点的PCA分析允许鉴别在将保护部位归类为跨越两个种族的年轻组和年老组中贡献最大差异的位点。使用PCA载荷来选择这些可变探针,施加到第一组分上的0.030载荷的截断值导致捕获年老组之间的最大变异性的322个探针。
保护部位的内在年龄预测值
从鉴别数据集鉴别以捕获内在年龄变化的322个CpG位点用于构建内在年龄模型,其中将与Horvath研究中使用的相同的弹性网用于10组n/10大小的训练数据集(对9个数据集进行训练和对1个数据集进行测试)。将这些重复10次,并且获得每次迭代的平均“准确度”以给出用于计算年龄的模型和针对每个探针的系数。
通过运行模型的几次迭代得到预测值列表。第一次迭代鉴别最佳预测值集。对于每个后续迭代,从训练集中排除来自先前迭代的所鉴别的预测值,以鉴别下一最佳预测值集。重复迭代,直到发现以ρ和误差裕度测量的预测准确度不如上述Horvath模型的预测准确度。
对于内在位点,进行了3次迭代。如表3所示,第一次鉴别了36个位点,第二次鉴别了53个位点,第三次鉴别了25个位点。
使用所得模型(其中从322个CpG位点的最终内在年龄列表中去除来自这3次迭代中的每一次迭代的位点)来估计来自测试数据集的保护样品的年龄。结果示于表4中。此外,计算得到的模型的防晒样品和日光暴露样品的平均年龄。结果示于表5中。表4和表5(斜体)中还显示了使用来自Horvath研究的353个位点预测防晒年龄的模型的准确度,以供参考。
表3–内在年龄评分的预测值集
表4–模型的准确度
根据表4中所示的准确度量度,与使用353个Horvath位点的模型(R2=0.93,误差=16.6岁)相比,包括在迭代1和迭代2中鉴别的位点的内在年龄的模型以更高或相等的准确度(R2=0.96,误差=5.7岁和R2=0.94,误差=12.6岁)和更好的误差执行。其余208个位点(其包括来自迭代3的25个位点)的测量准确度低于353个Horvath位点。因此,与Horvath模型相比,迭代1和迭代2的89个位点在预测内在年龄方面更佳。
表5–模型的平均年龄
预期来自日光暴露部位的样品的内在年龄会类似于来自防晒部位的样品的内在年龄。从表5可以看到,与Horvath模型相比,本研究的模型在日光暴露部位和防晒部位的平均年龄之间具有较小的差异。这表明本文描述的模型在预测内在年龄方面优于Horvath模型。
因此可以看到,使用选自表3中所示的迭代1和迭代2的CpG位点在确定皮肤内在年龄时提供更好的准确度。因此,本发明提供这些89个位点中>30个用于预测皮肤内在年龄。本发明还提供迭代2的53个位点作为优选组。本发明还提供迭代1的36个位点作为最优选组。
本发明的可选方案是,上述CpG位点也可以被替换,并且替代地使用最接近的基因。
表6提供了在迭代1和迭代2中鉴别的105个位点的注释(如Price等人,Epigenetics&Chromatin 2013,6:4,“Additional annotation enhances potential forbiologically-relevant analysis of the Illumina Infinium HumanMethylation450BeadChip array”中所述,使用Human Genome version HG19),包括最接近的基因名称。
表6–在迭代1和迭代2中鉴别的CpG位点的注释
Claims (9)
2.根据权利要求1所述的方法,其中使用来自所述组的>40个位点,更优选>45、>50、>55、>60、>65、>70、>75、>80、>85个位点,最优选所有89个位点。
6.根据权利要求5所述的试剂盒,其中所述引物或探针特异于生物样品中的>40个基因组DNA序列,更优选>45、>50、>55、>60、>65、>70、>75、>80、>85个基因组DNA序列,最优选所有89个基因组DNA序列。
7.根据权利要求5或6所述的试剂盒,其中所述引物或探针特异于皮肤样品中的基因组DNA序列。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP18177976 | 2018-06-15 | ||
EP18177976.0 | 2018-06-15 | ||
PCT/EP2019/065709 WO2019238935A1 (en) | 2018-06-15 | 2019-06-14 | Epigenetic method to estimate the intrinsic age of skin |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112218958A true CN112218958A (zh) | 2021-01-12 |
Family
ID=62684682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201980035763.0A Pending CN112218958A (zh) | 2018-06-15 | 2019-06-14 | 估计皮肤内在年龄的表观遗传学方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20210207201A1 (zh) |
EP (1) | EP3775275A1 (zh) |
CN (1) | CN112218958A (zh) |
MX (1) | MX2020012703A (zh) |
WO (1) | WO2019238935A1 (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150259742A1 (en) * | 2012-11-09 | 2015-09-17 | The Regents Of The University Of California | Methods for predicting age and identifying agents that induce or inhibit premature aging |
CN105765083A (zh) * | 2013-09-27 | 2016-07-13 | 加利福尼亚大学董事会 | 基于表观遗传学标记物来估计组织和细胞类型的年龄的方法 |
-
2019
- 2019-06-14 CN CN201980035763.0A patent/CN112218958A/zh active Pending
- 2019-06-14 WO PCT/EP2019/065709 patent/WO2019238935A1/en unknown
- 2019-06-14 EP EP19729570.2A patent/EP3775275A1/en active Pending
- 2019-06-14 MX MX2020012703A patent/MX2020012703A/es unknown
- 2019-06-15 US US17/057,771 patent/US20210207201A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150259742A1 (en) * | 2012-11-09 | 2015-09-17 | The Regents Of The University Of California | Methods for predicting age and identifying agents that induce or inhibit premature aging |
CN105765083A (zh) * | 2013-09-27 | 2016-07-13 | 加利福尼亚大学董事会 | 基于表观遗传学标记物来估计组织和细胞类型的年龄的方法 |
Non-Patent Citations (1)
Title |
---|
VANDIVER AR 等: "Age andsun exposure-related widespread genomic blocks of hypomethylation in nonmalignant skin", GENOME BIOL, vol. 16, no. 1 * |
Also Published As
Publication number | Publication date |
---|---|
EP3775275A1 (en) | 2021-02-17 |
WO2019238935A1 (en) | 2019-12-19 |
US20210207201A1 (en) | 2021-07-08 |
MX2020012703A (es) | 2021-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3899018B1 (en) | Cell-free dna end characteristics | |
CN108138233B (zh) | Dna混合物中组织的单倍型的甲基化模式分析 | |
US20200190568A1 (en) | Methods for detecting the age of biological samples using methylation markers | |
CA3180334A1 (en) | Methods for detection of donor-derived cell-free dna | |
CN108256292B (zh) | 一种拷贝数变异检测装置 | |
KR102334702B1 (ko) | 유방암 발생 위험의 평가 방법 | |
JP2015526101A (ja) | 染色体異常を検出する方法 | |
JP2023521308A (ja) | 合成トレーニングサンプルによるがん分類 | |
JP2020536585A (ja) | 乳癌発症リスクの評価法 | |
CN112204153A (zh) | 估计皮肤外在年龄的表观遗传学方法 | |
CN117079723B (zh) | 一种与肌萎缩侧索硬化症相关的生物标志物、诊断模型及其应用 | |
CN113201590B (zh) | 用于评估肝细胞癌早期复发风险的lncRNA、评估方法及装置 | |
CN112218958A (zh) | 估计皮肤内在年龄的表观遗传学方法 | |
US20180247019A1 (en) | Method for determining whether cells or cell groups are derived from same person, or unrelated persons, or parent and child, or persons in blood relationship | |
CN106119406B (zh) | 多发性肉芽肿血管炎及微小动脉炎的基因分型诊断试剂盒及使用方法 | |
EP3854881A1 (en) | Classifying subjects based on their biological response to uv irradiation | |
Nakabayashi | The Illumina Infinium methylation assay for genome-wide methylation analyses | |
CN115074439B (zh) | 一组nk/t细胞淋巴瘤预后相关基因、基因组预后模型及其用途 | |
US20240021267A1 (en) | Dynamically selecting sequencing subregions for cancer classification | |
US20240233872A9 (en) | Component mixture model for tissue identification in dna samples | |
US20240170099A1 (en) | Methylation-based age prediction as feature for cancer classification | |
US20220101947A1 (en) | Method for determining fetal fraction in maternal sample | |
KR20230072354A (ko) | 레티놀을 포함하는 처방의 자극 유무 예측 방법 및 시스템 | |
KR20230088246A (ko) | 비타민 c를 포함하는 처방의 피부 밝기 변화 예측 방법 및 시스템 | |
KR20220141658A (ko) | 피부색 판단용 유전자 다형성 마커 및 이의 용도 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20210907 Address after: Rotterdam Applicant after: Unilever Intellectual Property Holdings Ltd. Address before: Rotterdam Applicant before: Unilever Nederland B.V. |
|
TA01 | Transfer of patent application right |