Name: mapping alzheimer's disease variants to their target genes using computational analysis of chromatin configuration
Uploaded: 2020-01-09T13:00:00
Duration: 4 min 41 s
Description: Watch this Scientific Journal Video about Mapping Alzheimer's Disease Variants to Their Target Genes Using Computational Analysis of Chromatin Configuration at JoVE.com

这项工作得到了NIH授予R00MH113823（至H.W.）和R35GM128645（D.H.P.）、NARSAD青年调查员奖（授予H.W.）和西蒙斯基金会自闭症研究倡议（SFARI，至N.M.和H.W.）的SPARK赠款的支持。

1. 工作站设置
<ol><li>安装 R（3.5.0 版）和 RStudio 桌面。打开 RStudio。</li>
	<li>通过在 RStudio 中的控制台窗口中键入以下代码，在 R 中安装以下库。 如果（！BiocManager" % in% 行名（已安装.包（））） 安装.软件包（"生物管理"，存储库="https：//cran.r-project.org"） 生物管理器：：安装（"基因组学"） 生物经理：安装（"生物母体"） 生物经理：安装（"WGCNA"） 安装.软件包（"重塑"） 安装.软件包（"ggplot2"） 安装.软件包（"程序"） 安装.软件包（"gProfileR"） 安装.软件包（"整洁"） 安装.软件包（"ggpubr"）</li>
	<li>下载文件。 注： 在此协议中，所有文件都需要下载到 */工作目录。<ol>
<li>通过单击材料表中提供的链接下载以下文件。<ol>
<li>下载用于 AD 的精细映射可信 SNPs（Jansen 等人的第14号补充表 8 ）。 注： 分析之前，在 41588_2018_311_MOESM3_ESM.xlsx 中打开工作表 8，删除前三行，并将工作表另存为 Supplementary_Table_8_Jansen.txt，格式为 tab 分隔格式。</li>
			<li>从心理代码下载成人大脑中 10 kb 分辨率的 Hi-C 交互配置文件（如下所述anchored_chromatin_loops. 注： 此文件具有以下格式：染色体、TSS_start、TSS_end、Enhancer_start和Enhancer_end。在使用其他 Hi-C 数据集的情况下，此协议要求以高分辨率 （5~20 kb） 处理 Hi-C 数据集。</li>
			<li>从 PsychENCODE 下载单个单元格表达式数据集。 注：这些来自神经典型控制样本。</li>
			<li>从 BrainSpan 下载开发表达式数据集（如下所述为devExpr.rda）。 注：267666527是一个压缩文件，因此解压缩267666527以提取"columns_metadata.csv"、"expression_matrix.csv"和"rows_metadata.csv"以生成devExpr.rda（见第3节）。</li>
		</ol>
</li>
		<li>从Gencode版本19下载外音坐标（参见补充文件，描述为Gencode19_exon.bed和Gencode19_promoter.bed）。 注： 启动器定义为转录起始站点 （TSS） 上游的 2 kb。这些文件具有以下格式：染色体、开始、结束和基因。</li>
		<li>从生物市场下载基因注释文件（参见补充文件，如下所述为基因Anno.rda）。 注：此文件可用于匹配基于Ensembl基因ID和HUGO基因命名委员会（HGNC）符号的基因。</li>
	</ol>
</li>
</ol>2. 生成用于可信 SNPs 的 GRANGES 对象
<ol><li>通过在 RStudio 中的控制台窗口中键入以下代码在 R 中设置。 库（基因组学） 选项（字符串 AsFactors = F） setwd（"*/工作"） = 这是通往工作目录的路径。 criSNP = read.delim（"Supplementary_Table_8_Jansen.txt"，标题=T） 信元SNP = 信元SNP_credSNP$可信.因Causal="是"，]</li>
	<li>通过在 RStudio 中的控制台窗口中键入以下代码来创建 GRanges 对象。 信任 = GRanges（克里德SNP$Chr，IRanges（克里德SNP$bp，credSNP$bp），rsid_credSNP$SNP，P_credSNP$P） 保存（信任，文件="AD_credibleSNP.rda"）</li>
</ol>3. 位置映射
注： 对于每个步骤，在 RStudio 中的控制台窗口中键入相应的代码。
<ol><li>在 R 中设置。 选项（字符串AsFactors_F） 库（基因组学） 负载（"AD_credibleSNP.rda"） = （参见2）</li>
	<li>启动子/外置SNPs与基因的位置映射<ol>
