Genome-wide investigation of glycoside hydrolase 9 (GH9) gene family unveils implications in orchestrating the mastication trait of Citrus sinensis fruits.

Mastication trait of citrus significantly influences the fruit's overall quality and consumer preference. The accumulation of cellulose in fruits significantly impacts the mastication trait of citrus fruits, and the glycoside hydrolase 9 (GH9) family plays a crucial role in cellulose metabolism. In this study, we successfully identified 32 GH9 genes from the Citrus sinensis genome and subsequently conducted detailed bioinformatics analyses of the GH9 family. Additionally, we profiled the spatiotemporal expression patterns of CsGH9 genes across four distinct fruit tissue types and six crucial developmental stages of citrus fruits, leveraging transcriptome data. Parallel to this, we undertook a comparative analysis of transcriptome profiles and cellulose content among diverse fruit tissues spanning six developmental stages. Furthermore, to identify the pivotal genes involved in cellulose metabolism within the GH9 family during fruit maturity, we employed correlation analysis between cellulose content and gene expression in varying tissues across diverse citrus varieties. This analysis highlighted key genes such as CsGH9A2/6 and CsGH9B12/13/14/22. Collectively, this study provides an in-depth analysis of the GH9 gene family in citrus and offers novel molecular insights into the underlying mechanisms governing the mastication trait formation in citrus fruits.

Role and molecular mechanisms of HuangQiSiJunZi decoction for treating triple-negative breast cancer as explored via network pharmacology and bioinformatics analyses.

OBJECTIVE: In this study, we evaluated the molecular mechanisms of HuangQiSiJunZi Decoction (HQSJZD) for treating triple-negative breast cancer (TNBC) using network pharmacology and bioinformatics analyses. METHODS: Effective chemical components together with action targets of HQSJZD were selected based on the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP). Meanwhile, differentially expressed genes (DEGs) were extracted from TNBC sample data in The Cancer Genome Atlas (TCGA) database. Additionally, we built a protein-protein interaction (PPI) network and acquired hub genes. Gene Expression Omnibus(GEO) datasets were utilized to verify the accuracy of hub gene expression. Additionally, enrichment analyses were conducted on key genes. Furthermore, TNBC severity-related high-risk factors were screened through univariate together with multivariate Cox regressions; next, the logistic regression prediction model was built. Moreover, differential levels of 22 immune cell types in TNBC tissues compared with normal tissues were analyzed. The hub gene levels within pan-cancer and the human body were subsequently visualized and analyzed. Finally, quantitative PCR (RT-qPCR) was used to validate the correlation of the hub genes in TNBC cells. RESULTS: The study predicted 256 targets of active ingredients and 1791 DEGs in TNBC, and obtained 16 hub genes against TNBC. The prognostic signature based on FOS, MMP9, and PGR was independent in predicting survival. A total of seven types of immune cells, such as CD4 + memory T cells, showed a significant difference in infiltration (p < 0.05), and immune cells were related to the hub genes. The HPA database was adopted for hub gene analyses, and as determined, FOS was highly expressed in most human organs. The results of RT-qPCR validation for the FOS hub gene were consistent with those of bioinformatic analyses. CONCLUSION: HQSJZD might regulate the interleukin-17 and aging pathways via FOS genes to increase immune cell infiltration in TNBC tissues, and thus, may treat TNBC and improve the prognosis. The FOS genes are likely to be a new marker for TNBC.

SIRT1 and ZNF350 as novel biomarkers for osteoporosis: a bioinformatics analysis and experimental validation.

