Identification of upstream transcription factors (TFs) for expression signature genes in breast cancer

[Prokaryotic expression of a recombinant protein of adeno-associated virus capsid conserved regions and preparation of its polyclonal antibody]

OBJECTIVE

To express and purify the antigenic peptide of adeno-associated virus (AAV) capsid conserved regions in prokaryotic cells and prepare its rabbit polyclonal antibody.

METHODS

The DNA sequence encoding the conserved regions of AAV capsid protein was synthesized and cloned into the vector pET30a to obtain the plasmid pET30a-AAV-CR for prokaryotic expression and purification of the conserved peptides. Coomassie blue staining and Western blotting were used to identify the AAV conserved peptides. Japanese big ear white rabbits were immunized with AAV conserved region protein to prepare polyclonal antibody, with the rabbits injected with PBS as the control group. The antibody titer was determined with ELISA, and the performance of the antibody for recognizing capsid protein sequences of AAV1-AAV10 was assessed with Western blotting and immunofluorescence assay.

RESULTS

The plasmid pET30a-AAV-CR was successfully constructed, and a recombinant protein with a relative molecular mass of 17000 was obtained. The purified protein induced the production of antibodies against the conserved regions of AAV capsid in rabbits, and the titer of the purified antibodies reached 1:320 000. The antibodies were capable of recognizing a wide range of capsid protein sequences of AAV1-AAV10.

CONCLUSION

We successfully obtained the polyclonal antibodies against AAV capsid conserved region protein from rabbits, which facilitate future studies of AAV vector development and the biological functions of AAV.

ARID3A promotes the chemosensitivity of colon cancer by inhibiting AKR1C3

ARID3A is upregulated in colorectal cancer and can promote the proliferation and metastasis of cancer cells. However, patients with higher level of ARID3A have a better prognosis. This study aimed to uncover the mechanism by which ARID3A benefits the prognosis of colon cancer. Our results indicated that ARID3A upregulation enhanced the chemosensitivity of colon cancer cells to 5-Fluorouracil (5-FU), whereas ARID3A downregulation inhibited the chemosensitivity of colon cancer cells to 5-FU. Through database analysis, we found that AKR1C3, a drug resistance related gene, was the target of ARID3A. Moreover, AKR1C3 was downregulated in colon cancer tissues compared to normal tissues. Next, we assessed the interaction between AKR1C3 and ARID3A, and found that ARID3A inhibited the transcription of AKR1C3, leading to the downregulation of AKR1C3 in colon cancer cells. We also verified that AKR1C3 inhibited the chemosensitivity of colon cancer cells to 5-FU. Moreover, patients with higher ratio of ARID3A to AKR1C3 had a better prognosis. This study suggested that ARID3A promoted chemosensitivity of colon cancer cells by inhibiting AKR1C3 in colon cancer. The ratio of ARID3A to AKR1C3 is a good marker to predict the prognosis of colon cancer patients. This article is protected by copyright. All rights reserved.

Transcriptome profiling of lncRNA and co-expression network in the vaginal epithelial tissue of women with lubrication disorders

Background

Vaginal lubrication is a crucial physiological response that occurs at the beginning of sexual arousal. However, research on lubrication disorders (LD) is still in its infancy, and the role of long non-coding RNAs (lncRNAs) in LD remains unclear. This study aimed to explore the function of lncRNAs in the pathogenesis of vaginal LD.

Methods

The expression profiles of LD and normal control (NC) lncRNAs were examined using next-generation sequencing (NGS), and eight selected differentially expressed lncRNAs were verified by quantitative real-time PCR. We conducted GO annotation and KEGG pathway enrichment analyses to determine the principal functions of significantly deregulated genes. LncRNA-mRNA co-expression and protein-protein interaction (PPI) networks were constructed and the lncRNA transcription factors (TFs) were predicted.

Results

From the results, we identified 181,631 lncRNAs and 145,224 mRNAs in vaginal epithelial tissue. Subsequently, our preliminary judgment revealed a total of 499 up-regulated and 337 down-regulated lncRNAs in LD. The top three enriched GO items of the dysregulated lncRNAs included the following significant terms: “contractile fiber part,” “actin filament-based process,” and “contractile fiber”. The most enriched pathways were “cell-extracellular matrix interactions,” “muscle contraction,” “cell-cell communication,” and “cGMP-PKG signaling pathway”. Our results also showed that the lncRNA-mRNA co-expression network was a powerful platform for predicting lncRNA functions. We determined the three hub genes, ADCY5, CXCL12, and NMU, using PPI network construction and analysis. A total of 231 TFs were predicted with RHOXF1, SNAI2, ZNF354C and TBX15 were suspected to be involved in the mechanism of LD.

Conclusion

