Name: Mapping Alzheimers Disease Variants to Their Target Genes Using Computational Analysis of Chromatin Configuration (Video) | JoVE Journal
Uploaded: 2020-01-09T13:00:00.000+00:00
Duration: 4 min 41 s
Description: Discover how to computationally predict target genes of Alzheimer's disease GWAS risk variants using Hi-C chromatin interaction profiles. This protocol maps SNPs to genes and explores developmental expression patterns.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

This work was supported by the NIH grant R00MH113823 (to H.W.) and R35GM128645 (to D.H.P.), NARSAD Young Investigator Award (to H.W.), and SPARK grant from the Simons Foundation Autism Research Initiative (SFARI, to N.M. and H.W.).

1. Workstation Setup
<ol>
	<li>Install R (version 3.5.0) and RStudio Desktop. Open RStudio.</li>
	<li>Install the following libraries in R by typing the following code into the console window in RStudio. 
	if (!&#34;BiocManager&#34; %in% rownames(installed.packages())) 
	install.packages(&#34;BiocManager&#34;, repos=&#34;https://cran.r-project.org&#34;) 
	BiocManager::install(&#34;GenomicRanges&#34;) 
	BiocManager::install(&#34;biomaRt&#34;) 
	BiocManager::install(&#34;WGCNA&#34;) 
	install.packages(&#34;reshape&#34;) 
	install.packages(&#34;ggplot2&#34;) 
	install.packages(&#34;corrplot&#34;) 
	install.packages(&#34;gProfileR&#34;) 
	install.packages(&#34;tidyverse&#34;) 
	install.packages(&#34;ggpubr&#34;)</li>
	<li>Download files. 
	NOTE: In this protocol, all files are required to be downloaded to ~/work directory.
	<ol>
		<li>Download the following files by clicking the links provided in Table of Materials.
		<ol>
			<li>Download fine-mapped credible SNPs for AD (Supplementary Table 8 from Jansen et al.14). 
			NOTE: Before analysis, open sheet eight in 41588_2018_311_MOESM3_ESM.xlsx, remove the first three rows and save the sheet as Supplementary_Table_8_Jansen.txt with tab separated format.</li>
			<li>Download 10 kb resolution Hi-C interaction profiles in the adult brain from psychencode (described as Promoter-anchored_chromatin_loops.bed below). 
			NOTE: This file has the following format: chromosome, TSS_start, TSS_end, Enhancer_start, and Enhancer_end. In case other Hi-C datasets are used, this protocol requires Hi-C datasets processed at high resolution (5&#8722;20 kb).</li>
			<li>Download single cell expression datasets from the PsychENCODE. 
			NOTE: These are from neurotypical control samples.</li>
			<li>Download developmental expression datasets from the BrainSpan (described as devExpr.rda below). 
			NOTE: 267666527 is a zipped file, so unzip the 267666527 to extract &#34;columns_metadata.csv&#34;, &#34;expression_matrix.csv&#34;, and &#34;rows_metadata.csv&#34; to generate devExpr.rda (see section 3).</li>
		</ol>
		</li>
		<li>Download exonic coordinates (see Supplementary Files, described as Gencode19_exon.bed and Gencode19_promoter.bed below) from Gencode version 19. 
		NOTE: Promoters are defined as 2 kb upstream of transcription start site (TSS). These files have the following format: chromosome, start, end, and gene.</li>
		<li>Download gene annotation file (see Supplementary Files, described as geneAnno.rda below) from biomart. 
		NOTE: This file can be used to match genes based on Ensembl gene IDs and HUGO Gene Nomenclature Committee (HGNC) symbol.</li>
	</ol>
	</li>
</ol>
2. Generation of a GRanges Object for Credible SNPs
<ol>
	<li>Set up in R by typing the following code into the console window in RStudio. 
	library(GenomicRanges) 
	options(stringsAsFactors = F) 
	setwd(&#34;~/work&#34;) # This is the path to the working directory. 
	credSNP = read.delim(&#34;Supplementary_Table_8_Jansen.txt&#34;, header=T) 
	credSNP = credSNP[credSNP$Credible.Causal==&#34;Yes&#34;,]</li>
	<li>Make a GRanges object by typing the following code into the console window in RStudio. 
	credranges = GRanges(credSNP$Chr, IRanges(credSNP$bp, credSNP$bp), rsid=credSNP$SNP, P=credSNP$P) 
	save(credranges, file=&#34;AD_credibleSNP.rda&#34;)</li>
</ol>
3. Positional Mapping
NOTE: For each step, type the corresponding code into the console window in RStudio.
<ol>
	<li>Set up in R. 
	options(stringsAsFactors=F) 
	library(GenomicRanges) 
	load(&#34;AD_credibleSNP.rda&#34;) # (see 2)</li>
	<li>Positional mapping of promoter/exonic SNPs to genes
	<ol>
		<li>Load promoter and exonic region and generate a GRange object. 
		exon = read.table(&#34;Gencode19_exon.bed&#34;) 
		exonranges = GRanges(exon[,1],IRanges(exon[,2],exon[,3]),gene=exon[,4]) 
		promoter = read.table(&#34;Gencode19_promoter.bed&#34;) 
		promoterranges = GRanges(promoter[,1], IRanges(promoter[,2], promoter[,3]), gene=promoter[,4])</li>
		<li>Overlap credible SNPs with exonic regions. 
		olap = findOverlaps(credranges, exonranges) 
		credexon = credranges[queryHits(olap)] 
		mcols(credexon) = cbind(mcols(credexon), mcols(exonranges[subjectHits(olap)]))</li>
		<li>Overlap credible SNPs with promoter regions. 
		olap = findOverlaps(credranges, promoterranges) 
		credpromoter = credranges[queryHits(olap)] 
		mcols(credpromoter) = cbind(mcols(credpromoter), mcols(promoterranges[subjectHits(olap)]))</li>
	</ol>
	</li>
	<li>Link SNPs to their putative target genes using chromatin interactions.
	<ol>
		<li>Load Hi-C dataset and generate a GRange object. 
		hic = read.table(&#34;Promoter-anchored_chromatin_loops.bed &#34;, skip=1) 
		colnames(hic) = c(&#34;chr&#34;, &#34;TSS_start&#34;, &#34;TSS_end&#34;, &#34;Enhancer_start&#34;, &#34;Enhancer_end&#34;) 
		hicranges = GRanges(hic$chr, IRanges(hic$TSS_start, hic$TSS_end), enhancer=hic$Enhancer_start) 
		olap = findOverlaps(hicranges, promoterranges) 
		hicpromoter = hicranges[queryHits(olap)] 
		mcols(hicpromoter) = cbind(mcols(hicpromoter), mcols(promoterranges[subjectHits(olap)])) 
		hicenhancer = GRanges(seqnames(hicpromoter), IRanges(hicpromoter$enhancer, hicpromoter$enhancer+10000), gene=hicpromoter$gene)</li>
		<li>Overlap credible SNPs with Hi-C GRange object. 
		olap = findOverlaps(credranges, hicenhancer) 
		credhic = credranges[queryHits(olap)] 
		mcols(credhic) = cbind(mcols(credhic), mcols(hicenhancer[subjectHits(olap)]))</li>
	</ol>
	</li>
	<li>Compile AD candidate genes defined by positional mapping and chromatin interaction profiles. 
	### The resulting candidate genes for AD: 
	ADgenes = Reduce(union, list(credhic$gene, credexon$gene, credpromoter$gene)) 
	### to convert Ensembl Gene ID to HGNC symbol 
	load(&#34;geneAnno.rda&#34;) 
	ADhgnc = geneAnno1[match(ADgenes, geneAnno1$ensembl_gene_id), &#34;hgnc_symbol&#34;] 
	ADhgnc = ADhgnc[ADhgnc!=&#34;&#34;] 
	save(ADgenes, ADhgnc, file=&#34;ADgenes.rda&#34;) 
	write.table(ADhgnc, file=&#34;ADgenes.txt&#34;, row.names=F, col.names=F, quote=F, sep=&#34;\t&#34;)</li>
</ol>
4. Developmental Expression Trajectories
NOTE: For each step, type the corresponding code into the console window in RStudio.
<ol>
	<li>Set up in R. 
	library(reshape); library(ggplot2); library(GenomicRanges); library(biomaRt) 
	library(&#34;WGCNA&#34;) 
	options(stringsAsFactors=F)</li>
	<li>Process expression and meta data. 
	datExpr = read.csv(&#34;expression_matrix.csv&#34;, header = FALSE) 
	datExpr = datExpr[,-1] 
	datMeta = read.csv(&#34;columns_metadata.csv&#34;) 
	datProbes = read.csv(&#34;rows_metadata.csv&#34;) 
	datExpr = datExpr[datProbes$ensembl_gene_id!=&#34;&#34;,] 
	datProbes = datProbes[datProbes$ensembl_gene_id!=&#34;&#34;,] 
	datExpr.cr= collapseRows(datExpr, rowGroup = datProbes$ensembl_gene_id, rowID= rownames(datExpr)) 
	datExpr = datExpr.cr$datETcollapsed 
	gename = data.frame(datExpr.cr$group2row) 
	rownames(datExpr) = gename$group
	<ol>
		<li>Specify developmental stages. 
		datMeta$Unit = &#34;Postnatal&#34; 
		idx = grep(&#34;pcw&#34;, datMeta$age) 
		datMeta$Unit[idx] = &#34;Prenatal&#34; 
		idx = grep(&#34;yrs&#34;, datMeta$age) 
		datMeta$Unit[idx] = &#34;Postnatal&#34; 
		datMeta$Unit = factor(datMeta$Unit, levels=c(&#34;Prenatal&#34;, &#34;Postnatal&#34;))</li>
		<li>Select cortical regions. 
		datMeta$Region = &#34;SubCTX&#34; 
		r = c(&#34;A1C&#34;, &#34;STC&#34;, &#34;ITC&#34;, &#34;TCx&#34;, &#34;OFC&#34;, &#34;DFC&#34;, &#34;VFC&#34;, &#34;MFC&#34;, &#34;M1C&#34;, &#34;S1C&#34;, &#34;IPC&#34;, &#34;M1C-S1C&#34;, &#34;PCx&#34;, &#34;V1C&#34;, &#34;Ocx&#34;) 
		datMeta$Region[datMeta$structure_acronym %in% r] = &#34;CTX&#34; 
		datExpr = datExpr[,which(datMeta$Region==&#34;CTX&#34;)] 
		datMeta = datMeta[which(datMeta$Region==&#34;CTX&#34;),] 
		save(datExpr, datMeta, file=&#34;devExpr.rda&#34;)</li>
	</ol>
	</li>
	<li>Extract developmental expression profiles of AD risk genes. 
	load(&#34;ADgenes.rda&#34;) 
	exprdat = apply(datExpr[match(ADgenes, rownames(datExpr)),],2,mean,na.rm=T) 
	dat = data.frame(Region=datMeta$Region, Unit=datMeta$Unit, Expr=exprdat)</li>
	<li>Compare prenatal versus postnatal expression levels of AD risk genes. 
	pdf(file=&#34;developmental_expression.pdf&#34;) 
	ggplot(dat,aes(x=Unit, y=Expr, fill=Unit, alpha=Unit)) + ylab(&#34;Normalized expression&#34;) + geom_boxplot(outlier.size = NA) + ggtitle(&#34;Brain Expression&#34;) + xlab(&#34;&#34;) + scale_alpha_manual(values=c(0.2, 1)) + theme_classic() + theme(legend.position=&#34;na&#34;) 
	dev.off()</li>
</ol>
5. Cell-type Expression Profiles
NOTE: For each step, type the corresponding code into the console window in RStudio.
<ol>
	<li>Set up in R. 
	options(stringsAsFactors=F) 
	load(&#34;ADgenes.rda&#34;) 
	load(&#34;geneAnno.rda&#34;) 
	targetname = &#34;AD&#34; 
	targetgene = ADhgnc 
	cellexp = read.table(&#34;DER-20_Single_cell_expression_processed_TPM_backup.tsv&#34;,header=T,fill=T) 
	cellexp[1121,1] = cellexp[1120,1] 
	cellexp = cellexp[-1120,] 
	rownames(cellexp) = cellexp[,1] 
	cellexp = cellexp[,-1] 
	datExpr = scale(cellexp,center=T, scale=F) 
	datExpr = datExpr[,789:ncol(datExpr)]</li>
	<li>Extract cellular expression profiles of AD risk genes. 
	exprdat = apply(datExpr[match(targetgene, rownames(datExpr)),],2,mean,na.rm=T) 
	dat = data.frame(Group=targetname, cell=names(exprdat), Expr=exprdat) 
	dat$celltype = unlist(lapply(strsplit(dat$cell, split=&#34;[.]&#34;),&#39;[[&#39;,1)) 
	dat = dat[-grep(&#34;Ex|In&#34;,dat$celltype),] 
	dat$celltype = gsub(&#34;Dev&#34;,&#34;Fetal&#34;,dat$celltype) 
	dat$celltype = factor(dat$celltype, levels=c(&#34;Neurons&#34;,&#34;Astrocytes&#34;,&#34;Microglia&#34;,&#34;Endothelial&#34;, 
	Oligodendrocytes&#34;,&#34;OPC&#34;,&#34;Fetal&#34;)) 
	pdf(file=&#34;singlecell_expression_ADgenes.pdf&#34;) 
	ggplot(dat,aes(x=celltype, y=Expr, fill=celltype)) + 
	ylab(&#34;Normalized expression&#34;) + xlab(&#34;&#34;) + geom_violin() + theme(axis.text.x=element_text(angle = 90, hjust=1)) + theme(legend.position=&#34;none&#34;) + 
	ggtitle(paste0(&#34;Cellular expression profiles of AD risk genes&#34;)) 
	dev.off()</li>
</ol>
6. Gene Annotation Enrichment Analysis of AD Risk Genes
<ol>
	<li>Download and configure HOMER by typing the commands below in terminal. 
	mkdir homer 
	cd homer 
	wget http://homer.ucsd.edu/homer/configureHomer.pl 
	perl ./configureHomer.pl -install 
	perl ./configureHomer.pl -install human-p 
	perl ./configureHomer.pl -install human-o</li>
	<li>Run HOMER by typing the commands below in terminal. 
	export PATH=$PATH:~/work/homer/bin 
	findMotifs.pl ~/work/ADgenes.txt human ~/work/</li>
	<li>Plot the enriched terms by typing the following code into the console window in RStudio. 
	library(ggpubr) 
	options(stringsAsFactors=F) 
	pdf(&#34;GO_enrichment.pdf&#34;,width=15,height=8) 
	plot_barplot = function(dbname,name,color){ 
	input = read.delim(paste0(dbname,&#34;.txt&#34;),header=T) 
	input = input[,c(-1,-10,-11)] 
	input = unique(input) 
	input$FDR = p.adjust(exp(input$logP)) 
	input_sig = input[input$FDR &#60; 0.1,] 
	input_sig$FDR = -log10(input_sig$FDR) 
	input_sig = input_sig[order(input_sig$FDR),] 
	p = ggbarplot(input_sig, x = &#34;Term&#34;, y = &#34;FDR&#34;, fill = color, color = &#34;white&#34;, sort.val = &#34;asc&#34;, ylab = expression(-log[10](italic(FDR))), xlab = paste0(name,&#34; Terms&#34;), rotate = TRUE, label = paste0(input_sig$Target.Genes.in.Term,&#34;/&#34;,input_sig$Genes.in.Term), font.label = list(color = &#34;white&#34;, size = 9), lab.vjust = 0.5, lab.hjust = 1) 
	p = p+geom_hline(yintercept = -log10(0.05), linetype = 2, color = &#34;lightgray&#34;) 
	return(p) 
	} 
	p1 = plot_barplot(&#34;biological_process&#34;,&#34;GO Biological Process&#34;,&#34;#00AFBB&#34;) 
	p2 = plot_barplot(&#34;kegg&#34;,&#34;KEGG&#34;,&#34;#E7B800&#34;) 
	p3 = plot_barplot(&#34;reactome&#34;,&#34;Reactome&#34;,&#34;#FC4E07&#34;) 
	ggarrange(p1, p2, p3, labels = c(&#34;A&#34;, &#34;B&#34;, &#34;C&#34;), ncol = 2, nrow = 2) 
	dev.off()</li>
</ol>

<ol><li>Dekker, J., Misteli, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Long-Range+Chromatin+Interactions%5BTitle%5D)+AND+%22Cold+Spring+Harbor+Perspectives+in+Biology%22%5BJournal%5D)">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/pubmed?term=((The+long-range+interaction+landscape+of+gene+promoters%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">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/pubmed?term=((Enhancer+function%3A+mechanistic+and+genome-wide+insights+come+together%5BTitle%5D)+AND+%22Molecular+Cell%22%5BJournal%5D)">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/pubmed?term=((Exploring+the+three-dimensional+organization+of+genomes%3A+interpreting+chromatin+interaction+data%5BTitle%5D)+AND+%22Nature+Reviews+Genetics%22%5BJournal%5D)">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/pubmed?term=((Capture+Hi-C+reveals+novel+candidate+genes+and+complex+long-range+interactions+with+related+autoimmune+risk+loci%5BTitle%5D)+AND+%22Nature+Communications%22%5BJournal%5D)">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/pubmed?term=((Chromosome+conformation+elucidates+regulatory+relationships+in+developing+human+brain%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Chromosome conformation elucidates regulatory relationships in developing human brain.</a> Nature. 538 (7626), 523-527 (2016).</li>
<li>Jäger, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Capture+Hi-C+identifies+the+chromatin+interactome+of+colorectal+cancer+risk+loci%5BTitle%5D)+AND+%22Nature+Communications%22%5BJournal%5D)">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/pubmed?term=((Joint+genome-wide+association+study+of+progressive+supranuclear+palsy+identifies+novel+susceptibility+loci+and+genetic+correlation+to+neurodegenerative+diseases%5BTitle%5D)+AND+%22Molecular+Neurodegeneration%22%5BJournal%5D)">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/pubmed?term=((Comprehensive+functional+genomic+resource+and+integrative+model+for+the+adult+brain%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">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/pubmed?term=((Discovery+of+the+first+genome-wide+significant+risk+loci+for+attention+deficit%2Fhyperactivity+disorder%5BTitle%5D)+AND+%22Nature+Genetics%22%5BJournal%5D)">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/pubmed?term=((Identification+of+common+genetic+risk+variants+for+autism+spectrum+disorder%5BTitle%5D)+AND+%22Nature+Genetics%22%5BJournal%5D)">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/pubmed?term=((Genome+wide+meta-analysis+identifies+genomic+relationships%2C+novel+loci%2C+and+pleiotropic+mechanisms+across+eight+psychiatric+disorders%5BTitle%5D)+AND+%22bioRxiv%22%5BJournal%5D)">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/pubmed?term=((The+three-dimensional+landscape+of+the+genome+in+human+brain+tissue+unveils+regulatory+mechanisms+leading+to+schizophrenia+risk%5BTitle%5D)+AND+%22Schizophrenia+Research%22%5BJournal%5D)">The three-dimensional landscape of the genome in human brain tissue unveils regulatory mechanisms leading to schizophrenia risk.</a> Schizophrenia Research. , In press (2019).</li>
<li>Jansen, I. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome-wide+meta-analysis+identifies+new+loci+and+functional+pathways+influencing+Alzheimer%27s+disease+risk%5BTitle%5D)+AND+%22Nature+Genetics%22%5BJournal%5D)">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/pubmed?term=((Amyloid+%CE%B2+oligomers+in+Alzheimer%27s+disease+pathogenesis%2C+treatment%2C+and+diagnosis%5BTitle%5D)+AND+%22Acta+Neuropathologica%22%5BJournal%5D)">Amyloid β 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/pubmed?term=((Amyloid+%CE%B2+oligomers+%28A%CE%B2Os%29+in+Alzheimer%27s+disease%5BTitle%5D)+AND+%22Journal+of+Neural+Transmission%22%5BJournal%5D)">Amyloid β oligomers (Aβ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/pubmed?term=((Neuroinflammation+in+Alzheimer%27s+disease%5BTitle%5D)+AND+%22Lancet+Neurology%22%5BJournal%5D)">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/pubmed?term=((The+contribution+of+neuroinflammation+to+amyloid+toxicity+in+Alzheimer%27s+disease%5BTitle%5D)+AND+%22Journal+of+Neurochemistry%22%5BJournal%5D)">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/pubmed?term=((Microglia+in+Alzheimer%27s+disease%5BTitle%5D)+AND+%22The+Journal+of+Cell+Biology%22%5BJournal%5D)">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/pubmed?term=((Conserved+epigenomic+signals+in+mice+and+humans+reveal+immune+basis+of+Alzheimer%27s+disease%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">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/pubmed?term=((FINEMAP%3A+efficient+variable+selection+using+summary+data+from+genome-wide+association+studies%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">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/pubmed?term=((Identifying+causal+variants+at+loci+with+multiple+signals+of+association%5BTitle%5D)+AND+%22Genetics%22%5BJournal%5D)">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/pubmed?term=((Direct+Identification+of+Hundreds+of+Expression-Modulating+Variants+using+a+Multiplexed+Reporter+Assay%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">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/pubmed?term=((Genome-wide+quantitative+enhancer+activity+maps+identified+by+STARR-seq%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">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/pubmed?term=((The+human+genome+browser+at+UCSC%5BTitle%5D)+AND+%22Genome+Research%22%5BJournal%5D)">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/pubmed?term=((Spatio-temporal+transcriptome+of+the+human+brain%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Spatio-temporal transcriptome of the human brain.</a> Nature. 478 (7370), 483-489 (2011).</li>
</ol>

The authors have nothing to disclose.

Here we describe an analytic framework that can be used to functionally annotate GWS loci based on positional mapping and chromatin interactions. This process involves multiple steps (for more details see this review13). First, given that chromatin interaction profiles are highly cell-type specific, Hi-C data obtained from the appropriate cell/tissue types that best capture underlying biology of the disorder needs to be used. Given that AD is a neurodegenerative disorder, we used adult brain Hi-C data9 to annotate GWS loci. Second, each GWS locus often has up to hundreds of SNPs that are associated with the trait because of linkage disequilibrium (LD), so it is important to obtain putative causal (&#39;credible&#39;) SNPs by computationally predicting the causality through the use of fine-mapping algorithms21,22 or experimentally testing regulatory activities using high-throughput approaches such as massively parallel reporter assays (MPRA)23 or self-transcribing active regulatory region sequencing (STARR-seq)24. For the work described here, we used credible SNPs reported in Jansen et al.14. Third, promoter and exonic SNPs are annotated based on positional mapping. We used a simple positional mapping strategy in which SNPs were mapped to the genes when they overlapped with promoters (defined as 2 kb upstream of transcription start site) or exons. However, this approach can be further elaborated by assessing the functional consequences of exonic SNPs, such as whether the SNP induces nonsense mediated decay, missense variation, or nonsense variation. Fourth, chromatin interaction profiles from the appropriate tissue/cell type can be used to assign SNPs to their putative target genes based on physical proximity. We used interaction profiles anchored to promoters, but we can further refine or expand the interaction profiles by taking enhancer activities (guided by histone H3 K27 acetylation or chromatin accessibility) or exonic interactions into account. One important consideration in this process is to use consistent human genome build. For example, if genomic positions of summary statistics are not based on hg19 (i.e., hg18 or hg38), an appropriate version of the reference genome should be obtained or the summary statistics need to be converted to hg19 using liftover25.
We applied this framework to identify putative target genes for AD GWAS, assigning 284 SNPs to 112 AD risk genes. Using developmental expression profiles26 and cell-type specific expression profiles9, we then demonstrated that this gene set was consistent with what is known about AD pathology, revealing the cell types (microglia), biological functions (immune response and amyloid beta), and elevated risk upon age.
While we presented a framework that delineates potential target genes of AD and its underlying biology, it is of note that Hi-C based annotation can be expanded to annotate any non-coding variation. As more whole-genome sequencing data becomes available and our understanding about the non-coding rare variation grows, Hi-C will provide a key resource for interpretation of disease-associated genetic variants. A compendium of Hi-C resources obtained from multiple tissue and cell types will be therefore critical to facilitating a wide application of this framework to garner biological insights into various human traits and disease.

Genome-wide association studies (GWAS) have played a pivotal role in revealing the genetic basis of a range of human traits and diseases. This large-scale genotyping has uncovered thousands of genomic variants associated with phenotypes ranging from height to schizophrenia risk. However, despite the enormous success of GWAS in identifying disease and trait associated loci, a mechanistic understanding of how these variants contribute to phenotype has been challenging because most phenotype associated variants reside in the non-coding fraction of the human genome. Since these variants often overlap with predicted regulatory elements, they are likely to alter transcriptional control of a nearby gene. However, non-coding loci can influence transcription of genes at linear distances exceeding one megabase, making the genes affected by each variant difficult to identify. Three-dimensional (3D) chromatin structure plays an important role in mediating connections between distant regulatory loci and gene promoters and can be used to identify genes affected by phenotype associated single-nucleotide polymorphisms (SNPs).
Gene regulation is mediated by a complex process, which involves enhancer activation and chromatin loop formation that physically connect enhancers to gene promoters to which the transcriptional machinery can be directed1,2,3. Because chromatin loops often span several hundred kilobases (kb), detailed maps of 3D chromatin architecture are required to decipher gene regulatory mechanisms. Multiple chromatin conformation capture technologies have been invented to identify the 3D chromatin architecture4. Among these technologies, Hi-C provides the most comprehensive architecture, as it captures genome-wide 3D chromatin interaction profiles. Hi-C datasets have been quickly adapted to interpret non-coding genome-wide significant (GWS) loci5,6,7,8,9,10,11,12,13, as it can link non-coding variants to their putative target genes based on chromatin interaction profiles.
In this article, we outline a protocol to computationally predict putative target genes of GWAS risk variants using chromatin interaction profiles. We apply this protocol to map AD GWS loci14 to their target genes using Hi-C datasets in the adult human brain9. The resulting AD risk genes are characterized by other functional genomic datasets that include single cell transcriptomic and developmental expression profiles.

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

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

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

Genome-wide association studies (GWAS) have successfully identified hundreds of genomic loci that are associated with human traits and disease. However, because the majority of the genome-wide significant (GWS) loci fall onto the non-coding genome, the functional impact of many remain unknown. Three-dimensional chromatin interactions identified by Hi-C or its derivatives can provide useful tools to annotate these loci by linking non-coding variants to their actionable genes. Here, we outline a protocol to map GWAS non-coding variants to their putative genes using Alzheimer&#39;s disease (AD) GWAS and Hi-C datasets from human adult brain tissue. Putative causal single-nucleotide polymorphisms (SNPs) are identified by application of fine-mapping algorithms. SNPs are then mapped to their putative target genes using enhancer-promoter interactions based on Hi-C. The resulting gene set represents AD risk genes, as they are potentially regulated by AD risk variants. To garner further biological insights into molecular mechanisms underlying AD, we characterize AD risk genes using developmental brain expression data and brain single-cell expression profiles. This protocol can be expanded to any GWAS and Hi-C datasets to identify putative target genes and molecular mechanisms underlying various human traits and diseases.

The process described here was applied to a set of 800 credible SNPs that were defined by the original study14. Positional mapping revealed that 103 SNPs overlapped with promoters (43 unique genes) and 42 SNPs overlapped with exons (27 unique genes). After positional mapping, 84% (669) SNPs remained unannotated. Using Hi-C datasets in the adult brain, we were able to link an additional 208 SNPs to 64 genes based on physical proximity. In total, we mapped 284 AD credible SNPs to 112 AD risk genes (Figure 1A). AD risk genes were associated with amyloid precursor proteins, amyloid-beta formation, and immune response, reflecting the known biology of AD15,16,17,18 (Figure 1B-D). Developmental expression profiles of AD risk genes showed marked postnatal enrichment, indicative of the age-associated elevated risk of AD (Figure 2A). Finally, AD risk genes were highly expressed in microglia, primary immune cells in the brain (Figure 2B). This is in agreement with the recurrent findings that AD has a strong immune basis and microglia are the central player in AD pathogenesis14,19,20.
<img alt="Alzheimer's genetic study; KEGG, Reactome, GO analysis; SNP distribution; bar charts, gene terms." class="xfigimg" src="/files/ftp_upload/60428/60428fig1b.jpg" data-altupdatedat="1788282455000" data-alt="Figure 1"> 
Figure 1: Defining putative target genes of AD GWS loci. (A) Credible SNPs derived from the top 29 AD loci were categorized into promoter SNPs, exonic SNPs, and unannotated non-coding SNPs. Promoter and exonic SNPs were directly assigned to their target genes by positional mapping, while chromatin interaction profiles in the adult brain were additionally used to map SNPs based on physical interactions. (B-D) Enrichment of GO (B), KEGG (C), and Reactome (D) terms in AD risk genes was performed using HOMER as described in protocol section 6. The x axis represents the false discovery rate (FDR) corrected -log10 (P-value). Enriched terms with FDR &lt; 0.1 were plotted. Grey vertical lines represent FDR = 0.05. APP amyloid precursor protein. Numerator, the number of AD risk genes represented in each term; denominator, the number of genes in each term. <a href="https://www.jove.com/files/ftp_upload/60428/60428fig1alarge.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Gene expression analysis; box plot, violin plot; prenatal vs postnatal brain, AD gene profiles." class="xfigimg" src="/files/ftp_upload/60428/60428fig2.jpg" data-altupdatedat="1788282455000" data-alt="Figure 2"> 
Figure 2: Characterization of AD risk genes. (A) AD risk genes are highly expressed in the postnatal cortex compared to the prenatal cortex. (B) Violin plots depicting distributions of gene expression values (normalized expression) in different cell types from the cortex. These results show that AD risk genes are highly expressed in microglia, consistent with previous studies14. <a href="https://www.jove.com/files/ftp_upload/60428/60428fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Supplementary File 1. <a href="https://www.jove.com/files/ftp_upload/60428/Gencode19_exon.bed" target="_blank">Please click here to view this file (Right click to download).</a>
Supplementary File 2. <a href="https://www.jove.com/files/ftp_upload/60428/Gencode19_promoter.bed" target="_blank">Please click here to view this file (Right click to download).</a>
Supplementary File 3. <a href="https://www.jove.com/files/ftp_upload/60428/geneAnno2.rda" target="_blank">Please click here to view this file (Right click to download).</a>

We present a protocol to identify functional implications of non-coding variants identified by genome-wide association studies (GWAS) using three-dimensional chromatin interactions.

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

mapping-alzheimer-s-disease-variants-to-their-target-genes-using

Research

JoVE Journal

Genetics

Here we are presenting a chromatin immunoprecipitation (ChIP) procedure for genome-wide location analysis of protein isoforms that differ in a histone-binding domain. We are applying it to ChIP-Seq analysis to identify the targets of the KDM5A/JARID1A/RBP2 histone demethylase.

Genome-wide Analysis using ChIP to Identify Isoform-specific Gene Targets

We describe herein an assay by coupling DNA adenine methyltransferase identification (DamID) to high throughput sequencing (DamID-seq). This improved method provides a higher resolution and a wider dynamic range, and allows analyzing DamID-seq data in conjunction with other high throughput sequencing data such as ChIP-seq, RNA-seq, etc.

DamID-seq: Genome-wide Mapping of Protein-DNA Interactions by High Throughput Sequencing of Adenine-methylated DNA Fragments

We present a strategic plan and protocol for identifying non-coding genetic variants affecting transcription factor (TF) DNA binding. A detailed experimental protocol is provided for electrophoretic mobility shift assay (EMSA) and DNA affinity precipitation assay (DAPA) analysis of genotype-dependent TF DNA binding.

Screening for Functional Non-coding Genetic Variants Using Electrophoretic Mobility Shift Assay (EMSA) and DNA-affinity Precipitation Assay (DAPA)

This work describes an optimized methyl-CpG-binding domain (MBD) sequencing protocol and a computational pipeline to identify differentially methylated CpG-rich regions in chronic lymphocytic leukemia (CLL) patients.

Comprehensive DNA Methylation Analysis Using a Methyl-CpG-binding Domain Capture-based Method in Chronic Lymphocytic Leukemia Patients

Here, we describe an innovative method to determine complement receptor 1 (CR1) length polymorphisms for use in several applications, particularly the assessment of susceptibility to diseases such as Alzheimer's disease (AD). This method could be useful to better understand the role of CR1 isoforms in the pathogenesis of AD.

High-resolution Melting PCR for Complement Receptor 1 Length Polymorphism Genotyping: An Innovative Tool for Alzheimer's Disease Gene Susceptibility Assessment

Assay for Transposase-Accessible Chromatin coupled with high-throughput sequencing (ATAC-seq) is a genome-wide method to uncover accessible chromatin. This is a step-by-step ATAC-seq protocol, from molecular to the final computational analysis, optimized for human lymphocytes (Th1/Th2). This protocol can be adopted by researchers without prior experience in next-generation sequencing methods.

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

Here, we describe an optimized high-throughput ChIP-sequencing protocol and computational analyses pipeline for the determination of genome-wide chromatin state patterns from frozen tumor tissues and cell lines.

An Integrated Platform for Genome-wide Mapping of Chromatin States Using High-throughput ChIP-sequencing in Tumor Tissues

Targeted next-generation sequencing is a time- and cost-efficient approach that is becoming increasingly popular in both disease research and clinical diagnostics. The protocol described here presents the complex workflow required for sequencing and the bioinformatics process used to identify genetic variants that contribute to disease.

Targeted Next-generation Sequencing and Bioinformatics Pipeline to Evaluate Genetic Determinants of Constitutional Disease

This protocol describes the Capture Hi-C method used to characterize the 3D organization of megabased-sized targeted genomic regions at high-resolution, including boundaries of topologically associating domains (TADs) and long-range chromatin interactions between regulatory and other DNA sequence elements.

Deciphering High-Resolution 3D Chromatin Organization via Capture Hi-C

A protocol for mapping the three-dimensional genome organization with nucleosome resolution using the genome-wide chromosome conformation capture method Micro-C-XL is presented here.

Mapping Mammalian 3D Genome Interactions with Micro-C-XL

This study analyzed single-nuclei transcriptomes of thirty-three individuals with Alzheimer's disease (AD), revealing sex-specific DEGs in glial cells. Functional enrichment analysis highlighted synaptic, neural, and hormone-related pathways. Key genes, namely NLGN4Y and its regulators, were identified, and potential therapeutic candidates for gender-specific AD were proposed.

Biomarker Identification for Gender Specificity of Alzheimer's Disease Based on the Glial Transcriptome Profiles

Discover how to computationally predict target genes of Alzheimer's disease GWAS risk variants using Hi-C chromatin interaction profiles. This protocol maps SNPs to genes and explores developmental expression patterns.

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

Explore More Videos