A novel correlation between ATP5A1 gene expression and progression of human clear cell renal cell carcinoma identified by co‑expression analysis

ATP5A1 as a potential prognostic biomarker in clear-cell renal cell carcinoma

Background

Clear cell renal cell carcinoma (ccRCC), a malignant neoplasm originating in the renal tubules, is characterized by extended treatment durations and suboptimal therapeutic outcomes in clinical settings. Adenosine triphosphate (ATP) synthase F1 subunit α (ATP5A1), a subunit of mitochondrial ATP synthase, is integral to the energy metabolism of specific tumors. While prior research has established a link between ATP5A1 expression and malignancies, its precise function and clinical significance in ccRCC are yet to be elucidated. The study aims to investigate the role of ATP5A1 in ccRCC and to explore the underlying molecular mechanisms.

Methods

The RNA sequencing data from ccRCC and corresponding adjacent tissues were analyzed through The Cancer Genome Atlas to evaluate their diagnostic and prognostic implications. ATP5A1 expression in ccRCC was validated using the Human Protein Atlas database. The role of ATP5A1 in ccRCC was further characterized through a series of assays, including wound healing, transwell invasion, cell counting kit-8 proliferation, and flow cytometry.

Results

ATP5A1 expression levels were elevated across 17 tumor types while being notably downregulated in 15 others, including ccRCC, esophageal carcinoma, and colon adenocarcinoma. Compared to 293 cells and adjacent normal kidney tissues, renal cancer cells and tissues exhibited a significant reduction in ATP5A1 expression. An inverse relationship was observed between ATP5A1 expression and both the clinical stage and histological grade of ccRCC, yet it is positively associated with improved prognosis. Silencing ATP5A1 expression enhanced the malignant biological properties of ccRCC, while its upregulation inhibited these effects. Furthermore, ATP5A1 knockdown activated the Wnt/β-catenin signaling pathway, whereas its overexpression resulted in pathway suppression.

Conclusions

Collectively, this study indicates that ATP5A1 may serve as a potential biomarker for the diagnosis, prognosis, and therapeutic targeting of ccRCC.

BRG1 promotes liver cancer cell proliferation and metastasis by enhancing mitochondrial function and ATP5A1 synthesis through TOMM40

Hepatocellular carcinoma (HCC) is one of the most lethal malignant tumors worldwide. Brahma-related gene 1 (BRG1), as a catalytic ATPase, is a major regulator of gene expression and is known to mutate and overexpress in HCC. The purpose of this study was to investigate the mechanism of action of BRG1 in HCC cells. In our study, BRG1 was silenced or overexpressed in human HCC cell lines. Transwell and wound healing assays were used to analyze cell invasiveness and migration. Mitochondrial membrane potential (MMP) and mitochondrial permeability transition pore (mPTP) detection were used to evaluate mitochondrial function in HCC cells. Colony formation and cell apoptosis assays were used to evaluate the effect of BRG1/TOMM40/ATP5A1 on HCC cell proliferation and apoptosis/death. Immunocytochemistry (ICC), immunofluorescence (IF) staining and western blot analysis were used to determine the effect of BRG1 on TOMM40, ATP5A1 pathway in HCC cells. As a result, knockdown of BRG1 significantly inhibited cell proliferation and invasion, promoted apoptosis in HCC cells, whereas BRG1 overexpression reversed the above effects. Overexpression of BRG1 can up-regulate MMP level, inhibit mPTP opening and activate TOMM40, ATP5A1 expression. Our results suggest that BRG1, as an oncogene, promotes HCC progression by regulating TOMM40 affecting mitochondrial function and ATP5A1 synthesis. Targeting BRG1 may represent a new and effective way to prevent HCC development.

Volumetric body composition analysis of the Cancer Genome Atlas reveals novel body composition traits and molecular markers Associated with Renal Carcinoma outcomes

Clinically, the body mass index remains the most frequently used metric of overall obesity, although it is flawed by its inability to account for different adipose (i.e., visceral, subcutaneous, and inter/intramuscular) compartments, as well as muscle mass. Numerous prior studies have demonstrated linkages between specific adipose or muscle compartments to outcomes of multiple diseases. Although there are no universally accepted standards for body composition measurement, many studies use a single slice at the L3 vertebral level. In this study, we use computed tomography (CT) studies from patients in The Cancer Genome Atlas (TCGA) to compare current L3-based techniques with volumetric techniques, demonstrating potential limitations with level-based approaches for assessing outcomes. In addition, we identify gene expression signatures in normal kidney that correlate with fat and muscle body composition traits that can be used to predict sex-specific outcomes in renal cell carcinoma.