In this study, we constructed the lncRNA-mRNA co-expression network, predicted the lncRNA TFs, and comprehensively analyzed lncRNA expression profiles in LD, providing a basis for future studies on LD clinical biomarkers and therapeutic targets. Further research is also needed to fully determine lncRNA’s role in LD development.

Expression and prognostic values of ARID family members in breast cancer

The ARID family is a superfamily of 15 members containing a domain that interacts with AT-rich DNA elements. However, the expression and prognostic roles of each ARID in breast cancer are still elusive.

We used the TCGA and Kaplan-Meier plotter databases to assess the expression and prognostic values of ARID mRNA levels in breast cancer respectively. In the present study, 6 members were significantly lower in tumor tissues than those in the normal tissues, while 6 members were significantly higher. Further assessment of ARID expression in breast cancer with different molecular subtypes, 3 members were significantly higher in no-luminal molecular subtype than those in the luminal molecular subtype, and 6 members were significantly higher. In regard to prognostic values, high expression of ARID1A, ARID2, ARID3B, ARID4A, ARID5A, ARID5B, JARID1A were associated with favorable outcome, while ARID4B and JARID1B were correlated to a worse outcome. We further analyzed the prognostic value of ARID in different intrinsic subtypes and clinicopathological features of breast cancer. We found many meaningful ARID family biomarkers in breast cancer. The relevant results will expound the role of ARID in breast cancer and may further provide new insight to explore the ARID-targeting reagents for treating breast cancer patients.

Systems biology comprehensive analysis on breast cancer for identification of key gene modules and genes associated with TNM-based clinical stages

Breast cancer (BC), as one of the leading causes of death among women, comprises several subtypes with controversial and poor prognosis. Considering the TNM (tumor, lymph node, metastasis) based classification for staging of breast cancer, it is essential to diagnose the disease at early stages. The present study aims to take advantage of the systems biology approach on genome wide gene expression profiling datasets to identify the potential biomarkers involved at stage I, stage II, stage III, and stage IV as well as in the integrated group. Three HER2-negative breast cancer microarray datasets were retrieved from the GEO database, including normal, stage I, stage II, stage III, and stage IV samples. Additionally, one dataset was also extracted to test the developed predictive models trained on the three datasets. The analysis of gene expression profiles to identify differentially expressed genes (DEGs) was performed after preprocessing and normalization of data. Then, statistically significant prioritized DEGs were used to construct protein-protein interaction networks for the stages for module analysis and biomarker identification. Furthermore, the prioritized DEGs were used to determine the involved GO enrichment and KEGG signaling pathways at various stages of the breast cancer. The recurrence survival rate analysis of the identified gene biomarkers was conducted based on Kaplan-Meier methodology. Furthermore, the identified genes were validated not only by using several classification models but also through screening the experimental literature reports on the target genes. Fourteen (21 genes), nine (17 genes), eight (10 genes), four (7 genes), and six (8 genes) gene modules (total of 53 unique genes out of 63 genes with involving those with the same connectivity degree) were identified for stage I, stage II, stage III, stage IV, and the integrated group. Moreover, SMC4, FN1, FOS, JUN, and KIF11 and RACGAP1 genes with the highest connectivity degrees were in module 1 for abovementioned stages, respectively. The biological processes, cellular components, and molecular functions were demonstrated for outcomes of GO analysis and KEGG pathway assessment. Additionally, the Kaplan-Meier analysis revealed that 33 genes were found to be significant while considering the recurrence-free survival rate as an alternative to overall survival rate. Furthermore, the machine learning calcification models show good performance on the determined biomarkers. Moreover, the literature reports have confirmed all of the identified gene biomarkers for breast cancer. According to the literature evidence, the identified hub genes are highly correlated with HER2-negative breast cancer. The 53-mRNA signature might be a potential gene set for TNM based stages as well as possible therapeutics with potentially good performance in predicting and managing recurrence-free survival rates at stages I, II, III, and IV as well as in the integrated group. Moreover, the identified genes for the TNM-based stages can also be used as mRNA profile signatures to determine the current stage of the breast cancer.

Dysregulated lncRNA-miRNA-mRNA Network Reveals Patient Survival-Associated Modules and RNA Binding Proteins in Invasive Breast Carcinoma

