A Meta-Analysis Based Method for Prioritizing Candidate Genes Involved in a Pre-specific Function

Machine Learning-Driven discovery of immunogenic cell Death-Related biomarkers and molecular classification for diabetic ulcers

In this study, we redefine the diagnostic landscape of diabetic ulcers (DUs), a major diabetes complication. Our research uncovers new biomarkers linked to immunogenic cell death (ICD) in DUs by utilizing RNA-sequencing data of Gene Expression Omnibus (GEO) analysis combined with a comprehensive database interrogation. Employing a random forest algorithm, we have developed a diagnostic model that demonstrates improved accuracy in distinguishing DUs from normal tissue, with satisfactory results from ROC analysis. Beyond mere diagnosis, our model categorizes DUs into novel molecular classifications, which may enhance our comprehension of their underlying pathophysiology. This study bridges the gap between molecular insights and clinical practice. It sets the stage for transformative strategies in DUs management, marking a significant step forward in personalized medicine for diabetic patients.

Metabolism-related biomarkers, molecular classification, and immune infiltration in diabetic ulcers with validation

Diabetes mellitus (DM) can lead to diabetic ulcers (DUs), which are the most severe complications. Due to the need for more accurate patient classifications and diagnostic models, treatment and management strategies for DU patients still need improvement. The difficulty of diabetic wound healing is caused closely related to biological metabolism and immune chemotaxis reaction dysfunction. Therefore, the purpose of our study is to identify metabolic biomarkers in patients with DU and construct a molecular subtype-specific prognostic model that is highly accurate and robust. RNA-sequencing data for DU samples were obtained from the Gene Expression Omnibus (GEO) database. DU patients and normal individuals were compared regarding the expression of metabolism-related genes (MRGs). Then, a novel diagnostic model based on MRGs was constructed with the random forest algorithm, and classification performance was evaluated utilizing receiver operating characteristic (ROC) analysis. The biological functions of MRGs-based subtypes were investigated using consensus clustering analysis. A principal component analysis (PCA) was conducted to determine whether MRGs could distinguish between subtypes. We also examined the correlation between MRGs and immune infiltration. Lastly, qRT-PCR was utilized to validate the expression of the hub MRGs with clinical validations and animal experimentations. Firstly, 8 metabolism-related hub genes were obtained by random forest algorithm, which could distinguish the DUs from normal samples validated by the ROC curves. Secondly, DU samples could be consensus clustered into three molecular classifications by MRGs, verified by PCA analysis. Thirdly, associations between MRGs and immune infiltration were confirmed, with LYN and Type 1 helper cell significantly positively correlated; RHOH and TGF-β family remarkably negatively correlated. Finally, clinical validations and animal experiments of DU skin tissue samples showed that the expressions of metabolic hub genes in the DU groups were considerably upregulated, including GLDC, GALNT6, RHOH, XDH, MMP12, KLK6, LYN, and CFB. The current study proposed an auxiliary MRGs-based DUs model while proposing MRGs-based molecular clustering and confirmed the association with immune infiltration, facilitating the diagnosis and management of DU patients and designing individualized treatment plans.

An advanced systems biology framework of feature engineering for cold tolerance genes discovery from integrated omics and non-omics data in soybean

Soybean is sensitive to low temperatures during the crop growing season. An urgent demand for breeding cold-tolerant cultivars to alleviate the production loss is apparent to cope with this scenario. Cold-tolerant trait is a complex and quantitative trait controlled by multiple genes, environmental factors, and their interaction. In this study, we proposed an advanced systems biology framework of feature engineering for the discovery of cold tolerance genes (CTgenes) from integrated omics and non-omics (OnO) data in soybean. An integrative pipeline was introduced for feature selection and feature extraction from different layers in the integrated OnO data using data ensemble methods and the non-parameter random forest prioritization to minimize uncertainties and false positives for accuracy improvement of results. In total, 44, 143, and 45 CTgenes were identified in short-, mid-, and long-term cold treatment, respectively, from the corresponding gene-pool. These CTgenes outperformed the remaining genes, the random genes, and the other candidate genes identified by other approaches in an independent RNA-seq database. Furthermore, we applied pathway enrichment and crosstalk network analyses to uncover relevant physiological pathways with the discovery of underlying cold tolerance in hormone- and defense-related modules. Our CTgenes were validated by using 55 SNP genotype data of 56 soybean samples in cold tolerance experiments. This suggests that the CTgenes identified from our proposed systematic framework can effectively distinguish cold-resistant and cold-sensitive lines. It is an important advancement in the soybean cold-stress response. The proposed pipelines provide an alternative solution to biomarker discovery, module discovery, and sample classification underlying a particular trait in plants in a robust and efficient way.

Feature importance network reveals novel functional relationships between biological features in Arabidopsis thaliana

