Search NASA⌕ Search

SEARCH · Search NASA

Results for “Genetic prediction”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Intermediate Molecular Phenotypes to Identify Genetic Markers of Anthracycline-Induced Cardiotoxicity Risk

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Training Population Optimization for Genomic Selection in Miscanthus

Miscanthus is a perennial grass with potential for lignocellulosic ethanol production. To ensure its utility for this purpose, breeding efforts should focus on increasing genetic diversity of the nothospecies Miscanthus × giganteus (M×g) beyond the single clone used in many programs. Germplasm from the corresponding parental species M. sinensis (Msi) and M. sacchariflorus (Msa) could theoretically be used as training sets for genomic prediction of M×g clones with optimal genomic estimated breeding values for biofuel traits. To this end, we first showed that subpopulation structure makes a substantial contribution to the genomic selection (GS) prediction accuracies within a 538-member diversity panel of predominately Msi individuals and a 598-member diversity panels of Msa individuals. We then assessed the ability of these two diversity panels to train GS models that predict breeding values in an interspecific diploid 216-member M×g F2 panel. Low and negative prediction accuracies were observed when various subsets of the two diversity panels were used to train these GS models. To overcome the drawback of having only one interspecific M×g F2 panel available, we also evaluated prediction accuracies for traits simulated in 50 simulated interspecific M×g F2 panels derived from different sets of Msi and diploid Msa parents. The results revealed that genetic architectures with common causal mutations across Msi and Msa yielded the highest prediction accuracies. Ultimately, these results suggest that the ideal training set should contain the same causal mutations segregating within interspecific M×g populations, and thus efforts should be undertaken to ensure that individuals in the training and validation sets are as closely related as possible.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Identification of carbohydrate gene clusters obtained from in vitro fermentations as predictive biomarkers of prebiotic responses

Prebiotic fibers are non-digestible substrates that modulate the gut microbiome by promoting expansion of microbes having the genetic and physiological potential to utilize those molecules. Although several prebiotic substrates have been consistently shown to provide health benefits in human clinical trials, responder and non-responder phenotypes are often reported. These observations had led to interest in identifying, a priori, prebiotic responders and non-responders as a basis for personalized nutrition. In this study, we conducted in vitro fecal enrichments and applied shotgun metagenomics and machine learning tools to identify microbial gene signatures from adult subjects that could be used to predict prebiotic responders and non-responders. Using short chain fatty acids as a targeted response, we identified genetic features, consisting of carbohydrate active enzymes, transcription factors and sugar transporters, from metagenomic sequencing of in vitro fermentations for three prebiotic substrates: xylooligosacharides, fructooligosacharides, and inulin. A machine learning approach was then used to select substrate-specific gene signatures as predictive features. These features were found to be predictive for XOS responders with respect to SCFA production in an in vivo trial. Our results confirm the bifidogenic effect of commonly used prebiotic substrates along with inter-individual microbial responses towards these substrates. We successfully trained classifiers for the prediction of prebiotic responders towards XOS and inulin with robust accuracy (≥ AUC 0.9) and demonstrated its utility in a human feeding trial. Overall, the findings from this study highlight the practical implementation of pre-intervention targeted profiling of individual microbiomes to stratify responders and non-responders.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Empirical evidence that glucan-interacting amino acid side chains within the transmembrane channel collectively facilitate cellulose synthase function

The fundamental mechanism of cellulose synthesis is widely conserved across Kingdoms and depends on cellulose synthases, which are processive, dual-function, family 2 glycosyltransferases (GT-2). These enzymes polymerize glucose on the cytoplasmic side of the plasma membrane and export the glucan chain to the cell surface through an integral transmembrane (TM) channel. Structural studies of active plant cellulose synthases (CESAs) have revealed interactions between the nascent glucan chain and the side chains of polar, charged, and aromatic amino acid residues that line the TM channel. However, the functional consequences of modifying these side chains have not been tested in vivo in CESAs or other processive GT-2s. To test this, we used an established in vivo assay based on genetic complementation of CESA5 in the moss, Physcomitrium patens. For accurate prediction of glucan-interacting amino acid residues, we generated a complete homotrimeric molecular model of PpCESA5 using a combination of homology and de novo modeling. All-atom molecular dynamics-based analyses of contact metrics and interaction energy identified 23 amino acid residues with high propensity to interact with the nascent glucan chain within the TM channel or on the apoplastic surface of PpCESA5. Mutating any one of 18 of these amino acid residues to alanine, thereby removing their side chains, abolished or impaired CESA function, with the strongest effects observed upon the loss of charged amino acid side chains. This provides direct evidence to support the hypothesis that multiple amino acid residues collectively maintain a smooth energy landscape within the TM channel to facilitate glucan translocation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Transgenic Poplar Designed for Biofuels

