A High-Performance Computing Implementation of Iterative Random Forest for the Creation of Predictive Expression Networks

From Microbes to Microbiomes: Applications for Plant Health and Sustainable Agriculture

Plant-microbe interaction research has had a transformative trajectory, from individual microbial isolate studies to comprehensive analyses of plant microbiomes within the broader phytobiome framework. Acknowledging the indispensable role of plant microbiomes in shaping plant health, agriculture, and ecosystem resilience, we underscore the urgent need for sustainable crop production strategies in the face of contemporary challenges. We discuss how the synergies between advancements in 'omics technologies and artificial intelligence can help advance the profound potential of plant microbiomes. Furthermore, we propose a multifaceted approach encompassing translational considerations, transdisciplinary research initiatives, public-private partnerships, regulatory policy development, and pragmatic expectations for the practical application of plant microbiome knowledge across diverse agricultural landscapes. We advocate for strategic collaboration and intentional transdisciplinary efforts to unlock the benefits offered by plant microbiomes and address pressing global issues in food security. By emphasizing a nuanced understanding of plant microbiome complexities and fostering realistic expectations, we encourage the scientific community to navigate the transformative journey from discoveries in the laboratory to field applications. As companies specializing in agricultural microbes and microbiomes undergo shifts, we highlight the necessity of understanding how to approach sustainable agriculture with site-specific management solutions. While cautioning against overpromising, we underscore the excitement of exploring the many impacts of microbiome-plant interactions. We emphasize the importance of collaborative endeavors with societal partners to accelerate our collective capacity to harness the diverse and yet-to-be-discovered beneficial activities of plant microbiomes.

Learning epistatic polygenic phenotypes with Boolean interactions

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1, a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

Complex 33-beam simulated galactic cosmic radiation exposure impacts cognitive function and prefrontal cortex neurotransmitter networks in male mice

Astronauts will encounter extended exposure to galactic cosmic radiation (GCR) during deep space exploration, which could impair brain function. Here, we report that in male mice, acute or chronic GCR exposure did not modify reward sensitivity but did adversely affect attentional processes and increased reaction times. Potassium (K+)-stimulation in the prefrontal cortex (PFC) elevated dopamine (DA) but abolished temporal DA responsiveness after acute and chronic GCR exposure. Unlike acute GCR, chronic GCR increased levels of all other neurotransmitters, with differences evident between groups after higher K+-stimulation. Correlational and machine learning analysis showed that acute and chronic GCR exposure differentially reorganized the connection strength and causation of DA and other PFC neurotransmitter networks compared to controls which may explain space radiation-induced neurocognitive deficits.

Quantum biological insights into CRISPR-Cas9 sgRNA efficiency from explainable-AI driven feature engineering

CRISPR-Cas9 tools have transformed genetic manipulation capabilities in the laboratory. Empirical rules-of-thumb have been developed for only a narrow range of model organisms, and mechanistic underpinnings for sgRNA efficiency remain poorly understood. This work establishes a novel feature set and new public resource, produced with quantum chemical tensors, for interpreting and predicting sgRNA efficiency. Feature engineering for sgRNA efficiency is performed using an explainable-artificial intelligence model: iterative Random Forest (iRF). By encoding quantitative attributes of position-specific sequences for Escherichia coli sgRNAs, we identify important traits for sgRNA design in bacterial species. Additionally, we show that expanding positional encoding to quantum descriptors of base-pair, dimer, trimer, and tetramer sequences captures intricate interactions in local and neighboring nucleotides of the target DNA. These features highlight variation in CRISPR-Cas9 sgRNA dynamics between E. coli and H. sapiens genomes. These novel encodings of sgRNAs enhance our understanding of the elaborate quantum biological processes involved in CRISPR-Cas9 machinery.

Using iterative random forest to find geospatial environmental and Sociodemographic predictors of suicide attempts

Introduction

Despite a recent global decrease in suicide rates, death by suicide has increased in the United States. It is therefore imperative to identify the risk factors associated with suicide attempts to combat this growing epidemic. In this study, we aim to identify potential risk factors of suicide attempt using geospatial features in an Artificial intelligence framework.

Methods

