Identification and Analysis of Blood Gene Expression Signature for Osteoarthritis With Advanced Feature Selection Methods

Transcriptomic analysis of the HPT axis in a model of oligoasthenozoospermia induced by Adenine in rats

Male infertility is most commonly caused by oligozoospermia, and its pathogenesis is still poorly understood at the molecular level. This study used RNA sequencing (RNA-Seq) technology to identify candidate genes and regulatory pathways that regulate semen quality in the hypothalamic, pituitary, and testicular tissues of healthy rats and Adenine-induced oligozoospermia model rats. Semen quality testing and histological analysis of testicular tissues were performed on both groups of rats. We identified 627, 692, and 437 differentially expressed genes in the hypothalamus, pituitary gland, and testes, respectively. Functional analysis indicates that "neuronal projections," "positive regulation of hormone biosynthetic process," and "neuroactive ligand-receptor interaction pathways" are closely related to the hypothalamic-pituitary-testicular (HPT) axis hormone regulation and sperm production. Seven genes (Pomc, Rxfp1, Tac1, Npy, Insl3, Hsd3b3, Lhcgr) have been identified as key candidate genes responsible for regulating sperm quality within the HPT axis, potentially affecting rat reproductive function by influencing testicular development and testosterone secretion. These data provide a theoretical basis for further understanding the molecular mechanisms of reproductive performance in a rat model of oligoasthenozoospermia.

Peripheral Blood DNA Methylation-Based Machine Learning Models for Prediction of Knee Osteoarthritis Progression: Biologic Specimens and Data From the Osteoarthritis Initiative and Johnston County Osteoarthritis Project

OBJECTIVE

The lack of accurate biomarkers to predict knee osteoarthritis (OA) progression is a key unmet need in OA clinical research. The objective of this study was to develop baseline peripheral blood epigenetic biomarker models to predict knee OA progression.

METHODS

Genome-wide buffy coat DNA methylation patterns from 554 individuals from the Osteoarthritis Biomarkers Consortium (OABC) were determined using Illumina Infinium MethylationEPIC 850K arrays. Data were divided into model development and validation sets, and machine learning models were trained to classify future OA progression by knee pain, radiographic imaging, knee pain plus radiographic imaging, and any progression (pain, radiographic, or both). Parsimonious models using the top 13 CpG sites most frequently selected during development were tested on independent samples from participants in the Johnston County Osteoarthritis (JoCo OA) Project (n = 128) and a previously published Osteoarthritis Initiative (OAI) data set (n = 55).

RESULTS

Full models accurately classified future radiographic-only progression (mean ± SEM accuracy 87 ± 0.8%, area under the curve [AUC] 0.94 ± 0.004), pain-only progression (accuracy 89 ± 0.9%, AUC 0.97 ± 0.004), pain plus radiographic progression (accuracy 72 ± 0.7%, AUC 0.79 ± 0.006), and any progression (accuracy 78 ± 0.4%, AUC 0.86 ± 0.004). Pain-only and radiographic-only progressors were not distinguishable (mean ± SEM accuracy 58 ± 1%, AUC 0.62 ± 0.001). Parsimonious models showed similar performance and accurately classified future radiographic progressors in the OABC cohort and in both validation cohorts (mean ± SEM accuracy 80 ± 0.3%, AUC 0.88 ± 0.003 [using JoCo OA Project data], accuracy 80 ± 0.8%, AUC 0.89 ± 0.002 [using previous OAI data]).

CONCLUSION

Our data suggest that pain and structural progression share similar early systemic immune epigenotypes. Further studies should focus on evaluating the pathophysiologic consequences of differential DNA methylation and peripheral blood cell epigenotypes in individuals with knee OA.

CT- and MRI-Based 3D Reconstruction of Knee Joint to Assess Cartilage and Bone

