Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection.

Identification of host genes essential for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may reveal novel therapeutic targets and inform our understanding of coronavirus disease 2019 (COVID-19) pathogenesis. Here we performed genome-wide CRISPR screens in Vero-E6 cells with SARS-CoV-2, Middle East respiratory syndrome CoV (MERS-CoV), bat CoV HKU5 expressing the SARS-CoV-1 spike, and vesicular stomatitis virus (VSV) expressing the SARS-CoV-2 spike. We identified known SARS-CoV-2 host factors, including the receptor ACE2 and protease Cathepsin L. We additionally discovered pro-viral genes and pathways, including HMGB1 and the SWI/SNF chromatin remodeling complex, that are SARS lineage and pan-coronavirus specific, respectively. We show that HMGB1 regulates ACE2 expression and is critical for entry of SARS-CoV-2, SARS-CoV-1, and NL63. We also show that small-molecule antagonists of identified gene products inhibited SARS-CoV-2 infection in monkey and human cells, demonstrating the conserved role of these genetic hits across species. This identifies potential therapeutic targets for SARS-CoV-2 and reveals SARS lineage-specific and pan-CoV host factors that regulate susceptibility to highly pathogenic CoVs.

Recovering Gene Interactions from Single-Cell Data Using Data Diffusion.

Single-cell RNA sequencing technologies suffer from many sources of technical noise, including under-sampling of mRNA molecules, often termed "dropout," which can severely obscure important gene-gene relationships. To address this, we developed MAGIC (Markov affinity-based graph imputation of cells), a method that shares information across similar cells, via data diffusion, to denoise the cell count matrix and fill in missing transcripts. We validate MAGIC on several biological systems and find it effective at recovering gene-gene relationships and additional structures. Applied to the epithilial to mesenchymal transition, MAGIC reveals a phenotypic continuum, with the majority of cells residing in intermediate states that display stem-like signatures, and infers known and previously uncharacterized regulatory interactions, demonstrating that our approach can successfully uncover regulatory relations without perturbations.

Single-cell multiomics reveals persistence of HIV-1 in expanded cytotoxic T cell clones.

Understanding the drivers and markers of clonally expanding HIV-1-infected CD4+ T cells is essential for HIV-1 eradication. We used single-cell ECCITE-seq, which captures surface protein expression, cellular transcriptome, HIV-1 RNA, and TCR sequences within the same single cell to track clonal expansion dynamics in longitudinally archived samples from six HIV-1-infected individuals (during viremia and after suppressive antiretroviral therapy) and two uninfected individuals, in unstimulated conditions and after CMV and HIV-1 antigen stimulation. Despite antiretroviral therapy, persistent antigen and TNF responses shaped T cell clonal expansion. HIV-1 resided in Th1-polarized, antigen-responding T cells expressing BCL2 and SERPINB9 that may resist cell death. HIV-1 RNA+ T cell clones were larger in clone size, established during viremia, persistent after viral suppression, and enriched in GZMB+ cytotoxic effector memory Th1 cells. Targeting HIV-1-infected cytotoxic CD4+ T cells and drivers of clonal expansion provides another direction for HIV-1 eradication.

Immune dysregulation and autoreactivity correlate with disease severity in SARS-CoV-2-associated multisystem inflammatory syndrome in children.

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV-2 infection. We profiled MIS-C, adult COVID-19, and healthy pediatric and adult individuals using single-cell RNA sequencing, flow cytometry, antigen receptor repertoire analysis, and unbiased serum proteomics, which collectively identified a signature in MIS-C patients that correlated with disease severity. Despite having no evidence of active infection, MIS-C patients had elevated S100A-family alarmins and decreased antigen presentation signatures, indicative of myeloid dysfunction. MIS-C patients showed elevated expression of cytotoxicity genes in NK and CD8+ T cells and expansion of specific IgG-expressing plasmablasts. Clinically severe MIS-C patients displayed skewed memory T cell TCR repertoires and autoimmunity characterized by endothelium-reactive IgG. The alarmin, cytotoxicity, TCR repertoire, and plasmablast signatures we defined have potential for application in the clinic to better diagnose and potentially predict disease severity early in the course of MIS-C.

Distinguishing features of long COVID identified through immune profiling.