<li>加载启动子和外音区域，并生成 GRange 对象。 外扰子 = 读取表（"Gencode19_exon.bed"） 外向范围 = GRanges（exon_，1]，IRanges（exon_，2]，exon_，3_），基因[exon]，4]） 启动器 = 读取表（"Gencode19_promoter.bed"） 促进者 = GRanges（促进者，1]，IRanges（促进者，2]，促进者[，3]），基因[促进者]，4]）</li>
		<li>将可信的 SN 与外音区域重叠。 olap = 查找重叠（信任范围、外子范围） 信元 = 信元[查询Hits（olap）] mcols（克瑞克森） = cbind（mcols（克瑞克森），mcols（外兆[主题Hits（olap）]）</li>
		<li>将可信的 SN 与促进者区域重叠。 olap = 查找重叠（信任范围、促进器） 信任发起人 = 信任范围[查询Hits（olap）] mcols（信任者） = cbind（mcols（信任发起人），mcols（促进者[主题Hits（olap）]）</li>
	</ol>
</li>
	<li>使用染色质相互作用将SNPs与其假定目标基因联系起来。<ol>
<li>加载 Hi-C 数据集并生成 GRange 对象。 hic = read.table（"发起人anchored_chromatin_loops.bed"，跳过=1） 名字（hic）= c（"chr"，"TSS_start"、"TSS_end"、"Enhancer_start"、"Enhancer_end"） hicranges = GRanges（hic$chr，IRanges（hic$TSS_start，hic$TSS_end），增强器=hic$Enhancer_start） olap = 查找重叠（奇数、启动子） hic促进者 = hicranges_查询Hits（olap）] mcols（活动发起人） = cbind（mcols（活动发起人），mcols（促进器[主题Hits（olap）]） hic增强剂 = GRanges（seqnames（hicpromoter），IRanges（hicpromoter$增强剂，hic促进器$增强剂=10000），基因=hic促进者$基因）</li>
		<li>将可信的 SN 与 Hi-C GRange 对象重叠。 olap = 查找重叠（信任范围，希克增强器） 信元 = 信任范围[查询Hits（olap）] mcols（信任）= cbind（mcols（信任），mcols（希克增强剂[主题Hits（olap）]）</li>
	</ol>
</li>
	<li>编译由位置映射和染色质相互作用曲线定义的 AD 候选基因。 • AD 的候选基因： ADgenes = 减少（联合，列表（信不自他基因，信不自生自新$基因，可信基因）） • 将恩森布尔基因 ID 转换为 HGNC 符号 负载（"基因Anno.rda"） ADhgnc = 基因Anno1_匹配（ADgenes，基因Anno1$ensembl_gene_id），"hgnc_symbol"] ADhgnc = ADhgnc[ADhgnc！=""= 保存（ADgenes、ADhgnc、文件="ADgenes.rda"） write.table（ADhgnc，文件="ADgenes.txt"，行.names_F，col.names_F，报价=F，sep="\t"）</li>
</ol>4. 发展表达轨迹
注： 对于每个步骤，在 RStudio 中的控制台窗口中键入相应的代码。
<ol><li>在 R 中设置。 库（重塑）;库（ggplot2）;库（基因组学）;图书馆（生物硕士） 图书馆（"WGCNA"） 选项（字符串AsFactors_F）</li>
	<li>处理表达式和元数据。 datExpr = read.csv（"expression_matrix.csv"，标题 = FALSE） datExpr = datExpr_，-1] 达特梅塔 = read.csv（"columns_metadata.csv"） datProbes = read.csv（"rows_metadata.csv"） datExpr = datExpr_datProbes$ensembl_gene_id！""，"" datProbes = datProbes_datProbes$ensembl_gene_id！""，""，* datExpr.cr= 折叠行（datExpr，行组 = datProbes$ensembl_gene_id，行 ID = 行名（datExpr）） datExpr = datExpr.cr$datET 折叠 gename = 数据.帧（datExpr.cr$group2row） 行名（datExpr） = gename$组<ol>
<li>指定发展阶段。 达特梅塔 $单位 = "产后" idx = grep（"pcw"，达特梅塔美元年龄） 达特梅塔$单位[idx] = "产前" idx = grep（"年"，达特梅塔$年龄） 达特梅塔$单位[idx] = "产后" 达特梅塔元单位 = 系数（达美元$单位，水平=c（"产前"，"产后"））</li>
		<li>选择皮质区域。 达特梅塔$区域 = "子CTX" r = c（"A1C"、"STC"、"ITC"、"TCx"、"OFC"、"DFC"、"VFC"、"MFC"、"M1C"、"S1C"、"IPC"、"M1C-S1C"、"PCx"、"V1C"、"Ocx"） datMeta$区域_datMeta$structure_acronym %in% r= = "CTX" datExpr = datExpr_，其中（达特梅塔$Region="CTX"）* 达特梅塔 = 达特梅塔[其中（达特梅塔$Region="CTX"），* 保存（datExpr、datMeta、文件="devExpr.rda"）</li>
	</ol>
</li>
	<li>提取AD风险基因的发育表达图谱。 负载（"ADgenes.rda"） exprdat = 应用（datExpr_匹配（ADgenes、行名（datExpr））），2，均值，na.rm_T） dat = 数据.帧（区域_datMeta$Region，单位_datMeta$单位，Expr_exprdat）</li>
	<li>比较AD风险基因的产前和产后表达水平。 pdf（文件="developmental_expression.pdf"） ggplot（dat，aes（x=单位，y=Expr，填充=单位，alpha=单位））=ylab（"标准化表达式"）=geom_boxplot（异常值.大小= NA）= ggtitle（"大脑表达式"） = xlab（""）= scale_alpha_manual（值=c（0.2，1）=theme_classic（）= 主题（图例.) dev.off（）</li>
