Usecase ea
3. Workflow of GWAS-LD enrichment analysis
GWAS-LD enrichment analysis examines the number of TSSs overlapping with LD blocks on the genome or a given TSS ID list provided by the user (e.g., significant TSS sets from differential analysis). Here we introduce the workflow of this method.
1. Differential analysis is carried out with refTSS as the reference
- TSS tag counts are calculated against the refTSSv4 TSS reference using TSS sequencing data (e.g., CAGE, RAMPAGE, etc.). For instance, the `featureCounts` function in the Subread package can be useful for this process.
- Normalization and statistical testing were performed by standard packages such as DESeq2 and edgeR. For example, edgeR is used to calculate statistical significance between treated and control samples.
- The IDs of TSSs showing significant expression changes are extracted from the testing results. These TSS IDs are input into the refTSS website for GWAS-LD enrichment analysis. Here is an example file of TSS ID lists from CAGE analysis of neural crest cell differentiation experiments with iPS cells ([[1]]).
2. GWAS-LD enrichment analysis
From the menu, select refTSS-LD and paste the TSS IDs into the box.
3. Interpreting the results of the enrichment analysis
The table shows the results of the GWAS-LD analysis conducted using the example data. Results are sorted by order of FDR. The enriched TSS IDs are displayed in the "hits" column, allowing for inspection of which TSSs are concentrated in LD blocks grouped by traits.
4. Workflow of general annotation enrichment analysis
GO and KEGG pathway annotation enrichment analysis examines the over-representation of genes with biological annotations such as GO and KEGG pathways. Here we introduce the workflow of general annotation enrichment analysis using refTSS and external enrichment analysis tools.
1. General annotation enrichment analysis is carried out with refTSS as the reference
- TSS tag counts are calculated against the refTSSv4 TSS reference using TSS sequencing data (e.g., CAGE, RAMPAGE, etc.). The `featureCounts` function in the Subread package can be useful for this process.
- Normalization and statistical testing were performed by standard packages such as DESeq2 and edgeR. For example, edgeR is used for statistical significance calculations between treated and control samples.
- The IDs of TSSs showing significant expression changes are extracted from the differential expression analysis results.
- The TSS IDs are input into the ID conversion tool on the refTSS website.
- The extracted Gene Symbols are input into external enrichment analysis tools such as DAVID, MetaScape, and ClusterProfiler.