RERG is a novel ras-related, estrogen-regulated and growth-inhibitory gene in breast cancer

Construction of a novel radioresistance-related signature for prediction of prognosis, immune microenvironment and anti-tumour drug sensitivity in non-small cell lung cancer

BACKGROUND

Non-small cell lung cancer (NSCLC) is a fatal disease, and radioresistance is an important factor leading to treatment failure and disease progression. The objective of this research was to detect radioresistance-related genes (RRRGs) with prognostic value in NSCLC.

METHODS

The weighted gene coexpression network analysis (WGCNA) and differentially expressed genes (DEGs) analysis were performed to identify RRRGs using expression profiles from TCGA and GEO databases. The least absolute shrinkage and selection operator (LASSO) regression and random survival forest (RSF) were used to screen for prognostically relevant RRRGs. Multivariate Cox regression was used to construct a risk score model. Then, Immune landscape and drug sensitivity were evaluated. The biological functions exerted by the key gene LBH were verified by in vitro experiments.

RESULTS

Ninety-nine RRRGs were screened by intersecting the results of DEGs and WGCNA, then 11 hub RRRGs associated with survival were identified using machine learning algorithms (LASSO and RSF). Subsequently, an eight-gene (APOBEC3B, DOCK4, IER5L, LBH, LY6K, RERG, RMDN2 and TSPAN2) risk score model was established and demonstrated to be an independent prognostic factor in NSCLC on the basis of Cox regression analysis. The immune landscape and sensitivity to anti-tumour drugs showed significant disparities between patients categorized into different risk score subgroups. In vitro experiments indicated that overexpression of LBH enhanced the radiosensitivity of A549 cells, and knockdown LBH reversed the cytotoxicity induced by X-rays.

CONCLUSION

Our study developed an eight-gene risk score model with potential clinical value that can be adopted for choice of drug treatment and prognostic prediction. Its clinical routine use may assist clinicians in selecting more rational practices for individuals, which is important for improving the prognosis of NSCLC patients. These findings also provide references for the development of potential therapeutic targets.

Protein lipidation in the tumor microenvironment: enzymology, signaling pathways, and therapeutics

Protein lipidation is a pivotal post-translational modification that increases protein hydrophobicity and influences their function, localization, and interaction network. Emerging evidence has shown significant roles of lipidation in the tumor microenvironment (TME). However, a comprehensive review of this topic is lacking. In this review, we present an integrated and in-depth literature review of protein lipidation in the context of the TME. Specifically, we focus on three major lipidation modifications: S-prenylation, S-palmitoylation, and N-myristoylation. We emphasize how these modifications affect oncogenic signaling pathways and the complex interplay between tumor cells and the surrounding stromal and immune cells. Furthermore, we explore the therapeutic potential of targeting lipidation mechanisms in cancer treatment and discuss prospects for developing novel anticancer strategies that disrupt lipidation-dependent signaling pathways. By bridging protein lipidation with the dynamics of the TME, our review provides novel insights into the complex relationship between them that drives tumor initiation and progression.

MicroRNA-382 as a tumor suppressor during tumor progression

Despite the recent progresses in therapeutic and diagnostic methods, there is still a significantly high rate of mortality among cancer patients. One of the main reasons for the high mortality rate in cancer patients is late diagnosis, which leads to the failure of therapeutic strategies. Therefore, investigation of cancer biology can lead to the introduction of early diagnostic markers in these patients. MicroRNAs (miRNAs) play an important role in regulation of cellular processes associated with tumor progression. Due to the high stability of miRNAs in body fluids, these factors can be considered as the non-invasive tumor markers. Deregulation of miR-382 has been widely reported in different cancers. Therefore, in this review, we investigated the role of miR-382 during tumor development. It has shown that miR-382 has mainly a tumor suppressive, which inhibits the growth of tumor cells through the regulation of signaling pathways, RNA-binding proteins, and transcription factors. Therefore, miR-382 can be suggested as a diagnostic and therapeutic marker in cancer patients.

Perfluorooctanesulfonic acid (PFOS) induced cancer related DNA methylation alterations in human breast cells: A whole genome methylome study

DNA methylation plays a pivotal role in cancer. The ubiquitous contaminant perfluorooctanesulfonic acid (PFOS) has been epidemiologically associated with breast cancer, and can induce proliferation and malignant transformation of normal human breast epithelial cells (MCF-10A), but the information about its effect on DNA methylation is sparse. The aim of this study was to characterize the whole-genome methylome effects of PFOS in our breast cell model and compare the findings with previously demonstrated DNA methylation alterations in breast tumor tissues. The DNA methylation profile was assessed at single CpG resolution in MCF-10A cells treated with 1 μM PFOS for 72 h by using Enzymatic Methyl sequencing (EM-seq). We found 12,591 differentially methylated CpG-sites and 13,360 differentially methylated 100 bp tiles in the PFOS exposed breast cells. These differentially methylated regions (DMRs) overlapped with 2406 genes of which 494 were long non-coding RNA and 1841 protein coding genes. We identified 339 affected genes that have been shown to display altered DNA methylation in breast cancer tissue and several other genes related to cancer development. This includes hypermethylation of GACAT3, DELEC1, CASC2, LCIIAR, MUC16, SYNE1 and hypomethylation of TTN and KMT2C. DMRs were also found in estrogen receptor genes (ESR1, ESR2, ESRRG, ESRRB, GREB1) and estrogen responsive genes (GPER1, EEIG1, RERG). The gene ontology analysis revealed pathways related to cancer phenotypes such as cell adhesion and growth. These findings improve the understanding of PFOS's potential role in breast cancer and illustrate the value of whole-genome methylome analysis in uncovering mechanisms of chemical effects, identifying biomarker candidates, and strengthening epidemiological associations, potentially impacting risk assessment.

Genome-wide association study for growth traits with 1066 individuals in largemouth bass (Micropterus salmoides)

The largemouth bass is a native species of North America that was first introduced to mainland China in the 1980s. In recent years, it has been extensively farmed in China due to its high meat quality and broad adaptability. In this study, we collected growth trait data from 1,066 largemouth bass individuals across two populations. We generated an average of approximately 7× sequencing coverage for these fish using Illumina sequencers. From the samples, we identified 2,695,687 SNPs and retained 1,809,116 SNPs for further analysis after filtering. To estimate the number of genome-wide effective SNPs, we performed LD pruning with PLINK software and identified 77,935 SNPs. Our GWAS revealed 15 SNPs associated with six growth traits. We identified a total of 24 genes related to growth, with three genes-igf1, myf5, and myf6-directly associated with skeletal muscle development and growth, located near the leading SNP on chromosome 23. Other candidate genes are involved in the development of tissues and organs or other physiological processes. These findings provide a valuable set of SNPs and genes that could be useful for genetic breeding programs aimed at enhancing growth in largemouth bass.

DNA Hypomethylation Underlies Epigenetic Swapping between AGO1 and AGO1-V2 Isoforms in Tumors

Human tumors progress in part by accumulating epigenetic alterations, which include gains and losses of DNA methylation in different parts of the cancer cell genome. Recent work has revealed a link between these two opposite alterations by showing that DNA hypomethylation in tumors can induce the expression of transcripts that overlap downstream gene promoters and thereby induce their hypermethylation. Preliminary in silico evidence prompted us to investigate if this mechanism applies to the locus harboring AGO1, a gene that plays a central role in miRNA biogenesis and RNA interference. Inspection of public RNA-Seq datasets and RT-qPCR experiments show that an alternative transcript starting 13.4 kb upstream of AGO1 (AGO1-V2) is expressed specifically in testicular germ cells, and becomes aberrantly activated in different types of tumors, particularly in tumors of the esophagus, stomach, and lung. This expression pattern classifies AGO1-V2 into the group of "Cancer-Germline" (CG) genes. Analysis of transcriptomic and methylomic datasets provided evidence that transcriptional activation of AGO1-V2 depends on DNA demethylation of its promoter region. Western blot experiments revealed that AGO1-V2 encodes a shortened isoform of AGO1, corresponding to a truncation of 75 aa in the N-terminal domain, and which we therefore referred to as "∆NAGO1". Interestingly, significant correlations between hypomethylation/activation of AGO1-V2 and hypermethylation/repression of AGO1 were observed upon examination of tumor cell lines and tissue datasets. Overall, our study reveals the existence of a process of interdependent epigenetic alterations in the AGO1 locus, which promotes swapping between two AGO1 protein-coding mRNA isoforms in tumors.

TMO-Net: an explainable pretrained multi-omics model for multi-task learning in oncology

Cancer is a complex disease composing systemic alterations in multiple scales. In this study, we develop the Tumor Multi-Omics pre-trained Network (TMO-Net) that integrates multi-omics pan-cancer datasets for model pre-training, facilitating cross-omics interactions and enabling joint representation learning and incomplete omics inference. This model enhances multi-omics sample representation and empowers various downstream oncology tasks with incomplete multi-omics datasets. By employing interpretable learning, we characterize the contributions of distinct omics features to clinical outcomes. The TMO-Net model serves as a versatile framework for cross-modal multi-omics learning in oncology, paving the way for tumor omics-specific foundation models.

Molecular features of luminal breast cancer defined through spatial and single-cell transcriptomics

BACKGROUND

Intratumour heterogeneity is a hallmark of most solid tumours, including breast cancers. We applied spatial transcriptomics and single-cell RNA-sequencing on patient-derived xenografts (PDXs) to profile spatially resolved cell populations within oestrogen receptor-positive (ER+ ) breast cancer and to elucidate their importance in oestrogen-dependent tumour growth.

