Search

Myxozoans are a monophyletic taxon of approximately 2,400 described species of parasites from the phylum Cnidaria. The recent focus on their negative impacts on fisheries, on their evolution from free-living ancestors, and on their emergence into new fish host populations has stressed the critical need for genomic resources for this parasitic group. Here, we describe the genome assembly and annotation of Myxobolus rasmusseni, an emerging parasite of fathead minnows in Alberta, Canada. The assembly is 174.6 Mb in size, 68% of which is made up of repetitive elements, making it one of the most repetitive animal genomes sequenced to date. Through comparisons to other myxozoans, we show that widespread gene loss, a known phenomenon of this group of parasites, is consistent with closely related species. Additionally, we assembled the M. rasmusseni mitochondrial genome, which is nearly twice the size of the typical animal mitochondrial genome yet contains only five of the canonical mitochondrial protein-coding genes and open reading frames not found in other myxozoans. These results add to our understanding of the gene- and genome-level diversity observed in myxozoans., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Zoantharia is an order among the Hexacorallia (Anthozoa: Cnidaria), and includes at least 300 species. Previously reported genomes from scleractinian corals and actiniarian sea anemones have illuminated part of the hexacorallian diversification. However, little is known about zoantharian genomes and the early evolution of hexacorals. To explore genome evolution in this group of hexacorals, here, we report de novo genome assemblies of the zoantharians Palythoa mizigama (Pmiz) and Palythoa umbrosa (Pumb), both of which are members of the family Sphenopidae, and uniquely live in comparatively dark coral reef caves without symbiotic Symbiodiniaceae dinoflagellates. Draft genomes generated from ultra-low input PacBio sequencing totaled 373 and 319 Mbp for Pmiz and Pumb, respectively. Protein-coding genes were predicted in each genome, totaling 30,394 in Pmiz and 24,800 in Pumb, with each set having ∼93% BUSCO completeness. Comparative genomic analyses identified 3,036 conserved gene families, which were found in all analyzed hexacoral genomes. Some of the genes related to toxins, chitin degradation, and prostaglandin biosynthesis were expanded in these two Palythoa genomes and many of which aligned tandemly. Extensive gene family loss was not detected in the Palythoa lineage and five of ten putatively lost gene families likely had neuronal function, suggesting biased gene loss in Palythoa. In conclusion, our comparative analyses demonstrate evolutionary conservation of gene families in the Palythoa lineage from the common ancestor of hexacorals. Restricted loss of gene families may imply that lost neuronal functions were effective for environmental adaptation in these two Palythoa species., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The subfamily GH13_1 of alpha-amylases is typical of Fungi, but it is also found in some unicellular eukaryotes (e.g., Amoebozoa, choanoflagellates) and non-bilaterian Metazoa. Since a previous study in 2007, GH13_1 amylases were considered ancestral to the Unikonts, including animals, except Bilateria, such that it was thought to have been lost in the ancestor of this clade. The only alpha-amylases known to be present in Bilateria so far belong to the GH13_15 and 24 subfamilies (commonly called bilaterian alpha-amylases) and were likely acquired by horizontal transfer from a proteobacterium. The taxonomic scope of Eukaryota genomes in databases has been greatly increased ever since 2007. We have surveyed GH13_1 sequences in recent data from ca. 1600 bilaterian species, 60 non-bilaterian animals and also in unicellular eukaryotes. As expected, we found a number of those sequences in non-bilaterians: Anthozoa (Cnidaria) and in sponges, confirming the previous observations, but none in jellyfishes and in Ctenophora. Our main and unexpected finding is that such fungal (also called Dictyo-type) amylases were also consistently retrieved in several bilaterian phyla: hemichordates (deuterostomes), brachiopods and related phyla, some molluscs and some annelids (protostomes). We discuss evolutionary hypotheses possibly explaining the scattered distribution of GH13_1 across bilaterians, namely, the retention of the ancestral gene in those phyla only and/or horizontal transfers from non-bilaterian donors.

