The development of biocatalysts capable of fermenting xylose, a five-carbon sugar abundant in lignocellulosic biomass, is a key step to achieve a viable production of second-generation ethanol. In this work, a robust industrial strain of Saccharomyces cerevisiae was modified by the addition of essential genes for pentose metabolism. Subsequently, taken through cycles of adaptive evolution with selection for optimal xylose utilization, strains could efficiently convert xylose to ethanol with a yield of about 0.46 g ethanol/g xylose. Though evolved independently, two strains carried shared mutations: amplification of the xylose isomerase gene and inactivation of ISU1, a gene encoding a scaffold protein involved in the assembly of iron-sulfur clusters. In addition, one of evolved strains carried a mutation in SSK2, a member of MAPKKK signaling pathway. In validation experiments, mutating ISU1 or SSK2 improved the ability to metabolize xylose of yeast cells without adaptive evolution, suggesting that these genes are key players in a regulatory network for xylose fermentation. Furthermore, addition of iron ion to the growth media improved xylose fermentation even by non-evolved cells. Our results provide promising new targets for metabolic engineering of C5-yeasts and point to iron as a potential new additive for improvement of second-generation ethanol production.

Significance Individual bacteria that share identical genomes and growth environments can display substantial cell-to-cell differences in expression of stress-response genes and single-cell growth rates. This phenotypic heterogeneity can impact the survival of single cells facing sudden stress. However, the windows of time that cells spend in vulnerable or tolerant states are often unknown. We quantify the temporal expression of a suite of stress-response reporters, while simultaneously monitoring growth. We observe pulsatile expression across genes with a range of stress-response functions, finding that single-cell growth rates are often anticorrelated with reporter levels. These dynamic phenotypic differences have a concrete link to function, in which individual cells undergoing a pulse of elevated expression and slow growth are predisposed to survive antibiotic exposure.

The Saccharomyces cerevisiae strain JAY270/PE2 is a highly efficient biocatalyst used in the production of bioethanol from sugarcane feedstock. This strain is heterothallic and diploid, and its genome is characterized by abundant structural and nucleotide polymorphisms between homologous chromosomes. One of the reasons it is favored by many distilleries is that its cells do not normally aggregate, a trait that facilitates cell recycling during batch-fed fermentations. However, long-term propagation makes the yeast population vulnerable to the effects of genomic instability, which may trigger the appearance of undesirable phenotypes such as cellular aggregation. In pure cultures of JAY270, we identified the recurrent appearance of mutants displaying a mother-daughter cell separation defect resulting in rough colonies in agar media and fast sedimentation in liquid culture. We investigated the genetic basis of the colony morphology phenotype and found that JAY270 is heterozygous for a frameshift mutation in the ACE2 gene (ACE2/ace2-A7), which encodes a transcriptional regulator of mother-daughter cell separation. All spontaneous rough colony JAY270-derived isolates analyzed carried copy-neutral loss-of-heterozygosity (LOH) at the region of chromosome XII where ACE2 is located (ace2-A7/ace2-A7). We specifically measured LOH rates at the ACE2 locus, and at three additional chromosomal regions in JAY270 and in a conventional homozygous diploid laboratory strain. This direct comparison showed that LOH rates at all sites were quite similar between the two strain backgrounds. In this case study of genomic instability in an industrial strain, we showed that the JAY270 genome is dynamic and that structural changes to its chromosomes can lead to new phenotypes. However, our analysis also indicated that the inherent level of genomic instability in this industrial strain is normal relative to a laboratory strain. Our work provides an important frame of reference to contextualize the interpretation of instability processes observed in the complex genomes of industrial yeast strains.

Cells carrying deletions of genes encoding H-class ribonucleases display elevated rates of chromosome instability. The role of these enzymes is to remove RNA-DNA associations including persistent mRNA-DNA hybrids (R-loops) formed during transcription, and ribonucleotides incorporated into DNA during replication. RNases H1 and H2 can degrade the RNA component of R-loops, but only RNase H2 can initiate accurate ribonucleotide excision repair (RER). In order to examine the specific contributions of these activities to chromosome stability, we measured rates of loss-of-heterozygosity (LOH) in diploid Saccharomyces cerevisiae yeast strains carrying the rnh201-RED separation-of-function allele, encoding a version of RNase H2 that is RER-defective, but partly retains its other activity. The LOH rate in rnh201-RED was intermediate between RNH201 and rnh201Δ. In strains carrying a mutant version of DNA polymerase ε (pol2-M644G) that incorporates more ribonucleotides than normal, the LOH rate in rnh201-RED was as high as the rate measured in rnh201Δ. Topoisomerase 1 cleavage at sites of ribonucleotide incorporation has been recently shown to produce DNA double strand breaks. Accordingly, in both the POL2 and pol2-M644G backgrounds, the LOH elevation in rnh201-RED was suppressed by top1Δ. In contrast, in strains that incorporate fewer ribonucleotides (pol2-M644L) the LOH rate in rnh201-RED was low and independent of topoisomerase 1. These results suggest that both R-loop removal and RER contribute substantially to chromosome stability, and that their relative contributions may be variable across different regions of the genome. In this scenario, a prominent contribution of R-loop removal may be expected at highly transcribed regions, whereas RER may play a greater role at hotspots of ribonucleotide incorporation.

