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.

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.

How do bacterial endosymbionts work with so few genes?

The move from a free-living environment to a long-term residence inside a host eukaryotic cell has profound effects on bacterial function. While endosymbioses are found in many eukaryotes, from protists to plants to animals, the bacteria that form these host-beneficial relationships are even more diverse. Endosymbiont genomes can become radically smaller than their free-living relatives, and their few remaining genes show extreme compositional biases. The details of how these reduced and divergent gene sets work, and how they interact with their host cell, remain mysterious. This Unsolved Mystery reviews how genome reduction alters endosymbiont biology and highlights a "tipping point" where the loss of the ability to build a cell envelope coincides with a marked erosion of translation-related genes.

DeepReg: a deep learning hybrid model for predicting transcription factors in eukaryotic and prokaryotic genomes.

Deep learning models (DLMs) have gained importance in predicting, detecting, translating, and classifying a diversity of inputs. In bioinformatics, DLMs have been used to predict protein structures, transcription factor-binding sites, and promoters. In this work, we propose a hybrid model to identify transcription factors (TFs) among prokaryotic and eukaryotic protein sequences, named Deep Regulation (DeepReg) model. Two architectures were used in the DL model: a convolutional neural network (CNN), and a bidirectional long-short-term memory (BiLSTM). DeepReg reached a precision of 0.99, a recall of 0.97, and an F1-score of 0.98. The quality of our predictions, the bias-variance trade-off approach, and the characterization of new TF predictions were evaluated and compared against those produced by DeepTFactor, as well as against experimental data from three model organisms. Predictions based on our DLM tended to exhibit less variance and bias than those from DeepTFactor, thus increasing reliability and decreasing overfitting.

Blessing or curse: how the epigenetic resolution of host-transposable element conflicts shapes their evolutionary dynamics.

Transposable elements (TEs) are selfish genetic elements whose antagonistic interactions with hosts represent a common genetic conflict in eukaryotes. To resolve this conflict, hosts have widely adopted epigenetic silencing that deposits repressive marks at TEs. However, this mechanism is imperfect and fails to fully halt TE replication. Furthermore, TE epigenetic silencing can inadvertently spread repressive marks to adjacent functional sequences, a phenomenon considered a 'curse' of this conflict resolution. Here, we used forward simulations to explore how TE epigenetic silencing and its harmful side effects shape the evolutionary dynamics of TEs and their hosts. Our findings reveal that epigenetic silencing allows TEs and their hosts to stably coexist under a wide range of conditions, because the underlying molecular mechanisms give rise to copy-number dependency of the strength of TE silencing. Interestingly, contrary to intuitive expectations that TE epigenetic silencing should evolve to be as strong as possible, we found a selective benefit for modifier alleles that weaken TE silencing under biologically feasible conditions. These results reveal that the dual nature of TE epigenetic silencing, with both positive and negative effects, complicates its evolutionary trajectory and makes it challenging to determine whether TE epigenetic silencing is a 'blessing' or a 'curse'.

Advances in the mechanisms and applications of RNA silencing in crop protection.

RNA silencing (or RNA interference, RNAi) is a conserved mechanism for regulating gene expression in eukaryotes, which plays vital roles in plant development and response to biotic and abiotic stresses. The discovery of trans-kingdom RNAi and interspecies RNAi provides a theoretical basis for exploiting RNAi-based crop protection strategies. Here, we summarize the canonical RNAi mechanisms in plants and review representative studies associated with plant-pathogen interactions. Meanwhile, we also elaborate upon the principles of host-induced gene silencing, spray-induced gene silencing and microbe-induced gene silencing, and discuss their applications in crop protection, thereby providing help to establish novel RNAi-based crop protection strategies.

IPEV: identification of prokaryotic and eukaryotic virus-derived sequences in virome using deep learning.

BACKGROUND: The virome obtained through virus-like particle enrichment contains a mixture of prokaryotic and eukaryotic virus-derived fragments. Accurate identification and classification of these elements are crucial to understanding their roles and functions in microbial communities. However, the rapid mutation rates of viral genomes pose challenges in developing high-performance tools for classification, potentially limiting downstream analyses. FINDINGS: We present IPEV, a novel method to distinguish prokaryotic and eukaryotic viruses in viromes, with a 2-dimensional convolutional neural network combining trinucleotide pair relative distance and frequency. Cross-validation assessments of IPEV demonstrate its state-of-the-art precision, significantly improving the F1-score by approximately 22% on an independent test set compared to existing methods when query viruses share less than 30% sequence similarity with known viruses. Furthermore, IPEV outperforms other methods in accuracy on marine and gut virome samples based on annotations by sequence alignments. IPEV reduces runtime by at most 1,225 times compared to existing methods under the same computing configuration. We also utilized IPEV to analyze longitudinal samples and found that the gut virome exhibits a higher degree of temporal stability than previously observed in persistent personal viromes, providing novel insights into the resilience of the gut virome in individuals. CONCLUSIONS: IPEV is a high-performance, user-friendly tool that assists biologists in identifying and classifying prokaryotic and eukaryotic viruses within viromes. The tool is available at https://github.com/basehc/IPEV.