Background and Aims: Biological aspects of haustorial parasitism have significant effects on the configuration of the plastid genome. Approximately half the diversity of haustorial parasites belongs to the order Santalales, where a clearer picture of plastome evolution in relation to parasitism is starting to emerge. However, in previous studies of plastome evolution there is still a notable under-representation of members from non-parasitic and deep-branching hemiparasitic lineages, limiting evolutionary inference around the time of transition to a parasitic lifestyle. To expand taxon sampling relevant to this transition we therefore targeted three families of non-parasites (Erythropalaceae, Strombosiaceae, and Coulaceae), two families of root-feeding hemiparasites (Ximeniaceae and Olacaceae), and two families of uncertain parasitic status (Aptandraceae and Octoknemaceae). With data from these lineages we aimed to explore plastome evolution in relation to evolution of parasitism., Methods: From 29 new samples we sequenced and annotated plastomes and the nuclear ribosomal cistron. We examined phylogenetic patterns, plastome evolution, and patterns of relaxed or intensified selection in plastid genes. Available transcriptome data were analyzed to investigate potential transfer of infA to the nuclear genome., Results: Phylogenetic relationships indicate a single functional loss of all plastid ndh genes (ndhA-K) in a clade formed by confirmed parasites and Aptandraceae, and the loss coincides with major size and boundary shifts of the inverted repeat (IR) region. Depending on an autotrophic or heterotrophic lifestyle in Aptandraceae, plastome changes are either correlated with or predate evolution of parasitism. Phylogenetic patterns also indicate repeated loss of infA from the plastome, and based on presence of transcribed sequences with presequences corresponding to thylakoid luminal transit peptides, we infer that the genes were transferred to the nuclear genome., Conclusions: Except for the loss of the ndh complex, relatively few genes have been lost from the plastome in deep-branching root parasites in Santalales. Prior to loss of the ndh genes, they show signs of relaxed selection indicative of their dispensability. To firmly establish a potential correlation between ndh gene loss, plastome instability and evolution of parasitism, it is pertinent to refute or confirm a parasitic lifestyle all Santalales clades., (© The Author(s) 2024. Published by Oxford University Press on behalf of the Annals of Botany Company.)

The evolution of gene families is complex, involving gene-level evolutionary events such as gene duplication, horizontal gene transfer, and gene loss, and other processes such as incomplete lineage sorting (ILS). Because of this, topological differences often exist between gene trees and species trees. A number of models have been recently developed to explain these discrepancies, the most realistic of which attempts to consider both gene-level events and ILS. When unified in a single model, the interaction between ILS and gene-level events can cause polymorphism in gene copy number, which we refer to as copy number hemiplasy (CNH). In this paper, we extend the Wright-Fisher process to include duplications and losses over several species, and show that the probability of CNH for this process can be significant. We study how well two unified models-multilocus multispecies coalescent (MLMSC), which models CNH, and duplication, loss, and coalescence (DLCoal), which does not-approximate the Wright-Fisher process with duplication and loss. We then study the effect of CNH on gene family evolution by comparing MLMSC and DLCoal. We generate comparable gene trees under both models, showing significant differences in various summary statistics; most importantly, CNH reduces the number of gene copies greatly. If this is not taken into account, the traditional method of estimating duplication rates (by counting the number of gene copies) becomes inaccurate. The simulated gene trees are also used for species tree inference with the summary methods ASTRAL and ASTRAL-Pro, demonstrating that their accuracy, based on CNH-unaware simulations calibrated on real data, may have been overestimated., (© The Author(s) 2024. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For commercial re.use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site.for further information please contact journals.permissions@oup.com.)

