Increased glycolytic flux as an outcome of whole-genome duplication in yeast

Highly parallelized laboratory evolution of wine yeasts for enhanced metabolic phenotypes

Adaptive Laboratory Evolution (ALE) of microorganisms can improve the efficiency of sustainable industrial processes important to the global economy. However, stochasticity and genetic background effects often lead to suboptimal outcomes during laboratory evolution. Here we report an ALE platform to circumvent these shortcomings through parallelized clonal evolution at an unprecedented scale. Using this platform, we evolved 104 yeast populations in parallel from many strains for eight desired wine fermentation-related traits. Expansions of both ALE replicates and lineage numbers broadened the evolutionary search spectrum leading to improved wine yeasts unencumbered by unwanted side effects. At the genomic level, evolutionary gains in metabolic characteristics often coincided with distinct chromosome amplifications and the emergence of side-effect syndromes that were characteristic of each selection niche. Several high-performing ALE strains exhibited desired wine fermentation kinetics when tested in larger liquid cultures, supporting their suitability for application. More broadly, our high-throughput ALE platform opens opportunities for rapid optimization of microbes which otherwise could take many years to accomplish.

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.

Genome sequencing of Syzygium cumini (jamun) reveals adaptive evolution in secondary metabolism pathways associated with its medicinal properties

Syzygium cumini, also known as jambolan or jamun, is an evergreen tree widely known for its medicinal properties, fruits, and ornamental value. To understand the genomic and evolutionary basis of its medicinal properties, we sequenced S. cumini genome for the first time from the world's largest tree genus Syzygium using Oxford Nanopore and 10x Genomics sequencing technologies. We also sequenced and assembled the transcriptome of S. cumini in this study. The tetraploid and highly heterozygous draft genome of S. cumini had a total size of 709.9 Mbp with 61,195 coding genes. The phylogenetic position of S. cumini was established using a comprehensive genome-wide analysis including species from 18 Eudicot plant orders. The existence of neopolyploidy in S. cumini was evident from the higher number of coding genes and expanded gene families resulting from gene duplication events compared to the other two sequenced species from this genus. Comparative evolutionary analyses showed the adaptive evolution of genes involved in the phenylpropanoid-flavonoid (PF) biosynthesis pathway and other secondary metabolites biosynthesis such as terpenoid and alkaloid in S. cumini, along with genes involved in stress tolerance mechanisms, which was also supported by leaf transcriptome data generated in this study. The adaptive evolution of secondary metabolism pathways is associated with the wide range of pharmacological properties, specifically the anti-diabetic property, of this species conferred by the bioactive compounds that act as nutraceutical agents in modern medicine.

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.

Optimizing growth and biomass production of non-Saccharomyces wine yeast starters by overcoming sucrose consumption deficiency

The use of non-Saccharomyces yeasts as starters in winemaking has increased exponentially in the last years. For instance, non-conventional yeasts have proven useful for the improvement of the organoleptic profile and biocontrol. Active dry yeast starter production has been optimized for Saccharomyces cerevisiae, which may entail problems for the propagation of non-Saccharomyces yeasts. This work shows that the poor growth of Hanseniaspora vineae and Metschnikowia pulcherrima in molasses is related to a deficient sucrose consumption, linked to their low invertase activity. In order to address this issue, simple modifications to the cultivation media based hydrolysis and the reduction of sucrose concentration were performed. We performed biomass propagation simulations at a bench-top and bioreactor scale. The results show that cultivation in a hexose-based media improved biomass production in both species, as it solves their low invertase activity. The reduction in sugar concentration promoted a metabolic shift to a respiratory metabolism, which allowed a higher biomass yield, but did not improve total biomass production, due to the lower sugar availability. To evaluate the technological performance of these adaptations, we performed mixed grape juice fermentations with biomass produced in such conditions of M. pulcherrima and S. cerevisiae. The analysis of wines produced revealed that the different treatments we have tested did not have any negative impact on wine quality, further proving their applicability at an industrial level for the improvement of biomass production.

Chromosome-level reference genome of Tetrastigma hemsleyanum (Vitaceae) provides insights into genomic evolution and the biosynthesis of phenylpropanoids and flavonoids

Here, we present a high-quality chromosome-scale genome assembly (2.19 Gb) and annotation of Tetrastigma hemsleyanum, a perennial herbaceous liana native to subtropical China with diverse medicinal applications. Approximately 73% of the genome was comprised of transposable elements (TEs), of which long terminal repeat retrotransposons (LTR-RTs) were a predominant group (69% of the genome). The genome size increase of T. hemsleyanum (relative to Vitis species) was mostly due to the proliferation of LTR-RTs. Of the different modes of gene duplication identified, transposed duplication (TRD) and dispersed duplication (DSD) were the predominant ones. Genes, particularly those involved in the phenylpropanoid-flavonoid (PF) pathway and those associated with therapeutic properties and environmental stress resistance, were significantly amplified through recent tandem duplications. We dated the divergence of two intraspecific lineages in Southwest (SW) versus Central-South-East (CSE) China to the late Miocene (approximately 5.2 million years ago). Of those, the former showed more upregulated genes and metabolites. Based on resequencing data of 38 individuals representing both lineages, we identified various candidate genes related to 'response to stimulus' and 'biosynthetic process', including ThFLS11, which is putatively involved in flavonoid accumulation. Overall, this study provides abundant genomic resources for future evolutionary, ecological, and functional genomics studies in T. hemsleyanum and related species.

Engineering proton-coupled hexose uptake in Saccharomyces cerevisiae for improved ethanol yield

Background

In the yeast Saccharomyces cerevisiae, which is widely applied for industrial bioethanol production, uptake of hexoses is mediated by transporters with a facilitated diffusion mechanism. In anaerobic cultures, a higher ethanol yield can be achieved when transport of hexoses is proton-coupled, because of the lower net ATP yield of sugar dissimilation. In this study, the facilitated diffusion transport system for hexose sugars of S. cerevisiae was replaced by hexose–proton symport.

Results

Introduction of heterologous glucose– or fructose–proton symporters in an hxt⁰ yeast background strain (derived from CEN.PK2-1C) restored growth on the corresponding sugar under aerobic conditions. After applying an evolutionary engineering strategy to enable anaerobic growth, the hexose–proton symporter-expressing strains were grown in anaerobic, hexose-limited chemostats on synthetic defined medium, which showed that the biomass yield of the resulting strains was decreased by 44.0-47.6%, whereas the ethanol yield had increased by up to 17.2% (from 1.51 to 1.77 mol mol hexose⁻¹) compared to an isogenic strain expressing the hexose uniporter HXT5. To apply this strategy to increase the ethanol yield on sucrose, we constructed a platform strain in which all genes encoding hexose transporters, disaccharide transporters and disaccharide hydrolases were deleted, after which a combination of a glucose–proton symporter, fructose–proton symporter and extracellular invertase (SUC2) were introduced. After evolution, the resulting strain exhibited a 16.6% increased anaerobic ethanol yield (from 1.51 to 1.76 mol mol hexose equivalent⁻¹) and 46.6% decreased biomass yield on sucrose.

