---
title: '4DN FOF-CT: Chromatin Tracing Standard'
url: https://www.emergentmind.com/topics/4dn-fish-omics-format-chromatin-tracing-fof-ct
type: topic
---

# 4DN FOF-CT: Chromatin Tracing Standard

The 4DN FISH Omics Format–Chromatin Tracing (FOF-CT) is a community-developed standard for exchanging *processed* chromatin tracing data from multiplexed DNA FISH and FISH-omics experiments. It is designed for ball-and-stick chromatin traces and is implemented as a modular specification for tabular files (TXT/CSV/TSV) that encode per-bright-spot 3D positions, genomic coordinates, trace and cell identifiers, optional RNA-FISH and demultiplexing information, and quality and biological annotations, while remaining directly linkable to raw image data, microscope metadata, and 4DN/OME data infrastructures [2508.13255].

## 1. Scope, motivation, and conceptual domain

FOF-CT arose from a practical gap in spatial genomics data exchange. Before its introduction, laboratories commonly used bespoke formats such as MATLAB structs, spreadsheets, custom HDF5 layouts, and ad hoc CSV files for storing spot coordinates and traces. This impeded cross-study comparison, meta-analysis, and reuse, and often left chromatin tracing datasets inaccessible beyond figure panels. The standard was developed by the 4D Nucleome Imaging Working Group to address these limitations through a minimal mandatory core table and a set of optional but recommended companion tables aligned with OME, 4DN metadata schemas, and 4DN-BINA-OME-QUAREP microscopy metadata standards [2508.13255].

The format operates at the level of processed imaging results rather than raw image acquisition. Existing standards capture microscope configuration, acquisition metadata, and image organization, but do not specify how to represent chromatin traces as ordered sets of loci, how to bind those loci to genomic coordinates, how to connect them to biological annotations such as cell type or nuclear compartments, or how to store FISH-specific demultiplexing and quality-control information. FOF-CT fills this gap for DNA FISH-omics experiments that reconstruct chromatin folding by imaging arrays of genomic segments in single cells across multiple hybridization rounds [2508.13255].

Within the experimental space covered by the specification, current studies generally use one of two representations. In ball-and-stick methods, each targeted segment is represented by a single centroided bright spot per cell per hybridization round; this is the principal focus of FOF-CT. In volumetric methods, loci are represented as clouds of single-molecule localizations. FOF-CT does not yet define a full volumetric representation, but it already includes a Spot Demultiplexing table that stores per-localization events and is described as easy-to-extend to volumetric methods in future revisions [2508.13255].

The methods explicitly named as falling within the intended application space include Wang et al. 2016 chromatin tracing, ORCA, Hi-M, MINA, DNA seqFISH/seqFISH+, OligoFISSEQ, DNA-MERFISH, In situ Genomic Sequencing, and multi-omic extensions combining chromatin with RNA or protein measurements [2508.13255].

## 2. Modular schema and identifier architecture

FOF-CT is modular and comprises up to 12 interlinked tables. Every dataset must contain the DNA-Spot/Trace Core table; the remaining tables are optional or conditionally required depending on assay content and analysis outputs. All tables are plain text files with a header section and a column section. Header lines use three conventions: `#key: value` for human-readable metadata, `##key=value` for mandatory machine-readable metadata, and `#^key: value` for user-defined column descriptions. Relationships across tables are defined by unique identifiers, and the first column of each table is always a UID such as `Spot_ID`, `Trace_ID`, `Cell_ID`, `RNA_Spot_ID`, `Extra_Cell_ROI_ID`, or `Sub_Cell_ROI_ID` [2508.13255].

| Table | Status | Primary role |
|---|---|---|
| DNA-Spot/Trace Core | Mandatory | Per DNA bright spot |
| Spot Demultiplexing | Optional | Per localization event |
| Trace Data | Optional | Trace-level attributes |
| RNA Spot Data | Optional | RNA spot coordinates and IDs |
| Spot Quality | Optional | DNA spot QC |
| RNA_Spot Quality | Optional | RNA spot QC |
| Spot Biological Data | Optional | DNA spot biological annotations |
| RNA_Spot Biological Data | Optional | RNA spot biological annotations |
| Cell Data | Optional | Cell-level properties |
| Extra-Cell ROI Data | Optional | Tissue or organoid ROI context |
| Sub-Cell ROI Data | Optional | Subcellular structures |
| Cell/ROI Mapping | Optional | Explicit ROI boundaries |

