Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

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.

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.

Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

A mini-TGA protein modulates gene expression through heterogeneous association with transcription factors.

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.

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 over-promising, 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.

Protoplast fusion in Bacillus species produces frequent, unbiased, genome-wide homologous recombination.

In eukaryotes, fine-scale maps of meiotic recombination events have greatly advanced our understanding of the factors that affect genomic variation patterns and evolution of traits. However, in bacteria that lack natural systems for sexual reproduction, unbiased characterization of recombination landscapes has remained challenging due to variable rates of genetic exchange and influence of natural selection. Here, to overcome these limitations and to gain a genome-wide view on recombination, we crossed Bacillus strains with different genetic distances using protoplast fusion. The offspring displayed complex inheritance patterns with one of the parents consistently contributing the major part of the chromosome backbone and multiple unselected fragments originating from the second parent. Our results demonstrate that this bias was in part due to the action of restriction-modification systems, whereas genome features like GC content and local nucleotide identity did not affect distribution of recombination events around the chromosome. Furthermore, we found that recombination occurred uniformly across the genome without concentration into hotspots. Notably, our results show that species-level genetic distance did not affect genome-wide recombination. This study provides a new insight into the dynamics of recombination in bacteria and a platform for studying recombination patterns in diverse bacterial species.

Initiation of fatty acid biosynthesis in Pseudomonas putida KT2440.

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.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci.

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.

Potential Pathogenicity Determinants Identified from Structural Proteomics of SARS-CoV and SARS-CoV-2.

Despite SARS-CoV and SARS-CoV-2 being equipped with highly similar protein arsenals, the corresponding zoonoses have spread among humans at extremely different rates. The specific characteristics of these viruses that led to such distinct outcomes remain unclear. Here, we apply proteome-wide comparative structural analysis aiming to identify the unique molecular elements in the SARS-CoV-2 proteome that may explain the differing consequences. By combining protein modeling and molecular dynamics simulations, we suggest nonconservative substitutions in functional regions of the spike glycoprotein (S), nsp1, and nsp3 that are contributing to differences in virulence. Particularly, we explain why the substitutions at the receptor-binding domain of S affect the structure-dynamics behavior in complexes with putative host receptors. Conservation of functional protein regions within the two taxa is also noteworthy. We suggest that the highly conserved main protease, nsp5, of SARS-CoV and SARS-CoV-2 is part of their mechanism of circumventing the host interferon antiviral response. Overall, most substitutions occur on the protein surfaces and may be modulating their antigenic properties and interactions with other macromolecules. Our results imply that the striking difference in the pervasiveness of SARS-CoV-2 and SARS-CoV among humans seems to significantly derive from molecular features that modulate the efficiency of viral particles in entering the host cells and blocking the host immune response.

Atomic-level characterization of protein-protein association.

Despite the biological importance of protein-protein complexes, determining their structures and association mechanisms remains an outstanding challenge. Here, we report the results of atomic-level simulations in which we observed five protein-protein pairs repeatedly associate to, and dissociate from, their experimentally determined native complexes using a molecular dynamics (MD)-based sampling approach that does not make use of any prior structural information about the complexes. To study association mechanisms, we performed additional, conventional MD simulations, in which we observed numerous spontaneous association events. A shared feature of native association for these five structurally and functionally diverse protein systems was that if the proteins made contact far from the native interface, the native state was reached by dissociation and eventual reassociation near the native interface, rather than by extensive interfacial exploration while the proteins remained in contact. At the transition state (the conformational ensemble from which association to the native complex and dissociation are equally likely), the protein-protein interfaces were still highly hydrated, and no more than 20% of native contacts had formed.

Does proximal versus distal injury location of the medial ulnar collateral ligament of the elbow differentially impact elbow stability? An ultrasound-guided and robot-assisted biomechanical study.

