Global Differential Gene Expression Research Articles

Role of HRTPT in kidney proximal epithelial cell regeneration: Integrative differential expression and pathway analyses using microarray and scRNA-seq.

Damage to proximal tubules due to exposure to toxicants can lead to conditions such as acute kidney injury (AKI), chronic kidney disease (CKD) and ultimately end‐stage renal failure (ESRF). Studies have shown that kidney proximal epithelial cells can regenerate particularly after acute injury. In the previous study, we utilized an immortalized in vitro model of human renal proximal tubule epithelial cells, RPTEC/TERT1, to isolate HRTPT cell line that co‐expresses stem cell markers CD133 and CD24, and HREC24T cell line that expresses only CD24. HRTPT cells showed most of the key characteristics of stem/progenitor cells; however, HREC24T cells did not show any of these characteristics. The goal of this study was to further characterize and understand the global gene expression differences, upregulated pathways and gene interaction using scRNA‐seq in HRTPT cells. Affymetrix microarray analysis identified common gene sets and pathways specific to HRTPT and HREC24T cells analysed using DAVID, Reactome and Ingenuity software. Gene sets of HRTPT cells, in comparison with publicly available data set for CD133+ infant kidney, urine‐derived renal progenitor cells and human kidney‐derived epithelial proximal tubule cells showed substantial similarity in organization and interactions of the apical membrane. Single‐cell analysis of HRTPT cells identified unique gene clusters associated with CD133 and the 92 common gene sets from three data sets. In conclusion, the gene expression analysis identified a unique gene set for HRTPT cells and narrowed the co‐expressed gene set compared with other human renal–derived cell lines expressing CD133, which may provide deeper understanding in their role as progenitor/stem cells that participate in renal repair.

Abstract 27: Unravelling Cell-specific Interactions At The Preeclamptic Maternal-foetal Interface From Early To Late Pregnancy

Preeclamptic syndrome arrises in the fetal part of the placenta (villi). In this study, we analyse placental development by single nuclei RNA-sequencing in early and term pregnancy and draw conclusions about pathological processes in preeclampsia (PE) that originate early in gestation. We profiled the transcriptome of 101,067 nuclei obtained from a total of 12 pregnancies, spanning early, term and PE doners. Using unsupervised computational approaches, we identified 12 and 16 different cell types and states in decidua and villi, respectively. Our comprehensively identified catalogue of cell types and states aligns well with the previous single cell studies. We identified different subpopulations of syncytiotrophoblast and GATA3+/GREM2+ trophoblast stem cells (TSC) in villi. Through gestation, gene expression in cell populations from the matrisome or vascular environments show dynamic expression reflecting vascular development associated with spiral artery remodelling and concordant decidual stroma reorganisation. Global differential gene expression analysis shows that trophoblast cell types are most dysregulated in PE. Cell-cell communication analysis revealed important dysregulation between villi and decidual cell types. The secretory signalling characteristic of this trophoblastic disease may be used for early biomarker screening. Overall, this study paves the way to a deeper understanding of the early pathophysiology of PE. Figure 1: Villi (v) and decidua (d) cell clusters from early, late control and preeclampsia (PE) villi and decidua visualised as a UMAP. Datasets were integrated separately for each tissue and merged for embedding.

Abstract 1590: An experimental tumor model resistant to anti-PD-1 is sensitive to local mRNA immunotherapy

Abstract Immune checkpoint blockade elicits durable anti-cancer responses in the clinic, however a large proportion of patients do not benefit from treatment. In an effort to better understand acquired resistance to checkpoint blockade, we generated a mouse tumor model exhibiting in vivo resistance to anti-PD-1 antibody treatment. MC38 tumors acquired resistance to PD-1 blockade following serial in vivo passaging. Lack of sensitivity to PD-1 blockade was not attributed to dysregulation of PD-L1 or β2M expression, as both were expressed at similar levels in parental and resistant cells. Similarly, IFNγ signaling and antigen processing and presentation pathways were functional in both parental and resistant cell lines. RNA-sequencing revealed substantial differences in global gene expression, with tumors resistant to anti-PD-1 displaying a marked reduction in expression of immune-related genes relative to parental MC38 tumors. Indeed, resistant tumors exhibited reduced immune infiltration across multiple cell types, including T and NK cells. Pathway analysis revealed activation of TGFβ and Notch signaling in anti-PD-1 resistant tumors, and activation of these pathways was associated with poorer survival in human cancer patients. While pharmacological inhibition of TGFβ and Notch in combination with PD-1 blockade decelerated tumor growth, local administration of a mixture of mRNAs encoding immune-stimulatory cytokines (IL12sc, GM-CSF, IL15s and IFNα) potently induced regression of resistant tumors and resulted in complete tumor remission. Overall, this study describes a novel anti-PD-1 resistant mouse tumor model and underscores the role of two well-defined signaling pathways in response to immune checkpoint blockade. Furthermore, our data highlights the potential of intratumoral mRNA therapy in overcoming acquired resistance to PD-1 blockade. Citation Format: Marie Bernardo, Tatiana Tolstykh, Yu-an Zhang, Dinesh S. Bangari, Hui Cao, Kerstin A. Heyl, Joon Sang Lee, Natalia Malkova, Jack Pollard, Hui Qu, Fangxian Sun, Dmitri Wiederschain, Timothy R. Wagenaar. An experimental tumor model resistant to anti-PD-1 is sensitive to local mRNA immunotherapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 1590.