HSPD1 Supports Osteosarcoma Progression through Stabilizing ATP5A1 and thus Activation of AKT/mTOR Signaling

Malignant transformation is concomitant with excessive activation of stress response pathways. Heat shock proteins (HSPs) are stress-inducible proteins that play a role in folding and processing proteins, contributing to the non-oncogene addiction of stressed tumor cells. However, the detailed role of the HSP family in osteosarcoma has not been investigated. Bulk and single-cell transcriptomic data from the GEO and TARGET databases were used to identify HSPs associated with prognosis in osteosarcoma patients. The expression level of HSPD1 was markedly increased in osteosarcoma, correlating with a negative prognosis. Through in vitro and in vivo experiments, we systematically identified HSPD1 as an important contributor to the regulation of proliferation, metastasis, and apoptosis in osteosarcoma by promoting the epithelial-mesenchymal transition (EMT) and activating AKT/mTOR signaling. Subsequently, ATP5A1 was determined as a potential target of HSPD1 using immunoprecipitation followed by mass spectrometry. Mechanistically, HSPD1 may interact with ATP5A1 to reduce the K48-linked ubiquitination and degradation of ATP5A1, which ultimately activates the AKT/mTOR pathway to ensure osteosarcoma progression and EMT process. These findings expand the potential mechanisms by which HSPD1 exerts biological effects and provide strong evidence for its inclusion as a potential therapeutic target in osteosarcoma.

[High expression of ATP5A1 in gastric carcinoma is correlated with a poor prognosis and enhanced glucose metabolism in tumor cells]

OBJECTIVE

To analyze the expression level of ATP5A1 in gastric carcinoma and its influence on the prognosis of the patients and glucose metabolism in the tumor cells.

METHODS

We retrospectively analyzed the data of 115 patients undergoing radical resection of gastric carcinoma in our hospital from February, 2013 to November, 2016. ATP5A1 expression in the surgical specimens were detected using immunohistochemistry, and the long-term prognosis of the patients with high (n=58) and low ATP5A1 expression (n=57) were analyzed. In gastric carcinoma MGC803 cells, the effects of lentivirus-mediated ATP5A1 knockdown or overexpression on glucose metabolism were investigated. We also observed the growth and glucose metabolism of xenografts derived from MGC803 cells with ATP5A1 knockdown or overexpression in nude mice.

RESULTS

ATP5A1 was significantly overexpressed in gastric carcinoma tissues in close correlation with blood CEA and CA19-9 levels, pathological grade, T stage and N stage (P < 0.05). ATP5A1 overexpression was an independent risk factor for a significantly lowered 5-year survival rate of patients with gastric carcinoma (P < 0.05). ROC curve analysis demonstrated the predictive value of high ATP5A1 expression for the patients'prognosis (P < 0.001). In MGC803 cells, ATP5A1 overexpression significantly upregulated cellular glucose uptake and lactate production and increased the protein levels of HK2, PFK1, and LDHA (P < 0.05), while ATP5A1 knockdown produced the opposite changes (P < 0.05). In the tumor-bearing mice, overexpression of ATP5A1 increased glucose metabolism of the tumor cells and promoted tumor growth (P < 0.05). Overexpression of ATP5A1 promoted the expressions of p-JNK and p-JUN in MGC803 cells (P < 0.05), and the JNK inhibitor SP600125 significantly inhibited the enhancement of cellular glucose metabolism induced by ATP5A1 overexpression (P < 0.05).

CONCLUSION

High ATP5A1 expression in gastric cancer is associated a poor long-term prognosis of the patients, and its effect is mediated at least partly by promoting glucose metabolism of the cells through the JNK/JUN pathway.

Identification of serum exosomal metabolomic and proteomic profiles for remote ischemic preconditioning

BACKGROUND

Remote ischemic preconditioning (RIPC) refers to a brief episode of exposure to potential adverse stimulation and prevents injury during subsequent exposure. RIPC has been shown to increase tolerance to ischemic injury and improve cerebral perfusion status. Exosomes have a variety of activities, such as remodeling the extracellular matrix and transmitting signals to other cells. This study aimed to investigate the potential molecular mechanism of RIPC-mediated neuroprotection.

