Amoebozoan testate amoebae illuminate the diversity of heterotrophs and the complexity of ecosystems throughout geological time.

Heterotrophic protists are vital in Earth's ecosystems, influencing carbon and nutrient cycles and occupying key positions in food webs as microbial predators. Fossils and molecular data suggest the emergence of predatory microeukaryotes and the transition to a eukaryote-rich marine environment by 800 million years ago (Ma). Neoproterozoic vase-shaped microfossils (VSMs) linked to Arcellinida testate amoebae represent the oldest evidence of heterotrophic microeukaryotes. This study explores the phylogenetic relationship and divergence times of modern Arcellinida and related taxa using a relaxed molecular clock approach. We estimate the origin of nodes leading to extant members of the Arcellinida Order to have happened during the latest Mesoproterozoic and Neoproterozoic (1054 to 661 Ma), while the divergence of extant infraorders postdates the Silurian. Our results demonstrate that at least one major heterotrophic eukaryote lineage originated during the Neoproterozoic. A putative radiation of eukaryotic groups (e.g., Arcellinida) during the early-Neoproterozoic sustained by favorable ecological and environmental conditions may have contributed to eukaryotic life endurance during the Cryogenian severe ice ages. Moreover, we infer that Arcellinida most likely already inhabited terrestrial habitats during the Neoproterozoic, coexisting with terrestrial Fungi and green algae, before land plant radiation. The most recent extant Arcellinida groups diverged during the Silurian Period, alongside other taxa within Fungi and flowering plants. These findings shed light on heterotrophic microeukaryotes' evolutionary history and ecological significance in Earth's ecosystems, using testate amoebae as a proxy.

Mitochondrion-to-nucleus communication mediated by RNA export: a survey of potential mechanisms and players across eukaryotes.

The nucleus interacts with the other organelles to perform essential functions of the eukaryotic cell. Mitochondria have their own genome and communicate back to the nucleus in what is known as mitochondrial retrograde response. Information is transferred to the nucleus in many ways, leading to wide-ranging changes in nuclear gene expression and culminating with changes in metabolic, regulatory or stress-related pathways. RNAs are emerging molecules involved in this signalling. RNAs encode precise information and are involved in highly target-specific signalling, through a wide range of processes known as RNA interference. RNA-mediated mitochondrial retrograde response requires these molecules to exit the mitochondrion, a process that is still mostly unknown. We suggest that the proteins/complexes translocases of the inner membrane, polynucleotide phosphorylase, mitochondrial permeability transition pore, and the subunits of oxidative phosphorylation complexes may be responsible for RNA export.

Impact of steroid biosynthesis on the aerobic adaptation of eukaryotes.

Steroids are indispensable components of the eukaryotic cellular membrane and the acquisition of steroid biosynthesis was a key factor that enabled the evolution of eukaryotes. The polycyclic carbon structures of steroids can be preserved in sedimentary rocks as chemical fossils for billions of years and thus provide invaluable clues to trace eukaryotic evolution from the distant past. Steroid biosynthesis consists of (1) the production of protosteroids and (2) the subsequent modifications toward "modern-type" steroids such as cholesterol and stigmasterol. While protosteroid biosynthesis requires only two genes for the cyclization of squalene, complete modification of protosteroids involves ~10 additional genes. Eukaryotes universally possess at least some of those additional genes and thus produce modern-type steroids as major final products. The geological biomarker records suggest a prolonged period of solely protosteroid production in the mid-Proterozoic before the advent of modern-type steroids in the Neoproterozoic. It has been proposed that mid-Proterozoic protosteroids were produced by hypothetical stem-group eukaryotes that presumably possessed genes only for protosteroid production, even though in modern environments protosteroid production as a final product is found exclusively in bacteria. The host identity of mid-Proterozoic steroid producers is crucial for understanding the early evolution of eukaryotes. In this perspective, we discuss how geological biomarker data and genetic data complement each other and potentially provide a more coherent scenario for the evolution of steroids and associated early eukaryotes. We further discuss the potential impacts that steroids had on the evolution of aerobic metabolism in eukaryotes, which may have been an important factor for the eventual ecological dominance of eukaryotes in many modern environments.

LSTrAP-denovo: Automated Generation of Transcriptome Atlases for Eukaryotic Species Without Genomes.