Post-acute infection syndromes may develop after acute viral disease1. Infection with SARS-CoV-2 can result in the development of a post-acute infection syndrome known as long COVID. Individuals with long COVID frequently report unremitting fatigue, post-exertional malaise, and a variety of cognitive and autonomic dysfunctions2-4. However, the biological processes that are associated with the development and persistence of these symptoms are unclear. Here 275 individuals with or without long COVID were enrolled in a cross-sectional study that included multidimensional immune phenotyping and unbiased machine learning methods to identify biological features associated with long COVID. Marked differences were noted in circulating myeloid and lymphocyte populations relative to the matched controls, as well as evidence of exaggerated humoral responses directed against SARS-CoV-2 among participants with long COVID. Furthermore, higher antibody responses directed against non-SARS-CoV-2 viral pathogens were observed among individuals with long COVID, particularly Epstein-Barr virus. Levels of soluble immune mediators and hormones varied among groups, with cortisol levels being lower among participants with long COVID. Integration of immune phenotyping data into unbiased machine learning models identified the key features that are most strongly associated with long COVID status. Collectively, these findings may help to guide future studies into the pathobiology of long COVID and help with developing relevant biomarkers.

Causal identification of single-cell experimental perturbation effects with CINEMA-OT.

Recent advancements in single-cell technologies allow characterization of experimental perturbations at single-cell resolution. While methods have been developed to analyze such experiments, the application of a strict causal framework has not yet been explored for the inference of treatment effects at the single-cell level. Here we present a causal-inference-based approach to single-cell perturbation analysis, termed CINEMA-OT (causal independent effect module attribution + optimal transport). CINEMA-OT separates confounding sources of variation from perturbation effects to obtain an optimal transport matching that reflects counterfactual cell pairs. These cell pairs represent causal perturbation responses permitting a number of novel analyses, such as individual treatment-effect analysis, response clustering, attribution analysis, and synergy analysis. We benchmark CINEMA-OT on an array of treatment-effect estimation tasks for several simulated and real datasets and show that it outperforms other single-cell perturbation analysis methods. Finally, we perform CINEMA-OT analysis of two newly generated datasets: (1) rhinovirus and cigarette-smoke-exposed airway organoids, and (2) combinatorial cytokine stimulation of immune cells. In these experiments, CINEMA-OT reveals potential mechanisms by which cigarette-smoke exposure dulls the airway antiviral response, as well as the logic that governs chemokine secretion and peripheral immune cell recruitment.

Anatomical Diversity of the Adult Corticospinal Tract Revealed by Single-Cell Transcriptional Profiling.

The corticospinal tract (CST) forms a central part of the voluntary motor apparatus in all mammals. Thus, injury, disease, and subsequent degeneration within this pathway result in chronic irreversible functional deficits. Current strategies to repair the damaged CST are suboptimal in part because of underexplored molecular heterogeneity within the adult tract. Here, we combine spinal retrograde CST tracing with single-cell RNA sequencing (scRNAseq) in adult male and female mice to index corticospinal neuron (CSN) subtypes that differentially innervate the forelimb and hindlimb. We exploit publicly available datasets to confer anatomic specialization among CSNs and show that CSNs segregate not only along the forelimb and hindlimb axis but also by supraspinal axon collateralization. These anatomically defined transcriptional data allow us to use machine learning tools to build classifiers that discriminate between CSNs and cortical layer 2/3 and nonspinally terminating layer 5 neurons in M1 and separately identify limb-specific CSNs. Using these tools, CSN subtypes can be differentially identified to study postnatal patterning of the CST in vivo, leveraged to screen for novel limb-specific axon growth survival and growth activators in vitro, and ultimately exploited to repair the damaged CST after injury and disease.SIGNIFICANCE STATEMENT Therapeutic interventions designed to repair the damaged CST after spinal cord injury have remained functionally suboptimal in part because of an incomplete understanding of the molecular heterogeneity among subclasses of CSNs. Here, we combine spinal retrograde labeling with scRNAseq and annotate a CSN index by the termination pattern of their primary axon in the cervical or lumbar spinal cord and supraspinal collateral terminal fields. Using machine learning we have confirmed the veracity of our CSN gene lists to train classifiers to identify CSNs among all classes of neurons in primary motor cortex to study the development, patterning, homeostasis, and response to injury and disease, and ultimately target streamlined repair strategies to this critical motor pathway.

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium identifies target cells, alterations in gene expression, and cell state changes.

