Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “causal analysis”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Paleo-Megadroughts and Abrupt Climate Changes in the Speleothem Records. Final report

This project is motivated by the speleothem isotope records in Asia, which show regional responses in the hydrologic cycle to different climate forcings. Speleothem isotopic records are typically interpreted in terms of local precipitation variations or monsoon intensity. Our study demonstrates that non-local processes also play an important role. We started this project to understand the regional difference in speleothem isotopic composition between the Last Glacial Maximum (LGM) and the present-day. The record in Southwest China showed greater depletion during the LGM compared to those in East China. Our modeling and analysis showed that speleothems record, in addition, large scale changes in atmospheric circulation and moisture transport and their subsequent impact on precipitation. We developed an algorithm to partition total precipitation according to their formation dynamics, namely into frontal and non-frontal precipitation, and showed that the two have different trends and hence different causal mechanism. We then focused our subsequent attention on circulation impacts on precipitation changes. We applied a machine learning algorithm to detect rainbands in the ERA-Interim reanalysis product, and showed that the seasonal migrations of the rainbands are tied to the seasonal migrations of the jet stream, in particular the northerlies of the jet meanders. These northerlies, in turn, are partly topographic Rossby waves excited as the upstream westerlies impinge on the Tibetan Plateau. The seasonal variations of these upstream westerlies thus contribute to the seasonal movements of the rainbands and regional precipitation changes. Our analysis of the modern precipitation isotope record further confirms the importance of jet stream changes in the isotopic variations and shows that isotope-enriched years have reduced summer seasonality, with less pronounced northward migration of the jet.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On completeness and generalized symmetries in quantum field theory

We review a notion of completeness in QFT arising from the analysis of basic properties of the set of operator algebras attached to regions. In words, this completeness asserts that the physical observable algebras produced by local degrees of freedom are the maximal ones compatible with causality. We elaborate on equivalent statements to this completeness principle such as the non-existence of generalized symmetries and the uniqueness of the net of algebras. We clarify that for non-complete theories, the existence of generalized symmetries is unavoidable and further that they always come in dual pairs with precisely the same “size”. Moreover, the dual symmetries are always broken together, be it explicitly or effectively. Finally, we comment on several issues raised in recent literature, such as the relationship between completeness and modular invariance, dense sets of charges, and absence of generalized symmetries in the bulk of holographic theories.

Astronomy & Astrophysics↗
DOE OSTI2024

Stability analysis of WEST L-mode discharges with improved confinement from boron powder injection

WEST L-mode plasmas with dominant electron heating and no core torque source have observed improvements in confinement during boron (B) powder injection. These results are reminiscent of previous powder injection experiments on other devices and gaseous impurity seeding experiments on WEST. During powder injection, the stored energy increased up to 25% due to enhanced ion and electron heat and particle confinement. The improvements in confinement were not indicative of an L-H transition. To identify the dominant mechanisms and the causality chain behind these improvements in confinement, we employ interpretative modeling using METIS, predictive integrated modeling using a high-fidelity plasma simulator (HFPS), and stand-alone gyrokinetic simulations using quasi-linear gyrokinetic code. Interpretative modeling with METIS allowed for the estimation of missing data while maintaining good overall consistency with experiment. These results provided the initial conditions for fully predictive flux driven simulations using the HFPS. From these simulations, quasi-linear gyrokinetic analysis was performed at ρ = 0.5 and ρ = 0.65. Collisionality was found to be a strong candidate for the turbulence suppression mechanism at ρ = 0.5, while a combination of collisionality and the T e /T i ratio was found to be the likely mechanism at ρ = 0.65. The results additionally suggested that increased Z eff (through main ion dilution) could play a role in the improved confinement, but this could not be confirmed due to a lack of experimental measurements. The modeling framework established here can now be used to evaluate and exploit a variety of future powder injection experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Big-Data-Driven Geo-Spatiotemporal Correlation Analysis between Precursor Pollen and Influenza and its Implication to Novel Coronavirus Outbreak

