Unified Transcriptomic Signature of Arbuscular Mycorrhiza Colonization in Roots of Medicago truncatula by Integration of Machine Learning, Promoter Analysis, and Direct Merging Meta-Analysis

Abstract

Plant root symbiosis with Arbuscular mycorrhizal (AM) fungi improves uptake of water and mineral nutrients, improving plant development under stressful conditions. Unraveling the unified transcriptomic signature of a successful colonization provides a better understanding of symbiosis. We developed a framework for finding the transcriptomic signature of Arbuscular mycorrhiza colonization and its regulating transcription factors in roots of Medicago truncatula. Expression profiles of roots in response to AM species were collected from four separate studies and were combined by direct merging meta-analysis. Batch effect, the major concern in expression meta-analysis, was reduced by three normalization steps: Robust Multi-array Average algorithm, Z-standardization, and quartiling normalization. Then, expression profile of 33685 genes in 18 root samples of Medicago as numerical features, as well as study ID and Arbuscular mycorrhiza type as categorical features, were mined by seven models: RELIEF, UNCERTAINTY, GINI INDEX, Chi Squared, RULE, INFO GAIN, and INFO GAIN RATIO. In total, 73 genes selected by machine learning models were up-regulated in response to AM (Z-value difference > 0.5). Feature weighting models also documented that this signature is independent from study (batch) effect. The AM inoculation signature obtained was able to differentiate efficiently between AM inoculated and non-inoculated samples. The AP2 domain class transcription factor, GRAS family transcription factors, and cyclin-dependent kinase were among the highly expressed meta-genes identified in the signature. We found high correspondence between the AM colonization signature obtained in this study and independent RNA-seq experiments on AM colonization, validating the repeatability of the colonization signature. Promoter analysis of upregulated genes in the transcriptomic signature led to the key regulators of AM colonization, including the essential transcription factors for endosymbiosis establishment and development such as NF-YA factors. The approach developed in this study offers three distinct novel features: (I) it improves direct merging meta-analysis by integrating supervised machine learning models and normalization steps to reduce study-specific batch effects; (II) seven attribute weighting models assessed the suitability of each gene for the transcriptomic signature which contributes to robustness of the signature (III) the approach is justifiable, easy to apply, and useful in practice. Our integrative framework of meta-analysis, promoter analysis, and machine learning provides a foundation to reveal the transcriptomic signature and regulatory circuits governing Arbuscular mycorrhizal symbiosis and is transferable to the other biological settings.

Prediction of Alzheimer’s disease using blood gene expression data

“…First, the ComBat method, which was mainly used for internal and external validations, was used to remove batch effects among these datasets 27 , and then DEGs were selected using the "lmFit" function in the limma package. Second, the scaling & quartiling method by Mohammadi-Dehcheshmeh et al 71 was used to remove batch effects further, and then DEGs were selected using the "lmFit" function. Third, DEGs were computed for each dataset via the moderate t-test, and then rankings of p-values were used to curate meta-DEGs among three datasets 72 .…”

A Transcription Regulatory Sequence in the 5′ Untranslated Region of SARS-CoV-2 Is Vital for Virus Replication with an Altered Evolutionary Pattern against Human Inhibitory MicroRNAs

“…The availability of 212,346 (at 21 November 2020) SARS-CoV-2 genomic sequences, including 159,057 sequences of the full genome and high coverage, in GISAID ( ) [ 30 ] and NCBI provides the chance of pattern recognition in 5′UTR sequences of SARS-CoV-2, particularly against host microRNA inhibitory machinery, by machine learning models. Models and statistics such as decision tree classification based on association rule mining and deep learning [ 31 , 32 , 33 , 34 , 35 , 36 ] that have been used for eukaryotic promoter and UTR analysis, can be examined for UTR analysis of SARS-CoV-2. It should be noted that some of the sequences that have been deposited in GISAID as full genomes have incomplete or low-quality sequencing in 3′UTR and 5′UTR regions.…”

Machine learning compensates fold-change method and highlights oxidative phosphorylation in the brain transcriptome of Alzheimer’s disease

“…Although we did not combine datasets in this study, appropriate methods used for reducing the batch effect and differences between experiments 52 should be applied when combing datasets in future studies. We also noticed that random forest analysis dominated the identified gene features, indicating that future similar studies might focus on random forest first.…”