Conventional models of genome evolution are centered around the principle that mutations form independently of each other and build up slowly over time. We characterized the occurrence of bursts of genome-wide loss-of-heterozygosity (LOH) in Saccharomyces cerevisiae, providing support for an additional nonindependent and faster mode of mutation accumulation. We initially characterized a yeast clone isolated for carrying an LOH event at a specific chromosome site, and surprisingly found that it also carried multiple unselected rearrangements elsewhere in its genome. Whole-genome analysis of over 100 additional clones selected for carrying primary LOH tracts revealed that they too contained unselected structural alterations more often than control clones obtained without any selection. We also measured the rates of coincident LOH at two different chromosomes and found that double LOH formed at rates 14- to 150-fold higher than expected if the two underlying single LOH events occurred independently of each other. These results were consistent across different strain backgrounds and in mutants incapable of entering meiosis. Our results indicate that a subset of mitotic cells within a population can experience discrete episodes of systemic genomic instability, when the entire genome becomes vulnerable and multiple chromosomal alterations can form over a narrow time window. They are reminiscent of early reports from the classic yeast genetics literature, as well as recent studies in humans, both in cancer and genomic disorder contexts. The experimental model we describe provides a system to further dissect the fundamental biological processes responsible for punctuated bursts of structural genomic variation.

Cell-to-cell heterogeneity in gene expression and growth can have critical functional consequences, such as determining whether individual bacteria survive or die following stress. However, the timescales of the dynamics that underlie this heterogeneity are often unknown. This is a critical piece of information because different phenotypic outcomes can arise from gradual versus rapid changes in expression and growth. Using single-cell time-lapse microscopy, we conducted detailed measurements of gene expression and growth over many generations using 15 reporters in Escherichia coli, focusing on genes related to stress response. In constant environmental conditions without stress, we found many examples of pulsatile gene expression, suggesting that this may be a widespread dynamic property. Pulse lengths often exceeded the cell cycle, leading to multi-generational correlations. Temporal properties of the pulses, such as their frequency, duration, and amplitude, varied widely across the reporters. Single-cell growth rates were often anti-correlated with gene expression, with changes in growth typically preceding changes in expression. These dynamics and the timescales of the pulses have functional consequences, which we demonstrate by measuring single-cell survival after challenging cells with the antibiotic ciprofloxacin. Together, this work reveals that pulsatile expression dynamics are widespread in E. coli stress response networks. Further, these dynamic patterns of gene expression are closely linked with growth and have important consequences for survival.

No abstract

Abundant genomic heterozygosity can be found in wild strains of the budding yeast Saccharomyces cerevisiae isolated from industrial and clinical environments. The extent to which heterozygosity influences the phenotypes of these isolates is not fully understood. One such case is the PE-2/JAY270 strain, a natural hybrid widely adopted by sugarcane bioethanol distilleries for its ability to thrive under harsh biotic and abiotic stresses during industrial scale fermentation, however, it is not known whether or how the heterozygous configuration of the JAY270 genome contributes to its many desirable traits. In this study, we took a step toward exploring this question by conducting an initial functional characterization of JAY270’s heteroalleles. We manipulated the abundance and distribution of heterozygous alleles through inbreeding and targeted uniparental disomy (UPD). Unique combinations of homozygous alleles in each inbred strain revealed wide phenotypic variation for at least two important industrial traits: Heat stress tolerance and competitive growth. Quantitative trait loci analyses allowed the identification of broad genomic regions where genetic polymorphisms potentially impacted these traits, and there was no overlap between the loci associated with each. In addition, we adapted an approach to induce bidirectional UPD of three targeted pairs of chromosomes (IV, XIV, and XV), while heterozygosity was maintained elsewhere in the genome. In most cases UPD led to detectable phenotypic alterations, often in opposite directions between the two homozygous haplotypes in each UPD pair. Our results showed that both widespread and regional homozygosity could uncover cryptic phenotypic variation supported by the heteroalleles residing in the JAY270 genome. Interestingly, we characterized multiple examples of inbred and UPD strains that displayed heat tolerance or competitive growth phenotypes that were superior to their heterozygous parent. However, we propose that homozygosity for those regions may be associated with a decrease in overall fitness in the complex and dynamic distillery environment, and that may have contributed to slowing down the erosion of heterozygosity from the JAY270 genome. This study also laid a foundation for approaches that can be expanded to the identification of specific alleles of interest for industrial applications in this and other hybrid yeast strains.