The transcription factor encyclopedia

FOXO3 polymorphisms influence the risk and prognosis of rhabdomyosarcoma in children

Background

Rhabdomyosarcoma(RMS) is the most common soft tissue sarcoma in children and single nucleotide polymorphisms(SNPs) in certain genes influence risk of RMS. Although FOXO3 had been reported in multiple cancers including RMS, the role of FOXO3 polymorphisms in RMS remains unclear. In this case-control study, we evaluated the association of FOXO3 SNPs with RMS risk and prognosis in children.

Methods

Four FOXO3 SNPs(rs17069665 A>G, rs4946936 T>C, rs4945816 C>T and rs9400241 C>A) were genotyped in 110 RMS cases and 359 controls. The associations between FOXO3 polymorphisms and RMS risk were determined by odds ratios(ORs) with 95% confidence intervals(CIs). The associations of rs17069665 and rs4946936 with overall survival in RMS children were estimated using the Kaplan-Meier method and log-rank test. Functional analysis in silico was performed to estimate the probability that rs17069665 and rs4946936 might influence the regulation of FOXO3.

Results

We found that rs17069665 (GG vs. AA+AG, adjusted OR=2.96; 95%CI [1.10-3.32]; P=0.010) and rs4946936 (TC+CC vs. TT, adjusted OR=0.48; 95%CI [0.25-0.90]; P=0.023) were related to the increased and decreased RMS risk, respectively. Besides, rs17069665(P<0.001) and rs4946936(P<0.001) were associated with decreased and increased overall survival in RMS patients, respectively. Functional analysis showed that rs17069665 and rs4946936 might influence the transcription and expression of FOXO3 via altering the bindings to MYC, CTCF, and/or RELA.

Conclusions

This study revealed that FOXO3 polymorphisms influence the RMS susceptibility and prognosis in children, and might altered the expression of FOXO3. FOXO3 polymorphism was suggested as a biomarker for RMS susceptibility and prognosis.

Identifying Reproducible Transcription Regulator Coexpression Patterns with Single Cell Transcriptomics

The proliferation of single cell transcriptomics has potentiated our ability to unveil patterns that reflect dynamic cellular processes, rather than cell type compositional effects that emerge from bulk tissue samples. In this study, we leverage a broad collection of single cell RNA-seq data to identify the gene partners whose expression is most coordinated with each human and mouse transcription regulator (TR). We assembled 120 human and 103 mouse scRNA-seq datasets from the literature (>28 millions cells), constructing a single cell coexpression network for each. We aimed to understand the consistency of TR coexpression profiles across a broad sampling of biological contexts, rather than examine the preservation of context-specific signals. Our workflow therefore explicitly prioritizes the patterns that are most reproducible across cell types. Towards this goal, we characterize the similarity of each TR's coexpression within and across species. We create single cell coexpression rankings for each TR, demonstrating that this aggregated information recovers literature curated targets on par with ChIP-seq data. We then combine the coexpression and ChIP-seq information to identify candidate regulatory interactions supported across methods and species. Finally, we highlight interactions for the important neural TR ASCL1 to demonstrate how our compiled information can be adopted for community use.

Genome-Wide Identification and Expression Analysis of the SBP-Box Gene Family in Loquat Fruit Development

The loquat (Eriobotrya japonica L.) is a special evergreen tree, and its fruit is of high medical and health value as well as having stable market demand around the world. In recent years, research on the accumulation of nutrients in loquat fruit, such as carotenoids, flavonoids, and terpenoids, has become a hotspot. The SBP-box gene family encodes transcription factors involved in plant growth and development. However, there has been no report on the SBP-box gene family in the loquat genome and their functions in carotenoid biosynthesis and fruit ripening. In this study, we identified 28 EjSBP genes in the loquat genome, which were unevenly distributed on 12 chromosomes. We also systematically investigated the phylogenetic relationship, collinearity, gene structure, conserved motifs, and cis-elements of EjSBP proteins. Most EjSBP genes showed high expression in the root, stem, leaf, and inflorescence, while only five EjSBP genes were highly expressed in the fruit. Gene expression analysis revealed eight differentially expressed EjSBP genes between yellow- and white-fleshed fruits, suggesting that the EjSBP genes play important roles in loquat fruit development at the breaker stage. Notably, EjSBP01 and EjSBP19 exhibited completely opposite expression patterns between white- and yellow-fleshed fruits during fruit development, and showed a close relationship with SlCnr involved in carotenoid biosynthesis and fruit ripening, indicating that these two genes may participate in the synthesis and accumulation of carotenoids in loquat fruit. In summary, this study provides comprehensive information about the SBP-box gene family in the loquat, and identified two EjSBP genes as candidates involved in carotenoid synthesis and accumulation during loquat fruit development.

Targeted splicing therapy: new strategies for colorectal cancer

RNA splicing is the process of forming mature mRNA, which is an essential phase necessary for gene expression and controls many aspects of cell proliferation, survival, and differentiation. Abnormal gene-splicing events are closely related to the development of tumors, and the generation of oncogenic isoform in splicing can promote tumor progression. As a main process of tumor-specific splicing variants, alternative splicing (AS) can promote tumor progression by increasing the production of oncogenic splicing isoforms and/or reducing the production of normal splicing isoforms. This is the focus of current research on the regulation of aberrant tumor splicing. So far, AS has been found to be associated with various aspects of tumor biology, including cell proliferation and invasion, resistance to apoptosis, and sensitivity to different chemotherapeutic drugs. This article will review the abnormal splicing events in colorectal cancer (CRC), especially the tumor-associated splicing variants arising from AS, aiming to offer an insight into CRC-targeted splicing therapy.

The APC/C E3 ligase subunit ANAPC11 mediates FOXO3 protein degradation to promote cell proliferation and lymph node metastasis in urothelial bladder cancer

Urothelial bladder cancer (UBC) is one of the most prevalent malignancies worldwide, with striking tumor heterogeneity. Elucidating the molecular mechanisms that can be exploited for the treatment of aggressive UBC is a particularly relevant goal. Protein ubiquitination is a critical post-translational modification (PTM) that mediates the degradation of target protein via the proteasome. However, the roles of aberrant protein ubiquitination in UBC development and the underlying mechanisms by which it drives tumor progression remain unclear. In this study, taking advantage of clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein (Cas) 9 technology, we identified the ubiquitin E3 ligase ANAPC11, a critical subunit of the anaphase-promoting complex/cyclosome (APC/C), as a potential oncogenic molecule in UBC cells. Our clinical analysis showed that elevated expression of ANAPC11 was significantly correlated with high T stage, positive lymph node (LN) metastasis, and poor outcomes in UBC patients. By employing a series of in vitro experiments, we demonstrated that ANAPC11 enhanced the proliferation and invasiveness of UBC cells, while knockout of ANAPC11 inhibited the growth and LN metastasis of UBC cells in vivo. By conducting immunoprecipitation coupled with mass spectrometry, we confirmed that ANAPC11 increased the ubiquitination level of the Forkhead transcription factor FOXO3. The resulting decrease in FOXO3 protein stability led to the downregulation of the cell cycle regulator p21 and decreased expression of GULP1, a downstream effector of androgen receptor signaling. Taken together, these findings indicated that ANAPC11 plays an oncogenic role in UBC by modulating FOXO3 protein degradation. The ANAPC11-FOXO3 regulatory axis might serve as a novel therapeutic target for UBC.

Identification, characterization, and prognosis investigation of pivotal genes shared in different stages of breast cancer

One of the leading causes of death (20.1 per 100,000 women per year), breast cancer is the most prevalent cancer in females. Statistically, 95% of breast cancer are categorized as adenocarcinomas, and 55% of all patients may go into invasive phases; however, it can be successfully treated in approximately 70-80% of cases if diagnosed in the nascent stages. The emergence of breast tumor cells which are intensely resistant to conventional therapies, along with the high rate of metastasis occurrence, has highlighted the importance of finding novel strategies and treatments. One of the most advantageous schemes to alleviate this complication is to identify the common differentially expressed genes (DEGs) among primary and metastatic cancerous cells to use resultants for designing new therapeutic agents which are able to target both primary and metastatic breast tumor cells. In this study, the gene expression dataset with accession number GSE55715 was analyzed containing two primary tumor samples, three bone-metastatic samples, and three normal samples to distinguish the up- and down regulated genes in each stage compared to normal cells as control. In the next step, the common upregulated genes between the two experimental groups were detected by Venny online tool. Moreover, gene ontology, functions and pathways, gene-targeting microRNA, and influential metabolites were determined using EnrichR 2021 GO, KEGG pathways miRTarbase 2017, and HMDB 2021, respectively. Furthermore, elicited from STRING protein-protein interaction networks were imported to Cytoscape software to identify the hub genes. Then, identified hub genes were checked to validate the study using oncological databases. The results of the present article disclosed 1263 critical common DEGs (573 upregulated + 690 downregulated), including 35 hub genes that can be broadly used as new targets for cancer treatment and as biomarkers for cancer detection by evaluation of expression level. Besides, this study opens a new horizon to reveal unknown aspects of cancer signaling pathways by providing raw data evoked from in silico experiments. This study's outcomes can also be widely utilized in further lab research since it contains diverse information on common DEGs of varied stages and metastases of breast cancer, their functions, structures, interactions, and associations.

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

A novel HIF1α-STIL-FOXM1 axis regulates tumor metastasis

BACKGROUND

Metastasis is the major cause of morbidity and mortality in cancer that involves in multiple steps including epithelial-mesenchymal transition (EMT) process. Centrosome is an organelle that functions as the major microtubule organizing center (MTOC), and centrosome abnormalities are commonly correlated with tumor aggressiveness. However, the conclusive mechanisms indicating specific centrosomal proteins participated in tumor progression and metastasis remain largely unknown.

METHODS

The expression levels of centriolar/centrosomal genes in various types of cancers were first examined by in silico analysis of the data derived from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and European Bioinformatics Institute (EBI) datasets. The expression of STIL (SCL/TAL1-interrupting locus) protein in clinical specimens was further assessed by Immunohistochemistry (IHC) analysis and the oncogenic roles of STIL in tumorigenesis were analyzed using in vitro and in vivo assays, including cell migration, invasion, xenograft tumor formation, and metastasis assays. The transcriptome differences between low- and high-STIL expression cells were analyzed by RNA-seq to uncover candidate genes involved in oncogenic pathways. The quantitative polymerase chain reaction (qPCR) and reporter assays were performed to confirm the results. The chromatin immunoprecipitation (ChIP)-qPCR assay was applied to demonstrate the binding of transcriptional factors to the promoter.

RESULTS

The expression of STIL shows the most significant increase in lung and various other types of cancers, and is highly associated with patients' survival rate. Depletion of STIL inhibits tumor growth and metastasis. Interestingly, excess STIL activates the EMT pathway, and subsequently enhances cancer cell migration and invasion. Importantly, we reveal an unexpected role of STIL in tumor metastasis. A subset of STIL translocate into nucleus and associate with FOXM1 (Forkhead box protein M1) to promote tumor metastasis and stemness via FOXM1-mediated downstream target genes. Furthermore, we demonstrate that hypoxia-inducible factor 1α (HIF1α) directly binds to the STIL promoter and upregulates STIL expression under hypoxic condition.

CONCLUSIONS

Our findings indicate that STIL promotes tumor metastasis through the HIF1α-STIL-FOXM1 axis, and highlight the importance of STIL as a promising therapeutic target for lung cancer treatment.

CES2 sustains HNF4α expression to promote pancreatic adenocarcinoma progression through an epoxide hydrolase-dependent regulatory loop

Objective

Intra-tumoral expression of the serine hydrolase carboxylesterase 2 (CES2) contributes to the activation of the pro-drug irinotecan in pancreatic ductal adenocarcinoma (PDAC). Given other potential roles of CES2, we assessed its regulation, downstream effects, and contribution to tumor development in PDAC.

Methods

Association between the mRNA expression of CES2 in pancreatic tumors and overall survival was assessed using The Cancer Genome Atlas. Cell viability, clonogenic, and anchorage-independent growth assays as well as an orthotopic mouse model of PDAC were used to evaluate the biological relevance of CES2 in pancreatic cancer. CES2-driven metabolic changes were determined by untargeted and targeted metabolomic analyses.

Results

Elevated tumoral CES2 mRNA expression was a statistically significant predictor of poor overall survival in PDAC patients. Knockdown of CES2 in PDAC cells reduced cell viability, clonogenic capacity, and anchorage-independent growth in vitro and attenuated tumor growth in an orthotopic mouse model of PDAC. Mechanistically, CES2 was found to promote the catabolism of phospholipids resulting in HNF4α activation through a soluble epoxide hydrolase (sEH)-dependent pathway. Targeting of CES2 via siRNA or small molecule inhibitors attenuated HNF4α protein expression and reduced gene expression of classical/progenitor markers and increased basal-like markers. Targeting of the CES2-sEH-HNF4α axis using small molecule inhibitors of CES2 or sEH reduced cell viability.

Conclusions

We establish a novel regulatory loop between CES2 and HNF4α to sustain the progenitor subtype and promote PDAC progression and highlight the potential utility of CES2 or sEH inhibitors for the treatment of PDAC as part of non-irinotecan-containing regimens.

Effects of tobacco compound 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) on the expression of epigenetically regulated genes in lung carcinogenesis

Cigarette smoking is a major cause of lung cancer. Although tobacco smoking-induced genotoxicity has been well established, there is apparent lack of abundance functional epigenetic effects reported On cigarette smoke-induced lung carcinogenesis. The aim of this study was to determine effects of intratracheal administration of 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) utilizing target gene expression DNA methylation patterns in lung tissues of mice following twice weekly for 8 weeks treatment. An unbiased approach where genomic regions was undertaken to assess early methylation changes within mouse pulmonary tissues. A methylated-CpG island recovery assay (MIRA) was performed to map the DNA methylome in lung tissues, with the position of methylated DNA determined using a Genome Analyzer (MIRA-SEQ). Alterations in epigenetic-regulated target genes were confirmed with quantitative reverse transcription-PCR, which revealed 35 differentially hypermethylated genes including Cdkn1C, Hsf4, Hnf1a, Cdx1, and Hoxa5 and 30 differentially hypomethylated genes including Ddx4, Piwi1, Mdm2, and Pce1 in NNK-exposed lung tissue compared with controls. The main pathway of these genes for mediating biological information was analyzed using the Kyoto Encyclopedia of Genes and Genomes database. Among them, Rssf1 and Mdm2 were closely associated with NNK-induced lung carcinogenesis. Taken together, our data provide valuable resources for detecting cigarette smoke-induced lung carcinogenesis.

Experiment level curation of transcriptional regulatory interactions in neurodevelopment

To facilitate the development of large-scale transcriptional regulatory networks (TRNs) that may enable in-silico analyses of disease mechanisms, a reliable catalogue of experimentally verified direct transcriptional regulatory interactions (DTRIs) is needed for training and validation. There has been a long history of using low-throughput experiments to validate single DTRIs. Therefore, we reason that a reliable set of DTRIs could be produced by curating the published literature for such evidence. In our survey of previous curation efforts, we identified the lack of details about the quantity and the types of experimental evidence to be a major gap, despite the theoretical importance of such details for the identification of bona fide DTRIs. We developed a curation protocol to inspect the published literature for support of DTRIs at the experiment level, focusing on genes important to the development of the mammalian nervous system. We sought to record three types of low-throughput experiments: Transcription factor (TF) perturbation, TF-DNA binding, and TF-reporter assays. Using this protocol, we examined a total of 1,310 papers to assemble a collection of 1,499 unique DTRIs, involving 251 TFs and 825 target genes, many of which were not reported in any other DTRI resource. The majority of DTRIs (965; 64%) were supported by two or more types of experimental evidence and 27% were supported by all three. Of the DTRIs with all three types of evidence, 170 had been tested using primary tissues or cells and 44 had been tested directly in the central nervous system. We used our resource to document research biases among reports towards a small number of well-studied TFs. To demonstrate a use case for this resource, we compared our curation to a previously published high-throughput perturbation screen and found significant enrichment of the curated targets among genes differentially expressed in the developing brain in response to Pax6 deletion. This study demonstrates a proof-of-concept for the assembly of a high resolution DTRI resource to support the development of large-scale TRNs.

The capacity to computationally reconstruct gene regulatory networks using large-scale biological data is currently limited by the absence of a high confidence set of one-to-one regulatory interactions. Given the lengthy history of using small scale experimental assays to investigate individual interactions, we reason that a reliable collection of gene regulatory interactions could be compiled by systematically inspecting the published literature. To this end, we developed a curation protocol to examine and record evidence of regulatory interactions at the individual experiment level. Focusing on the area of brain development, we applied our pipeline to 1,310 publications. We identified 3,601 individual experiments, providing detailed information about 1,499 regulatory interactions. Many of these interactions have verified activity specifically in the embryonic brain. By capturing reports of regulatory interactions at this level of detail, we equip the users with more granular information than other similar resources, enabling more informed assessments of reliability.

Identification of novel genes and pathways regulated by the orphan nuclear receptor COUP-TFII in mouse MA-10 Leydig cells†

In males, Leydig cells are the main producers of testosterone and insulin-like 3 (INSL3), two hormones essential for sex differentiation and reproductive functions. Chicken ovalbumin upstream promoter-transcription factors I (COUP-TFI/NR2F1) and COUP-TFII (NR2F2) belong to the steroid/thyroid hormone nuclear receptor superfamily of transcription factors. In the testis, COUP-TFII is expressed and plays a role in the differentiation of cells committed to give rise to fully functional steroidogenic adult Leydig cells. Steroid production has also been shown to be diminished in COUP-TFII-depleted Leydig cells, indicating an important functional role in steroidogenesis. Until now, only a handful of target genes have been identified for COUP-TFII in Leydig cells. To provide new information into the mechanism of action of COUP-TFII in Leydig cells, we performed microarray analyses of COUP-TFII-depleted MA-10 Leydig cells. We identified 262 differentially expressed genes in COUP-TFII-depleted MA-10 cells. Many of the differentially expressed genes are known to be involved in lipid biosynthesis, lipid metabolism, male gonad development, and steroidogenesis. We validated the microarray data for a subset of the modulated genes by RT-qPCR. Downregulated genes included Hsd3b1, Cyp11a1, Prlr, Shp/Nr0b2, Fdx1, Scarb1, Inha and Gsta3. Finally, analysis of the Gsta3 and Inha gene promoters showed that at least two of the downregulated genes are potentially new direct targets for COUP-TFII. These data provide new evidence that further strengthens the important nature of COUP-TFII in steroidogenesis, androgen homeostasis, cellular defense, and differentiation in mouse Leydig cells.

Transcription at a Distance in the Budding Yeast Saccharomyces cerevisiae

Proper transcriptional regulation depends on the collaboration of multiple layers of control simultaneously. Cells tightly balance cellular resources and integrate various signaling inputs to maintain homeostasis during growth, development and stressors, among other signals. Many eukaryotes, including the budding yeast Saccharomyces cerevisiae, exhibit a non-random distribution of functionally related genes throughout their genomes. This arrangement coordinates the transcription of genes that are found in clusters, and can occur over long distances. In this work, we review the current literature pertaining to gene regulation at a distance in budding yeast.

Nanoparticle-based fluorescence probe for detection of NF-κB transcription factor in single cell via steric hindrance

A novel nanoparticle-based fluorescence probe was developed for NF-κB transcription factor detection and in situ imaging via steric hindrance. The probe contains gold nanoparticles (AuNPs) to quench fluorescence, and nucleic acids immobilized on the surface of AuNPs to output fluorescence. In the basal state, Cy5 labeled DNA1 folds its long chain into a hairpin structure and quenches fluorescence by forcing the Cy5 fluorophore close to the surface of AuNPs. After the probe enters the cell, the NF-κB transcription factor can bind to the κB site in the DNA duplex of the nucleic acids. The steric hindrance caused by NF-κB leads to the extension of the long chain of DNA1 and the removal of the Cy5 fluorophore from the surface of AuNPs, thereby restoring the fluorescence of the probe. By measuring NF-κB in cell lysis in vitro, the probe obtains a detection limit of 0.38 nM and the linear range from 0.5 to 16 nM. Repeated measurements showed the recovery in the cell nuclear extract was between 93.38 and 109.32%, with relative standard deviation less than 5%. By monitoring the sub-localization of the Cy5 fluorophore in single cell, the probe system can effectively distinguish active NF-κB (nucleus) and inactive NF-κB (cytoplasm) through in situ imaging. The well-designed probe will make up for the shortcomings of the existing technology, and reveal the regulatory role of transcription factors in many disease processes.

Dynamics of the transcriptional landscape during human fetal testis and ovary development

STUDY QUESTION

Which transcriptional program triggers sex differentiation in bipotential gonads and downstream cellular events governing fetal testis and ovary development in humans?

SUMMARY ANSWER

The characterization of a dynamically regulated protein-coding and non-coding transcriptional landscape in developing human gonads of both sexes highlights a large number of potential key regulators that show an early sexually dimorphic expression pattern.

WHAT IS KNOWN ALREADY

Gonadal sex differentiation is orchestrated by a sexually dimorphic gene expression program in XX and XY developing fetal gonads. A comprehensive characterization of its non-coding counterpart offers promising perspectives for deciphering the molecular events underpinning gonad development and for a complete understanding of the etiology of disorders of sex development in humans.

STUDY DESIGN, SIZE, DURATION

To further investigate the protein-coding and non-coding transcriptional landscape during gonad differentiation, we used RNA-sequencing (RNA-seq) and characterized the RNA content of human fetal testis (N = 24) and ovaries (N = 24) from 6 to 17 postconceptional week (PCW), a key period in sex determination and gonad development.

PARTICIPANTS/MATERIALS, SETTING, METHODS

First trimester fetuses (6-12 PCW) and second trimester fetuses (13-14 and 17 PCW) were obtained from legally induced normally progressing terminations of pregnancy. Total RNA was extracted from whole human fetal gonads and sequenced as paired-end 2 × 50 base reads. Resulting sequences were mapped to the human genome, allowing for the assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

This RNA-seq analysis of human fetal testes and ovaries at seven key developmental stages led to the reconstruction of 22 080 transcripts differentially expressed during testicular and/or ovarian development. In addition to 8935 transcripts displaying sex-independent differential expression during gonad development, the comparison of testes and ovaries enabled the discrimination of 13 145 transcripts that show a sexually dimorphic expression profile. The latter include 1479 transcripts differentially expressed as early as 6 PCW, including 39 transcription factors, 40 long non-coding RNAs and 20 novel genes. Despite the use of stringent filtration criteria (expression cut-off of at least 1 fragment per kilobase of exon model per million reads mapped, fold change of at least 2 and false discovery rate adjusted P values of less than <1%), the possibility of assembly artifacts and of false-positive differentially expressed transcripts cannot be fully ruled out.

LARGE-SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI GEO under accession number GSE116278.

LIMITATIONS, REASONS FOR CAUTION

The intrinsic nature of this bulk analysis, i.e. the sequencing of transcripts from whole gonads, does not allow direct identification of the cellular origin(s) of the transcripts characterized. Potential cellular dilution effects (e.g. as a result of distinct proliferation rates in XX and XY gonads) may account for a few of the expression profiles identified as being sexually dimorphic. Finally, transcriptome alterations that would result from exposure to pre-abortive drugs cannot be completely excluded. Although we demonstrated the high quality of the sorted cell populations used for experimental validations using quantitative RT-PCR, it cannot be totally excluded that some germline expression may correspond to cell contamination by, for example, macrophages.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of 1000 protein-coding and non-coding candidate genes showing an early, sexually dimorphic, expression pattern that have not previously been associated with sex differentiation. Collectively, these results increase our understanding of gonad development in humans, and contribute significantly to the identification of new candidate genes involved in fetal gonad differentiation. The results also provide a unique resource that may improve our understanding of the fetal origin of testicular and ovarian dysgenesis syndromes, including cryptorchidism and testicular cancers.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the French National Institute of Health and Medical Research (Inserm), the University of Rennes 1, the French School of Public Health (EHESP), the Swiss National Science Foundation [SNF n° CRS115_171007 to B.J.], the French National Research Agency [ANR n° 16-CE14-0017-02 and n° 18-CE14-0038-02 to F.C.], the Medical Research Council [MR/L010011/1 to P.A.F.] and the European Community's Seventh Framework Programme (FP7/2007-2013) [under grant agreement no 212885 to P.A.F.] and from the European Union's Horizon 2020 Research and Innovation Programme [under grant agreement no 825100 to P.A.F. and S.M.G.]. There are no competing interests related to this study.

Transcriptomic dynamics changes related to anthocyanin accumulation in the fleshy roots of carmine radish (Raphanus sativus L.) characterized using RNA-Seq

Carmine radish is famous for containing a natural red pigment (red radish pigment). However, the expression of anthocyanin biosynthesis-related genes during the dynamic development stages of the fleshy roots in carmine radish has not been fully investigated. Here, based on HPLC quantification of anthocyanin levels from our previous study, young fleshy roots of the carmine radish “Hongxin 1” obtained at the dynamic development stages of fleshy roots (seedling stage (SS), initial expansion (IE), full expansion (FE), bolting stage (BS), initial flowering stage (IFS), full bloom stage (FBS) and podding stage (PS)) were used for RNA-Seq. Approximately 126 comodulated DEGs related to anthocyanin biosynthesis (common DEGs in the dynamic growth stages of fleshy roots in carmine radish) were identified, from which most DEGs appeared to be likely to participate in anthocyanin biosynthesis, including two transcription factors, RsMYB and RsRZFP. In addition, some related proteins, e.g., RsCHS, RsDFR, RsANS, RsF′3H, RsF3GGT1, Rs3AT1, RsGSTF12, RsUFGT78D2 and RsUDGT-75C1, were found as candidate contributors to the regulatory mechanism of anthocyanin synthesis in the fleshy roots of carmine radish. In addition, 11 putative DEGs related to anthocyanin synthesis were evaluated by qRT-PCR via the (2-ΔΔCT) method; the Pearson correlation analysis indicated excellent concordance between the RNA-Seq and qRT-PCR results. Furthermore, GO enrichment analysis showed that “anthocyanin-containing compound biosynthetic process” and “anthocyanin-containing compound metabolic process” were commonly overrepresented in the dynamic growth stages of fleshy roots after the initial expansion stage. Moreover, five significantly enriched pathways were identified among the DEGs in the dynamic growth stages of fleshy roots in carmine radish, namely, flavonoid biosynthesis, flavone and flavonol biosynthesis, diterpenoid biosynthesis, anthocyanin biosynthesis, and benzoxazinoid biosynthesis. In conclusion, these results will expand our understanding of the complex molecular mechanisms of anthocyanin biosynthesis in the fleshy roots of carmine radish and the putative candidate genes involved in this process.

FOXO transcription factor family in cancer and metastasis

Forkhead box O (FOXO) transcription factors regulate diverse biological processes, affecting development, metabolism, stem cell maintenance and longevity. They have also been increasingly recognised as tumour suppressors through their ability to regulate genes essential for cell proliferation, cell death, senescence, angiogenesis, cell migration and metastasis. Mechanistically, FOXO proteins serve as key connection points to allow diverse proliferative, nutrient and stress signals to converge and integrate with distinct gene networks to control cell fate, metabolism and cancer development. In consequence, deregulation of FOXO expression and function can promote genetic disorders, metabolic diseases, deregulated ageing and cancer. Metastasis is the process by which cancer cells spread from the primary tumour often via the bloodstream or the lymphatic system and is the major cause of cancer death. The regulation and deregulation of FOXO transcription factors occur predominantly at the post-transcriptional and post-translational levels mediated by regulatory non-coding RNAs, their interactions with other protein partners and co-factors and a combination of post-translational modifications (PTMs), including phosphorylation, acetylation, methylation and ubiquitination. This review discusses the role and regulation of FOXO proteins in tumour initiation and progression, with a particular emphasis on cancer metastasis. An understanding of how signalling networks integrate with the FOXO transcription factors to modulate their developmental, metabolic and tumour-suppressive functions in normal tissues and in cancer will offer a new perspective on tumorigenesis and metastasis, and open up therapeutic opportunities for malignant diseases.

PAX8, an Emerging Player in Ovarian Cancer

Ovarian Cancer is one of the most lethal and widespread gynecological malignancies. It is the seventh leading cause of all cancer deaths worldwide. High-Grade Serous Cancer (HGSC), the most commonly occurring subtype, alone contributes to 70% of all ovarian cancer deaths. This is mainly attributed to the complete lack of symptoms during the early stages of the disease and absence of an early diagnostic marker.PAX8 is emerging as an important histological marker for most of the epithelial ovarian cancers, as it is expressed in about 90% of malignant ovarian cancers, specifically in HGSC. PAX8 is a member of the Paired-Box gene family (PAX1-9) of transcription factors whose expression is tightly controlled temporally and spatially. The PAX genes are well known for their role in embryonic development and their expression continues to persist in some adult tissues. PAX8 is required for the normal development of Müllerian duct that includes Fallopian tube, uterus, cervix, and upper part of vagina. In adults, it is expressed in the Fallopian tube and uterine epithelium and not in the ovarian epithelium. Considering the recent studies that predict the events preceding the tumorigenesis of HGSC from the Fallopian tube, PAX8 appears to have an important role in the development of ovarian cancer.In this chapter, we review some of the published findings to highlight the significance of PAX8 as an important marker and an emerging player in the pathogenesis of ovarian cancer. We also discuss regarding the future perspectives of PAX8 wherein it could contribute to the betterment of ovarian cancer diagnosis and treatment.

Hepatocyte nuclear factor 4α and cancer-related cell signaling pathways: a promising insight into cancer treatment

Hepatocyte nuclear factor 4α (HNF4α), a member of the nuclear receptor superfamily, is described as a protein that binds to the promoters of specific genes. It controls the expression of functional genes and is also involved in the regulation of numerous cellular processes. A large number of studies have demonstrated that HNF4α is involved in many human malignancies. Abnormal expression of HNF4α is emerging as a critical factor in cancer cell proliferation, apoptosis, invasion, dedifferentiation, and metastasis. In this review, we present emerging insights into the roles of HNF4α in the occurrence, progression, and treatment of cancer; reveal various mechanisms of HNF4α in cancer (e.g., the Wnt/β-catenin, nuclear factor-κB, signal transducer and activator of transcription 3, and transforming growth factor β signaling pathways); and highlight potential clinical uses of HNF4α as a biomarker and therapeutic target for cancer.

Knowledge extraction for assisted curation of summaries of bacterial transcription factor properties

Transcription factors (TFs) play a main role in transcriptional regulation of bacteria, as they regulate transcription of the genetic information encoded in DNA. Thus, the curation of the properties of these regulatory proteins is essential for a better understanding of transcriptional regulation. However, traditional manual curation of article collections to compile descriptions of TF properties takes significant time and effort due to the overwhelming amount of biomedical literature, which increases every day. The development of automatic approaches for knowledge extraction to assist curation is therefore critical. Here, we show an effective approach for knowledge extraction to assist curation of summaries describing bacterial TF properties based on an automatic text summarization strategy. We were able to recover automatically a median 77% of the knowledge contained in manual summaries describing properties of 177 TFs of Escherichia coli K-12 by processing 5961 scientific articles. For 71% of the TFs, our approach extracted new knowledge that can be used to expand manual descriptions. Furthermore, as we trained our predictive model with manual summaries of E. coli, we also generated summaries for 185 TFs of Salmonella enterica serovar Typhimurium from 3498 articles. According to the manual curation of 10 of these Salmonella typhimurium summaries, 96% of their sentences contained relevant knowledge. Our results demonstrate the feasibility to assist manual curation to expand manual summaries with new knowledge automatically extracted and to create new summaries of bacteria for which these curation efforts do not exist.

Database URL: The automatic summaries of the TFs of E. coli and Salmonella and the automatic summarizer are available in GitHub (https://github.com/laigen-unam/tf-properties-summarizer.git).

PAX6 mutation alters circadian rhythm and β cell function in mice without affecting glucose tolerance

The transcription factor PAX6 is involved in the development of the eye and pancreatic islets, besides being associated with sleep–wake cycles. Here, we investigated a point mutation in the RED subdomain of PAX6, previously described in a human patient, to present a comprehensive study of a homozygous Pax6 mutation in the context of adult mammalian metabolism and circadian rhythm. Pax6^Leca2 mice lack appropriate retinal structures for light perception and do not display normal daily rhythmic changes in energy metabolism. Despite β cell dysfunction and decreased insulin secretion, mutant mice have normal glucose tolerance. This is associated with reduced hepatic glucose production possibly due to altered circadian variation in expression of clock and metabolic genes, thereby evading hyperglycemia. Hence, our findings show that while the RED subdomain is important for β cell functional maturity, the Leca2 mutation impacts peripheral metabolism via loss of circadian rhythm, thus revealing pleiotropic effects of PAX6.

Nirav Chhabra et al. characterize adult mice carrying a homozygous mutation in Pax6 that was identified in a patient with foveal hypoplasia. They find that the Pax6 point mutation has pleiotropic effects, including defects in the mouse retinal structures, loss of the optic nerve, changes in energy metabolism and circadian rhythms, and dysregulation of genes expressed in the pancreas.

FOXO3 gene polymorphisms influence the risk of acute lymphoblastic leukemia in Chinese children

Acute lymphoblastic leukemia (ALL) is the most frequently diagnosed cancer in children and single-nucleotide polymorphisms (SNPs) in certain genes influence risk of ALL. Although FOXO3 had been demonstrated to be involved leukemia, the role of FOXO3 polymorphisms was still not clear. In the present study, we explored the association of FOXO3 SNPs with ALL risk in Chinese children. We genotyped four polymorphisms (rs17069665 A>G, rs4945816 T>C, rs4946936 C>T, and rs9400241 A>C) of FOXO3 in 425 ALL cases and 1339 health controls. The associations were estimated by odds ratios (ORs) with their 95% confidence intervals (CIs). Further analyses were performed to explore associations of rs17069665 and rs9400241 with ALL susceptibility in terms of age, gender, immunophenotype, minimal residual disease (MRD), and other clinical characteristics. We found rs17069665 related to the increased ALL risk (OR = 1.76; 95% CI = 1.02-3.04), rs9400241 related to decreased ALL risk (OR = 0.80; 95% CI = 0.64-0.99). The effects of rs17069665 on ALL risk were more predominant in males and children < 10 years, and patients with lower rates of platelet or neutrophil. As for rs9400241, the effects were more predominant in children < 10 years, and in patients with pre B ALL, positive MRD, anemia, or hepatomegaly. In conclusion, FOXO3 gene polymorphisms influence the risk of ALL in children and might be a potential biomarker for ALL susceptibility.

Genome-wide analysis of transcription factors related to anthocyanin biosynthesis in carmine radish (Raphanus sativus L.) fleshy roots

Carmine radish produced in Chongqing is famous for containing a natural red pigment (red radish pigment). However, the anthocyanin biosynthesis transcriptome and the expression of anthocyanin biosynthesis-related genes in carmine radish have not been fully investigated. Uncovering the mechanism of anthocyanin biosynthesis in the 'Hongxin 1' carmine radish cultivar has become a dominant research topic in this field. In this study, a local carmine radish cultivar named 'Hongxin 1' containing a highly natural red pigment was used to analyze transcription factors (TFs) related to anthocyanin biosynthesis during the dynamic development of fleshy roots. Based on RNA sequencing data, a total of 1,747 TFs in 64 TF families were identified according to their DNA-binding domains. Of those, approximately 71 differentially expressed transcription factors (DETFs) were commonly detected in any one stage compared with roots in the seedling stage (SS_root). Moreover, 26 transcripts of DETFs targeted by 74 miRNAs belonging to 25 miRNA families were identified, including MYB, WRKY, bHLH, ERF, GRAS, NF-YA, C2H2-Dof, and HD-ZIP. Finally, eight DETF transcripts belonging to the C2C2-Dof, bHLH and ERF families and their eight corresponding miRNAs were selected for qRT-PCR to verify their functions related to anthocyanin biosynthesis during the development of carmine radish fleshy roots. Finally, we propose a putative miRNA-target regulatory model associated with anthocyanin biosynthesis in carmine radish. Our findings suggest that sucrose synthase might act as an important regulator to modulate anthocyanin biosynthesis in carmine radish by inducing several miRNAs (miR165a-5p, miR172b, miR827a, miR166g and miR1432-5p) targeting different ERFs than candidate miRNAs in the traditional WMBW complex in biological processes.

Nuclear receptor HNF4α performs a tumor suppressor function in prostate cancer via its induction of p21-driven cellular senescence

Hepatocyte nuclear factor 4α (HNF4α, NR2A1) is a highly conserved member of the nuclear receptor superfamily. Recent advances reveal that it is a key transcriptional regulator of genes, broadly involved in xenobiotic and drug metabolism and also cancers of gastrointestinal tract. However, the exact functional roles of HNF4α in prostate cancer progression are still not fully understood. In this study, we determined the functional significance of HNF4α in prostate cancer. Our results showed that HNF4α exhibited a reduced expression pattern in clinical prostate cancer tissues, prostate cancer cell lines and xenograft model of castration-relapse prostate cancer. Stable HNF4α knockdown not only could promote cell proliferation and suppress doxorubicin (Dox)-induced cellular senescence in prostate cancer cells, but also confer resistance to paclitaxel treatment and enhance colony formation capacity and in vivo tumorigenicity of prostate cancer cells. On the contrary, ectopic overexpression of HNF4α could significantly inhibit the cell proliferation of prostate cancer cells, induce cell-cycle arrest at G₂/M phase and trigger the cellular senescence in prostate cancer cells by activation of p21 signal pathway in a p53-independent manner via its direct transactivation of CDKN1A. Together, our results show that HNF4α performs a tumor suppressor function in prostate cancer via a mechanism of p21-driven cellular senescence.

HNF4α-Deficient Fatty Liver Provides a Permissive Environment for Sex-Independent Hepatocellular Carcinoma

The incidence of hepatocellular carcinoma (HCC) is on the rise worldwide. Although the incidence of HCC in males is considerably higher than in females, the projected rates of HCC incidence are increasing for both sexes. A recently appreciated risk factor for HCC is the growing problem of nonalcoholic fatty liver disease, which is usually associated with obesity and the metabolic syndrome. In this study, we showed that under conditions of fatty liver, female mice were more likely to develop HCC than expected from previous models. Using an inducible knockout model of the tumor-suppressive isoform of hepatocyte nuclear factor 4 alpha ("P1-HNF4α") in the liver in combination with prolonged high fat (HF) diet, we found that HCC developed equally in male and female mice as early as 38 weeks of age. Similar sex-independent HCC occurred in the "STAM" model of mice, in which severe hyperglycemia and HF feeding results in rapid hepatic lipid deposition, fibrosis, and ultimately HCC. In both sexes, reduced P1-HNF4α activity, which also occurs under chronic HF diet feeding, increased hepatic lipid deposition and produced a greatly augmented circadian rhythm in IL6, a factor previously linked with higher HCC incidence in males. Loss of HNF4α combined with HF feeding induced epithelial-mesenchymal transition in an IL6-dependent manner. Collectively, these data provide a mechanism-based working hypothesis that could explain the rising incidence of aggressive HCC. SIGNIFICANCE: This study provides a mechanism for the growing incidence of hepatocellular carcinoma in both men and women, which is linked to nonalcoholic fatty liver disease.

HOCOMOCO: towards a complete collection of transcription factor binding models for human and mouse via large-scale ChIP-Seq analysis

We present a major update of the HOCOMOCO collection that consists of patterns describing DNA binding specificities for human and mouse transcription factors. In this release, we profited from a nearly doubled volume of published in vivo experiments on transcription factor (TF) binding to expand the repertoire of binding models, replace low-quality models previously based on in vitro data only and cover more than a hundred TFs with previously unknown binding specificities. This was achieved by systematic motif discovery from more than five thousand ChIP-Seq experiments uniformly processed within the BioUML framework with several ChIP-Seq peak calling tools and aggregated in the GTRD database. HOCOMOCO v11 contains binding models for 453 mouse and 680 human transcription factors and includes 1302 mononucleotide and 576 dinucleotide position weight matrices, which describe primary binding preferences of each transcription factor and reliable alternative binding specificities. An interactive interface and bulk downloads are available on the web: http://hocomoco.autosome.ru and http://www.cbrc.kaust.edu.sa/hocomoco11. In this release, we complement HOCOMOCO by MoLoTool (Motif Location Toolbox, http://molotool.autosome.ru) that applies HOCOMOCO models for visualization of binding sites in short DNA sequences.

Benchmark and integration of resources for the estimation of human transcription factor activities

The prediction of transcription factor (TF) activities from the gene expression of their targets (i.e., TF regulon) is becoming a widely used approach to characterize the functional status of transcriptional regulatory circuits. Several strategies and data sets have been proposed to link the target genes likely regulated by a TF, each one providing a different level of evidence. The most established ones are (1) manually curated repositories, (2) interactions derived from ChIP-seq binding data, (3) in silico prediction of TF binding on gene promoters, and (4) reverse-engineered regulons from large gene expression data sets. However, it is not known how these different sources of regulons affect the TF activity estimations and, thereby, downstream analysis and interpretation. Here we compared the accuracy and biases of these strategies to define human TF regulons by means of their ability to predict changes in TF activities in three reference benchmark data sets. We assembled a collection of TF-target interactions for 1541 human TFs and evaluated how different molecular and regulatory properties of the TFs, such as the DNA-binding domain, specificities, or mode of interaction with the chromatin, affect the predictions of TF activity. We assessed their coverage and found little overlap on the regulons derived from each strategy and better performance by literature-curated information followed by ChIP-seq data. We provide an integrated resource of all TF-target interactions derived through these strategies, with confidence scores, as a resource for enhanced prediction of TF activities.

PropaNet: Time-Varying Condition-Specific Transcriptional Network Construction by Network Propagation

Transcription factor (TF) has a significant influence on the state of a cell by regulating multiple down-stream genes. Thus, experimental and computational biologists have made great efforts to construct TF gene networks for regulatory interactions between TFs and their target genes. Now, an important research question is how to utilize TF networks to investigate the response of a plant to stress at the transcription control level using time-series transcriptome data. In this article, we present a new computational network, PropaNet, to investigate dynamics of TF networks from time-series transcriptome data using two state-of-the-art network analysis techniques, influence maximization and network propagation. PropaNet uses the influence maximization technique to produce a ranked list of TFs, in the order of TF that explains differentially expressed genes (DEGs) better at each time point. Then, a network propagation technique is used to select a group of TFs that explains DEGs best as a whole. For the analysis of Arabidopsis time series datasets from AtGenExpress, we used PlantRegMap as a template TF network and performed PropaNet analysis to investigate transcriptional dynamics of Arabidopsis under cold and heat stress. The time varying TF networks showed that Arabidopsis responded to cold and heat stress quite differently. For cold stress, bHLH and bZIP type TFs were the first responding TFs and the cold signal influenced histone variants, various genes involved in cell architecture, osmosis and restructuring of cells. However, the consequences of plants under heat stress were up-regulation of genes related to accelerating differentiation and starting re-differentiation. In terms of energy metabolism, plants under heat stress show elevated metabolic process and resulting in an exhausted status. We believe that PropaNet will be useful for the construction of condition-specific time-varying TF network for time-series data analysis in response to stress. PropaNet is available at http://biohealth.snu.ac.kr/software/PropaNet.

Retinoids induce antagonism between FOXO3A and FOXM1 transcription factors in human oral squamous cell carcinoma (OSCC) cells

To gain a greater understanding of oral squamous cell carcinoma (OSCC) we investigated the actions of all-trans-retinoic acid (RA; a retinoid), bexarotene (a pan-RXR agonist), and forkhead box (FOX) transcription factors in human OSCC-derived cell lines. RA and bexarotene have been shown to limit several oncogenic pathways in many cell types. FOXO proteins typically are associated with tumor suppressive activities, whereas FOXM1 acts as an oncogene when overexpressed in several cancers. RA and/or bexarotene increased the transcript levels of FOXO1, FOXO3A, and TRAIL receptors; reduced the transcript levels of FOXM1, Aurora kinase B (AURKB), and vascular endothelial growth factor A (VEGFA); and decreased the proliferation of OSCC-derived cell lines. Also, RA and/or bexarotene influenced the recruitment of FOXO3A and FOXM1 to target genes. Additionally, FOXM1 depletion reduced cell proliferation, decreased transcript levels of downstream targets of FOXM1, and increased transcript levels of TRAIL receptors. Overexpression of FOXO3A decreased proliferation and increased binding of histone deacetylases (HDACs) 1 and 2 at the FOXM1, AURKB, and VEGFA promoters. This research suggests novel influences of the drugs RA and bexarotene on the expression of FOXM1 and FOXO3A in transcriptional regulatory pathways of human OSCC.

Transcriptomic profile of two developmental stages of the cestode parasite Mesocestoides corti

Cestode development involves complex morphological and physiological changes. Here, we performed a differential expression analysis of gene transcripts between two developmental stages of the model cestode Mesocestoides corti. A RNA-seq-based approach was used to compare the transcriptomes of the tetrathyridium (larval, TT) and strobilated worm (ST) stages of the parasite. We found 19,053 transcripts, from which ∼45% were complete matches to genes previously annotated in the available M. corti draft genome sequence, ∼24% were considered novel isoforms, and ∼24% were considered potential novel transcripts. Stage-specific transcripts were found for both TTs (66) and STs (136), along with shared transcripts significantly overrepresented in one stage (342 in TTs, and 559 in STs). Differential expression and Gene Ontology term enrichment analyzes provided evidence of upregulation of different sets of transcripts associated with 'cytoskeleton', 'metabolism' and 'oxidation-reduction' processes in each stage, suggesting functional involvement of the corresponding genes with stage-specific features. Transcripts and processes enriched in the TT reflect typical larval processes that occur with the parasite in the intermediate host, such as asexual reproduction and budding, as well as active migration from the peritoneum to the liver and vice versa. In STs, transcripts associated with 'development', 'cell growth', and 'morphogenesis' were enriched, along with processes related to sexual reproduction, represented by the upregulation of numerous transcription factors, protein kinases, and histones. Overall, our results contributed to significantly increase the knowledge on the M. corti gene repertoire and expression profile in two developmental stages. Functional implications for the biology of larval and adult cestode parasites and for host-parasite interactions are discussed.

Deciphering the Molecular Mechanisms Sustaining the Estrogenic Activity of the Two Major Dietary Compounds Zearalenone and Apigenin in ER-Positive Breast Cancer Cell Lines

The flavone apigenin and the mycotoxin zearalenone are two major compounds found in the human diet which bind estrogen receptors (ERs), and therefore influence ER activity. However, the underlying mechanisms are not well known. To unravel the molecular mechanisms that could explain the differential effect of zearalenone and apigenin on ER-positive breast cancer cell proliferation, gene-reporter assays, chromatin immunoprecipitation (ChIP) experiments, proliferation assays and transcriptomic analysis were performed. We found that zearalenone and apigenin transactivated ERs and promoted the expression of estradiol (E2)-responsive genes. However, zearalenone clearly enhanced cellular proliferation, while apigenin appeared to be antiestrogenic in the presence of E2 in both ER-positive breast cancer cell lines, MCF-7 and T47D. The transcriptomic analysis showed that both compounds regulate gene expression in the same way, but with differences in intensity. Two major sets of genes were identified; one set was linked to the cell cycle and the other set was linked to stress response and growth arrest. Our results show that the transcription dynamics in gene regulation induced by apigenin were somehow different with zearalenone and E2 and may explain the differential effect of these compounds on the phenotype of the breast cancer cell. Together, our results confirmed the potential health benefit effect of apigenin, while zearalenone appeared to be a true endocrine-disrupting compound.

The forkhead-box family of transcription factors: key molecular players in colorectal cancer pathogenesis

Colorectal cancer (CRC) is the third most commonly occurring cancer worldwide and the fourth most frequent cause of death having an oncological origin. It has been found that transcription factors (TF) dysregulation, leading to the significant expression modifications of genes, is a widely distributed phenomenon regarding human malignant neoplasias. These changes are key determinants regarding tumour’s behaviour as they contribute to cell differentiation/proliferation, migration and metastasis, as well as resistance to chemotherapeutic agents. The forkhead box (FOX) transcription factor family consists of an evolutionarily conserved group of transcriptional regulators engaged in numerous functions during development and adult life. Their dysfunction has been associated with human diseases. Several FOX gene subgroup transcriptional disturbances, affecting numerous complex molecular cascades, have been linked to a wide range of cancer types highlighting their potential usefulness as molecular biomarkers. At least 14 FOX subgroups have been related to CRC pathogenesis, thereby underlining their role for diagnosis, prognosis and treatment purposes.

This manuscript aims to provide, for the first time, a comprehensive review of FOX genes’ roles during CRC pathogenesis. The molecular and functional characteristics of most relevant FOX molecules (FOXO, FOXM1, FOXP3) have been described within the context of CRC biology, including their usefulness regarding diagnosis and prognosis. Potential CRC therapeutics (including genome-editing approaches) involving FOX regulation have also been included. Taken together, the information provided here should enable a better understanding of FOX genes’ function in CRC pathogenesis for basic science researchers and clinicians.

Biomolecular Assemblies: Moving from Observation to Predictive Design

Biomolecular assembly is a key driving force in nearly all life processes, providing structure, information storage, and communication within cells and at the whole organism level. These assembly processes rely on precise interactions between functional groups on nucleic acids, proteins, carbohydrates, and small molecules, and can be fine-tuned to span a range of time, length, and complexity scales. Recognizing the power of these motifs, researchers have sought to emulate and engineer biomolecular assemblies in the laboratory, with goals ranging from modulating cellular function to the creation of new polymeric materials. In most cases, engineering efforts are inspired or informed by understanding the structure and properties of naturally occurring assemblies, which has in turn fueled the development of predictive models that enable computational design of novel assemblies. This Review will focus on selected examples of protein assemblies, highlighting the story arc from initial discovery of an assembly, through initial engineering attempts, toward the ultimate goal of predictive design. The aim of this Review is to highlight areas where significant progress has been made, as well as to outline remaining challenges, as solving these challenges will be the key that unlocks the full power of biomolecules for advances in technology and medicine.

Insulin signaling and reduced glucocorticoid receptor activity attenuate postprandial gene expression in liver

Hepatic circadian gene transcription is tightly coupled to feeding behavior, which has a profound impact on metabolic disorders associated with diet-induced obesity. Here, we describe a genomics approach to uncover mechanisms controlling hepatic postprandial gene expression. Combined transcriptomic and cistromic analysis identified hundreds of circadian-regulated genes and enhancers controlled by feeding. Postprandial suppression of enhancer activity was associated with reduced glucocorticoid receptor (GR) and Forkhead box O1 (FOXO1) occupancy of chromatin correlating with reduced serum corticosterone levels and increased serum insulin levels. Despite substantial co-occupancy of feeding-regulated enhancers by GR and FOXO1, selective disruption of corticosteroid and/or insulin signaling resulted in dysregulation of specific postprandial regulated gene programs. In combination, these signaling pathways operate a major part of the genes suppressed by feeding. Importantly, the feeding response was disrupted in diet-induced obese animals, which was associated with dysregulation of several corticosteroid- and insulin-regulated genes, providing mechanistic insights to dysregulated circadian gene transcription associated with obesity.

The liver is an essential organ regulating metabolic homeostasis in response to fluctuations of metabolites induced by daily rhythms of food intake. Homeostasis is maintained by precise dynamic regulation of signaling pathways controlling a wealth of enzymatic reactions involving lipid, bile acid, amino acid and glucose synthesis, storage, and redistribution in hepatocytes. Precise temporal expression of hepatic enzymes is crucial for metabolic homeostasis; a major part of circadian hepatic protein expression is regulated by precisely timed gene transcription. Here, we use a genomics approach to identify genes and regulatory regions of the genome involved in feeding-regulated gene expression. We find that transcription factors acting downstream of glucocorticoid and insulin signaling are enriched at regulatory regions repressed by feeding. Importantly, insulin and glucocorticoid signaling operate to cooperatively control the majority of feeding-mediated gene repression, and these signaling pathways are dysregulated in diet-induced obesity impacting dynamic hepatic gene expression.

Aberrant expression of alternative isoforms of transcription factors in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most prevalent malignancies worldwide and the second leading cause of death among all cancer types. Deregulation of the networks of tissue-specific transcription factors (TFs) observed in HCC leads to profound changes in the hepatic transcriptional program that facilitates tumor progression. In addition, recent reports suggest that substantial aberrations in the production of TF isoforms occur in HCC. In vitro experiments have identified distinct isoform-specific regulatory functions and related biological effects of liver-specific TFs that are implicated in carcinogenesis, which may be relevant for tumor progression and clinical outcome. This study reviews available data on the expression of isoforms of liver-specific and ubiquitous TFs in the liver and HCC and their effects, including HNF4α, C/EBPs, p73 and TCF7L2, and indicates that assessment of the ratio of isoforms and targeting specific TF variants may be beneficial for the prognosis and treatment of HCC.

A hnRNP K⁻AR-Related Signature Reflects Progression toward Castration-Resistant Prostate Cancer

The major challenge in castration-resistant prostate cancer (CRPC) remains the ability to predict the clinical responses to improve patient selection for appropriate treatments. The finding that androgen deprivation therapy (ADT) induces alterations in the androgen receptor (AR) transcriptional program by AR coregulators activity in a context-dependent manner, offers the opportunity for identifying signatures discriminating different clinical states of prostate cancer (PCa) progression. Gel electrophoretic analyses combined with western blot showed that, in androgen-dependent PCa and CRPC in vitro models, the subcellular distribution of spliced and serine-phosphorylated heterogeneous nuclear ribonucleoprotein K (hnRNP K) isoforms can be associated with different AR activities. Using mass spectrometry and bioinformatic analyses, we showed that the protein sets of androgen-dependent (LNCaP) and ADT-resistant cell lines (PDB and MDB) co-immunoprecipitated with hnRNP K varied depending on the cell type, unravelling a dynamic relationship between hnRNP K and AR during PCa progression to CRPC. By comparing the interactome of LNCaP, PDB, and MDB cell lines, we identified 51 proteins differentially interacting with hnRNP K, among which KLK3, SORD, SPON2, IMPDH2, ACTN4, ATP1B1, HSPB1, and KHDRBS1 were associated with AR and differentially expressed in normal and tumor human prostate tissues. This hnRNP K–AR-related signature, associated with androgen sensitivity and PCa progression, may help clinicians to better manage patients with CRPC.

The multisystemic functions of FOXD1 in development and disease

Transcription factors (TFs) participate in a wide range of cellular processes due to their inherent function as essential regulatory proteins. Their dysfunction has been linked to numerous human diseases. The forkhead box (FOX) family of TFs belongs to the "winged helix" superfamily, consisting of proteins sharing a related winged helix-turn-helix DNA-binding motif. FOX genes have been extensively present during vertebrates and invertebrates' evolution, participating in numerous molecular cascades and biological functions, such as embryonic development and organogenesis, cell cycle regulation, metabolism control, stem cell niche maintenance, signal transduction, and many others. FOXD1, a forkhead TF, has been related to different key biological processes such as kidney and retina development and embryo implantation. FOXD1 dysfunction has been linked to different pathologies, thereby constituting a diagnostic biomarker and a promising target for future therapies. This paper aims to present, for the first time, a comprehensive review of FOXD1's role in mouse development and human disease. Molecular, structural, and functional aspects of FOXD1 are presented in light of physiological and pathogenic conditions, including its role in human disease aetiology, such as cancer and recurrent pregnancy loss. Taken together, the information given here should enable a better understanding of FOXD1 function for basic science researchers and clinicians.

RNA cytosine methylation and methyltransferases mediate chromatin organization and 5-azacytidine response and resistance in leukaemia

The roles of RNA 5-methylcytosine (RNA:m⁵C) and RNA:m⁵C methyltransferases (RCMTs) in lineage-associated chromatin organization and drug response/resistance are unclear. Here we demonstrate that the RCMTs, namely NSUN3 and DNMT2, directly bind hnRNPK, a conserved RNA-binding protein. hnRNPK interacts with the lineage-determining transcription factors (TFs), GATA1 and SPI1/PU.1, and with CDK9/P-TEFb to recruit RNA-polymerase-II at nascent RNA, leading to formation of 5-Azacitidine (5-AZA)-sensitive chromatin structure. In contrast, NSUN1 binds BRD4 and RNA-polymerase-II to form an active chromatin structure that is insensitive to 5-AZA, but hypersensitive to the BRD4 inhibitor JQ1 and to the downregulation of NSUN1 by siRNAs. Both 5-AZA-resistant leukaemia cell lines and clinically 5-AZA-resistant myelodysplastic syndrome and acute myeloid leukaemia specimens have a significant increase in RNA:m⁵C and NSUN1-/BRD4-associated active chromatin. This study reveals novel RNA:m⁵C/RCMT-mediated chromatin structures that modulate 5-AZA response/resistance in leukaemia cells, and hence provides a new insight into treatment of leukaemia.

Resistance to chemotherapy is a serious issue that can be influenced by RNA epigenetics and chromatin structure. Here, the authors show in leukaemia cells that RNA 5-methylcytosine (RNA:m⁵C) and RNA:m⁵C methyltransferases (RCMTs) mediate chromatin structures that can modulate 5-Azacitidine response and resistance.

Coordination of signalling networks and tumorigenic properties by ABL in glioblastoma cells

The cytoplasmic tyrosine kinase ABL exerts positive or negative effects in solid tumours according to the cellular context, thus functioning as a “switch modulator”. The therapeutic effects of drugs targeting a set of signals encompassing ABL have been explored in several solid tumours. However, the net contribution of ABL inhibition by these agents remains elusive as these drugs also act on other signalling components. Here, using glioblastoma (GBM) as a cellular paradigm, we report that ABL inhibition exacerbates mesenchymal features as highlighted by down-regulation of epithelial markers and up-regulation of mesenchymal markers. Cells with permanent ABL inhibition exhibit enhanced motility and invasive capabilities, while proliferation and tumorigenic properties are reduced. Intriguingly, permanent ABL inhibition also interferes with GBM neurosphere formation and with expression of stemness markers in sphere-cultured GBM cells. Furthermore, we show that the molecular and biological characteristics of GBM cells with impaired ABL are reversible by restoring ABL levels, thus uncovering a remarkable plasticity of GBM cells to ABL threshold. A phospho-signalling screen revealed that loss of tumorigenic and self-renewal properties in GBM cells under permanent ABL inhibition coincide with drastic changes in the expression and/or phosphorylation levels of multiple signalling components. Our findings identify ABL as a crucial player for migration, invasion, proliferation, tumorigenic, and stem-cell like properties of GBM cells. Taken together, this work supports the notion that the oncogenic role of ABL in GBM cells is associated with its capability to coordinate a signalling setting that determines tumorigenic and stem-cell like properties.

Transcriptomic signatures of cellular and humoral immune responses in older adults after seasonal influenza vaccination identified by data-driven clustering

PBMC transcriptomes after influenza vaccination contain valuable information about factors affecting vaccine responses. However, distilling meaningful knowledge out of these complex datasets is often difficult and requires advanced data mining algorithms. We investigated the use of the data-driven Weighted Gene Correlation Network Analysis (WGCNA) gene clustering method to identify vaccine response-related genes in PBMC transcriptomic datasets collected from 138 healthy older adults (ages 50-74) before and after 2010-2011 seasonal trivalent influenza vaccination. WGCNA separated the 14,197 gene dataset into 15 gene clusters based on observed gene expression patterns across subjects. Eight clusters were strongly enriched for genes involved in specific immune cell types and processes, including B cells, T cells, monocytes, platelets, NK cells, cytotoxic T cells, and antiviral signaling. Examination of gene cluster membership identified signatures of cellular and humoral responses to seasonal influenza vaccination, as well as pre-existing cellular immunity. The results of this study illustrate the utility of this publically available analysis methodology and highlight genes previously associated with influenza vaccine responses (e.g., CAMK4, CD19), genes with functions not previously identified in vaccine responses (e.g., SPON2, MATK, CST7), and previously uncharacterized genes (e.g. CORO1C, C8orf83) likely related to influenza vaccine-induced immunity due to their expression patterns.

FOXM1 promotes invasion and migration of colorectal cancer cells partially dependent on HSPA5 transactivation

In this study, to investigate whether endoplastic reticulum (ER) stress correlated with FOXM1 in colorectal cancer, we analysed the mRNA levels of FOXM1 and ER stress markers HSPA5 and spliced XBP1 by qRT-PCR. FOXM1 mRNA levels were found to positively correlate with HSPA5 in colorectal cancer. However, no significant correlation between FOXM1 and spliced XBP1 mRNA levels was found. Theses results suggested the positive correlation between FOXM1 and HSPA5 in colorectal cancer was not associated with ER stress. Next, we provided evidences that FOXM1 promoted HSPA5 transcription by directly binding to and stimulating HSPA5 promoter. Moreover, a FOXM1-binding site mapped between -1019 and -1012 bp of the proximal HSPA5 promoter was identified. In addition, we found that enhancement of cell migration and invasion by FOXM1 was significantly attenuated by depletion of HSPA5 in colorectal cancer cell. Furthermore, FOXM1 triggered colorectal cancer cell migration and invasion was involved in activities of cell-surface HSPA5. Lastly, our results suggested FOXM1 facilitated the activities and expressions of MMP2 and 9 associated with cell-surface HSPA5 in colorectal cancer cells. Moreover, statistically significant positive correlations between FOXM1 and MMP2 mRNA expression, between HSPA5 and MMP2 were found in colorectal cancer tissue specimens. Together, our results suggested that FOXM1-HSPA5 signaling might be considered as a novel molecular target for designing novel therapeutic regimen to control colorectal cancer metastasis and progression.

Systematic target function annotation of human transcription factors

Background

Transcription factors (TFs), the key players in transcriptional regulation, have attracted great experimental attention, yet the functions of most human TFs remain poorly understood. Recent capabilities in genome-wide protein binding profiling have stimulated systematic studies of the hierarchical organization of human gene regulatory network and DNA-binding specificity of TFs, shedding light on combinatorial gene regulation. We show here that these data also enable a systematic annotation of the biological functions and functional diversity of TFs.

Result

We compiled a human gene regulatory network for 384 TFs covering the 146,096 TF–target gene (TF–TG) relationships, extracted from over 850 ChIP-seq experiments as well as the literature. By integrating this network of TF–TF and TF–TG relationships with 3715 functional concepts from six sources of gene function annotations, we obtained over 9000 confident functional annotations for 279 TFs. We observe extensive connectivity between TFs and Mendelian diseases, GWAS phenotypes, and pharmacogenetic pathways. Further, we show that TFs link apparently unrelated functions, even when the two functions do not share common genes. Finally, we analyze the pleiotropic functions of TFs and suggest that the increased number of upstream regulators contributes to the functional pleiotropy of TFs.

Conclusion

Our computational approach is complementary to focused experimental studies on TF functions, and the resulting knowledge can guide experimental design for the discovery of unknown roles of TFs in human disease and drug response.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0469-0) contains supplementary material, which is available to authorized users.

TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions

Transcription factors (TFs) are major trans-acting factors in transcriptional regulation. Therefore, elucidating TF-target interactions is a key step toward understanding the regulatory circuitry underlying complex traits such as human diseases. We previously published a reference TF-target interaction database for humans-TRRUST (Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining)-which was constructed using sentence-based text mining, followed by manual curation. Here, we present TRRUST v2 (www.grnpedia.org/trrust) with a significant improvement from the previous version, including a significantly increased size of the database consisting of 8444 regulatory interactions for 800 TFs in humans. More importantly, TRRUST v2 also contains a database for TF-target interactions in mice, including 6552 TF-target interactions for 828 mouse TFs. TRRUST v2 is also substantially more comprehensive and less biased than other TF-target interaction databases. We also improved the web interface, which now enables prioritization of key TFs for a physiological condition depicted by a set of user-input transcriptional responsive genes. With the significant expansion in the database size and inclusion of the new web tool for TF prioritization, we believe that TRRUST v2 will be a versatile database for the study of the transcriptional regulation involved in human diseases.

The FOXO3-FOXM1 axis: A key cancer drug target and a modulator of cancer drug resistance

The FOXO3 and FOXM1 forkhead box transcription factors, functioning downstream of the essential PI3K-Akt, Ras-ERK and JNK/p38MAPK signalling cascades, are crucial for cell proliferation, differentiation, cell survival, senescence, DNA damage repair and cell cycle control. The development of resistance to both conventional and newly emerged molecularly targeted therapies is a major challenge confronting current cancer treatment in the clinic. Intriguingly, the mechanisms of resistance to ‘classical’ cytotoxic chemotherapeutics and to molecularly targeted therapies are invariably linked to deregulated signalling through the FOXO3 and FOXM1 transcription factors. This is owing to the involvement of FOXO3 and FOXM1 in the regulation of genes linked to crucial drug action-related cellular processes, including stem cell renewal, DNA repair, cell survival, drug efflux, and deregulated mitosis. A better understanding of the mechanisms regulating the FOXO3-FOXM1 axis, as well as their downstream transcriptional targets and functions, may render these proteins reliable and early diagnostic/prognostic factors as well as crucial therapeutic targets for cancer treatment and importantly, for overcoming chemotherapeutic drug resistance.

Evaluation of Sox2 binding affinities for distinct DNA patterns using steered molecular dynamics simulation

Transcription factors (TFs) are gene expression regulators that bind to DNA in a sequence‐specific manner and determine the functional characteristics of the gene. It is worthwhile to study the unique characteristics of such specific TF‐binding pattern in DNA. Sox2 recognizes a 6‐ to 7‐base pair consensus DNA sequence; the central four bases of the binding site are highly conserved, whereas the two to three flanking bases are variable. Here, we attempted to analyze the binding affinity and specificity of the Sox2 protein for distinct DNA sequence patterns via steered molecular dynamics, in which a pulling force is employed to dissociate Sox2 from Sox2–DNA during simulation to study the behavior of a complex under nonequilibrium conditions. The simulation results revealed that the first two stacking bases of the binding pattern have an exclusive impact on the binding affinity, with the corresponding mutant complexes showing greater binding and longer dissociation time than the experimental complexes do. In contrast, mutation of the conserved bases tends to reduce the affinity, and mutation of the complete conserved region disrupts the binding. It might pave the way to identify the most likely binding pattern recognized by Sox2 based on the affinity of each configuration. The α2‐helix of Sox2 was found to be the key player in the Sox2–DNA association. The characterization of Sox2's binding patterns for the target genes in the genome helps in understanding of its regulatory functions.

Transcription Factors Responding to Pb Stress in Maize

Pb can damage the physiological function of human organs by entering the human body via food-chain enrichment. Revealing the mechanisms of maize tolerance to Pb is critical for preventing this. In this study, a Pb-tolerant maize inbred line, 178, was used to analyse transcription factors (TFs) expressed under Pb stress based on RNA sequencing data. A total of 464 genes expressed in control check (CK) or Pb treatment samples were annotated as TFs. Among them, 262 differentially expressed transcription factors (DETs) were identified that responded to Pb treatment. Furthermore, the DETs were classified into 4 classes according to their expression patterns, and 17, 12 and 2 DETs were significantly annotated to plant hormone signal transduction, basal transcription factors and base excision repair, respectively. Seventeen DETs were found to participate in the plant hormone signal transduction pathway, where basic leucine zippers (bZIPs) were the most significantly enriched TFs, with 12 members involved. We further obtained 5 Arabidopsis transfer DNA (T-DNA) mutants for 6 of the maize bZIPs, among which the mutants atbzip20 and atbzip47, representing ZmbZIP54 and ZmbZIP107, showed obviously inhibited growth of roots and above-ground parts, compared with wild type. Five highly Pb-tolerant and 5 highly Pb-sensitive in maize lines were subjected to DNA polymorphism and expression level analysis of ZmbZIP54 and ZmbZIP107. The results suggested that differences in bZIPs expression partially accounted for the differences in Pb-tolerance among the maize lines. Our results contribute to the understanding of the molecular regulation mechanisms of TFs in maize under Pb stress.

Cancer cell line specific co-factors modulate the FOXM1 cistrome

ChIP-seq has been commonly applied to identify genomic occupation of transcription factors (TFs) in a context-specific manner. It is generally assumed that a TF should have similar binding patterns in cells from the same or closely related tissues. Surprisingly, this assumption has not been carefully examined. To this end, we systematically compared the genomic binding of the cell cycle regulator FOXM1 in eight cell lines from seven different human tissues at binding signal, peaks and target genes levels. We found that FOXM1 binding in ER-positive breast cancer cell line MCF-7 are distinct comparing to those in not only other non-breast cell lines, but also MDA-MB-231, ER-negative breast cancer cell line. However, binding sites in MDA-MB-231 and non-breast cell lines were highly consistent. The recruitment of estrogen receptor alpha (ERα) caused the unique FOXM1 binding patterns in MCF-7. Moreover, the activity of FOXM1 in MCF-7 reflects the regulatory functions of ERα, while in MDA-MB-231 and non-breast cell lines, FOXM1 activities regulate cell proliferation. Our results suggest that tissue similarity, in some specific contexts, does not hold precedence over TF-cofactors interactions in determining transcriptional states and that the genomic binding of a TF can be dramatically affected by a particular co-factor under certain conditions.

Hypoxia-Inducible Factor 2α Mutation-Related Paragangliomas Classify as Discrete Pseudohypoxic Subcluster

Recently, activating mutations of the hypoxia-inducible factor 2α gene (HIF2A/EPAS1) have been recognized to predispose to multiple paragangliomas (PGLs) and duodenal somatostatinomas associated with polycythemia, and ocular abnormalities. Previously, mutations in the SDHA/B/C/D, SDHAF2, VHL, FH, PHD1, and PHD2 genes have been associated with HIF activation and the development of pseudohypoxic (cluster-1) PGLs. These tumors overlap in terms of tumor location, syndromic presentation, and noradrenergic phenotype to a certain extent. However, they also differ especially by clinical outcome and by presence of other tumors or abnormalities. In the present study, we aimed to establish additional molecular differences between HIF2A and non-HIF2A pseudohypoxic PGLs. RNA expression patterns of HIF2A PGLs (n = 6) from 2 patients were compared with normal adrenal medullas (n = 8) and other hereditary pseudohypoxic PGLs (VHL: n = 13, SDHB: n = 15, and SDHD: n = 14). Unsupervised hierarchical clustering showed that HIF2A PGLs made up a separate cluster from other pseudohypoxic PGLs. Significance analysis of microarray yielded 875 differentially expressed genes between HIF2A and other pseudohypoxic PGLs after normalization to adrenal medulla (false discovery rate 0.01). Prediction analysis of microarray allowed correct classification of all HIF2A samples based on as little as three genes (TRHDE, LRRC63, IGSF10; error rate: 0.02). Genes with the highest expression difference between normal medulla and HIF2A PGLs were selected for confirmatory quantitative reverse transcriptase polymerase chain reaction. In conclusion, HIF2A PGLs show a characteristic expression signature that separates them from non-HIF2A pseudohypoxic PGLs. Unexpectedly, the most significantly differentially expressed genes have not been previously described as HIF target genes.

Glyceollins trigger anti-proliferative effects through estradiol-dependent and independent pathways in breast cancer cells

Background

Estrogen receptors (ER) α and β are found in both women and men in many tissues, where they have different functions, including having roles in cell proliferation and differentiation of the reproductive tract. In addition to estradiol (E2), a natural hormone, numerous compounds are able to bind ERs and modulate their activities. Among these compounds, phytoestrogens such as isoflavones, which are found in plants, are promising therapeutics for several pathologies. Glyceollins are second metabolites of isoflavones that are mainly produced in soybean in response to an elicitor. They have potentially therapeutic actions in breast cancer by reducing the proliferation of cancer cells. However, the molecular mechanisms driving these effects remain elusive.

Methods

First, to determine the proliferative or anti-proliferative effects of glyceollins, in vivo and in vitro approaches were used. The length of epithelial duct in mammary gland as well as uterotrophy after treatment by E2 and glyceollins and their effect on proliferation of different breast cell line were assessed. Secondly, the ability of glyceollin to activate ER was assessed by luciferase assay. Finally, to unravel molecular mechanisms involved by glyceollins, transcriptomic analysis was performed on MCF-7 breast cancer cells.

Results

In this study, we show that synthetic versions of glyceollin I and II exert anti-proliferative effects in vivo in mouse mammary glands and in vitro in different ER-positive and ER-negative breast cell lines. Using transcriptomic analysis, we produce for the first time an integrated view of gene regulation in response to glyceollins and reveal that these phytochemicals act through at least two major pathways. One pathway involving FOXM1 and ERα is directly linked to proliferation. The other involves the HIF family and reveals that stress is a potential factor in the anti-proliferative effects of glyceollins due to its role in increasing the expression of REDD1, an mTORC1 inhibitor.

Conclusion

Overall, our study clearly shows that glyceollins exert anti-proliferative effects by reducing the expression of genes encoding cell cycle and mitosis-associated factors and biomarkers overexpressed in cancers and by increasing the expression of growth arrest-related genes. These results reinforce the therapeutic potential of glyceollins for breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12964-017-0182-1) contains supplementary material, which is available to authorized users.