Exploring the potential biological function of GRK2 in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is a malignant tumor of the gastrointestinal tract with high morbidity and mortality. There is growing evidence that GRK2 plays a key role in the development and progression of several human cancers. However, the role and potential mechanisms of GRK2 in colon cancer (COAD) are unclear. METHODS: The expression data of GRK2 was downloaded from The Cancer Genome Atlas database (TCGA). Variation in GRK2 was explored based on the cBioPortal database. The TIMER and TISCH2 databases were used to analyse the relationship between GRK2 expression and tumor immune microenvironment (TME). A log-rank test was used to compare the prognosis of high and low expression of GRK2 groups. Detecting the effect of GRK2 on cell cycle and apoptosis induced by 5-Fluorouracil (5-FU) through the flow cytometry and detection of apoptosis-related molecules by Western blot. RESULTS: We demonstrated that GRK2 has a potential oncogenic role. GRK2 expression was upregulated in COAD, which predicted poorer overall survival in COAD patients. The cellular assays showed that GRK2 plays a role in the growth and proliferation of colon cancer cells, also the expression of GRK2 have relationship with the sensitivity of 5-FU and cell cycle progression. CONCLUSIONS: Our results suggest that high GRK2 expression is closely associated with the development of tumor and affects the 5-FU sensitivity.

Fanconi Anemia Complementary Group A (FANCA) Facilitates the Occurrence and Progression of Liver Hepatocellular Carcinoma.

BACKGROUND: Liver hepatocellular carcinoma (LIHC) is a serious liver disease worldwide, and its pathogenesis is complicated. AIMS: This study investigated the potential role of FANCA in the advancement and prognosis of LIHC. METHODS: Public databases, quantitative reverse transcription polymerase chain reaction (qRT-PCR), western blot (WB) and immunohistochemistry (IHC) were employed to measure FANCA expression between tumor and normal samples. The relationship between FANCA expression and prognosis of LIHC patients were examined. Functional enrichment of FANCA-related genes was performed. Furthermore, univariate and multivariate analyses were conducted to determine the independent prognosis value of FANCA in LIHC. Finally, influence of FANCA knockout on the proliferation, migration, and invasion of HepG2 cell was validated with cloning formation, CCK8, and Transwell assays. RESULTS: Expression analysis presented that FANCA had high expression level in LIHC tissues and cells. Receiver operating characteristic (ROC) curve analysis showed that FANCA was of great diagnosis value in LIHC. Clinicopathological analysis revealed that FANCA was significantly greater expressed in the advanced stage than in the early stage of LIHC. Univariate, multivariate, and Kaplan-Meier survival analysis confirmed that high expression of FANCA was strongly associated with poor survival of LIHC patients. In addition, high level of FANCA in LIHC showed a negative association with immunoinfiltrated B cells, T cells, and stromal scores. Moreover, Knockout of FANCA significantly inhibited HepG2 cell proliferative activity, migration, and invasion ability. CONCLUSIONS: Our data revealed that high level of FANCA was closely associated with LIHC malignant progression, suggesting its potential utility as a diagnostic, predictive indicator, and therapeutic target.

Ferroptosis-related genes LPCAT3 and PGD are potential diagnostic biomarkers for osteoarthritis.

BACKGROUND: Osteoarthritis (OA) is the most common chronic joint disease and how ferroptosis contributes to OA has garnered much attention recently. Bioinformatics promoted the discovery of ferroptosis-related biomarkers for OA. But since OA is a whole-joint disease, sensitive biomarkers for OA are still limited. We herein focused on subchondral bone, a joint component often-ignored by existing bioinformatic reports, to identify ferroptosis-related diagnostic biomarkers for OA. METHOD: Microarray datasets GSE51588 and GSE55457 were downloaded from Gene Expression Omnibus database. Ferroptosis-related differential expression genes (Ferr-DEGs) between OA and normal samples were identified and their functional enrichment was analyzed. Common genes for OA diagnosis were selected from Ferr-DEGs using the combination of SVM-RFE, LASSO regression, and RandomForest machine learning algorithms. Common genes' diagnostic value was verified by receiver operating characteristic (ROC) curve and their association with immune infiltration was analyzed by CIBERSORT. Finally, candidate gene's expression was verified in chondrocytes from OA patients and in an in vitro IL-1ß-induced OA model, by RT-PCR. RESULTS: Two ferroptosis-related genes, LPCAT3 and PGD, were identified as OA diagnostic biomarkers and confirmed by ROC diagnostic test. The association of LPCAT3 and PGD with the infiltration of several types of immune cells was identified. The decreased expression of LPCAT3 and PGD was both confirmed in OA chondrocytes and IL-1ß-induced OA condition. CONCLUSIONS: We identified ferroptosis-related genes LPCAT3 and PGD as potential diagnostic biomarkers for OA, which may offer insight into the role of ferroptosis in OA and provides useful information for the diagnosis and treatment of OA.