We use iterative Random Forest, an explainable artificial intelligence method, to predict suicide attempts using data from the Million Veteran Program. This cohort incorporated 405,540 patients with 391,409 controls and 14,131 attempts. Our predictive model incorporates multiple climatic features at ZIP-code-level geospatial resolution. We additionally consider demographic features from the American Community Survey as well as the number of firearms and alcohol vendors per 10,000 people to assess the contributions of proximal environment, access to means, and restraint decrease to suicide attempts. In total 1,784 features were included in the predictive model.

Results

Our results show that geographic areas with higher concentrations of married males living with spouses are predictive of lower rates of suicide attempts, whereas geographic areas where males are more likely to live alone and to rent housing are predictive of higher rates of suicide attempts. We also identified climatic features that were associated with suicide attempt risk by age group. Additionally, we observed that firearms and alcohol vendors were associated with increased risk for suicide attempts irrespective of the age group examined, but that their effects were small in comparison to the top features.

Discussion

Taken together, our findings highlight the importance of social determinants and environmental factors in understanding suicide risk among veterans.

A primer on artificial intelligence in plant digital phenomics: embarking on the data to insights journey

Artificial intelligence (AI) has emerged as a fundamental component of global agricultural research that is poised to impact on many aspects of plant science. In digital phenomics, AI is capable of learning intricate structure and patterns in large datasets. We provide a perspective and primer on AI applications to phenome research. We propose a novel human-centric explainable AI (X-AI) system architecture consisting of data architecture, technology infrastructure, and AI architecture design. We clarify the difference between post hoc models and 'interpretable by design' models. We include guidance for effectively using an interpretable by design model in phenomic analysis. We also provide directions to sources of tools and resources for making data analytics increasingly accessible. This primer is accompanied by an interactive online tutorial.

Genetics of varicose veins reveals polygenic architecture and genetic overlap with arterial and venous disease

Varicose veins represent a common cause of cardiovascular morbidity, with limited available medical therapies. Although varicose veins are heritable and epidemiologic studies have identified several candidate varicose vein risk factors, the molecular and genetic basis remains uncertain. Here we analyzed the contribution of common genetic variants to varicose veins using data from the Veterans Affairs Million Veteran Program and four other large biobanks. Among 49,765 individuals with varicose veins and 1,334,301 disease-free controls, we identified 139 risk loci. We identified genetic overlap between varicose veins, other vascular diseases and dozens of anthropometric factors. Using Mendelian randomization, we prioritized therapeutic targets via integration of proteomic and transcriptomic data. Finally, topological enrichment analyses confirmed the biologic roles of endothelial shear flow disruption, inflammation, vascular remodeling and angiogenesis. These findings may facilitate future efforts to develop nonsurgical therapies for varicose veins.

Interpretable machine learning for genomics

High-throughput technologies such as next-generation sequencing allow biologists to observe cell function with unprecedented resolution, but the resulting datasets are too large and complicated for humans to understand without the aid of advanced statistical methods. Machine learning (ML) algorithms, which are designed to automatically find patterns in data, are well suited to this task. Yet these models are often so complex as to be opaque, leaving researchers with few clues about underlying mechanisms. Interpretable machine learning (iML) is a burgeoning subdiscipline of computational statistics devoted to making the predictions of ML models more intelligible to end users. This article is a gentle and critical introduction to iML, with an emphasis on genomic applications. I define relevant concepts, motivate leading methodologies, and provide a simple typology of existing approaches. I survey recent examples of iML in genomics, demonstrating how such techniques are increasingly integrated into research workflows. I argue that iML solutions are required to realize the promise of precision medicine. However, several open challenges remain. I examine the limitations of current state-of-the-art tools and propose a number of directions for future research. While the horizon for iML in genomics is wide and bright, continued progress requires close collaboration across disciplines.

Concrete Strength Prediction Using Different Machine Learning Processes: Effect of Slag, Fly Ash and Superplasticizer

Blast furnace slag (BFS) and fly ash (FA), as mining-associated solid wastes with good pozzolanic effects, can be combined with superplasticizer to prepare concrete with less cement utilization. Considering the important influence of strength on concrete design, random forest (RF) and particle swarm optimization (PSO) methods were combined to construct a prediction model and carry out hyper-parameter tuning in this study. Principal component analysis (PCA) was used to reduce the dimension of input features. The correlation coefficient (R), the explanatory variance score (EVS), the mean absolute error (MAE) and the mean square error (MSE) were used to evaluate the performance of the model. R = 0.954, EVS = 0.901, MAE = 3.746, and MSE = 27.535 of the optimal RF-PSO model on the testing set indicated the high generalization ability. After PCA dimensionality reduction, the R value decreased from 0.954 to 0.88, which was not necessary for the current dataset. Sensitivity analysis showed that cement was the most important feature, followed by water, superplasticizer, fine aggregate, BFS, coarse aggregate and FA, which was beneficial to the design of concrete schemes in practical projects. The method proposed in this study for estimation of the compressive strength of BFS-FA-superplasticizer concrete fills the research gap and has potential engineering application value.