Although studies of many respiratory viruses and pollens are often framed by both seasonal and health related perspectives, pollen has yet to be extensively examined as an important covariate to seasonal respiratory viruses (SRVs) in any context, including a causal one. This study contributes to those goals through an investigation of SRVs and pollen counts at selected regions across the Western Hemisphere. Two complementary decadal-scaled geospatial profiles were developed. One laterally spanned the US and was anchored by detailed pollen information for Albuquerque, New Mexico. The other straddled the equator to include Fortaleza, Brazil. We found that the geospatial and climatological patterns of pollen advancement and decline across the US every year presented a statistically significant correlation to the subsequent emergence and decline of SRVs. Other significant covariates included winds, temperatures, and atmospheric moisture. Our study indicates that areas of the US with lower geostrophic wind baselines are typically areas of persistently higher and earlier influenza like illness (ILI) cases. In addition to that continental- scaled contrast, many sites indicated seasonal highs of geostrophic winds and ILI which were closely aligned. These observations suggest extensive scale-dependent connectivity of viruses to geostrophic circulation. Pollen emergence and its own scale-dependent circulation may contribute to the geospatial and seasonal patterns of ILI. We explore some uncertainties associated with this investigation, and consider the possibility that in a temperate climate, following a Spring pollen emergence, a resulting increase in pollen triggered human Immunoglobulin E (IgE) antibodies may suppress ILIs for several months.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Herbaspirillum rubrisubalbicans as a Phytopathogenic Model to Study the Immune System of Sorghum bicolor

Herbaspirillum rubrisubalbicans is the causal agent of red stripe disease (RSD) and mottle stripe disease of sorghum and sugarcane, respectively. In all, 63 genotypes of Sorghum bicolor were inoculated with H. rubrisubalbicans, with 59 showing RSD symptoms. Quantitative trait loci (QTL) analysis in a recombinant inbred line (RIL) population identified several QTL associated with variation in resistance to RSD. RNA sequencing analysis identified a number of genes whose transcript levels were differentially regulated during H. rubrisubalbicans infection. Among those genes that responded to H. rubrisubalbicans inoculation were many involved in plant–pathogen interactions such as leucine-rich repeat receptors, mitogen-activated protein kinase 1, calcium-binding proteins, transcriptional factors (ethylene-responsive element binding factor), and callose synthase. Pretreatment of sorghum leaves with the pathogen-associated molecular pattern (PAMP) molecules flg22 and chitooctaose provided protection against subsequent challenge with the pathogen, suggesting that PAMP-triggered immunity plays an important role in the sorghum immunity response. These data present baseline information for the use of the genetically tractable H. rubrisubalbicans–sorghum pathosystem for the study of innate immunity and disease resistance in this important grain and bioenergy crop. Information gained from the use of this system is likely to be informative for other monocots, including those more intractable for experimental study (e.g., sugarcane).

Biochemistry & Molecular Biology↗
DOE OSTI2021

Limiting masses and radii of neutron stars and their implications

We combine the equation of state of dense matter up to twice nuclear saturation density 𝑛 sat obtained using chiral effective field theory (𝜒⁢EFT) and recent observations of neutron stars to gain insights about the high-density matter encountered in their cores. A key element in our study is the recent Bayesian analysis of correlated EFT truncation errors based on order-by-order calculations up to next-to-next-to-next-to-leading order in the 𝜒⁢EFT expansion. We refine the bounds on the maximum mass imposed by causality at high densities and provide stringent limits on the maximum and minimum radii of ∼1.4M ⊙ and ∼2.0M ⊙ stars. Including 𝜒⁢EFT predictions from 𝑛sat to 2𝑛 sat reduces the permitted ranges of the radius of a ~1.4M ⊙ star, 𝑅 1.4 , by ∼3.5km. If observations indicate 𝑅 1.4 < 11.2km, then our study implies that either the squared speed of sound 𝑐$^{2}_{𝑠}$ > 1/2 for densities above 2𝑛 sat or that 𝜒⁢EFT breaks down below 2𝑛 sat . Here, we also comment on the nature of the secondary compact object in GW190814 with mass ≃ 2.6M ⊙ and discuss the implications of massive neutron stars > 2.1⁢M ⊙ ⁡(2.6M ⊙ ) in future radio and gravitational-wave searches. Some form of strongly interacting matter with 𝑐$^{2}_{𝑠}$ > 0.35⁢(0.55) must be realized in the cores of such massive neutron stars. In the absence of phase transitions below 2𝑛 sat , the small tidal deformability inferred from GW170817 lends support for the relatively small pressure predicted by 𝜒⁢EFT for the baryon density 𝑛 B in the range 1–2𝑛 sat . Together they imply that the rapid stiffening required to support a high maximum mass should occur only when 𝑛 B ≳ 1.5–1.8 𝑛 sat .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A Portable Wave Tank and Wave Energy Converter for Engineering Dissemination and Outreach

