Epigenomic characterization of latent HIV infection identifies latency regulating transcription factors.

Transcriptional silencing of HIV in CD4 T cells generates a reservoir of latently infected cells that can reseed infection after interruption of therapy. As such, these cells represent the principal barrier to curing HIV infection, but little is known about their characteristics. To further our understanding of the molecular mechanisms of latency, we characterized a primary cell model of HIV latency in which infected cells adopt heterogeneous transcriptional fates. In this model, we observed that latency is a stable, heritable state that is transmitted through cell division. Using Assay of Transposon-Accessible Chromatin sequencing (ATACseq) we found that latently infected cells exhibit greatly reduced proviral accessibility, indicating the presence of chromatin-based structural barriers to viral gene expression. By quantifying the activity of host cell transcription factors, we observe elevated activity of Forkhead and Kruppel-like factor transcription factors (TFs), and reduced activity of AP-1, RUNX and GATA TFs in latently infected cells. Interestingly, latency reversing agents with different mechanisms of action caused distinct patterns of chromatin reopening across the provirus. We observe that binding sites for the chromatin insulator CTCF are highly enriched in the differentially open chromatin of infected CD4 T cells. Furthermore, depletion of CTCF inhibited HIV latency, identifying this factor as playing a key role in the initiation or enforcement of latency. These data indicate that HIV latency develops preferentially in cells with a distinct pattern of TF activity that promotes a closed proviral structure and inhibits viral gene expression. Furthermore, these findings identify CTCF as a novel regulator of HIV latency.

Pre-treatment differential correlation of gene expression and response to topical steroids in eosinophilic esophagitis.

Few predictors of response to topical corticosteroid (tCS) treatment have been identified in eosinophilic esophagitis (EoE). We aimed to determine whether baseline gene expression predicts histologic response to tCS treatment for EoE. We analyzed prospectively collected samples from incident EoE cases who were treated with tCS for 8 weeks in a development cohort (prospective study) or in an independent validation cohort (clinical trial). Whole transcriptome RNA expression was determined from a baseline (pre-treatment) RNA-later preserved esophageal biopsy. Baseline expression was compared between histologic responders (<15 eos/hpf) and non-responders (≥15 eos/hpf), and differential correlation was used to assess baseline gene expression by response status. In 87 EoE cases analyzed in the development set, there were no differentially expressed genes associated with treatment response (at false discovery rate = 0.1). However, differential correlation identified a module of 22 genes with statistically significantly high pairwise correlation in non-responders (mean correlation coefficient = 0.7) compared to low correlation in responders (coefficient = 0.3). When this 22-gene module was applied to the 89 EoE cases in the independent cohort, it was not validated to predict tCS response at the 15 eos/hpf threshold (mean correlation coefficient = 0.32 in responders and 0.25 in nonresponders). Exploration of other thresholds also did not validate any modules. Though we identified a 22 gene differential correlation module measured pre-treatment that was strongly associated with subsequent histologic response to tCS in EoE, this was not validated in an independent population. Alternative methods to predict steroid response should be explored.

microRNA-146a-5p association with the cardiometabolic disease risk factor TMAO.

Trimethylamine-N-oxide (TMAO), a microbial choline metabolism byproduct that is processed in the liver and excreted into circulation, is associated with increased atherosclerotic lesion formation and cardiovascular disease risk. Genetic regulators of TMAO levels are largely unknown. In the present study, we used 288 mice from a genetically heterogeneous mouse population [Diversity Outbred (DO)] to determine hepatic microRNA associations with TMAO in the context of an atherogenic diet. We also validated findings in two additional animal models of atherosclerosis: liver-specific insulin receptor knockout mice fed a chow diet (LIRKO) and African green monkeys fed high-fat/high-cholesterol diet. Small RNA-sequencing analysis in DO mice, LIRKO mice, and African green monkeys identified only one hepatic microRNA (miR-146a-5p) that is aberrantly expressed across all three models. Moreover, miR-146a-5p levels are associated with circulating TMAO after atherogenic diet in each of these models. We also performed high-resolution genetic mapping and identified a novel quantitative trait locus on Chromosome 12 for TMAO levels. This interval includes two genes, Numb and Dlst, which are inversely correlated with both miR-146a and TMAO and are predicted targets of miR-146a. Both of these genes have been validated as direct targets of miR-146a, though in other cellular contexts. This is the first report to our knowledge of a link between miR-146 and TMAO. Our findings suggest that miR-146-5p, as well as one or more genes at the Chromosome 12 QTL (possibly Numb or Dlst), is strongly linked to TMAO levels and likely involved in the control of atherosclerosis.