Conclusions

This study provides a proof-of-concept for the replacement of the endogenous hexose transporters of S. cerevisiae by hexose-proton symport, and the concomitant decrease in ATP yield, to greatly improve the anaerobic yield of ethanol on sugar. Moreover, the sugar-negative platform strain constructed in this study acts as a valuable starting point for future studies on sugar transport or development of cell factories requiring specific sugar transport mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02145-7.

Genome doubling enabled the expansion of yeast vesicle traffic pathways

Vesicle budding and fusion in eukaryotes depend on a suite of protein types, such as Arfs, Rabs, coats and SNAREs. Distinct paralogs of these proteins act at distinct intracellular locations, suggesting a link between gene duplication and the expansion of vesicle traffic pathways. Genome doubling, a common source of paralogous genes in fungi, provides an ideal setting in which to explore this link. Here we trace the fates of paralog doublets derived from the 100-Ma-old hybridization event that gave rise to the whole genome duplication clade of budding yeast. We find that paralog doublets involved in specific vesicle traffic functions and pathways are convergently retained across the entire clade. Vesicle coats and adaptors involved in secretory and early-endocytic pathways are retained as doublets, at rates several-fold higher than expected by chance. Proteins involved in later endocytic steps and intra-Golgi traffic, including the entire set of multi-subunit and coiled-coil tethers, have reverted to singletons. These patterns demonstrate that selection has acted to expand and diversify the yeast vesicle traffic apparatus, across species and time.

A chromosome-level genome assembly of an alpine plant Crucihimalaya lasiocarpa provides insights into high-altitude adaptation

It remains largely unknown how plants adapt to high-altitude habitats. Crucihimalaya (Brassicaceae) is an alpine genus occurring in the Qinghai-Tibet Plateau characterized by cold temperatures and strong ultraviolet radiation. Here, we generated a chromosome-level genome for C. lasiocarpa with a total size of 255.8 Mb and a scaffold N50 size of 31.9 Mb. We first examined the karyotype origin of this species and found that the karyotype of five chromosomes resembled the ancestral karyotype of the Brassicaceae family, while the other three showed strong chromosomal structural variations. In combination with the rough genome sequence of another congener (C. himalaica), we found that the significantly expanded gene families and positively selected genes involved in alpine adaptation have occurred since the origin of this genus. Our new findings provide valuable information for the chromosomal karyotype evolution of Brassicaceae and investigations of high-altitude environment adaptation of the genus.

Neofunctionalization of Glycolytic Enzymes: An Evolutionary Route to Plant Parasitism in the Oomycete Phytophthora nicotianae

Oomycetes, of the genus Phytophthora, comprise of some of the most devastating plant pathogens. Parasitism of Phytophthora results from evolution from an autotrophic ancestor and adaptation to a wide range of environments, involving metabolic adaptation. Sequence mining showed that Phytophthora spp. display an unusual repertoire of glycolytic enzymes, made of multigene families and enzyme replacements. To investigate the impact of these gene duplications on the biology of Phytophthora and, eventually, identify novel functions associated to gene expansion, we focused our study on the first glycolytic step on P. nicotianae, a broad host range pathogen. We reveal that this step is committed by a set of three glucokinase types that differ by their structure, enzymatic properties, and evolutionary histories. In addition, they are expressed differentially during the P. nicotianae life cycle, including plant infection. Last, we show that there is a strong association between the expression of a glucokinase member in planta and extent of plant infection. Together, these results suggest that metabolic adaptation is a component of the processes underlying evolution of parasitism in Phytophthora, which may possibly involve the neofunctionalization of metabolic enzymes.

Kinetic Modeling of Saccharomyces cerevisiae Central Carbon Metabolism: Achievements, Limitations, and Opportunities

Central carbon metabolism comprises the metabolic pathways in the cell that process nutrients into energy, building blocks and byproducts. To unravel the regulation of this network upon glucose perturbation, several metabolic models have been developed for the microorganism Saccharomyces cerevisiae. These dynamic representations have focused on glycolysis and answered multiple research questions, but no commonly applicable model has been presented. This review systematically evaluates the literature to describe the current advances, limitations, and opportunities. Different kinetic models have unraveled key kinetic glycolytic mechanisms. Nevertheless, some uncertainties regarding model topology and parameter values still limit the application to specific cases. Progressive improvements in experimental measurement technologies as well as advances in computational tools create new opportunities to further extend the model scale. Notably, models need to be made more complex to consider the multiple layers of glycolytic regulation and external physiological variables regulating the bioprocess, opening new possibilities for extrapolation and validation. Finally, the onset of new data representative of individual cells will cause these models to evolve from depicting an average cell in an industrial fermenter, to characterizing the heterogeneity of the population, opening new and unseen possibilities for industrial fermentation improvement.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Adaptations in metabolism and protein translation give rise to the Crabtree effect in yeast

Aerobic fermentation, also referred to as the Crabtree effect in yeast, is a well-studied phenomenon that allows many eukaryal cells to attain higher growth rates at high glucose availability. Not all yeasts exhibit the Crabtree effect, and it is not known why Crabtree-negative yeasts can grow at rates comparable to Crabtree-positive yeasts. Here, we quantitatively compared two Crabtree-positive yeasts, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and two Crabtree-negative yeasts, Kluyveromyces marxianus and Scheffersomyces stipitis, cultivated under glucose excess conditions. Combining physiological and proteome quantification with genome-scale metabolic modeling, we found that the two groups differ in energy metabolism and translation efficiency. In Crabtree-positive yeasts, the central carbon metabolism flux and proteome allocation favor a glucose utilization strategy minimizing proteome cost as proteins translation parameters, including ribosomal content and/or efficiency, are lower. Crabtree-negative yeasts, however, use a strategy of maximizing ATP yield, accompanied by higher protein translation parameters. Our analyses provide insight into the underlying reasons for the Crabtree effect, demonstrating a coupling to adaptations in both metabolism and protein translation.

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome-scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence-level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.

Fermenting Futures: an artistic view on yeast biotechnology