Wave energy converters are a nascent energy generation technology that harnesses the power in ocean waves. To assist in communicating both fundamental and complex concepts of wave energy, a small-scale portable wave tank and wave energy converter have been developed. The system has been designed using commercial off-the-shelf components, and all design hardware and software are openly available for replication. This project builds on prior research conducted at Sandia National Laboratories, particularly in the areas of WEC device design and control systems. By showcasing the principles of causal feedback control and innovative device design, SIWEED not only serves as a practical demonstration tool but also enhances the educational experience for users. This paper presents the detailed system design of this tool. Furthermore, via testing and analysis, we demonstrate the basic functionality of the system.

educational↗
DOE OSTI2023

Plant defensin MtDef4–derived antifungal peptide with multiple modes of action and potential as a bio–inspired fungicide

Chemical fungicides have been instrumental in protecting crops from fungal diseases. However, increasing fungal resistance to many of the single-site chemical fungicides calls for the development of new antifungal agents with novel modes of action (MoA). The sequence-divergent cysteine-rich antifungal defensins with multisite MoA are promising starting templates for design of novel peptide-based fungicides. Here, we experimentally tested such a set of 17-amino-acid peptides containing the γ-core motif of the antifungal plant defensin MtDef4. These designed peptides exhibited antifungal properties different from those of MtDef4. Focused analysis of a lead peptide, GMA4CG_V6, showed that it was a random coil in solution with little or no secondary structure elements. Additionally, it exhibited potent cation-tolerant antifungal activity against the plant fungal pathogen Botrytis cinerea, the causal agent of grey mould disease in fruits and vegetables. Its multisite MoA involved localization predominantly to the plasma membrane, permeabilization of the plasma membrane, rapid internalization into the vacuole and cytoplasm, and affinity for the bioactive phosphoinositides phosphatidylinositol 3-phosphate (PI3P), PI4P, and PI5P. The sequence motif RRRW was identified as a major determinant of the antifungal activity of this peptide. While topical spray application of GMA4CG_V6 on Nicotiana benthamiana and tomato plants provided preventive and curative suppression of grey mould disease symptoms, the peptide was not internalized into plant cells. Our findings open the possibility that truncated and modified defensin-derived peptides containing the γ-core sequence could serve as promising candidates for further development of bio-inspired fungicides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subseasonal Representation and Predictability of North American Weather Regimes Using Cluster Analysis

Abstract This study focuses on assessing the representation and predictability of North American weather regimes, which are persistent large-scale atmospheric patterns, in a set of initialized subseasonal reforecasts created using the Community Earth System Model, version 2 (CESM2). The k -means clustering was used to extract four key North American (10°–70°N, 150°–40°W) weather regimes within ERA5 reanalysis, which were used to interpret CESM2 subseasonal forecast performance. Results show that CESM2 can recreate the climatology of the four main North American weather regimes with skill but exhibits biases during later lead times with overoccurrence of the West Coast high regime and underoccurrence of the Greenland high and Alaskan ridge regimes. Overall, the West Coast high and Pacific trough regimes exhibited higher predictability within CESM2, partly related to El Niño. Despite biases, several reforecasts were skillful and exhibited high predictability during later lead times, which could be partly attributed to skillful representation of the atmosphere from the tropics to extratropics upstream of North America. The high predictability at the subseasonal time scale of these case-study examples was manifested as an “ensemble realignment,” in which most ensemble members agreed on a prediction despite ensemble trajectory dispersion during earlier lead times. Weather regimes were also shown to project distinct temperature and precipitation anomalies across North America that largely agree with observational products. This study further demonstrates that unsupervised learning methods can be used to uncover sources and limits of subseasonal predictability, along with systematic biases present in numerical prediction systems. Significance Statement North American weather regimes are large-scale atmospheric patterns that can persist for several days. Their skillful subseasonal (2 weeks or greater) prediction can provide valuable lead time to prepare for temperature and precipitation anomalies that can stress energy and water resources. The purpose of this study was to assess the climatological representation and subseasonal predictability of four key North American weather regimes using a research subseasonal prediction system and clustering analysis. We found that the Pacific trough and West Coast high regimes exhibited higher predictability than other regimes and that skillful representation of conditions across the tropics and extratropics can increase predictability during later lead times. Future work will quantify causal pathways associated with high predictability.

58 GEOSCIENCES↗
DOE OSTI2024

Decision-making based on Markov decision process in integrated artificial reasoning framework—Part I: Theory

