Gene Expression Omnibus Accession Research Articles

MaGENEgerZ: An Efficient Artificial Intelligence-Based Framework Can Extract More Expressed Genes and Biological Insights Underlying Breast Cancer Drug Response Mechanism

Understanding breast cancer drug response mechanisms can play a crucial role in improving treatment outcomes and survival rates. Existing bioinformatics-based approaches are far from perfect and do not adopt computational methods based on advanced artificial intelligence concepts. Therefore, we introduce a novel computational framework based on an efficient support vector machine (esvm) working as follows: First, we downloaded and processed three gene expression datasets related to breast cancer responding and non-responding to treatments from the gene expression omnibus (GEO) according to the following GEO accession numbers: GSE130787, GSE140494, and GSE196093. Our method esvm is formulated as a constrained optimization problem in its dual form as a function of λ. We recover the importance of each gene as a function of λ, y, and x. Then, we select p genes out of n, which are provided as input to enrichment analysis tools, Enrichr and Metascape. Compared to existing baseline methods, including deep learning, results demonstrate the superiority and efficiency of esvm, achieving high-performance results and having more expressed genes in well-established breast cancer cell lines, including MD-MB231, MCF7, and HS578T. Moreover, esvm is able to identify (1) various drugs, including clinically approved ones (e.g., tamoxifen and erlotinib); (2) seventy-four unique genes (including tumor suppression genes such as TP53 and BRCA1); and (3) thirty-six unique TFs (including SP1 and RELA). These results have been reported to be linked to breast cancer drug response mechanisms, progression, and metastasizing. Our method is available publicly on the maGENEgerZ web server.

Exploring personalized treatment for cardiac graft rejection based on a four-archetype analysis model and bioinformatics analysis

Heart transplantation is the gold standard for treating patients with advanced heart failure. Although improvements in immunosuppressive therapies have significantly reduced the frequency of cardiac graft rejection, the incidences of T cell-mediated rejection (TCMR) and antibody-mediated rejection remain almost unchanged. A four-archetype analysis (4AA) model, developed by Philip F. Halloran, illustrated this problem well. It provided a new dimension to improve the accuracy of diagnoses and an independent system for recalibrating the histology guidelines. However, this model was based on the invasive method of endocardial biopsy, which undoubtedly increased the postoperative risk of heart transplant patients. Currently, little is known regarding the associated genes and specific functions of the different phenotypes. We performed bioinformatics analysis (using machine-learning methods and the WGCNA algorithm) to screen for hub-specific genes related to different phenotypes, based Gene Expression Omnibus accession number GSE124897. More immune cell infiltration was observed with the ABMR, TCMR, and injury phenotypes than with the stable phenotype. Hub-specific genes for each of the four archetypes were verified successfully using an external test set (accession number GSE2596). Logistic-regression models based on TCMR-specific hub genes and common hub genes were constructed with accurate diagnostic utility (area under the curve > 0.95). RELA, NFKB1, and SOX14 were identified as transcription factors important for TCMR/injury phenotypes and common genes, respectively. Additionally, 11 Food and Drug Administration-approved drugs were chosen from the DrugBank Database for each four-archetype model. Tyrosine kinase inhibitors may be a promising new option for transplant rejection treatment. KRAS signaling in cardiac transplant rejection is worth further investigation. Our results showed that heart transplant rejection subtypes can be accurately diagnosed by detecting expression of the corresponding specific genes, thereby enabling precise treatment or medication.

Involvement of retinoic acid‑inducible gene‑I in radiation‑induced senescence of human umbilical vein endothelial cells.