Evaluating the Performance of Random Forest and Iterative Random Forest Based Methods when Applied to Gene Expression Data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Random Forest can be replaced in this process by iterative Random Forest (iRF), which performs variable selection and boosting. Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3 (RF-LOOP). We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

Provable Boolean interaction recovery from tree ensemble obtained via random forests

SignificanceRandom Forests (RFs) are among the most successful machine-learning algorithms in terms of prediction accuracy. In many domain problems, however, the primary goal is not prediction, but to understand the data-generation process-in particular, finding important features and feature interactions. There exists strong empirical evidence that RF-based methods-in particular, iterative RF (iRF)-are very successful in terms of detecting feature interactions. In this work, we propose a biologically motivated, Boolean interaction model. Using this model, we complement the existing empirical evidence with theoretical evidence for the ability of iRF-type methods to select desirable interactions. Our theoretical analysis also yields deeper insights into the general interaction selection mechanism of decision-tree algorithms and the importance of feature subsampling.

Murine Gut Microbiome Meta-analysis Reveals Alterations in Carbohydrate Metabolism in Response to Aging

Compositional and functional alterations to the gut microbiota during aging are hypothesized to potentially impact our health. Thus, determining aging-specific gut microbiome alterations is critical for developing microbiome-based strategies to improve health and promote longevity in the elderly. In this study, we performed a meta-analysis of publicly available 16S rRNA gene sequencing data from studies investigating the effect of aging on the gut microbiome in mice. Aging reproducibly increased gut microbial alpha diversity and shifted the microbial community structure in mice. We applied the bioinformatic tool PICRUSt2 to predict microbial metagenome function and established a random forest classifier to differentiate between microbial communities from young and old hosts and to identify aging-specific metabolic features. In independent validation data sets, this classifier achieved an area under the receiver operating characteristic curve (AUC) of 0.75 to 0.97 in differentiating microbiomes from young and old hosts. We found that 50% of the most important predicted aging-specific metabolic features were involved in carbohydrate metabolism. Furthermore, fecal short-chain fatty acid (SCFA) concentrations were significantly decreased in old mice, and the expression of the SCFA receptor Gpr41 in the colon was significantly correlated with the relative abundances of gut microbes and microbial carbohydrate metabolic pathways. In conclusion, this study identified aging-specific alterations in the composition and function of the gut microbiome and revealed a potential relationship between aging, microbial carbohydrate metabolism, fecal SCFA, and colonic Gpr41 expression. IMPORTANCE Aging-associated microbial alteration is hypothesized to play an important role in host health and longevity. However, investigations regarding specific gut microbes or microbial functional alterations associated with aging have had inconsistent results. We performed a meta-analysis across 5 independent studies to investigate the effect of aging on the gut microbiome in mice. Our analysis revealed that aging increased gut microbial alpha diversity and shifted the microbial community structure. To determine if we could reliably differentiate the gut microbiomes from young and old hosts, we established a random forest classifier based on predicted metagenome function and validated its performance against independent data sets. Alterations in microbial carbohydrate metabolism and decreased fecal short-chain fatty acid (SCFA) concentrations were key features of aging and correlated with host colonic expression of the SCFA receptor Gpr41. This study advances our understanding of the impact of aging on the gut microbiome and proposes a hypothesis that alterations in gut microbiota-derived SCFA-host GPR41 signaling are a feature of aging.

NF-κB perturbation reveals unique immunomodulatory functions in Prx1+ fibroblasts that promote development of atopic dermatitis

