Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae

Functional similarities and differences among subunits of the nascent polypeptide-associated complex (NAC) of Saccharomyces cerevisiae

Protein factors bind ribosomes near the tunnel exit, facilitating protein trafficking and folding. In eukaryotes, the heterodimeric nascent polypeptide-associated complex (NAC) is the most abundant-equimolar to ribosomes. Saccharomyces cerevisiae has a minor β-type subunit (Nacβ2) in addition to abundant Nacβ1, and therefore two NAC heterodimers, α/β1 and α/β12. The additional beta NAC gene arose at the time of the whole genome duplication that occurred in the S. cerevisiae lineage. Nacβ2 has been implicated in regulating the fate of messenger RNA encoding ribosomal protein Rpl4 during translation via its interaction with the Caf130 subunit of the regulatory CCR4-Not complex. We found that Nacβ2 residues just C-terminal to the globular domain are required for its interaction with Caf130 and its negative effect on the growth of cells lacking Acl4, the specialized chaperone for Rpl4. Substitution of these Nacβ2 residues at homologous positions in Nacβ1 results in a chimeric protein that interacts with Caf130 and slows the growth of ∆acl4 cells lacking Nacβ2. Furthermore, alteration of residues in the N-terminus of Nacβ2 or chimeric Nacβ1 previously shown to affect ribosome binding overcomes the growth defect of ∆acl4. Our results are consistent with a model in which Nacβ2's ribosome association per se or its precise positioning is necessary for productive recruitment of CCR4-Not via its interaction with the Caf130 subunit to drive Rpl4 messenger RNA degradation.

Functional genomic analysis of genes important for Candida albicans fitness in diverse environmental conditions

Fungal pathogens such as Candida albicans pose a significant threat to human health with limited treatment options available. One strategy to expand the therapeutic target space is to identify genes important for pathogen growth in host-relevant environments. Here, we leverage a pooled functional genomic screening strategy to identify genes important for fitness of C. albicans in diverse conditions. We identify an essential gene with no known Saccharomyces cerevisiae homolog, C1_09670C, and demonstrate that it encodes subunit 3 of replication factor A (Rfa3). Furthermore, we apply computational analyses to identify functionally coherent gene clusters and predict gene function. Through this approach, we predict the cell-cycle-associated function of C3_06880W, a previously uncharacterized gene required for fitness specifically at elevated temperatures, and follow-up assays confirm that C3_06880W encodes Iml3, a component of the C. albicans kinetochore with roles in virulence in vivo. Overall, this work reveals insights into the vulnerabilities of C. albicans.

Species-tree topology impacts the inference of ancient whole-genome duplications across the angiosperm phylogeny

PREMISE

The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.

METHODS

We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.

RESULTS

The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.

CONCLUSIONS

Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Whole-genome duplication in an algal symbiont bolsters coral heat tolerance

The algal endosymbiont Durusdinium trenchii enhances the resilience of coral reefs under thermal stress. D. trenchii can live freely or in endosymbiosis, and the analysis of genetic markers suggests that this species has undergone whole-genome duplication (WGD). However, the evolutionary mechanisms that underpin the thermotolerance of this species are largely unknown. Here, we present genome assemblies for two D. trenchii isolates, confirm WGD in these taxa, and examine how selection has shaped the duplicated genome regions using gene expression data. We assess how the free-living versus endosymbiotic lifestyles have contributed to the retention and divergence of duplicated genes, and how these processes have enhanced the thermotolerance of D. trenchii. Our combined results suggest that lifestyle is the driver of post-WGD evolution in D. trenchii, with the free-living phase being the most important, followed by endosymbiosis. Adaptations to both lifestyles likely enabled D. trenchii to provide enhanced thermal stress protection to the host coral.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

N-terminal signals in the SNX-BAR paralogs Vps5 and Vin1 guide endosomal coat complex formation

Endosomal coats incorporate membrane-binding subunits such as sorting nexin (SNX) proteins. The Saccharomyces cerevisiae SNX-BAR paralogs Vin1 and Vps5 are respective subunits of the endosomal VINE and retromer complexes whose dimerizing BAR domains are required for complex assembly and membrane association. However, a degree of promiscuity is predicted for yeast BAR-BAR pairings, and recent work has implicated the unstructured N-terminal domains of Vin1 and Vps5 in coat formation. Here, we map N-terminal signals in both SNX-BAR paralogs that contribute to the assembly and function of two distinct endosomal coats in vivo. Whereas Vin1 leverages a polybasic region and adjacent hydrophobic motif to bind Vrl1 and form VINE, the N-terminus of Vps5 interacts with the retromer subunit Vps29 at two sites, including a conserved hydrophobic pocket in Vps29 that engages other accessory proteins in humans. We also examined the sole isoform of Vps5 from the milk yeast Kluyveromyces lactis and found that ancestral yeasts may have used a nested N-terminal signal to form both VINE and retromer. Our results suggest that the specific assembly of Vps5-family SNX-BAR coats depends on inputs from unique N-terminal sequence features in addition to BAR domain coupling, expanding our understanding of endosomal coat biology.

