Viral reprogramming of host transcription initiation.

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

Transcriptomic analysis identifies B-lymphocyte kinase as a therapeutic target for desmoplastic small round cell tumor cancer stem cell-like cells.

Desmoplastic small round cell tumor (DSRCT) is an aggressive pediatric cancer caused by the EWSR1-WT1 fusion oncoprotein. The tumor is refractory to treatment with a 5-year survival rate of only 15-25%, necessitating the development of novel therapeutics, especially those able to target chemoresistant subpopulations. Novel in vitro cancer stem cell-like (CSC-like) culture conditions increase the expression of stemness markers (SOX2, NANOG) and reduce DSRCT cell line susceptibility to chemotherapy while maintaining the ability of DSRCT cells to form xenografts. To gain insights into this chemoresistant model, RNA-seq was performed to elucidate transcriptional alterations between DSRCT cells grown in CSC-like spheres and normal 2-dimensional adherent state. Commonly upregulated and downregulated genes were identified and utilized in pathway analysis revealing upregulation of pathways related to chromatin assembly and disassembly and downregulation of pathways including cell junction assembly and extracellular matrix organization. Alterations in chromatin assembly suggest a role for epigenetics in the DSRCT CSC-like state, which was further investigated with ATAC-seq, identifying over 10,000 differentially accessible peaks, including 4444 sphere accessible peaks and 6,120 adherent accessible peaks. Accessible regions were associated with higher gene expression, including increased accessibility of the CSC marker SOX2 in CSC-like culture conditions. These analyses were further utilized to identify potential CSC therapeutic targets, leading to the identification of B-lymphocyte kinase (BLK) as a CSC-enriched, EWSR1-WT1-regulated, druggable target. BLK inhibition and knockdown reduced CSC-like properties, including abrogation of tumorsphere formation and stemness marker expression. Importantly, BLK knockdown reduced DSRCT CSC-like cell chemoresistance, making its inhibition a promising target for future combination therapy.

Loss of feedback regulation between FAM3B and androgen receptor driving prostate cancer progression.

BACKGROUND: Although the fusion of the transmembrane serine protease 2 gene (TMPRSS2) with the erythroblast transformation-specific-related gene (ERG), or TMPRSS2-ERG, occurs frequently in prostate cancer, its impact on clinical outcomes remains controversial. Roughly half of TMPRSS2-ERG fusions occur through intrachromosomal deletion of interstitial genes and the remainder via insertional chromosomal rearrangements. Because prostate cancers with deletion-derived TMPRSS2-ERG fusions are more aggressive than those with insertional fusions, we investigated the impact of interstitial gene loss on prostate cancer progression. METHODS: We conducted an unbiased analysis of transcriptome data from large collections of prostate cancer samples and employed diverse in vitro and in vivo models combined with genetic approaches to characterize the interstitial gene loss that imposes the most important impact on clinical outcome. RESULTS: This analysis identified FAM3B as the top-ranked interstitial gene whose loss is associated with a poor prognosis. The association between FAM3B loss and poor clinical outcome extended to fusion-negative prostate cancers where FAM3B downregulation occurred through epigenetic imprinting. Importantly, FAM3B loss drives disease progression in prostate cancer. FAM3B acts as an intermediator of a self-governing androgen receptor feedback loop. Specifically, androgen receptor upregulates FAM3B expression by binding to an intronic enhancer to induce an enhancer RNA and facilitate enhancer-promoter looping. FAM3B, in turn, attenuates androgen receptor signaling. CONCLUSION: Loss of FAM3B in prostate cancer, whether through the TMPRSS2-ERG translocation or epigenetic imprinting, causes an exit from this autoregulatory loop to unleash androgen receptor activity and prostate cancer progression. These findings establish FAM3B loss as a new driver of prostate cancer progression and support the utility of FAM3B loss as a biomarker to better define aggressive prostate cancer.

Maternal Western diet programs cardiometabolic dysfunction and hypothalamic inflammation via epigenetic mechanisms predominantly in the male offspring.

