What are housekeeping genes?

The concept of “housekeeping gene” has been used for four decades but remains loosely defined. Housekeeping genes are commonly described as “essential for cellular existence regardless of their specific function in the tissue or organism”, and “stably expressed irrespective of tissue type, developmental stage, cell cycle state, or external signal”. However, experimental support for the tenet that gene essentiality is linked to stable expression across cell types, conditions, and organisms has been limited. Here we use genome-scale functional genomic screens together with bulk and single-cell sequencing technologies to test this link and optimize a quantitative and experimentally validated definition of housekeeping gene. Using the optimized definition, we identify, characterize, and provide as resources, housekeeping gene lists extracted from several human datasets, and 10 other animal species that include primates, chicken, and C. elegans. We find that stably expressed genes are not necessarily essential, and that the individual genes that are essential and stably expressed can considerably differ across organisms; yet the pathways enriched among these genes are conserved. Further, the level of conservation of housekeeping genes across the analyzed organisms captures their taxonomic groups, showing evolutionary relevance for our definition. Therefore, we present a quantitative and experimentally supported definition of housekeeping genes that can contribute to better understanding of their unique biological and evolutionary characteristics.     

WDR5 Expression Is Prognostic of Breast Cancer Outcome

WDR5 is a core component of the human mixed lineage leukemia-2 complex, which plays central roles in ER positive tumour cells and is a major driver of androgen-dependent prostate cancer cell proliferation. Given the similarities between breast and prostate cancers, we explore the potential prognostic value of WDR5 gene expression on breast cancer survival. Our findings reveal that WDR5 over-expression is associated with poor breast cancer clinical outcome in three gene expression data sets and BreastMark. The eQTL analysis reveals 130 trans-eQTL SNPs whose genes mapped with statistical significance are significantly associated with patient survival. These genes together with WDR5 are enriched with “cellular development, gene expression, cell cycle” signallings. Knocking down WDR5 in MCF7 dramatically decreases cell viability, but does not alter tumour cell response to doxorubicin. Our study reveals the prognostic value of WDR5 expression in breast cancer which is under long-range regulation of genes involved in cell cycle, and anthracycline could be coupled with treatments targeting WDR5 once such a regimen is available.     

Coordination between cell proliferation and apoptosis after DNA damage in Drosophila

Exposure to genotoxic stress promotes cell cycle arrest and DNA repair or apoptosis. These “life” or “death” cell fate decisions often rely on the activity of the tumor suppressor gene p53. Therefore, the precise regulation of p53 is essential to maintain tissue homeostasis and to prevent cancer development. However, how cell cycle progression has an impact on p53 cell fate decision-making is mostly unknown. In this work, we demonstrate that Drosophila p53 proapoptotic activity can be impacted by the G2/M kinase Cdk1. We find that cell cycle arrested or endocycle-induced cells are refractory to ionizing radiation-induced apoptosis. We show that p53 binding to the regulatory elements of the proapoptotic genes and its ability to activate their expression is compromised in experimentally arrested cells. Our results indicate that p53 genetically and physically interacts with Cdk1 and that p53 proapoptotic role is regulated by the cell cycle status of the cell. We propose a model in which cell cycle progression and p53 proapoptotic activity are molecularly connected to coordinate the appropriate response after DNA damage.Subject terms: Cell biology, Development, Gene regulation, Molecular biology     

Expression Quantitative Trait Loci Are Highly Sensitive to Cellular Differentiation State

Genetical genomics is a strategy for mapping gene expression variation to expression quantitative trait loci (eQTLs). We performed a genetical genomics experiment in four functionally distinct but developmentally closely related hematopoietic cell populations isolated from the BXD panel of recombinant inbred mouse strains. This analysis allowed us to analyze eQTL robustness/sensitivity across different cellular differentiation states. Although we identified a large number (365) of “static” eQTLs that were consistently active in all four cell types, we found a much larger number (1,283) of “dynamic” eQTLs showing cell-type–dependence. Of these, 140, 45, 531, and 295 were preferentially active in stem, progenitor, erythroid, and myeloid cells, respectively. A detailed investigation of those dynamic eQTLs showed that in many cases the eQTL specificity was associated with expression changes in the target gene. We found no evidence for target genes that were regulated by distinct eQTLs in different cell types, suggesting that large-scale changes within functional regulatory networks are uncommon. Our results demonstrate that heritable differences in gene expression are highly sensitive to the developmental stage of the cell population under study. Therefore, future genetical genomics studies should aim at studying multiple well-defined and highly purified cell types in order to construct as comprehensive a picture of the changing functional regulatory relationships as possible.     

