Gene Loss Predictably Drives Evolutionary Adaptation

Evolution of ion channels in cetaceans: a natural experiment in the tree of life

Cetaceans represent a natural experiment within the tree of life in which a lineage changed from terrestrial to aquatic habitats. This shift involved phenotypic modifications, representing an opportunity to explore the genetic bases of phenotypic diversity. Among the different molecular systems that maintain cellular homeostasis, ion channels are crucial for the proper physiological functioning of all living species. This study aims to explore the evolution of ion channels during the evolutionary history of cetaceans. To do so, we created a bioinformatic pipeline to annotate the repertoire of ion channels in the genome of the species included in our sampling. Our main results show that cetaceans have, on average, fewer protein-coding genes and a higher percentage of annotated ion channels than non-cetacean mammals. Signals of positive selection were detected in ion channels related to the heart, locomotion, visual and neurological phenotypes. Interestingly, we predict that the NaV1.5 ion channel of most toothed whales (odontocetes) is sensitive to tetrodotoxin, similar to NaV1.7, given the presence of tyrosine instead of cysteine, in a specific position of the ion channel. Finally, the gene turnover rate of the cetacean crown group is more than three times faster than that of non-cetacean mammals.

The evolutionary novelty of insect defensins: from bacterial killing to toxin neutralization

Insect host defense comprises two complementary dimensions, microbial killing-mediated resistance and microbial toxin neutralization-mediated resilience, both jointly providing protection against pathogen infections. Insect defensins are a class of effectors of innate immunity primarily responsible for resistance to Gram-positive bacteria. Here, we report a newly originated gene from an ancestral defensin via genetic deletion following gene duplication in Drosophila virilis, which confers an enhanced resilience to Gram-positive bacterial infection. This gene encodes an 18-mer arginine-rich peptide (termed DvirARP) with differences from its parent gene in its pattern of expression, structure and function. DvirARP specifically expresses in D. virilis female adults with a constitutive manner. It adopts a novel fold with a 310 helix and a two CXC motif-containing loop stabilized by two disulfide bridges. DvirARP exhibits no activity on the majority of microorganisms tested and only a weak activity against two Gram-positive bacteria. DvirARP knockout flies are viable and have no obvious defect in reproductivity but they are more susceptible to the DvirARP-resistant Staphylococcus aureus infection than the wild type files, which can be attributable to its ability in neutralization of the S. aureus secreted toxins. Phylogenetic distribution analysis reveals that DvirARP is restrictedly present in the Drosophila subgenus, but independent deletion variations also occur in defensins from the Sophophora subgenus, in support of the evolvability of this class of immune effectors. Our work illustrates for the first time how a duplicate resistance-mediated gene evolves an ability to increase the resilience of a subset of Drosophila species against bacterial infection.

Screening and In Silico Analyses of the Yeast Saccharomyces cerevisiae Σ1278b Bank Mutants Using Citral as a Natural Antimicrobial

The antimicrobial function of citral, one of the main compounds of the essential oils (EO) of the Citrus genus, and widely used by the food industry toward spoilage yeast, was previously proven. In this study, the possible mode of action of citral against yeast cells was evaluated by using a global deletome approach. Firstly, the suitability of Saccharomyces cerevisiae Σ1278b to serve as model yeast was assessed by determining its sensitivity to citral (MIC = 0.5 μL/mL). Subsequently, the complete library of Σ1278b haploid mutants deleted in 4019 non-essential genes was screened to identify potential molecular targets of citral. Finally, the deleted genes in the 590 mutants showing increased citral resistance was analyzed with an in-silico approach (Gene Ontology). The significantly enriched GO Terms were "cytoplasm", "vacuole", and "mitochondrion" (cellular components); "catalytic activity" (molecular function); "pseudohyphal growth" (biological process). For molecular function, resistant mutants were grouped into thiosulfate sulfur transferase activity, transferase activity, and oxidoreductase activity; for cellular components, resistant mutants were grouped as: cytoplasm, intracellular organelle, membrane-bounded organelle, mitochondrion, organelle membrane, and vacuole; and finally, with regard to biological process, deleted genes were grouped as: pseudohyphal growth, mitochondrion organization, lipid metabolic process, DNA recombination and repair, and proteolysis. Interestingly, many identified genes were associated with the cellular response to oxidative stress and ROS scavenging. These findings have important implications for the development of citral-based antimicrobials and the elucidation of its mechanism of action.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

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.

Loss-of-function mutations are main drivers of adaptations during short-term evolution

We noticed that during short-term experimental evolution and carcinogenesis, mutations causing gene inactivation (i.e., nonsense mutations or frameshifts) are frequent. Our meta-analysis of 65 experiments using modified dN/dS statistics indicated that nonsense mutations are adaptive in different experimental conditions and we empirically confirmed this prediction. Using yeast S. cerevisiae as a model we show that fixed or highly frequent gene loss-of-function mutations are almost exclusively adaptive in the majority of experiments.

The histone deacetylase Hda1 affects oxidative and osmotic stress response as well as mycoparasitic activity and secondary metabolite biosynthesis in Trichoderma atroviride