BioArt is a new discipline where artists employ materials and techniques of modern life sciences and create novel meanings of biology, often involving living organisms such as tissue culture, bacteria and yeasts, which may also be genetically engineered. The authors have engaged in a collaboration to develop 'Fermenting Futures', a project designed to explore the significance of yeast for early human history by enabling baking and brewing, all the way to industrial biotechnology and synthetic biology with their potential contributions to fight the climate change. Research in two of the authors' lab provides the materials and thematic lines for the artists to develop their installations. The two main pieces reflect on fermentation as a metabolic trait of baker's yeast and its enormous transformational power for human society, and on the application of synthetic biology to enable yeast to grow and produce materials from carbon dioxide. The role of BioArt to support public engagement and science dissemination is discussed, highlighting the importance of collaborations of scientists and artists on equal terms, as showcased here.

Truth in wine yeast

Evolutionary history and early association with anthropogenic environments have made Saccharomyces cerevisiae the quintessential wine yeast. This species typically dominates any spontaneous wine fermentation and, until recently, virtually all commercially available wine starters belonged to this species. The Crabtree effect, and the ability to grow under fully anaerobic conditions, contribute decisively to their dominance in this environment. But not all strains of Saccharomyces cerevisiae are equally suitable as starter cultures. In this article, we review the physiological and genetic characteristics of S. cerevisiae wine strains, as well as the biotic and abiotic factors that have shaped them through evolution. Limited genetic diversity of this group of yeasts could be a constraint to solving the new challenges of oenology. However, research in this field has for many years been providing tools to increase this diversity, from genetic engineering and classical genetic tools to the inclusion of other yeast species in the catalogues of wine yeasts. On occasion, these less conventional species may contribute to the generation of interspecific hybrids with S. cerevisiae. Thus, our knowledge about wine strains of S. cerevisiae and other wine yeasts is constantly expanding. Over the last decades, wine yeast research has been a pillar for the modernisation of oenology, and we can be confident that yeast biotechnology will keep contributing to solving any challenges, such as climate change, that we may face in the future.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Paralogs in the PKA Regulon Traveled Different Evolutionary Routes to Divergent Expression in Budding Yeast

Functional divergence of duplicate genes, or paralogs, is an important driver of novelty in evolution. In the model yeast Saccharomyces cerevisiae, there are 547 paralog gene pairs that survive from an interspecies Whole Genome Hybridization (WGH) that occurred ~100MYA. In this work, we report that ~1/6th (110) of these WGH paralogs pairs (or ohnologs) are differentially expressed with a striking pattern upon Protein Kinase A (PKA) inhibition. One member of each pair in this group has low basal expression that increases upon PKA inhibition, while the other has moderate and unchanging expression. For these genes, expression of orthologs upon PKA inhibition in the non-WGH species Kluyveromyces lactis and for PKA-related stresses in other budding yeasts shows unchanging expression, suggesting that lack of responsiveness to PKA was likely the typical ancestral phenotype prior to duplication. Promoter sequence analysis across related budding yeast species further revealed that the subsequent emergence of PKA-dependence took different evolutionary routes. In some examples, regulation by PKA and differential expression appears to have arisen following the WGH, while in others, regulation by PKA appears to have arisen in one of the two parental lineages prior to the WGH. More broadly, our results illustrate the unique opportunities presented by a WGH event for generating functional divergence by bringing together two parental lineages with separately evolved regulation into one species. We propose that functional divergence of two ohnologs can be facilitated through such regulatory divergence.

Modeling the energy metabolism in immune cells

In this review, we summarize and briefly discuss various approaches to modeling the metabolism in human immune cells, with a focus on energy metabolism. These approaches include metabolic reconstruction, elementary modes, and flux balance analysis, which are often subsumed under constraint-based modeling. Further approaches are evolutionary game theory and kinetic modeling. Many immune cells such as macrophages show the Warburg effect, meaning that glycolysis is upregulated upon activation. We outline a minimal model for explaining that effect using optimization. The effect of a confrontation with pathogen cells on immunometabolism is highlighted. Models describing the differences between M1 and M2 macrophages, ROS production in neutrophils, and tryptophan metabolism are discussed. Obstacles and future prospects are outlined.

A 37-amino acid loop in the Yarrowia lipolytica hexokinase impacts its activity and affinity and modulates gene expression

The oleaginous yeast Yarrowia lipolytica is a potent cell factory as it is able to use a wide variety of carbon sources to convert waste materials into value-added products. Nonetheless, there are still gaps in our understanding of its central carbon metabolism. Here we present an in-depth study of Y. lipolytica hexokinase (YlHxk1), a structurally unique protein. The greatest peculiarity of YlHxk1 is a 37-amino acid loop region, a structure not found in any other known hexokinases. By combining bioinformatic and experimental methods we showed that the loop in YlHxk1 is essential for activity of this protein and through that on growth of Y. lipolytica on glucose and fructose. We further proved that the loop in YlHxk1 hinders binding with trehalose 6-phosphate (T6P), a glycolysis inhibitor, as hexokinase with partial deletion of this region is 4.7-fold less sensitive to this molecule. We also found that YlHxk1 devoid of the loop causes strong repressive effect on lipase-encoding genes LIP2 and LIP8 and that the hexokinase overexpression in Y. lipolytica changes glycerol over glucose preference when cultivated in media containing both substrates.

Rapid radiation of angraecoids (Orchidaceae, Angraecinae) in tropical Africa characterised by multiple karyotypic shifts under major environmental instability

Angraecoid orchids present a remarkable diversity of chromosome numbers, which makes them a highly suitable system for exploring the impact of karyotypic changes on cladogenesis, diversification and morphological differentiation. We compiled an annotated cytotaxonomic checklist for 126 species of Angraecinae, which was utilised to reconstruct chromosomal evolution using a newly-produced, near-comprehensive phylogenetic tree that includes 245 angraecoid taxa. In tandem with this improved phylogenetic framework, using combined Bayesian, maximum likelihood and parsimony approaches on ITS-1 and five plastid markers, we propose a new cladistic nomenclature for the angraecoids, and we estimate a new timeframe for angraecoid radiation based on a secondary calibration, and calculate diversification rates using a Bayesian approach. Coincident divergence dates between clades with identical geographical distributions in the angraecoids and the pantropical orchid genus Bulbophyllum suggest that the same events may have intervened in the dispersal of these two epiphytic groups between Asia, continental Africa, Madagascar and the Neotropics. The major angraecoid lineages probably began to differentiate in the Middle Miocene, and most genera and species emerged respectively around the Late Miocene-Pliocene boundary and the Pleistocene. Ancestral state reconstruction using maximum likelihood estimation revealed an eventful karyotypic history dominated by descending dysploidy. Karyotypic shifts seem to have paralleled cladogenesis in continental tropical Africa, where approximately 90% of the species have descended from at least one inferred transition from n = 17-18 to n = 25 during the Middle Miocene Climatic Transition, followed by some clade-specific descending and ascending dysploidy from the Late Miocene to the Pleistocene. Conversely, detected polyploidy is restricted to a few species lineages mostly originating during the Pleistocene. No increases in net diversification could be related to chromosome number changes, and the apparent net diversification was found to be highest in Madagascar, where karyotypic stasis predominates. Finally, shifts in chromosome number appear to have paralleled the evolution of rostellum structure, leaflessness, and conspicuous changes in floral colour.

