xQTL NIAGADS

About

Overview

About the NIAGADS xQTL Portal

The NIAGADS xQTL Portal provides access to harmonized molecular quantitative trait locus (QTL) results generated from human brain and related multi-omic datasets. The portal is designed to help researchers explore how genetic variation may influence molecular activity across genes, variants, genomic regions, molecular assays, brain regions, and cell-type contexts.

This resource is intended for researchers who want to move from a genetic finding, such as a GWAS locus or variant, to interpretable functional genomic evidence. Users can search by gene, variant, or genomic region; compare associations across QTL types and datasets; view selected results in genomic context; and download summary statistics and metadata for downstream analysis.

The portal integrates xQTL results from major aging-brain and Alzheimer’s disease-relevant studies, including ROSMAP, MSBB, Knight-ADRC, MiGA, and related single-nucleus resources. Results are organized into a consistent framework to make cross-dataset exploration easier while preserving the context needed for careful interpretation.

Team and contributors

This resource is maintained by NIAGADS and the ADSP Functional Genomics Consortium xQTL team in collaboration with contributing investigators, data generators, analysts, and infrastructure teams. We are grateful to the study participants and their families, whose contributions made these datasets possible.

Read the paper

This work is described in-depth in a medRxiv manuscript:
A TAD-informed aging-brain xQTL atlas of multi-modal and cell-type-resolved regulatory variation. DOI: 10.64898/2026.05.21.26353713

What is an xQTL?

An xQTL, or molecular quantitative trait locus, is a statistical association between a genetic variant and a measured molecular trait. In this portal, the “x” refers to different molecular phenotypes, including gene expression, RNA splicing, DNA methylation, histone acetylation, protein abundance, and single-nucleus cell-type-specific expression.

For example, a variant associated with Alzheimer’s disease risk may also be associated with altered expression of a nearby gene, altered splicing of a transcript, differential methylation at a CpG site, altered histone acetylation at a regulatory element, or altered protein abundance. These associations can help prioritize candidate genes, molecular mechanisms, brain regions, and cell types for follow-up study.

An xQTL association should not be interpreted as proof of causality by itself. Association signals can reflect linkage disequilibrium, shared regulatory architecture, correlated molecular traits, study-specific effects, or unmeasured confounding. Portal results are best used for hypothesis generation and should be integrated with GWAS fine-mapping, colocalization, functional annotations, experimental evidence, and biological knowledge.

What can I do with this portal?

The portal supports several common research tasks:

  • Search a GWAS variant to determine whether that variant, or variants in the same region, are associated with molecular traits in human brain datasets.
  • Search a gene to identify variants associated with that gene or related molecular features across datasets, tissues, cell types, and QTL types.
  • Search a genomic region to examine local regulatory evidence across multiple omics types.
  • Compare results across QTL types, including eQTLs, sQTLs, mQTLs, haQTLs, pQTLs, and single-nucleus eQTLs.
  • View selected results in a genome browser to inspect genomic context, nearby genes, regulatory elements, and QTL tracks.
  • Download significant associations, fine-mapped associations, metadata, and larger QTL files for downstream analysis.

The portal is designed to make xQTL results easier to discover, but interpretation still requires attention to dataset, assay, molecular feature, brain region, cell type, ancestry, statistical model, and multiple-testing correction.

Data in the portal

Current release summary

The current release includes harmonized xQTL results across multiple human brain and related datasets.

Category Current release
Brain regions 14
Cell types 7
xQTL types 6
Samples 17,566
Significant xQTL associations 300M+ FDR < 0.05
HMT-significant xQTL associations 83M+
All tested associations 25B+ available through NIAGADS DSS

These counts summarize the scope of the release. Exact sample sizes and available molecular traits may differ by dataset, assay, QTL type, brain region, and cell-type context.

Study design at a glance

The portal combines results from multiple studies and molecular assays (see study-level summary table below). Because each association is context-specific, users should interpret results with the study design in mind.