METHODS

Two PDXs of 'ER-high' breast cancers with opposite oestrogen-mediated growth responses were investigated: oestrogen-suppressed GS3 (80-100% ER) and oestrogen-dependent SC31 (40-90% ER) models. The observation was validated via single-cell analyses on an 'ER-low' PDX, GS1 (5% ER). The results from our spatial and single-cell analyses were further supported by a public ER+ breast cancer single-cell dataset and protein-based dual immunohistochemistry (IHC) of SC31 examining important luminal cancer markers (i.e., ER, progesterone receptor and Ki67). The translational implication of our findings was assessed by clinical outcome analyses on publicly available cohorts.

RESULTS

Our space-gene-function study revealed four spatially distinct compartments within ER+ breast cancers. These compartments showed functional diversity (oestrogen-responsive, proliferative, hypoxia-induced and inflammation-related). The 'proliferative' population, rather than the 'oestrogen-responsive' compartment, was crucial for oestrogen-dependent tumour growth, leading to the acquisition of luminal B-like features. The cells expressing typical oestrogen-responsive genes like PGR were not directly linked to oestrogen-dependent proliferation. Dual IHC analyses demonstrated the distinct contribution of the Ki67+ proliferative cells toward oestrogen-mediated growth and their response to a CDK4/6 inhibitor. The gene signatures derived from the proliferative, hypoxia-induced and inflammation-related compartments were significantly correlated with worse clinical outcomes, while patients with the oestrogen-responsive signature showed better prognoses, suggesting that this compartment would not be directly associated with oestrogen-dependent tumour progression.

CONCLUSIONS

Our study identified the gene signature in our 'proliferative' compartment as an important determinant of luminal cancer subtypes. This 'proliferative' cell population is a causative feature of luminal B breast cancer, contributing toward its aggressive behaviours.

microRNA-382 as a tumor suppressor? Roles in tumorigenesis and clinical significance

MicroRNAs (miRNAs) are small single-stranded RNAs belonging to a class of non-coding RNAs with an average length of 18-22 nucleotides. Although not able to encode any protein, miRNAs are vastly studied and found to play role in various human physiologic as well as pathological conditions. A huge number of miRNAs have been identified in human cells whose expression is straightly regulated with crucial biological functions, while this number is constantly increasing. miRNAs are particularly studied in cancers, where they either can act with oncogenic function (oncomiRs) or tumor-suppressors role (referred as tumor-suppressor/oncorepressor miRNAs). miR-382 is a well-studied miRNA, which is revealed to play regulatory roles in physiological processes like osteogenic differentiation, hematopoietic stem cell differentiation and normal hematopoiesis, and liver progenitor cell differentiation. Notably, miR-382 deregulation is reported in pathologic conditions, such as renal fibrosis, muscular dystrophies, Rett syndrome, epidural fibrosis, atrial fibrillation, amelogenesis imperfecta, oxidative stress, human immunodeficiency virus (HIV) replication, and various types of cancers. The majority of oncogenesis studies have claimed miR-382 downregulation in cancers and suppressor impact on malignant phenotype of cancer cells in vitro and in vivo, while a few studies suggest opposite findings. Given the putative role of this miRNA in regulation of oncogenesis, assessment of miR-382 expression is suggested in a several clinical investigations as a prognostic/diagnostic biomarker for cancer patients. In this review, we have an overview to recent studies evaluated the role of miR-382 in oncogenesis as well as its clinical potential.

Multiparity increases the risk of diabetes by impairing the proliferative capacity of pancreatic β cells

Pregnancy imposes a substantial metabolic burden on women, but little is known about whether or how multiple pregnancies increase the risk of maternal postpartum diabetes. In this study, we assessed the metabolic impact of multiple pregnancies in humans and in a rodent model. Mice that underwent multiple pregnancies had increased adiposity, but their glucose tolerance was initially improved compared to those of age-matched virgin mice. Later, however, insulin resistance developed over time, but insulin secretory function and compensatory pancreatic β cell proliferation were impaired in multiparous mice. The β cells of multiparous mice exhibited aging features, including telomere shortening and increased expression of Cdkn2a. Single-cell RNA-seq analysis revealed that the β cells of multiparous mice exhibited upregulation of stress-related pathways and downregulation of cellular respiration- and oxidative phosphorylation-related pathways. In humans, women who delivered more than three times were more obese, and their plasma glucose concentrations were elevated compared to women who had delivered three or fewer times, as assessed at 2 months postpartum. The disposition index, which is a measure of the insulin secretory function of β cells, decreased when women with higher parity gained body weight after delivery. Taken together, our findings indicate that multiple pregnancies induce cellular stress and aging features in β cells, which impair their proliferative capacity to compensate for insulin resistance.

Missing link between tissue specific expressing pattern of ERβ and the clinical manifestations in LGBLEL

Purpose

Lacrimal gland benign lymphoepithelial lesion (LGBLEL) is an IgG4-related disease of unknown etiology with a risk for malignant transformation. Estrogen is considered to be related to LGBLEL onset.

Methods

Seventy-eight LGBLEL and 13 control clinical samples were collected and studied to determine the relationship between estrogen and its receptors and LGBLEL development.

Results

The serological analysis revealed no significant differences in the levels of three estrogens be-tween the LGBLEL and control groups. However, immunohistochemical analyses indicated that the expression levels of ERβ and its downstream receptor RERG were relatively lower in LGBLEL samples than in control samples, with higher expression in the lacrimal gland and lower expression in the lymphocyte infiltration region. However, low expression of ERα was detected. The transcriptome sequence analysis revealed upregulated genes associated with LGBLEL enriched in lymphocyte proliferation and activation function; downregulated genes were enriched in epithelial and vascular proliferation functions. The key genes and gene networks were further analyzed. Interactions between B cells and epithelial cells were analyzed due to the identified involvement of leukocyte subsets and epithelial cells. B cell proliferation was found to potentially contribute to lacrimal gland apoptosis.

Conclusion

Therefore, the tissue-heterogeneous expression pattern of ERβ is potentially related to the clinical manifestations and progression of LGBLEL, although further investigations are required to confirm this finding.

The genome of a vestimentiferan tubeworm (Ridgeia piscesae) provides insights into its adaptation to a deep-sea environment

BACKGROUND

Vestimentifera (Polychaeta, Siboglinidae) is a taxon of deep-sea worm-like animals living in deep-sea hydrothermal vents, cold seeps, and organic falls. The morphology and lifespan of Ridgeia piscesae, which is the only vestimentiferan tubeworm species found in the hydrothermal vents on the Juan de Fuca Ridge, vary greatly according to endemic environment. Recent analyses have revealed the genomic basis of adaptation in three vent- and seep-dwelling vestimentiferan tubeworms. However, the evolutionary history and mechanism of adaptation in R. piscesae, a unique species in the family Siboglinidae, remain to be investigated.

RESULT

We assembled a draft genome of R. piscesae collected at the Cathedral vent of the Juan de Fuca Ridge. Comparative genomic analysis showed that vent-dwelling tubeworms with a higher growth rate had smaller genome sizes than seep-dwelling tubeworms that grew much slower. A strong positive correlation between repeat content and genome size but not intron size and the number of protein-coding genes was identified in these deep-sea tubeworm species. Evolutionary analysis revealed that Ridgeia pachyptila and R. piscesae, the two tubeworm species that are endemic to hydrothermal vents of the eastern Pacific, started to diverge between 28.5 and 35 million years ago. Four genes involved in cell proliferation were found to be subject to positive selection in the genome of R. piscesae.

CONCLUSION

Ridgeia pachyptila and R. piscesae started to diverge after the formation of the Gorda/Juan de Fuca/Explorer ridge systems and the East Pacific Rise. The high growth rates of vent-dwelling tubeworms might be derived from their small genome sizes. Cell proliferation is important for regulating the growth rate in R. piscesae.

Identification and characterization of a proliferative cell population in estrogen receptor-positive metastatic breast cancer through spatial and single-cell transcriptomics

Background

Intratumor heterogeneity is a hallmark of most solid tumors, including breast cancers. We applied spatial transcriptomics and single-cell RNA-sequencing technologies to profile spatially resolved cell populations within estrogen receptor-positive (ER ⁺ ) metastatic breast cancers and elucidate their importance in estrogen-dependent tumor growth.

Methods

Spatial transcriptomics and single-cell RNA-sequencing were performed on two patient-derived xenografts (PDXs) of "ER-high" metastatic breast cancers with opposite estrogen-mediated growth responses: estrogen-suppressed GS3 (80-100% ER) and estrogen-stimulated SC31 (30-75% ER) models. The analyses included samples treated with and without 17β-estradiol. The findings were validated via scRNA-seq analyses on "ER-low" estrogen-accelerating PDX, GS1 (5% ER). The results from our spatial and single-cell analyses were further supported by the analysis of a publicly available single cell dataset and a protein-based dual immunohistochemical (IHC) evaluation using three important clinical markers [i.e., ER, progesterone receptor (PR), and Ki67]. The translational implication of these results was assessed by clinical outcome analyses on public breast cancer cohorts.

Results

