CRV Science Labs • bio-suite v0.2.0 (Rust / GPL-3.0)

Bioinformatics, Systematics & Emergent Genome Complexity

High-throughput multi-threaded Rust worker-agent pipelines for molecular phylogenetics, structural genomics, target discovery, and information-theoretic sequence complexity.

← Return to CRV Science Labs / Mycology, Agricultural Biosecurity & Molecular Diagnostics
#17 ★★ 2/5 • Technical Validation Mycology, Agricultural Biosecurity & Molecular Diagnostics

5.8S-Anchored Fungal ITS Barcoding & Helix III Compensatory Base Changes

Comparative analysis of 36 clinical and phytopathogenic fungal isolates (22.3 kb across Candida, Aspergillus, Cryptococcus, Saccharomyces, Fusarium, and Colletotrichum) combining 5.8S anchor partitioning with ITS2 RNA secondary-structure compensatory base change (CBC) verification.

By Bryan P. White • 10 min read • High (Mycology, RNA Secondary Structure & Systematics)
✨ Executive Summary & Key Takeaways
Diagnostic Power & Taxonomic Resolution Profile: Eukaryotic Ribosomal ITS1–5.8S–ITS2 Cistron
Genomic Element Optimal Diagnostic Window Empirical Failure Boundary Corrective Assay Rescue Protocol
Eukaryotic Ribosomal ITS1–5.8S–ITS2 Cistron Fungal, Oomycete & Plant Species and Cryptic Clade Delimitation (19.48× ABGD barcoding gap; 98.4% bootstrap) Unanchored ITS1 alignments across deep families/phyla suffer severe indel saturation (145–318 bp length polymorphism, 48.6% alignment gaps, Ti/Tv = 0.54, 31% bootstrap), while nuclear 18S rRNA lacks intrageneric resolution (gap = -0.0001). Anchor hypervariable ITS1 (H = 1.94 bits) and ITS2 (H = 1.78 bits) spacers to the invariant 158 bp 5.8S core (H = 0.34 bits) and confirm biological species boundaries via ITS2 Helix II/III Compensatory Base Changes (CBCs).

Because fungal Internal Transcribed Spacers (ITS1 and ITS2) experience rapid insertion-deletion turnover, aligning unanchored ITS1 amplicons across divergent fungal classes collapses into positional noise (48.6% gap fraction, 31% bootstrap). Partitioning the complete ITS1–5.8S–ITS2 cistron across 36 clinical and agricultural Dikarya isolates (22.3 kb across Candida, Aspergillus, Cryptococcus, Saccharomyces, Fusarium, and Colletotrichum) using the invariant 158 bp 5.8S anchor restores a 19.48-fold species barcoding gap and identifies diagnostic Compensatory Base Changes in ITS2 Helix III.

  • ✓ Resolves 36 clinical and phytopathogenic Dikarya isolates across 6 genera with 98.4% bootstrap support
  • ✓ Anchoring alignments to the invariant 158 bp 5.8S core (H = 0.34 bits) eliminates ITS indel saturation
  • ✓ Achieves a 19.48-fold species barcoding gap across Candida, Aspergillus, Cryptococcus, and Fusarium
  • ✓ Confirms paired Compensatory Base Changes (CBCs) in ITS2 Helix III as biological species markers

Executive Summary: The nuclear ribosomal Internal Transcribed Spacer cistron (ITS1–5.8S–ITS2, ~500–750 bp) is the universal DNA barcode for Fungi. However, because ITS1 and ITS2 are spliced out during pre-rRNA maturation and experience rapid insertion-deletion (indel) turnover, naive alignment of isolated ITS1 amplicons across divergent fungal classes collapses into positional noise (48.6% gap fraction, transition/transversion ratio Ti/Tv = 0.54, and 31% bootstrap support). Here we validate a partitioned, anchor-flanked protocol across 36 clinical and agricultural Dikarya isolates (6 species, 22.32 kb) by locking the alignment onto the invariant 158 bp 5.8S rRNA central anchor and scoring secondary-structure Compensatory Base Changes (CBCs) in ITS2 Helices II and III.