The mycoparasitic fungus Trichoderma atroviride is applied in agriculture as a biostimulant and biologic control agent against fungal pathogens that infest crop plants. Secondary metabolites are among the main agents determining the strength and progress of the mycoparasitic attack. However, expression of most secondary metabolism-associated genes requires specific cues, as they are silent under routine laboratory conditions due to their maintenance in an inactive heterochromatin state. Therefore, histone modifications are crucial for the regulation of secondary metabolism. Here, we functionally investigated the role of the class II histone deacetylase encoding gene hda1 of T. atroviride by targeted gene deletion, phenotypic characterization, and multi-omics approaches. Deletion of hda1 did not result in obvious phenotypic alterations but led to an enhanced inhibitory activity of secreted metabolites and reduced mycoparasitic abilities of T. atroviride against the plant-pathogenic fungi Botrytis cinerea and Rhizoctonia solani. The ∆hda1 mutants emitted altered amounts of four volatile organic compounds along their development, produced different metabolite profiles upon growth in liquid culture, and showed a higher susceptibility to oxidative and osmotic stress. Moreover, hda1 deletion affected the expression of several notable gene categories such as polyketide synthases, transcription factors, and genes involved in the HOG MAPK pathway.IMPORTANCEHistone deacetylases play crucial roles in regulating chromatin structure and gene transcription. To date, classical-Zn2+ dependent-fungal histone deacetylases are divided into two classes, of which each comprises orthologues of the two sub-groups Rpd3 and Hos2 and Hda1 and Hos3 of yeast, respectively. However, the role of these chromatin remodelers in mycoparasitic fungi is poorly understood. In this study, we provide evidence that Hda1, the class II histone deacetylases of the mycoparasitic fungus Trichoderma atroviride, regulates its mycoparasitic activity, secondary metabolite biosynthesis, and osmotic and oxidative stress tolerance. The function of Hda1 in regulating bioactive metabolite production and mycoparasitism reveals the importance of chromatin-dependent regulation in the ability of T. atroviride to successfully control fungal plant pathogens.

Deep Conservation and Unexpected Evolutionary History of Neighboring lncRNAs MALAT1 and NEAT1

Long non-coding RNAs (lncRNAs) have begun to receive overdue attention for their regulatory roles in gene expression and other cellular processes. Although most lncRNAs are lowly expressed and tissue-specific, notable exceptions include MALAT1 and its genomic neighbor NEAT1, two highly and ubiquitously expressed oncogenes with roles in transcriptional regulation and RNA splicing. Previous studies have suggested that NEAT1 is found only in mammals, while MALAT1 is present in all gnathostomes (jawed vertebrates) except birds. Here we show that these assertions are incomplete, likely due to the challenges associated with properly identifying these two lncRNAs. Using phylogenetic analysis and structure-aware annotation of publicly available genomic and RNA-seq coverage data, we show that NEAT1 is a common feature of tetrapod genomes except birds and squamates. Conversely, we identify MALAT1 in representative species of all major gnathostome clades, including birds. Our in-depth examination of MALAT1, NEAT1, and their genomic context in a wide range of vertebrate species allows us to reconstruct the series of events that led to the formation of the locus containing these genes in taxa from cartilaginous fish to mammals. This evolutionary history includes the independent loss of NEAT1 in birds and squamates, since NEAT1 is found in the closest living relatives of both clades (crocodilians and tuataras, respectively). These data clarify the origins and relationships of MALAT1 and NEAT1 and highlight an opportunity to study the change and continuity in lncRNA structure and function over deep evolutionary time.

Disruption of Leishmania flagellum attachment zone architecture causes flagellum loss

Leishmania are flagellated eukaryotic parasites that cause leishmaniasis and are closely related to the other kinetoplastid parasites such as Trypanosoma brucei. In all these parasites there is a cell membrane invagination at the base of the flagellum called the flagellar pocket, which is tightly associated with and sculpted by cytoskeletal structures including the flagellum attachment zone (FAZ). The FAZ is a complex interconnected structure linking the flagellum to the cell body and has critical roles in cell morphogenesis, function and pathogenicity. However, this structure varies dramatically in size and organisation between these different parasites, suggesting changes in protein localisation and function. Here, we screened the localisation and function of the Leishmania orthologues of T. brucei FAZ proteins identified in the genome-wide protein tagging project TrypTag. We identified 27 FAZ proteins and our deletion analysis showed that deletion of two FAZ proteins in the flagellum, FAZ27 and FAZ34 resulted in a reduction in cell body size, and flagellum loss in some cells. Furthermore, after null mutant generation, we observed distinct and reproducible changes to cell shape, demonstrating the ability of the parasite to adapt to morphological perturbations resulting from gene deletion. This process of adaptation has important implications for the study of Leishmania mutants.

Genome evolution in plants and the origins of innovation

Plant evolution has been characterised by a series of major novelties in their vegetative and reproductive traits that have led to greater complexity. Underpinning this diversification has been the evolution of the genome. When viewed at the scale of the plant kingdom, plant genome evolution has been punctuated by conspicuous instances of gene and whole-genome duplication, horizontal gene transfer and extensive gene loss. The periods of dynamic genome evolution often coincide with the evolution of key traits, demonstrating the coevolution of plant genomes and phenotypes at a macroevolutionary scale. Conventionally, plant complexity and diversity have been considered through the lens of gene duplication and the role of gene loss in plant evolution remains comparatively unexplored. However, in light of reductive evolution across multiple plant lineages, the association between gene loss and plant phenotypic diversity warrants greater attention.

Experimental evolution for cell biology

Evolutionary cell biology explores the origins, principles, and core functions of cellular features and regulatory networks through the lens of evolution. This emerging field relies heavily on comparative experiments and genomic analyses that focus exclusively on extant diversity and historical events, providing limited opportunities for experimental validation. In this opinion article, we explore the potential for experimental laboratory evolution to augment the evolutionary cell biology toolbox, drawing inspiration from recent studies that combine laboratory evolution with cell biological assays. Primarily focusing on approaches for single cells, we provide a generalizable template for adapting experimental evolution protocols to provide fresh insight into long-standing questions in cell biology.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Evaluation of PEN2-ATP6AP1 axis as an antiparasitic target for metformin based on phylogeny analysis and molecular docking

BACKGROUND

