Eukaryotic Biology Research Articles

Gene acquisition by giant transposons primes eukaryotes for rapid evolution via horizontal gene transfer.

Horizontal gene transfer (HGT) disseminates genetic information between species and is a powerful mechanism of adaptation. Yet, we know little about its underlying drivers in eukaryotes. Giant Starship transposons have been implicated as agents of fungal HGT, providing an unprecedented opportunity to reveal the evolutionary parameters behind this process. Here, we characterize the ssf gene cluster, which contributes to formaldehyde resistance, and use it to demonstrate how mobile element evolution shapes fungal adaptation. We found that ssf clusters have been acquired by various distantly related Starships, which each exhibit multiple instances of horizontal transfer across fungal species (at least nine events, including between different taxonomic orders). Many ssf clusters have subsequently integrated into their host's genome, illustrating how Starships shape the evolutionary trajectory of fungal hosts beyond any single transfer. Our results demonstrate the key role Starships play in mediating rapid and repeated adaptation via HGT, elevating the importance of mobile element evolution in eukaryotic biology.

Genetic Code Expansion Approaches to Decipher the Ubiquitin Code.

The covalent attachment of Ub (ubiquitin) to target proteins (ubiquitylation) represents one of the most versatile PTMs (post-translational modifications) in eukaryotic cells. Substrate modifications range from a single Ub moiety being attached to a target protein to complex Ub chains that can also contain Ubls (Ub-like proteins). Ubiquitylation plays pivotal roles in most aspects of eukaryotic biology, and cells dedicate an orchestrated arsenal of enzymes to install, translate, and reverse these modifications. The entirety of this complex system is coined the Ub code. Deciphering the Ub code is challenging due to the difficulty in reconstituting enzymatic machineries and generating defined Ub/Ubl-protein conjugates. This Review provides a comprehensive overview of recent advances in using GCE (genetic code expansion) techniques to study the Ub code. We highlight strategies to site-specifically ubiquitylate target proteins and discuss their advantages and disadvantages, as well as their various applications. Additionally, we review the potential of small chemical PTMs targeting Ub/Ubls and present GCE-based approaches to study this additional layer of complexity. Furthermore, we explore methods that rely on GCE to develop tools to probe interactors of the Ub system and offer insights into how future GCE-based tools could help unravel the complexity of the Ub code.

A functional screen for ubiquitin regulation identifies an E3 ligase secreted by Pseudomonas aeruginosa.

Ubiquitin signaling controls many aspects of eukaryotic biology, including targeted protein degradation and immune defense. Remarkably, invading bacterial pathogens have adapted secreted effector proteins that hijack host ubiquitination to gain control over host responses. These ubiquitin-targeted effectors can exhibit, for example, E3 ligase or deubiquitinase activities, often without any sequence or structural homology to eukaryotic ubiquitin regulators. Such convergence in function poses a challenge to the discovery of additional bacterial virulence factors that target ubiquitin. To overcome this, we have developed a workflow to harvest natively secreted bacterial effectors and functionally screen them for ubiquitin regulatory activities. After benchmarking this approach on diverse ligase and deubiquitinase activities from Salmonella Typhimurium, Enteropathogenic Escherichia coli, and Shigella flexneri, we applied it to the identification of a cryptic E3 ligase activity secreted by Pseudomonas aeruginosa. We identified an unreported P. aeruginosa E3 ligase, which we have termed Pseudomonas Ub ligase 1 (PUL-1), that resembles none of the other E3 ligases previously established in or outside of the eukaryotic system. Importantly, in an animal model of P. aeruginosa infection, PUL-1 ligase activity plays an important role in regulating virulence. Thus, our workflow for the functional identification of ubiquitin-targeted effector proteins carries promise for expanding our appreciation of how host ubiquitin regulation contributes to bacterial pathogenesis.

Measuring G protein activation by spectrally resolved imaging fluorescence fluctuation spectroscopy

