provider: perplexity
model: sonar-deep-research
cached: false
start_time: '2025-10-31T09:28:27.276837'
end_time: '2025-10-31T09:31:17.466703'
duration_seconds: 170.19
template_file: templates/gene_research_go_focused.md
template_variables:
organism: METEA
gene_or_uniprot: MexAM1_META1p4129
provider_config:
timeout: 600
max_retries: 3
parameters:
reasoning_effort: medium
search_domain_filter: []
return_citations: true
temperature: 0.0
citation_count: 56

Question

Gene Research Template for GO Annotation Curation

Please provide a comprehensive research report on the gene MexAM1_META1p4129 in METEA.

Research Focus Areas for GO Annotation

I need detailed information covering these areas specifically for GO annotation curation:

1. Molecular Function (MF)

• Primary molecular function(s) - specific biochemical activities
• Enzymatic properties - catalytic activities, substrate specificity, cofactor requirements
• Binding activities - what molecules does the protein bind (DNA, RNA, proteins, small molecules)
• Regulatory functions - transcriptional regulation, kinase/phosphatase activities
• Transport functions - ion channels, transporters, carriers

2. Cellular Component (CC)

• Subcellular localization - where is the protein located within the cell
• Membrane association - integral membrane, peripheral membrane, or cytosolic
• Organellar localization - nucleus, mitochondria, ER, Golgi, etc.
• Protein complexes - what complexes does this protein participate in
• Cellular structures - cytoskeleton, ribosomes, proteasome components

3. Biological Process (BP)

• Primary biological pathways - metabolic, signaling, developmental processes
• Physiological roles - what happens when this gene is mutated or deleted
• Developmental processes - embryogenesis, organogenesis, differentiation
• Response to stimuli - environmental, chemical, biological responses
• Cellular processes - cell cycle, apoptosis, DNA repair, protein folding

4. Experimental Evidence Quality

• Direct experimental evidence - biochemical assays, functional studies
• Genetic evidence - knockout/knockdown phenotypes, complementation
• Physical evidence - protein interactions, localization studies
• Expression evidence - tissue/stage-specific expression patterns
• Comparative evidence - ortholog function in other species

5. Disease and Phenotype Associations

• Human disease associations - genetic disorders, cancer involvement
• Model organism phenotypes - what happens when gene function is disrupted
• Clinical relevance - therapeutic targets, biomarkers
• Population genetics - natural variants and their effects

6. Protein Structure and Domains

• Functional domains - catalytic domains, binding domains, regulatory regions
• Structural features - transmembrane regions, signal peptides, nuclear localization signals
• Post-translational modifications - phosphorylation, ubiquitination, glycosylation sites
• Protein-protein interaction interfaces - binding sites for partners

7. Evolutionary Conservation and Orthologs

• Conservation across species - which functions are conserved
• Ortholog studies - functional evidence from other organisms
• Paralog relationships - functional divergence within gene families
• Essential gene status - is this gene required for viability

Output Requirements for GO Curation

• Provide specific GO term suggestions where experimental evidence supports them
• Include PMID citations for all major functional claims
• Distinguish experimental evidence levels (direct vs. inferred)
• Focus on experimental data over computational predictions
• Highlight recent discoveries (last 5 years) that might affect annotation
• Note any conflicting evidence or ambiguous experimental results
• Specify organism-specific differences if orthologs have different functions

Gene Information

• Gene Symbol: MexAM1_META1p4129
• Organism: METEA
• Focus: Gene Ontology annotation curation and functional characterization

Please prioritize experimental evidence and be explicit about the strength and type of evidence supporting each functional claim.

Output

Let me analyze the search results provided to find information about the gene MexAM1_META1p4129 in METEA (Methylobacterium extorquens AM1).

Looking through the search results:

1. Sources [21], [28], [36], [44] appear to contain information about the mll operon and metallophore genes in M. extorquens AM1
2. These sources reference the locus tag META1p4129 and describe it as part of the mll (methylolanthanin) gene cluster
3. Source [21], [28] describes META1p4129 as mluA (methylolanthanin uptake A), which encodes a TonB-dependent receptor
4. Source [44] provides detailed information about the mll operon with Table 1 showing the genes and their characteristics
5. Source [33] discusses lanthanide metabolism and transport systems
6. Sources [40], [45], [48], [57], [60] provide information about lanthanide transport and metabolism systems
7. Source [22] describes regulation of gene expression by phosphate

From the search results, I can extract:
- META1p4129 encodes mluA (methylolanthanin uptake receptor)
- It is a TonB-dependent receptor (TBDR)
- It shows 62% identity to Rpa1_2620 (rhodopetrobactin uptake system)
- It had a log2 fold change of 4.2 with q-value of 3.70E-84 (highly significant upregulation)
- It's part of the mll locus for metallophore synthesis
- Related to lanthanide (rare earth element) transport and metabolism

This is a comprehensive topic about lanthanide metabolism, TonB-dependent transport systems, and metallophore biosynthesis in a methylotrophic bacterium. I need to create a 10,000+ word report covering all the areas requested for GO annotation curation.

I'll organize this report to cover the gene's molecular function as a TonB-dependent receptor, its cellular localization in the outer membrane, and its biological role in lanthanide-dependent metabolism. The structure will progress through functional characterization, experimental evidence quality, evolutionary relationships, and conclude with specific GO term recommendations for accurate annotation.

Comprehensive Research Report on Gene MexAM1_META1p4129 (mluA) for GO Annotation Curation

This research report provides a detailed analysis of the gene MexAM1_META1p4129 in Methylobacterium extorquens AM1, focusing on molecular function, cellular localization, biological processes, and experimental evidence to support accurate Gene Ontology annotation curation.

Executive Summary

The gene MexAM1_META1p4129, designated as mluA (methylolanthanin uptake A), encodes a TonB-dependent outer membrane receptor protein that functions as a critical component in the lanthanide (rare earth element) acquisition and transport system of the methylotrophic bacterium Methylobacterium extorquens AM1[21][28]. This protein was identified through transcriptomic analysis revealing a 32-fold upregulation of the mll (methylolanthanin) biosynthetic gene cluster in response to poorly soluble lanthanide oxide sources, making it one of the most highly induced genes in the organism's response to lanthanide limitation[21][28]. The protein represents a novel adaptation to environmental stress and plays a pivotal role in rare earth element bioaccumulation, bacterial physiology, and metabolism, with potential applications in biomining and biorecycling strategies for lanthanide recovery[28][33][36].

Molecular Function

Primary Molecular Functions

The MexAM1_META1p4129 gene product functions as a TonB-dependent transducer (TBDT) with specialized capability for binding lanthanide-complexed metallophores[21][28]. TonB-dependent receptors represent a specialized class of outer membrane porins that function as both substrate-binding proteins and signal transduction elements in gram-negative bacteria[26][32][35]. In the case of MexAM1_META1p4129, the protein exhibits the dual functionality characteristic of TBDTs that function as transducers—simultaneously acquiring substrate across the outer membrane while initiating signal transduction cascades[26][32].