Our novel space-gene-function study revealed a "proliferative" cell population in addition to three major spatially distinct compartments within ER ⁺ metastatic breast cancers. These compartments showed functional diversity (i.e., estrogen-responsive, proliferative, hypoxia-induced, and inflammation-related). The "proliferative ( MKI67 ⁺ )" population, not "estrogen-responsive" compartment, was crucial for estrogen-dependent tumor growth, leading to the acquisition of luminal B features. The cells with induction of typical estrogen-responsive genes such as PGR were not directly linked to estrogen-dependent proliferation. Additionally, the dual IHC analyses demonstrated the distinct contribution of the Ki67 ⁺ proliferative cells toward estrogen-mediated growth and their response to palbociclib, a CDK4/6 inhibitor. The gene signatures developed from the proliferative, hypoxia-induced, and inflammation-related compartments were significantly correlated with worse clinical outcomes, while patients with the high estrogen-responsive scores showed better prognosis, confirming that the estrogen-responsive compartment would not be directly associated with estrogen-dependent tumor progression.

Conclusions

For the first time, our study elucidated a "proliferative" cell population distinctly distributed in ER ⁺ metastatic breast cancers. They contribute differently toward progression of these cancers, and the gene signature in the "proliferative" compartment is an important determinant of luminal cancer subtypes.

Absence of Scaffold Protein Tks4 Disrupts Several Signaling Pathways in Colon Cancer Cells

Tks4 is a large scaffold protein in the EGFR signal transduction pathway that is involved in several cellular processes, such as cellular motility, reactive oxygen species-dependent processes, and embryonic development. It is also implicated in a rare developmental disorder, Frank-ter Haar syndrome. Loss of Tks4 resulted in the induction of an EMT-like process, with increased motility and overexpression of EMT markers in colorectal carcinoma cells. In this work, we explored the broader effects of deletion of Tks4 on the gene expression pattern of HCT116 colorectal carcinoma cells by transcriptome sequencing of wild-type and Tks4 knockout (KO) cells. We identified several protein coding genes with altered mRNA levels in the Tks4 KO cell line, as well as a set of long non-coding RNAs, and confirmed these changes with quantitative PCR on a selected set of genes. Our results show a significant perturbation of gene expression upon the deletion of Tks4, suggesting the involvement of different signal transduction pathways over the well-known EGFR signaling.

The NFAT3/RERG Complex in Luminal Breast Cancers Is Required to Inhibit Cell Invasion and May Be Correlated With an Absence of Axillary Lymph Nodes Colonization

Luminal breast cancers represent 70% of newly diagnosed breast cancers per annum and have a relatively good prognosis compared with triple-negative breast cancers. Luminal tumors that are responsive to hormonal therapy are particularly associated with a favorable prognosis. Nonetheless, the absolute number of metastatic relapses in luminal cancers is larger than in triple-negative breast cancers. A better understanding of the biology of luminal cancers, control of metastases formation, and identification of predictive markers of their evolution are therefore still necessary. In this context, we previously disclosed the key role of NFAT3 in regulating luminal breast cancer invasion. We have now identified a specific inhibitory region, in the C-terminal part of NFAT3, required for the inhibition of invasion of the human luminal breast cancer cell line T-47D. Indeed, we showed that this 85 amino acid C-terminal region acts as a dominant negative form of NFAT3 and that its overexpression in the T-47D cell line led to increased cell invasion. Mechanistically, we have revealed that this region of NFAT3 interacts with the small Ras GTPase RERG (RAS like estrogen regulated growth inhibitor) and shown that RERG expression is required for NFAT3 to impede T-47D cell invasion. We have validated the association of NFAT3 with RERG in human luminal breast cancer tissues. We have shown an increase of the quantity of the NFAT3/RERG complexes in patients without axillary lymph node colonization and therefore proposed that the detection of this complex may be a non-invasive marker of axillary lymph node colonization.

Functional diversity in the RAS subfamily of small GTPases

RAS small GTPases regulate important signalling pathways and are notorious drivers of cancer development and progression. While most research to date has focused on understanding and addressing the oncogenic potential of three RAS oncogenes: HRAS, KRAS, and NRAS; the full RAS subfamily is composed of 35 related GTPases with diverse cellular functions. Most remain deeply understudied despite strong evolutionary conservation. Here, we highlight a group of 17 poorly characterized RAS GTPases that are frequently down-regulated in cancer and evidence suggests may function not as oncogenes, but as tumour suppressors. These GTPases remain largely enigmatic in terms of their cellular function, regulation, and interaction with effector proteins. They cluster within two families we designate as 'distal-RAS' (D-RAS; comprised of DIRAS, RASD, and RASL10) and 'CaaX-Less RAS' (CL-RAS; comprised of RGK, NKIRAS, RERG, and RASL11/12 GTPases). Evidence of a tumour suppressive role for many of these GTPases supports the premise that RAS subfamily proteins may collectively regulate cellular proliferation.

Identification of a novel prognosis-associated ceRNA network in lung adenocarcinoma via bioinformatics analysis

Background

Lung adenocarcinoma (LUAD) is the most common subtype of nonsmall-cell lung cancer (NSCLC) and has a high incidence rate and mortality. The survival of LUAD patients has increased with the development of targeted therapeutics, but the prognosis of these patients is still poor. Long noncoding RNAs (lncRNAs) play an important role in the occurrence and development of LUAD. The purpose of this study was to identify novel abnormally regulated lncRNA–microRNA (miRNA)–messenger RNA (mRNA) competing endogenous RNA (ceRNA) networks that may suggest new therapeutic targets for LUAD or relate to LUAD prognosis.

Methods

We used the SBC human ceRNA array V1.0 to screen for differentially expressed (DE) lncRNAs and mRNAs in four paired LUAD samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to annotate the DE lncRNAs and mRNAs. R bioinformatics packages, The Cancer Genome Atlas (TCGA) LUAD database, and Kaplan–Meier (KM) survival analysis tools were used to validate the microarray data and construct the lncRNA–miRNA–mRNA ceRNA regulatory network. Then, quantitative real-time PCR (qRT-PCR) was used to validate the DE lncRNAs in 7 LUAD cell lines.

Results

A total of 2819 DE lncRNAs and 2396 DE mRNAs (P < 0.05 and fold change ≥ 2 or ≤ 0.5) were identified in four paired LUAD tissue samples. In total, 255 of the DE lncRNAs were also identified in TCGA. The GO and KEGG analysis results suggested that the DE genes were most enriched in angiogenesis and cell proliferation, and were closely related to human cancers. Moreover, the differential expression of ENST00000609697, ENST00000602992, and NR_024321 was consistent with the microarray data, as determined by qRT-PCR validation in 7 LUAD cell lines; however, only ENST00000609697 was associated with the overall survival of LUAD patients (log-rank P = 0.029). Finally, through analysis of ENST00000609697 target genes, we identified the ENST00000609697–hsa-miR-6791-5p–RASL12 ceRNA network, which may play a tumor-suppressive role in LUAD.

Conclusion

ENST00000609697 was abnormally expressed in LUAD. Furthermore, downregulation of ENST00000609697 and its target gene RASL12 was associated with poor prognosis in LUAD. The ENST00000609697–hsa-miR-6791-5p–RASL12 axis may play a tumor-suppressive role. These results suggest new potential prognostic and therapeutic biomarkers for LUAD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12938-021-00952-x.

Single- and double-hit events in genes encoding immune targets before and after T cell-engaging antibody therapy in MM

T cell-engaging immunotherapies exert unprecedented single-agent activity in multiple myeloma (MM), thereby putting a yet unexplored selective pressure on the clonal architecture. In this study, we report on homozygous BCMA (TNFRSF17) gene deletion after BCMA-targeting T cell-redirecting bispecific antibody therapy in a heavily pretreated MM patient. Loss of BCMA protein expression persisted over subsequent relapses, with no response to treatment with anti-BCMA antibody drug conjugate. In light of the multiple alternative targets that are emerging in addition to BCMA, we extended our analyses to delineate a more complete picture of genetic alterations that may have an impact on immunotherapy targets in MM. We performed whole-genome sequencing and RNA sequencing in 100 MM patients (50 were newly diagnosed; 50 were relapsed/refractory) and identified a significant proportion of patients with aberrations in genes encoding immunotherapy targets; GPRC5D ranked first with 15% heterozygous deletions, followed by CD38 (10%), SDC1 (5%), and TNFRSF17 (4%). Notably, these heterozygous deletions did not lower the expression levels of respective genes, but they may represent a first hit that drives the acquisition of homozygous deletions and subsequent antigen-loss relapse upon targeted immunotherapy. In summary, we show preexisting vulnerability in genes encoding immunotargets before and homozygous deletions after T cell-engaging immunotherapy.

Transcriptional overlap links DNA hypomethylation with DNA hypermethylation at adjacent promoters in cancer

Tumor development involves alterations in DNA methylation patterns, which include both gains (hypermethylation) and losses (hypomethylation) in different genomic regions. The mechanisms underlying these two opposite, yet co-existing, alterations in tumors remain unclear. While studying the human MAGEA6/GABRA3 gene locus, we observed that DNA hypomethylation in tumor cells can lead to the activation of a long transcript (CT-GABRA3) that overlaps downstream promoters (GABRQ and GABRA3) and triggers their hypermethylation. Overlapped promoters displayed increases in H3K36me3, a histone mark deposited during transcriptional elongation and known to stimulate de novo DNA methylation. Consistent with such a processive mechanism, increases in H3K36me3 and DNA methylation were observed over the entire region covered by the CT-GABRA3 overlapping transcript. Importantly, experimental induction of CT-GABRA3 by depletion of DNMT1 DNA methyltransferase, resulted in a similar pattern of regional DNA hypermethylation. Bioinformatics analyses in lung cancer datasets identified other genomic loci displaying this process of coupled DNA hypo/hypermethylation, and some of these included tumor suppressor genes, e.g. RERG and PTPRO. Together, our work reveals that focal DNA hypomethylation in tumors can indirectly contribute to hypermethylation of nearby promoters through activation of overlapping transcription, and establishes therefore an unsuspected connection between these two opposite epigenetic alterations.

