Deep Research Report: GO:0005201 as Core Function of A0A2G9RZF1

Executive Judgment

Verdict: Over-annotated (Refuted)

The hypothesis that extracellular matrix structural constituent (GO:0005201) is a core function of A0A2G9RZF1 is refuted by three convergent and independent lines of evidence:

1. Semantic mismatch: GO:0005201 describes proteins that contribute to the structural integrity of the ECM (collagens, elastin, fibrillin) — not CUB-domain interaction modules that mediate protein-protein recognition. CUB domains are 110-residue β-sandwich folds functioning as modular binding surfaces (PMID: 21954942), fundamentally distinct from structural ECM scaffolds.

2. FunFam evidence chain failure: The PCPE-1 FunFam classification that motivated this annotation is doubly flawed. PCPE-1's experimentally validated function is peptidase activator activity (GO:0016504, IDA), not ECM structural support. Its single GO:0005201 annotation derives from a proteomics cataloging study (RCA evidence from PMID: 20551380), and individual CUB domains cannot perform PCPE-1's enhancing function — which requires cooperative CUB1+CUB2 binding with >1,000-fold higher affinity than single CUB domains (PMID: 19801683).

3. Gene prediction artifact: A0A2G9RZF1 is almost certainly a gene prediction fragment from a poorly assembled genome. Scaffold KV928989 encodes three consecutive ORFs whose combined domain architecture (CUB–EGF–CUB–CUB–EGF–CUB) matches the C-terminal region of a BMP-1/tolloid-like metalloprotease. A chromosome-level assembly published in 2024 (GCF_042186555.1) encodes complete BMP-1 (1,020 aa) and TLL1 (1,005 aa) — full-length tolloid-family metalloproteases that the 156-aa fragment represents a portion of.

The molecular function of A0A2G9RZF1 should be left unassigned. No GO annotation should be applied to this fragment. If the complete gene product were annotated, its molecular function would be metalloendopeptidase activity — not ECM structural support.

Summary

A0A2G9RZF1 is a 156-amino acid unreviewed (TrEMBL) protein from Aquarana catesbeiana (American bullfrog), predicted at protein existence level 4 (predicted) from whole-genome shotgun sequencing. Its sole recognizable feature is a single CUB domain (residues 31–147). The seed hypothesis proposed that this protein functions as an extracellular matrix structural constituent (GO:0005201), based on (a) the general extracellular localization of CUB-domain proteins, and (b) FunFam classification to PCPE-1 (Procollagen C-endopeptidase enhancer 1), a protein detected in ECM proteomics datasets.

This investigation systematically evaluated three critical questions: (1) Is GO:0005201 the correct term for what CUB domains do? (2) Does the PCPE-1 FunFam classification justify GO:0005201? (3) Is A0A2G9RZF1 a real, complete gene product? The answer to all three is no. CUB domains are protein-protein interaction modules, not structural ECM components. PCPE-1 is experimentally a peptidase activator, and its single GO:0005201 annotation derives from a low-confidence computational analysis. Most importantly, A0A2G9RZF1 appears to be a gene prediction fragment from a fragmented genome assembly, representing one CUB domain from a much larger (~1,000 aa) tolloid-family metalloprotease.

These findings were confirmed across three iterations of investigation, incorporating analysis of 32 primary literature papers, genomic context examination, and cross-assembly validation. The conclusion is robust: GO:0005201 is inappropriate for this protein at every level of analysis.

Key Findings

Finding 1: GO:0005201 Is Semantically Inappropriate for CUB-Domain Proteins

GO:0005201 (extracellular matrix structural constituent) is defined as "The action of a molecule that contributes to the structural integrity of the extracellular matrix." Examination of QuickGO annotations reveals this term is overwhelmingly applied to bona fide structural ECM proteins: collagens (COL1A1, COL4A1-6, COL5A2), elastin (ELN), fibrillin, and similar molecules that physically constitute the ECM scaffold. These are large, repetitive structural proteins whose presence is necessary for the mechanical and organizational properties of the matrix.

CUB domains, by contrast, are 110-residue protein motifs that adopt a β-sandwich fold and mediate protein-protein interactions. As established in the authoritative review by Gaboriaud et al. (PMID: 21954942): CUB domains are "110-residue protein motifs exhibiting a β-sandwich fold and mediating protein-protein interactions in various extracellular proteins." They function as modular binding surfaces — recognition and interaction elements — not as structural building blocks of the ECM. Annotating a single-CUB-domain protein with GO:0005201 conflates a protein-protein interaction module with a structural matrix component, which are fundamentally different molecular functions.

