Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Horizontal gene transfer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stable hypermutators revealed by the genomic landscape of genes involved in genome stability among yeast species

Mutator phenotypes are short-lived due to the rapid accumulation of deleterious mutations. Yet, recent observations reveal that certain fungi can undergo prolonged accelerated evolution after losing genes involved in DNA repair. Here, we surveyed 1,154 yeast genomes representing nearly all known yeast species of the subphylum Saccharomycotina (phylum Ascomycota) to examine the relationship between reduced gene repertoires broadly associated with genome stability functions (e.g., DNA repair, cell cycle) and elevated evolutionary rates. We identified three distantly related lineages—encompassing 12% of species—that had both the most streamlined sets of genes involved in genome stability (specifically DNA repair) and the highest evolutionary rates in the entire subphylum. Two of these “faster-evolving lineages” (FELs)—a subclade within the order Pichiales and the Wickerhamiella/Starmerella (W/S) clade (order Dipodascales)—are described here for the first time, while the third corresponds to a previously documented Hanseniaspora FEL. Examination of genome stability gene repertoires revealed a set of genes predominantly absent in these three FELs, suggesting a potential role in the observed acceleration of evolutionary rates. In the W/S clade, genomic signatures are consistent with a substantial mutational burden, including pronounced A|T bias and endogenous DNA damage. Interestingly, we found that the W/S clade also contains DNA repair genes possibly acquired through horizontal gene transfer, including a photolyase of bacterial origin. These findings highlight how hypermutators can persist across macroevolutionary timescales, potentially linked to the loss of genes related with genome stability, with horizontal gene transfer as a possible avenue for partial functional compensation.

DNA repair↗
DOE OSTI2021

Engineering organisms resistant to viruses and horizontally transferred genetic elements

Organisms resistant to horizontal gene transfer (HGT), and compositions and methods of use thereof are provided. The organisms are typically genomically recoded organisms (GRO), typically cells, having a genome in which at least one endogenous codon has been eliminated by reassignment of the codon to a synonymous or non-synonymous codon. The GRO typically include a recombinant expression construct for expression of at least one element of a ribosomal rescue pathway, typically lacking the eliminated codon. Typically, the disclosed cells are resistant to completed transfer and/or expression of a horizontally transferred genetic element (HTGE) from another organism compared a corresponding cell having a genome wherein the eliminated codon has not been eliminated. In some embodiments, the organism from which the HTGE is being transferred is a bacterium or a virus.

Isaacs, Farren J.↗
DOE OSTI2022

Phototrophy and carbon fixation in Chlorobi postdate the rise of oxygen

While most productivity on the surface of the Earth today is fueled by oxygenic photosynthesis, for much of Earth history it is thought that anoxygenic photosynthesis—using compounds like ferrous iron or sulfide as electron donors—drove most global carbon fixation. Anoxygenic photosynthesis is still performed by diverse bacteria in niche environments today. Of these, the Chlorobi (formerly green sulfur bacteria) are often interpreted as being particularly ancient and are frequently proposed to have fueled the biosphere during late Archean and early Paleoproterozoic time before the rise of oxygenic photosynthesis. Here, we perform comparative genomic, phylogenetic, and molecular clock analyses to determine the antiquity of the Chlorobi and their characteristic phenotypes. We show that contrary to common assumptions, the Chlorobi clade is relatively young, with anoxygenic phototrophy, carbon fixation via the rTCA pathway, and iron oxidation all significantly postdating the rise of oxygen ~2.3 billion years ago. The Chlorobi therefore could not have fueled the Archean biosphere, but instead represent a relatively young radiation of organisms which likely acquired the capacity for anoxygenic photosynthesis and other traits via horizontal gene transfer sometime after the evolution of oxygenic Cyanobacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Classification of bacterial plasmid and chromosome derived sequences using machine learning

Plasmids are important genetic elements that facilitate horizonal gene transfer between bacteria and contribute to the spread of virulence and antimicrobial resistance. Most bacterial genome sequences in the public archives exist in draft form with many contigs, making it difficult to determine if a contig is of chromosomal or plasmid origin. Using a training set of contigs comprising 10,584 chromosomes and 10,654 plasmids from the PATRIC database, we evaluated several machine learning models including random forest, logistic regression, XGBoost, and a neural network for their ability to classify chromosomal and plasmid sequences using nucleotide k-mers as features. Based on the methods tested, a neural network model that used nucleotide 6-mers as features that was trained on randomly selected chromosomal and plasmid subsequences 5kb in length achieved the best performance, outperforming existing out-of-the-box methods, with an average accuracy of 89.38% ± 2.16% over a 10-fold cross validation. The model accuracy can be improved to 92.08% by using a voting strategy when classifying holdout sequences. In both plasmids and chromosomes, subsequences encoding functions involved in horizontal gene transfer—including hypothetical proteins, transporters, phage, mobile elements, and CRISPR elements—were most likely to be misclassified by the model. This study provides a straightforward approach for identifying plasmid-encoding sequences in short read assemblies without the need for sequence alignment-based tools.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Novel Toxin Biosynthetic Gene Cluster in Harmful Algal Bloom-Causing Heteroscytonema crispum : Insights into the Origins of Paralytic Shellfish Toxins