There are currently limited Food and Drug Administration (FDA)-approved drugs and vaccines for the treatment or prevention of Coronavirus Disease 2019 (COVID-19). Enhanced understanding of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and pathogenesis is critical for the development of therapeutics. To provide insight into viral replication, cell tropism, and host-viral interactions of SARS-CoV-2, we performed single-cell (sc) RNA sequencing (RNA-seq) of experimentally infected human bronchial epithelial cells (HBECs) in air-liquid interface (ALI) cultures over a time course. This revealed novel polyadenylated viral transcripts and highlighted ciliated cells as a major target at the onset of infection, which we confirmed by electron and immunofluorescence microscopy. Over the course of infection, the cell tropism of SARS-CoV-2 expands to other epithelial cell types including basal and club cells. Infection induces cell-intrinsic expression of type I and type III interferons (IFNs) and interleukin (IL)-6 but not IL-1. This results in expression of interferon-stimulated genes (ISGs) in both infected and bystander cells. This provides a detailed characterization of genes, cell types, and cell state changes associated with SARS-CoV-2 infection in the human airway.

Uncovering axes of variation among single-cell cancer specimens.

While several tools have been developed to map axes of variation among individual cells, no analogous approaches exist for identifying axes of variation among multicellular biospecimens profiled at single-cell resolution. For this purpose, we developed 'phenotypic earth mover's distance' (PhEMD). PhEMD is a general method for embedding a 'manifold of manifolds', in which each datapoint in the higher-level manifold (of biospecimens) represents a collection of points that span a lower-level manifold (of cells). We apply PhEMD to a newly generated drug-screen dataset and demonstrate that PhEMD uncovers axes of cell subpopulational variation among a large set of perturbation conditions. Moreover, we show that PhEMD can be used to infer the phenotypes of biospecimens not directly profiled. Applied to clinical datasets, PhEMD generates a map of the patient-state space that highlights sources of patient-to-patient variation. PhEMD is scalable, compatible with leading batch-effect correction techniques and generalizable to multiple experimental designs.

Exploring single-cell data with deep multitasking neural networks.

It is currently challenging to analyze single-cell data consisting of many cells and samples, and to address variations arising from batch effects and different sample preparations. For this purpose, we present SAUCIE, a deep neural network that combines parallelization and scalability offered by neural networks, with the deep representation of data that can be learned by them to perform many single-cell data analysis tasks. Our regularizations (penalties) render features learned in hidden layers of the neural network interpretable. On large, multi-patient datasets, SAUCIE's various hidden layers contain denoised and batch-corrected data, a low-dimensional visualization and unsupervised clustering, as well as other information that can be used to explore the data. We analyze a 180-sample dataset consisting of 11 million T cells from dengue patients in India, measured with mass cytometry. SAUCIE can batch correct and identify cluster-based signatures of acute dengue infection and create a patient manifold, stratifying immune response to dengue.

The association between preoperative body composition and aerobic fitness in patients scheduled for colorectal surgery.

AIM: Although cardiopulmonary exercise testing (CPET) is considered the gold standard, a preoperative abdominal CT scan might also provide information concerning preoperative aerobic fitness for risk assessment. This study aimed to investigate the association between preoperative CT-scan-derived body composition variables and preoperative CPET variables of aerobic fitness in colorectal surgery. METHOD: In this retrospective cohort study, CT images at level L3 were analysed for skeletal muscle mass, skeletal muscle radiation attenuation, visceral adipose tissue (VAT) mass and subcutaneous adipose tissue mass. Regression analyses were performed to investigate the relation between CT-scan-derived body composition variables, CPET-derived aerobic fitness and other preoperative patient-related variables. Logistic regression analysis was performed to predict a preoperative anaerobic threshold (AT) ≤ 11.1 ml/kg/min as cut-off for having a high risk for postoperative complications. RESULTS: Data from 78 patients (45 men; mean [SD] age 74.5 [6.4 years]) were analysed. A correlation coefficient of 0.55 was observed between absolute AT and skeletal muscle mass index. Absolute AT (R2 of 51.1%) was lower in patients with a lower skeletal muscle mass index, together with higher age, lower body mass and higher American Society of Anesthesiologists (ASA) score. Higher ASA score (odds ratio 5.64; P = 0.033) and higher VAT mass (odds ratio 1.02; P = 0.036) were associated with an increased risk of an AT ≤ 11.1 ml/kg/min. CONCLUSION: Body composition variables from the preoperative CT scan were moderately associated with preoperative CPET-derived aerobic fitness. Higher ASA score and higher VAT mass were associated with an increased risk of an AT ≤ 11.1 ml/kg/min.