OBJECTIVE: Maternal exposure during pregnancy is a strong determinant of offspring health outcomes. Such exposure induces changes in the offspring epigenome resulting in gene expression and functional changes. In this study, we investigated the effect of maternal Western hypercaloric diet (HCD) programming during the perinatal period on neuronal plasticity and cardiometabolic health in adult offspring. METHODS: C57BL/6J dams were fed HCD for 1 month prior to mating with regular diet (RD) sires and kept on the same diet throughout pregnancy and lactation. At weaning, offspring were maintained on either HCD or RD for 3 months resulting in 4 treatment groups that underwent cardiometabolic assessments. DNA and RNA were extracted from the hypothalamus to perform whole genome methylation, mRNA, and miRNA sequencing followed by bioinformatic analyses. RESULTS: Maternal programming resulted in male-specific hypertension and hyperglycemia, with both males and females showing increased sympathetic tone to the vasculature. Surprisingly, programmed male offspring fed HCD in adulthood exhibited lower glucose levels, less insulin resistance, and leptin levels compared to non-programmed HCD-fed male mice. Hypothalamic genes involved in inflammation and type 2 diabetes were targeted by differentially expressed miRNA, while genes involved in glial and astrocytic differentiation were differentially methylated in programmed male offspring. These data were supported by our findings of astrogliosis, microgliosis and increased microglial activation in programmed males in the paraventricular nucleus (PVN). Programming induced a protective effect in male mice fed HCD in adulthood, resulting in lower protein levels of hypothalamic TGFß2, NF-κB2, NF-κBp65, Ser-pIRS1, and GLP1R compared to non-programmed HCD-fed males. Although TGFß2 was upregulated in male mice exposed to HCD pre- or post-natally, only blockade of the brain TGFß receptor in RD-HCD mice improved glucose tolerance and a trend to weight loss. CONCLUSIONS: Our study shows that maternal HCD programs neuronal plasticity in the offspring and results in male-specific hypertension and hyperglycemia associated with hypothalamic inflammation in mechanisms and pathways distinct from post-natal HCD exposure. Together, our data unmask a compensatory role of HCD programming, likely via priming of metabolic pathways to handle excess nutrients in a more efficient way.

A novel role of TRIM28 B box domain in L1 retrotransposition and ORF2p-mediated cDNA synthesis.

The long interspersed element 1 (LINE-1 or L1) integration is affected by many cellular factors through various mechanisms. Some of these factors are required for L1 amplification, while others either suppress or enhance specific steps during L1 propagation. Previously, TRIM28 has been identified to suppress transposable elements, including L1 expression via its canonical role in chromatin remodeling. Here, we report that TRIM28 through its B box domain increases L1 retrotransposition and facilitates shorter cDNA and L1 insert generation in cultured cells. Consistent with the latter, we observe that tumor specific L1 inserts are shorter in endometrial, ovarian, and prostate tumors with higher TRIM28 mRNA expression than in those with lower TRIM28 expression. We determine that three amino acids in the B box domain that are involved in TRIM28 multimerization are critical for its effect on both L1 retrotransposition and cDNA synthesis. We provide evidence that B boxes from the other two members in the Class VI TRIM proteins, TRIM24 and TRIM33, also increase L1 retrotransposition. Our findings could lead to a better understanding of the host/L1 evolutionary arms race in the germline and their interplay during tumorigenesis.

SpliceTools, a suite of downstream RNA splicing analysis tools to investigate mechanisms and impact of alternative splicing.

As a fundamental aspect of normal cell signaling and disease states, there is great interest in determining alternative splicing (AS) changes in physiologic, pathologic, and pharmacologic settings. High throughput RNA sequencing and specialized software to detect AS has greatly enhanced our ability to determine transcriptome-wide splicing changes. Despite the richness of this data, deriving meaning from sometimes thousands of AS events is a substantial bottleneck for most investigators. We present SpliceTools, a suite of data processing modules that arms investigators with the ability to quickly produce summary statistics, mechanistic insights, and functional significance of AS changes through command line or through an online user interface. Utilizing RNA-seq datasets for 186 RNA binding protein knockdowns, nonsense mediated RNA decay inhibition, and pharmacologic splicing inhibition, we illustrate the utility of SpliceTools to distinguish splicing disruption from regulated transcript isoform changes, we show the broad transcriptome footprint of the pharmacologic splicing inhibitor, indisulam, we illustrate the utility in uncovering mechanistic underpinnings of splicing inhibition, we identify predicted neo-epitopes in pharmacologic splicing inhibition, and we show the impact of splicing alterations induced by indisulam on cell cycle progression. Together, SpliceTools puts rapid and easy downstream analysis at the fingertips of any investigator studying AS.

Metabolic pathways enriched according to ERG status are associated with biochemical recurrence in Hispanic/Latino patients with prostate cancer.