Mlg1, a yeast acyltransferase located in ER membranes associated with mitochondria (MAMs), is involved in de novo synthesis and remodelling of phospholipids

In cells, phospholipids contain acyl chains of variable lengths and saturation, features that affect their functions. Their de novo synthesis in the endoplasmic reticulum takes place via the cytidine diphosphate diacylglycerol (CDP-DAG) and Kennedy pathways, which are conserved in eukaryotes. PA is a key intermediate for all phospholipids (PI, PIPs, PS, PE, PC, PG and CL). The de novo synthesis of PA occurs by acylation of glycerophosphate leading to the synthesis of 1-acyl lysoPA and subsequent acylation of 1-acyl lysoPA at the sn-2 position. Using membranes from Escherichia coli overexpressing MLG1, we showed that the yeast gene MLG1 encodes an acyltransferase, leading specifically to the synthesis of PA from 1-acyl lysoPA. Moreover, after their de novo synthesis, phospholipids can be remodelled by acyl exchange with one and/or two acyl chains exchanged at the sn-1 and/or sn-2 position. Based on shotgun lipidomics of the reference and mlg1Δ strains, as well as biochemical assays for acyltransferase activities, we identified an additional remodelling activity for Mlg1p, namely, incorporation of palmitic acid into the sn-1 position of PS and PE. By using confocal microscopy and subcellular fractionation, we also found that this acyltransferase is located in ER membranes associated with mitochondria, a finding that highlights the importance of these organelles in the global cellular metabolism of lipids.

Genome engineering on size reduction and complexity simplification: A review

BACKGROUND

Genome simplification is an important topic in the field of life sciences that has attracted attention from its conception to the present day. It can help uncover the essential components of the genome and, in turn, shed light on the underlying operating principles of complex biological systems. This has made it a central focus of both basic and applied research in the life sciences. With the recent advancements in related technologies and our increasing knowledge of the genome, now is an opportune time to delve into this topic.

AIM OF REVIEW

Our review investigates the progress of genome simplification from two perspectives: genome size reduction and complexity simplification. In addition, we provide insights into the future development trends of genome simplification.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Reducing genome size requires eliminating non-essential elements as much as possible. This process has been facilitated by advances in genome manipulation and synthesis techniques. However, we still need a better and clearer understanding of living systems to reduce genome complexity. As there is a lack of quantitative and clearly defined standards for this task, we have opted to approach the topic from various perspectives and present our findings accordingly.

Dynamic protein composition of Saccharomyces cerevisiae ribosomes in response to multiple stress conditions reflects alterations in translation activity

Ribosomes, intercellular macromolecules responsible for translation in the cell, are composed of RNAs and proteins. While rRNA makes the scaffold of the ribosome and directs the catalytic steps of protein synthesis, ribosomal proteins play a role in the assembly of the subunits and are essential for the proper structure and function of the ribosome. To date researchers identified heterogeneous ribosomes in different developmental and growth stages. We hypothesized that under stress conditions the heterogeneity of the ribosomes may provide means to prepare the cells for quick recovery. Therefore the aim of the study was the identification of heterogeneity of ribosomal proteins within the ribosomes in response to eleven stress conditions in Saccharomyces cerevisiae, by means of a liquid chromatography/high resolution mass spectrometry (LC-HRMS) and translation activity tests. Out of the total of 74 distinct ribosomal proteins identified in the study 14 small ribosomal subunit (RPS) and 8 large ribosomal subunit (RPL) proteins displayed statistically significant differential abundances within the ribosomes under stress. Additionally, significant alterations in the ratios of 7 ribosomal paralog proteins were observed. Accordingly, the translational activity of yeast ribosomes was altered after UV exposure, during sugar starvation, cold shock, high salt, anaerobic conditions, and amino acid starvation.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Chromosome-level Subgenome-aware de novo Assembly of Saccharomyces bayanus Provides Insight into Genome Divergence after Hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaption, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitated adaptation, we completed chromosome-level de novo assemblies of all 16 pairs chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380 (IFO11022), using Nanopore MinION long-read sequencing. Characterization of S. bayanus subgenomes and comparative analysis with the genomes of its parent species, S. uvarum and S. eubayanus, provide several new insights into understanding genome evolution after a relatively recent hybridization. For instance, multiple recombination events between the two subgenomes have been observed in each chromosome, followed by loss of heterozygosity (LOH) in most chromosomes in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of S. cerevisiae, such as those involved in the maltose metabolism. In addition, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus. The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaption to maltose brewing environments and mitigated the maladaptive effect of hybridization.

Long-lasting redundant gnrh1/3 expression in GnRH neurons enabled apparent switching of paralog usage during evolution