A phenomapping-derived tool to personalize the selection of anatomical vs. functional testing in evaluating chest pain (ASSIST).

AIMS: Coronary artery disease is frequently diagnosed following evaluation of stable chest pain with anatomical or functional testing. A more granular understanding of patient phenotypes that benefit from either strategy may enable personalized testing. METHODS AND RESULTS: Using participant-level data from 9572 patients undergoing anatomical (n = 4734) vs. functional (n = 4838) testing in the PROMISE (PROspective Multicenter Imaging Study for Evaluation of Chest Pain) trial, we created a topological representation of the study population based on 57 pre-randomization variables. Within each patient's 5% topological neighbourhood, Cox regression models provided individual patient-centred hazard ratios for major adverse cardiovascular events and revealed marked heterogeneity across the phenomap [median 1.11 (10th to 90th percentile: 0.52-2.61]), suggestive of distinct phenotypic neighbourhoods favouring anatomical or functional testing. Based on this risk phenomap, we employed an extreme gradient boosting algorithm in 80% of the PROMISE population to predict the personalized benefit of anatomical vs. functional testing using 12 model-derived, routinely collected variables and created a decision support tool named ASSIST (Anatomical vs. Stress teSting decIsion Support Tool). In both the remaining 20% of PROMISE and an external validation set consisting of patients from SCOT-HEART (Scottish COmputed Tomography of the HEART Trial) undergoing anatomical-first vs. functional-first assessment, the testing strategy recommended by ASSIST was associated with a significantly lower incidence of each study's primary endpoint (P = 0.0024 and P = 0.0321 for interaction, respectively), as well as a harmonized endpoint of all-cause mortality or non-fatal myocardial infarction (P = 0.0309 and P < 0.0001 for interaction, respectively). CONCLUSION: We propose a novel phenomapping-derived decision support tool to standardize the selection of anatomical vs. functional testing in the evaluation of stable chest pain, validated in two large and geographically diverse clinical trial populations.

Large-scale mapping of gene regulatory logic reveals context-dependent repression by transcriptional activators.

Transcription factors (TFs) are key mediators that propagate extracellular and intracellular signals through to changes in gene expression profiles. However, the rules by which promoters decode the amount of active TF into target gene expression are not well understood. To determine the mapping between promoter DNA sequence, TF concentration, and gene expression output, we have conducted in budding yeast a large-scale measurement of the activity of thousands of designed promoters at six different levels of TF. We observe that maximum promoter activity is determined by TF concentration and not by the number of binding sites. Surprisingly, the addition of an activator site often reduces expression. A thermodynamic model that incorporates competition between neighboring binding sites for a local pool of TF molecules explains this behavior and accurately predicts both absolute expression and the amount by which addition of a site increases or reduces expression. Taken together, our findings support a model in which neighboring binding sites interact competitively when TF is limiting but otherwise act additively.

Novel Toxin-Antitoxin Module SlvT-SlvA Regulates Megaplasmid Stability and Incites Solvent Tolerance in Pseudomonas putida S12.

Pseudomonas putida S12 is highly tolerant of organic solvents in saturating concentrations, rendering this microorganism suitable for the industrial production of various aromatic compounds. Previous studies revealed that P. putida S12 contains the single-copy 583-kbp megaplasmid pTTS12. pTTS12 carries several important operons and gene clusters facilitating P. putida S12 survival and growth in the presence of toxic compounds or other environmental stresses. We wished to revisit and further scrutinize the role of pTTS12 in conferring solvent tolerance. To this end, we cured the megaplasmid from P. putida S12 and conclusively confirmed that the SrpABC efflux pump is the major determinant of solvent tolerance on the megaplasmid pTTS12. In addition, we identified a novel toxin-antitoxin module (proposed gene names slvT and slvA, respectively) encoded on pTTS12 which contributes to the solvent tolerance phenotype and is important for conferring stability to the megaplasmid. Chromosomal introduction of the srp operon in combination with the slvAT gene pair created a solvent tolerance phenotype in non-solvent-tolerant strains, such as P. putida KT2440, Escherichia coli TG1, and E. coli BL21(DE3).IMPORTANCE Sustainable alternatives for high-value chemicals can be achieved by using renewable feedstocks in bacterial biocatalysis. However, during the bioproduction of such chemicals and biopolymers, aromatic compounds that function as products, substrates, or intermediates in the production process may exert toxicity to microbial host cells and limit the production yield. Therefore, solvent tolerance is a highly preferable trait for microbial hosts in the biobased production of aromatic chemicals and biopolymers. In this study, we revisit the essential role of megaplasmid pTTS12 from solvent-tolerant Pseudomonas putida S12 for molecular adaptation to an organic solvent. In addition to the solvent extrusion pump (SrpABC), we identified a novel toxin-antitoxin module (SlvAT) which contributes to short-term tolerance in moderate solvent concentrations, as well as to the stability of pTTS12. These two gene clusters were successfully expressed in non-solvent-tolerant strains of P. putida and Escherichia coli strains to confer and enhance solvent tolerance.