The redbanded stink bug, Piezodorus guildinii (Westwood) (Hemiptera: Pentatomidae), is a significant soybean pest in the Americas, which inflicts more physical damage on soybean than other native stink bugs. Studies suggest that its heightened impact is attributed to the aggressive digestive properties of its saliva. Despite its agricultural importance, the factors driving its greater ability to degrade plant tissues have remained unexplored in a genomic evolutionary context. In this study, we hypothesized that lineage-specific gene family expansions have increased the copy number of digestive genes expressed in the salivary glands. To investigate this, we annotated a previously published genome assembly of the redbanded stink bug, performed a comparative genomic analysis on 11 hemipteran species, and reconstructed patterns of gene duplication, gain, and loss in the redbanded stink bug. We also performed RNA-seq on the redbanded stink bug's salivary tissues, along with the rest of the body without salivary glands. We identified hundreds of differentially expressed salivary genes, including a subset lost in other stink bug lineages, but retained and expressed in the redbanded stink bug's salivary glands. These genes were significantly enriched with protein families involved in proteolysis, potentially explaining the redbanded stink bug's heightened damage to soybeans. Contrary to our hypothesis, we found no support for an enrichment of duplicated digestive genes that are also differentially expressed in the salivary glands of the redbanded stink bug. Nonetheless, these results provide insight into the evolution of this important crop pest, establishing a link between its genomic history and its agriculturally important physiology., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Transglutaminases (TGMs) cross-link proteins by introducing covalent bonds between glutamine and lysine residues. These cross-links are essential for epithelial cornification which enables tetrapods to live on land. Here, we investigated which evolutionary adaptations of vertebrates were associated with specific changes in the family of TGM genes. We determined the catalog of TGMs in the main clades of vertebrates, performed a comprehensive phylogenetic analysis of TGMs, and localized the distribution of selected TGMs in tissues. Our data suggest that TGM1 is the phylogenetically oldest epithelial TGM, with orthologs being expressed in the cornified teeth of the lamprey, a basal vertebrate. Gene duplications led to the origin of TGM10 in stem vertebrates, the origin of TGM2 in jawed vertebrates, and an increasing number of epithelium-associated TGM genes in the lineage leading to terrestrial vertebrates. TGM9 is expressed in the epithelial egg tooth, and its evolutionary origin in stem amniotes coincided with the evolution of embryonic development in eggs that are surrounded by a protective shell. Conversely, viviparous mammals have lost both the epithelial egg tooth and TGM9. TGM3 and TGM6 evolved as regulators of cornification in hair follicles and underwent pseudogenization upon the evolutionary loss of hair in cetaceans. Taken together, this study reveals the gain and loss of vertebrate TGM genes in association with the evolution of cornified skin appendages and suggests an important role of TGM9 in the evolution of amniotes., Competing Interests: Conflict of interest: The authors declare no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Comparative approaches have revealed both divergent and convergent paths to achieving shared developmental outcomes. Thus, only through assembling multiple case studies can we understand biological principles. Yet, despite appreciating the conservation-or lack thereof-of developmental networks, the conservation of epigenetic mechanisms regulating these networks is poorly understood. The nematode Pristionchus pacificus has emerged as a model system of plasticity and epigenetic regulation as it exhibits a bacterivorous or omnivorous morph depending on its environment. Here, we determined the "epigenetic toolkit" available to P. pacificus as a resource for future functional work on plasticity, and as a comparison with Caenorhabditis elegans to investigate the conservation of epigenetic mechanisms. Broadly, we observed a similar cast of genes with putative epigenetic function between C. elegans and P. pacificus. However, we also found striking differences. Most notably, the histone methyltransferase complex PRC2 appears to be missing in P. pacificus. We described the deletion/pseudogenization of the PRC2 genes mes-2 and mes-6 and concluded that both were lost in the last common ancestor of P. pacificus and a related species P. arcanus. Interestingly, we observed the enzymatic product of PRC2 (H3K27me3) by mass spectrometry and immunofluorescence, suggesting that a currently unknown methyltransferase has been co-opted for heterochromatin silencing. Altogether, we have provided an inventory of epigenetic genes in P. pacificus to compare with C. elegans. This inventory will enable reverse-genetic experiments related to plasticity and has revealed the first loss of PRC2 in a multicellular organism., Competing Interests: Conflicts of interest: The author(s) declare no conflicts of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here, we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation-mediated by the trans-regulator Spt10-is ancient, likely emerging between 320 and 380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage, which coincided with the loss of 1 copy of its paralogous core histone genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora fast-evolving lineage., Competing Interests: Conflicts of interest JDB is a Founder and Director of CDI Labs, Inc., a Founder of and consultant to Neochromosome, Inc., a Founder, SAB member of, and consultant to ReOpen Diagnostics, LLC, and serves or served on the Scientific Advisory Board of Logomix, Inc., Modern Meadow, Inc., Rome Therapeutics, Inc., Sample6, Inc., Sangamo, Inc., Tessera Therapeutics, Inc., and the Wyss Institute. JLS is a scientific advisor for WittGen Biotechnologies. JLS is an advisor for ForensisGroup, Inc., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