BACKGROUND: The role of ERG-status molecular subtyping in prognosis of prostate cancer (PCa) is still under debate. In this study, we identified differentially expressed genes (DEGs) according to ERG-status to explore their enriched pathways and implications in prognosis in Hispanic/Latino PCa patients. METHODS: RNA from 78 Hispanic PCa tissues from radical prostatectomies (RP) were used for RNA-sequencing. ERGhigh /ERGlow tumor groups were determined based on the 1.5-fold change median expression in non-tumor samples. DEGs with a False Discovery Rate (FDR) < 0.01 and a fold change >2 were identified between ERGhigh and ERGlow tumors and submitted to enrichment analysis in MetaCore. Survival and association analyses were performed to evaluate biochemical recurrence (BCR)-free survival. RESULTS: The identification of 150 DEGs between ERGhigh and ERGlow tumors revealed clustering of most of the non-BCR cases (60%) into de ERGhigh group and most of the BCR cases (60.8%) in ERGlow group. Kaplan-Meier survival curves showed a worst BCR-free survival for ERGlow patients, and a significant reduced risk of BCR was observed for ERGhigh cases (OR = 0.29 (95%CI, 0.10-0.8)). Enrichment pathway analysis identified metabolic-related pathways, such as the renin-angiotensin system and angiotensin maturation system, the linoleic acid metabolism, and polyamines metabolism in these ERG groups. CONCLUSIONS: ERGlow tumor cases were associated with poor BCR-free survival in our Hispanic/Latino patients, with metabolism-related pathways altered in the BCR progression. IMPACT: Our findings suggest the need to dissect the role of diet, metabolism, and lifestyle as risk factors for more aggressive PCa subtypes.

Apremilast in Recalcitrant Cutaneous Dermatomyositis: A Nonrandomized Controlled Trial.

Importance: Cutaneous disease in dermatomyositis has no standardized treatment approach and so presents a challenging task for patients and clinicians. Objective: To study the efficacy and safety of apremilast as an add-on therapy in patients with recalcitrant cutaneous dermatomyositis. Design, Setting, and Participants: This phase 2a, open-label, single-arm nonrandomized controlled trial was conducted at a single center from June 2018 to June 2021. Participants were 8 patients with recalcitrant cutaneous dermatomyositis, defined by a cutaneous disease activity severity index (CDASI) score greater than 5 despite treatment with steroids, steroid-sparing agents, or both. Data were analyzed from June 2018 to June 2021. Interventions: Apremilast 30 mg orally twice daily was added to ongoing treatment regimens. Main Outcomes and Measures: The primary outcome was the overall response rate (ORR) at 3 months. Key secondary outcomes were the safety and toxicity of apremilast and the durability of response at 6 months. The CDASI, muscle score, dermatology life quality index (DLQI), and depression assessments were performed at baseline and regularly until month 7. Skin biopsies were performed at baseline and 3 months after apremilast (defined as 3 months into active apremilast therapy) and tested for gene expression profiling and immunohistochemical stains. Adverse events were assessed using the Common Terminology Criteria for Adverse Events version 5.0. Results: Among 8 patients with recalcitrant cutaneous dermatomyositis (all women; mean [SD] age, 54 [15.9] years), a response was found at 3 months after apremilast among 7 patients (ORR, 87.5%). The mean (SD) decrease in CDASI was 12.9 (6.3) points at 3 months (P < .001). Apremilast was well tolerated, with no grade 3 or higher adverse events. Sequencing of RNA was performed on skin biopsies taken from 7 patients at baseline and at 3 months after therapy. Appropriate negative (ie, no primary antibody) and positive (ie, tonsil and spleen) controls were stained in parallel with each set of slides studied. Of 39 076 expressed genes, there were 195 whose expression changed 2-fold or more at P < .01 (123 downregulated and 72 upregulated genes). Gene set enrichment analysis identified 13 pathways in which apremilast was associated with downregulated expression, notably signal transducers and activators of transcription 1 (STAT1), STAT3, interleukin 4 (IL-4), IL-6, IL-12, IL-23, interferon γ (IFNγ), and tumor necrosis factor α (TNFα) pathways. In immunohistochemical staining, there was a mean (SD) decrease in phosphorylation levels STAT1 (22.3% [28.3%] positive cells) and STAT3 (13.4% [11.6%] positive cells) at the protein level, a downstream signaling pathway for the downregulated cytokines. Conclusions and Relevance: These findings suggest that apremilast was a safe and efficacious add-on treatment in recalcitrant dermatomyositis, with an overall response rate of 87.5% and associations with downregulation of multiple inflammatory pathways. Trial Registration: ClinicalTrials.gov Identifier: NCT03529955.

