BRIEF RESEARCH REPORT article

Front. Mar. Sci., 29 August 2025

Sec. Marine Molecular Biology and Ecology

Volume 12 - 2025 | https://doi.org/10.3389/fmars.2025.1661535

Determining the taxonomic and functional profile of marine bacterial copper systems involved in marine early copper surface colonization

  • 1. Department of Biological Sciences, University of Maryland Baltimore County, Baltimore, MD, United States

  • 2. Precise Systems, Inc, Lexington Park, MD, United States

  • 3. Chemistry Division, U.S. Naval Research Laboratory, Washington, DC, United States

Abstract

Copper-based antimicrobial paints are frequently used to inhibit the biofouling of marine vessels. However, some bacterial species can overcome this copper challenge and colonize the surfaces. The early adherent bacterial population of marine vessels plays an important role because of its ability to produce extracellular polymeric substances (EPSs), forming a thin layer of organic matter that traps nutrients from the water and protects other colonizers by blocking the toxic antifouling (AF) coatings. It is of interest to study the factors that drive the initial colonization of copper surfaces. We used a metagenomic sequencing approach to analyze the microbial diversity and potential functional abilities of early biofilm communities on copper surfaces and discovered enriched copper-specific functional traits in early colonizers compared to the bacterial community of the surrounding seawater. The functional analysis of early biofilm colonizers on copper surfaces in marine coastal environments provides insights into molecular mechanisms that support biofilm formation on copper surfaces. The taxa associated with copper resistance traits were found to dominate initial microbial communities on copper surfaces, allowing these organisms to survive on copper surfaces despite copper toxicity. Our analysis reveals the dominance of the genera Allomuricauda and Ruegeria, carrying several copper resistance genes, as the early colonizers of copper surfaces.

Introduction

In the marine environment, biofilms provide shelter and sustenance to the microorganisms, which receive protection from environmental changes, substrates for growth at relatively higher concentrations than found in the planktonic environment, and conditions favorable to genetic exchange and survivability (). Biofilm formation creates significant challenges in the marine industry (), causing an increase in drag forces that require additional shaft power (). Microbial attachment to vessel surfaces leads to higher maintenance costs and efforts, and risks introducing previously unknown biomes into new ecosystems (). The colonization of submerged surfaces is affected by various factors, including environmental characteristics, substrate types, and the biological properties of colonizers (; ).

Biofouling starts with the adhesion of bacterial cells and the formation of biofilms; these bacterial biofilms alter the surface and enable subsequent attachment of organisms such as diatoms, algae, and barnacles (). To limit biofouling, copper is widely used as a biocide due to its antimicrobial properties (). However, various microbial species are able to colonize copper-treated surfaces. The early colonizing population plays a crucial role in this process, as it can create a barrier between the interaction of toxic biocides and subsequent colonizers (; ).

Previous studies of microbial colonization on specific substrates have demonstrated that substrate type significantly influences the composition of biofilm communities (; ). On copper surfaces, microbial communities may develop resilience to copper toxicity through increased extracellular matrix production, copper resistance and detoxification processes, and the activation of stress-response pathways (). However, there are significant knowledge gaps in terms of the mechanisms employed by early colonizers to establish themselves on copper surfaces. Identifying these mechanisms is critical to understanding how biofilms initiate and evolve under the selective pressures of copper toxicity.

Prior research on early marine copper surface colonization has mainly focused on profiling microbial communities within the biofilm environment (; ; ; ). The functional research into copper-associated biofilm is limited to predictive functional approaches that use 16S rRNA sequencing (). Functional studies of copper-associated biofilm have primarily focused on long-term seawater exposure. Studies on marine biofilms on copper surfaces demonstrate that long-term seawater exposure drives the selection of specific functional genes, especially those involved in membrane transport ().

Many bacteria have evolved distinct regulatory mechanisms to cope with high copper concentrations. These mechanisms include resistance-nodulation-cell-division (RND) efflux pumps, P-type ATPases, metallothioneins, cation diffusion facilitator (CDF) family copper transporters, and oxidation of copper ions (; ; ; ).

