Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genetic models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

Abstract Background Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity. Results COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer ( Ostrinia nubilalis ). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity. Conclusions The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

GEPA: Reflective Prompt Evolution Can Outperform Reinforcement Learning

Large language models (LLMs) are increasingly adapted to downstream tasks via reinforcement learning (RL) methods like Group Relative Policy Optimization (GRPO), which often require thousands of rollouts to learn new tasks. We argue that the interpretable nature of language often provides a much richer learning medium for LLMs, compared to policy gradients derived from sparse, scalar rewards. To test this, we introduce GEPA (Genetic-Pareto), a prompt optimizer that thoroughly incorporates natural language reflection to learn high-level rules from trial and error. Given any AI system containing one or more LLM prompts, GEPA samples trajectories (e.g., reasoning, tool calls, and tool outputs) and reflects on them in natural language to diagnose problems, propose and test prompt updates, and combine complementary lessons from the Pareto frontier of its own attempts. As a result of GEPA's design, it can often turn even just a few rollouts into a large quality gain. Across six tasks, GEPA outperforms GRPO by 6% on average and by up to 20%, while using up to 35x fewer rollouts. GEPA also outperforms the leading prompt optimizer, MIPROv2, by over 10% (e.g., +12% accuracy on AIME-2025), and demonstrates promising results as an inference-time search strategy for code optimization. We release our code at https://github.com/gepa-ai/gepa.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Combining genome-wide association studies and expression quantitative trait nucleotide mapping with molecular and genetic validations to identify transcriptional networks regulating drought tolerance in Populus

Objectives: (i). To deploy a large-scale experimental drought trial for up to 1000 unique genotypes of Populus equipping the sites with controlled irrigation and drought treatments that are fully automated and monitored. FULLY COMPLETED (ii) To test the hypothesis that a suite of traits identified for drought tolerance in P. nigra can be measured in drought and control treatments in the wide germplasm collection of P. trichocarpa. FULLY COMPLETED (iii) To use established and novel GWAS model approaches to identify gene loci linked to drought tolerance traits on interest in P. trichocarpa. FULLY COMPLETED (iv) To undertake comparative analysis of GWAS results for drought tolerance traits in P. nigra and P. trichocarpa. PARTIALLY COMPLETED – remains active (v) Using RNAseq in P. trichocarpa, in droughted and control treatments to identify cis- and trans-regulated eQTN. FULLY COMPLETED (vi) Validate up to 50 cis-QTNs, from network hubs using transient protoplast assays. FULLY COMPLETED (vii) To establish Agrobacterium-based gene editing protocols in Populus. FULLY COMPLETED (viii) To utilize early leads from previous research to investigate at least 6 candidate genes for drought tolerance in Populus. FULLY COMPLETED (ix) To validate up to 20 candidate genes for drought tolerance in P. trichocarpa refined from the long-list tested in the transient assays for cis-acting hub gene targets. PARTIALLY COMPLETED- remains active.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

High-throughput platform for yeast morphological profiling predicts the targets of bioactive compounds

Abstract Morphological profiling is an omics-based approach for predicting intracellular targets of chemical compounds in which the dose-dependent morphological changes induced by the compound are systematically compared to the morphological changes in gene-deleted cells. In this study, we developed a reliable high-throughput (HT) platform for yeast morphological profiling using drug-hypersensitive strains to minimize compound use, HT microscopy to speed up data generation and analysis, and a generalized linear model to predict targets with high reliability. We first conducted a proof-of-concept study using six compounds with known targets: bortezomib, hydroxyurea, methyl methanesulfonate, benomyl, tunicamycin, and echinocandin B. Then we applied our platform to predict the mechanism of action of a novel diferulate-derived compound, poacidiene. Morphological profiling of poacidiene implied that it affects the DNA damage response, which genetic analysis confirmed. Furthermore, we found that poacidiene inhibits the growth of phytopathogenic fungi, implying applications as an effective antifungal agent. Thus, our platform is a new whole-cell target prediction tool for drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Dynamic nitrogen fixation in an aerobic endophyte of Populus