DNA methylation association with stage progression of head and neck squamous cell carcinoma

Head and Neck Squamous Cell Carcinoma (HNSCC) is the sixth most common cancer worldwide, which accounts for approximately 6% of all cases and is responsible for an estimated 2% of all cancer deaths. Despite progress in the treatment of squamous cell carcinomas, survival rates remain low. It is a fact that epigenetic modifications have numerous associations with biological processes and complex diseases such as cancer. Hence, a more systematic approach is needed to provide potential screening targets and have an effective therapy method. This study developed a workflow to analyze HM450 methylation arrays with mRNA expression profiles that identified novel signatures of epigenetic regulators for tumor progression. We identified differentially expressed genes and differentially methylated regions and the correlation between associated genes to identify epigenetic modifications underlying regulation roles. We have taken the differentiation direction of expressions into account during the integration of gene expression and DNA methylation modification to detect epigenetic regulators of core genes of tumor-stage progression. Enrichment analysis of selected key genes provides better insight into their functionality. Thus, we have investigated gene copy number alteration and mutations to filter differentially expressed genes, including some members of the fibroblast growth factor family and cyclin-dependent kinase inhibitor family with other potential known regulators. Our analysis has revealed the list of 61 commercial methylation probes positively correlated with 31 differentially expressed genes, which can be associated with HNSC metastasis stages. Most of these genes have already reported potential epigenetic regulators, and their role in cancer progression was studied. We suggest these selected probes of DNA methylation as potential targets of the epigenetic regulators in revealed genes that have displayed significant genetic and epigenetic modification behavior during cancer stage progression and tumor metastasis.

Association of RERG Expression with Female Survival Advantage in Malignant Pleural Mesothelioma

Sex differences in incidence, prognosis, and treatment response have been described for many cancers. In malignant pleural mesothelioma (MPM), a lethal disease associated with asbestos exposure, men outnumber women 4 to 1, but women consistently live longer than men following surgery-based therapy. This study investigated whether tumor expression of genes associated with estrogen signaling could potentially explain observed survival differences. Two microarray datasets of MPM tumors were analyzed to discover estrogen-related genes associated with survival. A validation cohort of MPM tumors was selected to balance the numbers of men and women and control for competing prognostic influences. The RAS like estrogen regulated growth inhibitor (RERG) gene was identified as the most differentially-expressed estrogen-related gene in these tumors and predicted prognosis in discovery datasets. In the sex-matched validation cohort, low RERG expression was significantly associated with increased risk of death among women. No association between RERG expression and survival was found among men, and no relationship between estrogen receptor protein or gene expression and survival was found for either sex. Additional investigations are needed to elucidate the molecular mechanisms underlying this association and its sex specificity.

The RIT1 C-terminus associates with lipid bilayers via charge complementarity

RIT1 is a member of the Ras superfamily of small GTPases involved in regulation of cellular signaling. Mutations to RIT1 are involved in cancer and developmental disorders. Like many Ras subfamily members, RIT1 is localized to the plasma membrane. However, RIT1 lacks the C-terminal prenylation that helps many other subfamily members adhere to cellular membranes. We used molecular dynamics simulations to examine the mechanisms by which the C-terminal peptide (CTP) of RIT1 associates with lipid bilayers. We show that the CTP is unstructured and that its membrane interactions depend on lipid composition. While a 12-residue region of the CTP binds strongly to anionic bilayers containing phosphatidylserine lipids, the CTP termini fray from the membrane allowing for accommodation of the RIT1 globular domain at the membrane-water interface.

Auto-antibodies against apolipoprotein A-1 block cancer cells proliferation and induce apoptosis

Auto-antibodies against apoA-1 (anti-apoA-1 IgGs) have been identified as important actors of atherosclerosis development through pro-inflammatory and pro-atherogenic properties and to also induce apoptosis in tumoral neuronal and lymphocyte derived cell lines through unknown mechanisms. The purpose of this study was to explore the cellular pathways involved in tumoral cell survival modulated by anti-apoA-1 antibodies. We observed that anti-apoA-1 antibodies induce growth arrest (in G2/M phase) and cell apoptosis through caspase 3 activation, accompanied by a selective p53 phosphorylation on serine 15. RNA sequencing indicated that anti-apoA-1 IgGs affect the expression of more than 950 genes belonging to five major groups of genes and respectively involved in i) cell proliferation inhibition, ii) p53 stabilisation and regulation, iii) apoptosis regulation, iv) inflammation regulation, and v) oxidative stress. In conclusion, anti-apoA-1 antibodies seem to have a role in blocking tumoral cell proliferation and survival, by activating a major tumor suppressor protein and by modulating the inflammatory and oxidative stress response. Further investigations are needed to explore a possible anti-cancer therapeutic approach of these antibodies in very specific and circumscribed conditions.

Identification of prognosis-related genes in the tumor microenvironment of stomach adenocarcinoma by TCGA and GEO datasets

Accumulating evidence has demonstrated that tumor microenvironment (TME) plays a crucial role in stomach adenocarcinoma (STAD) development, progression, prognosis and immunotherapeutic responses. How the genes in TME interact and behave is extremely crucial for tumor investigation. In the present study, we used gene expression data of STAD available from TCGA and GEO datasets to infer tumor purity using ESTIMATE algorithms, and predicted the associations between tumor purity and clinical features and clinical outcomes. Next, we calculated the differentially expressed genes (DEGs) from the comparisons of immune and stromal scores, and postulated key biological processes and pathways that the DEGs mainly involved in. Then, we analyzed the prognostic values of DEGs in TCGA dataset, and validated the results by GEO dataset. Finally, we used CIBERSORT computational algorithm to estimate the 22 tumor infiltrating immune cells (TIICs) subsets in STAD tissues. We found that stromal and immune scores were significantly correlated with STAD subtypes, clinical stages, Helicobacter polyri infection, and stromal scores could predict the clinical outcomes in STAD patients. Moreover, we screened 307 common DEGs in TCGA and GSE51105 datasets. In the prognosis analyses, we only found OGN, JAM2, RERG, OLFML2B, and ADAMTS1 genes were significantly associated with overall survival in TCGA and GSE84437 datasets, and these genes were correlated with the fractions of T cells, B cells, macrophages, monocytes, NK cells and DC cells, respectively. Our comprehensive analyses for transcriptional data not only improved the understanding of characteristics of TME, but also provided the targets for individual therapy in STAD.

Advances in targeted therapy mainly based on signal pathways for nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is a malignant epithelial carcinoma of the head and neck region which mainly distributes in southern China and Southeast Asia and has a crucial association with the Epstein-Barr virus. Based on epidemiological data, both incidence and mortality of NPC have significantly declined in recent decades grounded on the improvement of living standard and medical level in an endemic region, in particular, with the clinical use of individualized chemotherapy and intensity-modulated radiotherapy (IMRT) which profoundly contributes to the cure rate of NPC patients. To tackle the challenges including local recurrence and distant metastasis in the current NPC treatment, we discussed the implication of using targeted therapy against critical molecules in various signal pathways, and how they synergize with chemoradiotherapy in the NPC treatment. Combination treatment including targeted therapy and IMRT or concurrent chemoradiotherapy is presumably to be future options, which may reduce radiation or chemotherapy toxicities and open new avenues for the improvement of the expected functional outcome for patients with advanced NPC.

RASSF effectors couple diverse RAS subfamily GTPases to the Hippo pathway

Small guanosine triphosphatases (GTPases) of the RAS superfamily signal by directly binding to multiple downstream effector proteins. Effectors are defined by a folded RAS-association (RA) domain that binds exclusively to GTP-loaded (activated) RAS, but the binding specificities of most RA domains toward more than 160 RAS superfamily GTPases have not been characterized. Ten RA domain family (RASSF) proteins comprise the largest group of related effectors and are proposed to couple RAS to the proapoptotic Hippo pathway. Here, we showed that RASSF1-6 formed complexes with the Hippo kinase ortholog MST1, whereas RASSF7-10 formed oligomers with the p53-regulating effectors ASPP1 and ASPP2. Moreover, only RASSF5 bound directly to activated HRAS and KRAS, and RASSFs did not augment apoptotic induction downstream of RAS oncoproteins. Structural modeling revealed that expansion of the RASSF effector family in vertebrates included amino acid substitutions to key residues that direct GTPase-binding specificity. We demonstrated that the tumor suppressor RASSF1A formed complexes with the RAS-related GTPases GEM, REM1, REM2, and the enigmatic RASL12. Furthermore, interactions between RASSFs and RAS GTPases blocked YAP1 nuclear localization. Thus, these simple scaffolds link the activation of diverse RAS family small G proteins to Hippo or p53 regulation.

NF-κB-mediated miR-650 plays oncogenic roles and activates AKT/ERK/NF-κB pathways by targeting RERG in glioma cells

PURPOSE