METHODS

Sixty adult male military personnel participants were divided into the control group (n = 30) and the RIPC group (n = 30). We analyzed the differential metabolites and proteins in the serum exosomes of RIPC participants and control subjects.

RESULTS

Eighty-seven differentially expressed serum exosomal metabolites were found between the RIPC and control groups, which were enriched in pathways related to tyrosine metabolism, sphingolipid metabolism, serotonergic synapses, and multiple neurodegeneration diseases. In addition, there were 75 differentially expressed exosomal proteins between RIPC participants and controls, which involved the regulation of insulin-like growth factor (IGF) transport, neutrophil degranulation, vesicle-mediated transport, etc. Furthermore, we found differentially expressed theobromine, cyclo gly-pro, hemopexin (HPX), and apolipoprotein A1 (ApoA1), which are associated with neuroprotective benefits in ischemia/reperfusion injury. In addition, five potential metabolite biomarkers, including ethyl salicylate, ethionamide, piperic acid, 2, 6-di-tert-butyl-4-hydroxymethylphenol and zerumbone, that separated RIPC from control individuals were identified.

CONCLUSION

Our data suggest that serum exosomal metabolites are promising biomarkers for RIPC, and our results provide a rich dataset and framework for future analyses of cerebral ischemia‒reperfusion injury under ischemia/reperfusion conditions.

An in vitro approach to understand contribution of kidney cells to human urinary extracellular vesicles

Extracellular vesicles (EV) are membranous particles secreted by all cells and found in body fluids. Established EV contents include a variety of RNA species, proteins, lipids and metabolites that are considered to reflect the physiological status of their parental cells. However, to date, little is known about cell-type enriched EV cargo in complex EV mixtures, especially in urine. To test whether EV secretion from distinct human kidney cells in culture differ and can recapitulate findings in normal urine, we comprehensively analysed EV components, (particularly miRNAs, long RNAs and protein) from conditionally immortalised human kidney cell lines (podocyte, glomerular endothelial, mesangial and proximal tubular cells) and compared to EV secreted in human urine. EV from cell culture media derived from immortalised kidney cells were isolated by hydrostatic filtration dialysis (HFD) and characterised by electron microscopy (EM), nanoparticle tracking analysis (NTA) and Western blotting (WB). RNA was isolated from EV and subjected to miRNA and RNA sequencing and proteins were profiled by tandem mass tag proteomics. Representative sets of EV miRNAs, RNAs and proteins were detected in each cell type and compared to human urinary EV isolates (uEV), EV cargo database, kidney biopsy bulk RNA sequencing and proteomics, and single-cell transcriptomics. This revealed that a high proportion of the in vitro EV signatures were also found in in vivo datasets. Thus, highlighting the robustness of our in vitro model and showing that this approach enables the dissection of cell type specific EV cargo in biofluids and the potential identification of cell-type specific EV biomarkers of kidney disease.

Identification of miRNA-mRNA Regulatory Networks Associated with Diabetic Retinopathy using Bioinformatics Analysis

INTRODUCTION

Diabetic retinopathy (DR) is a major complication of diabetes and a leading cause of visual loss. This study aimed to explore biomarkers for DR that may provide additional reference to DR pathogenesis and development.

METHODS

The differentially expressed genes (DEGs) between the DR and control samples in the GSE53257 dataset were identified. Logistics analyses were performed to identify DR-associated miRNAs and genes, and correlation analysis was performed to determine the correlation between them in GSE160306.

RESULTS

A total of 114 DEGs in DR were identified in GSE53257. Three genes, including ATP5A1 (down), DAUFV2 (down), and OXA1L (down), were differentially expressed between DR and control samples in GSE160306. Univariate logistics analysis identified that ATP5A1 (OR=0.007, p = 1.40E-02), NDUFV2 (OR = 0.003, p = 6.40E-03), and OXA1L (OR = 0.093, p = 3.08E-02) were DR-associated genes. ATP5A1 and OXA1L were regulated by multiple miRNAs, of which hsa-let- 7b-5p (OR = 26.071, p = 4.40E-03) and hsa-miR-31-5p (OR = 4.188, p = 5.09E-02) were related to DR. ATP5A1 and OXA1L were closely correlated with each other in DR.