Despite the abundance of species with transcriptomic data, a significant number of species still lack sequenced genomes, making it difficult to study gene function and expression in these organisms. While de novo transcriptome assembly can be used to assemble protein-coding transcripts from RNA-sequencing (RNA-seq) data, the datasets used often only feature samples of arbitrarily selected or similar experimental conditions, which might fail to capture condition-specific transcripts. We developed the Large-Scale Transcriptome Assembly Pipeline for de novo assembled transcripts (LSTrAP-denovo) to automatically generate transcriptome atlases of eukaryotic species. Specifically, given an NCBI TaxID, LSTrAP-denovo can (1) filter undesirable RNA-seq accessions based on read data, (2) select RNA-seq accessions via unsupervised machine learning to construct a sample-balanced dataset for download, (3) assemble transcripts via over-assembly, (4) functionally annotate coding sequences (CDS) from assembled transcripts and (5) generate transcriptome atlases in the form of expression matrices for downstream transcriptomic analyses. LSTrAP-denovo is easy to implement, written in Python, and is freely available at https://github.com/pengkenlim/LSTrAP-denovo/.

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes.

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

Temporal variations of bacterial and eukaryotic community in coastal waters-implications for aquaculture.

Despite increased attention to the aquaculture environment, there is still a lack of understanding regarding the significance of water quality. To address this knowledge gap, this study utilized high-throughput sequencing of 16S rRNA and 18S rRNA to examine microbial communities (bacteria and eukaryotes) in coastal water over different months through long-term observations. The goal was to explore interaction patterns in the microbial community and identify potential pathogenic bacteria and red tide organisms. The results revealed significant differences in composition, diversity, and richness of bacterial and eukaryotic operational taxonomic units (OTUs) across various months. Principal coordinate analysis (PCoA) demonstrated distinct temporal variations in bacterial and eukaryotic communities, with significant differences (P = 0.001) among four groups: F (January-April), M (May), S (June-September), and T (October-December). Moreover, a strong association was observed between microbial communities and months, with most OTUs showing a distinct temporal preference. The Kruskal-Wallis test (P < 0.05) indicated significant differences in dominant bacterial and eukaryotic taxa among months, with each group exhibiting unique dominant taxa, including potential pathogenic bacteria and red tide organisms. These findings emphasize the importance of monitoring changes in potentially harmful microorganisms in aquaculture. Network analysis highlighted positive correlations between bacteria and eukaryotes, with bacteria playing a key role in network interactions. The key bacterial genera associated with other microorganisms varied significantly (P < 0.05) across different groups. In summary, this study deepens the understanding of aquaculture water quality and offers valuable insights for maintaining healthy aquaculture practices. KEY POINTS: • Bacterial and eukaryotic communities displayed distinct temporal variations. • Different months exhibited unique potential pathogenic bacteria and red tide organisms. • Bacteria are key taxonomic taxa involved in microbial network interactions.

Nitrate reduction to ammonium: a phylogenetic, physiological, and genetic aspects in Prokaryotes and eukaryotes.

The microbe-mediated conversion of nitrate (NO3-) to ammonium (NH4+) in the nitrogen cycle has strong implications for soil health and crop productivity. The role of prokaryotes, eukaryotes and their phylogeny, physiology, and genetic regulations are essential for understanding the ecological significance of this empirical process. Several prokaryotes (bacteria and archaea), and a few eukaryotes (fungi and algae) are reported as NO3- reducers under certain conditions. This process involves enzymatic reactions which has been catalysed by nitrate reductases, nitrite reductases, and NH4+-assimilating enzymes. Earlier reports emphasised that single-cell prokaryotic or eukaryotic organisms are responsible for this process, which portrayed a prominent gap. Therefore, this study revisits the similarities and uniqueness of mechanism behind NO3- -reduction to NH4+ in both prokaryotes and eukaryotes. Moreover, phylogenetic, physiological, and genetic regulation also shed light on the evolutionary connections between two systems which could help us to better explain the NO3--reduction mechanisms over time. Reports also revealed that certain transcription factors like NtrC/NtrB and Nit2 have shown a major role in coordinating the expression of NO3- assimilation genes in response to NO3- availability. Overall, this review provides a comprehensive information about the complex fermentative and respiratory dissimilatory nitrate reduction to ammonium (DNRA) processes. Uncovering the complexity of this process across various organisms may further give insight into sustainable nitrogen management practices and might contribute to addressing global environmental challenges.

