uht-tooling 0.1.2__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- uht_tooling-0.1.2/PKG-INFO +271 -0
- uht_tooling-0.1.2/README.md +244 -0
- uht_tooling-0.1.2/pyproject.toml +39 -0
- uht_tooling-0.1.2/setup.cfg +4 -0
- uht_tooling-0.1.2/src/uht_tooling/__init__.py +10 -0
- uht_tooling-0.1.2/src/uht_tooling/cli.py +368 -0
- uht_tooling-0.1.2/src/uht_tooling/models/__init__.py +0 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/__init__.py +0 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/design_gibson.py +368 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/design_slim.py +402 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/gui.py +595 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/mut_rate.py +2480 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/mutation_caller.py +432 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/nextera_designer.py +199 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/profile_inserts.py +441 -0
- uht_tooling-0.1.2/src/uht_tooling/workflows/umi_hunter.py +412 -0
- uht_tooling-0.1.2/src/uht_tooling.egg-info/PKG-INFO +271 -0
- uht_tooling-0.1.2/src/uht_tooling.egg-info/SOURCES.txt +20 -0
- uht_tooling-0.1.2/src/uht_tooling.egg-info/dependency_links.txt +1 -0
- uht_tooling-0.1.2/src/uht_tooling.egg-info/entry_points.txt +2 -0
- uht_tooling-0.1.2/src/uht_tooling.egg-info/requires.txt +20 -0
- uht_tooling-0.1.2/src/uht_tooling.egg-info/top_level.txt +1 -0
|
@@ -0,0 +1,271 @@
|
|
|
1
|
+
Metadata-Version: 2.4
|
|
2
|
+
Name: uht-tooling
|
|
3
|
+
Version: 0.1.2
|
|
4
|
+
Summary: Tooling for ultra-high throughput screening workflows.
|
|
5
|
+
Author: Matt115A
|
|
6
|
+
License: MIT
|
|
7
|
+
Requires-Python: >=3.8
|
|
8
|
+
Description-Content-Type: text/markdown
|
|
9
|
+
Requires-Dist: biopython==1.85
|
|
10
|
+
Requires-Dist: fuzzywuzzy==0.18.0
|
|
11
|
+
Requires-Dist: matplotlib==3.10.7
|
|
12
|
+
Requires-Dist: pandas==2.3.3
|
|
13
|
+
Requires-Dist: python-Levenshtein==0.27.3
|
|
14
|
+
Requires-Dist: pyyaml==6.0.3
|
|
15
|
+
Requires-Dist: pysam==0.23.3
|
|
16
|
+
Requires-Dist: scipy==1.15.3
|
|
17
|
+
Requires-Dist: seaborn==0.13.2
|
|
18
|
+
Requires-Dist: tabulate==0.9.0
|
|
19
|
+
Requires-Dist: tqdm==4.67.1
|
|
20
|
+
Requires-Dist: typer==0.20.0
|
|
21
|
+
Provides-Extra: gui
|
|
22
|
+
Requires-Dist: gradio==5.49.1; extra == "gui"
|
|
23
|
+
Provides-Extra: dev
|
|
24
|
+
Requires-Dist: pytest==9.0.0; extra == "dev"
|
|
25
|
+
Requires-Dist: black==25.9.0; extra == "dev"
|
|
26
|
+
Requires-Dist: ruff==0.14.4; extra == "dev"
|
|
27
|
+
|
|
28
|
+
# uht-tooling
|
|
29
|
+
|
|
30
|
+
Automation helpers for ultra-high-throughput molecular biology workflows. The package ships both a Typer-based CLI and an optional Gradio GUI that wrap the same workflow code paths.
|
|
31
|
+
|
|
32
|
+
---
|
|
33
|
+
|
|
34
|
+
## Installation
|
|
35
|
+
|
|
36
|
+
### Quick install (recommended, easiest file maintainance)
|
|
37
|
+
```bash
|
|
38
|
+
python -m pip install "uht-tooling[gui]"
|
|
39
|
+
```
|
|
40
|
+
|
|
41
|
+
This installs the core workflows plus the optional GUI dependencies (Gradio, pandas). Omit the `[gui]` extras if you only need the CLI:
|
|
42
|
+
|
|
43
|
+
```bash
|
|
44
|
+
python -m pip install uht-tooling
|
|
45
|
+
```
|
|
46
|
+
|
|
47
|
+
### Development install
|
|
48
|
+
```bash
|
|
49
|
+
git clone https://github.com/Matt115A/uht-tooling.git
|
|
50
|
+
cd uht-tooling
|
|
51
|
+
python -m pip install -e ".[gui,dev]"
|
|
52
|
+
```
|
|
53
|
+
|
|
54
|
+
The editable install exposes the latest sources, while the `dev` extras add linting and test tooling.
|
|
55
|
+
|
|
56
|
+
---
|
|
57
|
+
|
|
58
|
+
## Directory layout
|
|
59
|
+
|
|
60
|
+
- Reference inputs live under `data/<workflow>/`.
|
|
61
|
+
- Outputs (CSV, FASTA, plots, logs) are written to `results/<workflow>/`.
|
|
62
|
+
- All workflows log to `results/<workflow>/run.log` for reproducibility and debugging.
|
|
63
|
+
|
|
64
|
+
---
|
|
65
|
+
|
|
66
|
+
## Command-line interface
|
|
67
|
+
|
|
68
|
+
The CLI is exposed as the `uht-tooling` executable. List the available commands:
|
|
69
|
+
|
|
70
|
+
```bash
|
|
71
|
+
uht-tooling --help
|
|
72
|
+
```
|
|
73
|
+
|
|
74
|
+
Each command mirrors a workflow module. Common entry points:
|
|
75
|
+
|
|
76
|
+
| Command | Purpose |
|
|
77
|
+
| --- | --- |
|
|
78
|
+
| `uht-tooling nextera-primers` | Generate Nextera XT primer pairs from a binding-region CSV. |
|
|
79
|
+
| `uht-tooling design-slim` | Design SLIM mutagenesis primers from FASTA/CSV inputs. |
|
|
80
|
+
| `uht-tooling design-gibson` | Produce Gibson mutagenesis primers and assembly plans. |
|
|
81
|
+
| `uht-tooling mutation-caller` | Summarise amino-acid substitutions from long-read FASTQ files. |
|
|
82
|
+
| `uht-tooling umi-hunter` | Cluster UMIs and call consensus alleles. |
|
|
83
|
+
| `uht-tooling ep-library-profile` | Measure mutation rates without UMIs. |
|
|
84
|
+
| `uht-tooling profile-inserts` | Extract inserts defined by probe pairs. |
|
|
85
|
+
|
|
86
|
+
Each command provides detailed help, including option descriptions and expected file formats:
|
|
87
|
+
|
|
88
|
+
```bash
|
|
89
|
+
uht-tooling mutation-caller --help
|
|
90
|
+
```
|
|
91
|
+
|
|
92
|
+
You can pass multiple FASTQ paths using repeated `--fastq` options or glob patterns. Optional `--log-path` flags redirect logs if you prefer a location outside the default results directory.
|
|
93
|
+
|
|
94
|
+
---
|
|
95
|
+
|
|
96
|
+
## Workflow reference
|
|
97
|
+
|
|
98
|
+
### Nextera XT primer design
|
|
99
|
+
|
|
100
|
+
1. Prepare `data/nextera_designer/nextera_designer.csv` with a `binding_region` column. Row 1 should contain the forward region, row 2 the reverse region, both in 5'→3' orientation.
|
|
101
|
+
2. Optional: supply a YAML overrides file for index lists/prefixes via `--config`.
|
|
102
|
+
3. Run:
|
|
103
|
+
```bash
|
|
104
|
+
uht-tooling nextera-primers \
|
|
105
|
+
--binding-csv data/nextera_designer/nextera_designer.csv \
|
|
106
|
+
--output-csv results/nextera_designer/nextera_xt_primers.csv
|
|
107
|
+
```
|
|
108
|
+
4. Primer CSVs will be written to `results/nextera_designer/`, accompanied by a log file.
|
|
109
|
+
|
|
110
|
+
The helper is preloaded with twelve i5 and twelve i7 indices, enabling up to 144 unique amplicons. Downstream lab workflow suggestions (qPCR monitoring, SPRIselect cleanup) remain unchanged from earlier releases.
|
|
111
|
+
|
|
112
|
+
#### Wet-lab workflow notes
|
|
113
|
+
|
|
114
|
+
- Perform the initial amplification with an i5/i7 primer pair and monitor a small aliquot by qPCR. Cap thermocycling early so you only generate ~10% of the theoretical yield—this minimizes amplification bias.
|
|
115
|
+
- Purify the product with SPRIselect beads at approximately a 0.65:1 bead:DNA volume ratio to remove residual primers and short fragments.
|
|
116
|
+
- Confirm primer removal using electrophoresis (e.g., BioAnalyzer DNA chip) before moving to sequencing prep.
|
|
117
|
+
|
|
118
|
+
### SLIM primer design
|
|
119
|
+
|
|
120
|
+
- Inputs:
|
|
121
|
+
- `data/design_slim/slim_template_gene.fasta`
|
|
122
|
+
- `data/design_slim/slim_context.fasta`
|
|
123
|
+
- `data/design_slim/slim_target_mutations.csv` (single `mutations` column)
|
|
124
|
+
- Run:
|
|
125
|
+
```bash
|
|
126
|
+
uht-tooling design-slim \
|
|
127
|
+
--gene-fasta data/design_slim/slim_template_gene.fasta \
|
|
128
|
+
--context-fasta data/design_slim/slim_context.fasta \
|
|
129
|
+
--mutations-csv data/design_slim/slim_target_mutations.csv \
|
|
130
|
+
--output-dir results/design_slim/
|
|
131
|
+
```
|
|
132
|
+
- Output: `results/design_slim/SLIM_primers.csv` plus logs.
|
|
133
|
+
|
|
134
|
+
Mutation nomenclature examples:
|
|
135
|
+
- `A123G` (substitution)
|
|
136
|
+
- `T241Del` (deletion)
|
|
137
|
+
- `T241TS` (insert Ser after Thr241)
|
|
138
|
+
- `L46GP` (replace Leu46 with Gly-Pro)
|
|
139
|
+
|
|
140
|
+
#### Experimental blueprint
|
|
141
|
+
|
|
142
|
+
- Hands-on time is approximately three hours (excluding protein purification), with mutant protein obtainable in roughly three days.
|
|
143
|
+
- Conduct two PCRs per mutant set: (A) long forward with short reverse and (B) long reverse with short forward.
|
|
144
|
+
- Combine 10 µL from each PCR with 10 µL H-buffer (150 mM Tris pH 8, 400 mM NaCl, 60 mM EDTA) for a 30 µL annealing reaction: 99 °C for 3 min, then two cycles of 65 °C for 5 min followed by 30 °C for 15 min, hold at 4 °C.
|
|
145
|
+
- Transform directly into NEB 5-alpha or BL21 (DE3) cells without additional cleanup. The protocol has been validated for simultaneous introduction of dozens of mutations.
|
|
146
|
+
|
|
147
|
+
### Gibson assembly primers
|
|
148
|
+
|
|
149
|
+
- Inputs mirror the SLIM workflow but use `data/design_gibson/`.
|
|
150
|
+
- Link sub-mutations with `+` to specify multi-mutation assemblies (e.g., `A123G+T150A`).
|
|
151
|
+
- Run:
|
|
152
|
+
```bash
|
|
153
|
+
uht-tooling design-gibson \
|
|
154
|
+
--gene-fasta data/design_gibson/gibson_template_gene.fasta \
|
|
155
|
+
--context-fasta data/design_gibson/gibson_context.fasta \
|
|
156
|
+
--mutations-csv data/design_gibson/gibson_target_mutations.csv \
|
|
157
|
+
--output-dir results/design_gibson/
|
|
158
|
+
```
|
|
159
|
+
- Outputs include primer sets and an assembly-plan CSV.
|
|
160
|
+
|
|
161
|
+
If mutations fall within overlapping primer windows, design sequential reactions to avoid excessive primer reuse.
|
|
162
|
+
|
|
163
|
+
### Mutation caller (no UMIs)
|
|
164
|
+
|
|
165
|
+
1. Supply:
|
|
166
|
+
- `data/mutation_caller/mutation_caller_template.fasta`
|
|
167
|
+
- `data/mutation_caller/mutation_caller.csv` with `gene_flanks` and `gene_min_max` columns (two rows each).
|
|
168
|
+
- One or more FASTQ files via `--fastq`.
|
|
169
|
+
2. Run:
|
|
170
|
+
```bash
|
|
171
|
+
uht-tooling mutation-caller \
|
|
172
|
+
--template-fasta data/mutation_caller/mutation_caller_template.fasta \
|
|
173
|
+
--flanks-csv data/mutation_caller/mutation_caller.csv \
|
|
174
|
+
--fastq data/mutation_caller/*.fastq.gz \
|
|
175
|
+
--output-dir results/mutation_caller/ \
|
|
176
|
+
--threshold 10
|
|
177
|
+
```
|
|
178
|
+
3. Outputs: per-sample subdirectories with substitution summaries, co-occurrence matrices, and logs.
|
|
179
|
+
|
|
180
|
+
### UMI Hunter
|
|
181
|
+
|
|
182
|
+
- Inputs: `data/umi_hunter/template.fasta`, `data/umi_hunter/umi_hunter.csv`, and FASTQ reads.
|
|
183
|
+
- Command:
|
|
184
|
+
```bash
|
|
185
|
+
uht-tooling umi-hunter \
|
|
186
|
+
--template-fasta data/umi_hunter/template.fasta \
|
|
187
|
+
--config-csv data/umi_hunter/umi_hunter.csv \
|
|
188
|
+
--fastq data/umi_hunter/*.fastq.gz \
|
|
189
|
+
--output-dir results/umi_hunter/
|
|
190
|
+
```
|
|
191
|
+
- Tunable parameters include `--umi-identity-threshold` and `--consensus-mutation-threshold`.
|
|
192
|
+
|
|
193
|
+
### Profile inserts
|
|
194
|
+
|
|
195
|
+
- Prepare `data/profile_inserts/sample_probes.csv` with `upstream` and `downstream` columns.
|
|
196
|
+
- Run:
|
|
197
|
+
```bash
|
|
198
|
+
uht-tooling profile-inserts \
|
|
199
|
+
--probes-csv data/profile_inserts/sample_probes.csv \
|
|
200
|
+
--fastq data/profile_inserts/*.fastq.gz \
|
|
201
|
+
--output-dir results/profile_inserts/
|
|
202
|
+
```
|
|
203
|
+
- Outputs: extracted insert FASTA files, QC plots, metrics, and logs. Adjust fuzzy matching strictness via `--min-ratio`.
|
|
204
|
+
|
|
205
|
+
### EP library profiler (no UMIs)
|
|
206
|
+
|
|
207
|
+
- Inputs:
|
|
208
|
+
- `data/ep-library-profile/region_of_interest.fasta`
|
|
209
|
+
- `data/ep-library-profile/plasmid.fasta`
|
|
210
|
+
- FASTQ inputs (`--fastq` accepts multiple files)
|
|
211
|
+
- Run:
|
|
212
|
+
```bash
|
|
213
|
+
uht-tooling ep-library-profile \
|
|
214
|
+
--region-fasta data/ep-library-profile/region_of_interest.fasta \
|
|
215
|
+
--plasmid-fasta data/ep-library-profile/plasmid.fasta \
|
|
216
|
+
--fastq data/ep-library-profile/*.fastq.gz \
|
|
217
|
+
--output-dir results/ep-library-profile/
|
|
218
|
+
```
|
|
219
|
+
- Output bundle includes per-sample directories and a master summary TSV.
|
|
220
|
+
|
|
221
|
+
---
|
|
222
|
+
|
|
223
|
+
## GUI quick start (optional)
|
|
224
|
+
|
|
225
|
+
The Gradio GUI wraps the same workflows with upload widgets and result previews. Launch it directly:
|
|
226
|
+
|
|
227
|
+
```bash
|
|
228
|
+
python -m uht_tooling.workflows.gui
|
|
229
|
+
```
|
|
230
|
+
|
|
231
|
+
Key points:
|
|
232
|
+
- The server binds to `http://127.0.0.1:7860` by default and falls back to an available port if 7860 is busy. Copy http://127.0.0.1:7860 into your browser.
|
|
233
|
+
- Temporary working directories are created under the system temp folder and cleaned automatically.
|
|
234
|
+
- Output archives (ZIP files) mirror the directory structure produced by the CLI.
|
|
235
|
+
|
|
236
|
+
### Tabs and capabilities
|
|
237
|
+
|
|
238
|
+
1. **Nextera XT** – forward/reverse primer inputs with CSV preview.
|
|
239
|
+
2. **SLIM** – template/context FASTA text areas plus mutation list.
|
|
240
|
+
3. **Gibson** – multi-mutation support using `+` syntax.
|
|
241
|
+
4. **Mutation Caller** – upload FASTQ, template FASTA, and configuration CSV.
|
|
242
|
+
5. **UMI Hunter** – long-read UMI clustering with configurable thresholds.
|
|
243
|
+
6. **Profile Inserts** – probe CSV and multiple FASTQ uploads.
|
|
244
|
+
7. **EP Library Profile** – FASTQ uploads plus plasmid and region FASTA inputs.
|
|
245
|
+
|
|
246
|
+
### Workflow tips
|
|
247
|
+
|
|
248
|
+
- For large FASTQ datasets, the CLI remains the most efficient option (especially for automation or batch processing).
|
|
249
|
+
- Use the command-line flag `--share` in `python -m uht_tooling.workflows.gui` if you need to expose the GUI outside localhost.
|
|
250
|
+
|
|
251
|
+
### Troubleshooting
|
|
252
|
+
|
|
253
|
+
- **Port already bound:** the launcher automatically selects the next free port and logs the chosen URL.
|
|
254
|
+
- **Missing dependency:** ensure you installed with `pip install "uht-tooling[gui]"`.
|
|
255
|
+
- **Stopping the server:** press `Ctrl+C` in the terminal session running the GUI.
|
|
256
|
+
|
|
257
|
+
---
|
|
258
|
+
|
|
259
|
+
## Logging
|
|
260
|
+
|
|
261
|
+
Every workflow configures logging to the destination output directory. Inspect `run.log` for command echoes, parameter choices, and any warnings produced during execution. When providing bug reports, include this log file along with input metadata to streamline triage.
|
|
262
|
+
|
|
263
|
+
---
|
|
264
|
+
|
|
265
|
+
## Roadmap
|
|
266
|
+
|
|
267
|
+
- Replace deprecated Biopython command-line wrappers with native subprocess implementations.
|
|
268
|
+
- Expand CLI coverage to any remaining legacy scripts that are still invoked via `make`.
|
|
269
|
+
- Add documentation for automation pipelines and integrate continuous integration tests.
|
|
270
|
+
|
|
271
|
+
Contributions in the form of bug reports, pull requests, or feature suggestions are welcome. File issues on GitHub with clear reproduction steps and sample data when possible.
|
|
@@ -0,0 +1,244 @@
|
|
|
1
|
+
# uht-tooling
|
|
2
|
+
|
|
3
|
+
Automation helpers for ultra-high-throughput molecular biology workflows. The package ships both a Typer-based CLI and an optional Gradio GUI that wrap the same workflow code paths.
|
|
4
|
+
|
|
5
|
+
---
|
|
6
|
+
|
|
7
|
+
## Installation
|
|
8
|
+
|
|
9
|
+
### Quick install (recommended, easiest file maintainance)
|
|
10
|
+
```bash
|
|
11
|
+
python -m pip install "uht-tooling[gui]"
|
|
12
|
+
```
|
|
13
|
+
|
|
14
|
+
This installs the core workflows plus the optional GUI dependencies (Gradio, pandas). Omit the `[gui]` extras if you only need the CLI:
|
|
15
|
+
|
|
16
|
+
```bash
|
|
17
|
+
python -m pip install uht-tooling
|
|
18
|
+
```
|
|
19
|
+
|
|
20
|
+
### Development install
|
|
21
|
+
```bash
|
|
22
|
+
git clone https://github.com/Matt115A/uht-tooling.git
|
|
23
|
+
cd uht-tooling
|
|
24
|
+
python -m pip install -e ".[gui,dev]"
|
|
25
|
+
```
|
|
26
|
+
|
|
27
|
+
The editable install exposes the latest sources, while the `dev` extras add linting and test tooling.
|
|
28
|
+
|
|
29
|
+
---
|
|
30
|
+
|
|
31
|
+
## Directory layout
|
|
32
|
+
|
|
33
|
+
- Reference inputs live under `data/<workflow>/`.
|
|
34
|
+
- Outputs (CSV, FASTA, plots, logs) are written to `results/<workflow>/`.
|
|
35
|
+
- All workflows log to `results/<workflow>/run.log` for reproducibility and debugging.
|
|
36
|
+
|
|
37
|
+
---
|
|
38
|
+
|
|
39
|
+
## Command-line interface
|
|
40
|
+
|
|
41
|
+
The CLI is exposed as the `uht-tooling` executable. List the available commands:
|
|
42
|
+
|
|
43
|
+
```bash
|
|
44
|
+
uht-tooling --help
|
|
45
|
+
```
|
|
46
|
+
|
|
47
|
+
Each command mirrors a workflow module. Common entry points:
|
|
48
|
+
|
|
49
|
+
| Command | Purpose |
|
|
50
|
+
| --- | --- |
|
|
51
|
+
| `uht-tooling nextera-primers` | Generate Nextera XT primer pairs from a binding-region CSV. |
|
|
52
|
+
| `uht-tooling design-slim` | Design SLIM mutagenesis primers from FASTA/CSV inputs. |
|
|
53
|
+
| `uht-tooling design-gibson` | Produce Gibson mutagenesis primers and assembly plans. |
|
|
54
|
+
| `uht-tooling mutation-caller` | Summarise amino-acid substitutions from long-read FASTQ files. |
|
|
55
|
+
| `uht-tooling umi-hunter` | Cluster UMIs and call consensus alleles. |
|
|
56
|
+
| `uht-tooling ep-library-profile` | Measure mutation rates without UMIs. |
|
|
57
|
+
| `uht-tooling profile-inserts` | Extract inserts defined by probe pairs. |
|
|
58
|
+
|
|
59
|
+
Each command provides detailed help, including option descriptions and expected file formats:
|
|
60
|
+
|
|
61
|
+
```bash
|
|
62
|
+
uht-tooling mutation-caller --help
|
|
63
|
+
```
|
|
64
|
+
|
|
65
|
+
You can pass multiple FASTQ paths using repeated `--fastq` options or glob patterns. Optional `--log-path` flags redirect logs if you prefer a location outside the default results directory.
|
|
66
|
+
|
|
67
|
+
---
|
|
68
|
+
|
|
69
|
+
## Workflow reference
|
|
70
|
+
|
|
71
|
+
### Nextera XT primer design
|
|
72
|
+
|
|
73
|
+
1. Prepare `data/nextera_designer/nextera_designer.csv` with a `binding_region` column. Row 1 should contain the forward region, row 2 the reverse region, both in 5'→3' orientation.
|
|
74
|
+
2. Optional: supply a YAML overrides file for index lists/prefixes via `--config`.
|
|
75
|
+
3. Run:
|
|
76
|
+
```bash
|
|
77
|
+
uht-tooling nextera-primers \
|
|
78
|
+
--binding-csv data/nextera_designer/nextera_designer.csv \
|
|
79
|
+
--output-csv results/nextera_designer/nextera_xt_primers.csv
|
|
80
|
+
```
|
|
81
|
+
4. Primer CSVs will be written to `results/nextera_designer/`, accompanied by a log file.
|
|
82
|
+
|
|
83
|
+
The helper is preloaded with twelve i5 and twelve i7 indices, enabling up to 144 unique amplicons. Downstream lab workflow suggestions (qPCR monitoring, SPRIselect cleanup) remain unchanged from earlier releases.
|
|
84
|
+
|
|
85
|
+
#### Wet-lab workflow notes
|
|
86
|
+
|
|
87
|
+
- Perform the initial amplification with an i5/i7 primer pair and monitor a small aliquot by qPCR. Cap thermocycling early so you only generate ~10% of the theoretical yield—this minimizes amplification bias.
|
|
88
|
+
- Purify the product with SPRIselect beads at approximately a 0.65:1 bead:DNA volume ratio to remove residual primers and short fragments.
|
|
89
|
+
- Confirm primer removal using electrophoresis (e.g., BioAnalyzer DNA chip) before moving to sequencing prep.
|
|
90
|
+
|
|
91
|
+
### SLIM primer design
|
|
92
|
+
|
|
93
|
+
- Inputs:
|
|
94
|
+
- `data/design_slim/slim_template_gene.fasta`
|
|
95
|
+
- `data/design_slim/slim_context.fasta`
|
|
96
|
+
- `data/design_slim/slim_target_mutations.csv` (single `mutations` column)
|
|
97
|
+
- Run:
|
|
98
|
+
```bash
|
|
99
|
+
uht-tooling design-slim \
|
|
100
|
+
--gene-fasta data/design_slim/slim_template_gene.fasta \
|
|
101
|
+
--context-fasta data/design_slim/slim_context.fasta \
|
|
102
|
+
--mutations-csv data/design_slim/slim_target_mutations.csv \
|
|
103
|
+
--output-dir results/design_slim/
|
|
104
|
+
```
|
|
105
|
+
- Output: `results/design_slim/SLIM_primers.csv` plus logs.
|
|
106
|
+
|
|
107
|
+
Mutation nomenclature examples:
|
|
108
|
+
- `A123G` (substitution)
|
|
109
|
+
- `T241Del` (deletion)
|
|
110
|
+
- `T241TS` (insert Ser after Thr241)
|
|
111
|
+
- `L46GP` (replace Leu46 with Gly-Pro)
|
|
112
|
+
|
|
113
|
+
#### Experimental blueprint
|
|
114
|
+
|
|
115
|
+
- Hands-on time is approximately three hours (excluding protein purification), with mutant protein obtainable in roughly three days.
|
|
116
|
+
- Conduct two PCRs per mutant set: (A) long forward with short reverse and (B) long reverse with short forward.
|
|
117
|
+
- Combine 10 µL from each PCR with 10 µL H-buffer (150 mM Tris pH 8, 400 mM NaCl, 60 mM EDTA) for a 30 µL annealing reaction: 99 °C for 3 min, then two cycles of 65 °C for 5 min followed by 30 °C for 15 min, hold at 4 °C.
|
|
118
|
+
- Transform directly into NEB 5-alpha or BL21 (DE3) cells without additional cleanup. The protocol has been validated for simultaneous introduction of dozens of mutations.
|
|
119
|
+
|
|
120
|
+
### Gibson assembly primers
|
|
121
|
+
|
|
122
|
+
- Inputs mirror the SLIM workflow but use `data/design_gibson/`.
|
|
123
|
+
- Link sub-mutations with `+` to specify multi-mutation assemblies (e.g., `A123G+T150A`).
|
|
124
|
+
- Run:
|
|
125
|
+
```bash
|
|
126
|
+
uht-tooling design-gibson \
|
|
127
|
+
--gene-fasta data/design_gibson/gibson_template_gene.fasta \
|
|
128
|
+
--context-fasta data/design_gibson/gibson_context.fasta \
|
|
129
|
+
--mutations-csv data/design_gibson/gibson_target_mutations.csv \
|
|
130
|
+
--output-dir results/design_gibson/
|
|
131
|
+
```
|
|
132
|
+
- Outputs include primer sets and an assembly-plan CSV.
|
|
133
|
+
|
|
134
|
+
If mutations fall within overlapping primer windows, design sequential reactions to avoid excessive primer reuse.
|
|
135
|
+
|
|
136
|
+
### Mutation caller (no UMIs)
|
|
137
|
+
|
|
138
|
+
1. Supply:
|
|
139
|
+
- `data/mutation_caller/mutation_caller_template.fasta`
|
|
140
|
+
- `data/mutation_caller/mutation_caller.csv` with `gene_flanks` and `gene_min_max` columns (two rows each).
|
|
141
|
+
- One or more FASTQ files via `--fastq`.
|
|
142
|
+
2. Run:
|
|
143
|
+
```bash
|
|
144
|
+
uht-tooling mutation-caller \
|
|
145
|
+
--template-fasta data/mutation_caller/mutation_caller_template.fasta \
|
|
146
|
+
--flanks-csv data/mutation_caller/mutation_caller.csv \
|
|
147
|
+
--fastq data/mutation_caller/*.fastq.gz \
|
|
148
|
+
--output-dir results/mutation_caller/ \
|
|
149
|
+
--threshold 10
|
|
150
|
+
```
|
|
151
|
+
3. Outputs: per-sample subdirectories with substitution summaries, co-occurrence matrices, and logs.
|
|
152
|
+
|
|
153
|
+
### UMI Hunter
|
|
154
|
+
|
|
155
|
+
- Inputs: `data/umi_hunter/template.fasta`, `data/umi_hunter/umi_hunter.csv`, and FASTQ reads.
|
|
156
|
+
- Command:
|
|
157
|
+
```bash
|
|
158
|
+
uht-tooling umi-hunter \
|
|
159
|
+
--template-fasta data/umi_hunter/template.fasta \
|
|
160
|
+
--config-csv data/umi_hunter/umi_hunter.csv \
|
|
161
|
+
--fastq data/umi_hunter/*.fastq.gz \
|
|
162
|
+
--output-dir results/umi_hunter/
|
|
163
|
+
```
|
|
164
|
+
- Tunable parameters include `--umi-identity-threshold` and `--consensus-mutation-threshold`.
|
|
165
|
+
|
|
166
|
+
### Profile inserts
|
|
167
|
+
|
|
168
|
+
- Prepare `data/profile_inserts/sample_probes.csv` with `upstream` and `downstream` columns.
|
|
169
|
+
- Run:
|
|
170
|
+
```bash
|
|
171
|
+
uht-tooling profile-inserts \
|
|
172
|
+
--probes-csv data/profile_inserts/sample_probes.csv \
|
|
173
|
+
--fastq data/profile_inserts/*.fastq.gz \
|
|
174
|
+
--output-dir results/profile_inserts/
|
|
175
|
+
```
|
|
176
|
+
- Outputs: extracted insert FASTA files, QC plots, metrics, and logs. Adjust fuzzy matching strictness via `--min-ratio`.
|
|
177
|
+
|
|
178
|
+
### EP library profiler (no UMIs)
|
|
179
|
+
|
|
180
|
+
- Inputs:
|
|
181
|
+
- `data/ep-library-profile/region_of_interest.fasta`
|
|
182
|
+
- `data/ep-library-profile/plasmid.fasta`
|
|
183
|
+
- FASTQ inputs (`--fastq` accepts multiple files)
|
|
184
|
+
- Run:
|
|
185
|
+
```bash
|
|
186
|
+
uht-tooling ep-library-profile \
|
|
187
|
+
--region-fasta data/ep-library-profile/region_of_interest.fasta \
|
|
188
|
+
--plasmid-fasta data/ep-library-profile/plasmid.fasta \
|
|
189
|
+
--fastq data/ep-library-profile/*.fastq.gz \
|
|
190
|
+
--output-dir results/ep-library-profile/
|
|
191
|
+
```
|
|
192
|
+
- Output bundle includes per-sample directories and a master summary TSV.
|
|
193
|
+
|
|
194
|
+
---
|
|
195
|
+
|
|
196
|
+
## GUI quick start (optional)
|
|
197
|
+
|
|
198
|
+
The Gradio GUI wraps the same workflows with upload widgets and result previews. Launch it directly:
|
|
199
|
+
|
|
200
|
+
```bash
|
|
201
|
+
python -m uht_tooling.workflows.gui
|
|
202
|
+
```
|
|
203
|
+
|
|
204
|
+
Key points:
|
|
205
|
+
- The server binds to `http://127.0.0.1:7860` by default and falls back to an available port if 7860 is busy. Copy http://127.0.0.1:7860 into your browser.
|
|
206
|
+
- Temporary working directories are created under the system temp folder and cleaned automatically.
|
|
207
|
+
- Output archives (ZIP files) mirror the directory structure produced by the CLI.
|
|
208
|
+
|
|
209
|
+
### Tabs and capabilities
|
|
210
|
+
|
|
211
|
+
1. **Nextera XT** – forward/reverse primer inputs with CSV preview.
|
|
212
|
+
2. **SLIM** – template/context FASTA text areas plus mutation list.
|
|
213
|
+
3. **Gibson** – multi-mutation support using `+` syntax.
|
|
214
|
+
4. **Mutation Caller** – upload FASTQ, template FASTA, and configuration CSV.
|
|
215
|
+
5. **UMI Hunter** – long-read UMI clustering with configurable thresholds.
|
|
216
|
+
6. **Profile Inserts** – probe CSV and multiple FASTQ uploads.
|
|
217
|
+
7. **EP Library Profile** – FASTQ uploads plus plasmid and region FASTA inputs.
|
|
218
|
+
|
|
219
|
+
### Workflow tips
|
|
220
|
+
|
|
221
|
+
- For large FASTQ datasets, the CLI remains the most efficient option (especially for automation or batch processing).
|
|
222
|
+
- Use the command-line flag `--share` in `python -m uht_tooling.workflows.gui` if you need to expose the GUI outside localhost.
|
|
223
|
+
|
|
224
|
+
### Troubleshooting
|
|
225
|
+
|
|
226
|
+
- **Port already bound:** the launcher automatically selects the next free port and logs the chosen URL.
|
|
227
|
+
- **Missing dependency:** ensure you installed with `pip install "uht-tooling[gui]"`.
|
|
228
|
+
- **Stopping the server:** press `Ctrl+C` in the terminal session running the GUI.
|
|
229
|
+
|
|
230
|
+
---
|
|
231
|
+
|
|
232
|
+
## Logging
|
|
233
|
+
|
|
234
|
+
Every workflow configures logging to the destination output directory. Inspect `run.log` for command echoes, parameter choices, and any warnings produced during execution. When providing bug reports, include this log file along with input metadata to streamline triage.
|
|
235
|
+
|
|
236
|
+
---
|
|
237
|
+
|
|
238
|
+
## Roadmap
|
|
239
|
+
|
|
240
|
+
- Replace deprecated Biopython command-line wrappers with native subprocess implementations.
|
|
241
|
+
- Expand CLI coverage to any remaining legacy scripts that are still invoked via `make`.
|
|
242
|
+
- Add documentation for automation pipelines and integrate continuous integration tests.
|
|
243
|
+
|
|
244
|
+
Contributions in the form of bug reports, pull requests, or feature suggestions are welcome. File issues on GitHub with clear reproduction steps and sample data when possible.
|
|
@@ -0,0 +1,39 @@
|
|
|
1
|
+
[build-system]
|
|
2
|
+
requires = ["setuptools>=68", "wheel"]
|
|
3
|
+
build-backend = "setuptools.build_meta"
|
|
4
|
+
|
|
5
|
+
[project]
|
|
6
|
+
name = "uht-tooling"
|
|
7
|
+
version = "0.1.2"
|
|
8
|
+
description = "Tooling for ultra-high throughput screening workflows."
|
|
9
|
+
readme = "README.md"
|
|
10
|
+
requires-python = ">=3.8"
|
|
11
|
+
license = { text = "MIT" }
|
|
12
|
+
authors = [{ name = "Matt115A" }]
|
|
13
|
+
dependencies = [
|
|
14
|
+
"biopython==1.85",
|
|
15
|
+
"fuzzywuzzy==0.18.0",
|
|
16
|
+
"matplotlib==3.10.7",
|
|
17
|
+
"pandas==2.3.3",
|
|
18
|
+
"python-Levenshtein==0.27.3",
|
|
19
|
+
"pyyaml==6.0.3",
|
|
20
|
+
"pysam==0.23.3",
|
|
21
|
+
"scipy==1.15.3",
|
|
22
|
+
"seaborn==0.13.2",
|
|
23
|
+
"tabulate==0.9.0",
|
|
24
|
+
"tqdm==4.67.1",
|
|
25
|
+
"typer==0.20.0",
|
|
26
|
+
]
|
|
27
|
+
|
|
28
|
+
[project.optional-dependencies]
|
|
29
|
+
gui = ["gradio==5.49.1"]
|
|
30
|
+
dev = ["pytest==9.0.0", "black==25.9.0", "ruff==0.14.4"]
|
|
31
|
+
|
|
32
|
+
[project.scripts]
|
|
33
|
+
uht-tooling = "uht_tooling.cli:main"
|
|
34
|
+
|
|
35
|
+
[tool.setuptools]
|
|
36
|
+
package-dir = { "" = "src" }
|
|
37
|
+
|
|
38
|
+
[tool.setuptools.packages.find]
|
|
39
|
+
where = ["src"]
|