Expressed subtype of paralogous genes in functionally homologous cells sometimes show differences across species, the reasons for which have not been explained. The present study examined hypophysiotropic gonadotropin-releasing hormone (GnRH) neurons in vertebrates to investigate this mechanism. These neurons express either gnrh1 or gnrh3 paralogs, depending on the species, and apparent switching of the expressed paralogs in them occurred at least four times in vertebrate evolution. First, we found redundant expression of gnrh1 and gnrh3 in a single neuron in piranha and hypothesized that it may represent an ancestral GnRH system. Moreover, the gnrh1/gnrh3 enhancer of piranha induced reporter RFP/GFP co-expression in a single hypophysiotropic GnRH neuron in both zebrafish and medaka, whose GnRH neurons only express either gnrh3 or gnrh1. Thus, we propose that redundant expression of gnrh1/3 of relatively recent common ancestors may be the key to apparent switching of the paralog usage among present-day species.

Anti-Cariogenic Effects of S. cerevisiae and S. boulardii in S. mutans-C. albicans Cross-Kingdom In Vitro Models

Despite the well-documented health benefits of the probiotic Saccharomyces, its application in oral health has not been comprehensively assessed. Dental caries is a transmissible disease initiated by acid production of cariogenic bacteria and yeast, such as Streptococcus mutans and Candida albicans, on tooth enamel and followed by subsequent enamel demineralization. Here, we investigated the effect of two Saccharomyces strains (Saccharomyces boulardii and Saccharomyces cerevisiae) on S. mutans-C. albicans cross-kingdom interactions using a cariogenic planktonic model. Viable cells, pH changes, and gene expression were measured. S. cerevisiae and S. boulardii inhibited the growth of C. albicans in dual- and multi-species conditions at 4, 6, and 20 h. Saccharomyces also inhibited C. albicans hyphal formation. Furthermore, Saccharomyces reduced the acidity of the culture medium, which usually plummeted below pH 5 when S. mutans and C. albicans were present in the model. The presence of Saccharomyces maintained the culture medium above 6 even after overnight incubation, demonstrating a protective potential against dental enamel demineralization. S. boulardii significantly down-regulated S. mutans atpD and eno gene expression. Overall, our results shed light on a new promising candidate, Saccharomyces, for dental caries prevention due to its potential to create a less cariogenic environment marked by a neutral pH and reduced growth of C. albicans.

Comparative genomics reveals the dynamic evolutionary history of cement protein genes of barnacles from intertidal to deep-sea hydrothermal vents

Thoracican barnacles are a diverse group of marine organisms for which the availability of genome assemblies is currently limited. In this study, we sequenced the genomes of two neolepadoid species (Ashinkailepas kermadecensis, Imbricaverruca yamaguchii) from hydrothermal vents, in addition to two intertidal species. Genome sizes ranged from 481 to 1054 Mb, with repetitive sequence contents of 21.2% to 50.7%. Concordance rates of orthologs and heterozygosity rates were between 82.4% and 91.7% and between 1.0% and 2.1%, respectively, indicating high genetic diversity and heterozygosity. Based on phylogenomic analyses, we revised the nomenclature of cement genes encoding cement proteins that are not homologous to any known proteins. The major cement gene, CP100A, was found in all thoracican species, including vent-associated neolepadoids, and was hypothesised to be essential for thoracican settlement. Duplicated genes, CP100B and CP100C, were found only in balanids, suggesting potential functional redundancy or acquisition of new functions associated with the calcareous base. An ancestor of CP52 genes was duplicated dynamically among lepadids, pollicipedids with multiple copies on a single scaffold, and balanids with multiple sequential repeats of the conserved regions, but no CP52 genes were found in neolepadoids, providing insights into cement gene evolution among thoracican lineages. This study enhances our understanding of the adhesion mechanisms of thoracicans in underwater environments. The newly sequenced genomes provide opportunities for studying their evolution and ecology, shedding light on their adaptation to diverse marine environments, and contributing to our knowledge of barnacle biology with valuable genomic resources for further studies in this field.

Interplay between acetylation and ubiquitination of imitation switch chromatin remodeler Isw1 confers multidrug resistance in Cryptococcus neoformans

Cryptococcus neoformans poses a threat to human health, but anticryptococcal therapy is hampered by the emergence of drug resistance, whose underlying mechanisms remain poorly understood. Herein, we discovered that Isw1, an imitation switch chromatin remodeling ATPase, functions as a master modulator of genes responsible for in vivo and in vitro multidrug resistance in C. neoformans. Cells with the disrupted ISW1 gene exhibited profound resistance to multiple antifungal drugs. Mass spectrometry analysis revealed that Isw1 is both acetylated and ubiquitinated, suggesting that an interplay between these two modification events exists to govern Isw1 function. Mutagenesis studies of acetylation and ubiquitination sites revealed that the acetylation status of Isw1K97 coordinates with its ubiquitination processes at Isw1K113 and Isw1K441 through modulating the interaction between Isw1 and Cdc4, an E3 ligase. Additionally, clinical isolates of C. neoformans overexpressing the degradation-resistant ISW1K97Q allele showed impaired drug-resistant phenotypes. Collectively, our studies revealed a sophisticated acetylation-Isw1-ubiquitination regulation axis that controls multidrug resistance in C. neoformans.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants

KEY MESSAGE

The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

Neo-functionalization in Saccharomyces cerevisiae: a novel Nrg1-Rtg3 chimeric transcriptional modulator is essential to maintain mitochondrial DNA integrity