Metformin (Met), the first-line drug used in the treatment for type 2 diabetes mellitus, is effective against a variety of parasites. However, the molecular target of Met at clinical dose against various parasites remains unclear. Recently, low-dose Met (clinical dose) has been reported to directly bind PEN2 (presenilin enhancer protein 2) and initiate the lysosomal glucose-sensing pathway for AMPK activation via ATP6AP1 (V-type proton ATPase subunit S1), rather than perturbing AMP/ATP levels.

METHODS

To explore the possibility of PEN2-ATP6AP1 axis as a drug target of Met for the treatment of parasitic diseases, we identified and characterized orthologs of PEN2 and ATP6AP1 genes in parasites, by constructing phylogenetic trees, analyzing protein sequences and predicting interactions between Met and parasite PEN2.

RESULTS

The results showed that PEN2 and ATP6AP1 genes are only found together in a few of parasite species in the cestoda and nematoda groups. Indicated by molecular simulation, Met might function by interacting with PEN2 on V37/W38/E5 (Trichinella spiralis) with similar binding energy, and on F35/S39 (Caenorhabditis elegans) with higher binding energy, comparing to human PEN2. Hence, these results indicated that only the T. spiralis PEN2-ATP6AP1 axis has the potential to be the direct target of low-concentration Met. Together with contribution of host cells including immune cells in vivo, T. spiralis PEN2-ATP6AP1 axis might play roles in reducing parasite load at low-concentration Met. However, the mechanisms of low-concentration Met on other parasitic infections might be mainly achieved by regulating host cells, rather than directly targeting PEN2-ATP6AP1 axis.

CONCLUSIONS

These findings revealed the potential mechanisms by which Met treats various parasitic diseases, and shed new light on the development of antiparasitic drugs.

Whole genome sequencing of Canadian Saccharomyces cerevisiae strains isolated from spontaneous wine fermentations reveals a new Pacific West Coast Wine clade

Vineyards in wine regions around the world are reservoirs of yeast with oenological potential. Saccharomyces cerevisiae ferments grape sugars to ethanol and generates flavor and aroma compounds in wine. Wineries place a high-value on identifying yeast native to their region to develop a region-specific wine program. Commercial wine strains are genetically very similar due to a population bottleneck and in-breeding compared to the diversity of S. cerevisiae from the wild and other industrial processes. We have isolated and microsatellite-typed hundreds of S. cerevisiae strains from spontaneous fermentations of grapes from the Okanagan Valley wine region in British Columbia, Canada. We chose 75 S. cerevisiae strains, based on our microsatellite clustering data, for whole genome sequencing using Illumina paired-end reads. Phylogenetic analysis shows that British Columbian S. cerevisiae strains cluster into 4 clades: Wine/European, Transpacific Oak, Beer 1/Mixed Origin, and a new clade that we have designated as Pacific West Coast Wine. The Pacific West Coast Wine clade has high nucleotide diversity and shares genomic characteristics with wild North American oak strains but also has gene flow from Wine/European and Ecuadorian clades. We analyzed gene copy number variations to find evidence of domestication and found that strains in the Wine/European and Pacific West Coast Wine clades have gene copy number variation reflective of adaptations to the wine-making environment. The "wine circle/Region B", a cluster of 5 genes acquired by horizontal gene transfer into the genome of commercial wine strains is also present in the majority of the British Columbian strains in the Wine/European clade but in a minority of the Pacific West Coast Wine clade strains. Previous studies have shown that S. cerevisiae strains isolated from Mediterranean Oak trees may be the living ancestors of European wine yeast strains. This study is the first to isolate S. cerevisiae strains with genetic similarity to nonvineyard North American Oak strains from spontaneous wine fermentations.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.

Highly complete long-read genomes reveal pangenomic variation underlying yeast phenotypic diversity

Understanding the genetic causes of trait variation is a primary goal of genetic research. One way that individuals can vary genetically is through variable pangenomic genes: genes that are only present in some individuals in a population. The presence or absence of entire genes could have large effects on trait variation. However, variable pangenomic genes can be missed in standard genotyping workflows, owing to reliance on aligning short-read sequencing to reference genomes. A popular method for studying the genetic basis of trait variation is linkage mapping, which identifies quantitative trait loci (QTLs), regions of the genome that harbor causative genetic variants. Large-scale linkage mapping in the budding yeast Saccharomyces cerevisiae has found thousands of QTLs affecting myriad yeast phenotypes. To enable the resolution of QTLs caused by variable pangenomic genes, we used long-read sequencing to generate highly complete de novo genome assemblies of 16 diverse yeast isolates. With these assemblies, we resolved QTLs for growth on maltose, sucrose, raffinose, and oxidative stress to specific genes that are absent from the reference genome but present in the broader yeast population at appreciable frequency. Copies of genes also duplicate onto chromosomes where they are absent in the reference genome, and we found that these copies generate additional QTLs whose resolution requires pangenome characterization. Our findings show the need for highly complete genome assemblies to identify the genetic basis of trait variation.

Integrated analysis of SR-like protein kinases Sky1 and Sky2 links signaling networks with transcriptional regulation in Candida albicans