SCIFER: approach for analysis of LINE-1 mRNA expression in single cells at a single locus resolution.

BACKGROUND: Endogenous expression of L1 mRNA is the first step in an L1-initiated mutagenesis event. However, the contribution of individual cell types to patterns of organ-specific L1 mRNA expression remains poorly understood, especially at single-locus resolution. We introduce a method to quantify expression of mobile elements at the single-locus resolution in scRNA-Seq datasets called Single Cell Implementation to Find Expressed Retrotransposons (SCIFER). SCIFER aligns scRNA-Seq reads uniquely to the genome and extracts alignments from single cells by cell-specific barcodes. In contrast to the alignment performed using default parameters, this alignment strategy increases accuracy of L1 locus identification by retaining only reads that are uniquely mapped to individual L1 loci. L1 loci expressed in single cells are unambiguously identified using a list of L1 loci manually validated to be expressed in bulk RNA-Seq datasets generated from the same cell line or organ. RESULTS: Validation of SCIFER using MCF7 cells determined technical parameters needed for optimal detection of L1 expression in single cells. We show that unsupervised analysis of L1 expression in single cells exponentially inflates both the levels of L1 expression and the number of expressed L1 loci. Application of SCIFER to analysis of scRNA-Seq datasets generated from mouse and human testes identified that mouse Round Spermatids and human Spermatogonia, Spermatocytes, and Round Spermatids express the highest levels of L1 mRNA. Our analysis also determined that similar to mice, human testes from unrelated individuals share as much as 80% of expressed L1 loci. Additionally, SCIFER determined that individual mouse cells co-express different L1 sub-families and different families of transposable elements, experimentally validating their co-existence in the same cell. CONCLUSIONS: SCIFER detects mRNA expression of individual L1 loci in single cells. It is compatible with scRNA-Seq datasets prepared using traditional sequencing methods. Validated using a human cancer cell line, SCIFER analysis of mouse and human testes identified key cell types supporting L1 expression in these species. This will further our understanding of differences and similarities in endogenous L1 mRNA expression patterns in mice and humans.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells.

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of 'transitional' L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Ebola virus delta peptide is an enterotoxin.

During the 2013-2016 West African (WA) Ebola virus (EBOV) outbreak, severe gastrointestinal symptoms were common in patients and associated with poor outcome. Delta peptide is a conserved product of post-translational processing of the abundant EBOV soluble glycoprotein (sGP). The murine ligated ileal loop model was used to demonstrate that delta peptide is a potent enterotoxin. Dramatic intestinal fluid accumulation follows injection of biologically relevant amounts of delta peptide into ileal loops, along with gross alteration of villous architecture and loss of goblet cells. Transcriptomic analyses show that delta peptide triggers damage response and cell survival pathways and downregulates expression of transporters and exchangers. Induction of diarrhea by delta peptide occurs via cellular damage and regulation of genes that encode proteins involved in fluid secretion. While distinct differences exist between the ileal loop murine model and EBOV infection in humans, these results suggest that delta peptide may contribute to EBOV-induced gastrointestinal pathology.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation.

Introduction: B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells. Methods: We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection). Results: Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing. Discussion: Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

Elevated ATGL in colon cancer cells and cancer stem cells promotes metabolic and tumorigenic reprogramming reinforced by obesity.

Obesity is a worldwide epidemic associated with increased risk and progression of colon cancer. Here, we aimed to determine the role of adipose triglyceride lipase (ATGL), responsible for intracellular lipid droplet (LD) utilization, in obesity-driven colonic tumorigenesis. In local colon cancer patients, significantly increased ATGL levels in tumor tissue, compared to controls, were augmented in obese individuals. Elevated ATGL levels in human colon cancer cells (CCC) relative to non-transformed were augmented by an obesity mediator, oleic acid (OA). In CCC and colonospheres, enriched in colon cancer stem cells (CCSC), inhibition of ATGL prevented LDs utilization and inhibited OA-stimulated growth through retinoblastoma-mediated cell cycle arrest. Further, transcriptomic analysis of CCC, with inhibited ATGL, revealed targeted pathways driving tumorigenesis, and high-fat-diet obesity facilitated tumorigenic pathways. Inhibition of ATGL in colonospheres revealed targeted pathways in human colonic tumor crypt base cells (enriched in CCSC) derived from colon cancer patients. In CCC and colonospheres, we validated selected transcripts targeted by ATGL inhibition, some with emerging roles in colonic tumorigeneses (ATG2B, PCK2, PGAM1, SPTLC2, IGFBP1, and ABCC3) and others with established roles (MYC and MUC2). These findings demonstrate obesity-promoted, ATGL-mediated colonic tumorigenesis and establish the therapeutic significance of ATGL in obesity-reinforced colon cancer progression.

