Relationship of micro-RNA, mRNA and eIF Expression in Tamoxifen-Adapted MCF-7 Breast Cancer Cells: Impact of miR-1972 on Gene Expression, Proliferation and Migration.

BACKGROUND: Tamoxifen-adapted MCF-7-Tam cells represent an in-vitro model for acquired tamoxifen resistance, which is still a problem in clinics. We here investigated the correlation of microRNA-, mRNA- and eukaryotic initiation factors (eIFs) expression in this model. METHODS: MicroRNA- and gene expression were analyzed by nCounter and qRT-PCR technology; eIFs by Western blotting. Protein translation mode was determined using a reporter gene assay. Cells were transfected with a miR-1972-mimic. RESULTS: miR-181b-5p,-3p and miR-455-5p were up-, miR-375, and miR-1972 down-regulated and are significant in survival analysis. About 5% of the predicted target genes were significantly altered. Pathway enrichment analysis suggested a contribution of the FoxO1 pathway. The ratio of polio-IRES driven to cap-dependent protein translation shifted towards cap-dependent initiation. Protein expression of eIF2A, -4G, -4H and -6 decreased, whereas eIF3H was higher in MCF-7-Tam. Significant correlations between tamoxifen-regulated miRNAs and eIFs were found in representative breast cancer cell lines. Transfection with a miR-1972-mimic reverses tamoxifen-induced expression for a subset of genes and increased proliferation in MCF-7, but reduced proliferation in MCF-7-Tam, especially in the presence of 4OH-tamoxifen. Migration was inhibited in MCF-7-Tam cells. Translation mode remained unaffected. CONCLUSIONS: miR-1972 contributes to the orchestration of gene-expression and physiological consequences of tamoxifen adaption.

Importance of RpoD- and Non-RpoD-Dependent Expression of Horizontally Acquired Genes in Cupriavidus metallidurans.

The genome of the metal-resistant, hydrogen-oxidizing bacterium Cupriavidus metallidurans contains a large number of horizontally acquired plasmids and genomic islands that were integrated into its chromosome or chromid. For the C. metallidurans CH34 wild-type strain growing under nonchallenging conditions, 5,763 transcriptional starting sequences (TSSs) were determined. Using a custom-built motif discovery software based on hidden Markov models, patterns upstream of the TSSs were identified. The pattern TTGACA, -35.6 ± 1.6 bp upstream of the TSSs, in combination with a TATAAT sequence 15.8 ± 1.4 bp upstream occurred frequently, especially upstream of the TSSs for 48 housekeeping genes, and these were assigned to promoters used by RNA polymerase containing the main housekeeping sigma factor RpoD. From patterns upstream of the housekeeping genes, a score for RpoD-dependent promoters in C. metallidurans was derived and applied to all 5,763 TSSs. Among these, 2,572 TSSs could be associated with RpoD with high probability, 373 with low probability, and 2,818 with no probability. In a detailed analysis of horizontally acquired genes involved in metal resistance and not involved in this process, the TSSs responsible for the expression of these genes under nonchallenging conditions were assigned to RpoD- or non-RpoD-dependent promoters. RpoD-dependent promoters occurred frequently in horizontally acquired metal resistance and other determinants, which should allow their initial expression in a new host. However, other sigma factors and sense/antisense effects also contribute-maybe to mold in subsequent adaptation steps the assimilated gene into the regulatory network of the cell. IMPORTANCE In their natural environment, bacteria are constantly acquiring genes by horizontal gene transfer. To be of any benefit, these genes should be expressed. We show here that the main housekeeping sigma factor RpoD plays an important role in the expression of horizontally acquired genes in the metal-resistant hydrogen-oxidizing bacterium C. metallidurans. By conservation of the RpoD recognition consensus sequence, a newly arriving gene has a high probability to be expressed in the new host cell. In addition to integrons and genes travelling together with that for their sigma factor, conservation of the RpoD consensus sequence may be an important contributor to the overall evolutionary success of horizontal gene transfer in bacteria. Using C. metallidurans as an example, this publication sheds some light on the fate and function of horizontally acquired genes in bacteria.