It has been predicted that the highly degenerate mammalian Y chromosome will be lost eventually. Indeed, Y was lost in the Ryukyu spiny rat Tokudaia osimensis, but the fate of the formerly Y-linked genes is not completely known. We looked for all 12 ancestrally Y-linked genes in a draft T. osimensis genome sequence. Zfy1, Zfy2, Kdm5d, Eif2s3y, Usp9y, Uty, and Ddx3y are putatively functional and are now located on the X chromosome, whereas Rbmy, Uba1y, Ssty1, Ssty2, and Sry are missing or pseudogenized. Tissue expressions of the mouse orthologs of the retained genes are significantly broader/higher than those of the lost genes, suggesting that the destinies of the formerly Y-linked genes are related to their original expressions. Interestingly, patterns of gene retention/loss are significantly more similar than by chance across four rodent lineages where Y has been independently lost, indicating a level of certainty in the fate of Y-linked genes even when the chromosome is gone., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Amphibians are ideal for studying visual system evolution because their biphasic (aquatic and terrestrial) life history and ecological diversity expose them to a broad range of visual conditions. Here, we evaluate signatures of selection on visual opsin genes across Neotropical anurans and focus on three diurnal clades that are well-known for the concurrence of conspicuous colors and chemical defense (i.e., aposematism): poison frogs (Dendrobatidae), Harlequin toads (Bufonidae: Atelopus), and pumpkin toadlets (Brachycephalidae: Brachycephalus). We found evidence of positive selection on 44 amino acid sites in LWS, SWS1, SWS2, and RH1 opsin genes, of which one in LWS and two in RH1 have been previously identified as spectral tuning sites in other vertebrates. Given that anurans have mostly nocturnal habits, the patterns of selection revealed new sites that might be important in spectral tuning for frogs, potentially for adaptation to diurnal habits and for color-based intraspecific communication. Furthermore, we provide evidence that SWS2, normally expressed in rod cells in frogs and some salamanders, has likely been lost in the ancestor of Dendrobatidae, suggesting that under low-light levels, dendrobatids have inferior wavelength discrimination compared to other frogs. This loss might follow the origin of diurnal activity in dendrobatids and could have implications for their behavior. Our analyses show that assessments of opsin diversification in across taxa could expand our understanding of the role of sensory system evolution in ecological adaptation., Competing Interests: Conflict of interest: This publication is based in part on work by D.C.C. while serving at the National Science Foundation. Any opinion, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation or United States government., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Gene loss is a prevalent source of genetic variation in genome evolution. Calling loss events effectively and efficiently is a critical step for systematically characterizing their functional and phylogenetic profiles genome wide. Here, we developed a novel pipeline integrating orthologous inference and genome alignment. Interestingly, we identified 33 gene loss events that give rise to evolutionarily novel long noncoding RNAs (lncRNAs) that show distinct expression features and could be associated with various functions related to growth, development, immunity, and reproduction, suggesting loss relics as a potential source of functional lncRNAs in humans. Our data also demonstrated that the rates of protein gene loss are variable among different lineages with distinct functional biases., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The origin of microbial mercury methylation has long been a mystery. Here, we employed genome-resolved phylogenetic analyses to decipher the evolution of the mercury-methylating gene, hgcAB, constrain the ancestral origin of the hgc operon, and explain the distribution of hgc in Bacteria and Archaea. We infer the extent to which vertical inheritance and horizontal gene transfer have influenced the evolution of mercury methylators and hypothesize that evolution of this trait bestowed the ability to produce an antimicrobial compound (MeHg+) on a potentially resource-limited early Earth. We speculate that, in response, the evolution of MeHg+-detoxifying alkylmercury lyase (encoded by merB) reduced a selective advantage for mercury methylators and resulted in widespread loss of hgc in Bacteria and Archaea., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Change in gene family size has been shown to facilitate adaptation to different selective pressures. This includes gene duplication to increase dosage or diversification of enzymatic substrates and gene deletion due to relaxed selection. We recently found that the PON1 gene, an enzyme with arylesterase and lactonase activity, was lost repeatedly in different aquatic mammalian lineages, suggesting that the PON gene family is responsive to environmental change. We further investigated if these fluctuations in gene family size were restricted to mammals and approximately when this gene family was expanded within mammals. Using 112 metazoan protein models, we explored the evolutionary history of the PON family to characterize the dynamic evolution of this gene family. We found that there have been multiple, independent expansion events in tardigrades, cephalochordates, and echinoderms. In addition, there have been partial gene loss events in monotremes and sea cucumbers and what appears to be complete loss in arthropods, urochordates, platyhelminths, ctenophores, and placozoans. In addition, we show the mammalian expansion to three PON paralogs occurred in the ancestor of all mammals after the divergence of sauropsida but before the divergence of monotremes from therians. We also provide evidence of a novel PON expansion within the brushtail possum. In the face of repeated expansions and deletions in the context of changing environments, we suggest a range of selective pressures, including pathogen infection and mitigation of oxidative damage, are likely influencing the diversification of this dynamic gene family across metazoa., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Gene repertoire turnover is a characteristic of genome evolution. However, we lack well-replicated analyses of presence/absence patterns associated with different selection contexts. Here, we study ∼100 transcriptome assemblies across Macrostomum, a genus of simultaneously hermaphroditic flatworms exhibiting multiple convergent shifts in mating strategy and associated reproductive morphologies. Many species mate reciprocally, with partners donating and receiving sperm at the same time. Other species convergently evolved to mate by hypodermic injection of sperm into the partner. We find that for orthologous transcripts annotated as expressed in the body region containing the testes, sequences from hypodermically inseminating species diverge more rapidly from the model species, Macrostomum lignano, and have a lower probability of being observed in other species. For other annotation categories, simpler models with a constant rate of similarity decay with increasing genetic distance from M. lignano match the observed patterns well. Thus, faster rates of sequence evolution for hypodermically inseminating species in testis-region genes result in higher rates of homology detection failure, yielding a signal of rapid evolution in sequence presence/absence patterns. Our results highlight the utility of considering appropriate null models for unobserved genes, as well as associating patterns of gene presence/absence with replicated evolutionary events in a phylogenetic context., (© 2022 The Authors. Evolution published by Wiley Periodicals LLC on behalf of The Society for the Study of Evolution.)