For the observation of human joint cartilage, X-ray, computed tomography (CT) or magnetic resonance imaging (MRI) are the main diagnostic tools to evaluate pathologies or traumas. The current work introduces a set of novel measurements and 3D features based on MRI and CT data of the knee joint, used to reconstruct bone and cartilages and to assess cartilage condition from a new perspective. Forty-seven subjects presenting a degenerative disease, a traumatic injury or no symptoms or trauma were recruited in this study and scanned using CT and MRI. Using medical imaging software, the bone and cartilage of the knee joint were segmented and 3D reconstructed. Several features such as cartilage density, volume and surface were extracted. Moreover, an investigation was carried out on the distribution of cartilage thickness and curvature analysis to identify new markers of cartilage condition. All the extracted features were used with advanced statistics tools and machine learning to test the ability of our model to predict cartilage conditions. This work is a first step towards the development of a new gold standard of cartilage assessment based on 3D measurements.

Prognostic Effect of Inflammatory Genes on Stage I-III Colorectal Cancer-Integrative Analysis of TCGA Data

Simple Summary

Research interest in the role of inflammation in the progression and prognosis of colorectal cancer (CRC) is growing. In this study, we evaluated the expression and DNA methylation levels of inflammation-related genes in CRC tissues using the TCGA-COREAD dataset by integratively combining multi-omics features using machine learning. Statistical analysis was additionally performed to allow for interpretable, understandable, and clinically practical results. An integrative model combining expression, methylation, and clinical features had the highest performance. In multivariate analysis, the methylation levels of CEP250, RAB21, and TNPO3 were significantly associated with overall survival. Our study results implicate the importance of integrating expression and methylation information along with clinical information in the prediction of survival. CEP250, RAB21, and TNPO3 in the prediction model might have a crucial role in CRC prognosis and further improve our understanding of potential mechanisms linking inflammatory reactions and CRC progression.

Abstract

Background inflammatory status indicators have been reported as prognostic biomarkers of colorectal cancer (CRC). However, since inflammatory interactions with the colon involve various modes of action, the biological mechanism linking inflammation and CRC prognosis has not been fully elucidated. We comprehensively evaluated the predictive roles of the expression and methylation levels of inflammation-related genes for CRC prognosis and their pathophysiological associations. Method. An integrative analysis of 247 patients with stage I-III CRC from The Cancer Genome Atlas was conducted. Lasso-penalized Cox proportional hazards regression (Lasso-Cox) and statistical Cox proportional hazard regression (CPH) were used for the analysis. Results. Models to predict overall survival were designed with respective combinations of clinical variables, including age, sex, stage, gene expression, and methylation. An integrative model combining expression, methylation, and clinical features performed better (median C-index = 0.756) than the model with clinical features alone (median C-index = 0.726). Based on multivariate CPH with features from the best model, the methylation levels of CEP250, RAB21, and TNPO3 were significantly associated with overall survival. They did not share any biological process in functional networks. The 5-year survival rate was 29.8% in the low methylation group of CEP250 and 79.1% in the high methylation group (p < 0.001). Conclusion. Our study results implicate the importance of integrating expression and methylation information along with clinical information in the prediction of survival. CEP250, RAB21, and TNPO3 in the prediction model might have a crucial role in CRC prognosis and further improve our understanding of potential mechanisms linking inflammatory reactions and CRC progression.

A four-genes based diagnostic signature for osteoarthritis

Osteoarthritis (OA) is a primary leading cause of pain and disability. However, some cases are diagnosed at the later stage which delayed the timely treatment. This study aims to identify effective diagnostic signature for OA. The mRNA profile GSE48566 including 106 blood samples of OA patients and 33 blood samples of healthy individuals was downloaded from Gene Expression Omnibus (GEO) database. The potential OA-related genes were screened by weighted gene co-expression network analysis (WGCNA). Gene ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to reveal the functions or pathways of OA-related genes using the clusterProfiler function package of R software. Key genes significantly involved in OA progression were further screened by protein-protein interaction (PPI) network. The logistic regression model and the random forest model were conducted by bringing into optimal genes selected by stepwise regression analysis, and fivefold cross validation method was used to determine their reliability. A total of 146 genes, existed in three modules and might be associated with the occurrence of OA, were screened. 15 genes were screened from the PPI network and four genes, including CCR6, CLEC7A, IL18 and SRSF2, were further optimized. Finally, a logistic regression model and a random forest model were conducted by bringing into four optimal genes, and could reliably separate OA patients from healthy subjects. Our study established two effective diagnostic models based on CCR6, CLEC7A, IL18 and SRSF2, which could reliably separate OA patients from healthy subjects.