Fungal infections are a major global health burden where Candida albicans is among the most common fungal pathogen in humans and is a common cause of invasive candidiasis. Fungal phenotypes, such as those related to morphology, proliferation and virulence are mainly driven by gene expression, which is primarily regulated by kinase signaling cascades. Serine-arginine (SR) protein kinases are highly conserved among eukaryotes and are involved in major transcriptional processes in human and S. cerevisiae. Candida albicans harbors two SR protein kinases, while Sky2 is important for metabolic adaptation, Sky1 has similar functions as in S. cerevisiae. To investigate the role of these SR kinases for the regulation of transcriptional responses in C. albicans, we performed RNA sequencing of sky1Δ and sky2Δ and integrated a comprehensive phosphoproteome dataset of these mutants. Using a Systems Biology approach, we study transcriptional regulation in the context of kinase signaling networks. Transcriptomic enrichment analysis indicates that pathways involved in the regulation of gene expression are downregulated and mitochondrial processes are upregulated in sky1Δ. In sky2Δ, primarily metabolic processes are affected, especially for arginine, and we observed that arginine-induced hyphae formation is impaired in sky2Δ. In addition, our analysis identifies several transcription factors as potential drivers of the transcriptional response. Among these, a core set is shared between both kinase knockouts, but it appears to regulate different subsets of target genes. To elucidate these diverse regulatory patterns, we created network modules by integrating the data of site-specific protein phosphorylation and gene expression with kinase-substrate predictions and protein-protein interactions. These integrated signaling modules reveal shared parts but also highlight specific patterns characteristic for each kinase. Interestingly, the modules contain many proteins involved in fungal morphogenesis and stress response. Accordingly, experimental phenotyping shows a higher resistance to Hygromycin B for sky1Δ. Thus, our study demonstrates that a combination of computational approaches with integration of experimental data can offer a new systems biological perspective on the complex network of signaling and transcription. With that, the investigation of the interface between signaling and transcriptional regulation in C. albicans provides a deeper insight into how cellular mechanisms can shape the phenotype.

Gene loss during a transition to multicellularity

Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited novel interactions, suggesting that gene loss could play a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest gene loss, in addition to gene gain and co-option, may be important for the evolution developmental complexity.

Experimental approaches to study evolutionary cell biology using yeasts

The past century has witnessed tremendous advances in understanding how cells function. Nevertheless, how cellular processes have evolved is still poorly understood. Many studies have highlighted surprising molecular diversity in how cells from diverse species execute the same processes, and advances in comparative genomics are likely to reveal much more molecular diversity than was believed possible until recently. Extant cells remain therefore the product of an evolutionary history that we vastly ignore. Evolutionary cell biology has emerged as a discipline aiming to address this knowledge gap by combining evolutionary, molecular, and cellular biology thinking. Recent studies have shown how even essential molecular processes, such as DNA replication, can undergo fast adaptive evolution under certain laboratory conditions. These developments open new lines of research where the evolution of cellular processes can be investigated experimentally. Yeasts naturally find themselves at the forefront of this research line. Not only do they allow the observation of fast evolutionary adaptation, but they also provide numerous genomic, synthetic, and cellular biology tools already developed by a large community. Here we propose that yeasts can serve as an "evolutionary cell lab" to test hypotheses, principles, and ideas in evolutionary cell biology. We discuss various experimental approaches available for this purpose, and how biology at large can benefit from them.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

The ecological adaptation of the unparalleled plastome character evolution in slipper orchids

Plastomes may have undergone adaptive evolution in the process of plant adaptation to diverse environments, whereby species may differ in plastome characters. Cypripedioideae successfully colonized distinct environments and could be an ideal group for studying the interspecific variation and adaptive evolution of plastomes. Comparative study of plastomes, ancestral state reconstruction, phylogenetic-based analysis, ecological niche modelling, and selective pressure analysis were conducted to reveal the evolutionary patterns of plastomes in Cypripedioideae and their relationship with environmental factors. The plastomes of the three evolved genera had reduced plastome size, increased GC content, and compacted gene content compared to the basal group. Variations in plastome size and GC content are proved to have clear relationships with climate regions. Furthermore, ecological niche modelling revealed that temperature and water factors are important climatic factors contributing to the distributional difference which is directly correlated with the climate regions. The temperature-sensitive genes ndh genes, infA, and rpl20 were found to be either lost/pseudogenized or under positive selection in the evolved groups. Unparalleled plastome character variations were discovered in slipper orchids. Our study indicates that variations in plastome characters have adaptive consequences and that temperature and water factors are important climatic factors that affect plastome evolution. This research highlights the expectation that plants can facilitate adaptation to different environmental conditions with the changes in plastome and has added critical insight for understanding the process of plastome evolution in plants.

Divergent evolutionary trajectories of bryophytes and tracheophytes from a complex common ancestor of land plants

The origin of plants and their colonization of land fundamentally transformed the terrestrial environment. Here we elucidate the basis of this formative episode in Earth history through patterns of lineage, gene and genome evolution. We use new fossil calibrations, a relative clade age calibration (informed by horizontal gene transfer) and new phylogenomic methods for mapping gene family origins. Distinct rooting strategies resolve tracheophytes (vascular plants) and bryophytes (non-vascular plants) as monophyletic sister groups that diverged during the Cambrian, 515-494 million years ago. The embryophyte stem is characterized by a burst of gene innovation, while bryophytes subsequently experienced an equally dramatic episode of reductive genome evolution in which they lost genes associated with the elaboration of vasculature and the stomatal complex. Overall, our analyses reveal that extant tracheophytes and bryophytes are both highly derived from a more complex ancestral land plant. Understanding the origin of land plants requires tracing character evolution across a diversity of modern lineages.

The epitranscriptome in ageing and stress resistance: A systematic review

