Genetic background effects in quantitative genetics: gene-by-system interactions

Candida albicans natural diversity: a resource to dissect fungal commensalism and pathogenesis

Candida albicans is a ubiquitous fungus of humans. It is not only a component of the oral and intestinal microbiota of most healthy adults but also a major cause of mucosal disorders and life-threatening disseminated infections. Until recently, research on the biology and pathogenesis of the fungus was largely based on a single clinical isolate. We review investigations that have started to dissect a diverse set of C. albicans strains. Using different approaches to leverage the species' phenotypic and/or genetic diversity, these studies illuminate the wide range of interactions between fungus and host. While connecting genetic variants to phenotypes of interest remains challenging, research on C. albicans' natural diversity is central to understand fungal commensalism and pathogenesis.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

Functional Dichotomy for a Hyphal Repressor in Candida albicans

Nrg1 is a repressor of hypha formation and hypha-associated gene expression in the fungal pathogen Candida albicans. It has been well studied in the genetic background of the type strain SC5314. Here, we tested Nrg1 function in four other diverse clinical isolates through an analysis of nrg1Δ/Δ mutants, with SC5314 included as a control. In three strains, nrg1Δ/Δ mutants unexpectedly produced aberrant hyphae under inducing conditions, as assayed by microscopic observation and endothelial cell damage. The nrg1Δ/Δ mutant of strain P57055 had the most severe defect. We examined gene expression features under hypha-inducing conditions by RNA-sequencing (RNA-Seq) for the SC5314 and P57055 backgrounds. The SC5314 nrg1Δ/Δ mutant expressed six hypha-associated genes at reduced levels compared with wild-type SC5314. The P57055 nrg1Δ/Δ mutant expressed 17 hypha-associated genes at reduced levels compared with wild-type P57055, including IRF1, RAS2, and ECE1. These findings indicate that Nrg1 has a positive role in hypha-associated gene expression and that this role is magnified in strain P57055. Remarkably, the same hypha-associated genes affected by the nrg1Δ/Δ mutation in strain P57055 were also naturally expressed at lower levels in wild-type P57055 than those in wild-type SC5314. Our results suggest that strain P57055 is defective in a pathway that acts in parallel with Nrg1 to upregulate the expression of several hypha-associated genes. IMPORTANCE Hypha formation is a central virulence trait of the fungal pathogen Candida albicans. Control of hypha formation has been studied in detail in the type strain but not in other diverse C. albicans clinical isolates. Here, we show that the hyphal repressor Nrg1 has an unexpected positive role in hypha formation and hypha-associated gene expression, as revealed by the sensitized P57055 strain background. Our findings indicate that reliance on a single type strain limits understanding of gene function and illustrate that strain diversity is a valuable resource for C. albicans molecular genetic analysis.

Tra1 controls the transcriptional landscape of the aging cell

Gene expression undergoes considerable changes during the aging process. The mechanisms regulating the transcriptional response to cellular aging remain poorly understood. Here, we employ the budding yeast Saccharomyces cerevisiae to better understand how organisms adapt their transcriptome to promote longevity. Chronological lifespan assays in yeast measure the survival of nondividing cells at stationary phase over time, providing insights into the aging process of postmitotic cells. Tra1 is an essential component of both the yeast Spt-Ada-Gcn5 acetyltransferase/Spt-Ada-Gcn5 acetyltransferase-like and nucleosome acetyltransferase of H4 complexes, where it recruits these complexes to acetylate histones at targeted promoters. Importantly, Tra1 regulates the transcriptional response to multiple stresses. To evaluate the role of Tra1 in chronological aging, we took advantage of a previously characterized mutant allele that carries mutations in the TRA1 PI3K domain (tra1Q3). We found that loss of functions associated with tra1Q3 sensitizes cells to growth media acidification and shortens lifespan. Transcriptional profiling reveals that genes differentially regulated by Tra1 during the aging process are enriched for components of the response to stress. Notably, expression of catalases (CTA1, CTT1) involved in hydrogen peroxide detoxification decreases in chronologically aged tra1Q3 cells. Consequently, they display increased sensitivity to oxidative stress. tra1Q3 cells are unable to grow on glycerol indicating a defect in mitochondria function. Aged tra1Q3 cells also display reduced expression of peroxisomal genes, exhibit decreased numbers of peroxisomes, and cannot grow on media containing oleate. Thus, Tra1 emerges as an important regulator of longevity in yeast via multiple mechanisms.

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions

The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

Parallel evolution between genomic segments of seasonal human influenza viruses reveals RNA-RNA relationships

The influenza A virus (IAV) genome consists of eight negative-sense viral RNA (vRNA) segments that are selectively assembled into progeny virus particles through RNA-RNA interactions. To explore putative intersegmental RNA-RNA relationships, we quantified similarity between phylogenetic trees comprising each vRNA segment from seasonal human IAV. Intersegmental tree similarity differed between subtype and lineage. While intersegmental relationships were largely conserved over time in H3N2 viruses, they diverged in H1N1 strains isolated before and after the 2009 pandemic. Surprisingly, intersegmental relationships were not driven solely by protein sequence, suggesting that IAV evolution could also be driven by RNA-RNA interactions. Finally, we used confocal microscopy to determine that colocalization of highly coevolved vRNA segments is enriched over other assembly intermediates at the nuclear periphery during productive viral infection. This study illustrates how putative RNA interactions underlying selective assembly of IAV can be interrogated with phylogenetics.

The viruses responsible for influenza evolve rapidly during infection. Changes typically emerge in two key ways: through random mutations in the genetic sequence of the virus, or by reassortment. Reassortment can occur when two or more strains infect the same cell. Once in a cell, viral particles ‘open up’ to release their genetic material so it can make copies of itself using the cell’s machinery. The new copies of the genetic material of the virus are used to make new viral particles, which then envelop the genetic material and are released from the cell to infect other cells. If several strains of a virus infect the same cell, a new viral particle may pick up genetic segments from each of the infecting strains, creating a new strain via reassortment.