eHSCPr discriminating the cell identity involved in endothelial to hematopoietic transition

MOTIVATION

Hematopoietic stem cells (HSCs) give rise to all blood cells and play a vital role throughout the whole lifespan through their pluripotency and self-renewal properties. Accurately identifying the stages of early HSCs is extremely important, as it may open up new prospects for extracorporeal blood research. Existing experimental techniques for identifying the early stages of HSCs development are time-consuming and expensive. Machine learning has shown its excellence in massive single-cell data processing and it is desirable to develop related computational models as good complements to experimental techniques.

RESULTS

In this study, we presented a novel predictor called eHSCPr specifically for predicting the early stages of HSCs development. To reveal the distinct genes at each developmental stage of HSCs, we compared F-score with three state-of-art differential gene selection methods (limma, DESeq2, edgeR) and evaluated their performance. F-score captured the more critical surface markers of endothelial cells and hematopoietic cells, and the area under receiver operating characteristic curve (ROC) value was 0.987. Based on SVM, the 10-fold cross-validation accuracy of eHSCpr in the independent dataset and the training dataset reached 94.84% and 94.19%, respectively. Importantly, we performed transcription analysis on the F-score gene set, which indeed further enriched the signal markers of HSCs development stages. eHSCPr can be a powerful tool for predicting early stages of HSCs development, facilitating hypothesis-driven experimental design and providing crucial clues for the in vitro blood regeneration studies.

AVAILABILITY

http://bioinfor.imu.edu.cn/ehscpr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Methylation Pattern for Knee and Hip Osteoarthritis

Osteoarthritis is one of the most prevalent chronic joint diseases for middle-aged and elderly people. But in recent years, the number of young people suffering from the disease increases quickly. It is known that osteoarthritis is a common degenerative disease caused by the combination and interaction of many factors such as natural and environmental factors. DNA methylations reflect the effects of environmental factors. Several researches on DNA methylation at specific genes in OA cartilage indicated the great potential roles of DNA methylation in OA. To systematically investigate the methylation pattern in knee and hip osteoarthritis, we analyzed the methylation profiles in cartilage of 16 OA hip samples, 19 control hip samples and 62 OA knee samples. 12 discriminative methylation sites were identified using advanced minimal Redundancy Maximal Relevance (mRMR) and Incremental Feature Selection (IFS) methods. The SVM classifier of these 12 methylation sites from genes like MEIS1, GABRG3, RXRA, and EN1, can perfectly classify the OA hip samples, control hip samples and OA knee samples evaluated with LOOCV (Leave-One Out-Cross Validation). These 12 methylation sites can not only serve as biomarker, but also provide underlying mechanism of OA.

Machine Learning of Single-Cell Transcriptome Highly Identifies mRNA Signature by Comparing F-Score Selection with DGE Analysis

Human preimplantation development is a complex process involving dramatic changes in transcriptional architecture. For a better understanding of their time-spatial development, it is indispensable to identify key genes. Although the single-cell RNA sequencing (RNA-seq) techniques could provide detailed clustering signatures, the identification of decisive factors remains difficult. Additionally, it requires high experimental cost and a long experimental period. Thus, it is highly desired to develop computational methods for identifying effective genes of development signature. In this study, we first developed a predictor called EmPredictor to identify developmental stages of human preimplantation embryogenesis. First, we compared the F-score of feature selection algorithms with differential gene expression (DGE) analysis to find specific signatures of the development stage. In addition, by training the support vector machine (SVM), four types of signature subsets were comprehensively discussed. The prediction results showed that a feature subset with 1,881 genes from the F-score algorithm obtained the best predictive performance, which achieved the highest accuracy of 93.3% on the cross-validation set. Further function enrichment demonstrated that the gene set selected by the feature selection method was involved in more development-related pathways and cell fate determination biomarkers. This indicates that the F-score algorithm should be preferentially proposed for detecting key genes of multi-period data in mammalian early development.

Identification and Analysis of Dysfunctional Genes and Pathways in CD8+ T Cells of Non-Small Cell Lung Cancer Based on RNA Sequencing