CONCLUSION

The hsa-miR-31-5p-ATP5A1 and hsa-let-7b-5p-OXA1L axes might play novel and important roles in the pathogenesis and development of DR.

Tissue-Based Proteomic Profiling in Patients with Hyperplasia and Endometrial Cancer

Uterine cancers are among the most prevalent gynecological malignancies, and endometrial cancer (EC) is the most common in this group. This study used tissue-based proteomic profiling analysis in patients with endometrial cancer and hyperplasia, and control patients. Conventional 2D gel electrophoresis, followed by a mass spectrometry approach with bioinformatics, including a network pathway analysis pipeline, was used to identify differentially expressed proteins and associated metabolic pathways between the study groups. Thirty-six patients (twelve with endometrial cancer, twelve with hyperplasia, and twelve controls) were enrolled in this study. The mean age of the participants was 46–75 years. Eighty-seven proteins were significantly differentially expressed between the study groups, of which fifty-three were significantly differentially regulated (twenty-eight upregulated and twenty-five downregulated) in the tissue samples of EC patients compared to the control (Ctrl). Furthermore, 26 proteins were significantly dysregulated (8 upregulated and 18 downregulated) in tissue samples of hyperplasia (HY) patients compared to Ctrl. Thirty-two proteins (nineteen upregulated and thirteen downregulated) including desmin, peptidyl prolyl cis-trans isomerase A, and zinc finger protein 844 were downregulated in the EC group compared to the HY group. Additionally, fructose bisphosphate aldolase A, alpha enolase, and keratin type 1 cytoskeletal 10 were upregulated in the EC group compared to those in the HY group. The proteins identified in this study were known to regulate cellular processes (36%), followed by biological regulation (16%). Ingenuity pathway analysis found that proteins that are differentially expressed between EC and HY are linked to AKT, ACTA2, and other signaling pathways. The panels of protein markers identified in this study could be used as potential biomarkers for distinguishing between EC and HY and early diagnosis and progression of EC from hyperplasia and normal patients.

Targeting of Ubiquitin E3 Ligase RNF5 as a Novel Therapeutic Strategy in Neuroectodermal Tumors

RNF5, an endoplasmic reticulum (ER) E3 ubiquitin ligase, participates to the ER-associated protein degradation guaranteeing the protein homeostasis. Depending on tumor model tested, RNF5 exerts pro- or anti-tumor activity. The aim of this study was to elucidate the controversial role of RNF5 in neuroblastoma and melanoma, two neuroectodermal tumors of infancy and adulthood, respectively. RNF5 gene levels are evaluated in publicly available datasets reporting the gene expression profile of melanoma and neuroblastoma primary tumors at diagnosis. The therapeutic effect of Analog-1, an RNF5 pharmacological activator, was investigated on in vitro and in vivo neuroblastoma and melanoma models. In both neuroblastoma and melanoma patients the high expression of RNF5 correlated with a better prognostic outcome. Treatment of neuroblastoma and melanoma cell lines with Analog-1 reduced cell viability by impairing the glutamine availability and energy metabolism through inhibition of F1Fo ATP-synthase activity. This latter event led to a marked increase in oxidative stress, which, in turn, caused cell death. Similarly, neuroblastoma- and melanoma-bearing mice treated with Analog-1 showed a significant delay of tumor growth in comparison to those treated with vehicle only. These findings validate RNF5 as an innovative drug target and support the development of Analog-1 in early phase clinical trials for neuroblastoma and melanoma patients.

Defueling the cancer: ATP synthase as an emerging target in cancer therapy

Reprogramming of cellular metabolism is a hallmark of cancer. Mitochondrial ATP synthase (MAS) produces most of the ATP that drives the cell. High expression of the MAS-composing proteins is found during cancer and is linked to a poor prognosis in glioblastoma, ovarian cancer, prostate cancer, breast cancer, and clear cell renal cell carcinoma. Cell surface-expressed ATP synthase, translocated from mitochondrion to cell membrane, involves the angiogenesis, tumorigenesis, and metastasis of cancer. ATP synthase has therefore been considered a therapeutic target. We review recent various ATP synthase inhibitors that suppress tumor growth and are being tested for the clinic.