Unraveling the impact of climatic warming and wetting on eukaryotic microbial diversity and assembly mechanisms: A 10-year case study in Lake Bosten, NW China.

Over the last six decades, northwest China has undergone a significant climatic shift from "warm-dry" to "warm-wet", profoundly impacting the structures and functions of lake ecosystem across the region. However, the influences of this climatic transition on the diversity patterns, co-occurrence network, and assembly processes of eukaryotic microbial communities in lake ecosystem, along with the underlying mechanisms, remain largely unexplored. To bridge this knowledge gap, our study focused on Lake Bosten, the largest inland freshwater body in China, conducting a comprehensive analysis. Firstly, we examined the dynamics of key water quality parameters in the lake based on long-term monitoring data (1992-2022). Subsequently, we collected 93 water samples spanning two distinctive periods: low water level (WL) and high total dissolved solids (TDS) (PerWLTDS; 2010-2011; attributed to "warm-dry" climate), and high WL and low TDS (PerTDSWL; 2021-2022; associated with "warm-wet" climate). Eukaryotic microorganisms were further investigated using 18S rRNA gene sequencing and various statistical methods. Our findings revealed that climatic warming and wetting significantly increased eukaryotic microbial α-diversity (all Wilcox. test: P<0.05), while simultaneously reducing ß-diversity (all Wilcox. test: P<0.001) and network complexity. Through the two sampling periods, assembly mechanisms of eukaryotic microorganisms were predominantly influenced by dispersal limitation (DL) and drift (DR) within stochastic processes, alongside homogeneous selection (HoS) within deterministic processes. WL played a mediating role in eukaryotic microbial DL and HoS processes in the PerTDSWL, whereas water quality and α-diversity influenced the DL process in the PerWLTDS. Collectively, these results underscore the direct and indirect impacts of "warm-wet" conditions on the eukaryotic microorganisms within Lake Bosten. This study provides valuable insights into the evolutionary dynamics of lake ecosystems under such climatic conditions and aids in predicting the ecological ramifications of global climatic changes.

Planktonic eukaryotes in the Chesapeake Bay: contrasting responses of abundant and rare taxa to estuarine gradients.

Phytoplankton are important drivers of aquatic ecosystem function and environmental health. Their community compositions and distributions are directly impacted by environmental processes and human activities, including in the largest estuary in North America, the Chesapeake Bay. It is crucial to uncover how planktonic eukaryotes play fundamental roles as primary producers and trophic links and sustain estuarine ecosystems. In this study, we investigated the detailed community structure and spatiotemporal variations of planktonic eukaryotes in the Chesapeake Bay across space and time for three consecutive years. A clear seasonal and spatial shift of total, abundant, and rare planktonic eukaryotes was evident, and the pattern recurred interannually. Multiple harmful algal species have been identified in the Bay with varied distribution patterns, such as Karlodinium, Heterosigma akashiwo, Protoperidinium sp., etc. Compared to abundant taxa, rare subcommunities were more sensitive to environmental disturbance in terms of richness, diversity, and distribution. The combined effects of temporal variation (13.3%), nutrient availability (10.0%), and spatial gradients (8.8%) structured the distribution of eukaryotic microbial communities in the Bay. Similar spatiotemporal patterns between planktonic prokaryotes and eukaryotes suggest common mechanisms of adjustment, replacement, and species interaction for planktonic microbiomes under strong estuarine gradients. To our best knowledge, this work represents the first systematic study on planktonic eukaryotes in the Bay. A comprehensive view of the distribution of planktonic microbiomes and their interactions with environmental processes is critical in understanding the underlying microbial mechanisms involved in maintaining the stability, function, and environmental health of estuarine ecosystems. IMPORTANCE: Deep sequencing analysis of planktonic eukaryotes in the Chesapeake Bay reveals high community diversity with many newly recognized phytoplankton taxa. The Chesapeake Bay planktonic eukaryotes show distinct seasonal and spatial variability, with recurring annual patterns of total, abundant, and rare groups. Rare taxa mainly contribute to eukaryotic diversity compared to abundant groups, and they are more sensitive to spatiotemporal variations and environmental filtering. Temporal variations, nutrient availability, and spatial gradients significantly affect the distribution of eukaryotic microbial communities. Similar spatiotemporal patterns in prokaryotes and eukaryotes suggest common mechanisms of adjustment, substitution, and species interactions in planktonic microbiomes under strong estuarine gradients. Interannually recurring patterns demonstrate that diverse eukaryotic taxa have well adapted to the estuarine environment with a long residence time. Further investigations of how human activities impact estuarine planktonic eukaryotes are critical in understanding their essential ecosystem roles and in maintaining environmental safety and public health.