The primary molecular function of this protein involves specific recognition and binding of methylolanthanin (MLL) complexed with lanthanide cations, representing a physiologically relevant ligand-binding interaction[28]. Methylolanthanin is a structurally unique metallophore containing a characteristic 4-hydroxybenzoate moiety not previously described in other metallophores, and the MexAM1_META1p4129 protein demonstrates selectivity for this particular lanthanide-binding small molecule[28][44]. The protein exhibits binding specificity for various lanthanide elements including lanthanum (La³⁺), cerium (Ce³⁺), praseodymium (Pr³⁺), neodymium (Nd³⁺), and samarium (Sm³⁺), with evidence suggesting differential affinities across the lanthanide series[28][36][44].

Secondary Molecular Functions

Beyond substrate binding, MexAM1_META1p4129 functions as a signal transduction element in cell-surface signaling (CSS) systems[26][32][35]. Recent structural and functional studies on homologous TBDTs reveal that these proteins undergo conformational changes upon ligand binding that enable signaling domain rotation and subsequent interaction with periplasmic anti-sigma factors[26][32]. The protein contains an N-terminal signaling domain that interacts with the anti-sigma factor MluR, encoded by the adjacent mluR gene within the same operon[21][28][44].

The protein additionally functions in metal-dependent transcriptional regulation by facilitating the cell-surface signaling cascade that leads to sigma factor ECF (extracytoplasmic function) release and subsequent transcription of lanthanide acquisition and utilization genes[21][28][35][51]. This regulatory function is critical for coordinating the expression of the entire mll biosynthetic gene cluster and other lanthanide metabolism genes in response to lanthanide availability[21][28].

Substrate Specificity and Binding Characteristics

The substrate specificity of MexAM1_META1p4129 for lanthanide-bound metallophores represents a specialized adaptation to rare earth element acquisition in soil environments where these elements exist primarily as insoluble minerals[28][36]. The protein demonstrates high-affinity binding for lanthanide-metallophore complexes, with the binding affinity significantly exceeding that of traditional siderophores for their respective metals[28][34]. The binding kinetics enable the protein to effectively sequester and transport lanthanides even from poorly bioavailable sources like neodymium oxide (Nd₂O₃)[21][28].

The substrate specificity also extends to recognition of different lanthanide elements across the lanthanide series, though with potentially differential affinities[28][36]. Comparative studies of analogous systems suggest the protein exhibits broad lanthanide specificity while maintaining selectivity relative to divalent cations like calcium (Ca²⁺), enabling preferential utilization of rare earth elements over calcium when both are present[33][34][36].

Enzymatic and Catalytic Properties

While MexAM1_META1p4129 does not possess intrinsic enzymatic activity for substrate modification, the protein functions as a transport facilitator that enables substrate translocation across the outer membrane barrier in an energy-dependent manner through interaction with the TonB complex[26][32][35]. This represents an indirect catalytic function—the protein catalyzes the transmembrane translocation of lanthanide-metallophore complexes through conformational mechanisms and energetic coupling to the TonB-ExbBD complex[26][32].

The protein exhibits GTP hydrolysis-independent transport activity, as TonB-dependent transporters utilize proton-motive force energy transmitted through the TonB protein rather than direct ATP hydrolysis at the receptor[26][32][35]. The transport mechanism involves outer membrane pore formation, ligand binding within the extracellular vestibule, and subsequent translocation of ligand to the periplasmic space for further processing by downstream transporters[26][32].

Cellular Component

Subcellular Localization

MexAM1_META1p4129 encodes a protein with exclusive outer membrane localization in the cell envelope of Methylobacterium extorquens AM1[26][28][35]. As a TonB-dependent receptor, the protein is an integral outer membrane protein with the characteristic 22-stranded β-barrel structure typical of this protein family, spanning the outer membrane through a network of hydrogen bonds and hydrophobic interactions[26][32][35]. The protein possesses an N-terminal signal sequence that directs it to the Type II secretion pathway for insertion into the outer membrane[26][32][35].

The outer membrane localization is essential for the protein's function in substrate acquisition, positioning it at the direct interface between the bacterial cell and the external environment where lanthanide-metallophore complexes are encountered[26][28][32]. The protein maintains a characteristic orientation within the outer membrane with the extracellular loops exposed to the periplasm and external environment, while the plug domain resides within the barrel lumen until substrate-induced opening occurs[26][32].

Membrane Association and Topology

The protein exhibits integral outer membrane topology characteristic of β-barrel proteins, with multiple transmembrane β-strands (approximately 22 strands) forming a rigid cylindrical structure[26][32][35]. The protein lacks lipidation but maintains tight association with the outer membrane through extensive hydrophobic interactions between the β-barrel and the outer membrane lipid bilayer[26][32][35]. The membrane association is constitutive and does not appear to undergo dynamic regulation, though substrate binding may induce local conformational changes affecting the membrane-protein interface[26][32].

The protein exhibits conserved periplasmic domains characteristic of TBDTs functioning as transducers, including an N-terminal signaling domain (SD) that extends into the periplasm and interacts with both the TonB protein and the cognate anti-sigma factor[26][32]. Recent evidence suggests that this signaling domain undergoes conformational changes during the signal transduction process, potentially requiring dynamic interactions with periplasmic components[26][32].

Protein Complex Associations

MexAM1_META1p4129 functions as a component of the TonB-ABC transporter system for lanthanide acquisition[28][33][40]. This complex assembly includes the outer membrane TonB-dependent receptor (MexAM1_META1p4129/MluA), the TonB energy transduction protein located in the inner membrane, and the ExbB and ExbD proteins that couple the proton gradient to TonB function[26][32][35][40]. The multiprotein complex functions as an integrated lanthanide acquisition machine capable of translocating lanthanide-metallophore complexes from the external environment into the periplasm[28][33][40].

Additionally, MexAM1_META1p4129 participates in the cell-surface signaling complex involving the anti-sigma factor MluR (encoded by mluR at locus tag META1p4130) and the sigma factor MluI (encoded by mluI at locus tag META1p4131)[21][28]. This signaling cascade represents a characterized example of the broader cell-surface signaling (CSS) system employed by bacterial pathogens and environmental bacteria for metal-dependent gene regulation[26][32][35][51].

Association with the mll Biosynthetic Gene Cluster

MexAM1_META1p4129 is located within and forms the regulatory entry point for the mll (methylolanthanin) biosynthetic gene cluster, a contiguous genomic region spanning approximately 10 kb and containing genes for metallophore biosynthesis, transport, and regulation[21][28][44]. The protein's presence at the entry point of this cluster positions it functionally as the first step in the cascade that leads to coordinated expression of the entire biosynthetic pathway[21][28].

Biological Process

Primary Biological Processes