Estimation of Gene Regulatory Networks from Cancer Transcriptomics Data

Cancer is a genetic disease in which multiple genes are perturbed. Thus, information about the regulatory relationships between genes is necessary for the identification of biomarkers and therapeutic targets. In this review, methods for inference of gene regulatory networks (GRNs) from transcriptomics data that are used in cancer research are introduced. The methods are classified into three categories according to the analysis model. The first category includes methods that use pair-wise measures between genes, including correlation coefficient and mutual information. The second category includes methods that determine the genetic regulatory relationship using multivariate measures, which consider the expression profiles of all genes concurrently. The third category includes methods using supervised and integrative approaches. The supervised approach estimates the regulatory relationship using a supervised learning method that constructs a regression or classification model for predicting whether there is a regulatory relationship between genes with input data of gene expression profiles and class labels of prior biological knowledge. The integrative method is an expansion of the supervised method and uses more data and biological knowledge for predicting the regulatory relationship. Furthermore, simulation and experimental validation of the estimated GRNs are also discussed in this review. This review identified that most GRN inference methods are not specific for cancer transcriptome data, and such methods are required for better understanding of cancer pathophysiology. In addition, more systematic methods for validation of the estimated GRNs need to be developed in the context of cancer biology.

Identification of TYROBP and FCER1G as Key Genes with Prognostic Value in Clear Cell Renal Cell Carcinoma by Bioinformatics Analysis

The involvement of aberrantly expressed genes in the pathogenesis and progression of various human malignancies has been widely reported, including clear cell renal cell carcinoma (ccRCC). This study aimed to identify potential crucial genes in ccRCC and further investigate the role of these genes in ccRCC prognosis. Three gene expression profiles (GSE3, GSE6344 and GSE53000) were downloaded from GEO database. GEO2R was performed to identify the differentially expressed genes (DEGs) between ccRCC and normal samples. GO analysis and KEGG pathway enrichment analysis were applied for the function analysis. The DEGs were mapped into the PPI network, then the hub genes were identified and verified using the ONCOMINE database. Kaplan-Meier plotter was used to evaluate of the prognostic value of the identified hub genes. A total of 113 DEGs were identified from the three gene expression profiles, including 64 up-regulated genes and 69 down-regulated genes. DEGs were observed to be enriched in biological processes related to the progress and pathogenesis of human cancers. According to PPI network, 5 hub genes were collected, including TYROBP, C1QB, ITGB2, CD53 and FCER1G. Among them, CD53 was newly identified, and Kaplan-Meier survival curves suggested that high expression of CD53 was significantly associated with poor survival in ccRCC patients (log-rank P < 0.01). The present results may provide new insight into the understanding of molecular mechanisms and the clinical prognosis of ccRCC.

ATP5A1 Participates in Transcriptional and Posttranscriptional Regulation of Cancer-Associated Genes by Modulating Their Expression and Alternative Splicing Profiles in HeLa Cells

Background: Aberrant expression and alternative splicing of oncogenes are the driving events in tumor initiation and development. But how these events are regulated in cancer cells is largely unknown. Functions of ATP5A1, an important mitochondrial ATP synthase gene, in transcriptional and posttranscriptional regulation were explored in this study. Methods: ATP5A1 was overexpressed using plasmid-transformed HeLa cells, and its influence on cell apoptosis and proliferation is evaluated. Transcriptome sequencing was then performed using RNA-seq to study the changes in gene expression and regulation of alternative splicing events. Validation of the implicated genes was achieved using RT-qPCR analysis. Results: It was found that ATP5A1 could significantly promote cellular apoptosis, but it had no influence on cell proliferation. ATP5A1 overexpression significantly increased the expression levels of genes associated with the innate immune response, angiogenesis, and collagen catabolic processes. This included enrichment of MMP2 and MMP19. It was also found that ATP5A1 could interfere with the alternative splicing of hundreds of genes associated with glucose homeostasis, HIF-1 signaling activation, and several pathways associated with cancers. Eight ATP5A1-regulated differentially expressed genes and 3 genes altered by splicing were selected and validated using RT-qPCR analysis. Conclusions: In summary, we illustrate the regulatory functions of ATP5A1 on the transcriptome of HeLa cells by exploring its influence on gene expression and alternative splicing. The results suggest that ATP5A1 may play an important regulatory role in cervical cancer cells by regulating expression and alternative splicing of cancer-associated genes. This study provides novel insights into the current understanding of the mechanisms of ATP5A1 on carcinogenesis and cancer progression.

