Engineering topics
Publications and source records attributed to Jacobson, Daniel.
An interactive map showing Populus trichocarpa GWAS population structure estimated by ADMIXTURE (k=3, selected as optimal from k=2-11). Sampling locations are colored by their predominant ancestry proportion among the three inferred populations and geographic origins are searchable by genotype or river system using the search bar.
Abiotic stresses like heat and drought pose a major threat to forest ecosystems and are anticipated to increase due to anthropogenic climate change. One mechanism by which plants cope with stress
This narrative presents a demonstration of RWRtools methods to the important bioenergy feedstock Panicum virgatum (Switchgrass) by leveraging networks from the model plant Arabidopsis thaliana. We used a genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with Switchgrass shoot biomass. We then assigned these SNPs to Switchgrass genes, and using Phytozome (Goodstein et al., 2012), assigned these 38 Switchgrass GWAS genes to 32 Arabidopsis orthologs. These 32 genes are the starting set of genes (i.e. the seed gene set) for this demo narrative.
Computational and data-enabled science and engineering are revolutionizing advances throughout science and society, at all scales of computing. For example, teams in the U.S. Department of Energy’s Exascale Computing Project have been tackling new frontiers in modeling, simulation, and analysis by exploiting unprecedented exascale computing capabilities—building an advanced software ecosystem that supports next-generation applications and addresses disruptive changes in computer architectures. However, concerns are growing about the productivity of the developers of scientific software. Members of the Interoperable Design of Extreme-scale Application Software project serve as catalysts to address these challenges through fostering software communities, incubating and curating methodologies and resources, and disseminating knowledge to advance developer productivity and software sustainability. This article discusses how these synergistic activities are advancing scientific discovery—mitigating technical risks by building a firmer foundation for reproducible, sustainable science at all scales of computing, from laptops to clusters to exascale and beyond.
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. Here, 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.
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.
Not Available
As the US Department of Energy (DOE) computing facilities began deploying petascale systems in 2008, DOE was already setting its sights on exascale. In that year, DARPA published a report on the feasibility of reaching exascale. The report authors identified several key challenges in the pursuit of exascale including power, memory, concurrency, and resiliency. That report informed the DOE's computing strategy for reaching exascale. With the deployment of Oak Ridge National Laboratory's Frontier supercomputer, we have officially entered the exascale era. In this paper, we discuss Frontier's architecture, how it addresses those challenges, and describe some early application results from Oak Ridge Leadership Computing Facility's Center of Excellence and the Exascale Computing Project.
Plant phenotyping is typically a time-consuming and expensive endeavor, requiring large groups of researchers to meticulously measure biologically relevant plant traits, and is the main bottleneck in understanding plant adaptation and the genetic architecture underlying complex traits at population scale. In this work, we address these challenges by leveraging few-shot learning with convolutional neural networks to segment the leaf body and visible venation of 2,906 Populus trichocarpa leaf images obtained in the field. In contrast to previous methods, our approach (a) does not require experimental or image preprocessing, (b) uses the raw RGB images at full resolution, and (c) requires very few samples for training (e.g., just 8 images for vein segmentation). Traits relating to leaf morphology and vein topology are extracted from the resulting segmentations using traditional open-source image-processing tools, validated using real-world physical measurements, and used to conduct a genome-wide association study to identify genes controlling the traits. In this way, the current work is designed to provide the plant phenotyping community with (a) methods for fast and accurate image-based feature extraction that require minimal training data and (b) a new population-scale dataset, including 68 different leaf phenotypes, for domain scientists and machine learning researchers. All of the few-shot learning code, data, and results are made publicly available.
Plant microbiomes are assembled and modified through a complex milieu of biotic and abiotic factors. Despite dynamic and fluctuating contributing variables, specific host metabolites are consistently identified as important mediators of microbial interactions. We combine information from a large-scale metatranscriptomic dataset from natural poplar trees and experimental genetic manipulation assays in seedlings of the model plant Arabidopsis thaliana to converge on a conserved role for transport of the plant metabolite myo-inositol in mediating host-microbe interactions. While microbial catabolism of this compound has been linked to increased host colonization, we identify bacterial phenotypes that occur in both catabolism-dependent and -independent manners, suggesting that myo-inositol may additionally serve as a eukaryotic-derived signaling molecule to modulate microbial activities. Our data suggest host control of this compound and resulting microbial behavior are important mechanisms at play surrounding the host metabolite myo-inositol.
In this work we identify changes in high-resolution zones across the globe linked by environmental similarity that have implications for agriculture, bioenergy, and zoonosis. We refine exhaustive vector comparison methods with improved similarity metrics as well as provide multiple methods of amalgamation across 744 months of climatic data. The results of the vector comparison are captured as networks which are analyzed using static and longitudinal comparison methods to reveal locations around the globe experiencing dramatic changes in abiotic stress. Specifically we (i) incorporate updated similarity scores and provide a comparison between similarity metrics, (ii) implement a new feature for resource optimization, (iii) compare an agglomerative view to a longitudinal view, (iv) compare across 2-way and 3-way vector comparisons, (v) implement a new form of analysis, and (vi) demonstrate biological applications and discuss implications across a diverse set of species distributions by detecting changes that affect their habitats. Species of interest are related to agriculture (e.g., coffee, wine, chocolate), bioenergy (e.g., poplar, switchgrass, pennycress), as well as those living in zones of concern for zoonotic spillover that may lead to pandemics (e.g., eucalyptus, flying foxes).
Current methods for defining SARS-CoV-2 lineages ignore the vast majority of the SARS-CoV-2 genome. We develop and apply an exhaustive vector comparison method that directly compares all known SARS-CoV-2 genome sequences to produce novel lineage classifications. We utilize data-driven models that (i) accurately capture the complex interactions across the set of all known SARS-CoV-2 genomes, (ii) scale to leadership-class computing systems, and (iii) enable tracking how such strains evolve geospatially over time. We show that during the height of the original Omicron surge, countries across Europe, Asia, and the Americas had a spatially asynchronous distribution of Omicron sub-strains. Moreover, neighboring countries were often dominated by either different clusters of the same variant or different variants altogether throughout the pandemic. Analyses of this kind may suggest a different pattern of epidemiological risk than was understood from conventional data, as well as produce actionable insights and transform our ability to prepare for and respond to current and future biological threats.
Substance use disorder is a debilitating clinical condition in which behavioral dependence results from biological, environmental, genetic, and psychosocial factors. An epidemic surrounding the use and abuse of opioids is ravaging the world. While considerable efforts have explored the social drivers of addiction, a deeper understanding of biological causes and genetic vulnerabilities, preventative interventions, and effective treatments, have all proven elusive. This perspective article aims to remind readers that addictive natural compounds such as cocaine, nicotine, cathinone, or morphine, evolved as defensive metabolites to deter insect herbivory. The molecular mechanisms underlying motivational seeking and learning/reward show remarkable conservation since their early emergence in bilateral metazoans. An extended coevolutionary arms race subsequently weaponized these compounds into disruptors of learning, motivation, and incentivized attention. When plant chemical defenses attack insect physiology, humans are rendered susceptible due to strong conservation in the underlying molecular machinery. This perspective addresses the paradox that opiates were shaped to target insect neuropharmacology, even though this taxon appears to lack the recognized opioid receptor clade of mammals. We argue that the link is to be found in the allatostatin receptor, a basal ortholog of opioid receptors. Moreover, preliminary evidence indicates that morphine reduces Drosophila feeding and locomotion, concordant with a purported role as a defensive compound reducing herbivory. An implementation via allatostatin-mediated mechanisms is likely. This research argues for a broader heuristic perspective of substance abuse and a recognition of the evolutionary constraints that have likely shaped the biological drivers of opioid sensitivity and of its behavioral targets.
Identification of proteins dysregulated by COVID-19 infection is critically important for better understanding of its pathophysiology, building prognostic models, and identifying new targets. Plasma proteomic profiling of 4,301 proteins was performed in two independent datasets and tested for the association for three COVID-19 outcomes (infection, ventilation, and death). We identified 1,449 proteins consistently associated in both datasets with any of these three outcomes. We subsequently created highly accurate models that distinctively predict infection, ventilation, and death. These proteins were enriched in specific biological processes including cytokine signaling, Alzheimer’s disease, and coronary artery disease. Mendelian randomization and gene network analyses identified eight causal proteins and 141 highly connected hub proteins including 35 with known drug targets. Our findings provide distinctive prognostic biomarkers for two severe COVID-19 outcomes, reveal their relationship to Alzheimer’s disease and coronary artery disease, and identify potential therapeutic targets for COVID-19 outcomes.
Oat-based value-added products have increased their value as healthy foodstuff. Fusarium head blight (FHB) infections and the mycotoxins accumulated to the oat seeds, however, pose a challenge to oat production. The FHB infections are predicted to become more prevalent in the future changing climates and under more limited use of fungicides. Both these factors increase the pressure for breeding new resistant cultivars. Until now, however, genetic links in oats against FHB infection have been difficult to identify. Therefore, there is a great need for more effective breeding efforts, including improved phenotyping methods allowing time series analysis and the identification of molecular markers during disease progression. To these ends, dissected spikelets of several oat genotypes with different resistance profiles were studied by image-based methods during disease progression by Fusarium culmorum or F. langsethiae species. The chlorophyll fluorescence of each pixel in the spikelets was recorded after inoculation by the two Fusarium spp., and the progression of the infections was analyzed by calculating the mean maximum quantum yield of PSII (F v /F m ) values for each spikelet. The recorded values were (i) the change in the photosynthetically active area of the spikelet as percentage of its initial size, and (ii) the mean of F v /F m values of all fluorescent pixels per spikelet post inoculation, both indicative of the progression of the FHB disease. The disease progression was successfully monitored, and different stages of the infection could be defined along the time series. The data also confirmed the differential rate of disease progression by the two FHB causal agents. In addition, oat varieties with variable responses to the infections were indicated.
In this work, we use few-shot learning to segment the body and vein architecture of P. trichocarpa leaves from high-resolution scans obtained in the UC Davis common garden. Leaf and vein segmentation are formulated as separate tasks, in which convolutional neural networks (CNNs) are used to iteratively expand partial segmentations until reaching stopping criteria. Our leaf and vein segmentation approaches use just 50 and 8 manually traced images for training, respectively, and are applied to a set of 2,634 top and bottom leaf scans. We show that both methods achieve high segmentation accuracy, in some cases exceeding even human-level segmentation. The leaf and vein segmentations are subsequently used to extract 68 morphological traits using traditional open-source image processing tools, which are validated using real-world physical measurements. For a biological perspective, we perform a genome-wide association study using the vein density trait to discover novel genetic architectures associated with multiple physiological processes relating to leaf development and function. In addition to sharing all of the few-shot learning code (see https://github.com/jlager/few-shot-leaf-segmentation), we are releasing all images, manual segmentations, model predictions, 68 extracted leaf phenotypes, and a new set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes. The data folder includes all images, ground truth segmentations, predicted segmentations, and extracted leaf traits. All images encode the sample ID in the file name by indicating the treatment, block, row, position, and leaf side, respectively. For example, the file, C_1_1_2_bot.jpeg, indicates the control treatment, block 1, row 1, position 2, and the bottom side of the leaf. Tabulated results include position IDs as well as the corresponding genotype IDs. The images folder includes the 2,906 high-resolution leaf scans taken in the field. The leaf_masks folder includes 50 ground truth segmentations used for training the leaf tracing algorithm. The leaf_preds folder includes the 2,906 predicted segmentations from the leaf tracing algorithm. The vein_masks folder includes 8 ground truth segmentations used for training the vein growing algorithm. The vein_preds folder includes the 1,453 predicted segmentations from the vein growing algorithm. The vein_probs folder includes the 1,453 predicted probability maps from the vein growing algorithm before thresholding. The genomes folder includes the set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes with a README file detailing the steps taken. The results folder includes: (i) raw values of the 68 predicted leaf traits in digital_traits.tsv, (ii) manually measured values of petiole length and width in manual_traits.tsv, (iii) thin plate spline (TPS) adjusted values of the vein density trait in vein_density_tps_adj.tsv, (iv) best linear unbiased prediction (BLUP) adjusted values of the vein density trait in vein_density_blups.tsv, and (v) GWAS results for the vein density trait, including chromosome positions and corresponding P values, in gwas_results.csv.
The ever-increasing volumes of scientific data combined with sophisticated techniques for extracting information from them have led to the increasing popularity of ensemble workflows which are a collection of runs of individual workflows. A traditional approach followed by scientists to run ensembles is to rely on simple scripts to execute different runs and manage resources. This approach is not scalable and is error-prone, thereby motivating the development of workflow management systems that specialize in executing ensembles on HPC clusters. However, when the size of both the ensemble and the target system reach extreme scales, existing workflow management systems face new challenges that hamper their efficient execution. In this paper, we describe our experience scaling an ensemble workflow from the computational biology domain from the early design stages to the execution at extreme scale on Summit, a leadership class supercomputer at the Oak Ridge National Laboratory. We discuss challenges that arise when scaling ensembles to several million runs on thousands of HPC nodes. We identify challenges with composition of the ensemble itself, its execution at large scale, post-processing of the generated data, and scalability of the file system. Based on the experience acquired, we develop a generic vision of the capabilities and abstractions to add to existing workflow management systems to enable the execution of ensemble workflows at extreme scales. We believe that the understanding of these fundamental challenges will help application teams along with workflow system developers with designing the next generation of infrastructure for composing and executing extreme-scale ensemble workflows.