BACKGROUND: Osteoporosis (OP) is characterized by bone mass decrease and bone tissue microarchitectural deterioration in bone tissue. This study identified potential biomarkers for early diagnosis of OP and elucidated the mechanism of OP. METHODS: Gene expression profiles were downloaded from Gene Expression Omnibus (GEO) for the GSE56814 dataset. A gene co-expression network was constructed using weighted gene co-expression network analysis (WGCNA) to identify key modules associated with healthy and OP samples. Functional enrichment analysis was conducted using the R clusterProfiler package for modules to construct the transcriptional regulatory factor networks. We used the "ggpubr" package in R to screen for differentially expressed genes between the two samples. Gene set variation analysis (GSVA) was employed to further validate hub gene expression levels between normal and OP samples using RT-PCR and immunofluorescence to evaluate the potential biological changes in various samples. RESULTS: There was a distinction between the normal and OP conditions based on the preserved significant module. A total of 100 genes with the highest MM scores were considered key genes. Functional enrichment analysis suggested that the top 10 biological processes, cellular component and molecular functions were enriched. The Toll-like receptor signaling pathway, TNF signaling pathway, PI3K-Akt signaling pathway, osteoclast differentiation, JAK-STAT signaling pathway, and chemokine signaling pathway were identified by Kyoto Encyclopedia of Genes and Genomes pathway analysis. SIRT1 and ZNF350 were identified by Wilcoxon algorithm as hub differentially expressed transcriptional regulatory factors that promote OP progression by affecting oxidative phosphorylation, apoptosis, PI3K-Akt-mTOR signaling, and p53 pathway. According to RT-PCR and immunostaining results, SIRT1 and ZNF350 levels were significantly higher in OP samples than in normal samples. CONCLUSION: SIRT1 and ZNF350 are important transcriptional regulatory factors for the pathogenesis of OP and may be novel biomarkers for OP treatment.

Integrating bioinformatic resources to identify characteristics of rheumatoid arthritis-related usual interstitial pneumonia.

BACKGROUND: Rheumatoid arthritis (RA) is often accompanied by a common extra-articular manifestation known as RA-related usual interstitial pneumonia (RA-UIP), which is associated with a poor prognosis. However, the mechanism remains unclear. To identify potential mechanisms, we conducted bioinformatics analysis based on high-throughput sequencing of the Gene Expression Omnibus (GEO) database. RESULTS: Weighted gene co-expression network analysis (WGCNA) analysis identified 2 RA-positive related modules and 4 idiopathic pulmonary fibrosis (IPF)-positive related modules. A total of 553 overlapped differentially expressed genes (DEG) were obtained, of which 144 in the above modules were further analyzed. The biological process of "oxidative phosphorylation" was found to be the most relevant with both RA and IPF. Additionally, 498 up-regulated genes in lung tissues of RA-UIP were screened out and enriched by 7 clusters, of which 3 were closely related to immune regulation. The analysis of immune infiltration showed a characteristic distribution of peripheral immune cells in RA-UIP, compared with IPF-UIP in lung tissues. CONCLUSIONS: These results describe the complex molecular and functional landscape of RA-UIP, which will help illustrate the molecular pathological mechanism of RA-UIP and identify new biomarkers and therapeutic targets for RA-UIP in the future.

gjSOX9 Cloning, Expression, and Comparison with gjSOXs Family Members in Gekko japonicus.

SOX9 plays a crucial role in the male reproductive system, brain, and kidneys. In this study, we firstly analyzed the complete cDNA sequence and expression patterns for SOX9 from Gekko japonicus SOX9 (gjSOX9), carried out bioinformatic analyses of physiochemical properties, structure, and phylogenetic evolution, and compared these with other members of the gjSOX family. The results indicate that gjSOX9 cDNA comprises 1895 bp with a 1482 bp ORF encoding 494aa. gjSOX9 was not only expressed in various adult tissues but also exhibited a special spatiotemporal expression pattern in gonad tissues. gjSOX9 was predicted to be a hydrophilic nucleoprotein with a characteristic HMG-Box harboring a newly identified unique sequence, "YKYQPRRR", only present in SOXE members. Among the 20 SOX9 orthologs, gjSOX9 shares the closest genetic relationships with Eublepharis macularius SOX9, Sphacrodactylus townsendi SOX9, and Hemicordylus capensis SOX9. gjSOX9 and gjSOX10 possessed identical physicochemical properties and subcellular locations and were tightly clustered with gjSOX8 in the SOXE group. Sixteen gjSOX family members were divided into six groups: SOXB, C, D, E, F, and H with gjSOX8, 9, and 10 in SOXE among 150 SOX homologs. Collectively, the available data in this study not only facilitate a deep exploration of the functions and molecular regulation mechanisms of the gjSOX9 and gjSOX families in G. japonicus but also contribute to basic research regarding the origin and evolution of SOX9 homologs or even sex-determination mode in reptiles.

Identification and Characterization of the BBX Gene Family in Bambusa pervariabilis × Dendrocalamopsis grandis and Their Potential Role under Adverse Environmental Stresses.

