Suv39h-catalyzed H3K9me3 is critical for euchromatic genome organization and the maintenance of gene transcription.

H3K9me3-dependent heterochromatin is critical for the silencing of repeat-rich pericentromeric regions and also has key roles in repressing lineage-inappropriate protein-coding genes in differentiation and development. Here, we investigate the molecular consequences of heterochromatin loss in cells deficient in both SUV39H1 and SUV39H2 (Suv39DKO), the major mammalian histone methyltransferase enzymes that catalyze heterochromatic H3K9me3 deposition. We reveal a paradoxical repression of protein-coding genes in Suv39DKO cells, with these differentially expressed genes principally in euchromatic (Tn5-accessible, H3K4me3- and H3K27ac-marked) rather than heterochromatic (H3K9me3-marked) or polycomb (H3K27me3-marked) regions. Examination of the three-dimensional (3D) nucleome reveals that transcriptomic dysregulation occurs in euchromatic regions close to the nuclear periphery in 3D space. Moreover, this transcriptomic dysregulation is highly correlated with altered 3D genome organization in Suv39DKO cells. Together, our results suggest that the nuclear lamina-tethering of Suv39-dependent H3K9me3 domains provides an essential scaffold to support euchromatic genome organization and the maintenance of gene transcription for healthy cellular function.

Increasing histone acetylation improves sociability and restores learning and memory in KAT6B-haploinsufficient mice.

Mutations in genes encoding chromatin modifiers are enriched among mutations causing intellectual disability. The continuing development of the brain postnatally, coupled with the inherent reversibility of chromatin modifications, may afford an opportunity for therapeutic intervention following a genetic diagnosis. Development of treatments requires an understanding of protein function and models of the disease. Here, we provide a mouse model of Say-Barber-Biesecker-Young-Simpson syndrome (SBBYSS) (OMIM 603736) and demonstrate proof-of-principle efficacy of postnatal treatment. SBBYSS results from heterozygous mutations in the KAT6B (MYST4/MORF/QFK) gene and is characterized by intellectual disability and autism-like behaviors. Using human cells carrying SBBYSS-specific KAT6B mutations and Kat6b heterozygous mice (Kat6b+/-), we showed that KAT6B deficiency caused a reduction in histone H3 lysine 9 acetylation. Kat6b+/- mice displayed learning, memory, and social deficits, mirroring SBBYSS individuals. Treatment with a histone deacetylase inhibitor, valproic acid, or an acetyl donor, acetyl-carnitine (ALCAR), elevated histone acetylation levels in the human cells with SBBYSS mutations and in brain and blood cells of Kat6b+/- mice and partially reversed gene expression changes in Kat6b+/- cortical neurons. Both compounds improved sociability in Kat6b+/- mice, and ALCAR treatment restored learning and memory. These data suggest that a subset of SBBYSS individuals may benefit from postnatal therapeutic interventions.

ING4 and ING5 are essential for histone H3 lysine 14 acetylation and epicardial cell lineage development.

Inhibitor of growth 4 and 5 (ING4, ING5) are structurally similar chromatin-binding proteins in the KAT6A, KAT6B and KAT7 histone acetyltransferase protein complexes. Heterozygous mutations in the KAT6A or KAT6B gene cause human disorders with cardiac defects, but the contribution of their chromatin-adaptor proteins to development is unknown. We found that Ing5-/- mice had isolated cardiac ventricular septal defects. Ing4-/-Ing5-/- embryos failed to undergo chorioallantoic fusion and arrested in development at embryonic day 8.5, displaying loss of histone H3 lysine 14 acetylation, reduction in H3 lysine 23 acetylation levels and reduced developmental gene expression. Embryonic day 12.5 Ing4+/-Ing5-/- hearts showed a paucity of epicardial cells and epicardium-derived cells, failure of myocardium compaction, and coronary vasculature defects, accompanied by reduced expression of epicardium genes. Cell adhesion gene expression and proepicardium outgrowth were defective in the ING4- and ING5-deficient state. Our findings suggest that ING4 and ING5 are essential for heart development and promote epicardium and epicardium-derived cell fates and imply mutation of the human ING5 gene as a possible cause of isolated ventricular septal defects.

The histone acetyltransferase KAT6B is required for hematopoietic stem cell development and function.