Dataset Assay Brain region or tissue Cell type, if applicable Sample size QTL types available
Knight-ADRC RNA-seq, DNA methylation array, SomaScan protein aptamers PC Bulk 354-419 eQTL, mQTL, pQTL
MSBB RNA-seq, DNA methylation array, tandem mass tag mass spectroscopy FC BA 10, FC BA 22, FC BA 36, FC BA 44 Bulk 184-274 eQTL, mQTL, pQTL
ROSMAP RNA-seq, DNA methylation array, anti-H3K9ac ChIP-seq, tandem mass tag mass spectroscopy AC, DLPFC, PCC Bulk, monocyte 416-806 eQTL, haQTL, mQTL, pQTL, sQTL
MiGA RNA-seq MTG, STG, SVZ, THA mic 47-66 eQTL
CUIMC1 snRNA-seq DLPFC ast, exc, inh, mic, oli, opc 418-419 snuc-eQTL
MIT snRNA-seq DLPFC ast, exc, inh, mic, oli, opc 377-387 snuc-eQTL
CUIMC1,2+MIT snRNA-seq DLPFC ast, exc, inh, mic, oli, opc 733-737 snuc-eQTL

Users should consult dataset metadata before comparing effects across studies. Differences in assay type, sample composition, brain region, ancestry, disease status, cell-type composition, covariates, and molecular feature definitions can influence whether a signal is detected and how it should be interpreted.

xQTL types in this Atlas

The xQTL Atlas includes multiple classes of molecular QTLs. Each QTL type tests for associations between genetic variants and a specific molecular measurement.

QTL type Molecular trait Typical feature being tested Interpretation
eQTL Bulk gene expression Gene-level expression A variant is associated with expression level of a gene in bulk tissue.
sQTL RNA splicing Splicing event A variant is associated with splicing event.
mQTL DNA methylation methylation probe A variant is associated with methylation level at CG site
haQTL Histone acetylation Histone acetylation peak A variant is associated with histone acetylation levels.
pQTL Protein abundance Protein-level feature A variant is associated with measured protein abundance.
snuc-eQTL Single-nucleus gene expression Gene expression within a cell type A variant is associated with cell-type-specific expression in single-nucleus data.

Not all QTL types use the same feature identifiers. For eQTLs and snuc-eQTLs, the molecular feature is often a gene. For sQTLs, mQTLs, haQTLs, and pQTLs, the molecular feature may be a splice event, CpG site, histone acetylation peak, protein, peptide, or other assay-specific feature. Users should consult the result columns and metadata to determine exactly what molecular trait was tested.

Studies and cohorts in this Atlas

This Atlas includes 7 cohorts, all described in the table below. Note that the three single-nucleus studies are all ROSMAP derived, but have slight differences in individuals included.

Cohort Study inclusion overview Citation
Knight-ADRC The Charles F. and Joanne Knight Alzheimer’s Disease Research Center (Knight-ADRC) Memory and Aging Project cohort at Washington University School of Medicine in St. Louis enrolled individuals aged 65 years and older who demonstrated either no memory impairment or mild dementia at study entry. Yang, C., et al. Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders. Nat Neurosci 24, 1302-1312 (2021).
MiGA The Microglia Genomic Atlas (MiGA) cohort was constructed from postmortem human brain samples obtained from the Netherlands Brain Bank (NBB) and the Neuropathology Brain Bank and Research CoRE at Mount Sinai Hospital, New York Lopes, K.P., et al. Genetic analysis of the human microglial transcriptome across brain regions, aging and disease pathologies. Nat Genet 54, 4-17 (2022).
MSBB The Mount Sinai/JJ Peters VA Medical Center Brain Bank (MSBB) cohort comprises postmortem brain tissue from individuals representing the full spectrum of Alzheimer’s disease (AD) severity, from preclinical to advanced stages. Wang, M., et al. The Mount Sinai cohort of large-scale genomic, transcriptomic and proteomic data in Alzheimer's disease. Sci Data 5, 180185 (2018).
ROSMAP The Religious Orders Study (ROS) and Rush Memory and Aging Project (MAP), collectively known as ROSMAP, is comprised of older individuals without known dementia at enrollment who underwent annual standardized clinical assessments of global cognitive function and domain-specific cognitive performance. De Jager, P.L., et al. A multi-omic atlas of the human frontal cortex for aging and Alzheimer's disease research. Sci Data 5, 180142 (2018).
ROSMAP snRNA-seq sets: CUIMC1, MIT, and CUIMC1,2+MIT Discovery (CUIMC1), replication (MIT), and combined analysis (CUIMC1,2+MIT) sets include cortical cell subsets isolated from individuals in the ROSMAP cohort. These datasets were assembled at the Columbia University Irving Medical Center (CUIMC) in New York City, NY and Massachusetts Institute of Technology (MIT) in Cambridge, MA. Comandate-Lou, Natacha, et al. PLXNBI and other signaling drives a pathologic astrocyte state contributing to cognitive decline in Alzheimer's Disease. bioRxiv (2025).