Zinc finger protein (ZFP) transcription factors play a pivotal role in regulating plant growth, development, and response to biotic and abiotic stresses. Although extensively characterized in model organisms, these genes have yet to be reported in bamboo plants, and their expression information is lacking. Therefore, we identified 21 B-box (BBX) genes from a transcriptome analysis of Bambusa pervariabilis × Dendrocalamopsis grandis. Consequently, multiple sequence alignments and an analysis of conserved motifs showed that they all had highly similar structures. The BBX genes were divided into four subgroups according to their phylogenetic relationships and conserved domains. A GO analysis predicted multiple functions of the BBX genes in photomorphogenesis, metabolic processes, and biological regulation. We assessed the expression profiles of 21 BBX genes via qRT-PCR under different adversity conditions. Among them, eight genes were significantly up-regulated under water deficit stress (BBX4, BBX10, BBX11, BBX14, BBX15, BBX16, BBX17, and BBX21), nine under salt stress (BBX2, BBX3, BBX7, BBX9, BBX10, BBX12, BBX15, BBX16, and BBX21), twelve under cold stress (BBX1, BBX2, BBX4, BBX7, BBX10, BBX12, BBX14, BBX15, BBX17, BBX18, BBX19, and BBX21), and twelve under pathogen infestation stress (BBX1, BBX2, BBX4, BBX7, BBX10, BBX12, BBX14, BBX15, BBX17, BBX18, BBX19, and BBX21). Three genes (BBX10, BBX15, and BBX21) were significantly up-regulated under both biotic and abiotic stresses. These results suggest that the BBX gene family is integral to plant growth, development, and response to multivariate stresses. In conclusion, we have comprehensively analyzed the BDBBX genes under various adversity stress conditions, thus providing valuable information for further functional studies of this gene family.

The molecular mechanisms and environmental effects of modified clay control algal blooms in aquacultural water.

Modified clay (MC) technology is an effective method for controlling harmful algal blooms (HABs). Based on field experience, a bloom does not continue after treatment with MC, even though the residual HAB biomass accounts for 20-30% of the initial biomass. Laboratory studies using unialgal cultures have found that MC could inhibit the growth of the residual algal cells to prevent HABs. Nevertheless, the phytoplankton in field waters is diverse. Therefore, unclassified complex mechanisms may exist. To illustrate the molecular mechanisms through which MC controls HABs in the field and verify the previous laboratory findings, a series of experiments and bioinformatics analyses were conducted using bloom waters from aquacultural ponds. The results showed that a 72.29% removal efficiency of algal biomass could effectively control blooms. The metatranscriptomic results revealed that the number of downregulated genes (131,546) was greater than that of upregulated genes (24,318) at 3 h after MC addition. Among these genes, several genes related to DNA replication were downregulated; however, genes involved in DNA repair were upregulated. Metabolism-related pathways were the most significantly upregulated (q < 0.05), including photosynthesis and oxidative phosphorylation. The results also showed that MC reduced most of the biomass of the dominant phytoplankton species, likely by removing apical dominance, which increased the diversity and stability of the phytoplankton community. In addition to reducing the pathogenic bacterial density, MC reduced the concentrations of PO43- (96.22%) and SiO32- (66.77%), thus improving the aquaculture water quality, altering the phytoplankton community structure (the proportion of Diatomea decreased, and that of Chlorophyta increased), and inhibiting phytoplankton growth. These effects hindered the rapid development of large phytoplankton biomasses and allowed the community structure to remain stable, reducing HAB threats. This study illustrates the molecular mechanisms through which MC controls HABs in the field and provides a scientific method for removing HABs in aquacultural waters.

Exploring the potential mechanism of emetine against coronavirus disease 2019 combined with lung adenocarcinoma: bioinformatics and molecular simulation analyses.

