Search NASA⌕ Search

SEARCH · Search NASA

Results for “disequilibrium”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Temperature dependence of uranium and thorium partitioning in igneous zircons

Zircon is a key mineral in geochronology because of its chemical and physical durability and tendency to incorporate radioactive trace elements such as U and Th. Quantifying the partitioning of the actinide elements is critical to constrain initial non-secular equilibrium amounts of 230 Th in zircon. An excess or deficit of 206 Pb will be produced from such an initial excess/deficit of 230 Th from the secular equilibrium condition, which influences the calculated 206 Pb/ 238 U age (Schärer, 1984, Mattinson, 1973). However, there is no standard way to calculate Th/U partitioning ratios when applying age corrections to young igneous zircon, making uncertainties hard to estimate. To better understand U and Th partitioning of zircon/melt that is the cornerstone of the secular disequilibrium correction, zircon was synthesized in one-atmosphere experiments using basaltic andesite, andesite, and rhyolite starting materials, doped with Zr, U, and Th. Different experimental temperatures and oxygen fugacity conditions (ΔQFM−4 to ΔQFM+4) were explored to examine their effects on U and Th partitioning. In addition, we specifically quantify the effects of sector zoning, fractional crystallization, and melt composition on U and Th partitioning. By combining experimental and natural zircon data, we find that temperature has the primary control on the partitioning of U and Th in the zircon and calibrate an inverse relationship between these partition coefficients and crystallization temperatures. The calibrated equation can be applied to the 230 Th correction for an improvement in the accuracy of Th-corrected 206 Pb/ 238 U dates when the zircon crystallization temperature is known.

206Pb/238U dating↗
DOE OSTI2022

Manganese Oxidation States in Volcanic Soils across Annual Rainfall Gradients

Manganese (Mn) exists as Mn(II), Mn(III), or Mn(IV) in soils, and the Mn oxidation state controls the roles of Mn in numerous environmental processes. However, the variations of Mn oxidation states with climate remain unknown. Here we determined the Mn oxidation states in highly weathered bulk volcanic soils (primary minerals free) across two rainfall gradients covering mean annual precipitation (MAP) of 0.25–5 m in the Hawaiian Islands. With increasing MAP, the soil redox conditions generally shifted from oxic to suboxic and to anoxic despite fluctuating at each site; concurrently, the proportions of Mn(IV) and Mn(II) decreased and increased, respectively. Mn(III) was low at both low and high MAP, but accumulated substantially, up to 80% of total Mn, in soils with prevalent suboxic conditions at intermediate MAP. Mn(III) was likely hosted in Mn(III,IV) and iron(III) oxides or complexed with organic matter, and its distribution among these hosts varied with soil redox potentials and soil pH. Soil redox conditions and rainfall-driven leaching jointly controlled exchangeable Mn(II) in soils, with its concentration peaking at intermediate MAP. The Mn redox chemistry was at disequilibrium, with the oxidation states correlating with long-term average soil redox potentials better than with soil pH. The soil redox conditions likely fluctuated between oxic and anoxic conditions more frequently at intermediate than at low and high MAP, creating biogeochemical hot spots where Mn, Fe, and other redox-sensitive elements may be actively cycled.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Matrix Approach to Land Carbon Cycle Modeling

Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchanges among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Coupled Transport, Reactivity, and Mechanics in Fractured Shale Caprocks

Abstract Shales are low‐permeability caprocks that confine fluid, such as CO 2 , nuclear waste, and hydrogen, in storage formations. Stress‐induced fractures in shale caprocks provide pathways for fluid to leak and potentially contaminate fresh water aquifers. Fractured shales are also increasingly considered as resources for CO 2 sequestration, enhanced geothermal, and unconventional energy recovery. Injecting reactive fluids into shales introduces chemical disequilibrium, causing an onset of a series of dissolution, precipitation, and fines mobilization mechanisms. The reactions have rapid kinetics and significant impact on porosity and permeability; consequently, flow and storage properties of caprocks. While previous research has explored the separate effects of these reactions, this study aims to uncover their simultaneous occurrence and collective influence. This study unveils these highly coupled transport and reactivity mechanisms by tracking and visualizing the reaction‐induced alterations in the matrix, microcracks, and fractures of shales over time. We conducted brine injection experiments sequentially at pH 4 and 2 in a naturally fractured Wolfcamp shale sample while simultaneously imaging the dynamic processes using X‐ray computed tomography (CT). CT images are validated by finer resolution images obtained using micro‐CT and scanning electron microscopy. We also tracked the sample permeability and fluid chemistry using brine permeability and inductively coupled plasma mass spectrometry, respectively. Findings show that fluid primarily flowed through fractures, dissolving reactive minerals and mobilizing fines on fracture surfaces. Dissolution of fracture asperities under confining stress resulted in the closing of fractures. Clogging in narrow fracture pathways, caused by fines accumulation, diverted fluid flow into matrix pores.

