Sq vs. K Taxonomic Classification

• Sq taxonomic classification is defined by aligned sequence models that require global alignment to leverage position-specific evolutionary signals.
• K classification models use k-mer count vectors as a bag-of-words, offering flexibility for highly variable or unalignable regions.
• Bayesian nonparametric frameworks like BayesANT integrate both paradigms, balancing statistical power, novelty detection, and computational efficiency.

The Sq (aligned sequence) and K (k-mer count) taxonomic classification frameworks constitute two principal paradigms for modeling molecular sequence data in computational taxonomy. Their distinction is foundational, directly determining the likelihood function, prior structure, inference algorithms, and suitability to varying molecular datasets and application contexts. Theoretical analyses and empirical work—including those of BayesANT (Zito et al., 2022), and broader k-mer–based classifiers—demonstrate not only the diverse implementation requirements but also notable differences in calibration, interpretability, statistical power, and robustness to sequence variability.

1. Formal Definitions: Sq and K Representations

The Sq approach, or "aligned-sequence model," represents each DNA or RNA sequence as a vector of nucleotides at $p$ globally aligned sites. For a sequence $X=X_1,\ldots,X_p$ with $X_s\in\{A,C,G,T,-\}$, the underlying generative model assumes

$$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$

where $\theta_{v,s,\cdot}$ is a categorical distribution at site $s$ for taxon $v$.

By contrast, the K approach represents each sequence by its vector of observed $\kappa$-mer counts over the $4^\kappa$ possible canonical substrings of length $\kappa$, denoted $X=X_1,\ldots,X_p$0. The associated multinomial likelihood is

$X=X_1,\ldots,X_p$1

with no explicit positional information; $X=X_1,\ldots,X_p$2 parameterizes the expected k-mer composition within taxon $X=X_1,\ldots,X_p$3.

2. Bayesian Nonparametric Framework: BayesANT and Issues of Taxonomic Novelty

BayesANT (Zito et al., 2022) instantiates both Sq and K kernels within a Bayesian nonparametric hierarchy, employing Pitman–Yor species sampling priors at each internal branch of the taxonomy. For every parent node $X=X_1,\ldots,X_p$4 at rank $X=X_1,\ldots,X_p$5,

• The Pitman–Yor process hyperparameters $X=X_1,\ldots,X_p$6 control the probability assigned to previously seen versus new child taxa (see allocation rules below):

$X=X_1,\ldots,X_p$7

where $X=X_1,\ldots,X_p$8 is the number of sequences assigned so far to child $X=X_1,\ldots,X_p$9 under parent $X_s\in\{A,C,G,T,-\}$0, and $X_s\in\{A,C,G,T,-\}$1 is the current number of distinct children. This prior guarantees every query is assigned a nontrivial posterior probability of belonging to a novel taxon at each rank.

In both Sq and K settings, Dirichlet-multinomial conjugacy allows analytic marginalization over $X_s\in\{A,C,G,T,-\}$2, yielding closed forms for predictive probabilities at both known and novel leaves.

3. Likelihoods, Posteriors, and Predictive Inference

In the Sq model, with Dirichlet priors $X_s\in\{A,C,G,T,-\}$3, the marginal likelihood for a new sequence at leaf $X_s\in\{A,C,G,T,-\}$4 is

$X_s\in\{A,C,G,T,-\}$5

where $X_s\in\{A,C,G,T,-\}$6 aggregates training counts at site $X_s\in\{A,C,G,T,-\}$7 in $X_s\in\{A,C,G,T,-\}$8.

For the K model, with Dirichlet priors on the k-mer probability vector $X_s\in\{A,C,G,T,-\}$9, the posterior predictive for a new sequence is

$$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$0

where $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$1 is the observed count of k-mer $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$2 in $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$3's assigned training sequences.

At prediction time, BayesANT computes the posterior for each branch, calibrates via exponentiation ($$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$4, $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$5), and aggregates over leaves at each rank.

4. Comparative Methodological and Practical Considerations

Alignment Dependency

• Sq methods require globally aligned sequences (e.g., COI barcodes 658 bp in FinBOL), delivering site-specific priors and leveraging positionally informative substitutions.
• K methods make no alignment assumptions; k-mer counts or compositions serve as a bag-of-words over local motifs, critical in highly variable or unalignable regions (e.g., fungal ITS).

Statistical Power and Parameterization