Impact of Natural Genetic Variation on Gene Expression Dynamics

DNA sequence variation causes changes in gene expression, which in turn has profound effects on cellular states. These variations affect tissue development and may ultimately lead to pathological phenotypes. A genetic locus containing a sequence variation that affects gene expression is called an “expression quantitative trait locus” (eQTL). Whereas the impact of cellular context on expression levels in general is well established, a lot less is known about the cell-state specificity of eQTL. Previous studies differed with respect to how “dynamic eQTL” were defined. Here, we propose a unified framework distinguishing static, conditional and dynamic eQTL and suggest strategies for mapping these eQTL classes. Further, we introduce a new approach to simultaneously infer eQTL from different cell types. By using murine mRNA expression data from four stages of hematopoiesis and 14 related cellular traits, we demonstrate that static, conditional and dynamic eQTL, although derived from the same expression data, represent functionally distinct types of eQTL. While static eQTL affect generic cellular processes, non-static eQTL are more often involved in hematopoiesis and immune response. Our analysis revealed substantial effects of individual genetic variation on cell type-specific expression regulation. Among a total number of 3,941 eQTL we detected 2,729 static eQTL, 1,187 eQTL were conditionally active in one or several cell types, and 70 eQTL affected expression changes during cell type transitions. We also found evidence for feedback control mechanisms reverting the effect of an eQTL specifically in certain cell types. Loci correlated with hematological traits were enriched for conditional eQTL, thus, demonstrating the importance of conditional eQTL for understanding molecular mechanisms underlying physiological trait variation. The classification proposed here has the potential to streamline and unify future analysis of conditional and dynamic eQTL as well as many other kinds of QTL data.     

Robust Inference of Cell-to-Cell Expression Variations from Single- and K-Cell Profiling

Quantifying heterogeneity in gene expression among single cells can reveal information inaccessible to cell-population averaged measurements. However, the expression level of many genes in single cells fall below the detection limit of even the most sensitive technologies currently available. One proposed approach to overcome this challenge is to measure random pools of k cells (e.g., 10) to increase sensitivity, followed by computational “deconvolution” of cellular heterogeneity parameters (CHPs), such as the biological variance of single-cell expression levels. Existing approaches infer CHPs using either single-cell or k-cell data alone, and typically within a single population of cells. However, integrating both single- and k-cell data may reap additional benefits, and quantifying differences in CHPs across cell populations or conditions could reveal novel biological information. Here we present a Bayesian approach that can utilize single-cell, k-cell, or both simultaneously to infer CHPs within a single condition or their differences across two conditions. Using simulated as well as experimentally generated single- and k-cell data, we found situations where each data type would offer advantages, but using both together can improve precision and better reconcile CHP information contained in single- and k-cell data. We illustrate the utility of our approach by applying it to jointly generated single- and k-cell data to reveal CHP differences in several key inflammatory genes between resting and inflammatory cytokine-activated human macrophages, delineating differences in the distribution of ‘ON’ versus ‘OFF’ cells and in continuous variation of expression level among cells. Our approach thus offers a practical and robust framework to assess and compare cellular heterogeneity within and across biological conditions using modern multiplexed technologies.     

Tyrosine Phosphorylation of Dbl Regulates GTPase Signaling

Rho GTPases are molecular “switches” that cycle between “on” (GTP-bound) and “off” (GDP-bound) states and regulate numerous cellular activities such as gene expression, protein synthesis, cytoskeletal rearrangements, and metabolic responses. Dysregulation of GTPases is a key feature of many diseases, especially cancers. Guanine nucleotide exchange factors (GEFs) of the Dbl family are activated by mitogenic cell surface receptors and activate the Rho family GTPases Cdc42, Rac1, and RhoA. The molecular mechanisms that regulate GEFs from the Dbl family are poorly understood. Our studies reveal that Dbl is phosphorylated on tyrosine residues upon stimulation by growth factors and that this event is critical for the regulated activation of the GEF. These findings uncover a novel layer of complexity in the physiological regulation of this protein.     