58 GEOSCIENCES↗
DOE OSTI2024

A high internal heat flux and large core in a warm Neptune exoplanet

Interactions between exoplanetary atmospheres and internal properties have long been proposed to be drivers of the inflation mechanisms of gaseous planets and apparent atmospheric chemical disequilibrium conditions. However, transmission spectra of exoplanets have been limited in their ability to observationally confirm these theories owing to the limited wavelength coverage of the Hubble Space Telescope (HST) and inferences of single molecules, mostly H 2 O. Here, in this work, we present the panchromatic transmission spectrum of the approximately 750 K, low-density, Neptune-sized exoplanet WASP-107b using a combination of HST Wide Field Camera 3 (WFC3) and JWST Near-Infrared Camera (NIRCam) and Mid-Infrared Instrument (MIRI). From this spectrum, we detect spectroscopic features resulting from H 2 O (21σ), CH 4 (5σ), CO (7σ), CO 2 (29σ), SO 2 (9σ) and NH 3 (6σ). The presence of these molecules enables constraints on the atmospheric metal enrichment (M/H is 10–18× solar), vertical mixing strength (log 10 K zz = 8.4–9.0 cm 2 s -1 ) and internal temperature (>345 K). The high internal temperature is suggestive of tidally driven inflation4 acting on a Neptune-like internal structure, which can naturally explain the large radius and low density of the planet. These findings suggest that eccentricity-driven tidal heating is a critical process governing atmospheric chemistry and interior-structure inferences for most of the cool (<1,000 K) super-Earth-to-Saturn-mass exoplanet population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

Population structure limits the use of genomic data for predicting phenotypes and managing genetic resources in forest trees

There is overwhelming evidence that forest trees are locally adapted to climate. Thus, genecological models based on population phenotypes have been used to measure local adaptation, infer genetic maladaptation to climate, and guide assisted migration. However, instead of phenotypes, there is increasing interest in using genomic data for gene resource management. We used whole-genome resequencing and common-garden experiments to understand the genetic architecture of adaptive traits in black cottonwood. We studied the potential of using genome-wide association studies (GWAS) and genomic prediction to detect causal loci, identify climate-adapted phenotypes, and inform gene resource management. We analyzed population structure by partitioning phenotypic and genomic (single-nucleotide polymorphism) variation among 840 genotypes collected from 91 stands along 16 rivers. Most phenotypic variation (60 to 81%) occurred among populations and was strongly associated with climate. Population phenotypes were predicted well using genomic data (e.g., predictive abilityr> 0.9) but almost as well using climate or geography (r> 0.8). In contrast, genomic prediction within populations was poor (r< 0.2). We identified many GWAS associations among populations, but most appeared to be spurious based on pooled within-population analyses. Hierarchical partitioning of linkage disequilibrium and haplotype sharing suggested that within-population genomic prediction and GWAS were poor because allele frequencies of causal loci and linked markers differed among populations. Given the urgent need to conserve natural populations and ecosystems, our results suggest that climate variables alone can be used to predict population phenotypes, delineate seed zones and deployment zones, and guide assisted migration.

Science & Technology - Other Topics↗
DOE OSTI2025

MINE: maximally informative next experiment—toward a new GWAS experimental design and methodology

