Extensive proteome and functional genomic profiling of variability between genetically identical human B-lymphoblastoid cells.

In life-science research isogenic B-lymphoblastoid cell lines (LCLs) are widely known and preferred for their genetic stability - they are often used for studying mutations for example, where genetic stability is crucial. We have shown previously that phenotypic variability can be observed in isogenic B-lymphoblastoid cell lines. Isogenic LCLs present well-defined phenotypic differences on various levels, for example on the gene expression level or the chromatin level. Based on our investigations, the phenotypic variability of the isogenic LCLs is accompanied by certain genetic variation too. We have developed a compendium of LCL datasets that present the phenotypic and genetic variability of five isogenic LCLs from a multiomic perspective. In this paper, we present additional datasets generated with Next Generation Sequencing techniques to provide genomic and transcriptomic profiles (WGS, RNA-seq, single cell RNA-seq), protein-DNA interactions (ChIP-seq), together with mass spectrometry and flow cytometry datasets to monitor the changes in the proteome. We are sharing these datasets with the scientific community according to the FAIR principles for further investigations.

Potential Role of the Antidepressants Fluoxetine and Fluvoxamine in the Treatment of COVID-19.

Mapping non-canonical cellular pathways affected by approved medications can accelerate drug repurposing efforts, which are crucial in situations with a global impact such as the COVID-19 pandemic. Fluoxetine and fluvoxamine are well-established and widely-used antidepressive agents that act as serotonin reuptake inhibitors (SSRI-s). Interestingly, these drugs have been reported earlier to act as lysosomotropic agents, inhibitors of acid sphingomyelinase in the lysosomes, and as ligands of sigma-1 receptors, mechanisms that might be used to fight severe outcomes of COVID-19. In certain cases, these drugs were administered for selected COVID-19 patients because of their antidepressive effects, while in other cases, clinical studies were performed to assess the effect of these drugs on treating COVID-19 patients. Clinical studies produced promising data that encourage the further investigation of fluoxetine and fluvoxamine regarding their use in COVID-19. In this review, we summarize experimental data and the results of the performed clinical studies. We also provide an overview of previous knowledge on the tissue distribution of these drugs and by integrating this information with the published experimental results, we highlight the real opportunity of using these drugs in our fight against COVID-19.

Histological Grade of Endometrioid Endometrial Cancer and Relapse Risk Can Be Predicted with Machine Learning from Gene Expression Data.

The tumor grade of endometrioid endometrial cancer is used as an independent marker of prognosis and a key component in clinical decision making. It is reported that between grades 1 and 3, however, the intermediate grade 2 carries limited information; thus, patients with grade 2 tumors are at risk of both under- and overtreatment. We used RNA-sequencing data from the TCGA project and machine learning to develop a model which can correctly classify grade 1 and grade 3 samples. We used the trained model on grade 2 patients to subdivide them into low-risk and high-risk groups. With iterative retraining, we selected the most relevant 12 transcripts to build a simplified model without losing accuracy. Both models had a high AUC of 0.93. In both cases, there was a significant difference in the relapse-free survivals of the newly identified grade 2 subgroups. Both models could identify grade 2 patients that have a higher risk of relapse. Our approach overcomes the subjective components of the histological evaluation. The developed method can be automated to perform a prescreening of the samples before a final decision is made by pathologists. Our translational approach based on machine learning methods could allow for better therapeutic planning for grade 2 endometrial cancer patients.

A framework to assess the quality and impact of bioinformatics training across ELIXIR.

ELIXIR is a pan-European intergovernmental organisation for life science that aims to coordinate bioinformatics resources in a single infrastructure across Europe; bioinformatics training is central to its strategy, which aims to develop a training community that spans all ELIXIR member states. In an evidence-based approach for strengthening bioinformatics training programmes across Europe, the ELIXIR Training Platform, led by the ELIXIR EXCELERATE Quality and Impact Assessment Subtask in collaboration with the ELIXIR Training Coordinators Group, has implemented an assessment strategy to measure quality and impact of its entire training portfolio. Here, we present ELIXIR's framework for assessing training quality and impact, which includes the following: specifying assessment aims, determining what data to collect in order to address these aims, and our strategy for centralised data collection to allow for ELIXIR-wide analyses. In addition, we present an overview of the ELIXIR training data collected over the past 4 years. We highlight the importance of a coordinated and consistent data collection approach and the relevance of defining specific metrics and answer scales for consortium-wide analyses as well as for comparison of data across iterations of the same course.