Biological nitrogen fixation by microbial diazotrophs can contribute significantly to nitrogen availability in non-nodulating plant species. In this study of molecular mechanisms and gene expression relating to biological nitrogen fixation, the aerobic nitrogen-fixing endophyte Burkholderia vietnamiensis, strain WPB, isolated from Populus trichocarpa served as a model for endophyte–poplar interactions. Nitrogen-fixing activity was observed to be dynamic on nitrogen-free medium with a subset of colonies growing to form robust, raised globular like structures. Secondary ion mass spectrometry (NanoSIMS) confirmed that N-fixation was uneven within the population. A fluorescent transcriptional reporter (GFP) revealed that the nitrogenase subunit nifH is not uniformly expressed across genetically identical colonies of WPB and that only ~11% of the population was actively expressing the nifH gene. Higher nifH gene expression was observed in clustered cells through monitoring individual bacterial cells using single-molecule fluorescence in situ hybridization. Through 15 N 2 enrichment, we identified key nitrogenous metabolites and proteins synthesized by WPB and employed targeted metabolomics in active and inactive populations. We cocultivated WPB Pnif-GFP with poplar within a RhizoChip, a synthetic soil habitat, which enabled direct imaging of microbial nifH expression within root epidermal cells. We observed that nifH expression is localized to the root elongation zone where the strain forms a unique physical interaction with the root cells. This work employed comprehensive experimentation to identify novel mechanisms regulating both biological nitrogen fixation and beneficial plant–endophyte interactions.

15 N-tracking metabolomics↗
DOE OSTI2024

An emerging multi-omic understanding of the genetics of opioid addiction

Opioid misuse, addiction, and associated overdose deaths remain global public health crises. Despite the tremendous need for pharmacological treatments, current options are limited in number, use, and effectiveness. Fundamental leaps forward in our understanding of the biology driving opioid addiction are needed to guide development of more effective medication-assisted therapies. This Review focuses on the omics-identified biological features associated with opioid addiction. Recent GWAS have begun to identify robust genetic associations, including variants in OPRM1, FURIN, and the gene cluster SCAI/PPP6C/RABEPK. An increasing number of omics studies of postmortem human brain tissue examining biological features (e.g., histone modification and gene expression) across different brain regions have identified broad gene dysregulation associated with overdose death among opioid misusers. Drawn together by meta-analysis and multi-omic systems biology, and informed by model organism studies, key biological pathways enriched for opioid addiction–associated genes are emerging, which include specific receptors (e.g., GABAB receptors, GPCR, and Trk) linked to signaling pathways (e.g., Trk, ERK/MAPK, orexin) that are associated with synaptic plasticity and neuronal signaling. Studies leveraging the agnostic discovery power of omics and placing it within the context of functional neurobiology will propel us toward much-needed, field-changing breakthroughs, including identification of actionable targets for drug development to treat this devastating brain disease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Cellular response of keratinocytes to the entry and accumulation of nanoplastic particles

Plastic accumulation in the environment is rapidly increasing, and nanoplastics (NP), byproducts of environmental weathering of bulk plastic waste, pose a significant public health risk. Particles may enter the human body through many possible routes such as ingestion, inhalation, and skin absorption. However, studies on NP penetration and accumulation in human skin are limited. Loss or reduction of the keratinized skin barrier may enhance the skin penetration of NPs. The present study investigated the entry of NPs into a human skin system modeling skin with compromised barrier functions and cellular responses to the intracellular accumulations of NPs. Two in vitro models were employed to simulate human skin lacking keratinized barriers. The first model was an ex vivo human skin culture with the keratinized dermal layer (stratum corneum) removed. The second model was a 3D keratinocyte/dermal fibroblast cell co-culture model with stratified keratinocytes on the top and a monolayer of skin fibroblast cells co-cultured at the bottom. The penetration and accumulation of the NPs in different cell types were observed using fluorescent microscopy, confocal microscopy, and cryogenic electron microscopy (cryo-EM). The cellular responses of keratinocytes and dermal fibroblast cells to stress induced by NPs stress were measured. The genetic regulatory pathway of keratinocytes to the intracellular NPs was identified using transcript analyses and KEGG pathway analysis. The cellular uptake of NPs by skin cells was confirmed by imaging analyses. Transepidermal transport and penetration of NPs through the skin epidermis were observed. According to the gene expression and pathway analyses, an IL-17 signaling pathway was identified as the trigger for cellular responses to internal NP accumulation in the keratinocytes. The transepidermal NPs were also found in co-cultured dermal fibroblast cells and resulted in a large-scale transition from fibroblast cells to myofibroblast cells with enhanced production of α-smooth muscle actin and pro-Collagen Ia. The upregulation of inflammatory factors and cell activation may result in skin inflammation and ultimately trigger immune responses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The essential Rhodobacter sphaeroides CenKR two-component system regulates cell division and envelope biosynthesis