Members of the genus Populus (i.e., cottonwood, hybrid poplar) represent a promising source of lignocellulosic biomass for biofuels. However, one of the major factors negatively affecting poplar’s efficient conversion to biofuel is the inherent recalcitrance to enzymatic saccharification due to cell wall components such as lignin. To this effect, there have been efforts to modify gene expression to reduce biomass recalcitrance by changing cell wall properties. Furthermore, we review recent genetic modifications of poplar that led to change cell wall properties and the resulting effects on subsequent pretreatment efficacy and saccharification. Although genetic engineering’s impacts on cell wall properties are not fully predictable, recent studies have shown promising improvement in the biological conversion of transgenic poplar to biofuels.

09 BIOMASS FUELS↗
DOE OSTI2025

Intraspecific Reaction Norm Variation Controls the Eco-Evolutionary Consequences of Environmental Change

As environmental change accelerates globally, understanding concurrent organismal, species, and community responses is increasingly vital. Here, we examine these collective responses by incorporating genotype-specific thermal reaction norms into an eco-evolutionary predator-prey model, allowing us to track simultaneous phenotypic, ecological, and evolutionary responses to environmental change within ecological communities. We show that the reaction norms expressed by genotypes within a population determine how a community switches between different eco-evolutionary outcomes with changes in temperature. We identify how different components of phenotypic variation in thermal reaction norms—environmental (E), additive environmental and genetic (E + G), and gene-by-environment interactions (G × E)—influence eco-evolutionary dynamics and outcomes as temperature changes. Furthermore, our findings underscore how complex eco-evolutionary responses to environmental change ultimately emerge from variation in reaction norms among genotypes, offering new mechanistic insights into environmental impacts on adaptation, the maintenance of phenotypic and genetic variation, and ecological stability, which is crucial for understanding and predicting eco-evolutionary effects of rapid environmental change in the future.

Eco-evolutionary↗
DOE OSTI2024

Distinct function of Chlamydomonas CTRA-CTR transporters in Cu assimilation and intracellular mobilization

Abstract Successful acclimation to copper (Cu) deficiency involves a fine balance between Cu import and export. In the green alga Chlamydomonas reinhardtii, Cu import is dependent on a transcription factor, Copper Response Regulator 1 (CRR1), responsible for activating genes in Cu-deficient cells. Among CRR1 target genes are two Cu transporters belonging to the CTR/COPT gene family (CTR1 and CTR2) and a related soluble protein (CTR3). The ancestor of these green algal proteins was likely acquired from an ancient chytrid and contained conserved cysteine-rich domains (named the CTR-associated domains, CTRA) that are predicted to be involved in Cu acquisition. We show by reverse genetics that Chlamydomonas CTR1 and CTR2 are canonical Cu importers albeit with distinct affinities, while loss of CTR3 did not result in an observable phenotype under the conditions tested. Mutation of CTR1, but not CTR2, recapitulates the poor growth of crr1 in Cu-deficient medium, consistent with a dominant role for CTR1 in high-affinity Cu(I) uptake. On the other hand, the overaccumulation of Cu(I) (20 times the quota) in zinc (Zn) deficiency depends on CRR1 and both CTR1 and CTR2. CRR1-dependent activation of CTR gene expression needed for Cu over-accumulation can be bypassed by the provision of excess Cu in the growth medium. Over-accumulated Cu is sequestered into the acidocalcisome but can become remobilized by restoring Zn nutrition. This mobilization is also CRR1-dependent, and requires activation of CTR2 expression, again distinguishing CTR2 from CTR1 and consistent with the lower substrate affinity of CTR2. One sentence summary Regulation of Cu uptake and sequestration by members of the CTR family of proteins in Chlamydomonas.

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

A susceptibility gene signature for ERBB2-driven mammary tumour development and metastasis in collaborative cross mice