Mitochondrial copper and phosphate transporter specificity was defined early in the evolution of eukaryotes

The mitochondrial carrier family protein SLC25A3 transports both copper and phosphate in mammals, yet in Saccharomyces cerevisiae the transport of these substrates is partitioned across two paralogs: PIC2 and MIR1. To understand the ancestral state of copper and phosphate transport in mitochondria, we explored the evolutionary relationships of PIC2 and MIR1 orthologs across the eukaryotic tree of life. Phylogenetic analyses revealed that PIC2-like and MIR1-like orthologs are present in all major eukaryotic supergroups, indicating an ancient gene duplication created these paralogs. To link this phylogenetic signal to protein function, we used structural modeling and site-directed mutagenesis to identify residues involved in copper and phosphate transport. Based on these analyses, we generated an L175A variant of mouse SLC25A3 that retains the ability to transport copper but not phosphate. This work highlights the utility of using an evolutionary framework to uncover amino acids involved in substrate recognition by mitochondrial carrier family proteins.

Mitochondrial protein phosphorylation in yeast revisited

Protein phosphorylation is one of the best-known post-translational modifications occurring in all domains of life. In eukaryotes, protein phosphorylation affects all cellular compartments including mitochondria. High-throughput techniques of mass spectrometry combined with cell fractionation and biochemical methods yielded thousands of phospho-sites on hundreds of mitochondrial proteins. We have compiled the information on mitochondrial protein kinases and phosphatases and their substrates in Saccharomyces cerevisiae and provide the current state-of-the-art overview of mitochondrial protein phosphorylation in this model eukaryote. Using several examples, we describe emerging features of the yeast mitochondrial phosphoproteome and present challenges lying ahead in this exciting field.

Parallel Concerted Evolution of Ribosomal Protein Genes in Fungi and Its Adaptive Significance

Ribosomal protein (RP) genes encode structural components of ribosomes, the cellular machinery for protein synthesis. A single functional copy has been maintained in most of 78–80 RP families in animals due to evolutionary constraints imposed by gene dosage balance. Some fungal species have maintained duplicate copies in most RP families. The mechanisms by which the RP genes were duplicated and maintained and their functional significance are poorly understood. To address these questions, we identified all RP genes from 295 fungi and inferred the timing and nature of gene duplication events for all RP families. We found that massive duplications of RP genes have independently occurred by different mechanisms in three distantly related lineages: budding yeasts, fission yeasts, and Mucoromycota. The RP gene duplicates in budding yeasts and Mucoromycota were mainly created by whole genome duplication events. However, duplicate RP genes in fission yeasts were likely generated by retroposition, which is unexpected considering their dosage sensitivity. The sequences of most RP paralogs have been homogenized by repeated gene conversion in each species, demonstrating parallel concerted evolution, which might have facilitated the retention of their duplicates. Transcriptomic data suggest that the duplication and retention of RP genes increased their transcript abundance. Physiological data indicate that increased ribosome biogenesis allowed these organisms to rapidly consume sugars through fermentation while maintaining high growth rates, providing selective advantages to these species in sugar-rich environments.

Performance, genomic rearrangements, and signatures of adaptive evolution: Lessons from fermentative yeasts

The capacity of some yeasts to extract energy from single sugars, generating CO2 and ethanol (=fermentation), even in the presence of oxygen, is known as the Crabtree effect. This phenomenon represents an important adaptation as it allowed the utilization of the ecological niche given by modern fruits, an abundant source of food that emerged in the terrestrial environment in the Cretaceous. However, identifying the evolutionary events that triggered fermentative capacity in Crabtree-positive species is challenging, as microorganisms do not leave fossil evidence. Thus, key innovations should be inferred based only on traits measured under culture conditions. Here, we reanalyzed data from a common garden experiment where several proxies of fermentative capacity were recorded in Crabtree-positive and Crabtree-negative species, representing yeast phylogenetic diversity. In particular, we applied the "lasso-OU" algorithm which detects points of adaptive shifts, using traits that are proxies of fermentative performance. We tested whether multiple events or a single event explains the actual fermentative capacity of yeasts. According to the lasso-OU procedure, evolutionary changes in the three proxies of fermentative capacity that we considered (i.e., glycerol production, ethanol yield, and respiratory quotient) are consistent with a single evolutionary episode (a whole-genomic duplication, WGD), instead of a series of small genomic rearrangements. Thus, the WGD appears as the key event behind the diversification of fermentative yeasts, which by increasing gene dosage, and maximized their capacity of energy extraction for exploiting the new ecological niche provided by single sugars.

The lasting after-effects of an ancient polyploidy on the genomes of teleosts

The ancestor of most teleost fishes underwent a whole-genome duplication event three hundred million years ago. Despite its antiquity, the effects of this event are evident both in the structure of teleost genomes and in how the surviving duplicated genes still operate to drive form and function. I inferred a set of shared syntenic regions that survive from the teleost genome duplication (TGD) using eight teleost genomes and the outgroup gar genome (which lacks the TGD). I then phylogenetically modeled the TGD's resolution via shared and independent gene losses and applied a new simulation-based statistical test for the presence of bias toward the preservation of genes from one parental subgenome. On the basis of that test, I argue that the TGD was likely an allopolyploidy. I find that duplicate genes surviving from this duplication in zebrafish are less likely to function in early embryo development than are genes that have returned to single copy at some point in this species' history. The tissues these ohnologs are expressed in, as well as their biological functions, lend support to recent suggestions that the TGD was the source of a morphological innovation in the structure of the teleost retina. Surviving duplicates also appear less likely to be essential than singletons, despite the fact that their single-copy orthologs in mouse are no less essential than other genes.

Advances in yeast alcoholic fermentations for the production of bioethanol, beer and wine