Systems genetics identifies a co-regulated module of liver microRNAs associated with plasma LDL cholesterol in murine diet-induced dyslipidemia.

Chronically altered levels of circulating lipids, termed dyslipidemia, is a significant risk factor for a number of metabolic and cardiovascular morbidities. MicroRNAs (miRNAs) have emerged as important regulators of lipid balance, have been implicated in dyslipidemia, and have been proposed as candidate therapeutic targets in lipid-related disorders including atherosclerosis. A major limitation of most murine studies of miRNAs in lipid metabolic disorders is that they have been performed in just one (or very few) inbred strains, such as C57BL/6. Moreover, although individual miRNAs have been associated with lipid phenotypes, it is well understood that miRNAs likely work together in functional modules. To address these limitations, we implemented a systems genetics strategy using the Diversity Outbred (DO) mouse population. Specifically, we performed gene and miRNA expression profiling in the livers from ~300 genetically distinct DO mice after 18 wk on either a high-fat/high-cholesterol diet or a high-protein diet. Large-scale correlative analysis of these data with a wide range of cardio-metabolic end points revealed a co-regulated module of miRNAs significantly associated with circulating low-density lipoprotein cholesterol (LDL-C) levels. The hubs of this module were identified as miR-199a, miR-181b, miR-27a, miR-21_-_1, and miR-24. In sum, we demonstrate that a high-fat/high-cholesterol diet robustly rewires the miRNA regulatory network, and we identify a small group of co-regulated miRNAs that may exert coordinated effects to control circulating LDL-C.

Aminobisphosphonates reactivate the latent reservoir in people living with HIV-1.

Antiretroviral therapy (ART) is not curative due to the existence of cellular reservoirs of latent HIV-1 that persist during therapy. Current research efforts to cure HIV-1 infection include "shock and kill" strategies to disrupt latency using small molecules or latency-reversing agents (LRAs) to induce expression of HIV-1 enabling cytotoxic immune cells to eliminate infected cells. The modest success of current LRAs urges the field to identify novel drugs with increased clinical efficacy. Aminobisphosphonates (N-BPs) that include pamidronate, zoledronate, or alendronate, are the first-line treatment of bone-related diseases including osteoporosis and bone malignancies. Here, we show the use of N-BPs as a novel class of LRA: we found in ex vivo assays using primary cells from ART-suppressed people living with HIV-1 that N-BPs induce HIV-1 from latency to levels that are comparable to the T cell activator phytohemagglutinin (PHA). RNA sequencing and mechanistic data suggested that reactivation may occur through activation of the activator protein 1 signaling pathway. Stored samples from a prior clinical trial aimed at analyzing the effect of alendronate on bone mineral density, provided further evidence of alendronate-mediated latency reversal and activation of immune effector cells. Decay of the reservoir measured by IPDA was however not detected. Our results demonstrate the novel use of N-BPs to reverse HIV-1 latency while inducing immune effector functions. This preliminary evidence merits further investigation in a controlled clinical setting possibly in combination with therapeutic vaccination.

Aminobisphosphonates reactivate the latent reservoir in people living with HIV-1.