Background: Deeper insights into ERBB2-driven cancers are essential to develop new treatment approaches for ERBB2+ breast cancers (BCs). We employed the Collaborative Cross (CC) mouse model to unearth genetic factors underpinning Erbb2-driven mammary tumour development and metastasis. Methods: 732 F1 hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains were monitored for mammary tumour phenotypes. GWAS pinpointed SNPs that influence various tumour phenotypes. Multivariate analyses and models were used to construct the polygenic score and to develop a mouse tumour susceptibility gene signature (mTSGS), where the corresponding human ortholog was identified and designated as hTSGS. The importance and clinical value of hTSGS in human BC was evaluated using public datasets, encompassing TCGA, METABRIC, GSE96058, and I-SPY2 cohorts. The predictive power of mTSGS for response to chemotherapy was validated in vivo using genetically diverse MMTV-Erbb2 mice. Findings: Distinct variances in tumour onset, multiplicity, and metastatic patterns were observed in F1-hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains. Besides lung metastasis, liver and kidney metastases emerged in specific CC strains. GWAS identified specific SNPs significantly associated with tumour onset, multiplicity, lung metastasis, and liver metastasis. Multivariate analyses flagged SNPs in 20 genes (Stx6, Ramp1, Traf3ip1, Nckap5, Pfkfb2, Trmt1l, Rprd1b, Rer1, Sepsecs, Rhobtb1, Tsen15, Abcc3, Arid5b, Tnr, Dock2, Tti1, Fam81a, Oxr1, Plxna2, and Tbc1d31) independently tied to various tumour characteristics, designated as a mTSGS. hTSGS scores (hTSGSS) based on their transcriptional level showed prognostic values, superseding clinical factors and PAM50 subtype across multiple human BC cohorts, and predicted pathological complete response independent of and superior to MammaPrint score in I-SPY2 study. The power of mTSGS score for predicting chemotherapy response was further validated in an in vivo mouse MMTV-Erbb2 model, showing that, like findings in human patients, mouse tumours with low mTSGS scores were most likely to respond to treatment. Interpretation: Our investigation has unveiled many new genes predisposing individuals to ERBB2-driven cancer. Translational findings indicate that hTSGS holds promise as a biomarker for refining treatment strategies for patients with BC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Genetic program activity delineates risk, relapse, and therapy responsiveness in multiple myeloma

Despite recent advancements in the treatment of multiple myeloma (MM), nearly all patients ultimately relapse and many become refractory to multiple lines of therapies. Therefore, we not only need the ability to predict which patients are at high risk for disease progression but also a means to understand the mechanisms underlying their risk. Here, we report a transcriptional regulatory network (TRN) for MM inferred from cross-sectional multi-omics data from 881 patients that predicts how 124 chromosomal abnormalities and somatic mutations causally perturb 392 transcription regulators of 8549 genes to manifest in distinct clinical phenotypes and outcomes. We identified 141 genetic programs whose activity profiles stratify patients into 25 distinct transcriptional states and proved to be more predictive of outcomes than did mutations. The coherence of these programs and accuracy of our network-based risk prediction was validated in two independent datasets. We observed subtype-specific vulnerabilities to interventions with existing drugs and revealed plausible mechanisms for relapse, including the establishment of an immunosuppressive microenvironment. Investigation of the t(4;14) clinical subtype using the TRN revealed that 16% of these patients exhibit an extreme-risk combination of genetic programs (median progression-free survival of 5 months) that create a distinct phenotype with targetable genes and pathways.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Genome-wide functional screens enable the prediction of high activity CRISPR-Cas9 and -Cas12a guides in Yarrowia lipolytica

Abstract Genome-wide functional genetic screens have been successful in discovering genotype-phenotype relationships and in engineering new phenotypes. While broadly applied in mammalian cell lines and in E. coli , use in non-conventional microorganisms has been limited, in part, due to the inability to accurately design high activity CRISPR guides in such species. Here, we develop an experimental-computational approach to sgRNA design that is specific to an organism of choice, in this case the oleaginous yeast Yarrowia lipolytica . A negative selection screen in the absence of non-homologous end-joining, the dominant DNA repair mechanism, was used to generate single guide RNA (sgRNA) activity profiles for both SpCas9 and LbCas12a. This genome-wide data served as input to a deep learning algorithm, DeepGuide, that is able to accurately predict guide activity. DeepGuide uses unsupervised learning to obtain a compressed representation of the genome, followed by supervised learning to map sgRNA sequence, genomic context, and epigenetic features with guide activity. Experimental validation, both genome-wide and with a subset of selected genes, confirms DeepGuide’s ability to accurately predict high activity sgRNAs. DeepGuide provides an organism specific predictor of CRISPR guide activity that with retraining could be applied to other fungal species, prokaryotes, and other non-conventional organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Bridging Experiment and Theory to Reveal Compounds in K–Zn(Cd)–Bi Systems