Abstract The computational methodology of Genome Wide Association Studies (GWAS) currently has several limitations: (i) the number of observations (rows) on a quantitative trait tends to be smaller than the number of single nucleotide polymorphisms (SNPs) (columns) in the design matrix; (ii) each SNP is usually modeled separately, failing to acknowledge interaction between each other (ie epistasis); (iii) there is implicit linkage disequilibrium (LD) between neighboring SNPs due to their linkage. To overcome these issues, we developed a tool that uses ensemble methods to fit mixed linear models to GWAS data, and these ensemble methods include the development of a new experimental design approach in GWAS, which uses the resultant models and data to select the next informative experiment over time. This new adaptive and staged approach for GWAS experimental design was developed and tested in a 3 yr adaptive model-guided discovery experiment against a fixed classical design. In Sorghum bicolor a total of 79, 86, and 78 accessions were tested in years 1, 2, and 3, respectively out of 343 accessions available in the Bioenergy Association Panel (BAP) each identified for 232,303 SNPs, 1 every 2–3 kb in the genomes. We demonstrated the feasibility of MINE enacted with 8 people in the field per year over 3 yr vs in 1 large classical design enacted with 20 people in 1 yr. The MINE results for chromosomal regions identified controlling dry weight were confirmed against results from previous sorghum GWAS experiments and 1 large classical design for the BAP panel.

Genetics & Heredity↗
DOE OSTI2021

Genomic Signatures of a Major Adaptive Event in the Pathogenic Fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A genomic region associated with iteroparous spawning phenology is linked with age‐at‐maturity in female steelhead trout

Abstract Age‐at‐maturity and iteroparity are two life history variations of steelhead trout ( Oncorhynchus mykiss ) that are believed to increase population resilience and stability. While repeat‐spawning individuals are thought to have historically made up a substantial portion of the reproductive population in the Columbia River and the majority of females still attempt outmigration as kelts, return rates of repeat‐spawner are low throughout the basin and below 1% for the furthest migrating stocks. Notably, outmigrating adults exhibit variation in rematuration phenology, displaying either “consecutive” (reproduce immediately the following season) or “skip” (delay spawning for future seasons) spawning patterns. Here, we use low coverage whole genome sequencing of consecutive versus skip spawning female Columbia River steelhead from two populations to test for genomic differences between these two iteroparous phenotypes. We identified genomic regions on several chromosomes which were associated with the phenology of iteroparity, including a region on chromosome 25 containing two genes, estradiol receptor beta ( ER β) and glycoprotein hormone beta‐5 ( GPHB5 ), which, in mammals, are estrogen‐sensitive and expressed in reproductive tissues. Allele frequencies in this ER β/ GPHB5 region differed among female steelhead of different age at maturity, but not males. These genes also shared an island of linkage disequilibrium with the SIX6 gene, 600Kbp away on the same chromosome, a region of known association with age‐at‐maturity. These observations contribute to growing evidence that age‐at‐maturity and the phenology of iteroparity are determined by overlapping physiological processes and genetic pathways.

Willis, Stuart↗
DOE OSTI2024

Pressure Transducer Measurement Variability in Deep Wells Screened Across the Water Table

Abstract Automated water level measurements collected using vented pressure transducers in deep wells screened across the water table may exhibit a greater response to barometric pressure changes than the true water level. The cause was hypothesized to be disequilibrium in barometric pressure between the wellbores and land surface due to air exchange with the deep vadose zone. In this study, vented and nonvented pressure transducers were installed and operated simultaneously in two deep wells screened across the water table. A vent tube open to the atmosphere at land surface allowed for barometric compensation of the vented transducers. Two nonvented transducers were installed in each well, one submerged in the water and one above the water surface. The difference in readings allowed for barometric compensation. Manual measurements were also collected. It was confirmed that measurements from the vented transducers exhibited greater variability in response to barometric pressure changes than the nonvented transducers and manual measurements. Comparison of the downhole barometric pressure measurements to values from a nearby meteorology station showed the response in the wells to changes in barometric pressure was time‐lagged and attenuated. Thus, the reference pressure from land surface supplied to the vented transducers was not representative of the air pressure within the wells. This caused fluctuations of the transducer readings in response to barometric pressure changes to be greater than the true water level change. This issue can be resolved by the use of nonvented pressure transducers.

McDonald, John P.↗
DOE OSTI2024

Providing biological context for GWAS results using eQTL regulatory and co‐expression networks in Populus