Contrasting Stage Dependent Requirements for the Histone Methyltransferase EZH2 during the Evolution of Acute Myeloid Leukemia

Epigenetic regulators function in a context-dependent manner and are commonly mutated in tumors. Activating mutations and overexpression of the histone-methyltransferase, enhancer of zeste, homolog 2 (EZH2), occur in lymphoma and other malignancies, whilst loss-of-function mutations are predicted in myeloid malignancies. We study this apparent contradiction, examining the importance of cellular context and disease stage for Ezh2 loss, genetically deleting Ezh2 during the induction or maintenance of a single malignancy, Acute Myeloid Leukemia (AML).We initially focused on the role of Ezh2 in AML maintenance, using conditional Ezh2fl/fl mice. Following in vitro retroviral immortalization of haematopoietic stem and progenitor cells (HSPCs) from these with AML-fusion oncogenes: MLL-AF9 or AML1-ETO9a (reflecting two separate models of AML clinically associated with good and poor prognosis, markedly differing in mechanism of leukemic transformation and representing AML heterogeneity), these were then transduced with p-babe-Cre-puro vector (or empty vector) to delete Ezh2 . Strikingly, Ezh2 loss completely abolished cellular growth and colony formation, whilst empty vector transduced cells maintained their clonogenic output. To interrogate Ezh2 as a target in AML maintenance in vivo, we generated secondary MLL-AF9 or AML1-ETO9a leukemias through transplantation of bulk leukemia derived from Ezh2fl/fl ; Mx1-Cre primary leukemic mice. Following engraftment, deletion of Ezh2 was achieved by pIpC administration. Loss of Ezh2 in vivo markedly impaired both secondary MLL-AF9 and AML1-ETO9a leukemia progression. These data suggest that Ezh2 is an oncogene across multiple disparate subtypes of established AML and propose it as a possible therapeutic target in this context .Addressing this hypothesis, we examined the role of pharmacological inhibition of EZH2 in AML, utilizing GSK343, an S-Adenosyl Methionine competitive EZH2 inhibitor. Through colony formation and functional assays, we demonstrated in vitro sensitivities of MLL-AF9 transformed murine BM HSPCs, primary murine Ezh2+/+ MLL-AF9 tumours , the AML1-ETO fusion driven human cell line Kasumi-1 and primary leukemic cells isolated from patients with various AML genotypes to GSK343 at similar IC50 (~10µM) across all cells. In contrast, minimal effects of GSK343 were seen on normal primary CD34+ HSCs clonogenic function. Together these data provide evidence of therapeutic potential, validating EZH2 as a target in AML.To characterise the role of Ezh2 during AML induction in vivo, we utilized Ezh2fl/fl; Mx1-Cre mice to generate Ezh2-/- mice and compared them to Ezh2+/+ controls. HSPCs from these were retrovirally transduced with MLL-AF9 or AML1-ETO9a and transplanted into lethally irradiated wild-type recipients, and leukemias allowed to develop. Remarkably in the context of disease induction, and in stark contrast to the maintenance experiments, Ezh2 functioned as a tumor suppressor gene, with its loss significantly accelerating leukemia generation across both subtypes of AML (compared to control).To interrogate the molecular mechanisms of these findings we performed integrated genomic analysis with differential global gene expression analysis via RNA-sequencing (in Ezh2+/+ vs Ezh2-/- non-transformed HSPCs, Ezh2+/+ vs Ezh2-/- AML1-ETO9a and Ezh2+/+ vs Ezh2-/- MLL-AF9 leukemias) and ChIP-sequencing of the histone modifications H3K27me3 and H3K4me3. These revealed that during disease induction, Ezh2 represses a subset of bivalent promoters that resolve towards gene activation upon Ezh2 loss, inducing a number of genes including the transcription factor Plag1 (with known oncogenic function in CBFB-MYH11 murine leukemogenesis). Of note, overexpression of Plag1 phenocopied Ezh2 loss i n vivo to significantly accelerate MLL-AF9 leukemia induction.In addition, we demonstrated virtual mutual exclusivity of genes de-repressed in established MLL-AF9 leukemias following pharmacological Ezh2 inhibition compared to those de-repressed in the induction of Ezh2-/- MLL-AF9 leukemias. This de-repression of entirely separate gene programmes explains the starkly contrasting phenotypic outcomes of targeting Ezh2 in AML induction and maintenance, and further provides reassurance that EZH2 inhibition in established AML is a safe strategy. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Comparative Transcriptomic Analysis of Spermatozoa From High- and Low-Fertile Crossbred Bulls: Implications for Fertility Prediction.

