Genome-wide signatures of adaptation to extreme environments in red algae

Two high-quality Prototheca zopfii genomes provide new insights into their evolution as obligate algal heterotrophs and their pathogenicity

The majority of the nearly 10,000 described species of green algae are photoautotrophs; however, some species have lost their ability to photosynthesize and become obligate heterotrophs that rely on parasitism for survival. Two high-quality genomes of the heterotrophic algae Prototheca zopfii Pz20 and Pz23 were obtained using short- and long-read genomic as well as transcriptomic data. The genome sizes were 31.2 Mb and 31.3 Mb, respectively, and contig N50 values of 1.99 Mb and 1.26 Mb. Although P. zopfii maintained its plastid genome, the transition to heterotrophy led to a reduction in both plastid and nuclear genome size, including the loss of photosynthesis-related genes from both the nuclear and plastid genomes and the elimination of genes encoding for carotenoid oxygenase and pheophorbide an oxygenase. The loss of genes, including basic leucine-zipper (bZIP) transcription factors, flavin adenine dinucleotide-linked oxidase, and helicase, could have played a role in the transmission of autotrophy to heterotrophs and in the processes of abiotic stress resistance and pathogenicity. A total of 66 (1.37%) and 73 (1.49%) genes were identified as potential horizontal gene transfer events in the two P. zopfii genomes, respectively. Genes for malate synthase and isocitrate lyase, which are horizontally transferred from bacteria, may play a pivotal role in carbon and nitrogen metabolism as well as the pathogenicity of Prototheca and non-photosynthetic organisms. The two high-quality P. zopfii genomes provide new insights into their evolution as obligate heterotrophs and pathogenicity.

IMPORTANCE

The genus Prototheca, characterized by its heterotrophic nature and pathogenicity, serves as an exemplary model for investigating pathobiology. The limited understanding of the protothecosis infectious disease is attributed to the lack of genomic resources. Using HiFi long-read sequencing, both nuclear and plastid genomes were generated for two strains of P. zopfii. The findings revealed a concurrent reduction in both plastid and nuclear genome size, accompanied by the loss of genes associated with photosynthesis, carotenoid oxygenase, basic leucine-zipper (bZIP) transcription factors, and others. The analysis of horizontal gene transfer revealed the presence of 1.37% and 1.49% bacterial genes, including malate synthase and isocitrate lyase, which play crucial roles in carbon and nitrogen metabolism, as well as pathogenicity and obligate heterotrophy. The two high-quality P. zopfii genomes represent valuable resources for investigating their adaptation and evolution as obligate heterotrophs, as well as for developing future prevention and treatment strategies against protothecosis.

Methylmercury exposure of the sponge O. lobularis induces strong tissue and cell defects

Mediterranean marine biota suffers from various anthropogenic threats. Among them, pollutants such as mercury (Hg) represent important environmental issues that are exacerbated by bioaccumulation and bioamplification along food webs via its organic form, monomethylmercury (MMHg). To date, very little is known regarding the impact of mercury on Porifera and the few available studies have been exclusively focused on Demospongiae. This work studies the effect of MMHgCl at different biological levels of Oscarella lobularis (Porifera, Homoscleromorpha). Bioaccumulation assays show that MMHgCl significantly accumulated in sponge tissues after a 96-h exposure to 0.1 μg L-1. Toxicity assays (LC5096h) show a sensibility that depends on life-stage (adult vs bud). Additionally, we show that the exposure to 1 μg L-1 MMHgCl negatively impacts the epithelial integrity and the regeneration process in buds, as shown by the loss of cell-cell contacts and the alteration of osculum morphogenesis. For the first time in a sponge, a whole set of genes classically involved in metal detoxification and in antioxidant response were identified. Significant changes in catalase, superoxide dismutase and nuclear factor (erythroid-derived 2)-like 2 expressions in exposed juveniles were measured. Such an integrative approach from the physiological to the molecular scales on a non-model organism expands our knowledge concerning sensitivity and toxicity mechanisms induced by MMHg in Porifera, raising new questions regarding the possible defences used by marine sponges.

Mechanisms of adaptation and evolution in Toxoplasma gondii