This paper presents a decision-making framework based on an integrated artificial reasoning framework and Markov decision process (MDP). The integrated artificial reasoning framework provides a physics-based approach that converts system information into state transition models, and the analysis result will be represented by the transition probabilities that can be used with an MDP to find a traceable and explainable optimal pathway. A dynamic Bayesian network (DBN) is well suited for representing the structure of an MDP. The causality information among process variables (or among subsystems) is mathematically represented in a DBN by the conditional probabilities of the node’s states provided different probabilities of the parent node’s states. To define node states in a physically understandable manner, we used multilevel flow modeling (MFM). An MFM follows the fundamental energy and mass conservation laws and supports the selection of process variables that represent the system of interest so that causal relations among process variables are properly captured. An MFM-based DBN supports developing state transition models in an MDP to capture the effect of process variables of system having physical relations. The operators of the target system can capture stochastic system dynamics as multiple subsystem state transitions based on their physical relations and uncertainties coming from component degradation or random failures. We analyzed a simplified exemplary system to illustrate an optimal operational policy using the suggested approach.

Markov decision process↗
DOE OSTI2023

Data driven pathway analysis and forecast of global warming and sea level rise

Climate change is a critical issue of our time, and its causes, pathways, and forecasts remain a topic of broader discussion. In this paper, we present a novel data driven pathway analysis framework to identify the key processes behind mean global temperature and sea level rise, and to forecast the magnitude of their increase from the present to 2100. Based on historical data and dynamic statistical modeling alone, we have established the causal pathways that connect increasing greenhouse gas emissions to increasing global mean temperature and sea level, with its intermediate links encompassing humidity, sea ice coverage, and glacier mass, but not for sunspot numbers. Our results indicate that if no action is taken to curb anthropogenic greenhouse gas emissions, the global average temperature would rise to an estimated 3.28 °C (2.46–4.10 °C) above its pre-industrial level while the global sea level would be an estimated 573 mm (474–671 mm) above its 2021 mean by 2100. However, if countries adhere to the greenhouse gas emission regulations outlined in the 2021 United Nations Conference on Climate Change (COP26), the rise in global temperature would lessen to an average increase of 1.88 °C (1.43–2.33 °C) above its pre-industrial level, albeit still higher than the targeted 1.5 °C, while the sea level increase would reduce to 449 mm (389–509 mm) above its 2021 mean by 2100.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Identifying human failure events (HFEs) for external hazard probabilistic risk assessment

In recent years, several advancements in nuclear power plant (NPP) probabilistic risk assessment (PRA) have been driven by increased understanding of external hazards, plant response, and uncertainties. However, major sources of uncertainty associated with external hazard PRA remain. One important source is how risk-significant human actions that are carried out to enable plant response and recovery from natural hazards cause the close coupling of physical impacts on plants and overall plant risk during these hazard events. This makes human reliability and human-plant interactions important elements to consider in resolving PRA gaps in external hazards. One of the challenges in considering human response in external hazard probabilistic risk assessment (XHPRA) is that most existing human reliability analysis (HRA) models were not developed for assessing actions outside the control room (termed ex-control room actions) and hazard response. To support this new scope, HRA models will need to be developed or modified to support identification of human activities, causal factors, and uncertainties inherent in external hazard response, thereby providing insights regarding event timing and physical event conditions as they relate to human performance. In this study, there are two main objectives: (1) evaluate the applicability of an existing cognitive-based HRA method, Phoenix, to ex-control room actions, and (2) identify sources of uncertainty to be characterized or reduced in order to make this method suitable for XHPRA. The first step of such work is performed by assessing the suitability of existing HRA methods to support identifying human failure events (HFEs) for human response to flooding hazards. These HFEs are human actions or inactions that are involved in human responses to flooding hazards and could contribute to the loss of a critical function for the plant in the scenario being examined. Here, in this work, decomposition analyses using the cognitive-based Phoenix HRA model are used to identify HFEs. The Phoenix method was found to be suitable for analyzing ex-control room actions as well as identifying specific HFEs and underlying crew failure modes (CFMs). However, the method's suitability for use in ex-control room actions would benefit from expanding the available CFMs to accommodate a larger variety of physical and communication tasks.

42 ENGINEERING↗
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Hazard-Based Model of Activity Generation Using Vehicle Trajectory Data

Travel behavior modeling suffers from lack of longitudinal data due to the high cost of procurement and difficulty with ensuring data quality. In the recent years ubiquitous data intercepted from social media, cell detail record, GPS traces are increasingly being explored to analyze human mobility. However, these data lack key demographic information which limits their use for developing causal models of traveler behavior. The primary contribution of the paper lies in using vehicle trajectory data from Ann Arbor along with multiple other data sources for developing a hazard-based activity generation model. The different category of data used in the analysis include (i) vehicle trajectory data (VTD), (ii) American Census Survey (ACS), (iii) land use data, (iv) household travel survey (HTS) data, and (v) smart location data from Environmental Protection Agency (EPA). Home and work are identified as the two most frequently visited locations from the VTD and validated against land use categories of the region. Next a random utility-based model is developed using HTS data to identify other activity purposes as a function of the polynomials of distance from home, activity duration and departure time. The developed model is then applied to identify the non-home activity purposes in the VTD. Demographic and accessibility information are appended to the VTD from ACS and smart location databases respectively. The fmal database is then used to develop a hazard-based activity generation model. (C) 2020 The Authors. Published by Elsevier B.V.