BACKGROUND: The location (proximal vs. distal) of elbow medial ulnar collateral ligament (MUCL) tears impacts clinical outcomes of nonoperative treatment. The purposes of our study were to (1) determine whether selective releases of the MUCL could be performed under ultrasound (US) guidance without disrupting overlying soft tissues, (2) assess the difference in medial elbow stability for proximal and distal releases of the MUCL using stress US and a robotic testing device, and (3) elucidate the flexion angle that resulted in the greatest amount of medial elbow laxity after MUCL injury. METHODS: Sixteen paired, fresh-frozen elbow specimens were used. Valgus laxity was evaluated with both US and robotic-assisted measurements before and after selective MUCL releases. A percutaneous US-guided technique was used to perform proximal MUCL releases in 8 elbows and to perform distal MUCL releases in their matched pairs. The robot was used to determine the elbow flexion angle at which the maximum valgus displacement occurred for both proximally and distally released specimens. Open dissection was then performed to assess the accuracy of the percutaneous releases. RESULTS: Percutaneous US-guided releases were successfully performed in 15 of 16 specimens. The proximal release resulted in greater valgus angle displacement (11° ± 2°) than the distal release (8° ± 2°) between flexion angles of 30° and 70° (P < .0001 at 30°, P < .0001 at 40°, P = .001 at 50°, P = .005 at 60°, and P = .020 at 70°). Valgus displacement between release locations did not reach the level of statistical significance between 80° and 120° (P = .051 at 80°, P = .131 at 90°, P = .245 at 100°, P = .400 at 110°, and P = .532 at 120°). When we compared the values for the mean increase in US delta gap (stressed - supported state) from before to after MUCL release, the proximally released elbows had larger increases than the distally released elbows (5.0 mm proximal vs. 3.7 mm distal, P = .032). After MUCL release, maximum mean valgus displacement occurred at 49° of flexion. CONCLUSIONS: US-guided selective releases of the MUCL can be performed reliably without violating the overlying musculature. Valgus instability is not of greater magnitude for distal releases when compared with proximal releases. This findings suggests there must be alternative factors to explain the difference in clinical prognosis between distal and proximal tears. The observed flexion angle for maximum valgus laxity could have important implications for elbow positioning during US or fluoroscopic stress examination, as well as surgical repair or reconstruction of the MUCL.

Pharmacoepidemiology, Machine Learning, and COVID-19: An Intent-to-Treat Analysis of Hydroxychloroquine, With or Without Azithromycin, and COVID-19 Outcomes Among Hospitalized US Veterans.

Hydroxychloroquine (HCQ) was proposed as an early therapy for coronavirus disease 2019 (COVID-19) after in vitro studies indicated possible benefit. Previous in vivo observational studies have presented conflicting results, though recent randomized clinical trials have reported no benefit from HCQ among patients hospitalized with COVID-19. We examined the effects of HCQ alone and in combination with azithromycin in a hospitalized population of US veterans with COVID-19, using a propensity score-adjusted survival analysis with imputation of missing data. According to electronic health record data from the US Department of Veterans Affairs health care system, 64,055 US Veterans were tested for the virus that causes COVID-19 between March 1, 2020 and April 30, 2020. Of the 7,193 veterans who tested positive, 2,809 were hospitalized, and 657 individuals were prescribed HCQ within the first 48-hours of hospitalization for the treatment of COVID-19. There was no apparent benefit associated with HCQ receipt, alone or in combination with azithromycin, and there was an increased risk of intubation when HCQ was used in combination with azithromycin (hazard ratio = 1.55; 95% confidence interval: 1.07, 2.24). In conclusion, we assessed the effectiveness of HCQ with or without azithromycin in treatment of patients hospitalized with COVID-19, using a national sample of the US veteran population. Using rigorous study design and analytic methods to reduce confounding and bias, we found no evidence of a survival benefit from the administration of HCQ.

Characterization of aromatic acid/proton symporters in Pseudomonas putida KT2440 toward efficient microbial conversion of lignin-related aromatics.