Bacterial two-component systems (TCSs) often function through the detection of an extracytoplasmic stimulus and the transduction of a signal by a transmembrane sensory histidine kinase. This kinase then initiates a series of reversible phosphorylation modifications to regulate the activity of a cognate, cytoplasmic response regulator as a transcription factor. Several TCSs have been implicated in the regulation of cell cycle dynamics, cell envelope integrity, or cell wall development in Escherichia coli and other well-studied Gram-negative model organisms. However, many α-proteobacteria lack homologs to these regulators, so an understanding of how α-proteobacteria orchestrate extracytoplasmic events is lacking. In this work we identify an essential TCS, CenKR ( C ell en velope K inase and R egulator), in the α-proteobacterium Rhodobacter sphaeroides and show that modulation of its activity results in major morphological changes. Using genetic and biochemical approaches, we dissect the requirements for the phosphotransfer event between CenK and CenR, use this information to manipulate the activity of this TCS in vivo , and identify genes that are directly and indirectly controlled by CenKR in Rb . sphaeroides . Combining ChIP-seq and RNA-seq, we show that the CenKR TCS plays a direct role in maintenance of the cell envelope, regulates the expression of subunits of the Tol-Pal outer membrane division complex, and indirectly modulates the expression of peptidoglycan biosynthetic genes. CenKR represents the first TCS reported to directly control the expression of Tol-Pal machinery genes in Gram-negative bacteria, and we predict that homologs of this TCS serve a similar function in other closely related organisms. We propose that Rb . sphaeroides genes of unknown function that are directly regulated by CenKR play unknown roles in cell envelope biosynthesis, assembly, and/or remodeling in this and other α-proteobacteria.

Lakey, Bryan D. (ORCID:0000000332795242)↗
DOE OSTI2020

NREL's Capabilities in Algae-Based Chemicals, Polymers, and Fuels

NREL is advancing research and development (R&D) to maximize photosynthetic carbon capture into biomass for conversion to a broad portfolio of valuable products. Areas of focus include: (1) Carbon capture improvement through photosynthesis engineering; (2) Biotechnology and genetic toolbox development for broad species portfolio; (3) Conversion and fractionation of algal biomass to create high-value products and fuels; (4) Biological seaweed conversion to volatile fatty acids as fuel and product feedstocks; (5) Novel polyurethane production synthesis from fully renewable algae-based feedstocks; (6) Advanced analytical characterization of biomass feedstocks; and (7) Process techno-economic and life cycle modeling.

algae↗
DOE OSTI2025

A conserved chaperone protein is required for the formation of a noncanonical type VI secretion system spike tip complex