Lung cancer, the most common of malignant tumors, is typically of the non-small cell (NSCLC) type. T-cell-based immunotherapies are a promising and powerful approach to treating NSCLCs. To characterize the CD8⁺ T cells of non-small cell lung cancer, we re-analyzed the published RNA-Seq gene expression profiles of 36 CD8⁺ T cell isolated from tumor (TIL) samples and 32 adjacent uninvolved lung (NTIL) samples. With an advanced Monte Carlo method of feature selection, we identified the CD8⁺ TIL specific expression patterns. These patterns revealed the key dysfunctional genes and pathways in CD8⁺ TIL and shed light on the molecular mechanisms of immunity and use of immunotherapy.

Integrative Analysis of Methylation and Gene Expression in Lung Adenocarcinoma and Squamous Cell Lung Carcinoma

Lung cancer is a highly prevalent type of cancer with a poor 5-year survival rate of about 4-17%. Eighty percent lung cancer belongs to non-small-cell lung cancer (NSCLC). For a long time, the treatment of NSCLC has been mostly guided by tumor stage, and there has been no significant difference between the therapy strategy of lung adenocarcinoma (LUAD) and squamous cell lung carcinoma (SCLC), the two major subtypes of NSCLC. In recent years, important molecular differences between LUAD and SCLC are increasingly identified, indicating that targeted therapy will be more and more histologically specific in the future. To investigate the LUAD and SCLC difference on multi-omics scale, we analyzed the methylation and gene expression data together. With the Boruta method to remove irrelevant features and the MCFS (Monte Carlo Feature Selection) method to identify the significantly important features, we identified 113 key methylation features and 23 key gene expression features. HNF1B and TP63 were found to be dysfunctional on both methylation and gene expression levels. The experimentally determined interaction network suggested that TP63 may play an important role in connecting methylation genes and expression genes. Many of the discovered signature genes have been supported by literature. Our results may provide directions of precision diagnosis and therapy of LUAD and SCLC.

The Functional Effects of Key Driver KRAS Mutations on Gene Expression in Lung Cancer

Lung cancer is a common malignant cancer. Kirsten rat sarcoma oncogene (KRAS) mutations have been considered as a key driver for lung cancers. KRAS p.G12C mutations were most predominant in NSCLC which was comprised about 11–16% of lung adenocarcinomas (p.G12C accounts for 45–50% of mutant KRAS). But it is still not clear how the KRAS mutation triggers lung cancers. To study the molecular mechanisms of KRAS mutation in lung cancer. We analyzed the gene expression profiles of 156 KRAS mutation samples and other negative samples with two stage feature selection approach: (1) minimal Redundancy Maximal Relevance (mRMR) and (2) Incremental Feature Selection (IFS). At last, 41 predictive genes for KRAS mutation were identified and a KRAS mutation predictor was constructed. Its leave one out cross validation MCC was 0.879. Our results were helpful for understanding the roles of KRAS mutation in lung cancer.

A pilot study of peripheral blood DNA methylation models as predictors of knee osteoarthritis radiographic progression: data from the Osteoarthritis Initiative (OAI)

Knee osteoarthritis (OA) is a leading cause of chronic disability worldwide, but no diagnostic or prognostic biomarkers are available. Increasing evidence supports epigenetic dysregulation as a contributor to OA pathogenesis. In this pilot study, we investigated epigenetic patterns in peripheral blood mononuclear cells (PBMCs) as models to predict future radiographic progression in OA patients enrolled in the longitudinal Osteoarthritis Initiative (OAI) study. PBMC DNA was analyzed from baseline OAI visits in 58 future radiographic progressors (joint space narrowing at 24 months, sustained at 48 months) compared to 58 non-progressors. DNA methylation was quantified via Illumina microarrays and beta- and M-values were used to generate linear classification models. Data were randomly split into a 60% development and 40% validation subsets, models developed and tested, and cross-validated in a total of 40 cycles. M-value based models outperformed beta-value based models (ROC-AUC 0.81 ± 0.01 vs. 0.73 ± 0.02, mean ± SEM, comparison p = 0.002), with a mean classification accuracy of 73 ± 1% (mean ± SEM) for M- and 69 ± 1% for beta-based models. Adjusting for covariates did not significantly alter model performance. Our findings suggest that PBMC DNA methylation-based models may be useful as biomarkers of OA progression and warrant additional evaluation in larger patient cohorts.

