Single cell functional genomics reveals the importance of mitochondria in cell-to-cell phenotypic variation

Microfluidic single-cell measurements of oxidative stress as a function of cell cycle position

Single-cell measurements routinely demonstrate high levels of variation between cells, but fewer studies provide insight into the analytical and biological sources of this variation. This is particularly true of chemical cytometry, in which individual cells are lysed and their contents separated, compared to more established single-cell measurements of the genome and transcriptome. To characterize population-level variation and its sources, we analyzed oxidative stress levels in 1278 individual Dictyostelium discoideum cells as a function of exogenous stress level and cell cycle position. Cells were exposed to varying levels of oxidative stress via singlet oxygen generation using the photosensitizer Rose Bengal. Single-cell data reproduced the dose-response observed in ensemble measurements by CE-LIF, superimposed with high levels of heterogeneity. Through experiments and data analysis, we explored possible biological sources of this heterogeneity. No trend was observed between population variation and oxidative stress level, but cell cycle position was a major contributor to heterogeneity in oxidative stress. Cells synchronized to the same stage of cell division were less heterogeneous than unsynchronized cells (RSD of 37-51% vs 93%), and mitotic cells had higher levels of reactive oxygen species than interphase cells. While past research has proposed changes in cell size during the cell cycle as a source of biological noise, the measurements presented here use an internal standard to normalize for effects of cell volume, suggesting a more complex contribution of cell cycle to heterogeneity of oxidative stress.

Cell cycle-dependent Cu uptake explained the heterogenous responses of Chlamydomonas to Cu exposure

Growing evidence suggested that microorganisms exhibited heterogeneous sensitivity to toxicants, but their underlying mechanisms remain largely unknown. The asynchronous cell cycle progression in natural population implies the connection between cell cycle and heterogeneity. Here, the heterogenous responses of Chlamydomonas reinhardtii upon Cu stress were confirmed with the aid of a fluorometric probe for imaging Cu(I), implying the connection with cell cycle. Our results further indicated that the increase of labile Cu(I) was related to the cell division, leading to the fluctuation of labile Cu(I) with diurnal cycle and cell cycle, respectively. However, lack of Cu mainly influenced the cell division. We demonstrated that G2/M phase was the critical stage requiring high Cu quota during cell division. Specifically, algae at G2/M phase required 10-fold of Cu quota compared with that at G1 phase, which was related to the mitochondrial replication. Eventually, the heterogeneous Cu uptake ability of algae at different cell phases led to the heterogeneous responses to Cu exposure. Overall, Cu could influence the cell cycle through mediating the cell division, and in turn algae at different cell phases exhibited different Cu sensitivities. This study firstly uncovered the underlying mechanisms of heterogeneous Cu sensitivity for phytoplankton, which could help to evaluate the potential ecological risks of Cu.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

Proteasome-dependent truncation of the negative heterochromatin regulator Epe1 mediates antifungal resistance

Epe1 histone demethylase restricts H3K9-methylation-dependent heterochromatin, preventing it from spreading over, and silencing, gene-containing regions in fission yeast. External stress induces an adaptive response allowing heterochromatin island formation that confers resistance on surviving wild-type lineages. Here we investigate the mechanism by which Epe1 is regulated in response to stress. Exposure to caffeine or antifungals results in Epe1 ubiquitylation and proteasome-dependent removal of the N-terminal 150 residues from Epe1, generating truncated Epe1 (tEpe1) which accumulates in the cytoplasm. Constitutive tEpe1 expression increases H3K9 methylation over several chromosomal regions, reducing expression of underlying genes and enhancing resistance. Reciprocally, constitutive non-cleavable Epe1 expression decreases resistance. tEpe1-mediated resistance requires a functional JmjC demethylase domain. Moreover, caffeine-induced Epe1-to-tEpe1 cleavage is dependent on an intact cell integrity MAP kinase stress signaling pathway, mutations in which alter resistance. Thus, environmental changes elicit a mechanism that curtails the function of this key epigenetic modifier, allowing heterochromatin to reprogram gene expression, thereby bestowing resistance to some cells within a population. H3K9me-heterochromatin components are conserved in human and crop-plant fungal pathogens for which a limited number of antifungals exist. Our findings reveal how transient heterochromatin-dependent antifungal resistant epimutations develop and thus inform on how they might be countered.

Robustness: linking strain design to viable bioprocesses