BACKGROUND: Patients with lung adenocarcinoma (LUAD) may be more predisposed to coronavirus disease 2019 (COVID-19) and have a poorer prognosis. Currently, there is still a lack of effective anti-LUAD/COVID-19 drugs. Thus, this study aimed to screen for an effective anti-LUAD/COVID-19 drug and explore the potential mechanisms. METHODS: Firstly, we performed differentially expressed gene (DEG) analysis on LUAD transcriptome profiling data in The Cancer Genome Atlas (TCGA), where intersections with COVID-19-related genes were screened out. Then, we conducted Cox proportional hazards analyses on these LUAD/COVID-19 DEGs to construct a risk score. Next, LUAD/COVID-19 DEGs were uploaded on Connectivity Map to obtain drugs for anti-LUAD/COVID-19. Finally, we used network pharmacology, molecular docking, and molecular dynamics (MD) simulation to explore the drug's therapeutic targets and potential mechanisms for anti-LUAD/COVID-19. RESULTS: We identified 230 LUAD/COVID-19 DEGs and constructed a risk score containing 7 genes (BTK, CCL20, FURIN, LDHA, TRPA1, ZIC5, and SDK1) that could classify LUAD patients into two risk groups. Then, we screened emetine as an effective drug for anti-LUAD/COVID-19. Network pharmacology analyses identified 6 potential targets (IL6, DPP4, MIF, PRF1, SERPING1, and SLC6A4) for emetine in anti-LUAD/COVID-19. Molecular docking and MD simulation analyses showed that emetine exhibited excellent binding capacities to DDP4 and the main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). CONCLUSIONS: This study found that emetine may inhibit the entry and replication of SARS-CoV-2 and enhance tumor immunity by bounding to DDP4 and Mpro.

Identification and validation of the N6-methyladenosine RNA methylation regulator ZC3H13 as a novel prognostic marker and potential target for hepatocellular carcinoma.

N6-methyladenosine (m6A) RNA methylation has been implicated in various malignancies. This study aimed to identify prognostic signature based on m6A methylation regulators for hepatocellular carcinoma (HCC) and provide candidate targets for HCC treatment. In this study, the expression levels, prognostic values, correlation with tumor grades and genetic variations of m6A-related genes in HCC were evaluated using bioinformatics analyses. Interestingly, the results show that methyltransferase zinc finger CCCH-type containing 13 (ZC3H13) was expressed at a significantly low level in HCC. Survival outcome analysis suggested that significant correlations existed between ZC3H13 downregulation and poor overall survival (OS) and poor recurrence-free survival (RFS) in HCC patients. Therefore, ZC3H13 was chosen for further experimental validation. The expression of ZC3H13 in HCC cell lines was investigated by western blotting. Knockdown of ZC3H13 significantly enhanced the migration and invasion of HCC cells, as demonstrated by wound healing and transwell assays. Moreover, upregulating ZC3H13 repressed the growth of xenograft tumors in vivo. Functional and pathway enrichment analyses indicated that ZC3H13 might be involved in transcriptional dysregulation or the JAK-STAT signaling pathway in cancer. Additionally, ZC3H13 expression was significantly correlated with lymphocytes and immunomodulators. Therefore, ZC3H13 is a promising candidate as a novel biomarker and therapeutic target for HCC.

Silencing oncogene cell division cycle associated 5 induces apoptosis and G1 phase arrest of non-small cell lung cancer cells via p53-p21 signaling pathway.

BACKGROUNDS: As a regulator of cell cycle, cell division cycle-associated 5 (CDCA5) is involved in the progression of various malignant tumors. However, the potential relationship between CDCA5 and lung cancer has not been reported. METHODS: In our study, we analyzed the expression of CDCA5 in a variety of malignant tumors, performed Kaplan-Meier survival analysis of lung adenocarcinoma (LUAD), explored the potential relationship between CDCA5 expression and clinicopathological characteristics, assessed the predictive capability of at different stages of clinicopathological characteristics, revealed the enriched functions and signaling pathways among LUAD paitents with high CDCA5 expression, and investigated the correlation between PD-1, PD-L1, and CDCA5 through bioinformatics analyses. Subsequently, we performed quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) and western blotting (WB) to demonstrate that CDCA5 mediates the p53-p21 pathway and regulates the cell cycle. RESULT: CDCA5 is probably involved in the occurrence and development of NSCLC, and function as a reliable biomarker for predicting the survival outcomes of patients with early stage of patients with LUAD. Furthermore, CDCA5 may be a promising indicator of immunotherapy efficacy. In addition, silencing the expression of CDCA5 significantly increased the proportion of apoptotic NSCLC cells, and caused NSCLC cells to be arrested in the G1 phase. CONSLUSION: In conclusion, CDCA5 regulated the cell cycle of NSCLC cells by mediating the p53-p21 signaling pathway, participating in the development and progression of NSCLC patients.

Identification of GINS1 as a therapeutic target in the cancer patients infected with COVID-19: a bioinformatics and system biology approach.

