Search NASA⌕ Search

SEARCH · Search NASA

Results for “genetic testing”

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 91 records · Page 5

Intraspecific variability in plant and soil chemical properties in a common garden plantation of the energy crop Populus

Optimizing crops for synergistic soil carbon (C) sequestration can enhance CO 2 removal in food and bioenergy production systems. Yet, in bioenergy systems, we lack an understanding of how intraspecies variation in plant traits correlates with variation in soil biogeochemistry. This knowledge gap is exacerbated by both the heterogeneity and difficulty of measuring belowground traits. Here, we provide initial observations of C and nutrients in soil and root and stem tissues from a common garden field site of diverse, natural variant, Populus trichocarpa genotypes—established for aboveground biomass-to-biofuels research. Our goal was to explore the value of such field sites for evaluating genotype-specific effects on soil C, which ultimately informs the potential for optimizing bioenergy systems for both aboveground productivity and belowground C storage. To do this, we investigated variation in chemical traits at the scale of individual trees and genotypes and we explored correlations among stem, root, and soil samples. We observed substantial variation in soil chemical properties at the scale of individual trees and specific genotypes. While correlations among elements were observed both within and among sample types (soil, stem, root), above-belowground correlations were generally poor. We did not observe genotype-specific patterns in soil C in the top 10 cm, but we did observe genotype associations with soil acid-base chemistry (soil pH and base cations) and bulk density. Finally, a specific phenotype of interest (high vs low lignin) was unrelated to soil biogeochemistry. Our pilot study supports the usefulness of decade-old, genetically-variable, Populus bioenergy field test plots for understanding plant genotype effects on soil properties. Finally, this study contributes to the advancement of sampling methods and baseline data for Populus systems in the Pacific Northwest, USA. Further species- and region-specific efforts will enhance C predictability across scales in bioenergy systems and, ultimately, accelerate the identification of genotypes that optimize yield and carbon storage.

54 ENVIRONMENTAL SCIENCES↗

Data from: "Lab and Field Warming Similarly Advance Germination Date and Limit Germination Rate for High and Low Elevation Provenances of Two Widespread Subalpine Conifers"

This archive contains data used to draw conclusions in “Lab and Field Warming Similarly Advance Germination Date and Limit Germination Rate for High and Low Elevation Provenances of Two Widespread Subalpine Conifers”, by Kueppers et al. 2017. Cone collection and field experiments were conducted on Niwot Ridge, Colorado, with the latter completed within the Alpine Treeline Warming Experiment (ATWE). Laboratory germination experiments were conducted at the University of California Merced. Geospatial files are also provided for additional geographical context.There are four data file formats in this package: comma-separated values (.csv), Microsoft Excel (.xlsx), keyhole markup language (.kml), and ESRI shapefile (.shp). The latter two are geospatial file formats. .csv files can be opened using any simple text editor program such as TextEdit (Mac) and Notepad (Windows), and .xlsx files can be opened using Microsoft Excel or Google Sheets. The .kml file contains coordinate data, and is compatible with Google Earth and Google Maps. ESRI shapefiles (.shp) are compatible with any geospatial software able to read the file type, such as QGIS and ESRI’s ArcGIS Suite. This data user’s guide is available in as a .pdf, and can be opened using software such as Adobe Acrobat Reader.-----------------------------------------------------------------------------------------------------------------------------------Accurately predicting upslope shifts in subalpine tree ranges with warming requires understanding how future forest populations will be affected by climate change, as these are the seed sources for new tree line and alpine populations. Early life history stages are particularly sensitive to climate and are also influenced by genetic variation among populations. We tested the climate sensitivity of germination and initial development for two widely distributed subalpine conifers, using controlled-environment growth chambers with one temperature regime from subalpine forest in the Colorado Rocky Mountains and one 5 °C warmer, and two soil moisture levels. We tracked germination rate and timing, rate of seedling development, and seedling morphology for two seed provenances separated by ~300 m elevation. Warming advanced germination timing and initial seedling development by a total of ~2 weeks, advances comparable to mean differences between provenances. Advances were similar for both provenances and species; however, warming reduced the overall germination rate, as did low soil moisture, only for Picea engelmannii. A three-year field warming and watering experiment planted with the same species and provenances yielded responses qualitatively consistent with the lab trials. Together these experiments indicate that in a warmer, drier climate, P. engelmannii germination, and thus regeneration, could decline, which could lead to declining subalpine forest populations, while Pinus flexilis forest populations could remain robust as a seed source for upslope range shifts.