Yeasts have a long-standing relationship with humankind that has widened in recent years to encompass production of diverse foods, beverages, fuels and medicines. Here, key advances in the field of yeast fermentation applied to alcohol production, which represents the predominant product of industrial biotechnology, will be presented. More specifically, we have selected industries focused in producing bioethanol, beer and wine. In these bioprocesses, yeasts from the genus Saccharomyces are still the main players, with Saccharomyces cerevisiae recognized as the preeminent industrial ethanologen. However, the growing demand for new products has opened the door to diverse yeasts, including non-Saccharomyces strains. Furthermore, the development of synthetic media that successfully simulate industrial fermentation medium will be discussed along with a general overview of yeast fermentation modeling.

Porcine Fetal-Derived Fibroblasts Alter Gene Expression and Mitochondria to Compensate for Hypoxic Stress During Culture

The Warburg effect is characterized by decreased mitochondrial oxidative phosphorylation and increased glycolytic flux in adequate oxygen. The preimplantation embryo has been described to have characteristics of the Warburg effect, including similar changes in gene expression and mitochondria, which are more rudimentary in appearance. We hypothesized hypoxia would facilitate anaerobic glycolysis in fibroblasts thereby promoting gene expression and media metabolite production reflecting the Warburg effect hallmarks in early embryos. Additionally, we speculated that hypoxia would induce a rudimentary small mitochondrial phenotype observed in several cell types evidenced to demonstrate the Warburg effect. While many have examined the role hypoxia plays in pathological conditions, few studies have investigated changes in primary cells which could be used in somatic cell nuclear transfer. We found that cells grown in 1.25% O₂ had normal cell viability and more, but smaller mitochondria. Several hypoxia-inducible genes were identified, including seven genes for glycolytic enzymes. In conditioned media from hypoxic cells, the quantities of gluconolactone, cytosine, and uric acid were decreased indicating higher consumption than control cells. These results indicate that fibroblasts alter gene expression and mitochondria to compensate for hypoxic stress and maintain viability. Furthermore, the metabolic changes observed, making them more similar to preimplantation embryos, could be facilitating nuclear reprogramming making these cells more amendable to future use in somatic cell nuclear transfer.

Functional link between mitochondria and Rnr3, the minor catalytic subunit of yeast ribonucleotide reductase

Ribonucleotide reductase (RNR) is an essential holoenzyme required for de novo synthesis of dNTPs. The Saccharomyces cerevisiae genome encodes for two catalytic subunits, Rnr1 and Rnr3. While Rnr1 is required for DNA replication and DNA damage repair, the function(s) of Rnr3 is unknown. Here, we show that carbon source, an essential nutrient, impacts Rnr1 and Rnr3 abundance: Non-fermentable carbon sources or limiting concentrations of glucose down regulate Rnr1 and induce Rnr3 expression. Oppositely, abundant glucose induces Rnr1 expression and down regulates Rnr3. The carbon source dependent regulation of Rnr3 is mediated by Mec1, the budding yeast ATM/ATR checkpoint response kinase. Unexpectedly, this regulation is independent of all currently known components of the Mec1 DNA damage response network, including Rad53, Dun1, and Tel1, implicating a novel Mec1 signalling axis. rnr3Δ leads to growth defects under respiratory conditions and rescues temperature sensitivity conferred by the absence of Tom6, a component of the mitochondrial TOM (translocase of outer membrane) complex responsible for mitochondrial protein import. Together, these results unveil involvement of Rnr3 in mitochondrial functions and Mec1 in mediating the carbon source dependent regulation of Rnr3.

Harmonizing semantic annotations for computational models in biology

Life science researchers use computational models to articulate and test hypotheses about the behavior of biological systems. Semantic annotation is a critical component for enhancing the interoperability and reusability of such models as well as for the integration of the data needed for model parameterization and validation. Encoded as machine-readable links to knowledge resource terms, semantic annotations describe the computational or biological meaning of what models and data represent. These annotations help researchers find and repurpose models, accelerate model composition and enable knowledge integration across model repositories and experimental data stores. However, realizing the potential benefits of semantic annotation requires the development of model annotation standards that adhere to a community-based annotation protocol. Without such standards, tool developers must account for a variety of annotation formats and approaches, a situation that can become prohibitively cumbersome and which can defeat the purpose of linking model elements to controlled knowledge resource terms. Currently, no consensus protocol for semantic annotation exists among the larger biological modeling community. Here, we report on the landscape of current annotation practices among the COmputational Modeling in BIology NEtwork community and provide a set of recommendations for building a consensus approach to semantic annotation.

Dominance of wine Saccharomyces cerevisiae strains over S. kudriavzevii in industrial fermentation competitions is related to an acceleration of nutrient uptake and utilization

Grape must is a sugar-rich habitat for a complex microbiota which is replaced by Saccharomyces cerevisiae strains during the first fermentation stages. Interest on yeast competitive interactions has recently been propelled due to the use of alternative yeasts in the wine industry to respond to new market demands. The main issue resides in the persistence of these yeasts due to the specific competitive activity of S. cerevisiae. To gather deeper knowledge of the molecular mechanisms involved, we performed a comparative transcriptomic analysis during fermentation carried out by a wine S. cerevisiae strain and a strain representative of the cryophilic S. kudriavzevii, which exhibits high genetic and physiological similarities to S. cerevisiae, but also differences of biotechnological interest. In this study, we report that transcriptomic response to the presence of a competitor is stronger in S. cerevisiae than in S. kudriavzevii. Our results demonstrate that a wine S. cerevisiae industrial strain accelerates nutrient uptake and utilization to outcompete the co-inoculated yeast, and that this process requires cell-to-cell contact to occur. Finally, we propose that this competitive phenotype evolved recently, during the adaptation of S. cerevisiae to man-manipulated fermentative environments, since a non-wine S. cerevisiae strain, isolated from a North American oak, showed a remarkable low response to competition.

Minor Isozymes Tailor Yeast Metabolism to Carbon Availability

Gene duplication is one of the main evolutionary paths to new protein function. Typically, duplicated genes either accumulate mutations and degrade into pseudogenes or are retained and diverge in function. Some duplicated genes, however, show long-term persistence without apparently acquiring new function. An important class of isozymes consists of those that catalyze the same reaction in the same compartment, where knockout of one isozyme causes no known functional defect. Here we present an approach to assigning specific functional roles to seemingly redundant isozymes. First, gene expression data are analyzed computationally to identify conditions under which isozyme expression diverges. Then, knockouts are compared under those conditions. This approach revealed that the expression of many yeast isozymes diverges in response to carbon availability and that carbon source manipulations can induce fitness phenotypes for seemingly redundant isozymes. A driver of these fitness phenotypes is differential allosteric enzyme regulation, indicating isozyme divergence to achieve more-optimal control of metabolism.