Finding 2: PCPE-1's Validated Function Is Peptidase Activator Activity, Not ECM Structural Support

The FunFam classification that linked A0A2G9RZF1 to PCPE-1 was the primary justification for considering GO:0005201. However, examining PCPE-1's actual experimental annotations reveals a critical mismatch. Human PCPE-1 (UniProt Q15113) has been experimentally characterized with the following GO molecular functions supported by IDA (Inferred from Direct Assay) evidence from PMID: 12393877:

• Collagen binding (GO:0005518)
• Heparin binding (GO:0008201)
• Peptidase activator activity (GO:0016504)

PCPE-1's single GO:0005201 annotation is supported only by RCA (Reviewed Computational Analysis) evidence from a proteomics study of human aortic ECM (PMID: 20551380), classified by the BHF-UCL curation group. This is a computational classification based on finding the protein in an ECM proteomics dataset — it does not demonstrate that PCPE-1 structurally contributes to the ECM.

As stated by Berry et al. (PMID: 30078642): "Procollagen C-proteinase enhancer-1 (PCPE-1) is a secreted protein that specifically accelerates proteolytic release of the C-propeptides from fibrillar procollagens, a crucial step in fibril assembly." This is a regulatory/catalytic enhancement function, fundamentally different from structural ECM support. Even if A0A2G9RZF1 were a genuine PCPE-1 ortholog (which it is not), GO:0016504 (peptidase activator activity), not GO:0005201, would be the appropriate molecular function term.

Finding 3: A Single CUB Domain Cannot Perform PCPE-1 Function

Even setting aside the semantic mismatch with GO:0005201, transferring any PCPE-1 function to A0A2G9RZF1 is unjustified because a single CUB domain is insufficient. PCPE-1 is a 449-amino acid protein with three functional domains: CUB1 (37–149), CUB2 (159–273), and NTR (318–437). Both CUB domains and the NTR domain contribute to PCPE-1's biological activities.

Critical experimental evidence from Blanc et al. (PMID: 19801683) demonstrated conclusively that "only those containing both CUB1 and CUB2 were capable of enhancing BMP-1 activity and binding to a mini-procollagen substrate with nanomolar affinity. Both these properties were lost by individual CUB domains, which had dissociation constants at least three orders of magnitude higher." This >1,000-fold loss of binding affinity means that a single CUB domain fragment cannot perform any of PCPE-1's characterized functions. The FunFam classification, while computationally valid at the domain level, does not support functional transfer to a protein containing only one of the required domains.

Bourhis et al. (PMID: 17446170) further confirmed that "Procollagen C-proteinase enhancers (PCPE-1 and -2) are extracellular glycoproteins that can stimulate the C-terminal processing of fibrillar procollagens by tolloid proteinases such as bone morphogenetic protein-1. They consist of two CUB domains (CUB1 and -2) that alone account for PCPE-enhancing activity and one C-terminal NTR domain." A0A2G9RZF1 at 156 amino acids possesses none of this required architecture.

Finding 4: A0A2G9RZF1 Is a Gene Prediction Fragment of a BMP-1/Tolloid-Like Metalloprotease

The most decisive finding emerged from examining the genomic context. Scaffold KV928989 (134 kb) from the RCv2.1 assembly encodes three consecutive open reading frames:

The combined domain modules from these three ORFs (CUB – EGF – CUB – CUB – EGF – CUB) match the C-terminal region of a tolloid-like (mTLD) protease architecture. The neighboring ORF A0A2G9RZI6 is explicitly marked as "Fragment" in UniProt and classified by PANTHER as BMP-1. The scaffold contains multiple assembly gaps. No complete (>700 aa) BMP-1/tolloid-like protein was found elsewhere in the bullfrog proteome from this assembly, consistent with the gene being split across gaps.

Finding 5: CUB Domains in Tolloid Proteases Function as Substrate-Recognition Modules

