14.3 Feature Footprinting

In addition to footprinting motifs, ArchR enables footprinting of any user-defined feature set. To illustrate this capability, we will use the plotFootprints() function to create a TSS insertion profile (introduced previously during the section on data quality control). A TSS insertion profile is just a specialized sub-case of footprinting.

As discussed in the previous section, footprinting is performed using group coverage files which are derived from pseudo-bulk replicates. We originally created these in a previous chapter to perform peak calling. If you haven’t already added group coverages to your ArchRProject, lets do that now.

projHeme5 <- addGroupCoverages(ArchRProj = projHeme5, groupBy = "Clusters2")

We create TSS insertion profiles without normalization for Tn5 bias. The main difference from our previous analyses is that we specify flank = 2000 to extend these footprints 2000 bp on either side of each TSS.

seTSS <- getFootprints(
  ArchRProj = projHeme5, 
  positions = GRangesList(TSS = getTSS(projHeme5)), 
  groupBy = "Clusters2",
  flank = 2000
)

## ArchR logging to : ArchRLogs/ArchR-getFootprints-10b7a47a82723-Date-2020-04-15_Time-11-17-52.log
## If there is an issue, please report to github with logFile!
## 2020-04-15 11:17:52 : Computing Kmer Bias Table, 0.012 mins elapsed.
## 2020-04-15 11:18:24 : Finished Computing Kmer Tables, 0.535 mins elapsed.
## 2020-04-15 11:18:24 : Computing Footprints, 0.547 mins elapsed.
## 2020-04-15 11:18:36 : Computing Footprints Bias, 0.734 mins elapsed.
## 2020-04-15 11:18:43 : Summarizing Footprints, 0.852 mins elapsed.

We can then plot the TSS insertion profiles for each cell group using plotFootprints().

plotFootprints(
  seFoot = seTSS,
  ArchRProj = projHeme5, 
  normMethod = "None",
  plotName = "TSS-No-Normalization",
  addDOC = FALSE,
  flank = 2000,
  flankNorm = 100
)

## ArchR logging to : ArchRLogs/ArchR-plotFootprints-10b7a1b7faa1-Date-2020-04-15_Time-11-18-46.log
## If there is an issue, please report to github with logFile!
## 2020-04-15 11:18:46 : Plotting Footprint : TSS (1 of 1), 0.007 mins elapsed.
## Normalizing by flanking regions
## NormMethod = None
## ArchR logging successful to : ArchRLogs/ArchR-plotFootprints-10b7a1b7faa1-Date-2020-04-15_Time-11-18-46.log