Organ-, sex- and age-dependent patterns of endogenous L1 mRNA expression at a single locus resolution.

Expression of L1 mRNA, the first step in the L1 copy-and-paste amplification cycle, is a prerequisite for L1-associated genomic instability. We used a reported stringent bioinformatics method to parse L1 mRNA transcripts and measure the level of L1 mRNA expressed in mouse and rat organs at a locus-specific resolution. This analysis determined that mRNA expression of L1 loci in rodents exhibits striking organ specificity with less than 0.8% of loci shared between organs of the same organism. This organ specificity in L1 mRNA expression is preserved in male and female mice and across age groups. We discovered notable differences in L1 mRNA expression between sexes with only 5% of expressed L1 loci shared between male and female mice. Moreover, we report that the levels of total L1 mRNA expression and the number and spectrum of expressed L1 loci fluctuate with age as independent variables, demonstrating different patterns in different organs and sexes. Overall, our comparisons between organs and sexes and across ages ranging from 2 to 22 months establish previously unforeseen dynamic changes in L1 mRNA expression in vivo. These findings establish the beginning of an atlas of endogenous L1 mRNA expression across a broad range of biological variables that will guide future studies.

Transcriptome analysis reveals sexual disparities in gene expression in rat brain microvessels.

Sex is an important determinant of brain microvessels (MVs) function and susceptibility to cerebrovascular and neurological diseases, but underlying mechanisms are unclear. Using high throughput RNA sequencing analysis, we examined differentially expressed (DE) genes in brain MVs from young, male, and female rats. Bioinformatics analysis of the 23,786 identified genes indicates that 298 (1.2%) genes were DE using False Discovery Rate criteria (FDR; p < 0.05), of which 119 (40%) and 179 (60%) genes were abundantly expressed in male and female MVs, respectively. Nucleic acid binding, enzyme modulator, and transcription factor were the top three DE genes, which were more highly expressed in male than female MVs. Synthesis of glycosylphosphatidylinositol (GPI), biosynthesis of GPI-anchored proteins, steroid and cholesterol synthesis, were the top three significantly enriched canonical pathways in male MVs. In contrast, respiratory chain, ribosome, and 3 ́-UTR-mediated translational regulation were the top three enriched canonical pathways in female MVs. Different gene functions of MVs were validated by proteomic analysis and western blotting. Our novel findings reveal major sex disparities in gene expression and canonical pathways of MVs and these differences provide a foundation to study the underlying mechanisms and consequences of sex-dependent differences in cerebrovascular and other neurological diseases.

Transcriptional signatures of Zika virus infection in astrocytes.

Astrocytes are an early and important target of Zika virus (ZIKV) infection in the developing brain, but the impacts of infection on astrocyte function remain controversial. Given that nonhuman primate (NHP) models of ZIKV infection replicate aspects of neurologic disease seen in human infections, we cultured primary astrocytes from the brain tissue of infant rhesus macaques and then infected the cells with Asian or African lineage ZIKV to identify transcriptional patterns associated with infection in these cells. The African lineage virus appeared to have greater infectivity and promote stronger antiviral signaling, but infection by either strain ultimately produced typical virus response patterns. Both viruses induced hypoxic stress, but the Asian lineage strain additionally had an effect on metabolic and lipid biosynthesis pathways. Together, these findings describe an NHP astrocyte model that may be used to assess transcriptional signatures following ZIKV infection.

SON inhibits megakaryocytic differentiation via repressing RUNX1 and the megakaryocytic gene expression program in acute megakaryoblastic leukemia.