Toxoplasma has high host flexibility, infecting all nucleated cells of mammals and birds. This implies that during its infective process the parasite must constantly adapt to different environmental situations, which in turn leads to modifications in its metabolism, regulation of gene transcription, translation of mRNAs and stage specific factors. There are conserved pathways that support these adaptations, which we aim to elucidate in this review. We begin by exploring the widespread epigenetic mechanisms and transcription regulators, continue with the supportive role of Heat Shock Proteins (Hsp), the translation regulation, stress granules, and finish with the emergence of contingency genes in highly variable genomic domains, such as subtelomeres. Within epigenetics, the discovery of a new histone variant of the H2B family (H2B.Z), contributing to T. gondii virulence and differentiation, but also gene expression regulation and its association with the metabolic state of the parasite, is highlighted. Associated with the regulation of gene expression are transcription factors (TFs). An overview of the main findings on TF and development is presented. We also emphasize the role of Hsp90 and Tgj1 in T. gondii metabolic fitness and the regulation of protein translation. Translation regulation is also highlighted as a mechanism for adaptation to conditions encountered by the parasite as well as stress granules containing mRNA and proteins generated in the extracellular tachyzoite. Another important aspect in evolution and adaptability are the subtelomeres because of their high variability and gene duplication rate. Toxoplasma possess multigene families of membrane proteins and contingency genes that are associated with different metabolic stresses. Among them parasite differentiation and environmental stresses stand out, including those that lead tachyzoite to bradyzoite conversion. Finally, we are interested in positioning protozoa as valuable evolution models, focusing on research related to the Extended Evolutionary Synthesis, based on models recently generated, such as extracellular adaptation and ex vivo cyst recrudescence.

Accumulation and bio-oxidation of arsenite mediated by thermoacidophilic Cyanidiales: innate potential biomaterials toward arsenic remediation

Addressing geogenic and anthropogenic arsenic (As) pollution is critical for environmental health. This study explored arsenite [As(III)] removal using Cyanidiales, particularly Cyanidium caldarium (Cc) and Galdieria partita (Gp), under acidic to neutral pH, and determined As(III) detoxification mechanisms in relation to As speciation and protein secondary structure in Cyanidiales. Regarding As(III) sorption amounts, Cc outperformed Gp, reaching 83.2 mg g-1 of removal at pH 5.0. Wherein, 23.5 % of sorbed As on Cc presented as arsenate [As(V)] complexation with polysaccharides, alongside other predominant species including As(III)-cysteine (41.2 %) and As(III)-polysaccharides (35.3 %) complexes. This suggested that As(III) was directly transported into cells, rather than As(V). Coupled with the formation of As(III)-cysteine complexes within cells, these mechanisms may be key to efficiently accumulating As(III) in Cyanidiales during the 6-h incubation. These results highlight the potential of Cyanidiales for sustainable As(III) remediation and provide new insights into managing As(III) toxicity.

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

Effects of benzophenone-3 and its metabolites on the marine diatom Chaetoceros neogracilis: Underlying mechanisms and environmental implications

The wide application of benzophenones (BPs), such as benzophenone-3 (BP3), as an ingredient in sunscreens, cosmetics, coatings, and plastics, has led to their global contamination in aquatic environments. Using the marine diatom Chaetoceros neogracilis as a model, this study assessed the toxic effects and mechanisms of BP3 and its two major metabolites (BP8 and BP1). The results showed that BP3 exhibited higher toxicity on C. neogracilis than BP8 and BP1, with their 72-h median effective concentrations being 0.4, 0.8 and 4 mg/L, respectively. Photosynthesis efficiencies were significantly reduced after exposure to environmentally relevant concentrations of the three benzophenones, while cell viability, membrane integrity, membrane potential, and metabolic activities could be further impaired at their higher concentrations. Comparative transcriptomic analysis, followed by gene ontology and KEGG pathway enrichment analyses unraveled that all the three tested benzophenones disrupted photosynthesis and nitrogen metabolism of the diatom through alteration of similar pathways. The toxic effect of BP3 was also attributable to its unique inhibitory effects on eukaryotic ribosome biosynthesis and DNA replication. Taken together, our findings underscore that benzophenones may pose a significant threat to photosynthesis, oxygen production, primary productivity, carbon fixation, and the nitrogen cycle of diatom in coastal waters worldwide.

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Dollo Parsimony Overestimates Ancestral Gene Content Reconstructions

Ancestral reconstruction is a widely used technique that has been applied to understand the evolutionary history of gain and loss of gene families. Ancestral gene content can be reconstructed via different phylogenetic methods, but many current and previous studies employ Dollo parsimony. We hypothesize that Dollo parsimony is not appropriate for ancestral gene content reconstruction inferences based on sequence homology, as Dollo parsimony is derived from the assumption that a complex character cannot be regained. This premise does not accurately model molecular sequence evolution, in which false orthology can result from sequence convergence or lateral gene transfer. The aim of this study is to test Dollo parsimony's suitability for ancestral gene content reconstruction and to compare its inferences with a maximum likelihood-based approach that allows a gene family to be gained more than once within a tree. We first compared the performance of the two approaches on a series of artificial data sets each of 5,000 genes that were simulated according to a spectrum of evolutionary rates without gene gain or loss, so that inferred deviations from the true gene count would arise only from errors in orthology inference and ancestral reconstruction. Next, we reconstructed protein domain evolution on a phylogeny representing known eukaryotic diversity. We observed that Dollo parsimony produced numerous ancestral gene content overestimations, especially at nodes closer to the root of the tree. These observations led us to the conclusion that, confirming our hypothesis, Dollo parsimony is not an appropriate method for ancestral reconstruction studies based on sequence homology.