Several factors are known to affect the success of the reassortment process. For example, if the new strain acquires a genetic defect that hinders its replication cycle, it is likely to die out quickly. Other times, this trading of genetic information can create a strain that is more resistant to the human immune system, allowing it to sweep across the globe and cause a deadly pandemic. However, a key part of the reassortment process that still remains unclear is how genome segments from two different influenza strains recognize each other before merging together to create hybrid daughter viruses.

To explore this further, Jones et al. used a technique called fluorescence microscopy. They found that genome segments that evolved along similar paths were more likely to cluster in the same area inside infected cells, and therefore, more likely to be reassorted together into a new strain during assembly of daughter viruses. This suggests that assembly may guide the evolutionary path taken by individual genomic segments. Jones et al. also looked at the evolution of different genome segments collected from patients suffering from seasonal influenza, and found that these segments had a distinct evolutionary path to those in pandemic-causing strains.

This research provides new insights into the role of reassortment in the evolution of influenza viruses during infection. In particular, it suggests that how the genome segments interact with one another may have a previously unknown and important role in guiding this evolution. These insights could be used to predict future reassortment events based on evolutionary relationships between influenza virus genomic segments, and may in the future be used as part of risk assessment tools to predict the emergence of new pandemic strains.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

Stress modulation as a means to improve yeasts for lignocellulose bioconversion

The second-generation (2G) fermentation environment for lignocellulose conversion presents unique challenges to the fermentative organism that do not necessarily exist in other industrial fermentations. While extreme osmotic, heat, and nutrient starvation stresses are observed in sugar- and starch-based fermentation environments, additional pre-treatment-derived inhibitor stress, potentially exacerbated by stresses such as pH and product tolerance, exist in the 2G environment. Furthermore, in a consolidated bioprocessing (CBP) context, the organism is also challenged to secrete enzymes that may themselves lead to unfolded protein response and other stresses. This review will discuss responses of the yeast Saccharomyces cerevisiae to 2G-specific stresses and stress modulation strategies that can be followed to improve yeasts for this application. We also explore published -omics data and discuss relevant rational engineering, reverse engineering, and adaptation strategies, with the view of identifying genes or alleles that will make positive contributions to the overall robustness of 2G industrial strains. KEYPOINTS: • Stress tolerance is a key driver to successful application of yeast strains in biorefineries. • A wealth of data regarding stress responses has been gained through omics studies. • Integration of this knowledge could inform engineering of fit for purpose strains.

Homeostatic plasticity fails at the intersection of autism-gene mutations and a novel class of common genetic modifiers

We identify a set of common phenotypic modifiers that interact with five independent autism gene orthologs (RIMS1, CHD8, CHD2, WDFY3, ASH1L) causing a common failure of presynaptic homeostatic plasticity (PHP) in Drosophila. Heterozygous null mutations in each autism gene are demonstrated to have normal baseline neurotransmission and PHP. However, PHP is sensitized and rendered prone to failure. A subsequent electrophysiology-based genetic screen identifies the first known heterozygous mutations that commonly genetically interact with multiple ASD gene orthologs, causing PHP to fail. Two phenotypic modifiers identified in the screen, PDPK1 and PPP2R5D, are characterized. Finally, transcriptomic, ultrastructural and electrophysiological analyses define one mechanism by which PHP fails; an unexpected, maladaptive up-regulation of CREG, a conserved, neuronally expressed, stress response gene and a novel repressor of PHP. Thus, we define a novel genetic landscape by which diverse, unrelated autism risk genes may converge to commonly affect the robustness of synaptic transmission.

Circuit diversification in a biofilm regulatory network

Genotype-phenotype relationships can vary extensively among members of a species. One cause of this variation is circuit diversification, the alteration of gene regulatory relationships among members of a species. Circuit diversification is thought to be a starting point for the circuit divergence or rewiring that occurs during speciation. How widespread is circuit diversification? Here we address this question with the fungal pathogen Candida albicans, which forms biofilms rich in distinctive hyphal cells as a prelude to infection. Our understanding of the biofilm/hyphal regulatory network comes primarily from studies of one clinical isolate, strain SC5314, and its marked derivatives. We used CRISPR-based methods to create mutations of four key biofilm transcription factor genes–BCR1, UME6, BRG1, and EFG1 –in SC5314 and four additional clinical isolates. Phenotypic analysis revealed that mutations in BCR1 or UME6 have variable impact across strains, while mutations in BRG1 or EFG1 had uniformly severe impact. Gene expression, sampled with Nanostring probes and examined comprehensively for EFG1 via RNA-Seq, indicates that regulatory relationships are highly variable among isolates. Our results suggest that genotype-phenotype relationships vary in this strain panel in part because of differences in control of BRG1 by BCR1, a hypothesis that is supported through engineered constitutive expression of BRG1. Overall, the data show that circuit diversification is the rule, not the exception, in this biofilm/hyphal regulatory network.

Much of what we know about microbial pathogens is derived from in-depth analysis of one or a few standard laboratory strains. This statement is especially true for the fungal pathogen Candida albicans, because most studies have centered on strain SC5314 and its genetically marked derivatives. Here we examine the functional impact of mutations of four key biofilm regulators across five different clinical isolates. We observe that functional impact of the mutations, based on biological phenotypes and gene expression effects, varies extensively among the isolates. Our results support the idea that gene function should be validated with multiple strain isolates. In addition, our results indicate that a core regulatory network, which comprises regulatory relationships common to multiple isolates, may be enriched for functionally relevant genes.