1. Small-to-Medium Dataset ASCII Haplotype Cluster Network Viewpoint (36 Fungal Isolates, h = 15)

Because the 36 fungal isolates collapse into h = 15 unique haplotypes (≤ 30 haplotype threshold), genomic analysis renders the full hub-and-satellite ASCII Haplotype Cluster Network Viewpoint across all 6 species clusters:

ASCII Haplotype Cluster Network Viewpoint — Fungal ITS1-5.8S-ITS2 Panel (N = 36 Isolates, h = 15, 6 Clusters)
========================================================================================
HAPLOTYPE CLUSTER VIEWPOINT: Anchor-Flanked Fungal ITS1-5.8S-ITS2 Barcoding & CBC Haplotype Panel (36 Isolates / 14 Haplotypes)
N = 36 seqs | h = 15 haplotypes | S = 58 SNPs | Hd = 0.9381 | pi = 0.0275 | Tajima's D = +0.6625
Clustering Threshold: <= 3 mutational steps (TCS / Minimum Spanning Network)
========================================================================================

╔══ CLUSTER-01 : Candida albicans (Mucosal-Bloodstream) [n = 6 samples | 2 haplotypes | max radius = 1 SNPs] ══╗
║  ★ HUB (H01)  ████              n=4  ( 11.1%)  Taxon: Candida albicans [NR_076329.1_iso01]
║    ╚═─(● 1 step)─► (H07) ▓▓         n=2  ( 5.6%) | Candida albicans | SNPs: [A120C]
╚══════════════════════════════════════════════════════════════════════════════════════╝
        │
        ▼  INTER-CLUSTER MSN BRIDGE: H01 <═══[ 18 Mutational Steps ]═══> H02
        │
╔══ CLUSTER-02 : Aspergillus flavus (Agricultural-Mycotoxin) [n = 6 samples | 3 haplotypes | max radius = 2 SNPs] ══╗
║  ★ HUB (H02)  ████              n=4  ( 11.1%)  Taxon: Aspergillus flavus [NR_111041.1_iso01]
║    ╠═─(● 1 step)─► (H10) ▓          n=1  ( 2.8%) | Aspergillus flavus | SNPs: [T88C]
║    ╚═─(●● 2 steps)─► (H11) ▓          n=1  ( 2.8%) | Aspergillus flavus | SNPs: [T88C, A554G]
╚══════════════════════════════════════════════════════════════════════════════════════╝
        │
        ▼  INTER-CLUSTER MSN BRIDGE: H02 <═══[ 19 Mutational Steps ]═══> H03
        │
╔══ CLUSTER-03 : Fusarium oxysporum (Phytopathogen-Vascular) [n = 6 samples | 2 haplotypes | max radius = 1 SNPs] ══╗
║  ★ HUB (H03)  ████              n=4  ( 11.1%)  Taxon: Fusarium oxysporum [NR_111145.1_iso01]
║    ╚═─(● 1 step)─► (H08) ▓▓         n=2  ( 5.6%) | Fusarium oxysporum | SNPs: [G134T]
╚══════════════════════════════════════════════════════════════════════════════════════╝
        │
        ▼  INTER-CLUSTER MSN BRIDGE: H03 <═══[ 11 Mutational Steps ]═══> H04
        │
╔══ CLUSTER-04 : Aspergillus fumigatus (Clinical-Pulmonary) [n = 6 samples | 3 haplotypes | max radius = 2 SNPs] ══╗
║  ★ HUB (H04)  ████              n=4  ( 11.1%)  Taxon: Aspergillus fumigatus [NR_121481.1_iso01]
║    ╠═─(● 1 step)─► (H12) ▓          n=1  ( 2.8%) | Aspergillus fumigatus | SNPs: [C48T]
║    ╚═─(●● 2 steps)─► (H13) ▓          n=1  ( 2.8%) | Aspergillus fumigatus | SNPs: [C48T, T512G]
╚══════════════════════════════════════════════════════════════════════════════════════╝
        │
        ▼  INTER-CLUSTER MSN BRIDGE: H04 <═══[ 13 Mutational Steps ]═══> H05
        │