Pseudomonas putida KT2440 (hereafter KT2440) is a well-studied platform bacterium for the production of industrially valuable chemicals from heterogeneous mixtures of aromatic compounds obtained from lignin depolymerization. KT2440 can grow on lignin-related monomers, such as ferulate (FA), 4-coumarate (4CA), vanillate (VA), 4-hydroxybenzoate (4HBA), and protocatechuate (PCA). Genes associated with their catabolism are known, but knowledge about the uptake systems remains limited. In this work, we studied the KT2440 transporters of lignin-related monomers and their substrate selectivity. Based on the inhibition by protonophores, we focused on five genes encoding aromatic acid/H+ symporter family transporters categorized into major facilitator superfamily that uses the proton motive force. The mutants of PP_1376 (pcaK) and PP_3349 (hcnK) exhibited significantly reduced growth on PCA/4HBA and FA/4CA, respectively, while no change was observed on VA for any of the five gene mutants. At pH 9.0, the conversion of these compounds by hcnK mutant (FA/4CA) and vanK mutant (VA) was dramatically reduced, revealing that these transporters are crucial for the uptake of the anionic substrates at high pH. Uptake assays using 14C-labeled substrates in Escherichia coli and biosensor-based assays confirmed that PcaK, HcnK, and VanK have ability to take up PCA, FA/4CA, and VA/PCA, respectively. Additionally, analyses of the predicted protein structures suggest that the size and hydropathic properties of the substrate-binding sites of these transporters determine their substrate preferences. Overall, this study reveals that at physiological pH, PcaK and HcnK have a major role in the uptake of PCA/4HBA and FA/4CA, respectively, and VanK is a VA/PCA transporter. This information can contribute to the engineering of strains for the efficient conversion of lignin-related monomers to value-added chemicals.

Isoflurane anesthesia disrupts the cortical metabolome.

Identifying similarities and differences in the brain metabolome during different states of consciousness has broad relevance for neuroscience and state-dependent autonomic function. This study focused on the prefrontal cortex (PFC) as a brain region known to modulate states of consciousness. Anesthesia was used as a tool to eliminate wakefulness. Untargeted metabolomic analyses were performed on microdialysis samples obtained from mouse PFC during wakefulness and during isoflurane anesthesia. Analyses detected 2,153 molecules, 91 of which could be identified. Analytes were grouped as detected during both wakefulness and anesthesia (n = 61) and as unique to wakefulness (n = 23) or anesthesia (n = 7). Data were analyzed using univariate and multivariate approaches. Relative to wakefulness, during anesthesia there was a significant (q < 0.0001) fourfold change in 21 metabolites. During anesthesia 11 of these 21 molecules decreased and 10 increased. The Kyoto Encyclopedia of Genes and Genomes database was used to relate behavioral state-specific changes in the metabolome to metabolic pathways. Relative to wakefulness, most of the amino acids and analogs measured were significantly decreased during isoflurane anesthesia. Nucleosides and analogs were significantly increased during anesthesia. Molecules associated with carbohydrate metabolism, maintenance of lipid membranes, and normal cell functions were significantly decreased during anesthesia. Significant state-specific changes were also discovered among molecules comprising lipids and fatty acids, monosaccharides, and organic acids. Considered together, these molecules regulate point-to-point transmission, volume conduction, and cellular metabolism. The results identify a novel ensemble of candidate molecules in PFC as putative modulators of wakefulness and the loss of wakefulness.NEW & NOTEWORTHY The loss of wakefulness caused by a single concentration of isoflurane significantly altered levels of interrelated metabolites in the prefrontal cortex. The results support the interpretation that states of consciousness reflect dynamic interactions among cortical neuronal networks involving a humbling number of molecules that comprise the brain metabolome.

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.

Linking crop traits to transcriptome differences in a progeny population of tetraploid potato.