Up-regulated cytotrophoblast DOCK4 contributes to over-invasion in placenta accreta spectrum

The syndrome of cytotrophoblast invasion beyond the normal boundary (in the superficial myometrium) is collectively termed placenta accreta spectrum. The incidence of this condition is rising. However, little is known about the underlying molecular changes. Global transcriptomic profiling of cytotrophoblasts isolated from these cases, as compared to gestational age-matched controls, revealed numerous changes in gene expression involving diverse pathways, including cell signaling, migration, and immune functions. DOCK4 was the most highly up-regulated mRNA in the cases. Mutations in this gene are mechanistically linked to cancer progression. Overexpression of DOCK4 in primary cytotrophoblasts increased their invasiveness. This study provides molecular insights into the pathways driving placenta accreta spectrum and suggests numerous future directions.

In humans, a subset of placental cytotrophoblasts (CTBs) invades the uterus and its vasculature, anchoring the pregnancy and ensuring adequate blood flow to the fetus. Appropriate depth is critical. Shallow invasion increases the risk of pregnancy complications, e.g., severe preeclampsia. Overly deep invasion, the hallmark of placenta accreta spectrum (PAS), increases the risk of preterm delivery, hemorrhage, and death. Previously a rare condition, the incidence of PAS has increased to 1:731 pregnancies, likely due to the rise in uterine surgeries (e.g., Cesarean sections). CTBs track along scars deep into the myometrium and beyond. Here we compared the global gene expression patterns of CTBs from PAS cases to gestational age-matched control cells that invaded to the normal depth from preterm birth (PTB) deliveries. The messenger RNA (mRNA) encoding the guanine nucleotide exchange factor, DOCK4, mutations of which promote cancer cell invasion and angiogenesis, was the most highly up-regulated molecule in PAS samples. Overexpression of DOCK4 increased CTB invasiveness, consistent with the PAS phenotype. Also, this analysis identified other genes with significantly altered expression in this disorder, potential biomarkers. These data suggest that CTBs from PAS cases up-regulate a cancer-like proinvasion mechanism, suggesting molecular as well as phenotypic similarities in the two pathologies.

Integrative Analysis of Siglec-15 mRNA in Human Cancers Based on Data Mining

Objective: Cancer is expected to be the leading cause of death worldwide within the 21st century and is the single most important obstacle to extending life expectancy. Unfortunately, the most effective approach to combating cancers remains a complex and unsolved problem. Siglec-15 is a member of the Siglec family and plays a conserved regulatory role in the immune system of vertebrates. Previous studies on Siglec-15 have focused on its function in osteoclast regulation. The purpose of this study was to explore the significance of Siglec-15 mRNA in human cancer mainly based on information obtained from online databases. Method: Data were collected from several online databases. Serial analysis of gene expression (SAGE) and Virtual Northern, UALCAN Database Analysis, Catalog of Somatic Mutations in Cancer (COSMIC) analysis, the cBio cancer genomics portal, Cancer Regulome tools and data, Kaplan-Meier Plotter Analysis and the UCSC Xena website were used to analyze the data. Results: Compared with normal tissues, Siglec-15 up-regulation was widely observed in tuomrs. Differences in Siglec-15 expression were associated with different prognoses. Siglec-15 mutations are widely observed in tumors and interact with different genes in different cancer types. Conclusion: Siglec-15 is a potential target for the expansion of cancer immunotherapy.

Identification of ATP synthase α subunit as a new maternal factor capable of protecting zebrafish embryos from bacterial infection