</ol>5. 单元格类型的表达式配置文件
注： 对于每个步骤，在 RStudio 中的控制台窗口中键入相应的代码。
<ol><li>在 R 中设置。 选项（字符串AsFactors_F） 负载（"ADgenes.rda"） 负载（"基因Anno.rda"） 目标名称 = "AD" 目标基因 = ADhgnc 单元格exp = read.表（"DER-20_Single_cell_expression_processed_TPM_backup.tsv"，标题=T，填充=T） 单元格exp[1121，1] = 细胞exp[1120，1] 细胞外号 = 细胞exp_-1120，* 行名（单元格） = 单元格exp_，1] 细胞外号 = 细胞exp_，-1] datExpr = 刻度（cellexp，center=T，刻度=F） datExpr = datExpr_，789：ncol（datExpr）]</li>
	<li>提取AD风险基因的细胞表达谱。 exprdat = 应用（datExpr_匹配（目标基因、行名（datExpr））），2，均值，na.rm_T） dat = 数据.帧（组=目标名称，单元格名称（exprdat），Expr_exprdat） dat $celltype = 未列表（l 应用（strsplit（dat$cell，拆分=[.]），''''，1）） dat = dat_-grep（"Ex_在"，dat$细胞类型），* dat $细胞型 = gsub（"Dev"，"胎儿"，dat$细胞型） dat$细胞类型 = 因子（dat$细胞型，水平=c（"神经元"，"星形细胞"，"微胶"，"内皮"， "奥托根德罗茨"，"OPC"，"胎儿"）） pdf（文件="singlecell_expression_ADgenes.pdf"） ggplot（dat，aes（x=细胞类型，y=Expr，填充=细胞类型）））* ylab（"标准化表达式"） = xlab（"） = geom_violin（） = 主题（轴.text.x=element_text（角度 = 90，hjust =1））） = 主题（图例"无"） | ggtitle（粘贴0（"AD风险基因的细胞表达特征"）） dev.off（）</li>
</ol>6. AD风险基因的基因注释扩充分析
<ol><li>通过在终端中键入以下命令下载并配置 HOMER。 姆克迪尔·荷马 cd 主 wget http://homer.ucsd.edu/homer/configureHomer.pl perl ./configureHomer.pl -安装 perl ./configureHomer.pl -安装人-p perl ./configureHomer.pl -安装人-o</li>
	<li>通过在终端中键入以下命令来运行 HOMER。 导出 PATH=$PATH：*/工作/主站/箱 findMotifs.pl [/工作/ADgenes.txt 人类 /工作/</li>
	<li>通过在 RStudio 中的控制台窗口中键入以下代码来绘制丰富术语。 库（格普布） 选项（字符串AsFactors_F） pdf（"GO_enrichment.pdf"，宽度=15，高度=8） plot_barplot = 函数（名称、名称、颜色）* 输入 = read.delim（粘贴0（dbname，".txt"），标题=T） 输入 = 输入，c（-1，-10，-11）* 输入 = 唯一（输入） 输入$FDR = p.调整（exp（输入$logP）） input_sig = 输入 [输入 $FDR < 0.1，] input_sig$FDR = -log10（input_sig$FDR） input_sig = input_sig[订单（input_sig$FDR），] p = ggbarplot（input_sig，x = "术语"，y = "FDR"，填充 = 颜色，颜色 = "白色"， sort.val = "asc"，ylab = 表达式（-log=10]（斜体（FDR）），xlab = 粘贴0（名称，术语），旋转 = TRUE，标签 = 粘贴0（input_sig$Target.Genes.in.term，"/"input_sig$Genes.in.term），字体.label = 列表（颜色 = "白色"，大小 = 9），lab.vjust = 0.5，lab.5 p = p= geom_hline（yintercept = -log10（0.05），线型 = 2，颜色 ="浅灰色"） 返回（p） } p1 = plot_barplot（"biological_process"，"GO生物过程"，"#00AFBB"） p2 = plot_barplot（"kegg"，"KEGG"，"#E7B800"） p3 = plot_barplot（"反应"，"反应"，"#FC4E07"） gg排列（p1，p2，p3，标签 = c（"A"，"B"，"C"），ncol = 2，nrow = 2） dev.off（）</li>
</ol>

<ol>
	<li>Dekker, J., Misteli, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Long-Range+Chromatin+Interactions.">Long-Range Chromatin Interactions.</a> Cold Spring Harbor Perspectives in Biology. 7 (10), a019356 (2015).</li><li>Sanyal, A., Lajoie, B. R., Jain, G., Dekker, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+long-range+interaction+landscape+of+gene+promoters.">The long-range interaction landscape of gene promoters.</a> Nature. 489 (7414), 109-113 (2012).</li><li>Plank, J. L., Dean, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Enhancer+function:+mechanistic+and+genome-wide+insights+come+together.">Enhancer function: mechanistic and genome-wide insights come together.</a> Molecular Cell. 55 (1), 5-14 (2014).</li><li>Dekker, J., Marti-Renom, M. A., Mirny, L. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Exploring+the+three-dimensional+organization+of+genomes:+interpreting+chromatin+interaction+data.">Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data.</a> Nature Reviews Genetics. 14 (6), 390-403 (2013).</li><li>Martin, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Capture+Hi-C+reveals+novel+candidate+genes+and+complex+long-range+interactions+with+related+autoimmune+risk+loci.">Capture Hi-C reveals novel candidate genes and complex long-range interactions with related autoimmune risk loci.</a> Nature Communications. 6, 10069 (2015).</li><li>Won, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Chromosome+conformation+elucidates+regulatory+relationships+in+developing+human+brain.">Chromosome conformation elucidates regulatory relationships in developing human brain.</a> Nature. 538 (7626), 523-527 (2016).</li><li>J&#228;ger, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Capture+Hi-C+identifies+the+chromatin+interactome+of+colorectal+cancer+risk+loci.">Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci.</a> Nature Communications. 6, 6178 (2015).</li><li>Chen, J. A. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Joint+genome-wide+association+study+of+progressive+supranuclear+palsy+identifies+novel+susceptibility+loci+and+genetic+correlation+to+neurodegenerative+diseases.">Joint genome-wide association study of progressive supranuclear palsy identifies novel susceptibility loci and genetic correlation to neurodegenerative diseases.</a> Molecular Neurodegeneration. 13 (1), 41 (2018).</li><li>Wang, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comprehensive+functional+genomic+resource+and+integrative+model+for+the+adult+brain.">Comprehensive functional genomic resource and integrative model for the adult brain.</a> Science. 362 (6420), eaat8464 (2018).</li><li>Demontis, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Discovery+of+the+first+genome-wide+significant+risk+loci+for+attention+deficit/hyperactivity+disorder.">Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder.</a> Nature Genetics. 51 (1), 63-75 (2019).</li><li>Grove, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Identification+of+common+genetic+risk+variants+for+autism+spectrum+disorder.">Identification of common genetic risk variants for autism spectrum disorder.</a> Nature Genetics. 51 (3), 431-444 (2019).</li><li>Lee, P. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome+wide+meta-analysis+identifies+genomic+relationships,+novel+loci,+and+pleiotropic+mechanisms+across+eight+psychiatric+disorders.">Genome wide meta-analysis identifies genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders.</a> bioRxiv. , 528117 (2019).</li><li>Mah, W., Won, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+three-dimensional+landscape+of+the+genome+in+human+brain+tissue+unveils+regulatory+mechanisms+leading+to+schizophrenia+risk.">The three-dimensional landscape of the genome in human brain tissue unveils regulatory mechanisms leading to schizophrenia risk.</a> Schizophrenia Research. , (2019).</li><li>Jansen, I. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome-wide+meta-analysis+identifies+new+loci+and+functional+pathways+influencing+Alzheimer's+disease+risk.">Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer's disease risk.</a> Nature Genetics. 51 (3), 404-413 (2019).</li><li>Viola, K. L., Klein, W. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Amyloid+&#946;+oligomers+in+Alzheimer's+disease+pathogenesis,+treatment,+and+diagnosis.">Amyloid &#946; oligomers in Alzheimer's disease pathogenesis, treatment, and diagnosis.</a> Acta Neuropathologica. 129 (2), 183-206 (2015).</li><li>Mroczko, B., Groblewska, M., Litman-Zawadzka, A., Kornhuber, J., Lewczuk, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Amyloid+&#946;+oligomers+(A&#946;Os)+in+Alzheimer's+disease.">Amyloid &#946; oligomers (A&#946;Os) in Alzheimer's disease.</a> Journal of Neural Transmission. 125 (2), 177-191 (2018).</li><li>Heneka, M. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Neuroinflammation+in+Alzheimer's+disease.">Neuroinflammation in Alzheimer's disease.</a> Lancet Neurology. 14 (4), 388-405 (2015).</li><li>Minter, M. R., Taylor, J. M., Crack, P. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+contribution+of+neuroinflammation+to+amyloid+toxicity+in+Alzheimer's+disease.">The contribution of neuroinflammation to amyloid toxicity in Alzheimer's disease.</a> Journal of Neurochemistry. 136 (3), 457-474 (2016).</li><li>Hansen, D. V., Hanson, J. E., Sheng, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Microglia+in+Alzheimer's+disease.">Microglia in Alzheimer's disease.</a> The Journal of Cell Biology. 217 (2), 459-472 (2018).</li><li>Gjoneska, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Conserved+epigenomic+signals+in+mice+and+humans+reveal+immune+basis+of+Alzheimer's+disease.">Conserved epigenomic signals in mice and humans reveal immune basis of Alzheimer's disease.</a> Nature. 518 (7539), 365-369 (2015).</li><li>Benner, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=FINEMAP:+efficient+variable+selection+using+summary+data+from+genome-wide+association+studies.">FINEMAP: efficient variable selection using summary data from genome-wide association studies.</a> Bioinformatics. 32 (10), 1493-1501 (2016).</li><li>Hormozdiari, F., Kostem, E., Kang, E. Y., Pasaniuc, B., Eskin, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Identifying+causal+variants+at+loci+with+multiple+signals+of+association.">Identifying causal variants at loci with multiple signals of association.</a> Genetics. 198 (2), 497-508 (2014).</li><li>Tewhey, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Direct+Identification+of+Hundreds+of+Expression-Modulating+Variants+using+a+Multiplexed+Reporter+Assay.">Direct Identification of Hundreds of Expression-Modulating Variants using a Multiplexed Reporter Assay.</a> Cell. 165 (6), 1519-1529 (2016).</li><li>Arnold, C. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome-wide+quantitative+enhancer+activity+maps+identified+by+STARR-seq.">Genome-wide quantitative enhancer activity maps identified by STARR-seq.</a> Science. 339 (6123), 1074-1077 (2013).</li><li>Kent, W. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+human+genome+browser+at+UCSC.">The human genome browser at UCSC.</a> Genome Research. 12 (6), 996-1006 (2002).</li><li>Kang, H. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Spatio-temporal+transcriptome+of+the+human+brain.">Spatio-temporal transcriptome of the human brain.</a> Nature. 478 (7370), 483-489 (2011).</li></ol>

作者没有什么可透露的。

在这里，我们描述了一个分析框架，它可用于基于位置映射和染色质相互作用对GWS位点进行功能性说明。此过程涉及多个步骤（有关更多详细信息，请参阅此评论13）。首先，鉴于染色质相互作用的轮廓是高度细胞类型的特异性，Hi-C数据从适当的细胞/组织类型获得，最好地捕获疾病的基础生物学需要使用。鉴于AD是一种神经退行性疾病，我们使用成人大脑Hi-C数据9来对GWS位点进行分文。其次，由于联动不平衡 （LD），每个 GWS 位点通常具有多达数百个与特征相关的 SNPs，因此，通过使用精细映射算法21、22或使用高通量方法（如大规模并行报告器 （MPRA）23或自转录活动区域排序）计算预测因果关系，获得假定因果（"可信"）SNPs 非常重要。STARR-seq）24.对于这里描述的工作，我们使用了Jansen等人14中报道的可信的SNPs。第三，启动子和外音SNPs是基于位置映射进行加法。我们使用一个简单的位置映射策略，其中 SNPs 在与启动器（定义为转录起始位上游的 2 kb）或外源子重叠时映射到基因。但是，可以通过评估外音 SNPs 的功能后果来进一步阐述此方法，例如 SNP 是否导致无意义的中介衰变、误解变异或无意义变异。第四，来自相应组织/细胞类型的染色质相互作用曲线可用于根据物理接近性为其假定目标基因分配SNPs。我们使用锚定到启动器的交互配置文件，但我们可以通过将增强剂活动（由组蛋白 H3 K27 乙酰化或染色质可访问性引导）或外音相互作用考虑进一步优化或扩展交互配置文件。在这个过程中的一个重要考虑因素是使用一致的人类基因组构建。例如，如果汇总统计的基因组位置不基于 hg19（即 hg18 或 hg38），则应获得参考基因组的适当版本，或者需要使用提升25将汇总统计数据转换为 hg19。
我们应用这个框架来识别AD GWAS的假定目标基因，将284个SNPs分配给112个AD风险基因。使用发育表达谱26和细胞类型特定表达谱9，我们随后证明，该基因集与已知的AD病理学一致，揭示了细胞类型（微胶质）、生物功能（免疫反应和淀粉样蛋白β），以及年龄风险升高。
虽然我们提出了一个框架，其中标定了AD的潜在目标基因及其基础生物学，但值得注意的是，基于Hi-C的注释可以展开，以注释任何非编码变异。随着更多的全基因组测序数据可用，我们对非编码罕见变异的理解也越来越多，Hi-C将为解释疾病相关基因变异提供关键资源。因此，从多种组织和细胞类型中获得的Hi-C资源简编对于促进广泛应用这一框架以收集对各种人类特征和疾病的生物见解至关重要。

全基因组关联研究（GWAS）在揭示一系列人类特征和疾病的遗传基础方面发挥了关键作用。这种大规模的基因分型已经发现了数千种与表型相关的基因组变异，范围从身高到精神分裂症风险。然而，尽管GWAS在识别疾病和特征相关位点方面取得了巨大成功，但机械学对这些变异如何促成表型的理解一直具有挑战性，因为大多数表型相关变异都存在于非编码中人类基因组的一小部分。由于这些变异通常与预测的调控元素重叠，它们可能会改变附近基因的转录控制。然而，非编码位点可以影响基因的转录，线性距离超过一兆碱基，使受每个变异影响的基因难以识别。三维（3D）染色质结构在调解远距离调控位点和基因启动子之间的连接方面起着重要作用，可用于识别受表型相关单核苷酸多态性（SNPs）影响的基因。
基因调控是由一个复杂的过程，其中涉及增强剂激活和染色质循环形成，物理连接增强剂的基因启动子，转录机制可以定向1，2，3。由于染色质环通常跨越几百千基（kb），需要3D染色质结构的详细地图来破译基因调控机制。已发明多种染色质构象捕获技术，以识别3D染色质架构4。在这些技术中，Hi-C 提供了最全面的架构，因为它捕获了全基因组的 3D 染色质相互作用配置文件。Hi-C数据集已经迅速被调整，以解释非编码全基因组显著（GWS）位点5，6，7，8，9，10，11，12，13，因为它可以将非编码变异链接到其假定的目标基因基于染色质相互作用配置文件。
在本文中，我们概述了使用染色质相互作用谱计算预测GWAS风险变异的假定目标基因的协议。我们应用这个协议来映射ADGWS位点14到他们的目标基因使用Hi-C数据集在成人人脑9。由此产生的AD风险基因由其他功能基因组数据集为特征，包括单细胞转录组和发育表达谱。

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>10 kb resolution Hi-C interaction profiles in the adult brain from psychencode</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="http://adult.psychencode.org/Datasets/Integrative/Promoter-anchored_chromatin_loops.bed">http://adult.psychencode.org/</a></td>
		</tr>
		<tr>
			<td>Developmental expression datasets</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="http://www.brainspan.org/api/v2/well_known_file_download/267666527">http://www.brainspan.org/</a></td>
		</tr>
		<tr>
			<td>Fine-mapped credible SNPs for AD (Supplementary Table 8 from Jansen et al.14)</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="https://static-content.springer.com/esm/art%3A10.1038%2Fs41588-018-0311-9/MediaObjects/41588_2018_311_MOESM3_ESM.xlsx">https://static-content.springer.com/</a></td>
		</tr>
		<tr>
			<td>HOMER</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="http://homer.ucsd.edu/homer/configureHomer.pl">http://homer.ucsd.edu/</a></td>
		</tr>
		<tr>
			<td>R (version 3.5.0)</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="https://www.r-project.org/">https://www.r-project.org/</a></td>
		</tr>
		<tr>
			<td>RStudio Desktop</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="https://www.rstudio.com/products/rstudio/download/">https://www.rstudio.com/</a></td>
		</tr>
		<tr>
			<td>Single cell expression datasets</td>
			<td></td>
			<td></td>
			<td><a target="_blank" href="http://adult.psychencode.org/Datasets/Derived/SC_Decomp/DER-20_Single_cell_expression_processed_TPM_backup.tsv">http://adult.psychencode.org/</a></td>
		</tr>
	</tbody>
</table>

mapping alzheimer's disease variants to their target genes using computational analysis of chromatin configuration

全基因组关联研究（GWAS）已经成功地确定了数百个与人类特征和疾病相关的基因组位点。然而，由于大多数全基因组显著（GWS）位点都属于非编码基因组，许多功能影响仍不得而知。由Hi-C或其衍生物识别的三维染色质相互作用可以通过将非编码变异与其可操作基因联系起来，为这些位点进行说明提供有用的工具。在这里，我们概述了一个协议，以映射GWAS非编码变体，以映射到他们的假定基因使用阿尔茨海默氏病（AD）GWAS和Hi-C数据集从人类成人脑组织。通过应用精细映射算法来识别假定因果单核苷酸多态性 （SNPs）。然后，使用基于 Hi-C 的增强剂-促进者相互作用，将 SNPs 映射到其假定目标基因。生成的基因集代表 AD 风险基因，因为它们可能受 AD 风险变异的调控。为了进一步深入了解AD背后的分子机制，我们使用发育性大脑表达数据和大脑单细胞表达图谱来描述AD风险基因。该协议可以扩展到任何GWAS和Hi-C数据集，以识别各种人类特征和疾病背后的假定靶基因和分子机制。

此处描述的过程适用于原始研究14定义的一组800个可信的SNPs。位置图显示，103个SNPs与启动子重叠（43个唯一基因），42个SNPs与外显子重叠（27个唯一基因）。位置映射后，84% （669） SNPs 仍未被指示。在成人大脑中使用Hi-C数据集，我们能够根据物理接近性将另外208个SNPs与64个基因联系起来。我们总共将 284 个 AD 可信 SNPs 映射到 112 个 AD 风险基因（图 1A）。AD风险基因与淀粉样蛋白前体蛋白、淀粉样蛋白-β形成和免疫反应有关，反映了AD15、16、17、18的已知生物学（图1B-D）。AD风险基因的发育表达特征显示产后显著富集，表明与年龄相关的AD风险升高（图2A）。最后，AD风险基因在微胶质、脑原免疫细胞中高度表达（图2B）。这与复发性发现一致，即AD具有很强的免疫基础，微胶质是AD发病机制14、19、20的中心作用者。
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/60428/60428fig1b.jpg"> 图1：定义ADGWS位点的假定目标基因。（A） 从前 29 个 AD 位点派生的可信 SNPs 被分类为启动子SNPs、外音SNPs和未指定非编码的SNPs。发起人和外音SNPs通过位置映射直接分配给其目标基因，而成人大脑中的色母相互作用曲线则另外用于根据物理相互作用绘制SNPs。（B-D）AD风险基因中的GO（ B）、KEGG（C）和反应素（D）项的富集使用协议第6节所述。x 轴表示已更正的错误发现率 （FDR） -log10（P 值）。绘制了 FDR < 0.1 的丰富术语。灰色垂直线表示 FDR = 0.05。APP淀粉样蛋白前体蛋白。分子，每个术语中表示的AD风险基因的数量;分母，每个术语中的基因数量。<a href="https://www.jove.com/files/ftp_upload/60428/60428fig1alarge.jpg" target="_blank">请点击此处查看此图的较大版本。</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/60428/60428fig2.jpg"> 图2：AD风险基因的表征。（A） 与产前皮层相比， AD 风险基因在产后皮层中高度表达.（B） 小提琴图描绘了不同细胞类型中不同细胞类型的基因表达值（规范化表达）的分布。这些结果表明，AD风险基因在微胶质中高度表达，与以前的研究14一致。<a href="https://www.jove.com/files/ftp_upload/60428/60428fig2large.jpg" target="_blank">请点击此处查看此图的较大版本。</a>
补充文件1。<a href="https://www.jove.com/files/ftp_upload/60428/Gencode19_exon.bed" target="_blank">请点击此处查看此文件（右键单击下载）。</a>
补充文件2。<a href="https://www.jove.com/files/ftp_upload/60428/Gencode19_promoter.bed" target="_blank">请点击此处查看此文件（右键单击下载）。</a>
补充文件3.<a href="https://www.jove.com/files/ftp_upload/60428/geneAnno.rda" target="_blank">请点击此处查看此文件（右键单击下载）。</a>

Watch this Scientific Journal Video about Mapping Alzheimer's Disease Variants to Their Target Genes Using Computational Analysis of Chromatin Configuration at JoVE.com

Mapping Alzheimer's Disease Variants to Their Target Genes Using Computational Analysis of Chromatin Configuration

我们提出了一个协议，用于识别使用三维染色质相互作用的全基因组关联研究（GWAS）识别的非编码变异的功能影响。

利用染色质配置的计算分析，将阿尔茨海默氏病变异体映射到其目标基因

Genome-wide association studies (GWAS) have successfully identified hundreds of genomic loci that are associated with human traits and disease. However, ...

虽然GWAS已经成功地识别了与人类特征和疾病相关的基因组区域，但这些风险变异的生物影响尚不清楚。在这里，我们概述了一个协议，使用染色质相互作用配置文件计算预测GWAS风险变异的预计目标基因。通常，识别风险基因是了解疾病机制的第一步，并允许正常的治疗方法。我们希望，这项工作的结果最终能够导致诊断和治疗阿尔茨海默氏症的最终策略。这项技术的主要优点是，通过使用3D色度接触频率，我们可以识别受阿尔茨海默病影响的基因风险方差，即使他们是数千甚至数百万的碱基对。尝试此协议时，熟悉 R 或 X 对系统至关重要，因为用户需要对系统执行整个协议。若要执行此计算协议，请参阅文本手稿或屏幕上的代码。首先，在 R 中设置，为可信的单核苷酸多态性或 SNPS 生成 G 范围对象。对于位置映射，在 R 中设置，然后加载启动子和外向区域并生成 G 范围对象。将可信的 SNPS 与外向区域和启动区域重叠。要使用 Chromaton 相互作用将 SNPS 链接到其生成的目标基因，请加载 Hi C 数据集并生成 G 范围对象。将可信的 SNPS 与 Hi C G 范围对象重叠。并编译AD候选基因，由位置映射和色度体相互作用配置文件定义。接下来，探索发展轨迹。在 R 中设置并处理表达式元数据。指定发展阶段并选择皮质区域。提取AD风险基因的发育表达特征，比较产前与产后表达水平。通过 R 中设置并提取 AD 风险的细胞表达式配置文件来调查细胞类型表达式配置文件。最后，对AD风险基因进行基因注释富集分析。下载并配置荷马。然后运行荷马，并绘制与R工作室丰富的术语。使用此过程对一组 800 个可信的 SNPs 进行了调查。位置映射显示，103 个 SNP 与启动子重叠，42 个 SNP 与 Exons 重叠，而 84% 的 SNP 未通知。在成人大脑中使用高C数据集，另外208个SNPs与64个基于物理接近的基因相关。总共有284个AD可信的SNPs被映射到112个AD风险基因。AD风险基因与淀粉样蛋白、淀粉样蛋白形成和免疫反应有关，这反映了该疾病的已知生物学。AD风险基因的发育表达特征显示，产后富集显著，表明与年龄相关的疾病风险升高。最后，这些基因在微胶质中高度表达，这是大脑中的主要免疫细胞，支持了AD具有很强的免疫基础的经常性发现。在这里，我们使用来自脑组织的Hi-C数据来分析阿尔茨海默病风险差异的生物学影响。然而，将此方法应用于另一个 GWAS 研究，相关组织中新的 Hi-C 数据水平至关重要。这些结果可以通过基于脆性的技术、增强剂报告器检测或与其他功能基因组数据集（如殷拓）相交来进一步研究和验证。在这里，我们识别了几十个阿尔茨海默病风险基因，我们期望这些基因的识别可以帮助我们了解它们以前在阿尔茨海默病中未知的角色。

Research

JoVE Journal

Genetics

Mapping RNA-RNA Interactions Globally Using Biotinylated Psoralen

使用生物素化补骨脂素全球定位RNA-RNA相互作用

Knowing how RNAs interact with themselves and with others is key to understanding RNA based gene regulation in the cell. While examples of RNA-RNA ...

科研

遗传学

Sequential Salt Extractions for the Analysis of Bulk Chromatin Binding Properties of Chromatin Modifying Complexes

序贯盐萃取法分析染色质修饰配合物的体染色质结合特性

Elucidation of the binding properties of chromatin-targeting proteins can be very challenging due to the complex nature of chromatin and the heterogeneous ...

Analysis of the c-KIT Ligand Promoter Using Chromatin Immunoprecipitation

使用染色质免疫沉淀分析c-KIT配体启动子

Multiple cellular processes, including DNA replication and repair, DNA recombination, and gene expression, require interactions between proteins and DNA. ...

Retroviral Scanning: Mapping MLV Integration Sites to Define Cell-specific Regulatory Regions

逆转录病毒扫描：映射MLV集成站点以定义细胞特异性监管区域

Moloney murine leukemia (MLV) virus-based retroviral vectors integrate predominantly in acetylated enhancers and promoters. For this reason, mLV ...

Mapping Genome-wide Accessible Chromatin in Primary Human T Lymphocytes by ATAC-Seq

ATAC 基因组在原发性人类 T 淋巴细胞中的染色体全易接近染色质

Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a method used for the identification of open (accessible) regions ...

Immunostaining for DNA Modifications: Computational Analysis of Confocal Images

DNA 修饰的免疫: 共焦图像的计算分析

For several decades, 5-methylcytosine (5mC) has been thought to be the only DNA modification with a functional significance in metazoans. The discovery of ...

Using Sniper-Cas9 to Minimize Off-target Effects of CRISPR-Cas9 Without the Loss of On-target Activity Via Directed Evolution

利用 niper-cas9 在不损失目标活动的情况下通过定向进化, 最大限度地减少 crispr-cas9 的非目标效应

The development of clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (Cas9) into therapeutic modalities requires the ...

An In Vitro Protocol for Evaluating MicroRNA Levels, Functions, and Associated Target Genes in Tumor Cells

用于评估肿瘤细胞中的微RNA水平、功能和相关靶基因的体外协议

MicroRNAs (miRNAs) are small regulatory RNAs which are recognized to modulate numerous intracellular signaling pathways in several diseases including ...

In Vivo Functional Study of Disease-associated Rare Human Variants Using Drosophila

In Vivo Functional Study of Disease-associated Rare Human Variants Using Drosophila

在Vivo功能研究的疾病相关的稀有人类变异使用果蝇

Advances in sequencing technology have made whole-genome and whole-exome datasets more accessible for both clinical diagnosis and cutting-edge human ...