BACKGROUND: Potato is the third most consumed crop in the world. Breeding for traits such as yield, product quality and pathogen resistance are main priorities. Identifying molecular signatures of these and other important traits is important in future breeding efforts. In this study, a progeny population from a cross between a breeding line, SW93-1015, and a cultivar, Désirée, was studied by trait analysis and RNA-seq in order to develop understanding of segregating traits at the molecular level and identify transcripts with expressional correlation to these traits. Transcript markers with predictive value for field performance applicable under controlled environments would be of great value for plant breeding. RESULTS: A total of 34 progeny lines from SW93-1015 and Désirée were phenotyped for 17 different traits in a field in Nordic climate conditions and controlled climate settings. A master transcriptome was constructed with all 34 progeny lines and the parents through a de novo assembly of RNA-seq reads. Gene expression data obtained in a controlled environment from the 34 lines was correlated to traits by different similarity indices, including Pearson and Spearman, as well as DUO, which calculates the co-occurrence between high and low values for gene expression and trait. Our study linked transcripts to traits such as yield, growth rate, high laying tubers, late and tuber blight, tuber greening and early flowering. We found several transcripts associated to late blight resistance and transcripts encoding receptors were associated to Dickeya solani susceptibility. Transcript levels of a UBX-domain protein was negatively associated to yield and a GLABRA2 expression modulator was negatively associated to growth rate. CONCLUSION: In our study, we identify 100's of transcripts, putatively linked based on expression with 17 traits of potato, representing both well-known and novel associations. This approach can be used to link the transcriptome to traits. We explore the possibility of associating the level of transcript expression from controlled, optimal environments to traits in a progeny population with different methods introducing the application of DUO for the first time on transcriptome data. We verify the expression pattern for five of the putative transcript markers in another progeny population.

Plasticity of Phymatotrichopsis omnivora infection strategies is dependent on host and nonhost plant responses.

Necrotrophic fungi constitute the largest group of plant fungal pathogens that cause heavy crop losses worldwide. Phymatotrichopsis omnivora is a broad host, soil-borne necrotrophic fungal pathogen that infects over 2,000 dicotyledonous plants. The molecular basis of such broad host range is unknown. We conducted cell biology and transcriptomic studies in Medicago truncatula (susceptible), Brachypodium distachyon (resistant/nonhost), and Arabidopsis thaliana (partially resistant) to understand P. omnivora virulence mechanisms. We performed defence gene analysis, gene enrichments, and correlational network studies during key infection stages. We identified that P. omnivora infects the susceptible plant as a traditional necrotroph. However, it infects the partially resistant plant as a hemi-biotroph triggering salicylic acid-mediated defence pathways in the plant. Further, the infection strategy in partially resistant plants is determined by the host responses during early infection stages. Mutant analyses in A. thaliana established the role of small peptides PEP1 and PEP2 in defence against P. omnivora. The resistant/nonhost B. distachyon triggered stress responses involving sugars and aromatic acids. Bdwat1 mutant analysis identified the role of cell walls in defence. This is the first report that describes the plasticity in infection strategies of P. omnivora providing insights into broad host range.

Evolutionary reinforcement learning of dynamical large deviations.

We show how to bound and calculate the likelihood of dynamical large deviations using evolutionary reinforcement learning. An agent, a stochastic model, propagates a continuous-time Monte Carlo trajectory and receives a reward conditioned upon the values of certain path-extensive quantities. Evolution produces progressively fitter agents, potentially allowing the calculation of a piece of a large-deviation rate function for a particular model and path-extensive quantity. For models with small state spaces, the evolutionary process acts directly on rates, and for models with large state spaces, the process acts on the weights of a neural network that parameterizes the model's rates. This approach shows how path-extensive physics problems can be considered within a framework widely used in machine learning.

Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits.