How good are global DNA-based environmental surveys for detecting all protist diversity? Arcellinida as an example of biased representation.

Metabarcoding approaches targeting microeukaryotes have deeply changed our vision of protist environmental diversity. The public repository EukBank consists of 18S v4 metabarcodes from 12,672 samples worldwide. To estimate how far this database provides a reasonable overview of all eukaryotic diversity, we used Arcellinida (lobose testate amoebae) as a case study. We hypothesised that (1) this approach would allow the discovery of unexpected diversity, but also that (2) some groups would be underrepresented because of primer/sequencing biases. Most of the Arcellinida sequences appeared in freshwater and soil, but their abundance and diversity appeared underrepresented. Moreover, 84% of ASVs belonged to the suborder Phryganellina, a supposedly species-poor clade, whereas the best-documented suborder (Glutinoconcha, 600 described species) was only marginally represented. We explored some possible causes of these biases. Mismatches in the primer-binding site seem to play a minor role. Excessive length of the target region could explain some of these biases, but not all. There must be some other unknown factors involved. Altogether, while metabarcoding based on ribosomal genes remains a good first approach to document microbial eukaryotic clades, alternative approaches based on other genes or sequencing techniques must be considered for an unbiased picture of the diversity of some groups.

How mitochondrial cristae illuminate the important role of oxygen during eukaryogenesis.

Inner membranes of mitochondria are extensively folded, forming cristae. The observed overall correlation between efficient eukaryotic ATP generation and the area of internal mitochondrial inner membranes both in unicellular organisms and metazoan tissues seems to explain why they evolved. However, the crucial use of molecular oxygen (O2) as final acceptor of the electron transport chain is still not sufficiently appreciated. O2 was an essential prerequisite for cristae development during early eukaryogenesis and could be the factor allowing cristae retention upon loss of mitochondrial ATP generation. Here I analyze illuminating bacterial and unicellular eukaryotic examples. I also discuss formative influences of intracellular O2 consumption on the evolution of the last eukaryotic common ancestor (LECA). These considerations bring about an explanation for the many genes coming from other organisms than the archaeon and bacterium merging at the start of eukaryogenesis.

Compatibility of termination mechanisms in bacterial transcription with inference on eukaryotic models.

Transcription termination has evolved to proceed through diverse mechanisms. For several classes of terminators, multiple models have been debatably proposed. Recent single-molecule studies on bacterial terminators have resolved several long-standing controversies. First, termination mode or outcome is twofold rather than single. RNA is released alone before DNA or together with DNA from RNA polymerase (RNAP), i.e. with RNA release for termination, RNAP retains on or dissociates off DNA, respectively. The concomitant release, described in textbooks, results in one-step decomposition of transcription complexes, and this 'decomposing termination' prevails at ρ factor-dependent terminators. Contrastingly, the sequential release was recently discovered abundantly from RNA hairpin-dependent intrinsic terminations. RNA-only release allows RNAP to diffuse on DNA in both directions and recycle for reinitiation. This 'recycling termination' enables one-dimensional reinitiation, which would be more expeditious than three-dimensional reinitiation by RNAP dissociated at decomposing termination. Second, while both recycling and decomposing terminations occur at a hairpin-dependent terminator, four termination mechanisms compatibly operate at a ρ-dependent terminator with ρ in alternative modes and even intrinsically without ρ. RNA-bound catch-up ρ mediates recycling termination first and decomposing termination later, while RNAP-prebound stand-by ρ invokes only decomposing termination slowly. Without ρ, decomposing termination occurs slightly and sluggishly. These four mechanisms operate on distinct timescales, providing orderly fail-safes. The stand-by mechanism is benefited by terminational pause prolongation and modulated by accompanying riboswitches more greatly than the catch-up mechanisms. Conclusively, any mechanism alone is insufficient to perfect termination, and multiple mechanisms operate compatibly to achieve maximum possible efficiency under separate controls.