Crossbred bulls produced by crossing Bos taurus and Bos indicus suffer with high incidence of infertility/subfertility problems; however, the etiology remains poorly understood. The uncertain predictability and the inability of semen evaluation techniques to maintain constant correlation with fertility demand for alternate methods for bull fertility prediction. Therefore, in this study, the global differential gene expression between high- and low-fertile crossbred bull sperm was assessed using a high-throughput RNA sequencing technique with the aim to identify transcripts associated with crossbred bull fertility. Crossbred bull sperm contained transcripts for 13,563 genes, in which 2,093 were unique to high-fertile and 5,454 were unique to low-fertile bulls. After normalization of data, a total of 776 transcripts were detected, in which 84 and 168 transcripts were unique to high-fertile and low-fertile bulls, respectively. A total of 176 transcripts were upregulated (fold change > 1) and 209 were downregulated (<1) in low-fertile bulls. Gene ontology analysis identified that the sperm transcripts involved in the oxidative phosphorylation pathway and biological process such as multicellular organism development, spermatogenesis, and in utero embryonic development were downregulated in low-fertile crossbred bull sperm. Sperm transcripts upregulated and unique to low-fertile bulls were majorly involved in translation (biological process) and ribosomal pathway. With the use of RT-qPCR, selected sperm transcripts (n = 12) were validated in crossbred bulls (n = 12) with different fertility ratings and found that the transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes was significantly (p < 0.05) lower in low-fertile bulls than high-fertile bulls and was positively (p < 0.05) correlated with conception rate. It is inferred that impaired oxidative phosphorylation could be the predominant reason for low fertility in crossbred bulls and that transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes could serve as potential biomarkers for fertility in crossbred bulls.

Temporal Transcriptome Analysis Reveals Dynamic Gene Expression Patterns Driving β-Cell Maturation.

Newly differentiated pancreatic β cells lack proper insulin secretion profiles of mature functional β cells. The global gene expression differences between paired immature and mature β cells have been studied, but the dynamics of transcriptional events, correlating with temporal development of glucose-stimulated insulin secretion (GSIS), remain to be fully defined. This aspect is important to identify which genes and pathways are necessary for β-cell development or for maturation, as defective insulin secretion is linked with diseases such as diabetes. In this study, we assayed through RNA sequencing the global gene expression across six β-cell developmental stages in mice, spanning from β-cell progenitor to mature β cells. A computational pipeline then selected genes differentially expressed with respect to progenitors and clustered them into groups with distinct temporal patterns associated with biological functions and pathways. These patterns were finally correlated with experimental GSIS, calcium influx, and insulin granule formation data. Gene expression temporal profiling revealed the timing of important biological processes across β-cell maturation, such as the deregulation of β-cell developmental pathways and the activation of molecular machineries for vesicle biosynthesis and transport, signal transduction of transmembrane receptors, and glucose-induced Ca2+ influx, which were established over a week before β-cell maturation completes. In particular, β cells developed robust insulin secretion at high glucose several days after birth, coincident with the establishment of glucose-induced calcium influx. Yet the neonatal β cells displayed high basal insulin secretion, which decreased to the low levels found in mature β cells only a week later. Different genes associated with calcium-mediated processes, whose alterations are linked with insulin resistance and deregulation of glucose homeostasis, showed increased expression across β-cell stages, in accordance with the temporal acquisition of proper GSIS. Our temporal gene expression pattern analysis provided a comprehensive database of the underlying molecular components and biological mechanisms driving β-cell maturation at different temporal stages, which are fundamental for better control of the in vitro production of functional β cells from human embryonic stem/induced pluripotent cell for transplantation-based type 1 diabetes therapy.

Abstract PO048: Loss of intra-tumoral RIG-I immune signaling is a potential microbiome-mediated mechanism underlying poor anti-tumor immunity and immunotherapy resistance in mucosal melanoma