54 ENVIRONMENTAL SCIENCES↗

Altered ß-cell Prohormone Processing and Secretion in Type 1 Diabetes

Analysis of data from clinical cohorts and more recently from human pancreatic tissue, indicate that defects in prohormone processing are an early and persistent component of type 1 diabetes pathogenesis. In this Perspectives article, we review the current state of knowledge of alterations in islet prohormone expression and processing in type 1 diabetes, and consider the clinical impact of these findings. Lingering questions, including pathologic etiologies and consequences of altered prohormone expression and secretion in type 1 diabetes, and the natural history of circulating prohormone production in health and disease are considered. Finally, key steps required to move forward in this area are outlined, including longitudinal testing of relevant clinical populations, studies that probe the genetics of altered prohormone processing, the need for combined functional and histologic testing of human pancreatic tissues, continued interrogation of the intersection between prohormone processing and autoimmunity, and optimal assays or approaches for analysis. Successful resolution of these questions may offer the potential to use altered prohormone processing as a pathogenic anchor that can be used to cluster different endotypes enveloping a large set of variables, and inform therapeutic strategies aimed at personalized intervention in the natural history of type 1 diabetes.

Islets, Insulin, Proinsulin, Islet Amyloid Polypep↗

Nucleic Acid-Based Sensing Techniques for Diagnostics and Surveillance of Influenza

Influenza virus poses a threat to global health by causing seasonal outbreaks as well as three pandemics in the 20th century. In humans, disease is primarily caused by influenza A and B viruses, while influenza C virus causes mild disease mostly in children. Influenza D is an emerging virus found in cattle and pigs. To mitigate the morbidity and mortality associated with influenza, rapid and accurate diagnostic tests need to be deployed. However, the high genetic diversity displayed by influenza viruses presents a challenge to the development of a robust diagnostic test. Nucleic acid-based tests are more accurate than rapid antigen tests for influenza and are therefore better candidates to be used in both diagnostic and surveillance applications. We review various nucleic acid-based techniques that have been applied towards the detection of influenza viruses in order to evaluate their utility as both diagnostic and surveillance tools. We discuss both traditional as well as novel methods to detect influenza viruses by covering techniques that require nucleic acid amplification or direct detection of viral RNA as well as comparing advantages and limitations for each method. There has been substantial progress in the development of nucleic acid-based sensing techniques for the detection of influenza virus. However, there is still an urgent need for a rapid and reliable influenza diagnostic test that can be used at point-of-care in order to enhance responsiveness to both seasonal and pandemic influenza outbreaks.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

Ecological connectivity and habitat loss shape patterns of genetic diversity in a threatened salamander

Context The maintenance of genetic diversity is essential for preserving adaptive potential in populations, yet it is increasingly threatened by landscape alteration. The field of landscape genetics offers a framework for assessing how patch-level landscape conditions, modeled at multiple scales, influence genetic diversity. Objectives We sought to assess how local environmental features and connectivity influence genetic diversity across 74 four-toed salamander (Hemidactylium scutatum) breeding wetlands in the southeastern United States. Methods Using next-generation sequencing data and hierarchical Bayesian models, we examined genome-wide heterozygosity in relation to local landscape features and ecological connectivity. We also assessed the scale of effect of landscape features and tested for temporal lag effects. Results Genetic diversity was lower in wetlands with higher levels of historic deforestation and lower connectivity. An interaction between deforestation and connectivity indicated that deforestation had stronger negative effects in isolated wetlands but weaker effects in well-connected wetlands. Accounting for scale of effect and temporal lags was critical for detecting these relationships. Conclusions Our analyses highlight the importance of assessing the spatial scale (scale of effect) and temporal lag of landscape features to detect key drivers of genetic diversity. In line with population genetic theory, our results indicate that the genetic consequences of habitat loss do not affect populations uniformly and are most severe in isolated populations where gene flow cannot buffer against loss of diversity. Altogether, we highlight the importance of considering the interaction of habitat loss and connectivity in conservation genetic management.