Caused by both eukaryotic dinoflagellates and prokaryotic cyanobacteria, harmful algal blooms are events of severe ecological, economic, and public health consequence, and their incidence has become more common of late. Despite coordinated research efforts to identify and characterize the genomes of harmful algal bloom-causing organisms, the genomic basis and evolutionary origins of paralytic shellfish toxins produced by harmful algal blooms remain at best incomplete. The paralytic shellfish toxin saxitoxin has an especially complex genomic architecture and enigmatic phylogenetic distribution, spanning dinoflagellates and multiple cyanobacterial genera. Using filtration and extraction techniques to target the desired cyanobacteria from nonaxenic culture, coupled with a combination of short- and long-read sequencing, we generated a reference-quality hybrid genome assembly for Heteroscytonema crispum UTEX LB 1556, a freshwater, paralytic shellfish toxin-producing cyanobacterium thought to have the largest known genome in its phylum. We report a complete, novel biosynthetic gene cluster for the paralytic shellfish toxin saxitoxin. Leveraging this biosynthetic gene cluster, we find support for the hypothesis that paralytic shellfish toxin production has appeared in divergent Cyanobacteria lineages through widespread and repeated horizontal gene transfer. This work demonstrates the utility of long-read sequencing and metagenomic assembly toward advancing our understanding of paralytic shellfish toxin biosynthetic gene cluster diversity and suggests a mechanism for the origin of paralytic shellfish toxin biosynthetic genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Synthetic overlapping genes stabilize genetic systems

Overlapping genes—wherein two different proteins are translated from alternative reading frames of the same DNA sequence—provide a means to stabilize an engineered gene by directly linking its evolutionary fate with that of an overlapping gene. However, creating overlapping gene pairs is challenging, as it requires redesigning both protein products to accommodate overlap constraints. Here, we present a new “overlapping, alternate-frame insertion” (OAFI) method for creating synthetic overlapping genes by inserting an “inner” gene, encoded in an alternate frame, into a flexible region of an “outer” gene. Using OAFI, we create new overlapping gene pairs of genetic reporters and bacterial toxins within an antibiotic resistance gene. We show that both the inner and outer genes retain function despite redesign, with translation of the inner gene influenced by its overlap position in the outer gene. Importantly, we show that, despite these inner gene sequences not contributing to outer gene function, selection for the outer gene alters the permitted inactivating mutations in the inner gene, and that overlapping toxins can restrict horizontal gene transfer of the antibiotic resistance gene. Overall, OAFI offers a versatile tool for synthetic biology, expanding the applications of overlapping genes in gene stabilization and biocontainment.

Biological and medical sciences↗
DOE OSTI2022

Nramp: Deprive and conquer?

Solute carriers 11 (Slc11) evolved from bacterial permease (MntH) to eukaryotic antibacterial defense (Nramp) while continuously mediating proton (H+)-dependent manganese (Mn2+) import. Also, Nramp horizontal gene transfer (HGT) toward bacteria led to mntH polyphyly. Prior demonstration that evolutionary rate-shifts distinguishing Slc11 from outgroup carriers dictate catalytic specificity suggested that resolving Slc11 family tree may provide a function-aware phylogenetic framework. Hence, MntH C (MC) subgroups resulted from HGTs of prototype Nramp (pNs) parologs while archetype Nramp (aNs) correlated with phagocytosis. PHI-Blast based taxonomic profiling confirmed MntH B phylogroup is confined to anaerobic bacteria vs. MntH A (MA)’s broad distribution; suggested niche-related spread of MC subgroups; established that MA-variant MH, which carries ‘eukaryotic signature’ marks, predominates in archaea. Slc11 phylogeny shows MH is sister to Nramp. Site-specific analysis of Slc11 charge network known to interact with the protonmotive force demonstrates sequential rate-shifts that recapitulate Slc11 evolution. 3D mapping of similarly coevolved sites across Slc11 hydrophobic core revealed successive targeting of discrete areas. The data imply that pN HGT could advantage recipient bacteria for H+-dependent Mn2+ acquisition and Alphafold 3D models suggest conformational divergence among MC subgroups. It is proposed that Slc11 originated as a bacterial stress resistance function allowing Mn2+-dependent persistence in conditions adverse for growth, and that archaeal MH could contribute to eukaryogenesis as a Mn2+ sequestering defense perhaps favoring intracellular growth-competent bacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Catabolic pathway acquisition by rhizosphere bacteria readily enables growth with a root exudate component but does not affect root colonization

