Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical and quantitative genetics”

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.

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↗

Predictive models of the genetic bases underlying budding yeast fitness in multiple environments

Abstract The ability of organisms to adapt and survive depends on the effects of genes and the environment on fitness. However, the multigenic nature of fitness and genotype-by-environment interactions hinder our understanding of the genetic basis of fitness. Here, we established fitness prediction models for 35 environments using machine learning and existing fitness data and different genetic variant types for a Saccharomyces cerevisiae population. Models revealed that the predictive ability of genetic variants varied across environments, with copy number variants explaining the majority of fitness variation in most cases. Model interpretation showed that different variant types identified distinct gene sets associated with predictive variants. These gene sets were significantly enriched in experimentally validated genes affecting fitness in only a subset of environments, indicating that many genes influencing fitness remain unexplored. Notably, non-experimentally validated genes were more important than validated ones for fitness predictions. Gene contributions to predictions were both isolate- and environment-dependent, pointing to gene-by-gene and gene-by-environment interactions. Furthermore, models uncovered experimentally validated and novel candidate genetic interactions for a well-characterized stress, the fungicide benomyl. These findings highlight the feasibility of identifying the genetic basis of fitness by using different genetic variant types and offer novel targets for future functional analysis.

DNA copy number variations↗

Unraveling the Genetics of Two Key Biomass Traits that Differentiate Upland and Lowland Tetraploid Switchgrass Ecotypes, Colonization by Mycorrhizal Fungi and Frost Tolerance

The overall objectives of the proposal were to dissect the genetics of the ability of switchgrass to interact with arbuscular mycorrhizal fungi (AMF) and to tolerate subzero temperatures, and to investigate the effect of AMF colonization on freezing tolerance and biomass production. In order to achieve this, we generated genetic maps in an F2 mapping population derived from a cross between the lowland genotype AP13 and the upland genotype VS16. In addition to providing a framework for identification of quantitative trait loci (QTL), the genetic maps were used to validate and improve the switchgrass genome assembly. We mapped phenotypic quantitative trait loci (QTL) for AMF colonization and cold tolerance in the F 2 population. We demonstrated that identification of QTL for AMF colonization was achieved only if AMF genus-specific colonization levels were used as phenotypes. No QTL were obtained when the overall level of AMF colonization was used as a trait. We conducted RNASeq on a subset of the F2 progeny following cold acclimation, and developed statistical methods to conduct expression QTL (eQTL) analyses in a polyploid organism. A preliminary analysis of the results identified eQTL for a range of genes that were differentially regulated in the lowland AP13 compared to the upland VS16 under cold acclimation. Once sequencing of the transcriptomes of a temporal replicate of the mapping population that was cold-acclimated using the same conditions has been completed, a final eQTL analysis will be conducted.

09 BIOMASS FUELS↗

Genome-Scale Transcription-Translation Mapping Reveals Features of Zymomonas mobilis Transcription Units and Promoters

Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis . This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

59 BASIC BIOLOGICAL SCIENCES↗

GWAS supported by computer vision identifies large numbers of candidate regulators of in planta regeneration in Populus trichocarpa

Plant regeneration is an important dimension of plant propagation and a key step in the production of transgenic plants. However, regeneration capacity varies widely among genotypes and species, the molecular basis of which is largely unknown. Association mapping methods such as genome-wide association studies (GWAS) have long demonstrated abilities to help uncover the genetic basis of trait variation in plants; however, the performance of these methods depends on the accuracy and scale of phenotyping. To enable a large-scale GWAS of in planta callus and shoot regeneration in the model tree Populus, we developed a phenomics workflow involving semantic segmentation to quantify regenerating plant tissues over time. We found that the resulting statistics were of highly non-normal distributions, and thus employed transformations or permutations to avoid violating assumptions of linear models used in GWAS. We report over 200 statistically supported quantitative trait loci (QTLs), with genes encompassing or near to top QTLs including regulators of cell adhesion, stress signaling, and hormone signaling pathways, as well as other diverse functions. Our results encourage models of hormonal signaling during plant regeneration to consider keystone roles of stress-related signaling (e.g. involving jasmonates and salicylic acid), in addition to the auxin and cytokinin pathways commonly considered. The putative regulatory genes and biological processes we identified provide new insights into the biological complexity of plant regeneration, and may serve as new reagents for improving regeneration and transformation of recalcitrant genotypes and species.

