gene duplication plays a key role in evolution by providing dna templates for evolutionary innovation , while preventing interference with the cellular function of the original genes ( 13 ) . immediately after gene duplication , the resulting paralog pairs are usually identical and therefore functionally redundant . unless duplication confers a selective advantage , either via gene dosage effects or via mutational acquisition of modified or new functions , duplicated genes will eventually be pseudogenized and/or lost from the genome ( 2 , 3 ) . gene duplication played a key role in the evolutionary history of the yeast saccharomyces cerevisiae , an intensively investigated model in eukaryotic evolutionary biology . a whole - genome duplication event ( wgd ) in an ancestor of s. cerevisiae , ca . despite the long time interval that separates current s. cerevisiae strains from the wgd , many surviving paralog pairs still exhibit a substantial degree of functional redundancy , as indicated by neutral or weak phenotypes of single - paralog deletion mutants ( 69 ) . however , the selective pressures that caused the long - term retention of these functionally redundant paralogs remain elusive ( 1 , 10 ) . the embden - meyerhof - parnas ( emp ) pathway , the main route for oxidation of glucose to pyruvate in all eukaryotes and in many other organisms , is among the most slowly evolving metabolic pathways ( 11 , 12 ) . in all taxa , paralog families are found for the structural genes that encode the emp enzymes ( 1214 ) . under conditions of oxygen limitation and/or sugar excess , s. cerevisiae couples the emp pathway to the fermentative production of ethanol via pyruvate decarboxylase and nad - dependent alcohol dehydrogenase ( 15 , 16 ) . in this article , we use the term  glycolysis  to indicate the set of 12 enzyme - catalyzed reactions in yeast that convert intracellular glucose to ethanol . in s. cerevisiae , no fewer than 8 of the 12 enzyme reactions in glycolysis are represented by multiple paralogous genes ( fig . 1 ) . this incidence represents a significant overrepresentation ( p = 1.9  10 , based on hypergeometric distribution analysis ) relative to the ca . the wgd event and resulting duplication of genes involved in central metabolism have been implicated in the appearance of several key physiological characteristics of s. cerevisiae . in particular , the duplication of glycolytic genes has been proposed to have contributed to the strong tendency of s. cerevisiae to produce ethanol under aerobic conditions ( crabtree effect ) and its high glycolytic capacity ( 1820 ) . however , the impact of reducing the number of glycolytic paralogs on these and other physiological characteristics of s. cerevisiae has not been systematically explored . in all paralogous gene sets in yeast glycolysis , with the notable exception of the phosphofructokinase gene , gene expression and gene deletion studies support the definition of a single major paralog and one to four minor paralogs ( fig . all paralogs have retained their original catalytic function , although their context - dependent expression profiles differ ( fig . deletion of minor paralogs for individual glycolytic enzymes has minor effects on enzyme activities in cell extracts and on the specific growth rate under standard laboratory conditions ( usually shake flask cultivation on yeast extract , peptone , and glucose ) ( 2125 ) ( fig . several hypotheses have been forwarded to explain the neutral effect of paralog deletion ( 6 , 2631 ) . these include experimental limitations , such as the poor sensitivity of fitness screens ( 32 ) and the narrow range of cultivation condition tested ( generally complex medium ) ( 33 ) . additionally , analysis of deletion mutants in which paralogs that encode a single glycolytic enzyme are inactivated can not reveal synergistic effects of the minor paralogs of different glycolytic enzymes . such synergistic effects might , for example , arise from the well - documented phenomenon of distribution of metabolic control over multiple enzymes in metabolic pathways ( 34 , 35 ) or from other regulatory or catalytic interactions . the following eight enzymes in yeast glycolysis are encoded by parologous genes : hexokinase / glucokinase ( hxk ) , phosphofructokinase ( pfk ) , glyceraldehyde-3-phosphate dehydrogenase ( gapdh ) , phosphoglycerate mutase ( gpm ) , enolase ( eno ) , pyruvate kinase ( pyk ) , pyruvate decarboxylase ( pdc ) , and alcohol dehydrogenase ( adh ) . percentages of gene and protein similarity ( gesy and prsy , respectively ) between paralogs , the type of duplication event ( duev ) , whole - genome duplication ( wgd ) or other small scale duplication ( ssd ) , and the proposed fate ( pseudogenization [ p ] , subfunctionalization [ s ] , or neofunctionalization [ n ] ) of the different paralogs are indicated in the tables in each panel . separation between wgd and small - scale duplications pre - wgd ( ssd * ) and post - wgd ( ssd ) was based on the information in the yeast gene order browser ( ygob ; http://ygob.ucd.ie/ ) ( 119 ) . scatter plots show the comparisons between expression levels of the different glycolytic paralogs in s. cerevisiae measured under 170 different conditions ( 75 ) ( see table s6 in the supplemental material ) . except for pfk , where both paralogs pfk1 and pfk2 were considered paralogs with equivalent contributions , the expression levels of the major glycolytic paralogs ( labeled in red ) are indicated by a red line . bar graphs display in vitro enzyme activities of glycolytic enzymes in mutants carrying individual and multiple glycolytic gene knockouts . values are presented as percentages of in vitro enzyme activities measured in reference strains ( rs [ denoted by an orange dotted line ] ) as reported in the literature ( 2125 ) . insight into the importance , under laboratory conditions , of the minor glycolytic paralogs is essential for understanding an apparent genetic redundancy in a key ubiquitous metabolic pathway in an important model organism and industrial platform . in addition , analysis of the extent of genetic redundancy in central metabolism is highly relevant for the complete redesign and construction of entirely synthetic yeast genomes , as pursued in the synthetic yeast 2.0 initiative ( 36 ) . moreover , if complexity in yeast glycolysis can be significantly reduced by elimination of  redundant  isoenzymes , this could eliminate uncertainties and thereby facilitate the formulation and validation of mathematical models that describe the kinetics of this key metabolic pathway ( 37 , 38 ) . the goals of the present study are to experimentally explore genetic redundancy in yeast glycolysis by cumulative deletion of minor paralogs and to provide a new experimental platform for fundamental yeast research by constructing a yeast strain with a functional  minimal glycolysis  ( mg ) . to this end , we deleted 13 minor paralogs , leaving only the 14 major paralogs for the s. cerevisiae glycolytic pathway . the cumulative impact of deletion of all minor paralogs was investigated by two complementary approaches . a first , quantitative analysis focused on the impact on glycolytic flux under a number of controlled cultivation conditions that in wild - type strains result in different glycolytic fluxes . these quantitative growth studies were combined with transcriptome , enzyme activity , and intracellular metabolite assays to capture potential small phenotypic effects . a second , semiquantitative characterization explored the phenotype of the  minimal glycolysis  ( mg ) strain under a wide array of experimental conditions to identify potential context - dependent phenotypes . plasmid propagation and isolation were performed with chemically competent escherichia coli dh5 ( z - competent transformation kit ; zymo research , orange , ca ) cultivated in lysogeny broth ( lb ) medium ( 39 , 40 ) supplemented with 100 mg liter ampicillin ( lbamp ) when required . all yeast strains are derived from the cen.pk family ( 4143 ) and are listed in table s3 in the supplemental material . all strains were stored at 80c in 1-ml aliquots of 30% glycerol and the appropriate medium . cen.pk102-12a was selected as the parental strain for the minimal glycolysis ( mg ) strain . the order of gene deletions that led to the final mg strain imx372 was glk1 , hxk1 , tdh1 , tdh2 , gpm2 , gpm3 , eno1 , pyk2 , pdc5 , pdc6 , adh2 , adh5 , and adh4 . the genes glk1 , hxk1 , tdh1 , and tdh2 were deleted using the auxotrophic and dominant markers sphis5 ( 44 , 45 ) , klleu2 ( 45 , 46 ) , kanmx ( 47 , 48 ) , and hphnt1 ( 49 , 50 ) , respectively . gpm2 to adh4 deletions were performed using a strategy of selection and counterselection with the klura3/5-fluoroorotic acid ( 5-foa ) system ( 51 ) for the recovery of the marker module . the dominant marker modules kanmx and hphnt1 were removed using deletion cassettes containing klura3 and amdsym , respectively . these markers were sequentially removed as reported previously ( 54 ) . to obtain a prototrophic strain , a cassette containing the marker module klura3 was integrated in the tdh1 locus ; this generated the mg strain . all integrative cassettes were constructed using phusion hot start ii high - fidelity polymerase ( thermo scientific , landsmeer , the netherlands ) following the manufacturer 's recommendations and with the primer pairs and plasmids listed in tables s4 and s5 in the supplemental material , respectively , as the templates . correct integrations , deletions , and marker excision were confirmed by pcr using dreamtaq polymerase ( thermo scientific ) and following the manufacturer 's recommendations . genomic dna that served as the template for pcr was obtained by extraction with 0.05 n naoh directly from single colonies or by purification using the yeastar genomic dna kit ( zymo research , orange , ca ) following the manufacturer 's recommendations . all pcr products were loaded on gels containing 1% ( wt / vol ) agarose ( thermo scientific ) and 1 tris - acetate - edta ( tae ) buffer ( thermo scientific ) . integrative cassettes were gel purified using the zymoclean gel dna recovery kit ( zymo research ) . yeast strain transformations were performed with the lithium acetate protocol as previously described ( 55 ) . plasmids were isolated from e. coli with the genelute plasmid miniprep kit ( sigma - aldrich , st . complex and nonselective media for growth rate determination and propagation contained 10 g liter yeast extract , 20 g liter peptone , and 20 g liter glucose ( ypd ) . when selection was required , ypd was supplemented with 200 mg liter g418 ( ypd+g418 ) or 200 mg liter hygromycin ( ypd+hyg ) . synthetic medium ( sm ) containing 3 g liter kh2po4 , 0.5 g liter mgso47h2o , 5 g liter ( nh4)2so4 , 1 ml liter of a trace element solution , and 1 ml liter of a vitamin solution as previously described ( 56 ) . sm was supplemented with 20 g liter glucose ( smg ) for propagation , growth rate determination , and batch cultures . when auxotrophic strains were cultivated , smg was supplemented with 150 mg liter uracil , 125 mg liter histidine , and/or 500 mg liter leucine , according to strain auxotrophy ( 57 ) . selection of strains lacking the klura3 marker was done in smg supplemented with 150 mg liter uracil and 1 g liter 5-fluoroorotic acid ( smg5-foa ) . unless otherwise stated , liquid cultures were performed in 500-ml shake flasks with a working volume of 100 ml and incubated at 30c and 250 rpm , the targeted starting optical density at 660 nm ( od660 ) , measured using a libra s11 spectrophotometer ( biocrom , cambridge , united kingdom ) , being 0.5 . unless otherwise stated , cultures on solid media were incubated for 3 days at 30c . growth of the mg strain and of the prototrophic control strain cen.pk113-7d was performed on agar plates containing ypd or sm plus 2% glucose ( smg ) , sm plus 2% galactose ( smgal ) , sm plus 2% maltose ( smmal ) , sm plus 2% sucrose ( smsuc ) , sm plus 3% ( vol / vol ) ethanol ( smetoh ) , sm plus 6% glycerol ( smglyc ) , sm plus 5% glucose ( smhg ) , sm plus 0.5% glucose ( smlg ) , smg plus 2.5 , 5 , or 10 mm licl ( smg - lit ) , smg plus 200 or 500 mm nacl ( smgsa ) , smg plus 50 , 100 , or 200 m cdcl2 ( smgcad ) , smg plus 1 , 1.5 , or 2 m sorbitol ( smgsor ) , smg at ph 4 , 6 , or 7.5 ( smgph ) , or smg plus 1 mm h2o2 ( smgox ) . spot plates were inoculated with 10 to 10 cells obtained by serial dilutions of liquid cultures grown to the exponential phase on smg . all plates were incubated at 30c for 3 days , except for plates containing smglyc , which were incubated for 6 days . growth rate determination was performed by od660 measurement of cultures grown in 500-ml shake flasks containing 100 ml of ypd , smg , or smetoh . cultures for growth rate measurements were inoculated with precultures grown until the late exponential phase under identical conditions ( i.e. , the same medium and temperature ) . glucose / ethanol / glucose switches were started by inoculation of smg shake flasks with a preculture grown overnight on the same medium . at the mid - exponential phase ( od660 of ca . 4 ) samples were taken , washed twice with sterile demineralized water , and used to inoculate shake flasks containing 100 ml of smetoh . when reaching an od660 of 3 , a sample from the culture was washed twice with sterile demineralized water and used to inoculate a 500-ml shake flask containing 100 ml smg . growth was then monitored until the late exponential phase . to evaluate tolerance of mg to oxidative stress , cells grown to the exponential phase in shake flasks containing smg were transferred to a shake flask containing smg supplemented with 4.5 mm h2o2 . samples were then taken every 20 min for 1 h , diluted , and plated on smg plates to reach approximately 100 cells per plate . aerobic and anaerobic batch cultures were performed in 2-liter laboratory fermenters ( applikon , schiedam , the netherlands ) with a 1-liter working volume . antifoam emulsion c ( sigma ) was added to a concentration of 0.2 g liter from a 20% ( vol / vol ) solution autoclaved separately ( 121c ) . prior to inoculation , glucose was added to a final concentration of 20 g liter from a sterile ( 110c ) 50% glucose solution and 1 ml of a vitamin solution ( 56 ) . anaerobic cultures were supplemented with 0.01 g liter ergosterol and 0.42 g liter tween 80 dissolved in ethanol as previously described ( 56 ) . compressed air or gaseous nitrogen ( quality , 4.5 ; > 99.995% vol n2 , < 10 volumes per million [ vpm ] o2 pollution ) ( linde gas , schiedam , the netherlands ) for aerobic and anaerobic cultures , respectively , was sparged to the bioreactor at a rate of 500 ml min via an ion science saga digital flow meter ( cambridge , united kingdom ) . precultures were started by inoculation of smg shake flasks with 1 stock vial of the appropriate strain . 16 h of incubation , these cultures were used to inoculate new smg shake flasks . the starting od660 was tightly controlled to reach an od660 of 1 after 8 h of incubation . these exponentially growing cultures were washed twice with demineralized water ( spinning at 5,000 rpm for 4 min ) and resuspended in 100 ml demineralized water . this was the cell suspension that was used to inoculate fermenters with a starting od660 of 0.1 . the concentration of extracellular metabolite in culture supernatants was measured by high - performance liquid chromatography ( hplc ) analysis using a aminex hpx-87h ion - exchange column operated at 60c with 5 mm h2so4 as the mobile phase at a flow rate of 0.6 ml min . agarose plugs for the different strains were prepared using the contour - clamped homogenous electric field ( chef ) yeast genomic dna plugs kit ( bio - rad , richmond , ca ) , following the manufacturer 's recommendations , and used for chef electrophoresis . mg and cen.pk113-7d vacuoles were stained with the red fluorescent dye fm4 - 64 ( excitation / emission , 515/640 nm ) ( thermo fisher scientific ) following the manufacturer 's recommendations . yeast cells and vacuoles were visualized with an imager - z1 microscope equipped with an axiocam mr camera , an ec plan - neofluar 100/1.3 oil ph3 m27 objective , and the filter set bp 535/25 , ft 580 , and lp 590 ( carl - zeiss , oberkochen , germany ) . samples equivalent to 62.5 mg of dry weight biomass taken at the mid - exponential phase ( ca . 10 h after inoculation ) of aerobic batch cultures were used to obtain cell extracts as previously described ( 58 ) . measurement of the activity of glycolytic enzymes was carried out as previously described ( 59 ) , except for phosphofructokinase , the activity of which was determined as described in reference 60 . enzyme activities are expressed as micromoles of substrate per minute per milligram of protein or units per milligram of protein . protein concentrations in cell extracts were determined according to reference 61 with bovine serum albumin as a standard . samples ( 1.2 ml each ) were taken from aerobic batch cultures using a rapid sampling setup ( 62 ) and placed directly into vials containing 5 ml 80% methanol precooled at 40c . the samples were washed with precooled 80% methanol , and the extraction was performed with boiling ethanol as previously described ( 63 ) . c - labeled cell extract was used as an internal standard for metabolite quantification ( 64 ) . intracellular glucose ( gluc ) , glucose-6-phosphate ( gluc-6-p ) , fructose-6-phosphate ( fruc-6-p ) , dihydroxyacetone phosphate ( dhap ) , glyceraldehyde-3-phosphate ( gap ) , 3-phosphoglycerate ( 3pg ) , 2-phosphoglycerate ( 2pg ) , and pyruvate ( pyr ) were measured by gas chromatography - mass spectrometry ( gc - ms ) according to reference 65 . intracellular fructose-1,6-bisphosphate ( fruc-1,6-bp ) and phosphoenolpyruvate ( pep ) were measured by liquid chromatography - mass spectrometry ( lc - ms ) according to reference 66 . samples for transcriptome sequencing ( rna - seq ) were obtained from aerobic batch cultures at the mid - exponential phase of growth on glucose ( ca . sampling , rapid quenching in liquid nitrogen , and rna extraction were performed as previously described ( 68 ) . sequencing was performed with the illumina hiseq 2500 and carried out by baseclear ( leiden , the netherlands ) . a data set of 51-bp single reads of at least 1 the data were aligned to the reference using the burrow - wheeler alignment tool bwa ( 69 , 70 ) . gene expression levels were estimated using fragments per kilobase per million ( fpkm ) values by the cufflinks software ( 71 ) . to identify differential gene expression between strains cen.pk113-7d and imx372 , rna - seq data comparison was performed and statistically assessed using cuffdiff ( 71 ) . high - quality genomic dna of the strains cen.pk102-12a and imx372 ( mg ) was obtained using the qiagen 100/g kit ( qiagen , hilden , germany ) following the manufacturer 's recommendations . libraries of 300-bp inserts were constructed and paired end sequenced ( 100-bp reads ) using an illlumina hiseq 2500 sequencer ( baseclear bv , leiden , the netherlands ) . a minimum data quantity of 950 mb was generated , representing a minimum 80-fold coverage . the sequence reads were mapped onto cen.pk113-7d genome ( 42 ) using the burrows - wheeler alignment tool bwa and further processed using samtools ( 69 , 70 ) . default settings were used , except that the maximum read depth was set to 400 ( d400 ) . to minimize false - positive mutation calls , custom perl scripts were used for further mutation filtering as follows : ( i ) mutation calls containing ambiguous bases in mapping consensus were filtered out , ( ii ) only the single - nucleotide variations with a quality of at least 20 were kept ( with variant quality defined as the phred - scaled probability that the mutation call is incorrect [ 72 , 73 ] ) , ( iii ) mutations with a depth of coverage below 10 were discarded , and ( iv ) the mutations found in cen.pk102-12a were subtracted from the list sequence . eventually , the single - nucleotide variations were physically positioned and functionally annotated according to the cen.pk113-7d sequence annotation . the rna - seq data generated in this study have been submitted to the genome expression omnibus database and assigned accession no . the sequence data generated in this study are searchable at ncbi  entrez ( http://www.ncbi.nlm.nih.gov/ ) under bioproject prjna269221 . to enable construction of a  minimal glycolysis  ( mg ) yeast strain , the first goal was to identify the major paralogs that should be retained in the final strain design . this assessment was based on information from the literature on phenotypes of relevant deletion mutants and on nonglycolytic roles (  moonlighting functions  ) ( 74 ) of glycolytic isoenzymes . furthermore , a compendium of s. cerevisiae transcriptome data obtained under a wide range of controlled cultivation conditions ( 75 ) was used to compare expression profiles of glycolytic paralogs . major paralogs were elected based on the following , nonexclusive criteria ( fig . 1 ; see fig . s1 in the supplemental material ) : ( i ) the highest transcript level over a range of growth conditions , ( ii ) the most extensive loss of enzyme activity in cell extracts upon deletion , ( iii ) moonlighting functions with a strong impact on specific growth rate or robustness ( hxk2 , eno1 , and eno2 ) ( 76 , 77 ) , and ( iv ) the strongest decrease in specific growth rate upon deletion . based on these criteria , 11 of the 27 glycolytic genes in s. cerevisiae were assessed to be major paralogs : hxk2 , pgi , fba1 , tpi1 , tdh3 , gpm1 , eno2 , pgi1 , pyk1 ( cdc19 ) , pdc1 , and adh1 , while 13 minor paralogs ( i.e. , hxk1 , glk1 , tdh1 , tdh2 , gpm2 , gpm3 , eno1 , pyk2 , pdc5 , pdc6 , adh2 , adh4 , and adh5 ) were selected for gene deletion . in s. cerevisiae , the pkf1 and pfk2 paralogs have evolved into subunits of a hetero - octameric phosphofructokinase ( 78 ) . since deletion of either gene substantially decreases fitness ( 7981 ) ( see fig . s1 in the supplemental material ) , both were retained in the minimal glycolysis design . adh3 encodes a mitochondrial alcohol dehydrogenase involved in an anaerobic redox shuttle across the mitochondrial inner membrane ( 82 ) . to prevent reduced growth rates under anaerobic conditions , sequential deletion of the 13 minor paralogs in a haploid s. cerevisiae strain belonging to the cen.pk family ( 4143 ) yielded the prototrophic strain imx372 , which we will refer to as the  mg resequencing of the genome of the mg strain confirmed the correct deletion of all 13 minor glycolytic genes . although transformation of s. cerevisiae can be mutagenic ( 83 , 84 ) , only 20 single - nucleotide differences were detected in the mg strain relative to its ancestor strain , cen.pk102-12a . eleven of these differences occurred within open reading frames ( orfs ) , of which 10 resulted in amino acid changes ( see table s1 in the supplemental material ) . whole - genome sequencing and karyotyping indicated the duplication of two short sections of chromosome iii ( 16.5 kbp , from ycr019w to ycr027c ) and chromosome v ( 19.1 kbp , from yer093c - a to yer104w ) ( see fig . these regions do not carry genes involved in central carbon metabolism , and no interchromosomal rearrangements were observed . to explore the physiological impact of the simultaneous deletion of all 13 minor glycolytic paralogs , we quantitatively compared specific growth rates , product formation , and gene expression in the mg strain and in a congenic , prototrophic reference strain with a full complement of glycolytic genes . in these studies , a synthetic , chemically defined medium was used , since complex media do not enable cells to express their full genetic potential ( 33 ) and complicate quantitative physiological analysis . in aerobic , glucose - grown bioreactor batch cultures , specific rates of growth , substrate consumption , and product formation were not significantly affected by deletion of all 13 minor glycolytic paralogs ( fig . the only exception to this observation concerned acetate production , which was slightly higher in the mg strain . also after the diauxic shift , where the glycolytic flux changes direction as the aerobic yeast cultures consumed the ethanol and acetate produced during the initial growth phase on glucose , the biomass formation and substrate consumption rates of the two strains were virtually identical . in anaerobic yeast cultures , the absence of oxidative phosphorylation makes glycolysis the only pathway for energy conservation . also under these more demanding conditions , the specific growth rate and metabolic fluxes of the mg strain did not differ significantly from those of the congenic reference strain ( fig . biomass production and extracellular metabolite profiles from aerobic and anaerobic batch cultures in bioreactors of the minimal glycolysis ( mg ) strain and a congenic reference strain ( indicated by  c  ) . shown are biomass and extracellular metabolite profiles from aerobic ( a ) and anaerobic ( b ) controlled batch cultures of the minimal glycolysis strain ( mg ) ( open circles ) and the prototrophic reference strain cen.pk113-7d ( closed circles ) . the data shown in the graphs are the average and average deviation of the mean from two independent cultures for each strain . the specific growth rate (  [ per hour ] ) and biomass - specific rates of glucose consumption ( qs ) , ethanol production ( qetoh ) , glycerol production ( qglyc ) , and acetate production ( qace ) ( all expressed in millimoles per gram dry weight per hour ) represent the average and standard deviation of data from four independent cultures for each strain . * , statistically significant ( by two - tailed t test assuming the same variance in the populations ) differences between the two tested strains ( p = 0.03 and p = 0.02 for  and qace , respectively , in aerobic cultures ) . to further compare growth of the mg strain in aerobic bioreactor batch cultures , its glycolytic enzyme activities in cell extracts , intracellular metabolite concentrations , and transcriptome remarkably , no significant differences were observed in the transcript levels of any of the major glycolytic paralogs or in vitro glycolytic enzyme activities ( fig . 3 ) . with the exception of slightly higher intracellular concentrations of acetaldehyde , concentrations of glycolytic intermediates in the mg strain did not significantly differ from those in the reference strain . to explore potential impacts of the deletion of 13 minor glycolytic paralogs outside glycolysis , genome - wide transcript levels of the mg strain and the reference strain were compared . as few as 17 genes showed a significantly different transcript level ( see table s2 in the supplemental material ) , 12 of which were located on the duplicated regions on chromosomes iii and v. an in - depth , comprehensive analysis of the mg strain in glucose - grown bioreactor batch cultures therefore failed to identify substantial impacts of the minor glycolytic paralogs on fluxes , intracellular metabolite concentrations , or gene expression . in vitro enzyme activities and intracellular metabolite concentrations in aerobic batch cultures in bioreactors of the minimal glycolysis ( mg ) strain and a congenic reference strain . ( a ) average values of four independent culture replicates of the in vitro activities for the glycolytic enzymes in the mg strain ( white bars ) and the prototrophic reference strain cen.pk113-7d ( black bars ) . hxk , hexokinase / glucokinase ; pgi , phosphoglucose isomerase ; pfk , phosphofructokinase ; fbpa , fructose - bisphosphate aldolase ; tpi , triose - phosphate isomerase ; gapdh , glyceraldehyde-3-phosphate dehydrogenase ; pgk , phosphoglycerate kinase ; gpm , phosphoglycerate mutase ; eno , enolase ; pyk , pyruvate kinase ; pdc , pyruvate decarboxylase ; adh , alcohol dehydrogenase . ( b ) intracellular glycolytic metabolite profiles of the mg strain ( open circles ) and of cen.pk113-7d ( close circles ) from aerobic batch cultures . average values from two independent culture replicates are shown , and the average deviations of the mean are indicated by error bars . the vertical orange dotted line indicates the time at which glucose was depleted . if during evolution of s. cerevisiae , its glycolytic paralogs have evolved different roles through sub- or neofunctionalization , their deletion may only cause an observable phenotype under specific growth conditions . the growth rates of the mg and reference strains were therefore compared under a wide range of selected growth conditions . during fast growth in shake flasks on complex ( yeast extract - peptone - glucose [ ypd ] ) medium ( fig . 4 ) , the glycolytic pathway operates at rates closer to its maximum capacity than during growth at 30c on synthetic medium ( 60 , 85 ) . however , no difference in specific growth rates between the mg and reference strains was observed under these conditions . the mg strain also showed the same growth rate as the reference strain at high temperature ( 37c ) ( fig . growth of the minimal glycolysis ( mg ) strain and a congenic reference strain under different growth conditions in shake flasks . ( a ) specific growth rates of s. cerevisiae strains mg ( white bars ) and cen.pk113-7d ( reference strain , black bars ) . smg , synthetic medium with glucose ; smetoh , synthetic medium with ethanol as the carbon source ; ypd , complex medium with glucose as the carbon source . the labels ( b ) growth , measured as change in the culture 's optical density at 660 nm ( od660 ) , of the mg strain ( open circles ) and reference strain ( closed circles ) during carbon source switches . all data represent the average and average deviation of the mean from two independent culture replicates . absence of the minor glycolytic paralogs did not affect growth of the mg strain on ethanol ( fig . 2a and 4a ) . these results supported the conclusion from the aerobic bioreactor batch cultures that efficient gluconeogenesis does not require any of the minor paralogs ( as has previously been proposed for adh2 ) . while minor glycolytic paralogs apparently do not contribute to maximum specific growth rates on individual carbon sources , they might be involved in substrate transitions . for instance , pyk2 , the glucose - repressed and fructose-1,6-bisphosphate - insensitive paralog of pyk1 ( 21 ) , has been proposed to prevent futile cycling resulting from simultaneous operation of glycolysis and gluconeogenesis during transitions between fermentable and nonfermentable carbon sources ( 87 , 88 ) . however , no significant impact of the deletion of minor glycolytic paralogs was observed during transitions from glucose to ethanol and back to glucose ( fig . 4b ) . to investigate possible phenotypes of the mg strain under a wider range of environmental conditions , growth on solid medium 5p to w ) , and at high concentrations of lithium and sodium ions , to which strains from the cen.pk lineage are hypersensitive ( 89 ) ( fig . 5h to l ) . only growth with glycerol as a carbon source indicated a slightly reduced growth rate of the mg strain . however , this phenotype was not observed during growth on liquid medium with glycerol ( see fig . s3 in the supplemental material ) and may therefore have been caused by a contaminating substrate in the agar used for the plate experiments . growth of the minimal glycolysis ( mg ) strain and a congenic reference strain on different solid media . serial dilutions of cell suspensions of the mg strain and of the reference strain cen.pk113-7d were plated on agar media with 24 different compositions . the conditions included osmotic stress ( a to c ) , different phs ( d to f ) , oxidative stress ( g ) , three different salts ( h to o ) , six different carbon sources ( p to w ) , and complex medium ( x ) . with the exception of panel x ( ypd ) , all plates contained synthetic medium ( sm ) . all plates contained 2% glucose , with the exception of plates with different carbon sources , which contained galactose ( 2% ) , maltose ( 2% ) , sucrose ( 2% ) , ethanol ( 3% , vol / vol ) , or glycerol ( 6% ) and plates with low ( 0.03 m [ 5% ] ) and high ( 0.3 m [ 50% ] ) concentrations of glucose . colonies of reference and mg strains are denoted by black bands and orange bands , respectively , at the bottom of each panel . additional growth conditions tested on solid media were chosen based on information from the literature ( 74 , 90 ) on sub- or neofunctionalization of specific deleted glycolytic paralogs . in addition to their enolase activity , eno1 and eno2 have both been implicated in vacuolar assembly , but in the absence of eno1 , eno2 is able to support both functions ( 76 ) . in the vacuole and the atpase associated with it , disruption causes a deleterious phenotype in s. cerevisiae when cultivated under alkaline conditions ( 91 ) . as expected , the mg strain grew normally at ph 7.5 , thereby indicating the absence of major vacuolar malfunction ( fig . tdh3 and tdh2 , the two minor isoenzymes of glyceraldehyde-3-phosphate dehydrogenase , have been proposed to protect s. cerevisiae against oxidative stress because of their different levels of thiolation ( 92 ) . however , exposure to oxidative stress by growth on hydrogen peroxide - containing plates ( fig . 5 g ) showed that the absence of tdh2 did not affect the oxidative stress resistance of the mg strain . sub- and neofunctionalization have resulted in a different glucose - dependent regulation of glk1/hxk1 and pyk2 and in different regulatory properties of the proteins that they encode from those encoded by hxk2 and pyk1 , respectively ( 21 , 93 ) . however , growth of the mg strain at low , intermediate , and high glucose concentrations was not impaired relative to that of the reference strain ( fig . finally , pdc6 , which encodes a pyruvate decarboxylase isoenzyme with a low cysteine and methionine content , is specifically induced under sulfur limitation ( 94 ) . cultivation in the presence of cadmium increases abundance of pdc6 ( 9496 ) ; since cells require a high level of glutathione production for detoxification , sulfur amino acids are then used for this process and are less available for protein synthesis . despite the absence of pdc6 , the mg strain grew as well as the reference strain in the presence of cadmium . deletion of 13 of the 27 glycolytic paralogs in s. cerevisiae yielded a  minimal glycolysis  ( mg ) yeast strain whose most spectacular characteristic was the absence of any pronounced phenotype under a wide range of laboratory growth conditions . one of the hypotheses that has been proposed to explain retention of functionally redundant paralogs during evolution is a contribution to gene dosage and , thereby , to the capacity of the pathway or process in which they operate ( 97 , 98 ) . the high glycolytic rates in anaerobic cultures of the mg strain ( 18 mmol glucose per g of dry biomass per h ) did not significantly differ from those of wild - type strains ( 99 ) ( fig . 2 ) . this result argues against gene dosage effects as a means for fixing minor glycolytic paralogs in the yeast genome . instead , our observations indicate that duplication of glycolytic genes is not a prerequisite for achieving the high glycolytic fluxes and fermentative capacities that are characteristic of s. cerevisiae and essential for many of its industrial applications ( 100 , 101 ) . it might be argued that the gene dosage hypothesis does apply to phosphofructokinase , as deletion of either pfk1 or pfk2 substantially reduces enzyme activity and fitness ( 7981 ) ( see fig . however , the presence of pfk1 and pfk2 in a hetero - octameric complex can also be seen as a case of neofunctionalization , in which two paralogs have been fixed by an acquired mutual dependency . the near - wild - type growth kinetics of the mg strain in aerobic and anaerobic cultures is difficult to reconcile with the hypothesis that duplication of glycolytic genes during the wgd event played a major role in increasing its glycolytic capacity or in causing the phenomenon of aerobic fermentation ( the crabtree effect ) ( 20 ) . our inability to identify a phenotype after deletion of all minor paralogs of glycolytic genes in s. cerevisiae does not imply that such a phenotype does not exist . the range of conditions tested represents an infinitesimal fraction of the environmental conditions to which s. cerevisiae may have been exposed to in its evolutionary history . the absence of a clear phenotype under standard laboratory conditions makes the mg strain an even more interesting platform for future high - throughput studies to investigate its phenotype under more conditions . other characteristic features of s. cerevisiae like tolerance to high ethanol concentrations or growth at high gravity are interesting conditions to test , as well as the robustness of the mg strain under nonstandard conditions , such as sporulation , starvation , severe calorie restriction , and dynamic nutrient supply regimens . however , guessing the environmental factors that have conferred an evolutionary advantage to strains carrying gene duplications presents a formidable challenge , more particularly because the real niche for s. cerevisiae is still a matter of debate , to the extent that s. cerevisiae may be considered a generalist that does not necessarily favor a specific type of environment ( 102 ) . for further functional analysis studies of the minor glycolytic paralogs , the set of intermediate strains used for construction of the mg strains ( see table s3 in the supplemental material ) can be used to rapidly identify which minor paralog or paralogs contribute to any newly identified phenotypes . moreover , the sensitivity of fitness analyses can be improved  for example , by competitive cultivation of wild - type and mg strains in mixed cultures . in s. cerevisiae , overrepresentation of paralogous gene sets is not unique for glycolysis ( 23 of 27 glycolytic genes are paralogs , corresponding to 85% ) but also occurs for metabolic genes in general . of all metabolic genes , 44% are paralogs ( 98 ) , while this percentage is only 5% for the entire yeast genome ( 5 , 17 ) . currently , a large research effort is under way to enable design and assembly of entirely synthetic yeast genomes ( 36 ) . for the design of compact , functional synthetic genomes , it will be highly interesting to investigate whether our results on yeast glycolysis can be extrapolated to minor paralogs of genes encoding key enzymes in other central metabolic pathways . at present , predictions of the outcome of such experiments the clear impact of  minor paralogs  of transaldolase and transketolase on pentose fermentation kinetics by engineered s. cerevisiae ( 103 ) is a clear example that , at least under some conditions , minor paralogs can make a clear contribution to metabolic flux . the deletion of all minor glycolytic paralogs had a surprisingly small effect on the yeast transcriptome in aerobic , glucose - grown batch cultures . this result is in marked contrast with several studies in which deletion of individual major glycolytic paralogs led to transcriptional upregulation of its minor paralogs . upregulation of minor paralogs is the upregulation of pdc5 upon deletion of pdc1 ( 104 ) . such an asymmetric cross - regulation is consistent with the backup theory for fixation of duplicated genes ( 17 ) , which postulates that ( minor ) paralogs provide a buffer against deleterious mutations . although upregulation of minor paralogs and even repair of a null mutation in a major glycolytic paralog by recombination with a minor paralog have been demonstrated ( pdc1/pdc5 ) ( 104 ) , the general evolutionary significance of the backup theory is still a matter of debate ( 98 , 105 ) . the selective advantages responsible for the overrepresentation of paralogous genes among the structural genes of yeast glycolysis therefore remain an intriguing conundrum . as the mg strain was constructed before the advent of crispr - cas9 technology ( 106 ) , its construction involved 16 consecutive rounds of transformation and 11 marker recycling steps . a previous study ( 107 ) in which all 20 hexose transporter ( hxt ) genes in s. cerevisiae were deleted by repeated rounds of transformation and marker recycling resulted in massive genomic rearrangements , which were recently explained from repeated use of the loxp system for marker recovery ( 120 ) . repeated use of the loxp system enables recombination across loxp sites that are left in the genome after marker removal ( 108 ) . strain is mainly used as an excellent platform for functional analysis of individual transporter strains ( 107 , 109117 ) , in which role its chromosomal rearrangements are hardly relevant . while the mg strain similarly offers an interesting platform for functional analysis of ( heterologous ) glycolytic genes by one - step gene replacement , we designed and constructed it with the specific aim of building a robust platform for quantitative studies in systems biology . by using ura3 and amdsym ( 53 ) as counterselectable , recyclable marker genes , furthermore , interference of auxotrophies in the interpretation of quantitative growth studies ( 57 ) was avoided by making the mg strain prototrophic . the availability of a well - defined yeast platform with a minimal complement of glycolytic enzymes should provide clear advantages for quantitative modeling of the kinetics and regulation of glycolysis , as it eliminates the intrinsic uncertainties caused by the simultaneous , context - dependent expression of different isoenzymes . in view of the important role of glucose transport in the kinetics of yeast glycolysis , we are currently constructing mg variants with a single hexose transporter . in addition to providing a relevant test bed for mathematical modeling of glycolysis , the mg strain provides an interesting , simplified starting point for laboratory studies on yeast glycolysis . previous studies have indicated that long - term laboratory evolution can have a major impact on glycolytic genes and their expression ( 59 , 118 ) . comparison of evolution of the mg strain with that of strains that carry a full complement of glycolytic genes provides an interesting starting point for studies on the impact of gene duplication on evolutionary flexibility of a ubiquitous central metabolic pathway .
as a result of ancestral whole - genome and small - scale duplication events , the genomes of saccharomyces cerevisiae and many eukaryotes still contain a substantial fraction of duplicated genes . in all investigated organisms , metabolic pathways , and more particularly glycolysis , 
 are specifically enriched for functionally redundant paralogs . in ancestors of the saccharomyces lineage , the duplication of glycolytic genes is purported to have played an important role leading to s. cerevisiae 's current lifestyle favoring fermentative metabolism even in the presence of oxygen and characterized by a high glycolytic capacity . in modern s. cerevisiae strains , 
 the 12 glycolytic reactions leading to the biochemical conversion from glucose to ethanol are encoded by 27 paralogs . in order to experimentally explore the physiological role of this genetic redundancy , 
 a yeast strain with a minimal set of 14 paralogs was constructed ( the  minimal glycolysis  [ mg ] strain ) . remarkably , a combination of a quantitative systems approach and semiquantitative analysis in a wide array of growth environments revealed the absence of a phenotypic response to the cumulative deletion of 13 glycolytic paralogs . 
 this observation indicates that duplication of glycolytic genes is not a prerequisite for achieving the high glycolytic fluxes and fermentative capacities that are characteristic of s. cerevisiae and essential for many of its industrial applications and argues against gene dosage effects as a means of fixing minor glycolytic paralogs in the yeast genome . 
 the mg strain was carefully designed and constructed to provide a robust prototrophic platform for quantitative studies and has been made available to the scientific community .