Identification of plants' functional counterpart of the metazoan mediator of DNA Damage checkpoint 1

Induction of DNA damage triggers rapid phosphorylation of the histone H2A.X (γH2A.X). In animals, mediator of DNA damage checkpoint 1 (MDC1) binds γH2A.X through a tandem BRCA1 carboxyl-terminal (tBRCT) domain and mediates recruitment of downstream effectors of DNA damage response (DDR). However, readers of this modification in plants have remained elusive. We show that from the Arabidopsis BRCT domain proteome, BCP1-4 proteins with tBRCT domains are involved in DDR. Through its tBRCT domain BCP4 binds γH2A.X in vitro and localizes to DNA damage-induced foci in an H2A.X-dependent manner. BCP4 also contains a domain that interacts directly with NBS1 and thus acts as a functional counterpart of MDC1. We also show that BCP1, that contains two tBRCT domains, co-localizes with γH2A.X but it does not bind γH2A.X suggesting functional similarity with human PAXIP1. A phylogenetic analysis supports that PAXIP1 and MDC1 in metazoa and their plant counterparts evolved independently from common ancestors with tBRCT domains. Collectively, our study reveals missing components and provides mechanistic and evolutionary insights into plant DDR.

The structure of PSI-LHCI from Cyanidium caldarium provides evolutionary insights into conservation and diversity of red-lineage LHCs

Light-harvesting complexes (LHCs) are diversified among photosynthetic organisms, and the structure of the photosystem I-LHC (PSI-LHCI) supercomplex has been shown to be variable depending on the species of organisms. However, the structural and evolutionary correlations of red-lineage LHCs are unknown. Here, we determined a 1.92-Å resolution cryoelectron microscopic structure of a PSI-LHCI supercomplex isolated from the red alga Cyanidium caldarium RK-1 (NIES-2137), which is an important taxon in the Cyanidiophyceae. We subsequently investigated the correlations of PSI-LHCIs from different organisms through structural comparisons and phylogenetic analysis. The PSI-LHCI structure obtained shows five LHCI subunits surrounding a PSI-monomer core. The five LHCIs are composed of two Lhcr1s, two Lhcr2s, and one Lhcr3. Phylogenetic analysis of LHCs bound to PSI in the red-lineage algae showed clear orthology of LHCs between C. caldarium and Cyanidioschyzon merolae, whereas no orthologous relationships were found between C. caldarium Lhcr1-3 and LHCs in other red-lineage PSI-LHCI structures. These findings provide evolutionary insights into conservation and diversity of red-lineage LHCs associated with PSI.

Rapid development of increased neonicotinoid tolerance in non-target freshwater amphipods

The comprehensive assessment of the long-term impacts of constant exposure to pollutants on wildlife populations remains a relatively unexplored area of ecological risk assessment. Empirical evidence to suggest that multigenerational exposure affects the susceptibility of organisms is scarce, and the underlying mechanisms in the natural environment have yet to be fully understood. In this study, we first examined the arthropod candidate species, Gammarus roeselii that - unlike closely related species - commonly occurs in many contaminated river systems of Central Europe. This makes it a suitable study organism to investigate the development of tolerances and phenotypic adaptations along pollution gradients. In a 96-h acute toxicity assay with the neonicotinoid thiacloprid, we indeed observed a successive increase in tolerance in populations coming from contaminated regions. This was accompanied by a certain phenotypic change, with increased investment into reproduction. To address the question of whether these changes are plastic or emerged from longer lasting evolutionary processes, we conducted a multigeneration experiment in the second part of our study. Here, we used closely-related Hyalella azteca and pre-exposed them for multiple generations to sublethal concentrations of thiacloprid in a semi-static design (one week renewal of media containing 0.1 or 1.0 µg/L thiacloprid). The pre-exposed individuals were then used in acute toxicity assays to see how quickly such adaptive responses can develop. Over only two generations, the tolerance to the neonicotinoid almost doubled, suggesting developmental plasticity as a plausible mechanism for the rapid adaptive response to strong selection factors such as neonicotinoid insecticides. It remains to be discovered whether the plasticity of rapidly developed tolerance is species-specific and explains why closely related species - which may not have comparable adaptive response capabilities - disappear in polluted habitats. Overall, our findings highlight the neglected role of developmental plasticity during short- and long-term exposure of natural populations to pollution. Moreover, our results show that even pollutant levels seven times lower than concentrations found in the study region have a clear impact on the developmental trajectories of non-target species.