In 2012, the threshold radiation dose (0.5 Gy) for cardiovascular and cerebrovascular diseases was revised, and this threshold dose may be exceeded during procedures involving radiation such as interventional radiology. Therefore, in addition to regulating radiation dose, it is necessary to develop strategies to prevent and mitigate the development of cardiovascular disease. Cellular senescence is irreversible arrest of cell proliferation. Although cellular senescence is one of the mechanisms for suppressing cancer, it also has adverse effects. For example, senescence of vascular endothelial cells is involved in development of vascular disorders. However, the mechanisms underlying induction of cellular senescence are not fully understood. Therefore, the present study explored the factors involved in the radiation-induced senescence in human umbilical vein endothelial cells (HUVECs). The present study reanalyzed the gene expression data of senescent normal human endothelial cells and fibroblast after irradiation (NCBI Gene Expression Omnibus accession no. GSE130727) and microarray data of HUVECs 24 h after irradiation (NCBI Gene Expression Omnibus accession no. GSE76484). Numerous genes related to viral infection and inflammation were upregulated in radiation-induced senescent cells. In addition, the gene group involved in the retinoic acid-inducible gene-I (RIG-I)-like receptor (RLR) signaling pathway, which plays an important role to induce anti-viral response, was altered in irradiated HUVECs. Therefore, to investigate the involvement of RIG-I and melanoma differentiation-associated gene 5 (MDA5), which are RLRs, in radiation-induced senescence of HUVECs, the protein expression of RIG-I and MDA5 and the activity of senescence-associated β-galactosidase (SA-β-gal), a representative senescence marker, were analyzed. Of note, knockdown of RIG-I in HUVECs significantly decreased radiation-increased proportion of cells with high SA-β-gal activity (i.e., senescent cells), whereas this phenomenon was not observed in MDA5-knockdown cells. Taken together, the present results suggested that RIG-I, but not MDA5, was associated with radiation-induced senescence in HUVECs.

Formylated Peptide Receptor-1-Mediated Gut Inflammation as a Therapeutic Target in Inflammatory Bowel Disease.

Formylated peptide receptor (FPR)-1 is a G-coupled receptor that senses foreign bacterial and host-derived mitochondrial formylated peptides (FPs), leading to innate immune system activation. We sought to investigate the role of FPR1-mediated inflammation and its potential as a therapeutic target in inflammatory bowel disease (IBD). We characterized FPR1 gene and protein expression in 8 human IBD (~1000 patients) datasets with analysis on disease subtype, mucosal inflammation, and drug response. We performed in vivo dextran-sulfate sodium (DSS) colitis in C57/BL6 FPR1 knockout mice. In ex vivo studies, we studied the role of mitochondrial FPs and pharmacological blockade of FPR1 using cyclosporin H in human peripheral blood neutrophils. Finally, we assess mitochondrial FPs as a potential mechanistic biomarker in the blood and stools of patients with IBD. Detailed in silico analysis in human intestinal biopsies showed that FPR1 is highly expressed in IBD (n = 207 IBD vs 67 non-IBD controls, P < .001), and highly correlated with gut inflammation in ulcerative colitis (UC) and Crohn's disease (CD) (both P < .001). FPR1 receptor is predominantly expressed in leukocytes, and we showed significantly higher FPR1+ve neutrophils in inflamed gut tissue section in IBD (17 CD and 24 UC; both P < .001). Further analysis in 6 independent IBD (data available under Gene Expression Omnibus accession numbers GSE59071, GSE206285, GSE73661, GSE16879, GSE92415, and GSE235970) showed an association with active gut inflammation and treatment resistance to infliximab, ustekinumab, and vedolizumab. FPR1 gene deletion is protective in murine DSS colitis with lower gut neutrophil inflammation. In the human ex vivo neutrophil system, mitochondrial FP, nicotinamide adenine dinucleotide dehydrogenase subunit-6 (ND6) is a potent activator of neutrophils resulting in higher CD62L shedding, CD63 expression, reactive oxygen species production, and chemotactic capacity; these effects are inhibited by cyclosporin H. We screened for mitochondrial ND6 in IBD (n = 54) using ELISA and detected ND6 in stools with median values of 2.2 gg/mL (interquartile range [IQR] 0.0-4.99; range 0-53.3) but not in blood. Stool ND6 levels, however, were not significantly correlated with paired stool calprotectin, C-reactive protein, and clinical IBD activity. Our data suggest that FPR1-mediated neutrophilic inflammation is a tractable target in IBD; however, further work is required to clarify the clinical utility of mitochondrial FPs as a potential mechanistic marker for future stratification.

Baricitinib protects mice from sepsis-induced cardiac dysfunction and multiple-organ failure.