epiGBS2: Improvements and evaluation of highly multiplexed, epiGBS-based reduced representation bisulfite sequencing.

Several reduced-representation bisulfite sequencing methods have been developed in recent years to determine cytosine methylation de novo in nonmodel species. Here, we present epiGBS2, a laboratory protocol based on epiGBS with a revised and user-friendly bioinformatics pipeline for a wide range of species with or without a reference genome. epiGBS2 is cost- and time-efficient and the computational workflow is designed in a user-friendly and reproducible manner. The library protocol allows a flexible choice of restriction enzymes and a double digest. The bioinformatics pipeline was integrated in the Snakemake workflow management system, which makes the pipeline easy to execute and modular, and parameter settings for important computational steps flexible. We implemented bismark for alignment and methylation analysis and we preprocessed alignment files by double masking to enable single nucleotide polymorphism calling with Freebayes (epiFreebayes). The performance of several critical steps in epiGBS2 was evaluated against baseline data sets from Arabidopsis thaliana and great tit (Parus major), which confirmed its overall good performance. We provide a detailed description of the laboratory protocol and an extensive manual of the bioinformatics pipeline, which is publicly accessible on github (https://github.com/nioo-knaw/epiGBS2) and zenodo (https://doi.org/10.5281/zenodo.4764652).

Transcriptional Analysis of Total CD8+ T Cells and CD8+CD45RA- Memory T Cells From Young and Old Healthy Blood Donors.

Memory CD8+ T cells accumulate with aging, while the naïve T cell compartment decreases, leading to an increased susceptibility to infections and a decreased vaccine efficiency. To get deeper insights into the underlying mechanisms, this study aims to determine the age-dependent expression profile of total versus memory CD8+ T cells from young and old donors. Total CD8+ and CD8+CD45RA- memory T cells isolated from young (<30 years) and old (>60 years) donors were stimulated with anti-CD3 and anti-CD28 antibodies for 48h before analyzing the cytokine secretion and activation markers by flow cytometry and changes in the expression profiles using RNA sequencing. Gene ontology (GO) term enrichment analyses were performed for up-regulated and uniquely expressed transcripts identified in the T cell populations of both age groups. Total and memory CD8+ T cells from old donors expressed significantly higher CD25 levels and have an increased cytokine secretion. While approximately 1,500 up-regulated transcripts were identified in all groups, CD8+CD45RA- memory T cells of old donors had approximately 500 more uniquely expressed transcripts. Four GO terms related to the JAK-STAT pathway were identified for up-regulated transcripts in the total CD8+ T cells of old donors, whereas CD8+CD45RA- memory T cells GO terms related to adjacent pathways, like JNK and MAPK/ERK, were found. Additionally, the unique transcripts of CD8+CD45RA- memory T cells of old donors were related to the JNK, MAPK and IL-12 pathways. For both T cell populations of the old donors, cytokine and JAK-STAT pathway transcripts were up-regulated. Thus, an age-dependent effect was observed on the transcriptomes of total and memory CD8+ T cells. The CD8+ CD45RA- memory T cells from old donors maintained the increased cytokine secretion of the total CD8+ T cell population and the increased JAK-STAT pathway transcripts, which have an impact on inflammation and senescence.

Acute Effects of an Inorganic Phosphorus Additive on Mineral Metabolism and Cardiometabolic Risk Factors in Healthy Subjects.

CONTEXT: Hyperphosphatemia and high levels of fibroblast growth factor 23 (FGF23) are risk factors for cardiovascular events in patients with chronic kidney diseases. However, the impact of an inorganic phosphorus additive in healthy people is largely unknown. OBJECTIVE: We aimed to investigate the acute effect of excessive dietary phosphorus administered as sodium dihydrogen phosphate on the postprandial levels of Pi and FGF23 and the response to food. METHODS: This study was a double-blind placebo-controlled crossover study with 29 healthy male and female participants from the general community who were administered a single dose of either 700 mg phosphorus (NaH2PO4) or a sodium-adjusted placebo in combination with a test meal. Postprandial plasma levels of Pi and FGF23 were measured. RESULTS: Compared with placebo, oral phosphorus increased the plasma Pi level, which remained elevated during the ensuing 8 hours (at 480 minutes: 1.31 vs 1.16 mmol/l; P < 0.001), increased urinary Pi (iAUC0-480 789 vs 95 mmol/mmol; P < 0.001), reduced tubular Pi reabsorption (iAUC0-480 -31.5 vs -6.2; P < 0.001), decreased urinary calcium (iAUC0-240 30.6 vs 53.0 mmol/mmol; P = 0.009), and stimulated the release of parathyroid hormone (iAUC0-480 2212 vs 768 ng/l; P < 0.001). However, the FGF23 levels did not change. Postprandial levels of glucose, insulin, and lipids were not substantially affected by phosphorus vs placebo. CONCLUSION: An oral phosphorus load can induce elevated postprandial levels of circulating Pi for hours in healthy subjects, despite rapid homeostatic counterreactions. FGF23 levels and the postprandial response to food were not affected.

SESN2 Knockdown Increases Betulinic Acid-Induced Radiosensitivity of Hypoxic Breast Cancer Cells.

Betulinic acid (BA) is a natural compound well known for its anti-inflammatory, anti-viral, anti-bacterial, anti-malarial effects and anti-tumor properties. Its enhanced cytotoxicity in tumor cells and induction of cell death in various cancer entities qualifies BA as an interesting candidate for novel treatment concepts. Our analyses showed enhanced cytotoxicity and radiosensitization under hypoxic conditions in human breast cancer cells. So far, the underlying mechanisms are unknown. Therefore, we investigated the BA-treated human breast cancer cell lines MDA-MB-231 and MCF-7 under normoxic and hypoxic conditions based on microarray technology. Hypoxia and BA regulated a variety of genes in both breast cancer cell lines. KEGG pathway analysis identified an enrichment of the p53 pathway in MCF-7 cells (wtp53) under hypoxia. In MDA-MB-231 cells (mtp53) an additional BA incubation was required to activate the p53 signaling pathway. Fourteen down-regulated and up-regulated genes of the p53 pathway were selected for further validation via qRT-PCR in a panel of five breast cancer cell lines. The stress-induced gene Sestrin-2 (SESN2) was identified as one of the most strongly up-regulated genes after BA treatment. Knockdown of SESN2 enhanced BA-induced ROS production, DNA damage, radiosensitivity and reduced autophagy in breast cancer cells. Our results identified SESN2 as an important target to enhance the radiobiological and anti-tumor effects of BA on breast cancer cells.

Detection of SARS-CoV-2 Derived Small RNAs and Changes in Circulating Small RNAs Associated with COVID-19.

Cleavage of double-stranded RNA is described as an evolutionary conserved host defense mechanism against viral infection. Small RNAs are the product and triggers of post transcriptional gene silencing events. Up until now, the relevance of this mechanism for SARS-CoV-2-directed immune responses remains elusive. Herein, we used high throughput sequencing to profile the plasma of active and convalescent COVID-19 patients for the presence of small circulating RNAs. The existence of SARS-CoV-2 derived small RNAs in plasma samples of mild and severe COVID-19 cases is described. Clusters of high siRNA abundance were discovered, homologous to the nsp2 3'-end and nsp4 virus sequence. Four virus-derived small RNA sequences have the size of human miRNAs, and a target search revealed candidate genes associated with ageusia and long COVID symptoms. These virus-derived small RNAs were detectable also after recovery from the disease. The additional analysis of circulating human miRNAs revealed differentially abundant miRNAs, discriminating mild from severe cases. A total of 29 miRNAs were reduced or absent in severe cases. Several of these are associated with JAK-STAT response and cytokine storm.

tSFM 1.0: tRNA Structure-Function Mapper.

MOTIVATION: Structure-conditioned information statistics have proven useful to predict and visualize tRNA Class-Informative Features (CIFs) and their evolutionary divergences. Although permutation P-values can quantify the significance of CIF divergences between two taxa, their naive Monte Carlo approximation is slow and inaccurate. The Peaks-over-Threshold approach of Knijnenburg et al. (2009) promises improvements to both speed and accuracy of permutation P-values, but has no publicly available API. RESULTS: We present tRNA Structure-Function Mapper (tSFM) v1.0, an open-source, multi-threaded application that efficiently computes, visualizes and assesses significance of single- and paired-site CIFs and their evolutionary divergences for any RNA, protein, gene or genomic element sequence family. Multiple estimators of permutation P-values for CIF evolutionary divergences are provided along with confidence intervals. tSFM is implemented in Python 3 with compiled C extensions and is freely available through GitHub (https://github.com/tlawrence3/tSFM) and PyPI. AVAILABILITY AND IMPLEMENTATION: The data underlying this article are available on GitHub at https://github.com/tlawrence3/tSFM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Taxonomic analysis of metagenomic data with kASA.

The taxonomic analysis of sequencing data has become important in many areas of life sciences. However, currently available tools for that purpose either consume large amounts of RAM or yield insufficient quality and robustness. Here, we present kASA, a k-mer based tool capable of identifying and profiling metagenomic DNA or protein sequences with high computational efficiency and a user-definable memory footprint. We ensure both high sensitivity and precision by using an amino acid-like encoding of k-mers together with a range of multiple k's. Custom algorithms and data structures optimized for external memory storage enable a full-scale taxonomic analysis without compromise on laptop, desktop, and HPCC.

Identifying High-Risk Triple-Negative Breast Cancer Patients by Molecular Subtyping.

INTRODUCTION: Triple-negative breast cancer (TNBC) is considered the most aggressive type of breast cancer (BC) with limited options for therapy. TNBC is a heterogeneous disease and tumors have been classified into TNBC subtypes using gene expression profiling to distinguish basal-like 1, basal-like 2, immunomodulatory, mesenchymal, mesenchymal stem-like, luminal androgen receptor (LAR), and one nonclassifiable group (called unstable). OBJECTIVES: The aim of this study was to verify the clinical relevance of molecular subtyping of TNBCs to improve the individual indication of systemic therapy. PATIENTS AND METHODS: Molecular subtyping was performed in 124 (82%) of 152 TNBC tumors that were obtained from a prospective, multicenter cohort including 1,270 histopathologically confirmed invasive, nonmetastatic BCs (NCT01592825). Treatment was guideline-based. TNBC subtypes were correlated with recurrence-free interval (RFI) and overall survival (OS) after 5 years of observation. RESULTS: Using PAM50 analysis, 87% of the tumors were typed as basal with an inferior clinical outcome compared to patients with nonbasal tumors. Using the TNBCtype-6 classifier, we identified 23 (15%) of TNBCs as LAR subtype. After standard adjuvant or neoadjuvant chemotherapy, patients with LAR subtype showed the most events for 5-year RFI (66.7 vs. 80.6%) and the poorest probability of 5-year OS (60.0 vs. 84.4%) compared to patients with non-LAR disease (RFI: adjusted hazard ratio [aHR] = 1.87, 95% confidence interval [CI] 0.69-5.05, p = 0.211; OS: aHR = 2.74, 95% CI 1.06-7.10, p = 0.037). CONCLUSION: Molecular analysis and subtyping of TNBC may be relevant to identify patients with LAR subtype. These cancers seem to be less sensitive to conventional chemotherapy, and new treatment options, including androgen receptor-blocking agents and immune checkpoint inhibitors, have to be explored.

Current Understanding of the HIF-1-Dependent Metabolism in Oral Squamous Cell Carcinoma.

Oral squamous cell carcinoma (OSCC) is the 10th most frequent human malignancy and is thus a global burden. Despite some progress in diagnosis and therapy, patients' overall survival rate, between 40 and 55%, has stagnated over the last four decades. Since the tumor node metastasis (TNM) system is not precise enough to predict the disease outcome, additive factors for diagnosis, prognosis, prediction and therapy resistance are urgently needed for OSCC. One promising candidate is the hypoxia inducible factor-1 (HIF-1), which functions as an early regulator of tumor aggressiveness and is a key promoter of energy adaptation. Other parameters comprise the composition of the tumor microenvironment, which determines the availability of nutrients and oxygen. In our opinion, these general processes are linked in the pathogenesis of OSCC. Based on this assumption, the review will summarize the major features of the HIF system-induced activities, its target proteins and related pathways of nutrient utilization and metabolism that are essential for the initiation, progression and therapeutic stratification of OSCC.

Insights gained from a comprehensive all-against-all transcription factor binding motif benchmarking study.

BACKGROUND: Positional weight matrix (PWM) is a de facto standard model to describe transcription factor (TF) DNA binding specificities. PWMs inferred from in vivo or in vitro data are stored in many databases and used in a plethora of biological applications. This calls for comprehensive benchmarking of public PWM models with large experimental reference sets. RESULTS: Here we report results from all-against-all benchmarking of PWM models for DNA binding sites of human TFs on a large compilation of in vitro (HT-SELEX, PBM) and in vivo (ChIP-seq) binding data. We observe that the best performing PWM for a given TF often belongs to another TF, usually from the same family. Occasionally, binding specificity is correlated with the structural class of the DNA binding domain, indicated by good cross-family performance measures. Benchmarking-based selection of family-representative motifs is more effective than motif clustering-based approaches. Overall, there is good agreement between in vitro and in vivo performance measures. However, for some in vivo experiments, the best performing PWM is assigned to an unrelated TF, indicating a binding mode involving protein-protein cooperativity. CONCLUSIONS: In an all-against-all setting, we compute more than 18 million performance measure values for different PWM-experiment combinations and offer these results as a public resource to the research community. The benchmarking protocols are provided via a web interface and as docker images. The methods and results from this study may help others make better use of public TF specificity models, as well as public TF binding data sets.

Fold-Change-Specific Enrichment Analysis (FSEA): Quantification of Transcriptional Response Magnitude for Functional Gene Groups.

Gene expression profiling data contains more information than is routinely extracted with standard approaches. Here we present Fold-Change-Specific Enrichment Analysis (FSEA), a new method for functional annotation of differentially expressed genes from transcriptome data with respect to their fold changes. FSEA identifies Gene Ontology (GO) terms, which are shared by the group of genes with a similar magnitude of response, and assesses these changes. GO terms found by FSEA are fold-change-specifically (e.g., weakly, moderately, or strongly) affected by a stimulus under investigation. We demonstrate that many responses to abiotic factors, mutations, treatments, and diseases occur in a fold-change-specific manner. FSEA analyses suggest that there are two prevailing responses of functionally-related gene groups, either weak or strong. Notably, some of the fold-change-specific GO terms are invisible by classical algorithms for functional gene enrichment, Singular Enrichment Analysis (SEA), and Gene Set Enrichment Analysis (GSEA). These are GO terms not enriched compared to the genome background but strictly regulated by a factor within specific fold-change intervals. FSEA analysis of a cancer-related transcriptome suggested that the gene groups with a tightly coordinated response can be the valuable source to search for possible regulators, markers, and therapeutic targets in oncogenic processes. Availability and Implementation: FSEA is implemented as the FoldGO Bioconductor R package and a web-server.

Performance of Mapping Approaches for Whole-Genome Bisulfite Sequencing Data in Crop Plants.

DNA methylation is involved in many different biological processes in the development and well-being of crop plants such as transposon activation, heterosis, environment-dependent transcriptome plasticity, aging, and many diseases. Whole-genome bisulfite sequencing is an excellent technology for detecting and quantifying DNA methylation patterns in a wide variety of species, but optimized data analysis pipelines exist only for a small number of species and are missing for many important crop plants. This is especially important as most existing benchmark studies have been performed on mammals with hardly any repetitive elements and without CHG and CHH methylation. Pipelines for the analysis of whole-genome bisulfite sequencing data usually consists of four steps: read trimming, read mapping, quantification of methylation levels, and prediction of differentially methylated regions (DMRs). Here we focus on read mapping, which is challenging because un-methylated cytosines are transformed to uracil during bisulfite treatment and to thymine during the subsequent polymerase chain reaction, and read mappers must be capable of dealing with this cytosine/thymine polymorphism. Several read mappers have been developed over the last years, with different strengths and weaknesses, but their performances have not been critically evaluated. Here, we compare eight read mappers: Bismark, BismarkBwt2, BSMAP, BS-Seeker2, Bwameth, GEM3, Segemehl, and GSNAP to assess the impact of the read-mapping results on the prediction of DMRs. We used simulated data generated from the genomes of Arabidopsis thaliana, Brassica napus, Glycine max, Solanum tuberosum, and Zea mays, monitored the effects of the bisulfite conversion rate, the sequencing error rate, the maximum number of allowed mismatches, as well as the genome structure and size, and calculated precision, number of uniquely mapped reads, distribution of the mapped reads, run time, and memory consumption as features for benchmarking the eight read mappers mentioned above. Furthermore, we validated our findings using real-world data of Glycine max and showed the influence of the mapping step on DMR calling in WGBS pipelines. We found that the conversion rate had only a minor impact on the mapping quality and the number of uniquely mapped reads, whereas the error rate and the maximum number of allowed mismatches had a strong impact and leads to differences of the performance of the eight read mappers. In conclusion, we recommend BSMAP which needs the shortest run time and yields the highest precision, and Bismark which requires the smallest amount of memory and yields precision and high numbers of uniquely mapped reads.

Selective egg cell polyspermy bypasses the triploid block.

Polyploidization, the increase in genome copies, is considered a major driving force for speciation. We have recently provided the first direct in planta evidence for polyspermy induced polyploidization. Capitalizing on a novel sco1-based polyspermy assay, we here show that polyspermy can selectively polyploidize the egg cell, while rendering the genome size of the ploidy-sensitive central cell unaffected. This unprecedented result indicates that polyspermy can bypass the triploid block, which is an established postzygotic polyploidization barrier. In fact, we here show that most polyspermy-derived seeds are insensitive to the triploid block suppressor admetos. The robustness of polyspermy-derived plants is evidenced by the first transcript profiling of triparental plants and our observation that these idiosyncratic organisms segregate tetraploid offspring within a single generation. Polyspermy-derived triparental plants are thus comparable to triploids recovered from interploidy crosses. Our results expand current polyploidization concepts and have important implications for plant breeding.

Ever since Darwin published his most famous book on the theory of evolution, scientists have sought to identify the mechanisms that drive the formation of new species. This is especially true for plant biologists who have long been fascinated by the extraordinary diversity of flowering plants.Many species of flowering plant first evolved after a dramatic increase in the DNA content of an individual plant, a process termed polyploidization. Most explanations for polyploidization involve a pollen grain making sperm that mistakenly contain two sets of chromosomes rather than one. Yet, it is difficult to reconcile this explanation with an important aspect of plant reproduction ­ the so-called "triploid block".Fertilization in flowering plants is more complicated than in animals. While one sperm fertilizes the egg cell to make the plant embryo, a second sperm from the same pollen grain must fertilize another cell to form the endosperm, the tissue that will nourish the embryo as it develops. This means that sperm with twice the normal number of chromosomes would affect the DNA content of both the embryo and the endosperm. Yet, an endosperm that receives extra paternal DNA typically halts the development of the seed via a process known as the triploid block, meaning it was not clear how often this process would actually result in a polyploid plant.In 2017, researchers reported that plants can, on rare occasions, generate polyploid offspring via a different route: the fertilization of one egg with two sperm rather than one. Now, Mao et al. ­ who include several researchers involved in the 2017 study ­ show that this process, termed "polyspermy", can introduce extra copies of DNA into just the egg cell, meaning it can bypass the triploid block of the endosperm.The experiments involved a model plant called Arabidopsis, and a screen of over 55,000 seeds identified about a dozen with embryos that had three parents, one mother and two fathers. Notably, most of these three-parent embryos developed in seeds that contained endosperm with the regular number of chromosomes and hence escaped the triploid block.These new results show that polyspermy provides plants with a means to essentially sneak extra copies of DNA 'behind the back' of the DNA-sensitive endosperm and into the next generation. They also give new insight in how polyploidization may have shaped the evolution of flowering plants and have important implications for agriculture where the breeding of new "hybrid" crops has often been limited by incompatibilities in the endosperm.

Urban areas as hotspots for bees and pollination but not a panacea for all insects.

Urbanisation is an important global driver of biodiversity change, negatively impacting some species groups whilst providing opportunities for others. Yet its impact on ecosystem services is poorly investigated. Here, using a replicated experimental design, we test how Central European cities impact flying insects and the ecosystem service of pollination. City sites have lower insect species richness, particularly of Diptera and Lepidoptera, than neighbouring rural sites. In contrast, Hymenoptera, especially bees, show higher species richness and flower visitation rates in cities, where our experimentally derived measure of pollination is correspondingly higher. As well as revealing facets of biodiversity (e.g. phylogenetic diversity) that correlate well with pollination, we also find that ecotones in insect-friendly green cover surrounding both urban and rural sites boost pollination. Appropriately managed cities could enhance the conservation of Hymenoptera and thereby act as hotspots for pollination services that bees provide to wild flowers and crops grown in urban settings.

A single ChIP-seq dataset is sufficient for comprehensive analysis of motifs co-occurrence with MCOT package.

Recognition of composite elements consisting of two transcription factor binding sites gets behind the studies of tissue-, stage- and condition-specific transcription. Genome-wide data on transcription factor binding generated with ChIP-seq method facilitate an identification of composite elements, but the existing bioinformatics tools either require ChIP-seq datasets for both partner transcription factors, or omit composite elements with motifs overlapping. Here we present an universal Motifs Co-Occurrence Tool (MCOT) that retrieves maximum information about overrepresented composite elements from a single ChIP-seq dataset. This includes homo- and heterotypic composite elements of four mutual orientations of motifs, separated with a spacer or overlapping, even if recognition of motifs within composite element requires various stringencies. Analysis of 52 ChIP-seq datasets for 18 human transcription factors confirmed that for over 60% of analyzed datasets and transcription factors predicted co-occurrence of motifs implied experimentally proven protein-protein interaction of respecting transcription factors. Analysis of 164 ChIP-seq datasets for 57 mammalian transcription factors showed that abundance of predicted composite elements with an overlap of motifs compared to those with a spacer more than doubled; and they had 1.5-fold increase of asymmetrical pairs of motifs with one more conservative 'leading' motif and another one 'guided'.

Causes and Consequences of A Glutamine Induced Normoxic HIF1 Activity for the Tumor Metabolism.

The transcription factor hypoxia-inducible factor 1 (HIF1) is the crucial regulator of genes that are involved in metabolism under hypoxic conditions, but information regarding the transcriptional activity of HIF1 in normoxic metabolism is limited. Different tumor cells were treated under normoxic and hypoxic conditions with various drugs that affect cellular metabolism. HIF1α was silenced by siRNA in normoxic/hypoxic tumor cells, before RNA sequencing and bioinformatics analyses were performed while using the breast cancer cell line MDA-MB-231 as a model. Differentially expressed genes were further analyzed and validated by qPCR, while the activity of the metabolites was determined by enzyme assays. Under normoxic conditions, HIF1 activity was significantly increased by (i) glutamine metabolism, which was associated with the release of ammonium, and it was decreased by (ii) acetylation via acetyl CoA synthetase (ACSS2) or ATP citrate lyase (ACLY), respectively, and (iii) the presence of L-ascorbic acid, citrate, or acetyl-CoA. Interestingly, acetylsalicylic acid, ibuprofen, L-ascorbic acid, and citrate each significantly destabilized HIF1α only under normoxia. The results from the deep sequence analyses indicated that, in HIF1-siRNA silenced MDA-MB-231 cells, 231 genes under normoxia and 1384 genes under hypoxia were transcriptionally significant deregulated in a HIF1-dependent manner. Focusing on glycolysis genes, it was confirmed that HIF1 significantly regulated six normoxic and 16 hypoxic glycolysis-associated gene transcripts. However, the results from the targeted metabolome analyses revealed that HIF1 activity affected neither the consumption of glucose nor the release of ammonium or lactate; however, it significantly inhibited the release of the amino acid alanine. This study comprehensively investigated, for the first time, how normoxic HIF1 is stabilized, and it analyzed the possible function of normoxic HIF1 in the transcriptome and metabolic processes of tumor cells in a breast cancer cell model. Furthermore, these data imply that HIF1 compensates for the metabolic outcomes of glutaminolysis and, subsequently, the Warburg effect might be a direct consequence of the altered amino acid metabolism in tumor cells.

Prediction of regulatory targets of alternative isoforms of the epidermal growth factor receptor in a glioblastoma cell line.

BACKGROUND: The epidermal growth factor receptor (EGFR) is a major regulator of proliferation in tumor cells. Elevated expression levels of EGFR are associated with prognosis and clinical outcomes of patients in a variety of tumor types. There are at least four splice variants of the mRNA encoding four protein isoforms of EGFR in humans, named I through IV. EGFR isoform I is the full-length protein, whereas isoforms II-IV are shorter protein isoforms. Nevertheless, all EGFR isoforms bind the epidermal growth factor (EGF). Although EGFR is an essential target of long-established and successful tumor therapeutics, the exact function and biomarker potential of alternative EGFR isoforms II-IV are unclear, motivating more in-depth analyses. Hence, we analyzed transcriptome data from glioblastoma cell line SF767 to predict target genes regulated by EGFR isoforms II-IV, but not by EGFR isoform I nor other receptors such as HER2, HER3, or HER4. RESULTS: We analyzed the differential expression of potential target genes in a glioblastoma cell line in two nested RNAi experimental conditions and one negative control, contrasting expression with EGF stimulation against expression without EGF stimulation. In one RNAi experiment, we selectively knocked down EGFR splice variant I, while in the other we knocked down all four EGFR splice variants, so the associated effects of EGFR II-IV knock-down can only be inferred indirectly. For this type of nested experimental design, we developed a two-step bioinformatics approach based on the Bayesian Information Criterion for predicting putative target genes of EGFR isoforms II-IV. Finally, we experimentally validated a set of six putative target genes, and we found that qPCR validations confirmed the predictions in all cases. CONCLUSIONS: By performing RNAi experiments for three poorly investigated EGFR isoforms, we were able to successfully predict 1140 putative target genes specifically regulated by EGFR isoforms II-IV using the developed Bayesian Gene Selection Criterion (BGSC) approach. This approach is easily utilizable for the analysis of data of other nested experimental designs, and we provide an implementation in R that is easily adaptable to similar data or experimental designs together with all raw datasets used in this study in the BGSC repository, https://github.com/GrosseLab/BGSC .