Type VI secretion systems (T6SSs) are dynamic protein nanomachines found in Gram-negative bacteria that deliver toxic effector proteins into target cells in a contact-dependent manner. Prior to secretion, many T6SS effector proteins require chaperones and/or accessory proteins for proper loading onto the structural components of the T6SS apparatus. However, despite their established importance, the precise molecular function of several T6SS accessory protein families remains unclear. In this study, we set out to characterize the DUF2169 family of T6SS accessory proteins. Using gene co-occurrence analyses, we find that DUF2169-encoding genes strictly co-occur with genes encoding T6SS spike complexes formed by valine-glycine repeat protein G (VgrG) and DUF4150 domains. Although structurally similar to Pro-Ala-Ala-Arg (PAAR) domains, “PAAR-like” DUF4150 domains lack PAAR motifs and instead contain a conserved PIPY motif, leading us to designate them PIPY domains. Next, we present both genetic and biochemical evidence that PIPY domains require a cognate DUF2169 protein to form a functional T6SS spike complex with VgrG. This contrasts with canonical PAAR proteins, which bind VgrG on their own to form functional spike complexes. By solving the first crystal structure of a DUF2169 protein, we show that this T6SS accessory protein adopts a novel protein fold. Furthermore, biophysical and structural modeling data suggest that DUF2169 contains a dynamic loop that physically interacts with a hydrophobic patch on the surface of its cognate PIPY domain. Based on these findings, we propose a model whereby DUF2169 proteins function as molecular chaperones that maintain VgrG–PIPY spike complexes in a secretion-competent state prior to their export by the T6SS apparatus.

DUF2169↗
DOE OSTI2025

Breaking the barrier of human-annotated training data for machine learning-aided plant research using aerial imagery

Machine learning (ML) can accelerate biological research. However, the adoption of such tools to facilitate phenotyping based on sensor data has been limited by (i) the need for a large amount of human-annotated training data for each context in which the tool is used and (ii) phenotypes varying across contexts defined in terms of genetics and environment. This is a major bottleneck because acquiring training data is generally costly and time-consuming. This study demonstrates how a ML approach can address these challenges by minimizing the amount of human supervision needed for tool building. A case study was performed to compare ML approaches that examine images collected by an uncrewed aerial vehicle to determine the presence/absence of panicles (i.e. “heading”) across thousands of field plots containing genetically diverse breeding populations of 2 Miscanthus species. Automated analysis of aerial imagery enabled the identification of heading approximately 9 times faster than in-field visual inspection by humans. Leveraging an Efficiently Supervised Generative Adversarial Network (ESGAN) learning strategy reduced the requirement for human-annotated data by 1 to 2 orders of magnitude compared to traditional, fully supervised learning approaches. The ESGAN model learned the salient features of the data set by using thousands of unlabeled images to inform the discriminative ability of a classifier so that it required minimal human-labeled training data. This method can accelerate the phenotyping of heading date as a measure of flowering time in Miscanthus across diverse contexts (e.g. in multistate trials) and opens avenues to promote the broad adoption of ML tools.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

pixelvar79/ESGAN-Flowering-Detection-paper

Machine learning (ML) can accelerate biological research. However, the adoption of such tools to facilitate phenotyping based on sensor data has been limited by (i) the need for a large amount of human-annotated training data for each context in which the tool is used and (ii) phenotypes varying across contexts defined in terms of genetics and environment. This is a major bottleneck because acquiring training data is generally costly and time-consuming. This study demonstrates how a ML approach can address these challenges by minimizing the amount of human supervision needed for tool building. A case study was performed to compare ML approaches that examine images collected by an uncrewed aerial vehicle to determine the presence/absence of panicles (i.e. “heading”) across thousands of field plots containing genetically diverse breeding populations of 2 Miscanthus species. Automated analysis of aerial imagery enabled the identification of heading approximately 9 times faster than in-field visual inspection by humans. Leveraging an Efficiently Supervised Generative Adversarial Network (ESGAN) learning strategy reduced the requirement for human-annotated data by 1 to 2 orders of magnitude compared to traditional, fully supervised learning approaches. The ESGAN model learned the salient features of the data set by using thousands of unlabeled images to inform the discriminative ability of a classifier so that it required minimal human-labeled training data. This method can accelerate the phenotyping of heading date as a measure of flowering time in Miscanthus across diverse contexts (e.g. in multistate trials) and opens avenues to promote the broad adoption of ML tools.

Varela, Sebastian↗
DOE OSTI2022

Artificial intelligence based analysis of nanoindentation load–displacement data using a genetic algorithm