59 BASIC BIOLOGICAL SCIENCES↗

Automated segmentation of soft X-ray tomography: Native cellular structure with submicron resolution at high-throughput for whole-cell quantitative imaging in yeast

Soft X-ray tomography (SXT) is an invaluable tool for quantitatively analyzing cellular structures at suboptical isotropic resolution. However, it has traditionally depended on manual segmentation, limiting its scalability for large datasets. Here, we leverage a deep learning-based autosegmentation pipeline to segment and label cellular structures in hundreds of cells across three Saccharomyces cerevisiae strains. This task-based pipeline uses manual iterative refinement to improve segmentation accuracy for key structures, including the cell body, nucleus, vacuole, and lipid droplets, enabling high-throughput and precise phenotypic analysis. Using this approach, we quantitatively compared the three-dimensional (3D) whole-cell morphometric characteristics of wild-type, VPH1-GFP, and vac14 strains, uncovering detailed strain-specific cell and organelle size and shape variations. We show the utility of SXT data for precise 3D curvature analysis of entire organelles and cells and detection of fine morphological features using surface meshes. Our approach facilitates comparative analyses with high spatial precision and statistical throughput, uncovering subtle morphological features at the single-cell and population level. This workflow significantly enhances our ability to characterize cell anatomy and supports scalable studies on the mesoscale, with applications in investigating cellular architecture, organelle biology, and genetic research across diverse biological contexts.