If A0A2G9RZF1's CUB domain is part of a tolloid-like protease, its function would be substrate recognition and presentation — not ECM structural support. Lee et al. (PMID: 18664565) demonstrated in Xenopus Xolloid that "the first and second CUB domains bind Chordin and present it to the protease domain." This substrate-recognition function is the canonical role of CUB domains in the tolloid/BMP-1 protease family and is fundamentally distinct from structural ECM support.

Additional evidence from Drosophila tolloid (PMID: 25642644) showed that N-terminal CUB domains interact with Collagen IV to enhance Tolloid activity toward its substrate Sog, while C-terminal CUB domains mediate Sog interaction. This bipartite CUB domain function (ECM anchoring + substrate binding) fine-tunes protease activity but is entirely regulatory, not structural.

Finding 6: Chromosome-Level Assembly Confirms Gene Fragmentation

The definitive confirmation came from the chromosome-level assembly GCF_042186555.1 (ASM4218655v1, 2024), which encodes:

• Complete BMP-1: XP_073478370.1, 1,020 aa, 16 exons on LG03
• Complete TLL1: XP_073462190.1, 1,005 aa, 21 exons on LG01

These are full-length tolloid-family metalloproteases with the expected multi-domain architecture. The original RCv2.1 assembly had a scaffold N50 of only 39,368 bp compared to 691,824,178 bp in the chromosome-level assembly — a ~17,000-fold improvement in contiguity. The three fragmented ORFs on scaffold KV928989 (156 + 228 + 250 = 634 aa combined) represent portions of one of these ~1,000 aa genes, with the N-terminal protease domain and additional domains falling in assembly gaps.

This finding is critical for curation: the UniProt proteome (UP000228934) still references the 2017 draft assembly. When updated to the chromosome-level assembly, these fragment entries should be superseded by complete gene models.

Evidence Matrix

GO Curation Implications

Recommended Action: Do NOT Annotate with GO:0005201; Leave MF Unassigned

GO:0005201 (extracellular matrix structural constituent) should not be applied to A0A2G9RZF1 under any evidence code. The term is semantically incorrect for a CUB-domain protein, the evidence chain through PCPE-1 FunFam classification does not support it even for PCPE-1 itself, and the protein is a gene prediction artifact.

GO Decision Table:

Key Distinction

Even if A0A2G9RZF1 were a standalone protein (which genomic evidence refutes), the appropriate MF for a single-CUB protein would relate to protein binding in the extracellular space — never GO:0005201. The seed hypothesis itself correctly notes that "the short length (156 aa) and absence of additional functional domains preclude confident functional assignment." This assessment is confirmed and strengthened by our analysis.

Mechanistic Scope

Direct Gene-Product Activity Under Test

The hypothesis tests whether A0A2G9RZF1 directly functions as a structural component of the extracellular matrix — physically contributing to ECM architecture by being incorporated into the matrix scaffold, analogous to collagens, proteoglycans, or elastin.

What the CUB Domain Actually Does

CUB domains are non-catalytic protein-protein interaction modules. In the tolloid-family protease context (the most likely identity for the complete gene), CUB domains function as:

1. Substrate-recognition modules — binding target proteins (e.g., Chordin, procollagen C-propeptides) and presenting them to the catalytic protease domain (PMID: 18664565)
2. ECM-anchoring modules — interacting with Collagen IV to localize the protease (PMID: 25642644)

Neither function constitutes "structural ECM support."

Separation of Activities

Conflicts and Alternatives

1. FunFam Classification to PCPE-1

The CATH FunFam classification (2.60.120.290:FF:000005) groups A0A2G9RZF1 with PCPE-1. While computationally valid at the domain-fold level, this classification does not imply functional equivalence. PCPE-1 requires CUB1 + CUB2 + NTR for its characterized activities. The FunFam grouping reflects structural similarity of the CUB fold, not functional annotation transfer. CUB domains are found across functionally diverse proteins: complement factors (C1r, C1s), endocytic receptors (cubilin — 27 CUB domains; PMID: 30295181), spermadhesins, developmental proteases, neurotransmitter receptor modulators (PMID: 21093502), and innate immune receptors.

2. PANTHER OVOCHYMASE-RELATED Classification

PANTHER classifies A0A2G9RZF1 under PTHR24251 (OVOCHYMASE-RELATED). Ovochymase is a Xenopus egg polyprotein with multiple CUB domains interspersed between serine protease domains (PMID: 10500163). This family-level classification is broad and includes proteins with diverse functions (metalloproteases, serine proteases, enhancers). It should not be used to infer a specific molecular function for a single-CUB fragment.