In Saccharomyces cerevisiae, the transcriptional repressor Nrg1 (Negative Regulator of Glucose-repressed genes) and the β-Zip transcription factor Rtg3 (ReTroGrade regulation) mediate glucose repression and signalling from the mitochondria to the nucleus, respectively. Here, we show a novel function of these two proteins, in which alanine promotes the formation of a chimeric Nrg1/Rtg3 regulator that represses the ALT2 gene (encoding an alanine transaminase paralog of unknown function). An NRG1/NRG2 paralogous pair, resulting from a post-wide genome small-scale duplication event, is present in the Saccharomyces genus. Neo-functionalization of only one paralog resulted in the ability of Nrg1 to interact with Rtg3. Both nrg1Δ and rtg3Δ single mutant strains were unable to use ethanol and showed a typical petite (small) phenotype on glucose. Neither of the wild-type genes complemented the petite phenotype, suggesting irreversible mitochondrial DNA damage in these mutants. Neither nrg1Δ nor rtg3Δ mutant strains expressed genes encoded by any of the five polycistronic units transcribed from mitochondrial DNA in S. cerevisiae. This, and the direct measurement of the mitochondrial DNA gene complement, confirmed that irreversible damage of the mitochondrial DNA occurred in both mutant strains, which is consistent with the essential role of the chimeric Nrg1/Rtg3 regulator in mitochondrial DNA maintenance.

Impact of PpSpi1, a glycosylphosphatidylinositol-anchored cell wall glycoprotein, on cell wall defects of N-glycosylation-engineered Pichia pastoris

IMPORTANCE

Engineering of biological pathways in various microorganisms is a promising direction for biotechnology. Since the existing microbial cells have evolved over a long period of time, any artificial engineering may cause some unexpected and harmful effects on them. Systematically studying and evaluating these engineered strains are very important and necessary. In order to produce therapeutic proteins with human-like N-glycan structures, much progress has been achieved toward the humanization of N-glycosylation pathways in yeasts. The properties of a P. pastoris strain with humanized N-glycosylation machinery were carefully evaluated in this study. Our work has identified a key glycoprotein (PpSpi1) responsible for the poor growth and morphological defects of this glycoengineered strain. Overexpression of PpSpi1 could significantly rescue the growth defect of the glycoengineered P. pastoris and facilitate its future industrial applications.

Comparing Mitochondrial Activity, Oxidative Stress Tolerance, and Longevity of Thirteen Ascomycota Yeast Species

Aging is characterized by a number of hallmarks including loss of mitochondrial homeostasis and decay in stress tolerance, among others. Unicellular eukaryotes have been widely used to study chronological aging. As a general trait, calorie restriction and activation of mitochondrial respiration has been proposed to contribute to an elongated lifespan. Most aging-related studies have been conducted with the Crabtree-positive yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, and with deletion collections deriving from these conventional yeast models. We have performed an unbiased characterization of longevity using thirteen fungi species, including S. cerevisiae and S. pombe, covering a wide range of the Ascomycota clade. We have determined their mitochondrial activity by oxygen consumption, complex IV activity, and mitochondrial redox potential, and the results derived from these three methodologies are highly overlapping. We have phenotypically compared the lifespans of the thirteen species and their capacity to tolerate oxidative stress. Longevity and elevated tolerance to hydrogen peroxide are correlated in some but not all yeasts. Mitochondrial activity per se cannot anticipate the length of the lifespan. We have classified the strains in four groups, with members of group 1 (Kluyveromyces lactis, Saccharomyces bayanus and Lodderomyces elongisporus) displaying high mitochondrial activity, elevated resistance to oxidative stress, and elongated lifespan.

Paralogous translation factors target distinct mRNAs to differentially regulate tolerance to oxidative stress in yeast

Translation initiation factor 4G (eIF4G) is an integral component of the eIF4F complex which is key to translation initiation for most eukaryotic mRNAs. Many eIF4G isoforms have been described in diverse eukaryotic organisms but we currently have a poor understanding of their functional roles and whether they regulate translation in an mRNA specific manner. The yeast Saccharomyces cerevisiae expresses two eIF4G isoforms, eIF4G1 and eIF4G2, that have previously been considered as functionally redundant with any phenotypic differences arising due to alteration in eIF4G expression levels. Using homogenic strains that express eIF4G1 or eIF4G2 as the sole eIF4G isoforms at comparable expression levels to total eIF4G, we show that eIF4G1 is specifically required to mediate the translational response to oxidative stress. eIF4G1 binds the mRNA cap and remains associated with actively translating ribosomes during oxidative stress conditions and we use quantitative proteomics to show that eIF4G1 promotes oxidative stress-specific proteome changes. eIF4G1, but not eIF4G2, binds the Slf1 LARP protein which appears to mediate the eIF4G1-dependent translational response to oxidative stress. We show similar isoform specific roles for eIF4G in human cells suggesting convergent evolution of multiple eIF4G isoforms offers significant advantages especially where translation must continue under stress conditions.

Streamlining N-terminally anchored yeast surface display via structural insights into S. cerevisiae Pir proteins