Summary Our study utilized genome‐wide association studies (GWAS) to link nucleotide variants to traits in Populus trichocarpa , a species with rapid linkage disequilibrium decay. The aim was to overcome the challenge of interpreting statistical associations at individual loci without sufficient biological context, which often leads to reliance solely on gene annotations from unrelated model organisms. We employed an integrative approach that included GWAS targeting multiple traits using three individual techniques for lignocellulose phenotyping, expression quantitative trait loci (eQTL) analysis to construct transcriptional regulatory networks around each candidate locus and co‐expression analysis to provide biological context for these networks, using lignocellulose biosynthesis in Populus trichocarpa as a case study. The research identified three candidate genes potentially involved in lignocellulose formation, including one previously recognized gene (Potri.005G116800/VND1, a critical regulator of secondary cell wall formation) and two genes (Potri.012G130000/AtSAP9 and Potri.004G202900/BIC1) with newly identified putative roles in lignocellulose biosynthesis. Our integrative approach offers a framework for providing biological context to loci associated with trait variation, facilitating the discovery of new genes and regulatory networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Nonphotochemical quenching kinetics GWAS in sorghum identifies genes that may play conserved roles in maize and Arabidopsis thaliana photoprotection

SUMMARY Photosynthetic organisms must cope with rapid fluctuations in light intensity. Nonphotochemical quenching (NPQ) enables the dissipation of excess light energy as heat under high light conditions, whereas its relaxation under low light maximizes photosynthetic productivity. We quantified variation in NPQ kinetics across a large sorghum ( Sorghum bicolor ) association panel in four environments, uncovering significant genetic control for NPQ. A genome‐wide association study (GWAS) confidently identified three unique regions in the sorghum genome associated with NPQ and suggestive associations in an additional 61 regions. We detected strong signals from the sorghum ortholog of Arabidopsis thaliana Suppressor Of Variegation 3 ( SVR3 ) involved in plastid–nucleus signaling. By integrating GWAS results for NPQ across maize ( Zea mays ) and sorghum‐association panels, we identified a second gene, Non‐yellowing 1 ( NYE1 ), originally studied by Gregor Mendel in pea ( Pisum sativum ) and involved in the degradation of photosynthetic pigments in light‐harvesting complexes. Analysis of nye1 insertion alleles in A. thaliana confirmed the effect of this gene on NPQ kinetics in eudicots. We extended our comparative genomics GWAS framework across the entire maize and sorghum genomes, identifying four additional loci involved in NPQ kinetics. These results provide a baseline for increasing the accuracy and speed of candidate gene identification for GWAS in species with high linkage disequilibrium.

Plant Sciences↗
DOE OSTI2021

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Compositionally and density stratified igneous terrain in Jezero crater, Mars

Before Perseverance, Jezero crater’s floor was variably hypothesized to have a lacustrine, lava, volcanic airfall, or aeolian origin. SuperCam observations in the first 286 Mars days on Mars revealed a volcanic and intrusive terrain with compositional and density stratification. The dominant lithology along the traverse is basaltic, with plagioclase enrichment in stratigraphically higher locations. Stratigraphically lower, layered rocks are richer in normative pyroxene. The lowest observed unit has the highest inferred density and is olivine-rich with coarse (1.5 millimeters) euhedral, relatively unweathered grains, suggesting a cumulate origin. This is the first martian cumulate and shows similarities to martian meteorites, which also express olivine disequilibrium. Alteration materials including carbonates, sulfates, perchlorates, hydrated silicates, and iron oxides are pervasive but low in abundance, suggesting relatively brief lacustrine conditions. Orbital observations link the Jezero floor lithology to the broader Nili-Syrtis region, suggesting that density-driven compositional stratification is a regional characteristic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Frequency of Sex: Population Genomics Reveals Differences in Recombination and Population Structure of the Aflatoxin-Producing Fungus Aspergillus flavus

The apparent rarity of sex in many fungal species has raised questions about how much sex is needed to purge deleterious mutations and how differences in frequency of sex impact fungal evolution. We sought to determine how differences in the extent of recombination between populations of Aspergillus flavus impact the evolution of genes associated with the synthesis of aflatoxin, a notoriously potent carcinogen. We sequenced the genomes of, and quantified aflatoxin production in, 94 isolates of A. flavus sampled from seven states in eastern and central latitudinal transects of the United States. The overall population is subdivided into three genetically differentiated populations (A, B, and C) that differ greatly in their extent of recombination, diversity, and aflatoxin-producing ability. Estimates of the number of recombination events and linkage disequilibrium decay suggest relatively frequent sex only in population A. Population B is sympatric with population A but produces significantly less aflatoxin and is the only population where the inability of nonaflatoxigenic isolates to produce aflatoxin was explained by multiple gene deletions. Population expansion evident in population B suggests a recent introduction or range expansion. Population C is largely nonaflatoxigenic and restricted mainly to northern sampling locations through restricted migration and/or selection. Despite differences in the number and type of mutations in the aflatoxin gene cluster, codon optimization and site frequency differences in synonymous and nonsynonymous mutations suggest that low levels of recombination in some A. flavus populations are sufficient to purge deleterious mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A Review of Variables Impacting the Indoor Inhalation Radon Equilibrium Factor (FEQ)