Abstract Cutaneous melanoma (CM) has one of the highest response rates across cancers to immune checkpoint blockade (ICB) therapies, but mucosal melanoma (MM) responds poorly to ICB. MM patients have very few effective treatments options, making it critical to understand and improve anti-tumor immunity in these patients. MM is a rare subtype of melanoma, accounting for approximately 1% of all malignant melanoma diagnoses. MM occurs on mucosal surfaces such as nasal, vaginal, and anorectal, which provide a highly unique tumor microenvironment compared to CM arising in the skin. We compared the tumor immune microenvironment between CM and MM by multiplex immunofluorescence and discovered MM has decreased frequency of infiltrating immune cells. Using RNA sequencing, we investigated intra-tumoral differences in global gene expression between CM and MM that might explain the lack of immune cell infiltration and ICB response in MM. Compared to CM, MM had decreased expression of numerous innate immune genes in the RIG-I pathogen-sensing pathway. The RIG-I pathway produces an inflammatory response to destroy foreign pathogens, but is also important for ICB sensitivity and anti-tumor immunity in CM. Studies in CM demonstrate the RIG-I pathway enhances antigen presentation and promotes a favorable immune cell profile in the tumor microenvironment. Interestingly, the tumor microbiome, which is vastly different between CM and MM, has been shown to regulate the RIG-I pathway. We measured total fungal and bacterial load in MM and CM primary tumors and found increased levels of both bacteria and fungi in MM compared to CM. This raises the possibility that microbiome-mediated suppression of the RIG-I pathway underlies poor anti-tumor immunity and ICB response in MM, and that re-activation of this pathway may be a novel therapeutic strategy for overcoming ICB resistance in MM patients. Direct RIG-I agonists are still in clinical development, but FDA-approved hypomethylating agents, particularly 5’ aza-deoxycytidine (decitabine), can activate the RIG-I pathway and improve anti-tumor immunity by increasing the expression of RIG-I pathway activating genes, endogenous RNA retroviruses, tumor antigens, and natural killer (NK) cell ligands. We treated MM cell lines with decitabine and observed strong induction of RIG-I pathway genes, NK ligands, and tumor antigens. In conclusion, we have identified the loss of intra-tumoral RIG-I signaling as a potential microbiome-mediated mechanism underlying the poor tumor immune microenvironment and ICB response in MM. Combining decitabine, which is FDA-approved and now available orally, with ICB represents a novel treatment strategy for these difficult to treat MM patients that can be rapidly translated into the clinic. Citation Format: Morgan MacBeth, Richard Tobin, Robert Van Gulick, Martin D. McCarter, William A. Robinson, Kasey L. Couts. Loss of intra-tumoral RIG-I immune signaling is a potential microbiome-mediated mechanism underlying poor anti-tumor immunity and immunotherapy resistance in mucosal melanoma [abstract]. In: Proceedings of the AACR Virtual Special Conference on the Evolving Tumor Microenvironment in Cancer Progression: Mechanisms and Emerging Therapeutic Opportunities; in association with the Tumor Microenvironment (TME) Working Group; 2021 Jan 11-12. Philadelphia (PA): AACR; Cancer Res 2021;81(5 Suppl):Abstract nr PO048.

Transcriptomics of Long-Term Meditation Practice: Evidence for Prevention or Reversal of Stress Effects Harmful to Health.

Background and Objectives: Stress can overload adaptive mechanisms, leading to epigenetic effects harmful to health. Research on the reversal of these effects is in its infancy. Early results suggest some meditation techniques have health benefits that grow with repeated practice. This study focused on possible transcriptomic effects of 38 years of twice-daily Transcendental Meditation® (TM®) practice. Materials and Methods: First, using Illumina® BeadChip microarray technology, differences in global gene expression in peripheral blood mononuclear cells (PBMCs) were sought between healthy practitioners and tightly matched controls (n = 12, age 65). Second, these microarray results were verified on a subset of genes using quantitative polymerase chain reaction (qPCR) and were validated using qPCR in larger TM and control groups (n = 45, age 63). Bioinformatics investigation employed Ingenuity® Pathway Analysis (IPA®), DAVID, Genomatix, and R packages. Results: The 200 genes and loci found to meet strict criteria for differential expression in the microarray experiment showed contrasting patterns of expression that distinguished the two groups. Differential expression relating to immune function and energy efficiency were most apparent. In the TM group, relative to the control, all 49 genes associated with inflammation were downregulated, while genes associated with antiviral and antibody components of the defense response were upregulated. The largest expression differences were shown by six genes related to erythrocyte function that appeared to reflect a condition of lower energy efficiency in the control group. Results supporting these gene expression differences were obtained with qPCR-measured expression both in the well-matched microarray groups and in the larger, less well-matched groups. Conclusions: These findings are consistent with predictions based on results from earlier randomized trials of meditation and may provide evidence for stress-related molecular mechanisms underlying reductions in anxiety, post-traumatic stress disorder (PTSD), cardiovascular disease (CVD), and other chronic disorders and diseases.

An experimental model of anti-PD-1 resistance exhibits activation of TGFß and Notch pathways and is sensitive to local mRNA immunotherapy