Surface display co-opts yeast's innate ability to embellish its cell wall with mannoproteins, thus converting the yeast's outer surface into a growing and self-sustaining catalyst. However, the efficient toolbox for converting the enzyme of interest into its surface-displayed isoform is currently lacking, especially if the isoform needs to be anchored to the cell wall near the isoform's N-terminus, e.g., through a short GPI-independent protein anchor. Aiming to advance such N-terminally anchored surface display, we employed in silico and machine-learning strategies to study the 3D structure, function, genomic organisation, and evolution of the Pir protein family, whose members evolved to covalently attach themselves near their N-terminus to the β-1,3-glucan of the cell wall. Through the newly-gained insights, we rationally engineered 14 S. cerevisiae Hsp150 (Pir2)-based fusion proteins. We quantified their performance, uncovering guidelines for efficient yeast surface display while developing a construct that promoted a 2.5-fold more efficient display of a reporter protein than the full-length Hsp150. Moreover, we developed a Pir-tag, i.e., a peptide spanning only 4.5 kDa but promoting as efficient surface display of a reporter protein as the full-length Hsp150. These constructs fortify the existing surface display toolbox, allowing for a prompt and routine refitting of intracellular proteins into their N-terminally anchored isoforms.

Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

Microtubule search-and-capture model evaluates the effect of chromosomal volume conservation on spindle assembly during mitosis

Variation in the chromosome numbers can arise from the erroneous mitosis or fusion and fission of chromosomes. While the mitotic errors lead to an increase or decrease in the overall chromosomal substance in the daughter cells, fission and fusion keep this conserved. Variations in chromosome numbers are assumed to be a crucial driver of speciation. For example, the members of the muntjac species are known to have very different karyotypes with the chromosome numbers varying from 2n=70+3B in the brown brocket deer to 2n=46 in the Chinese muntjac and 2n=6/7 in the Indian muntjac. The chromosomal content in the nucleus of these closely related mammals is roughly the same and various chromosome fusion and fission pathways have been suggested as the evolution process of these karyotypes. Similar trends can also be found in lepidoptera and yeast species which show a wide variation of chromosome numbers. The effect of chromosome number variation on the spindle assembly time and accuracy is still not properly addressed. We computationally investigate the effect of conservation of the total chromosomal substance on the spindle assembly during prometaphase. Our results suggest that chromosomal fusion pathways aid the microtubule-driven search and capture of the kinetochore in cells with monocentric chromosomes. We further report a comparative analysis of the site and percentage of amphitelic captures, dependence on cell shape, and position of the kinetochore in respect to chromosomal volume partitioning.

The evolution of anaerobic growth in Saccharomycotina yeasts

Humans rely on the ability of budding yeasts to grow without oxygen in industrial scale fermentations that produce beverages, foods, and biofuels. Oxygen is deeply woven into the energy metabolism and biosynthetic capabilities of budding yeasts. While diverse ecological habitats may provide wide varieties of different carbon and nitrogen sources for yeasts to utilize, there is no direct substitute for molecular oxygen, only a range of availability. Understanding how a small subset of budding yeasts evolved the ability to grow without oxygen could expand the set of useful species in industrial scale fermentations as well as provide insight into the cryptic field of yeast ecology. However, we still do not yet appreciate the full breadth of species that can growth without oxygen, what genes underlie this adaptation, and how these genes have evolved.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Mslar: Microbial synthetic lethal and rescue database

Synthetic lethality (SL) occurs when mutations in two genes together lead to cell or organism death, while a single mutation in either gene does not have a significant impact. This concept can also be extended to three or more genes for SL. Computational and experimental methods have been developed to predict and verify SL gene pairs, especially for yeast and Escherichia coli. However, there is currently a lack of a specialized platform to collect microbial SL gene pairs. Therefore, we designed a synthetic interaction database for microbial genetics that collects 13,313 SL and 2,994 Synthetic Rescue (SR) gene pairs that are reported in the literature, as well as 86,981 putative SL pairs got through homologous transfer method in 281 bacterial genomes. Our database website provides multiple functions such as search, browse, visualization, and Blast. Based on the SL interaction data in the S. cerevisiae, we review the issue of duplications' essentiality and observed that the duplicated genes and singletons have a similar ratio of being essential when we consider both individual and SL. The Microbial Synthetic Lethal and Rescue Database (Mslar) is expected to be a useful reference resource for researchers interested in the SL and SR genes of microorganisms. Mslar is open freely to everyone and available on the web at http://guolab.whu.edu.cn/Mslar/.

Direct tests of cytochrome c and c1 functions in the electron transport chain of malaria parasites

