ABSTRACT
While the proliferation of data-driven omics technologies has continued to accelerate, methods of identifying relationships among large-scale changes from omics experiments have stagnated. It is therefore imperative to develop methods that can identify key mechanisms among one or more omics experiments in order to advance biological discovery. To solve this problem, here we describe the network-based algorithm MENTOR - Multiplex Embedding of Networks for Team-Based Omics Research. We demonstrate MENTOR's utility as a supervised learning approach to successfully partition a gene set containing multiple ontological functions into their respective functions. Subsequently, we used MENTOR as an unsupervised learning approach to identify important biological functions pertaining to the host genetic architectures in Populus trichocarpa associated with microbial abundance of multiple taxa. Moreover, as open source software designed with scientific teams in mind, we demonstrate the ability to use the output of MENTOR to facilitate distributed interpretation of omics experiments.
ABSTRACT
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.
Subject(s)
Agriculture , Microbiota , Plants , Microbiota/physiology , Plants/microbiology , Crops, Agricultural/microbiologyABSTRACT
Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects the quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids had a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well-characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 126 independent genetic loci, 69 of which are new. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level and cognitive traits. Integration of the genome-wide association studies findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, ß-blockers and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.
Subject(s)
Chronic Pain , Veterans , Adult , Humans , Chronic Pain/drug therapy , Chronic Pain/genetics , Genome-Wide Association Study/methods , Pain Measurement , Quality of Life , Genetic Predisposition to Disease , Polymorphism, Single Nucleotide/geneticsABSTRACT
Herbicide-resistant weeds are increasingly a problem in crop fields when exposed to similar chemistry over time. To avoid future yield losses, identifying herbicidal chemistry needs to be accelerated. We screened 50,000 small molecules using a liquid-handling robot and light microscopy focusing on pre-emergent herbicides in the family of cellulose biosynthesis inhibitors. Through phenotypic, chemical, genetic, and in silico methods we uncovered 6-{[4-(2-fluorophenyl)-1-piperazinyl]methyl}-N-(2-methoxy-5-methylphenyl)-1,3,5-triazine-2,4-diamine (fluopipamine). Symptomologies support fluopipamine as a putative antagonist of cellulose synthase enzyme 1 (CESA1) from Arabidopsis (Arabidopsis thaliana). Ectopic lignification, inhibition of etiolation, phenotypes including loss of anisotropic cellular expansion, swollen roots, and live cell imaging link fluopipamine to cellulose biosynthesis inhibition. Radiolabeled glucose incorporation of cellulose decreased in short-duration experiments when seedlings were incubated in fluopipamine. To elucidate the mechanism, ethylmethanesulfonate mutagenized M2 seedlings were screened for fluopipamine resistance. Two loci of genetic resistance were linked to CESA1. In silico docking of fluopipamine, quinoxyphen, and flupoxam against various CESA1 mutations suggests that an alternative binding site at the interface between CESA proteins is necessary to preserve cellulose polymerization in compound presence. These data uncovered potential fundamental mechanisms of cellulose biosynthesis in plants along with feasible leads for herbicidal uses.
Subject(s)
Arabidopsis Proteins , Arabidopsis , Herbicides , Arabidopsis/metabolism , Arabidopsis Proteins/metabolism , Cellulose/chemistry , Cell Wall/metabolism , Glucosyltransferases/chemistry , Seedlings/metabolism , Herbicides/pharmacology , Herbicides/metabolismABSTRACT
Fungal specialized metabolites are a major source of beneficial compounds that are routinely isolated, characterized, and manufactured as pharmaceuticals, agrochemical agents, and industrial chemicals. The production of these metabolites is encoded by biosynthetic gene clusters that are often silent under standard growth conditions. There are limited resources for characterizing the direct link between abiotic stimuli and metabolite production. Herein, we introduce a network analysis-based, data-driven algorithm comprising two routes to characterize the production of specialized fungal metabolites triggered by different exogenous compounds: the direct route and the auxiliary route. Both routes elucidate the influence of treatments on the production of specialized metabolites from experimental data. The direct route determines known and putative metabolites induced by treatments and provides additional insight over traditional comparison methods. The auxiliary route is specific for discovering unknown analytes, and further identification can be curated through online bioinformatic resources. We validated our algorithm by applying chitooligosaccharides and lipids at two different temperatures to the fungal pathogen Aspergillus fumigatus. After liquid chromatography-mass spectrometry quantification of significantly produced analytes, we used network centrality measures to rank the treatments' ability to elucidate these analytes and confirmed their identity through fragmentation patterns or in silico spiking with commercially available standards. Later, we examined the transcriptional regulation of these metabolites through real-time quantitative polymerase chain reaction. Our data-driven techniques can complement existing metabolomic network analysis by providing an approach to track the influence of any exogenous stimuli on metabolite production. Our experimental-based algorithm can overcome the bottlenecks in elucidating novel fungal compounds used in drug discovery.
ABSTRACT
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.
Subject(s)
CRISPR-Cas Systems , RNA, Guide, CRISPR-Cas Systems , Artificial Intelligence , DNA , Escherichia coli/genetics , Gene Editing , HumansABSTRACT
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.
ABSTRACT
Metabolite genome-wide association studies (mGWASs) are increasingly used to discover the genetic basis of target phenotypes in plants such as Populus trichocarpa, a biofuel feedstock and model woody plant species. Despite their growing importance in plant genetics and metabolomics, few mGWASs are experimentally validated. Here, we present a functional genomics workflow for validating mGWAS-predicted enzyme-substrate relationships. We focus on uridine diphosphate-glycosyltransferases (UGTs), a large family of enzymes that catalyze sugar transfer to a variety of plant secondary metabolites involved in defense, signaling, and lignification. Glycosylation influences physiological roles, localization within cells and tissues, and metabolic fates of these metabolites. UGTs have substantially expanded in P. trichocarpa, presenting a challenge for large-scale characterization. Using a high-throughput assay, we produced substrate acceptance profiles for 40 previously uncharacterized candidate enzymes. Assays confirmed 10 of 13 leaf mGWAS associations, and a focused metabolite screen demonstrated varying levels of substrate specificity among UGTs. A substrate binding model case study of UGT-23 rationalized observed enzyme activities and mGWAS associations, including glycosylation of trichocarpinene to produce trichocarpin, a major higher-order salicylate in P. trichocarpa. We identified UGTs putatively involved in lignan, flavonoid, salicylate, and phytohormone metabolism, with potential implications for cell wall biosynthesis, nitrogen uptake, and biotic and abiotic stress response that determine sustainable biomass crop production. Our results provide new support for in silico analyses and evidence-based guidance for in vivo functional characterization.
ABSTRACT
The intracellular cholesterol transporter NPC1 functions in late endosomes and lysosomes to efflux unesterified cholesterol, and its deficiency causes Niemann-Pick disease Type C, an autosomal recessive lysosomal disorder characterized by progressive neurodegeneration and early death. Here, we use single-nucleus RNA-seq on the forebrain of Npc1-/- mice at P16 to identify cell types and pathways affected early in pathogenesis. Our analysis uncovers significant transcriptional changes in the oligodendrocyte lineage during developmental myelination, accompanied by diminished maturation of myelinating oligodendrocytes. We identify upregulation of genes associated with neurogenesis and synapse formation in Npc1-/- oligodendrocyte lineage cells, reflecting diminished gene silencing by H3K27me3. Npc1-/- oligodendrocyte progenitor cells reproduce impaired maturation in vitro, and this phenotype is rescued by treatment with GSK-J4, a small molecule inhibitor of H3K27 demethylases. Moreover, mobilizing stored cholesterol in Npc1-/- mice by a single administration of 2-hydroxypropyl-ß-cyclodextrin at P7 rescues myelination, epigenetic marks, and oligodendrocyte gene expression. Our findings highlight an important role for NPC1 in oligodendrocyte lineage maturation and epigenetic regulation, and identify potential targets for therapeutic intervention.
Subject(s)
Niemann-Pick Disease, Type C , Animals , Mice , Cell Lineage , Cholesterol/metabolism , Epigenesis, Genetic , Membrane Transport Proteins/metabolism , Niemann-Pick Disease, Type C/genetics , Niemann-Pick Disease, Type C/metabolism , Oligodendroglia/metabolismABSTRACT
Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.
Subject(s)
Depressive Disorder, Major , Veterans , Humans , Suicidal Ideation , Veterans/psychology , Genome-Wide Association Study , Depressive Disorder, Major/genetics , Suicide, Attempted/psychology , Risk FactorsABSTRACT
Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids played a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 125 independent genetic loci, 82 of which are novel. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level, and cognitive traits. Integration of the GWAS findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, beta-blockers, and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.
ABSTRACT
For plants, distinguishing between mutualistic and pathogenic microbes is a matter of survival. All microbes contain microbe-associated molecular patterns (MAMPs) that are perceived by plant pattern recognition receptors (PRRs). Lysin motif receptor-like kinases (LysM-RLKs) are PRRs attuned for binding and triggering a response to specific MAMPs, including chitin oligomers (COs) in fungi, lipo-chitooligosaccharides (LCOs), which are produced by mycorrhizal fungi and nitrogen-fixing rhizobial bacteria, and peptidoglycan in bacteria. The identification and characterization of LysM-RLKs in candidate bioenergy crops including Populus are limited compared to other model plant species, thus inhibiting our ability to both understand and engineer microbe-mediated gains in plant productivity. As such, we performed a sequence analysis of LysM-RLKs in the Populus genome and predicted their function based on phylogenetic analysis with known LysM-RLKs. Then, using predictive models, molecular dynamics simulations, and comparative structural analysis with previously characterized CO and LCO plant receptors, we identified probable ligand-binding sites in Populus LysM-RLKs. Using several machine learning models, we predicted remarkably consistent binding affinity rankings of Populus proteins to CO. In addition, we used a modified Random Walk with Restart network-topology based approach to identify a subset of Populus LysM-RLKs that are functionally related and propose a corresponding signal transduction cascade. Our findings provide the first look into the role of LysM-RLKs in Populus-microbe interactions and establish a crucial jumping-off point for future research efforts to understand specificity and redundancy in microbial perception mechanisms.
ABSTRACT
Deciphering the mechanisms of bacterial fatty acid biosynthesis is crucial for both the engineering of bacterial hosts to produce fatty acid-derived molecules and the development of new antibiotics. However, gaps in our understanding of the initiation of fatty acid biosynthesis remain. Here, we demonstrate that the industrially relevant microbe Pseudomonas putida KT2440 contains three distinct pathways to initiate fatty acid biosynthesis. The first two routes employ conventional ß-ketoacyl-ACP synthase III enzymes, FabH1 and FabH2, that accept short- and medium-chain-length acyl-CoAs, respectively. The third route utilizes a malonyl-ACP decarboxylase enzyme, MadB. A combination of exhaustive in vivo alanine-scanning mutagenesis, in vitro biochemical characterization, X-ray crystallography, and computational modeling elucidate the presumptive mechanism of malonyl-ACP decarboxylation via MadB. Given that functional homologs of MadB are widespread throughout domain Bacteria, this ubiquitous alternative fatty acid initiation pathway provides new opportunities to target a range of biotechnology and biomedical applications.
Subject(s)
3-Oxoacyl-(Acyl-Carrier-Protein) Synthase , Pseudomonas putida , Pseudomonas putida/genetics , Pseudomonas putida/metabolism , 3-Oxoacyl-(Acyl-Carrier-Protein) Synthase/genetics , Mutagenesis , Fatty AcidsABSTRACT
TGA (TGACG-binding) transcription factors, which bind their target DNA through a conserved basic region leucine zipper (bZIP) domain, are vital regulators of gene expression in salicylic acid (SA)-mediated plant immunity. Here, we investigated the role of StTGA2.1, a potato (Solanum tuberosum) TGA lacking the full bZIP, which we named a mini-TGA. Such truncated proteins have been widely assigned as loss-of-function mutants. We, however, confirmed that StTGA2.1 overexpression compensates for SA-deficiency, indicating a distinct mechanism of action compared with model plant species. To understand the underlying mechanisms, we showed that StTGA2.1 can physically interact with StTGA2.2 and StTGA2.3, while its interaction with DNA was not detected. We investigated the changes in transcriptional regulation due to StTGA2.1 overexpression, identifying direct and indirect target genes. Using in planta transactivation assays, we confirmed that StTGA2.1 interacts with StTGA2.3 to activate StPRX07, a member of class III peroxidases (StPRX), which are known to play role in immune response. Finally, via structural modeling and molecular dynamics simulations, we hypothesized that the compact molecular architecture of StTGA2.1 distorts DNA conformation upon heterodimer binding to enable transcriptional activation. This study demonstrates how protein truncation can lead to distinct functions and that such events should be studied carefully in other protein families.
Subject(s)
Basic-Leucine Zipper Transcription Factors , Transcription Factors , Transcription Factors/genetics , Transcription Factors/metabolism , Transcriptional Activation , Gene Expression , Basic-Leucine Zipper Transcription Factors/genetics , Basic-Leucine Zipper Transcription Factors/metabolism , Gene Expression Regulation, PlantABSTRACT
Importance: Suicide is a leading cause of death; however, the molecular genetic basis of suicidal thoughts and behaviors (SITB) remains unknown. Objective: To identify novel, replicable genomic risk loci for SITB. Design, Setting, and Participants: This genome-wide association study included 633â¯778 US military veterans with and without SITB, as identified through electronic health records. GWAS was performed separately by ancestry, controlling for sex, age, and genetic substructure. Cross-ancestry risk loci were identified through meta-analysis. Study enrollment began in 2011 and is ongoing. Data were analyzed from November 2021 to August 2022. Main Outcome and Measures: SITB. Results: A total of 633â¯778 US military veterans were included in the analysis (57â¯152 [9%] female; 121â¯118 [19.1%] African ancestry, 8285 [1.3%] Asian ancestry, 452â¯767 [71.4%] European ancestry, and 51â¯608 [8.1%] Hispanic ancestry), including 121â¯211 individuals with SITB (19.1%). Meta-analysis identified more than 200 GWS (P < 5 × 10-8) cross-ancestry risk single-nucleotide variants for SITB concentrated in 7 regions on chromosomes 2, 6, 9, 11, 14, 16, and 18. Top single-nucleotide variants were largely intronic in nature; 5 were independently replicated in ISGC, including rs6557168 in ESR1, rs12808482 in DRD2, rs77641763 in EXD3, rs10671545 in DCC, and rs36006172 in TRAF3. Associations for FBXL19 and AC018880.2 were not replicated. Gene-based analyses implicated 24 additional GWS cross-ancestry risk genes, including FURIN, TSNARE1, and the NCAM1-TTC12-ANKK1-DRD2 gene cluster. Cross-ancestry enrichment analyses revealed significant enrichment for expression in brain and pituitary tissue, synapse and ubiquitination processes, amphetamine addiction, parathyroid hormone synthesis, axon guidance, and dopaminergic pathways. Seven other unique European ancestry-specific GWS loci were identified, 2 of which (POM121L2 and METTL15/LINC02758) were replicated. Two additional GWS ancestry-specific loci were identified within the African ancestry (PET112/GATB) and Hispanic ancestry (intergenic locus on chromosome 4) subsets, both of which were replicated. No GWS loci were identified within the Asian ancestry subset; however, significant enrichment was observed for axon guidance, cyclic adenosine monophosphate signaling, focal adhesion, glutamatergic synapse, and oxytocin signaling pathways across all ancestries. Within the European ancestry subset, genetic correlations (r > 0.75) were observed between the SITB phenotype and a suicide attempt-only phenotype, depression, and posttraumatic stress disorder. Additionally, polygenic risk score analyses revealed that the Million Veteran Program polygenic risk score had nominally significant main effects in 2 independent samples of veterans of European and African ancestry. Conclusions and Relevance: The findings of this analysis may advance understanding of the molecular genetic basis of SITB and provide evidence for ESR1, DRD2, TRAF3, and DCC as cross-ancestry candidate risk genes. More work is needed to replicate these findings and to determine if and how these genes might impact clinical care.
Subject(s)
Veterans , Humans , Female , Male , Suicidal Ideation , Genome-Wide Association Study , TNF Receptor-Associated Factor 3/genetics , Genetic Loci/genetics , Nucleotides , Polymorphism, Single Nucleotide/genetics , Genetic Predisposition to Disease/genetics , Proteins , Protein Serine-Threonine Kinases/geneticsABSTRACT
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.
ABSTRACT
Lipo-chitooligosaccharides (LCOs) are historically known for their role as microbial-derived signaling molecules that shape plant symbiosis with beneficial rhizobia or mycorrhizal fungi. Recent studies showing that LCOs are widespread across the fungal kingdom have raised questions about the ecological function of these compounds in organisms that do not form symbiotic relationships with plants. To elucidate the ecological function of these compounds, we investigate the metabolomic response of the ubiquitous human pathogen Aspergillus fumigatus to LCOs. Our metabolomics data revealed that exogenous application of various types of LCOs to A. fumigatus resulted in significant shifts in the fungal metabolic profile, with marked changes in the production of specialized metabolites known to mediate ecological interactions. Using network analyses, we identify specific types of LCOs with the most significant effect on the abundance of known metabolites. Extracts of several LCO-induced metabolic profiles significantly impact the growth rates of diverse bacterial species. These findings suggest that LCOs may play an important role in the competitive dynamics of non-plant-symbiotic fungi and bacteria. This study identifies specific metabolomic profiles induced by these ubiquitously produced chemicals and creates a foundation for future studies into the potential roles of LCOs as modulators of interkingdom competition. IMPORTANCE The activation of silent biosynthetic gene clusters (BGC) for the identification and characterization of novel fungal secondary metabolites is a perpetual motion in natural product discoveries. Here, we demonstrated that one of the best-studied symbiosis signaling compounds, lipo-chitooligosaccharides (LCOs), play a role in activating some of these BGCs, resulting in the production of known, putative, and unknown metabolites with biological activities. This collection of metabolites induced by LCOs differentially modulate bacterial growth, while the LCO standards do not convey the same effect. These findings create a paradigm shift showing that LCOs have a more prominent role outside of host recognition of symbiotic microbes. Importantly, our work demonstrates that fungi use LCOs to produce a variety of metabolites with biological activity, which can be a potential source of bio-stimulants, pesticides, or pharmaceuticals.
Subject(s)
Chitosan , Mycorrhizae , Humans , Chitin , Chitosan/pharmacology , Oligosaccharides/pharmacologyABSTRACT
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.
ABSTRACT
To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.