3. Could A0A2G9RZF1 Be a Genuine Single-CUB Protein?

This alternative is unlikely but not absolutely excludable. Single CUB-domain proteins exist (e.g., spermadhesins PSP-I/PSP-II in pig seminal plasma; PMID: 16819821). However, arguments against this interpretation are strong:

• Protein existence level 4 (predicted, no transcript/protein evidence)
• TrEMBL (unreviewed) status
• Genomic context: three consecutive ORFs with complementary domains on a gap-containing scaffold
• No complete tolloid gene found elsewhere in the same assembly
• Chromosome-level assembly encodes full-length versions (1,005–1,020 aa)
• ~17,000-fold improvement in assembly contiguity resolved the fragmentation

Even if A0A2G9RZF1 were standalone, its function would be extracellular protein binding — still NOT GO:0005201.

4. Paralog/Family Confusion Risk

The PTHR24251 family includes ovochymase (egg protease with 5 CUB domains), BMP-1/tolloid metalloproteases, PCPE-1/2 enhancers, and various uncharacterized CUB-containing proteins. Homology at the domain level does not predict function at the protein level. The FunFam classification likely reflects the general CUB fold similarity rather than specific PCPE-1 functional identity.

Knowledge Gaps

Discriminating Tests

Computational (Immediately Feasible)

1. Scaffold-to-chromosome alignment: BLAST or minimap2 alignment of scaffold KV928989 against GCF_042186555.1 would definitively identify which complete gene (BMP-1 or TLL1) corresponds to the A0A2G9RZF1 locus. Cost: minutes of computation.

2. PCPE-1 ortholog search: Reciprocal best BLAST of human PCPE-1 (Q15113) against the chromosome-level bullfrog proteome would identify whether a genuine PCPE-1 ortholog exists, independent of A0A2G9RZF1.

3. RNA-seq validation: Search SRA/ENA for A. catesbeiana transcriptome data and map to the chromosome-level assembly to confirm expression of the complete tolloid gene.

4. AlphaFold/ESMFold: Structure prediction for A0A2G9RZF1 to assess whether it folds into a stable CUB domain with intact calcium-binding site — informative about whether even the fragment is structurally viable.

Experimental (Definitive)

1. BMP-1 enhancer activity assay: Test whether the isolated CUB domain can enhance BMP-1 activity on procollagens. Expected result: negative, based on PMID: 19801683.

2. SPR/ITC binding assays: Test binding of the isolated CUB domain to known tolloid substrates (Chordin, procollagen C-propeptide) and ECM structural proteins. Would directly distinguish substrate-recognition from structural roles.

3. Full-length gene cloning: Clone the complete BMP-1/TLL1 gene from A. catesbeiana cDNA and characterize its enzymatic activity.

Curation Leads

Lead 1: Remove GO:0005201 from Core Functions (HIGH PRIORITY)

Action: Do not annotate A0A2G9RZF1 with GO:0005201. Remove this term from consideration as a core function.

Rationale: Three independent evidence lines refute this annotation (semantic mismatch, PCPE-1 function mismatch, gene fragmentation). No experimental evidence supports it.

References to verify: - PMID: 19801683 — "Out of all the forms tested, only those containing both CUB1 and CUB2 were capable of enhancing BMP-1 activity" - PMID: 21954942 — "CUB domains are 110-residue protein motifs exhibiting a β-sandwich fold and mediating protein-protein interactions" - UniProt Q15113 GO annotations — PCPE-1 IDA terms are GO:0005518, GO:0008201, GO:0016504 (not GO:0005201)

Confidence: High.

Lead 2: Flag A0A2G9RZF1 as Gene Prediction Fragment (HIGH PRIORITY)

Action: Flag A0A2G9RZF1, A0A2G9RZH1, and A0A2G9RZI6 as probable fragments of a single BMP-1/tolloid-like gene. Do not annotate fragments with function-level GO terms.

Supporting evidence: - Three consecutive ORFs on scaffold KV928989 with complementary tolloid-family domains - Scaffold N50 of 39,368 bp in original assembly - Chromosome-level assembly encodes complete BMP-1 (1,020 aa) and TLL1 (1,005 aa) - A0A2G9RZI6 explicitly marked as "Fragment" and classified as BMP-1