Horizontal gene transfer (HGT) is a fundamental evolutionary process that plays a key role in bacterial evolution. The likelihood of a successful transfer event is expected to depend on the precise balance of costs and benefits resulting from pathway acquisition. Most experimental analyses of HGT have focused on phenotypes that have large fitness benefits under appropriate selective conditions, such as antibiotic resistance. However, many examples of HGT involve phenotypes that are predicted to provide smaller benefits, such as the ability to catabolize additional carbon sources. We have experimentally simulated the consequences of one such HGT event in the laboratory, studying the effects of transferring a pathway for catabolism of the plant-derived aromatic compound salicyl alcohol between rhizosphere isolates from the Pseudomonas genus. We find that pathway acquisition enables rapid catabolism of salicyl alcohol with only minor disruptions to the existing metabolic and regulatory networks of the new host. However, this new catabolic potential does not confer a measurable fitness advantage during competitive growth in the rhizosphere. We conclude that the phenotype of salicyl alcohol catabolism is readily transferable but is selectively neutral under environmentally relevant conditions. We propose that this condition is common and that HGT of many pathways will be self-limiting because the selective benefits are small.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. In this work, we demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Horizontal transfer of chromosomal DNA mediated by an integrative and conjugative element generates frequent localized recombination in Novosphingobium aromaticivorans

Horizontal gene transfer is an important evolutionary process by which DNA is exchanged between cells that are physically co-located but not direct evolutionary descendants. Horizontal transfer of highly divergent DNA is relatively easy to detect and can produce major phenotypic changes, exemplified by the acquisition of antibiotic resistance determinants. However, transfer of high-identity DNA, for example, between strains of the same species, is likely to be more frequent, harder to detect, and highly impactful in aggregate. In this work, we demonstrate that soil isolates of the alphaproteobacterium Novosphingobium aromaticivorans can exchange chromosomal DNA, leading to multiple unselected recombination events spanning approximately 10% of the chromosome. Chromosomal recombination was directional and more efficient near an integrative and conjugative element (ICE), and required a relaxase found in the ICE. Recombination could not be observed in strains from closely related Novosphingobium species. In combination, these results suggest that ICE-mediated recombination can efficiently recombine DNA within N. aromaticivorans, increasing the adaptive potential of the species while also enforcing species boundaries through preferential intraspecific recombination.

Allemann, Marco [ORNL]↗
DOE OSTI2023

Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea

Sulfate-coupled anaerobic oxidation of methane (AOM) is performed by multicellular consortia of anaerobic methanotrophic archaea (ANME) in obligate syntrophic partnership with sulfate-reducing bacteria (SRB). Diverse ANME and SRB clades co-associate but the physiological basis for their adaptation and diversification is not well understood. In this work, we used comparative metagenomics and phylogenetics to investigate the metabolic adaptation among the 4 main syntrophic SRB clades (HotSeep-1, Seep-SRB2, Seep-SRB1a, and Seep-SRB1g) and identified features associated with their syntrophic lifestyle that distinguish them from their non-syntrophic evolutionary neighbors in the phylum Desulfobacterota. We show that the protein complexes involved in direct interspecies electron transfer (DIET) from ANME to the SRB outer membrane are conserved between the syntrophic lineages. In contrast, the proteins involved in electron transfer within the SRB inner membrane differ between clades, indicative of convergent evolution in the adaptation to a syntrophic lifestyle. Our analysis suggests that in most cases, this adaptation likely occurred after the acquisition of the DIET complexes in an ancestral clade and involve horizontal gene transfers within pathways for electron transfer (CbcBA) and biofilm formation (Pel). We also provide evidence for unique adaptations within syntrophic SRB clades, which vary depending on the archaeal partner. Among the most widespread syntrophic SRB, Seep-SRB1a, subclades that specifically partner ANME-2a are missing the cobalamin synthesis pathway, suggestive of nutritional dependency on its partner, while closely related Seep-SRB1a partners of ANME-2c lack nutritional auxotrophies. Our work provides insight into the features associated with DIET-based syntrophy and the adaptation of SRB towards it.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Pore confinement enhances but surface adhesion reduces bacterial cell-to-cell conjugation