This study investigates the facile hydride synthesis method guided by theoretical predictions to explore the K–T–Bi (T = Zn, Cd) phase spaces. Using an adaptive genetic algorithm (AGA) and density functional theory (DFT), candidate compositions are identified for experimental validation via a facile hydrides route, permitting experimental screening of K–Zn–Bi and “empty” K–Cd–Bi systems. The previously reported KZnBi and KZn 2 Bi 2 are synthesized alongside newly discovered KCdBi and KCd 2 Bi 2 . While the AGA and DFT predict the stability of these compounds, structural predictions align with the experiment only for KZnBi and KZn 2 Bi 2 . Single-crystal X-ray structure refinements confirm that KZnBi and KZn 2 Bi 2 adopt the hexagonal ZrBeSi- and tetragonal ThCr 2 Si 2 -structure types, respectively. KCdBi has tetragonal PbClF-structure type and KCd 2 Bi 2 belongs to the ThCr 2 Si 2 -structure type. A trend based on the ratio of the metal ionic radii allows to rationalize variation in the structure types within the ATBi family (A = Li–Cs), correctly identifying KCdBi as isostructural to NaZnBi. Thermal stability studied by high-temperature powder X-ray diffraction reveals that Zn-containing compounds melt at higher temperatures (821 K for KZn 2 Bi 2 ) than Cd-containing KCd 2 Bi 2 (635 K). This study highlights the efficacy of combining rapid synthesis techniques with predictive modeling, though structural predictions show some limitations in accuracy.

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

FluxRETAP: a REaction TArget Prioritization genome-scale modeling technique for selecting genetic targets

MOTIVATION: Metabolic engineering is rapidly evolving as a result of new advances in synthetic biology tools and automation platforms that enable high throughput strain construction, as well as the development of machine learning tools (ML) for biology. However, selecting genetic engineering targets that effectively guide the metabolic engineering process is still challenging. ML can provide predictive power for synthetic biology, but current technical limitations prevent the independent use of ML approaches without previous biological knowledge. RESULTS: Here, we present FluxRETAP, a simple and computationally inexpensive method that leverages the prior mechanistic knowledge embedded in genome-scale models for suggesting targets for genetic overexpression, downregulation or deletion, with the final goal of increasing the production of a desired metabolite. This method can provide a list of desirable engineering targets that can be combined with current ML pipelines. FluxRETAP captured 100% of reaction targets experimentally verified to improve Escherichia coli isoprenol production, 50% of targets that experimentally improved taxadiene production in E. coli and ∼60% of genetic targets from a verified minimal constrained cut-set in Pseudomonas putida, while providing additional high priority targets that could be tested. Overall, FluxRETAP is an efficient algorithm for identifying a prioritized list of testable genetic and reaction targets. AVAILABILITY AND IMPLEMENTATION: FluxRETAP is implemented in python and released under the creative commons license. The implementation and code are freely available at: https://github.com/JBEI/FluxRETAP.

Czajka, Jeffrey J↗
DOE OSTI2025

Tutorial: Machine-Learning-Based CREASE-2D Analysis of 2D SAXS Profiles to Characterize Anisotropic Nanostructures in Soft Materials

We present a tutorial to guide users on how to extend the Computational Reverse Engineering Analysis of Scattering Experiments-2D (CREASE-2D) framework to interpret their experimental two-dimensional small-angle scattering (SAS) data from soft materials (e.g., polymers, peptide amphiphiles, biomolecular fibrils). Unlike most traditional SAS analysis approaches, which typically rely on azimuthally averaged onedimensional (1D) profiles, CREASE-2D utilizes the complete 2D scattering profile to reveal information about anisotropy in the structure. In past applications, CREASE has provided insights into complex structural features, including the cross-sectional shapes of assembled nanostructures and dispersity in these features, which are difficult to discern with existing analytical models. While (1D- ) CREASE has been applied to SANS and SAXS data, this tutorial shares the steps for implementing CREASE-2D using an example of a dipeptide solution system, for which we have SAXS data. We present details for these steps involved in using CREASE-2D to interpret SAXS profiles: how to preprocess SAXS data, define relevant structural features, generate three-dimensional real-space structures for specific values of these features, train a machine learning (ML) surrogate model to predict scattering profiles for given structural features, and optimize these features using genetic algorithms (GA). Then, we use these steps to interpret complex 2DSAXS data collected from dipeptide solutions that, in microscopy images, exhibit nanoscale structures that could be elliptical tubes/ flat tapes/cylinders or a combination of these cross sections. Open-source codes, computational hardware, and software requirements, as well as the strengths and limitations of this protocol, are also presented. We expect researchers working with (soft) biomaterials, peptide amphiphiles, amphiphilic polymer solutions, polymer nanocomposites, and blends of particles/polymers will find this CREASE-2D method and this tutorial of use.

