While processing 10x Chromium Single Cell Multiome ATAC + Gene Expression sequencing data via ‘cellranger-arc count’ pipeline, Feature linkages analysis is performed as pairs of genomic features, such as peaks and genes, that have significant correlation in signals across cells. Because it provides a basis for inferring enhancer-gene targeting relationships and constructing transcriptional networks. The features with strong linkages are considered to be “co-expressed” and enriched for a shared regulatory mechanism.
Partek Flow provides the opportunity to our users to explore the linkage relationships among different features including peaks and genes, peaks and peaks, and genes and genes. A tab-delimited file containing information of feature linkages inferred from Flow Cell Ranger - ATAC task will be loaded into Integrative Genome Viewer (IGV)[1] for exploration if Feature linkage analysis task has been completed successfully.
Running Feature linkage analysis
To run Feature linkage analysis task (Figure 1),
- Click one datanode that has both features of ATAC and gene expression;
- Click the Feature linkage analysis task under Peak analysis section in the toolbox;
- Click the Finish button to complete the submission.
Figure 1. Feature analysis task in Flow
There will be no inputs needed if the FASTQ is converted to counts matrix within Flow. However, if users processed the FASTQ files outside of Partek, and imported the counts matrix into Flow later. The feature_linkage.bedpe file in outs/analysis/feature_linkage from Cell Ranger output will be needed for each sample (Figure 2) to complete the analysis.
Figure 2. Specify the .bedpe file for project started with counts matrix in Flow.
Task report
In the task report, each pair will generate a list of peaks displayed in a table (Figure 6). Use the drop down menu next to Peaks detected for... to select the pair.
Figure 3. Peaks report on each IP vs control pair
In the report table, each row is a region of a peak and includes the following information:
- Absolute summit: base pair location of peak summit
- Pileup: pileup height at peak summit
- -log10(pvalue): negative log10 pvalue for the peak summit
- Fold enrichment: fold enrichment for the peak summit against random Poisson distribution with local lambda
- -log10(qvalue): negative log10 qvalue at peak summit
- a peak name generated by the MACS2 algorithm
Click the browse to peak button (
) to invoke chromosome view and zoom into that location.
Click the Download button at the lower-right corner to download the peaks in a text file.
References
- Zhang Y, Liu T, et al. Model-based Analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
Additional Assistance
If you need additional assistance, please visit our support page to submit a help ticket or find phone numbers for regional support.
| Your Rating: |
|
Results: |
|
2 | rates |
Overview
Content Tools