The activation of heterotrimeric G proteins through G-protein-coupled receptors (GPCRs) is a ubiquitous signaling mechanism in eukaryotic biology. The three principal molecular components of this cascade are the GPCR, Gα subunit, and Gβγ subunit. Measurement of interactions between these components and their downstream effectors in live cells is paramount to understanding how cells fine-tune their physiology in response to many external stimuli. Multicolor fluorescence fluctuation spectroscopy (FFS) approaches allow the sensitive detection of heteromeric interactions by using spectrally distinct fluorophores to label biomolecules of interest. We considered three imaging FFS approaches to measuring molecular interactions from the signals produced by a spectrally resolved confocal microscopy: raster spectral image correlation spectroscopy (RSICS), spectral spatial cumulant analysis, and native resolution spatial cumulant analysis. We characterized these approaches using simulation and experiments on heteromers with known stoichiometries. We found that RSICS had the best sensitivity for measuring heteromeric interactions and employed it to measure G protein complexes. As measured by RSICS, interactions between the G protein subunits Gαi1 and Gβ1γ2 were sensitive to the stimulation of two GPCRs, the D2 dopamine receptor and the α-2A adrenergic receptor. Interactions between GPCRs and G proteins were not detectable above background, supporting a collisional model of GPCR/G protein interactions in contrast to a preassembly model where strong interactions would be present. These data are uniquely available by this FFS framework, which is appropriate for not only multiplexed measurements of G protein biology but any dynamic protein complexes in the cell.

Atomic Tailoring of Ubiquitin Side Chains Influences E2‐Mediated Ubiquitin Chain Formation

AbstractUbiquitin (Ub) is a small, highly conserved protein essential for eukaryotic biology, and is unique in its formation of polyubiquitin chains by conjugation to one of its seven lysine side chains. Here we report that atomic tailoring of Ub side chains – i. e. the insertion, deletion, or replacement of specific atoms – has significant and unexpected consequences on the enzymatic conjugation of Ub oligomers by isopeptide bond formation mediated by E2 conjugating enzymes. These studies employed chemical synthesis and ligation methods to prepare numerous specifically tailored Ub monomers on multi‐milligram scales. While some modifications including N‐terminal acylation and methionine replacement did not affect protein folding or Ub chain formation with Ube2 K, other modifications had a pronounced effect of oligomerization with Ubc13/Mms2. We observed that Ala46Hse mutation obliterates the ability of this Ub monomer to accept another Ub at Lys63 in Ubc13‐mediated conjugations. Exhaustive replacement of all seven lysines with shorter surrogates Orn, Dab, or Dap essentially blocks Ub chain formation, and in the case of Dap, precludes proper folding of the Ub protein.

The translocation activity of Rad54 reduces crossover outcomes during homologous recombination.

Homologous recombination (HR) is a template-based DNA double-strand break repair pathway that requires the selection of an appropriate DNA sequence to facilitate repair. Selection occurs during a homology search that must be executed rapidly and with high fidelity. Failure to efficiently perform the homology search can result in complex intermediates that generate genomic rearrangements, a hallmark of human cancers. Rad54 is an ATP dependent DNA motor protein that functions during the homology search by regulating the recombinase Rad51. How this regulation reduces genomic exchanges is currently unknown. To better understand how Rad54 can reduce these outcomes, we evaluated several amino acid mutations in Rad54 that were identified in the COSMIC database. COSMIC is a collection of amino acid mutations identified in human cancers. These substitutions led to reduced Rad54 function and the discovery of a conserved motif in Rad54. Through genetic, biochemicaland single-molecule approaches, we show that disruption of this motif leads to failure in stabilizing early strand invasion intermediates, causing increased crossovers between homologous chromosomes. Our study also suggests that the translocation rate of Rad54 is a determinant in balancing genetic exchange. The latch domain's conservation implies an interaction likely fundamental to eukaryotic biology.

Evaluating protein prenylation of human and viral CaaX sequences using a humanized yeast system.