The mitochondrial electron transport chain (ETC) of Plasmodium malaria parasites is a major antimalarial drug target, but critical cytochrome (cyt) functions remain unstudied and enigmatic. Parasites express two distinct cyt c homologs (c and c-2) with unusually sparse sequence identity and uncertain fitness contributions. P. falciparum cyt c-2 is the most divergent eukaryotic cyt c homolog currently known and has sequence features predicted to be incompatible with canonical ETC function. We tagged both cyt c homologs and the related cyt c1 for inducible knockdown. Translational repression of cyt c and cyt c1 was lethal to parasites, which died from ETC dysfunction and impaired ubiquinone recycling. In contrast, cyt c-2 knockdown or knockout had little impact on blood-stage growth, indicating that parasites rely fully on the more conserved cyt c for ETC function. Biochemical and structural studies revealed that both cyt c and c-2 are hemylated by holocytochrome c synthase, but UV-vis absorbance and EPR spectra strongly suggest that cyt c-2 has an unusually open active site in which heme is stably coordinated by only a single axial amino acid ligand and can bind exogenous small molecules. These studies provide a direct dissection of cytochrome functions in the ETC of malaria parasites and identify a highly divergent Plasmodium cytochrome c with molecular adaptations that defy a conserved role in eukaryotic evolution.

Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex

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.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole-genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3500 instances where only one of two paralogous proteins undergoes posttranslational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially modified pairs of paralogous proteins. We found that the most common modifications-phosphorylation, ubiquitylation, and acylation but not N-glycosylation-occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility but did mirror known differences in kinase-substrate interactions. Thus, differences in posttranslational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large-scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

Functional genome annotation and transcriptome analysis of Pseudozyma hubeiensis BOT-O, an oleaginous yeast that utilizes glucose and xylose at equal rates

Pseudozyma hubeiensis is a basidiomycete yeast that has the highly desirable traits for lignocellulose valorisation of being equally efficient at utilization of glucose and xylose, and capable of their co-utilization. The species has previously mainly been studied for its capacity to produce secreted biosurfactants in the form of mannosylerythritol lipids, but it is also an oleaginous species capable of accumulating high levels of triacylglycerol storage lipids during nutrient starvation. In this study, we aimed to further characterize the oleaginous nature of P. hubeiensis by evaluating metabolism and gene expression responses during storage lipid formation conditions with glucose or xylose as a carbon source. The genome of the recently isolated P. hubeiensis BOT-O strain was sequenced using MinION long-read sequencing and resulted in the most contiguous P. hubeiensis assembly to date with 18.95 Mb in 31 contigs. Using transcriptome data as experimental support, we generated the first mRNA-supported P. hubeiensis genome annotation and identified 6540 genes. 80% of the predicted genes were assigned functional annotations based on protein homology to other yeasts. Based on the annotation, key metabolic pathways in BOT-O were reconstructed, including pathways for storage lipids, mannosylerythritol lipids and xylose assimilation. BOT-O was confirmed to consume glucose and xylose at equal rates, but during mixed glucose-xylose cultivation glucose was found to be taken up faster. Differential expression analysis revealed that only a total of 122 genes were significantly differentially expressed at a cut-off of |log₂ fold change| ≥ 2 when comparing cultivation on xylose with glucose, during exponential growth and during nitrogen-starvation. Of these 122 genes, a core-set of 24 genes was identified that were differentially expressed at all time points. Nitrogen-starvation resulted in a larger transcriptional effect, with a total of 1179 genes with significant expression changes at the designated fold change cut-off compared with exponential growth on either glucose or xylose.

Reconstruction of hundreds of reference ancestral genomes across the eukaryotic kingdom

Ancestral sequence reconstruction is a fundamental aspect of molecular evolution studies and can trace small-scale sequence modifications through the evolution of genomes and species. In contrast, fine-grained reconstructions of ancestral genome organizations are still in their infancy, limiting our ability to draw comprehensive views of genome and karyotype evolution. Here we reconstruct the detailed gene contents and organizations of 624 ancestral vertebrate, plant, fungi, metazoan and protist genomes, 183 of which are near-complete chromosomal gene order reconstructions. Reconstructed ancestral genomes are similar to their descendants in terms of gene content as expected and agree precisely with reference cytogenetic and in silico reconstructions when available. By comparing successive ancestral genomes along the phylogenetic tree, we estimate the intra- and interchromosomal rearrangement history of all major vertebrate clades at high resolution. This freely available resource introduces the possibility to follow evolutionary processes at genomic scales in chronological order, across multiple clades and without relying on a single extant species as reference.

SCORPiOs, a Novel Method to Reconstruct Gene Phylogenies in the Context of a Known WGD Event

Phylogenetic gene trees recapitulate the evolutionary history of genes across species, forming an essential framework for comparative genomic studies. In particular, within the context of whole-genome duplications (WGDs), they serve as a basis to investigate patterns of duplicate gene retention and loss, timing of genome rediploidization, and, more generally, to explore the functional consequences of the duplication in descending species. Yet, despite ever more sophisticated models to describe the evolution of gene sequences, building accurate gene trees remains a challenge in ancient polyploid taxons. WGDs generate complex gene families with many duplicated copies and recurrent gene losses, which complicate this task even more. Here, we describe how to use SCORPiOs, a novel method that leverages synteny conservation to provide more accurate phylogenies in the presence of a known WGD event.

Evolution of pathogenicity-associated genes in Rhizoctonia solani AG1-IA by genome duplication and transposon-mediated gene function alterations

BACKGROUND