• The Sq approach parameterizes $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$6 probabilities per leaf, where $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$7 is the sequence length. This yields highly parsimonious models in well-aligned marker regions.
• The K approach grows exponentially in $$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$8 ($$P(X\mid\theta_v) = \prod_{s=1}^p \prod_{g\in\{A,C,G,T\}} \theta_{v,s,g}^{\mathbb{I}(X_s=g)},$$9-dimensional multinomial). Large $\theta_{v,s,\cdot}$0 leads to sparsity and risk of overfitting, while small $\theta_{v,s,\cdot}$1 (e.g., 2–5) offers a compromise between expressivity and sample complexity.

Predictive Performance and Discovery of Novel Taxa

Empirical results in (Zito et al., 2022) show, for the ornithopter segment of the Canadian barcode library:

• With the aligned Sq kernel, species-level identification reaches 85.2% (random split, S1) or 70.6% (half novel, S2), and BayesANT recognizes 77.9–93.8% of truly novel species.
• The m-2 (sitewise 2-mer) kernel yields indistinguishable performance ($\theta_{v,s,\cdot}$20.5% difference), indicating only moderate gains from per-site higher-order motifs over pure position-independence.
• The K kernel achieves similar performance for small $\theta_{v,s,\cdot}$3 (not explicitly benchmarked in (Zito et al., 2022)), and is standard when alignment is inapplicable.

A plausible implication is that for loci permitting reliable global alignment (e.g., animal COI), the Sq kernel is preferable, offering compactness and leveraging sitewise phylogenetic signal. For taxonomic groups and markers where alignment is infeasible (ITS, hypervariable regions), K is the default approach.

5. Extensions and Practical Implementations

Mixture Models, Trees, and Machine Learning Integration

Beyond BayesANT, k-mer–based (K) taxonomic classification underpins methods such as MetaPalette (Koslicki et al., 2016), which employs "palette painting"—modeling samples as sparse mixtures over reference and hypothetical k-mer palettes, exploiting $\theta_{v,s,\cdot}$4 as large as $\theta_{v,s,\cdot}$5–$\theta_{v,s,\cdot}$6 to achieve high strain specificity. Classification is then solved by penalized nonnegative least squares, recovering abundances of both known and novel representatives.

Resource-efficient implementations, such as those evaluated in (Fuhl et al., 2023), further process normalized k-mer frequency vectors using standard machine learning classifiers—subspace k-NN, shallow neural nets, bagged decision trees—yielding genus-level precision to $\theta_{v,s,\cdot}$7 with sub-millisecond per-read runtimes and memory footprints down to a few MB, far below alignment-based alternatives.

Hybrid and Compression-Aware Models

MEM-based approaches (sequence-based or Sq) such as those in (Draesslerová et al., 2024) extend the FM-index paradigm to arbitrarily long exact matches (not restricted to fixed k), allowing more flexible, alignment-free sequence comparison. Hybrid indices employing lossy compression (e.g., KATKA kernels, minimizer digests) achieve substantial reductions in storage and indexing time at minimal accuracy costs, narrowing the gap between pure K-mer and full-sequence Sq approaches.

6. Scope, Limitations, and Selection Guidelines

• Sq is preferred when: robust global alignment is possible, sitewise evolutionary signals are informative, and the repertoire of taxa is relatively stable (i.e., low likelihood of unclassifiable reads).
• K is necessary when: alignment is infeasible or not meaningful, sequences display high indel or motif variability, or when computational/algorithmic simplicity is needed for very large-scale datasets.
• For both approaches, Bayesian nonparametric priors (Pitman–Yor or Dirichlet processes) enable principled discovery and quantification of novel taxa at every rank, avoiding negative transfer or overconfidence due to hard assignment.
• Parameter choices (alignment reference, $\theta_{v,s,\cdot}$8, Dirichlet concentration, feature preprocessing) crucially impact both statistical accuracy and computational trade-offs.

7. Summary Table: Sq vs. K Classification Paradigms

Aspect	Sq (Aligned sequence)	K (k-mer count)
Data structure	Positionwise nucleotide	k-mer frequency/count vector
Alignment required	Yes	No
Feature dimension	$\theta_{v,s,\cdot}$9 ($s$0=length)	$s$1
Parametric model	Product-multinomial	Multinomial
Sensitive to indels	Yes	No
Performance (COI)	85.2%–70.6% (species)	Comparable (for small $s$2)
Novelty detection	Supported (BayesANT)	Supported (BayesANT)
Computational cost	Moderate (if $s$3 moderate)	Grows rapidly with $s$4

Selection between Sq and K should be grounded in the molecular properties of the marker loci, the need for alignment, and the trade-off between computational resources and classification accuracy as demonstrated in empirical benchmarking studies (Zito et al., 2022, Koslicki et al., 2016, Fuhl et al., 2023, Draesslerová et al., 2024).