Microbial cell factories are becoming increasingly popular for the sustainable production of various chemicals. Metabolic engineering has led to the design of advanced cell factories; however, their long-term yield, titer, and productivity falter when scaled up and subjected to industrial conditions. This limitation arises from a lack of robustness - the ability to maintain a constant phenotype despite the perturbations of such processes. This review describes predictable and stochastic industrial perturbations as well as state-of-the-art technologies to counter process variability. Moreover, we distinguish robustness from tolerance and discuss the potential of single-cell studies for improving system robustness. Finally, we highlight ways of achieving consistent and comparable quantification of robustness that can guide the selection of strains for industrial bioprocesses.

A conserved expression signature predicts growth rate and reveals cell & lineage-specific differences

Isogenic cells cultured together show heterogeneity in their proliferation rate. To determine the differences between fast and slow-proliferating cells, we developed a method to sort cells by proliferation rate, and performed RNA-seq on slow and fast proliferating subpopulations of pluripotent mouse embryonic stem cells (mESCs) and mouse fibroblasts. We found that slowly proliferating mESCs have a more naïve pluripotent character. We identified an evolutionarily conserved proliferation-correlated transcriptomic signature that is common to all eukaryotes: fast cells have higher expression of genes for protein synthesis and protein degradation. This signature accurately predicted growth rate in yeast and cancer cells, and identified lineage-specific proliferation dynamics during development, using C. elegans scRNA-seq data. In contrast, sorting by mitochondria membrane potential revealed a highly cell-type specific mitochondria-state related transcriptome. mESCs with hyperpolarized mitochondria are fast proliferating, while the opposite is true for fibroblasts. The mitochondrial electron transport chain inhibitor antimycin affected slow and fast subpopulations differently. While a major transcriptional-signature associated with cell-to-cell heterogeneity in proliferation is conserved, the metabolic and energetic dependency of cell proliferation is cell-type specific.

By performing RNA sequencing on cells sorted by their proliferation rate, this study identifies a gene expression signature capable of predicting proliferation rates in diverse eukaryotic cell types and species. This signature, applied to single-cell RNA sequencing data from embryos of the roundworm C. elegans, reveals lineage-specific proliferation differences during development. In contrast to the universality of the proliferation signature, mitochondria and metabolism related genes show a high degree of cell-type specificity; mouse pluripotent stem cells (mESCs) and differentiated cells (fibroblasts) exhibit opposite relations between mitochondria state and proliferation. Furthermore, we identified a slow proliferating subpopulation of mESCs with higher expression of pluripotency genes. Finally, we show that fast and slow proliferating subpopulations are differentially sensitive to mitochondria inhibitory drugs in different cell types.

Highlights:

A FACS-based method to determine the transcriptomes of fast and slow proliferating subpopulations.

A universal proliferation-correlated transcriptional signature indicates high protein synthesis and degradation in fast proliferating cells across cell types and species.

Applied to scRNA-seq, the expression signature predicts the global proliferation slowdown during C. elegans development.

Mitochondria membrane potential predicts proliferation rate in a cell-type specific manner, with ETC complex III inhibitor having distinct effects on fibroblasts vs mESCs.

A rapidly reversible mutation generates subclonal genetic diversity and unstable drug resistance

Most genetic changes have negligible reversion rates. As most mutations that confer resistance to an adverse condition (e.g., drug treatment) also confer a growth defect in its absence, it is challenging for cells to genetically adapt to transient environmental changes. Here, we identify a set of rapidly reversible drug-resistance mutations in Schizosaccharomyces pombe that are caused by microhomology-mediated tandem duplication (MTD) and reversion back to the wild-type sequence. Using 10,000× coverage whole-genome sequencing, we identify nearly 6,000 subclonal MTDs in a single clonal population and determine, using machine learning, how MTD frequency is encoded in the genome. We find that sequences with the highest-predicted MTD rates tend to generate insertions that maintain the correct reading frame, suggesting that MTD formation has shaped the evolution of coding sequences. Our study reveals a common mechanism of reversible genetic variation that is beneficial for adaptation to environmental fluctuations and facilitates evolutionary divergence.

Insights on the Control of Yeast Single-Cell Growth Variability by Members of the Trehalose Phosphate Synthase (TPS) Complex