Rhizoctonia solani is a polyphagous fungal pathogen that causes diseases in crops. The fungal strains are classified into anastomosis groups (AGs); however, genomic complexity, diversification into the AGs and the evolution of pathogenicity-associated genes remain poorly understood.

RESULTS

We report a recent whole-genome duplication and sequential segmental duplications in AG1-IA strains of R. solani. Transposable element (TE) clusters have caused loss of synteny in the duplicated blocks and introduced differential structural alterations in the functional domains of several pathogenicity-associated paralogous gene pairs. We demonstrate that the TE-mediated structural variations in a glycosyl hydrolase domain and a GMC oxidoreductase domain in two paralogous pairs affect the pathogenicity of R. solani. Furthermore, to investigate the association of TEs with the natural selection and evolution of pathogenicity, we sequenced the genomes of forty-two rice field isolates of R. solani AG1-IA. The genomic regions with high population mutation rates and with the lowest nucleotide diversity are enriched with TEs. Genetic diversity analysis predicted the genes that are most likely under diversifying and purifying selections. We present evidence that a smaller variant of a glucosamine phosphate N-acetyltransferase (GNAT) protein, predicted to be under purifying selection, and an LPMP_AA9 domain-containing protein, predicted to be under diversifying selection, are important for the successful pathogenesis of R. solani in rice as well as tomato.

CONCLUSIONS

Our study has unravelled whole-genome duplication, TE-mediated neofunctionalization of genes and evolution of pathogenicity traits in R. solani AG1-IA. The pathogenicity-associated genes identified during the study can serve as novel targets for disease control.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.

Direct Tests of Cytochrome Function in the Electron Transport Chain of Malaria Parasites

The mitochondrial electron transport chain (ETC) of Plasmodium malaria parasites is a major antimalarial drug target, but critical cytochrome functions remain unstudied and enigmatic. Parasites express two distinct cyt c homologs ( c and c -2) with unusually sparse sequence identity and uncertain fitness contributions. P. falciparum cyt c -2 is the most divergent eukaryotic cyt c homolog currently known and has sequence features predicted to be incompatible with canonical ETC function. We tagged both cyt c homologs and the related cyt c ₁ for inducible knockdown. Translational repression of cyt c and cyt c ₁ was lethal to parasites, which died from ETC dysfunction and impaired ubiquinone recycling. In contrast, cyt c -2 knockdown or knock-out had little impact on blood-stage growth, indicating that parasites rely fully on the more conserved cyt c for ETC function. Biochemical and structural studies revealed that both cyt c and c -2 are hemylated by holocytochrome c synthase, but UV-vis absorbance and EPR spectra strongly suggest that cyt c -2 has an unusually open active site in which heme is stably coordinated by only a single axial amino-acid ligand and can bind exogenous small molecules. These studies provide a direct dissection of cytochrome functions in the ETC of malaria parasites and identify a highly divergent Plasmodium cytochrome c with molecular adaptations that defy a conserved role in eukaryotic evolution.

SIGNIFICANCE STATEMENT

Mitochondria are critical organelles in eukaryotic cells that drive oxidative metabolism. The mitochondrion of Plasmodium malaria parasites is a major drug target that has many differences from human cells and remains poorly studied. One key difference from humans is that malaria parasites express two cytochrome c proteins that differ significantly from each other and play untested and uncertain roles in the mitochondrial electron transport chain (ETC). Our study revealed that one cyt c is essential for ETC function and parasite viability while the second, more divergent protein has unusual structural and biochemical properties and is not required for growth of blood-stage parasites. This work elucidates key biochemical properties and evolutionary differences in the mitochondrial ETC of malaria parasites.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

WGDI: A user-friendly toolkit for evolutionary analyses of whole-genome duplications and ancestral karyotypes

