docs: rename AxoMEME to HyphAeon in model card
Browse files
README.md
CHANGED
|
@@ -1,6 +1,6 @@
|
|
| 1 |
---
|
| 2 |
license: mit
|
| 3 |
-
library_name:
|
| 4 |
tags:
|
| 5 |
- biology
|
| 6 |
- bioinformatics
|
|
@@ -11,19 +11,19 @@ tags:
|
|
| 11 |
- onnx
|
| 12 |
---
|
| 13 |
|
| 14 |
-
#
|
| 15 |
|
| 16 |
**Ultra-fast neural inference of episodic positive selection in molecular sequences.**
|
| 17 |
|
| 18 |
-
|
| 19 |
arbitrary codon alignments and phylogenetic trees. It replaces numerical Maximum Likelihood
|
| 20 |
Estimation with an axial geometric transformer using Continuous 4D Tree Rotary Position
|
| 21 |
Embeddings (Tree-RoPE), reaching a 100x-1,000x speedup over HyPhy MEME and PAML CodeML while
|
| 22 |
matching or exceeding their empirical statistical power at a strictly controlled false
|
| 23 |
positive rate.
|
| 24 |
|
| 25 |
-
- **Repository:** https://github.com/veg/
|
| 26 |
-
- **Architecture notes:** [ARCHITECTURE.md](https://github.com/veg/
|
| 27 |
- **Developed by:** Sergei L. Kosakovsky Pond and collaborators
|
| 28 |
- **License:** MIT
|
| 29 |
|
|
@@ -33,13 +33,13 @@ The same model weights power three analyses, exposed as CLI subcommands:
|
|
| 33 |
|
| 34 |
| Command | Analysis | Output |
|
| 35 |
| :--- | :--- | :--- |
|
| 36 |
-
| `
|
| 37 |
| `epistasis` | **Epistatic sector mining (ESSM)**: phylogenetic branch attribution, inter-site co-selection networks, in-silico Selection Deep Mutational Scanning, and multi-scale epistatic sectors. | co-selection edges, sectors, DMS sweep (JSON / GraphML / CSV) |
|
| 38 |
| `phenotype` | **Directional phenotype-genotype association (PhyloWAS)**: associates sites with a trait across the tree and extracts Phenotype-Associated Residue Signatures (PARS). Supports curated phenotype presets (marine, echolocation, high-altitude, longevity, …) or a user trait table. | per-site trait association, p-values, FDR |
|
| 39 |
|
| 40 |
`epistasis` and `phenotype` are **derived analyses built on the selection model** (its per-site drive and
|
| 41 |
branch-level attributions); they do not require separate weights. The results below and the variants
|
| 42 |
-
concern the core `
|
| 43 |
|
| 44 |
## Model variants
|
| 45 |
|
|
@@ -76,27 +76,27 @@ standard MEME null mixture, `0.5 * delta(0) + 0.5 * chi^2(1)`.
|
|
| 76 |
Install the package, which pulls these weights automatically on first use:
|
| 77 |
|
| 78 |
```bash
|
| 79 |
-
pip install git+https://github.com/veg/
|
| 80 |
```
|
| 81 |
|
| 82 |
### Command line
|
| 83 |
|
| 84 |
```bash
|
| 85 |
-
|
| 86 |
--alignment my_genes.fasta \
|
| 87 |
--tree my_genes.nwk \
|
| 88 |
--output results.json \
|
| 89 |
--csv results.csv
|
| 90 |
```
|
| 91 |
|
| 92 |
-
If `--tree` is omitted,
|
| 93 |
|
| 94 |
### Python
|
| 95 |
|
| 96 |
```python
|
| 97 |
-
from
|
| 98 |
|
| 99 |
-
model = PhyloAxialTransformer.from_pretrained("datamonkey/
|
| 100 |
codons, aas, dists, mds, invariable, taxa, num_sites = load_alignment_and_tree(
|
| 101 |
"my_genes.fasta", "my_genes.nwk"
|
| 102 |
)
|
|
@@ -107,7 +107,7 @@ To load the **viral fine-tuned** variant explicitly:
|
|
| 107 |
```python
|
| 108 |
import torch
|
| 109 |
from safetensors.torch import load_file
|
| 110 |
-
from
|
| 111 |
|
| 112 |
model = PhyloAxialTransformer(embed_dim=384, num_layers=6, num_heads=12,
|
| 113 |
window_size=1, num_thresholds=16)
|
|
@@ -133,7 +133,7 @@ flag. Sites that are invariable across the alignment are assigned an LRT of 0.
|
|
| 133 |
Taking HyPhy MEME (p <= 0.10, asymptotic LRT >= 4.605) as ground truth across 84 empirical
|
| 134 |
datasets from 9 independent literature studies (43,302 codons, up to 476 taxa):
|
| 135 |
|
| 136 |
-
| Literature study & system | Datasets | Codons | ROC-AUC | PR-AUC | PPV | FPR | Spearman rho | HyPhy MLE |
|
| 137 |
| :--- | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: |
|
| 138 |
| Abdul et al. (2018), SMC5/6 complex | 9 | 7,073 | 0.990 | 0.752 | 100.0% | 0.00% | 0.874 | 561.0 s | 3.84 s | 171.3x |
|
| 139 |
| Nisson et al. (2025), CCDC137 (HIV Vpr) | 1 | 290 | 0.940 | 0.650 | 75.0% | 0.35% | 0.833 | 80.0 s | 0.21 s | 380.9x |
|
|
@@ -193,7 +193,7 @@ so browser or cross-language ports should account for this.
|
|
| 193 |
|
| 194 |
## Limitations and intended use
|
| 195 |
|
| 196 |
-
|
| 197 |
sites under episodic diversifying selection.
|
| 198 |
|
| 199 |
- It is trained to approximate HyPhy MEME, so it inherits that method's assumptions and is
|
|
@@ -215,10 +215,10 @@ sites under episodic diversifying selection.
|
|
| 215 |
A publication is in preparation. Until it appears, cite the repository:
|
| 216 |
|
| 217 |
```bibtex
|
| 218 |
-
@software{
|
| 219 |
-
title = {
|
| 220 |
author = {Kosakovsky Pond, Sergei L. and collaborators},
|
| 221 |
-
url = {https://github.com/veg/
|
| 222 |
year = {2026}
|
| 223 |
}
|
| 224 |
```
|
|
|
|
| 1 |
---
|
| 2 |
license: mit
|
| 3 |
+
library_name: hyphaeon
|
| 4 |
tags:
|
| 5 |
- biology
|
| 6 |
- bioinformatics
|
|
|
|
| 11 |
- onnx
|
| 12 |
---
|
| 13 |
|
| 14 |
+
# HyphAeon
|
| 15 |
|
| 16 |
**Ultra-fast neural inference of episodic positive selection in molecular sequences.**
|
| 17 |
|
| 18 |
+
HyphAeon detects site-level episodic diversifying positive selection (dN/dS > 1) across
|
| 19 |
arbitrary codon alignments and phylogenetic trees. It replaces numerical Maximum Likelihood
|
| 20 |
Estimation with an axial geometric transformer using Continuous 4D Tree Rotary Position
|
| 21 |
Embeddings (Tree-RoPE), reaching a 100x-1,000x speedup over HyPhy MEME and PAML CodeML while
|
| 22 |
matching or exceeding their empirical statistical power at a strictly controlled false
|
| 23 |
positive rate.
|
| 24 |
|
| 25 |
+
- **Repository:** https://github.com/veg/HyphAeon
|
| 26 |
+
- **Architecture notes:** [ARCHITECTURE.md](https://github.com/veg/HyphAeon/blob/main/ARCHITECTURE.md)
|
| 27 |
- **Developed by:** Sergei L. Kosakovsky Pond and collaborators
|
| 28 |
- **License:** MIT
|
| 29 |
|
|
|
|
| 33 |
|
| 34 |
| Command | Analysis | Output |
|
| 35 |
| :--- | :--- | :--- |
|
| 36 |
+
| `meme` | **Episodic positive selection**: per-site MEME LRT surrogate (the core model). | per-site LRT + p-value |
|
| 37 |
| `epistasis` | **Epistatic sector mining (ESSM)**: phylogenetic branch attribution, inter-site co-selection networks, in-silico Selection Deep Mutational Scanning, and multi-scale epistatic sectors. | co-selection edges, sectors, DMS sweep (JSON / GraphML / CSV) |
|
| 38 |
| `phenotype` | **Directional phenotype-genotype association (PhyloWAS)**: associates sites with a trait across the tree and extracts Phenotype-Associated Residue Signatures (PARS). Supports curated phenotype presets (marine, echolocation, high-altitude, longevity, …) or a user trait table. | per-site trait association, p-values, FDR |
|
| 39 |
|
| 40 |
`epistasis` and `phenotype` are **derived analyses built on the selection model** (its per-site drive and
|
| 41 |
branch-level attributions); they do not require separate weights. The results below and the variants
|
| 42 |
+
concern the core `meme` model (aliases: `predict`, `site-selection`), which the other two consume.
|
| 43 |
|
| 44 |
## Model variants
|
| 45 |
|
|
|
|
| 76 |
Install the package, which pulls these weights automatically on first use:
|
| 77 |
|
| 78 |
```bash
|
| 79 |
+
pip install git+https://github.com/veg/HyphAeon.git
|
| 80 |
```
|
| 81 |
|
| 82 |
### Command line
|
| 83 |
|
| 84 |
```bash
|
| 85 |
+
hyphaeon meme \
|
| 86 |
--alignment my_genes.fasta \
|
| 87 |
--tree my_genes.nwk \
|
| 88 |
--output results.json \
|
| 89 |
--csv results.csv
|
| 90 |
```
|
| 91 |
|
| 92 |
+
If `--tree` is omitted, HyphAeon extracts an embedded tree from a NEXUS or FASTA alignment.
|
| 93 |
|
| 94 |
### Python
|
| 95 |
|
| 96 |
```python
|
| 97 |
+
from hyphaeon import PhyloAxialTransformer, load_alignment_and_tree
|
| 98 |
|
| 99 |
+
model = PhyloAxialTransformer.from_pretrained("datamonkey/hyphaeon").eval()
|
| 100 |
codons, aas, dists, mds, invariable, taxa, num_sites = load_alignment_and_tree(
|
| 101 |
"my_genes.fasta", "my_genes.nwk"
|
| 102 |
)
|
|
|
|
| 107 |
```python
|
| 108 |
import torch
|
| 109 |
from safetensors.torch import load_file
|
| 110 |
+
from hyphaeon import PhyloAxialTransformer
|
| 111 |
|
| 112 |
model = PhyloAxialTransformer(embed_dim=384, num_layers=6, num_heads=12,
|
| 113 |
window_size=1, num_thresholds=16)
|
|
|
|
| 133 |
Taking HyPhy MEME (p <= 0.10, asymptotic LRT >= 4.605) as ground truth across 84 empirical
|
| 134 |
datasets from 9 independent literature studies (43,302 codons, up to 476 taxa):
|
| 135 |
|
| 136 |
+
| Literature study & system | Datasets | Codons | ROC-AUC | PR-AUC | PPV | FPR | Spearman rho | HyPhy MLE | HyphAeon CPU | Speedup |
|
| 137 |
| :--- | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: | :---: |
|
| 138 |
| Abdul et al. (2018), SMC5/6 complex | 9 | 7,073 | 0.990 | 0.752 | 100.0% | 0.00% | 0.874 | 561.0 s | 3.84 s | 171.3x |
|
| 139 |
| Nisson et al. (2025), CCDC137 (HIV Vpr) | 1 | 290 | 0.940 | 0.650 | 75.0% | 0.35% | 0.833 | 80.0 s | 0.21 s | 380.9x |
|
|
|
|
| 193 |
|
| 194 |
## Limitations and intended use
|
| 195 |
|
| 196 |
+
HyphAeon is a research tool for comparative molecular evolution, intended as a fast screen for
|
| 197 |
sites under episodic diversifying selection.
|
| 198 |
|
| 199 |
- It is trained to approximate HyPhy MEME, so it inherits that method's assumptions and is
|
|
|
|
| 215 |
A publication is in preparation. Until it appears, cite the repository:
|
| 216 |
|
| 217 |
```bibtex
|
| 218 |
+
@software{hyphaeon,
|
| 219 |
+
title = {HyphAeon: Ultra-Fast Neural Inference of Episodic Positive Selection in Molecular Sequences},
|
| 220 |
author = {Kosakovsky Pond, Sergei L. and collaborators},
|
| 221 |
+
url = {https://github.com/veg/HyphAeon},
|
| 222 |
year = {2026}
|
| 223 |
}
|
| 224 |
```
|