BACKGROUND: Coronavirus disease 2019 (COVID-19) caused a series of biological changes in cancer patients which have rendered the original treatment ineffective and increased the difficulty of clinical treatment. However, the clinical treatment for cancer patients infected with COVID-19 is currently unavailable. Since bioinformatics is an effective method to understand undiscovered biological functions, pharmacological targets, and therapeutic mechanisms. The aim of this study was to investigate the influence of COVID-19 infection in cancer patients and to search the potential treatments. METHODS: Firstly, we obtained the COVID-19-associated genes from seven databases and analyzed the cancer pathogenic genes from Gene Expression Omnibus (GEO) databases, respectively. The Cancer/COVID-19-associated genes were shown by Venn analyses. Moreover, we demonstrated the signaling pathways and biological functions of pathogenic genes in Cancer/COVID-19. RESULTS: We identified that Go-Ichi-Ni-San complex subunit 1 (GINS1) is the potential therapeutic target in Cancer/COVID-19 by GEPIA. The high expression of GINS1 was not only promoting the development of cancers but also affecting their prognosis. Furthermore, eight potential compounds of Cancer/COVID-19 were identified from CMap and molecular docking analysis. CONCLUSION: We revealed the GINS1 is a potential therapeutic target in cancer patients infected with COVID-19 for the first time, as COVID-19 will be a severe and prolonged pandemic. However, the findings have not been verified actually cancer patients infected with COVID-19, and further studies are needed to demonstrate the functions of GINS1 and the clinical treatment of the compounds.

Expression profiles and functional prediction of long non-coding RNAs LINC01133, ZEB1-AS1 and ABHD11-AS1 in the luminal subtype of breast cancer.

BACKGROUND: Luminal breast cancer (BC) is the most frequent subtype accounting for more than 70% of BC. LncRNAs, a class of non-coding RNAs with more than 200 nucleotides, are involved in a variety of cellular processes and biological functions. Abberant expression is related to the development of various cancers, such as breast cancer. LINC01133, ZEB1-AS1, and ABHD11-AS1 were reported to be dysregulated in different cancers. However, their expression level in luminal BC remains poorly known. The aim of the present study was to evaluate the potential roles of these lncRNAs in BC, especially in luminal subtypes. METHODS: A comprehensive analysis was performed using the Lnc2Cancer database to identify novel cancer-associated lncRNA candidates. After conducting a literature review, three novel lncRNAs named LINC01133, ZEB1-AS1, and ABHD11-AS1 were chosen as target genes of the present study. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to evaluate the expression level of the mentioned lncRNAs in both luminal BC tissues and cell lines. Then, the correlation of the three mentioned lncRNAs expression with clinicopathological characteristics of the patients was studied. Moreover, several datasets were used to discover the potential roles and functions of LINC01133, ZEB1-AS1 and ABHD11-AS1 in luminal subtype of BC. RESULTS: According to the qRT-PCR assay, the expression levels of LINC01133 and ZEB1-AS1 were decreased in luminal BC tissues and cell lines. On the other hand, ABHD11-AS1 was upregulated in the above-mentioned samples. The expression levels of LINC01133, ZEB1-AS1, and ABHD11-AS1 were not associated with any of the clinical features. Also, the results obtained from the bioinformatics analyses were consistent with qRT-PCR data. Functional annotation of the co-expressed genes with the target lncRNAs, protein-protein interactions and significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways across luminal BC were also obtained using bioinformatics analysis. CONCLUSIONS: Taken together, our findings disclosed the dysregulation of LINC01133, ZEB1-AS1, and ABHD11-AS1 in luminal BC. It was revealed that LINC01133 and ZEB1-AS1 expression was significantly downregulated in luminal BC tissues and cell lines, while ABHD11-AS1 was upregulated considerably in the mentioned tissues and cell lines. Also, bioinformatics and systems biology analyses have helped to identify the possible role of these lncRNAs in luminal BC. However, further analysis is needed to confirm the current findings.

The impact of PSRC1 overexpression on gene and transcript expression profiling in the livers of ApoE-/- mice fed a high-fat diet.