Reference to verify: PMID: 18664565 — CUB domains in tolloid = substrate recognition modules, not structural ECM

Confidence: High.

Lead 3: Confirm Rejection of Metalloendopeptidase Activity (MODERATE)

Action: Confirm the seed hypothesis's recommendation against annotating metalloendopeptidase activity (GO:0004222) for this fragment. The protease domain is absent from A0A2G9RZF1.

Note: If the complete gene product is identified as BMP-1 or TLL1, GO:0004222 would be appropriate for the full-length protein — but not for this CUB-domain-only fragment.

Confidence: High.

Lead 4: Retain GO:0005576 as CC Annotation (MODERATE)

Action: Retain GO:0005576 (extracellular region) as a Cellular Component annotation with IEA-level evidence, based on the near-universal extracellular localization of CUB-domain proteins.

Confidence: Moderate.

Lead 5: Leave Molecular Function Unassigned (HIGH PRIORITY)

Action: Do not assign any specific MF term. The function is genuinely unknown for this fragment. The original seed hypothesis statement — "The precise molecular function of A0A2G9RZF1 is unknown" — is the most accurate assessment and should be retained.

Confidence: High.

Lead 6: Reference Complete Gene Models for Future Curation (INFORMATIONAL)

Action: When the UniProt proteome is updated to the chromosome-level assembly (GCF_042186555.1), A0A2G9RZF1 should be superseded by: - BMP-1: XP_073478370.1 (1,020 aa, 16 exons, LG03) — NCBI Gene 141133120 - TLL1: XP_073462190.1 (1,005 aa, 21 exons, LG01) — NCBI Gene 141113146

These complete gene products would carry appropriate metalloendopeptidase and substrate-binding annotations.

Mechanistic Model

The following diagram summarizes the evidence architecture:

SEED HYPOTHESIS CHAIN (refuted at each step):

  A0A2G9RZF1 ──FunFam──> PCPE-1 ──RCA──> GO:0005201     
  (156 aa,        ↓         ↓              ↓              
   1 CUB)     FOLD MATCH  WEAK RCA    SEMANTIC MISMATCH   
                  ↓         ↓              ↓              
              Single CUB  IDA terms    GO:0005201 = 
              cannot do    are NOT      collagens,        
              PCPE-1 work  GO:0005201   elastin           

ACTUAL IDENTITY (supported):

  Scaffold KV928989 (fragmented assembly, N50=39kb):
  ┌─────────────────────────────────────────────────┐
  │ ...gap... [CUB] ...gap... [EGF-CUB] ...gap...  │
  │           ↑                ↑                    │
  │     A0A2G9RZF1      A0A2G9RZH1                 │
  │                                                 │
  │ [CUB-EGF-CUB] ...gap...                        │
  │       ↑                                         │
  │  A0A2G9RZI6 (PANTHER: BMP-1, "Fragment")        │
  └─────────────────────────────────────────────────┘
                        ↓ resolves to
  Chromosome-level assembly (N50=692Mb):
  ┌─────────────────────────────────────────────────┐
  │ BMP-1 (1020 aa): Protease-CUB-EGF-CUB-CUB-EGF │
  │ or                                              │
  │ TLL1 (1005 aa): Protease-CUB-EGF-CUB-CUB-EGF  │
  └─────────────────────────────────────────────────┘

Interpretation: A0A2G9RZF1 is one CUB domain from a tolloid-family metalloprotease that was incorrectly predicted as a standalone gene due to assembly fragmentation. The CUB domain's function within the complete protein would be substrate recognition — binding target proteins and presenting them to the protease domain for cleavage. This is a regulatory/binding function within a larger enzyme, not an ECM structural role.

Evidence Base: Key Literature

Core Evidence Against GO:0005201

Gaboriaud et al. (2011) — Structure and properties of the Ca(2+)-binding CUB domain, a widespread ligand-recognition unit involved in major biological functions. PMID: 21954942 Comprehensive review establishing CUB domains as protein-protein interaction modules with β-sandwich fold. Key evidence: CUB domains mediate ligand recognition across diverse biological systems (complement, development, endocytosis, reproduction) — they are interaction modules, not structural components.