The DNA-Spot/Trace Core table is the heart of the format. Each row corresponds to a DNA bright spot in a given cell and typically includes `Spot_ID`, `Trace_ID`, `X`, `Y`, `Z`, `Chrom`, `Chrom_Start`, `Chrom_End`, and `Cell_ID`, with optional `Sub_Cell_ROI_ID`, `Extra_Cell_ROI_ID`, fluorescence-channel information, and additional biological annotations. The Spot Demultiplexing table stores localization-level records such as `Loc_ID`, `Spot_ID`, `X`, `Y`, `Z`, and `Fluor`, allowing reconstruction of per-spot localization clouds and advanced QC. Trace Data stores attributes of a whole chromatin trace, RNA tables mirror the DNA tables for multiplexed RNA FISH, and Cell/ROI tables place traces within cellular, subcellular, tissue, or organoid context [2508.13255].

This organization makes the schema simultaneously minimal and extensible. A dataset can remain compact when only spot-level positions are available, yet it can also capture trace-level state, RNA co-assays, nuclear landmarks, segmentation geometry, and per-localization provenance without breaking identifier consistency [2508.13255].

## 3. Representation of chromatin traces and derived geometric objects

In ball-and-stick chromatin tracing, a trace is the reconstructed 3D path of a specific contiguous genomic region in a single cell and often per allele. Formally, a trace \(t\) is an ordered list
$$
t = \left( (g_1, \mathbf{x}_1), (g_2, \mathbf{x}_2), \dots, (g_N, \mathbf{x}_N) \right),
$$
where each genomic segment \(g_i = (\text{Chrom}_i, \text{Start}_i, \text{End}_i)\) is paired with a spatial coordinate \(\mathbf{x}_i = (x_i, y_i, z_i)\) [2508.13255].

FOF-CT encodes this representation implicitly through the core table. All rows sharing the same `Trace_ID` belong to one trace, and genomic ordering is obtained from `Chrom_Start` and `Chrom_End` or a locus index if present. Trace reconstruction therefore consists of selecting rows with a common `Trace_ID`, sorting them by genomic order, and extracting the \((X,Y,Z)\) coordinates. Coordinate units and frame are supplied in header metadata such as `##coordinate_unit=µm`, and the coordinates are defined with respect to the microscopy detection system, typically after drift and chromatic correction [2508.13255].

The format adopts BED-like genomic coordinates through `Chrom`, `Chrom_Start`, and `Chrom_End`. For modified genomes, it requires a description of the genomic modification, a supplemental VCF file specifying variants, names of inserted or deleted fragments, and targeted-segment coordinates relative to those modifications. This establishes a direct bridge from traced loci to reference and variant genome annotations [2508.13255].

Once traces are reconstructed, standard geometric objects follow directly. For a trace with coordinates \(\mathbf{x}_i\), the single-cell pairwise distance matrix is
$$
d_{ij} = \| \mathbf{x}_i - \mathbf{x}_j \|_2
= \sqrt{ (x_i - x_j)^2 + (y_i - y_j)^2 + (z_i - z_j)^2 }.
$$
A contact map can be defined by thresholding distances at \(d_0\), for example \(250\,\mathrm{nm}\),
$$
C_{ij}^{(c,t)} =
\begin{cases}
1 & \text{if } d_{ij} < d_0 \\
0 & \text{otherwise}.
\end{cases}
$$
Across cells, one obtains mean distance matrices \(\bar d_{ij}\) and mean contact probabilities \(P_{ij}\). The radius of gyration is likewise computed from coordinates as
$$
R_g^2 = \frac{1}{N}\sum_{i=1}^{N}\|\mathbf{x}_i-\mathbf{x}_{\mathrm{cm}}\|^2,
$$
with \(\mathbf{x}_{\mathrm{cm}}\) the trace centroid [2508.13255].

Quality and uncertainty are handled by auxiliary tables rather than a single imposed uncertainty field. Spot Quality includes quantities such as `Centroid_Intensity` and `X_Drift`, `Y_Drift`, `Z_Drift`, and optional metrics such as fit residuals, PSF width, chromatic aberration corrections, and localization uncertainty. Spot Demultiplexing can support localization-cloud-derived uncertainty estimates through the distribution of per-spot localizations. This architecture preserves the processed coordinates needed for downstream modeling while keeping assay-specific QC explicit [2508.13255].

## 4. Population heterogeneity, FISH–Hi-C reconciliation, and model-aware annotation