In this study, we used a shotgun metagenomics approach to identify early bacterial colonizers of copper surfaces and functional gene profiles from biofilm and seawater communities, with a focus on copper resistance pathways. Our analysis revealed that taxa contributing copper resistance traits dominate the early stage of biofilm. This suggests the significance of early colonizers in establishing the necessary conditions to mitigate copper toxicity and facilitate the adherence of other fouling organisms. Our work establishes the fundamental groundwork for future research comparing functional signatures from various copper-related surfaces, thereby developing strategies against marine biofouling.

Methods

Sample collection and processing

The copper panels (90–10 Cu-Ni alloy, 4 X 4 inches, 0.25 inch) were deployed for 7 days at Key West, Florida (October 2023). The copper panels were suspended in an ambient seawater, continuous flow-through trough (Supplementary Figure 1). At the time of deployment, the seawater temperature was 25°C, pH was 8.02, salinity was 36 and DO was 7.5 mg/L. These parameters remained relatively stable until sampling. A thin layer of biofilm was observed and scraped off the surface after 7 days of immersion. The seawater samples from the same area were collected by filtering 1 liter of seawater through a sterile 0.22 µm membrane filter. After collection, the biofilm and seawater were suspended in a nucleic acid stabilizer, DNA/RNA shield (Zymo Research). Samples were immediately frozen at -20°C until shipped.

DNA extraction and quality control

DNA from biofilm was extracted using a DNeasy Power biofilm kit (Qiagen), and seawater samples were extracted using a DNeasy Blood and Tissue kit (Qiagen). After extraction, DNA was quantified using a nanodrop spectrophotometer, and the amplification of bacterial 16S rRNA gene fragments was validated by polymerase chain reaction (PCR) using universal PCR primers ().

Quality control, library preparation, and sequencing

DNA concentration and fragment size were quantified using the Agilent Femto Pulse capillary system. The sequencing libraries were prepared following the HiFi SMRTbell library preparation protocol. The fragment lengths of the prepared libraries were validated using the Femto Pulse system, revealing that the fragment lengths of water samples were smaller than those of biofilm samples. To account for this size difference, water and biofilm libraries were processed at a 1:5 ratio, respectively, and sequencing was performed on a single SMRT cell. Libraries were sequenced on an 8M SMRTcell on a Sequel II platform using circular consensus sequencing protocol. After sequencing, the reads were demultiplexed to obtain high-quality HiFi FASTQ reads for all samples.

Taxonomic annotation of HiFi reads

The taxonomic classification and functional annotation of the reads were performed using established tools compiled by Pacific Biosciences in their GitHub repository for long-read metagenomic workflows (). The reads were analyzed utilizing the Diamond and MEGAN analysis pipeline ().

A local database of NCBI NR protein sequences was created by using the Diamond alignment tool (). The HiFi sequence reads were aligned with the compiled protein database using the Diamond BLASTX option. For this analysis, we used the default parameters provided by the Diamond tool, selecting the -top5 hits parameter. An alignment was considered significant if it was within the top 5 percent of the high-scoring hit range. Given the long-read data, we included the frameshift-aware and range-culling parameters for Diamond alignment.

The aligned sequences were processed for taxonomic binning using the MEGAN 7 software (). MEGAN utilizes NCBI reference taxonomy to generate and assign taxonomy to reads by applying an interval union Lowest Common Ancestor (LCA) parameter. It determines the taxon assignment for each read by analyzing sequence conservation throughout the taxonomic tree. We labeled a node using the total number of aligned bases, considering the length of long reads, as the alignment coverage can vary significantly among different reads ().

MEGAN default parameters were used for the analysis, except for the minimum support parameter, which was set to 0.01 (which provides optimal tradeoff between precision and recall for HiFi metagenomics data) (), and abundance, set to be reported as the number of aligned bases. Thus, a taxon was reported if at least 0.01% of the aligned bases were assigned to that node. If the minimum support value was not met, MEGAN assigned the reads to its lowest common ancestor where the support parameter value was satisfied. The total number of aligned bases was reported as count data. The relative abundance at each taxonomic category was calculated by normalizing counts to the read counts from the water samples (minimum number of counts among all samples).

The taxonomic assignments were stored as read-match archive (RMA) files, which were used for interactive analysis and visualization in the MEGAN user interface. Comparative analyses of aligned base counts at different taxonomic levels were performed to estimate the relative abundance of taxa across the two environments.