Genome-wide association studies (GWAS) have great promise for identifying the loci that contribute to adaptive variation, but the complex genetic architecture of many quantitative traits presents a substantial challenge. We measured 14 morphological and physiological traits and identified single nucleotide polymorphism (SNP)-phenotype associations in a Populus trichocarpa population distributed from California, USA to British Columbia, Canada. We used whole-genome resequencing data of 882 trees with more than 6.78 million SNPs, coupled with multitrait association to detect polymorphisms with potentially pleiotropic effects. Candidate genes were validated with functional data. Broad-sense heritability (H2 ) ranged from 0.30 to 0.56 for morphological traits and 0.08 to 0.36 for physiological traits. In total, 4 and 20 gene models were detected using the single-trait and multitrait association methods, respectively. Several of these associations were corroborated by additional lines of evidence, including co-expression networks, metabolite analyses, and direct confirmation of gene function through RNAi. Multitrait association identified many more significant associations than single-trait association, potentially revealing pleiotropic effects of individual genes. This approach can be particularly useful for challenging physiological traits such as water-use efficiency or complex traits such as leaf morphology, for which we were able to identify credible candidate genes by combining multitrait association with gene co-expression and co-methylation data.

Genus-Wide Assessment of Lignocellulose Utilization in the Extremely Thermophilic Genus Caldicellulosiruptor by Genomic, Pangenomic, and Metagenomic Analyses.

Metagenomic data from Obsidian Pool (Yellowstone National Park, USA) and 13 genome sequences were used to reassess genus-wide biodiversity for the extremely thermophilic Caldicellulosiruptor The updated core genome contains 1,401 ortholog groups (average genome size for 13 species = 2,516 genes). The pangenome, which remains open with a revised total of 3,493 ortholog groups, encodes a variety of multidomain glycoside hydrolases (GHs). These include three cellulases with GH48 domains that are colocated in the glucan degradation locus (GDL) and are specific determinants for microcrystalline cellulose utilization. Three recently sequenced species, Caldicellulosiruptor sp. strain Rt8.B8 (renamed here Caldicellulosiruptor morganii), Thermoanaerobacter cellulolyticus strain NA10 (renamed here Caldicellulosiruptor naganoensis), and Caldicellulosiruptor sp. strain Wai35.B1 (renamed here Caldicellulosiruptor danielii), degraded Avicel and lignocellulose (switchgrass). C. morganii was more efficient than Caldicellulosiruptor bescii in this regard and differed from the other 12 species examined, both based on genome content and organization and in the specific domain features of conserved GHs. Metagenomic analysis of lignocellulose-enriched samples from Obsidian Pool revealed limited new information on genus biodiversity. Enrichments yielded genomic signatures closely related to that of Caldicellulosiruptor obsidiansis, but there was also evidence for other thermophilic fermentative anaerobes (Caldanaerobacter, Fervidobacterium, Caloramator, and Clostridium). One enrichment, containing 89.8% Caldicellulosiruptor and 9.7% Caloramator, had a capacity for switchgrass solubilization comparable to that of C. bescii These results refine the known biodiversity of Caldicellulosiruptor and indicate that microcrystalline cellulose degradation at temperatures above 70°C, based on current information, is limited to certain members of this genus that produce GH48 domain-containing enzymes.IMPORTANCE The genus Caldicellulosiruptor contains the most thermophilic bacteria capable of lignocellulose deconstruction, which are promising candidates for consolidated bioprocessing for the production of biofuels and bio-based chemicals. The focus here is on the extant capability of this genus for plant biomass degradation and the extent to which this can be inferred from the core and pangenomes, based on analysis of 13 species and metagenomic sequence information from environmental samples. Key to microcrystalline hydrolysis is the content of the glucan degradation locus (GDL), a set of genes encoding glycoside hydrolases (GHs), several of which have GH48 and family 3 carbohydrate binding module domains, that function as primary cellulases. Resolving the relationship between the GDL and lignocellulose degradation will inform efforts to identify more prolific members of the genus and to develop metabolic engineering strategies to improve this characteristic.