Forty Years of Inferential Methods in the Journals of the Society for Molecular Biology and Evolution.

We are launching a series to celebrate the 40th anniversary of the first issue of Molecular Biology and Evolution. In 2024, we will publish virtual issues containing selected papers published in the Society for Molecular Biology and Evolution journals, Molecular Biology and Evolution and Genome Biology and Evolution. Each virtual issue will be accompanied by a perspective that highlights the historic and contemporary contributions of our journals to a specific topic in molecular evolution. This perspective, the first in the series, presents an account of the broad array of methods that have been published in the Society for Molecular Biology and Evolution journals, including methods to infer phylogenies, to test hypotheses in a phylogenetic framework, and to infer population genetic processes. We also mention many of the software implementations that make methods tractable for empiricists. In short, the Society for Molecular Biology and Evolution community has much to celebrate after four decades of publishing high-quality science including numerous important inferential methods.

Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events.

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.

Foraminifera as a model of the extensive variability in genome dynamics among eukaryotes.

Knowledge of eukaryotic life cycles and associated genome dynamics stems largely from research on animals, plants, and a small number of "model" (i.e., easily cultivable) lineages. This skewed sampling results in an underappreciation of the variability among the many microeukaryotic lineages, which represent the bulk of eukaryotic biodiversity. The range of complex nuclear transformations that exists within lineages of microbial eukaryotes challenges the textbook understanding of genome and nuclear cycles. Here, we look in-depth at Foraminifera, an ancient (∼600 million-year-old) lineage widely studied as proxies in paleoceanography and environmental biomonitoring. We demonstrate that Foraminifera challenge the "rules" of life cycles developed largely from studies of plants and animals. To this end, we synthesize data on foraminiferal life cycles, focusing on extensive endoreplication within individuals (i.e., single cells), the unusual nuclear process called Zerfall, and the separation of germline and somatic function into distinct nuclei (i.e., heterokaryosis). These processes highlight complexities within lineages and expand our understanding of the dynamics of eukaryotic genomes.

Taxon-rich transcriptomics supports higher-level phylogeny and major evolutionary trends in Foraminifera.

Foraminifera, classified in the supergroup Rhizaria, are a common and highly diverse group of mainly marine protists. Despite their evolutionary and ecological importance, only limited genomic data (one partial genome and nine transcriptomic datasets) have been published for this group. Foraminiferal molecular phylogeny is largely based on 18S rRNA gene sequence analysis. However, due to highly variable evolutionary rates of substitution in ribosomal genes plus the existence of intragenomic variation at this locus, the relationships between and within foraminiferal classes remain uncertain. We analyze transcriptomic data from 28 species, adding 19 new species to the previously published dataset, including members of the strongly under-represented class Monothalamea. A phylogenomic reconstruction of Rhizaria, rooted with alveolates and stramenopiles, based on 199 genes and 68 species supports the monophyly of Foraminifera and their sister relationship to Polycystinea. The phylogenomic tree of Foraminifera is very similar to the 18S rRNA tree, with the paraphyletic single-chambered monothalamids giving rise to the multi-chambered Tubothalamea and Globothalamea. Within the Monothalamea, our analyses confirm the monophyly of the giant, deep-sea xenophyophores that branch within clade C and indicate the basal position of monothalamous clades D and E. The multi-chambered Globothalamea are monophyletic and comprise the paraphyletic Textulariida and monophyletic Rotaliida. Our phylogenomic analyses support major evolutionary trends of Foraminifera revealed by ribosomal phylogenies and reinforce their current higher-level classification.

Macronuclear development in ciliates, with a focus on nuclear architecture.

Ciliates are defined by the presence of dimorphic nuclei as they have both a somatic macronucleus and germline micronucleus within each individual cell. The size and structure of both germline micronuclei and somatic macronuclei vary tremendously among ciliates. Except just after conjugation (i.e. the nuclear exchange in their life cycle), the germline micronucleus is transcriptionally inactive and contains canonical chromosomes that will be inherited between generations. In contrast, the transcriptionally active macronucleus contains chromosomes that vary in size in different classes of ciliates, with some lineages having extensively fragmented gene-sized somatic chromosomes while others contain longer multigene chromosomes. Here, we describe the variation in somatic macronuclear architecture in lineages sampled across the ciliate tree of life, specifically focusing on lineages with extensively fragmented chromosomes (e.g. the classes Phyllopharyngea and Spirotrichea). Further, we synthesize information from the literature on the development of ciliate macronuclei, focusing on changes in nuclear architecture throughout life cycles. These data highlight the tremendous diversity among ciliate nuclear cycles, extend our understanding of patterns of genome evolution, and provide insight into different germline and somatic nuclear features (e.g. nuclear structure and development) among eukaryotes.

Epigenetic influences of mobile genetic elements on ciliate genome architecture and evolution.

Mobile genetic elements (MGEs) are transient genetic material that can move either within a single organism's genome or between individuals or species. While historically considered "junk" DNA (i.e., deleterious or at best neutral), more recent studies reveal the potential adaptive advantages MGEs provide in lineages across the tree of life. Ciliates, a group of single-celled microbial eukaryotes characterized by nuclear dimorphism, exemplify how epigenetic influences from MGEs shape genome architecture and patterns of molecular evolution. Ciliate nuclear dimorphism may have evolved as a response to transposon invasion and ciliates have since co-opted transposons to carry out programmed DNA deletion. Another example of the effect of MGEs is in providing mechanisms for lateral gene transfer (LGT) from bacteria, which introduces genetic diversity and, in several cases, may drive ecological specialization in ciliates. As a third example, the integration of viral DNA, likely through transduction, provides new genetic materials and can change the way host cells defend themselves against other viral pathogens. We argue that the acquisition of MGEs through non-Mendelian patterns of inheritance, coupled with their effects on ciliate genome architecture and persistence throughout evolutionary history, exemplify how the transmission of mobile elements should be considered a mechanism of transgenerational epigenetic inheritance.

