A successful FFPE-CUTAC reaction produces a DNA library of approximately 200 base pair fragments, corresponding to adapter-ligated FFPE-CUTAC fragments. This product forms a peak on the capillary electrophoresis trace, distinct from the primer-dimer region below 100 base pairs and high-molecular-weight background above 800 base pairs (Figure 2). Libraries with a visible 200 bp peak but low overall signal are often still suitable for sequencing and produce sufficient mapped reads for genome-wide analysis. Lane 3, which showed very low yield by capillary electrophoresis, nonetheless produced 18.5 million raw reads (Table 1). Lane 5, derived from a 90-week-old mouse colon with a larger tissue input, produced a stronger electrophoresis signal and higher mapped read fraction, illustrating the range of outcomes across input conditions. A library showing no detectable peak above the primer-dimer region, or one dominated by a high-molecular-weight smear without a discrete product, should be considered potentially failed; however, sequencing may still be attempted and can be successful if there were higher-yield reactions in the same experiment.
Samples were pooled according to the measured library concentrations, targeting a sequencing depth of 10–20 million fragments per sample, and subsequently submitted for sequencing on a high-throughput sequencer in dual-index paired-end mode. Sequencing reads from a successful library mapped predominantly to the reference genome, with the fraction of reads aligning to the nuclear genome typically exceeding 80% and only a small minority of reads mapping to mitochondrial DNA or to the Escherichia coli genome (Table 1). The latter derives from residual E. coli DNA carried over during pAg-Tn5 production and, because the same volume of enzyme is added to each reaction within an experiment, can serve as a proxy for sequencing normalization across samples. After alignment, mapped fragments accumulated as sharp peaks at the transcription start sites of highly expressed genes, consistent with enrichment of RNA Polymerase II at active regulatory elements rather than background (Figure 3).
Table 1 shows sequencing statistics for replicate samples. A reaction was considered unsuccessful if the corresponding sequencing data showed a low fraction of genome-mapped reads (mm10_Fraction_Mapped below approximately 0.5), a chrM_Fraction above approximately 0.1, the absence of enrichment at known active regulatory elements, or a fragment-size distribution dominated by low- or high-molecular-weight material. Libraries with Estimated Library Size (ELS) values below approximately 4 million may have limited complexity and yield few unique fragments, even with deeper sequencing. Browser tracks for these samples showed peaks in reads at the highly abundant long non-coding RNAs Malat1 and Neat1.
Further quality control can be performed post-sequencing using computational tools from open source packages to visualize and evaluate signal quality. UCSC bigWigMerge11 can be used to combine normalized coverage tracks into a single merged track. The deepTools suite12 provides a range of QC tools, including plotCorrelation for replicate correlation matrices, plotPCA for principal component analysis across samples, computeMatrix and plotHeatmap for visualizing signal enrichment around annotated features like transcription start sites, and bamCompare for generating ratio tracks between conditions. Fractions of Reads in Peaks (FRiP) scores, calculated with tools such as BEDTools13 intersect or MACS214, provide a measure of library specificity. Finally, fragment size distribution plots generated from mapped paired-end reads should confirm that fragments in the 50-200 base pair range are enriched.

Figure 1: Overview of the protocol for FFPE-CUTAC. The workflow for Formalin-Fixed Paraffin-Embedded Cleavage Under Targeted Accessible Chromatin (FFPE-CUTAC) begins with deparaffinizing, followed by de-crosslinking of formalin-fixed, paraffin-embedded (FFPE) tissues at high heat. The protein of interest is targeted using a primary antibody and a secondary antibody, which is then conjugated to protein A-Tn5 transposase. The slide is incubated with magnesium to activate the Tn5, then the tissue is scraped into a PCR tube. The cells are incubated with Proteinase K to release DNA, then amplified and sequenced. Abbreviations: FFPE = formalin-fixed, paraffin-embedded; CUTAC = Cleavage Under Targeted Accessible Chromatin; pAg-Tn5 = protein A/G-Tn5 transposase; PCR = polymerase chain reaction; TLPK = thermolabile proteinase K. Please click here to view a larger version of this figure.

Figure 2: Representative capillary electrophoresis image. Gel analysis of libraries generated after PCR amplification and assessed by capillary electrophoresis before high-throughput sequencing. The electronic ladder on the left provides the reference for estimating fragment size, with successful reactions generally yielding a visible band at approximately 200 base pairs. Lanes 1 through 4 represent technical replicates of a 12-week-old mouse colon, with 1 having 50% of tissue input, 2 having 30%, 3, having 10%, and 4 having the remaining tissue. The library in lane 5 came from a larger tissue section of a 90-week-old mouse colon for comparison. Although occasional reactions showed little visible yield (see lane 3), they often produced sufficient sequencing for efficient genome mapping. Reactions with strong bands (see lane 5) tend to have a higher mapped fraction of reads and a larger library. Abbreviations: PCR = polymerase chain reaction; bp = base pairs. Please click here to view a larger version of this figure.

Figure 3: Genome track visualization. RNA Polymerase II (RNAPII) enrichment at promoters for highly abundant nuclear long non-coding RNA Neat1 and Malat1. Peaks form at the promoters of active genes because there is high RNAPII signal. Abbreviation: RNAPII = RNA polymerase II. Please click here to view a larger version of this figure.
| Sample | Conc. [pg/µL] | Raw_Reads | Clipped_Reads | mm10_fragments | chrM_fraction | ecoli_Fragments | mm10_Fraction_Mapped | Fraction_Duplication | estimated_Library_Size |
| 1 | 50.4 | 17022392 | 16944511 | 14106793 | 0.034 | 13157 | 0.833 | 0.423147 | 12158636 |
| 2 | 40.1 | 17652810 | 17563896 | 14618178 | 0.033 | 13414 | 0.832 | 0.513631 | 9082386 |
| 3 | 3.85 | 18471771 | 18354167 | 15240815 | 0.033 | 14311 | 0.83 | 0.755529 | 3812279 |
| 4 | No peak recorded | 17241161 | 17143339 | 14222035 | 0.033 | 12015 | 0.83 | 0.562674 | 7377301 |
| 5 | 390 | 9269872 | 9233822 | 8533086 | 0.017 | 7856 | 0.924 | 0.180875 | 26006109 |
Table 1: Representative sequencing statistics. Raw_Reads is the total number of sequencing reads produced before any processing steps. Clipped_Reads is the number of reads retained after adapter sequences are trimmed off. Mm10_Fragments is the number of read pairs that map to the mm10 mouse reference genome. ChrM_Fraction is the proportion of total mapped fragments that align to the mitochondrial chromosome. Ecoli_Fragents is the number of read pairs mapped to the Escherichia coli reference genome, generally carried along with the bacterially-produced pAg-Tn5 protein. Residual E. coli fragments associated with pAg-Tn5 production can be used as a proxy for normalization across samples processed with the same amount of enzyme. Mm10_Fraction_Mapped is the proportion of reads that map successfully to the mm10 mouse reference genome. Fraction_Duplication is duplicate reads, inferred by the MarkDuplicates (Picard) package. ELS is the estimated library size, calculated by MarkDuplicates (Picard). Abbreviations: Conc. = concentration; pg/µL = picograms per microliter; mm10 = mouse reference genome assembly mm10; chrM = mitochondrial chromosome; E. coli = Escherichia coli; ELS = estimated library size.