NR2F2 Orphan Nuclear Receptor is Involved in Estrogen Receptor Alpha-Mediated Transcriptional Regulation in Luminal A Breast Cancer Cells.

Nuclear Receptor Subfamily 2 Group F Member 2 (NR2F2) is a member of the steroid/thyroid hormone receptor superfamily with a crucial role in organogenesis, angiogenesis, cardiovascular development and tumorigenesis. However, there is limited knowledge about the cistrome and transcriptome of NR2F2 in breast cancer. In this study, we mapped the regulatory mechanism by NR2F2 using functional genomic methods. To investigate the clinical significance of NR2F2 in breast cancer, The Cancer Genome Atlas (TCGA) data were used. These results show that a high NR2F2 is associated with better survival of a specific subset of patients, namely those with luminal A breast cancer. Therefore, genome-wide NR2F2 and estrogen receptor alpha (ERα) binding sites were mapped in luminal A breast cancer cells using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), revealing that most NR2F2 overlap with ERα that are co-occupied by forkhead box A1 (FOXA1) and GATA binding protein 3 (GATA3) in active enhancer regions. NR2F2 overlaps with highly frequent ERα chromatin interactions, which are essential for the formation of ERα-bound super-enhancers. In the process of the transcriptome profiling of NR2F2-depleted breast cancer cells such differentially expressed genes have been identified that are involved in endocrine therapy resistance and are also ERα target genes. Overall, these findings demonstrate that the NR2F2 nuclear receptor has a key role in ERα-mediated transcription and it can offer a potential therapeutic target in patients with luminal A breast cancer.

Alternatively Constructed Estrogen Receptor Alpha-Driven Super-Enhancers Result in Similar Gene Expression in Breast and Endometrial Cell Lines.

Super-enhancers (SEs) are clusters of highly active enhancers, regulating cell type-specific and disease-related genes, including oncogenes. The individual regulatory regions within SEs might be simultaneously bound by different transcription factors (TFs) and co-regulators, which together establish a chromatin environment conducting to effective transcription. While cells with distinct TF profiles can have different functions, how different cells control overlapping genetic programs remains a question. In this paper, we show that the construction of estrogen receptor alpha-driven SEs is tissue-specific, both collaborating TFs and the active SE components greatly differ between human breast cancer-derived MCF-7 and endometrial cancer-derived Ishikawa cells; nonetheless, SEs common to both cell lines have similar transcriptional outputs. These results delineate that despite the existence of a combinatorial code allowing alternative SE construction, a single master regulator might be able to determine the overall activity of SEs.

COUP-TFII is a modulator of cell-type-specific genetic programs based on genomic localization maps.

Chicken ovalbumin upstream promoter transcription factor II (COUP-TFII) is a member of the steroid/thyroid hormone receptor superfamily, but its ligand has not yet been identified. Little is known about the role of the COUP-TFII nuclear receptor in cancer cells. In this study, we mapped the cistrome of COUP-TFII in three different cancer cells, namely breast cancer cells (MCF-7), myelogenous leukaemia cells (K562) and liver cancer cells (HepG2) using publicly available ChIP-seq data. Our results show that COUP-TFII co-localises with master transcription factors (TFs) in a cell-specific manner such as estrogen receptor alpha in MCF-7, hepatocyte nuclear factor alpha in HepG2, and GATA-binding factor in K562, while the shared, non-specific COUP-TFII binding sites are co-occupied by CTCF. We identified chromatin environments for these COUP-TFII and master TF co-bound sites together with COUP-TFII and CTCF co-bound sites. Our results show that COUP-TFII and master TF co-bound sites are marked with active enhancer specific histone modifications (H3K27ac and H3K4me1), while COUP-TFII and CTCF co-bound sites reveal active promoter specific histone marks (H3K27ac and H3K4me3). These results describe the genomic context and role of COUP-TFII in the cell-type specific transcriptional programs. Furthermore, we report that the VEGFA gene regulated by shared COUP-TFII and CTCF co-bound regulatory elements is involved in long-range looping in a cell-type-independent manner. These findings provide a genomic insight into the regulation and angiogenic role of COUP-TFII.

Classification of different types of estrogen receptor alpha binding sites in MCF-7 cells.