activity generation↗
DOE OSTI2024

XA21-mediated resistance to Xanthomonas oryzae pv. oryzae is dose dependent

The rice receptor kinase XA21 confers broad-spectrum resistance to Xanthomonas oryzae pv. oryzae ( Xoo ), the causal agent of rice bacterial blight disease. To investigate the relationship between the expression level of XA21 and resulting resistance, we generated independent HA-XA21 transgenic rice lines accumulating the XA21 immune receptor fused with an HA epitope tag. Whole-genome sequence analysis identified the T-DNA insertion sites in sixteen independent T0 events. Further, through quantification of the HA-XA21 protein and assessment of the resistance to Xoo strain PXO99 in six independent transgenic lines, we observed that XA21-mediated resistance is dose dependent. In contrast, based on the four agronomic traits quantified in these experiments, yield is unlikely to be affected by the expression level of HA-XA21 . These findings extend our knowledge of XA21-mediated defense and contribute to the growing number of well-defined genomic landing pads in the rice genome that can be targeted for gene insertion without compromising yield.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Image-based time series analysis to establish differential disease progression for two Fusarium head blight pathogens in oat spikelets with variable resistance

Oat-based value-added products have increased their value as healthy foodstuff. Fusarium head blight (FHB) infections and the mycotoxins accumulated to the oat seeds, however, pose a challenge to oat production. The FHB infections are predicted to become more prevalent in the future changing climates and under more limited use of fungicides. Both these factors increase the pressure for breeding new resistant cultivars. Until now, however, genetic links in oats against FHB infection have been difficult to identify. Therefore, there is a great need for more effective breeding efforts, including improved phenotyping methods allowing time series analysis and the identification of molecular markers during disease progression. To these ends, dissected spikelets of several oat genotypes with different resistance profiles were studied by image-based methods during disease progression by Fusarium culmorum or F. langsethiae species. The chlorophyll fluorescence of each pixel in the spikelets was recorded after inoculation by the two Fusarium spp., and the progression of the infections was analyzed by calculating the mean maximum quantum yield of PSII (F v /F m ) values for each spikelet. The recorded values were (i) the change in the photosynthetically active area of the spikelet as percentage of its initial size, and (ii) the mean of F v /F m values of all fluorescent pixels per spikelet post inoculation, both indicative of the progression of the FHB disease. The disease progression was successfully monitored, and different stages of the infection could be defined along the time series. The data also confirmed the differential rate of disease progression by the two FHB causal agents. In addition, oat varieties with variable responses to the infections were indicated.

oats↗
DOE OSTI2023

Operation Optimization using Reinforcement Learning with Integrated Artificial Reasoning Framework

In large and complex systems, operational decision-making requires a systematic analysis with a vast amount of data from both process parameters and component status monitoring. In this paper, we present an integrated artificial reasoning approach for system state transition models that can help operational decision-making with explainable and traceable reasoning. The integrated artificial reasoning framework is a physics-based approach of defining the system structure in a Bayesian network, so we leveraged it in a Markov decision process (MDP) for finding optimal operational solutions. In our proposed framework, the MDP is implemented on a dynamic Bayesian network (DBN), which represents causalities in a system. The multilevel flow modeling was utilized in order to extract these causalities in a more efficient and objective manner. Since multilevel flow modeling is based on the fundamental energy and mass conservation laws, the target system is decomposed into several mass, energy, and information structures, which serve as the basis for a DBN. The MDP consists of the processes of finding a solution for the Bellman equation, which can be derived from the conditional probability equations of the constructed DBN. System operators can capture stochastic system dynamics as multiple subsystem state transitions based on their physical relations and uncertainties coming from the component degradation process or random failures. We analyzed a simplified example system to illustrate finding an optimal operational policy with this approach.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 hours post-fertilization after exposure to eight retene concentrations. The concentrations were selected to produce effects ranging from no phenotype to mortality and malformations in 100% of animals at 5 days post-fertilization. We identified a concentration-response relationship between retene teratogenicity and differential gene expression in both number of DEGs and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of AHR-dependent teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-ß (TGF-ß) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module 7 genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module 7 genes are largely responsible for retene toxicity.

Toxin, Zebrafish, Retene, Transcriptomics, network↗