A Four Amino Acid Metabolism-Associated Genes (AMGs) Signature for Predicting Overall Survival Outcomes and Immunotherapeutic Efficacy in Hepatocellular Carcinoma.

Metabolites are important indicators of cancer and mutations in genes involved in amino acid metabolism may influence tumorigenesis. Immunotherapy is an effective cancer treatment option; however, its relationship with amino acid metabolism has not been reported. In this study, RNA-seq data for 371 liver cancer patients were acquired from TCGA and used as the training set. Data for 231 liver cancer patients were obtained from ICGC and used as the validation set to establish a gene signature for predicting liver cancer overall survival outcomes and immunotherapeutic responses. Four reliable groups based on 132 amino acid metabolism-related DEGs were obtained by consistent clustering of 371 HCC patients and a four-gene signature for prediction of liver cancer survival outcomes was developed. Our data show that in different clinical groups, the overall survival outcomes in the high-risk group were markedly low relative to the low-risk group. Univariate and multivariate analyses revealed that the characteristics of the 4-gene signature were independent prognostic factors for liver cancer. The ROC curve revealed that the risk characteristic is an efficient predictor for 1-, 2-, and 3-year HCC survival outcomes. The GSVA and KEGG pathway analyses revealed that high-risk score tumors were associated with all aspects of the degree of malignancy in liver cancer. There were more mutant genes and greater immune infiltrations in the high-risk groups. Assessment of the three immunotherapeutic cohorts established that low-risk score patients significantly benefited from immunotherapy. Then, we established a prognostic nomogram based on the TCGA cohort. In conclusion, the 4-gene signature is a reliable diagnostic marker and predictor for immunotherapeutic efficacy.

MeGATAs, functional generalists in interactions between cassava growth and development, and abiotic stresses.

The proteins with DNA-binding preference to the consensus DNA sequence (A/T) GATA (A/G) belong to a GATA transcription factor family, with a wide array of biological processes in plants. Cassava (Manihot esculenta) is an important food crop with high production of starch in storage roots. Little was however known about cassava GATA domain-containing genes (MeGATAs). Thirty-six MeGATAs, MeGATA1 to MeGATA36, were found in this study. Some MeGATAs showed a collinear relationship with orthologous genes of Arabidopsis, poplar and potato, rice, maize and sorghum. Eight MeGATA-encoded proteins (MeGATAs) analysed were all localized in the nucleus. Some MeGATAs had potentials of binding ligands and/or enzyme activity. One pair of tandem-duplicated MeGATA17-MeGATA18 and 30 pairs of whole genome-duplicated MeGATAs were found. Fourteen MeGATAs showed low or no expression in the tissues. Nine analysed MeGATAs showed expression responses to abiotic stresses and exogenous phytohormones. Three groups of MeGATA protein interactions were found. Fifty-three miRNAs which can target 18 MeGATAs were identified. Eight MeGATAs were found to target other 292 cassava genes, which were directed to radial pattern formation and phyllome development by gene ontology enrichment, and autophagy by Kyoto Encyclopaedia of Genes and Genomes enrichment. These data suggest that MeGATAs are functional generalists in interactions between cassava growth and development, abiotic stresses and starch metabolism.

A cancer stem cell associated gene signature for predicting overall survival of hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is the most prevalent type of primary liver cancer characterized by high mortality and morbidity rate. The lack of effective treatments and the high frequency of recurrence lead to poor prognosis of patients with HCC. Therefore, it is important to develop robust prediction tools for predicting the prognosis of HCC. Recent studies have shown that cancer stem cells (CSC) participate in HCC progression. The aim of this study was to explore the prognostic value of CSC-related genes and establish a prediction model based on data from The Cancer Genome Atlas (TCGA) database. In this study, 475 CSC-related genes were obtained from the Molecular Signature Database and 160 differentially expressed CSC-related genes in HCC patients were identified using the limma R package in the TCGA database. A total of 79 CSC-related genes were found to be associated with overall survival (OS). Using the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regressions, a 3-gene signature (RAB10, TCOF1, and PSMD14) was constructed. Receiver operating characteristic (ROC) curves and Kaplan-Meier survival curves were constructed to test the prediction performance of the signature. Performance of the signature was validated using the International Cancer Genome Consortium (ICGC) dataset. In addition, immune feature and functional enrichment analyses were carried out to explore the underlying mechanisms. Moreover, a co-expression network was constructed using the weighted gene correlation network analysis (WGCNA) method to select genes significantly associated with risk scores in HCC in the TCGA dataset. The SGO2 gene was found to be significantly associated with risk scores of HCC. In vitro experiments revealed that it can promote HCC cell proliferation. Therefore, SGO2 may be a potential therapeutic target for HCC treatment. The constructed nomogram can help clinicians make decisions about HCC treatment.

