
Citation: | Solomon E. Owumi, Uche O. Arunsi, Moses T. Otunla, Imisioluwa O. Oluwasuji. Exposure to lead and dietary furan intake aggravates hypothalamus-pituitary-testicular axis toxicity in chronic experimental rats[J]. The Journal of Biomedical Research, 2023, 37(2): 100-114. DOI: 10.7555/JBR.36.20220108 |
Nowadays, colorectal cancer (CRC) is one of the most deadly cancers and almost 900 000 CRC-related deaths were reported each year in the world[1]. With the understanding of pathophysiology of the disease, different treatment options to improve the survival rates of CRC patients have been developed in the world. The 5-year survival rate of CRC patients was >90% when the patients were diagnosed at early stages [2]. However, due to lacking early detection methods, many CRC patients were diagnosed at an advanced stage or in the metastasis status. And the 5-year survival rate for those diagnosed with metastasis was at approximately 12%[3]. Recently, a new kind of analysis method has been used to identify the differential expression genes between CRC and normal tissues based on the high-throughput sequencing platforms, such as microarrays. This is a promising tool with extensive clinical applications, including molecular diagnosis, prognosis prediction, new drug targets discovery, etc[4–6]. Furthermore, microarray assay combining bioinformatics analysis made it possible to analyze the gene expression on mRNA level in CRC progression. For example, several studies have used this method to identify key genes in CRC development through comparing with normal samples, and showed that the key genes were involved in different signal pathways, biological processes, and molecular functions[7–12].
However, with a relatively limited degree of overlap, we still can not find reliable biomarkers or drug targets. Therefore, the discovery of novel biomarkers for early detection and prognosis prediction of CRC is urgently required.
In the present study, we targeted to find key genes to develop novel biomarkers or drug targets for CRC. Therefore, we chose four Gene Expression Omnibus (GEO) datasets, GSE113513, GSE21510, GSE44076, and GSE32323, and used bioinformatics methods to screen the significant differentially expressed genes (DEGs) between CRC tissues and normal tissues. Gene ontology (GO) and KEGG Pathway analyses were used to find the biological roles of these DEGs through DAVID database. Furthermore, the PPI network of DEGs was constructed and key modules or hub genes were selected with Molecular Complex Detection (MCODE) plugin of Cytoscape software. And the clinical significance was validated by GEPIA database. Finally, small active candidate molecules were identified to develop new drugs through Connectivity Map (CMap) database. In brief, we found four hub genes involved in CRC, which may provide novel potential biomarkers for CRC prognosis, and two potential candidate drugs for CRC.
To explore the differential gene expression profiles between CRC and normal tissues, we searched the NCBI-GEO database to collect enough and adequate tissues. A total of 4 GEO datasets were selected, including GSE113513, GSE21510, GSE44076, and GSE32323. These mRNA profiles were based on platform GPL15207 (GSE113513), GPL570 (GSE21510 and GSE32323), and GPL13667 (GSE44076). A total of 253 CRC samples and 203 normal samples were chosen for this study, including 14 pairs of cancer and normal samples in GSE113513, 124 CRC samples and 24 normal samples in GSE21510, 98 pairs of cancer and normal samples plus 50 healthy donor tissues in GSE44076, and 17 pairs of cancer and normal samples in GSE32323.
To identify the DEGs, we used the NCBI-GEO2R online tool to analyze these datasets. Subsequently, adjusted P-value <0.05 and |log 2(fold change)| >1 were set as the cutoff criteria to screen the significant DEGs of each dataset. Finally, Venn diagrams were performed to get the overlap significant DEGs of the 4 datasets.
To find the biological functional roles of DEGs, GO and KEGG pathway analyses were performed through DAVID database. Significant results of molecular function (MF), biological process (BP), cellular component (CC), and biological pathways were selected with P-value <0.05.
The DEGs profiles were submitted to STRING database for exploring their potential interactions. The interactions with a combined score >0.4 were considered significant. Subsequently, the interaction files were downloaded and imported into Cytoscape software to construct the PPI network. The MCODE plugin was used to find key modules of the whole PPI network with a degree cutoff=2, node score cutoff=0.2, K-core=2, and max depth=100. The hub genes were then selected with connectivity degree >10. Furthermore, KEGG pathway analyses of the significant modules were performed with P-value <0.05.
To verify the hub genes we found, we used GEPIA database to analyze their expression and clinical prognostic information in 270 CRC patients. And the survival curve, stage analysis and box plot were performed to show the clinical implications of hub genes.
To find potential small active molecules to develop new drugs for treating CRC, we uploaded DEGs probe profiles into the CMap database. This database can help to predict small molecules that induce or reverse gene expression signature with a score from −1 to 1. And the molecules which value from 0 closer to −1 were functioned as reversing the cancer cell status.
Analyzed with the GEO2R online tool, a total of 1763, 4411, 2428, and 2276 DEGs were extracted from GSE113513, GSE21510, GSE44076, and GSE32323, respectively, using adjusted P-value <0.05 and |log 2(fold change)| >1 as cutoff criteria. The volcano plots of DEGs in each dataset were shown in Fig. 1A . And the Venn diagrams showed that 865 overlap DEGs were identified from these four datasets, including 374 significantly upregulated genes and 491 downregulated genes (Fig. 1B and Supplementary Table 1 , available online).
To explore the biological functional roles of the overlap DEGs, GO and KEGG analyses were performed on DAVID database. And the top 20 terms were listed in the charts (Fig. 2A –D and Supplementary Table 2 , available online). The GO analysis results consist of three functional categories, including BP, CC, and MF. In the BP group, DEGs were mainly enriched in cell proliferation (Fig. 2A ). In the CC group, DEGs were enriched in cytoplasm (Fig. 2B ). And in the MF group, DEGs were enriched in protein binding (Fig. 2C ). KEGG pathway analysis showed that DEGs were enriched in metabolic pathways, pathways in cancer and cell cycle (Fig. 2D ). The details of the top 20 terms were listed in Supplementary Table 2 .
Using the STRING online database and Cytoscape software, a total of 865 DEGs were filtered into the PPI network complex, containing 863 nodes and 5817 edges (Fig. 3 ). Based on degree scores using the MCODE plugin, two key modules were detected from the whole PPI network complex. Module 1 contained 61 nodes and 1648 edges, and DEGs were enriched in cell cycle, oocyte meiosis, progesterone-mediated oocyte maturation, DNA replication and p53 signaling pathway (Fig. 4A and B ). Module 2 had 55 nodes and 625 edges, and these DEGs were enriched in chemokine signaling pathway, ribosome biogenesis in eukaryotes, cytokine-cytokine receptor interaction, pathways in cancer, purine metabolism, RNA polymerase, retrograde endocannabinoid signaling, TNF signaling pathway, legionellosis, regulation of lipolysis in adipocytes, NOD-like receptor signaling pathway, cytosolic DNA-sensing pathway and gastric acid secretion (Fig. 4C and D ). Additionally, the top 20 hub genes, CDK1, CCNB1, MYC, CCNA2, MAD2L1, AURKA, TOP2A, CDC6, UBE2C, CHEK1, RRM2, BUB1B, TTK, TRIP13, TPX2, BUB1, NCAPG, KIF2C, KIF23, and MCM4 were identified with higher degrees of connectivity. These hub genes were enriched in cell cycle, progesterone-mediated oocyte maturation, oocyte meiosis, p53 signaling pathway, and HTLV-I infection (Fig. 4E and F ).
To validate the hub genes we got from this study, we uploaded the hub genes list into GEPIA database and explored the correlation between hub genes expression and the clinical characteristics of CRC. It was found that HMMR, PAICS, ETFDH, and SCG2 were significant DEGs in 270 CRC samples from GEPIA (Fig. 5A ). And these four genes could represent the important prognostic biomarkers for predicting the survival of CRC patients (Fig. 5B ). Meanwhile, PAICS and SCG2 were related to the stages of CRC progression (Fig. 5C ). The summaries of four hub genes were shown in Table 1 .
Gene | Summary | Microarray datasets [P-value, log2(fold change)] | |||
GSE113513 | GSE21510 | GSE44076 | GSE32323 | ||
HMMR | The protein encoded by this gene is involved in cell motility. It is expressed in breast tissue and together with other proteins, it forms a complex with BRCA1 and BRCA2, thus is potentially associated with higher risk of breast cancer. Alternatively spliced transcript variants encoding different isoforms have been noted for this gene. | 6.82e−04, 1.11 | 3.84e−35, 3.22 | 5.32e−22, 1.20 | 2.63e−04, 1.63 |
PAICS | This gene encodes a bifunctional enzyme containing phosphoribosylaminoimidazole carboxylase activity in its N-terminal region and phosphoribosylaminoimidazole succinocarboxamide synthetase in its C-terminal region. It catalyzes steps 6 and 7 of purine biosynthesis. The gene is closely linked and divergently transcribed with a locus that encodes an enzyme in the same pathway, and transcription of the two genes is coordinately regulated. The human genome contains several pseudogenes of this gene. Multiple transcript variants encoding different isoforms have been found for this gene. | 4.77e−08, 1.31 | 5.71e−47, 2.48 | 2.14e−55, 1.49 | 1.38e−06, 1.73 |
ETFDH | This gene encodes a component of the electron-transfer system in mitochondria and is essential for electron transfer from a number of mitochondrial flavin-containing dehydrogenases to the main respiratory chain. Mutations in this gene are associated with glutaric acidemia. Alternatively spliced transcript variants that encode distinct isoforms have been observed. | 1.78e−08, −1.47 | 1.67e−27, −1.92 | 1.80e−71, −1.85 | 3.14e−08, −1.78 |
SCG2 | The protein encoded by this gene is a member of the chromogranin/secretogranin family of neuroendocrine secretory proteins. Studies in rodents suggest that the full-length protein, secretogranin II, is involved in the packaging or sorting of peptide hormones and neuropeptides into secretory vesicles. The full-length protein is cleaved to produce the active peptide secretoneurin, which exerts chemotaxic effects on specific cell types, and EM66, whose function is unknown. | 2.37e−07, −2.50 | 1.13e−17, −1.86 | 4.83e−25, −1.90 | 2.91e−07, −2.28 |
To search candidate small molecules for developing potential drugs to treat CRC, we uploaded DEGs probe profiles into the CMap database. And the predicted results were download and filtered with enrichment score <0 and P-value <0.05. The results were shown in Table 2 . And Fig. 5D listed the top 20 small molecules with their enrichment scores and P-values. Therefore, these small molecules may be the targets to develop new drugs or therapies of CRC. Among these molecules, Blebbistatin and Sulconazole may be selected for new clinical trials.
CMap name | Count | Enrichment | P-value | CMap name | Count | Enrichment | P-value | |
DL-thiorphan | 2 | −0.976 | 0.001 27 | 8-azaguanine | 4 | −0.644 | 0.038 11 | |
Quinostatin | 2 | −0.878 | 0.029 90 | Acepromazine | 4 | −0.640 | 0.040 08 | |
1,4-chrysenequinone | 2 | −0.866 | 0.035 75 | Doxazosin | 4 | −0.639 | 0.040 48 | |
Menadione | 2 | −0.852 | 0.043 24 | Mefloquine | 5 | −0.637 | 0.016 28 | |
Blebbistatin | 2 | −0.844 | 0.048 59 | Trifluridine | 4 | −0.637 | 0.041 87 | |
trazodone | 3 | −0.813 | 0.013 20 | Clomipramine | 4 | −0.636 | 0.042 51 | |
Thioguanosine | 4 | −0.809 | 0.002 59 | Trioxysalen | 4 | −0.633 | 0.044 24 | |
Sulconazole | 4 | −0.805 | 0.002 82 | Corbadrine | 4 | −0.629 | 0.046 73 | |
0297417-0002B | 3 | −0.792 | 0.018 27 | Etofenamate | 4 | −0.624 | 0.049 27 | |
Etacrynic acid | 3 | −0.776 | 0.022 99 | Phthalylsulfathiazole | 5 | −0.616 | 0.022 71 | |
Latamoxef | 3 | −0.759 | 0.028 76 | Oxetacaine | 5 | −0.610 | 0.025 23 | |
GW-8510 | 4 | −0.757 | 0.007 12 | Amiodarone | 5 | −0.601 | 0.029 00 | |
Triflusal | 3 | −0.755 | 0.030 11 | Meticrane | 5 | −0.599 | 0.029 62 | |
Piperidolate | 3 | −0.743 | 0.034 63 | Vorinostat | 12 | −0.596 | 0.000 18 | |
Daunorubicin | 4 | −0.733 | 0.010 19 | Tranylcypromine | 5 | −0.592 | 0.032 64 | |
Repaglinide | 4 | −0.727 | 0.011 42 | Bromocriptine | 5 | −0.573 | 0.042 46 | |
Camptothecin | 3 | −0.725 | 0.042 85 | Zimeldine | 5 | −0.572 | 0.043 38 | |
Norcyclobenzaprine | 4 | −0.723 | 0.012 02 | Clemizole | 5 | −0.568 | 0.045 48 | |
Ronidazole | 3 | −0.711 | 0.049 19 | Astemizole | 5 | −0.567 | 0.046 14 | |
Ellipticine | 4 | −0.708 | 0.015 14 | Meclozine | 5 | −0.565 | 0.047 08 | |
Gliclazide | 4 | −0.707 | 0.015 28 | Procaine | 5 | −0.565 | 0.047 46 | |
Phenoxybenzamine | 4 | −0.704 | 0.015 76 | Cloperastine | 6 | −0.551 | 0.031 56 | |
0175029-0000 | 6 | −0.701 | 0.001 63 | Famprofazone | 6 | −0.543 | 0.036 21 | |
Tyloxapol | 4 | −0.694 | 0.018 52 | Dipyridamole | 6 | −0.537 | 0.039 63 | |
Skimmianine | 4 | −0.690 | 0.020 01 | Promazine | 6 | −0.523 | 0.047 85 | |
Bisacodyl | 4 | −0.690 | 0.020 03 | Trifluoperazine | 16 | −0.504 | 0.000 20 | |
Resveratrol | 9 | −0.684 | 0.000 04 | Prochlorperazine | 16 | −0.434 | 0.002 99 | |
Medrysone | 6 | −0.680 | 0.002 92 | LY-294002 | 61 | −0.428 | 0.000 00 | |
Bepridil | 4 | −0.674 | 0.025 22 | Trichostatin A | 182 | −0.425 | 0.000 00 | |
Pyrvinium | 6 | −0.672 | 0.003 34 | Fluphenazine | 18 | −0.406 | 0.003 80 | |
Nortriptyline | 4 | −0.669 | 0.027 07 | Thioridazine | 20 | −0.388 | 0.003 27 | |
Apigenin | 4 | −0.667 | 0.027 65 | Alpha-estradiol | 16 | −0.357 | 0.025 35 | |
Methylergometrine | 4 | −0.664 | 0.028 71 | Geldanamycin | 15 | −0.349 | 0.039 26 | |
Bufexamac | 4 | −0.660 | 0.030 10 | Sirolimus | 44 | −0.335 | 0.000 08 | |
Prestwick-1084 | 4 | −0.656 | 0.032 41 | Chlorpromazine | 19 | −0.329 | 0.024 17 | |
Deptropine | 4 | −0.651 | 0.034 95 | Wortmannin | 18 | −0.314 | 0.046 91 | |
Dextromethorphan | 4 | −0.647 | 0.036 86 | Tanespimycin | 62 | −0.311 | 0.000 00 | |
Protriptyline | 4 | −0.646 | 0.037 36 | Fulvestrant | 40 | −0.225 | 0.029 99 |
In summary, we chose GSE113513, GSE21510, GSE44076, and GSE32323 GEO datasets and found 865 significant DEGs between CRC tissues and normal tissues. Subsequently, the biological roles of these DEGs were confirmed with enrichment pathway analysis. Furthermore, the four hub genes, HMMR, PAICS, ETFDH, and SCG2 were identified as important prognostic biomarkers for predicting the survival of CRC patients based on the GEPIA database. Finally, blebbistatin and sulconazole were picked out to develop new drugs through CMap database (Fig. 6 ).
In our study, we chose four GEO datasets and used bioinformatics methods to get 865 DEGs (374 upregulated and 491 downregulated). KEGG pathway analysis showed that the key modules were mainly metabolic pathways, pathways in cancer, cell cycle, purine metabolism, pancreatic secretion, thyroid hormone signaling pathway and Wnt signaling pathway. The PPI network was constructed including 863 nodes and 5817 edges. The four hub genes, HMMR, PAICS, ETFDH, and SCG2 were remarkably related to the prognosis of patients. Furthermore, two small molecules, blebbistatin and sulconazole, also have been identified as potential candidates to develop new drugs.
Recently, findings about DEGs or molecular biomarkers of CRC have been increasingly reported. Based on integrated analysis of GSE32323, GSE74602, and GSE113513 datasets, and TCGA databases, CCL19, CXCL1, CXCL5, CXCL11, CXCL12, GNG4, INSL5, NMU, PYY, and SST were identified as hub genes. And 9 genes including SLC4A4, NFE2L3, GLDN, PCOLCE2, TIMP1, CCL28, SCGB2A1, AXIN2, and MMP1 was related to predicting overall survivals of CRC patients[8]. Moreover, TOP2A, MAD2L1, CCNB1, CHEK1, CDC6, and UBE2C were indicated as hub genes, and TOP2A, MAD2L1, CDC6, and CHEK1 may serve as prognostic biomarkers in CRC[10]. In addition, CEACAM7, SLC4A4, GCG, and CLCA1 genes were associated with unfavorable prognosis in CRC[11]. According to analysis of GEO datasets and survival analysis by GEPIA database, AURKA, CCNB1, CCNF, and EXO1 were significantly associated with longer overall survival. Moreover, CMap predicted that DL-thiorphan, repaglinide, MS-275, and quinostatin have the potential to treat CRC[9]. In this study, we have identified four hub genes as new potential biomarkers to predict the prognosis of CRC patients and two new small molecules. Further studies are needed to develop new drugs to treat CRC.
Several studies have reported that these hub genes play important roles in cancer development. For instance, HMMR expression level was remarkably correlated with the progression and prognosis of breast cancer[13], bladder cancer[14], prostate cancer[15–16], lung cancer[17–19], hepatocellular carcinoma (HCC)[20–22], and gastric cancer[23]. Furthermore, HMMR was confirmed to maintain its oncogenic properties and resistance to chemotherapy through activating TGF-β/Smad-2 signaling pathway[24]. And HMMR was highly expressed in glioblastoma and related to support the self-renewal and tumorigenic potential of glioblastoma stem cells[25]. PAICS was also upregulated in several kinds of cancer tissues and it promotes cancer cells proliferation, migration, and invasion[26–31]. The expression level of ETFDH was found significantly decreased in HCC tissues, and this low expression was related to poor overall survival in patients[32]. However, the role of SCG2 in cancer remains unclear.
In the present study, HMMR, PAICS, ETFDH, and SCG2 were significantly up or down regulated in CRC tissues compared with those in normal samples, and the survival rate of CRC patients was positively correlated with the expression of these genes. Besides, several small molecules with potential therapeutic efficacy were identified through bioinformatics analyses, including blebbistatin and sulconazole. Blebbistatin has been reported to inhibit cell migration and invasiveness of pancreatic adenocarcinoma[31], and decrease spreading and migration of breast cancer cells[33]. Moreover, blebbistatin has shown its antitumorigenic properties in HCC cells[34]. Another small molecule, sulconazole, also inhibited the proliferation and formation of breast cancer stem cells through blocking the NF-κB/IL-8 signaling pathway[35]. Although these two molecules have significant antitumor activity, their specific roles in CRC development need to be further clarified.
This work was supported by grants from the National Natural Science Foundation of China (No. 81672748 and No. 81871936).
CLC number: X592; R114, Document code: A
The authors reported no conflict of interests.
[1] |
Skakkebæk NE, Lindahl-Jacobsen R, Levine H, et al. Environmental factors in declining human fertility[J]. Nat Rev Endocrinol, 2022, 18(3): 139–157. doi: 10.1038/s41574-021-00598-8
|
[2] |
Owumi SE, Otunla MT, Arunsi UO, et al. 3-Indolepropionic acid upturned male reproductive function by reducing oxido-inflammatory responses and apoptosis along the hypothalamic-pituitary-gonadal axis of adult rats exposed to chlorpyrifos[J]. Toxicology, 2021, 463: 152996. doi: 10.1016/j.tox.2021.152996
|
[3] |
Bisconti M, Simon JF, Grassi S, et al. Influence of risk factors for male infertility on sperm protein composition[J]. Int J Mol Sci, 2021, 22(23): 13164. doi: 10.3390/ijms222313164
|
[4] |
Calogero AE, Fiore M, Giacone F, et al. Exposure to multiple metals/metalloids and human semen quality: a cross-sectional study[J]. Ecotoxicol Environ Saf, 2021, 215: 112165. doi: 10.1016/j.ecoenv.2021.112165
|
[5] |
Wrzecińska M, Kowalczyk A, Cwynar P, et al. Disorders of the Reproductive Health of Cattle as a Response to Exposure to Toxic Metals[J]. Biology, 2021, 10(9): 882. doi: 10.3390/biology10090882
|
[6] |
Rahman Z, Singh VP. The relative impact of toxic heavy metals (THMs) (arsenic (As), cadmium (Cd), chromium (Cr)(VI), mercury (Hg), and lead (Pb)) on the total environment: an overview[J]. Environ Monit Assess, 2019, 191(7): 419. doi: 10.1007/s10661-019-7528-7
|
[7] |
Charkiewicz AE, Backstrand JR. Lead toxicity and pollution in Poland[J]. Int J Environ Res Public Health, 2020, 17(12): 4385. doi: 10.3390/ijerph17124385
|
[8] |
Advisory Committee on Childhood Lead Poisoning Prevention. Recommendations for blood lead screening of young children enrolled in Medicaid: targeting a group at high risk[J]. MMWR Recomm Rep, 2000, 49(RR-14): 1–13. https://pubmed.ncbi.nlm.nih.gov/11147551/
|
[9] |
Vorvolakos T, Arseniou S, Samakouri M. There is no safe threshold for lead exposure: α literature review[J]. Psychiatriki, 2016, 27(3): 204–214. doi: 10.22365/jpsych.2016.273.204
|
[10] |
Rocha A, Trujillo KA. Neurotoxicity of low-level lead exposure: history, mechanisms of action, and behavioral effects in humans and preclinical models[J]. Neurotoxicology, 2019, 73: 58–80. doi: 10.1016/j.neuro.2019.02.021
|
[11] |
Kumar S. Occupational and environmental exposure to lead and reproductive health impairment: an overview[J]. Indian J Occup Environ Med, 2018, 22(3): 128–137. doi: 10.4103/ijoem.IJOEM_126_18
|
[12] |
Abdel-Emam RA, Ahmed EA. Ameliorative effect of L-carnitine on chronic lead-induced reproductive toxicity in male rats[J]. Vet Med Sci, 2021, 7(4): 1426–1435. doi: 10.1002/vms3.473
|
[13] |
Seok YJ, Her JY, Kim YG, et al. Furan in thermally processed foods - a review[J]. Toxicol Res, 2015, 31(3): 241–253. doi: 10.5487/TR.2015.31.3.241
|
[14] |
European Food Safety Authority. Update on furan levels in food from monitoring years 2004–2010 and exposure assessment[J]. EFSA J, 2011, 9(9): 2347. doi: 10.2903/j.efsa.2011.2347
|
[15] |
EFSA Panel on Contaminants in the Food Chain, Knutsen HK, Alexander J, et al. Risks for public health related to the presence of furan and methylfurans in food[J]. EFSA J, 2017, 15(10): e05005. doi: 10.2903/j.efsa.2017.5005
|
[16] |
Faul F, Erdfelder E, Lang AG, et al. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences[J]. Behav Res Methods, 2007, 39(2): 175–191. doi: 10.3758/BF03193146
|
[17] |
Owumi SE, Adedara IA, Farombi EO, et al. Protocatechuic acid modulates reproductive dysfunction linked to furan exposure in rats[J]. Toxicology, 2020, 442: 152556. doi: 10.1016/j.tox.2020.152556
|
[18] |
Kataba A, Botha TL, Nakayama SMM, et al. Environmentally relevant lead (Pb) water concentration induce toxicity in zebrafish (Danio rerio) larvae[J]. Comp Biochem Physiol C Toxicol Pharmacol, 2022, 252: 109215. doi: 10.1016/j.cbpc.2021.109215
|
[19] |
Cressey D. Best way to kill lab animals sought[J]. Nature, 2013, 500(7461): 130–131. doi: 10.1038/500130a
|
[20] |
Zemjanis R. Collection and evaluation of semen[M]//Zemjanis R. Diagnostic and Therapeutic Techniques in Animal Reproduction. 2nd ed. Baltimore: Williams & Wilkins, 1970.
|
[21] |
World Health Organization. WHO laboratory manual for the examination and processing of human semen[EB/OL]. [2019-09-01]. https://www.who.int/publications/i/item/9789240030787.
|
[22] |
Wells ME, Awa OA. New technique for assessing Acrosomal characteristics of spermatozoa[J]. J Dairy Sci, 1970, 53(2): 227–232. doi: 10.3168/jds.S0022-0302(70)86184-7
|
[23] |
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding[J]. Anal Biochem, 1976, 72(1–2): 248–254. doi: 10.1006/abio.1976.9999
|
[24] |
Malymy M, Horecker BL. Alkaline phosphatase[M]//Willis AW. Methods in Enzymology. New York: Academic Press, 1966: 639–642.
|
[25] |
Vanha‐Perttula T, Nikkanen V. Acid phosphatases of the rat testis in experimental conditions[J]. Eur J Endocrinol, 1973, 72(2): 376–390. doi: 10.1530/acta.0.0720376
|
[26] |
Wolf BHM, Weening RS, Schutgens RBH, et al. Detection of glucose-6-phosphate dehydrogenase deficiency in erythrocytes: a spectrophotometric assay and a fluorescent spot test compared with a cytochemical method[J]. Clin Chim Acta, 1987, 168(2): 129–136. doi: 10.1016/0009-8981(87)90281-6
|
[27] |
Vassault A. Lactate dehydrogenase. UV-method with pyruvate and NADH[M]//Bergmeyer HU. Methods of Enzymatic Analysis. 3rd ed. New York: Plenum, 1993: 118–125.
|
[28] |
Misra HP, Fridovich I. The role of superoxide anion in the autoxidation of epinephrine and a simple assay for superoxide dismutase[J]. J Biol Chem, 1972, 247(10): 3170–3175. doi: 10.1016/S0021-9258(19)45228-9
|
[29] |
Ellman GL. Tissue sulfhydryl groups[J]. Arch Biochem Biophys, 1959, 82(1): 70–77. doi: 10.1016/0003-9861(59)90090-6
|
[30] |
Clairborne A. Catalase activity[M]. Boca Raton: CRC Press, 1995: 283–284.
|
[31] |
Habig WH, Pabst MJ, Jakoby WB. Glutathione S-transferases: the first enzymatic step in mercapturic acid formation[J]. J Biol Chem, 1974, 249(22): 7130–7139. doi: 10.1016/S0021-9258(19)42083-8
|
[32] |
Rotruck JT, Pope AL, Ganther HE, et al. Selenium: biochemical role as a component of glutathione peroxidase[J]. Science, 1973, 179(4073): 588–590. doi: 10.1126/science.179.4073.588
|
[33] |
Jollow DJ, Mitchell JR, Zampaglione N, et al. Bromobenzene-induced liver necrosis. Protective role of glutathione and evidence for 3, 4-bromobenzene oxide as the hepatotoxic metabolite[J]. Pharmacology, 1974, 11(3): 151–169. doi: 10.1159/000136485
|
[34] |
Bergmeyer HI, Gawehn K. Methods of enzymatic analysis[M]. 2nd ed. Amsterdam: Elsevier, 1974: 521–522.
|
[35] |
Ohkawa H, Ohishi N, Yagi K. Assay for lipid peroxides in animal tissues by Thiobarbituric acid reaction[J]. Anal Biochem, 1979, 95(2): 351–358. doi: 10.1016/0003-2697(79)90738-3
|
[36] |
Owumi SE, Dim UJ. Manganese suppresses oxidative stress, inflammation and caspase-3 activation in rats exposed to chlorpyrifos[J]. Toxicol Rep, 2019, 6: 202–209. doi: 10.1016/j.toxrep.2019.02.007
|
[37] |
Green LC, Wagner DA, Glogowski J, et al. Analysis of nitrate, nitrite, and [15N]nitrate in biological fluids[J]. Anal Biochem, 1982, 126(1): 131–138. doi: 10.1016/0003-2697(82)90118-X
|
[38] |
Granell S, Gironella M, Bulbena O, et al. Heparin mobilizes xanthine oxidase and induces lung inflammation in acute pancreatitis[J]. Crit Care Med, 2003, 31(2): 525–530. doi: 10.1097/01.CCM.0000049948.64660.06
|
[39] |
Bancroft JD, Gamble M. Theory and practice of histological techniques[M]. 6th ed. Amsterdam: Elsevier, 2007.
|
[40] |
Owumi SE, Irozuru CE, Arunsi UO, et al. Caffeic acid mitigates aflatoxin B1-mediated toxicity in the male rat reproductive system by modulating inflammatory and apoptotic responses, testicular function, and the redox-regulatory systems[J]. J Food Biochem, 2022, 46(5): e14090. doi: 10.1111/jfbc.14090
|
[41] |
Teo CH, Soga T, Parhar I. Lithium chloride enhances serotonin induced calcium activity in EGFP-GnIH neurons[J]. Sci Rep, 2020, 10(1): 13876. doi: 10.1038/s41598-020-70710-x
|
[42] |
Parhar IS, Ogawa S, Ubuka T. Reproductive neuroendocrine pathways of social behavior[J]. Front Endocrinol, 2016, 7: 28. doi: 10.3389/fendo.2016.00028
|
[43] |
Peruquetti RL, Taboga SR, Azeredo-Oliveira MT. Expression of acid phosphatase in the seminiferous epithelium of vertebrates[J]. Genet Mol Res, 2010, 9(2): 620–628. doi: 10.4238/vol9-2gmr730
|
[44] |
Gurel H, Baspinar N, Akalin PP, et al. Erythrocyte and spermatozoa glucose-6-phosphate dehydrogenase activity in merino rams: an experimental study[J]. Int J Reprod Biomed, 2018, 16(6): 373–378. doi: 10.29252/ijrm.16.6.373
|
[45] |
Alibawi FNAA, Al-Morshidy SY, Alhuweizi AG. The alkaline phosphatase levels in the seminal plasma and sperms of sub-fertile patients and Normospermic men[C]//Proceedings of 2012 International Conference on Applied Life Sciences. Konya, Turkey. ISALS, 2012: 217–222.
|
[46] |
Owumi SE, Anaikor RA, Arunsi UO, et al. Chlorogenic acid co-administration abates tamoxifen-mediated reproductive toxicities in male rats: An experimental approach[J]. J Food Biochem, 2021, 45(2): e13615. doi: 10.1111/jfbc.13615
|
[47] |
Bidanchi RM, Lalrindika L, Khushboo M, et al. Antioxidative, anti-inflammatory and anti-apoptotic action of Ellagic acid against lead acetate induced testicular and hepato-renal oxidative damages and pathophysiological changes in male Long Evans rats[J]. Environ Pollut, 2022, 302: 119048. doi: 10.1016/j.envpol.2022.119048
|
[48] |
Yuan Y, Wang Z, Nan B, et al. Salidroside alleviates liver inflammation in furan-induced mice by regulating oxidative stress and endoplasmic reticulum stress[J]. Toxicology, 2021, 461: 152905. doi: 10.1016/j.tox.2021.152905
|
[49] |
Owumi SE, Adedara IA, Otunla MT, et al. Influence of furan and lead co-exposure at environmentally relevant concentrations on neurobehavioral performance, redox-regulatory system and apoptotic responses in rats[J]. Environ Toxicol Pharmacol, 2022, 97: 104011. doi: 10.1016/j.etap.2022.104011
|
[1] | Dongmei Zhu, Lingli Luo, Hanjie Zeng, Zheng Zhang, Min Huang, Suming Zhou. Knockdown of 11β-hydroxysteroid dehydrogenase type 1 alleviates LPS-induced myocardial dysfunction through the AMPK/SIRT1/PGC-1α pathway[J]. The Journal of Biomedical Research, 2023, 37(4): 290-301. DOI: 10.7555/JBR.36.20220212 |
[2] | Wang Lei, Zhang Kun, Ding Xin, Wang Yan, Bai Hui, Yang Qing, Ben Jingjing, Zhang Hanwen, Li Xiaoyu, Chen Qi, Zhu Xudong. Fucoidan antagonizes diet-induced obesity and inflammation in mice[J]. The Journal of Biomedical Research, 2021, 35(3): 197-205. DOI: 10.7555/JBR.34.20200153 |
[3] | Venakatesh Archana G., Mathew Johann J., Coleman Scott, Yang Longqiu, Liu Geoffrey L., Li Marilyn M., Liu Henry. Effects of milrinone on inflammatory response-related gene expressions in cultured rat cardiomyocytes[J]. The Journal of Biomedical Research, 2019, 33(4): 258-263. DOI: 10.7555/JBR.32.20170085 |
[4] | Alireza Shirpoor, Farzaneh Hosseini Gharalari, Yousef Rasmi, Elaheh Heshmati. Ginger extract attenuates ethanol-induced pulmonary histological changes and oxidative stress in rats[J]. The Journal of Biomedical Research, 2017, 31(6): 521-527. DOI: 10.7555/JBR.31.20160151 |
[5] | Qin Wei, Yeping Bian, Fuchao Yu, Qiang Zhang, Guanghao Zhang, Yang Li, Songsong Song, Xiaomei Ren, Jiayi Tong. Chronic intermittent hypoxia induces cardiac inflammation and dysfunction in a rat obstructive sleep apnea model[J]. The Journal of Biomedical Research, 2016, 30(6): 490-495. DOI: 10.7555/JBR.30.20160110 |
[6] | Iyaswamy Ashok, Rathinasamy Sheeladevi, Dapkupar Wankhar. Acute effect of aspartame-induced oxidative stress in Wistar albino rat brain[J]. The Journal of Biomedical Research, 2015, 29(5): 390-396. DOI: 10.7555/JBR.28.20120118 |
[7] | Guoliang Meng, Jing Wang, Yujiao Xiao, Wenli Bai, Liping Xie, Liyang Shan, Philip K Moore, Yong Ji. GYY4137 protects against myocardial ischemia and reperfusion injury by attenuating oxidative stress and apoptosis in rats[J]. The Journal of Biomedical Research, 2015, 29(3): 203-213. DOI: 10.7555/JBR.28.20140037 |
[8] | Venkatashivam Shiva Kumar, Anuchandra Ramchandra Rajmane, Mohammad Adil, Amit Dattatraya Kandhare, Pinaki Ghosh, Subhash Laxman Bodhankar. Naringin ameliorates acetic acid induced colitis through modulation of endogenous oxido-nitrosative balance and DNA damage in rats[J]. The Journal of Biomedical Research, 2014, 28(2): 132-145. DOI: 10.7555/JBR.27.20120082 |
[9] | Ze Guo, Bingrong Zhou, Wei Li, Xuejun Sun, Dan Luo. Hydrogen-rich saline protects against ultraviolet B radiation injury in rats[J]. The Journal of Biomedical Research, 2012, 26(5): 365-371. DOI: 10.7555/JBR.26.20110037 |
[10] | Bo Cui, Stewart P. Johnson, Nancy Bullock, Francis Ali-Osman, Darell D. Bigner, Henry S. Friedman. Decoupling of DNA damage response signaling from DNA damages underlies temozolomide resistance in glioblastoma cells[J]. The Journal of Biomedical Research, 2010, 24(6): 424-435. DOI: 10.1016/S1674-8301(10)60057-7 |
1. | Liu J, Yang H, Li P, et al. Microarray analysis points to LMNB1 and JUN as potential target genes for predicting metastasis promotion by etoposide in colorectal cancer. Sci Rep, 2024, 14(1): 23661. DOI:10.1038/s41598-024-72674-8 |
2. | Horaira MA, Islam MA, Kibria MK, et al. Bioinformatics screening of colorectal-cancer causing molecular signatures through gene expression profiles to discover therapeutic targets and candidate agents. BMC Med Genomics, 2023, 16(1): 64. DOI:10.1186/s12920-023-01488-w |
3. | Tang YP, Yin YX, Xie MZ, et al. Systematic Analysis of the Clinical Significance of Hyaluronan-Mediated Motility Receptor in Colorectal Cancer. Front Mol Biosci, 2021, 8: 733271. DOI:10.3389/fmolb.2021.733271 |
Gene | Summary | Microarray datasets [P-value, log2(fold change)] | |||
GSE113513 | GSE21510 | GSE44076 | GSE32323 | ||
HMMR | The protein encoded by this gene is involved in cell motility. It is expressed in breast tissue and together with other proteins, it forms a complex with BRCA1 and BRCA2, thus is potentially associated with higher risk of breast cancer. Alternatively spliced transcript variants encoding different isoforms have been noted for this gene. | 6.82e−04, 1.11 | 3.84e−35, 3.22 | 5.32e−22, 1.20 | 2.63e−04, 1.63 |
PAICS | This gene encodes a bifunctional enzyme containing phosphoribosylaminoimidazole carboxylase activity in its N-terminal region and phosphoribosylaminoimidazole succinocarboxamide synthetase in its C-terminal region. It catalyzes steps 6 and 7 of purine biosynthesis. The gene is closely linked and divergently transcribed with a locus that encodes an enzyme in the same pathway, and transcription of the two genes is coordinately regulated. The human genome contains several pseudogenes of this gene. Multiple transcript variants encoding different isoforms have been found for this gene. | 4.77e−08, 1.31 | 5.71e−47, 2.48 | 2.14e−55, 1.49 | 1.38e−06, 1.73 |
ETFDH | This gene encodes a component of the electron-transfer system in mitochondria and is essential for electron transfer from a number of mitochondrial flavin-containing dehydrogenases to the main respiratory chain. Mutations in this gene are associated with glutaric acidemia. Alternatively spliced transcript variants that encode distinct isoforms have been observed. | 1.78e−08, −1.47 | 1.67e−27, −1.92 | 1.80e−71, −1.85 | 3.14e−08, −1.78 |
SCG2 | The protein encoded by this gene is a member of the chromogranin/secretogranin family of neuroendocrine secretory proteins. Studies in rodents suggest that the full-length protein, secretogranin II, is involved in the packaging or sorting of peptide hormones and neuropeptides into secretory vesicles. The full-length protein is cleaved to produce the active peptide secretoneurin, which exerts chemotaxic effects on specific cell types, and EM66, whose function is unknown. | 2.37e−07, −2.50 | 1.13e−17, −1.86 | 4.83e−25, −1.90 | 2.91e−07, −2.28 |
CMap name | Count | Enrichment | P-value | CMap name | Count | Enrichment | P-value | |
DL-thiorphan | 2 | −0.976 | 0.001 27 | 8-azaguanine | 4 | −0.644 | 0.038 11 | |
Quinostatin | 2 | −0.878 | 0.029 90 | Acepromazine | 4 | −0.640 | 0.040 08 | |
1,4-chrysenequinone | 2 | −0.866 | 0.035 75 | Doxazosin | 4 | −0.639 | 0.040 48 | |
Menadione | 2 | −0.852 | 0.043 24 | Mefloquine | 5 | −0.637 | 0.016 28 | |
Blebbistatin | 2 | −0.844 | 0.048 59 | Trifluridine | 4 | −0.637 | 0.041 87 | |
trazodone | 3 | −0.813 | 0.013 20 | Clomipramine | 4 | −0.636 | 0.042 51 | |
Thioguanosine | 4 | −0.809 | 0.002 59 | Trioxysalen | 4 | −0.633 | 0.044 24 | |
Sulconazole | 4 | −0.805 | 0.002 82 | Corbadrine | 4 | −0.629 | 0.046 73 | |
0297417-0002B | 3 | −0.792 | 0.018 27 | Etofenamate | 4 | −0.624 | 0.049 27 | |
Etacrynic acid | 3 | −0.776 | 0.022 99 | Phthalylsulfathiazole | 5 | −0.616 | 0.022 71 | |
Latamoxef | 3 | −0.759 | 0.028 76 | Oxetacaine | 5 | −0.610 | 0.025 23 | |
GW-8510 | 4 | −0.757 | 0.007 12 | Amiodarone | 5 | −0.601 | 0.029 00 | |
Triflusal | 3 | −0.755 | 0.030 11 | Meticrane | 5 | −0.599 | 0.029 62 | |
Piperidolate | 3 | −0.743 | 0.034 63 | Vorinostat | 12 | −0.596 | 0.000 18 | |
Daunorubicin | 4 | −0.733 | 0.010 19 | Tranylcypromine | 5 | −0.592 | 0.032 64 | |
Repaglinide | 4 | −0.727 | 0.011 42 | Bromocriptine | 5 | −0.573 | 0.042 46 | |
Camptothecin | 3 | −0.725 | 0.042 85 | Zimeldine | 5 | −0.572 | 0.043 38 | |
Norcyclobenzaprine | 4 | −0.723 | 0.012 02 | Clemizole | 5 | −0.568 | 0.045 48 | |
Ronidazole | 3 | −0.711 | 0.049 19 | Astemizole | 5 | −0.567 | 0.046 14 | |
Ellipticine | 4 | −0.708 | 0.015 14 | Meclozine | 5 | −0.565 | 0.047 08 | |
Gliclazide | 4 | −0.707 | 0.015 28 | Procaine | 5 | −0.565 | 0.047 46 | |
Phenoxybenzamine | 4 | −0.704 | 0.015 76 | Cloperastine | 6 | −0.551 | 0.031 56 | |
0175029-0000 | 6 | −0.701 | 0.001 63 | Famprofazone | 6 | −0.543 | 0.036 21 | |
Tyloxapol | 4 | −0.694 | 0.018 52 | Dipyridamole | 6 | −0.537 | 0.039 63 | |
Skimmianine | 4 | −0.690 | 0.020 01 | Promazine | 6 | −0.523 | 0.047 85 | |
Bisacodyl | 4 | −0.690 | 0.020 03 | Trifluoperazine | 16 | −0.504 | 0.000 20 | |
Resveratrol | 9 | −0.684 | 0.000 04 | Prochlorperazine | 16 | −0.434 | 0.002 99 | |
Medrysone | 6 | −0.680 | 0.002 92 | LY-294002 | 61 | −0.428 | 0.000 00 | |
Bepridil | 4 | −0.674 | 0.025 22 | Trichostatin A | 182 | −0.425 | 0.000 00 | |
Pyrvinium | 6 | −0.672 | 0.003 34 | Fluphenazine | 18 | −0.406 | 0.003 80 | |
Nortriptyline | 4 | −0.669 | 0.027 07 | Thioridazine | 20 | −0.388 | 0.003 27 | |
Apigenin | 4 | −0.667 | 0.027 65 | Alpha-estradiol | 16 | −0.357 | 0.025 35 | |
Methylergometrine | 4 | −0.664 | 0.028 71 | Geldanamycin | 15 | −0.349 | 0.039 26 | |
Bufexamac | 4 | −0.660 | 0.030 10 | Sirolimus | 44 | −0.335 | 0.000 08 | |
Prestwick-1084 | 4 | −0.656 | 0.032 41 | Chlorpromazine | 19 | −0.329 | 0.024 17 | |
Deptropine | 4 | −0.651 | 0.034 95 | Wortmannin | 18 | −0.314 | 0.046 91 | |
Dextromethorphan | 4 | −0.647 | 0.036 86 | Tanespimycin | 62 | −0.311 | 0.000 00 | |
Protriptyline | 4 | −0.646 | 0.037 36 | Fulvestrant | 40 | −0.225 | 0.029 99 |