Chen, Jianhua [Lawrence Berkeley National Laborato↗

Influence of ontogenetic development, temperature, and pCO2 on otolith calcium carbonate polymorph composition in sturgeons

Changes to calcium carbonate (CaCO 3 ) biomineralization in aquatic organisms is among the many predicted effects of climate change. Because otolith (hearing/orientation structures in fish) CaCO 3 precipitation and polymorph composition are controlled by genetic and environmental factors, climate change may be predicted to affect the phenotypic plasticity of otoliths. We examined precipitation of otolith polymorphs (aragonite, vaterite, calcite) during early life history in two species of sturgeon, Lake Sturgeon, (Acipenser fulvescens) and White Sturgeon (A. transmontanus), using quantitative X-ray microdiffraction. Both species showed similar fluctuations in otolith polymorphs with a significant shift in the proportions of vaterite and aragonite in sagittal otoliths coinciding with the transition to fully exogenous feeding. We also examined the effect of the environment on otolith morphology and polymorph composition during early life history in Lake Sturgeon larvae reared in varying temperature (16/22 °C) and pCO 2 (1000/2500 µatm) environments for 5 months. Fish raised in elevated temperature had significantly increased otolith size and precipitation of large single calcite crystals. Interestingly, pCO 2 had no statistically significant effect on size or polymorph composition of otoliths despite blood pH exhibiting a mild alkalosis, which is contrary to what has been observed in several studies on marine fishes. These results suggest climate change may influence otolith polymorph composition during early life history in Lake Sturgeon.

54 ENVIRONMENTAL SCIENCES↗

Top-down proteomics

Proteoforms arising from posttranslational modifications, genetic polymorphisms, and RNA splice variants, play a pivotal role as the key drivers in biology. Thus, a comprehensive understanding of proteoforms is essential for unraveling the intricacies of biological systems and bridging the gap between genotype and phenotype. By analyzing whole proteins without digestion, top-down proteomics (TDP) provides a holistic view of the proteome and presents a next-generation approach for deciphering protein function, uncovering disease mechanisms, and advancing precision medicine. This Primer embarks on a journey into the world of TDP by encapsulating its historical context, underlying principles, recent advances, and an outlook on the future of TDP. The experimental section navigates instrumentation, sample preparation, intact protein separation, tandem mass spectrometry techniques, and data collection. Results decipher raw data, visualize intact protein spectra, unravel data analysis, and explain proteoform identification, characterization, and quantitation, as well as statistical analysis. Various applications of TDP spanning the human proteoform project, biomedical, biopharmaceutical, and clinical applications are described. These are complemented by discussions on measurement reproducibility, limitations, and a forward-looking perspective outlining uncharted waters where the field can advance, and potential exciting future applications of TDP.

Roberts, David S.↗

Comparative differential cuproproteomes of Rhodobacter capsulatus reveal novel copper homeostasis related proteins

Abstract Copper (Cu) is an essential, but toxic, micronutrient for living organisms and cells have developed sophisticated response mechanisms towards both the lack and the excess of Cu in their environments. In this study, we achieved a global view of Cu-responsive changes in the prokaryotic model organism Rhodobacter capsulatus using label-free quantitative differential proteomics. Semi-aerobically grown cells under heterotrophic conditions in minimal medium (∼0.3 μM Cu) were compared with cells supplemented with either 5 μM Cu or with 5 mM of the Cu-chelator bathocuproine sulfonate. Mass spectrometry based bottom-up proteomics of unfractionated cell lysates identified 2430 of the 3632 putative proteins encoded by the genome, producing a robust proteome dataset for R. capsulatus. Use of biological and technical replicates for each growth condition yielded high reproducibility and reliable quantification for 1926 of the identified proteins. Comparison of cells grown under Cu-excess or Cu-depleted conditions to those grown under minimal Cu-sufficient conditions revealed that 75 proteins exhibited statistically significant (p < 0.05) abundance changes, ranging from 2- to 300-fold. A subset of the highly Cu-responsive proteins was orthogonally probed using molecular genetics, validating that several of them were indeed involved in cellular Cu homeostasis.

59 BASIC BIOLOGICAL SCIENCES↗

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Merging Fungal and Bacterial Community Profiles via an Internal Control

Integrated measurements of fungi and bacteria are critical to understand how interactions between these taxa drive key processes in ecosystems ranging from soils to animal guts. High-throughput amplicon sequencing is commonly used to census microbiomes, but the genetic markers targeted for fungi and bacteria (typically ribosomal regions) are domain-specific so profiling must be performed separately, obscuring relationships between these groups. To solve this problem, we developed a spike-in method with an internal control (IC) construct containing primer sites commonly used for bacterial and fungal taxonomic profiling. The internal control offers several advantages: estimation of absolute abundances, estimation of fungal to bacterial ratios (F:B), integration of bacterial and fungal profiles for holistic community analysis, and lower costs compared to other quantitation methods. To validate the IC as a scaling method, we compared IC-derived measures of F:B to measures from quantitative PCR (qPCR) using a commercial mock community (the ZymoBiomic Microbial Community DNA Standard II, containing two fungi and eight bacteria) and complex environmental samples. For both the mock community and the environmental samples, the IC produced F:B values that were statistically consistent with qPCR. Merging the environmental fungal and bacterial profiles based on the IC-derived F:B values revealed new relationships among samples in terms of community similarity. As a result, this IC method is the first spike-in method to employ a single construct for cross-domain amplicon sequencing, offering more reliable measurements.

59 BASIC BIOLOGICAL SCIENCES↗

Quantitative differentiation of benign and misfolded glaucoma-causing myocilin variants on the basis of protein thermal stability

Accurate predictions of the pathogenicity of mutations associated with genetic diseases are key to the success of precision medicine. Inherited missense mutations in the myocilin (MYOC) gene, within its olfactomedin (OLF) domain, constitute the strongest genetic link to primary open-angle glaucoma via a toxic gain of function, and thus MYOC is an attractive precision-medicine target. However, not all mutations in MYOC cause glaucoma, and common variants are expected to be neutral polymorphisms. The Genome Aggregation Database (gnomAD) lists ~100 missense variants documented within OLF, all of which are relatively rare (allele frequency <0.001%) and nearly all are of unknown pathogenicity. To distinguish disease-causing OLF variants from benign OLF variants, we first characterized the most prevalent population-based variants using a suite of cellular and biophysical assays, and identified two variants with features of aggregation-prone familial disease variants. Next, we considered all available biochemical and clinical data to demonstrate that pathogenic and benign variants can be differentiated statistically based on a single metric: the thermal stability of OLF. Our results motivate genotyping MYOC in patients for clinical monitoring of this widespread, painless and irreversible ocular disease.

59 BASIC BIOLOGICAL SCIENCES↗

Knowledge graph-aided Bayesian active learning for top- K genetic interaction discovery

In silico methods for predicting the effects of multi-gene perturbations hold great promise for advancing functional genomics, computational drug discovery, and disease modeling. However, the development of these predictive algorithms for mammalian systems has been hampered by limited datasets and high experimental costs. In this study, we present a Bayesian active learning framework designed to discover pairwise host gene knockdowns that effectively inhibit viral proliferation in an in vitro HIV-1 infection model. Our method leverages a biological knowledge graph as side information and employs a computationally efficient batch diversification approach. We evaluated this framework using a dataset of viral load measurements obtained from multi-day dual-gene depletion experiments, encompassing all possible pairwise knockdowns of over 350 host genes associated with HIV infection. We demonstrate that our framework rapidly identifies the most effective gene knockdown pairs for reducing viral load. Furthermore, we show that incorporating side information enhances performance during the early stages of active learning (low data regime), while our batch diversification strategy significantly boosts performance in later stages (high data regime). This framework is general and can be adapted to explore gene interactions in other contexts, such as synthetic lethality prediction and mapping epistatic effects across quantitative trait loci.

Computational biology and bioinformatics↗

Expandable Li Percolation Network: The Effects of Site Distortion in Cation-Disordered Rock-Salt Cathode Material

Cation-disordered rock-salt (DRX) materials receive intensive attention as a new class of cathode candidates for high-capacity lithium-ion batteries (LIBs). Unlike traditional layered cathode materials, DRX materials have a three-dimensional (3D) percolation network for Li + transportation. The disordered structure poses a grand challenge to a thorough understanding of the percolation network due to its multiscale complexity. In this work, we introduce the large supercell modeling for DRX material Li 1.16 Ti 0.37 Ni 0.37 Nb 0.10 O 2 (LTNNO) via the reverse Monte Carlo (RMC) method combined with neutron total scattering. Here, through a quantitative statistical analysis of the material’s local atomic environment, we experimentally verified the existence of short-range ordering (SRO) and uncovered an element-dependent behavior of transition metal (TM) site distortion. A displacement from the original octahedral site for Ti 4+ cations is pervasive throughout the DRX lattice. Density functional theory (DFT) calculations revealed that site distortions quantified by the centroid offsets could alter the migration barrier for Li + diffusion through the tetrahedral channels, which can expand the previously proposed theoretical percolating network of Li. The estimated accessible Li content is highly consistent with the observed charging capacity. The newly developed characterization method here uncovers the expandable nature of the Li percolation network in DRX materials, which may provide valuable guidelines for the design of superior DRX materials.

25 ENERGY STORAGE↗

An interactive machine learning platform for analyzing multi-particle coincidence data from cold target recoil ion momentum spectroscopy

We present SCULPT (Supervised Clustering and Uncovering Latent Patterns with Training), a comprehensive software platform for analyzing tabulated high-dimensional multi-particle coincidence data from Cold Target Recoil Ion Momentum Spectroscopy (COLTRIMS) experiments. The software addresses critical challenges in modern momentum spectroscopy by integrating advanced machine learning techniques with physics-informed analysis in an interactive web-based environment. SCULPT implements uniform manifold approximation and projection for non-linear dimensionality reduction to reveal correlations in high-dimensional data. We also discuss potential extensions to deep autoencoders for feature learning and genetic programming for automated discovery of physically meaningful observables. A novel adaptive confidence scoring system provides quantitative reliability assessments by evaluating user-selected clustering quality metrics with predefined weights that reflect each metric’s robustness. The platform features configurable molecular profiles for different experimental systems, interactive visualization with selection tools, and comprehensive data filtering capabilities. Utilizing a subset of SCULPT’s capabilities, we analyze photo-double-ionization data measured using the COLTRIMS method for three-body dissociation of the D 2 O molecule, revealing distinct fragmentation channels and their correlations with physics parameters. The software’s modular architecture and web-based implementation make it accessible to the broader atomic and molecular physics community, significantly reducing the time required for complex multi-dimensional analyses. This opens the door to finding and isolating rare events exhibiting non-linear correlations on the fly during experimental measurements, which can help steer exploration and improve the efficiency of experiments.

Artificial neural networks↗