Glioma is the most common cancer in the central nervous system and has a high mortality rate. Despite advances that have been made in the treatment of glioma, its prognosis still remains poor. Dysregulation of miRNAs has been reported in many cancers, including glioma. Here, we set out to assess the role of miR-650 in glioma, including its diagnostic and therapeutic potential.

METHODS

miR-650 and RAS-like estrogen-regulated growth inhibitor (RERG) expression levels were analyzed using qRT-PCR in primary glioma tissues and cell lines. Cell Counting Kit-8, 5-ethynyl-2'-deoxyuridine, colony formation, Western blotting, scratch wound healing, Transwell, adhesion, autophagy, immunofluorescence, luciferase reporter, electrophoretic mobility shift, tumor xenograft and flow cytometry assays were employed to investigate the mechanisms underlying the effect of miR-650 and RERG on glioma development.

RESULTS

miR-650 was found to be up-regulated in glioma tissues and cell lines compared to non-cancerous brain tissues and neural progenitor cells, respectively. We also found that miR-650 promoted cell proliferation, migration and invasion in glioma cells, and enhanced glioma tumor formation and growth in vivo. We identified and validated RERG as a direct target of miR-650. RERG was shown to act as a tumor suppressor in glioma cells, and its suppressor roles were rescued by miR-650. We found that nuclear factor (NF)-κB bound to the promoter of miR-650 and enhanced its expression. PH domain and leucine rich repeat protein phosphatase 2 (PHLPP2), as a co-factor of the RERG/PHLPP2 complex, mediated miR-650-induced activation of the protein kinase B/extracellular-signal-regulated kinase/NF-κB signaling pathways.

CONCLUSIONS

Our data revealed novel functional roles for miR-650 in glioma development and may provide new avenues for future clinical applications.

Overexpression of Long Non-Coding RNA FGF14-AS2 Inhibits Colorectal Cancer Proliferation Via the RERG/Ras/ERK Signaling by Sponging microRNA-1288-3p

Colorectal cancer remains one of most common cancer types with poor prognosis globally. Recent years, numerous studies depicted pivotal roles of lncRNAs in colorectal cancer progression. This study aimed to investigate the role of FGF14-AS2 in colorectal cancer development. FGF14-AS2 was found as a significantly downregulated lncRNA in TCGA dataset. Via RT-qPCR, we confirmed the downregulation of FGF14-AS2 in collected colorectal carcinoma samples. Transfection of plasmid containing full length of FGF14-AS2 repressed cell proliferation and induced elevation of cell apoptosis in colorectal cancer cells. In addition, FGF14-AS2 overexpression inactivated MAPK/ERK signaling in cells. Bioinformatic analysis and subsequent cell-based assays showed that FGF14-AS2 sponging miR-1288-3p, an oncogenic miRNA in colorectal cancer. RERG, the regulator of Ras/ERK pathway, was predicted and verified as target gene of miR-1288. Via downregulation of miR-1288, FGF14-AS2 elevated RERG expression in colorectal cancer cells. Rescue assays indicated that FGF14-AS2 relied on regulation of RERG to control cell proliferation and apoptosis in colorectal cancer. Taken together, the current study demonstrated FGF14-AS2 as a regulator of colorectal cancer development via downregulation of miR-1288-3p and inactivation of Ras/ERK signaling.

Selection of viral variants during persistent infection of insectivorous bat cells with Middle East respiratory syndrome coronavirus

Coronaviruses that cause severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) are speculated to have originated in bats. The mechanisms by which these viruses are maintained in individuals or populations of reservoir bats remain an enigma. Mathematical models have predicted long-term persistent infection with low levels of periodic shedding as a likely route for virus maintenance and spillover from bats. In this study, we tested the hypothesis that bat cells and MERS coronavirus (CoV) can co-exist in vitro. To test our hypothesis, we established a long-term coronavirus infection model of bat cells that are persistently infected with MERS-CoV. We infected cells from Eptesicus fuscus with MERS-CoV and maintained them in culture for at least 126 days. We characterized the persistently infected cells by detecting virus particles, protein and transcripts. Basal levels of type I interferon in the long-term infected bat cells were higher, relative to uninfected cells, and disrupting the interferon response in persistently infected bat cells increased virus replication. By sequencing the whole genome of MERS-CoV from persistently infected bat cells, we identified that bat cells repeatedly selected for viral variants that contained mutations in the viral open reading frame 5 (ORF5) protein. Furthermore, bat cells that were persistently infected with ΔORF5 MERS-CoV were resistant to superinfection by wildtype virus, likely due to reduced levels of the virus receptor, dipeptidyl peptidase 4 (DPP4) and higher basal levels of interferon in these cells. In summary, our study provides evidence for a model of coronavirus persistence in bats, along with the establishment of a unique persistently infected cell culture model to study MERS-CoV-bat interactions.

Multiple Loci Control Eyespot Number Variation on the Hindwings of Bicyclus anynana Butterflies

The underlying genetic changes that regulate the appearance and disappearance of repeated traits, or serial homologs, remain poorly understood. One hypothesis is that variation in genomic regions flanking master regulatory genes, also known as input-output genes, controls variation in trait number, making the locus of evolution almost predictable. Another hypothesis implicates genetic variation in up- or downstream loci of master control genes. Here, we use the butterfly Bicyclus anynana, a species that exhibits natural variation in eyespot number on the dorsal hindwing, to test these two hypotheses. We first estimated the heritability of dorsal hindwing eyespot number by breeding multiple butterfly families differing in eyespot number and regressing eyespot numbers of offspring on midparent values. We then estimated the number and identity of independent genetic loci contributing to eyespot number variation by performing a genome-wide association study with restriction site-associated DNA sequencing from multiple individuals varying in number of eyespots sampled across a freely breeding laboratory population. We found that dorsal hindwing eyespot number has a moderately high heritability of ∼0.50 and is characterized by a polygenic architecture. Previously identified genomic regions involved in eyespot development, and novel ones, display high association with dorsal hindwing eyespot number, suggesting that homolog number variation is likely determined by regulatory changes at multiple loci that build the trait, and not by variation at single master regulators or input-output genes.

Infant Formula Feeding Changes the Proliferative Status in Piglet Neonatal Mammary Glands Independently of Estrogen Signaling

BACKGROUND

Soy infant formula contains isoflavones, which are able to bind to and activate estrogen receptor (ER) pathways. The mammary gland is sensitive to estrogens, raising concern that the use of soy formulas may promote premature development.

OBJECTIVE

We aimed to determine if soy formula feeding increases mammary gland proliferation and differentiation in comparison to other infant postnatal diets.

METHODS

White-Dutch Landrace piglets aged 2 d received either sow milk (Sow), or were provided milk formula (Milk), soy formula (Soy), milk formula supplemented with 17-beta-estradiol (2 mg/(kg·d); M + E2), or milk formula supplemented with genistein (84 mg/L of diet; M + G) until day 21. Mammary gland proliferation and differentiation was assessed by histology, and real-time RT-PCR confirmation of differentially expressed genes identified by microarray analysis.

RESULTS

Mammary terminal end bud numbers were 19-31% greater in the Milk, Soy, and M + G groups relative to the Sow and M + E2, P <0.05. Microarray analysis identified differentially expressed genes between each formula-fed group relative to the Sow (±1.7-fold, P <0.05). Real-time RT-PCR confirmed 2- to 4-fold increases in mRNA transcripts of genes involved in cell proliferation, insulin-like growth factor 1 (IGF1), fibroblast growth factor 10 (FGF10), and fibroblast growth factor 18 (FGF18), in all groups relative to the Sow, P <0.05. In contrast, genes involved in cell differentiation and ductal morphogenesis, angiotensin II receptor type 2 (AGTR2), microtubule associated protein 1b (MAP1B), and kinesin family member 26b (KIF26B), were significantly upregulated by 2-, 4-, and 13-fold, respectively, in the M + E2 group. Additionally, mRNA expression of ER-specific gene targets, progesterone receptor (PGR), was increased by 12-fold, and amphiregulin (AREG) and Ras-like estrogen regulated growth inhibitor (RERG) expression by 1.5-fold in the M + E2 group, P <0.05. In the soy and M + G groups, mRNA expressions of fatty acid synthesis genes were increased 2- to 4-fold.

CONCLUSIONS

Our data indicate soy formula feeding does not promote ER-signaling in the piglet mammary gland. Infant formula feeding (milk- or soy-based) may initiate proliferative pathways independently of estrogenic signaling.

miR-532-5p promotes breast cancer proliferation and migration by targeting RERG

Aberrant expression of microRNAs (miRNAs/miRs) mediates the initiation and progression of breast cancer. Therefore, it is important to investigate the molecular mechanisms of miRNAs and their effects on breast cancer progression. In the present study, miR-532-5p was highly expressed in breast cancer tissues compared with normal tissues. In addition, expression of ras-related and estrogen-regulated growth inhibitor (RERG), a tumor suppressor in breast cancer, was negatively correlated with miR-532-5p expression. Inhibition of miR-532-5p significantly elevated RERG at both mRNA and protein levels and inactivated the mitogen-activated protein kinase (MAPK)/ERK signaling pathway. Overexpression of miR-532-5p decreased RERG expression and activated the MAPK/ERK signaling in breast cancer cell line MDA-MB-231. Bioinformatic analysis indicated that RERG 3'-untraslated region contained a putative binding site for miR-532-5p. Dual luciferase assay further validated RERG as a target gene of miR-532-5p. Notably, downregulation of miR-532-5p inhibited MDA-MB-231 cell proliferation and migration, which was partially attenuated upon RERG knockdown. In conclusion, the current study revealed an oncogenic role of miR-532-5p in breast cancer cells via direct targeting of RERG expression.