ABSTRACT Immune checkpoint blockade elicits durable anti-cancer responses in the clinic, however a large proportion of patients do not benefit from treatment. Several mechanisms of innate and acquired resistance to checkpoint blockade have been defined and include mutations of MHC I and IFNγ signaling pathways. However, such mutations occur in a low frequency of patients and additional mechanisms have yet to be elucidated. In an effort to better understand acquired resistance to checkpoint blockade, we generated a mouse tumor model exhibiting in vivo resistance to anti-PD-1 antibody treatment. MC38 tumors acquired resistance to PD-1 blockade following serial in vivo passaging. Lack of sensitivity to PD-1 blockade was not attributed to dysregulation of PD-L1 or β2M expression, as both were expressed at similar levels in parental and resistant cells. Similarly, IFNγ signaling and antigen processing and presentation pathways were functional in both parental and resistant cell lines. Unbiased gene expression analysis was used to further characterize potential resistance mechanisms. RNA-sequencing revealed substantial differences in global gene expression, with tumors resistant to anti-PD-1 displaying a marked reduction in expression of immune-related genes relative to parental MC38 tumors. Indeed, resistant tumors exhibited reduced immune infiltration across multiple cell types, including T and NK cells. Pathway analysis revealed activation of TGFβ and Notch signaling in anti-PD-1 resistant tumors, and activation of these pathways was associated with poorer survival in human cancer patients. While pharmacological inhibition of TGFβ and Notch in combination with PD-1 blockade decelerated tumor growth, a local mRNA-based immunotherapy potently induced regression of resistant tumors, resulting in complete tumor remission, and resensitized tumors to treatment with anti-PD-1. Overall, this study describes a novel anti-PD-1 resistant mouse tumor model and underscores the role of two well-defined signaling pathways in response to immune checkpoint blockade. Furthermore, our data highlights the potential of intratumoral mRNA therapy in overcoming acquired resistance to PD-1 blockade.

Increased expression of interferon regulated and antiviral response genes in CD31+/CD102+ lung microvascular endothelial cells from systemic sclerosis patients with end-stage interstitial lung disease.

Systemic sclerosis (SSc) is characterised by severe fibroproliferative vasculopathy, fibrosis in skin and multiple internal organs, and humoral, cellular and innate immunity abnormalities. Vascular alterations are the earliest and most severe SSc manifestations, however, the mechanisms responsible have remained elusive. To investigate the molecular abnormalities involved in SSc-vasculopathy we examined global gene expression differences between highly purified lung microvascular endothelial cells (MVECs) from patients with SSc-interstitial lung disease (SSc-ILD) and normal lung MVECs. MVECs were isolated from fresh transplanted lungs from patients with SSc-ILD. Sequential CD31 and CD102 immunopurification was performed to obtain highly purified CD31+/CD102+ lung MVECs. Global gene expression analysis was successfully performed in CD31+/CD102+ MVEC from two SSc-ILD patients and from two normal lungs. RT-PCR, Western blots, and indirect immunofluorescence validated the gene expression results. Numerous interferon-regulated genes (IRGs) including IFI44, IFI44L, IFI6, IFIH1, IFIT1, ISG-15, BST-2/Tetherin, and RSAD2/Viperin, genes encoding innate immunity antiviral responses (OAS1, OAS2, OAS3, OASL) and antiviral MX1 and MX2 proteins, and mesenchymal cell-specific genes were significantly overexpressed in CD31+/CD102+ SSc-ILD lung MVECs. Highly purified CD31+/CD102+ MVECs from lungs from SSc patients with end stage SSc-ILD displayed remarkable overexpression of numerous IRGs and of genes encoding antiviral innate immune response and antiviral proteins. These observations suggest that interferon-induced and antiviral response proteins may participate in the pathogenesis of SSc vasculopathy and SSc-ILD. The CD31+/CD102+ lung MVECs from SSc-ILD also showed elevated expression of mesenchymal cell-specific genes confirming the presence of endothelial to mesenchymal transition in SSc-ILD.

Comparative Transcriptomics of Shiga Toxin-Producing and Commensal Escherichia coli and Cytokine Responses in Colonic Epithelial Cell Culture Infections.

Ingestion of Shiga toxin-producing Escherichia coli (STEC) can result in a range of illness severity from asymptomatic to hemorrhagic colitis and death; thus risk assessment of STEC strains for human pathogenicity is important in the area of food safety. Illness severity depends in part on the combination of virulence genes carried in the genome, which can vary between strains even of identical serotype. To better understand how core genes are regulated differently among strains and to identify possible novel STEC virulence gene candidates that could be added to the risk assessment repertoire, we used comparative transcriptomics to investigate global gene expression differences between two STEC strains associated with severe illness and a commensal E. coli strain during in vitro intestinal epithelial cell (IEC) infections. Additionally, we compared a wide array of concomitant cytokine levels produced by the IECs. The cytokine expression levels were examined for a pattern representing STEC pathogenicity; however, while one STEC strain appeared to elicit a proinflammatory response, infection by the other strain produced a pattern comparable to the commensal E. coli. This result may be explained by the significant differences in gene content and expression observed between the STEC strains. RNA-Seq analysis revealed considerable disparity in expression of genes in the arginine and tryptophan biosynthesis/import pathways between the STEC strains and the commensal E. coli strain, highlighting the important role some amino acids play in STEC colonization and survival. Contrasting differential expression patterns were observed for genes involved in respiration among the three strains suggesting that metabolic diversity is a strategy utilized to compete with resident microflora for successful colonization. Similar temporal expression results for known and putative virulence genes were observed in the STEC strains, revealing strategies used for survival prior to and after initial adherence to IECs. Additionally, three genes encoding hypothetical proteins located in mobile genetic elements were, after interrogation of a large set of E. coli genomes, determined to likely represent novel STEC virulence factors.