Asgard archaeal selenoproteome reveals a roadmap for the archaea-to-eukaryote transition of selenocysteine incorporation machinery.

Selenocysteine (Sec) is encoded by the UGA codon that normally functions as a stop signal and is specifically incorporated into selenoproteins via a unique recoding mechanism. The translational recoding of UGA as Sec is directed by an unusual RNA structure, the SECIS element. Although archaea and eukaryotes adopt similar Sec encoding machinery, the SECIS elements have no similarities to each other with regard to sequence and structure. We analyzed >400 Asgard archaeal genomes to examine the occurrence of both Sec encoding system and selenoproteins in this archaeal superphylum, the closest prokaryotic relatives of eukaryotes. A comprehensive map of Sec utilization trait has been generated, providing the most detailed understanding of the use of this nonstandard amino acid in Asgard archaea so far. By characterizing the selenoproteomes of all organisms, several selenoprotein-rich phyla and species were identified. Most Asgard archaeal selenoprotein genes possess eukaryotic SECIS-like structures with varying degrees of diversity. Moreover, euryarchaeal SECIS elements might originate from Asgard archaeal SECIS elements via lateral gene transfer, indicating a complex and dynamic scenario of the evolution of SECIS element within archaea. Finally, a roadmap for the transition of eukaryotic SECIS elements from archaea was proposed, and selenophosphate synthetase may serve as a potential intermediate for the generation of ancestral eukaryotic SECIS element. Our results offer new insights into a deeper understanding of the evolution of Sec insertion machinery.

EUKARYOME: the rRNA gene reference database for identification of all eukaryotes.

Molecular identification of micro- and macroorganisms based on nuclear markers has revolutionized our understanding of their taxonomy, phylogeny and ecology. Today, research on the diversity of eukaryotes in global ecosystems heavily relies on nuclear ribosomal RNA (rRNA) markers. Here, we present the research community-curated reference database EUKARYOME for nuclear ribosomal 18S rRNA, internal transcribed spacer (ITS) and 28S rRNA markers for all eukaryotes, including metazoans (animals), protists, fungi and plants. It is particularly useful for the identification of arbuscular mycorrhizal fungi as it bridges the four commonly used molecular markers-ITS1, ITS2, 18S V4-V5 and 28S D1-D2 subregions. The key benefits of this database over other annotated reference sequence databases are that it is not restricted to certain taxonomic groups and it includes all rRNA markers. EUKARYOME also offers a number of reference long-read sequences that are derived from (meta)genomic and (meta)barcoding-a unique feature that can be used for taxonomic identification and chimera control of third-generation, long-read, high-throughput sequencing data. Taxonomic assignments of rRNA genes in the database are verified based on phylogenetic approaches. The reference datasets are available in multiple formats from the project homepage, http://www.eukaryome.org.

GeneMark-ETP significantly improves the accuracy of automatic annotation of large eukaryotic genomes.

Large-scale genomic initiatives, such as the Earth BioGenome Project, require efficient methods for eukaryotic genome annotation. Here we present an automatic gene finder, GeneMark-ETP, integrating genomic-, transcriptomic-, and protein-derived evidence that has been developed with a focus on large plant and animal genomes. GeneMark-ETP first identifies genomic loci where extrinsic data are sufficient for making gene predictions with "high confidence." The genes situated in the genomic space between the high-confidence genes are predicted in the next stage. The set of high-confidence genes serves as an initial training set for the statistical model. Further on, the model parameters are iteratively updated in the rounds of gene prediction and parameter re-estimation. Upon reaching convergence, GeneMark-ETP makes the final predictions and delivers the whole complement of predicted genes. GeneMark-ETP outperforms gene finders using a single type of extrinsic evidence. Comparisons with gene finders MAKER2 and TSEBRA, those that use both transcript- and protein-derived extrinsic evidence, show that GeneMark-ETP delivers state-of-the-art gene-prediction accuracy, with the margin of outperforming existing approaches increasing in its application to larger and more complex eukaryotic genomes.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments.

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

mRNA context and translation factors determine decoding in alternative nuclear genetic codes.