Blanc et al. (2009) — Strong cooperativity and loose geometry between CUB domains are the basis for procollagen C-proteinase enhancer activity. PMID: 19801683 Demonstrated that individual CUB domains lose >1,000-fold binding affinity. Critical for ruling out functional transfer from PCPE-1 to single-CUB proteins. This is the strongest single piece of evidence against the PCPE-1-based inference chain.

Bourhis et al. (2007) — Insights into how CUB domains can exert specific functions while sharing a common fold. PMID: 17446170 Confirmed the CUB1-CUB2-NTR architecture requirement for PCPE function. Identified PCPE-specific residues in CUB1 that are necessary but insufficient without CUB2 cooperativity.

Evidence for Tolloid Fragment Identity

Lee et al. (2009) — Molecular determinants of Xolloid action in vivo. PMID: 18664565 Demonstrated CUB domain function in amphibian tolloid proteases: substrate recognition and presentation to the protease domain. Directly relevant as Xenopus is the closest well-characterized amphibian model.

Winstanley et al. (2015) — Synthetic enzyme-substrate tethering obviates the Tolloid-ECM interaction during Drosophila BMP gradient formation. PMID: 25642644 Showed bipartite CUB domain function in tolloid: ECM anchoring (N-terminal CUBs) and substrate binding (C-terminal CUBs). Confirms CUB domains in tolloids serve regulatory/binding, not structural, roles.

PCPE-1 Function Characterization

Berry et al. (2018) — Structural Basis for the Acceleration of Procollagen Processing by Procollagen C-Proteinase Enhancer-1. PMID: 30078642) Crystal structure of PCPE-1 revealing its mechanism as a regulatory accelerator of procollagen processing — not an ECM structural protein.

Salza et al. (2014) — Extended interaction network of procollagen C-proteinase enhancer-1 in the extracellular matrix. PMID: 24117177 Identified 17 binding partners of PCPE-1, confirming its role as an interaction hub, not structural scaffold. CUB1CUB2 fragment inhibits angiogenesis — a regulatory function.

Limitations

1. No direct experimental evidence exists for A0A2G9RZF1. All conclusions are based on computational analysis, domain architecture reasoning, genomic context, and inference from well-characterized homologs in other species.

2. The fragment-to-gene assignment is inferential. We have not performed the BLAST alignment of KV928989 against the chromosome-level assembly that would definitively identify which gene (BMP-1 or TLL1) the fragment belongs to.

3. Assembly-based reasoning has inherent uncertainty. While the evidence is strong that A0A2G9RZF1 is a fragment, we cannot exclude unusual genomic rearrangements or lineage-specific gene fission events in bullfrog, though these would be extraordinary claims requiring extraordinary evidence.

4. Literature is from model organisms. CUB domain function, PCPE-1 biochemistry, and tolloid protease characterization are primarily from human, mouse, Xenopus, and Drosophila. Direct evidence from Aquarana catesbeiana is absent.

5. The 32 papers reviewed focused on PCPE-1 and tolloid biology. We did not exhaustively survey all possible functions of isolated CUB domains in non-model amphibians, though no evidence from the reviewed literature suggested an ECM structural role for any CUB-domain protein.

Proposed Follow-up Actions

Immediate (Computational)

1. Scaffold alignment: BLAST scaffold KV928989 against the chromosome-level assembly to identify the exact gene correspondence (BMP-1 or TLL1).
2. Proteome refresh: Request or monitor UniProt proteome update for A. catesbeiana to the chromosome-level assembly.

Short-term (Curation)

1. Remove GO:0005201: Remove this term from consideration for A0A2G9RZF1 and any derived annotation pipelines.
2. Flag as fragment: Mark A0A2G9RZF1 as a gene prediction fragment in the curation system.
3. No MF annotation: Leave molecular function unassigned on this entry.

Medium-term (Analysis)

1. Annotate complete genes: Apply appropriate annotations (GO:0004222 metalloendopeptidase activity, GO:0006508 proteolysis, GO:0005576 extracellular region) to the complete BMP-1 (XP_073478370.1) and TLL1 (XP_073462190.1) from the chromosome-level assembly.
2. PCPE-1 ortholog identification: Search the chromosome-level proteome for genuine PCPE-1/PCOLCE1 orthologs and annotate appropriately with GO:0016504.

Report generated by autonomous scientific discovery agent across 3 investigation iterations. 32 papers reviewed, 6 findings confirmed. All conclusions are computational and require curator verification.