With radon and its daughter products estimated as the second leading cause of lung cancer in the United States, it is imperative to understand their relative equilibrium inside commercial, community, and residential dwellings. The radon indoor inhalation fractional equilibrium factor (F eq ) quantifies the disequilibrium between radon and its progeny in indoor air, and recent advances have shown how air exchange rates (ACH) influence F eq . These numerically derived ACH-dependent F eq values are incorporated into the U.S. EPA's Radon Vapor Intrusion Screening Level (RVISL) calculator, which assists risk assessors in evaluating radon exposure. To advance the risk assessment science of actinon (Rn-219), thoron (Rn-220), and radon (Rn-222), the impact of variables such as indoor aerosol concentration and composition, outdoor air quality, household-specific characteristics, and environmental/meteorological conditions on the F eq must be examined. The primary objective of this research is to analyze these additional variables to determine the usefulness of incorporating such adjustment factors into the RVISL calculator and to identify areas of future research. Studies regarding the influence of these parameters are presented along with recommendations regarding the adjustment of the numerically derived F eq value. For example, elevated indoor aerosol concentrations, such as those originating from outdoor PM 2.5 or cigarette smoke, increase the abundance of accumulation- mode particles indoors, which in turn raises F eq values by facilitating the attachment of radon progeny to these aerosols. Smoking increases both the bronchial dose and F eq , while regions with high smog levels demonstrate the impact of regional air quality on F eq . In contrast, air cleaning systems and purifiers have been shown to reduce the F eq , suggesting that these systems could help mitigate radon exposure. Additionally, higher F eq values are typically observed during winter when ventilation rates are lower. This paper presents adjustment factors that may be applied to the RVISL F eq , emphasizing the need for further research to refine these variables and ensure accurate risk assessments in diverse environments. Applying these adjustment factors will minimize calculator over- and underestimations, providing a more accurate representation of the real-world risk associated with radon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Genome-Wide Association Study Reveals Complex Genetic Architecture of Cadmium and Mercury Accumulation and Tolerance Traits in Medicago truncatula

Heavy metals are an increasing problem due to contamination from human sources that and can enter the food chain by being taken up by plants. Understanding the genetic basis of accumulation and tolerance in plants is important for reducing the uptake of toxic metals in crops and crop relatives, as well as for removing heavy metals from soils by means of phytoremediation. Following exposure of Medicago truncatula seedlings to cadmium (Cd) and mercury (Hg), we conducted a genome-wide association study using relative root growth (RRG) and leaf accumulation measurements. Cd and Hg accumulation and RRG had heritability ranging 0.44 – 0.72 indicating high genetic diversity for these traits. The Cd and Hg trait associations were broadly distributed throughout the genome, indicated the traits are polygenic and involve several quantitative loci. For all traits, candidate genes included several membrane associated ATP-binding cassette transporters, P-type ATPase transporters, oxidative stress response genes, and stress related UDP-glycosyltransferases. The P-type ATPase transporters and ATP-binding cassette protein-families have roles in vacuole transport of heavy metals, and our findings support their wide use in physiological plant responses to heavy metals and abiotic stresses. We also found associations between Cd RRG with the genes CAX3 and PDR3, two linked adjacent genes, and leaf accumulation of Hg associated with the genes NRAMP6 and CAX9. When plant genotypes with the most extreme phenotypes were compared, we found significant divergence in genomic regions using population genomics methods that contained metal transport and stress response gene ontologies. Several of these genomic regions show high linkage disequilibrium (LD) among candidate genes suggesting they have evolved together. Minor allele frequency (MAF) and effect size of the most significant SNPs was negatively correlated with large effect alleles being most rare. This is consistent with purifying selection against alleles that increase toxicity and abiotic stress. Conversely, the alleles with large affect that had higher frequencies that were associated with the exclusion of Cd and Hg. Overall, macroevolutionary conservation of heavy metal and stress response genes is important for improvement of forage crops by harnessing wild genetic variants in gene banks such as the Medicago HapMap collection.

59 BASIC BIOLOGICAL SCIENCES↗