Evidence of whole-genome duplications (WGDs) and subsequent karyotype changes has been detected in most major lineages of living organisms on Earth. To clarify the complex resulting multi-layered patterns of gene collinearity in genome analyses, there is a need for convenient and accurate toolkits. To meet this need, we developed WGDI (Whole-Genome Duplication Integrated analysis), a Python-based command-line tool that facilitates comprehensive analysis of recursive polyploidization events and cross-species genome alignments. WGDI supports three main workflows (polyploid inference, hierarchical inference of genomic homology, and ancestral chromosome karyotyping) that can improve the detection of WGD and characterization of WGD-related events based on high-quality chromosome-level genomes. Significantly, it can extract complete synteny blocks and facilitate reconstruction of detailed karyotype evolution. This toolkit is freely available at GitHub (https://github.com/SunPengChuan/wgdi). As an example of its application, WGDI convincingly clarified karyotype evolution in Aquilegia coerulea and Vitis vinifera following WGDs and rejected the hypothesis that Aquilegia contributed as a parental lineage to the allopolyploid origin of core dicots.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Systematic analysis of Histidine photosphoto transfer gene family in cotton and functional characterization in response to salt and around tolerance

BACKGROUND

Phosphorylation regulated by the two-component system (TCS) is a very important approach signal transduction in most of living organisms. Histidine phosphotransfer (HP) is one of the important members of the TCS system. Members of the HP gene family have implications in plant stresses tolerance and have been deeply studied in several crops. However, upland cotton is still lacking with complete systematic examination of the HP gene family.

RESULTS

A total of 103 HP gene family members were identified. Multiple sequence alignment and phylogeny of HPs distributed them into 7 clades that contain the highly conserved amino acid residue "XHQXKGSSXS", similar to the Arabidopsis HP protein. Gene duplication relationship showed the expansion of HP gene family being subjected with whole-genome duplication (WGD) in cotton. Varying expression profiles of HPs illustrates their multiple roles under altering environments particularly the abiotic stresses. Analysis is of transcriptome data signifies the important roles played by HP genes against abiotic stresses. Moreover, protein regulatory network analysis and VIGS mediated functional approaches of two HP genes (GhHP23 and GhHP27) supports their predictor roles in salt and drought stress tolerance.

CONCLUSIONS

This study provides new bases for systematic examination of HP genes in upland cotton, which formulated the genetic makeup for their future survey and examination of their potential use in cotton production.

Hybridization order is not the driving factor behind biases in duplicate gene losses among the hexaploid Solanaceae

We model the post-hexaploidy evolution of four genomes from the Solanaceae, a group of flowering plants comprising tomatoes, potatoes and their relatives. The hexaploidy that these genomes descend from occurred through two sequential allopolyploidy events and was marked by the unequal losses of duplicated genes from the different progenitor subgenomes. In contrast with the hexaploid Brassiceae (broccoli and its relatives), where the subgenome with the most surviving genes arrived last in the hexaploidy, among the Solanaceae the most preserved subgenome descends from one of the original two tetraploid progenitors. In fact, the last-arriving subgenome in these plants actually has the fewest surviving genes in the modern genomes. We explore whether the distribution of repetitive elements (REs) in these genomes can explain the biases in gene losses, but while the signals we find are broadly consistent with a role for high RE density in driving gene losses, the REs turn over so quickly that little signal of the RE condition at the time of paleopolyploidy is extant in the modern genomes.

Distinct chromosomal “niches” in the genome of Saccharomyces cerevisiae provide the background for genomic innovation and shape the fate of gene duplicates

Nearly one third of Saccharomyces cerevisiae protein coding sequences correspond to duplicate genes, equally split between small-scale duplicates (SSD) and whole-genome duplicates (WGD). While duplicate genes have distinct properties compared to singletons, to date, there has been no systematic analysis of their positional preferences. In this work, we show that SSD and WGD genes are organized in distinct gene clusters that occupy different genomic regions, with SSD being more peripheral and WGD more centrally positioned close to centromeric chromatin. Duplicate gene clusters differ from the rest of the genome in terms of gene size and spacing, gene expression variability and regulatory complexity, properties that are also shared by singleton genes residing within them. Singletons within duplicate gene clusters have longer promoters, more complex structure and a higher number of protein–protein interactions. Particular chromatin architectures appear to be important for gene evolution, as we find SSD gene-pair co-expression to be strongly associated with the similarity of nucleosome positioning patterns. We propose that specific regions of the yeast genome provide a favourable environment for the generation and maintenance of small-scale gene duplicates, segregating them from WGD-enriched genomic domains. Our findings provide a valuable framework linking genomic innovation with positional genomic preferences.

Sharing the wealth: The versatility of proteins targeted to peroxisomes and other organelles

Peroxisomes are eukaryotic organelles with critical functions in cellular energy and lipid metabolism. Depending on the organism, cell type, and developmental stage, they are involved in numerous other metabolic and regulatory pathways. Many peroxisomal functions require factors also relevant to other cellular compartments. Here, we review proteins shared by peroxisomes and at least one different site within the cell. We discuss the mechanisms to achieve dual targeting, their regulation, and functional consequences. Characterization of dual targeting is fundamental to understand how peroxisomes are integrated into the metabolic and regulatory circuits of eukaryotic cells.

Current Ethanol Production Requirements for the Yeast Saccharomyces cerevisiae

An increase in global energy demand has caused oil prices to reach record levels in recent times. High oil prices together with concerns over CO₂ emissions have resulted in renewed interest in renewable energy. Nowadays, ethanol is the principal renewable biofuel. However, the industrial need for increased productivity, wider substrate range utilization, and the production of novel compounds leads to renewed interest in further extending the use of current industrial strains by exploiting the immense, and still unknown, potential of natural yeast strains. This review seeks to answer the following questions: (a) which characteristics should S. cerevisiae have for the current production of first- and second-generation ethanol? (b) Why are alcohol-tolerance and thermo-tolerance characteristics required? (c) Which genes are related to these characteristics? (d) What are the advances that can be achieved with the isolation of new organisms from the environment?

Can Introns Stabilize Gene Duplication?

Gene duplication is considered one of the most important events that determine the evolution of genomes. However, the neo-duplication condition of a given gene is particularly unstable due to recombination events. Several mechanisms have been proposed to justify this step. In this "opinion article" we propose a role for intron sequences in stabilizing gene duplication by limiting and reducing the identity of the gene sequence between the two duplicated copies. A review of the topic and a detailed hypothesis are presented.