As the habitats of bacteria, soil pore network and surface properties control the distribution, adhesion, and motility of bacteria in soils. These physical processes in turn influence bacterial accesses to nutrients and bacterial interactions. Our understanding on the pore- and surface-mediated bacterial interactions is currently limited. In this research, we evaluated the effects of soil pore confinement and surface adhesion on conjugation-based bacterial interactions. The interaction was measured by plasmid transfer between donor and recipient cells within the population of soil bacterium Pseudomonas putida. We found that the presence of porous sand media led to a net increase in conjugation frequency compared to sand-free liquid control. The increase is attributed to the facilitated effect of pore confinement on the collision of bacteria within pores. In contrast, bacterial adhesion to sand surfaces under elevated ionic strength conditions decreased the conjugation frequency as a result of mobility reduction on the surface. Such collision and adhesion mechanisms jointly drive the conjugation as a function of pore and surface properties of porous media. Finally, these results provide valuable insights into the roles of soil pores and surfaces in regulating horizontal gene transfer, an essential cell-to-cell interaction sustaining key processes of soil ecology and health.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Characterization of the thermophilic xylanase Fsa02490Xyn from the hyperthermophile Fervidibacter sacchari belonging to glycoside hydrolase family 10

Fervidibacter sacchari is an aerobic hyperthermophile belonging to the phylum Armatimonadota that degrades a variety of polysaccharides. Its genome encodes 117 enzymes with one or more annotated glycoside hydrolase (GH) domain, but the roles of these putative GHs in polysaccharide catabolism are poorly defined. Here, we describe one F. sacchari enzyme encoding a GH10 domain, Fsa02490Xyn, that was previously shown to be active on Miscanthus, oat β-glucan, and beech-wood xylan, with optimal activity at 90-100 °C. We show that Fsa02490Xyn is also active on birch-wood xylan and gellan gum. The pH range on beech-wood xylan was 4.5 to 9.5 (pHopt 7.0-8.0). Fsa024940Xyn had a Km of 2.375 mm, Vmax of 1250 μm·min-1, and kcat/Km of 1.259 × 104 s-1·m-1 when using a para-nitrophenyl-?-xylobioside assay. A phylogenetic analysis of GH10 family enzymes revealed a large clade of enzymes from diverse members of the class Fervidibacteria, including Fsa02490Xyn and a second enzyme from F. sacchari, with apparent horizontal gene transfer within Fervidibacteria and between Fervidibacteria and thermophilic Bacillota. This study establishes Fsa02490Xyn as a hyperthermophilic GH10 enzyme with endo-β-1,4-xylanase activity and identifies a large clade of homologous GH10 enzymes within the class Fervidibacteria. Impact statement The depolymerization of xylan at high temperatures is important because this process limits the degradation of polysaccharides in nature and the synthesis of biofuels from plant wastes. Our study is also important because F. sacchari is one of only a few cultivated members of the Armatimonadota, which are polysaccharide-degradation specialists.

Armatimonadota↗
DOE OSTI2025

Comparative genomics of Aspergillus nidulans and section Nidulantes

Aspergillus nidulans is an important model organism for eukaryotic biology and the reference for the section Nidulantes in comparative studies. In this study, we de novo sequenced the genomes of 25 species of this section. Whole-genome phylogeny of 34 Aspergillus species and Penicillium chrysogenum clarifies the position of clades inside section Nidulantes. Comparative genomics reveals a high genetic diversity between species with 684 up to 2433 unique protein families. Furthermore, we categorized 2118 secondary metabolite gene clusters (SMGC) into 603 families across Aspergilli, with at least 40 % of the families shared between Nidulantes species. Genetic dereplication of SMGC and subsequent synteny analysis provides evidence for horizontal gene transfer of a SMGC. Proteins that have been investigated in A. nidulans as well as its SMGC families are generally present in the section Nidulantes, supporting its role as model organism. The set of genes encoding plant biomass-related CAZymes is highly conserved in section Nidulantes, while there is remarkable diversity of organization of MAT-loci both within and between the different clades. This study provides a deeper understanding of the genomic conservation and diversity of this section and supports the position of A. nidulans as a reference species for cell biology.

