Publications
2024
-
(2024) Nano Today. 54, 102084. Abstract
Two-dimensional (2D) transition metal dichalcogenides (TMDs), such as molybdenum disulfide (MoS2) and tungsten disulfide (WS2), have attracted considerable interest in biomedicine due to their unique combination of physicochemical properties. The effect of nanomaterials on immune cells and their biodistribution are critical aspects of their clinical translation. However, understanding the interactions of these emerging 2D nanomaterials with the complex pool of peripheral blood mononuclear cells (PBMCs) at the single-cell level and their in situ localization in the main organs still needs to be discovered, preventing their translation in medical settings. Here, we report in-depth immune profiling of water-based and defect-free 2D formulations of MoS2 and WS2 through the simultaneous label-free tracking of their immune cell interactions, both ex vivo in human PBMCs and in vivo in mice by high-dimensional analytical approaches, as well as their biodistribution. For comparison, we studied graphene, the hitherto most explored 2D material for biomedical applications. First, we assessed the impact at the protein and gene level by multiplex protein arrays and RNA sequencing, demonstrating a very modest effect of MoS2 and WS2 on immune cell functionality compared to graphene. Then, a single-cell view of the effects of MoS2 and WS2 on 16 primary human immune cell types in terms of viability and functionality was obtained by single-cell mass cytometry by time-of-flight (CyTOF). We explored over 30 markers looking at multiple cell parameters. Finally, we present evidence that MoS2 and WS2 are visible, without the need for labeling, at the single-cell and tissue level by CyTOF, imaging mass cytometry, and multiplexed ion beam imaging by time-of-flight (MIBI-TOF). In particular, MoS2 and WS2 could be detected in the molybdenum (95Mo) and tungsten (180–186W) channels, respectively, which are not used for commercial mass cytometry tags, allowing for the simultaneous interrogation of a wide variety of biological parameters ex vivo and in vivo following intravenous administration of the TMDs. Indeed, we demonstrated the accumulation of TMDs in the main organs by MIBI-TOF and they could also be identified in specific immune cell subsets by CyTOF. Among the two TMDs studied, WS2 exhibited the highest brightness and signal intensity in all the cell subpopulations and tissues analyzed. In conclusion, we identified TMDs as immune-compatible nanoplatforms, traceable at the single-cell and tissue (sub-organ) levels, thus opening up new perspectives for their exploration in biomedicine.
2023
-
(2023) Proceedings of the National Academy of Sciences. 120, 44, e230851112. Abstract
The immune system is a complex network of cells with critical functions in health and disease. However, a comprehensive census of the cells comprising the immune system is lacking. Here, we estimated the abundance of the primary immune cell types throughout all tissues in the human body. We conducted a literature survey and integrated data from multiplexed imaging and methylome-based deconvolution. We also considered cellular mass to determine the distribution of immune cells in terms of both number and total mass. Our results indicate that the immune system of a reference 73 kg man consists of 1.8 x 1012 cells (95% CI 1.5-2.3 x 1012), weighing 1.2 kg (95% CI 0.8-1.9). Lymphocytes constitute 40% of the total number of immune cells and 15% of the mass and are mainly located in the lymph nodes and spleen. Neutrophils account for similar proportions of both the number and total mass of immune cells, with most neutrophils residing in the bone marrow. Macrophages, present in most tissues, account for 10% of immune cells but contribute nearly 50% of the total cellular mass due to their large size. The quantification of immune cells within the human body presented here can serve to understand the immune function better and facilitate quantitative modeling of this vital system.
-
(2023) Nature Communications. 14, 4302. Abstract
Multiplexed imaging enables measurement of multiple proteins in situ, offering an unprecedented opportunity to chart various cell types and states in tissues. However, cell classification, the task of identifying the type of individual cells, remains challenging, labor-intensive, and limiting to throughput. Here, we present CellSighter, a deep-learning based pipeline to accelerate cell classification in multiplexed images. Given a small training set of expert-labeled images, CellSighter outputs the label probabilities for all cells in new images. CellSighter achieves over 80% accuracy for major cell types across imaging platforms, which approaches inter-observer concordance. Ablation studies and simulations show that CellSighter is able to generalize its training data and learn features of protein expression levels, as well as spatial features such as subcellular expression patterns. CellSighter’s design reduces overfitting, and it can be trained with only thousands or even hundreds of labeled examples. CellSighter also outputs a prediction confidence, allowing downstream experts control over the results. Altogether, CellSighter drastically reduces hands-on time for cell classification in multiplexed images, while improving accuracy and consistency across datasets.
-
(2023) Nature. 619, 7970, p. 595-605 Abstract
Beginning in the first trimester, fetally derived extravillous trophoblasts (EVTs) invade the uterus and remodel its spiral arteries, transforming them into large, dilated blood vessels. Several mechanisms have been proposed to explain how EVTs coordinate with the maternal decidua to promote a tissue microenvironment conducive to spiral artery remodelling (SAR)1–3. However, it remains a matter of debate regarding which immune and stromal cells participate in these interactions and how this evolves with respect to gestational age. Here we used a multiomics approach, combining the strengths of spatial proteomics and transcriptomics, to construct a spatiotemporal atlas of the human maternal–fetal interface in the first half of pregnancy. We used multiplexed ion beam imaging by time-of-flight and a 37-plex antibody panel to analyse around 500,000 cells and 588 arteries within intact decidua from 66 individuals between 6 and 20 weeks of gestation, integrating this dataset with co-registered transcriptomics profiles. Gestational age substantially influenced the frequency of maternal immune and stromal cells, with tolerogenic subsets expressing CD206, CD163, TIM-3, galectin-9 and IDO-1 becoming increasingly enriched and colocalized at later time points. By contrast, SAR progression preferentially correlated with EVT invasion and was transcriptionally defined by 78 gene ontology pathways exhibiting distinct monotonic and biphasic trends. Last, we developed an integrated model of SAR whereby invasion is accompanied by the upregulation of pro-angiogenic, immunoregulatory EVT programmes that promote interactions with the vascular endothelium while avoiding the activation of maternal immune cells.
-
(2023) Cancer Cell. 41, 3, p. 404-420 Abstract
The tumor microenvironment (TME) is composed of many different cellular and acellular components that together drive tumor growth, invasion, metastasis, and response to therapies. Increasing realization of the significance of the TME in cancer biology has shifted cancer research from a cancer-centric model to one that considers the TME as a whole. Recent technological advancements in spatial profiling methodologies provide a systematic view and illuminate the physical localization of the components of the TME. In this review, we provide an overview of major spatial profiling technologies. We present the types of information that can be extracted from these data and describe their applications, findings and challenges in cancer research. Finally, we provide a future perspective of how spatial profiling could be integrated into cancer research to improve patient diagnosis, prognosis, stratification to treatment and development of novel therapeutics.
2022
-
(2022) Advanced materials (Weinheim). 34, 45, 2205154. Abstract
There is a critical unmet need to detect and image 2D materials within single cells and tissues while surveying a high degree of information from single cells. Here, a versatile multiplexed label‐free single‐cell detection strategy is proposed based on single‐cell mass cytometry by time‐of‐flight (CyTOF) and ion‐beam imaging by time‐of‐flight (MIBI‐TOF). This strategy, “Label‐free sINgle‐cell tracKing of 2D matErials by mass cytometry and MIBI‐TOF Design” (LINKED), enables nanomaterial detection and simultaneous measurement of multiple cell and tissue features. As a proof of concept, a set of 2D materials, transition metal carbides, nitrides, and carbonitrides (MXenes), is selected to ensure mass detection within the cytometry range while avoiding overlap with more than 70 currently available tags, each able to survey multiple biological parameters. First, their detection and quantification in 15 primary human immune cell subpopulations are demonstrated. Together with the detection, mass cytometry is used to capture several biological aspects of MXenes, such as their biocompatibility and cytokine production after their uptake. Through enzymatic labeling, MXenes’ mediation of cell–cell interactions is simultaneously evaluated. In vivo biodistribution experiments using a mixture of MXenes in mice confirm the versatility of the detection strategy and reveal MXene accumulation in the liver, blood, spleen, lungs, and relative immune cell subtypes. Finally, MIBI‐TOF is applied to detect MXenes in different organs revealing their spatial distribution. The label‐free detection of 2D materials by mass cytometry at the single‐cell level, on multiple cell subpopulations and in multiple organs simultaneously, will enable exciting new opportunities in biomedicine.
“Label‐free sINgle‐cell tracKing of 2D matErials by mass cytometry and MIBI‐TOF Design” (LINKED) is a multiplexed label‐free single‐cell detection strategy enabling nanomaterial tracking and simultaneous measurement of multiple cell and tissue features. MXenes are designed to be detected within the mass cytometry range. Their detection and quantification are demonstrated in 15 human immune cell subpopulations and in different organs. -
(2022) Molecular Systems Biology. 18, 5, e10726. Abstract
Twenty years ago, molecular biology transitioned from predominantly studying genes as isolated elements to viewing them as part of complex modules, giving rise to the field of systems biology. This transition was made possible by technological advances that allowed to simultaneously measure the expression levels of thousands of genes in a single experiment and drove a shift toward analyses identifying gene sets, modules, and pathways involved in a biological process of interest. Today we are excitingly facing a similar turning point in cell biology, where single-cell technologies have enabled us to approach cells as cellular modules.
-
(2022) Nature biotechnology. 40, 4, p. 555-565 Abstract
A principal challenge in the analysis of tissue imaging data is cell segmentation—the task of identifying the precise boundary of every cell in an image. To address this problem we constructed TissueNet, a dataset for training segmentation models that contains more than 1 million manually labeled cells, an order of magnitude more than all previously published segmentation training datasets. We used TissueNet to train Mesmer, a deep-learning-enabled segmentation algorithm. We demonstrated that Mesmer is more accurate than previous methods, generalizes to the full diversity of tissue types and imaging platforms in TissueNet, and achieves human-level performance. Mesmer enabled the automated extraction of key cellular features, such as subcellular localization of protein signal, which was challenging with previous approaches. We then adapted Mesmer to harness cell lineage information in highly multiplexed datasets and used this enhanced version to quantify cell morphology changes during human gestation. All code, data and models are released as a community resource.
-
(2022) Nature Immunology. 23, 2, p. 318-329 Abstract
Tuberculosis (TB) in humans is characterized by formation of immune-rich granulomas in infected tissues, the architecture and composition of which are thought to affect disease outcome. However, our understanding of the spatial relationships that control human granulomas is limited. Here, we used multiplexed ion beam imaging by time of flight (MIBI-TOF) to image 37 proteins in tissues from patients with active TB. We constructed a comprehensive atlas that maps 19 cell subsets across 8 spatial microenvironments. This atlas shows an IFN-γ-depleted microenvironment enriched for TGF-β, regulatory T cells and IDO1+ PD-L1+ myeloid cells. In a further transcriptomic meta-analysis of peripheral blood from patients with TB, immunoregulatory trends mirror those identified by granuloma imaging. Notably, PD-L1 expression is associated with progression to active TB and treatment response. These data indicate that in TB granulomas, there are local spatially coordinated immunoregulatory programs with systemic manifestations that define active TB.
-
(2022) Cell. 185, 2, p. 299-310.e18 Abstract
Ductal carcinoma in situ (DCIS) is a pre-invasive lesion that is thought to be a precursor to invasive breast cancer (IBC). To understand the changes in the tumor microenvironment (TME) accompanying transition to IBC, we used multiplexed ion beam imaging by time of flight (MIBI-TOF) and a 37-plex antibody staining panel to interrogate 79 clinically annotated surgical resections using machine learning tools for cell segmentation, pixel-based clustering, and object morphometrics. Comparison of normal breast with patient matched DCIS and IBC revealed coordinated transitions between four TME states that were delineated based on the location and function of myoepithelium, fibroblasts, and immune cells. Surprisingly, myoepithelial disruption was more advanced in DCIS patients that did not develop IBC, suggesting this process could be protective against recurrence. Taken together, this HTAN Breast PreCancer Atlas study offers insight into drivers of IBC relapse and emphasizes the importance of the TME in regulating these processes.
-
(2022) Single-Cell Protein Analysis. p. 147-156 Abstract
Multiplexed Ion Beam Imaging by Time of Flight (MIBI-TOF) enables high-dimensional imaging in situ of clinical specimens at single-cell resolution. In MIBI-TOF, tissue sections are stained with dozens of metal-labeled antibodies, whose abundance and location are read by secondary ionization mass spectrometry. The result is a multi-dimensional image, depicting sub-cellular expression and localization for dozens of distinct proteins in situ. Here, we describe the staining and imaging procedures of a MIBI-TOF experiment.
2021
-
(2021) Annals of the New York Academy of Sciences. 1506, 1, p. 74-97 Abstract
Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state—transcriptional, epigenetic, and other characteristics—can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17–19, 2021, experts in single cell biology met virtually for the Keystone eSymposium “Single Cell Biology” to discuss advances both in single cell applications and technologies.
-
(2021) Cancer Cell. 39, 8, p. 1015-1017 Abstract
Tumor heterogeneity was traditionally considered in the genetic terms, but it has now been broadened into many more facets. These facets represent a challenge in our understanding of cancer etiology but also provide opportunity for us to understand prognosis and therapy response.
-
(2021) Communications Biology. 4, 852. Abstract
Triple-negative breast cancer (TNBC), the poorest-prognosis breast cancer subtype, lacks clinically approved biomarkers for patient risk stratification, treatment management, and immunotherapies. Prior literature has shown that interrogation of the tumor-immune microenvironment (TIME) may be a promising approach for the discovery of novel biomarkers that can fill these gaps. Recent developments in high-dimensional tissue imaging technology, such as multiplexed ion beam imaging (MIBI), provide spatial context to protein expression in the TIME, opening doors for in-depth characterization of cellular processes. We developed a computational pipeline for the robust examination of the TIME using MIBI. We discover that profiling the functional proteins involved in cell-to-cell interactions in the TIME predicts recurrence and overall survival in TNBC. The interactions between CD45RO and Beta Catenin and CD45RO and HLA-DR were the most relevant for patient stratification. We demonstrated the clinical relevance of the immunoregulatory proteins PD-1, PD-L1, IDO, and Lag3 by tying their interactions to recurrence and survival. Multivariate analysis revealed that our methods provide additional prognostic information compared to clinical variables. Our novel computational pipeline produces interpretable results, and is generalizable to other cancer types.
-
(2021) Frontiers in Immunology. 12, 652631. Abstract
Multiplex imaging technologies are now routinely capable of measuring more than 40 antibody-labeled parameters in single cells. However, lateral spillage of signals in densely packed tissues presents an obstacle to the assignment of high-dimensional spatial features to individual cells for accurate cell-type annotation. We devised a method to correct for lateral spillage of cell surface markers between adjacent cells termed REinforcement Dynamic Spillover EliminAtion (REDSEA). The use of REDSEA decreased contaminating signals from neighboring cells. It improved the recovery of marker signals across both isotopic (i.e., Multiplexed Ion Beam Imaging) and immunofluorescent (i.e., Cyclic Immunofluorescence) multiplexed images resulting in a marked improvement in cell-type classification.
-
(2021) PLoS Computational Biology. 17, 4, e1008887. Abstract
Mass Based Imaging (MBI) technologies such as Multiplexed Ion Beam Imaging by time of flight (MIBI-TOF) and Imaging Mass Cytometry (IMC) allow for the simultaneous measurement of the expression levels of 40 or more proteins in biological tissue, providing insight into cellular phenotypes and organization in situ. Imaging artifacts, resulting from the sample, assay or instrumentation complicate downstream analyses and require correction by domain experts. Here, we present MBI Analysis User Interface (MAUI), a series of graphical user interfaces that facilitate this data pre-processing, including the removal of channel crosstalk, noise and antibody aggregates. Our software streamlines these steps and accelerates processing by enabling real-time and interactive parameter tuning across multiple images.
2020
-
(2020) BioRxiv. Abstract
Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis that is distinctly characterized by granuloma formation within infected tissues. Granulomas are dynamic and organized immune cell aggregates that limit dissemination, but can also hinder bacterial clearance. Consequently, outcome in TB is influenced by how granuloma structure and composition shift the balance between these two functions. To date, our understanding of what factors drive granuloma function in humans is limited. With this in mind, we used Multiplexed Ion Beam Imaging by Time-of-Flight (MIBI-TOF) to profile 37 proteins in tissues from thirteen patients with active TB disease from the U.S. and South Africa. With this dataset, we constructed a comprehensive tissue atlas where the lineage, functional state, and spatial distribution of 19 unique cell subsets were mapped onto eight phenotypically-distinct granuloma microenvironments. This work revealed an immunosuppressed microenvironment specific to TB granulomas with spatially coordinated co-expression of IDO1 and PD-L1 by myeloid cells and proliferating regulatory T cells. Interestingly, this microenvironment lacked markers consistent with T-cell activation, supporting a myeloid-mediated mechanism of immune suppression. We observed similar trends in gene expression of immunoregulatory proteins in a confirmatory transcriptomic analysis of peripheral blood collected from over 1500 individuals with latent or active TB infection and healthy controls across 29 cohorts spanning 14 countries. Notably, PD-L1 gene expression was found to correlate with TB progression and treatment response, supporting its potential use as a blood-based biomarker. Taken together, this study serves as a framework for leveraging independent cohorts and complementary methodologies to understand how local and systemic immune responses are linked in human health and disease.
-
(2020) Journal of Computational Biology. 8, p. 1204-1218 Abstract
Recent in situ multiplexed profiling techniques provide insight into microenvironment formation, maintenance, and transformation through a lens of localized cellular phenotype distribution. In this article, we introduce a method for recovering signatures of microenvironments from such data. We use topic models to identify characteristic cell types overrepresented in neighborhoods that serve as proxies for microenvironment. Furthermore, by assuming spatial coherence among neighboring microenvironments our model limits the number of parameters that need to be learned and permits data-driven decisions about the size of cellular neighborhoods. We apply this method to uncover anatomically known structures in mouse spleen-identifying distinct population of spleen B cells that are defined by their characteristic neighborhoods. Next, we apply the method to a dataset of triple-negative breast cancer tumors from 41 patients to study the structure of tumor-immune boundary. We uncover previously reported tumor-immune microenvironment near the tumor-immune boundary enriched for immune cells with high Indoleamine-pyrrole 2,3-dioxygenase (IDO) and Programmed death ligand 1 (PD-L1) and a novel, immunosuppressed, microenvironment-enriched for cells expressing CD45 and FoxP3.
-
(2020) Nature Reviews Cancer. 1, 2, p. 156-157 Abstract
Tumor heterogeneity remains an obstacle to effective clinical management of breast cancer. Two new studies use high-dimensional imaging of single-cell protein expression in situ in clinical samples to link genomic alterations to multi-cellular features of the tumor microenvironment and reveal breast-cancer phenotypes associated with clinical outcome.
2019
-
(2019) Science Advances. 5, 10, 5851. Abstract
Understanding tissue structure and function requires tools that quantify the expression of multiple proteins while preserving spatial information. Here, we describe MIBI-TOF (multiplexed ion beam imaging by time of flight), an instrument that uses bright ion sources and orthogonal time-of-flight mass spectrometry to image metal-tagged antibodies at subcellular resolution in clinical tissue sections. We demonstrate quantitative, full periodic table coverage across a five-log dynamic range, imaging 36 labeled antibodies simultaneously with histochemical stains and endogenous elements. We image fields of view up to 800 mu m x 800 mu m at resolutions down to 260 nm with sensitivities approaching single-molecule detection. We leverage these properties to interrogate intrapatient heterogeneity in tumor organization in triple-negative breast cancer, revealing regional variability in tumor cell phenotypes in contrast to a structured immune response. Given its versatility and sample back-compatibility, MIBI-TOF is positioned to leverage existing annotated, archival tissue cohorts to explore emerging questions in cancer, immunology, and neurobiology.
-
(2019) Nature Communications. 10, 68. Abstract
Steady-state protein abundance is set by four rates: transcription, translation, mRNA decay and protein decay. A given protein abundance can be obtained from infinitely many combinations of these rates. This raises the question of whether the natural rates for each gene result from historical accidents, or are there rules that give certain combinations a selective advantage? We address this question using high-throughput measurements in rapidly growing cells from diverse organisms to find that about half of the rate combinations do not exist: genes that combine high transcription with low translation are strongly depleted. This depletion is due to a trade-off between precision and economy: high transcription decreases stochastic fluctuations but increases transcription costs. Our theory quantitatively explains which rate combinations are missing, and predicts the curvature of the fitness function for each gene. It may guide the design of gene circuits with desired expression levels and noise.
2018
-
(2018) Cell. 174, 6, p. 1373-1387 Abstract
The immune system is critical in modulating cancer progression, but knowledge of immune composition, phenotype, and interactions with tumor is limited. We used multiplexed ion beam imaging by time-of-flight (MIBI-TOF) to simultaneously quantify in situ expression of 36 proteins covering identity, function, and immune regulation at sub-cellular resolution in 41 triple-negative breast cancer patients. Multi-step processing, including deep-learning-based segmentation, revealed variability in the composition of tumor-immune populations across individuals, reconciled by overall immune infiltration and enriched co-occurrence of immune subpopulations and checkpoint expression. Spatial enrichment analysis showed immune mixed and compartmentalized tumors, coinciding with expression of PD1, PD-L1, and IDO in a cell-type- and location-specific manner. Ordered immune structures along the tumor-immune border were associated with compartmentalization and linked to survival. These data demonstrate organization in the tumor-immune microenvironment that is structured in cellular composition, spatial arrangement, and regulatory-protein expression and provide a framework to apply multiplexed imaging to immune oncology.
2016
-
(2016) Cell. 166, 5, p. 1282-1294.e18 Abstract
Data of gene expression levels across individuals, cell types, and disease states is expanding, yet our understanding of how expression levels impact phenotype is limited. Here, we present a massively parallel system for assaying the effect of gene expression levels on fitness in Saccharomyces cerevisiae by systematically altering the expression level of 100 genes at 100 distinct levels spanning a 500-fold range at high resolution. We show that the relationship between expression levels and growth is gene and environment specific and provides information on the function, stoichiometry, and interactions of genes. Wild-type expression levels in some conditions are not optimal for growth, and genes whose fitness is greatly affected by small changes in expression level tend to exhibit lower cell-to-cell variability in expression. Our study addresses a fundamental gap in understanding the functional significance of gene expression regulation and offers a framework for evaluating the phenotypic effects of expression variation.
-
(2016) PLoS ONE. 11, 4, 0153344. Abstract
Most proteins show changes in level across growth conditions. Many of these changes seem to be coordinated with the specific growth rate rather than the growth environment or the protein function. Although cellular growth rates, gene expression levels and gene regulation have been at the center of biological research for decades, there are only a few models giving a base line prediction of the dependence of the proteome fraction occupied by a gene with the specific growth rate. We present a simple model that predicts a widely coordinated increase in the fraction of many proteins out of the proteome, proportionally with the growth rate. The model reveals how passive redistribution of resources, due to active regulation of only a few proteins, can have proteome wide effects that are quantitatively predictable. Our model provides a potential explanation for why and how such a coordinated response of a large fraction of the proteome to the specific growth rate arises under different environmental conditions. The simplicity of our model can also be useful by serving as a baseline null hypothesis in the search for active regulation. We exemplify the usage of the model by analyzing the relationship between growth rate and proteome composition for the model microorganism E.coli as reflected in recent proteomics data sets spanning various growth conditions. We find that the fraction out of the proteome of a large number of proteins, and from different cellular processes, increases proportionally with the growth rate. Notably, ribosomal proteins, which have been previously reported to increase in fraction with growth rate, are only a small part of this group of proteins. We suggest that, although the fractions of many proteins change with the growth rate, such changes may be partially driven by a global effect, not necessarily requiring specific cellular control mechanisms.
2015
-
(2015) Genome Research. 25, p. 1893-1902 Abstract
Genetically identical cells exposed to the same environment display variability in gene expression (noise), with important consequences for the fidelity of cellular regulation and biological function. Although population average gene expression is tightly coupled to growth rate, the effects of changes in environmental conditions on expression variability are not known. Here, we measure the single-cell expression distributions of approximately 900 Saccharomyces cerevisiæ promoters across four environmental conditions using flow cytometry, and find that gene expression noise is tightly coupled to the environment and is generally higher at lower growth rates. Nutrient-poor conditions, which support lower growth rates, display elevated levels of noise for most promoters, regardless of their specific expression values. We present a simple model of noise in expression that results from having an asynchronous population, with cells at different cell-cycle stages, and with different partitioning of the cells between the stages at different growth rates. This model predicts non-monotonic global changes in noise at different growth rates as well as overall higher variability in expression for cell-cycle-regulated genes in all conditions. The consistency between this model and our data, as well as with noise measurements of cells growing in a chemostat at well-defined growth rates, suggests that cell-cycle heterogeneity is a major contributor to gene expression noise. Finally, we identify gene and promoter features that play a role in gene expression noise across conditions. Our results show the existence of growth-related global changes in gene expression noise and suggest their potential phenotypic implications.
2014
-
(2014) Genome Research. 24, 10, p. 1698-1706 Abstract
Genetically identical cells exhibit large variability (noise) in gene expression, with important consequences for cellular function. Although the amount of noise decreases with and is thus partly determined by the mean expression level, the extent to which different promoter sequences can deviate away from this trend is not fully known. Here, we present a high-throughput method for measuring promoter-driven noise for thousands of designed synthetic promoters in parallel. We use it to investigate how promoters encode different noise levels and find that the noise levels of promoters with similar mean expression levels can vary more than one order of magnitude, with nucleosome-disfavoring sequences resulting in lower noise and more transcription factor binding sites resulting in higher noise. We propose a kinetic model of gene expression that takes into account the nonspecific DNA binding and one-dimensional sliding along the DNA, which occurs when transcription factors search for their target sites. We show that this assumption can improve the prediction of the mean-independent component of expression noise for our designed promoter sequences, suggesting that a transcription factor target search may affect gene expression noise. Consistent with our findings in designed promoters, we find that binding-site multiplicity in native promoters is associated with higher expression noise. Overall, our results demonstrate that small changes in promoter DNA sequence can tune noise levels in a manner that is predictable and partly decoupled from effects on the mean expression levels. These insights may assist in designing promoters with desired noise levels.
2013
-
(2013) Molecular Systems Biology. 9, 701. Abstract
Most genes change expression levels across conditions, but it is unclear which of these changes represents specific regulation and what determines their quantitative degree. Here, we accurately measured activities of ∼900 S. cerevisiae and ∼1800 E. coli promoters using fluorescent reporters. We show that in both organisms 60-90% of promoters change their expression between conditions by a constant global scaling factor that depends only on the conditions and not on the promoter's identity. Quantifying such global effects allows precise characterization of specific regulation - promoters deviating from the global scale line. These are organized into few functionally related groups that also adhere to scale lines and preserve their relative activities across conditions. Thus, only several scaling factors suffice to accurately describe genome-wide expression profiles across conditions. We present a parameter-free passive resource allocation model that quantitatively accounts for the global scaling factors. It suggests that many changes in expression across conditions result from global effects and not specific regulation, and provides means for quantitative interpretation of expression profiles.
-
Sequence features of yeast and human core promoters that are predictive of maximal promoter activity(2013) Nucleic Acids Research. 41, 11, p. 5569-5581 Abstract
The core promoter is the region in which RNA polymerase II is recruited to the DNA and acts to initiate transcription, but the extent to which the core promoter sequence determines promoter activity levels is largely unknown. Here, we identified several base content and k-mer sequence features of the yeast core promoter sequence that are highly predictive of maximal promoter activity. These features are mainly located in the region 75 bp upstream and 50 bp downstream of the main transcription start site, and their associations hold for both constitutively active promoters and promoters that are induced or repressed in specific conditions. Our results unravel several architectural features of yeast core promoters and suggest that the yeast core promoter sequence downstream of the TATA box (or of similar sequences involved in recruitment of the pre-initiation complex) is a major determinant of maximal promoter activity. We further show that human core promoters also contain features that are indicative of maximal promoter activity; thus, our results emphasize the important role of the core promoter sequence in transcriptional regulation.
-
(2013) PLoS Computational Biology. 9, 3, e1002934. Abstract
A full understanding of gene regulation requires an understanding of the contributions that the various regulatory regions have on gene expression. Although it is well established that sequences downstream of the main promoter can affect expression, our understanding of the scale of this effect and how it is encoded in the DNA is limited. Here, to measure the effect of native S. cerevisiae 3′ end sequences on expression, we constructed a library of 85 fluorescent reporter strains that differ only in their 3′ end region. Notably, despite being driven by the same strong promoter, our library spans a continuous twelve-fold range of expression values. These measurements correlate with endogenous mRNA levels, suggesting that the 3′ end contributes to constitutive differences in mRNA levels. We used deep sequencing to map the 3′UTR ends of our strains and show that determination of polyadenylation sites is intrinsic to the local 3′ end sequence. Polyadenylation mapping was followed by sequence analysis, we found that increased A/T content upstream of the main polyadenylation site correlates with higher expression, both in the library and genome-wide, suggesting that native genes differ by the encoded efficiency of 3′ end processing. Finally, we use single cells fluorescence measurements, in different promoter activation levels, to show that 3′ end sequences modulate protein expression dynamics differently than promoters, by predominantly affecting the size of protein production bursts as opposed to the frequency at which these bursts occur. Altogether, our results lead to a more complete understanding of gene regulation by demonstrating that 3′ end regions have a unique and sequence dependent effect on gene expression.
-
(2013) GENOME BIOLOGY. 14, 11, Abstract
A new study exploits the time-dependence of formaldehyde cross-linking in the commonly used chromatin immunoprecipitation (ChIP) assay to infer the on and off rates for site-specific chromatin interactions.
2012
-
(2012) Nature Genetics. 44, 7, p. 743-750 Abstract
Understanding how precise control of gene expression is specified within regulatory DNA sequences is a key challenge with far-reaching implications. Many studies have focused on the regulatory role of transcription factor-binding sites. Here, we explore the transcriptional effects of different elements, nucleosome-disfavoring sequences and, specifically, poly(dA:dT) tracts that are highly prevalent in eukaryotic promoters. By measuring promoter activity for a large-scale promoter library, designed with systematic manipulations to the properties and spatial arrangement of poly(dA:dT) tracts, we show that these tracts significantly and causally affect transcription. We show that manipulating these elements offers a general genetic mechanism, applicable to promoters regulated by different transcription factors, for tuning expression in a predictable manner, with resolution that can be even finer than that attained by altering transcription factor sites. Overall, our results advance the understanding of the regulatory code and suggest a potential mechanism by which promoters yielding prespecified expression patterns can be designed.
-
(2012) Nature Biotechnology. 30, 6, p. 521-+ Abstract
Despite extensive research, our understanding of the rules according to which cis-regulatory sequences are converted into gene expression is limited. We devised a method for obtaining parallel, highly accurate gene expression measurements from thousands of designed promoters and applied it to measure the effect of systematic changes in the location, number, orientation, affinity and organization of transcription-factor binding sites and nucleosome-disfavoring sequences. Our analyses reveal a clear relationship between expression and binding-site multiplicity, as well as dependencies of expression on the distance between transcription-factor binding sites and gene starts which are transcription-factor specific, including a striking similar to 10-bp periodic relationship between gene expression and binding-site location. We show how this approach can measure transcription-factor sequence specificities and the sensitivity of transcription-factor sites to the surrounding sequence context, and compare the activity of 75 yeast transcription factors. Our method can be used to study both cis and trans effects of genotype on transcriptional, post-transcriptional and translational control.
2011
-
(2011) Genome Research. 21, 12, p. 2114-2128 Abstract
Coordinate regulation of ribosomal protein (RP) genes is key for controlling cell growth. In yeast, it is unclear how this regulation achieves the required equimolar amounts of the different RP components, given that some RP genes exist in duplicate copies, while others have only one copy. Here, we tested whether the solution to this challenge is partly encoded within the DNA sequence of the RP promoters, by fusing 110 different RP promoters to a fluorescent gene reporter, allowing us to robustly detect differences in their promoter activities that are as small as ∼10%. We found that single-copy RP promoters have significantly higher activities, suggesting that proper RP stoichiometry is indeed partly encoded within the RP promoters. Notably, we also partially uncovered how this regulation is encoded by finding that RP promoters with higher activity have more nucleosome-disfavoring sequences and characteristic spatial organizations of these sequences and of binding sites for key RP regulators. Mutations in these elements result in a significant decrease of RP promoter activity. Thus, our results suggest that intrinsic (DNA-dependent) nucleosome organization may be a key mechanism by which genomes encode biologically meaningful promoter activities. Our approach can readily be applied to uncover how transcriptional programs of other promoters are encoded.
2010
-
(2010) Trends in Genetics. 26, 8, p. 335-340 Abstract
The recently discovered prokaryotic immune system known as CRISPR (clustered regularly interspaced short palindromic repeats) is based on small RNAs ('spacers') that restrict phage and plasmid infection. It has been hypothesized that CRISPRs can also regulate self gene expression by utilizing spacers that target self genes. By analyzing CRISPRs from 330 organisms we found that one in every 250 spacers is self-targeting, and that such self-targeting occurs in 18% of all CRISPR-bearing organisms. However, complete lack of conservation across species, combined with abundance of degraded repeats near self-targeting spacers, suggests that self-targeting is a form of autoimmunity rather than a regulatory mechanism. We propose that accidental incorporation of self nucleic acids by CRISPR can incur an autoimmune fitness cost, and this could explain the abundance of degraded CRISPR systems across prokaryotes.