STAT3 and GR Cooperate to Drive Gene Expression and Growth of Basal-Like Triple-Negative Breast Cancer.

Breast cancers are divided into subtypes with different prognoses and treatment responses based on global differences in gene expression. Luminal breast cancer gene expression and proliferation are driven by estrogen receptor alpha, and targeting this transcription factor is the most effective therapy for this subtype. By contrast, it remains unclear which transcription factors drive the gene expression signature that defines basal-like triple-negative breast cancer, and there are no targeted therapies approved to treat this aggressive subtype. In this study, we utilized integrated genomic analysis of DNA methylation, chromatin accessibility, transcription factor binding, and gene expression in large collections of breast cancer cell lines and patient tumors to identify transcription factors responsible for the basal-like gene expression program. Glucocorticoid receptor (GR) and STAT3 bind to the same genomic regulatory regions, which were specifically open and unmethylated in basal-like breast cancer. These transcription factors cooperated to regulate expression of hundreds of genes in the basal-like gene expression signature, which were associated with poor prognosis. Combination treatment with small-molecule inhibitors of both transcription factors resulted in synergistic decreases in cell growth in cell lines and patient-derived organoid models. This study demonstrates that GR and STAT3 cooperate to regulate the basal-like breast cancer gene expression program and provides the basis for improved therapy for basal-like triple-negative breast cancer through rational combination of STAT3 and GR inhibitors. SIGNIFICANCE: This study demonstrates that GR and STAT3 cooperate to activate the canonical gene expression signature of basal-like triple-negative breast cancer and that combination treatment with STAT3 and GR inhibitors could provide synergistic therapeutic efficacy.

Multiple Known Mechanisms and a Possible Role of an Enhanced Immune System in Bt-Resistance in a Field Population of the Bollworm, Helicoverpa zea: Differences in Gene Expression with RNAseq.

Several different agricultural insect pests have developed field resistance to Bt (Bacillus thuringiensis) proteins (ex. Cry1Ac, Cry1F, etc.) expressed in crops, including corn and cotton. In the bollworm, Helicoverpa zea, resistance levels are increasing; recent reports in 2019 show up to 1000-fold levels of resistance to Cry1Ac, a major insecticidal protein in Bt-crops. A common method to analyze global differences in gene expression is RNA-seq. This technique was used to measure differences in global gene expression between a Bt-susceptible and Bt-resistant strain of the bollworm, where the differences in susceptibility to Cry1Ac insecticidal proteins were 100-fold. We found expected gene expression differences based on our current understanding of the Bt mode of action, including increased expression of proteases (trypsins and serine proteases) and reduced expression of Bt-interacting receptors (aminopeptidases and cadherins) in resistant bollworms. We also found additional expression differences for transcripts that were not previously investigated, i.e., transcripts from three immune pathways-Jak/STAT, Toll, and IMD. Immune pathway receptors (ex. PGRPs) and the IMD pathway demonstrated the highest differences in expression. Our analysis suggested that multiple mechanisms are involved in the development of Bt-resistance, including potentially unrecognized pathways.

Thermotolerant coral symbionts modulate heat stress-responsive genes in their hosts.

Some corals may become more resistant to bleaching by shuffling their Symbiodiniaceae communities toward thermally tolerant species, and manipulations to boost the abundance of these symbionts in corals may increase resilience in warming oceans. However, the thermotolerant symbiont Durusdinium trenchii may reduce growth and fecundity in Caribbean corals, and these tradeoffs need to be better understood as this symbiont spreads through the region. We sought to understand how D.trenchii modulates coral gene expression by manipulating symbiont communities in Montastraea cavernosa to produce replicate ramets containing D.trenchii together with paired ramets of these same genets (n=3) containing Cladocopium C3 symbionts. We then examined differences in global gene expression between corals hosting Durusdinium and Cladocopium under control temperatures, and in response to short-term heat stress. We identified numerous transcriptional differences associated with symbiont identity, which explained 2%-14% of the transcriptional variance. Corals with D.trenchii upregulated genes related to translation, ribosomal structure and biogenesis, and downregulated genes related to extracellular structures, and carbohydrate and lipid transport and metabolism, relative to corals with Cladocopium. Unexpectedly, these changes were similar to those observed in Cladocopium-dominated corals in response to heat stress, suggesting that thermotolerant D.trenchii may cause corals to increase expression of heat stress-responsive genes, explaining both the increased heat tolerance and the associated energetic tradeoffs in corals containing D.trenchii. These findings provide insight into the ecological changes occurring on contemporary coral reefs in response to climate change, and the diverse ways in which different symbionts modulate emergent phenotypes of their hosts.

Transcriptome analysis of testis reveals the effects of developmental exposure to bisphenol a or 17α-ethinylestradiol in medaka (Oryzias latipes)