Theobald, Sebastian [Technical University of Denma↗
DOE OSTI2022

Virioplankton assemblages from challenger deep, the deepest place in the oceans

Hadal ocean biosphere, that is, the deepest part of the world’s oceans, harbors a unique microbial community, suggesting a potential uncovered co-occurring virioplankton assemblage. Herein, we reveal the unique virioplankton assemblages of the Challenger Deep, comprising 95,813 non-redundant viral contigs from the surface to the hadal zone. Almost all of the dominant viral contigs in the hadal zone were unclassified, potentially related to Alteromonadales and Oceanospirillales. 2,586 viral auxiliary metabolic genes from 132 different KEGG orthologous groups were mainly related to the carbon, nitrogen, sulfur, and arsenic metabolism. Lysogenic viral production and integrase genes were augmented in the hadal zone, suggesting the prevalence of viral lysogenic life strategy. Abundant rve genes in the hadal zone, which function as transposase in the caudoviruses, further suggest the prevalence of viral-mediated horizontal gene transfer. This study provides fundamental insights into the virioplankton assemblages of the hadal zone, reinforcing the necessity of incorporating virioplankton into the hadal biogeochemical cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Expanding the genomic encyclopedia of Actinobacteria with 824 isolate reference genomes

The phylum Actinobacteria includes important human pathogens like Mycobacterium tuberculosis and Corynebacterium diphtheriae and renowned producers of secondary metabolites of commercial interest, yet only a small part of its diversity is represented by sequenced genomes. Here, we present 824 actinobacterial isolate genomes in the context of a phylum-wide analysis of 6,700 genomes including public isolates and metagenome-assembled genomes (MAGs). We estimate that only 30%–50% of projected actinobacterial phylogenetic diversity possesses genomic representation via isolates and MAGs. A comparison of gene functions reveals novel determinants of host-microbe interaction as well as environment-specific adaptations such as potential antimicrobial peptides. We identify plasmids and prophages across isolates and uncover extensive prophage diversity structured mainly by host taxonomy. Analysis of >80,000 biosynthetic gene clusters reveals that horizontal gene transfer and gene loss shape secondary metabolite repertoire across taxa. Our observations illustrate the essential role of and need for high-quality isolate genome sequences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Plant microbiomes harbor potential to promote nutrient turnover in impoverished substrates of a Brazilian biodiversity hotspot

The substrates of the Brazilian campos rupestres, a grassland ecosystem, have extremely low concentrations of phosphorus and nitrogen, imposing restrictions to plant growth. Despite that, this ecosystem harbors almost 15% of the Brazilian plant diversity, raising the question of how plants acquire nutrients in such a harsh environment. Here, we set out to uncover the taxonomic profile, the compositional and functional differences and similarities, and the nutrient turnover potential of microbial communities associated with two plant species of the campos rupestres-dominant family Velloziaceae that grow over distinct substrates (soil and rock). Using amplicon sequencing data, we show that, despite the pronounced composition differentiation, the plant-associated soil and rock communities share a core of highly efficient colonizers that tend to be highly abundant and is enriched in 21 bacterial families. Functional investigation of metagenomes and 522 metagenome-assembled genomes revealed that the microorganisms found associated to plant roots are enriched in genes involved in organic compound intake, and phosphorus and nitrogen turnover. We show that potential for phosphorus transport, mineralization, and solubilization are mostly found within bacterial families of the shared microbiome, such as Xanthobacteraceae and Bryobacteraceae. We also detected the full repertoire of nitrogen cycle-related genes and discovered a lineage of Isosphaeraceae that acquired nitrogen-fixing potential via horizontal gene transfer and might be also involved in nitrification via a metabolic handoff association with Binataceae. We highlight that plant-associated microbial populations in the campos rupestres harbor a genetic repertoire with potential to increase nutrient availability and that the microbiomes of biodiversity hotspots can reveal novel mechanisms of nutrient turnover.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Patterns and ecological drivers of viral communities in acid mine drainage sediments across Southern China

Recent advances in environmental genomics have provided unprecedented opportunities for the investigation of viruses in natural settings. Yet, our knowledge of viral biogeographic patterns and the corresponding drivers is still limited. Here, we perform metagenomic deep sequencing on 90 acid mine drainage (AMD) sediments sampled across Southern China and examine the biogeography of viruses in this extreme environment. The results demonstrate that prokaryotic communities dictate viral taxonomic and functional diversity, abundance and structure, whereas other factors especially latitude and mean annual temperature also impact viral populations and functions. In silico predictions highlight lineage-specific virus-host abundance ratios and richness-dependent virus-host interaction structure. Further functional analyses reveal important roles of environmental conditions and horizontal gene transfers in shaping viral auxiliary metabolic genes potentially involved in phosphorus assimilation. Our findings underscore the importance of both abiotic and biotic factors in predicting the taxonomic and functional biogeographic dynamics of viruses in the AMD sediments.

59 BASIC BIOLOGICAL SCIENCES↗