Estrogen Receptor alpha (ERα) is a ligand-activated transcription factor and it has a prominent role in both physiological and pathological processes of the reproductive system. ERα has been investigated extensively in breast cancer and the MCF-7 breast-cancer-derived cell line is a widely used model for the study of its behavior. In this paper we provide a systematic catalog of the possible scenarios of binding to more than 80,000 ERα transcription factor binding sites based on the mechanism of ERα binding to DNA (upon both vehicle and estradiol (E2) treatment). A key feature of the estrogen-driven genetic programs is the presence or absence of the specific response element referred to as the estrogen response element (ERE). While ERα-driven super-enhancers are key components of estrogen-dependent genetic programs, three additional classes of enhancers could be identified: one with the presence of ERE where the ERα bound to the DNA prior of E2-treatment, one where the E2 was required for ERα binding even in the presence of ERE, and one where the ERα binding is established through the response elements of the collaborating factors. Our results suggest that different scenarios of ERα binding result in different genetic programs.

Correction: Lyophilized human cells stored at room temperature preserve multiple RNA species at excellent quality for RNA sequencing.

[This corrects the article DOI: 10.18632/oncotarget.25764.].

Lyophilized human cells stored at room temperature preserve multiple RNA species at excellent quality for RNA sequencing.

Biobanks operating at ambient temperatures would dramatically reduce the costs associated with standard cryogenic storage. In the present study, we used lyophilization to stabilize unfractionated human cells in a dried state at room temperature and tested the yield and integrity of the isolated RNA by microfluidic electrophoresis, RT-qPCR and RNA sequencing. RNA yields and integrity measures were not reduced for lyophilized cells (unstored, stored for two weeks or stored for two months) compared to their paired controls. The abundance of the selected mRNAs with various expression levels, as well as enhancer-associated RNAs and cancer biomarker long non-coding RNAs (MALAT1, GAS5 and TUG1), were not significantly different between the two groups as assessed by RT-qPCR. RNA sequencing data of three lyophilized samples stored for two weeks at room temperature revealed a high degree of similarity with their paired controls in terms of the RNA biotype distribution, cumulative gene diversity, gene body read coverage and per base mismatch rate. Among the 28 differentially expressed genes transcriptional regulators, as well as certain transcript properties suggestive of a residual active decay mechanism were enriched. Our study suggests that freeze-drying of human cells is a suitable alternative for the long-term stabilization of total RNA in whole human cells for routine diagnostics and high-throughput biomedical research.

MicroRNA expression profiles identify disease-specific alterations in systemic lupus erythematosus and primary Sjögren's syndrome.

The discovery of microRNAs (miRNAs) and their critical role in genetic control opened new avenues in understanding of various biological processes including immune cell lineage commitment, differentiation, proliferation and apoptosis. However, a given miRNA may have hundreds of different mRNA targets and a target might be regulated by multiple miRNAs, thus the characterisation of dysregulated miRNA expression profiles could give a better insight into the development of immunological disturbances in autoimmune diseases. The aim of our study was to examine the changes in miRNA expression profiles in patients with systemic lupus erythematosus (SLE) and primary Sjögren's syndrome (pSS). Eight SLE patients, 8 pSS patients and 7 healthy subjects were enrolled in the investigation. MiRNAs were isolated from peripheral blood mononuclear cells, and expression patterns were determined with Illumina next-generation sequencing technology. Since the immunopathogenesis of pSS and SLE encompasses pronounced B cell hyperactivity along with specific autoantibody production, we paid a special attention on the association between miRNA expression levels and altered peripheral B cell distribution. In SLE patients 135, while in pSS patients 26 miRNAs showed altered expression. Interestingly, the 25 miRNAs including miR-146a, miR-16 and miR-21, which were over-expressed in pSS patients, were found to be elevated in SLE group, as well. On the contrary, we observed the down-regulation of miR-150-5p, which is a novel and unique finding in pSS. Levels of several miRNAs over-expressed in SLE, were not changed in pSS, such as miR-148a-3p, miR-152, miR-155, miR-223, miR-224, miR-326 and miR-342. Expression levels of miR-223-5p, miR-150-5p, miR-155-5p and miR-342-3p, which miRNAs are potentially linked to B cell functions, showed associations with the B cell proportions within peripheral blood mononuclear cells. The observed differences in miRNA expression profiles and the better understanding of immune regulatory mechanisms of miRNAs may help to elucidate the pathogenesis of SLE and pSS.