The genetic code is a set of instructions that determine how the information in our genetic material is translated into amino acids. In general, it is universal for all organisms, from viruses and bacteria to humans. However, in the last few decades, exceptions to this rule have been identified both in pro- and eukaryotes. In this review, we discuss the 16 described alternative eukaryotic nuclear genetic codes and observe theories of their appearance in evolution. We consider possible molecular mechanisms that allow codon reassignment. Most reassignments in nuclear genetic codes are observed for stop codons. Moreover, in several organisms, stop codons can simultaneously encode amino acids and serve as termination signals. In this case, the meaning of the codon is determined by the additional factors besides the triplets. A comprehensive review of various non-standard coding events in the nuclear genomes provides a new insight into the translation mechanism in eukaryotes.

Eukaryotic genomic data uncover an extensive host range of mirusviruses.

A recent marine metagenomic study has revealed the existence of a novel group of viruses designated mirusviruses, which are proposed to form an evolutionary link between two realms of double-stranded DNA viruses, Varidnaviria and Duplodnaviria. Metagenomic data suggest that mirusviruses infect microeukaryotes in the photic layer of the ocean, but their host range remains largely unknown. In this study, we investigated the presence of mirusvirus marker genes in 1,901 publicly available eukaryotic genome assemblies, mainly derived from unicellular eukaryotes, to identify potential hosts of mirusviruses. Mirusvirus marker sequences were identified in 915 assemblies spanning 227 genera across eight supergroups of eukaryotes. The habitats of the putative mirusvirus hosts included not only marine but also other diverse environments. Among the major capsid protein (MCP) signals in the genome assemblies, we identified 85 sequences that showed high sequence and structural similarities to reference mirusvirus MCPs. A phylogenetic analysis of these sequences revealed their distant evolutionary relationships with the seven previously reported mirusvirus clades. Most of the scaffolds with these MCP sequences encoded multiple mirusvirus homologs, suggesting that mirusviral infection contributes to the alteration of the host genome. We also identified three circular mirusviral genomes within the genomic data of the oil-producing thraustochytrid Schizochytrium sp. and the endolithic green alga Ostreobium quekettii. Overall, mirusviruses probably infect a wide spectrum of eukaryotes and are more diverse than previously reported.

Connecting Species-Specific Extents of Genome Reduction in Mitochondria and Plastids.

Mitochondria and plastids have both dramatically reduced their genomes since the endosymbiotic events that created them. The similarities and differences in the evolution of the two organelle genome types have been the target of discussion and investigation for decades. Ongoing work has suggested that similar mechanisms may modulate the reductive evolution of the two organelles in a given species, but quantitative data and statistical analyses exploring this picture remain limited outside of some specific cases like parasitism. Here, we use cross-eukaryote organelle genome data to explore evidence for coevolution of mitochondrial and plastid genome reduction. Controlling for differences between clades and pseudoreplication due to relatedness, we find that extents of mtDNA and ptDNA gene retention are related to each other across taxa, in a generally positive correlation that appears to differ quantitatively across eukaryotes, for example, between algal and nonalgal species. We find limited evidence for coevolution of specific mtDNA and ptDNA gene pairs, suggesting that the similarities between the two organelle types may be due mainly to independent responses to consistent evolutionary drivers.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes.

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Natural history of eukaryotic DNA viruses with double jelly-roll major capsid proteins.

The phylum Preplasmiviricota (kingdom Bamfordvirae, realm Varidnaviria) is a broad assemblage of diverse viruses with comparatively short double-stranded DNA genomes (<50 kbp) that produce icosahedral capsids built from double jelly-roll major capsid proteins. Preplasmiviricots infect hosts from all cellular domains, testifying to their ancient origin, and, in particular, are associated with six of the seven supergroups of eukaryotes. Preplasmiviricots comprise four major groups of viruses, namely, polintons, polinton-like viruses (PLVs), virophages, and adenovirids. We used protein structure modeling and analysis to show that protein-primed DNA polymerases (pPolBs) of polintons, virophages, and cytoplasmic linear plasmids encompass an N-terminal domain homologous to the terminal proteins (TPs) of prokaryotic PRD1-like tectivirids and eukaryotic adenovirids that are involved in protein-primed replication initiation, followed by a viral ovarian tumor-like cysteine deubiquitinylase (vOTU) domain. The vOTU domain is likely responsible for the cleavage of the TP from the large pPolB polypeptide and is inactivated in adenovirids, in which TP is a separate protein. Many PLVs and transpovirons encode a distinct derivative of polinton-like pPolB that retains the TP, vOTU, and pPolB polymerization palm domains but lacks the exonuclease domain and instead contains a superfamily 1 helicase domain. Analysis of the presence/absence and inactivation of the vOTU domains and replacement of pPolB with other DNA polymerases in eukaryotic preplasmiviricots enabled us to outline a complete scenario for their origin and evolution.