Newly designed foraminifera primers identify habitat-specific lineages through metabarcoding analyses.

Foraminifera include diverse shell-building lineages found in a wide array of aquatic habitats from the deep-sea to intertidal zones to brackish and freshwater ecosystems. Recent estimates of morphological and molecular foraminifera diversity have increased the knowledge of foraminiferal diversity, which is critical as these lineages are used as bioindicators of past and present environmental perturbation. However, a comparative analysis of foraminiferal biodiversity between their major habitats (freshwater, brackish, intertidal, and marine) is underexplored, particularly using molecular tools. Here, we present a metabarcoding survey of foraminiferal diversity across different ecosystems using newly designed foraminifera-specific primers that target the hypervariable regions of the foraminifera SSU-rRNA gene (~250-300 bp long). We tested these primer sets on four foraminifera species and then across several environments: the intertidal zone, coastal ecosystems, and freshwater vernal pools. We retrieved 655 operational taxonomic units (OTUs); the majority of which are undetermined taxa that have no closely matching sequences in the reference database. Furthermore, we identified 163 OTUs with distinct habitat preferences. Most of the observed OTUs belonged to lineages of single-chambered foraminifera, including poorly explored freshwater foraminifera which encompass a clade of Reticulomyxa-like forms. Our pilot study provides the community with an additional set of newly designed and taxon-specific primers to elucidate foraminiferal diversity across different habitats.

Opinion: Genetic Conflict With Mobile Elements Drives Eukaryotic Genome Evolution, and Perhaps Also Eukaryogenesis.

Through analyses of diverse microeukaryotes, we have previously argued that eukaryotic genomes are dynamic systems that rely on epigenetic mechanisms to distinguish germline (i.e., DNA to be inherited) from soma (i.e., DNA that undergoes polyploidization, genome rearrangement, etc.), even in the context of a single nucleus. Here, we extend these arguments by including two well-documented observations: (1) eukaryotic genomes interact frequently with mobile genetic elements (MGEs) like viruses and transposable elements (TEs), creating genetic conflict, and (2) epigenetic mechanisms regulate MGEs. Synthesis of these ideas leads to the hypothesis that genetic conflict with MGEs contributed to the evolution of a dynamic eukaryotic genome in the last eukaryotic common ancestor (LECA), and may have contributed to eukaryogenesis (i.e., may have been a driver in the evolution of FECA, the first eukaryotic common ancestor). Sex (i.e., meiosis) may have evolved within the context of the development of germline-soma distinctions in LECA, as this process resets the germline genome by regulating/eliminating somatic (i.e., polyploid, rearranged) genetic material. Our synthesis of these ideas expands on hypotheses of the origin of eukaryotes by integrating the roles of MGEs and epigenetics.

High Diversity of Testate Amoebae (Amoebozoa, Arcellinida) Detected by HTS Analyses in a New England Fen using Newly Designed Taxon-specific Primers.

Testate (shell-building) amoebae, such as the Arcellinida (Amoebozoa), are useful bioindicators for climate change. Though past work has relied on morphological analyses to characterize Arcellinida diversity, genetic analyses revealed the presence of multiple cryptic species underlying morphospecies. Here, we design and deploy Arcellinida-specific primers for the SSU-rDNA gene to assess the community composition on the molecular level in a pilot study of two samplings from a New England fen: (1) 36-cm horizontal transects and vertical cores; and (2) 26-m horizontal transects fractioned into four size classes (2-10, 10-35, 35-100, and 100-300 µm). Analyses of these data show the following: (1) a considerable genetic diversity within Arcellinida, much of which comes from morphospecies lacking sequences on GenBank; (2) communities characterized by DNA (i.e. active + quiescent) are distinct from those characterized by RNA (i.e. active, indicator of biomass); (3) active communities on the surface tend to be more similar to one another than to core communities, despite considerable heterogeneity; and (4) analyses of communities fractioned by size find some lineages (OTUs) that are abundant in disjunct size categories, suggesting the possibility of life-history stages. Together, these data demonstrate the potential of these primers to elucidate the diversity of Arcellinida communities in diverse habitats.

De novo Sequencing, Assembly, and Annotation of the Transcriptome for the Free-Living Testate Amoeba Arcella intermedia.

Arcella, a diverse understudied genus of testate amoebae is a member of Tubulinea in Amoebozoa group. Transcriptomes are a powerful tool for characterization of these organisms as they are an efficient way of characterizing the protein-coding potential of the genome. In this work, we employed both single-cell and clonal populations transcriptomics to create a reference transcriptome for Arcella. We compared our results with annotations of Dictyostelium discoideum, a model Amoebozoan. We assembled a pool of 38 Arcella intermedia transcriptomes, which after filtering are composed of a total of 14,712 translated proteins. There are GO categories enriched in Arcella including mainly intracellular signal transduction pathways; we also used KEGG to annotate 11,546 contigs, which also have similar distribution to Dictyostelium. A large portion of data is still impossible to assign to a gene family, probably due to a combination of lineage-specific genes, incomplete sequences in the transcriptome and rapidly evolved genes. Some absences in pathways could also be related to low expression of these genes. We provide a reference database for Arcella, and we highlight the emergence of the need for further gene discovery in Arcella.