Hemidactylium scutatum↗

Evolutionary Design of a Phased Array Antenna Element

We present an evolved S-band phased array antenna element design that meets the requirements of NASA's TDRS-C communications satellite scheduled for launch early next decade. The original specification called for two types of elements, one for receive only and one for transmit/receive. We were able to evolve a single element design that meets both specifications thereby simplifying the antenna and reducing testing and integration costs. The highest performance antenna found using a genetic algorithm and stochastic hill-climbing has been fabricated and tested. Laboratory results are largely consistent with simulation. Researchers have been investigating evolutionary antenna design and optimization since the early 1990s, and the field has grown in recent years its computer speed has increased and electromagnetic simulators have improved. Many antenna types have been investigated, including wire antennas, antenna arrays and quadrifilar helical antennas. In particular, our laboratory evolved a wire antenna design for NASA's Space Technology 5 (ST5) spacecraft. This antenna has been fabricated, tested, and is scheduled for launch on the three spacecraft in 2006.

Globus, Al↗

DropSynth 2.0: high-fidelity multiplexed gene synthesis in emulsions

Abstract Multiplexed assays allow functional testing of large synthetic libraries of genetic elements, but are limited by the designability, length, fidelity and scale of the input DNA. Here, we improve DropSynth, a low-cost, multiplexed method that builds gene libraries by compartmentalizing and assembling microarray-derived oligonucleotides in vortexed emulsions. By optimizing enzyme choice, adding enzymatic error correction and increasing scale, we show that DropSynth can build thousands of gene-length fragments at >20% fidelity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient population representation with more genetic markers increases performance of a steelhead (Oncorhynchus mykiss) genetic stock identification baseline

Abstract Genetic stock identification (GSI) is an important fisheries management tool to identify the origin of fish harvested in mixed stock fisheries. Periodic updates of genetic baselines can improve performance via the addition of unsampled or under‐sampled populations and the inclusion of more informative markers. We used a combination of baselines to evaluate how population representation, marker number, and marker type affected the performance and accuracy of genetic stock assignments (self‐assignment, bias, and holdout group tests) for steelhead ( Oncorhynchus mykiss ) in the Snake River basin. First, we compared the performance of an existing genetic baseline with a newly developed one which had a reduced number of individuals from more populations using the same set of markers. Self‐assignment rates were significantly higher ( p < 0.001; +5.4%) for the older, larger baseline, bias did not differ significantly between the two, but there was a significant improvement in performance for the new baseline in holdout results ( p < 0.001; mean increase of 25.0%). Second, we compared the performance of the new baseline with increased numbers of genetic markers (~2x increase of single‐nucleotide polymorphisms; SNPs) for the same set of baseline individuals. In this comparison, results produced significantly higher rates of self‐assignment ( p < 0.001; +9.7%) but neither bias nor leave‐one‐out were significantly affected. Third, we compared 334 SNPs versus opportunistically discovered microhaplotypes from the same amplicons for the new baseline, and showed the latter produced significantly higher rates of self‐assignment ( p < 0.01; +2.6%), similar bias, but slightly lower holdout performance (−0.1%). Combined, we show the performance of genetic baselines can be improved via representative and efficient sampling, that increased marker number consistently improved performance over the original baseline, and that opportunistic discovery of microhaplotypes can lead to small improvements in GSI performance.

Hargrove, John S.↗