SEED subsystems of copper with MEGAN

The assignment of reads to functional categories was determined by their best-scoring match to a reference protein with a known function. SEED-based functional annotations were used to extract reads associated with copper. Copper-associated genes were identified as a part of the membrane transport subsystem (including sub-subsystems: copper transport & blue copper proteins, the copper transport system, and copper uptake system), and the stress response, virulence subsystem (sub-systems: copper homeostasis) of SEED. These copper-associated reads were used for copper gene-specific homology search.

Copper gene-specific functional profiling with protein homology search

Copper-associated genes and their regulators were identified by literature review and were used as the proteins of interest for a homology search (; ). Profile Hidden Markov models (HMMs) of these proteins of interest were obtained from the PFAM and TIGRFAM (; ). These HMMs were compiled into a local database using the HMMER suite.

The HMMER esl-translate tool was utilized to obtain the six-frame translation of the individual copper-associated reads. Translated reads were scanned against the compiled HMM database to identify protein homologs using the hmmscan tool. The copper-associated proteins of interest and their respective protein families are shown in Table 1. The search process was automated using a Python script. The hit count for copper-associated proteins was evaluated as a function of limiting e-value in 8 representative bacterial species to define an e-value threshold in a previous study (). The hits with e-values below 10E-30 were considered significant and thus classified as sequence homologs. The copper-associated protein hits from biofilm and water samples were analyzed to identify copper gene abundance levels. The abundance of specific copper genes across two samples was normalized with the same scaling factor as taxonomic abundance to account for the differences in the total sequence between biofilm and water.

Table 1

HMM Protein IDSystemSpeciesNCBI AccessionFunction
TIGR02044CueRE. coliNP_415020.1Copper efflux regulator
TIGR01480CueOE. coliNP_414665.1Multicopper oxidase
PcoAE. coliANH09828.1
CopAPseudomonas syringaeAQX42270.1
TIGR01511.1CopAE. coliNP_415017.1Copper efflux pump
TIGR01386CusSE. coliNP_415102.1two-component regulatory system
PcoSE. coli PlasmidCAA58530.1
CopSP. syringaeQBZ81958.1
TIGR01387CusRE. coliNP_415103.1two-component regulatory system
PcoRE. coli PlasmidCAA58529.1
CopRP. syringaeQBZ81959.1
TIGR00914.1CusAE. coliNP_415107.1transenvelope copper efflux system
TIGR01730.1CusBE. coliNP_415106.1transenvelope copper efflux system
TIGR01845CusCE. coliNP_415104.1transenvelope copper efflux system
PF11604CusFE. coliNP_415105.1copper chaperone
PF05275CopBP. syringaeAQX42189.1plasmid-borne copper resistance
PcoBE. coliANH09778.1
PF04234CopCP. syringaeAQX42188.1plasmid-borne copper resistance
PcoCE. coliANH09779.1
PF05425CopDP. syringaeAQX42268.1plasmid-borne copper resistance
PcoDE. coliANH09780.1
TIG02698CopYEnterococcus hiraeCAA86835.1CopY family regulator
PF02583CsoRM. tuberculosisNP_215482.1CsoR/RcnR family regulator
PF00440ComRE. coliWP_001336528TetR-like regulator
TIGR02937CorEMyxococcus xanthusABF90558.1copper-responsive ECF σ factor
PF00403CopZP. syringaeAQX41994.1copper chaperone

Copper-associated protein models used for homology search.

The Wilcoxon rank-sum test method enabled comparison between the overall distributions of copper gene abundance within biofilm samples and water samples (Non-parametric distribution verified by Shapiro-Wilk test p<0.05). A binomial test was used to detect copper genes, which showed statistically significant count differences between the two samples. The Benjamini-Hochberg correction method with a 5% false discovery rate threshold was applied to adjust p-values accounting for multiple hypothesis testing. The statistical analyses and data visualization processes were performed using R with the ggplot2 package for graphical representation and data manipulation.

Taxonomic assignment of copper genes

The reads that mapped to copper systems in SEED classification were assigned taxonomy with the DIAMOND and MEGAN pipeline to identify the bacterial genera carrying copper resistance genes. The taxa associated with the copper gene profile in the biofilm were then compared to the overall bacterial profile to assess the role of these bacterial genera carrying copper-resistant genes in the biofilm environment.