Breast cancer is the most common cancer in women, but few biomarkers are effective in clinic. Previous studies have shown the important roles of non-coding RNAs in diagnosis, prognosis, and therapy selection for breast cancer and have suggested the significance of integrating molecules at different levels to interpret the mechanism of breast cancer. Here, we collected transcriptome data including long non-coding RNA (lncRNA), microRNA (miRNA), and mRNA for ~1,200 samples, including 1079 invasive breast carcinoma samples and 104 normal samples, from The Cancer Genome Atlas (TCGA) project. We identified differentially expressed lncRNAs, miRNAs, and mRNAs that distinguished invasive carcinoma samples from normal samples. We further constructed an integrated dysregulated network consisting of differentially expressed lncRNAs, miRNAs, and mRNAs and found housekeeping and cancer-related functions. Moreover, 58 RNA binding proteins (RBPs) involved in biological processes that are essential to maintain cell survival were found in the dysregulated network, and 10 were correlated with overall survival. In addition, we identified two modules that stratify patients into high- and low-risk subgroups. The expression patterns of these two modules were significantly different in invasive carcinoma versus normal samples, and some molecules were high-confidence biomarkers of breast cancer. Together, these data demonstrated an important clinical application for improving outcome prediction for invasive breast cancers.

Identification of long non-coding RNA and circular RNA in mice after intra-tracheal instillation with fine particulate matter

BACKGROUND

Fine particulate matter (PM_2.5) exposure has been proved to be associated with respiratory diseases in epidemiological studies, but the underlying mechanisms are not clear. One of the most important mechanisms involved is inflammation. Non-coding RNAs are proposed to play crucial roles in epigenetic modulation and post-transcriptional regulation. Identification of non-coding RNAs can show us the new insight into the molecular toxicity of PM_2.5.

MATERIALS AND METHODS

Intra-tracheal instillation of saline or PM_2.5 was performed in BALB/c Mice once a week for consecutive eight weeks. Genomewide transcriptome profiling of coding genes, long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) in mice lung were done by ribosomal RNA-depleted RNA sequencing. Lung histological alternations were observed in haematoxylin and eosin (HE) staining sections. The expressions of pro-inflammatory cytokines and Nucleotide-binding oligomerization domain-like receptor protein 3 (NLRP3) inflammasome were quantified by qRT-PCR、ELISA and Western blot.

RESULTS

1873 coding genes, 885 lncRNAs and 142 circRNAs were differentially expressed in lung tissues of the saline and PM_2.5 exposed mice. The upregulated expressions of lncRNA NONMMUT065867, lncRNA NONMMUT064312, lncRNA NONMMUT018123 and the downregulated expressions of circRNA CBT15_circR_1011, circRNA mm9_circ_005915 were identified by qRT-PCR in PM_2.5 group. The pulmonary inflammation score was higher in PM_2.5 group. What's more, the expressions of pro-inflammatory cytokines and NLRP3 inflammasome were upregulated in PM_2.5 exposed mice.

CONCLUSION

PM_2.5 causes lung inflammation and increases the expression of NLRP3 inflammasome. The identified novel lncRNAs and circRNAs may paly important role in the development of lung inflammation caused by PM_2.5.

Emerging Role of SOX Proteins in Breast Cancer Development and Maintenance

The SOX genes encode a family of more than 20 transcription factors that are critical regulators of embryogenesis and developmental processes and, when aberrantly expressed, have been shown to contribute to tumor development and progression in both an oncogenic and tumor suppressive role. Increasing evidence demonstrates that the SOX proteins play essential roles in multiple cellular processes that mediate or contribute to oncogenic transformation and tumor progression. In the context of breast cancer, SOX proteins function both as oncogenes and tumor suppressors and have been shown to be associated with tumor stage and grade and poor prognosis. Experimental evidence demonstrates that a subset of SOX proteins regulate critical aspects of breast cancer biology including cancer stemness and multiple signaling pathways leading to altered cell proliferation, survival, and tumor development; EMT, cell migration and metastasis; as well as other tumor associated characteristics. This review will summarize the role of SOX family members as important mediators of tumorigenesis in breast cancer, with an emphasis on the triple negative or basal-like subtype of breast cancer, as well as examine the therapeutic potential of these genes and their downstream targets.

hsa-miR-301a- and SOX10-dependent miRNA-TF-mRNA regulatory circuits in breast cancer

Breast cancer is the most common cancer among women and the molecular pathways that play main roles in breast cancer regulation are still not completely understood. MicroRNAs (miRNAs) and transcription factors (TFs) are important regulators of gene expression. It is important to unravel the relation of TFs, miRNAs, and their targets within regulatory networks to clarify the processes that cause breast cancer and the progression of it. In this study, mRNA and miRNA expression studies including breast tumors and normal samples were extracted from the GEO microarray database. Two independent mRNA studies and a miRNA study were selected and reanalyzed. Differentially expressed (DE) mRNAs and miRNAs between breast tumor and normal samples were listed by using BRBArray Tools. CircuitsDB2 analysis conducted with DE miRNAs and mRNAs resulted in 3 significant circuits that are SOX10- and hsamiR-301a-dependent. The following significant circuits were characterized and validated bioinformatically by using web-based tools: SOX10→hsa-miR-301a→HOXA3, SOX10→hsa-miR-301a→KIT, and SOX10→hsa-miR-301a→NFIB. It can be concluded that regulatory motifs involving miRNAs and TFs may be useful for understanding breast cancer regulation and for predicting new biomarkers.