A high incidence of acute megakaryoblastic leukemia (AMKL) in Down syndrome patients implies that chromosome 21 genes have a pivotal role in AMKL development, but the functional contribution of individual genes remains elusive. Here, we report that SON, a chromosome 21-encoded DNA- and RNA-binding protein, inhibits megakaryocytic differentiation by suppressing RUNX1 and the megakaryocytic gene expression program. As megakaryocytic progenitors differentiate, SON expression is drastically reduced, with mature megakaryocytes having the lowest levels. In contrast, AMKL cells express an aberrantly high level of SON, and knockdown of SON induced the onset of megakaryocytic differentiation in AMKL cell lines. Genome-wide transcriptome analyses revealed that SON knockdown turns on the expression of pro-megakaryocytic genes while reducing erythroid gene expression. Mechanistically, SON represses RUNX1 expression by directly binding to the proximal promoter and two enhancer regions, the known +23 kb enhancer and the novel +139 kb enhancer, at the RUNX1 locus to suppress H3K4 methylation. In addition, SON represses the expression of the AP-1 complex subunits JUN, JUNB, and FOSB which are required for late megakaryocytic gene expression. Our findings define SON as a negative regulator of RUNX1 and megakaryocytic differentiation, implicating SON overexpression in impaired differentiation during AMKL development.

Somatic mutations in the DNA repairome in prostate cancers in African Americans and Caucasians.

Most hereditary tumors show aberrations in DNA repair genes or their regulators. In contrast, only a minority of sporadic tumors show alterations in these genes. As a result, genomic instability is currently considered an enhancer of tumorigenesis rather than an obligatory event in this process. However, tumor heterogeneity presents a significant technical challenge for most cancer genomics studies performed at less than 100× mean resolution depth. To address the importance of genomic instability in prostate carcinogenesis and tumor progression, we performed ultrahigh depth exome sequencing of 124 DNA damage repair/response (repairome) genes in 63 tumors and matched normal tissue samples in African Americans and Caucasians. The average sequence depth was 712-fold for DNA isolated from normal tissue and 368-fold for FFPE tumors. We identified 671 somatic mutations in tumors from African Americans and 762 somatic mutations in tumors in Caucasians. The most frequently mutated DNA repairome genes were EXO1, ATR, POLQ, NEIL3, ERCC6, BRCA2, BRCA1, XPC, JAG1, RPA1, POLE, ATM, and LIG1 in African American men, and POLQ, NEIL3, POLB, BRCA2, EXO1, ERCC6, ATR, RBBP8, BRCA1, ATM, JAG1, XPC, and POLE in Caucasians. We found that 89% of tumors had at least one mutation in nucleotide excision repair pathway genes in African Americans, whereas >40% of tumors had mutations in base excision repair pathway genes in Caucasians. We further identified a marginal increase in mutation rate in tumors in African Americans with increasing age. Tumors in Caucasians did not show a correlation with age, but a progressive increase in the mutation rate was observed at higher Gleason scores. Our data reveal significant differences in the molecular signatures in the DNA repairome in prostate cancer between African Americans and Caucasians. These data also have substantial implications regarding the well-known health disparities in prostate cancer, such as the higher mortality in African Americans than Caucasians.

Bacterial TLR4 and NOD2 signaling linked to reduced mitochondrial energy function in active inflammatory bowel disease.

Inflammatory bowel disease (IBD) has been linked to active signaling with bacterial components and reduced mitochondrial ATP production; however, synergism between both of these disease characteristics remains unclear. We aimed to determine in human IBD transcriptomes the link between a transcriptional signature unique to intestinal cells (ICs) with reduced mitochondrial ATP production (Mito-0) and bacteria triggered signaling using a bioinformatics approach. We generated an IC Mito-0 panel comprised of 199 differentially expressed (DE) transcripts mediated by reduced mitochondrial ATP function (DEGseq, log2 fold-change > |2|, p < .001). Transcripts from this panel were involved in diverse biological functions including regulation of mitochondrial energy (lower ATP), extracellular matrix, cell-cell contact, cytoskeleton, growth, metabolism, and inflammation. Next, unsupervised hierarchical clustering showed that the Mito-0 panel distinctly separated inflamed IBD from non-inflamed transcriptomes, which was also supported by principal component analysis (PCA) revealing distinct variation between sample types based on presence of the Mito-0 signature (PCA, p = 8.77e-09). Utilizing three independent IBD cohorts, we validated that 60 novel transcripts from the Mito-0 panel were significantly increased in inflamed tissue. Subsequently, KEGG generated bacterial TLR4 and NOD2 transcriptional signatures strongly associated with inflamed IBD transcriptomes and with the Mito-0 signature as determined by Spearman's analysis (coefficient of correlation, r = 0.92, p < .05). Herein, using a comprehensive analysis we demonstrated existence of an axis between bacteria triggered signaling and reduced mitochondrial energy function. Furthermore, we identified and validated novel transcripts within this axis as potential drivers and therapeutic targets for human IBD.