Rhizobiales are well-known plant-root nitrogen-fixing symbionts, but the functions of insect-associated Rhizobiales are poorly understood. We obtained genomes of three strains associated with Acromyrmex leaf-cutting ants and show that, in spite of being extracellular gut symbionts, they lost all pathways for essential amino acid biosynthesis, making them fully dependent on their hosts. Comparison with 54 Rhizobiales genomes showed that all insect-associated Rhizobiales lost the ability to fix nitrogen and that the Acromyrmex symbionts had exceptionally also lost the urease genes. However, the Acromyrmex strains share biosynthesis pathways for riboflavin vitamin, queuosine and a wide range of antioxidant enzymes likely to be beneficial for the ant fungus-farming symbiosis. We infer that the Rhizobiales symbionts catabolize excess of fungus-garden-derived arginine to urea, supplementing complementary Mollicutes symbionts that turn arginine into ammonia and infer that these combined symbiont activities stabilize the fungus-farming mutualism. Similar to the Mollicutes symbionts, the Rhizobiales species have fully functional CRISPR/Cas and R-M phage defenses, suggesting that these symbionts are important enough for the ant hosts to have precluded the evolution of metabolically cheaper defenseless strains., (© The Author(s) 2022. Published by Oxford University Press on behalf of FEMS.)

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Chromatin remodelers play a fundamental role in the assembly of chromatin, regulation of transcription, and DNA repair. Biochemical and functional characterizations of the CHD family of chromatin remodelers from a variety of model organisms have shown that these remodelers participate in a wide range of activities. However, because the evolutionary history of CHD homologs is unclear, it is difficult to predict which of these activities are broadly conserved and which have evolved more recently in individual eukaryotic lineages. Here, we performed a comprehensive phylogenetic analysis of 8,042 CHD homologs from 1,894 species to create a model for the evolution of this family across eukaryotes with a particular focus on the timing of duplications that gave rise to the diverse copies observed in plants, animals, and fungi. Our analysis confirms that the three major subfamilies of CHD remodelers originated in the eukaryotic last common ancestor, and subsequent losses occurred independently in different lineages. Improved taxon sampling identified several subfamilies of CHD remodelers in plants that were absent or highly divergent in the model plant Arabidopsis thaliana. Whereas the timing of CHD subfamily expansions in vertebrates corresponds to whole genome duplication events, the mechanisms underlying CHD diversification in land plants appear more complicated. Analysis of protein domains reveals that CHD remodeler diversification has been accompanied by distinct transitions in domain architecture, contributing to the functional differences observed between these remodelers. This study demonstrates the importance of proper taxon sampling when studying ancient evolutionary events to prevent misinterpretation of subsequent lineage-specific changes and provides an evolutionary framework for functional and comparative analysis of this critical chromatin remodeler family across eukaryotes., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Fungal species of the Ceratocystidaceae grow on their host plants using a variety of different lifestyles, from saprophytic to highly pathogenic. Although many genomes of fungi in the Ceratocystidaceae are publicly available, it is not known how the genes that encode catechol dioxygenases (CDOs), enzymes involved in the degradation of phenolic plant defense compounds, differ among members of the Ceratocystidaceae. The aim of this study was therefore to identify and characterize the genes encoding CDOs in the genomes of Ceratocystidaceae representatives. We found that genes encoding CDOs are more abundant in pathogenic necrotrophic species of the Ceratocystidaceae and less abundant in saprophytic species. The loss of the CDO genes and the associated 3-oxoadipate catabolic pathway appears to have occurred in a lineage-specific manner. Taken together, this study revealed a positive association between CDO gene copy number and fungal lifestyle in Ceratocystidaceae representatives., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Species tree inference from gene family trees is becoming increasingly popular because it can account for discordance between the species tree and the corresponding gene family trees. In particular, methods that can account for multiple-copy gene families exhibit potential to leverage paralogy as informative signal. At present, there does not exist any widely adopted inference method for this purpose. Here, we present SpeciesRax, the first maximum likelihood method that can infer a rooted species tree from a set of gene family trees and can account for gene duplication, loss, and transfer events. By explicitly modeling events by which gene trees can depart from the species tree, SpeciesRax leverages the phylogenetic rooting signal in gene trees. SpeciesRax infers species tree branch lengths in units of expected substitutions per site and branch support values via paralogy-aware quartets extracted from the gene family trees. Using both empirical and simulated data sets we show that SpeciesRax is at least as accurate as the best competing methods while being one order of magnitude faster on large data sets at the same time. We used SpeciesRax to infer a biologically plausible rooted phylogeny of the vertebrates comprising 188 species from 31,612 gene families in 1 h using 40 cores. SpeciesRax is available under GNU GPL at https://github.com/BenoitMorel/GeneRax and on BioConda., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Most organisms possess time-keeping devices called circadian clocks. At the molecular level, circadian clocks consist of transcription-translation feedback loops (TTFLs). Although some components of the negative TTFL are conserved across the animals, important differences exist between typical models, such as mouse and the fruit fly. In Drosophila, the key components are PERIOD (PER) and TIMELESS (TIM-d) proteins, whereas the mammalian clock relies on PER and CRYPTOCHROME (CRY-m). Importantly, how the clock has maintained functionality during evolutionary transitions between different states remains elusive. Therefore, we systematically described the circadian clock gene setup in major bilaterian lineages and identified marked lineage-specific differences in their clock constitution. Then we performed a thorough functional analysis of the linden bug Pyrrhocoris apterus, an insect species comprising features characteristic of both the Drosophila and the mammalian clocks. Unexpectedly, the knockout of timeless-d, a gene essential for the clock ticking in Drosophila, did not compromise rhythmicity in P. apterus, it only accelerated its pace. Furthermore, silencing timeless-m, the ancestral timeless type ubiquitously present across animals, resulted in a mild gradual loss of rhythmicity, supporting its possible participation in the linden bug clock, which is consistent with timeless-m role suggested by research on mammalian models. The dispensability of timeless-d in P. apterus allows drawing a scenario in which the clock has remained functional at each step of transition from an ancestral state to the TIM-d-independent PER + CRY-m system operating in extant vertebrates, including humans., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Acidic chitinase (Chia) digests the chitin of insects in the omnivorous stomach and the chitinase activity in carnivorous Chia is significantly lower than that of the omnivorous enzyme. However, mechanistic and evolutionary insights into the functional changes in Chia remain unclear. Here we show that a noninsect-based diet has caused structural and functional changes in Chia during the course of evolution in Carnivora. By creating mouse-dog chimeric Chia proteins and modifying the amino acid sequences, we revealed that F214L and A216G substitutions led to the dog enzyme activation. In 31 Carnivora, Chia was present as a pseudogene with stop codons in the open reading frame (ORF) region. Importantly, the Chia proteins of skunk, meerkat, mongoose, and hyena, which are insect-eating species, showed high chitinolytic activity. The cat Chia pseudogene product was still inactive even after ORF restoration. However, the enzyme was activated by matching the number and position of Cys residues to an active form and by introducing five meerkat Chia residues. Mutations affecting the Chia conformation and activity after pseudogenization have accumulated in the common ancestor of Felidae due to functional constraints. Evolutionary analysis indicates that Chia genes are under relaxed selective constraint in species with noninsect-based diets except for Canidae. These results suggest that there are two types of inactivating processes in Carnivora and that dietary changes affect the structure and activity of Chia., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Molecular genetic data have recently been incorporated in attempts to reconstruct the ecology of the ancestral snake, though this has been limited by a paucity of data for one of the two main extant snake taxa, the highly fossorial Scolecophidia. Here we present and analyze vision genes from the first eye-transcriptomic and genome-wide data for Scolecophidia, for Anilios bicolor, and A. bituberculatus, respectively. We also present immunohistochemistry data for retinal anatomy and visual opsin-gene expression in Anilios. Analyzed in the context of 19 lepidosaurian genomes and 12 eye transcriptomes, the new genome-wide and transcriptomic data provide evidence for a much more reduced visual system in Anilios than in non-scolecophidian (=alethinophidian) snakes and in lizards. In Anilios, there is no evidence of the presence of 7 of the 12 genes associated with alethinophidian photopic (cone) phototransduction. This indicates extensive gene loss and many of these candidate gene losses occur also in highly fossorial mammals with reduced vision. Although recent phylogenetic studies have found evidence for scolecophidian paraphyly, the loss in Anilios of visual genes that are present in alethinophidians implies that the ancestral snake had a better-developed visual system than is known for any extant scolecophidian., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Colletotrichum species are known to engage in unique sexual behaviors that differ significantly from the mating strategies of other filamentous ascomycete species. For example, most ascomycete fungi require the expression of both the MAT1-1-1 and MAT1-2-1 genes to induce sexual reproduction. In contrast, all isolates of Colletotrichum harbor only the MAT1-2-1 gene and yet, are capable of recognizing suitable mating partners and producing sexual progeny. The molecular mechanisms contributing to mating types and behaviors in Colletotrichum are, however, unknown. A comparative genomics approach analyzing 35 genomes, representing 31 Colletotrichum species and two Verticillium species, was used to elucidate a putative molecular mechanism underlying the unique sexual behaviors observed in Colletotrichum species. The existence of only the MAT1-2 idiomorph was confirmed across all species included in this study. Comparisons of the loci harboring the two mating pheromones and their cognate receptors revealed interesting patterns of gene presence and absence. The results showed that these genes have been lost multiple, independent times over the evolutionary history of this genus. These losses indicate that the pheromone pathway no longer plays an active role in mating type determination, suggesting an undiscovered mechanism by which mating partner recognition is controlled in these species. This further suggests that there has been a redirection of the underlying genetic mechanisms that regulate sexual development in Colletotrichum species. This research thus provides a foundation from which further interrogation of this topic can take place., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Arylalkylamine N-acetyltransferase (AANAT) plays a crucial role in synchronizing internal biological functions to circadian and circannual changes. Generally speaking, only one copy of AANAT gene has been found in mammals, however, three independent duplications of this gene were detected in several cetartiodactyl lineages (i.e., Suidae, Hippopotamidae, and Pecora), which originated in the middle Eocene, a geological period characterized with the increased climate seasonality. Lineage-specific expansions of AANAT and the associated functional enhancement in these lineages strongly suggest an improvement in regulating photoperiodic response to adapt to seasonal climate changes. In contrast, independent inactivating mutations or deletions of the AANAT locus were identified in the four pineal-deficient clades (cetaceans, sirenians, xenarthrans, and pangolins). Loss of AANAT function in cetaceans and sirenians could disrupt the sleep-promoting effects of pineal melatonin, which might contribute to increasing wakefulness, adapting these clades to underwater sleep. The absence of AANAT and pineal glands in xenarthrans and pangolins may be associated with their body temperature maintenance. The present work demonstrates a far more complex and intriguing evolutionary pattern and functional diversity of mammalian AANAT genes than previously thought and provides further evidence for understanding AANAT evolution as driven by rhythmic adaptations in mammals., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Insights into the evolution of ancestral complexes and pathways are generally achieved through careful and time-intensive manual analysis often using phylogenetic profiles of the constituent proteins. This manual analysis limits the possibility of including more protein-complex components, repeating the analyses for updated genome sets or expanding the analyses to larger scales. Automated orthology inference should allow such large-scale analyses, but substantial differences between orthologous groups generated by different approaches are observed. We evaluate orthology methods for their ability to recapitulate a number of observations that have been made with regard to genome evolution in eukaryotes. Specifically, we investigate phylogenetic profile similarity (co-occurrence of complexes), the last eukaryotic common ancestor's gene content, pervasiveness of gene loss and the overlap with manually determined orthologous groups. Moreover, we compare the inferred orthologies to each other. We find that most orthology methods reconstruct a large last eukaryotic common ancestor, with substantial gene loss, and can predict interacting proteins reasonably well when applying phylogenetic co-occurrence. At the same time, derived orthologous groups show imperfect overlap with manually curated orthologous groups. There is no strong indication of which orthology method performs better than another on individual or all of these aspects. Counterintuitively, despite the orthology methods behaving similarly regarding large-scale evaluation, the obtained orthologous groups differ vastly from one another. Availability and implementation The data and code underlying this article are available in github and/or upon reasonable request to the corresponding author: https://github.com/ESDeutekom/ComparingOrthologies., (© The Author(s) 2020. Published by Oxford University Press.)