Modifications of RNA, collectively called the "epitranscriptome", might provide novel biomarkers and innovative targets for interventions in geroscience but are just beginning to be studied in the context of ageing and stress resistance. RNA modifications modulate gene expression by affecting translation initiation and speed, miRNA binding, RNA stability, and RNA degradation. Nonetheless, the precise underlying molecular mechanisms and physiological consequences of most alterations of the epitranscriptome are still only poorly understood. We here systematically review different types of modifications of rRNA, tRNA and mRNA, the methodology to analyze them, current challenges in the field, and human disease associations. Furthermore, we compiled evidence for a connection between individual enzymes, which install RNA modifications, and lifespan in yeast, worm and fly. We also included resistance to different stressors and competitive fitness as search criteria for genes potentially relevant to ageing. Promising candidates identified by this approach include RCM1/NSUN5, RRP8, and F33A8.4/ZCCHC4 that introduce base methylations in rRNA, the methyltransferases DNMT2 and TRM9/ALKBH8, as well as factors involved in the thiolation or A to I editing in tRNA, and finally the m⁶A machinery for mRNA.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Saccharomyces cerevisiae as a Model System for Eukaryotic Cell Biology, from Cell Cycle Control to DNA Damage Response

The yeast Saccharomyces cerevisiae has been used for bread making and beer brewing for thousands of years. In addition, its ease of manipulation, well-annotated genome, expansive molecular toolbox, and its strong conservation of basic eukaryotic biology also make it a prime model for eukaryotic cell biology and genetics. In this review, we discuss the characteristics that made yeast such an extensively used model organism and specifically focus on the DNA damage response pathway as a prime example of how research in S. cerevisiae helped elucidate a highly conserved biological process. In addition, we also highlight differences in the DNA damage response of S. cerevisiae and humans and discuss the challenges of using S. cerevisiae as a model system.

Actin-Related Protein 4 and Linker Histone Sustain Yeast Replicative Ageing

Ageing is accompanied by dramatic changes in chromatin structure organization and genome function. Two essential components of chromatin, the linker histone Hho1p and actin-related protein 4 (Arp4p), have been shown to physically interact in Saccharomyces cerevisiae cells, thus maintaining chromatin dynamics and function, as well as genome stability and cellular morphology. Disrupting this interaction has been proven to influence the stability of the yeast genome and the way cells respond to stress during chronological ageing. It has also been proven that the abrogated interaction between these two chromatin proteins elicited premature ageing phenotypes. Alterations in chromatin compaction have also been associated with replicative ageing, though the main players are not well recognized. Based on this knowledge, here, we examine how the interaction between Hho1p and Arp4p impacts the ageing of mitotically active yeast cells. For this purpose, two sets of strains were used—haploids (WT(n), arp4, hho1Δ and arp4 hho1Δ) and their heterozygous diploid counterparts (WT(2n), ARP4/arp4, HHO1/hho1Δ and ARP4 HHO1/arp4 hho1Δ)—for the performance of extensive morphological and physiological analyses during replicative ageing. These analyses included a comparative examination of the yeast cells’ chromatin structure, proliferative and reproductive potential, and resilience to stress, as well as polysome profiles and chemical composition. The results demonstrated that the haploid chromatin mutants arp4 and arp4 hho1Δ demonstrated a significant reduction in replicative and total lifespan. These findings lead to the conclusion that the importance of a healthy interaction between Arp4p and Hho1p in replicative ageing is significant. This is proof of the concomitant importance of Hho1p and Arp4p in chronological and replicative ageing.

Evolutionary consequences of genomic deletions and insertions in the woolly mammoth genome

Woolly mammoths had a set of adaptations that enabled them to thrive in the Arctic environment. Many mammoth-specific single nucleotide polymorphisms (SNPs) responsible for unique mammoth traits have been previously identified from ancient genomes. However, a multitude of other genetic variants likely contributed to woolly mammoth evolution. In this study, we sequenced two woolly mammoth genomes and combined these with previously sequenced mammoth and elephant genomes to conduct a survey of mammoth-specific deletions and indels. We find that deletions are highly enriched in non-coding regions, suggesting selection against structural variants that affect protein sequences. Nonetheless, at least 87 woolly mammoth genes contain deletions or indels that modify the coding sequence, including genes involved in skeletal morphology and hair growth. These results suggest that deletions and indels contributed to the unique phenotypic adaptations of the woolly mammoth, and were potentially critical to surviving in its natural environment.

•

Two new high-quality woolly mammoth genomes have been generated

•

A new method was used to identify deletions and insertions in woolly mammoths

•

At least 87 genes have been affected by deletions or indels in the mammoth lineage

•

Genes involved in skeletal morphology and hair growth are affected by deletions

Chromosomal-level reference genome assembly of the North American wolverine (Gulo gulo luscus): a resource for conservation genomics

We report a chromosomal-level genome assembly of a male North American wolverine (Gulo gulo luscus) from the Kugluktuk region of Nunavut, Canada. The genome was assembled directly from long-reads, comprising: 758 contigs with a contig N50 of 36.6 Mb; contig L50 of 20; base count of 2.39 Gb; and a near complete representation (99.98%) of the BUSCO 5.2.2 set of 9,226 genes. A presumptive chromosomal-level assembly was generated by scaffolding against two chromosomal-level Mustelidae reference genomes, the ermine and the Eurasian river otter, to derive a final scaffold N50 of 144.0 Mb and a scaffold L50 of 7. We annotated a comprehensive set of genes that have been associated with models of aggressive behavior, a trait which the wolverine is purported to have in the popular literature. To support an integrated, genomics-based wildlife management strategy at a time of environmental disruption from climate change, we annotated the principal genes of the innate immune system to provide a resource to study the wolverine’s susceptibility to new infectious and parasitic diseases. As a resource, we annotated genes involved in the modality of infection by the coronaviruses, an important class of viral pathogens of growing concern as shown by the recent spillover infections by severe acute respiratory syndrome coronavirus-2 to naïve wildlife. Tabulation of heterozygous single nucleotide variants in our specimen revealed a heterozygosity level of 0.065%, indicating a relatively diverse genetic pool that would serve as a baseline for the genomics-based conservation of the wolverine, a rare cold-adapted carnivore now under threat.