In this work, we developed an automated tool, Nanoindentation Neo package for the analysis of nanoindentation load–displacement curves using a Genetic Algorithm (GA) applied to the Oliver-Pharr method (Oliver et al.,1992). For some materials, such as polycrystalline isotropic graphites, Least Squares Fitting (LSF) of the unload curve can produce unrealistic fit parameters. These graphites exhibit sharply peaked unloading curves not easily fit using the LSF, which tends to overestimate the indenter tip geometry parameter. To tackle this problem, we extended our general materials characterization tool Neo for EXAFS analysis (Terry et al., 2021) to fit nanoindentation data. Nanoindentation Neo automatically processes and analyzes nanoindentation data with minimal user input while producing meaningful fit parameters. GA, a robust metaheuristic method, begins with a population of temporary solutions using model parameters called chromosomes; from these we evaluate a fitness value for each solution, and select the best solutions to mix with random solutions producing the next generation. A mutation operator then modifies existing solutions by random perturbations, and the optimal solution is selected. We tested the GA method using Silica and Al reference standards. We fit samples of graphite and a high entropy alloy (HEA) consisting of BCC and FCC phases.

42 ENGINEERING↗
DOE OSTI2021

Systematic discovery of pseudomonad genetic factors involved in sensitivity to tailocins

Tailocins are bactericidal protein complexes produced by a wide variety of bacteria that kill closely related strains and may play a role in microbial community structure. Thanks to their high specificity, tailocins have been proposed as precision antibacterial agents for therapeutic applications. Compared to tailed phages, with whom they share an evolutionary and morphological relationship, bacterially produced tailocins kill their host upon production but producing strains display resistance to self-intoxication. Though lipopolysaccharide (LPS) has been shown to act as a receptor for tailocins, the breadth of factors involved in tailocin sensitivity, and the mechanisms behind resistance to self-intoxication, remain unclear. Here, we employed genome-wide screens in four non-model pseudomonads to identify mutants with altered fitness in the presence of tailocins produced by closely related pseudomonads. Our mutant screens identified O-antigen composition and display as most important in defining sensitivity to our tailocins. In addition, the screens suggest LPS thinning as a mechanism by which resistant strains can become more sensitive to tailocins. Furthermore, we validate many of these novel findings, and extend these observations of tailocin sensitivity to 130 genome-sequenced pseudomonads. This work offers insights into tailocin–bacteria interactions, informing the potential use of tailocins in microbiome manipulation and antibacterial therapy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Abstract Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Genotoxicity evaluation using primary hepatocytes isolated from rhesus macaque ( Macaca mulatta )

Non-human primates (NHPs) have played a vital role in fundamental, pre-clinical, and translational studies because of their high physiological and genetic similarity to humans. We report a method to isolate primary hepatocytes from the livers of rhesus macaques (Macaca mulatta) after in situ whole liver perfusion. Isolated primary macaque hepatocytes (PMHs) were treated with various compounds known to have different pathways of genotoxicity/carcinogenicity and the resulting DNA damage was evaluated using the high-throughput CometChip assay. The comet data were quantified using benchmark dose (BMD) modeling and the BMD50 values for treatments of PMHs were compared with those generated from primary human hepatocytes (PHHs) in our previous study. The results showed that despite varying CYP450 enzyme activities, PMHs had the same sensitivity and specificity as PHHs in detecting four indirect-acting (i.e., requiring metabolic activation) and seven direct-acting genotoxicants/carcinogens, as well as five non-carcinogens that are negative or equivocal for genotoxicity in vivo. The BMD50 estimates and their confidence intervals revealed species differences for DNA damage potency, especially for direct-acting compounds. The present study provides a practical method for maximizing the use of animal tissues by isolating primary hepatocytes from NHPs. Our data support the use of PMHs as a reliable surrogate of PHHs for evaluating the genotoxic hazards of chemical substances for humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Gene loss during a transition to multicellularity

Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited novel interactions, suggesting that gene loss could play a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest gene loss, in addition to gene gain and co-option, may be important for the evolution developmental complexity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

54 ENVIRONMENTAL SCIENCES↗