Single-cell variability of growth is a biological phenomenon that has attracted growing interest in recent years. Important progress has been made in the knowledge of the origin of cell-to-cell heterogeneity of growth, especially in microbial cells. To better understand the origins of such heterogeneity at the single-cell level, we developed a new methodological pipeline that coupled cytometry-based cell sorting with automatized microscopy and image analysis to score the growth rate of thousands of single cells. This allowed investigating the influence of the initial amount of proteins of interest on the subsequent growth of the microcolony. As a preliminary step to validate this experimental setup, we referred to previous findings in yeast where the expression level of Tsl1, a member of the Trehalose Phosphate Synthase (TPS) complex, negatively correlated with cell division rate. We unfortunately could not find any influence of the initial TSL1 expression level on the growth rate of the microcolonies. We also analyzed the effect of the natural variations of trehalose-6-phosphate synthase (TPS1) expression on cell-to-cell growth heterogeneity, but we did not find any correlation. However, due to the already known altered growth of the tps1Δ mutants, we tested this strain at the single-cell level on a permissive carbon source. This mutant showed an outstanding lack of reproducibility of growth rate distributions as compared to the wild-type strain, with variable proportions of non-growing cells between cultivations and more heterogeneous microcolonies in terms of individual growth rates. Interestingly, this variable behavior at the single-cell level was reminiscent to the high variability that is also stochastically suffered at the population level when cultivating this tps1Δ strain, even when using controlled bioreactors.

Intrinsic Exercise Capacity and Mitochondrial DNA Lead to Opposing Vascular-Associated Risks

Exercise capacity is a strong predictor of all-cause morbidity and mortality in humans. However, the associated hemodynamic traits that link this valuable indicator to its subsequent disease risks are numerable. Additionally, exercise capacity has a substantial heritable component and genome-wide screening indicates a vast amount of nuclear and mitochondrial DNA (mtDNA) markers are significantly associated with traits of physical performance. A long-term selection experiment in rats confirms a divide for cardiovascular risks between low- and high-capacity runners (LCR and HCR, respectively), equipping us with a preclinical animal model to uncover new mechanisms. Here, we evaluated the LCR and HCR rat model system for differences in vascular function at the arterial resistance level. Consistent with the known divide between health and disease, we observed that LCR rats present with resistance artery and perivascular adipose tissue dysfunction compared to HCR rats that mimic qualities important for health, including improved vascular relaxation. Uniquely, we show by generating conplastic strains, which LCR males with mtDNA of female HCR (LCR-mt^HCR/Tol) present with improved vascular function. Conversely, HCR-mt^LCR/Tol rats displayed indices for cardiac dysfunction. The outcome of this study suggests that the interplay between the nuclear genome and the maternally inherited mitochondrial genome with high intrinsic exercise capacity is a significant factor for improved vascular physiology, and animal models developed on an interaction between nuclear and mtDNA are valuable new tools for probing vascular risk factors in the offspring.

Towards Accurate Genotype-Phenotype Correlations in the CYP2D6 Gene

Establishing accurate and large-scale genotype-phenotype correlations and predictions of individual response to pharmacological treatments are two of the holy grails of Personalized Medicine. These tasks are challenging and require an integrated knowledge of the complex processes that regulate gene expression and, ultimately, protein functionality in vivo, the effects of mutations/polymorphisms and the different sources of interindividual phenotypic variability. A remarkable example of our advances in these challenging tasks is the highly polymorphic CYP2D6 gene, which encodes a cytochrome P450 enzyme involved in the metabolization of many of the most marketed drugs (including SARS-Cov-2 therapies such as hydroxychloroquine). Since the introduction of simple activity scores (AS) over 10 years ago, its ability to establish genotype-phenotype correlations on the drug metabolizing capacity of this enzyme in human population has provided lessons that will help to improve this type of score for this, and likely many other human genes and proteins. Multidisciplinary research emerges as the best approach to incorporate additional concepts to refine and improve such functional/activity scores for the CYP2D6 gene, as well as for many other human genes associated with simple and complex genetic diseases.

A Landscape of Murine Long Non-Coding RNAs Reveals the Leading Transcriptome Alterations in Adipose Tissue during Aging

Aging is an inevitable process that involves profound physiological changes. Long non-coding RNAs (lncRNAs) are emerging as important regulators in various biological processes but are not systemically studied in aging. To provide an organism-wide lncRNA landscape during aging, we conduct comprehensive RNA sequencing (RNA-seq) analyses across the mouse lifespan. Of the 1,675 aging-regulated lncRNAs (AR-lncRNAs) identified, the majority are connected to inflammation-related biological pathways. AR-lncRNAs exhibit high tissue specificity; conversely, those with higher tissue specificity are preferentially regulated during aging. White adipose tissue (WAT) displays the highest number of AR-lncRNAs and develops the most dynamic crosstalk between AR-lncRNA and AR-mRNA during aging. An adipose-enriched AR-lncRNA, lnc-adipoAR1, is negatively correlated with aging, and knocking it down inhibits adipogenesis, phenocopying the compromised adipogenic capacity of aged fat. Our works together reveal AR-lncRNAs as essential components in aging and suggest that although each tissue ages in a distinct manner, WAT is a leading contributor to aging-related health decline.