A central issue for chromatin tracing is that population-average contact assays and imaging-derived distance measurements are not related by a simple monotonic rule once cell populations are heterogeneous. In the Generalized Rouse Model for Chromosomes (GRMC), a homogeneous population yields a one-to-one relation between mean distance \(\langle R_{mn}\rangle\) and contact probability \(P_{mn}\), with the asymptotic scaling
$$
\langle R_{mn}\rangle \approx r_c\,P_{mn}^{-1/3}
$$
for small \(P_{mn}\). Under that expectation, \(P_{mn} < P_{kl}\) implies \(\langle R_{mn}\rangle > \langle R_{kl}\rangle\). The FISH–Hi-C paradox is the empirical observation that some locus pairs violate this ordering: a pair can have higher Hi-C contact frequency but also a larger mean FISH distance [1904.09485].

The proposed resolution is conformational heterogeneity. In a two-subpopulation mixture with fractions \(\eta\) and \(1-\eta\),
$$
\langle R_{mn}\rangle = \eta\langle R_{1,mn}\rangle + (1-\eta)\langle R_{2,mn}\rangle,
$$
and
$$
P_{mn} = \eta P_{1,mn} + (1-\eta)P_{2,mn}.
$$
A pair that is tightly looped in only a subset of cells can therefore have a larger overall contact probability and simultaneously a larger ensemble-average distance than a more uniformly folded comparison pair. Using high-resolution Hi-C on human chromosome 5 and FISH data for eight locus pairs, the authors fit experimental FISH CDFs, inferred \(\eta\), \(\langle R_1\rangle\), and \(\langle R_2\rangle\), and obtained a Pearson correlation of 0.87 between model-derived relative contact frequencies and Hi-C, with robustness for reasonable choices of \(r_c\) such as \(10\)–\(20\,\mathrm{nm}\) [1904.09485].

The same study reports that CTCF loops have \(\langle R_1\rangle \approx 0.25\)–\(0.35\,\mu\mathrm{m}\) almost independent of genomic separation, whereas peak-control pairs have larger distances of approximately \(0.47\)–\(0.67\,\mu\mathrm{m}\). For 212 locus pairs from optical imaging, reconstructed distributions \(P(\langle R\rangle)\) were frequently broad and multi-modal, and the coefficient of variation \(\mathrm{CV}=\sigma(\langle R\rangle)/\mu(\langle R\rangle)\) was frequently \(>0.5\), indicating pervasive heterogeneity across scales from local loops to larger chromosome conformations [1904.09485].

Within a FOF-CT context, this suggests that per-cell coordinates alone are often insufficient when the objective is integration with Hi-C or polymer modeling. The full distribution of distances, empirical CDFs, contact thresholds \(r_c\), sample sizes, and mixture-model quantities such as \(\eta\), \(\langle R_1\rangle\), \(\langle R_2\rangle\), or even continuous \(P(\langle R\rangle)\) reconstructions are analytically relevant metadata. A plausible implication is that FOF-CT-linked workflows are especially valuable when they preserve not only point estimates such as mean or median distances, but also the heterogeneity structure needed to interpret population-average contact maps [1904.09485].

## 5. Spectral loop detection and single-cell loop annotation

Multiplexed super-resolution FISH data support a second class of model-aware analysis in which loops are inferred from full 3D trajectories rather than from pairwise proximity alone. In the spectral framework applied to MERFISH data from Bintu et al. on chromosome 21 in HCT116 cells, each conformation is an ordered sequence of 83 3D points spanning \(34.6\,\mathrm{Mb}\) to \(37.1\,\mathrm{Mb}\) at approximately \(30\,\mathrm{kb}\) resolution, with wild-type and auxin-treated conditions and roughly \(20{,}000\) configurations used to build \(\Lambda\)-plots [2310.15701].

The method treats the polymer as a discrete 3D signal indexed by locus number and applies a discrete cosine transform. Its core observable is the log-spectral ratio
$$
\Lambda(x_n)=\log\left(\frac{\langle x_2^2\rangle^2}{\langle x_1^2\rangle\langle x_4^2\rangle}\right),
$$
where \(\langle x_p^2\rangle\) is the power spectral density at mode \(p\). Internal loops of unknown size and position are detected by scanning sub-chains defined by center \(\iota\) and length \(\eta\), producing a 2D \(\Lambda\)-plot over \((\iota,\eta)\). The theory gives a characteristic maximum at \(\mu_0 \approx 1.34767\) with \(\Lambda_{\max} \approx 2.55882\), so a detected maximum yields a loop midpoint \(\iota_0\), an inferred loop size \(\theta \approx \eta/\mu_0\), and corresponding anchor indices [2310.15701].