Sepsis is one of the major complications of surgery resulting in high morbidity and mortality, but there are no specific therapies for sepsis-induced organ dysfunction. Data obtained under Gene Expression Omnibus accession GSE131761 were re-analyzed and showed an increased gene expression of Janus Kinase 2 (JAK2) and Signal Transducer and Activator of Transcription 3 (STAT3) in the whole blood of post-operative septic patients. Based on these results, we hypothesized that JAK/STAT activation may contribute to the pathophysiology of septic shock and, hence, investigated the effects of baricitinib (JAK1/JAK2 inhibitor) on sepsis-induced cardiac dysfunction and multiple-organ failure (MOF). In a mouse model of post-trauma sepsis induced by midline laparotomy and cecal ligation and puncture (CLP), 10-week-old male (n=32) and female (n=32) C57BL/6 mice received baricitinib (1mg/kg; i.p.) or vehicle at 1h or 3h post-surgery. Cardiac function was assessed at 24h post-CLP by echocardiography in vivo, and the degree of MOF was analyzed by determination of biomarkers in the serum. The potential mechanism underlying both the cardiac dysfunction and the effect of baricitinib was analyzed by western blot analysis in the heart. Trauma and subsequent sepsis significantly depressed the cardiac function and induced multiple-organ failure, associated with an increase in the activation of JAK2/STAT3, NLRP3 inflammasome and NF- κβ pathways in the heart of both male and female animals. These pathways were inhibited by the administration of baricitinib post the onset of sepsis. Moreover, treatment with baricitinib at 1h or 3h post-CLP protected mice from sepsis-induced cardiac injury and multiple-organ failure. Thus, baricitinib may be repurposed for trauma-associated sepsis.

Differential gene expression and protein-protein interaction network profiling of sulfur mustard-exposed rabbit corneas employing RNA-seq data and bioinformatics tools