Types of associations in this Atlas

Following xQTL calling or identification of fine-mapped associations, filtering takes place to prioritize associations with statistical significance. These types of filters, and which ones are kept in the Atlas, are explained below.

Association name Description Track suffix Available in Search?
All All xQTL associations passing quality control. all No, NIAGADS DSS download required.
Benjamini-Hochberg significant xQTLs Derived from the all-significance set. Associations passing per-target Benjamini-Hochberg multiple testing correction (BH FDR <0.05). Synonymous with FDR-significant xQTLs. bh No, NIAGADS DSS download required.
Hierarchical multiple testing significant xQTLs Derived from the all-significance set. Associations passing the hierarchical multiple testing significance threshold, requiring molecular targets with target-level BH FDR <0.05 as well as other significance thresholds. Please see HMT details in the manuscript: Methods: Calling significant xQTL associations using multiple testing correction. hmt Yes
Single-context fine-mapping, all results Fine-mapping output for each context, including all reported fine-mapping results. For m/ha results (via fSuSiE), all records belong to a 95% credible set. For fine-mapping of other molecular traits (via SuSiE), there was no PIP or credible set requirement. scfmAll No, NIAGADS DSS download required.
Single-context fine-mapping, 95% credible set Fine-mapping output restricted to variants assigned to 95% credible sets. Additionally, variant-target pairs were required to have appeared in the corresponding HMT significant file and have a credible set PIP sum of QC’d records of at least 0.95. scfmCs95 Yes

Note: while fine-mapping results may be referred to by their corresponding xQTL-type, e.g., eQTL fine-mapping results, these fine-mapping results are not downstream of xQTL calling and are produced totally independently from the xQTL calling. This is only to indicate that the same molecular traits are being analyzed.

Why harmonization matters

Brain functional genomics datasets are often generated by different groups using different assays, file formats, metadata conventions, genome annotations, statistical models, and analysis pipelines. Without harmonization, it can be difficult to search, compare, or reuse results across datasets.

The NIAGADS xQTL Portal reduces this friction by organizing association results, metadata, annotations, and downloadable files into a consistent framework. Harmonization is intended to make results easier to find and compare, but it does not remove all study-specific differences.

All data provided in the xQTL atlas has been harmonized:

Harmonization item Description
Genome build GRCh38/hg38
Variant identifier format chr:pos:ref:alt
Allele convention non-reference/alternative allele is used as an effect allele, with beta/Z reflecting the effect of the non-reference/alternative allele
Gene annotation Ensembl v103
Molecular feature identifiers gene symbol, ENSG Ensembl gene ID, and target ID (cg ID for mQTLs, UNIPROT ID for pQTLs, LeafCutter splice-event IDs for sQTLs, Ensembl gene ID for eQTLs, and histone acetylation peak ID for haQTLs) annotations are provided for every molecular target
Coordinate system portal is using coordinates 1-based where relevant. Track files use 0-based BED convention
Target normalization Regardless of molecular trait, all outputs include the original QTL target, a mapped gene symbol, and a mapped Ensembl ID. Histone acetylation (haQTLs), DNA methylation (mQTLs), and splicing events (sQTLs) are mapped according to overlapped Gencode gene features as described in Methods.
Covariates age at death, sex, and post-mortem interval (also see Methods)
Multiple testing per-target Benjamini-Hochberg FDR correction and cross-target HMT correction
Fine-mapping SuSiE for e/s/pQTLs and fSuSiE for epigenetic traits (mQTL and haQTL) with 95% credible sets
Metadata standardization single table describing all outputs with output file local path, name, size, coverage, md5, tissue and cell type, QTL type, association/significance level, biosample term id (e.g., CL_0000576, UBERON_0001873), and source publication information