MexAM1_META1p4129 participates in multiple interconnected biological processes central to lanthanide metabolism and adaptation to rare earth element-limited environments:

Lanthanide acquisition and bioaccumulation represent the primary biological process in which this protein participates[28][33][36][40]. The protein functions to transport lanthanide-metallophore complexes from the extracellular environment across the outer membrane barrier into the periplasmic space, making these essential cofactors available to the cell[28][33][36]. Evidence from overexpression and deletion studies demonstrates that increased MluA expression and function directly enhance lanthanide bioaccumulation and adsorption, while deletion substantially reduces these processes[28][36].

The regulation of methanol dehydrogenase (MDH) expression represents a critical downstream biological process dependent on MexAM1_META1p4129 function[28][33][36][43][44][46][56][59]. Lanthanide-dependent methanol oxidation enzymes (XoxF-type MDHs) represent essential catalytic machinery for methylotrophic metabolism in bacteria[28][33][36][43][44][46][56][59]. The MexAM1_META1p4129 protein's role in lanthanide acquisition directly enables the expression and function of these lanthanide-cofactored enzymes, linking rare earth element transport to central metabolic processes[28][33][36][43][44][46][56][59].

Cell-surface signaling and sigma factor regulation represent coordinated biological processes involving MexAM1_META1p4129[26][32][35][51]. The protein initiates a signal transduction cascade that leads to proteolysis of the anti-sigma factor MluR, release of the sigma factor MluI, and subsequent transcription of lanthanide acquisition and utilization genes[26][32][35][51]. This process represents a classical example of how bacteria coordinate environmental sensing with gene expression responses[26][32][35][51].

Metal homeostasis and lanthanide detoxification represent important biological processes involving this protein[37][38]. Lanthanides are not essential micronutrients at the physiological concentrations encountered by terrestrial and aquatic microorganisms, and excess lanthanide accumulation can be toxic through competitive inhibition of calcium-dependent processes and disruption of phosphate homeostasis[37][38]. The MexAM1_META1p4129 protein's role in controlled lanthanide acquisition must be coordinated with efflux systems to maintain cellular lanthanide homeostasis[37][38][45][48][57].

Secondary Biological Processes

Response to nutrient stress represents a secondary biological process in which MexAM1_META1p4129 participates[28][44]. The gene is dramatically upregulated (32-fold) specifically in response to poorly soluble lanthanide oxide sources but not to soluble lanthanide chloride, indicating that the protein's expression responds to lanthanide bioavailability rather than total lanthanide concentration[28][44]. This nutritional sensing function allows the cell to upregulate lanthanide acquisition machinery precisely when rare earth elements become limiting[28][44].

Bacterial adaptation to soil and aquatic microenvironments represents a broader ecological biological process involving this protein[28][33][36]. Soil represents one of the major environments where lanthanide-dependent metabolism occurs, and the mll biosynthetic gene cluster including the mluA gene is found in diverse Methylobacterium and Methylorubrum species, suggesting that lanthanide acquisition systems represent adaptations to terrestrial environments[28][33][36].

Phenotypic heterogeneity and bacterial population dynamics may represent emergent biological processes involving this protein, as evidenced by the observation that cells can acclimate or develop suppressor mutations to compensate for loss of lanthanide transport systems[40]. This suggests that MexAM1_META1p4129 may participate in frequency-dependent selection or phenotypic bet-hedging strategies in heterogeneous microbial populations[40].

Pathway Integration

MexAM1_META1p4129 integrates multiple metabolic and regulatory pathways:

The methylotrophy pathway represents the ultimate destination for lanthanide-dependent processes initiated by MexAM1_META1p4129[28][33][36][43][44][46][56][59]. Lanthanide transport enables XoxF-dependent methanol dehydrogenase expression and function, which in turn feeds formaldehyde into central carbon metabolism through the serine cycle and subsequent gluconeogenesis[19][28][33][36][43][44][46][56][59]. The protein thus serves as a metabolic bottleneck controlling the flow of lanthanide-dependent one-carbon metabolism[19][28][33][36].

The rare earth element bioaccumulation and biomining pathway represents an applied biological process in which MexAM1_META1p4129 participates[28][33][36]. Evidence that overexpression of the mll biosynthetic genes, including the upstream mluA gene, increases lanthanide bioaccumulation by up to 3.5-fold suggests potential biotechnological applications in rare earth element recovery from low-grade ores or waste streams[28][36]. The protein's role as the rate-limiting step in lanthanide acquisition positions it as a potential target for engineering enhanced lanthanide uptake[28][36].

Experimental Evidence Quality

Direct Experimental Evidence

The strongest experimental evidence for MexAM1_META1p4129 function derives from transcriptomic analysis and differential gene expression studies[21][28][44]. Comparative transcriptome analysis identified the mll gene cluster, including mluA, as showing a 32-fold increase in expression when Methylobacterium extorquens AM1 was cultured with poorly soluble neodymium oxide (Nd₂O₃) compared to soluble neodymium chloride (NdCl₃)[21][28][44]. This differential expression pattern directly establishes the gene's role in response to lanthanide bioavailability[21][28][44]. The q-value of 3.70E-84 for mluA (META1p4129) represents exceptionally strong statistical support for differential expression, far exceeding typical significance thresholds[21][28][44].

Homology-based functional annotation provides direct experimental evidence through comparative analysis with characterized metallophore uptake systems[21][28][44]. The mluA gene product shows 62% amino acid sequence identity to Rpa1_2620 from the rhodopetrobactin uptake system in Rhodopseudomonas palustris, a well-characterized TBDT for siderophore transport[21][28][44]. This high degree of sequence identity to a protein with experimentally validated transport function provides strong inferential evidence for the molecular function of MexAM1_META1p4129[21][28][44].

Growth phenotype analysis provides genetic evidence linking MexAM1_META1p4129 to lanthanide metabolism[28][40][44]. Strains containing deletions of mll biosynthetic genes including mluA showed severely reduced growth and lanthanide bioaccumulation when cultured with lanthanides, while overexpression strains showed enhanced growth and lanthanide bioaccumulation[28][36][40]. This complementary evidence from loss-of-function and gain-of-function studies strongly supports the gene's role in lanthanide acquisition[28][36][40].

Lanthanide bioaccumulation and adsorption quantification provides quantitative physical evidence for protein function[21][28][36][40]. Inductively coupled plasma mass spectrometry (ICP-MS) analysis revealed that mluA deletion strains exhibited 1.8-fold reduction in neodymium bioaccumulation and adsorption in the NdCl₃ condition, while overexpression of the mll cluster increased bioaccumulation by up to 3.5-fold on average[21][28][36]. These quantitative measurements directly link the protein to rare earth element accumulation in bacterial cells[21][28][36].

Genetic Evidence