Results and discussion

Bacterial diversity analysis at phylum and genus taxonomic levels

Biofilm and water DNA were sequenced using PacBio circular consensus sequencing (CCS) protocol to obtain HiFi reads (Supplementary Table 1). Bacterial reads accounted for over 90% of the total sequences and were thus isolated from overall reads. To identify bacterial reads, sequences were first aligned to the NCBI database, which were then assigned taxonomy using MEGAN. Reads specific to bacteria were extracted for downstream analysis. Taxonomic composition within the two samples was assessed by calculating the relative abundance, which is defined as the proportion of aligned bases assigned to a specific taxon relative to the total bases in the dataset. The relative abundance of bacterial taxa was examined at the phylum and genus taxonomic levels.

Taxonomic analysis revealed distinct variations in microbial composition between the biofilm and seawater samples. At a phylum level, the biofilm and water samples were quite similar: Pseudomonadota dominated in both communities, followed by Bacteroidota. The overall distribution at a phylum level confirms previously described composition in marine biofilm environments, showing a greater abundance of phyla such as Planctomycetota, and Cyanobacteriota, known to be associated with organisms such as marine microalgae, macroalgae, sponges, protozoa, and other invertebrates (; ; ) (Supplementary Figure 2).

The biofilm microbial community mainly consisted of the genus Allomuricauda (11.09%) (phylum Bacteroidota), followed by members of Ruegeria (5.54%) (phylum Pseudomonadota). Both genera were also found in the seawater environment, but their abundance levels remained lower than those in the biofilm (Allomuricauda 2.52% and Ruegeria 1.78%). Additional dominant genera identified in the biofilm included Ulvibacterium (phylum Bacteroidota), Kordiimonas (phylum Pseudomonadota), Winogradskyella (phylum Bacteroidota), Congregibacter (phylum Pseudomonadota), Ekhidna (phylum Bacteroidota), Alkalimarinus (phylum Pseudomonadota), and Tateyamaria (phylum Pseudomonadota) (Figure 1). Some of these genera could not be detected in seawater at the used sequencing depth, suggesting they may be present at a relative abundance below the detection limit in the seawater.

Figure 1

In contrast, the dominant genus in seawater was Pelagibacter (phylum Pseudomonadota), representing approximately 10% of the community profile (Figure 1).

The genus Allomuricauda (recently reclassified as Flagellimonas) was previously identified as a key microbial genus associated with copper surface colonization (; ), and harbors copper-associated genes from both the Cus and Cop systems (). Allomuricauda is found in diverse marine habitats, demonstrating its ecological versatility and ability to survive under different environmental stress conditions (), including copper rich biofilm environment.

Ruegeria also exhibits copper-related adaptations, including a specific strain isolated from red algae, that harbors multiple chromosome-borne and plasmid-borne copper resistance operons, establishing its association with copper-enriched habitats (). Some species of Ruegeria, including R. pomeroyi DSS-3 (previously classified as Silicibacter pomeroyi) are moderate copiotrophs, as indicated by their genome analysis, thus suggesting higher abundance as seen in the biofilm environment (; ). Other dominant genera from the biofilm such as Kordimonas, Winogradskyella, and Congregibacter are known to secrete an EPS matrix, which supports their survival and adaptability in the biofilm (; ; ). In addition to copper genes, these traits may play a significant role in the dominance of these genera in the biofilm.

Previous studies have found that facultative anaerobic microbes dominate mature copper surface biofilm (). Our study on early-stage biofilm was dominated by representatives from primarily aerobic genera, including Allomuricauda, Ruegeria, Winogradskyella, Alkalimarinus Tateyamaria, and the obligate aerobe, Kordimonas. Other genera, such as Congregibacter and Ekhidna, are facultative anaerobes, allowing them to thrive in environments with oxygen gradients, such as biofilms. Our results highlight the dominance of taxa involved in aerobic metabolism in the initial biofilm environment.

Enrichment of copper gene panel

SEED-based functional annotations from MEGAN were used to extract reads associated with copper. The reads annotated to copper-associated genes in the SEED subsystems category were used to identify specific copper resistance proteins by homology search.