Transcriptome dynamics of long noncoding RNAs and transcription factors demarcate human neonatal, adult, and human mesenchymal stem cell-derived engineered cartilage

The engineering of a native-like articular cartilage (AC) is a long-standing objective that could serve the clinical needs of millions of patients suffering from osteoarthritis and cartilage injury. An incomplete understanding of the developmental stages of AC has contributed to limited success in this endeavor. Using next generation RNA sequencing, we have transcriptionally characterized two critical stages of AC development in humans-that is, immature neonatal and mature adult, as well as tissue-engineered cartilage derived from culture expanded human mesenchymal stem cells. We identified key transcription factors (TFs) and long noncoding RNAs (lncRNAs) as candidate drivers of the distinct phenotypes of these tissues. AGTR2, SCGB3A1, TFCP2L1, RORC, and TBX4 stand out as key TFs, whose expression may be capable of reprogramming engineered cartilage into a more expandable and neonatal-like cartilage primed for maturation into biomechanically competent cartilage. We also identified that the transcriptional profiles of many annotated but poorly studied lncRNAs were dramatically different between these cartilages, indicating that lncRNAs may also be playing significant roles in cartilage biology. Key neonatal-specific lncRNAs identified include AC092818.1, AC099560.1, and KC877982. Collectively, our results suggest that tissue-engineered cartilage can be optimized for future clinical applications by the specific expression of TFs and lncRNAs.

Identifying Methylation Pattern and Genes Associated with Breast Cancer Subtypes

Breast cancer is regarded worldwide as a severe human disease. Various genetic variations, including hereditary and somatic mutations, contribute to the initiation and progression of this disease. The diagnostic parameters of breast cancer are not limited to the conventional protein content and can include newly discovered genetic variants and even genetic modification patterns such as methylation and microRNA. In addition, breast cancer detection extends to detailed breast cancer stratifications to provide subtype-specific indications for further personalized treatment. One genome-wide expression–methylation quantitative trait loci analysis confirmed that different breast cancer subtypes have various methylation patterns. However, recognizing clinically applied (methylation) biomarkers is difficult due to the large number of differentially methylated genes. In this study, we attempted to re-screen a small group of functional biomarkers for the identification and distinction of different breast cancer subtypes with advanced machine learning methods. The findings may contribute to biomarker identification for different breast cancer subtypes and provide a new perspective for differential pathogenesis in breast cancer subtypes.

Integrated Analysis of Germline and Tumor DNA Identifies New Candidate Genes Involved in Familial Colorectal Cancer

Colorectal cancer (CRC) shows aggregation in some families but no alterations in the known hereditary CRC genes. We aimed to identify new candidate genes which are potentially involved in germline predisposition to familial CRC. An integrated analysis of germline and tumor whole-exome sequencing data was performed in 18 unrelated CRC families. Deleterious single nucleotide variants (SNV), short insertions and deletions (indels), copy number variants (CNVs) and loss of heterozygosity (LOH) were assessed as candidates for first germline or second somatic hits. Candidate tumor suppressor genes were selected when alterations were detected in both germline and somatic DNA, fulfilling Knudson's two-hit hypothesis. Somatic mutational profiling and signature analysis were also performed. A series of germline-somatic variant pairs were detected. In all cases, the first hit was presented as a rare SNV/indel, whereas the second hit was either a different SNV (3 genes) or LOH affecting the same gene (141 genes). BRCA2, BLM, ERCC2, RECQL, REV3L and RIF1 were among the most promising candidate genes for germline CRC predisposition. The identification of new candidate genes involved in familial CRC could be achieved by our integrated analysis. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

Genetic and genomic analyses of testicular hypoplasia in Nellore cattle

Reproductive performance is a key indicator of the long-term sustainability of any livestock production system. Testicular hypoplasia (TH) is a morphological and functional reproductive disorder that affects bulls around the world and consequently causes major economic losses due to reduced fertility rates. Despite the improvements in management practices to enhance performance of affected animals, the use of hypoplastic animals for reproduction might contribute to expand the prevalence of this disorder. The aim of this study was to identify genomic regions that are associated with TH in Nellore cattle by performing a genome-wide association study (GWAS) and functional analyses. Phenotypic and pedigree data from 47,563 animals and genotypes (500,689 Single Nucleotide Polymorphism, SNPs) from 265 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. The estimated breeding values (EBVs) were calculated by fitting a single-trait threshold animal model using a Bayesian approach. The SNP effects were estimated using the Bayes C method and de-regressed EBVs for TH as the response variable (pseudo-phenotype). The top-15 ranking windows (5-adjacent SNPs) that explained the highest proportion of variance were identified for further functional and biological network analyses. The posterior mean (95% highest posterior density) of the heritability for TH was 0.16 (0.08; 0.23). The most important genomic windows were located on BTA1, BTA3, BTA4, BTA5, BTA9, BTA22, BTA23, and BTA25. These windows explained together 22.69% of the total additive genetic variance for TH. Strong candidate genes associated with metabolism and synthesis of steroids, cell survival, spermatogenesis process and sperm motility were identified, which might play an important role in the expression of TH. Our findings contribute to a better biological understanding of TH and future characterization of causal variants might enable improved genomic prediction of this trait in beef cattle.

Site-specific differences in osteoblast phenotype, mechanical loading response and estrogen receptor-related gene expression

The osteoporosis-resistant nature of skull bones implies inherent differences exist between their cellular responses and those of other osteoporosis-susceptible skeletal sites. Phenotypic differences in calvarial and femoral osteoblastic responses to induction of osteogenesis, mechanical loading, estrogen, growth factor and cytokine stimulation were investigated. Primary rat calvarial and femoral adult male osteoblasts were cultured and osteoblastic mineralisation and maturation determined using Alizarin Red staining and expression of osteogenic marker genes assessed. Expression of known mechanically-responsive genes was compared between sites following loading of scaffold-seeded cells in a bioreactor. Cell proliferation and differentiation following growth factor and estrogen stimulation were also compared. Finally expression of estrogen receptors and associated genes during osteogenic differentiation were investigated. Calvarial osteoblasts exhibited delayed maturation (45d. vs 21d.) and produced less mineralised matrix than femoral osteoblasts when osteogenically induced. PDGF-BB and FGF2 both caused a selective increase in proliferation and decrease in osteoblastic differentiation of femoral osteoblasts. Mechanical stimulation resulted in the induction of the expression of Ccl2 and Anx2a selectively in femoral osteoblasts, but remained unchanged in calvarial cells. Estrogen receptor beta expression was selectively upregulated 2-fold in calvarial osteoblasts. Most interestingly, the estrogen responsive transcriptional repressor RERG was constitutively expressed at 1000-fold greater levels in calvarial compared with femoral osteoblasts. RERG expression in calvarial osteoblasts was down regulated during osteogenic induction whereas upregulation occurred in femoral osteoblasts. Bone cells of the skull are inherently different to those of the femur, and respond differentially to a range of stimuli. These site-specific differences may have important relevance in the development of strategies to tackle metabolic bone disorders.

White-nose syndrome is associated with increased replication of a naturally persisting coronaviruses in bats

Spillover of viruses from bats to other animals may be associated with increased contact between them, as well as increased shedding of viruses by bats. Here, we tested the prediction that little brown bats (Myotis lucifugus) co-infected with the M. lucifugus coronavirus (Myl-CoV) and with Pseudogymnoascus destructans (Pd), the fungus that causes bat white-nose syndrome (WNS), exhibit different disease severity, viral shedding and molecular responses than bats infected with only Myl-CoV or only P. destructans. We took advantage of the natural persistence of Myl-CoV in bats that were experimentally inoculated with P. destructans in a previous study. Here, we show that the intestines of virus-infected bats that were also infected with fungus contained on average 60-fold more viral RNA than bats with virus alone. Increased viral RNA in the intestines correlated with the severity of fungus-related pathology. Additionally, the intestines of bats infected with fungus exhibited different expression of mitogen-activated protein kinase pathway and cytokine related transcripts, irrespective of viral presence. Levels of coronavirus antibodies were also higher in fungal-infected bats. Our results suggest that the systemic effects of WNS may down-regulate anti-viral responses in bats persistently infected with M. lucifugus coronavirus and increase the potential of virus shedding.

Transcriptome analysis for identification of candidate genes related to sex determination and growth in Charybdis japonica

Charybdis japonica is an important cultured crab in China and it exhibits sex differences in their growth. Growth is an important economic trait that is controlled by many genes. In order to discover the growth-related regulatory mechanisms, whole-body transcriptomic sequencing and comparative analyses in different genders of C. japonica were conducted based on Illumina RNA-seq technology. In the present study, we obtained 62,879,204 and 60,226,334 clean reads in female and male libraries, respectively. 25,000,000 clean reads of every library were randomly selected and compared with Nt database to examine the possible contamination. Results showed that all clean reads were distributed among C. japonica or other species that were closely relative to this species, indicating no-pollution. De novo assembly was performed and a total of 32,543 and 44,174 unigenes were produced in female and male of C. japonica, respectively. Among all the unigenes, 12,591 and 14,455 unigenes of female and male crabs were annotated based on protein databases. Moreover, a total of 33,926 unigenes were found to contain ORFs and 52,839 SSRs were detected. The contrast between male and female C. japonica identifying 1939 unigenes were significantly differentially expressed. In addition, we specifically discussed some gene functions and pathways potentially associated with sex determination and growth. This is the first systematic report of whole transcriptome in C. japonica. The transcriptome information provides a basic resource for further studies on understanding the molecular basis of biological processes in C. japonica and other crustaceans.