Sulfur mustard (SM) ocular exposure severely damages the cornea and causes vision impairment. At present, no specific therapy exists to mitigate SM-induced corneal injury and vision loss. This study performed transcriptome profiling of naïve, SM-damaged, and SM-undamaged rabbit corneas using RNA-seq analysis and bioinformatic tools to gain a better mechanistic understanding and develop SM-specific medical countermeasures. The mRNA profiles of rabbit corneas 4 weeks post SM vapor exposure were generated using Illumina-NextSeq deep sequencing (Gene Expression Omnibus accession # GSE127708). The RNA sequences of naïve (n = 4), SM-damaged (n = 5), and SM-undamaged (n = 5) corneas were subjected to differential expression (DE) analysis after quality control profiling with FastQC. DE analysis was performed using HISAT2, StringTie, and DESeq2. The log2(FC)±2 and adjusted p˂0.05 were chosen to identify the most relevant genes. A total of 5930 differentially expressed genes (DEGs) (upregulated: 3196, downregulated: 2734) were found in SM-damaged corneas compared to naïve corneas, whereas SM-undamaged corneas showed 1884 DEGs (upregulated: 1029, downregulated: 855) compared to naïve corneas. DE profiling of SM-damaged corneas to SM-undamaged corneas revealed 985 genes (upregulated: 308, downregulated: 677). The DE profiles were subsequently subjected to signaling pathway enrichment, and protein‒protein interactions (PPIs) were analyzed. Pathway enrichment was performed for the genes associated with cellular apoptosis, death, adhesion, migration, differentiation, proliferation, extracellular matrix, and tumor necrosis factor production. To identify novel targets, we narrowed the pathway analysis to upregulated and downregulated genes associated with cell proliferation and differentiation, and PPI networks were developed. Furthermore, protein targets associated with cell differentiation and proliferation that may play vital roles in corneal fibrosis and wound healing post SM injury were identified.

Skeletal Muscle's Role in Prenatal Inter-organ Communication: A Phenogenomic Study with Qualitative Citation Analysis.

Gene targeting in mice allows for a complete elimination of skeletal (striated or voluntary) musculature in the body, from the beginning of its development, resulting in our ability to study the consequences of this ablation on other organs. Here I focus on the relationship between the muscle and lung, motor neurons, skeleton, and special senses. Since the inception of my independent laboratory, in 2000, with my team, we published more than 30 papers (and a book chapter), nearly 400 pages of data, on these specific relationships. Here I trace, using Web of Science, nearly 600 citations of this work, to understand its impact. The current report contains a summary of our work and its impact, NCBI's Gene Expression Omnibus accession numbers of all our microarray data, and three clear future directions doable by anyone using our publicly available data. Together, this effort furthers our understanding of inter-organ communication during prenatal development.

Uterine fluid microRNAs are dysregulated in women with recurrent implantation failure.

STUDY QUESTIONIs the composition of microRNAs (miRNAs) in uterine fluid (UF) of women with recurrent implantation failure (RIF) different from that of healthy fertile women?SUMMARY ANSWERThe composition of miRNAs in UF of women with RIF is different from that of healthy fertile women and the dysregulated miRNAs are associated with impaired endometrial receptivity and embryo implantation.WHAT IS KNOWN ALREADYIt has previously been demonstrated that the miRNAs secreted from endometrial cells into the UF contribute to the achievement of endometrial receptivity. Endometrial miRNAs are dysregulated in women with RIF.STUDY DESIGN, SIZE, DURATIONIn this descriptive laboratory case–control study, miRNA abundancy was compared between UF collected during implantation phase from healthy fertile women (n = 17) and women with RIF (n = 34), which was defined as three failed IVF cycles with high-quality embryos.PARTICIPANTS/MATERIALS, SETTING, METHODSRecruitment of study subjects and sampling of UF were performed at two university clinics in Stockholm, Sweden and Tartu, Estonia. The study participants monitored their menstrual cycles using an LH test kit. The UF samples were collected on Day LH + 7–9 by flushing with saline. Samples were processed for small RNA sequencing and mapped for miRNAs. The differential abundance of miRNAs in UF was compared between the two groups using differential expression analysis (DESeq2). Further downstream analyses, including miRNA target gene prediction (miRTarBase), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis (g:Profiler) and external validation using relevant published data, were performed on the dysregulated miRNAs. Two miRNAs were technically validated with quantitative real-time PCR (RT-PCR).MAIN RESULTS AND THE ROLE OF CHANCEAfter processing of the sequencing data, there were 15 samples in the healthy fertile group and 33 samples in the RIF group. We found 61 differentially abundant UF miRNAs (34 upregulated and 27 downregulated) in RIF compared to healthy women with a false discovery rate of <0.05 and a fold change (FC) of −2 or 2. When analyzed with published literature, we found that several of the differentially abundant miRNAs are expressed in endometrial epithelial cells and have been reported in endometrial extracellular vesicles and in association with endometrial receptivity and RIF. Their predicted target genes were further expressed both in the trophectodermal cells of blastocyst-stage embryos and endometrial mid-secretory epithelial cells, as assessed by publicly available single-cell transcriptome-sequencing studies. Pathway analysis further revealed that 25 pathways, having key roles in endometrial receptivity and implantation, were significantly enriched. Hsa-miR-486-5p (FC −20.32; P-value = 0.004) and hsa-miR-92b-3p (FC −9.72; P-value = 0.004) were successfully technically validated with RT-PCR.LARGE SCALE DATAThe data are available in Gene Expression Omnibus (GEO) at https://www.ncbi.nlm.nih.gov/geo/ with GEO accession number: GSE173289.LIMITATIONS, REASONS FOR CAUTIONThis is a descriptive study with a limited number of study participants. Moreover, the identified differentially abundant miRNAs should be validated in a larger study cohort, and the predicted miRNA target genes and enriched pathways in RIF need to be confirmed and further explored in vitro.WIDER IMPLICATIONS OF THE FINDINGSRIF is a major challenge in the current IVF setting with no diagnostic markers nor effective treatment options at hand. For the first time, total miRNAs have been extensively mapped in receptive phase UF of both healthy women with proven fertility and women diagnosed with RIF. Our observations shed further light on the molecular mechanisms behind RIF, with possible implications in future biomarker and clinical treatment studies.STUDY FUNDING/COMPETING INTEREST(S)This work was financially supported by the Swedish Research Council (2017-00932), a joint grant from Region Stockholm and Karolinska Institutet (ALF Medicine 2020, FoUI-954072), Estonian Research Council (PRG1076), Horizon 2020 innovation (ERIN, EU952516) and European Commission and Enterprise Estonia (EU48695). The authors have no competing interests to declare for the current study.

GeoCancerPrognosticDatasetsRetriever: a bioinformatics tool to easily identify cancer prognostic datasets on Gene Expression Omnibus (GEO).

Having multiple datasets is a key aspect of robust bioinformatics analyses, because it allows researchers to find possible confirmation of the discoveries made on multiple cohorts. For this purpose, Gene Expression Omnibus (GEO) can be a useful database, since it provides hundreds of thousands of microarray gene expression datasets freely available for download and usage. Despite this large availability, collecting prognostic datasets of a specific cancer type from GEO can be a long, time-consuming and energy-consuming activity for any bioinformatician, who needs to execute it manually by first performing a search on the GEO website and then by checking all the datasets found one by one. To solve this problem, we present here geoCancerPrognosticDatasetsRetriever, a Perl 5 application which reads a cancer type and a list of microarray platforms, searches for prognostic gene expression datasets of that cancer type and based on those platforms available on GEO, and returns the GEO accession codes of those datasets, if found. Our bioinformatics tool can easily generate in a few minutes a list of cancer prognostic datasets that otherwise would require numerous hours of manual work to any bioinformatician. geoCancerPrognosticDatasetsRetriever can handily retrieve multiple prognostic datasets of gene expression of any cancer type, laying the foundations for numerous bioinformatics studies and meta-analyses that can have a strong impact on oncology research. geoCancerPrognosticDatasetsRetriever is freely available under the GPLv2 license on the Comprehensive Perl Archive Network (CPAN) at https://metacpan.org/pod/App::geoCancerPrognosticDatasetsRetriever and on GitHub at https://github.com/AbbasAlameer/geoCancerPrognosticDatasetsRetriever. Supplementary data are available at Bioinformatics online.

Ribo-Seq and RNA-Seq of TMA46 ( DFRP1) and GIR2 ( DFRP2) knockout yeast strains.

In eukaryotes, stalled and collided ribosomes are recognized by several conserved multicomponent systems, which either block protein synthesis in situ and resolve the collision locally, or trigger a general stress response. Yeast ribosome-binding GTPases RBG1 (DRG1 in mammals) and RBG2 (DRG2) form two distinct heterodimers with TMA46 (DFRP1) and GIR2 (DFRP2), respectively, both involved in mRNA translation. Accumulated evidence suggests that the dimers play partially redundant roles in elongation processivity and resolution of ribosome stalling and collision events, as well as in the regulation of GCN1-mediated signaling involved in ribosome-associated quality control (RQC). They also genetically interact with SLH1 (ASCC3) helicase, a key component of RQC trigger (RQT) complex disassembling collided ribosomes. Here, we present RNA-Seq and ribosome profiling (Ribo-Seq) data from S. cerevisiae strains with individual deletions of the TMA46 and GIR2 genes. Raw RNA-Seq and Ribo-Seq data as well as gene-level readcounts are available in NCBI Gene Expression Omnibus (GEO) repository under GEO accession GSE185458 and GSE185286.

CXCL10 and its related key genes as potential biomarkers for psoriasis: Evidence from bioinformatics and real-time quantitative polymerase chain reaction.

Although several studies have attempted to investigate the etiology of and mechanism underlying psoriasis, the precise molecular mechanism remains unclear. Our study aimed to explore the molecular mechanism underlying psoriasis based on bioinformatics.GSE30999, GSE34248, GSE41662, and GSE50790 datasets were obtained from the Gene Expression Omnibus database. The Gene Expression Omnibus profiles were integrated to obtain differentially expressed genes in R software. Then a series of analyses was performed, such as Gene Ontology annotation, Kyoto Encyclopedia of Genes and Genomes pathway analysis, protein-protein interaction network analysis, among others. The key genes were obtained by CytoHubba, and validated by real-time quantitative polymerase chain reaction.A total of 359 differentially expressed genes were identified between 270 paired lesional and non-lesional skin groups. The common enriched pathways were nucleotide-binding and oligomerization domain-like receptor signaling pathway, and cytokine-cytokine receptor interaction. Seven key genes were identified, including CXCL1, ISG15, CXCL10, STAT1, OASL, IFIT1, and IFIT3. These key genes were validated as upregulated in the 4 datasets and M5-induced HaCaT cells.Our study identified 7 key genes, namely CXCL1, ISG15, CXCL10, STAT1, OASL, IFIT1, and IFIT3, and 2 mostly enriched pathways (nucleotide-binding and oligomerization domain-like receptor signaling pathway, and cytokine-cytokine receptor interaction) involved in psoriatic pathogenesis. More importantly, CXCL1, ISG15, STAT1, OASL, IFIT1, IFIT3, and especially CXCL10 may be potential biomarkers. Therefore, our findings may bring a new perspective to the molecular mechanism underlying psoriasis and suggest potential biomarkers.

Temperature-dependent effects of house fly proto-Y chromosomes on gene expression

Code, commands, and data to perform analyses of house fly RNA-seq data. House fly has a stable polygenic sex determination system. The male determining factor (Mdmd) is commonly found on the Y chromosome (Y^M) or the third chromosome (III^M). These proto-Y chromosomes are clinally distributed, with Y^M found most commonly in northern latitudes and III^M most commonly found at southern latitudes, hinting at possible genotype-by-temperature interactions that maintain the polymorphism. If this distribution is maintained by temperature-dependent selection pressures, we expect the fitness of III^M and Y^M flies to vary across developmental temperatures. These temperature-dependent effects could be driven by differential gene expression in III^M and Y^M males at different temperatures. Here, we performed RNA-seq experiments to study how genotype-by-temperature interactions affect gene expression in male houseflies carrying different proto-Y chromosomes. We raised a Y^M strain known as IsoCS and a III^M strain known as CSrab at 18°C and 29°C for two generations. We dissected 5 heads and 15-20 pairs of testes for each of three replicates of each genotype-by-temperature combination. We carried out RNA-seq on these tissues. Data are available from NCBI Gene Expression Omnibus accession GSE136188.

Toward Colorectal Cancer Biomarkers: The Role of Genetic Variation, Wnt Pathway, and Long Noncoding RNAs.

Colorectal cancer (CRC) is the third leading cause of death worldwide, comprising nearly 8% of cancer-related deaths per year. In South Korea, for example, CRC is the second most common cancer in men, and third in women. This study reports on the association of CRC with genetic variations in long noncoding RNAs, activators, and inhibitors of a cell proliferation pathway. Five normal colon mucosa tissue samples and their matched five-stage IV CRC samples were evaluated (dataset Gene Expression Omnibus accession: GSE50760). We identified more than 5000 differentially expressed genes (DEGs). The Wnt pathway had the greatest portion of DEGs, including activators, inhibitors, and associated long noncoding RNAs (lncRNAs), suggesting the importance of Wnt pathway in CRC. The following genes were aberrantly expressed: WIF1, SFRP4, CD82, WNT2, WNT3, WNT5A, HOTAIR, CRNDE, and UCA1. Notably, HOTAIR is known to silence WIF1, and WIF1 inhibits the Wnt ligands to negatively regulate the pathway. The lncRNA CRNDE positively regulates WNT5A, while UCA1 positively regulates WNT2 and WNT3. We note that HOTAIR was unable to silence WIF1. CRNDE and UCA1 were found to be upregulated, which may explain the high expression of the WIF1 targets. Furthermore, 10 single-nucleotide polymorphisms (SNPs) were identified in five of the candidate genes above. A possible novel SNP in CD82, chr11:44619242T > C, was predicted to introduce a ZBTB7A binding site. These SNPs are hypothesized to contribute to aberrant and discrepant regulation of the Wnt pathway in a context of CRC pathogenesis. These findings collectively inform future research on diagnostics and therapeutics innovation in CRC.

Single-cell RNA sequencing data suggest a role for angiotensin-converting enzyme 2 in kidney impairment in patients infected with 2019-novel coronavirus.

Single-cell RNA sequencing data suggest a role for angiotensin-converting enzyme 2 in kidney impairment in patients infected with 2019-novel coronavirus.

TWIST1 transcriptionally regulates glycolytic genes to promote the Warburg metabolism in pancreatic cancer

Reprogrammed glucose metabolism is essential for tumor initiation and development, especially for pancreatic ductal adenocarcinoma (PDAC). Most cancer cells rely on aerobic glycolysis, a phenomenon termed “the Warburg effect”, to support uncontrolled proliferation and evade apoptosis. However, the direct regulators of the Warburg effect remain areas of active investigation. In this study, we found that the highly conserved transcription factor, TWIST1, is a crucial regulator of aerobic glycolysis in PDAC. Genetic silencing of TWIST1 significantly inhibited the glycolytic phenotypes of PDAC cells as revealed by reduced glucose uptake, lactate production, and extracellular acidification rate, which can be restored by re-expression of siRNA-resistant TWIST1. Moreover, tamoxifen-inducible expression of TWIST1 promoted the Warburg metabolism of PDAC cells. Mechanistically, by luciferase reporter assay and chromatin immunoprecipitation experiment, we showed that TWIST1 can directly increase the expression of several glycolytic genes, including SLC2A1, HK2, ENO1, and PKM2. Of note, the transcriptional regulation by TWIST1 was not dependent on HIF1α or c-Myc. In The Cancer Genome Atlas and Gene Expression Omnibus accession GSE15471, we confirmed that TWIST1 was closely associated with the glycolysis pathway. Collectively, our findings indicate that TWIST1 is likely to act as important regulator of the Warburg effect in PDAC.

Differential expression of long non-coding RNA and mRNA in children with Henoch-Schönlein purpura nephritis.

Long non-coding RNAs (lncRNAs) serve an essential role in regulating immunological functions. However, their impact on Henoch-Schönlein purpura nephritis (HSPN), has remained elusive. The present study determined the expression of lncRNAs and mRNAs in the peripheral blood of 6 children with HSPN and recruited 4 healthy children for comparative study. High-throughput sequencing revealed outstanding differences in lncRNA and mRNA expression, which were verified through reverse transcription-quantitative polymerase chain reaction. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were used to investigate the associated biological functions and possible mechanisms of lncRNAs and mRNAs in HSPN. A total of 820 differentially expressed lncRNAs between the two groups were identified, of which 34 were upregulated and 786 were downregulated. Simultaneously, a total of 3,557 mRNAs were also identified to be differentially expressed, of which 1,232 were upregulated and 2,325 were downregulated. The results revealed that the expression of lncRNAs including ENST00000378432, ENST00000571370, uc001kfc.1 and uc010qna.2 was decreased in HSPN patients compared with that in healthy controls. These lncRNAs were associated with the p53 signaling pathway and apoptosis-associated genes (AKT2, tumor protein 53, phosphatase and tensin homolog and FAS). Further studies of those lncRNAs will be performed to elucidate their functions in apoptosis. Complete raw data files were deposited in the Gene Expression Omnibus (GEO) at National Center for Biotechnology information under the GEO accession no. GSE102114 (www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE102114).

In silico analysis of the potential mechanism of telocinobufagin on breast cancer MCF-7 cells

Backgrounds and aimsThe extractives from a ChanSu, traditional Chinese medicine, have been discovered to possess anti-inflammatory and tumor-suppressing abilities. However, the molecular mechanism of telocinobufagin, a compound extracted from ChanSu, on breast cancer cells has not been clarified. The aim of this study is to investigate the underlying mechanism of telocinobufagin on breast cancer cells. Methods and materialsThe differentially expressed genes after telocinobufagin treatment on breast cancer cells were searched and downloaded from Gene Expression Omnibus (GEO), ArrayExpress and literatures. Bioinformatics tools were applied to further explore the potential mechanism of telocinobufagin in breast cancer using the Kyoto Encyclopedia of genes and genomes (KEGG) pathway, Gene ontology (GO) enrichment, panther, and protein–protein interaction analyses. To better comprehend the role of telocinobufagin in breast cancer, we also queried the Connectivity Map using the gene expression profiles of telocinobufagin treatment. ResultsOne GEO accession (GSE85871) provided 1251 differentially expressed genes after telocinobufagin treatment on MCF-7 cells. The pathway of neuroactive ligand-receptor interaction, cell adhesion molecules (CAMs), intestinal immune network for IgA production, hematopoietic cell lineage and calcium signaling pathway were the key pathways from KEGG analysis. IGF1 and KSR1, owning to higher protein levels in breast cancer tissues, IGF1 and KSR1 could be the hub genes related to telocinobufagin treatment. It was indicated that the molecular mechanism of telocinobufagin resembled that of fenspiride. ConclusionsTelocinobufagin might regulate neuroactive ligand-receptor interaction pathway to exert its influences in breast cancer MCF-7 cells, and its molecular mechanism might share some similarities with fenspiride. This study only presented a comprehensive picture of the role of telocinobufagin in breast cancer MCF-7 cells using big data. However, more thorough and deeper researches are required to add to the validity of this study.

Author Correction: TET proteins safeguard bivalent promoters from de novo methylation in human embryonic stem cells.

In the version of this article initially published, in the Methods, the Gene Expression Omnibus accession code for H3K36me3 ChIP-seq data was incorrectly given asGSM1003585 instead of GSM733725. The error has been corrected in the HTML, PDF and print versions of the article.

Publisher Correction: TET proteins safeguard bivalent promoters from de novo methylation in human embryonic stem cells.

The version of the Supplementary Text and Figures file initially posted was missing Supplementary Tables 1-6 and the Supplementary Note and used incorrect versions of the supplementary figures.

Abstract 4569: Effect of smoking on tumor-infiltrating immune cell composition and prognosis in non-small cell lung cancer

Abstract Introduction: The immune system plays an important role in tumor progression and treatment response such as with new immune modulating therapies. It is unclear in lung cancer what the effect of smoking is upon intratumoral immune cell composition and their function. In our study publicly available expression data were used to find evidence for the influence of smoking on the composition of intratumoral immune cell fractions and their prognostic impact. Methods: We searched the Gene Expression Omnibus (GEO) for human non-small cell lung cancer (NSCLC) samples. Analysis was confined to samples hybridized to the Affymetrix HG-U133 Plus 2.0 (GEO accession number GPL570) platforms. Expression profiles were downloaded and curated for duplicates and corrupt files. Different immune cell fractions were estimated using CIBERSORT and the LM22 leukocyte signature matrix. We combined this information with available clinical data. We used Cox regression analyses to evaluate the prognostic impact of the different immune fractions and Mann-Whitney U tests to test significant differences. Results: Smokers with NSCLC had a significantly (p&amp;lt;0.05) higher percentage of plasma cells, CD8+ T-cells, active CD4+ T-cells, follicular helper cells, monocytes, M0 macrophages, resting dendritic cells, active mast cells and neutrophils compared to non-smokers. Smokers had a lower percentage of naïve B-cells, memory B-cells, resting CD4+ T-cells and resting mast cells. Non-smokers significantly benefit from a higher percentage of intratumoral memory B-cells in the immune infiltrate as compared to normal tissue (HR 0.92, 95% CI 0.85-1.00, p=0,049), while the smokers significantly benefit from higher intratumoral immune factions of resting CD4+ T-cells (HR 0.93, 95% CI 0.89-0.98, p=0.001). Higher percentages of regulatory T-cells (HR 1.23, 95% CI 1.03-1.47, p=0.025), neutrophils (HR 1.11, 95% CI 1.02-1.21, p=0.014) and follicular helper cells (HR 1.12, 95% CI 1.01-1.24, p=0.03) had a significant detrimental effect on survival of smokers. There was a borderline beneficial effect for smokers of the intratumoral NK cell faction (HR 0.61, 95% CI 0.36-1.02, p=0.057). For non-smokers, a detrimental effect was observed for the intratumoral plasma cell fraction (HR 1.04, 95% CI 1.00-1.09, p=0.073) and the naïve B-cell fraction (HR 1.17, 95%CI 0.99-1.39, p=0.075). Conclusion: In NSCLC significant different composition of intratumoral immune cells were observed between a smoking and non-smoking population. These differences were associated with differences in overall survival. Citation Format: Menno Tamminga, T. Jeroen N. Hiltermann, Ed Schuuring, Rudolf S. Fehrmann, Harry J. Groen. Effect of smoking on tumor-infiltrating immune cell composition and prognosis in non-small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4569. doi:10.1158/1538-7445.AM2017-4569