Ribotin: automated assembly and phasing of rDNA morphs.

MOTIVATION: The ribosomal DNA (rDNA) arrays are highly repetitive and homogenous regions which exist in all life. Due to their repetitiveness, current assembly methods do not fully assemble the rDNA arrays in humans and many other eukaryotes, and so variation within the rDNA arrays cannot be effectively studied. RESULTS: Here, we present the tool ribotin to assemble full length rDNA copies, or morphs. Ribotin uses a combination of highly accurate long reads and extremely long nanopore reads to resolve the variation between rDNA morphs. We show that ribotin successfully recovers the most abundant morphs in human and nonhuman genomes. We also find that genome wide consensus sequences of the rDNA arrays frequently produce a mosaic sequence that does not exist in the genome. AVAILABILITY AND IMPLEMENTATION: Ribotin is available on https://github.com/maickrau/ribotin and as a package on bioconda.

Poly(A) tale: From A to A; RNA polyadenylation in prokaryotes and eukaryotes.

Most eukaryotic mRNAs and different non-coding RNAs undergo a form of 3' end processing known as polyadenylation. Polyadenylation machinery is present in almost all organisms except few species. In bacteria, the machinery has evolved from PNPase, which adds heteropolymeric tails, to a poly(A)-specific polymerase. Differently, a complex machinery for accurate polyadenylation and several non-canonical poly(A) polymerases are developed in eukaryotes. The role of poly(A) tail has also evolved from serving as a degradative signal to a stabilizing modification that also regulates translation. In this review, we discuss poly(A) tail emergence in prokaryotes and its development into a stable, yet dynamic feature at the 3' end of mRNAs in eukaryotes. We also describe how appearance of novel poly(A) polymerases gives cells flexibility to shape poly(A) tail. We explain how poly(A) tail dynamics help regulate cognate RNA metabolism in a context-dependent manner, such as during oocyte maturation. Finally, we describe specific mRNAs in metazoans that bear stem-loops instead of poly(A) tails. We conclude with how recent discoveries about poly(A) tail can be applied to mRNA technology. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Processing > 3' End Processing RNA Turnover and Surveillance > Regulation of RNA Stability.

Cip1, a CDK regulator, determines heterothallic mating or homothallic selfing in a protist.

Mating type (sex) plays a crucial role in regulating sexual reproduction in most extant eukaryotes. One of the functions of mating types is ensuring self-incompatibility to some extent, thereby promoting genetic diversity. However, heterothallic mating is not always the best mating strategy. For example, in low-density populations or specific environments, such as parasitic ones, species may need to increase the ratio of potential mating partners. Consequently, many species allow homothallic selfing (i.e., self-fertility or intraclonal mating). Throughout the extensive evolutionary history of species, changes in environmental conditions have influenced mating strategies back and forth. However, the mechanisms through which mating-type recognition regulates sexual reproduction and the dynamics of mating strategy throughout evolution remain poorly understood. In this study, we show that the Cip1 protein is responsible for coupling sexual reproduction initiation to mating-type recognition in the protozoal eukaryote Tetrahymena thermophila. Deletion of the Cip1 protein leads to the loss of the selfing-avoidance function of mating-type recognition, resulting in selfing without mating-type recognition. Further experiments revealed that Cip1 is a regulatory subunit of the Cdk19-Cyc9 complex, which controls the initiation of sexual reproduction. These results reveal a mechanism that regulates the choice between mating and selfing. This mechanism also contributes to the debate about the ancestral state of sexual reproduction.

Phylogenomic position of eupelagonemids, abundant, and diverse deep-ocean heterotrophs.

Eupelagonemids, formerly known as Deep Sea Pelagic Diplonemids I (DSPD I), are among the most abundant and diverse heterotrophic protists in the deep ocean, but little else is known about their ecology, evolution, or biology in general. Originally recognized solely as a large clade of environmental ribosomal subunit RNA gene sequences (SSU rRNA), branching with a smaller sister group DSPD II, they were postulated to be diplonemids, a poorly studied branch of Euglenozoa. Although new diplonemids have been cultivated and studied in depth in recent years, the lack of cultured eupelagonemids has limited data to a handful of light micrographs, partial SSU rRNA gene sequences, a small number of genes from single amplified genomes, and only a single formal described species, Eupelagonema oceanica. To determine exactly where this clade goes in the tree of eukaryotes and begin to address the overall absence of biological information about this apparently ecologically important group, we conducted single-cell transcriptomics from two eupelagonemid cells. A SSU rRNA gene phylogeny shows that these two cells represent distinct subclades within eupelagonemids, each different from E. oceanica. Phylogenomic analysis based on a 125-gene matrix contrasts with the findings based on ecological survey data and shows eupelagonemids branch sister to the diplonemid subgroup Hemistasiidae.