Prenylated proteins are prevalent in eukaryotic biology (∼1-2% of proteins) and are associated with human disease, including cancer, premature aging and infections. Prenylated proteins with a C-terminal CaaX sequence are targeted by CaaX-type prenyltransferases and proteases. To aid investigations of these enzymes and their targets, we developed Saccharomyces cerevisiae strains that express these human enzymes instead of their yeast counterparts. These strains were developed in part to explore human prenyltransferase specificity because of findings that yeast FTase has expanded specificity for sequences deviating from the CaaX consensus (i.e. atypical sequence and length). The humanized yeast strains displayed robust prenyltransferase activity against CaaX sequences derived from human and pathogen proteins containing typical and atypical CaaX sequences. The system also recapitulated prenylation of heterologously expressed human proteins (i.e. HRas and DNAJA2). These results reveal that substrate specificity is conserved for yeast and human farnesyltransferases but is less conserved for type I geranylgeranyltransferases. These yeast systems can be easily adapted for investigating the prenylomes of other organisms and are valuable new tools for helping define the human prenylome, which includes physiologically important proteins for which the CaaX modification status is unknown.

Thermotolerance in S. cerevisiae as a model to study extracellular vesicle biology.

The budding yeast Saccharomyces cerevisiae is a proven model organism for elucidating conserved eukaryotic biology, but to date its extracellular vesicle (EV) biology is understudied. Here, we show yeast transmit information through the extracellular medium that increases survival when confronted with heat stress and demonstrate the EV-enriched samples mediate this thermotolerance transfer. These samples contain vesicle-like particles that are exosome-sized and disrupting exosome biogenesis by targeting endosomal sorting complexes required for transport (ESCRT) machinery inhibits thermotolerance transfer. We find that Bro1, the yeast ortholog of the human exosome biomarker ALIX, is present in EV samples, and use Bro1 tagged with green fluorescent protein (GFP) to track EV release and uptake by endocytosis. Proteomics analysis reveals that heat shock protein 70 (HSP70) family proteins are enriched in EV samples that provide thermotolerance. We confirm the presence of the HSP70 ortholog stress-seventy subunit A2 (Ssa2) in EV samples and find that mutant yeast cells lacking SSA2 produce EVs but they fail to transfer thermotolerance. We conclude that Ssa2 within exosomes shared between yeast cells contributes to thermotolerance. Through this work, we advance Saccharomyces cerevisiae as an emerging model organism for elucidating molecular details of eukaryotic EV biology and establish a role for exosomes in heat stress and proteostasis that seems to be evolutionarily conserved.

Expanding TheCellVision.org: a central repository for visualizing and mining high-content cell imaging projects.

We previously constructed TheCellVision.org, a central repository for visualizing and mining data from yeast high-content imaging projects. At its inception, TheCellVision.org housed two high-content screening (HCS) projects providing genome-scale protein abundance and localization information for the budding yeast Saccharomyces cerevisiae, as well as a comprehensive analysis of the morphology of its endocytic compartments upon systematic genetic perturbation of each yeast gene. Here, we report on the expansion of TheCellVision.org by the addition of two new HCS projects and the incorporation of new global functionalities. Specifically, TheCellVision.org now hosts images from the Cell Cycle Omics project, which describes genome-scale cell cycle-resolved dynamics in protein localization, protein concentration, gene expression, and translational efficiency in budding yeast. Moreover, it hosts PIFiA, a computational tool for image-based predictions of protein functional annotations. Across all its projects, TheCellVision.org now houses >800,000 microscopy images along with computational tools for exploring both the images and their associated datasets. Together with the newly added global functionalities, which include the ability to query genes in any of the hosted projects using either yeast or human gene names, TheCellVision.org provides an expanding resource for single-cell eukaryotic biology.

Mechanism of millisecond Lys48-linked poly-ubiquitin chain formation by cullin-RING ligases