Antiretroviral therapy (ART) is not curative due to the existence of cellular reservoirs of latent HIV-1 that persist during therapy. Current research efforts to cure HIV-1 infection include "shock and kill" strategies to disrupt latency using small molecules or latency-reversing agents (LRAs) to induce expression of HIV-1 enabling cytotoxic immune cells to eliminate infected cells. The modest success of current LRAs urges the field to identify novel drugs with increased clinical efficacy. Aminobisphosphonates (N-BPs) that include pamidronate, zoledronate, or alendronate, are the first-line treatment of bone-related diseases including osteoporosis and bone malignancies. Here, we show the use of N-BPs as a novel class of LRA: we found in ex vivo assays using primary cells from ART-suppressed people living with HIV-1 that N-BPs induce HIV-1 from latency to levels that are comparable to the T cell activator phytohemagglutinin (PHA). RNA sequencing and mechanistic data suggested that reactivation may occur through activation of the activator protein 1 signaling pathway. Stored samples from a prior clinical trial aimed at analyzing the effect of alendronate on bone mineral density, provided further evidence of alendronate-mediated latency reversal and activation of immune effector cells. Decay of the reservoir measured by IPDA was however not detected. Our results demonstrate the novel use of N-BPs to reverse HIV-1 latency while inducing immune effector functions. This preliminary evidence merits further investigation in a controlled clinical setting possibly in combination with therapeutic vaccination.

Stretching muscle cells induces transcriptional and splicing transitions and changes in SR proteins.

Alternative splicing is an RNA processing mechanism involved in skeletal muscle development and pathology. Muscular diseases exhibit splicing alterations and changes in mechanobiology leading us to investigate the interconnection between mechanical forces and RNA processing. We performed deep RNA-sequencing after stretching muscle cells. First, we uncovered transcriptional changes in genes encoding proteins involved in muscle function and transcription. Second, we observed that numerous mechanosensitive genes were part of the MAPK pathway which was activated in response to stretching. Third, we revealed that stretching skeletal muscle cells increased the proportion of alternatively spliced cassette exons and their inclusion. Fourth, we demonstrated that the serine and arginine-rich proteins exhibited stronger transcriptional changes than other RNA-binding proteins and that SRSF4 phosphorylation is mechanosensitive. Identifying SRSF4 as a mechanosensitive RNA-binding protein that might contribute to crosstalk between mechanotransduction, transcription, and splicing could potentially reveal novel insights into muscular diseases, particularly those with unknown etiologies.

Tumor Cell Extrinsic Synaptogyrin 3 Expression as a Diagnostic and Prognostic Biomarker in Head and Neck Cancer.

Over 70% of oropharyngeal head and neck squamous cell carcinoma (HNSC) cases in the United States are positive for human papillomavirus (HPV) yet biomarkers for stratifying oropharyngeal head and neck squamous cell carcinoma (HNSC) patient risk are limited. We used immunogenomics to identify differentially expressed genes in immune cells of HPV(+) and HPV(-) squamous carcinomas. Candidate genes were tested in clinical specimens using both quantitative RT-PCR and IHC and validated by IHC using the Carolina Head and Neck Cancer Study (CHANCE) tissue microarray of HNSC cases. We performed multiplex immunofluorescent staining to confirm expression within the immune cells of HPV(+) tumors, receiver operating characteristic (ROC) curve analyses, and assessed survival outcomes. The neuronal gene Synaptogyrin-3 (SYNGR3) is robustly expressed in immune cells of HPV(+) squamous cancers. Multiplex immunostaining and single cell RNA-seq analyses confirmed SYNGR3 expression in T cells, but also unexpectedly in B cells of HPV(+) tumors. ROC curve analyses revealed that combining SYNGR3 and p16 provides more sensitivity and specificity for HPV detection compared to p16 IHC alone. SYNGR3-high HNSC patients have significantly better prognosis with five-year OS and DSS rates of 60% and 71%, respectively. Moreover, combining p16 localization and SYNGR3 expression can further risk stratify HPV(+) patients such that high cytoplasmic, low nuclear p16 do significantly worse (Hazard Ratio, 8.6; P = 0.032) compared to patients with high cytoplasmic, high nuclear p16. SYNGR3 expression in T and B cells is associated with HPV status and enhanced survival outcomes of HNSC patients.

Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice.

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