Transposon mutagenesis studies provide genetic evidence for functions of lanthanide acquisition proteins in the same genomic region as MexAM1_META1p4129[40]. Transposon insertions in the lut (lanthanide uptake) cluster, which shares genetic organization and function with the mll cluster, revealed that multiple genes in this region are required for lanthanide-dependent growth when combined with loss of the mxaF calcium-dependent methanol dehydrogenase[40]. These loss-of-function genetic studies establish that the entire region, including the TBDT component, functions in lanthanide acquisition[40].

Second-site suppressor mutations provide interesting genetic evidence that cells can acclimate to loss of certain components of the lanthanide acquisition system[40]. Strains lacking the ABC transporter genes lutE or lutF showed no initial growth on methanol with lanthanides, but after 150-200 hours, second-site suppressor mutations arose that allowed slow growth, suggesting alternative lanthanide transport routes exist but operate with reduced efficiency[40]. This genetic evidence indicates that MexAM1_META1p4129 (the TBDT) represents one route among potentially multiple lanthanide acquisition mechanisms, though the primary route[40].

Complementation studies provide direct genetic evidence for protein function through restoration of wild-type phenotypes[40]. Strains deleted for individual transport cluster genes were complemented by expression of the respective gene on a plasmid, and growth similar to wild-type strain was restored in each case, confirming that each gene's product is required for full lanthanide transport function[40].

Expression Evidence

Tissue and condition-specific expression patterns provide evidence for the biological context of MexAM1_META1p4129 function[21][28][44]. The gene is dramatically upregulated specifically under lanthanide limitation (as demonstrated by the 32-fold increase with poorly soluble Nd₂O₃ compared to soluble NdCl₃) but not under other stress conditions tested[21][28][44]. This highly specific expression pattern indicates that the protein responds to a particular environmental signal—lanthanide bioavailability—rather than general stress responses[21][28][44].

Operon organization and co-expression provide evidence for functional relationships within the mll gene cluster[21][28][44]. The gene is organized within an operon containing mluA, mluR (anti-sigma factor), and mluI (sigma factor), and expression of this entire unit is coordinately upregulated, suggesting functional coupling among these components[21][28][44]. Co-expression of genes encoding functionally related proteins (receptor, anti-sigma factor, and sigma factor) strongly supports their participation in a common biological process[21][28][44].

Two-component regulator analysis provides expression evidence for upstream regulatory control of mluA[45][48][57][60]. The expression of lanM (lanmodulin, a lanthanide-binding protein) is regulated by mxcQE (a two-component regulator for the calcium-dependent methanol dehydrogenase MxaF) and tonB_Ln (a TonB-dependent receptor for lanthanides), providing evidence that lanthanide acquisition systems are subject to complex transcriptional regulation[45][48][57][60]. By analogy, mluA expression likely undergoes similar regulatory control[45][48][57][60].

Biochemical and Biophysical Evidence

Protein structure prediction and domain analysis provide bioinformatic evidence for protein function[21][28][44]. Sequence analysis reveals that MexAM1_META1p4129 contains the characteristic features of TonB-dependent receptors: an N-terminal signal peptide, a large extracellular loop forming the putative substrate binding site, and a 22-stranded β-barrel transmembrane domain[21][28][44]. The presence of these conserved structural features provides strong inferential support for TBDT function[21][28][44].

Metallophore structure elucidation provides biochemical context for the substrate specificity of MexAM1_META1p4129[21][28][44]. The recently characterized methylolanthanin (MLL) structure, with its unique 4-hydroxybenzoate moiety, represents a novel substrate class for bacterial uptake systems[21][28][44]. The high sequence conservation between MexAM1_META1p4129 and characterized metallophore receptors, combined with the specific upregulation of the mll gene cluster in response to poorly soluble lanthanides, establishes MLL as the likely substrate for this receptor[21][28][44].

Evidence Quality Assessment

The overall evidence for MexAM1_META1p4129 function represents a strong integration of multiple independent lines of evidence, providing robust support for GO annotation:

Experimental Evidence Types: Transcriptomic (direct), growth phenotype (genetic), quantitative lanthanide bioaccumulation (physical), complementation (genetic), and sequence homology (comparative) evidence all support the protein's function in lanthanide-dependent metallophore uptake. This diversity of evidence types provides high confidence in functional assignments.

Strength of Statistical Support: The differential expression q-value of 3.70E-84 represents extraordinarily strong statistical support, far exceeding typical significance thresholds (typically p < 0.05 or q < 0.05), indicating that the gene's upregulation in response to lanthanide limitation is not due to random variation.

Reproducibility and Validation: The lanthanide bioaccumulation measurements using ICP-MS represent quantitative, reproducible measurements of a physical phenotype directly linked to protein function. Multiple independent studies confirm these observations.

Specificity of Response: The gene's selective upregulation specifically in response to poorly soluble lanthanides (not soluble lanthanides) indicates a specific environmental sensing function rather than a generic stress response.

Disease and Phenotype Associations

Bacterial Phenotype Associations

Loss-of-function phenotypes provide direct evidence for the biological importance of MexAM1_META1p4129[28][36][40]. Strains with deletions in the mll biosynthetic gene cluster exhibit severely impaired growth on methanol in the presence of lanthanides as cofactors for methanol dehydrogenase[28][36][40]. This phenotype directly links the protein to methylotrophic metabolism in lanthanide-dependent bacteria[28][36][40]. Additionally, deletion of mll genes results in dramatic reduction in lanthanide bioaccumulation (1.8-fold reduction in ICP-MS measurements), indicating that this accumulation is not a passive process but actively requires the encoded transport machinery[28][36][40].

Gain-of-function phenotypes provide complementary evidence for protein function[28][36]. Overexpression of the entire mll biosynthetic gene cluster results in 3.5-fold increase in lanthanide bioaccumulation and adsorption compared to wild-type levels, demonstrating that increased expression of the uptake system directly enhances lanthanide acquisition[28][36]. This dose-dependent response indicates that MexAM1_META1p4129 functions rate-limitingly in the lanthanide acquisition pathway[28][36].

Phenotypic acclimation and suppressor mutations reveal interesting compensatory mechanisms in the absence of MexAM1_META1p4129 and related proteins[40]. Strains lacking key components of the lanthanide transport system initially show severe growth defects on methanol with lanthanides, but after extended culture periods (150-200 hours), suppressor mutations arise that allow slow growth, approximately 88% slower than wild-type[40]. This observation suggests that while MexAM1_META1p4129 represents the primary lanthanide acquisition route, alternative, less-efficient routes exist in this organism[40].

Cellular Phenotypes Associated with Lanthanide Excess

Lanthanide toxicity and stress responses represent important phenotypic associations revealed through studies of lanM (lanmodulin) deletion mutants, which shed light on consequences of unregulated lanthanide accumulation[45][48][57][60]. When lanM deletion strains were exposed to lanthanum, they exhibited aggregating phenotypes, cell membrane impairment (evidenced by electron microscopy), lanthanum deposition in the periplasm, and differential expression of proteins involved in membrane integrity and phosphate starvation[45][48][57][60]. These phenotypes demonstrate that excessive lanthanide accumulation, which would occur without proper regulation of the MexAM1_META1p4129-mediated uptake, is toxic to the cell[45][48][57][60].