The same framework derives loop frequency for mixed populations and complements the population-level spectral analysis with a neural-network classifier for single-cell loop presence. The network takes as input a normalized vector of symmetric distances around an inferred loop midpoint, uses a fully connected hidden layer with ReLU activation and a sigmoid output, and is trained on synthetic looped and non-looped random walks. The reported training sets comprise \(20{,}000\) looped and \(20{,}000\) non-looped samples, with \(5{,}000\) validation and \(2{,}000\) test samples per category, and the ideal synthetic test accuracy is approximately \(90\%\). Applied to MERFISH data, the method identified 14 loops in wild-type cells and 7 loops in auxin-treated cells, and yielded per-loop frequencies, per-cell loop probabilities, loop-positive and loop-negative subsets, average distance maps, subset-specific \(\Lambda\)-plots, ShRec3D-like mean configurations, and loop-loop correlation matrices showing hierarchical, flower-like, and anti-correlated patterns [2310.15701].

These analyses are not themselves part of the mandatory FOF-CT schema, but they map naturally onto FOF-CT-linked data products. The raw requirements of the spectral and neural-network workflow are ordered 3D coordinates, locus identity, genomic order, condition labels, and missing-data flags; all are native to the processed-data layer targeted by FOF-CT. A plausible extension is therefore to accompany core FOF-CT tables with derived spectral outputs such as \(\Lambda\)-plots, loop annotation tables, per-cell loop probabilities, and subset-specific summaries, especially in studies centered on multiplexed super-resolution FISH data [2310.15701].

## 6. Repositories, reuse, interoperability, and future directions

FOF-CT is designed for FAIR dissemination. Datasets are deposited in the 4D Nucleome Data Portal and the OME Image Data Resource (IDR), where they are linked to 4DN experiment and sample metadata, microscope hardware JSON files compliant with 4DN-BINA-OME-QUAREP standards, and raw images in OME-compatible formats. In IDR, FOF-CT-like analysis tables are queryable through a dedicated namespace, and `Spot_ID` can refer to an Image Shape identifier that is directly visualizable in the iViewer [2508.13255].

The scale of reuse is already substantial. Supplemental Table II in the FOF-CT report lists 212 public imaging datasets from 28 publications converted to the format. These include Wang et al. 2016, Bintu et al. 2018, Mateo et al. 2019, Su et al. 2020, Takei et al. 2021, ORCA-based datasets, Hi-M datasets, MINA datasets, and perturbative or disease-oriented studies. Reported use cases include cross-lab comparison of distance matrices, integration with Hi-C and other omics through genomic coordinates, multi-way contact analysis, polymer modeling, domain detection, loop detection, imputation, and reproducible exchange across tools such as pyHiM, traceratops, and CIMA [2508.13255].

Interoperability follows from several design choices. Genomic coordinates use BED conventions, variant handling uses VCF, microscopy metadata is aligned with OME and 4DN-BINA-OME-QUAREP standards, and standardized columns such as `Chrom`, `Cell_ID`, and `Trace_ID` allow search and indexing across repositories. Reusability is reinforced by rich header metadata describing lab, assay conditions, genome build, coordinate units, and analysis software versions, together with the explicit separation between raw imaging assets and processed FOF-CT tables [2508.13255].

Current limitations are also explicit. The present formulation focuses on ball-and-stick representations; volumetric data are only partially covered through Spot Demultiplexing. Metadata completeness can vary across legacy datasets, heterogeneity of imaging platforms complicates uniform representation of localization precision, and genome-scale datasets can make plain-text tables large even when compression and tooling remain practical. Future directions identified for the standard include full volumetric representations, integration of emerging methods such as LoopTrace, GOLD FISH, NOVA, Tn5-FISH, FRET-FISH, and perturb-tracing, stronger linkage to visualization platforms, and continued community extension of the specification and associated tooling [2508.13255].

In practice, FOF-CT functions as a processed-data exchange layer for single-cell chromatin trajectories. Its immediate strength is the standardized representation of spots, traces, RNA co-assays, QC, and spatial context. Its broader significance lies in enabling chromatin tracing datasets to be reused for heterogeneous-population modeling, cross-assay reconciliation, loop calling, imputation, and integrative 4D nucleome analysis without bespoke reverse engineering of every dataset [2508.13255].

Source: https://www.emergentmind.com/topics/4dn-fish-omics-format-chromatin-tracing-fof-ct