E3 ubiquitin ligases, in collaboration with E2 ubiquitin-conjugating enzymes, modify proteins with poly-ubiquitin chains. Cullin-RING ligase (CRL) E3s use Cdc34/UBE2R-family E2s to build Lys48-linked poly-ubiquitin chains to control an enormous swath of eukaryotic biology. Yet the molecular mechanisms underlying this exceptional linkage specificity and millisecond kinetics of poly-ubiquitylation remain unclear. Here we obtain cryogenic-electron microscopy (cryo-EM) structures that provide pertinent insight into how such poly-ubiquitin chains are forged. The CRL RING domain not only activates the E2-bound ubiquitin but also shapes the conformation of a distinctive UBE2R2 loop, positioning both the ubiquitin to be transferred and the substrate-linked acceptor ubiquitin within the active site. The structures also reveal how the ubiquitin-like protein NEDD8 uniquely activates CRLs during chain formation. NEDD8 releases the RING domain from the CRL, but unlike previous CRL–E2 structures, does not contact UBE2R2. These findings suggest how poly-ubiquitylation may be accomplished by many E2s and E3s.

Evaluating the Activity of Ras/Rap GTPases in Dictyostelium.

Ras and Rap small GTPases of the Ras superfamily act as molecular switches to control diverse cellular processes as part of different signaling pathways. Dictyostelium expresses several Ras and Rap proteins, and their study has and continues to greatly contribute to our understanding of their role in eukaryote biology. To study the activity of Ras and Rap proteins in Dictyostelium, several assays based on their interaction with the Ras binding domain of known eukaryotic Ras/Rap effectors have been developed and proved extremely useful to study their regulation and cellular roles. Here, we describe methods to assess Ras/Rap activity biochemically using a pull-down assay and through live-cell imaging using fluorescent reporters.

ATG9B is a tissue-specific homotrimeric lipid scramblase that can compensate for ATG9A

ABSTRACT Macroautophagy/autophagy is a fundamental aspect of eukaryotic biology, and the autophagy-related protein ATG9A is part of the core machinery facilitating this process. In addition to ATG9A vertebrates encode ATG9B, a poorly characterized paralog expressed in a subset of tissues. Herein, we characterize the structure of human ATG9B revealing the conserved homotrimeric quaternary structure and explore the conformational dynamics of the protein. Consistent with the experimental structure and computational chemistry, we establish that ATG9B is a functional lipid scramblase. We show that ATG9B can compensate for the absence of ATG9A in starvation-induced autophagy displaying similar subcellular trafficking and steady-state localization. Finally, we demonstrate that ATG9B can form a heteromeric complex with ATG2A. By establishing the molecular structure and function of ATG9B, our results inform the exploration of niche roles for autophagy machinery in more complex eukaryotes and reveal insights relevant across species. Abbreviation: ATG: autophagy related; CHS: cholesteryl hemisuccinate; cryo-EM: single-particle cryogenic electron microscopy; CTF: contrast transfer function: CTH: C- terminal α helix; FSC: fourier shell correlation; HDIR: HORMA domain interacting region; LMNG: lauryl maltose neopentyl glycol; MD: molecular dynamics simulations; MSA: multiple sequence alignment; NBD-PE: 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(7-nitro-2–1,3-benzoxadiazol-4-yl ammonium salt); POPC: palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine; RBG: repeating beta groove domain; RMSD: root mean square deviation; SEC: size-exclusion chromatography; TMH: transmembrane helix

An orthogonalized PYR1-based CID module with reprogrammable ligand-binding specificity

Plants sense abscisic acid (ABA) using chemical-induced dimerization (CID) modules, including the receptor PYR1 and HAB1, a phosphatase inhibited by ligand-activated PYR1. This system is unique because of the relative ease with which ligand recognition can be reprogrammed. To expand the PYR1 system, we designed an orthogonal ‘*’ module, which harbors a dimer interface salt bridge; X-ray crystallographic, biochemical and in vivo analyses confirm its orthogonality. We used this module to create PYR1*MANDI/HAB1* and PYR1*AZIN/HAB1*, which possess nanomolar sensitivities to their activating ligands mandipropamid and azinphos-ethyl. Experiments in Arabidopsis thaliana and Saccharomyces cerevisiae demonstrate the sensitive detection of banned organophosphate contaminants using living biosensors and the construction of multi-input/output genetic circuits. Our new modules enable ligand-programmable multi-channel CID systems for plant and eukaryotic synthetic biology that can empower new plant-based and microbe-based sensing modalities.