Release notes

V1.0 — First version of NIAGADS xQTL data released on 06/18/2026 (NIAGADS DSS NG00184).

Tutorials and how-to

How to search the portal

Users can search the portal by gene, variant, or genomic region.

Example searches include:

  • Gene: BIN1
  • Variant: rs6733839
  • Variant ID: chr2:127135234:C:T
  • Region: chr2:127048026-127107288

Search behavior should be interpreted using the following rules:

Search type Accepted input Notes
Gene search HGNC symbols Aliases are not currently supported
Variant search rsID or chr:pos:ref:alt GRCh38/hg38 coordinates must be provided
Region search chr:start-end Region coordinates must be provided using 1-based GRCh38/hg38 genomic build. Region size is limited to be at most 500,000 bp

Search results may include associations from multiple datasets, brain regions, cell types, QTL types, and molecular features. Users should filter or stratify results by context before drawing biological conclusions.

Important questions to consider when reviewing search results:

  • Is the queried variant itself associated with the molecular trait, or is the signal driven by nearby variants in linkage disequilibrium?
  • Is the molecular feature a gene, splice event, CpG site, histone acetylation peak, protein, or other assay-specific feature?
  • Which dataset, brain region, and cell-type context support the association?
  • Is the association significant after the appropriate multiple-testing correction?
  • Does the effect direction align with the biological hypothesis?
  • Is there fine-mapping evidence supporting a causal signal?

How to use the genome browser

The genome browser can be used to inspect xQTL results in genomic context. This is useful for understanding whether a signal overlaps nearby genes, regulatory annotations, chromatin features, or other QTL tracks.

When using the genome browser, users should consider:

  • Whether the displayed region contains multiple genes or regulatory elements.
  • Whether multiple QTL types support the same candidate gene or feature.
  • Whether association peaks align across datasets, tissues, or cell types.
  • Whether the queried variant is the strongest association signal or one of many linked variants.
  • Whether fine-mapped variants fall within plausible regulatory annotations.

Browser views are useful for visual interpretation, but downstream analyses should rely on downloaded summary statistics and metadata.

How to download data

The portal supports several download routes depending on your needs.

Goal Recommended download route
Download results for a specific gene, variant, or region Use the Download interface, and enter your variant/gene/region of interest to download its corresponding TSV.
Download significant associations for a query After running a search for your query, navigate to the Downloads sub-tab and click on 'Download all significant xQTL associations for <your query> (TSV)'.
Download fine-mapped associations for a query After running a search for your query, navigate to the Downloads sub-tab and click on 'Download fine-mapped xQTL associations for <your query> (TSV)'.
Download all significant results for a QTL type Use the NIAGADS DSS Open Access Data Portal and filter by dataset NG00184.
Download all tested associations Use NIAGADS DSS Open Access Data Portal, where full association files are available.
Download metadata Go to the Download tab and scroll down. Underneath Download xQTL track metadata you can find the complete metadata in JSON and TSV format.
Understand file contents Download and read the associated README file found in the NIAGADS DSS Open Access Data Portal by filtering the Portal table interface by README in the 'FILE NAME' field.

Downloaded summary files may include fields such as rsID, variantId, targetGene, #chrom, chromStart, chromEnd, Z, beta, p-value, FDR, and other summary statistics. Full tracks and larger association files are available through the NIAGADS DSS Open Access Data Portal.

Because full association files can be large, users should check file size, compression format, indexing, genome build, and README documentation before large-scale downstream analysis.

Understanding results

How to interpret a result

Each xQTL result represents a statistical association between a genetic variant and a molecular trait in a specific context. That context may include the dataset, QTL type, assay, brain region, cell type, molecular feature, statistical model, and multiple-testing correction.

A typical result should be interpreted using the following fields.

Field Meaning Interpretation notes
variantId Variant tested in the association model Variant IDs include chromosome, 1-based position (GRCh38/hg38) and ref:alt alleles.
rsID dbSNP identifier, when available Some variants may have no rsID if a variant is not found in dbSNP.
targetGene and molecular target ID Molecular trait being tested May be a gene, splice event, CpG site, peak, protein, or other feature.
QTL type Molecular phenotype being measured eQTL, sQTL, mQTL, haQTL, pQTL, or snuc-eQTL.
Dataset Source study or analysis dataset Important for sample composition and study design.
Brain region Brain region or tissue context Signals may be context/brain region-specific.
Cell type Cell type Especially important for single-nucleus results.
beta Estimated effect size beta for the non-reference/alternative allele.
Z Test statistic Z-score for the non-reference/alternative allele. Z-score sign is consistent with the sign for beta.
p-value Nominal association evidence Does not account for multiple testing.
FDR Benjamini-Hochberg False discovery rate per-target FDR computed for all variants tested for association with the target
HMT significance Hierarchical multiple-testing result Indicates significance under the hierarchical testing framework (more stringent cross-target correction, see Methods).
Fine-mapping fields SuSiE/fSuSiE 95% credible sets, posterior inclusion probability, conditional effect Prioritizes candidate regulatory variants

Effect direction and allele convention

All effect statistics (beta, Z-score) are reported for non-reference / alternative allele.

Positive beta values indicate that the effect allele is associated with higher measured abundance or activity of the molecular trait. Negative beta values indicate that the effect allele is associated with lower measured abundance or activity. The effect allele is defined as non-reference/alternative allele.

For example, if a variant has a positive eQTL beta for BIN1, that means the specified effect allele is associated with increased BIN1 expression in the reported dataset and context. If the same variant has a negative beta in another brain region or cell type, the difference may reflect tissue specificity, cell-type specificity, or sampling variation and assay differences.

Statistical significance

The portal reports statistical evidence such as p-values, FDR values, HMT significance, and fine-mapping results. These statistics answer related but distinct questions.

  • A p-value describes evidence for association in a specific test.
  • FDR describes significance after per-target multiple-testing correction
  • HMT significance reflects a hierarchical multiple-testing procedure.
  • Fine-mapping statistics prioritize variants that may explain an association signal within a locus.

FDR was calculated for each target separately, e.g. for each gene for eQTLs, each splicing event for sQTLs, or each methylation site for mQTLs.

Recommended interpretation workflow

A typical workflow for interpreting a GWAS locus is:

  1. Search the lead GWAS variant or associated genomic region.
  2. Identify significant xQTL associations in the region.
  3. Check whether associations point to one or more plausible genes or molecular features.
  4. Review the QTL type, dataset, brain region, cell type, and molecular assay.
  5. Examine effect direction and allele convention.
  6. Evaluate statistical support using p-value, FDR, HMT significance, and fine-mapping evidence.
  7. Compare evidence across datasets and molecular phenotypes.
  8. Use the genome browser to inspect regional context.
  9. Download summary statistics and metadata for formal downstream analysis.
  10. Integrate results with GWAS fine-mapping, LD, colocalization, functional annotation, and experimental evidence.

Limitations

The portal is designed to support interpretation and hypothesis generation, but users should keep several limitations in mind.

  • xQTL associations do not prove that a variant causally regulates a molecular trait.
  • A significant xQTL and a disease association in the same region may reflect different causal variants unless supported by colocalization or fine-mapping evidence.
  • Effect sizes may not be directly comparable across datasets, assays, QTL types, tissues, or cell types.
  • Signals can be influenced by sample size, ancestry, disease composition, technical covariates, brain region, postmortem factors, and cell-type composition.
  • Bulk-tissue QTLs may reflect mixtures of cell types.
  • Single-nucleus QTLs may be cell-type-specific but can have smaller sample sizes or different statistical power.
  • Variant identifiers, genome builds, and allele conventions must be checked before integrating portal results with external GWAS or QTL resources.
  • Multiple-testing correction should be interpreted according to the documented correction scope.
  • Fine-mapping results prioritize likely causal variants but do not provide experimental proof.

Users should integrate portal results with independent genetic, molecular, clinical, and experimental evidence before making strong claims about disease mechanisms.

File formats and metadata

The portal provides metadata in JSON and TSV formats, along with variant-, gene-, and region-based summary TSVs. Full association files and track-style files are available through NIAGADS DSS (NG00184).

Search result TSVs

The TSV obtained from the Downloads subtab after performing a Search has the core statistic results as well as some contextual metadata. These are explained below in order of appearance in these TSVs.

Column Description
variantId Variant represented by chromosome, genomic position, reference allele, and alternative allele.
targetGene Gene mapped to the molecular trait as a gene symbol.
xQTLtype Type of QTL analysis (eQTL, snuc-eQTL, haQTL, mQTL, pQTL, sQTL).
context Abbreviation for the brain region, tissue, or cell type context.
study Cohort of origin.
p Nominal p-value for the molecular trait and variant correlation, as calculated using the Student's t-test.*
FDR Multiple-testing adjusted value using Benjamini-Hochberg p-value correction.*
Z Association Z-score for xQTLs, as calculated by beta/se.*
beta Slope of the correlation between the molecular trait and the effect allele.*
se Standard error of the slope of the correlation between the molecular trait and the effect allele.*
target Molecular trait targeted by the variant.
chrom Chromosome of the variant in the selected association.
position 1-based genomic position of the variant in the selected association.
tssDistance Genomic position of associated variant minus the TSS of the target gene.
cisOrTad Whether the association could be detected using a traditional ±1Mb window (cis), or required use of TADB-enhanced (topologically associated domains and adjacent boundaries) testing windows (tad).
rsID dbSNP identifier of variant, when available.
PIP Fine-mapping results only: Posterior inclusion probability, for fine-mapping associations.
conditional_effect Fine-mapping results only: Effect size of the conditional model (e/p/s data) or interpolated effect of variant on the m/ha target.
cs95_size Fine-mapping results only: Number of unique variants in this credible set.
cs95_set_id Fine-mapping results only: Fine-mapping credible set assignment, if available.
context_long Full context definition for the brain region, tissue, or cell type.

* indicates statistic calculated as part of xQTL calling. All fine-mapping statistics are annotated with corresponding xQTL result, so p in a fine-mapping summary TSV refers to its corresponding xQTL association, not a p-value from the fine-mapping analysis.

Comprehensive Atlas metadata

The comprehensive metadata found at the bottom of Downloads has the following fields:

Column Description
Identifier The unique identifier for this file, indicating it's dataset ID, version, chromosome, variant type, and significance level respectively.
Data Source The cohort this file belongs to. One of: Knight-ADRC, MiGA, MSBB, ROSMAP, ROSMAP_CUIMC1, ROSMAP_CUIMC1_2_MIT, and ROSMAP_MIT.
File name The basename of the .bed.gz file.
Number of intervals Total number of unique associations. Note that this Atlas includes empty files.
bp covered Total genomic coverage by this file.
Output type The QTL type, variant type, and association/significance level of this file.
Genome build The reference genome for this file. For this Atlas this is always 'hg38'.
cell type The brain cell type or tissue of origin for this sample.
Biosample type Brief characterization of the isolated sample.
Biosamples term id Ontology term matching the sample's cell type.
Tissue category Broad tissue of origin of the cell type. For this Atlas this is always 'Brain'.
Assay The type of measurement or analysis used to produce this file.
File format Number of BED format fields followed by number of custom fields and the format name.
File size Size of file in bytes as obtained by stat -c %s
Release date The date this file was made public.
Processed file md5 md5 hash of this file.
Link out URL ADSP FunGen xQTL Atlas home page link.
Data Category Broad categorization of this file type. QTL or fine-mapping
Track Description Extra fields including study name and original study pubmed id.
system category The system of the body the Tissue Category belongs to. For this Atlas, always 'Nervous'.
Life Stage Stage of development for this source sample. For this Atlas, always 'Adult'.

Citation and help

Citation and acknowledgment

When using this portal or its data, please cite:

A TAD-informed aging-brain xQTL atlas of multi-modal and cell-type-resolved regulatory variation. DOI: 10.64898/2026.05.21.26353713

Users should also include the required NIAGADS dataset acknowledgment language in publication acknowledgment sections and follow the acknowledgment instructions associated with the relevant NIAGADS dataset record.

Contact

For questions, data issues, or documentation suggestions, please contact help@niagads.org.

When reporting an issue, please include the query used, the QTL type, dataset, browser URL or downloaded file name, and a short description of the unexpected result.