Genetic association study of intron variants in the forkhead box protein P3 gene in Chinese patients diagnosed with cervical cancer.

The aim of this study was to investigate the effects of forkhead box protein P3 (FOXP3) intron single nucleotide variants (SNVs) in high-risk human papilloma virus (HR-HPV) infection and cervical cancer (CC) malignant lesions. We performed FOXP3 genotyping in 350 patients with CC and 350 healthy controls using the ImLDR multiple single nucleotide polymorphism genotyping technology. The heterozygous mutation TC in rs2294021 decreased the risk of HR-HPV infection and CC malignant lesions (TC vs. TT: OR = 0.71, 95% CI = 0.51-0.99); the dominant model TC+CC and allele C in rs2294021 decreased the risk of CC malignant lesions (TC+CC vs. TT: OR = 0.69, 95% CI = 0.50-0.95; C vs. T: OR = 0.78, 95% CI = 0.63-0.97). The heterozygous mutation GA, dominant model GA+AA and allele A in rs3761549 also decreased the risk of HR-HPV infection and CC malignant lesions (GA vs. GG: OR = 0.70, 95% CI = 0.51-0.96; GA+AA vs. GG: OR = 0.69, 95% CI = 0.51-0.94; A vs. G: OR = 0.75, 95% CI = 0.58-0.96). Patients with CC and HR-HPV infection carrying rs2294021 TC and rs3761549 GA had lower expression of FOXP3 protein. Haplotype analysis revealed that T-C-A decreased the risk of HR-HPV infection. Furthermore, we found a significant association between immune cells infiltration and prognosis in patients with CC. Our findings demonstrated that rs2294021 and rs3761549 variants may protect against HR-HPV and CC malignant lesions by downregulating FOXP3 and that FOXP3 was associated with immune cells infiltration, which affected the prognosis of CC.

NLK is required for Ras/ERK/SRF/ELK signaling to tune skeletal muscle development by phosphorylating SRF and antagonizing the SRF/MKL pathway.

Serum response factor (SRF) regulates differentiation and proliferation by binding to RhoA-actin-activated MKL or Ras-MAPK-activated ELK transcriptional coactivators, but the molecular mechanisms responsible for SRF regulation remain unclear. Here, we show that Nemo-like kinase (NLK) is required for the promotion of SRF/ELK signaling in human and mouse cells. NLK was found to interact with and phosphorylate SRF at serine residues 101/103, which in turn enhanced the association between SRF and ELK. The enhanced affinity of SRF/ELK antagonized the SRF/MKL pathway and inhibited mouse myoblast differentiation in vitro. In a skeletal muscle-specific Nlk conditional knockout mouse model, forming muscle myofibers underwent hypertrophic growth, resulting in an increased muscle and body mass phenotype. We propose that both phosphorylation of SRF by NLK and phosphorylation of ELKs by MAPK are required for RAS/ELK signaling, confirming the importance of this ancient pathway and identifying an important role for NLK in modulating muscle development in vivo.

Analyses of open-access multi-omics data sets reveal genetic and expression characteristics of maize ZmCCT family genes.

Flowering in maize (Zea mays) is influenced by photoperiod. The CO, CO-like/COL and TOC1 (CCT) domain protein-encoding genes in maize, ZmCCTs, are particularly important for photoperiod sensitivity. However, little is known about CCT protein-encoding gene number across plant species or among maize inbred lines. Therefore, we analysed CCT protein-encoding gene number across plant species, and characterized ZmCCTs in different inbred lines, including structural variations (SVs), copy number variations (CNVs), expression under stresses, dark-dark (DD) and dark-light (DL) cycles, interaction network and associations with maize quantitative trait loci (QTLs) by referring to the latest v4 genome data of B73. Gene number varied greatly across plant species, more in polyploids than in diploids. The numbers of ZmCCTs identified were 58 in B73, 59 in W22, 48 in Mo17, and 57 in Huangzao4 for temperate maize inbred lines, and 68 in tropical maize inbred line SK. Some ZmCCTs underwent duplications and presented chromosome collinearity. Structural variations and CNVs were found but they had no germplasm specificity. Forty-two ZmCCTs responded to stresses. Expression of 37 ZmCCTs in embryonic leaves during seed germination of maize under DD and DL cycles was roughly divided into five patterns of uphill pattern, downhill-pattern, zigzag-pattern, └-pattern and ⅃-pattern, indicating some of them have a potential to perceive dark and/or dark-light transition. Thirty-three ZmCCTs were co-expressed with 218 other maize genes; and 24 ZmCCTs were associated with known QTLs. The data presented in this study will help inform further functions of ZmCCTs.

Dark response genes: a group of endogenous pendulum/timing players in maize?