Influence of genomic structural variations and nutritional conditions on the emergence of quorum sensing-dependent gene regulation defects in Burkholderia glumae

Bacteria often change their genetic and physiological traits to survive in harsh environments. To determine whether, in various strains of Burkholderia glumae, genomic diversity is associated with the ability to adapt to ever-changing environments, whole genomes of 44 isolates from different hosts and regions were analyzed. Whole-genome phylogenetic analysis of the 44 isolates revealed six clusters and two divisions. While all isolates possessed chromosomes 1 and 2, strains BGR80S and BGR81S had one chromosome resulting from the merging of the two chromosomes. Upon comparison of genomic structures to the prototype BGR1, inversions, deletions, and rearrangements were found within or between chromosomes 1 and/or 2 in the other isolates. When three isolates—BGR80S, BGR15S, and BGR21S, representing clusters III, IV, and VI, respectively—were grown in Luria-Bertani medium, spontaneous null mutations were identified in qsmR encoding a quorum-sensing master regulator. Six days after subculture, qsmR mutants were found at detectable frequencies in BGR15S and BGR21S, and reached approximately 40% at 8 days after subculture. However, the qsmR mutants appeared 2 days after subculture in BGR80S and dominated the population, reaching almost 80%. No qsmR mutant was detected at detectable frequency in BGR1 or BGR13S. The spontaneous qsmR mutants outcompeted their parental strains in the co-culture. Daily addition of glucose or casamino acids to the batch cultures of BGR80S delayed emergence of qsmR mutants and significantly reduced their incidence. These results indicate that spontaneous qsmR mutations are correlated with genomic structures and nutritional conditions.

Gene loss and compensatory evolution promotes the emergence of morphological novelties in budding yeast

Deleterious mutations are generally considered to be irrelevant for morphological evolution. However, they could be compensated by conditionally beneficial mutations, thereby providing access to new adaptive paths. Here we use high-dimensional phenotyping of laboratory-evolved budding yeast lineages to demonstrate that new cellular morphologies emerge exceptionally rapidly as a by-product of gene loss and subsequent compensatory evolution. Unexpectedly, the capacities for invasive growth, multicellular aggregation and biofilm formation also spontaneously evolve in response to gene loss. These multicellular phenotypes can be achieved by diverse mutational routes and without reactivating the canonical regulatory pathways. These ecologically and clinically relevant traits originate as pleiotropic side effects of compensatory evolution and have no obvious utility in the laboratory environment. The extent of morphological diversity in the evolved lineages is comparable to that of natural yeast isolates with diverse genetic backgrounds and lifestyles. Finally, we show that both the initial gene loss and subsequent compensatory mutations contribute to new morphologies, with their synergistic effects underlying specific morphological changes. We conclude that compensatory evolution is a previously unrecognized source of morphological diversity and phenotypic novelties.

Evolutionary Adaptation by Repetitive Long-Term Cultivation with Gradual Increase in Temperature for Acquiring Multi-Stress Tolerance and High Ethanol Productivity in Kluyveromyces marxianus DMKU 3-1042

During ethanol fermentation, yeast cells are exposed to various stresses that have negative effects on cell growth, cell survival, and fermentation ability. This study, therefore, aims to develop Kluyveromyces marxianus-adapted strains that are multi-stress tolerant and to increase ethanol production at high temperatures through a novel evolutionary adaptation procedure. K. marxianus DMKU 3-1042 was subjected to repetitive long-term cultivation with gradual increases in temperature (RLCGT), which exposed cells to various stresses, including high temperatures. In each cultivation step, 1% of the previous culture was inoculated into a medium containing 1% yeast extract, 2% peptone, and 2% glucose, and cultivation was performed under a shaking condition. Four adapted strains showed increased tolerance to ethanol, furfural, hydroxymethylfurfural, and vanillin, and they also showed higher production of ethanol in a medium containing 16% glucose at high temperatures. One showed stronger ethanol tolerance. Others had similar phenotypes, including acetic acid tolerance, though genome analysis revealed that they had different mutations. Based on genome and transcriptome analyses, we discuss possible mechanisms of stress tolerance in adapted strains. All adapted strains gained a useful capacity for ethanol fermentation at high temperatures and improved tolerance to multi-stress. This suggests that RLCGT is a simple and efficient procedure for the development of robust strains.

Convergent Loss of Chemoreceptors across Independent Origins of Slave-Making in Ants

The evolution of an obligate parasitic lifestyle often leads to the reduction of morphological and physiological traits, which may be accompanied by loss of genes and functions. Slave-making ants are social parasites that exploit the work force of closely related ant species for social behaviors such as brood care and foraging. Recent divergence between these social parasites and their hosts enables comparative studies of gene family evolution. We sequenced the genomes of eight ant species, representing three independent origins of ant slavery. During the evolution of eusociality, chemoreceptor genes multiplied due to the importance of chemical communication in insect societies. We investigated the evolutionary fate of these chemoreceptors and found that slave-making ant genomes harbored only half as many gustatory receptors as their hosts', potentially mirroring the outsourcing of foraging tasks to host workers. In addition, parasites had fewer odorant receptors and their loss shows striking patterns of convergence across independent origins of parasitism, in particular in orthologs often implicated in sociality like the 9-exon odorant receptors. These convergent losses represent a rare case of convergent molecular evolution at the level of individual genes. Thus, evolution can operate in a way that is both repeatable and reversible when independent ant lineages lose important social traits during the transition to a parasitic lifestyle.

Intraspecific Variation in the Rates of Mutations Causing Structural Variation in Daphnia magna

