Enhancer mutations modulate the severity of chemotherapy-induced myelosuppression.

Non-small cell lung cancer is often diagnosed at advanced stages, and many patients are still treated with classical chemotherapy. The unselective nature of chemotherapy often results in severe myelosuppression. Previous studies showed that protein-coding mutations could not fully explain the predisposition to myelosuppression. Here, we investigate the possible role of enhancer mutations in myelosuppression susceptibility. We produced transcriptome and promoter-interaction maps (using HiCap) of three blood stem-like cell lines treated with carboplatin or gemcitabine. Taking advantage of publicly available enhancer datasets, we validated HiCap results in silico and in living cells using epigenetic CRISPR technology. We also developed a network approach for interactome analysis and detection of differentially interacting genes. Differential interaction analysis provided additional information on relevant genes and pathways for myelosuppression compared with differential gene expression analysis at the bulk level. Moreover, we showed that enhancers of differentially interacting genes are highly enriched for variants associated with differing levels of myelosuppression. Altogether, our work represents a prominent example of integrative transcriptome and gene regulatory datasets analysis for the functional annotation of noncoding mutations.

SPP1/osteopontin: a driver of fibrosis and inflammation in degenerative ascending aortic aneurysm?

Degenerative ascending aortic aneurysm (AscAA) is a silent and potentially fatal disease characterized by excessive vascular inflammation and fibrosis. We aimed to characterize the cellular and molecular signature for the fibrotic type of endothelial mesenchymal transition (EndMT) that has previously been described in degenerative AscAA. Patients undergoing elective open-heart surgery for AscAA and/or aortic valve repair were recruited. Gene expression in the intima-media of the ascending aorta was measured in 22 patients with non-dilated and 24 with dilated aortas, and candidate genes were identified. Protein expression was assessed using immunohistochemistry. Interacting distal gene enhancer regions were identified using targeted chromosome conformation capture (HiCap) in untreated and LPS-treated THP1 cells, and the associated transcription factors were analyzed. Differential expression analysis identified SPP1 (osteopontin) as a key gene in the signature of fibrotic EndMT in patients with degenerative AscAA. The aortic intima-media expression of SPP1 correlated with the expression of inflammatory markers, the level of macrophage infiltration, and the aortic diameter. HiCap analysis, followed by transcription factor binding analysis, identified ETS1 as a potential regulator of SPP1 expression under inflammatory conditions. In conclusion, the present findings suggest that SPP1 may be involved in the development of the degenerative type of AscAA. KEY MESSAGES: In the original manuscript titled "SPP1/osteopontin, a driver of fibrosis and inflammation in degenerative ascending aortic aneurysm?" by David Freiholtz, Otto Bergman, Saliendra Pradhananga, Karin Lång, Flore-Anne Poujade, Carl Granath, Christian Olsson, Anders Franco-Cereceda, Pelin Sahlén, Per Eriksson, and Hanna M Björck, we present novel findings on regulatory factors on osteopontin (SPP1) expression in immune cells involved in degenerative ascending aortic aneurysms (AscAA). The central findings convey: SPP1 is a potential driver of the fibrotic endothelial-to-mesenchymal transition in AscAA. SPP1/osteopontin expression in AscAA is predominately by immune cells. ETS1 is a regulatory transcription factor of SPP1 expression in AscAA immune cells.

Current challenges in understanding the role of enhancers in disease.

Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.

Variants That Differentiate Wolf and Dog Populations Are Enriched in Regulatory Elements.

Research on the genetics of domestication most often focuses on the protein-coding exons. However, exons cover only a minor part (1-2%) of the canine genome, whereas functional mutations may be located also in regions beyond the exome, in regulatory regions. Therefore, a large proportion of phenotypical differences between dogs and wolves may remain genetically unexplained. In this study, we identified variants that have high allelic frequency differences (i.e., highly differentiated variants) between wolves and dogs across the canine genome and investigated the potential functionality. We found that the enrichment of highly differentiated variants was substantially higher in promoters than in exons and that such variants were enriched also in enhancers. Several enriched pathways were identified including oxytocin signaling, carbohydrate digestion and absorption, cancer risk, and facial and body features, many of which reflect phenotypes of potential importance during domestication, including phenotypes of the domestication syndrome. The results highlight the importance of regulatory mutations during dog domestication and motivate the functional annotation of the noncoding part of the canine genome.

Whole-genome sequencing and gene network modules predict gemcitabine/carboplatin-induced myelosuppression in non-small cell lung cancer patients.

Gemcitabine/carboplatin chemotherapy commonly induces myelosuppression, including neutropenia, leukopenia, and thrombocytopenia. Predicting patients at risk of these adverse drug reactions (ADRs) and adjusting treatments accordingly is a long-term goal of personalized medicine. This study used whole-genome sequencing (WGS) of blood samples from 96 gemcitabine/carboplatin-treated non-small cell lung cancer (NSCLC) patients and gene network modules for predicting myelosuppression. Association of genetic variants in PLINK found 4594, 5019, and 5066 autosomal SNVs/INDELs with p ≤ 1 × 10-3 for neutropenia, leukopenia, and thrombocytopenia, respectively. Based on the SNVs/INDELs we identified the toxicity module, consisting of 215 unique overlapping genes inferred from MCODE-generated gene network modules of 350, 345, and 313 genes, respectively. These module genes showed enrichment for differentially expressed genes in rat bone marrow, human bone marrow, and human cell lines exposed to carboplatin and gemcitabine (p < 0.05). Then using 80% of the patients as training data, random LASSO reduced the number of SNVs/INDELs in the toxicity module into a feasible prediction model consisting of 62 SNVs/INDELs that accurately predict both the training and the test (remaining 20%) data with high (CTCAE 3-4) and low (CTCAE 0-1) maximal myelosuppressive toxicity completely, with the receiver-operating characteristic (ROC) area under the curve (AUC) of 100%. The present study shows how WGS, gene network modules, and random LASSO can be used to develop a feasible and tested model for predicting myelosuppressive toxicity. Although the proposed model predicts myelosuppression in this study, further evaluation in other studies is required to determine its reproducibility, usability, and clinical effect.