Reversible crowdedness of pH-responsive and host–guest active polymersomes: Mimicking µm-sized cell structures

The structure–function characteristics of isolated artificial organelles (AOs) in protocells are mainly known, but there are few reports on clustered or aggregated AOs. To imitate µm-sized complex and heterogeneous cell structures, approaches are needed that enable reversible changes in the aggregation state of colloidal structures in response to chemical, biological, and external stimuli. To construct adaptive organelle-like or cell-like reorganization characteristics, we present an advanced crosslinking strategy to fabricate clustered polymersomes as a platform based on host–guest interactions between azobenzene-containing polymersomes (Azo-Psomes) and a β-cyclodextrin-modified polymer (β-CD polymer) as a crosslinker. First, the reversible (dis)assembly of clustered Azo-Psomes is carried out by the alternating input of crosslinker and adamantane-PEG3000 as a decrosslinker. Moreover, cluster size dependence is demonstrated by environmental pH. These offer the controlled fabrication of various homogeneous and heterogeneous Azo-Psomes structures, including the size regulation and visualization of clustered AOs through a fluorescent enzymatic cascade reaction. Finally, a temperature-sensitive crosslinking agent with β-CD units can promote the coaggregation of Azo-Psomes mediated by temperature changes. Overall, these (co–)clustered Azo-Psomes and their successful transformation in AOs may provide new features for modelling biological systems for eukaryotic cells and systems biology.

The RNA-binding protein RBP42 regulates cellular energy metabolism in mammalian-infective Trypanosoma brucei.

RNA-binding proteins (RBPs) are key players in coordinated post-transcriptional regulation of functionally related genes, defined as RNA regulons. RNA regulons play particularly critical roles in parasitic trypanosomes, which exhibit unregulated co-transcription of long unrelated gene arrays. In this report, we present a systematic analysis of an essential RBP, RBP42, in the mammalian-infective bloodstream form of African trypanosome and show that RBP42 is a key regulator of parasite's central carbon and energy metabolism. Using individual-nucleotide resolution UV cross-linking and immunoprecipitation to identify genome-wide RBP42-RNA interactions, we show that RBP42 preferentially binds within the coding region of mRNAs encoding core metabolic enzymes. Global quantitative transcriptomic and proteomic analyses reveal that loss of RBP42 reduces the abundance of target mRNA-encoded proteins, but not target mRNA, suggesting a positive translational regulatory role of RBP42. Significant changes in central carbon metabolic intermediates, following loss of RBP42, further support its critical role in cellular energy metabolism. Trypanosoma brucei infection, transmitted through the bite of blood-feeding tsetse flies, causes deadly diseases in humans and livestock. This disease, if left untreated, is almost always fatal. Existing therapies are toxic and difficult to administer. During T. brucei's lifecycle in two different host environments, the parasite progresses through distinctive life stages with major morphological and metabolic changes, requiring precise alteration of parasite gene expression program. In the absence of regulated transcription, post-transcriptional processes mediated by RNA-binding proteins play critical roles in T. brucei gene regulation. In this study, we show that the RNA-binding protein RBP42 plays crucial roles in cellular energy metabolic regulation of this important human pathogen. Metabolic dysregulation observed in RBP42 knockdown cells offers a breadth of potential interest to researchers studying parasite biology and can also impact research in general eukaryotic biology.

Classification of Features across Five CURE Networks Reveals Opportunities to Improve Course Design, Instruction, and Equity.