miR-135 promotes proliferation and stemness of oesophageal squamous cell carcinoma by targeting RERG

MicroRNA (miRNA) plays an important role in tumourigenesis and cancer development by regulating oncogenes or tumour suppressor that are implicated in cell cycle, cell mobility and even cell senescence. Due to the resistance to enzymes that could degrade nucleotides, miRNAs have been considered proper for diagnosis and prognosis evaluation of cancer. The present study was designed to investigate miRNA associated with ESCC and identified effective miRNAs, which could serve as biomarker or targets. We first performed miRNA profiling to identify a subset of dysregulated miRNAs in ESCC. miR-135, miR-451 and miR-186 were the most differentially expressed miRNAs. Subsequent RT-PCR validated that miR-135 was upregulated in ESCC cell lines TE2 and TE9, implying the promise as a prognostic and diagnostic marker. The Cox regression analysis suggests the correlation of miR-135 expression and tumour stages. Survival analysis demonstrated metastatic samples largely have higher miR-135 expression. Downregulation of miR-135 suppressed proliferation and invasion of TE2 and TE9 cell lines. Subsequent target prediction combined with functional enrichment analysis identified "Small GTPase superfamily" that are possibly targeted by miR-135, which offers candidates for further investigation. Finally, RERG was identified as a target of miR-135. Downregulation of RERG could inhibit the cell proliferation and sphere formation ability of TE2 and TE9. Taken together, miR-135 was proved to promote tumour development of ESCC, which promises the prospect of using miR-135 as a biomarker indicator in diagnosis and prognosis.

In Vivo Study of the Effects of ERβ on Apoptosis and Proliferation of Hormone-Independent Prostate Cancer Cell Lines PC-3M

OBJECTIVE

To evaluate the in vivo therapeutic effects of attenuated Salmonella carrying PCDNA3.1-ERβ plasmid in hormone-independent prostatic cancer in nude mice and to clarify the mechanism by which estrogen receptor β (ERβ) induces apoptosis and proliferation in prostatic cancer cells in mice.

METHODS

The orthotopic prostatic cancer models of mice were randomly divided as follows: MOCK group, treated with PBS, PQ group, treated with attenuated Salmonella alone, PQ-PCDNA3.1 group, treated with attenuated Salmonella carrying PCDNA3.1 plasmid, and PQ-PCDNA3.1-ERβ group, treated with the attenuated Salmonella carrying PCDNA3.1-ERβ plasmid. Then, 10 μl of the plasmid-containing solution, comprising 1 × 107 cfu of the bacteria, was administered via intranasal delivery to each group except the MOCK group. The experimental methods included flow cytometry and terminal deoxyribonucleotidyl transferase-mediated dUTP-digoxigenin nick end-labelling (TUNEL) assay, immunohistochemistry, and western blotting.

RESULTS

Compared with the MOCK, PQ, and PQ-PCDNA3.1 groups, the weights of tumors in the PQ-PCDNA3.1-ERβ group were significantly reduced. The results of flow cytometry and TUNEL assay revealed that the number of apoptotic cells in the PQ-PCDNA3.1-ERβ group significantly increased. Compared with PQ-PCDNA3.1 group, the protein expression levels of ERβ, Bad, p-caspase 9, p-caspase 3, and cleaved PARP in the PQ-PCDNA3.1-ERβ group were significantly increased, while the expression levels of Akt, p-Akt, and Bcl-xl were decreased (P < 0.05).

CONCLUSION

The attenuated Salmonella carrying PCDNA3.1-ERβ plasmid could inhibit the growth of orthotopic prostatic cancer in mice by increasing the expression of ERβ.

Estrogen receptor α dependent regulation of estrogen related receptor β and its role in cell cycle in breast cancer

BACKGROUND

Breast cancer (BC) is highly heterogeneous with ~ 60-70% of estrogen receptor positive BC patient's response to anti-hormone therapy. Estrogen receptors (ERs) play an important role in breast cancer progression and treatment. Estrogen related receptors (ERRs) are a group of nuclear receptors which belong to orphan nuclear receptors, which have sequence homology with ERs and share target genes. Here, we investigated the possible role and clinicopathological importance of ERRβ in breast cancer.

METHODS

Estrogen related receptor β (ERRβ) expression was examined using tissue microarray slides (TMA) of Breast Carcinoma patients with adjacent normal by immunohistochemistry and in breast cancer cell lines. In order to investigate whether ERRβ is a direct target of ERα, we investigated the expression of ERRβ in short hairpin ribonucleic acid knockdown of ERα breast cancer cells by western blot, qRT-PCR and RT-PCR. We further confirmed the binding of ERα by electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP), Re-ChIP and luciferase assays. Fluorescence-activated cell sorting analysis (FACS) was performed to elucidate the role of ERRβ in cell cycle regulation. A Kaplan-Meier Survival analysis of GEO dataset was performed to correlate the expression of ERRβ with survival in breast cancer patients.

RESULTS

Tissue microarray (TMA) analysis showed that ERRβ is significantly down-regulated in breast carcinoma tissue samples compared to adjacent normal. ER + ve breast tumors and cell lines showed a significant expression of ERRβ compared to ER-ve tumors and cell lines. Estrogen treatment significantly induced the expression of ERRβ and it was ERα dependent. Mechanistic analyses indicate that ERα directly targets ERRβ through estrogen response element and ERRβ also mediates cell cycle regulation through p18, p21^cip and cyclin D1 in breast cancer cells. Our results also showed the up-regulation of ERRβ promoter activity in ectopically co-expressed ERα and ERRβ breast cancer cell lines. Fluorescence-activated cell sorting analysis (FACS) showed increased G0/G1 phase cell population in ERRβ overexpressed MCF7 cells. Furthermore, ERRβ expression was inversely correlated with overall survival in breast cancer. Collectively our results suggest cell cycle and tumor suppressor role of ERRβ in breast cancer cells which provide a potential avenue to target ERRβ signaling pathway in breast cancer.

CONCLUSION

Our results indicate that ERRβ is a negative regulator of cell cycle and a possible tumor suppressor in breast cancer. ERRβ could be therapeutic target for the treatment of breast cancer.

High N-Acetyltransferase 1 Expression Is Associated with Estrogen Receptor Expression in Breast Tumors, but Is not Under Direct Regulation by Estradiol, 5α-androstane-3β,17β-Diol, or Dihydrotestosterone in Breast Cancer Cells

N-acetyltransferase 1 (NAT1) is an enzyme that metabolizes carcinogens, which suggests a potential role in breast carcinogenesis. High NAT1 expression in breast tumors is associated with estrogen receptor α (ERα+) and the luminal subtype. We report that NAT1 mRNA transcript, protein, and enzyme activity were higher in human breast tumors with high expression of ERα/ESR1 compared with normal breast tissue. There was a strong correlation between NATb promoter and NAT1 protein expression/enzyme activity. High NAT1 expression in tumors was not the result of adipocytes, as evidenced by low perilipin (PLIN) expression. ESR1, NAT1, and XBP1 expression were associated in tumor biopsies. Direct regulation of NAT1 transcription by estradiol (E2) was investigated in ERα (+) MCF-7 and T47D breast cancer cells. E2 did not increase NAT1 transcript expression but increased progesterone receptor expression in a dose-dependent manner. Likewise, NAT1 transcript levels were not increased by dihydrotestosterone (DHT) or 5α-androstane-3β, (3β-adiol) 17β-diol. Dithiothreitol increased levels of the activated, spliced XBP1 in ERα (+) MCF-7 and T47D breast cancer cells but did not affect NAT1 or ESR1 expression. We conclude that NAT1 expression is not directly regulated by E2, DHT, 3β-adiol, or dithiothreitol despite high NAT1 and ESR1 expression in luminal A breast cancer cells, suggesting that ESR1, XBP1, and NAT1 expression may share a common transcriptional network arising from the luminal epithelium associated with better survival in breast cancer. Clusters of high-expression genes, including NAT1, in breast tumors might serve as potential targets for novel therapeutic drug development.

Mediator complex (MED) 7: a biomarker associated with good prognosis in invasive breast cancer, especially ER+ luminal subtypes

Background

Mediator complex (MED) proteins have a key role in transcriptional regulation, some interacting with the oestrogen receptor (ER). Interrogation of the METABRIC cohort suggested that MED7 may regulate lymphovascular invasion (LVI). Thus MED7 expression was assessed in large breast cancer (BC) cohorts to determine clinicopathological significance.

Methods

MED7 gene expression was investigated in the METABRIC cohort (n = 1980) and externally validated using bc-GenExMiner v4.0. Immunohistochemical expression was assessed in the Nottingham primary BC series (n = 1280). Associations with clinicopathological variables and patient outcome were evaluated.

Results

High MED7 mRNA and protein expression was associated with good prognostic factors: low grade, smaller tumour size, good NPI, positive hormone receptor status (p < 0.001), and negative LVI (p = 0.04) status. Higher MED7 protein expression was associated with improved BC-specific survival within the whole cohort and ER+/luminal subgroup. Pooled MED7 gene expression data in the external validation cohort confirmed association with better survival, corroborating with the protein expression. On multivariate analysis, MED7 protein was independently predictive of longer BC-specific survival in the whole cohort and Luminal A subtype (p < 0.001).