Overexpression of a novel osteopetrosis-related gene CCDC154 suppresses cell proliferation by inducing G2/M arrest

Osteopetrosis, a disorder of skeletal bone, can cause death during childhood. We previously described a new spontaneous autosomal recessive osteopetrosis mouse mutant, “new toothless” (ntl). In this study, we reported for the first time the identification, cloning and characterization of the coiled-coil domain-containing 154 (CCDC154), a novel gene whose deletion of ~5 kb sequence including exons 1–6 was completely linked to the ntl mutant. The CCDC154 was conserved between mouse and human and is wildly expressed in mouse tissues. The cellular localization of CCDC154 was in the early endosomes. Overexpression of CCDC154 inhibited cell proliferation of HEK293 cells by inducing G₂/M arrest. CCDC154 also inhibited tumor cell growth, and the soft agar assay revealed a significant decrease of the colony size of Hela cells upon transfection of CCDC154. Our results indicate that CCDC154 is a novel osteopetrosis-related gene involved in cell cycle regulation and tumor suppression growth.     

A Novel TetR-Regulating Peptide Turns off rtTA-Mediated Activation of Gene Expression

Conditional regulation of gene expression is a powerful and indispensable method for analyzing gene function. The “Tet-On” system is a tool widely used for that purpose. Here, the transregulator rtTA mediates expression of a gene of interest after addition of the small molecule effector doxycycline. Although very effective in rapidly turning on gene expression, the system is hampered by the long half-life of doxycycline which makes shutting down gene expression rapidly very difficult to achieve. We isolated an rtTA-binding peptide by in vivo selection that acts as a doxycycline antagonist and leads to rapid and efficient shut down of rtTA-mediated reporter gene expression in a human cell line. This peptide represents the basis for novel effector molecules which complement the “Tet-system” by enabling the investigator to rapidly turn gene expression not just on at will, but now also off.     

O-GlcNAc and the Epigenetic Regulation of Gene Expression

O-GlcNAcylation is an abundant nutrient-driven modification linked to cellular signaling and regulation of gene expression. Utilizing precursors derived from metabolic flux, O-GlcNAc functions as a homeostatic regulator. The enzymes of O-GlcNAc cycling, OGT and O-GlcNAcase, act in mitochondria, the cytoplasm, and the nucleus in association with epigenetic “writers” and “erasers” of the histone code. Both O-GlcNAc and O-phosphate modify repeats within the RNA polymerase II C-terminal domain (CTD). By communicating with the histone and CTD codes, O-GlcNAc cycling provides a link between cellular metabolic status and the epigenetic machinery. Thus, O-GlcNAcylation is poised to influence trans-generational epigenetic inheritance.     

Determining the Variability of Lesion Size Measurements from CT Patient Data Sets Acquired under “No Change” Conditions

PURPOSE: To determine the variability of lesion size measurements in computed tomography data sets of patients imaged under a “no change” (“coffee break”) condition and to determine the impact of two reading paradigms on measurement variability. METHOD AND MATERIALS: Using data sets from 32 non-small cell lung cancer patients scanned twice within 15 minutes (“no change”), measurements were performed by five radiologists in two phases: (1) independent reading of each computed tomography dataset (timepoint): (2) a locked, sequential reading of datasets. Readers performed measurements using several sizing methods, including one-dimensional (1D) longest in-slice dimension and 3D semi-automated segmented volume. Change in size was estimated by comparing measurements performed on both timepoints for the same lesion, for each reader and each measurement method. For each reading paradigm, results were pooled across lesions, across readers, and across both readers and lesions, for each measurement method. RESULTS: The mean percent difference (± SD) when pooled across both readers and lesions for 1D and 3D measurements extracted from contours was 2.8 ± 22.2% and 23.4 ± 105.0%, respectively, for the independent reads. For the locked, sequential reads, the mean percent differences (± SD) reduced to 2.52 ± 14.2% and 7.4 ± 44.2% for the 1D and 3D measurements, respectively. CONCLUSION: Even under a “no change” condition between scans, there is variation in lesion size measurements due to repeat scans and variations in reader, lesion, and measurement method. This variation is reduced when using a locked, sequential reading paradigm compared to an independent reading paradigm.