A k-mer-Based Approach for Phylogenetic Classification of Taxa in Environmental Genomic Data

In the age of genome sequencing, whole-genome data is readily and frequently generated, leading to a wealth of new information that can be used to advance various fields of research. New approaches, such as alignment-free phylogenetic methods that utilize k-mer-based distance scoring, are becoming increasingly popular given their ability to rapidly generate phylogenetic information from whole-genome data. However, these methods have not yet been tested using environmental data, which often tends to be highly fragmented and incomplete. Here, we compare the results of one alignment-free approach (which utilizes the D2 statistic) to traditional multi-gene maximum likelihood trees in 3 algal groups that have high-quality genome data available. In addition, we simulate lower-quality, fragmented genome data using these algae to test method robustness to genome quality and completeness. Finally, we apply the alignment-free approach to environmental metagenome assembled genome data of unclassified Saccharibacteria and Trebouxiophyte algae, and single-cell amplified data from uncultured marine stramenopiles to demonstrate its utility with real datasets. We find that in all instances, the alignment-free method produces phylogenies that are comparable, and often more informative, than those created using the traditional multi-gene approach. The k-mer-based method performs well even when there are significant missing data that include marker genes traditionally used for tree reconstruction. Our results demonstrate the value of alignment-free approaches for classifying novel, often cryptic or rare, species, that may not be culturable or are difficult to access using single-cell methods, but fill important gaps in the tree of life.

The Polycomb repressive complex 2 deposits H3K27me3 and represses transposable elements in a broad range of eukaryotes

The mobility of transposable elements (TEs) contributes to evolution of genomes. Their uncontrolled activity causes genomic instability; therefore, expression of TEs is silenced by host genomes. TEs are marked with DNA and H3K9 methylation, which are associated with silencing in flowering plants, animals, and fungi. However, in distantly related groups of eukaryotes, TEs are marked by H3K27me3 deposited by the Polycomb repressive complex 2 (PRC2), an epigenetic mark associated with gene silencing in flowering plants and animals. The direct silencing of TEs by PRC2 has so far only been shown in one species of ciliates. To test if PRC2 silences TEs in a broader range of eukaryotes, we generated mutants with reduced PRC2 activity and analyzed the role of PRC2 in extant species along the lineage of Archaeplastida and in the diatom P. tricornutum. In this diatom and the red alga C. merolae, a greater proportion of TEs than genes were repressed by PRC2, whereas a greater proportion of genes than TEs were repressed by PRC2 in bryophytes. In flowering plants, TEs contained potential cis-elements recognized by transcription factors and associated with neighbor genes as transcriptional units repressed by PRC2. Thus, silencing of TEs by PRC2 is observed not only in Archaeplastida but also in diatoms and ciliates, suggesting that PRC2 deposited H3K27me3 to silence TEs in the last common ancestor of eukaryotes. We hypothesize that during the evolution of Archaeplastida, TE fragments marked with H3K27me3 were selected to shape transcriptional regulation, controlling networks of genes regulated by PRC2.

Red macroalgae in the genomic era

Rhodophyta (or red algae) are a diverse and species-rich group that forms one of three major lineages in the Archaeplastida, a eukaryotic supergroup whose plastids arose from a single primary endosymbiosis. Red algae are united by several features, such as relatively small intron-poor genomes and a lack of cytoskeletal structures associated with motility like flagella and centrioles, as well as a highly efficient photosynthetic capacity. Multicellular red algae (or macroalgae) are one of the earliest diverging eukaryotic lineages to have evolved complex multicellularity, yet despite their ecological, evolutionary, and commercial importance, they have remained a largely understudied group of organisms. Considering the increasing availability of red algal genome sequences, we present a broad overview of fundamental aspects of red macroalgal biology and posit on how this is expected to accelerate research in many domains of red algal biology in the coming years.

Origin of minicircular mitochondrial genomes in red algae