Genetic Architecture Modulates Diet-Induced Hepatic mRNA and miRNA Expression Profiles in Diversity Outbred Mice.

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Lastly, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

FGFR4 regulates tumor subtype differentiation in luminal breast cancer and metastatic disease.

Mechanisms driving tumor progression from less aggressive subtypes to more aggressive states represent key targets for therapy. We identified a subset of luminal A primary breast tumors that give rise to HER2-enriched (HER2E) subtype metastases, but remain clinically HER2 negative (cHER2-). By testing the unique genetic and transcriptomic features of these cases, we developed the hypothesis that FGFR4 likely participates in this subtype switching. To evaluate this, we developed 2 FGFR4 genomic signatures using a patient-derived xenograft (PDX) model treated with an FGFR4 inhibitor, which inhibited PDX growth in vivo. Bulk tumor gene expression analysis and single-cell RNA sequencing demonstrated that the inhibition of FGFR4 signaling caused molecular switching. In the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohort, FGFR4-induced and FGFR4-repressed signatures each predicted overall survival. Additionally, the FGFR4-induced signature was an independent prognostic factor beyond subtype and stage. Supervised analysis of 77 primary tumors with paired metastases revealed that the FGFR4-induced signature was significantly higher in luminal/ER+ tumor metastases compared with their primaries. Finally, multivariate analysis demonstrated that the FGFR4-induced signature also predicted site-specific metastasis for lung, liver, and brain, but not for bone or lymph nodes. These data identify a link between FGFR4-regulated genes and metastasis, suggesting treatment options for FGFR4-positive patients, whose high expression is not caused by mutation or amplification.

Adropin: An endocrine link between the biological clock and cholesterol homeostasis.

OBJECTIVE: Identify determinants of plasma adropin concentrations, a secreted peptide translated from the Energy Homeostasis Associated (ENHO) gene linked to metabolic control and vascular function. METHODS: Associations between plasma adropin concentrations, demographics (sex, age, BMI) and circulating biomarkers of lipid and glucose metabolism were assessed in plasma obtained after an overnight fast in humans. The regulation of adropin expression was then assessed in silico, in cultured human cells, and in animal models. RESULTS: In humans, plasma adropin concentrations are inversely related to atherogenic LDL-cholesterol (LDL-C) levels in men (n = 349), but not in women (n = 401). Analysis of hepatic Enho expression in male mice suggests control by the biological clock. Expression is rhythmic, peaking during maximal food consumption in the dark correlating with transcriptional activation by RORα/γ. The nadir in the light phase coincides with the rest phase and repression by Rev-erb. Plasma adropin concentrations in nonhuman primates (rhesus monkeys) also exhibit peaks coinciding with feeding times (07:00 h, 15:00 h). The ROR inverse agonists SR1001 and the 7-oxygenated sterols 7-ß-hydroxysterol and 7-ketocholesterol, or the Rev-erb agonist SR9009, suppress ENHO expression in cultured human HepG2 cells. Consumption of high-cholesterol diets suppress expression of the adropin transcript in mouse liver. However, adropin over expression does not prevent hypercholesterolemia resulting from a high cholesterol diet and/or LDL receptor mutations. CONCLUSIONS: In humans, associations between plasma adropin concentrations and LDL-C suggest a link with hepatic lipid metabolism. Mouse studies suggest that the relationship between adropin and cholesterol metabolism is unidirectional, and predominantly involves suppression of adropin expression by cholesterol and 7-oxygenated sterols. Sensing of fatty acids, cholesterol and oxysterols by the RORα/γ ligand-binding domain suggests a plausible functional link between adropin expression and cellular lipid metabolism. Furthermore, the nuclear receptors RORα/γ and Rev-erb may couple adropin synthesis with circadian rhythms in carbohydrate and lipid metabolism.