Course-based undergraduate research experiences (CUREs) are tools used to introduce students to authentic participation in science. Several specific CUREs have been shown to benefit students' interest and retention in the biological sciences. Nevertheless, CUREs vary greatly in terms of their context, methodology, and degree of research authenticity, so different types of CUREs may differently influence student outcomes. This programmatic diversity poses a challenge to educators who want to better understand which course components and features are reliably present in a CURE curriculum. To address these issues, we identified, catalogued, and classified 112 potential features of CUREs across the biosciences. To develop the list, we interviewed instructors experienced with teaching individual and large networked CUREs across a diversity of the biological disciplines, including: Squirrel-Net (field-based animal behavior), SEA-PHAGES (wet lab microbiology and computational microbiology), Tiny Earth (environmental and wet lab microbiology), PARE (environmental microbiology), and the Genomics Education Partnership (eukaryotic computational biology). Twenty-five interviewees contributed expert content in terms of CURE features and classification of those items into an organized list. The resulting list's categories encompasses student experiences with the following: (i) the scientific process; (ii) technical aspects of science; (iii) the professional development associated with research; and (iv) building scientific identity. The most striking insight was that CUREs vary widely in terms of which features they contain, since different CUREs will by necessity have different approaches to science and student involvement. We also identified several features commonly thought to be crucial to CUREs yet have ambiguous definitions. This ambiguity can potentially confound efforts to make CUREs research-authentic and aligned with the central goals of science. We disambiguate these terms and represent their varied meanings throughout the classification. We also provide instructor-friendly supplementary worksheets along with considerations for instructors interested in expanding their CURE course design, instruction, and equity.

Evolution and phylogenetic distribution of endo-α-mannosidase.

While glycans underlie many biological processes, such as protein folding, cell adhesion, and cell-cell recognition, deep evolution of glycosylation machinery remains an understudied topic. N-linked glycosylation is a conserved process in which mannosidases are key trimming enzymes. One of them is the glycoprotein endo-α-1,2-mannosidase which participates in the initial trimming of mannose moieties from an N-linked glycan inside the cis-Golgi. It is unique as the only endo-acting mannosidase found in this organelle. Relatively little is known about its origins and evolutionary history; so far it was reported to occur only in vertebrates. In this work, a taxon-rich bioinformatic survey to unravel the evolutionary history of this enzyme, including all major eukaryotic clades and a wide representation of animals, is presented. The endomannosidase was found to be more widely distributed in animals and other eukaryotes. The protein motif changes in context of the canonical animal enzyme were tracked. Additionally, the data show the two canonical vertebrate endomannosidase genes, MANEA and MANEAL, arose at the second round of the two vertebrate genome duplications and one more vertebrate paralog, CMANEAL, is uncovered. Finally, a framework where N-glycosylation co-evolved with complex multicellularity is described. A better understanding of the evolution of core glycosylation pathways is pivotal to understanding biology of eukaryotes in general, and the Golgi apparatus in particular. This systematic analysis of the endomannosidase evolution is one step toward this goal.

A high-throughput sensory assay for parasitic and free-living nematodes.

Sensory pathways first elucidated in Caenorhabditis elegans are conserved across free-living and parasitic nematodes, even though each species responds to a diverse array of compounds. Most nematode sensory assays are performed by tallying observations of worm behavior on two-dimensional planes using agarose plates. These assays have been successful in the study of volatile sensation but are poorly suited for investigation of water-soluble gustation or parasitic nematodes without a free-living stage. In contrast, gustatory assays tend to be tedious, often limited to the manipulation of a single individual at a time. We have designed a nematode sensory assay using a microfluidics device that allows for the study of gustation in a 96-well, three-dimensional environment. This device is suited for free-living worms and parasitic worms that spend their lives in an aqueous environment, and we have used it to show that ivermectin inhibits the gustatory ability of vector-borne parasitic nematodes. Insight box Nematodes are powerful model organisms for understanding the sensory biology of multicellular eukaryotes, and many parasitic species cause disease in humans. Simple sensory assays performed on agarose plates have been the bedrock for establishing the neuronal, genetic, and developmental foundations for many sensory modalities in nematodes. However, these classical assays are poorly suited for translational movement of many parasitic nematodes and the sensation of water-soluble molecules (gustation). We have designed a device for high-throughput nematode sensory assays in a gel matrix. This 'gustatory microplate' is amenable to several species and reveals novel responses by free-living and parasitic nematodes to cues and drugs.