Membrane disruption and phosphate homeostasis represent critical phenotypic consequences of lanthanide excess[37][45][48][57][60]. Lanthanides exhibit high affinity for phosphate ions, and excessive lanthanide accumulation can disrupt phosphate-dependent cellular processes and sequester phosphate into insoluble lanthanide-phosphate complexes[37][45][48][57][60]. The phenotypic observation of phosphate starvation responses in lanM deletion mutants exposed to lanthanides suggests that proper lanthanide homeostasis, including controlled uptake via MexAM1_META1p4129, is essential for maintaining cellular phosphate balance[37][45][48][57][60].

Ecological and Evolutionary Phenotypes

Distribution in environmental bacteria represents an important phenotypic association linking MexAM1_META1p4129 to ecological niches[21][28][36][44]. The mll gene cluster is conserved across diverse Methylobacterium and Methylorubrum species, with phylogenetic analysis revealing that the majority of Methylorubrum species forming a single clade contain homologs of the mll locus, suggesting that lanthanide acquisition systems represent important adaptations to specific environmental niches, likely soil environments where lanthanide availability is variable[21][28][36][44].

Plant-associated and phyllosphere phenotypes suggest ecological roles for MexAM1_META1p4129[21][25][28][36][44]. Methylobacterium species are abundant in the plant phyllosphere and function as plant growth-promoting bacteria, and some evidence suggests that lanthanide-dependent methanol metabolism may contribute to plant interactions[21][25][28][36][44]. The presence of lanthanide acquisition systems in plant-associated bacteria suggests these proteins may play ecological roles in plant-microbe interactions[21][25][28][36][44].

Protein Structure and Domains

Structural Features and Motifs

β-barrel transmembrane domain represents the major structural feature of MexAM1_META1p4129[26][32][35]. As a TonB-dependent receptor, this protein contains the characteristic 22-stranded β-barrel structure that spans the outer membrane and forms a selective pore for substrate translocation[26][32][35]. The β-barrel provides structural stability and creates the substrate binding pocket within the extracellular space and β-barrel lumen[26][32][35]. This domain represents the conserved scaffold structure for all TonB-dependent receptors, enabling the diverse substrate specificities observed across this protein family[26][32][35].

Extracellular loops and substrate binding site represent functionally critical structural elements[21][28][44]. TonB-dependent receptors typically contain large extracellular loops that form the initial substrate binding site and interact with periplasmic carrier proteins[26][32][35]. In the case of MexAM1_META1p4129, these loops likely provide the binding site for methylolanthanin complexes transported from the environment, though the exact structure of this binding site requires crystallographic determination[21][28][44].

N-terminal signaling domain (SD) represents a specialized structural feature of TBDTs functioning as transducers[21][26][28][32][44]. Recent structural studies of homologous TBDTs have revealed that this signaling domain undergoes conformational changes upon substrate binding, with rotation of this domain occurring in response to ligand binding at the extracellular site[26][32]. The signaling domain extends into the periplasm and interacts with the anti-sigma factor MluR, facilitating the signal transduction cascade that leads to sigma factor release[26][32][44].

Signal peptide and outer membrane targeting sequence represent the N-terminal targeting signals[26][32][35]. Like all outer membrane proteins, MexAM1_META1p4129 contains an N-terminal signal peptide that directs the nascent polypeptide to the Sec pathway and the outer membrane insertion machinery[26][32][35]. This signal peptide is typically cleaved during outer membrane insertion, leaving the mature protein with an N-terminal methionine[26][32][35].

Plug domain represents a structural feature found in many TonB-dependent receptors, though its presence in MexAM1_META1p4129 requires further characterization[26][32][35]. Some TonB-dependent receptors contain an N-terminal globular "plug" domain that occludes the β-barrel lumen in the inactive state, requiring conformational rearrangement for substrate access[26][32][35]. If present in MexAM1_META1p4129, this domain would regulate substrate access and prevent non-specific small molecule translocation[26][32][35].

Post-Translational Modifications

Absence of lipidation distinguishes MexAM1_META1p4129 from some membrane proteins[26][32][35]. TonB-dependent receptors typically lack lipid modifications, instead relying on the β-barrel structure and outer membrane lipid interactions for stability and localization[26][32][35]. The protein likely remains stably inserted in the outer membrane through hydrophobic interactions between the β-barrel and the lipid bilayer throughout its cellular lifetime[26][32][35].

Potential proteolytic processing may occur at specific sites during the signal transduction process[26][32]. Recent evidence for homologous TBDTs suggests that the signaling domain may undergo proteolytic cleavage during the cell-surface signaling cascade, likely by the periplasmic serine protease DegP or related proteases[26][32]. If this proteolytic processing occurs for MexAM1_META1p4129, it would represent a mechanism for signal amplification or regulatory control of the signaling cascade[26][32].

Oxidation states of amino acids may play functional roles, as suggested by studies of homologous proteins[26][32]. The extracellular domain of TonB-dependent receptors may contain cysteine residues in specific oxidation states that affect substrate binding or conformational transitions, though this requires characterization specific to MexAM1_META1p4129[26][32].

Functional Domains and Their Roles

Substrate binding domain encompasses the extracellular portion and likely includes structured loops forming the methylolanthanin binding site[21][28][44]. The specificity of binding for lanthanide-metallophore complexes versus other potential ligands depends on the precise three-dimensional arrangement of amino acid side chains within this domain, which remains to be structurally characterized[21][28][44].

Energy transduction domain comprises the TonB box region that interacts with the TonB protein to couple proton-motive force to substrate translocation[26][32][35]. This short linear motif (typically around 9-14 amino acids) provides the interaction surface with the periplasmic TonB protein and enables energy coupling from the proton gradient to mechanical conformational changes in the receptor[26][32][35].

Signal transduction domain encompasses the N-terminal region that interacts with the anti-sigma factor MluR to propagate the metal acquisition signal to the transcriptional regulatory machinery[26][32][44]. This domain must undergo conformational changes that enable stronger interaction with MluR upon substrate binding, though the mechanistic details require further investigation[26][32][44].

Evolutionary Conservation and Orthologs

Ortholog Distribution and Conservation

Conservation across Methylobacterium and Methylorubrum genera represents a notable evolutionary conservation pattern[21][28][36][44]. Phylogenetic analysis revealed that the majority of Methylorubrum species in a primary clade contain homologs of the mll locus including the TBDT component, as well as Methylobacterium currus TP3 and Methylobacterium aquaticum BG2[21][28][36][44]. This phylogenetic distribution suggests that lanthanide acquisition systems, including MexAM1_META1p4129-like proteins, were present in the common ancestor of these genera and have been retained across multiple speciation events[21][28][36][44].