The histone lysine acetyltransferase KAT6B (MYST4, MORF, QKF) is the target of recurrent chromosomal translocations causing hematological malignancies with poor prognosis. Using Kat6b germline deletion and overexpression in mice, we determined the role of KAT6B in the hematopoietic system. We found that KAT6B sustained the fetal hematopoietic stem cell pool but did not affect viability or differentiation. KAT6B was essential for normal levels of histone H3 lysine 9 (H3K9) acetylation but not for a previously proposed target, H3K23. Compound heterozygosity of Kat6b and the closely related gene, Kat6a, abolished hematopoietic reconstitution after transplantation. KAT6B and KAT6A cooperatively promoted transcription of genes regulating hematopoiesis, including the Hoxa cluster, Pbx1, Meis1, Gata family, Erg, and Flt3. In conclusion, we identified the hematopoietic processes requiring Kat6b and showed that KAT6B and KAT6A synergistically promoted HSC development, function, and transcription. Our findings are pertinent to current clinical trials testing KAT6A/B inhibitors as cancer therapeutics.

Systemic inflammatory response syndrome triggered by blood-borne pathogens induces prolonged dendritic cell paralysis and immunosuppression.

Blood-borne pathogens can cause systemic inflammatory response syndrome (SIRS) followed by protracted, potentially lethal immunosuppression. The mechanisms responsible for impaired immunity post-SIRS remain unclear. We show that SIRS triggered by pathogen mimics or malaria infection leads to functional paralysis of conventional dendritic cells (cDCs). Paralysis affects several generations of cDCs and impairs immunity for 3-4 weeks. Paralyzed cDCs display distinct transcriptomic and phenotypic signatures and show impaired capacity to capture and present antigens in vivo. They also display altered cytokine production patterns upon stimulation. The paralysis program is not initiated in the bone marrow but during final cDC differentiation in peripheral tissues under the influence of local secondary signals that persist after resolution of SIRS. Vaccination with monoclonal antibodies that target cDC receptors or blockade of transforming growth factor ß partially overcomes paralysis and immunosuppression. This work provides insights into the mechanisms of paralysis and describes strategies to restore immunocompetence post-SIRS.

Identification of aberrant luminal progenitors and mTORC1 as a potential breast cancer prevention target in BRCA2 mutation carriers.

Inheritance of a BRCA2 pathogenic variant conveys a substantial life-time risk of breast cancer. Identification of the cell(s)-of-origin of BRCA2-mutant breast cancer and targetable perturbations that contribute to transformation remains an unmet need for these individuals who frequently undergo prophylactic mastectomy. Using preneoplastic specimens from age-matched, premenopausal females, here we show broad dysregulation across the luminal compartment in BRCA2mut/+ tissue, including expansion of aberrant ERBB3lo luminal progenitor and mature cells, and the presence of atypical oestrogen receptor (ER)-positive lesions. Transcriptional profiling and functional assays revealed perturbed proteostasis and translation in ERBB3lo progenitors in BRCA2mut/+ breast tissue, independent of ageing. Similar molecular perturbations marked tumours bearing BRCA2-truncating mutations. ERBB3lo progenitors could generate both ER+ and ER- cells, potentially serving as cells-of-origin for ER-positive or triple-negative cancers. Short-term treatment with an mTORC1 inhibitor substantially curtailed tumorigenesis in a preclinical model of BRCA2-deficient breast cancer, thus uncovering a potential prevention strategy for BRCA2 mutation carriers.

Unraveling the diversity and functions of tissue-resident plasma cells.

Antibody-secreting plasma cells (PCs) are generated in secondary lymphoid organs but are reported to reside in an emerging range of anatomical sites. Analysis of the transcriptome of different tissue-resident (Tr)PC populations revealed that they each have their own transcriptional signature indicative of functional adaptation to the host tissue environment. In contrast to expectation, all TrPCs were extremely long-lived, regardless of their organ of residence, with longevity influenced by intrinsic factors like the immunoglobulin isotype. Analysis at single-cell resolution revealed that the bone marrow is unique in housing a compendium of PCs generated all over the body that retain aspects of the transcriptional program indicative of their tissue of origin. This study reveals that extreme longevity is an intrinsic property of TrPCs whose transcriptome is imprinted by signals received both at the site of induction and within the tissue of residence.

Dividing out quantification uncertainty allows efficient assessment of differential transcript expression with edgeR.