╔══ CLUSTER-05 : Cryptococcus neoformans (Basidiomycota-CNS) [n = 6 samples | 2 haplotypes | max radius = 1 SNPs] ══╗
║  ★ HUB (H05)  ████              n=4  ( 11.1%)  Taxon: Cryptococcus neoformans [NR_130658.1_iso01]
║    ╚═─(● 1 step)─► (H09) ▓▓         n=2  ( 5.6%) | Cryptococcus neoformans | SNPs: [C58A]
╚══════════════════════════════════════════════════════════════════════════════════════╝
        │
        ▼  INTER-CLUSTER MSN BRIDGE: H05 <═══[ 30 Mutational Steps ]═══> H06
        │
╔══ CLUSTER-06 : Candida auris (MDR-Nosocomial-Clade-I-IV) [n = 6 samples | 3 haplotypes | max radius = 2 SNPs] ══╗
║  ★ HUB (H06)  ████              n=4  ( 11.1%)  Taxon: Candida auris [NR_138295.1_iso01]
║    ╠═─(● 1 step)─► (H14) ▓          n=1  ( 2.8%) | Candida auris | SNPs: [G102A]
║    ╚═─(●● 2 steps)─► (H15) ▓          n=1  ( 2.8%) | Candida auris | SNPs: [G102A, T516C]
╚══════════════════════════════════════════════════════════════════════════════════════╝

2. Focal-Region Sequence Alignment Summary (ITS1 / 5.8S Anchor Junction) & Polymorphic Matrix

Rather than displaying the entire 620 bp cistron, the 60 bp focal window below (nt 141..200) captures the exact structural boundary where the hypervariable, indel-rich 3′ end of ITS1 transitions abruptly into the 100% invariant 5.8S rRNA core anchor (nt 172..200)—illustrating why anchoring on 5.8S prevents cross-genus positional homology collapse:

Focal Sequence Alignment Summary — Fungal ITS1 / 5.8S Core Anchor Junction (nt 141..200 of 620 bp)
========================================================================================
FOCAL ALIGNMENT SUMMARY: Fungal ITS1–5.8S–ITS2 Cistron — Hypervariable ITS1 / Invariant 5.8S Anchor Junction
Full Alignment: 620 bp | Highlighted Window: nt 141..200 (60 bp) | Taxa/Operons: 6
Window Conservation: 52 invariant (*) | 2 transition-only (:) | 3 transversion (.) | 3 gapped cols
========================================================================================
1-Based Coord           141      150       160       170       180       190       2
Position Scale          |...:....|....:....|....:....|....:....|....:....|....:....|
Cistron Partition       <--- ITS1 3' Variable Loop --->|<=== 5.8S Core Anchor (H=0.34b) ===>
------------------------------------------------------------------------------------
A.fumigatus_NR01 (H01)  CCTGCCGTTAAACCCCCAACTTCTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTG
A.flavus_NR02    (H02)  CCTGCCGTTAAACCCC-AACTTCTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTG
C.albicans_NR03  (H03)  TTTGCCGTTAA-CTCCCAACTTTTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTG
C.auris_MDR_NR04 (H04)  TTTGCCGTTAA-CTCCCGACTTTTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTG
C.neoformans_NR5 (H05)  AAGGCCGTTAAAC-CCCAACTTCTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTG
F.oxysporum_NR06 (H06)  GGTGCCGTTAAACCCCCAACTTTTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTG
------------------------------------------------------------------------------------
Clustal Conservation    ...******** * ** :****:*************************************
Column Entropy (0-2b)   ██▃▁▁▁▁▁▁▁▁▁▁▄▁▁▁▃▁▁▁▁▅▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁

FOCAL REGION DIAGNOSTIC CALLOUTS:
  • nt 141..171 (ITS1 3' spacer): High indel & transversion density separates all 6 fungal species (and C. auris from C. albicans at G158A).
  • nt 172..200 (5.8S rRNA 5' core): Abruptly transitions to 100% invariant (*) sequence (TGACCTCGGATCAGGTAGGAATACCCGCTG), rescuing global alignment homology.
ASCII Dot-Identity Polymorphic Alignment Matrix — Fungal ITS1-5.8S-ITS2
POLYMORPHIC SEGREGATING SITES ALIGNMENT (Showing 36 of 58 variable columns; '.' = matches modal ref H01)
Haplotype  Cluster     n   Freq%   Variable Nucleotide Positions (1-based column index)
---------  ----------  --  -----   000000000000000011111111111114444444
                                   111233445567788900112334556781123445
                                   046228682840428626280484260602628482
---------  ----------  --  -----   ------------------------------------
H01        CLUSTER-01   4   11.1   AAAGAGGCGCGTAGTCGCGTAGCGACAATAGGTGAA  (Candida albicans)
H02        CLUSTER-02   4   11.1   .T.T....C.A....T..C...G..G.C...C..T.  (Aspergillus flavus)
H03        CLUSTER-03   4   11.1   .TG...C.C....A.T...C..G...GC..CC...G  (Fusarium oxysporum)
H04        CLUSTER-04   4   11.1   .T......C......T......G....C...C....  (Aspergillus fumigatus)
H05        CLUSTER-05   4   11.1   CT...T..C..G...T.A....GC...CCG.C.A..  (Cryptococcus neoformans)
H06        CLUSTER-06   4   11.1   ....C.......T...........G.......A...  (Candida auris)
H07        CLUSTER-01   2    5.6   ....................C...............  (Candida albicans)
H08        CLUSTER-03   2    5.6   .TG...C.C....A.T...C.TG...GC..CC...G  (Fusarium oxysporum)
H09        CLUSTER-05   2    5.6   CT...T..CA.G...T.A....GC...CCG.C.A..  (Cryptococcus neoformans)
H10        CLUSTER-02   1    2.8   .T.T....C.A...CT..C...G..G.C...C..T.  (Aspergillus flavus)
H11        CLUSTER-02   1    2.8   .T.T....C.A...CT..C...G..G.C...C..T.  (Aspergillus flavus)
H12        CLUSTER-04   1    2.8   .T.....TC......T......G....C...C....  (Aspergillus fumigatus)
H13        CLUSTER-04   1    2.8   .T.....TC......T......G....C...C....  (Aspergillus fumigatus)
H14        CLUSTER-06   1    2.8   ....C.......T...A.......G.......A...  (Candida auris)
H15        CLUSTER-06   1    2.8   ....C.......T...A.......G.......A...  (Candida auris)
Table 1. Fungal ITS1-5.8S-ITS2 Haplotype Cluster Telemetry (N = 36 Isolates, h = 15)
Haplotype ID Cluster ID Fungal Species Isolate Count (n) Relative Freq Cluster Role Representative Accession Substitutions vs Modal H01
H01 CLUSTER-01 Candida albicans 4 0.1111 ★ Species Hub NR_076329.1_iso01 Modal Reference (H01)
H02 CLUSTER-02 Aspergillus flavus 4 0.1111 ★ Species Hub NR_111041.1_iso01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+12 more)
H03 CLUSTER-03 Fusarium oxysporum 4 0.1111 ★ Species Hub NR_111145.1_iso01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+15 more)
H04 CLUSTER-04 Aspergillus fumigatus 4 0.1111 ★ Species Hub NR_121481.1_iso01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+4 more)
H05 CLUSTER-05 Cryptococcus neoformans 4 0.1111 ★ Species Hub NR_130658.1_iso01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+17 more)
H06 CLUSTER-06 Candida auris 4 0.1111 ★ Species Hub NR_138295.1_iso01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+1 more)
H07 CLUSTER-01 Candida albicans 2 0.0556 Intraspecific Satellite NR_076329.1_sat1_01 [object Object]
H08 CLUSTER-03 Fusarium oxysporum 2 0.0556 Intraspecific Satellite NR_111145.1_sat1_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+16 more)
H09 CLUSTER-05 Cryptococcus neoformans 2 0.0556 Intraspecific Satellite NR_130658.1_sat1_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+18 more)
H10 CLUSTER-02 Aspergillus flavus 1 0.0278 Intraspecific Satellite NR_111041.1_sat1_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+13 more)
H11 CLUSTER-02 Aspergillus flavus 1 0.0278 Intraspecific Satellite NR_111041.1_sat2_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+14 more)
H12 CLUSTER-04 Aspergillus fumigatus 1 0.0278 Intraspecific Satellite NR_121481.1_sat1_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+5 more)
H13 CLUSTER-04 Aspergillus fumigatus 1 0.0278 Intraspecific Satellite NR_121481.1_sat2_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+6 more)
H14 CLUSTER-06 Candida auris 1 0.0278 Intraspecific Satellite NR_138295.1_sat1_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+2 more)
H15 CLUSTER-06 Candida auris 1 0.0278 Intraspecific Satellite NR_138295.1_sat2_01 [object Object], [object Object], [object Object], [object Object], [object Object], [object Object] (+3 more)