Isozymes are enzymes that differ in sequence but catalyze the same chemical reactions. Despite their apparent redundancy, isozymes are often retained over evolutionary time, suggesting that they contribute to fitness. We developed an unsupervised computational method for identifying environmental conditions under which isozymes are likely to make fitness contributions. This method analyzes published gene expression data to find specific experimental perturbations that induce differential isozyme expression. In yeast, we found that isozymes are strongly enriched in the pathways of central carbon metabolism and that many isozyme pairs show anticorrelated expression during the respirofermentative shift. Building on these observations, we assigned function to two minor central carbon isozymes, aconitase 2 (ACO2) and pyruvate kinase 2 (PYK2). ACO2 is expressed during fermentation and proves advantageous when glucose is limiting. PYK2 is expressed during respiration and proves advantageous for growth on three-carbon substrates. PYK2’s deletion can be rescued by expressing the major pyruvate kinase only if that enzyme carries mutations mirroring PYK2’s allosteric regulation. Thus, central carbon isozymes help to optimize allosteric metabolic regulation under a broad range of potential nutrient conditions while requiring only a small number of transcriptional states.

IMPORTANCE Gene duplication is one of the main evolutionary paths to new protein function. Typically, duplicated genes either accumulate mutations and degrade into pseudogenes or are retained and diverge in function. Some duplicated genes, however, show long-term persistence without apparently acquiring new function. An important class of isozymes consists of those that catalyze the same reaction in the same compartment, where knockout of one isozyme causes no known functional defect. Here we present an approach to assigning specific functional roles to seemingly redundant isozymes. First, gene expression data are analyzed computationally to identify conditions under which isozyme expression diverges. Then, knockouts are compared under those conditions. This approach revealed that the expression of many yeast isozymes diverges in response to carbon availability and that carbon source manipulations can induce fitness phenotypes for seemingly redundant isozymes. A driver of these fitness phenotypes is differential allosteric enzyme regulation, indicating isozyme divergence to achieve more-optimal control of metabolism.

Candida glabrata: A Lot More Than Meets the Eye

Candida glabrata is an opportunistic human fungal pathogen that causes superficial mucosal and life-threatening bloodstream infections in individuals with a compromised immune system. Evolutionarily, it is closer to the non-pathogenic yeast Saccharomyces cerevisiae than to the most prevalent Candida bloodstream pathogen, C. albicans. C. glabrata is a haploid budding yeast that predominantly reproduces clonally. In this review, we summarize interactions of C. glabrata with the host immune, epithelial and endothelial cells, and the ingenious strategies it deploys to acquire iron and phosphate from the external environment. We outline various attributes including cell surface-associated adhesins and aspartyl proteases, biofilm formation and stress response mechanisms, that contribute to the virulence of C. glabrata. We further discuss how, C. glabrata, despite lacking morphological switching and secreted proteolytic activity, is able to disarm macrophage, dampen the host inflammatory immune response and replicate intracellularly.

The Regulation of Cbf1 by PAS Kinase Is a Pivotal Control Point for Lipogenesis vs. Respiration in Saccharomyces cerevisiae

PAS kinase 1 (Psk1) is a key regulator of respiration in Saccharomyces cerevisiae. Herein the molecular mechanisms of this regulation are explored through the characterization of its substrate, Centromere binding factor 1 (Cbf1). CBF1-deficient yeast displayed a significant decrease in cellular respiration, while PAS kinase-deficient yeast, or yeast harboring a Cbf1 phosphosite mutant (T211A) displayed a significant increase. Transmission electron micrographs showed an increased number of mitochondria in PAS kinase-deficient yeast consistent with the increase in respiration. Although the CBF1-deficient yeast did not appear to have an altered number of mitochondria, a mitochondrial proteomics study revealed significant differences in the mitochondrial composition of CBF1-deficient yeast including altered Atp3 levels, a subunit of the mitochondrial F₁-ATP synthase complex. Both beta-galactosidase reporter assays and western blot analysis confirmed direct transcriptional control of ATP3 by Cbf1. In addition, we confirmed the regulation of yeast lipid genes LAC1 and LAG1 by Cbf1. The human homolog of Cbf1, Upstream transcription factor 1 (USF1), is also known to be involved in lipid biogenesis. Herein, we provide the first evidence for a role of USF1 in respiration since it appeared to complement Cbf1in vivo as determined by respiration phenotypes. In addition, we confirmed USF1 as a substrate of human PAS kinase (hPASK) in vitro. Combined, our data supports a model in which Cbf1/USF1 functions to partition glucose toward respiration and away from lipid biogenesis, while PAS kinase inhibits respiration in part through the inhibition of Cbf1/USF1.

A single Gal4-like transcription factor activates the Crabtree effect in Komagataella phaffii

The Crabtree phenotype defines whether a yeast can perform simultaneous respiration and fermentation under aerobic conditions at high growth rates. It provides Crabtree positive yeasts an evolutionary advantage of consuming glucose faster and producing ethanol to outcompete other microorganisms in sugar rich environments. While a number of genetic events are associated with the emergence of the Crabtree effect, its evolution remains unresolved. Here we show that overexpression of a single Gal4-like transcription factor is sufficient to convert Crabtree-negative Komagataella phaffii (Pichia pastoris) into a Crabtree positive yeast. Upregulation of the glycolytic genes and a significant increase in glucose uptake rate due to the overexpression of the Gal4-like transcription factor leads to an overflow metabolism, triggering both short-term and long-term Crabtree phenotypes. This indicates that a single genetic perturbation leading to overexpression of one gene may have been sufficient as the first molecular event towards respiro-fermentative metabolism in the course of yeast evolution.

Aerobic ethanol production, a phenomenon referred as Crabtree effect, allows yeast to outcompete other microorganisms in sugar rich environments. Here, the authors show that overexpression of a Gal4-like transcription factor can transform Komagataella phaffii from Crabtree effect negative to positive.

The Genetic Makeup and Expression of the Glycolytic and Fermentative Pathways Are Highly Conserved Within the Saccharomyces Genus