Differential expression analysis of RNA-seq is one of the most commonly performed bioinformatics analyses. Transcript-level quantifications are inherently more uncertain than gene-level read counts because of ambiguous assignment of sequence reads to transcripts. While sequence reads can usually be assigned unambiguously to a gene, reads are very often compatible with multiple transcripts for that gene, particularly for genes with many isoforms. Software tools designed for gene-level differential expression do not perform optimally on transcript counts because the read-to-transcript ambiguity (RTA) disrupts the mean-variance relationship normally observed for gene level RNA-seq data and interferes with the efficiency of the empirical Bayes dispersion estimation procedures. The pseudoaligners kallisto and Salmon provide bootstrap samples from which quantification uncertainty can be assessed. We show that the overdispersion arising from RTA can be elegantly estimated by fitting a quasi-Poisson model to the bootstrap counts for each transcript. The technical overdispersion arising from RTA can then be divided out of the transcript counts, leading to scaled counts that can be input for analysis by established gene-level software tools with full statistical efficiency. Comprehensive simulations and test data show that an edgeR analysis of the scaled counts is more powerful and efficient than previous differential transcript expression pipelines while providing correct control of the false discovery rate. Simulations explore a wide range of scenarios including the effects of paired vs single-end reads, different read lengths and different numbers of replicates.

Three-dimensional genome architecture coordinates key regulators of lineage specification in mammary epithelial cells.

Although lineage-specific genes have been identified in the mammary gland, little is known about the contribution of the 3D genome organization to gene regulation in the epithelium. Here, we describe the chromatin landscape of the three major epithelial subsets through integration of long- and short-range chromatin interactions, accessibility, histone modifications, and gene expression. While basal genes display exquisite lineage specificity via distal enhancers, luminal-specific genes show widespread promoter priming in basal cells. Cell specificity in luminal progenitors is largely mediated through extensive chromatin interactions with super-enhancers in gene-body regions in addition to interactions with polycomb silencer elements. Moreover, lineage-specific transcription factors appear to be controlled through cell-specific chromatin interactivity. Finally, chromatin accessibility rather than interactivity emerged as a defining feature of the activation of quiescent basal stem cells. This work provides a comprehensive resource for understanding the role of higher-order chromatin interactions in cell-fate specification and differentiation in the adult mouse mammary gland.

Unraveling the timeline of gene expression: A pseudotemporal trajectory analysis of single-cell RNA sequencing data.

Background: Single-cell RNA sequencing (scRNA-seq) technologies have rapidly developed in recent years. The droplet-based single cell platforms enable the profiling of gene expression in tens of thousands of cells per sample. The goal of a typical scRNA-seq analysis is to identify different cell subpopulations and their respective marker genes. Additionally, trajectory analysis can be used to infer the developmental or differentiation trajectories of cells. Methods: This article demonstrates a comprehensive workflow for performing trajectory inference and time course analysis on a multi-sample single-cell RNA-seq experiment of the mouse mammary gland. The workflow uses open-source R software packages and covers all steps of the analysis pipeline, including quality control, doublet prediction, normalization, integration, dimension reduction, cell clustering, trajectory inference, and pseudo-bulk time course analysis. Sample integration and cell clustering follows the Seurat pipeline while the trajectory inference is conducted using the monocle3 package. The pseudo-bulk time course analysis uses the quasi-likelihood framework of edgeR. Results: Cells are ordered and positioned along a pseudotime trajectory that represented a biological process of cell differentiation and development. The study successfully identified genes that were significantly associated with pseudotime in the mouse mammary gland. Conclusions: The demonstrated workflow provides a valuable resource for researchers conducting scRNA-seq analysis using open-source software packages. The study successfully demonstrated the usefulness of trajectory analysis for understanding the developmental or differentiation trajectories of cells. This analysis can be applied to various biological processes such as cell development or disease progression, and can help identify potential biomarkers or therapeutic targets.

Benchmarking long-read RNA-sequencing analysis tools using in silico mixtures.

The lack of benchmark data sets with inbuilt ground-truth makes it challenging to compare the performance of existing long-read isoform detection and differential expression analysis workflows. Here, we present a benchmark experiment using two human lung adenocarcinoma cell lines that were each profiled in triplicate together with synthetic, spliced, spike-in RNAs (sequins). Samples were deeply sequenced on both Illumina short-read and Oxford Nanopore Technologies long-read platforms. Alongside the ground-truth available via the sequins, we created in silico mixture samples to allow performance assessment in the absence of true positives or true negatives. Our results show that StringTie2 and bambu outperformed other tools from the six isoform detection tools tested, DESeq2, edgeR and limma-voom were best among the five differential transcript expression tools tested and there was no clear front-runner for performing differential transcript usage analysis between the five tools compared, which suggests further methods development is needed for this application.