We performed a HMMER homology search using the profile HMM of proteins involved in bacterial copper homeostasis to identify its sequence homologs from the sequencing reads.

Bacterial copper regulatory mechanism involves genes responsible for efflux transport, detoxification, and chaperone functions. Copper homeostasis is primarily maintained with the chromosome-borne Copper efflux (Cue) system, which includes genes involved in cytoplasmic efflux and copper oxidases, and Copper sensing (Cus) system that includes the periplasmic efflux complexes, and copper chaperone (). Several bacteria that can survive in copper-rich environments may harbor additional one or many copies of plasmid-encoded genes that confer copper resistance, such as Pco and Cop systems, which include genes encoding copper oxidases and membrane proteins (; ).

We found a statistically significant difference in copper-related gene abundance between biofilm and water samples (Wilcoxon Rank test, p < 0.05). The statistically significant difference in individual copper gene homologs between biofilm and water samples was confirmed using a binomial test (Figure 2B). Copper-related genes, including copA, cusB, and cusA exhibited statistically significant enrichment in the biofilm metagenome (Binomial test, adjusted p-value < 0.05) (Figures 2A, B).

Figure 2

Among all copper-related genes, the P-type ATPase copA, which exports Cu(I) from the cytoplasm to the periplasm, was the most abundant and exhibited the most significant difference between the two samples. Additionally, the biofilm showed significant enrichment of the periplasmic copper efflux system cusA and cusB genes, enabling proton-driven periplasmic copper ion translocation across the outer membrane (). This observed enrichment of the copper protein from biofilm metagenome highlights the critical role of bacterial copper efflux mechanisms within this environment.

Taxonomic assignment to copper-associated genes

The reads annotated as copper-associated genes from the SEED annotations of the metagenome were taxonomically classified using the Diamond and Megan pipeline to get insights into the bacteria responsible for copper genes in the biofilm environment. This was done by using Diamond and Megan taxonomic analysis on copper-associated reads.

Taxonomic analysis of the biofilm metagenome revealed key contributors to copper surface colonization. We found that copper-related genes found in the biofilm were predominantly linked to the bacterial genera Allomuricauda and Ruegeria (Figure 3). These taxa accounted for over 45% of the biofilm copper-related genes identified through metagenomic analysis. These two taxa also emerged as the most dominant taxa from the overall taxonomic classification. This highlights that the taxa contributing copper genes dominate the early biofilm environment on copper surfaces, suggesting that this allows these organisms to survive on copper surfaces despite copper toxicity and thereby supports the attachment of other fouling organisms.

Figure 3

Overall, our observation suggests the enrichment of copper proteins involved in efflux function in the biofilm (Figure 2). The copper toxicity mitigated by efflux via actively pumping the copper ions out of the cell may lead to an increased localized copper concentration near the biofilm and expose the nearby taxa to higher copper stress. Previous studies have shown spatial gradients of copper within the biofilm when biofilm is grown in the presence of copper. These differences in copper concentration across biofilm have been mainly attributed to the presence of biofilm matrix, which could affect the diffusion of copper within biofilm (; ). Considering the observed enrichment of efflux proteins in biofilm metagenomes, we hypothesize that copper efflux proteins may additionally lead to localized copper accumulation. Studies on bacterial multidrug efflux pumps have shown evidence of increased local concentration of antibiotics, thereby affecting growth of bacterial population in the surrounding microenvironment (). The biofilm environment thus may also experience higher localized copper concentrations because of active efflux mechanisms, thereby explaining the dominance of taxa contributing copper genes.

Conclusion

The genetic pathways encoded within microbes of the marine biofilm environment, and the biological processes actively expressed in marine biofilm habitat, referred to as the functional potential and the functional activity respectively, remain poorly understood. To address this gap, we explored the functional potential by characterizing the microbial diversity and copper specific gene profile of early-stage marine copper-associated biofilms. By comparing the functional profiles associated with copper genes of copper-substrate grown biofilms with those of ambient seawater, functional signatures unique to early colonizers were identified.

The biofilm community was enriched with genes associated with copper resistance, particularly those encoding efflux complexes. These included two efflux proteins: the ATPase-driven P-type copA and the RND transporter proteins of the cusABC system (; ). The CopA P-type ATPase exports Cu(I) ions from the cytoplasm to the periplasm (), while the RND efflux proteins mediate the proton-driven translocation of copper ions across the outer membrane via a periplasmic efflux system ().

