diff --git a/xqtl_flowchart.html b/xqtl_flowchart.html index 7402863c..f7e6e5ce 100644 --- a/xqtl_flowchart.html +++ b/xqtl_flowchart.html @@ -175,8 +175,8 @@

FunGen-xQTL Computational Protocol

+ data-out="mol-pheno,geno-prep,fine-map,gwas-integ,enrich" + data-tip="Outputs: ref genome & annotations → Molecular Phenotypes; variant annotations → Genotype Preprocessing; LD reference → Fine-mapping (RSS); LD + gene annotations → GWAS Integration; LD scores → Enrichment (S-LDSC)">
📚
Step 1
Reference Data
Reference genome, gene annotations, LD matrices, and TADB boundaries shared across all downstream steps
@@ -289,8 +289,8 @@

FunGen-xQTL Computational Protocol

+ data-out="cov-prep,qtl-assoc,fine-map" + data-tip="Inputs: variant annotations ← Reference Data Outputs: formatted genotypes → QTL Testing; genotype PCs → Covariate Preprocessing; individual genotypes → Fine-mapping (SuSiE/mvSuSiE/fSuSiE)">
🧪
Step 3a
Genotype Preprocessing
QC raw VCF, convert to PLINK, filter by HWE/missingness, compute kinship and genotype PCs
@@ -321,8 +321,8 @@

FunGen-xQTL Computational Protocol

+ data-out="cov-prep,qtl-assoc,fine-map" + data-tip="Inputs: raw phenotype matrices ← Molecular Phenotypes Outputs: phenotype BED files → QTL Testing; phenotype data → Covariate Preprocessing; individual phenotypes → Fine-mapping (SuSiE/mvSuSiE/fSuSiE)">
📊
Step 3b
Phenotype Preprocessing
Annotate genomic coordinates, impute missing values, and reformat phenotype matrices into BED files for QTL testing
@@ -347,8 +347,8 @@

FunGen-xQTL Computational Protocol

+ data-out="qtl-assoc,fine-map" + data-tip="Inputs: genotype PCs ← Genotype Preprocessing; phenotype data ← Phenotype Preprocessing Outputs: covariate files (PCs + hidden factors) → QTL Testing; covariates → Fine-mapping (SuSiE/mvSuSiE/fSuSiE)">
🔧
Step 3c
Covariate Preprocessing
Merge genotype PCs with phenotypic covariates, then estimate hidden confounders (PEER / Marchenko PCA)
@@ -374,8 +374,8 @@

FunGen-xQTL Computational Protocol

+ data-out="mash,gwas-integ" + data-tip="Inputs: formatted genotypes + phenotype BED + covariate files Outputs: per-context effect sizes & SEs → MASH; QTL sumstats → GWAS Integration (ColocBoost)">
Step 4
QTL Association Testing
Scan all variant–phenotype pairs (cis & trans) with TensorQTL, then apply hierarchical multiple-testing correction
@@ -432,9 +432,9 @@

FunGen-xQTL Computational Protocol

+ data-tip="Inputs: individual genotypes, phenotypes & covariates ← Step 3 (SuSiE/mvSuSiE/fSuSiE); mixture priors ← MASH (mvSuSiE); LD reference ← Reference Data (RSS) Outputs: credible sets + PIPs + TWAS weights → GWAS Integration, Enrichment, EMS">
🎯
Step 5b
Multiomics Regression / Fine-mapping
Pinpoint causal variants within QTL loci using Bayesian credible sets (SuSiE, mvSuSiE, RSS) and generate TWAS weights
@@ -473,9 +473,9 @@

FunGen-xQTL Computational Protocol

+ data-tip="Inputs: credible sets, PIPs, TWAS weights ← Fine-mapping; GWAS sumstats, LD reference & gene annotations ← Reference Data; QTL sumstats ← QTL Association (ColocBoost) Outputs: colocalization posteriors, gene associations → Enrichment">
🔗
Step 6
GWAS Integration
Link QTL signals to GWAS loci via colocalization (SuSiE-enloc, ColocBoost) and gene-level associations (TWAS / cTWAS)
@@ -489,7 +489,7 @@

FunGen-xQTL Computational Protocol

MWETWAS & cTWAS Gene-level trait association using expression weights - + MWEColocBoost Shared-variant discovery across molecular traits @@ -503,7 +503,7 @@

FunGen-xQTL Computational Protocol

📉
@@ -573,9 +573,15 @@

FunGen-xQTL Computational Protocol

['pheno-prep', 'qtl-assoc', 'phenotype BED files', '#0891b2'], ['cov-prep', 'qtl-assoc', 'covariate files', '#0f766e'], ['qtl-assoc', 'mash', 'per-context effect sizes', '#b45309'], - ['qtl-assoc', 'fine-map', 'nominal summary stats', '#b45309'], + ['geno-prep', 'fine-map', 'individual genotypes (SuSiE/mvSuSiE/fSuSiE)', '#1d4ed8'], + ['pheno-prep', 'fine-map', 'individual phenotypes (SuSiE/mvSuSiE/fSuSiE)', '#0891b2'], + ['cov-prep', 'fine-map', 'covariates (SuSiE/mvSuSiE/fSuSiE)', '#0f766e'], + ['ref-data', 'fine-map', 'LD reference (RSS)', '#7c3aed'], ['mash', 'fine-map', 'mixture priors (mvSuSiE)', '#1e40af'], ['fine-map', 'gwas-integ', 'credible sets + TWAS weights', '#be123c'], + ['ref-data', 'gwas-integ', 'GWAS sumstats + LD reference', '#7c3aed'], + ['qtl-assoc', 'gwas-integ', 'QTL sumstats (ColocBoost)', '#b45309'], + ['ref-data', 'enrich', 'LD scores (S-LDSC)', '#7c3aed'], ['fine-map', 'enrich', 'fine-mapped variants', '#be123c'], ['fine-map', 'ems', 'variant features', '#be123c'], ['gwas-integ', 'enrich', 'colocalization posteriors', '#7e22ce'],