Mutations that cause structural variation are important sources of genetic variation upon which other evolutionary forces can act, however, they are difficult to observe and therefore few direct estimates of their rate and spectrum are available. Understanding mutation rate evolution, however, requires adding to the limited number of species for which direct estimates are available, quantifying levels of intraspecific variation in mutation rates, and assessing whether rate estimates co-vary across types of mutation. Here, we report structural variation-causing mutation rates (svcMRs) for six categories of mutations (short insertions and deletions, long deletions and duplications, and deletions and duplications at copy number variable sites) from nine genotypes of Daphnia magna collected from three populations in Finland, Germany, and Israel using a mutation accumulation approach. Based on whole-genome sequence data and validated using simulations, we find svcMRs are high (two orders of magnitude higher than base substitution mutation rates measured in the same lineages), highly variable among populations, and uncorrelated across categories of mutation. Furthermore, to assess the impact of scvMRs on the genome, we calculated rates while adjusting for the lengths of events and ran simulations to determine if the mutations occur in genic regions more or less frequently than expected by chance. Our results pose a challenge to most prevailing theories aimed at explaining the evolution of the mutation rate, underscoring the importance of obtaining additional mutation rate estimates in more genotypes, for more types of mutation, in more species, in order to improve our future understanding of mutation rates, their variation, and their evolution.

Adaptations and Diversity of Antarctic Fishes: A Genomic Perspective

Antarctic notothenioid fishes are the classic example of vertebrate adaptive radiation in a marine environment. Notothenioids diversified from a single common ancestor ∼25 Mya to more than 140 species today, and they represent ∼90% of fish biomass on the continental shelf of Antarctica. As they diversified in the cold Southern Ocean, notothenioids evolved numerous traits, including osteopenia, anemia, cardiomegaly, dyslipidemia, and aglomerular kidneys, that are beneficial or tolerated in their environment but are pathological in humans. Thus, notothenioids are models for understanding adaptive radiations, physiological and biochemical adaptations to extreme environments, and genetic mechanisms of human disease. Since 2014, 16 notothenioid genomes have been published, which enable a first-pass holistic analysis of the notothenioid radiation and the genetic underpinnings of novel notothenioid traits. Here, we review the notothenioid radiation from a genomic perspective and integrate our insights with recent observations from other fish radiations. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Ploidy and recombination proficiency shape the evolutionary adaptation to constitutive DNA replication stress

In haploid budding yeast, evolutionary adaptation to constitutive DNA replication stress alters three genome maintenance modules: DNA replication, the DNA damage checkpoint, and sister chromatid cohesion. We asked how these trajectories depend on genomic features by comparing the adaptation in three strains: haploids, diploids, and recombination deficient haploids. In all three, adaptation happens within 1000 generations at rates that are correlated with the initial fitness defect of the ancestors. Mutations in individual genes are selected at different frequencies in populations with different genomic features, but the benefits these mutations confer are similar in the three strains, and combinations of these mutations reproduce the fitness gains of evolved populations. Despite the differences in the selected mutations, adaptation targets the same three functional modules in strains with different genomic features, revealing a common evolutionary response to constitutive DNA replication stress.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Functional or Vestigial? The Genomics of the Pineal Gland in Xenarthra

Vestigial organs are historical echoes of past phenotypes. Determining whether a specific organ constitutes a functional or vestigial structure can be a challenging task, given that distinct levels of atrophy may arise between and within lineages. The mammalian pineal gland, an endocrine organ involved in melatonin biorhythmicity, represents a classic example, often yielding contradicting anatomical observations. In Xenarthra (sloths, anteaters, and armadillos), a peculiar mammalian order, the presence of a distinct pineal organ was clearly observed in some species (i.e., Linnaeus's two-toed sloth), but undetected in other closely related species (i.e., brown-throated sloth). In the nine-banded armadillo, contradicting evidence supports either functional or vestigial scenarios. Thus, to untangle the physiological status of the pineal gland in Xenarthra, we used a genomic approach to investigate the evolution of the gene hub responsible for melatonin synthesis and signaling. We show that both synthesis and signaling compartments are eroded and were probably lost independently among Xenarthra orders. Additionally, by expanding our analysis to 157 mammal genomes, we offer a comprehensive view showing that species with very distinctive habitats and lifestyles have convergently evolved a similar phenotype: Cetacea, Pholidota, Dermoptera, Sirenia, and Xenarthra. Our findings suggest that the recurrent inactivation of melatonin genes correlates with pineal atrophy and endorses the use of genomic analyses to ascertain the physiological status of suspected vestigial structures.

Parallel in vivo experimental evolution reveals that increased stress resistance was key for the emergence of persistent tuberculosis bacilli

Pathogenomic evidence suggests that Mycobacterium tuberculosis (MTB) evolved from an environmental ancestor similar to Mycobacterium canettii, a rare human pathogen. Although the adaptations responsible for this transition are poorly characterized, the ability to persist in humans seems to be important. We set out to identify the adaptations contributing to the evolution of persistence in MTB. We performed an experimental evolution of eight M. canettii populations in mice; four populations were derived from the isolate STB-K (phylogenomically furthest from MTB) and four from STB-D (closest to MTB), which were monitored for 15 and 6 cycles, respectively. We selected M. canettii mutants with enhanced persistence in vivo compared with the parental strains, which were phenotypically closer to MTB. Genome sequencing of 140 mutants and complementation analysis revealed that mutations in two loci were responsible for enhanced persistence. Most of the tested mutants were more resistant than their parental strains to nitric oxide, an important effector of immunity. Modern MTB were similarly more resistant to nitric oxide than M. canettii. Our findings demonstrate phenotypic convergence during experimental evolution of M. canettii, which mirrors natural evolution of MTB. Furthermore, they indicate that the ability to withstand host-induced stresses was key for the emergence of persistent MTB.

Massive Gene Loss and Function Shuffling in Appendicularians Stretch the Boundaries of Chordate Wnt Family Evolution

Gene loss is a pervasive source of genetic variation that influences species evolvability, biodiversity and the innovation of evolutionary adaptations. To better understand the evolutionary patterns and impact of gene loss, here we investigate as a case study the evolution of the wingless (Wnt) family in the appendicularian tunicate Oikopleura dioica, an emergent EvoDevo model characterized by its proneness to lose genes among chordates. Genome survey and phylogenetic analyses reveal that only four of the thirteen Wnt subfamilies have survived in O. dioica—Wnt5, Wnt10, Wnt11, and Wnt16,—representing the minimal Wnt repertoire described in chordates. While the loss of Wnt4 and Wnt8 likely occurred in the last common ancestor of tunicates, representing therefore a synapomorphy of this subphylum, the rest of losses occurred during the evolution of appendicularians. This work provides the first complete Wnt developmental expression atlas in a tunicate and the first insights into the evolution of Wnt developmental functions in appendicularians. Our work highlights three main evolutionary patterns of gene loss: (1) conservation of ancestral Wnt expression domains not affected by gene losses; (2) function shuffling among Wnt paralogs accompanied by gene losses; and (3) extinction of Wnt expression in certain embryonic directly correlated with gene losses. Overall our work reveals that in contrast to “conservative” pattern of evolution of cephalochordates and vertebrates, O. dioica shows an even more radical “liberal” evolutionary pattern than that described ascidian tunicates, stretching the boundaries of the malleability of Wnt family evolution in chordates.

Genome-resolved metagenomics suggests a mutualistic relationship between Mycoplasma and salmonid hosts

Salmonids are important sources of protein for a large proportion of the human population. Mycoplasma species are a major constituent of the gut microbiota of salmonids, often representing the majority of microbiota. Despite the frequent reported dominance of salmonid-related Mycoplasma species, little is known about the phylogenomic placement, functions and potential evolutionary relationships with their salmonid hosts. In this study, we utilise 2.9 billion metagenomic reads generated from 12 samples from three different salmonid host species to I) characterise and curate the first metagenome-assembled genomes (MAGs) of Mycoplasma dominating the intestines of three different salmonid species, II) establish the phylogeny of these salmonid candidate Mycoplasma species, III) perform a comprehensive pangenomic analysis of Mycoplasma, IV) decipher the putative functionalities of the salmonid MAGs and reveal specific functions expected to benefit the host. Our data provide a basis for future studies examining the composition and function of the salmonid microbiota.