Identifying Methylation Pattern and Genes Associated with Breast Cancer Subtypes

Breast cancer is regarded worldwide as a severe human disease. Various genetic variations, including hereditary and somatic mutations, contribute to the initiation and progression of this disease. The diagnostic parameters of breast cancer are not limited to the conventional protein content and can include newly discovered genetic variants and even genetic modification patterns such as methylation and microRNA. In addition, breast cancer detection extends to detailed breast cancer stratifications to provide subtype-specific indications for further personalized treatment. One genome-wide expression–methylation quantitative trait loci analysis confirmed that different breast cancer subtypes have various methylation patterns. However, recognizing clinically applied (methylation) biomarkers is difficult due to the large number of differentially methylated genes. In this study, we attempted to re-screen a small group of functional biomarkers for the identification and distinction of different breast cancer subtypes with advanced machine learning methods. The findings may contribute to biomarker identification for different breast cancer subtypes and provide a new perspective for differential pathogenesis in breast cancer subtypes.

Expression Profile Analysis Identifies a Novel Five-Gene Signature to Improve Prognosis Prediction of Glioblastoma

Glioblastoma multiforme (GBM) is the most aggressive primary central nervous system malignant tumor. The median survival of GBM patients is 12–15 months, and the 5 years survival rate is less than 5%. More novel molecular biomarkers are still urgently required to elucidate the mechanisms or improve the prognosis of GBM. This study aimed to explore novel biomarkers for GBM prognosis prediction. The gene expression profiles from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets of GBM were downloaded. A total of 2241 overlapping differentially expressed genes (DEGs) were identified from TCGA and GSE7696 datasets. By univariate COX regression survival analysis, 292 survival-related genes were found among these DEGs (p < 0.05). Functional enrichment analysis was performed based on these survival-related genes. A five-gene signature (PTPRN, RGS14, G6PC3, IGFBP2, and TIMP4) was further selected by multivariable Cox regression analysis and a prognostic model of this five-gene signature was constructed. Based on this risk score system, patients in the high-risk group had significantly poorer survival results than those in the low-risk group. Moreover, with the assistance of GEPIA http://gepia.cancer-pku.cn/, all five genes were found to be differentially expressed in GBM tissues compared with normal brain tissues. Furthermore, the co-expression network of the five genes was constructed based on weighted gene co-expression network analysis (WGCNA). Finally, this five-gene signature was further validated in other datasets. In conclusion, our study identified five novel biomarkers that have potential in the prognosis prediction of GBM.

Evaluation of Plasma Extracellular Vesicle MicroRNA Signatures for Lung Adenocarcinoma and Granuloma With Monte-Carlo Feature Selection Method

Extracellular Vesicle (EV) is a compilation of secreted vesicles, including micro vesicles, large oncosomes, and exosomes. It can be used in non-invasive diagnosis. MicroRNAs (miRNAs) processed by exosomes can be detected by liquid biopsy. To objectively evaluate the discriminative ability of miRNAs from whole plasma, EV and EV-free plasma, we analyzed the miRNA expression profiles in whole plasma, EV and EV-free plasma of 10 lung adenocarcinoma and 9 granuloma patients. With Monte-Carlo feature selection method, the top discriminative miRNAs in whole plasma, EV and EV-free plasma were identified, and they were quite different. Using the Repeated Incremental Pruning to Produce Error Reduction (RIPPER) method, we learned the classification rules: in whole plasma, granuloma patients did not express hsa-miR-223-3p while the lung adenocarcinoma patients expressed hsa-miR-223-3p; in EV, the hsa-miR-23b-3p was highly expressed in granuloma patients but not lung adenocarcinoma patients; in EV-free plasma, hsa-miR-376a-3p was expressed in granuloma patients but barely expressed in lung adenocarcinoma patients. For prediction performance, whole plasma had the highest weighted accuracy and EV outperformed EV-free plasma. Our results suggested that EV can be used as lung cancer biomarker. However, since it is less stable and not easy to detect, there are still technological difficulties to overcome.