Avian genomes are small and lack some genes that are conserved in the genomes of most other vertebrates including nonavian sauropsids. One hypothesis stated that paralogs may provide biochemical or physiological compensation for certain gene losses; however, no functional evidence has been reported to date. By integrating evolutionary analysis, physiological genomics, and experimental gene interference, we clearly demonstrate functional compensation for gene loss. A large-scale phylogenetic analysis of over 1,400 SLC2 gene sequences identifies six new SLC2 genes from nonmammalian vertebrates and divides the SLC2 gene family into four classes. Vertebrates retain class III SLC2 genes but partially lack the more recent duplicates of classes I and II. Birds appear to have completely lost the SLC2A4 gene that encodes an important insulin-sensitive GLUT in mammals. We found strong evidence for positive selection, indicating that the N-termini of SLC2A4 and SLC2A12 have undergone diversifying selection in birds and mammals, and there is a significant correlation between SLC2A12 functionality and basal metabolic rates in endotherms. Physiological genomics have uncovered that SLC2A12 expression and allelic variants are associated with insulin sensitivity and blood glucose levels in wild birds. Functional tests have indicated that SLC2A12 abrogation causes hyperglycemia, insulin resistance, and high relative activity, thus increasing energy expenditures that resemble a diabetic phenotype. These analyses suggest that the SLC2A12 gene not only functionally compensates insulin response for SLC2A4 loss but also affects daily physical behavior and basal metabolic rate during bird evolution, highlighting that older genes retain a higher level of functional diversification., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Digitaria exilis, white fonio, is a minor but vital crop of West Africa that is valued for its resilience in hot, dry, and low-fertility environments and for the exceptional quality of its grain for human nutrition. Its success is hindered, however, by a low degree of plant breeding and improvement., Findings: We sequenced the fonio genome with long-read SMRT-cell technology, yielding a ∼761 Mb assembly in 3,329 contigs (N50, 1.73 Mb; L50, 126). The assembly approaches a high level of completion, with a BUSCO score of >99%. The fonio genome was found to be a tetraploid, with most of the genome retained as homoeologous duplications that differ overall by ∼4.3%, neglecting indels. The 2 genomes within fonio were found to have begun their independent divergence ∼3.1 million years ago. The repeat content (>49%) is fairly standard for a grass genome of this size, but the ratio of Gypsy to Copia long terminal repeat retrotransposons (∼6.7) was found to be exceptionally high. Several genes related to future improvement of the crop were identified including shattering, plant height, and grain size. Analysis of fonio population genetics, primarily in Mali, indicated that the crop has extensive genetic diversity that is largely partitioned across a north-south gradient coinciding with the Sahel and Sudan grassland domains., Conclusions: We provide a high-quality assembly, annotation, and diversity analysis for a vital African crop. The availability of this information should empower future research into further domestication and improvement of fonio., (© The Author(s) 2021. Published by Oxford University Press GigaScience.)