Chymotrypsin-like Elastase-1 Mediates Progressive Emphysema in Alpha-1 Antitrypsin Deficiency

Alpha-1 antitrypsin (AAT) deficiency is a rare disease affecting approximately 1 in 2000 White individuals with approximately 10% of these individuals developing AATD lung disease. This lung disease is marked by progressive alveolar loss despite the withdrawal of triggering agents such as cigarette smoke. Although the PiZZ genotype is the most common mutation, there are over 100 described variants making true population estimates difficult and AAT deficiency is typically diagnosed by reduced levels of AAT in the blood. Typically developing in the fourth and fifth decades of life, AAT-deficient lung disease is marked by progressive emphysema and is the fourth leading indication for lung transplantation. AAT augmentation therapy does not prevent disease progression making the development of new therapeutic approaches critical. Chymotrypsin-like elastase 1 (CELA1) is a serine protease synthesized and secreted by alveolar type 2 cells with a physiologic role in reducing postnatal lung elastance. At a molecular level, CELA1 binds and cleaves non-crosslinked, hydrophobic domains of tropoelastin, and its binding to lung elastin fibers is increased with strain—similar to other pancreatic elastases. CELA1 is neutralized by covalent binding with AAT, and Cela1 -/- mice were completely protected from emphysema in an antisense oligonucleotide model of AAT-deficient emphysema. This model, however, did not include any injury apart from administration of the antisense oligonucleotide with levels of emphysema exceeding that seen in mice with genetic ablation of 5 Serpina1 paralogues and subjected to tracheal lipopolysacharide or cigarette smoke. Here, we use this murine genetic model of AAT deficiency to test the role of the CELA1 gene in AAT-deficient emphysema using multiple models to show that CELA1 has a role in progressive airspace enlargement in AAT-deficiency independent of inflammation.

60 APPLIED LIFE SCIENCES↗

Detecting Arsenic Contamination Using Satellite Imagery and Machine Learning

Arsenic, a potent carcinogen and neurotoxin, affects over 200 million people globally. Current detection methods are laborious, expensive, and unscalable, being difficult to implement in developing regions and during crises such as COVID-19. This study attempts to determine if a relationship exists between soil’s hyperspectral data and arsenic concentration using NASA’s Hyperion satellite. It is the first arsenic study to use satellite-based hyperspectral data and apply a classification approach. Four regression machine learning models are tested to determine this correlation in soil with bare land cover. Raw data are converted to reflectance, problematic atmospheric influences are removed, characteristic wavelengths are selected, and four noise reduction algorithms are tested. The combination of data augmentation, Genetic Algorithm, Second Derivative Transformation, and Random Forest regression (R 2 =0.840 and normalized root mean squared error (re-scaled to [0,1]) = 0.122) shows strong correlation, performing better than past models despite using noisier satellite data (versus lab-processed samples). Three binary classification machine learning models are then applied to identify high-risk shrub-covered regions in ten U.S. states, achieving strong accuracy (=0.693) and F1-score (=0.728). Overall, these results suggest that such a methodology is practical and can provide a sustainable alternative to arsenic contamination detection.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Genetic interactions between polycystin-1 and Wwtr1 in osteoblasts define a novel mechanosensing mechanism regulating bone formation in mice