Our previous studies have confirmed that proline/serine-rich coiled-coil 1 (PSRC1) overexpression can regulate blood lipid levels and inhibit atherosclerosis (AS) development. In the current study, the gene and transcript expression profiles in the livers of ApoE-/- mice overexpressing PSRC1 were investigated. HiSeq X Ten RNA sequencing (RNA-seq) analysis was used to examine the differentially expressed genes (DEGs) and differentially expressed transcripts in the livers of PSRC1-overexpressing ApoE-/- and control mice. Then, Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on these DEGs and on long noncoding RNA (lncRNA) predicted target genes. A total of 1892 significant DEGs were identified: 1431 were upregulated (e.g., Cyp2a4, Obp2a, and Sertad4), and 461 were downregulated (e.g., Moxd1, Egr1, and Elovl3). In addition, 8184 significant differentially expressed transcripts were identified, 4908 of which were upregulated and 3276 of which were downregulated. Furthermore, 1106 significant differentially expressed lncRNAs were detected, 713 of which were upregulated and 393 of which were downregulated. Quantitative reverse transcription PCR (qRT-PCR) verified changes in 10 randomly selected DEGs. GO analyses showed that the DEGs and predicted lncRNA target genes were mostly enriched for actin binding and lipid metabolism. KEGG biological pathway analyses showed that the DEGs in the livers of PSRC1-overexpressing ApoE-/- mice were enriched in the mitogen-activated protein kinase (MAPK) pathway. These findings reveal that PSRC1 may affect liver actin polymerization and cholesterol metabolism-related genes or pathways. These mRNAs and lncRNAs may represent new biomarkers and targets for the diagnosis and therapy of lipid metabolism disturbance and AS.

Soil Microbial Community Responses After Amendment with Thermally Altered Pinus radiata Needles.

Post-fire litter layers are composed of leaves and woody debris that predominantly fall during or soon after the fire event. These layers are distinctly different to pre-fire litters due to their common origin and deposition time. However, heterogeneity can arise from the variable thermal conditions in the canopy during fire. Therefore, in this study, we used thermally altered pine needles (heated to 40 °C, 150 °C, 260 °C and 320 °C for 1 h) in a laboratory incubation study for 43 days. These samples were measured for respiration throughout and extracted for DNA at the experiment's end; soil ribosomal RNA was analysed using Illumina sequencing (16S and internal transcribed spacer amplicons). The addition of pine needles heated to 40 °C or 150 °C caused a substantial shift in community structure, decreased alpha diversity and significantly increased soil respiration relative to the control treatment. In contrast, pine needles heated to 260 °C or 320 °C had little effect on microbial community structure or soil respiration. These results indicate that highly thermally altered needles are not microbially decomposed during the first 43 days of exposure and therefore that biomass temperature may have significant effects on post-fire litter decomposition and carbon flux. This research outlines an important knowledge gap in forest fire responses that may affect post-fire carbon emission estimates.

Integrated bioinformatics analyses identified SCL3-induced regulatory network in Arabidopsis thaliana roots.

OBJECTIVES: To uncover key genes and pathways regulated by SCL3, a GRAS transcription factor, in the context of gibberellin (GA) in the roots of the model plant Arabidopsis thaliana. RESULTS: Gene expression profiles of ga1-3 mutant and ga1-3 and scl3 double mutant are considerably similar to each other, revealed by Principal Component Analysis (PCA). More than 400 significantly Differentially Expressed Genes (DEGs) among the Arabidopsis thaliana roots of ga1-3 mutant, ga1-3 and scl3 double mutant and GA loss/SCL3 gain mutant were uncovered by comprehensive bioinformatics analyses. Protein synthesis pathway, including RPL proteins, RPS proteins, etc., and flavonoid biosynthesis pathway, including TT4, F3H, TT5, CHIL, etc. were significantly increased when SCL3 expression was higher than normal by means of pathway enrichment analysis and protein-protein interaction analysis, which is further supported by comparison analyses between wild type samples and SCL3 overexpressed roots. CONCLUSION: Protein synthesis and flavonoid biosynthesis were regulated by SCL3 in the context of GA in Arabidopsis thaliana root system identified by comprehensive bioinformatic analyses.

Identification of MMP1 as a potential gene conferring erlotinib resistance in non-small cell lung cancer based on bioinformatics analyses.