Conclusions

MED7 is an important prognostic marker in BC, particularly in ER+luminal subtypes, associated with improved survival and warrants future functional analysis.

Mutation status coupled with RNA-sequencing data can efficiently identify important non-significantly mutated genes serving as diagnostic biomarkers of endometrial cancer

Background

Endometrial cancers (ECs) are one of the most common types of malignant tumor in females. Substantial efforts had been made to identify significantly mutated genes (SMGs) in ECs and use them as biomarkers for the classification of histological subtypes and the prediction of clinical outcomes. However, the impact of non-significantly mutated genes (non-SMGs), which may also play important roles in the prognosis of EC patients, has not been extensively studied. Therefore, it is essential for the discovery of biomarkers in ECs to further investigate the non-SMGs that were highly associated with clinical outcomes.

Results

For the 9681 non-SMGs reported by the mutation annotation pipeline, there were 1053, 1273 and 395 non-SMGs differentially expressed between the patient groups divided by the clinical endpoints of histological grade, histological type as well as the International Federation of Gynecology and Obstetrics (FIGO) stage of ECs, respectively. In the gene set enrichment analysis, the cancer-related pathways, namely neuroactive ligand-receptor interaction signaling pathway, cAMP signaling pathway and calcium signaling pathway, were significantly enriched with the differentially expressed non-SMGs for all the three endpoints. We further identified 23, 19 and 24 non-SMGs, which were highly associated with histological grade, histological type and FIGO stage, respectively, from the differentially expressed non-SMGs by using the variable combination population analysis (VCPA) approach and found that 69.6% (16/23), 78.9% (15/19) and 66.7% (16/24) of the identified non-SMGs had been previously reported to be correlated with cancers. In addition, the averaged areas under the receiver operating characteristic curve (AUCs) achieved by the predictive models with identified non-SMGs as predictors in predicting histological type, histological grade, and FIGO stage were 0.993, 0.961 and 0.832, respectively, which were superior to those achieved by the models with SMGs as features (averaged AUCs = 0.928, 0.864 and 0.535, resp.).

Conclusions

Besides the SMGs, the non-SMGs reported in the mutation annotation analysis may also involve the crucial genes that were highly associated with clinical outcomes. Combining the mutation status with the gene expression profiles can efficiently identify the cancer-related non-SMGs as predictors for cancer prognostic prediction and provide more supplemental candidates for the discovery of biomarkers.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1891-6) contains supplementary material, which is available to authorized users.

Genome-wide association study identified copy number variants associated with sporadic colorectal cancer risk

Background

Multiple single nucleotide polymorphisms (SNPs) have been associated with colorectal cancer (CRC) risk. The role of structural or copy number variants (CNV) in CRC, however, remained unclear. We investigated the role of CNVs in patients with sporadic CRC.

Methods

A genome-wide association study (GWAS) was performed on 1000 Singapore Chinese patients aged 50 years or more with no family history of CRC and 1000 ethnicity-matched, age-matched and gender-matched healthy controls using the Affymetrix SNP 6 platform. After 16 principal component corrections, univariate and multivariate segmentations followed by association testing were performed on 1830 samples that passed quality assurance tests.

Results

A rare CNV region (CNVR) at chromosome 14q11 (OR=1.92 (95% CI 1.59 to 2.32), p=2.7e-12) encompassing CHD8, and common CNVR at chromosomes 3q13.12 (OR=1.54 (95% CI 1.33 to 1.77), p=2.9e-9) and 12p12.3 (OR=1.69 (95% CI 1.41 to 2.01), p=2.8e-9) encompassing CD47 and RERG/ARHGDIB, respectively, were significantly associated with CRC risk. CNV loci were validated in an independent replication panel using an optimised copy number assay. Whole-genome expression data in matched tumours of a subset of cases demonstrated that copy number loss at CHD8 was significantly associated with dysregulation of several genes that perturb the Wnt, TP53 and inflammatory pathways.

Conclusions

A rare CNVR at 14q11 encompassing the chromatin modifier CHD8 was significantly associated with sporadic CRC risk. Copy number loss at CHD8 altered expressions of genes implicated in colorectal tumourigenesis.

MicroRNA-382-5p aggravates breast cancer progression by regulating the RERG/Ras/ERK signaling axis

Aberrant activation of the Ras/ERK pathway mediates breast cancer initiation and aggressiveness. Therefore, it is important to identify miRNAs that modulate the Ras/ERK pathway during breast carcinogenesis and progression. The Ras GTPase superfamily member RERG (Ras-related and estrogen-regulated growth inhibitor) acts as a tumor suppressor to reduce breast cancer cell proliferation and tumor formation and has been suggested to have a regulatory role in the Ras/ERK pathway. In this study, we found that RERG exerted its tumor suppressor role by attenuating the activation of Ras/ERK signaling effectors. Furthermore, we found that miR-382-5p directly targets and represses RERG to attenuate the inhibitory effects of RERG on the oncogenic Ras/ERK pathway. Thereby, miR-382-5p promoted breast cancer cell viability, clonogenicity, survival, migration, invasion and in vivo tumorigenesis/metastasis. In clinical interpretation, miR-382-5p expression was negatively correlated with RERG expression, and it also significantly functioned as an independent oncomiR for the higher incidence and poorer prognosis of breast cancer. This novel connection highlights new diagnostic and prognostic roles for miR-382-5p and RERG in breast cancer.

Comprehensive Analysis of the Unfolded Protein Response in Breast Cancer Subtypes

Purpose

Triple-negative breast cancers (TNBCs) are associated with a worse prognosis and patients with TNBC have fewer therapeutic options than patients with non-TNBC. Recently, the IRE1α-XBP1 branch of the unfolded protein response (UPR) was implicated in TNBC prognosis on the basis of a relatively small patient population, suggesting the diagnostic and therapeutic value of this pathway in TNBCs. In addition, the IRE1α-XBP1 and hypoxia-induced factor 1 α (HIF1α) pathways have been identified as interacting partners in TNBC, suggesting a novel mechanism of regulation. To comprehensively evaluate and validate these findings, we investigated the relative activities and relevance to patient survival of the UPR and HIF1α pathways in different breast cancer subtypes in large populations of patients.

Materials and Methods

We performed a comprehensive analysis of gene expression and survival data from large cohorts of patients with breast cancer. The patients were stratified based on the average expression of the UPR or HIF1α gene signatures.

Results

We identified a strong positive association between the XBP1 gene signature and estrogen receptor-positive status or the HIF1α gene signature, as well as the predictive value of the XBP1 gene signature for survival of patients who are estrogen receptor negative, or have TNBC or HER2⁺. In contrast, another important UPR branch, the ATF4/CHOP pathway, lacks prognostic value in breast cancer in general. Activity of the HIF1α pathway is correlated with patient survival in all the subtypes evaluated.

Conclusion

These findings clarify the relevance of the UPR pathways in different breast cancer subtypes and underscore the potential therapeutic importance of the IRE1α-XBP1 branch in breast cancer treatment.

RERG suppresses cell proliferation, migration and angiogenesis through ERK/NF-κB signaling pathway in nasopharyngeal carcinoma

Background

Nasopharyngeal carcinoma (NPC) is a malignancy of the head and neck that is prevalent in Southeast Asia and southern China. Recent studies in epigenetics suggest that DNA methylation plays a pivotal role in the onset and progression of cancer. Combining the methyl-DNA binding domain capture technique and cDNA microarray analysis, we identified a unique hypermethylated gene, RERG (Ras-like estrogen-regulated growth inhibitor), that was down-regulated in NPC tissues. RERG is a tumor suppressor gene that was first reported in breast cancer. However, the functions of RERG are largely unknown in other tumor types.

Methods

RERG expression was assessed in human subjects (NPC primary tissues and non-cancer tissues) and cell lines (NPC cell lines and an immortalized epithelial cell line NP460). Further, we investigated the methylation rate of RERG in both human subject and cell lines. 5-Aza-2’-deoxycytidine (Aza) or combined with trichostatin A (TSA) were treated to three NPC cell lines (HK1, C666-1 and HK1_EBV). In addition, the role of RERG in NPC cells and its underlying mechanisms were explored by overexpression of RERG in NPC cell lines.

Results

RERG was significantly down-regulated in NPC cancer nests compared to normal nasopharyngeal epithelium cells. Furthermore, the RERG promoter was frequently methylated in NPC tissues and cell lines. The RERG methylation rate yielded an area under the curve (AUC) of receiver operating characteristic (ROC) curve was 0.897 (95%CI: 0.818–0.976). The down-regulation of RERG was restored in NPC cells treated with Aza and TSA. In addition, ectopic expression of RERG in NPC cell lines resulted in a significant suppression of cell proliferation, clonogenicity, migration and invasion. RERG-overexpressing cells showed significantly slower growth and less angiogenesis in tumor xenografts in nude mice. RERG suppressed the ERK/NF-κB signaling pathway and inhibited tumor growth and angiogenesis with down-regulation of MMPs and IL8 in tumors of nude mouse xenografts.

Conclusions

Our results suggest that RERG is frequently silenced by promoter CpG methylation in NPC, and acts as a functional tumor suppressor by suppressing the ERK/NF-κB signaling pathway. These findings support the potential use of RERG as a novel molecular target in NPC therapy.

Electronic supplementary material

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