Endocrine disrupting chemicals (EDCs) can induce abnormalities in organisms via alteration of molecular pathways and subsequent disruption of endocrine functions. Bisphenol A (BPA) and 17α-ethinylestradiol (EE2) are ubiquitous EDCs in the environment. Many aquatic organisms, including fish, are often exposed to varying concentrations of BPA and EE2 throughout their lifespan. Both BPA and EE2 can activate estrogenic signaling pathways and cause adverse effects on reproduction via alteration of pathways associated with steroidogenesis. However, transcriptional pathways that are affected by chronic exposure to these two ubiquitous environmental estrogens during embryonic, larval, and juvenile stages are not clearly understood. In the present study, we examined transcriptional alterations in the testis of medaka fish (Oryzias latipes) chronically exposed to a low concentration of BPA or EE2. Medaka were exposed to BPA (10 μg/L) or EE2 (0.01 μg/L) from 8 h post-fertilization (as embryos) to adulthood 50 days post fertilization (dpf), and transcriptional alterations in the testis were examined by RNA sequencing (RNA-seq). Transcriptomic profiling revealed 651 differentially expressed genes (DEGs) between BPA-exposed and control testes, while 1475 DEGs were found between EE2-exposed and control testes. Gene ontology (GO) analysis showed a significant enrichment of “intracellular receptor signaling pathway”, “response to steroid hormone” and “hormone-mediated signaling pathway” in the BPA-induced DEGs, and of “cilium organization”, “microtubule-based process” and “organelle assembly” in the EE2-induced DEGs. Pathway analysis showed significant enrichment of “integrin signaling pathway” in both treatment groups, and of “cadherin signaling pathway”, “Alzheimer disease-presenilin pathway” in EE2-induced DEGs. Single nucleotide polymorphism (SNP) and insertion-deletion (Indel) analysis found no significant differences in mutation rates with either BPA or EE2 treatments. Taken together, global gene expression differences in testes of medaka during early stages of gametogenesis were responsive to chronic BPA and EE2 exposure.

Differential expression of lung adenocarcinoma transcriptome with signature of tobacco exposure

Smoking accounts for almost 80–90% of lung cancer cases, which is also the most frequent cause of cancer-related deaths in humans. With over 60 carcinogens in tobacco smoke, cells dividing at the time of carcinogen exposure are at particular risk of neoplasia. The present study aimed to investigate global gene expression differences in lung adenocarcinoma (LUAD) tumour samples of current smokers and non-smokers, in an attempt to elucidate biological mechanisms underlying divergent smoking effects. Current and non-smoker tumour samples were analysed using bioinformatics tools, examining differences in molecular drivers of cancer initiation and progression, as well as evaluating the effect of smoking and sex on epithelial mesenchymal transition (EMT). As a result, we identified 1150 differentially expressed genes showing visible differences in the expression profiles between the smoking subgroups. The genes were primarily involved in cell cycle, DNA replication, DNA repair, VEGF, GnRH, ErbB and T cell receptor signalling pathways. Our results show that smoking clearly affected E2F transcriptional activity and DNA repair pathways including mismatch repair, base excision repair and homologous recombination. We observed that sex could modify the effects of PLA2G2A and PRG4 in LUAD tumour samples, whereas sex and smoking status might possibly have a biological effect on the EMT-related genes: HEY2, OLFM1, SFRP1 and STRAP. We also identified potential epigenetic changes smoking solely might have on EMT-related genes, which may serve as potential diagnostic and prognostic biomarkers for LUAD patients.

IFN-I Independent Antiviral Immune Response to Vesicular Stomatitis Virus Challenge in Mouse Brain.

Type I interferon (IFN-I) plays a pivotal role during viral infection response in the central nervous system (CNS). The IFN-I can orchestrate and regulate most of the innate immune gene expression and myeloid cell dynamics following a noncytopathic virus infection. However, the role of IFN-I in the CNS against viral encephalitis is not entirely clear. Here we have implemented the combination of global differential gene expression profiling followed by bioinformatics analysis to decipher the CNS immune response in the presence and absence of the IFN-I signaling. We observed that vesicular stomatitis virus (VSV) infection induced 281 gene changes in wild-type (WT) mice primarily associated with IFN-I signaling. This was accompanied by an increase in antiviral response through leukocyte vascular patrolling and leukocyte influx along with the expression of potent antiviral factors. Surprisingly, in the absence of the IFN-I signaling (IFNAR−/− mice), a significantly higher (1357) number of genes showed differential expression compared to the WT mice. Critical candidates such as IFN-γ, CCL5, CXCL10, and IRF1, which are responsible for the recruitment of the patrolling leukocytes, are also upregulated in the absence of IFN-I signaling. The computational network analysis suggests the presence of the IFN-I independent pathway that compensates for the lack of IFN-I signaling in the brain. The analysis shows that TNF-α is connected maximally to the networked candidates, thus emerging as a key regulator of gene expression and recruitment of myeloid cells to mount antiviral action. This pathway could potentiate IFN-γ release; thereby, synergistically activating IRF1-dependent ISG expression and antiviral response.

RNA-seq reveals downregulated osteochondral genes potentially related to tibia bacterial chondronecrosis with osteomyelitis in broilers