Venetoclax, alone and in combination with the BH3 mimetic S63845, depletes HIV-1 latently infected cells and delays rebound in humanized mice.

HIV-1 persists indefinitely in people living with HIV (PLWH) on antiretroviral therapy (ART). If ART is stopped, the virus rapidly rebounds from long-lived latently infected cells. Using a humanized mouse model of HIV-1 infection and CD4+ T cells from PLWH on ART, we investigate whether antagonizing host pro-survival proteins can prime latent cells to die and facilitate HIV-1 clearance. Venetoclax, a pro-apoptotic inhibitor of Bcl-2, depletes total and intact HIV-1 DNA in CD4+ T cells from PLWH ex vivo. This venetoclax-sensitive population is enriched for cells with transcriptionally higher levels of pro-apoptotic BH3-only proteins. Furthermore, venetoclax delays viral rebound in a mouse model of persistent HIV-1 infection, and the combination of venetoclax with the Mcl-1 inhibitor S63845 achieves a longer delay in rebound compared with either intervention alone. Thus, selective inhibition of pro-survival proteins can induce death of HIV-1-infected cells that persist on ART, extending time to viral rebound.

Modeling group heteroscedasticity in single-cell RNA-seq pseudo-bulk data.

Group heteroscedasticity is commonly observed in pseudo-bulk single-cell RNA-seq datasets and its presence can hamper the detection of differentially expressed genes. Since most bulk RNA-seq methods assume equal group variances, we introduce two new approaches that account for heteroscedastic groups, namely voomByGroup and voomWithQualityWeights using a blocked design (voomQWB). Compared to current gold-standard methods that do not account for group heteroscedasticity, we show results from simulations and various experiments that demonstrate the superior performance of voomByGroup and voomQWB in terms of error control and power when group variances in pseudo-bulk single-cell RNA-seq data are unequal.

Neither random nor censored: estimating intensity-dependent probabilities for missing values in label-free proteomics.

MOTIVATION: Mass spectrometry proteomics is a powerful tool in biomedical research but its usefulness is limited by the frequent occurrence of missing values in peptides that cannot be reliably quantified (detected) for particular samples. Many analysis strategies have been proposed for missing values where the discussion often focuses on distinguishing whether values are missing completely at random (MCAR), missing at random (MAR) or missing not at random (MNAR). RESULTS: Statistical models and algorithms are proposed for estimating the detection probabilities and for evaluating how much statistical information can or cannot be recovered from the missing value pattern. The probability that an intensity is detected is shown to be accurately modeled as a logit-linear function of the underlying intensity, showing that missing value process is intermediate between MAR and censoring. The detection probability asymptotes to 100% for high intensities, showing that missing values unrelated to intensity are rare. The rule applies globally to each dataset and is appropriate for both high and lowly expressed peptides. A probability model is developed that allows the distribution of unobserved intensities to be inferred from the observed values. The detection probability model is incorporated into a likelihood-based approach for assessing differential expression and successfully recovers statistical power compared to omitting the missing values from the analysis. In contrast, imputation methods are shown to perform poorly, either reducing statistical power or increasing the false discovery rate to unacceptable levels. AVAILABILITY AND IMPLEMENTATION: Data and code to reproduce the results shown in this article are available from https://mengbo-li.github.io/protDP/.

Overexpression of Lmo2 initiates T-lymphoblastic leukemia via impaired thymocyte competition.

Cell competition has recently emerged as an important tumor suppressor mechanism in the thymus that inhibits autonomous thymic maintenance. Here, we show that the oncogenic transcription factor Lmo2 causes autonomous thymic maintenance in transgenic mice by inhibiting early T cell differentiation. This autonomous thymic maintenance results in the development of self-renewing preleukemic stem cells (pre-LSCs) and subsequent leukemogenesis, both of which are profoundly inhibited by restoration of thymic competition or expression of the antiapoptotic factor BCL2. Genomic analyses revealed the presence of Notch1 mutations in pre-LSCs before subsequent loss of tumor suppressors promotes the transition to overt leukemogenesis. These studies demonstrate a critical role for impaired cell competition in the development of pre-LSCs in a transgenic mouse model of T cell acute lymphoblastic leukemia (T-ALL), implying that this process plays a role in the ontogeny of human T-ALL.