BACKGROUND: Non-small cell lung cancer (NSCLC) is the major type of lung cancer with high morbidity and poor prognosis. Erlotinib, an inhibitor of epidermal growth factor receptor (EGFR), has been clinically applied for NSCLC treatment. Nevertheless, the erlotinib acquired resistance of NSCLC occurs inevitably in recent years. METHODS: Through analyzing two microarray datasets, erlotinib resistant NSCLC cells microarray (GSE80344) and NSCLC tissue microarray (GSE19188), the differentially expressed genes (DEGs) were screened via R language. DEGs were then functionally annotated by Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, which up-regulated more than 2-folds in both datasets were further functionally analyzed by Oncomine, GeneMANIA, R2, Coremine, and FunRich. RESULTS: We found that matrix metalloproteinase 1 (MMP1) may confer the erlotinib therapeutic resistance in NSCLC. MMP1 highly expressed in erlotinib-resistant cells and NSCLC tissues, and it associated with poor overall survival. In addition, MMP1 may be associated with COPS5 and be involve in an increasing transcription factors HOXA9 and PBX1 in erlotinib resistance. CONCLUSIONS: Generally, these results demonstrated that MMP1 may play a crucial role in erlotinib resistance in NSCLC, and MMP1 could be a prognostic biomarker for erlotinib treatment.

Annexin A8 can serve as potential prognostic biomarker and therapeutic target for ovarian cancer: based on the comprehensive analysis of Annexins.

BACKGROUND: Annexins are involved in vesicle trafficking, cell proliferation and apoptosis, but their functional mechanisms in ovarian cancer remain unclear. In this study, we analyzed Annexins in ovarian cancer using different databases and selected Annexin A8 (ANXA8), which showed the greatest prognostic value, for subsequent validation in immunohistochemical (IHC) assays. METHODS: The mRNA expression levels, genetic variations, prognostic values and gene-gene interaction network of Annexins in ovarian cancer were analyzed using the Oncomine, Gene Expression Profiling Interactive Analysis (GEPIA), cBioPortal, Kaplan-Meier plotter and GeneMANIA database. ANXA8 was selected for analyzing the biological functions and pathways of its co-expressed genes, and its correlation with immune system responses via the Database for Annotation, Visualization, and Integrated Discovery (DAVID) and the TISIDB database, respectively. We validated the expression of ANXA8 in ovarian cancer via IHC assays and analyzed its correlation with clinicopathological parameters and prognosis. RESULTS: ANXA2/3/8/11 mRNA expression levels were significantly upregulated in ovarian cancer, and ANXA5/6/7 mRNA expression levels were significantly downregulated. Prognostic analysis suggested that significant correlations occurred between ANXA2/4/8/9 mRNA upregulation and poor overall survival, and between ANXA8/9/11 mRNA upregulation and poor progression-free survival in patients with ovarian serous tumors. Taken together, results suggested that ANXA8 was most closely associated with ovarian cancer tumorigenesis and progression. Further analyses indicated that ANXA8 may be involved in cell migration, cell adhesion, and vasculature development, as well as in the regulation of PI3K-Akt, focal adhesion, and proteoglycans. Additionally, ANXA8 expression was significantly correlated with lymphocytes and immunomodulators. The IHC results showed that ANXA8 expression was higher in the malignant tumor group than in the borderline and benign tumor groups and normal ovary group, and high ANXA8 expression was an independent risk factor for survival and prognosis of ovarian cancer patients (P = 0.013). CONCLUSIONS: Members of the Annexin family display varying degrees of abnormal expressions in ovarian cancer. ANXA8 was significantly highly expressed in ovarian cancer, and high ANXA8 expression was significantly correlated with poor prognosis. Therefore, ANXA8 is a high candidate as a novel biomarker and therapeutic target for ovarian cancer.

CircRNA expression profile and functional analysis in testicular tissue of patients with non-obstructive azoospermia.

BACKGROUND: Non-obstructive azoospermia (NOA) is a multifactorial disorder whose molecular basis remains largely unknown. Circular RNAs (CircRNAs), a novel class of endogenous RNAs, have been recognized to play important roles in many biological processes. However, little is known about the expression patterns and functions of circRNAs in human testes involved in NOA. METHODS: In this study, the testicular circRNA expression profile were explored in NOA patients and the controls by high-throughput circRNA microarray. Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to confirm the microarray data. Bioinformatics analyses including the circRNA/miRNA/mRNA interaction network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used to predict the functions of differentially expressed circRNAs. RESULTS: A total of 368 differentially down-regulated and 526 up-regulated circRNAs were detected in NOA patients. These findings have been verified by qRT-PCR on 6 selected circRNAs. Among these differentially expressed circRNAs, the hsa_circRNA_0023313 was obviously up-regulated in testicular tissue of NOA patients. The most likely potential target miRNA for hsa_circRNA_0023313 include hsa-miR-520d-3p, hsa-miR-373-3p, hsa-miR-372-3p, hsa-miR-302c-3p and hsa-miR-130b-5p. Function analysis indicated that hsa_circRNA_0023313 was ubiquitin-protein transferase activity and chromatin binding. KEGG analysis revealed that the top five pathways related to hsa_circRNA_0023313 were endocytosis, meiosis, FoxO signaling pathway, ubiquitin mediated proteolysis and AMPK signaling pathway. CONCLUSIONS: This is the first report that the testicular circRNA expression profile is altered in NOA patients indicating circRNAs might play important roles in regulating spermatogenesis and be potential biomarkers for the diagnosis and treatment of NOA.