BackgroundBacterial chondronecrosis with osteomyelitis (BCO) develops in the growth plate (GP) of the proximal femur and tibia and is initiated by damage to the less mineralized chondrocytes followed by colonization of opportunistic bacteria. This condition affects approximately 1% of all birds housed, being considered one of the major causes of lameness in fast growing broilers. Although several studies have been previously performed aiming to understand its pathogenesis, the molecular mechanisms involved with BCO remains to be elucidated. Therefore, this study aimed to generate a profile of global differential gene expression involved with BCO in the tibia of commercial broilers, through RNA sequencing analysis to identity genes and molecular pathways involved with BCO in chickens.ResultsOur data showed 192 differentially expressed (DE) genes: 63 upregulated and 129 downregulated in the GP of the tibia proximal epiphysis of BCO-affected broilers. Using all DE genes, six Biological Processes (BP) were associated with bone development (connective tissue development, cartilage development, skeletal system development, organ morphogenesis, system development and skeletal system morphogenesis). The analyses of the upregulated genes did not indicate any significant BP (FDR < 0.05). However, with the downregulated genes, the same BP were identified when using all DE genes in the analysis, with a total of 26 coding genes explaining BCO in the tibia: ACAN, ALDH1A2, CDH7, CHAD, CHADL, COL11A1, COMP, CSGALNACT1, CYR61, FRZB, GAL3ST1, HAPLN1, IHH, KIF26B, LECT1, LPPR1, PDE6B, RBP4A, SERINC5, SFRP1, SOX8, SOX9, TENM2, THBS1, UCHL1 and WFIKKN2. In addition, seven transcription factors were also associated to BCO: NFATC2, MAFB, HIF1A-ARNT, EWSR1-FLI1, NFIC, TCF3 and NF-KAPPAB.ConclusionsOur data show that osteochondral downregulated genes are potential molecular causes of BCO in broilers, and the bacterial process seems to be, in fact, a secondary condition. Sixteen genes responsible for bone and cartilage formation were downregulated in BCO-affected broilers being strong candidate genes to trigger this disorder.

Heat stress response in the closest algal relatives of land plants reveals conserved stress signaling circuits.

All land plants (embryophytes) share a common ancestor that likely evolved from a filamentous freshwater alga. Elucidating the transition from algae to embryophytes - and the eventual conquering of Earth's surface - is one of the most fundamental questions in plant evolutionary biology. Here, we investigated one of the organismal properties that might have enabled this transition: resistance to drastic temperature shifts. We explored the effect of heat stress in Mougeotia and Spirogyra, two representatives of Zygnematophyceae - the closest known algal sister lineage to land plants. Heat stress induced pronounced phenotypic alterations in their plastids, and high-performance liquid chromatography-tandem mass spectroscopy-based profiling of 565 transitions for the analysis of main central metabolites revealed significant shifts in 43 compounds. We also analyzed the global differential gene expression responses triggered by heat, generating 92.8Gbp of sequence data and assembling a combined set of 8905 well-expressed genes. Each organism had its own distinct gene expression profile; less than one-half of their shared genes showed concordant gene expression trends. We nevertheless detected common signature responses to heat such as elevated transcript levels for molecular chaperones, thylakoid components, and - corroborating our metabolomic data - amino acid metabolism. We also uncovered the heat-stress responsiveness of genes for phosphorelay-based signal transduction that links environmental cues, calcium signatures and plastid biology. Our data allow us to infer the molecular heat stress response that the earliest land plants might have used when facing the rapidly shifting temperature conditions of the terrestrial habitat.

Transcriptome Profile of Nicotinic Receptor-Linked Sensitization of Beta Amyloid Neurotoxicity

Understanding the specific gene changes underlying the prodromic stages of Alzheimer’s disease pathogenesis will aid the development of new, targeted therapeutic strategies for this neurodegenerative disorder. Here, we employed RNA-sequencing to analyze global differential gene expression in a defined model nerve cell line expressing α4β2 nicotinic receptors (nAChRs), high-affinity targets for beta amyloid (Aβ). The nAChR-expressing neuronal cells were treated with nanomolar Aβ1–42 to gain insights into the molecular mechanisms underlying Aβ-induced neurotoxicity in the presence of this sensitizing target receptor. We identified 15 genes (out of 15,336) that were differentially expressed upon receptor-linked Aβ treatment. Genes up-regulated with Aβ treatment were associated with calcium signaling and axonal vesicle transport (including the α4 nAChR subunit, the calcineurin regulator RCAN3, and KIF1C of the kinesin family). Downregulated genes were associated with metabolic, apoptotic or DNA repair pathways (including APBA3, PARP1 and RAB11). Validation of the differential expression was performed via qRT-PCR and immunoblot analysis in the defined model nerve cell line and primary mouse neurons. Further verification was performed using immunocytochemistry. In conclusion, we identified apparent changes in gene expression on Aβ treatment in the presence of the sensitizing nAChRs, linked to early-stage Aβ-induced neurotoxicity, which may represent novel therapeutic targets.