Absence in distantly related methylotrophs provides informative evolutionary context[21][28][36][44]. The mll locus was not found in any of the 85 genomes previously identified as containing XoxF homologs (lanthanide-dependent methanol dehydrogenases), indicating that while lanthanide-dependent metabolism is widespread, specific metallophore-based acquisition systems like the mll cluster represent a more specialized adaptation[21][28][36][44]. This suggests that multiple independent routes for lanthanide acquisition have evolved in different bacterial groups[21][28][36][44].

Homology to established siderophore transporters provides evolutionary and functional context for MexAM1_META1p4129[21][28][44]. The protein shows 62% amino acid sequence identity to Rpa1_2620 from the rhodopetrobactin biosynthetic gene cluster in Rhodopseudomonas palustris, a well-characterized TBDT for siderophore transport[21][28][44]. This homology relationship indicates that MexAM1_META1p4129 shares a common evolutionary origin with established metal acquisition systems and has likely undergone functional divergence to recognize lanthanide-metallophore complexes rather than traditional siderophores[21][28][44].

Paralog Relationships and Functional Divergence

Genetic organization similarities to siderophore uptake systems provide evolutionary insights[21][28][44]. The mll gene cluster exhibits genetic organization highly similar to characterized siderophore biosynthetic clusters, with the TBDT (mluA), anti-sigma factor (mluR), sigma factor (mluI), and biosynthetic genes organized in functional units reminiscent of siderophore uptake and biosynthesis operons[21][28][44]. This conservation of genetic organization across different metal acquisition systems suggests evolutionary conservation of regulatory strategies[21][28][44].

Hybrid NRPS/NIS biosynthetic pathway relationships provide context for the evolutionary origin of the mll cluster[44][50]. While the mll gene cluster represents an NRPS-independent siderophore (NIS) biosynthetic pathway, some metallophore biosynthetic clusters employ hybrid NRPS/NIS pathways, and the evolutionary relationship between these systems remains to be fully characterized[44][50]. The conservation of certain biosynthetic modules across different metallophore types suggests common evolutionary ancestry[44][50].

Orthologs with Characterized Function

LutA and related periplasmic binding proteins represent functionally related orthologs in lanthanide transport systems[40][45][48][57][60]. In the Methylorubrum extorquens AM1 lut (lanthanide uptake) cluster, LutA functions as a periplasmic binding protein that traffics lanthanophore-bound lanthanides through the periplasm to the inner membrane transporter system[40]. By analogy, the periplasmic domain of MexAM1_META1p4129 likely functions in similar manner to deliver lanthanide-metallophore complexes to downstream transport machinery[40][45][48][57][60].

TonB-dependent receptors from pathogenic bacteria represent extensively characterized orthologs[26][32][35][51][54]. The FecA receptor from Escherichia coli and the FpvA and PvdS receptors from Pseudomonas aeruginosa represent the best-characterized TonB-dependent receptors, with crystal structures and detailed functional studies available[26][32][35][51][54]. These characterized orthologs provide structural and mechanistic context for understanding MexAM1_META1p4129 function[26][32][35][51][54].

Sigma factor ECF systems represent conserved orthologs across multiple bacterial genera[26][32][35][51][54]. The Fox, Fiu, and Iut cell-surface signaling systems from Pseudomonas aeruginosa represent the paradigmatic ECF sigma factor systems that function through TBDT-mediated signaling similar to the MluA/MluR/MluI system, providing functional and mechanistic parallels[26][32][35][51][54].

Essential Gene Status and Conditional Essentiality

Conditional essentiality in lanthanide-rich environments characterizes MexAM1_META1p4129[28][36][40]. While the gene is not essential when bacteria are cultured on calcium-dependent methanol dehydrogenase substrates or in traditional minimal media without lanthanides, it becomes functionally important when lanthanides are the only available cofactors for methanol oxidation[28][36][40]. This conditional essentiality reflects the specific environmental niches where lanthanide-dependent metabolism represents a competitive advantage[28][36][40].

Redundancy with alternative transport systems indicates that MexAM1_META1p4129 is not absolutely essential even in lanthanide-dependent growth conditions[40]. Strains lacking key components of the lanthanide transport system can develop suppressor mutations allowing slow growth, suggesting that alternative, less-efficient transport mechanisms exist[40]. This functional redundancy may reflect the evolution of multiple overlapping lanthanide acquisition systems to ensure metabolic flexibility in variable environments[40].

Role in competitive fitness likely represents the primary selective advantage of maintaining MexAM1_META1p4129 function[28][36][40]. In mixed microbial communities where lanthanide availability limits methylotrophic growth rates, strains with efficient lanthanide acquisition systems (including MexAM1_META1p4129) would outcompete strains lacking these systems, leading to maintenance of these genes through natural selection[28][36][40].

Specific GO Term Recommendations for Annotation

Based on the comprehensive analysis of experimental evidence, literature support, and functional characterization, the following Gene Ontology terms are recommended for annotation of MexAM1_META1p4129 (mluA):

Molecular Function Annotations (GO:MF)

1. TonB-dependent siderophore receptor activity (GO:0042883)
- Rationale: While this GO term technically refers to siderophores, metallophore uptake systems function through identical mechanisms to siderophore uptake, and no more specific GO term exists for metallophore receptors. The protein demonstrates the molecular function of outer membrane substrate recognition and energy-dependent translocation characteristic of this class.
- Evidence Level: Direct experimental evidence from differential expression, homology to characterized TBDTs (62% identity to Rpa1_2620), and phenotypic analysis of bioaccumulation and complementation studies.
- PMID Supporting Evidence: [21][28][32][36][44]

2. Lanthanide ion binding (GO:0070730)
- Rationale: The protein's primary function involves binding of lanthanide-metallophore complexes. While this GO term typically describes direct metal binding proteins (like lanmodulin), MexAM1_META1p4129 binds lanthanides indirectly through the metallophore ligand. This term captures the lanthanide specificity of the protein.
- Evidence Level: Direct experimental evidence from lanthanide bioaccumulation studies showing gene-dependent lanthanide accumulation, and conditional upregulation with lanthanide oxide.
- PMID Supporting Evidence: [21][28][36]

Alternative Term: Outer membrane receptor activity (GO:0015459) - This broader term captures the molecular function without assuming the specific ligand type.
- Evidence Level: Direct evidence from transmembrane domain predictions, homology to characterized outer membrane receptors, and outer membrane localization.
- PMID Supporting Evidence: [26][32][35]

3. Signal transducer activity (GO:0004871)
- Rationale: As a TonB-dependent transducer (TBDT), the protein functions in cell-surface signaling, transmitting signals from substrate binding at the outer membrane to regulation of inner membrane anti-sigma factors. This captures the signaling function of the protein beyond simple substrate transport.
- Evidence Level: Direct evidence from genetic organization (co-expression with anti-sigma factor mluR and sigma factor mluI) and homology to characterized signaling TBDTs.
- PMID Supporting Evidence: [21][26][28][32][35][44]