Promoter anchored interaction landscape of THP-1 macrophages captures early immune response processes.

Macrophages are highly plastic immune cells with temporally distinct transcriptome changes upon lipopolysaccride (LPS) activation. However, to what extent transcriptome reprogramming is mediated via spatial chromatin looping is not well studied. We generated high resolution chromatin interaction maps for LPS-stimulated THP-1 macrophages (0 and 2 h) using capture Hi-C. Success of LPS stimulation was validated with transcriptome sequencing. Circa 2900 genes changed their interaction profile upon LPS stimulation and those gaining interactions were enriched for LPS response relevant processes, suggesting a substantial role for distal regulation. Immune and cardiovascular risk variants were enriched within the interacting regions, thereby providing insights into macrophage biology.

A Multi-Omics Approach to Liver Diseases: Integration of Single Nuclei Transcriptomics with Proteomics and HiCap Bulk Data in Human Liver.

The liver is the largest solid organ and a primary metabolic hub. In recent years, intact cell nuclei were used to perform single-nuclei RNA-seq (snRNA-seq) for tissues difficult to dissociate and for flash-frozen archived tissue samples to discover unknown and rare cell subpopulations. In this study, we performed snRNA-seq of a liver sample to identify subpopulations of cells based on nuclear transcriptomics. In 4282 single nuclei, we detected, on average, 1377 active genes and we identified seven major cell types. We integrated data from 94,286 distal interactions (p < 0.05) for 7682 promoters from a targeted chromosome conformation capture technique (HiCap) and mass spectrometry proteomics for the same liver sample. We observed a reasonable correlation between proteomics and in silico bulk snRNA-seq (r = 0.47) using tissue-independent gene-specific protein abundancy estimation factors. We specifically looked at genes of medical importance. The DPYD gene is involved in the pharmacogenetics of fluoropyrimidine toxicity and some of its variants are analyzed for clinical purposes. We identified a new putative polymorphic regulatory element, which may contribute to variation in toxicity. Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer and we investigated all known risk genes. We identified a complex regulatory landscape for the SLC2A2 gene with 16 candidate enhancers. Three of them harbor somatic motif breaking and other mutations in HCC in the Pan Cancer Analysis of Whole Genomes dataset and are candidates to contribute to malignancy. Our results highlight the potential of a multi-omics approach in the study of human diseases.

High-Resolution Regulatory Maps Connect Vascular Risk Variants to Disease-Related Pathways.

BACKGROUND: Genetic variant landscape of coronary artery disease is dominated by noncoding variants among which many occur within putative enhancers regulating the expression levels of relevant genes. It is crucial to assign the genetic variants to their correct genes both to gain insights into perturbed functions and better assess the risk of disease. METHODS: In this study, we generated high-resolution genomic interaction maps (≈750 bases) in aortic endothelial, smooth muscle cells and THP-1 (human leukemia monocytic cell line) macrophages stimulated with lipopolysaccharide using Hi-C coupled with sequence capture targeting 25 429 features, including variants associated with coronary artery disease. We also sequenced their transcriptomes and mapped putative enhancers using chromatin immunoprecipitation with an antibody against H3K27Ac. RESULTS: The regions interacting with promoters showed strong enrichment for enhancer elements and validated several previously known interactions and enhancers. We detected interactions for 727 risk variants obtained by genome-wide association studies and identified novel, as well as established genes and functions associated with cardiovascular diseases. We were able to assign potential target genes for additional 398 genome-wide association studies variants using haplotype information, thereby identifying additional relevant genes and functions. Importantly, we discovered that a subset of risk variants interact with multiple promoters and their expression levels were strongly correlated. CONCLUSIONS: In summary, we present a catalog of candidate genes regulated by coronary artery disease-related variants and think that it will be an invaluable resource to further the investigation of cardiovascular pathologies and disease.

Visualization and analysis of gene expression in tissue sections by spatial transcriptomics.

Analysis of the pattern of proteins or messengerRNAs (mRNAs) in histological tissue sections is a cornerstone in biomedical research and diagnostics. This typically involves the visualization of a few proteins or expressed genes at a time. We have devised a strategy, which we call "spatial transcriptomics," that allows visualization and quantitative analysis of the transcriptome with spatial resolution in individual tissue sections. By positioning histological sections on arrayed reverse transcription primers with unique positional barcodes, we demonstrate high-quality RNA-sequencing data with maintained two-dimensional positional information from the mouse brain and human breast cancer. Spatial transcriptomics provides quantitative gene expression data and visualization of the distribution of mRNAs within tissue sections and enables novel types of bioinformatics analyses, valuable in research and diagnostics.