The Pectobacterium pangenome, with a focus on Pectobacterium brasiliense, shows a robust core and extensive exchange of genes from a shared gene pool

BACKGROUND

Bacterial plant pathogens of the Pectobacterium genus are responsible for a wide spectrum of diseases in plants, including important crops such as potato, tomato, lettuce, and banana. Investigation of the genetic diversity underlying virulence and host specificity can be performed at genome level by using a comprehensive comparative approach called pangenomics. A pangenomic approach, using newly developed functionalities in PanTools, was applied to analyze the complex phylogeny of the Pectobacterium genus. We specifically used the pangenome to investigate genetic differences between virulent and avirulent strains of P. brasiliense, a potato blackleg causing species dominantly present in Western Europe.

RESULTS

Here we generated a multilevel pangenome for Pectobacterium, comprising 197 strains across 19 species, including type strains, with a focus on P. brasiliense. The extensive phylogenetic analysis of the Pectobacterium genus showed robust distinct clades, with most detail provided by 452,388 parsimony-informative single-nucleotide polymorphisms identified in single-copy orthologs. The average Pectobacterium genome consists of 47% core genes, 1% unique genes, and 52% accessory genes. Using the pangenome, we zoomed in on differences between virulent and avirulent P. brasiliense strains and identified 86 genes associated to virulent strains. We found that the organization of genes is highly structured and linked with gene conservation, function, and transcriptional orientation.

CONCLUSION

The pangenome analysis demonstrates that evolution in Pectobacteria is a highly dynamic process, including gene acquisitions partly in clusters, genome rearrangements, and loss of genes. Pectobacterium species are typically not characterized by a set of species-specific genes, but instead present themselves using new gene combinations from the shared gene pool. A multilevel pangenomic approach, fusing DNA, protein, biological function, taxonomic group, and phenotypes, facilitates studies in a flexible taxonomic context.

LI Detector: a framework for sensitive colony-based screens regardless of the distribution of fitness effects

Microbial growth characteristics have long been used to investigate fundamental questions of biology. Colony-based high-throughput screens enable parallel fitness estimation of thousands of individual strains using colony growth as a proxy for fitness. However, fitness estimation is complicated by spatial biases affecting colony growth, including uneven nutrient distribution, agar surface irregularities, and batch effects. Analytical methods that have been developed to correct for these spatial biases rely on the following assumptions: (1) that fitness effects are normally distributed, and (2) that most genetic perturbations lead to minor changes in fitness. Although reasonable for many applications, these assumptions are not always warranted and can limit the ability to detect small fitness effects. Beneficial fitness effects, in particular, are notoriously difficult to detect under these assumptions. Here, we developed the linear interpolation-based detector (LI Detector) framework to enable sensitive colony-based screening without making prior assumptions about the underlying distribution of fitness effects. The LI Detector uses a grid of reference colonies to assign a relative fitness value to every colony on the plate. We show that the LI Detector is effective in correcting for spatial biases and equally sensitive toward increase and decrease in fitness. LI Detector offers a tunable system that allows the user to identify small fitness effects with unprecedented sensitivity and specificity. LI Detector can be utilized to develop and refine gene-gene and gene-environment interaction networks of colony-forming organisms, including yeast, by increasing the range of fitness effects that can be reliably detected.