Meiosis: The silk moth and the elephant.

In most eukaryotes, balanced chromosome segregation at meiosis requires crossovers, but female Bombyx mori lack these structures. Instead, the synaptonemal complex is repurposed to compensate for this absence of crossovers, a remarkable example of exaptation.

Dynamic DNA N 6-adenine methylation (6mA) governs the encystment process, showcased in the unicellular eukaryote Pseudocohnilembus persalinus.

The formation of resting cysts commonly found in unicellular eukaryotes is a complex and highly regulated survival strategy against environmental stress that involves drastic physiological and biochemical changes. Although most studies have focused on the morphology and structure of cysts, little is known about the molecular mechanisms that control this process. Recent studies indicate that DNA N 6-adenine methylation (6mA) could be dynamically changing in response to external stimuli; however, its potential role in the regulation of cyst formation remains unknown. We used the ciliate Pseudocohnilembus persalinus, which can be easily induced to form cysts to investigate the dynamic pattern of 6mA in trophonts and cysts. Single-molecule real-time (SMRT) sequencing reveals high levels of 6mA in trophonts that decrease in cysts, along with a conversion of symmetric 6mA to asymmetric 6mA. Further analysis shows that 6mA, a mark of active transcription, is involved in altering the expression of encystment-related genes through changes in 6mA levels and 6mA symmetric-to-asymmetric conversion. Most importantly, we show that reducing 6mA levels by knocking down the DNA 6mA methyltransferase PpAMT1 accelerates cyst formation. Taken together, we characterize the genome-wide 6mA landscape in P. persalinus and provide insights into the role of 6mA in gene regulation under environmental stress in eukaryotes. We propose that 6mA acts as a mark of active transcription to regulate the encystment process along with symmetric-to-asymmetric conversion, providing important information for understanding the molecular response to environmental cues from the perspective of 6mA modification.

Modeling alternative translation initiation sites in plants reveals evolutionarily conserved cis-regulatory codes in eukaryotes.

mRNA translation relies on identifying translation initiation sites (TISs) in mRNAs. Alternative TISs are prevalent across plant transcriptomes, but the mechanisms for their recognition are unclear. Using ribosome profiling and machine learning, we developed models for predicting alternative TISs in the tomato (Solanum lycopersicum). Distinct feature sets were predictive of AUG and nonAUG TISs in 5' untranslated regions and coding sequences, including a novel CU-rich sequence that promoted plant TIS activity, a translational enhancer found across dicots and monocots, and humans and viruses. Our results elucidate the mechanistic and evolutionary basis of TIS recognition, whereby cis-regulatory RNA signatures affect start site selection. The TIS prediction model provides global estimates of TISs to discover neglected protein-coding genes across plant genomes. The prevalence of cis-regulatory signatures across plant species, humans, and viruses suggests their broad and critical roles in reprogramming the translational landscape.

When Paul Berg meets Donald Crothers: an achiral connection through protein biosynthesis.

Outliers in scientific observations are often ignored and mostly remain unreported. However, presenting them is always beneficial since they could reflect the actual anomalies that might open new avenues. Here, we describe two examples of the above that came out of the laboratories of two of the pioneers of nucleic acid research in the area of protein biosynthesis, Paul Berg and Donald Crothers. Their work on the identification of D-aminoacyl-tRNA deacylase (DTD) and 'Discriminator hypothesis', respectively, were hugely ahead of their time and were partly against the general paradigm at that time. In both of the above works, the smallest and the only achiral amino acid turned out to be an outlier as DTD can act weakly on glycine charged tRNAs with a unique discriminator base of 'Uracil'. This peculiar nature of glycine remained an enigma for nearly half a century. With a load of available information on the subject by the turn of the century, our work on 'chiral proofreading' mechanisms during protein biosynthesis serendipitously led us to revisit these findings. Here, we describe how we uncovered an unexpected connection between them that has implications for evolution of different eukaryotic life forms.