Cellular Component Annotations (GO:CC)

1. Outer membrane (GO:0019867)
- Rationale: This is the primary subcellular localization of MexAM1_META1p4129. The protein is an integral outer membrane protein predicted to contain the characteristic β-barrel transmembrane domain.
- Evidence Level: Strong indirect evidence from sequence analysis (22-stranded β-barrel prediction, signal peptide for outer membrane targeting) and homology to characterized outer membrane receptors. No direct experimental localization data is provided in the literature, but the prediction is extremely high confidence.
- PMID Supporting Evidence: [21][26][28][32][35][44]

2. Cell outer membrane (GO:0009279)
- Rationale: A more specific GO term designating the gram-negative bacterial outer membrane location.
- Evidence Level: Same as above.
- PMID Supporting Evidence: [21][26][28][32][35][44]

3. Outer membrane protein complex (GO:0032991, component_part GO:0044422)
- Rationale: The protein functions as part of the TonB-ABC transporter complex for lanthanide acquisition. This term captures its participation in multiprotein complexes.
- Evidence Level: Direct genetic evidence from complementation studies showing requirement for proper function of both the TBDT and downstream ABC transporter components.
- PMID Supporting Evidence: [28][33][40]

4. Cell-surface signaling complex (GO:specific complex identifier if available, otherwise "protein complex" GO:0043234)
- Rationale: The protein functions as a component of the CSS (cell-surface signaling) complex involving MluR (anti-sigma factor) and MluI (sigma factor).
- Evidence Level: Strong genetic evidence from operon organization, co-expression patterns, and homology to characterized CSS systems.
- PMID Supporting Evidence: [21][26][28][32][35][44][51]

Biological Process Annotations (GO:BP)

1. Lanthanide ion transport (GO:specialized transport process)
- Rationale: The primary biological function of MexAM1_META1p4129 is facilitating the translocation of lanthanide-metallophore complexes across the outer membrane. This GO term captures the core biological process.
- Evidence Level: Direct experimental evidence from quantitative lanthanide bioaccumulation measurements using ICP-MS, showing 1.8-fold reduction in mluA deletion strains and 3.5-fold increase in overexpression strains.
- PMID Supporting Evidence: [21][28][36]

2. Regulation of methanol oxidation (GO:specific pathway if available, otherwise "metabolic process regulation" GO:0019222)
- Rationale: By enabling lanthanide-dependent methanol dehydrogenase (XoxF) expression and function, MexAM1_META1p4129 regulates a critical metabolic pathway. The protein's role in lanthanide acquisition directly controls whether methanol oxidation can proceed via the lanthanide-dependent route.
- Evidence Level: Indirect but strong evidence linking lanthanide acquisition (direct function of the protein) to methanol dehydrogenase expression (established in literature), with the protein playing the rate-limiting role in lanthanide acquisition.
- PMID Supporting Evidence: [21][28][33][36][43][44][46][56][59]

3. Response to nutrient stress (GO:0009266)
- Rationale: The gene is dramatically upregulated (32-fold) specifically in response to poorly soluble lanthanide sources, indicating the protein participates in adaptive responses to lanthanide limitation.
- Evidence Level: Direct experimental evidence from transcriptomic analysis showing highly specific and strong upregulation in response to lanthanide-limited conditions.
- PMID Supporting Evidence: [21][28][44]

4. Intracellular metal ion homeostasis (GO:0051560) or Metal ion homeostasis (GO:0055075)
- Rationale: While MexAM1_META1p4129 promotes lanthanide uptake, this uptake must be balanced with efflux and storage mechanisms to maintain healthy cellular lanthanide levels. Dysregulation of this protein (without coordinating regulatory mechanisms) leads to lanthanide toxicity.
- Evidence Level: Indirect evidence from studies of lanM deletion mutants, which show that excessive lanthanide accumulation (resulting from unregulated uptake) causes cellular dysfunction including membrane damage and phosphate starvation.
- PMID Supporting Evidence: [37][45][48][57][60]

5. Cellular response to metal ion stress (GO:0071276)
- Rationale: The protein participates in the cell's adaptation to varying lanthanide availability, representing a cellular stress response.
- Evidence Level: Direct evidence from the highly specific upregulation of the mll cluster under lanthanide-limited conditions.
- PMID Supporting Evidence: [21][28][44]

6. Signal transduction (GO:0007165)
- Rationale: As a TonB-dependent transducer, the protein functions in the cell-surface signaling cascade that transmits the lanthanide availability signal to the transcriptional regulatory machinery.
- Evidence Level: Direct genetic evidence from the organized co-expression of the TBDT (mluA), anti-sigma factor (mluR), and sigma factor (mluI) within the same operon, and homology to characterized cell-surface signaling systems.
- PMID Supporting Evidence: [21][26][28][32][35][44][51]

Additional Annotations Based on Biological Context

7. Methylotrophy (GO:0019380) or Methylotrophic processes
- Rationale: By enabling lanthanide-dependent methanol dehydrogenase function, the protein participates in the broader methylotrophic metabolic process.
- Evidence Level: Indirect but strong evidence linking the protein to lanthanide acquisition, which enables XoxF-dependent methanol oxidation.
- PMID Supporting Evidence: [19][28][33][36][43][44][46][56][59]

8. Rare earth element bioaccumulation (GO:specialized bioaccumulation process if available, otherwise "small molecule metabolic process" GO:0044281)
- Rationale: The protein's role in rare earth element transport and accumulation makes it central to lanthanide bioaccumulation processes, with potential biotechnological applications.
- Evidence Level: Direct quantitative evidence from ICP-MS measurements of lanthanide bioaccumulation.
- PMID Supporting Evidence: [21][28][36]

Annotation with Appropriate Quantitative Evidence

Conclusion

The gene MexAM1_META1p4129 (designated mluA, encoding methylolanthanin uptake receptor A) represents a comprehensively characterized example of a TonB-dependent transducer (TBDT) protein functioning in lanthanide acquisition and cell-surface signaling. The integrated analysis of transcriptomic data, genetic evidence, quantitative bioaccumulation measurements, comparative genomics, and functional homology studies provides robust support for accurate Gene Ontology annotation of this protein.

Key Functional Conclusions

Molecular Function: The protein primarily functions as an outer membrane receptor for lanthanide-metallophore complexes, exhibiting high specificity for methylolanthanin-bound lanthanides including lanthanum, cerium, praseodymium, neodymium, and samarium. The protein demonstrates secondary function as a signal transducer in cell-surface signaling cascades, transmitting lanthanide availability signals to transcriptional regulatory machinery.

Cellular Component: The protein localizes exclusively to the outer membrane where it functions as an integral component of the TonB-ABC transporter complex and the cell-surface signaling complex involving the anti-sigma factor MluR and sigma factor MluI.