Characterisation of the OTU domain deubiquitinase complement of Toxoplasma gondii.

The phylum Apicomplexa contains several parasitic species of medical and agricultural importance. The ubiquitination machinery remains, for the most part, uncharacterised in apicomplexan parasites, despite the important roles that it plays in eukaryotic biology. Bioinformatic analysis of the ubiquitination machinery in apicomplexan parasites revealed an expanded ovarian tumour domain-containing (OTU) deubiquitinase (DUB) family in Toxoplasma, potentially reflecting functional importance in apicomplexan parasites. This study presents comprehensive characterisation of Toxoplasma OTU DUBs. AlphaFold-guided structural analysis not only confirmed functional orthologues found across eukaryotes, but also identified apicomplexan-specific enzymes, subsequently enabling discovery of a cryptic OTU DUB in Plasmodium species. Comprehensive biochemical characterisation of 11 Toxoplasma OTU DUBs revealed activity against ubiquitin- and NEDD8-based substrates and revealed ubiquitin linkage preferences for Lys6-, Lys11-, Lys48-, and Lys63-linked chain types. We show that accessory domains in Toxoplasma OTU DUBs impose linkage preferences, and in case of apicomplexan-specific TgOTU9, we discover a cryptic ubiquitin-binding domain that is essential for TgOTU9 activity. Using the auxin-inducible degron (AID) to generate knockdown parasite lines, TgOTUD6B was found to be important for Toxoplasma growth.

Giant Virus Infection Signatures Are Modulated by Euphotic Zone Depth Strata and Iron Regimes of the Subantarctic Southern Ocean.

Viruses can alter the abundance, evolution, and metabolism of microorganisms in the ocean, playing a key role in water column biogeochemistry and global carbon cycles. Large efforts to measure the contribution of eukaryotic microorganisms (e.g., protists) to the marine food web have been made, yet the in situ activities of the ecologically relevant viruses that infect these organisms are not well characterized. Viruses within the phylum Nucleocytoviricota ("giant viruses") are known to infect a diverse range of ecologically relevant marine protists, yet how these viruses are influenced by environmental conditions remains under-characterized. By employing metatranscriptomic analyses of in situ microbial communities along a temporal and depth-resolved gradient, we describe the diversity of giant viruses at the Southern Ocean Time Series (SOTS), a site within the subpolar Southern Ocean. Using a phylogeny-guided taxonomic assessment of detected giant virus genomes and metagenome-assembled genomes, we observed depth-dependent structuring of divergent giant virus families mirroring dynamic physicochemical gradients in the stratified euphotic zone. Analyses of transcribed metabolic genes from giant viruses suggest viral metabolic reprogramming of hosts from the surface to a 200-m depth. Lastly, using on-deck incubations reflecting a gradient of iron availability, we show that modulating iron regimes influences the activity of giant viruses in the field. Specifically, we show enhanced infection signatures of giant viruses under both iron-replete and iron-limited conditions. Collectively, these results expand our understanding of how the water column's vertical biogeography and chemical surroundings affect an important group of viruses within the Southern Ocean. IMPORTANCE The biology and ecology of marine microbial eukaryotes is known to be constrained by oceanic conditions. In contrast, how viruses that infect this important group of organisms respond to environmental change is less well known, despite viruses being recognized as key microbial community members. Here, we address this gap in our understanding by characterizing the diversity and activity of "giant" viruses within an important region in the sub-Antarctic Southern Ocean. Giant viruses are double-stranded DNA (dsDNA) viruses of the phylum Nucleocytoviricota and are known to infect a wide range of eukaryotic hosts. By employing a metatranscriptomics approach using both in situ samples and microcosm manipulations, we illuminated both the vertical biogeography and how changing iron availability affects this primarily uncultivated group of protist-infecting viruses. These results serve as a foundation for our understanding of how the open ocean water column structures the viral community, which can be used to guide models of the viral impact on marine and global biogeochemical cycling.