Engineering topics

Pavicic, Mirko

Publications and source records attributed to Pavicic, Mirko.

8 records

DOE OSTI2024

Analyses of GWAS signal using GRIN identify additional genes contributing to suicidal behavior

Genome-wide association studies (GWAS) identify genetic variants underlying complex traits but are limited by stringent genome-wide significance thresholds. We present GRIN (Gene set Refinement through Interacting Networks), which increases confidence in the expanded gene set by retaining genes strongly connected by biological networks when GWAS thresholds are relaxed. GRIN was validated on both simulated interrelated gene sets as well as multiple GWAS traits. From multiple GWAS summary statistics of suicide attempt, a complex phenotype, GRIN identified additional genes that replicated across independent cohorts and retained biologically interrelated genes despite a relaxed significance threshold. We present a conceptual model of how these retained genes interact through neurobiological pathways that may influence suicidal behavior, and identify existing drugs associated with these pathways that would not have been identified under traditional GWAS thresholds. We demonstrate GRIN’s utility in boosting GWAS results by increasing the number of true positive genes identified from GWAS results.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

RWRtoolkit Narrative

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.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Transcriptomic Analysis Reveals Novel Regulators of the Scots Pine Stilbene Pathway

Stilbenes accumulate in Scots pine heartwood where they have important roles in protecting wood from decaying fungi. They are also part of active defense responses, and their production is induced by different (a)biotic stressors. The specific transcriptional regulators as well as the enzyme responsible for activating the stilbene precursor cinnamate in the pathway are still unknown. UV-C radiation was the first discovered artificial stress activator of the pathway. Here, we describe a large-scale transcriptomic analysis of pine needles in response to UV-C and treatment with translational inhibitors, both activating the transcription of stilbene pathway genes. We used the data to identify putative candidates for the missing CoA ligase and for pathway regulators. We further showed that the pathway is transcriptionally activated by phosphatase inhibitor, ethylene and jasmonate treatments, as in grapevine, and that the stilbene synthase promoter retains its inducibility in some of the tested conditions in Arabidopsis, a species that normally does not synthesize stilbenes. Shared features between gymnosperm and angiosperm regulation and partially retained inducibility in Arabidopsis suggest that pathway regulation occurs not only via ancient stress-response pathway(s) but also via species-specific regulators. Understanding which genes control the biosynthesis of stilbenes in Scots pine aids breeding of more resistant trees.

4-Coumarate-CoA ligase↗

DOE OSTI2023

Few-Shot Learning Enables Population-Scale Analysis of Leaf Traits in Populus trichocarpa

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.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Image-based time series analysis to establish differential disease progression for two Fusarium head blight pathogens in oat spikelets with variable resistance

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.

oats↗

DOE OSTI2023

Supporting information for Few-Shot Learning Enables Population-Scale Analysis of Leaf Traits in Populus trichocarpa

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.

09 BIOMASS FUELS↗

DOE OSTI2021

Image-Based Methods to Score Fungal Pathogen Symptom Progression and Severity in Excised Arabidopsis Leaves

Image-based symptom scoring of plant diseases is a powerful tool for associating disease resistance with plant genotypes. Advancements in technology have enabled new imaging and image processing strategies for statistical analysis of time-course experiments. There are several tools available for analyzing symptoms on leaves and fruits of crop plants, but only a few are available for the model plant Arabidopsis thaliana (Arabidopsis). Arabidopsis and the model fungus Botrytis cinerea (Botrytis) comprise a potent model pathosystem for the identification of signaling pathways conferring immunity against this broad host-range necrotrophic fungus. Here, we present two strategies to assess severity and symptom progression of Botrytis infection over time in Arabidopsis leaves. Thus, a pixel classification strategy using color hue values from red-green-blue (RGB) images and a random forest algorithm was used to establish necrotic, chlorotic, and healthy leaf areas. Secondly, using chlorophyll fluorescence (ChlFl) imaging, the maximum quantum yield of photosystem II (Fv/Fm) was determined to define diseased areas and their proportion per total leaf area. Both RGB and ChlFl imaging strategies were employed to track disease progression over time. This has provided a robust and sensitive method for detecting sensitive or resistant genetic backgrounds. A full methodological workflow, from plant culture to data analysis, is described.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

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

Abstract 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 Asp 614 Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro 323 Leu 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.

59 BASIC BIOLOGICAL SCIENCES↗