Skin is composed of diverse cell populations that cooperatively maintain homeostasis. Up-regulation of the nuclear factor κB (NF-κB) pathway may lead to the development of chronic inflammatory disorders of the skin, but its role during the early events remains unclear. Through analysis of single-cell RNA sequencing data via iterative random forest leave one out prediction, an explainable artificial intelligence method, we identified an immunoregulatory role for a unique paired related homeobox-1 (Prx1)⁺ fibroblast subpopulation. Disruption of Ikkb-NF-κB under homeostatic conditions in these fibroblasts paradoxically induced skin inflammation due to the overexpression of C-C motif chemokine ligand 11 (CCL11; or eotaxin-1) characterized by eosinophil infiltration and a subsequent T_H2 immune response. Because the inflammatory phenotype resembled that seen in human atopic dermatitis (AD), we examined human AD skin samples and found that human AD fibroblasts also overexpressed CCL11 and that perturbation of Ikkb-NF-κB in primary human dermal fibroblasts up-regulated CCL11. Monoclonal antibody treatment against CCL11 was effective in reducing the eosinophilia and T_H2 inflammation in a mouse model. Together, the murine model and human AD specimens point to dysregulated Prx1⁺ fibroblasts as a previously unrecognized etiologic factor that may contribute to the pathogenesis of AD and suggest that targeting CCL11 may be a way to treat AD-like skin lesions.

Antiviral Strategies Against SARS-CoV-2: A Systems Biology Approach

The unprecedented scientific achievements in combating the COVID-19 pandemic reflect a global response informed by unprecedented access to data. We now have the ability to rapidly generate a diversity of information on an emerging pathogen and, by using high-performance computing and a systems biology approach, we can mine this wealth of information to understand the complexities of viral pathogenesis and contagion like never before. These efforts will aid in the development of vaccines, antiviral medications, and inform policymakers and clinicians. Here we detail computational protocols developed as SARS-CoV-2 began to spread across the globe. They include pathogen detection, comparative structural proteomics, evolutionary adaptation analysis via network and artificial intelligence methodologies, and multiomic integration. These protocols constitute a core framework on which to build a systems-level infrastructure that can be quickly brought to bear on future pathogens before they evolve into pandemic proportions.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

BACKGROUND

A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic.

RESULTS

Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus.

CONCLUSIONS

These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

Neurotransmitter networks in mouse prefrontal cortex are reconfigured by isoflurane anesthesia

This study quantified eight small-molecule neurotransmitters collected simultaneously from prefrontal cortex of C57BL/6J mice (n = 23) during wakefulness and during isoflurane anesthesia (1.3%). Using isoflurane anesthesia as an independent variable enabled evaluation of the hypothesis that isoflurane anesthesia differentially alters concentrations of multiple neurotransmitters and their interactions. Machine learning was applied to reveal higher order interactions among neurotransmitters. Using a between-subjects design, microdialysis was performed during wakefulness and during anesthesia. Concentrations (nM) of acetylcholine, adenosine, dopamine, GABA, glutamate, histamine, norepinephrine, and serotonin in the dialysis samples are reported (means ± SD). Relative to wakefulness, acetylcholine concentration was lower during isoflurane anesthesia (1.254 ± 1.118 vs. 0.401 ± 0.134, P = 0.009), and concentrations of adenosine (29.456 ± 29.756 vs. 101.321 ± 38.603, P < 0.001), dopamine (0.0578 ± 0.0384 vs. 0.113 ± 0.084, P = 0.036), and norepinephrine (0.126 ± 0.080 vs. 0.219 ± 0.066, P = 0.010) were higher during anesthesia. Isoflurane reconfigured neurotransmitter interactions in prefrontal cortex, and the state of isoflurane anesthesia was reliably predicted by prefrontal cortex concentrations of adenosine, norepinephrine, and acetylcholine. A novel finding to emerge from machine learning analyses is that neurotransmitter concentration profiles in mouse prefrontal cortex undergo functional reconfiguration during isoflurane anesthesia. Adenosine, norepinephrine, and acetylcholine showed high feature importance, supporting the interpretation that interactions among these three transmitters may play a key role in modulating levels of cortical and behavioral arousal.NEW & NOTEWORTHY This study discovered that interactions between neurotransmitters in mouse prefrontal cortex were altered during isoflurane anesthesia relative to wakefulness. Machine learning further demonstrated that, relative to wakefulness, higher order interactions among neurotransmitters were disrupted during isoflurane administration. These findings extend to the neurochemical domain the concept that anesthetic-induced loss of wakefulness results from a disruption of neural network connectivity.