The genera Allomuricauda and Ruegeria are dominant bacterial genera that contributed significant copper resistance genes in the biofilm metagenome, even though their presence was less abundant in the seawater metagenome. These two genera accounted for more than 45% of the copper-related genes in the early colonization profile, highlighting that the taxa contributing copper traits dominate early copper surface microbial communities.

While previous studies examined taxonomic characterization of the biofilm communities, we performed the first metagenomic functional analysis of early biofilm colonizers on copper surfaces in marine coastal environments, providing insights into molecular mechanisms that support biofilm formation on copper surfaces. Results from this study creates a foundation for subsequent studies to analyze functional signatures beyond copper systems associated with biofilm formation, that can support the development of new genetic strategies against marine biofouling. This study is limited to profiling the taxonomic composition of microbes and their functional potential of copper-associated biofilm and its surrounding seawater. In this approach, gene abundance has been used as a proxy for microbial activity within the community. It does not, however, determine which genes are actively expressed in the biofilm and seawater communities. As a result, this approach does not offer a clear insight into gene function. To fill this gap, our research must be broadened to incorporate RNA-level analyses that can reflect the active biological processes within the biofilm. Future studies should integrate transcriptomic and proteomic approaches, which can provide direct evidence of actively transcribed genes and quantify proteins. These insights can provide a comprehensive understanding of microbial interactions within the biofilm and can expand on microbial function, metabolic, and regulatory interactions. Additionally, these findings of functional enrichment of copper resistance genes in the copper-associated biofilm could be integrated along with studies characterizing genes related to biofilm formation and surface adhesion, like EPS production, to identify conserved functional activity that can aid in designing effective antifouling strategies against these pathways.

Statements

Data availability statement

The sequence data generated from this study have been deposited in the NCBI Sequence Read Archive under BioProject accession PRJNA1263017.

Author contributions

PG: Investigation, Formal Analysis, Methodology, Validation, Visualization, Data curation, Software, Writing – review & editing, Writing – original draft. VM: Formal Analysis, Methodology, Writing – review & editing. MF: Resources, Project administration, Writing – review & editing, Methodology, Funding acquisition. IE: Writing – review & editing, Software, Resources, Validation, Supervision, Methodology. KC: Project administration, Supervision, Methodology, Conceptualization, Validation, Funding acquisition, Writing – review & editing, Resources.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by UMBC internal funding. Portions of this work were funded through the U.S. Naval Research Laboratories under its base funds to the Nanoscience Institute, but the paper does not represent the U.S. Navy, the Office of Naval Research, or the U.S. Government.

Acknowledgments

We thank Olga Shevchenko and Bruce Brewster Kingham at the University of Delaware DNA Sequencing & Genotyping Center for their assistance in sample quality control validation and Pacbio HiFi sequencing. The hardware used in the computational studies is part of the UMBC High Performance Computing Facility (HPCF). The facility is supported by the U.S. National Science Foundation through the MRI program (grant nos. CNS-0821258, CNS-1228778, OAC-1726023, and CNS-1920079) and the SCREMS program (grant no. DMS-0821311), with additional substantial support from the University of Maryland, Baltimore County (UMBC).

Conflict of interest

Author VM was employed by Precise Systems, Inc.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2025.1661535/full#supplementary-material

References

Summary

Keywords

copper, surface colonization, taxonomic profiling, functional profiling, metagenomics, biofilm, seawater, long-read sequencing

Citation

Gautam P, Molina V, First M, Erill I and Cusick KD (2025) Determining the taxonomic and functional profile of marine bacterial copper systems involved in marine early copper surface colonization. Front. Mar. Sci. 12:1661535. doi: 10.3389/fmars.2025.1661535

Received

07 July 2025

Accepted

08 August 2025

Published

29 August 2025

Volume

12 - 2025

Edited by

Yu Zhang, Shanghai Jiao Tong University, China

Reviewed by

Khaled Mohammed Geba, Menoufia University, Egypt

Hafiz Zeshan Wadood, Lahore Garrison University, Pakistan

Updates

Copyright

*Correspondence: Kathleen D. Cusick,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics