Transcription Factor RREB1: from Target Genes towards Biological Functions

A Proteomics Approach Identifies RREB1 as a Crucial Molecular Target of Imidazo-Pyrazole Treatment in SKMEL-28 Melanoma Cells

Cutaneous melanoma is the most dangerous and deadly form of human skin malignancy. Despite its rarity, it accounts for a staggering 80% of deaths attributed to cutaneous cancers overall. Moreover, its final stages often exhibit resistance to drug treatments, resulting in unfavorable outcomes. Hence, ensuring access to novel and improved chemotherapeutic agents is imperative for patients grappling with this severe ailment. Pyrazole and its fused systems derived thereof are heteroaromatic moieties widely employed in medicinal chemistry to develop effective drugs for various therapeutic areas, including inflammation, pain, oxidation, pathogens, depression, and fever. In a previous study, we described the biochemical properties of a newly synthesized group of imidazo-pyrazole compounds. In this paper, to improve our knowledge of the pharmacological properties of these molecules, we conduct a differential proteomic analysis on a human melanoma cell line treated with one of these imidazo-pyrazole derivatives. Our results detail the changes to the SKMEL-28 cell line proteome induced by 24, 48, and 72 h of 3e imidazo-pyrazole treatment. Notably, we highlight the down-regulation of the Ras-responsive element binding protein 1 (RREB1), a member of the zinc finger transcription factors family involved in the tumorigenesis of melanoma. RREB1 is a downstream element of the MAPK pathway, and its activation is mediated by ERK1/2 through phosphorylation.

Association between MCU Gene Polymorphisms with Obesity: Findings from the All of Us Research Program

Obesity is a public health crisis, and its prevalence disproportionately affects African Americans in the United States. Dysregulation of organelle calcium homeostasis is associated with obesity. The mitochondrial calcium uniporter (MCU) complex is primarily responsible for mitochondrial calcium homeostasis. Obesity is a multifactorial disease in which genetic underpinnings such as single-nucleotide polymorphisms (SNPs) may contribute to disease progression. The objective of this study was to identify genetic variations of MCU with anthropometric measurements and obesity in the All of Us Research Program.

METHODS

We used an additive genetic model to assess the association between obesity traits (body mass index (BMI), waist and hip circumference) and selected MCU SNPs in 19,325 participants (3221 normal weight and 16,104 obese). Eleven common MCU SNPs with a minor allele frequency ≥ 5% were used for analysis.

RESULTS

We observed three MCU SNPs in self-reported Black/African American (B/AA) men, and six MCU SNPs in B/AA women associated with increased risk of obesity, whereas six MCU SNPs in White men, and nine MCU SNPs in White women were protective against obesity development.

CONCLUSIONS

This study found associations of MCU SNPs with obesity, providing evidence of a potential predictor of obesity susceptibility in B/AA adults.

Rreb1 is a key transcription factor in Sertoli cell maturation and function and spermatogenesis in mouse

Spermatogenesis is a developmental process driven by interactions between germ cells and Sertoli cells. This process depends on appropriate gene expression, which might be regulated by transcription factors. This study focused on Rreb1, a zinc finger transcription factor, and explored its function and molecular mechanisms in spermatogenesis in a mouse model. Our results showed that RREB1 was predominantly expressed in the Sertoli cells of the testis. The decreased expression of RREB1 following injection of siRNA caused impaired Sertoli cell development, which was characterized using a defective blood-testis barrier structure and decreased expression of Sertoli cell functional maturity markers; its essential trigger might be SMAD3 destabilization. The decreased expression of RREB1 in mature Sertoli cells influenced the cell structure and function, which resulted in abnormal spermatogenesis, manifested as oligoasthenoteratozoospermia, and we believe RREB1 plays this role by regulating the transcription of Fshr and Wt1. RREB1 has been reported to activate Fshr transcription, and we demonstrated that the knockdown of Rreb1 caused a reduction in follicle-stimulating hormone receptor (FSHR) in the testis, which could be the cause of the increased sperm malformation. Furthermore, we confirmed that RREB1 directly activates Wt1 promoter activity, and RREB1 downregulation induced the decreased expression of Wt1 and its downstream polarity-associated genes Par6b and E-cadherin, which caused increased germ-cell death and reduced sperm number and motility. In conclusion, RREB1 is a key transcription factor essential for Sertoli cell development and function and is required for normal spermatogenesis.

Modeling methamphetamine use disorder and relapse in animals: short- and long-term epigenetic, transcriptional., and biochemical consequences in the rat brain

Methamphetamine use disorder (MUD) is a neuropsychiatric disorder characterized by binge drug taking episodes, intervals of abstinence, and relapses to drug use even during treatment. MUD has been modeled in rodents and investigators are attempting to identify its molecular bases. Preclinical experiments have shown that different schedules of methamphetamine self-administration can cause diverse transcriptional changes in the dorsal striatum of Sprague-Dawley rats. In the present review, we present data on differentially expressed genes (DEGs) identified in the rat striatum following methamphetamine intake. These include genes involved in transcription regulation, potassium channel function, and neuroinflammation. We then use the striatal data to discuss the potential significance of the molecular changes induced by methamphetamine by reviewing concordant or discordant data from the literature. This review identified potential molecular targets for pharmacological interventions. Nevertheless, there is a need for more research on methamphetamine-induced transcriptional consequences in various brain regions. These data should provide a more detailed neuroanatomical map of methamphetamine-induced changes and should better inform therapeutic interventions against MUD.

Mechanistic insights into super-enhancer-driven genes as prognostic signatures in patients with glioblastoma

BACKGROUND

Glioblastoma (GBM) is one of the most common malignant brain tumors in adults and is characterized by high aggressiveness and rapid progression, poor treatment, high recurrence rate, and poor prognosis. Although super-enhancer (SE)-driven genes haven been recognized as prognostic markers for several cancers, whether it can be served as effective prognostic markers for patients with GBM has not been evaluated.

METHODS

We first combined histone modification data with transcriptome data to identify SE-driven genes associated with prognosis in patients with GBM. Second, we developed a SE-driven differentially expressed genes (SEDEGs) risk score prognostic model by univariate Cox analysis, KM survival analysis, multivariate Cox analysis and least absolute shrinkage and selection operator (LASSO) regression. Its reliability in predicting was verified by two external data sets. Third, through mutation analysis, immune infiltration, we explored the molecular mechanisms of prognostic genes. Next, Genomics of Drug Sensitivity in Cancer (GDSC) and the Connectivity Map (cMap) database were employed to assess different sensitivities to chemotherapeutic agents and small-molecule drug candidates between high- and low-risk patients. Finally, SEanalysis database was chosen to identify SE-driven transcription factors (TFs) regulating prognostic markers which will reveal a potential SE-driven transcriptional regulatory network.

RESULTS

First, we developed a 11-gene risk score prognostic model (NCF2, MTHFS, DUSP6, G6PC3, HOXB2, EN2, DLEU1, LBH, ZEB1-AS1, LINC01265, and AGAP2-AS1) selected from 1,154 SEDEGs, which is not only an independent prognostic factor for patients, but also can effectively predict the survival rate of patients. The model can effectively predict 1-, 2- and 3-year survival of patients and was validated in external Chinese Glioma Genome Atlas (CGGA) and Gene Expression Omnibus (GEO) datasets. Second, the risk score was positively correlated with the infiltration of regulatory T cell, CD4 memory activated T cell, activated NK cell, neutrophil, resting mast cell, M0 macrophage, and memory B cell. Third, we found that high-risk patients showed higher sensitivity than low-risk patients to both 27 chemotherapeutic agents and 4 small-molecule drug candidates which might benefit further precision therapy for GBM patients. Finally, 13 potential SE-driven TFs imply how SE regulates GBM patient's prognosis.

CONCLUSION

The SEDEG risk model not only helps to elucidate the impact of SEs on the course of GBM, but also provides a bright future for prognosis determination and choice of treatment for GBM patients.

Meta-analysis of integrated ChIP-seq and transcriptome data revealed genomic regions affected by estrogen receptor alpha in breast cancer

BACKGROUND

The largest group of patients with breast cancer are estrogen receptor-positive (ER+) type. The estrogen receptor acts as a transcription factor and triggers cell proliferation and differentiation. Hence, investigating ER-DNA interaction genomic regions can help identify genes directly regulated by ER and understand the mechanism of ER action in cancer progression.

METHODS

In the present study, we employed a workflow to do a meta-analysis of ChIP-seq data of ER+ cell lines stimulated with 10 nM and 100 nM of E2. All publicly available data sets were re-analyzed with the same platform. Then, the known and unknown batch effects were removed. Finally, the meta-analysis was performed to obtain meta-differentially bound sites in estrogen-treated MCF7 cell lines compared to vehicles (as control). Also, the meta-analysis results were compared with the results of T47D cell lines for more precision. Enrichment analyses were also employed to find the functional importance of common meta-differentially bound sites and associated genes among both cell lines.

RESULTS

Remarkably, POU5F1B, ZNF662, ZNF442, KIN, ZNF410, and SGSM2 transcription factors were recognized in the meta-analysis but not in individual studies. Enrichment of the meta-differentially bound sites resulted in the candidacy of pathways not previously reported in breast cancer. PCGF2, HNF1B, and ZBED6 transcription factors were also predicted through the enrichment analysis of associated genes. In addition, comparing the meta-analysis results of both ChIP-seq and RNA-seq data showed that many transcription factors affected by ER were up-regulated.

CONCLUSION

The meta-analysis of ChIP-seq data of estrogen-treated MCF7 cell line leads to the identification of new binding sites of ER that have not been previously reported. Also, enrichment of the meta-differentially bound sites and their associated genes revealed new terms and pathways involved in the development of breast cancer which should be examined in future in vitro and in vivo studies.

Does an apple a day keep away diseases? Evidence and mechanism of action

Apples and their products exemplify the recently reemphasized link between dietary fruit intake and the alleviation of human disease. Their consumption does indeed improve human health due to their high phytochemical content. To identify potentially relevant articles from clinical trials, some epidemiological studies and meta-analyses, and in vitro and in vivo studies (cell cultures and animal models), PubMed was searched from January 1, 2012, to May 15, 2022. This review summarized the potential effects of apple and apple products (juices, puree, pomace, dried apples, extracts rich in apple bioactives and single apple bioactives) on health. Apples and apple products have protective effects against cardiovascular diseases, cancer, as well as mild cognitive impairment and promote hair growth, healing of burn wounds, improve the oral environment, prevent niacin-induced skin flushing, promote the relief of UV-induced skin pigmentation, and improve the symptoms of atopic dermatitis as well as cedar hay fever among others. These effects are associated with various mechanisms, such as vascular endothelial protection, blood lipids lowering, anti-inflammatory, antioxidant, antiapoptotic, anti-invasion, and antimetastatic effects. Meanwhile, it has provided an important reference for the application and development of medicine, nutrition, and other fields.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Moderate-intensity endurance training improves late phase β-cell function in adults with type 2 diabetes

Physical activity is important for type 2 diabetes treatment, yet the underlying mechanisms for these beneficial effects of exercise are not fully understood. Here, we investigated the effects of exercise training on biphasic β-cell insulin secretory function, a key factor regulating blood glucose. Adults with type 2 diabetes (7F/3M, age 49 ± 5 years, BMI 30 ± 3 kg/m²) completed a 10-week moderate-intensity exercise program and multiple components of glucose homeostasis were measured. Training improved glycemic control, insulin sensitivity, and processing of proinsulin-to-insulin. Training increased late phase β-cell function by 38% (p = 0.01), which was correlated with changes in VO_2peak suggesting training response-dependent effects. Ras-Responsive Element Binding Protein 1 (RREB1) concentrations, a protein postulated to increase type 2 diabetes risk, were inversely correlated with increases in training-induced late-phase disposition index, consistent with an inhibitory role of RREB1 on insulin secretion. Moderate-intensity exercise training improves late-phase β-cell function and glycemic control in adults with type 2 diabetes.

Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions

BACKGROUND

The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci.

METHODS

We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci.

RESULTS

We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci.

CONCLUSIONS

Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

High-resolution Nanopore methylome-maps reveal random hyper-methylation at CpG-poor regions as driver of chemoresistance in leukemias

Aberrant DNA methylation at CpG dinucleotides is a cancer hallmark that is associated with the emergence of resistance to anti cancer treatment, though molecular mechanisms and biological significance remain elusive. Genome scale methylation maps by currently used methods are based on chemical modification of DNA and are best suited for analyses of methylation at CpG rich regions (CpG islands). We report the first high coverage whole-genome map in cancer using the long read nanopore technology, which allows simultaneous DNA-sequence and -methylation analyses on native DNA. We analyzed clonal epigenomic/genomic evolution in Acute Myeloid Leukemias (AMLs) at diagnosis and relapse, after chemotherapy. Long read sequencing coupled to a novel computational method allowed definition of differential methylation at unprecedented resolution, and showed that the relapse methylome is characterized by hypermethylation at both CpG islands and sparse CpGs regions. Most differentially methylated genes, however, were not differentially expressed nor enriched for chemoresistance genes. A small fraction of under-expressed and hyper-methylated genes at sparse CpGs, in the gene body, was significantly enriched in transcription factors (TFs). Remarkably, these few TFs supported large gene-regulatory networks including 50% of all differentially expressed genes in the relapsed AMLs and highly-enriched in chemoresistance genes. Notably, hypermethylated regions at sparse CpGs were poorly conserved in the relapsed AMLs, under-represented at their genomic positions and showed higher methylation entropy, as compared to CpG islands. Analyses of available datasets confirmed TF binding to their target genes and conservation of the same gene-regulatory networks in large patient cohorts. Relapsed AMLs carried few patient specific structural variants and DNA mutations, apparently not involved in drug resistance. Thus, drug resistance in AMLs can be mainly ascribed to the selection of random epigenetic alterations at sparse CpGs of a few transcription factors, which then induce reprogramming of the relapsing phenotype, independently of clonal genomic evolution.

Development and validation of a multivariable genotype-informed gestational diabetes prediction algorithm for clinical use in the Mexican population: insights into susceptibility mechanisms

INTRODUCTION

Gestational diabetes mellitus (GDM) is underdiagnosed in Mexico. Early GDM risk stratification through prediction modeling is expected to improve preventative care. We developed a GDM risk assessment model that integrates both genetic and clinical variables.

RESEARCH DESIGN AND METHODS

Data from pregnant Mexican women enrolled in the 'Cuido mi Embarazo' (CME) cohort were used for development (107 cases, 469 controls) and data from the 'Mónica Pretelini Sáenz' Maternal Perinatal Hospital (HMPMPS) cohort were used for external validation (32 cases, 199 controls). A 2-hour oral glucose tolerance test (OGTT) with 75 g glucose performed at 24-28 gestational weeks was used to diagnose GDM. A total of 114 single-nucleotide polymorphisms (SNPs) with reported predictive power were selected for evaluation. Blood samples collected during the OGTT were used for SNP analysis. The CME cohort was randomly divided into training (70% of the cohort) and testing datasets (30% of the cohort). The training dataset was divided into 10 groups, 9 to build the predictive model and 1 for validation. The model was further validated using the testing dataset and the HMPMPS cohort.

RESULTS

Nineteen attributes (14 SNPs and 5 clinical variables) were significantly associated with the outcome; 11 SNPs and 4 clinical variables were included in the GDM prediction regression model and applied to the training dataset. The algorithm was highly predictive, with an area under the curve (AUC) of 0.7507, 79% sensitivity, and 71% specificity and adequately powered to discriminate between cases and controls. On further validation, the training dataset and HMPMPS cohort had AUCs of 0.8256 and 0.8001, respectively.

CONCLUSIONS

We developed a predictive model using both genetic and clinical factors to identify Mexican women at risk of developing GDM. These findings may contribute to a greater understanding of metabolic functions that underlie elevated GDM risk and support personalized patient recommendations.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

A novel method to assess copy number variation in melanoma: Droplet digital PCR for precise quantitation of the RREB1 gene in formalin-fixed, paraffin-embedded melanocytic neoplasms, a proof-of-concept study

BACKGROUND

Melanocytic neoplasms can be challenging to diagnose. One well-established diagnostic aid is the detection of copy number variation (CNV) in a few key genetic loci using conventional methods such as fluorescence in situ hybridization (FISH) and chromosomal microarray (CMA). Droplet digital polymerase chain reaction (ddPCR) is a novel, cost-effective, rapid, and automated method to detect CNV.

METHODS

We perform the first investigation of ddPCR to assay Ras-responsive element-binding protein-1 (RREB1), the most common CNV in melanoma using formalin-fixed, paraffin-embedded (FFPE) melanocytic lesion samples; CMA data are used as the gold standard. Archival samples from 2013 to 2021 were analyzed, including 153 data points from 39 FFPE samples representing 34 patients. Benign, borderline, malignant, and metastatic melanocytic neoplasms were examined.

RESULTS

ddPCR showed a sensitivity and specificity of 93.8% and 95.7% using one reference gene, and 87.5% and 100% using a different reference gene for RREB1 gain detection.

CONCLUSIONS

Here we show that ddPCR can provide inexpensive, rapid, and robust data on the commonest copy number alteration in melanoma. Future development and validation could provide a useful ancillary tool in the diagnosis of challenging melanocytic lesions.

In silico identification and in vitro expression analysis of breast cancer-related m6A-SNPs

Research on m6A-associated SNPs (m6A-SNPs) has emerged recently due to their possible critical roles in many key biological processes. In this sense, several investigations have identified m6A-SNPs in different diseases. In order to gain a more complete understanding of the role that m6A-SNPs can play in breast cancer, we performed an in silico analysis to identify the m6A-SNPs associated with breast cancer and to evaluate their possible effects. For this purpose, we downloaded SNPs related to breast cancer and a list of m6A-SNPs from public databases in order to identify which ones appear in both. Subsequently, we assessed the identified m6A-SNPs in silico by expression quantitative trait loci (eQTL) analysis and differential gene expression analysis. We genotyped the m6A-SNPs found in the in silico analysis in 35 patients with breast cancer, and we carried out a gene expression analysis experimentally on those that showed differences. Our results identified 981 m6A-SNPs related to breast cancer. Four m6A-SNPs showed an eQTL effect and only three were in genes that presented an altered gene expression. When the three m6A-SNPs were evaluated in the tissue sample of our breast cancer patients, only the m6A-SNP rs76563149 located in ZNF354A gene presented differences in allele frequencies and a low gene expression in breast cancer tissues, especially in luminal B HER2+ subtype. Future investigations of these m6A-SNPs should expand the study in different ethnic groups and increase the sample sizes to test their association with breast cancer and elucidate their molecular function.

Ensemble machine learning identifies genetic loci associated with future worsening of disability in people with multiple sclerosis

Limited studies have been conducted to identify and validate multiple sclerosis (MS) genetic loci associated with disability progression. We aimed to identify MS genetic loci associated with worsening of disability over time, and to develop and validate ensemble genetic learning model(s) to identify people with MS (PwMS) at risk of future worsening. We examined associations of 208 previously established MS genetic loci with the risk of worsening of disability; we learned ensemble genetic decision rules and validated the predictions in an external dataset. We found 7 genetic loci (rs7731626: HR 0.92, P = 2.4 × 10^-5; rs12211604: HR 1.16, P = 3.2 × 10^-7; rs55858457: HR 0.93, P = 3.7 × 10^-7; rs10271373: HR 0.90, P = 1.1 × 10^-7; rs11256593: HR 1.13, P = 5.1 × 10^-57; rs12588969: HR = 1.10, P = 2.1 × 10^-10; rs1465697: HR 1.09, P = 1.7 × 10^-128) associated with risk worsening of disability; most of which were located near or tagged to 13 genomic regions enriched in peptide hormones and steroids biosynthesis pathways by positional and eQTL mapping. The derived ensembles produced a set of genetic decision rules that can be translated to provide additional prognostic values to existing clinical predictions, with the additional benefit of incorporating relevant genetic information into clinical decision making for PwMS. The present study extends our knowledge of MS progression genetics and provides the basis of future studies regarding the functional significance of the identified loci.

circXRCC5 foster gastric cancer growth and metastasis by the HNRNPC/circXRCC5/miR-655-3p/RREB1/UBA2 positive feedback loop

Gastric cancer (GC) is one of the most common malignancies, leading to millions of deaths each year. Here, we investigated the molecular mechanisms of GC, with a focus on circXRCC5/miR-655-3p/RREB1/UBA2 axis. circXRCC5 was identified in 62 paired cancer specimens and adjacent normal tissues by genome-wide bioinformatics analysis and verified by qRT-PCR and Sanger sequencing. Knockdown or exogenous expression of circXRCC5 was performed to validate the functional significance of circXRCC5 using both in vitro and in vivo assays, including CCK-8, colony formation, EdU incorporation, transwell system, as well as animal experiments. RNA immunoprecipitation, biotinylated RNA pull-down, ChIP, and dual-luciferase assays were employed to validate the regulatory network of circXRCC5/miR-655-3p/RREB1/UBA2. Frequently elevated circXRCC5 in GC tissues and cell lines was associated with poor prognosis of GC patients. Functionally, circXRCC5 overexpression facilitated GC cell proliferation, migration, and invasion, as well as promoted tumor growth and metastasis in vivo. Mechanistically, circXRCC5 served as a sponge of miR-655-3p to induce upregulation of RREB1. RREB1 was identified as a transcriptional activator of UBA2, thus contributing to GC tumorigenesis. Moreover, RNA binding protein (RBP) HNRNPC was proved to interact with circXRCC5 to promote circXRCC5 biogenesis. Collectively, circXRCC5 facilitates GC progression through the HNRNPC/circXRCC5/miR-655-3p/RREB1/UBA2 axis, which might bring novel therapeutic strategies for GC treatment.

Genomic insights into host and parasite interactions during intracellular infection by Toxoplasma gondii

To gain insights into the molecular interactions of an intracellular pathogen and its host cell, we studied the gene expression and chromatin states of human fibroblasts infected with the Apicomplexan parasite Toxoplasma gondii. We show a striking activation of host cell genes that regulate a number of cellular processes, some of which are protective of the host cell, others likely to be advantageous to the pathogen. The simultaneous capture of host and parasite genomic information allowed us to gain insights into the regulation of the T. gondii genome. We show how chromatin accessibility and transcriptional profiling together permit novel annotation of the parasite's genome, including more accurate mapping of known genes and the identification of new genes and cis-regulatory elements. Motif analysis reveals not only the known T. gondii AP2 transcription factor-binding site but also a previously-undiscovered candidate TATA box-containing motif at one-quarter of promoters. By inferring the transcription factor and upstream cell signaling responses involved in the host cell, we can use genomic information to gain insights into T. gondii's perturbation of host cell physiology. Our resulting model builds on previously-described human host cell signalling responses to T. gondii infection, linked to induction of specific transcription factors, some of which appear to be solely protective of the host cell, others of which appear to be co-opted by the pathogen to enhance its own survival.

Genetic polymorphisms of PKLR gene and their associations with milk production traits in Chinese Holstein cows

Our previous work had confirmed that pyruvate kinase L/R (PKLR) gene was expressed differently in different lactation periods of dairy cattle, and participated in lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, suggesting that PKLR is a candidate gene to affect milk production traits in dairy cattle. Here, we verified whether this gene has significant genetic association with milk yield and composition traits in a Chinese Holstein cow population. In total, we identified 21 single nucleotide polymorphisms (SNPs) by resequencing the entire coding region and partial flanking region of PKLR gene, in which, two SNPs were located in 5′ promoter region, two in 5′ untranslated region (UTR), three in introns, five in exons, six in 3′ UTR and three in 3′ flanking region. The single marker association analysis displayed that all SNPs were significantly associated with milk yield, fat and protein yields or protein percentage (p ≤ 0.0497). The haplotype block containing all the SNPs, predicted by Haploview, had a significant association with fat yield and protein percentage (p ≤ 0.0145). Further, four SNPs in 5′ regulatory region and eight SNPs in UTR and exon regions were predicted to change the transcription factor binding sites (TFBSs) and mRNA secondary structure, respectively, thus affecting the expression of PKLR, leading to changes in milk production phenotypes, suggesting that these SNPs might be the potential functional mutations for milk production traits in dairy cattle. In conclusion, we demonstrated that PKLR had significant genetic effects on milk production traits, and the SNPs with significant genetic effects could be used as candidate genetic markers for genomic selection (GS) in dairy cattle.

Anti-cancer therapeutic strategies based on HGF/MET, EpCAM, and tumor-stromal cross talk

As an intelligent disease, tumors apply several pathways to evade the immune system. It can use alternative routes to bypass intracellular signaling pathways, such as nuclear factor-κB (NF-κB), Wnt, and mitogen-activated protein (MAP)/phosphoinositide 3-kinase (PI3K)/mammalian target of rapamycin (mTOR). Therefore, these mechanisms lead to therapeutic resistance in cancer. Also, these pathways play important roles in the proliferation, survival, migration, and invasion of cells. In most cancers, these signaling pathways are overactivated, caused by mutation, overexpression, etc. Since numerous molecules share these signaling pathways, the identification of key molecules is crucial to achieve favorable consequences in cancer therapy. One of the key molecules is the mesenchymal-epithelial transition factor (MET; c-Met) and its ligand hepatocyte growth factor (HGF). Another molecule is the epithelial cell adhesion molecule (EpCAM), which its binding is hemophilic. Although both of them are involved in many physiologic processes (especially in embryonic stages), in some cancers, they are overexpressed on epithelial cells. Since they share intracellular pathways, targeting them simultaneously may inhibit substitute pathways that tumor uses to evade the immune system and resistant to therapeutic agents.

Regulatory mechanisms and potential medical applications of HNF1A-AS1 in cancers

Long noncoding RNAs (lncRNAs) are defined as a class of non-protein-coding RNAs that are longer than 200 nucleotides. Previous studies have shown that lncRNAs play a vital role in the progression of multiple diseases, which highlights their potential for medical applications. The lncRNA hepatocyte nuclear factor 1 homeobox A (HNF1A) antisense RNA 1 (HNF1A-AS1) is known to be abnormally expressed in multiple cancers. HNF1A-AS1 exerts its oncogenic roles through a variety of molecular mechanisms. Moreover, aberrant HNF1A-AS1 expression is associated with diverse clinical features in cancer patients. Therefore, HNF1A-AS1 is a promising biomarker for tumor diagnosis and prognosis and thus a potential candidate for tumor therapy. This review summarizes current studies on the role and the underlying mechanisms of HNF1A-AS1 various cancer types, including gastric cancer, liver cancer, glioma, lung cancer, colorectal cancer, breast cancer, bladder cancer, osteosarcoma, esophageal adenocarcinoma, hemangioma, oral squamous cell carcinoma, laryngeal squamous cell carcinoma, cervical cancer, as well as gastroenteropancreatic neuroendocrine neoplasms. We also describe the diagnostic, prognostic, and therapeutic value of HNF1A-AS1 for multiple cancer patients.

β-Arrestin2 Is Critically Involved in the Differential Regulation of Phosphosignaling Pathways by Thyrotropin-Releasing Hormone and Taltirelin

In recent years, thyrotropin-releasing hormone (TRH) and its analogs, including taltirelin (TAL), have demonstrated a range of effects on the central nervous system that represent potential therapeutic agents for the treatment of various neurological disorders, including neurodegenerative diseases. However, the molecular mechanisms of their actions remain poorly understood. In this study, we investigated phosphosignaling dynamics in pituitary GH1 cells affected by TRH and TAL and the putative role of β-arrestin2 in mediating these effects. Our results revealed widespread alterations in many phosphosignaling pathways involving signal transduction via small GTPases, MAP kinases, Ser/Thr- and Tyr-protein kinases, Wnt/β-catenin, and members of the Hippo pathway. The differential TRH- or TAL-induced phosphorylation of numerous proteins suggests that these ligands exhibit some degree of biased agonism at the TRH receptor. The different phosphorylation patterns induced by TRH or TAL in β-arrestin2-deficient cells suggest that the β-arrestin2 scaffold is a key factor determining phosphorylation events after TRH receptor activation. Our results suggest that compounds that modulate kinase and phosphatase activity can be considered as additional adjuvants to enhance the potential therapeutic value of TRH or TAL.

Identification of Novel Regulatory Regions Induced by Intrauterine Growth Restriction in Rat Islets

Intrauterine growth restriction (IUGR) leads to the development of type 2 diabetes in adulthood, and the permanent alterations in gene expression implicate an epigenetic mechanism. Using a rat model of IUGR, we performed TrueSeq-HELP Tagging to assess the association of DNA methylation changes and gene dysregulation in islets. We identified 511 differentially methylated regions (DMRs) and 4377 significantly altered single CpG sites. Integrating the methylome and our published transcriptome data sets resulted in the identification of pathways critical for islet function. The identified DMRs were enriched with transcription factor binding motifs, such as Elk1, Etv1, Foxa1, Foxa2, Pax7, Stat3, Hnf1, and AR. In silico analysis of 3-dimensional chromosomal interactions using human pancreas and islet Hi-C data sets identified interactions between 14 highly conserved DMRs and 35 genes with significant expression changes at an early age, many of which persisted in adult islets. In adult islets, there were far more interactions between DMRs and genes with significant expression changes identified with Hi-C, and most of them were critical to islet metabolism and insulin secretion. The methylome was integrated with our published genome-wide histone modification data sets from IUGR islets, resulting in further characterization of important regulatory regions of the genome altered by IUGR containing both significant changes in DNA methylation and specific histone marks. We identified novel regulatory regions in islets after exposure to IUGR, suggesting that epigenetic changes at key transcription factor binding motifs and other gene regulatory regions may contribute to gene dysregulation and an abnormal islet phenotype in IUGR rats.

Disease-Associated Risk Variants in ANRIL Are Associated with Tumor-Infiltrating Lymphocyte Presence in Primary Melanomas in the Population-Based GEM Study

BACKGROUND

Genome-wide association studies have reported that genetic variation at ANRIL (CDKN2B-AS1) is associated with risk of several chronic diseases including coronary artery disease, coronary artery calcification, myocardial infarction, and type 2 diabetes mellitus. ANRIL is located at the CDKN2A/B locus, which encodes multiple melanoma tumor suppressors. We investigated the association of these variants with melanoma prognostic characteristics.

METHODS

The Genes, Environment, and Melanoma Study enrolled 3,285 European origin participants with incident invasive primary melanoma. For each of ten disease-associated SNPs at or near ANRIL, we used linear and logistic regression modeling to estimate, respectively, the per allele mean changes in log of Breslow thickness and ORs for presence of ulceration and tumor-infiltrating lymphocytes (TIL). We also assessed effect modification by tumor NRAS/BRAF mutational status.

RESULTS

Rs518394, rs10965215, and rs564398 passed false discovery and were each associated (P ≤ 0.005) with TILs, although only rs564398 was independently associated (P = 0.0005) with TILs. Stratified by NRAS/BRAF mutational status, rs564398*A was significantly positively associated with TILs among NRAS/BRAF mutant, but not wild-type, cases. We did not find SNP associations with Breslow thickness or ulceration.

CONCLUSIONS

ANRIL rs564398 was associated with TIL presence in primary melanomas, and this association may be limited to NRAS/BRAF-mutant cases.

IMPACT

Pathways related to ANRIL variants warrant exploration in relationship to TILs in melanoma, especially given the impact of TILs on immunotherapy and survival.

Side effect prediction based on drug-induced gene expression profiles and random forest with iterative feature selection

One in every ten drug candidates fail in clinical trials mainly due to efficacy and safety related issues, despite in-depth preclinical testing. Even some of the approved drugs such as chemotherapeutics are notorious for their side effects that are burdensome on patients. In order to pave the way for new therapeutics with more tolerable side effects, the mechanisms underlying side effects need to be fully elucidated. In this work, we addressed the common side effects of chemotherapeutics, namely alopecia, diarrhea and edema. A strategy based on Random Forest algorithm unveiled an expression signature involving 40 genes that predicted these side effects with an accuracy of 89%. We further characterized the resulting signature and its association with the side effects using functional enrichment analysis and protein-protein interaction networks. This work contributes to the ongoing efforts in drug development for early identification of side effects to use the resources more effectively.

Superenhancer-transcription factor regulatory network in malignant tumors

Objective

This study aims to identify superenhancer (SE)-transcriptional factor (TF) regulatory network related to eight common malignant tumors based on ChIP-seq data modified by histone H3K27ac in the enhancer region of the SRA database.

Methods

H3K27ac ChIP-seq data of eight common malignant tumor samples were downloaded from the SRA database and subjected to comparison with the human reference genome hg19. TFs regulated by SEs were screened with HOMER software. Core regulatory circuitry (CRC) in malignant tumor samples was defined through CRCmapper software and validated by RNA-seq data in TCGA. The findings were substantiated in bladder cancer cell experiments.

Results

Different malignant tumors could be distinguished through the H3K27ac signal. After SE identification in eight common malignant tumor samples, 35 SE-regulated genes were defined as malignant tumor-specific. SE-regulated specific TFs effectively distinguished the types of malignant tumors. Finally, we obtained 60 CRC TFs, and SMAD3 exhibited a strong H3K27ac signal in eight common malignant tumor samples. In vitro experimental data verified the presence of a SE-TF regulatory network in bladder cancer, and SE-TF regulatory network enhanced the malignant phenotype of bladder cancer cells.

Conclusion

The SE-TF regulatory network with SMAD3 as the core TF may participate in the carcinogenesis of malignant tumors.

Knockdown of RREB1 inhibits cell proliferation via enhanced p16 expression in gastric cancer

Gastric cancer (GC) is the most common gastrointestinal malignancy worldwide. However, the molecular mechanisms of the progression of GC are not fully understood. Ras-responsive element binding protein 1 (RREB1) is an oncogene in many types of cancer that is involved in various biological processes, such as DNA damage repair, cell growth and proliferation, cell differentiation, fat development, and fasting glucose balance. In this study, we demonstrate the role of RREB1 in gastric cancer. First, by immunohistochemistry staining (IHC) and bioinformatics analysis, we demonstrated the expression of RREB1 in gastric cancer and paired normal gastric tissues. Then, we established RREB1 overexpression and knockdown cell lines via lentiviral transfection and detected cell proliferation by using MTT, colony-forming, cell cycle and apoptosis assays in vitro. We demonstrated the effect of RREB1 on cell proliferation in vivo by using a subcutaneous xenograft tumor model in nude mice. Finally, by using Western blotting and IHC, we demonstrated the possible mechanism by which RREB1 affects cell proliferation. The IHC and bioinformatics analyses demonstrated that RREB1 was highly expressed in gastric cancer and showed that RREB1-expressing patients had a larger tumor size and more lymphovascular invasion than RREB1-negative patients. Knockdown of RREB1 inhibited cell proliferation in vivo and in vitro. Knockdown of RREB1 enhanced p16 expression in vivo and in vitro, and p16 expression was negatively related to RREB1 in gastric cancer tissue. RREB1 was highly expressed in gastric cancer, and knockdown of RREB1 inhibited cell proliferation via enhanced p16 expression.

miR-26a enhances colorectal cancer cell growth by targeting RREB1 deacetylation to activate AKT-mediated glycolysis

We previously reported the inhibitory effects of microRNA-26a (miR-26a) on the conversion of pyruvate to acetyl coenzyme A in glucose metabolism by directly targeting pyruvate dehydrogenase protein X component in colorectal cancer (CRC) cells (Chen B et al., BMC Cancer 2014). Here, using microRNA in situ hybridization, we confirmed that miR-26a levels were elevated in 77 human CRC tissue samples and further investigated the key miR-26a-mediated metabolic regulation elements and signaling pathways in CRC cells through quantitative proteomic dissection combined with cancer cell biology and biochemical loss-of-function analysis. We found that AKT transcription signaling was a target pathway via miR-26a-mediated deacetylation modification of Ras-responsive element-binding protein 1 (RREB1) at the Lys-60 residue. miR-26a improved the deacetylation level of RREB1, thus contributing to RREB1 binding to the AKT1 promoter to activate AKT transcription and its related signaling pathway in glycolysis. Moreover, miR-26a promoted CRC tumorigenesis in CRC cells and subcutaneous xenograft mice. Thus, miR-26a is a key regulator of CRC tumorigenesis that mediates the deacetylation modification of RREB1 to enhance AKT1 transcription and downstream target gene expression in glycolysis for CRC growth.

Identifying molecular targets for reverse aging using integrated network analysis of transcriptomic and epigenomic changes during aging

Aging is associated with widespread physiological changes, including skeletal muscle weakening, neuron system degeneration, hair loss, and skin wrinkling. Previous studies have identified numerous molecular biomarkers involved in these changes, but their regulatory mechanisms and functional repercussions remain elusive. In this study, we conducted next-generation sequencing of DNA methylation and RNA sequencing of blood samples from 51 healthy adults between 20 and 74 years of age and identified aging-related epigenetic and transcriptomic biomarkers. We also identified candidate molecular targets that can reversely regulate the transcriptomic biomarkers of aging by reconstructing a gene regulatory network model and performing signal flow analysis. For validation, we screened public experimental data including gene expression profiles in response to thousands of chemical perturbagens. Despite insufficient data on the binding targets of perturbagens and their modes of action, curcumin, which reversely regulated the biomarkers in the experimental dataset, was found to bind and inhibit JUN, which was identified as a candidate target via signal flow analysis. Collectively, our results demonstrate the utility of a network model for integrative analysis of omics data, which can help elucidate inter-omics regulatory mechanisms and develop therapeutic strategies against aging.

The transcription factor Rreb1 regulates epithelial architecture, invasiveness, and vasculogenesis in early mouse embryos

Ras-responsive element-binding protein 1 (Rreb1) is a zinc-finger transcription factor acting downstream of RAS signaling. Rreb1 has been implicated in cancer and Noonan-like RASopathies. However, little is known about its role in mammalian non-disease states. Here, we show that Rreb1 is essential for mouse embryonic development. Loss of Rreb1 led to a reduction in the expression of vasculogenic factors, cardiovascular defects, and embryonic lethality. During gastrulation, the absence of Rreb1 also resulted in the upregulation of cytoskeleton-associated genes, a change in the organization of F-ACTIN and adherens junctions within the pluripotent epiblast, and perturbed epithelial architecture. Moreover, Rreb1 mutant cells ectopically exited the epiblast epithelium through the underlying basement membrane, paralleling cell behaviors observed during metastasis. Thus, disentangling the function of Rreb1 in development should shed light on its role in cancer and other diseases involving loss of epithelial integrity.

RREB1-MKL2 fusion in a spindle cell sinonasal sarcoma: biphenotypic sinonasal sarcoma or ectomesenchymal chondromyxoid tumor in an unusual site?

Biphenotypic sinonasal sarcoma (BSNS) is a rare, low grade spindle cell sarcoma, recently recognized in the WHO classification of head and neck tumors, which is characterized by a dual myogenic and neural differentiation and recurrent gene fusions, often involving PAX3-MAML3, and less commonly PAX3 fusions with other partners such as NCOA1, NCOA2, or WWTR1. Yet, in about 4% of tumors no gene rearrangements are identified. Herein, we describe a RREB1-MKL2 fusion in a BSNS lesion occurring in a 73-year-old female patient with a right maxillo-ethmoidal angle lesion. The polypoid, moderately cellular tumor with infiltrative submucosal growth was composed of fascicles of relatively bland spindle cells embedded in a loose collagenous matrix. The tumor cells showed moderate amounts of eosinophilic cytoplasm with indistinct borders and uniform, pale, ovoid to slender nuclei. The slowly proliferating neoplastic cells co-expressed smooth muscle actin and S100, and showed focal nuclear positivity for ß-catenin, while lacking staining for cytokeratins, desmin, myogenin, caldesmon, glial fibrillary acid protein, and SOX-10. Molecular analysis by targeted RNA-based next-generation sequencing identified an in-frame fusion between exon 8 of RREB1 and exon 11 of MKL2, a genetic event that was reported to be a molecular hallmark of ectomesenchymal chondromyxoid tumor. Gene rearrangements in both genes were independently verified by fluorescence in situ hybridization (FISH). To evaluate its recurrent potential an additional group of 15 fusion negative BSNS were tested for abnormalities in RREB1 and MKL2 genes by FISH, but no additional positive cases were identified.

Pan-cancer analysis of non-coding recurrent mutations and their possible involvement in cancer pathogenesis

Cancer-related mutations have been mainly identified in protein-coding regions. Recent studies have demonstrated that mutations in non-coding regions of the genome could also be a risk factor for cancer. However, the non-coding regions comprise 98% of the total length of the human genome and contain a huge number of mutations, making it difficult to interpret their impacts on pathogenesis of cancer. To comprehensively identify cancer-related non-coding mutations, we focused on recurrent mutations in non-coding regions using somatic mutation data from COSMIC and whole-genome sequencing data from The Cancer Genome Atlas (TCGA). We identified 21 574 recurrent mutations in non-coding regions that were shared by at least two different samples from both COSMIC and TCGA databases. Among them, 580 candidate cancer-related non-coding recurrent mutations were identified based on epigenomic and chromatin structure datasets. One of such mutation was located in RREB1 binding site that is thought to interact with TEAD1 promoter. Our results suggest that mutations may disrupt the binding of RREB1 to the candidate enhancer region and increase TEAD1 expression levels. Our findings demonstrate that non-coding recurrent mutations and coding mutations may contribute to the pathogenesis of cancer.

Genome-Wide Association Analyses of Fertility Traits in Beef Heifers

The ability of livestock to reproduce efficiently is critical to the sustainability of animal agriculture. Antral follicle count (AFC) and reproductive tract scores (RTS) can be used to estimate fertility in beef heifers, but the genetic mechanisms influencing variation in these measures are not well understood. Two genome-wide association studies (GWAS) were conducted to identify the significant loci associated with these traits. In total, 293 crossbred beef heifers were genotyped on the Bovine GGP 50K chip and genotypes were imputed to 836,121 markers. A GWAS was performed with the AFC phenotype for 217 heifers with a multi-locus mixed model, conducted using the year, age at time of sampling and principal component analysis groupings as the covariates. The RTS GWAS was performed with 289 heifers using an additive correlation/trend test comparing prepubertal to pubertal heifers. The loci on chromosomes 2, 3 and 23 were significant in the AFC GWAS and the loci on chromosomes 2, 8, 10 and 11 were significant in the RTS GWAS. The significant region on chromosome 2 was similar between both analyses. These regions contained genes associated with cell proliferation, transcription, apoptosis and development. This study proposes candidate genes for beef cattle fertility, although future research is needed to elucidate the precise mechanisms.

Praja2 suppresses the growth of gastric cancer by ubiquitylation of KSR1 and inhibiting MEK-ERK signal pathways

Gastric cancer (GC) is a common malignant tumor, which has a high incidence and fatality. Therefore, it is important to clarify the molecular mechanism of the occurrence and development for GC and to find more effective treatments and targeted drugs. In this study, we found that the kinase suppressor of Ras1 (KSR1) was increased in GC tissues and cell lines. Silencing of KSR1 inhibited the proliferation, migration and invasion of MKN-45 cells. E3 ligase Praja2 was downregulated in GC tissues and cell lines. In addition, praja2 promoted ubiquitylation of KSR1, but inhibited MEK-ERK signal pathways. Functional analysis indicated overexpression of praja2 inhibited the proliferation, migration and invasion of MKN-45 cells, while MG132 or FGF2 treatment removed the inhibitory effects of praja2 on GC progression. In vivo tumorigenesis experiments indicated praja2 inhibited tumor growth via KSR1-MEK-ERK axis. In conclusion, praja2 promoted the ubiquitylation and degradation of KSR1, which disturbed MEK- ERK signaling and inhibited GC progression. Our study might provide a novel target for GC clinical treatment.

Cancer and SOX proteins: New insight into their role in ovarian cancer progression/inhibition

Transcription factors are potential targets in disease therapy, particularly in cancer. This is due to the fact that transcription factors regulate a variety of cellular events, and their modulation has opened a new window in cancer therapy. Sex-determining region Y (SRY)-related high-mobility group (HMG) box (SOX) proteins are potential transcription factors that are involved in developmental processes such as embryogenesis. It has been reported that abnormal expression of SOX proteins is associated with development of different cancers, particularly ovarian cancer (OC). In the present review, our aim is to provide a mechanistic review of involvement of SOX members in OC. SOX members may suppress and/or promote aggressiveness and proliferation of OC cells. Clinical studies have also confirmed the potential of transcription factors as diagnostic and prognostic factors in OC. Notably, studies have demonstrated the relationship between SOX members and other molecular pathways such as ST6Ga1-I, PI3K, ERK and so on, leading to more complexity. Furthermore, SOX members can be affected by upstream mediators such as microRNAs, long non-coding RNAs, and so on. It is worth mentioning that the expression of each member of SOX proteins is corelated with different stages of OC. Furthermore, their expression determines the response of OC cells to chemotherapy. These topics are discussed in this review to shed some light on role of SOX transcription factors in OC.