Ribosome decision graphs for the representation of eukaryotic RNA translation complexity.

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, within both annotated protein-coding and noncoding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term ribosome decision graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the latter "translons." Nondeterministic events, such as initiation, reinitiation, selenocysteine insertion, or ribosomal frameshifting, are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions and for analyzing genetic variation and quantitative genome-wide data on translation for characterization of regulatory modulators of translation.

Dynamic evolution of the mTHF gene family associated with primary metabolism across life.

BACKGROUND: The folate cycle of one-carbon (C1) metabolism, which plays a central role in the biosynthesis of nucleotides and amino acids, demonstrates the significance of metabolic adaptation. We investigated the evolutionary history of the methylenetetrahydrofolate dehydrogenase (mTHF) gene family, one of the main drivers of the folate cycle, across life. RESULTS: Through comparative genomic and phylogenetic analyses, we found that several lineages of Archaea lacked domains vital for folate cycle function such as the mTHF catalytic and NAD(P)-binding domains of FolD. Within eukaryotes, the mTHF gene family diversified rapidly. For example, several duplications have been observed in lineages including the Amoebozoa, Opisthokonta, and Viridiplantae. In a common ancestor of Opisthokonta, FolD and FTHFS underwent fusion giving rise to the gene MTHFD1, possessing the domains of both genes. CONCLUSIONS: Our evolutionary reconstruction of the mTHF gene family associated with a primary metabolic pathway reveals dynamic evolution, including gene birth-and-death, gene fusion, and potential horizontal gene transfer events and/or amino acid convergence.

18S rDNA sequence-structure phylogeny of the eukaryotes simultaneously inferred from sequences and their individual secondary structures.

OBJECTIVE: The eukaryotic tree of life has been subject of numerous studies ever since the nineteenth century, with more supergroups and their sister relations being decoded in the last years. In this study, we reconstructed the phylogeny of eukaryotes using complete 18S rDNA sequences and their individual secondary structures simultaneously. After the sequence-structure data was encoded, it was automatically aligned and analyzed using sequence-only as well as sequence-structure approaches. We present overall neighbor-joining trees of 211 eukaryotes as well as the respective profile neighbor-joining trees, which helped to resolve the basal branching pattern. A manually chosen subset was further inspected using neighbor-joining, maximum parsimony, and maximum likelihood analyses. Additionally, the 75 and 100 percent consensus structures of the subset were predicted. RESULTS: All sequence-structure approaches show improvements compared to the respective sequence-only approaches: the average bootstrap support per node of the sequence-structure profile neighbor-joining analyses with 90.3, was higher than the average bootstrap support of the sequence-only profile neighbor-joining analysis with 73.9. Also, the subset analyses using sequence-structure data were better supported. Furthermore, more subgroups of the supergroups were recovered as monophyletic and sister group relations were much more comparable to results as obtained by multi-marker analyses.

Macroalgal eukaryotic microbiome composition indicates novel phylogenetic diversity and broad host spectrum of oomycete pathogens.

Seaweeds are important components of marine ecosystems with emerging potential in aquaculture and as sources of biofuel, food products and pharmacological compounds. However, an increasingly recognised threat to natural and industrial seaweed populations is infection with parasitic single-celled eukaryotes from the relatively understudied oomycete lineage. Here we examine the eukaryomes of diverse brown, red and green marine macroalgae collected from polar (Baffin Island), cold-temperate (Falkland Islands) and tropical (Ascension Island) locations, with a focus on oomycete and closely related diatom taxa. Using 18S rRNA gene amplicon sequencing, we show unexpected genetic and taxonomic diversity of the eukaryomes, a strong broad-brush association between eukaryome composition and geographic location, and some evidence of association between eukaryome structure and macroalgal phylogenetic relationships (phylosymbiosis). However, the oomycete fraction of the eukaryome showed disparate patterns of diversity and structure, highlighting much weaker association with geography and no evidence of phylosymbiosis. We present several novel haplotypes of the most common oomycete Eurychasma dicksonii and report for the first time a cosmopolitan distribution and absence of host specificity of this important pathogen. This indicates rich diversity in macroalgal oomycete pathogens and highlights that these pathogens may be generalist and highly adaptable to diverse environmental conditions.