Development and Validation of the Quick COVID-19 Severity Index: A Prognostic Tool for Early Clinical Decompensation.

STUDY OBJECTIVE: The goal of this study is to create a predictive, interpretable model of early hospital respiratory failure among emergency department (ED) patients admitted with coronavirus disease 2019 (COVID-19). METHODS: This was an observational, retrospective, cohort study from a 9-ED health system of admitted adult patients with severe acute respiratory syndrome coronavirus 2 (COVID-19) and an oxygen requirement less than or equal to 6 L/min. We sought to predict respiratory failure within 24 hours of admission as defined by oxygen requirement of greater than 10 L/min by low-flow device, high-flow device, noninvasive or invasive ventilation, or death. Predictive models were compared with the Elixhauser Comorbidity Index, quick Sequential [Sepsis-related] Organ Failure Assessment, and the CURB-65 pneumonia severity score. RESULTS: During the study period, from March 1 to April 27, 2020, 1,792 patients were admitted with COVID-19, 620 (35%) of whom had respiratory failure in the ED. Of the remaining 1,172 admitted patients, 144 (12.3%) met the composite endpoint within the first 24 hours of hospitalization. On the independent test cohort, both a novel bedside scoring system, the quick COVID-19 Severity Index (area under receiver operating characteristic curve mean 0.81 [95% confidence interval {CI} 0.73 to 0.89]), and a machine-learning model, the COVID-19 Severity Index (mean 0.76 [95% CI 0.65 to 0.86]), outperformed the Elixhauser mortality index (mean 0.61 [95% CI 0.51 to 0.70]), CURB-65 (0.50 [95% CI 0.40 to 0.60]), and quick Sequential [Sepsis-related] Organ Failure Assessment (0.59 [95% CI 0.50 to 0.68]). A low quick COVID-19 Severity Index score was associated with a less than 5% risk of respiratory decompensation in the validation cohort. CONCLUSION: A significant proportion of admitted COVID-19 patients progress to respiratory failure within 24 hours of admission. These events are accurately predicted with bedside respiratory examination findings within a simple scoring system.

The reproducibility of skeletal muscle signal intensity on routine magnetic resonance imaging in Crohn's disease.

BACKGROUND AND AIM: Myosteatosis is a prognostic factor in cancer and liver cirrhosis. It can be determined noninvasively using computed tomography or, as shown recently, by magnetic resonance (MR) imaging. The primary aim was to analyze the reproducibility of skeletal muscle signal intensity on routine MR-enterographies, as indicator of myosteatosis, in Crohn's disease (CD) and to explore the association between skeletal muscle signal intensity at diagnosis with time to intestinal resection. METHODS: CD patients undergoing MR-enterography within 6 months from diagnosis and having a maximum of 5 years follow-up were included. Skeletal muscle signal intensity was analyzed on T1-weighted fat-saturated post-contrast images. Intra-observer and inter-observer reproducibilities were assessed by intra-class correlation coefficient and Cohen's kappa. Intra-observer and inter-observer variabilities were determined by Pearson correlation coefficient and displayed by Bland-Altman plots. Time to intestinal resection was studied by Kaplan-Meier analysis. RESULTS: Median time between diagnosis and MR-enterography was 5 weeks (inter-quartile range 1-9) in 35 CD patients. Skeletal muscle signal intensity showed good intra-class correlation and substantial agreement (for intra-observer, intraclass correlation coefficient = 0.948, κ = 0.677; and inter-observer reproducibility, intraclass correlation coefficient = 0.858, κ = 0.622). Resection free survival was shorter in the low skeletal muscle signal intensity group (P = 0.037). CONCLUSION: Skeletal muscle signal intensity on routine MR-enterographies is reproducible and was associated with unfavorable disease outcome, indicating potential clinical relevance.