Stem cell plasticity, acetylation of H3K14, and de novo gene activation rely on KAT7.

In the conventional model of transcriptional activation, transcription factors bind to response elements and recruit co-factors, including histone acetyltransferases. Contrary to this model, we show that the histone acetyltransferase KAT7 (HBO1/MYST2) is required genome wide for histone H3 lysine 14 acetylation (H3K14ac). Examining neural stem cells, we find that KAT7 and H3K14ac are present not only at transcribed genes but also at inactive genes, intergenic regions, and in heterochromatin. KAT7 and H3K14ac were not required for the continued transcription of genes that were actively transcribed at the time of loss of KAT7 but indispensable for the activation of repressed genes. The absence of KAT7 abrogates neural stem cell plasticity, diverse differentiation pathways, and cerebral cortex development. Re-expression of KAT7 restored stem cell developmental potential. Overexpression of KAT7 enhanced neuron and oligodendrocyte differentiation. Our data suggest that KAT7 prepares chromatin for transcriptional activation and is a prerequisite for gene activation.

Survey of activation-induced genome architecture reveals a novel enhancer of Myc.

The transcription factor Myc is critically important in driving cell proliferation, a function that is frequently dysregulated in cancer. To avoid this dysregulation Myc is tightly controlled by numerous layers of regulation. One such layer is the use of distal regulatory enhancers to drive Myc expression. Here, using chromosome conformation capture to examine B cells of the immune system in the first hours after their activation, we reveal a previously unidentified enhancer of Myc. The interactivity of this enhancer coincides with a dramatic, but discrete, spike in Myc expression 3 h post-activation. However, genetic deletion of this region, has little impact on Myc expression, Myc protein level or in vitro and in vivo cell proliferation. Examination of the enhancer deleted regulatory landscape suggests that enhancer redundancy likely sustains Myc expression. This work highlights not only the importance of temporally examining enhancers, but also the complexity and dynamics of the regulation of critical genes such as Myc.

Activation of stably silenced genes by recruitment of a synthetic de-methylating module.

Stably silenced genes that display a high level of CpG dinucleotide methylation are refractory to the current generation of dCas9-based activation systems. To counter this, we create an improved activation system by coupling the catalytic domain of DNA demethylating enzyme TET1 with transcriptional activators (TETact). We show that TETact demethylation-coupled activation is able to induce transcription of suppressed genes, both individually and simultaneously in cells, and has utility across a number of cell types. Furthermore, we show that TETact can effectively reactivate embryonic haemoglobin genes in non-erythroid cells. We anticipate that TETact will expand the existing CRISPR toolbox and be valuable for functional studies, genetic screens and potential therapeutics.

Loss of TIP60 (KAT5) abolishes H2AZ lysine 7 acetylation and causes p53, INK4A, and ARF-independent cell cycle arrest.

Histone acetylation is essential for initiating and maintaining a permissive chromatin conformation and gene transcription. Dysregulation of histone acetylation can contribute to tumorigenesis and metastasis. Using inducible cre-recombinase and CRISPR/Cas9-mediated deletion, we investigated the roles of the histone lysine acetyltransferase TIP60 (KAT5/HTATIP) in human cells, mouse cells, and mouse embryos. We found that loss of TIP60 caused complete cell growth arrest. In the absence of TIP60, chromosomes failed to align in a metaphase plate during mitosis. In some TIP60 deleted cells, endoreplication occurred instead. In contrast, cell survival was not affected. Remarkably, the cell growth arrest caused by loss of TIP60 was independent of the tumor suppressors p53, INK4A and ARF. TIP60 was found to be essential for the acetylation of H2AZ, specifically at lysine 7. The mRNA levels of 6236 human and 8238 mouse genes, including many metabolism genes, were dependent on TIP60. Among the top 50 differentially expressed genes, over 90% were downregulated in cells lacking TIP60, supporting a role for TIP60 as a key co-activator of transcription. We propose a primary role of TIP60 in H2AZ lysine 7 acetylation and transcriptional activation, and that this fundamental role is essential for cell proliferation. Growth arrest independent of major tumor suppressors suggests TIP60 as a potential anti-cancer drug target.