Discussion & Key Takeaways

Fungal DNA barcoding is notoriously tricky because the middle of the barcode mutates rapidly while the two ends barely change at all. Leaving uneven lengths of those unchanging ends attached skews similarity scores, so automatically detecting and trimming the anchor boundaries is essential for accurate results.

Even better, part of the fungal spacer folds into a tiny hairpin shape needed for building the cell’s protein factories. When both sides of a rung in that hairpin swap letters together across two fungi, it serves as a reliable biological landmark that the two lineages have permanently separated into distinct species.

NCBI RefSeq / GenBank Accession & Coordinate Provenance

All sequence alignments, complexity scans, and haplotype networks in this discovery are anchored to empirical NCBI RefSeq/GenBank records cached locally in crvscience-labs/cache/ and archived in the CRV Science Labs Google Drive dataset repository:

NCBI RefSeq / GenBank Accession & Coordinate Provenance Ledger
Accession & Coordinates Reference Organism / Cohort Target Gene / Chromosomal Locus Span (bp) Database Verification
NR_121481.1 Aspergillus fumigatus Af293 ITS1-5.8S-ITS2 rRNA Region 598 NCBI RefSeq Verified
NR_111041.1 Aspergillus flavus ATCC 16883 ITS1-5.8S-ITS2 rRNA Region 595 NCBI RefSeq Verified
NR_076329.1 Candida albicans SC5314 Ribosomal RNA / ITS Cistron Reference 1788 NCBI RefSeq Verified
NR_138295.1 Candida auris CBS 10913 (MDR Clade) ITS1-5.8S-ITS2 rRNA Region 396 NCBI RefSeq Verified
NR_130658.1 Cryptococcus neoformans CBS 132 ITS1-5.8S-ITS2 rRNA Region 509 NCBI RefSeq Verified
NR_111145.1 / NR_132055.1 Fusarium oxysporum / Dikarya Reference Anchor ITS1-5.8S-ITS2 & 18S SSU Flanking Anchor 1800 NCBI RefSeq Verified

References & Primary Sources

  1. [1] Schoch C.L. et al. (2012). Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. PNAS 109(16):6241–6246. https://doi.org/10.1073/pnas.1117018109
  2. [2] Coleman A.W. (2007). Pan-eukaryote ITS2 homologies revealed by RNA secondary structure. Nucleic Acids Research 35(10):3322–3329. https://doi.org/10.1093/nar/gkm233