Molecular mechanisms transducing physical forces in the bone microenvironment to regulate bone mass are poorly understood. Here, we used mouse genetics, mechanical loading, and pharmacological approaches to test the possibility that polycystin-1 and Wwtr1 have interdependent mechanosensing functions in osteoblasts. We created and compared the skeletal phenotypes of control Pkd1 flox/+ ;Wwtr1 flox/+ , Pkd1 Oc-cKO , Wwtr1 Oc-cKO , and Pkd1/Wwtr1 Oc-cKO mice to investigate genetic interactions. Consistent with an interaction between polycystins and Wwtr1 in bone in vivo, Pkd1/Wwtr1 Oc-cKO mice exhibited greater reductions of BMD and periosteal MAR than either Wwtr1 Oc-cKO or Pkd1 Oc-cKO mice. Micro-CT 3D image analysis indicated that the reduction in bone mass was due to greater loss in both trabecular bone volume and cortical bone thickness in Pkd1/Wwtr1 Oc-cKO mice compared to either Pkd1 Oc-cKO or Wwtr1 Oc-cKO mice. Pkd1/Wwtr1 Oc-cKO mice also displayed additive reductions in mechanosensing and osteogenic gene expression profiles in bone compared to Pkd1 Oc-cKO or Wwtr1 Oc-cKO mice. Moreover, we found that Pkd1/Wwtr1 Oc-cKO mice exhibited impaired responses to tibia mechanical loading in vivo and attenuation of load-induced mechanosensing gene expression compared to control mice. Finally, control mice treated with a small molecule mechanomimetic, MS2 that activates the polycystin complex resulted in marked increases in femoral BMD and periosteal MAR compared to vehicle control. In contrast, Pkd1/Wwtr1 Oc-cKO mice were resistant to the anabolic effects of MS2. These findings suggest that PC1 and Wwtr1 form an anabolic mechanotransduction signaling complex that mediates mechanical loading responses and serves as a potential novel therapeutic target for treating osteoporosis.

60 APPLIED LIFE SCIENCES↗

BETO 2021 Peer Review - Cyanobacteria Photosynthetic Energy Platform (1.3.4.301)

Cyanobacteria Photosynthetic Energy Platform (WBS 1.3.4.301) The primary challenge in algal technology is low energy conversion efficiency, currently up to a few percentages from photon energy to biomass energy. Most incoming solar energy is lost in the processes of photosynthesis, through mechanisms that are not fully understood. There is large room for improvement. This project develops cyanobacteria genetic tools to improve photosynthetic efficiency and biomass productivity through manipulation of energy regulation mechanisms. Increasing algal biomass productivity by 20% translates to $2 reduction of MFSP. Genetic manipulation is an essential tool in understanding and improving algal productivity. While genetic manipulation is challenging for most algae, model cyanobacteria, with in-house mutant library and genetic engineering toolboxes, allow for rapid hypothesis testing and transfer of lessons to other cyanobacteria and eukaryotic algae. We developed Energome (energy-ome) concept to guide novel engineering strategies and unleash unused potential in photosynthesis. Energome consists of cellular energy management mechanisms. For example, we found that cyanobacteria use futile cycles around glycogen and sucrose to dissipate ATP. By genetic modification, the model cyanobacterium Synechocystis 6803 showed higher energy levels, increased photosynthetic capacity, and 20% improvement in biomass yield in simulated outdoor light conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

National Prevalence of Caprine Prion Protein Genetic Variability at Codons 146, 211, and 222 in Goat Herds in the United States

Scrapie is a neurodegenerative disease that impacts sheep and goats, characterized by gradual and progressive changes in neurological function. Recent research shows that the scrapie incubation period is significantly influenced by specific variations in amino acids within the prion protein gene (PRNP). The objective of this study was to estimate the national prevalence of caprine PRNP genetic variability at codons 146, 211, and 222 in goat populations across the United States. A total of 3052 blood, ear tissue, and brain tissue samples were collected from goats from 50 states. The participating states were categorized into four Veterinary Service (VS) district regions. The samples underwent DNA extraction, and the PRNP variants corresponding to codons 146, 211, and 222 were amplified and sequenced. The analysis of PRNP variants, when compared to the PRNP reference sequence, revealed seven alleles in twelve genotypes. The homozygous 146NN, 211RR, and 222QQ alleles, which have been linked to an increased risk of scrapie, were found to be the most prevalent among all the goats. The heterozygous 222QK, 211RQ, 146SD, 146ND, and 146NS alleles and the homozygous 222KK, 146SS, and 146DD alleles, known to be associated with reduced scrapie susceptibility and a prolonged incubation period after experimental challenge, were found in 1.098% (222QK), 2.33% (211RQ), 0.58% (146SD), 3.13% (146ND), 20.68% (146NS), 0.005% (222KK), 3.31% (146SS), and 0.67% (146DD) of goats, respectively. The 222QK allele was found most frequently in goats tested from the east (VS District 1, 1.59%) and southwest (VS District 4, 1.08%) regions, whereas the 211RQ allele was found most often in goats tested from the Midwest (VS District 2, 8.03%) and east (VS District 1, 6.53%) regions. The 146NS allele was found most frequently in goats tested from the northwest (VS District 3, 29.02%) and southwest (VS District 4, 20.69%) regions. Our results showed that the prevalence of less susceptible genotypes at PRNP codon 146 may be sufficient to use genetic susceptibility testing in some herds. This may reduce the number of goats removed as part of a herd clean-up plan and may promote the selective breeding goats for less susceptible alleles in high-risk herds at the national level.

