celA1 Resolved · high auto-curated

H37Rv Rv0062 · MTBC0 - · 380 aa · 65552–66694 (+) · RefSeq YP_177689.1

Annotation: from legacy to revised

Legacy (H37Rv / Mycobrowser)cellulase CelA
MTBC0 PGAP re-annotation
Revised (this work)Cellulase CelA. Pfam: Glyco_hydro_6 (PF01341.24).

Auto-curated: this verdict and function were generated by rules from PGAP + Pfam + Foldseek and have not been hand-reviewed.

Annotated on the H37Rv protein: this gene has no 1:1 ancestral MTBC0 anchor (PE/PPE, paralogue, IS element, or otherwise unanchored CDS).

Curated reference (UniProt)

UniProt Q79G13 TrEMBL · unreviewed · Evidence at protein level
UniProt nameGlucanase
EC (curated) EC 3.2.1.-

Functional vocabulary (eggNOG-mapper, orthology transfer)

COG category G Carbohydrate transport and metabolism
Preferred namecelA1
eggNOG descriptionBelongs to the glycosyl hydrolase family 6
Orthologous groupCOG5297
EC number EC 3.2.1.4
KEGG orthology K01179
KEGG pathways map00500, map01100
CAZy family GH5, GH9

Orthology-based transfer (eggNOG 5.0.2, diamond). EC/KO/GO/CAZy are computed annotations, not manual curation; cross-check against the primary literature before treating a specific reaction as established.

Conservation & selection (intra-MTBC, 145 209 strains)

pN/pS 0.344 · purifying
Polymorphic sites (≥ 0.1% of strains) 18 synonymous, 17 missense, 0 nonsense, 1 frameshift
Disruption 1 distinct premature-stop/frameshift site(s); most common in 0.34% of strains (487) · clonal

pN/pS from segregating SNPs (singletons removed) normalised by possible sites. Low pN/pS = purifying selection (a strong signal that a "hypothetical" is a real, constrained gene). A high pN/pS is ambiguous: relaxed constraint or positive selection (drug resistance, antigenic variation) inflate it; e.g. rpoB/katG/pncA score high here for resistance, not loss of function. A clonal disruption (one allele over a clade) suggests lineage pseudogenisation; a convergent one (many independent alleles) is typical of resistance loss-of-function.

Domains (Pfam, hmmscan --cut_ga)

PfamAccessioni-EvalueResiduesDescription
Glyco_hydro_6PF01341.24 1.0e-6593–378 Glycosyl hydrolases family 6

Functional interaction network (STRING v12, guilt-by-association)

Closest characterised functional partner: celA2b (Rv1090, (MTV017.43), len: 151 aa. Probable celA2b,second part of cellulase (endoglucanase), similar to C-terminus of others e.g. O08468 cell), high confidence from genomic context alone (score 942 excluding text-mining).

PartnerProductScoreNo text-miningChannels (≥400)
Rv1090 celA2b Rv1090, (MTV017.43), len: 151 aa. Probable celA2b,second part of cellulase (endoglucanase), similar to C-terminus of others e.g. O08468 cell 972 942 ctx cooccurence:441 database:900 textmining:547
Rv1089A celA2a Rv1089A, len: 34 aa. Probable celA2a, first part of cellulase (endoglucanase), similar to N-terminus of others. This region is a possible MT 903 904 database:900
Rv0186 bglS beta-glucosidase BglS 981 902 database:900 textmining:819
Rv0061c hyp hypothetical protein 520 520 ctx neighborhood:520
Rv1084 hyp hypothetical protein 536 517 experimental:465
Rv0518 hyp hypothetical protein 541 512 experimental:405
Rv0452 transcriptional regulator 482 482 ctx cooccurence:479
Rv0063 oxidoreductase 468 468 ctx neighborhood:461
Rv0517 acyltransferase 446 447 ctx cooccurence:427
Rv0067c transcriptional regulator 443 443 ctx cooccurence:443
Rv2418c octT hyp hypothetical protein 474 441 experimental:405
Rv1075c hyp hypothetical protein 472 439 experimental:405
Rv0575c oxidoreductase 425 426 ctx cooccurence:424
Rv0048c membrane protein 401 402
Rv0315 beta-1,3-glucanase 915 386 textmining:868

STRING combines evidence channels (neighborhood, fusion, cooccurrence, coexpression, experimental, database, text-mining) into a 0–1000 score. The ctx badge marks edges carried by the genomic-context channels (conserved neighborhood, fusion, phylogenetic co-occurrence), which are independent of orthology and structure and the strongest signal for an unknown gene. The no text-mining column recomputes the score from data alone, so a link that does not depend on the literature is visible. Association is a function hypothesis, not proof: corroborate with the operon context and the primary literature before assigning a function.

Evidence

  • Annotation from H37Rv (no MTBC0 1:1 anchor; H37Rv protein used): cellulase CelA
  • Pfam (hmmscan --cut_ga): Glyco_hydro_6 PF01341.24 (E=1e-65)
  • (auto-curated by rules from PGAP + Pfam + Foldseek; not hand-reviewed)

Sources

  • Ancestral sequence & coordinates: Harrison LB et al. (2024), An imputed ancestral reference genome for the MTBC, doi:10.1101/2023.09.07.556366
  • Product annotation: NCBI PGAP on MTBC0; legacy from H37Rv NC_000962.3 (RefSeq YP_177689.1)
  • Domains: Pfam-A via hmmscan --cut_ga — Glyco_hydro_6 (PF01341.24)
  • Sequence-level signal: ESM Atlas (EvolutionaryScale × BioHub) — exploratory
  • Controlled vocabulary: eggNOG-mapper 2.1.12 (Cantalapiedra et al. 2021, doi:10.1093/molbev/msab293), eggNOG 5.0 DB (Huerta-Cepas et al. 2019) — OG COG5297
  • Curated reference: UniProt Q79G13 (TrEMBL, unreviewed; Evidence at protein level)
  • Intra-MTBC selection: pN/pS and disruption from SPDI variants of 145 209 MTBC strains (this work, local collection vs H37Rv NC_000962.3)
  • Interaction network: STRING v12.0 (Szklarczyk et al. 2023, doi:10.1093/nar/gkac1000), taxon 83332, CC-BY 4.0 — 29 functional partner(s); context anchor celA2b
  • Primary literature: none located yet; annotation rests on the domain/homology sources above.

Ancestral MTBC0 protein sequence

>H37Rv|Rv0062|celA1
MTRRTGQRWRGTLPGRRPWTRPAPATCRRHLAFVELRHYFARVMSSAIGSVARWIVPLLGVAAVASIGVIADPVRVVRAPALILVDAANPLAGKPFYVDPASAAMVAARNANPPNAELTSVANTPQSYWLDQAFPPATVGGTVARYTGAAQAAGAMPVLTLYGIPHRDCGSYASGGFATGTDYRGWIDAVASGLGSSPATIIVEPDALAMADCLSPDQRQERFDLVRYAVDTLTRDPAAAVYVDAGHSRWLSAEAMAARLNDVGVGRARGFSLNVSNFYTTDEEIGYGEAISGLTNGSHYVIDTSRNGAGPAPDAPLNWCNPSGRALGAPPTTATAGAHADAYLWIKRPGESDGTCGRGEPQAGRFVSQYAIDLAHNAGQ