PHENOS: a high-throughput and flexible tool for microorganism growth phenotyping on solid media

Simple large-scale quantitative phenotyping and antimicrobial susceptibility testing with Q-PHAST

The characterization of antimicrobial susceptibility and other relevant phenotypes in large collections of microbial isolates is a common need across research and clinical microbiology laboratories. Robotization provides unprecedented throughput but involves costs that are prohibitive for the average laboratory. Here, using affordable materials and open-source software, we developed Q-PHAST (Quantitative PHenotyping and Antimicrobial Susceptibility Testing), a unique solution for cost-effective, large-scale phenotyping in a standard microbiology laboratory. Single colonies are grown in a deep 96-well master plate, from which diluted aliquots are used to generate 96 spots on different experimental plates containing solid medium with the substance and concentration of interest. These plates are incubated on inexpensive flatbed scanners that monitor the growth of each spot by obtaining images every 15 min. A simple, python-based software, which can be used via a graphical interface on various operating systems ( https://github.com/Gabaldonlab/Q-PHAST ), analyzes the images to infer growth, fitness (e.g., doubling rate) and susceptibility (e.g., minimum inhibitory concentration) measures. With <120 min of hands-on time per day for three consecutive days, ready-to-use results are obtained and presented in tables or graphs. This solution enables non-experts with limited resources to perform accurate quantitative phenotyping on hundreds of strains in parallel.

Impact of inter-species hybridisation on antifungal drug response in the Saccharomyces genus

BACKGROUND

Antifungal drug resistance presents one of the major concerns for global public health, and hybridization allows the development of high fitness organisms that can better survive in restrictive conditions or in presence of antifungal agents. Hence, understanding how allelic variation can influence antifungal susceptibility in hybrid organisms is important for the development of targeted treatments. Here, we exploited recent advances in multigenerational breeding of hemiascomycete hybrids to study the impact of hybridisation on antifungal resistance and identify quantitative trait loci responsible for the phenotype.

RESULTS

The offspring of Saccharomyces cerevisiae x S. kudriavzevii hybrids were screened in the presence of six antifungal drugs and revealed a broad phenotypic diversity across the progeny. QTL analysis was carried out comparing alleles between pools of high and low fitness offspring, identifying hybrid-specific genetic regions involved in resistance to fluconazole, micafungin and flucytosine. We found both drug specific and pleiotropic regions, including 41 blocks containing genes not previously associated with resistance phenotypes. We identified linked genes that influence the same trait, namely a hybrid specific 'super' QTL, and validated, via reciprocal hemizygosity analysis, two causal genes, BCK2 and DNF1. The co-location of genes with similar phenotypic impact supports the notion of an adaption process that limits the segregation of advantageous alleles via recombination.

CONCLUSIONS

This study demonstrates the value of QTL studies to elucidate the hybrid-specific mechanisms of antifungal susceptibility. We also show that an inter-species hybrid model system in the Saccharomyces background, can help to decipher the trajectory of antifungal drug resistance in pathogenic hybrid lineages.

MiniRead: A simple and inexpensive do-it-yourself device for multiple analyses of micro-organism growth kinetics

Fitness in micro-organisms can be proxied by growth parameters on different media and/or temperatures. This is achieved by measuring optical density at 600 nm using a spectrophotometer, which measures the effect of absorbance and side scattering due to turbidity of cells suspensions. However, when growth kinetics must be monitored in many 96-well plates at the same time, buying several 96-channel spectrophotometers is often beyond budgets. The MiniRead device presented here is a simple and inexpensive do-it-yourself 96-well temperature-controlled turbidimeter designed to measure the interception of white light via absorption or side scattering through liquid culture medium. Turbidity is automatically recorded in each well at regular time intervals for up to several days or weeks. Output tabulated text files are recorded into a micro-SD memory card to be easily transferred to a computer. We propose also an R package which allows (1) to compute the nonlinear calibration curves required to convert raw readings into cell concentration values, and (2) to analyze growth kinetics output files to automatically estimate proxies of growth parameters such as lag time, maximum growth rate, or cell concentration at the plateau.

Multi-Attribute Subset Selection enables prediction of representative phenotypes across microbial populations

The interpretation of complex biological datasets requires the identification of representative variables that describe the data without critical information loss. This is particularly important in the analysis of large phenotypic datasets (phenomics). Here we introduce Multi-Attribute Subset Selection (MASS), an algorithm which separates a matrix of phenotypes (e.g., yield across microbial species and environmental conditions) into predictor and response sets of conditions. Using mixed integer linear programming, MASS expresses the response conditions as a linear combination of the predictor conditions, while simultaneously searching for the optimally descriptive set of predictors. We apply the algorithm to three microbial datasets and identify environmental conditions that predict phenotypes under other conditions, providing biologically interpretable axes for strain discrimination. MASS could be used to reduce the number of experiments needed to identify species or to map their metabolic capabilities. The generality of the algorithm allows addressing subset selection problems in areas beyond biology.

A Normalization Protocol Reduces Edge Effect in High-Throughput Analyses of Hydroxyurea Hypersensitivity in Fission Yeast

Edge effect denotes better growth of microbial organisms situated at the edge of the solid agar media. Although the precise reason underlying edge effect is unresolved, it is generally attributed to greater nutrient availability with less competing neighbors at the edge. Nonetheless, edge effect constitutes an unavoidable confounding factor that results in misinterpretation of cell fitness, especially in high-throughput screening experiments widely employed for genome-wide investigation using microbial gene knockout or mutant libraries. Here, we visualize edge effect in high-throughput high-density pinning arrays and report a normalization approach based on colony growth rate to quantify drug (hydroxyurea)-hypersensitivity in fission yeast strains. This normalization procedure improved the accuracy of fitness measurement by compensating cell growth rate discrepancy at different locations on the plate and reducing false-positive and -negative frequencies. Our work thus provides a simple and coding-free solution for a struggling problem in robotics-based high-throughput screening experiments.

Beyond Saccharomyces pastorianus for modern lager brews: Exploring non-cerevisiae Saccharomyces hybrids with heterotic maltotriose consumption and novel aroma profile

Non-domesticated, wild Saccharomyces yeasts have promising characteristics for beer diversification, particularly when used in the generation of de novo interspecific hybrids. A major motivation for the current work was the question whether attractive novel Saccharomyces interspecific hybrids can be created for the production of exotic lager beers without using the genomic resources of the ale yeast Saccharomyces cerevisiae. Importantly, maltotriose utilization is an essential characteristic typically associated with domesticated ale/lager brewing strains. A high-throughput screening on nearly 200 strains representing all eight species of the Saccharomyces genus was conducted. Three Saccharomyces mikatae strains were able to aerobically grow on maltotriose as the sole carbon source, a trait until recently unidentified for this species. Our screening also confirmed the recently reported maltotriose utilization of the S. jurei strain D5095T. Remarkably, de novo hybrids between a maltotriose-utilizing S. mikatae or S. jurei strain and the maltotriose-negative Saccharomyces eubayanus strain CBS 12357T displayed heterosis and outperformed both parents with regard to aerobically utilizing maltotriose as the sole source of carbon. Indeed, the maximum specific growth rates on this sugar were comparable to the well-known industrial strain, Saccharomyces pastorianus CBS 1513. In lager brewing settings (oxygen-limited), the new hybrids were able to ferment maltose, while maltotriose was not metabolized. Favorable fruity esters were produced, demonstrating that the novel hybrids have the potential to add to the diversity of lager brewing.

Using colony size to measure fitness in Saccharomyces cerevisiae

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

Restoring fertility in yeast hybrids: Breeding and quantitative genetics of beneficial traits

Hybrids between species can harbor a combination of beneficial traits from each parent and may exhibit hybrid vigor, more readily adapting to new harsher environments. Interspecies hybrids are also sterile and therefore an evolutionary dead end unless fertility is restored, usually via auto-polyploidisation events. In the Saccharomyces genus, hybrids are readily found in nature and in industrial settings, where they have adapted to severe fermentative conditions. Due to their hybrid sterility, the development of new commercial yeast strains has so far been primarily conducted via selection methods rather than via further breeding. In this study, we overcame infertility by creating tetraploid intermediates of Saccharomyces interspecies hybrids to allow continuous multigenerational breeding. We incorporated nuclear and mitochondrial genetic diversity within each parental species, allowing for quantitative genetic analysis of traits exhibited by the hybrids and for nuclear-mitochondrial interactions to be assessed. Using pooled F12 generation segregants of different hybrids with extreme phenotype distributions, we identified quantitative trait loci (QTLs) for tolerance to high and low temperatures, high sugar concentration, high ethanol concentration, and acetic acid levels. We identified QTLs that are species specific, that are shared between species, as well as hybrid specific, in which the variants do not exhibit phenotypic differences in the original parental species. Moreover, we could distinguish between mitochondria-type-dependent and -independent traits. This study tackles the complexity of the genetic interactions and traits in hybrid species, bringing hybrids into the realm of full genetic analysis of diploid species, and paves the road for the biotechnological exploitation of yeast biodiversity.

Functional and transcriptional profiling of non-coding RNAs in yeast reveal context-dependent phenotypes and in trans effects on the protein regulatory network

Non-coding RNAs (ncRNAs), including the more recently identified Stable Unannotated Transcripts (SUTs) and Cryptic Unstable Transcripts (CUTs), are increasingly being shown to play pivotal roles in the transcriptional and post-transcriptional regulation of genes in eukaryotes. Here, we carried out a large-scale screening of ncRNAs in Saccharomyces cerevisiae, and provide evidence for SUT and CUT function. Phenotypic data on 372 ncRNA deletion strains in 23 different growth conditions were collected, identifying ncRNAs responsible for significant cellular fitness changes. Transcriptome profiles were assembled for 18 haploid ncRNA deletion mutants and 2 essential ncRNA heterozygous deletants. Guided by the resulting RNA-seq data we analysed the genome-wide dysregulation of protein coding genes and non-coding transcripts. Novel functional ncRNAs, SUT125, SUT126, SUT035 and SUT532 that act in trans by modulating transcription factors were identified. Furthermore, we described the impact of SUTs and CUTs in modulating coding gene expression in response to different environmental conditions, regulating important biological process such as respiration (SUT125, SUT126, SUT035, SUT432), steroid biosynthesis (CUT494, SUT053, SUT468) or rRNA processing (SUT075 and snR30). Overall, these data capture and integrate the regulatory and phenotypic network of ncRNAs and protein-coding genes, providing genome-wide evidence of the impact of ncRNAs on cellular homeostasis.

A quarter of the yeast genome comprises non-coding RNA molecules (ncRNAs), which do not translate into proteins but are involved in the regulation of gene expression. ncRNAs can affect nearby genes by physically interfering with their transcription (cis mode of action), or they interact with DNA, proteins or other RNAs to regulate the expression of distant genes (trans mode of action). Examples of cis-acting ncRNAs have been broadly described, however, genome-wide studies to identify functional trans-acting ncRNAs involved in global gene regulation are still lacking. Here, we used a ncRNA yeast deletion collection to score ncRNA impact on cellular function in different environmental conditions. A group of 20 ncRNA deletion mutants with broad fitness diversity were selected to investigate the ncRNA effect on the protein and ncRNA expression network. We showed a high correlation between altered phenotypes and global transcriptional changes, in an environmental dependent manner. We confirmed the trans acting regulation of ncRNAs in the genome and their role in altering the expression of transcription factors. These findings support the notion of the involvement of ncRNAs in fine tuning cellular expression via regulation of transcription factors, as an advantageous RNA-mediated mechanism that can be fast and cost-effective for the cells.

Genome-Wide Dynamic Evaluation of the UV-Induced DNA Damage Response

Genetic screens in Saccharomyces cerevisiae have allowed for the identification of many genes as sensors or effectors of DNA damage, typically by comparing the fitness of genetic mutants in the presence or absence of DNA-damaging treatments. However, these static screens overlook the dynamic nature of DNA damage response pathways, missing time-dependent or transient effects. Here, we examine gene dependencies in the dynamic response to ultraviolet radiation-induced DNA damage by integrating ultra-high-density arrays of 6144 diploid gene deletion mutants with high-frequency time-lapse imaging. We identify 494 ultraviolet radiation response genes which, in addition to recovering molecular pathways and protein complexes previously annotated to DNA damage repair, include components of the CCR4-NOT complex, tRNA wobble modification, autophagy, and, most unexpectedly, 153 nuclear-encoded mitochondrial genes. Notably, mitochondria-deficient strains present time-dependent insensitivity to ultraviolet radiation, posing impaired mitochondrial function as a protective factor in the ultraviolet radiation response.