Previous studies have shown that ATP synthase α subunit (ATP5A1) plays multiple roles, but our understanding of its biologic functions remains poor and incomprehensive. Here, we clearly demonstrated that zebrafish ATP5A1 was a newly characterized lipoteichoic acid (LTA)- and LPS-binding protein abundantly stored in the eggs and embryos of zebrafish. Zebrafish ATP5A1 acted not only as a pattern recognition receptor, capable of identifying LTA and LPS, but also as an effector molecule, capable of inhibiting the growth of both gram-positive and -negative bacteria. ATP5A1 could disrupt the bacterial membranes by a combined action of membrane depolarization and permeabilization. We also found that the N-terminal 65 residues were critical for the antibacterial activity of zebrafish ATP5A1. In particular, we showed that microinjection of exogenous recombinant (r)ATP5A1 into early embryos could promote their resistance against pathogenic Aeromonas hydrophila challenge, and this pathogen-resistant activity was markedly reduced by the coinjection of anti-ATP5A1 antibody or by the knockdown with morpholino for atp5a1 but not by the coinjection of anti-actin antibody. Moreover, each egg/embryo contains a sufficient amount of ATP5A1 in vivo to kill A. hydrophila. Furthermore, the N-terminal 65 residues 1-65 of ATP5A1 α subunit (rA_1-65) with in vitro antibacterial activity also promoted the resistance of embryos against A. hydrophila, but the N-terminal 69 residues 66-134 (rA_66-134) or C-terminal residues 135-551 (rA_135-551) of ATP5A1 α subunit without in vitro antibacterial activity did not. Finally, we showed that the antibacterial activity of the N-terminal 65 residues of ATP5A1 α subunit was conserved throughout animal evolution. Collectively, these results indicate that ATP5A1 is a novel maternal immunocompetent factor that can protect the early embryos of zebrafish from bacterial infection. This work also provides a new viewpoint for understanding the biologic roles of ATP5A1, which is ubiquitously present in animals.-Ni, S., Zhou, Y., Chen, Y., Du, X., Zhang, S. Identification of ATP synthase α subunit as a new maternal factor capable of protecting zebrafish embryos from bacterial infection.

PPWD1 is associated with the occurrence of postmenopausal osteoporosis as determined by weighted gene co‑expression network analysis

Postmenopausal osteoporosis (PMO) is the most common type of primary osteoporosis (OP), a systemic skeletal disease. Although many factors have been revealed to contribute to the occurrence of PMO, specific biomarkers for the early diagnosis and therapy of PMO are not available. In the present study, a weighted gene co-expression network analysis (WGCNA) was performed to screen gene modules associated with menopausal status. The turquoise module was verified as the clinically significant module, and 12 genes (NUP133, PSMD12, PPWD1, RBM8A, CRNKL1, PPP2R5C, RBM22, PIK3CB, SKIV2L2, PAPOLA, SRSF1 and COPS2) were identified as ‘real’ hub genes in both the protein-protein interaction (PPI) network and co-expression network. Furthermore, gene expression analysis by microarray in blood monocytes from pre- and post-menopausal women revealed an increase in the expression of these hub genes in postmenopausal women. However, only the expression of peptidylprolyl isomerase domain and WD repeat containing 1 (PPWD1) was correlated with bone mineral density (BMD) in postmenopausal women. In the validation set, a similar expression pattern of PPWD1 was revealed. Functional enrichment analysis revealed that the fatty acid metabolism pathway was significantly abundant in the samples that exhibited a higher expression of PPWD1. Collectively, PPWD1 is indicated as a potential diagnostic biomarker for the occurrence of PMO.

Weighted correlation network analysis of triple-negative breast cancer progression: Identifying specific modules and hub genes based on the GEO and TCGA database

Triple-negative breast cancer (TNBC) represents an aggressive malignancy of frequent high histologic grade with no effective specific targeted therapies. The present study aimed to identify specific modules and hub genes that may influence the progression of TNBC. The key words ‘breast cancer’ were used to search microarray datasets in the Gene Expression Omnibus and The Cancer Genome Atlas databases that included 5 datasets. A total of 11 co-expression modules were constructed based on the expression levels of 5,782 genes obtained from 456 patients with TNBC using the weighted correlation network analysis (WGCNA). The results demonstrated that the red module was significantly associated with relapse-free survival (RFS) in patients with TNBC [hazard ratio (HR)=0.381, 95% confidence interval (CI), 0.183–0.793; P=0.010]. The functional enrichment analysis revealed that the biological processes corresponding to the red module were ‘mRNA processing’, ‘histone lysine methylation’ and ‘regulation of TOR signaling’. In addition, Hedgehog signaling pathways were considered to serve a critical role in the development of this disease (P<0.001). A total of 12 hub genes were identified, of which α-thalassemia/mental retardation syndrome X-linked (ATRX) was significantly associated with RFS in patients with TNBC (HR=0.601; 95%CI, 0.376–0.960; P=0.033). The receiver operating characteristic curve indicated that ATRX could distinguish relapse from non-relapse in patients with TNBC (area under the curve=0.570; P=0.023). In conclusion, the present study demonstrated that ATRX was associated with TNBC progression, which suggested that ATRX may be involved in a recombination-mediated telomere maintenance mechanism.