The ability of the yeast Saccharomyces cerevisiae to convert glucose, even in the presence of oxygen, via glycolysis and the fermentative pathway to ethanol has played an important role in its domestication. Despite the extensive knowledge on these pathways in S. cerevisiae, relatively little is known about their genetic makeup in other industrially relevant Saccharomyces yeast species. In this study we explore the diversity of the glycolytic and fermentative pathways within the Saccharomyces genus using S. cerevisiae, Saccharomyces kudriavzevii, and Saccharomyces eubayanus as paradigms. Sequencing data revealed a highly conserved genetic makeup of the glycolytic and fermentative pathways in the three species in terms of number of paralogous genes. Although promoter regions were less conserved between the three species as compared to coding sequences, binding sites for Rap1, Gcr1 and Abf1, main transcriptional regulators of glycolytic and fermentative genes, were highly conserved. Transcriptome profiling of these three strains grown in aerobic batch cultivation in chemically defined medium with glucose as carbon source, revealed a remarkably similar expression of the glycolytic and fermentative genes across species, and the conserved classification of genes into major and minor paralogs. Furthermore, transplantation of the promoters of major paralogs of S. kudriavzevii and S. eubayanus into S. cerevisiae demonstrated not only the transferability of these promoters, but also the similarity of their strength and response to various environmental stimuli. The relatively low homology of S. kudriavzevii and S. eubayanus promoters to their S. cerevisiae relatives makes them very attractive alternatives for strain construction in S. cerevisiae, thereby expanding the S. cerevisiae molecular toolbox.

Integrating Networks, Phylogenomics, and Population Genomics for the Study of Polyploidy

Duplication events are regarded as sources of evolutionary novelty, but our understanding of general trends for the long-term trajectory of additional genomic material is still lacking. Organisms with a history of whole genome duplication (WGD) offer a unique opportunity to study potential trends in the context of gene retention and/or loss, gene and network dosage, and changes in gene expression. In this review, we discuss the prevalence of polyploidy across the tree of life, followed by an overview of studies investigating genome evolution and gene expression. We then provide an overview of methods in network biology, phylogenomics, and population genomics that are critical for advancing our understanding of evolution post-WGD, highlighting the need for models that can accommodate polyploids. Finally, we close with a brief note on the importance of random processes in the evolution of polyploids with respect to neutral versus selective forces, ancestral polymorphisms, and the formation of autopolyploids versus allopolyploids.

Expression properties exhibit correlated patterns with the fate of duplicated genes, their divergence, and transcriptional plasticity in Saccharomycotina

Gene duplication is an important source of novelties and genome complexity. What genes are preserved as duplicated through long evolutionary times can shape the evolution of innovations. Identifying factors that influence gene duplicability is therefore an important aim in evolutionary biology. Here, we show that in the yeast Saccharomyces cerevisiae the levels of gene expression correlate with gene duplicability, its divergence, and transcriptional plasticity. Genes that were highly expressed before duplication are more likely to be preserved as duplicates for longer evolutionary times and wider phylogenetic ranges than genes that were lowly expressed. Duplicates with higher expression levels exhibit greater divergence between their gene copies. Duplicates that exhibit higher expression divergence are those enriched for TATA-containing promoters. These duplicates also show transcriptional plasticity, which seems to be involved in the origin of adaptations to environmental stresses in yeast. While the expression properties of genes strongly affect their duplicability, divergence and transcriptional plasticity are enhanced after gene duplication. We conclude that highly expressed genes are more likely to be preserved as duplicates due to their promoter architectures, their greater tolerance to expression noise, and their ability to reduce the noise-plasticity conflict.

Heterosis Is a Systemic Property Emerging From Non-linear Genotype-Phenotype Relationships: Evidence From in Vitro Genetics and Computer Simulations

Heterosis, the superiority of hybrids over their parents for quantitative traits, represents a crucial issue in plant and animal breeding as well as evolutionary biology. Heterosis has given rise to countless genetic, genomic and molecular studies, but has rarely been investigated from the point of view of systems biology. We hypothesized that heterosis is an emergent property of living systems resulting from frequent concave relationships between genotypic variables and phenotypes, or between different phenotypic levels. We chose the enzyme-flux relationship as a model of the concave genotype-phenotype (GP) relationship, and showed that heterosis can be easily created in the laboratory. First, we reconstituted in vitro the upper part of glycolysis. We simulated genetic variability of enzyme activity by varying enzyme concentrations in test tubes. Mixing the content of "parental" tubes resulted in "hybrids," whose fluxes were compared to the parental fluxes. Frequent heterotic fluxes were observed, under conditions that were determined analytically and confirmed by computer simulation. Second, to test this model in a more realistic situation, we modeled the glycolysis/fermentation network in yeast by considering one input flux, glucose, and two output fluxes, glycerol and acetaldehyde. We simulated genetic variability by randomly drawing parental enzyme concentrations under various conditions, and computed the parental and hybrid fluxes using a system of differential equations. Again we found that a majority of hybrids exhibited positive heterosis for metabolic fluxes. Cases of negative heterosis were due to local convexity between certain enzyme concentrations and fluxes. In both approaches, heterosis was maximized when the parents were phenotypically close and when the distributions of parental enzyme concentrations were contrasted and constrained. These conclusions are not restricted to metabolic systems: they only depend on the concavity of the GP relationship, which is commonly observed at various levels of the phenotypic hierarchy, and could account for the pervasiveness of heterosis.

Natural fruits, flowers, honey, and honeybees harbor Helicobacter pylori-positive yeasts

BACKGROUND

For controlling Helicobacter pylori infection in humans, its environmental reservoir should be determined. In this study, yeast isolates from an isolated village in Iran were studied for the intracellular occurrence of H. pylori.

MATERIALS AND METHODS

In this study, yeasts were isolated from 29 samples, including oral swabs from villagers (n = 7), flowers and fruits (n = 6), honey and honeybees (n = 12) and miscellaneous samples (4). Yeasts were classified into 12 RFLP groups and identified by amplification of 26S rDNA and sequencing. DNA extracted from the yeast cells was examined for the presence of H. pylori using PCR.

RESULTS

Of the 29 yeasts, 27 were members of different genera of Ascomycete. H. pylori was detected in 5 of 9 Candida (55.5%), 4 of 5 Komagataella (80%), 3 of 4 Pichia (100%), 2 of 2 Cytobasidia (100%), 2 of 2 Hansenia (100%), 1 of 1 Meyerozyma (100%) and 2 of 3 not sequenced (66.6%) yeasts. Distribution of 19 of 29 (65.5%) H. pylori-positive yeasts within 4 groups was as follows: 1 of 7(14.3%) in oral swabs, 5 of 6 (83.3%) in flowers and fruits, 10 of 12 (83.3%) in honey and the bee group and 3 of 4 (75%) in miscellaneous.

CONCLUSIONS