Expression Levels of Warburg-Effect Related microRNAs Correlate with each Other and that of Histone Deacetylase Enzymes in Adult Hematological Malignancies with Emphasis on Acute Myeloid Leukemia.

Disruption of epigenetic regulation and characteristic metabolic alterations (known as the Warburg-effect) are well-known hallmarks of cancer. In our study we investigated the expression levels of microRNAs and histone deacetylase enzymes via RT-qPCR in bone marrow specimens of adult patients suffering from hematological malignancies (total cohort n = 40), especially acute myeloid leukemia (n = 27). The levels of the three examined Warburg-effect related microRNAs (miR-378*, miR-23b, miR-26a) positively correlated with each other and the oncogenic miR-155 and miR-125b, while negatively with the level of the tumorsuppressor miR-124. Significant relationships have been confirmed between the levels of SIRT6, HDAC4 and the microRNAs listed above. In NPM1-mutated AML (n = 6), the level of miR-125b was significantly lower than in the group of AML patients not carrying this mutation (n = 13) (p < 0.05). In M5 FAB type of AML (n = 5), the level of miR-124 was significantly higher compared to the M2 group (n = 7) (p < 0.05). In two cases of FAB M5 AML, the levels of SIRT6 and miR-26a increased during the first 4 weeks of treatment. In the total cohort, white blood cell count at the time of the diagnosis significantly correlated with the levels of HDAC4, SIRT6, miR-124 and miR-26a. Our results suggest that Warburg-effect related microRNAs may have important role in the pathogenesis of leukemia, and the potential oncogenic property of HDAC4 and SIRT6 cannot be excluded in hematological malignancies. Elevated level of miR-125b can contribute to adverse prognosis of AML without NPM1 mutation. The prevailment of the tumorsuppressor property of miR-124 may depend on the accompanying genetic alterations.

Improved transgene expression in doxycycline-inducible embryonic stem cells by repeated chemical selection or cell sorting.

Transgene-mediated programming is a preeminent strategy to direct cellular identity. To facilitate cell fate switching, lineage regulating genes must be efficiently and uniformly induced. However, gene expression is often heterogeneous in transgenic systems. Consistent with this notion, a non-uniform reporter gene expression was detected in our doxycycline (DOX)-regulated, murine embryonic stem (ES) cell clones. Interestingly, a significant fraction of cells within each clone failed to produce any reporter signals upon DOX treatment. We found that the majority of these non-responsive cells neither carry reporter transgene nor geneticin/G418 resistance. This observation suggested that our ES cell clones contained non-recombined cells that survived the G418 selection which was carried out during the establishment of these clones. We successfully eliminated most of these corrupted cells with repeated chemical (G418) selection, however, even after prolonged G418 treatments, a few cells remained non-responsive due to epigenetic silencing. We found that cell sorting has been the most efficient approach to select those cells which can uniformly and stably induce the integrated transgene in this ES cell based platform. Together, our data revealed that post-cloning chemical re-selection or cell sorting strongly facilitate the production of ES cell lines with a uniform transgene induction capacity.

Effect of Bleomycin Hydrolase Gene Polymorphism on Late Pulmonary Complications of Treatment for Hodgkin Lymphoma.

BACKGROUND: Bleomycin hydrolase (BLMH), an enzyme that inactivates bleomycin, may be a potential candidate that could influence pulmonary function in ABVD (doxorubicin, bleomycin, vinblastin, dacarbasine)-treated Hodgkin lymphoma (HL) patients. PATIENTS AND METHODS: We hypothesized that the BLMH gene SNP A1450G (rs1050565) influences BLMH activity and late pulmonary toxicity. St. George Respiratory Questionnaire, lung scintigraphy and spirometry were used to determine lung function. TaqMan genotyping assay was used to determine genotype distribution of 131 previously treated HL patients. RESULTS: Significantly more favorable results were seen in the wild-type A/A genotype group than those in the group containing the mutated allele: A/G+G/G in retrospective pulmonary tests of ABVD treated patients. CONCLUSION: Besides limitations of the current study, bleomycin pharmacokinetics should be further evaluated in patients with BLMH variations, hence identify those cases even in the frontline setting, where bleomycin should be omitted and replaced with targeted therapy.

STAT2/IRF9 directs a prolonged ISGF3-like transcriptional response and antiviral activity in the absence of STAT1.