LncRNA Expression Profile of Human Thoracic Aortic Dissection by High-Throughput Sequencing.

BACKGROUND/AIMS: In this study, the long non-coding RNA (lncRNA) expression profile in human thoracic aortic dissection (TAD), a highly lethal cardiovascular disease, was investigated. METHODS: Human TAD (n=3) and normal aortic tissues (NA) (n=3) were examined by high-throughput sequencing. Bioinformatics analyses were performed to predict the roles of aberrantly expressed lncRNAs. Quantitative real-time polymerase chain reaction (qRT-PCR) was applied to validate the results. RESULTS: A total of 269 lncRNAs (159 up-regulated and 110 down-regulated) and 2, 255 mRNAs (1 294 up-regulated and 961 down-regulated) were aberrantly expressed in human TAD (fold-change> 1.5, P< 0.05). QRT-PCR results of five dysregulated genes were consistent with HTS data. A lncRNA-mRNA coexpression analysis showed positive correlations between the up-regulated lncRNA (ENSG00000269936) and its adjacent up-regulated mRNA (MAP2K6, R=0.940, P< 0.01), and between the down-regulated lncRNA_1421 and its down-regulated mRNAs (FBLN5, R=0.950, P< 0.01; ACTA2, R=0.96, P< 0.01; TIMP3, R=0.96, P< 0.05). The lncRNA-miRNA-mRNA network indicated that the up-regulated lncRNA XIST and p21 had similar sequences targeted by has-miR-17-5p. The results of luciferase assay and fluorescence immuno-cytochemistry were consistent with that. And qRT-PCR results showed that lncRNA XIST and p21 were expressed at a higher level and has-miR-17-5p was expressed at a lower level in TAD than in NA. The predicted binding motifs of three up-regulated lncRNAs (ENSG00000248508, ENSG00000226530, and EG00000259719) were correlated with up-regulated RUNX1 (R=0.982, P< 0.001; R=0.967, P< 0.01; R=0.960, P< 0.01, respectively). CONCLUSIONS: Our study revealed a set of dysregulated lncRNAs and predicted their multiple potential functions in human TAD. These findings suggest that lncRNAs are novel potential therapeutic targets for human TAD.

Identifying and engineering a critical amino acid residue to enhance the catalytic efficiency of Pseudomonas sp. methyl parathion hydrolase.

Methyl parathion hydrolase (MPH) that hydrolyzes a wide range of organophosphorus pesticides can be used to remediate land polluted by the pesticides. Here, the catalytic efficiency of methyl parathion hydrolase from Pseudomonas sp. (WBC-3) was enhanced by searching and engineering a critical site far away from the binding pocket. In the first round, a four-site mutant with a modest increased catalytic efficiency (3.2-fold kcat/Km value of the wild type) was obtained with random mutagenesis. By splitting and re-combining the four substitutions in the mutant, the critical site S277, was identified to show the most significant effects of improving binding affinity and catalytic efficiency. With further site-saturation mutagenesis focused on the residue S277, another two substitutions were discovered to have even more significant decrease in Km (40.2 and 47.6 µM) and increased in kcat/Km values (9.5- and 10.3-fold of the wild type) compared to the original four-site mutant (3.0- and 3.2-fold). In the three-dimensional structure, residue S277 is located at a hinge region of a loop, which could act as a "lid" at the substrate entering to the binding pocket. This suggests that substitutions of residue S277 could affect substrate binding via conformational change in substrate entrance region. This work provides a valuable protocol combining random mutagenesis, site-saturation mutagenesis, structural and bioinformatics analyses to obtain mutants with high catalytic efficiency from a screening library of a modest size (3200 strains).