CREASE↗
DOE OSTI2023

Targeting tissues via dynamic human systems modeling in generative design

Drug discovery is a complex, costly process with high failure rates. A successful drug should bind to a target, be deliverable to an intended site of activity, and promote a desired pharmacological effect without causing toxicity. Typically, these factors are evaluated in series over the course of a pipeline where the number of candidates is sequentially whittled down from a very large initial pool. One promise of AI-driven discovery is the opportunity to evaluate multiple facets of drug performance in parallel. However, despite ML-driven advancements, current models for pharmacological property prediction are exclusively trained to predict molecular properties, ignoring important, dynamic biodistribution and bioactivity effects. Here, we present our progress towards incorporating quantitative systems physiology models into an AI-enabled drug design and molecular generation pipeline. Within a genetic algorithm, we include human-relevant physiologically based pharmacokinetic (PBPK) models. These PBPK models leverage properties that are predicted by a fine-tuned molecular language model. Together, these models will aid in capturing the mapping between molecules and therapeutic outcomes that is necessary to accelerate the drug discovery process.

Fox, Zach↗
DOE OSTI2021

Gut microbiome partially mediates and coordinates the effects of genetics on anxiety-like behavior in Collaborative Cross mice

Abstract Growing evidence suggests that the gut microbiome (GM) plays a critical role in health and disease. However, the contribution of GM to psychiatric disorders, especially anxiety, remains unclear. We used the Collaborative Cross (CC) mouse population-based model to identify anxiety associated host genetic and GM factors. Anxiety-like behavior of 445 mice across 30 CC strains was measured using the light/dark box assay and documented by video. A custom tracking system was developed to quantify seven anxiety-related phenotypes based on video. Mice were assigned to a low or high anxiety group by consensus clustering using seven anxiety-related phenotypes. Genome-wide association analysis (GWAS) identified 141 genes (264 SNPs) significantly enriched for anxiety and depression related functions. In the same CC cohort, we measured GM composition and identified five families that differ between high and low anxiety mice. Anxiety level was predicted with 79% accuracy and an AUC of 0.81. Mediation analyses revealed that the genetic contribution to anxiety was partially mediated by the GM. Our findings indicate that GM partially mediates and coordinates the effects of genetics on anxiety.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

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

Evaluating genetic and environmental controls on silicon accumulation in Populus trichocarpa leaves

Background and aims: Silicon (Si) accumulation and formation of non-crystalline silica particles (phytoliths) inside plant cells contribute to plant function, soil carbon storage, and Si reservoirs in soils. A better understanding of how genetic and environmental factors control Si accumulation may improve understanding of Si cycling between plants and soil and its impacts on other ecosystem processes. Methods: This study examined how genetic potential for Si uptake and silicification translated to Si concentrations and phytolith formation in Populus trichocarpa leaves. Leaves were collected from four different genotypes of cultivated poplar trees that were grown in a common garden and possess established differences in leaf-level gene expression related to Si transport and silicification. Scanning electron microscopy and subsequent image analysis were used to quantify Si concentrations and phytolith abundance and features. Results: Phytoliths were identified in nearly all leaves but were often absent at very low (< 0.5 wt.%) Si concentrations at the leaf surface. The proportion of leaf area occupied by phytoliths increased as Si content increased. Silicon content and phytolith abundance at the leaf surface did not vary across genotypes, but surface Si was weakly correlated with plant available Si in soil. Conclusion: Phytoliths were prevalent in P. trichocarpa leaves; however, the genetic potential for Si uptake and silicification at the leaf-level poorly predicted foliar Si and phytolith abundance. Uniformly low expression of Si uptake genes in the roots may have reduced differences across genotypes. Consequently, Si accumulation in leaves was likely more strongly influenced by soil Si bioavailability, phenology, or additional environmental factors.

Image analysis↗