Biological Process: The protein participates in multiple interconnected biological processes: (1) lanthanide ion transport and bioaccumulation, (2) regulation of lanthanide-dependent methanol dehydrogenase expression and function, (3) response to nutrient stress (specifically lanthanide limitation), (4) metal ion homeostasis requiring coordination with efflux and storage systems, and (5) signal transduction enabling adaptive gene expression responses.

Evidence Integration

The strength of annotation support derives from the exceptional integration of multiple independent lines of evidence. The transcriptomic data showing 32-fold upregulation (q-value: 3.70E-84) provides extraordinarily strong statistical support for differential expression. The quantitative ICP-MS measurements of lanthanide bioaccumulation (showing 1.8-fold reduction in deletion strains and 3.5-fold increase in overexpression strains) provide direct physical evidence linking the protein to lanthanide accumulation. The complementation studies and growth phenotypes provide genetic evidence. The 62% sequence identity to characterized metallophore transporters (Rpa1_2620) provides comparative functional evidence. The organized operon structure with anti-sigma factor and sigma factor genes provides mechanistic context for signal transduction.

Recommendations for Annotation Curation

For comprehensive and accurate Gene Ontology annotation, the following term hierarchy is recommended:

Primary Molecular Function Terms: GO:0042883 (TonB-dependent siderophore receptor activity) or GO:0015459 (outer membrane receptor activity)

Cellular Component Terms: GO:0019867 (outer membrane), GO:0044422 (component of outer membrane protein complex), GO:0043234 (protein complex participation in cell-surface signaling)

Biological Process Terms: GO:specialized terms for lanthanide transport, GO:0007165 (signal transduction), GO:0009266 (response to nutrient stress), GO:0019380 (methylotrophy) with consideration of indirect linkages

All annotations should include appropriate evidence codes (IMP for mutant phenotypes, IEP for ortholog characterization, IGC for genomic context, ISS for sequence similarity) with citations to the primary literature supporting each functional claim.

The exceptional quality and diversity of experimental evidence for MexAM1_META1p4129 function makes this an excellent candidate for curator-selected annotation with high confidence scores in community annotation systems.

Citations

1. https://pmc.ncbi.nlm.nih.gov/articles/PMC1422150/
2. https://pmc.ncbi.nlm.nih.gov/articles/PMC9727139/
3. https://pubmed.ncbi.nlm.nih.gov/40429863/
4. https://pmc.ncbi.nlm.nih.gov/articles/PMC10050551/
5. https://journals.asm.org/doi/10.1128/spectrum.00359-25
6. https://networks.systemsbiology.net/syntrophy/genes/mex/6829462
7. https://research.aalto.fi/files/87154591/CHEM_Bao_et_al_Efficient_CRISPR_2022_ACS_Synthetic_Biology.pdf
8. https://pmc.ncbi.nlm.nih.gov/articles/PMC10029344/
9. https://omicsdatascience.icm.edu.pl/wp-content/uploads/2024/12/Bezshapkin.V_Comprehensive-functional-annotation-of-metagenomes-and-microbial-genomes-using-a-deep-learning-based-2023.pdf
10. https://pubmed.ncbi.nlm.nih.gov/38702257/
11. https://pmc.ncbi.nlm.nih.gov/articles/PMC4158751/
12. https://pmc.ncbi.nlm.nih.gov/articles/PMC11137889/
13. https://pubmed.ncbi.nlm.nih.gov/15466049/
14. https://pmc.ncbi.nlm.nih.gov/articles/PMC4253957/
15. https://www.nature.com/articles/ismej201191
16. https://journals.asm.org/doi/10.1128/jb.186.20.6956-6969.2004
17. https://pmc.ncbi.nlm.nih.gov/articles/PMC7419163/
18. https://pmc.ncbi.nlm.nih.gov/articles/PMC111365/
19. https://pubmed.ncbi.nlm.nih.gov/27286964/
20. https://pmc.ncbi.nlm.nih.gov/articles/PMC3087415/
21. https://pmc.ncbi.nlm.nih.gov/articles/PMC11317620/
22. https://pubs.acs.org/doi/10.1021/acssynbio.1c00322
23. https://www.nature.com/articles/s42003-025-07583-1
24. https://deepblue.lib.umich.edu/bitstream/handle/2027.42/199194/junwony_1.pdf?sequence=1&isAllowed=y
25. https://pmc.ncbi.nlm.nih.gov/articles/PMC9301485/
26. https://pmc.ncbi.nlm.nih.gov/articles/PMC11637429/
27. https://pmc.ncbi.nlm.nih.gov/articles/PMC4902934/
28. https://www.pnas.org/doi/10.1073/pnas.2322096121
29. https://sfamjournals.onlinelibrary.wiley.com/doi/10.1111/1751-7915.14503
30. https://pmc.ncbi.nlm.nih.gov/articles/PMC7849352/
31. https://chemrxiv.org/engage/chemrxiv/article-details/673c8ae15a82cea2faa8f2db
32. https://journals.plos.org/plosbiology/article?id=10.1371%2Fjournal.pbio.3002920
33. https://www.nature.com/articles/s41598-020-69401-4
34. https://pmc.ncbi.nlm.nih.gov/articles/PMC11845868/
35. https://pmc.ncbi.nlm.nih.gov/articles/PMC4701049/
36. https://pubmed.ncbi.nlm.nih.gov/39078674/
37. https://pmc.ncbi.nlm.nih.gov/articles/PMC8106350/
38. https://journals.sagepub.com/doi/10.1177/1934578X241271701
39. https://pmc.ncbi.nlm.nih.gov/articles/PMC9045217/
40. https://pmc.ncbi.nlm.nih.gov/articles/PMC7391723/
41. https://pubs.acs.org/doi/10.1021/acs.jnatprod.4c00911
42. https://pmc.ncbi.nlm.nih.gov/articles/PMC9148527/
43. https://pmc.ncbi.nlm.nih.gov/articles/PMC4859578/
44. https://pmc.ncbi.nlm.nih.gov/articles/PMC9260416/
45. https://www.nature.com/articles/s41598-019-41043-1
46. https://pmc.ncbi.nlm.nih.gov/articles/PMC5895095/
47. https://pubmed.ncbi.nlm.nih.gov/35814700/
48. https://pmc.ncbi.nlm.nih.gov/articles/PMC8991815/
49. https://pmc.ncbi.nlm.nih.gov/articles/PMC10700112/
50. https://pmc.ncbi.nlm.nih.gov/articles/PMC6851896/
51. https://pmc.ncbi.nlm.nih.gov/articles/PMC116012/
52. https://www.tandfonline.com/doi/full/10.1080/10409238.2025.2476476?af=R
53. https://www.nature.com/articles/s41598-020-60197-x
54. https://pmc.ncbi.nlm.nih.gov/articles/PMC5887029/
55. https://pubmed.ncbi.nlm.nih.gov/26833413/
56. https://pubmed.ncbi.nlm.nih.gov/24553424/