Different genera of osmotolerant yeasts from flowers, fruits, honey, and honeybees contained H. pylori in their vacuole. High frequency of H. pylori-positive yeasts in these samples might be related to their high sugar content. Insects such as honeybees that facilitate transfer and easy access of these yeasts to nectars serve as the main reservoirs of these yeasts, playing an important role in their protection and dispersal. Accordingly, H. pylori inside these yeasts can be carried by honeybees to different sugar- and nutrient-rich environments. Sugar-rich environments and honeybees play an important role in distribution of H. pylori-positive yeasts in nature.

Preferential retention of genes from one parental genome after polyploidy illustrates the nature and scope of the genomic conflicts induced by hybridization

Polyploidy is increasingly seen as a driver of both evolutionary innovation and ecological success. One source of polyploid organisms' successes may be their origins in the merging and mixing of genomes from two different species (e.g., allopolyploidy). Using POInT (the Polyploid Orthology Inference Tool), we model the resolution of three allopolyploidy events, one from the bakers' yeast (Saccharomyces cerevisiae), one from the thale cress (Arabidopsis thaliana) and one from grasses including Sorghum bicolor. Analyzing a total of 21 genomes, we assign to every gene a probability for having come from each parental subgenome (i.e., derived from the diploid progenitor species), yielding orthologous segments across all genomes. Our model detects statistically robust evidence for the existence of biased fractionation in all three lineages, whereby genes from one of the two subgenomes were more likely to be lost than those from the other subgenome. We further find that a driver of this pattern of biased losses is the co-retention of genes from the same parental genome that share functional interactions. The pattern of biased fractionation after the Arabidopsis and grass allopolyploid events was surprisingly constant in time, with the same parental genome favored throughout the lineages' history. In strong contrast, the yeast allopolyploid event shows evidence of biased fractionation only immediately after the event, with balanced gene losses more recently. The rapid loss of functionally associated genes from a single subgenome is difficult to reconcile with the action of genetic drift and suggests that selection may favor the removal of specific duplicates. Coupled to the evidence for continuing, functionally-associated biased fractionation after the A. thaliana At-α event, we suggest that, after allopolyploidy, there are functional conflicts between interacting genes encoded in different subgenomes that are ultimately resolved through preferential duplicate loss.

Reciprocally Retained Genes in the Angiosperm Lineage Show the Hallmarks of Dosage Balance Sensitivity

In several organisms, particular functional categories of genes, such as regulatory and complex-forming genes, are preferentially retained after whole-genome multiplications but rarely duplicate through small-scale duplication, a pattern referred to as reciprocal retention. This peculiar duplication behavior is hypothesized to stem from constraints on the dosage balance between the genes concerned and their interaction context. However, the evidence for a relationship between reciprocal retention and dosage balance sensitivity remains fragmentary. Here, we identified which gene families are most strongly reciprocally retained in the angiosperm lineage and studied their functional and evolutionary characteristics. Reciprocally retained gene families exhibit stronger sequence divergence constraints and lower rates of functional and expression divergence than other gene families, suggesting that dosage balance sensitivity is a general characteristic of reciprocally retained genes. Gene families functioning in regulatory and signaling processes are much more strongly represented at the top of the reciprocal retention ranking than those functioning in multiprotein complexes, suggesting that regulatory imbalances may lead to stronger fitness effects than classical stoichiometric protein complex imbalances. Finally, reciprocally retained duplicates are often subject to dosage balance constraints for prolonged evolutionary times, which may have repercussions for the ease with which genome multiplications can engender evolutionary innovation.

The evolution, structure and function of the ray finned fish (Actinopterygii) glucocorticoid receptors

Basal ray-finned fish (Actinopterygii) possess a single glucocorticoid receptor (GR) and when compared to the lobe-finned vertebrate (Sarcopterygii) GR possess nine additional amino acids between the zinc-finger of the DNA binding domain. A whole genome duplication event which occurred between 320 and 350MYA in the teleost lineage following the split from the basal ray-finned fish resulted in 2 GRs: one GR group, GR1, has retained the 9 amino acids insert whereas the other group, GR2, has not. The exception to this is the zebrafish, that have lost one of the GRs, but they do possess 2 GRs with a splice variant that lacks the C-terminal portion of the GR to form GRβ which acts as a dominant-repressor of the wildtype GR. Another splice variant sees the basal ray-finned GR and teleost GR1 without the 9 amino acids insert. The molecular basis for GRs retention is beginning to be unravelled. In Pantadon buchholzi, rainbow trout, carp, marine and Japanese medaka GR2 is more sensitive to glucocorticoids (GC), thus potentially playing a more significant role in regulating gene expression at basal circulatory GC concentrations. However, this division in GC sensitivity is not seen in other species. The few studies to evaluate the significance of the 9 amino acid insert have shown that it affect maximal transactivational activity the extent to which is dependent on the number of glucocorticoid response elements (GREs) present in the reporter plasmid. The retention of these GRs would suggest there was an evolutionary advantage, which saw the development of a complex regulatory process to mediate the actions of the glucocorticoids.

Recent advances in understanding the roles of whole genome duplications in evolution

Ancient whole-genome duplications (WGDs)- paleopolyploidy events-are key to solving Darwin's 'abominable mystery' of how flowering plants evolved and radiated into a rich variety of species. The vertebrates also emerged from their invertebrate ancestors via two WGDs, and genomes of diverse gymnosperm trees, unicellular eukaryotes, invertebrates, fishes, amphibians and even a rodent carry evidence of lineage-specific WGDs. Modern polyploidy is common in eukaryotes, and it can be induced, enabling mechanisms and short-term cost-benefit assessments of polyploidy to be studied experimentally. However, the ancient WGDs can be reconstructed only by comparative genomics: these studies are difficult because the DNA duplicates have been through tens or hundreds of millions of years of gene losses, mutations, and chromosomal rearrangements that culminate in resolution of the polyploid genomes back into diploid ones (rediploidisation). Intriguing asymmetries in patterns of post-WGD gene loss and retention between duplicated sets of chromosomes have been discovered recently, and elaborations of signal transduction systems are lasting legacies from several WGDs. The data imply that simpler signalling pathways in the pre-WGD ancestors were converted via WGDs into multi-stranded parallelised networks. Genetic and biochemical studies in plants, yeasts and vertebrates suggest a paradigm in which different combinations of sister paralogues in the post-WGD regulatory networks are co-regulated under different conditions. In principle, such networks can respond to a wide array of environmental, sensory and hormonal stimuli and integrate them to generate phenotypic variety in cell types and behaviours. Patterns are also being discerned in how the post-WGD signalling networks are reconfigured in human cancers and neurological conditions. It is fascinating to unpick how ancient genomic events impact on complexity, variety and disease in modern life.