Understanding how the different cellular components are working together to form a living cell requires multidisciplinary approaches combining molecular and computational biology. Machine learning shows great potential in life sciences, as it can find novel relationships between biological features. Here, we constructed a dataset of 11,801 gene features for 31,522 Arabidopsis thaliana genes and developed a machine learning workflow to identify linked features. The detected linked features are visualised as a Feature Important Network (FIN), which can be mined to reveal a variety of novel biological insights pertaining to gene function. We demonstrate how FIN can be used to generate novel insights into gene function. To make this network easily accessible to the scientific community, we present the FINder database, available at finder.plant.tools.

Autophagy and Programmed Cell Death Are Critical Pathways in Jasmonic Acid Mediated Saline Stress Tolerance in Oryza sativa

Saline stress is the most limiting condition impacting the plant growth, development, and productivity. In this present study, jasmonic acid (JA) was used as a foliar spray on the rice seedlings grown under saline stress. Increase in photosynthetic pigments, anthocyanin, and total protein content was observed with JA treatment while NaCl showed reduction in biochemical constituents and enhanced antioxidant enzyme activity. The leaf cells of NaCl-treated seedlings accumulated more ROS and had more fragmented nuclei, whereas JA decreased the accumulation and fragmentation during saline stress. In NaCl treatment, gene expression analysis showed many fold upregulation in comparison with other treatments. The results suggest that JA acts as a promoter for growth, physiological, biochemical, and cellular contents, as well as ameliorate the effects of saline stress. The expression of genes demonstrated that saline stress may promote autophagy, which leads to autophagic cell death, and improve tolerance to saline stress in rice seedlings via the jasmonic acid signaling pathway. However, the mechanism by which jasmonate signaling induces autophagy and cell death is unknown and requires further exploration.

Status and prospects of genome-wide association studies in plants

Genome-wide association studies (GWAS) have developed into a powerful and ubiquitous tool for the investigation of complex traits. In large part, this was fueled by advances in genomic technology, enabling us to examine genome-wide genetic variants across diverse genetic materials. The development of the mixed model framework for GWAS dramatically reduced the number of false positives compared with naïve methods. Building on this foundation, many methods have since been developed to increase computational speed or improve statistical power in GWAS. These methods have allowed the detection of genomic variants associated with either traditional agronomic phenotypes or biochemical and molecular phenotypes. In turn, these associations enable applications in gene cloning and in accelerated crop breeding through marker assisted selection or genetic engineering. Current topics of investigation include rare-variant analysis, synthetic associations, optimizing the choice of GWAS model, and utilizing GWAS results to advance knowledge of biological processes. Ongoing research in these areas will facilitate further advances in GWAS methods and their applications.

Prioritization and Evaluation of Flooding Tolerance Genes in Soybean [Glycine max (L.) Merr.]

Soybean [Glycine max (L.) Merr.] is one of the most important legume crops abundant in edible protein and oil in the world. In recent years there has been increasingly more drastic weather caused by climate change, with flooding, drought, and unevenly distributed rainfall gradually increasing in terms of the frequency and intensity worldwide. Severe flooding has caused extensive losses to soybean production and there is an urgent need to breed strong soybean seeds with high flooding tolerance. The present study demonstrates bioinformatics big data mining and integration, meta-analysis, gene mapping, gene prioritization, and systems biology for identifying prioritized genes of flooding tolerance in soybean. A total of 83 flooding tolerance genes (FTgenes), according to the appropriate cut-off point, were prioritized from 36,705 test genes collected from multidimensional genomic features linking to soybean flooding tolerance. Several validation results using independent samples from SoyNet, genome-wide association study, SoyBase, GO database, and transcriptome databases all exhibited excellent agreement, suggesting these 83 FTgenes were significantly superior to others. These results provide valuable information and contribution to research on the varieties selection of soybean.

Transcriptome-Wide Annotation of m5C RNA Modifications Using Machine Learning

The emergence of epitranscriptome opened a new chapter in gene regulation. 5-methylcytosine (m⁵C), as an important post-transcriptional modification, has been identified to be involved in a variety of biological processes such as subcellular localization and translational fidelity. Though high-throughput experimental technologies have been developed and applied to profile m⁵C modifications under certain conditions, transcriptome-wide studies of m⁵C modifications are still hindered by the dynamic and reversible nature of m⁵C and the lack of computational prediction methods. In this study, we introduced PEA-m5C, a machine learning-based m⁵C predictor trained with features extracted from the flanking sequence of m⁵C modifications. PEA-m5C yielded an average AUC (area under the receiver operating characteristic) of 0.939 in 10-fold cross-validation experiments based on known Arabidopsis m⁵C modifications. A rigorous independent testing showed that PEA-m5C (Accuracy [Acc] = 0.835, Matthews correlation coefficient [MCC] = 0.688) is remarkably superior to the recently developed m⁵C predictor iRNAm5C-PseDNC (Acc = 0.665, MCC = 0.332). PEA-m5C has been applied to predict candidate m⁵C modifications in annotated Arabidopsis transcripts. Further analysis of these m⁵C candidates showed that 4nt downstream of the translational start site is the most frequently methylated position. PEA-m5C is freely available to academic users at: https://github.com/cma2015/PEA-m5C.