Raw report: aromatic_catabolism_network Feeds into: cofitness-network-architecture, condition-space-dimensionality, condition-specific-fitness, cross-species-fitness-transferability, evidence-triangulation-for-functional-annotation, experimental-prioritization-of-functional-dark-matter, fitness-module-detection-sensitivity, gene-essentiality, genomic-dispersal-functional-coupling, metabolic-model-gapfilling, multi-omics-integration, respiratory-capacity-and-nadh-load
Overview
This report defines a 51-gene support network for quinate catabolism in Acinetobacter baylyi ADP1 (acinetobacter-baylyi-adp1). The network surrounds the beta-ketoadipate-pathway and spans core aromatic degradation, Complex I / NADH dehydrogenase, iron acquisition, PQQ biosynthesis, and transcriptional regulation. The analysis combines growth phenotypes, flux-balance analysis (FBA, a constraint-based model that predicts metabolic fluxes by optimizing for growth), genomic organization, co-fitness correlations, and cross-species ortholog-transferred fitness data. 1
Key Findings
A 51-gene support network
The 51 quinate-specific genes organize into a coherent metabolic dependency network around the β-ketoadipate pathway. Co-fitness analysis assigns 44/51 genes (86%) to four functional subsystems — the core aromatic degradation pathway (8 genes), Complex I / NADH dehydrogenase (21 genes), iron acquisition (7 genes), and PQQ biosynthesis (2 genes) — plus 6 transcriptional regulators. The report’s wording counts the regulators within the 44-gene total, so the four subsystem counts alone do not make up the 44 assigned genes. A further 7 genes remain unassigned because of low co-fitness with all subsystems. The core pathway converts quinate through protocatechuate and β-ketoadipate to succinyl-CoA and acetyl-CoA, which enter the TCA cycle. 1
Co-fitness heatmap among 51 quinate-specific genes. From aromatic_catabolism_network.
The biochemical rationale has three parts. Quinate dehydrogenase (quiA) is a PQQ-dependent quinoprotein, which the report uses to explain the 2 PQQ-biosynthesis genes (pqqC, pqqD). Protocatechuate 3,4-dioxygenase (pcaGH) is a non-heme Fe²⁺-dependent ring-cleavage enzyme, which the report uses to explain the 7 iron-acquisition genes (siderophore biosynthesis, ExbD/TolR transport, ferrichrome receptor, TonB-dependent receptor). Succinyl-CoA and acetyl-CoA from β-ketoadipate are fully oxidized in the TCA cycle, generating NADH at multiple steps, and the report proposes that this requires complex-i for NADH reoxidation under high flux. 1
Complex I is the largest support subsystem and an FBA blind spot
Complex I, or NADH:ubiquinone oxidoreductase, accounts for 21/51 quinate-specific genes (41%) and is the largest support subsystem. FBA captures 1.76× higher Complex I flux on aromatic substrates than on the comparison condition, with fluxes of 0.55 versus 0.31, but predicts 0% essentiality. The model therefore detects increased respiratory demand without identifying Complex I as a bottleneck. 1
FBA-predicted flux heatmap for quinate-specific genes. From aromatic_catabolism_network.
A total of 30/51 quinate-specific genes have no FBA reaction mappings. These unmapped genes include cofactor-supply, iron-acquisition, regulatory, and newly identified Complex I-associated functions, showing that the model represents core metabolism more completely than the infrastructure required to support it. 1
The report interprets this discrepancy as a limitation of growth-optimizing linear programming: FBA can redistribute flux through alternative routes, whereas Complex I is a single multi-subunit complex whose disruption can eliminate complex function. This interpretation is supported by the observation that 10/13 Complex I operon subunits independently produce quinate-specific growth defects. 1
Genomic independence with metabolic coupling
The support subsystems occupy distinct chromosomal regions rather than a shared genomic neighborhood. The Complex I operon lies at 714–729 kb, the pca/qui pathway at 1,709–1,724 kb, PQQ biosynthesis at 2,461 kb, and iron-acquisition genes are scattered across 4 loci. No cross-category operons were identified, except within the aromatic pathway itself. 1
The Complex I operon contains 13 nuoA–N subunits on the same strand with <100 bp intergenic distances. The pca/qui region forms a 12-gene operon spanning pcaIJFBDCHG-quiABC plus transport genes; the report states that this operon prediction matches prior experimental transcriptional mapping of the pca/qui genes (Dal et al. 2005). Across the chromosome, 9 genomic clusters contain ≥2 quinate-specific genes, with mild overall clustering expressed as an observed/expected nearest-neighbor distance ratio of 0.89. 1
Co-fitness assigns previously unknown genes
Of 23 genes initially categorized as Other or Unknown, co-fitness assigns 16 to support subsystems with medium or high confidence. Two DUF-domain proteins, ACIAD3137 (UPF0234) and ACIAD2176 (DUF2280), correlate with Complex I genes at r > 0.98 and are candidate uncharacterized Complex I accessory factors. Within-category correlations are higher than between-category correlations, with mean r = 0.992 for Complex I and r = 0.961 for the aromatic pathway. 1
The co-fitness analysis also recovered pcaC, a 4-carboxymuconolactone decarboxylase that was initially miscategorized by keyword matching, and nominated two DUF proteins (ACIAD3137, ACIAD2176) as probable Complex I accessory factors. However, the analysis uses only 8 conditions and 8-dimensional growth vectors, limiting the resolution of gene-gene correlations; the 11 Complex I-associated assignments beyond the core nuo operon are based on phenotypic correlation and may represent indirect connections rather than physical association. 1
Complex I dependence tracks NADH-generating substrates
Ortholog-transferred fitness data from the kescience-fitnessbrowser contains 12,241 entries covering 2,005 genes and 13 conditions. Complex I orthologs have significantly worse fitness on aromatic conditions than on the comparison conditions, with mean fitness values of -1.35 versus -0.77 and Mann-Whitney p < 0.0001. 1
Per-condition analysis refines this result: the largest Complex I defects relative to background occur on acetate (-1.55) and succinate (-1.39), which are non-aromatic substrates that also generate high NADH flux through the TCA cycle. The report therefore reads the cross-species data as a high-TCA-flux dependency rather than a strictly aromatic one: Complex I is reported as dispensable on glucose and lactate but essential on substrates that generate high TCA-cycle flux (aromatics, acetate, succinate). The report attributes ADP1’s apparent quinate specificity of Complex I to an alternative NADH dehydrogenase, ndh 2, that likely compensates on simpler substrates. This compensation is an inferred hypothesis, not a result the report establishes. 1
The cross-species evidence is not definitive for ADP1 because the transferred data mixes organisms with different respiratory-chain architectures. Direct Complex I fitness measurements on aromatic substrates in a single organism would provide a stronger test of whether the dependency is caused by aromatic catabolism itself or by high NADH flux. 1
Relation to established pathway biology
The PQQ and iron dependencies are consistent with established biochemistry: PQQ-dependent quinate dehydrogenases require PQQ, and protocatechuate 3,4-dioxygenase is a non-heme Fe²⁺-dependent intradiol ring-cleavage enzyme. In ADP1, the growth-phenotype specificity of PQQ genes is consistent with prior transcriptomic evidence that 4/5 PQQ-biosynthesis genes are upregulated on quinate versus succinate. 1
The report presents the Complex I dependency as its novel finding. Bacterial respiratory chains can contain both proton-pumping Complex I (NDH-1) and non-pumping NDH-2; the proposed interpretation is that NDH-2 can substitute under lower NADH flux but cannot match Complex I capacity under high flux. 1
Caveats and Open Directions
The PQQ dependency is not exclusively aromatic: PQQ-biosynthesis genes also appear as glucose-specific in the adp1_deletion_phenotypes project, where they are associated with PQQ-dependent glucose dehydrogenase. 1
The report’s summary states that aromatic catabolism requires 51 genes, “over 8×” the core pathway, and that the support requirement exceeds the pathway itself “by 7:1”. Neither ratio is consistent with the report’s own counts of 51 total genes and 8 core-pathway genes. The intended denominators and the meaning of “support requirement” are unresolved, so these ratios should not be relied on. 1
The 8-condition co-fitness matrix provides approximately 5 independent dimensions, so additional conditions are needed to sharpen subsystem boundaries. The ortholog-transferred cross-species data is confounded by organism-specific respiratory architectures, and the non-core Complex I assignments rely on correlation rather than direct physical evidence. 1
The proposed NDH-2 compensation mechanism remains untested. The report proposes searching the ADP1 genome for NDH-2 and testing its deletion on quinate versus glucose; experimentally validating ACIAD3137 (UPF0234) and ACIAD2176 (DUF2280), which show r > 0.98 co-fitness with Complex I, by protein-protein interaction or co-purification studies; expanding the condition panel with benzoate, catechol, vanillate, iron limitation, and respiratory inhibitors; comparing Complex I retention across aromatic-degrading species using pangenome data; and adding PQQ biosynthesis, iron homeostasis, and respiratory-chain capacity constraints to the ADP1 FBA model. 1
Figures
- Quinate-specific genes by functional category (figures/support_network_categories.png)
- FBA-predicted flux heatmap for quinate-specific genes (figures/fba_flux_heatmap.png)
- Chromosome map of quinate-specific genes (figures/chromosome_map.png)
- Genomic clusters of quinate-specific genes (figures/gene_clusters.png)
- Co-fitness heatmap among 51 quinate-specific genes (figures/cofitness_heatmap.png)
- Unknown gene correlations with known subsystems (figures/unknown_assignments_heatmap.png)
- Within vs between category co-fitness (figures/cofitness_within_between.png)
- Cross-species fitness on aromatic vs non-aromatic conditions (figures/cross_species_fitness.png)
- Complex I fitness relative to background by condition (figures/complex_I_vs_background.png) 1
Cited Background Literature
The report’s reference list includes de Berardinis et al. (2008), “A complete collection of single-gene deletion mutants of Acinetobacter baylyi ADP1” (Molecular Systems Biology 4:174; PMID: 18319726); an article titled “Catabolite repression of aromatic compound degradation in Acinetobacter baylyi”; a work titled “The β-ketoadipate pathway”; and an article titled “Novel metabolic features in Acinetobacter baylyi ADP1 revealed by a multiomics approach”. These are cited background works, and their titles alone do not establish results within this report. 1
Slots Into
- cofitness-network-architecture — the 51-gene co-fitness network, subsystem assignments, r > 0.98 DUF candidates, and within-category correlations provide a condition-specific functional architecture. 1
- condition-specific-fitness — the quinate-specific phenotypes and cross-condition Complex I results refine the distinction between aromatic dependence and high-NADH-flux dependence. 1
- metabolic-model-gapfilling — the 0% predicted Complex I essentiality, 1.76× flux increase, and 30/51 unmapped genes identify cofactor, respiratory, and regulatory blind spots for FBA. 1
- gene-essentiality — the contrast between model-predicted 0% essentiality and observed defects for 10/13 Complex I subunits tests how complex-level constraints shape gene essentiality. 1
- multi-omics-integration — the report integrates fitness phenotypes, FBA, genomic organization, co-fitness, and prior transcriptomic evidence to explain aromatic-catabolism dependencies. 1
- respiratory-capacity-and-nadh-load — Complex I accounts for 21/51 quinate-specific genes. The largest cross-species Complex I defects, on acetate (-1.55) and succinate (-1.39), point to NADH load rather than aromatic chemistry, and NDH-2 compensation is proposed but untested. 1
- genomic-dispersal-functional-coupling — the support subsystems sit in distinct chromosomal regions with no cross-category operons, yet are metabolically coupled through PQQ, iron, and NADH demands. 1
- cross-species-fitness-transferability — ortholog-transferred Fitness Browser data (12,241 entries, 2,005 genes, 13 conditions) support a Complex I dependency but mix organisms with different respiratory-chain architectures. 1
- evidence-triangulation-for-functional-annotation — co-fitness assigned 16 of 23 Other/Unknown genes and corrected the keyword miscategorization of pcaC. 1
- experimental-prioritization-of-functional-dark-matter — ACIAD3137 (UPF0234) and ACIAD2176 (DUF2280), with r > 0.98 Complex I co-fitness, are candidate accessory factors awaiting experimental validation. 1
- condition-space-dimensionality — the 8-condition growth matrix provides only ~5 independent dimensions, limiting subsystem resolution. 1
- fitness-module-detection-sensitivity — co-fitness assignments derived from only 8 conditions have limited resolution and may capture indirect rather than physical associations. 1