Evidence is accumulating for the existence of a signal transducer and activator of transcription 2 (STAT2)/interferon regulatory factor 9 (IRF9)-dependent, STAT1-independent interferon alpha (IFNα) signalling pathway. However, no detailed insight exists into the genome-wide transcriptional regulation and the biological implications of STAT2/IRF9-dependent IFNα signalling as compared with interferon-stimulated gene factor 3 (ISGF3). In STAT1-defeicient U3C cells stably overexpressing human STAT2 (hST2-U3C) and STAT1-deficient murine embryonic fibroblast cells stably overexpressing mouse STAT2 (mST2-MS1KO) we observed that the IFNα-induced expression of 2'-5'-oligoadenylate synthase 2 (OAS2) and interferon-induced protein with tetratricopeptide repeats 1 (Ifit1) correlated with the kinetics of STAT2 phosphorylation, and the presence of a STAT2/IRF9 complex requiring STAT2 phosphorylation and the STAT2 transactivation domain. Subsequent microarray analysis of IFNα-treated wild-type (WT) and STAT1 KO cells overexpressing STAT2 extended our observations and identified ∼120 known antiviral ISRE-containing interferon-stimulated genes (ISGs) commonly up-regulated by STAT2/IRF9 and ISGF3. The STAT2/IRF9-directed expression profile of these IFN-stimulated genes (ISGs) was prolonged as compared with the early and transient response mediated by ISGF3. In addition, we identified a group of 'STAT2/IRF9-specific' ISGs, whose response to IFNα was ISGF3-independent. Finally, STAT2/IRF9 was able to trigger an antiviral response upon encephalomyocarditis virus (EMCV) and vesicular stomatitis Indiana virus (VSV). Our results further prove that IFNα-activated STAT2/IRF9 induces a prolonged ISGF3-like transcriptome and generates an antiviral response in the absence of STAT1. Moreover, the existence of 'STAT2/IRF9-specific' target genes predicts a novel role of STAT2 in IFNα signalling.

Hmgb1 can facilitate activation of the matrilin-1 gene promoter by Sox9 and L-Sox5/Sox6 in early steps of chondrogenesis.

The architectural high mobility group box 1 (Hmgb1) protein acts as both a nuclear and an extracellular regulator of various biological processes, including skeletogenesis. Here we report its contribution to the evolutionarily conserved, distinctive regulation of the matrilin-1 gene (Matn1) expression in amniotes. We previously demonstrated that uniquely assembled proximal promoter elements restrict Matn1 expression to specific growth plate cartilage zones by allowing varying doses of L-Sox5/Sox6 and Nfi proteins to fine-tune their Sox9-mediated transactivation. Here, we dissected the regulatory mechanisms underlying the activity of a conserved distal promoter element 1. We show that this element carries three Sox-binding sites, works as an enhancer in vivo, and allows promoter activation by the Sox5/6/9 chondrogenic trio. In early steps of chondrogenesis, declining Hmgb1 expression overlaps with the onset of Sox9 expression. Unlike repression in late steps, Hmgb1 overexpression in early chondrogenesis increases Matn1 promoter activation by the Sox trio, and forced Hmgb1 expression in COS-7 cells facilitates induction of Matn1 expression by the Sox trio. The conserved Matn1 control elements bind Hmgb1 and SOX9 with opposite efficiency in vitro. They show higher HMGB1 than SOX trio occupancy in established chondrogenic cell lines, and HMGB1 silencing greatly increases MATN1 and COL2A1 expression. Together, these data thus suggest a model whereby Hmgb1 helps recruit the Sox trio to the Matn1 promoter and thereby facilitates activation of the gene in early chondrogenesis. We anticipate that Hmgb1 may similarly affect transcription of other cartilage-specific genes.

Activation of retinoic acid receptor signaling coordinates lineage commitment of spontaneously differentiating mouse embryonic stem cells in embryoid bodies.

Retinoid signaling has been implicated in embryonic stem cell differentiation. Here we present a systematic analysis of gene expression changes in mouse embryonic stem cells (mESCs), during their spontaneous differentiation into embryoid bodies and the effect of all-trans retinoic acid (ATRA) on this process. We show that retinoic acid is present in the serum and is sufficient to activate retinoid signaling at a basal level in undifferentiated mESCs. This signal disappears during embryoid body formation. However exogenously added ATRA resets the spontaneous differentiation programs in embryoid bodies and initiates a distinct genetic program. These data suggest that retinoid signaling not only promotes a particular pathway but also acts as a context dependent general coordinator of the differentiation states in embryonic stem cells.