Comment on "Pressure enhancement in carbon nanopores: a major confinement effect" by Y. Long, J. C. Palmer, B. Coasne, M. Sliwinska-Bartkowiak and K. E. Gubbins, Phys. Chem. Chem. Phys., 2011, 13, 17163.

In an article from Phys. Chem. Chem. Phys., 2011, 13, 17163, it is claimed that the microscopic local pressure is very high in a certain family of physical systems and that this phenomenon explains previously reported phase equilibrium and chemical reaction equilibrium data. The results provided in the article are based on two arbitrary choices. Thus, the results are arbitrary, and the conclusions appear to be unjustified.

Low thoracic muscle radiation attenuation is associated with postoperative pneumonia following partial hepatectomy for colorectal metastasis.

BACKGROUND: Low skeletal muscle radiation attenuation (SM-RA) is indicative of myosteatosis and diminished muscle function. It is predictive of poor outcome following oncological surgery in several cancer types. Postoperative pneumonia is a known risk factor for increased postoperative mortality. We hypothesized that low SM-RA of the respiratory muscles at the 4th thoracic-vertebra (T4) is associated with postoperative pneumonia following liver surgery. METHODS: Postoperative pneumonia was identified using prospective infection control data. Computed tomography body composition analysis was performed at the L3-and T4 level to determine SM-RA. Body composition variables were corrected for confounders and related to postoperative pneumonia and admission time by multivariable logistic regression. RESULTS: Body composition analysis of 180 patients was performed. Twenty-one patients developed postoperative pneumonia (11.6%). Multivariable analysis showed that low T4 SM-RA as well as low L3 SM-RA were significantly associated with postoperative pneumonia (OR 3.65, 95% CI 1.41-9.49, p < 0.01) and (OR 3.22, 95% CI 1.20-8.61, p = 0.02, respectively). CONCLUSION: Low SM-RA at either the L3-or T4-level is associated with a higher risk of postoperative pneumonia following CLRM resection.

A novel human cell culture model to study visceral smooth muscle phenotypic modulation in health and disease.

Adaptation of the smooth muscle cell (SMC) phenotype is essential for homeostasis and is often involved in pathologies of visceral organs (e.g., uterus, bladder, gastrointestinal tract). In vitro studies of the behavior of visceral SMCs under (patho)-physiological conditions are hampered by a spontaneous, uncontrolled phenotypic modulation of visceral SMCs under regular tissue culture conditions. We aimed to develop a new visceral SMC culture model that allows controlled phenotypic modulation. Human uterine SMCs [ULTR and telomerase-immortalized human myometrial cells (hTERT-HM)] were grown to confluency and kept for up to 6 days on regular tissue culture surfaces or basement membrane (BM) matrix-coated surfaces in the presence of 0-10% serum. mRNA and protein expression and localization of SMC-specific phenotype markers and their transcriptional regulators were investigated by quantitative PCR, Western blotting, and immunofluorescence. Maintaining visceral SMCs confluent for 6 days increased α-smooth muscle actin (1.9-fold) and smooth muscle protein 22-α (3.1-fold), whereas smooth muscle myosin heavy chain was only slightly upregulated (1.3-fold). Culturing on a BM matrix-coated surface further increased these proteins and also markedly promoted mRNA expression of γ-smooth muscle actin (15.0-fold), smoothelin (3.5-fold), h-caldesmon (5.2-fold), serum response factor (7.6-fold), and myocardin (8.1-fold). Whereas additional serum deprivation only minimally affected contractile markers, platelet-derived growth factor-BB and transforming growth factor ß1 consistently reduced versus increased their expression. In conclusion, we present a simple and reproducible visceral SMC culture system that allows controlled phenotypic modulation toward both the synthetic and the contractile phenotype. This may greatly facilitate the identification of factors that drive visceral SMC phenotypic changes in health and disease.

Noise in gene expression is coupled to growth rate.

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 cerevisiae 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.