Melatonin enhances the apoptotic effects and modulates the changes in gene expression induced by docetaxel in MCF‑7 human breast cancer cells

Results from clinical trials and multiple in vivo and in vitro studies point to melatonin as a promising adjuvant molecule with many beneficial effects when concomitantly administered with chemotherapy. Melatonin palliates side‑effects and enhances the efficacy of chemotherapeutic agents. However, the mechanisms through which melatonin regulates molecular changes induced by chemotherapeutic agents remain largely unknown. In this study, we demonstrated that melatonin enhanced the anti-proliferative and apoptotic responses to low doses of docetaxel in breast cancer cells. Importantly, these effects were more potent when melatonin was added prior to docetaxel. Treatment with 1 µM docetaxel (equivalent to the therapeutic dosage) induced changes in gene expression profiles and melatonin modulated these changes. Specifically, docetaxel downregulated TP53, cyclin-dependent kinase inhibitor 1A (CDKN1A) and cadherin 13 (CDH13), and upregulated mucin 1 (MUC1), GATA binding protein 3 (GATA3) and c-MYC, whereas melatonin counteracted these effects. Melatonin further stimulated the expression of the pro-apoptotic BAD and BAX genes, and enhanced the inhibition of the anti-apoptotic gene BCL-2 induced by docetaxel. The findings of this study suggest that melatonin is a molecule with potential for use as an adjuvant in cancer chemotherapy, which may have implications for designing clinical trials using chemotherapeutic drugs in combination with melatonin.

Identification and functional analysis of SOX10 phosphorylation sites in melanoma

The transcription factor SOX10 plays an important role in vertebrate neural crest development, including the establishment and maintenance of the melanocyte lineage. SOX10 is also highly expressed in melanoma tumors, and SOX10 expression increases with tumor progression. The suppression of SOX10 in melanoma cells activates TGF-β signaling and can promote resistance to BRAF and MEK inhibitors. Since resistance to BRAF/MEK inhibitors is seen in the majority of melanoma patients, there is an immediate need to assess the underlying biology that mediates resistance and to identify new targets for combinatorial therapeutic approaches. Previously, we demonstrated that SOX10 protein is required for tumor initiation, maintenance and survival. Here, we present data that support phosphorylation as a mechanism employed by melanoma cells to tightly regulate SOX10 expression. Mass spectrometry identified eight phosphorylation sites contained within SOX10, three of which (S24, S45 and T240) were selected for further analysis based on their location within predicted MAPK/CDK binding motifs. SOX10 mutations were generated at these phosphorylation sites to assess their impact on SOX10 protein function in melanoma cells, including transcriptional activation on target promoters, subcellular localization, and stability. These data further our understanding of SOX10 protein regulation and provide critical information for identification of molecular pathways that modulate SOX10 protein levels in melanoma, with the ultimate goal of discovering novel targets for more effective combinatorial therapeutic approaches for melanoma patients.

The nuclear receptor ERβ engages AGO2 in regulation of gene transcription, RNA splicing and RISC loading

BACKGROUND

The RNA-binding protein Argonaute 2 (AGO2) is a key effector of RNA-silencing pathways It exerts a pivotal role in microRNA maturation and activity and can modulate chromatin remodeling, transcriptional gene regulation and RNA splicing. Estrogen receptor beta (ERβ) is endowed with oncosuppressive activities, antagonizing hormone-induced carcinogenesis and inhibiting growth and oncogenic functions in luminal-like breast cancers (BCs), where its expression correlates with a better prognosis of the disease.

RESULTS

Applying interaction proteomics coupled to mass spectrometry to characterize nuclear factors cooperating with ERβ in gene regulation, we identify AGO2 as a novel partner of ERβ in human BC cells. ERβ-AGO2 association was confirmed in vitro and in vivo in both the nucleus and cytoplasm and is shown to be RNA-mediated. ChIP-Seq demonstrates AGO2 association with a large number of ERβ binding sites, and total and nascent RNA-Seq in ERβ + vs ERβ - cells, and before and after AGO2 knock-down in ERβ + cells, reveals a widespread involvement of this factor in ERβ-mediated regulation of gene transcription rate and RNA splicing. Moreover, isolation and sequencing by RIP-Seq of ERβ-associated long and small RNAs in the cytoplasm suggests involvement of the nuclear receptor in RISC loading, indicating that it may also be able to directly control mRNA translation efficiency and stability.

CONCLUSIONS

These results demonstrate that AGO2 can act as a pleiotropic functional partner of ERβ, indicating that both factors are endowed with multiple roles in the control of key cellular functions.