Dollo's law posits that evolutionary losses are irreversible, thereby narrowing the potential paths of evolutionary change. While phenotypic reversals to ancestral states have been observed, little is known about their underlying genetic causes. The genomes of budding yeasts have been shaped by extensive reductive evolution, such as reduced genome sizes and the losses of metabolic capabilities. However, the extent and mechanisms of trait reacquisition after gene loss in yeasts have not been thoroughly studied. Here, through phylogenomic analyses, we reconstructed the evolutionary history of the yeast galactose utilization pathway and observed widespread and repeated losses of the ability to utilize galactose, which occurred concurrently with the losses of GALactose (GAL) utilization genes. Unexpectedly, we detected multiple galactose-utilizing lineages that were deeply embedded within clades that underwent ancient losses of galactose utilization. We show that at least two, and possibly three, lineages reacquired the GAL pathway via yeast-to-yeast horizontal gene transfer. Our results show how trait reacquisition can occur tens of millions of years after an initial loss via horizontal gene transfer from distant relatives. These findings demonstrate that the losses of complex traits and even whole pathways are not always evolutionary dead-ends, highlighting how reversals to ancestral states can occur., (© The Author(s) 2020. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Plectus murrayi is one of the most common and locally abundant invertebrates of continental Antarctic ecosystems. Because it is readily cultured on artificial medium in the laboratory and highly tolerant to an extremely harsh environment, P. murrayi is emerging as a model organism for understanding the evolutionary origin and maintenance of adaptive responses to multiple environmental stressors, including freezing and desiccation. The de novo assembled genome of P. murrayi contains 225.741 million base pairs and a total of 14,689 predicted genes. Compared to Caenorhabditis elegans, the architectural components of P. murrayi are characterized by a lower number of protein-coding genes, fewer transposable elements, but more exons, than closely related taxa from less harsh environments. We compared the transcriptomes of lab-reared P. murrayi with wild-caught P. murrayi and found genes involved in growth and cellular processing were up-regulated in lab-cultured P. murrayi, while a few genes associated with cellular metabolism and freeze tolerance were expressed at relatively lower levels. Preliminary comparative genomic and transcriptomic analyses suggest that the observed constraints on P. murrayi genome architecture and functional gene expression, including genome decay and intron retention, may be an adaptive response to persisting in a biotically simplified, yet consistently physically harsh environment., (© The Author(s) 2020. Published by Oxford University Press on behalf of Genetics Society of America.)