Analysis of gene co‑expression network reveals prognostic significance of CNFN in patients with head and neck cancer

In patients with head and neck cancer (HNC), lymph node (N) metastases are associated with cancer aggressiveness and poor prognosis. Identifying meaningful gene modules and representative biomarkers relevant to the N stage helps predict prognosis and reveal mechanisms underlying tumor progression. The present study used a step-wise approach for weighted gene co-expression network analysis (WGCNA). Dataset GSE65858 was subjected to WGCNA. RNA sequencing data of HNC downloaded from the Cancer Genome Atlas (TCGA) and dataset GSE39366 were utilized to validate the results. Following data preprocessing, 4,295 genes were screened, and blue and black modules associated with the N stage of HNC were identified. A total of 16 genes [keratinocyte differentiation associated protein, suprabasin, cornifelin (CNFN), small proline rich protein 1B, desmoglein 1 (DSG1), chromosome 10 open reading frame 99, keratin 16 pseudogene 3, gap junction protein β2, dermokine, LY6/PLAUR domain containing 3, transmembrane protein 79, phospholipase A2 group IVE, transglutaminase 5, potassium two pore domain channel subfamily K member 6, involucrin, kallikrein related peptidase 8] that had a negative association with the N-stage in the blue module, and two genes (structural maintenance of chromosomes 4 and mutS homolog 6) that had a positive association in the black module, were identified to be candidate hub genes. Following further validation in TCGA and dataset GSE65858, it was identified that CNFN and DSG1 were associated with the clinical stage of HNC. Survival analysis of CNFN and DSG1 was subsequently performed. Patients with increased expression of CNFN displayed better survival probability in dataset GSE65858 and TCGA. Therefore, CNFN was selected as the hub gene for further verification in the Gene Expression Profiling Interactive Analysis database. Finally, functional enrichment and gene set enrichment analyses were performed using datasets GSE65858 and GSE39366. Three gene sets, namely ‘P53 pathway’, ‘estrogen response early’ and ‘estrogen response late’, were enriched in the two datasets. In conclusion, CNFN, identified via the WGCNA algorithm, may contribute to the prediction of lymph node metastases and prognosis, probably by regulating the pathways associated with P53, and the early and late estrogen response.

A brief procedure for big data analysis of gene expression

There are a lot of biological and experimental data from genomics, proteomics, drug screening, medicinal chemistry, etc. A large amount of data must be analyzed by special methods of statistics, bioinformatics, and computer science. Big data analysis is an effective way to build scientific hypothesis and explore internal mechanism. Here, gene expression is taken as an example to illustrate the basic procedure of the big data analysis.

First chloroplast genomics study of Phoenix dactylifera (var. Naghal and Khanezi): A comparative analysis

Date palm (Phoenix dactylifera L.) is one of the oldest fruit crops in the arid regions of the Middle East. However, little information is available regarding its plastid genomes. In this study, we sequenced the chloroplast (cp) genomes of two economically important but genomically unexplored date palm cultivars of Phoenix dactylifera (var. Naghal and Khanezi). The data assembly and genome annotation revealed a typical quadripartite structure similar to Arecaceae, and the genome sizes of Naghal and Khanezi were 158,210 bp and 158,211 bp, respectively. Structurally, both cp genomes were comprised of four regions: a pair of inverted repeats (27,273 bp for Khanezi and for Naghal 27,272 bp), a large single-copy region (86,090 bp and 86,092 bp) and a small single-copy region (17,575 bp and 17,574 bp). Both genomes had 138 representative genes, whereas 227 and 229 randomly distributed microsatellites were also observed in Khanezi and Naghal, respectively. Phylogenetic analysis based on the whole cp genomes and 68 shared genes showed identical phylogenetic trees of Khanezi and Naghal forming clades with Khalas and Aseel cultivars, respectively. The current study showed detailed comparative cp genome analysis, which could be essential for broader population genetics and molecular studies of these four date palm cultivars.