Eukaryotic organelle genomes are generally of conserved size and gene content within phylogenetic groups. However, significant variation in genome structure may occur. Here, we report that the Stylonematophyceae red algae contain multipartite circular mitochondrial genomes (i.e., minicircles) which encode one or two genes bounded by a specific cassette and a conserved constant region. These minicircles are visualized using fluorescence microscope and scanning electron microscope, proving the circularity. Mitochondrial gene sets are reduced in these highly divergent mitogenomes. Newly generated chromosome-level nuclear genome assembly of Rhodosorus marinus reveals that most mitochondrial ribosomal subunit genes are transferred to the nuclear genome. Hetero-concatemers that resulted from recombination between minicircles and unique gene inventory that is responsible for mitochondrial genome stability may explain how the transition from typical mitochondrial genome to minicircles occurs. Our results offer inspiration on minicircular organelle genome formation and highlight an extreme case of mitochondrial gene inventory reduction.

Revised classification of the Cyanidiophyceae based on plastid genome data with descriptions of the Cavernulicolales ord. nov. and Galdieriales ord. nov. (Rhodophyta)

The Cyanidiophyceae, an extremophilic red algal class, is distributed worldwide in extreme environments. Species grow either in acidic hot environments or in dim light conditions (e.g., "cave Cyanidium"). The taxonomy and classification systems are currently based on morphological, eco-physiological, and molecular phylogenetic characters; however, previous phylogenetic results showed hidden diversity of the Cyanidiophyceae and suggested a revision of the classification system. To clarify phylogenetic relationships within this red algal class, we employ a phylogenomic approach based on 15 plastomes (10 new) and 15 mitogenomes (seven new). Our phylogenies show consistent relationships among four lineages (Galdieria, "cave Cyanidium", Cyanidium, and Cyanidioschyzon lineages). Each lineage is distinguished by organellar genome characteristics. The "cave Cyanidium" lineage is a distinct clade that diverged after the Galdieria clade but within a larger monophyletic clade that included the Cyanidium and Cyanidioschyzon lineages. Because the "cave Cyanidium" lineage is a mesophilic lineage that differs substantially from the other three thermoacidophilic lineages, we describe it as a new order (Cavernulicolales). Based on this evidence, we reclassified the Cyanidiophyceae into four orders: Cyanidiales, Cyanidioschyzonales, Cavernulicolales ord. nov., and Galdieriales ord. nov. The genetic distance among these four orders is comparable to, or greater than, the distances found between other red algal orders and subclasses. Three new genera (Cavernulicola, Gronococcus, Sciadococcus), five new species (Galdieria javensis, Galdieria phlegrea, Galdieria yellowstonensis, Gronococcus sybilensis, Sciadococcus taiwanensis), and a new nomenclatural combination (Cavernulicola chilensis) are proposed.

Algae obscura: The potential of rare species as model systems

Model organism research has provided invaluable knowledge about foundational biological principles. However, most of these studies have focused on species that are in high abundance, easy to cultivate in the lab, and represent only a small fraction of extant biodiversity. Here, we present three examples of rare algae with unusual features that we refer to as "algae obscura." The Cyanidiophyceae (Rhodophyta), Glaucophyta, and Paulinella (rhizarian) lineages have all transitioned out of obscurity to become models for fundamental evolutionary research. Insights have been gained into the prevalence and importance of eukaryotic horizontal gene transfer, early Earth microbial community dynamics, primary plastid endosymbiosis, and the origin of Archaeplastida. By reviewing the research that has come from the exploration of these organisms, we demonstrate that underappreciated algae have the potential to help us formulate, refine, and substantiate core hypotheses and that such organisms should be considered when establishing future model systems.

Eukaryotes contribute more than bacteria to the recovery of freshwater ecosystem functions under different drought durations

Global climate change mostly impacts river ecosystems by affecting microbial biodiversity and ecological functions. Considering the high functional redundancy of microorganisms, the unknown relationship between biodiversity and ecosystem functions obstructs river ecological research, especially under the influence of increasing weather extremes, such as in intermittent rivers and ephemeral streams (IRES). Herein, dry-wet alternation experiments were conducted in artificial stream channels for 25 and 90 days of drought, both followed by 20 days of rewetting. The dynamic recovery of microbial biodiversity and ecosystem functions (represented by ecosystem metabolism and denitrification rate) were determined to analyse biodiversity-ecosystem-function (BEF) relationships after different drought durations. There was a significant difference between bacterial and eukaryotic biodiversity recovery after drought. Eukaryotic biodiversity was more sensitive to drought duration than bacterial, and the eukaryotic network was more stable under dry-wet alternations. Based on the establishment of partial least squares path models, we found that eukaryotic biodiversity has a stronger effect on ecosystem functions than bacteria after long-term drought. Indeed, this work represents a significant step forward for further research on the ecosystem functions of IRES, especially emphasizing the importance of eukaryotic biodiversity in the BEF relationship.