Diversity and Postzygotic Evolution of the Mitochondrial Genome in Hybrids of Saccharomyces Species Isolated by Double Sterility Barrier

Eukaryotic species are reproductively isolated by sterility barriers that prevent interspecies fertilization (prezygotic sterility barrier) or the fertilization results in infertile offspring (postzygotic sterility barrier). The Saccharomyces species are isolated by postzygotic sterility barriers. Their allodiploid hybrids form no viable gametes (ascospores) and the viable ascospores of the allotetraploids cannot fertilize (conjugate). Our previous work revealed that this mechanism of reproductive isolation differs from those operating in plants and animals and we designated it double sterility barrier (the failure of homeologous chromosomes to pair and the repression of mating by mating-type heterozygosity). Other studies implicated nucleo-mitochondrial incompatibilities in the sterility of the Saccharomyces hybrids, a mechanism assumed to play a central role in the reproductive isolation of animal species. In this project the mitochondrial genomes of 50 cevarum (S. cerevisiae × S. uvarum) hybrids were analyzed. 62% had S. cerevisiae mitotypes, 4% had S. uvarum mitotypes, and 34% had recombinant mitotypes. All but one hybrid formed viable spores indicating that they had genomes larger than allodiploid. Most of these spores were sterile (no sporulation in the clone of vegetative descendants; a feature characteristic of allodiploids). But regardless of their mitotypes, most hybrids could also form fertile alloaneuploid spore clones at low frequencies upon the loss of the MAT-carrying chromosome of the S. uvarum subgenome during meiosis. Hence, the cevarum alloploid nuclear genome is compatible with both parental mitochondrial genomes as well as with their recombinants, and the sterility of the hybrids is maintained by the double sterility barrier (determined in the nuclear genome) rather than by nucleo-mitochondrial incompatibilities. During allotetraploid sporulation both the nuclear and the mitochondrial genomes of the hybrids could segregate but no correlation was observed between the sterility or the fertility of the spore clones and their mitotypes. Nucleo-mitochondrial incompatibility was manifested as respiration deficiency in certain meiotic segregants. As respiration is required for meiosis-sporulation but not for fertilization (conjugation), these segregants were deficient only in sporulation. Thus, the nucleo-mitochondrial incompatibility affects the sexual processes only indirectly through the inactivation of respiration and causes only partial sterility in certain segregant spore clones.

Allele-specific single-cell RNA sequencing reveals different architectures of intrinsic and extrinsic gene expression noises

Gene expression noise refers to the variation of the expression level of a gene among isogenic cells in the same environment, and has two sources: extrinsic noise arising from the disparity of the cell state and intrinsic noise arising from the stochastic process of gene expression in the same cell state. Due to the low throughput of the existing method for measuring the two noise components, the architectures of intrinsic and extrinsic expression noises remain elusive. Using allele-specific single-cell RNA sequencing, we here estimate the two noise components of 3975 genes in mouse fibroblast cells. Our analyses verify predicted influences of several factors such as the TATA-box and microRNA targeting on intrinsic or extrinsic noises and reveal gene function-associated noise trends implicating the action of natural selection. These findings unravel differential regulations, optimizations, and biological consequences of intrinsic and extrinsic noises and can aid the construction of desired synthetic circuits.

Single cell functional genomics reveals the importance of mitochondria in cell-to-cell phenotypic variation

Mutations frequently have outcomes that differ across individuals, even when these individuals are genetically identical and share a common environment. Moreover, individual microbial and mammalian cells can vary substantially in their proliferation rates, stress tolerance, and drug resistance, with important implications for the treatment of infections and cancer. To investigate the causes of cell-to-cell variation in proliferation, we used a high-throughput automated microscopy assay to quantify the impact of deleting >1500 genes in yeast. Mutations affecting mitochondria were particularly variable in their outcome. In both mutant and wild-type cells mitochondrial membrane potential - but not amount - varied substantially across individual cells and predicted cell-to-cell variation in proliferation, mutation outcome, stress tolerance, and resistance to a clinically used anti-fungal drug. These results suggest an important role for cell-to-cell variation in the state of an organelle in single cell phenotypic variation.