59 BASIC BIOLOGICAL SCIENCES↗

Drought increases microbial allocation to stress tolerance but with few tradeoffs among community-level traits

Climate change will increase soil drying, altering microbial communities via increasing water stress and decreasing resource availability. The responses of these microbial communities to changing environments could be governed by physiological tradeoffs between high yield, resource acquisition, and stress tolerance (YAS framework). We leveraged a unique field experiment that manipulates both drought and carbon availability across two years and three land uses, and we measured both physiological (with bioassays) and genetic (with metagenomics) microbial traits at the community level to test the following hypotheses: 1. Drought increases microbial allocation to stress tolerance functions, at both physiological and genetic levels. 2. Because microbes are resource-limited under drought, increased carbon will enable greater expression of stress tolerance. 3. All three key life history traits described in the YAS framework will trade off. Drought did increase microbial physiological investment in stress tolerance (measured via trehalose production), but we saw few other changes in microbial communities under drought. Adding carbon to plots increased resource acquisition (measured via enzyme activity and resource acquisition gene abundance) and stress tolerance (trehalose assay), but did so in both drought and average rainfall environments. Here, we found little evidence of trait tradeoffs, as a negative correlation between rRNA copy number and resource acquisition gene abundance was the only significant negative correlation between traits that we found that was consistent across years (for physiological or genetic traits). In summary, we found C addition, and to a lesser extent, drought, altered microbial community function and functional genes. However, resources did not alter drought response in a way that was consistent with theory of life history tradeoffs.

59 BASIC BIOLOGICAL SCIENCES↗

Association mapping across a multitude of traits collected in diverse environments in maize

Abstract Classical genetic studies have identified many cases of pleiotropy where mutations in individual genes alter many different phenotypes. Quantitative genetic studies of natural genetic variants frequently examine one or a few traits, limiting their potential to identify pleiotropic effects of natural genetic variants. Widely adopted community association panels have been employed by plant genetics communities to study the genetic basis of naturally occurring phenotypic variation in a wide range of traits. High-density genetic marker data—18M markers—from 2 partially overlapping maize association panels comprising 1,014 unique genotypes grown in field trials across at least 7 US states and scored for 162 distinct trait data sets enabled the identification of of 2,154 suggestive marker-trait associations and 697 confident associations in the maize genome using a resampling-based genome-wide association strategy. The precision of individual marker-trait associations was estimated to be 3 genes based on a reference set of genes with known phenotypes. Examples were observed of both genetic loci associated with variation in diverse traits (e.g., above-ground and below-ground traits), as well as individual loci associated with the same or similar traits across diverse environments. Many significant signals are located near genes whose functions were previously entirely unknown or estimated purely via functional data on homologs. This study demonstrates the potential of mining community association panel data using new higher-density genetic marker sets combined with resampling-based genome-wide association tests to develop testable hypotheses about gene functions, identify potential pleiotropic effects of natural genetic variants, and study genotype-by-environment interaction.

59 BASIC BIOLOGICAL SCIENCES↗

Learning epistatic polygenic phenotypes with Boolean interactions

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1 , a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

59 BASIC BIOLOGICAL SCIENCES↗