MAIN CONCLUSION: Maize has a set of dark response genes, expression of which is influenced by multiple factor and varies with maize inbred lines but without germplasm specificity. The response to photoperiod is a common biological issue across the species kingdoms. Dark is as important as light in photoperiod. However, further in-depth understanding of responses of maize (Zea mays) to light and dark transition under photoperiod is hindered due to the lack of understanding of dark response genes. With multiple public "-omic" datasets of temperate and tropical/subtropical maize, 16 maize dark response genes, ZmDRGs, were found and had rhythmic expression under dark and light-dark cycle. ZmDRGs 6-8 were tandemly duplicated. ZmDRGs 2, 13, and 14 had a chromosomal collinearity with other maize genes. ZmDRGs 1-11 and 13-16 had copy-number variations. ZmDRGs 2, 9, and 16 showed 5'-end sequence deletion mutations. Some ZmDRGs had chromatin interactions and underwent DNA methylation and/or m6A mRNA methylation. Chromosomal histones associated with 15 ZmDRGs were methylated and acetylated. ZmDRGs 1, 2, 4, 9, and 13 involved photoperiodic phenotypes. ZmDRG16 was within flowering-related QTLs. ZmDRGs 1, 3, and 6-11 were present in cis-acting expression QTLs (eQTLs). ZmDRGs 1, 4, 6-9, 11, 12, and 14-16 showed co-expression with other maize genes. Some of ZmDRG-encoded ZmDRGs showed obvious differences in abundance and phosphorylation. CONCLUSION: Sixteen ZmDRGs 1-16 are associated with the dark response of maize. In the process of post-domestication and/or breeding, the ZmDRGs undergo the changes without germplasm specificity, including epigenetic modifications, gene copy numbers, chromatin interactions, and deletion mutations. In addition to effects by these factors, ZmDRG expression is influenced by promoter elements, cis-acting eQTLs, and co-expression networks.

Decrypting tubby-like protein gene family of multiple functions in starch root crop cassava.

Tubby-like proteins (TLPs) are ubiquitous in eukaryotes and function in abiotic stress tolerance of some plants. Cassava (Manihot esculenta Crantz) is a high-yield starch root crop and has a high tolerance to poor soil conditions and abiotic stress. However, little is known about TLP gene characteristics and their expression in cassava. We identified cassava TLP genes, MeTLPs, and further analysed structure, duplication, chromosome localization and collinearity, cis-acting elements in the promoter regions and expression patterns of MeTLPs, and three-dimensional structure of the encoded proteins MeTLPs. In conclusion, there is a MeTLP family containing 13 members, which are grouped into A and C subfamilies. There are 11 pairs of MeTLPs that show the duplication which took place between 10.11 and 126.69 million years ago. Two MeTLPs 6 and 9 likely originate from one gene in an ancestral species, may be common ancestors for other MeTLPs and would most likely not be eligible for ubiquitin-related protein degradation because their corresponding proteins (MeTLPs 6 and 9) have no the F-box domain in the N-terminus. MeTLPs feature differences in the number from TLPs in wheat, apple, Arabidopsis, poplar and maize, and are highlighted by segmental duplication but more importantly by the chromosomal collinearity with potato StTLPs. MeTLPs are at least related to abiotic stress tolerance in cassava. However, the subtle differences in function among MeTLPs are predictable partly because of their differential expression profiles, which are coupled with various cis­acting elements existing in the promoter regions depending on genes.

Cassava AGPase genes and their encoded proteins are different from those of other plants.

MAIN CONCLUSION: Cassava AGPase and AGPase genes have some unique characteristics. ADP-glucose pyrophosphorylase (AGPase) is a rate-limiting enzyme for starch synthesis. In this study, cassava AGPase genes (MeAGP) were analyzed based on six cultivars and one wild species. A total of seven MeAGPs was identified, including four encoding AGPase large subunits (MeAGPLs 1, 2, 3 and 4) and three encoding AGPase small subunits (MeAGPSs 1, 2 and 3). The copy number of MeAGPs varied in cassava germplasm materials. There were 14 introns for MeAGPLs 1, 2 and 3, 13 introns for MeAGPL4, and 8 introns for other three MeAGPSs. Multiple conservative amino acid sequence motifs were found in the MeAGPs. There were differences in amino acids at binding sites of substrates and regulators among different MeAGP subunits and between MeAGPs and a potato AGPase small subunit (1YP2:B). MeAGPs were all located in chloroplasts. MeAGP expression was not only associated with gene copy number and types/combinations, regions and levels of the DNA methylation but was also affected by environmental factors with the involvement of various transcription factors in multiple regulation networks and in various cis-elements in the gene promoter regions. The MeAGP activity also changed with environmental conditions and had potential differences among the subunits. Taken together, MeAGPs differ in number from those of Arabidopsis, potato, maize, banana, sweet potato, and tomato.