SEARCH · Engineering Papers

Results for “genetic models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

NASA NTRS1999

GH/IGF-I Transgene Expression on Muscle Homeostasis

We propose to test the hypothesis that the growth hormone/ insulin like growth factor-I axis through autocrine/paracrine mechanisms may provide long term muscle homeostasis under conditions of prolonged weightlessness. As a key alternative to hormone replacement therapy, ectopic production of hGH, growth hormone releasing hormone (GHRH), and IGF-I will be studied for its potential on muscle mass impact in transgenic mice under simulated microgravity. Expression of either hGH or IGF-I would provide a chronic source of a growth-promoting protein whose biosynthesis or secretion is shut down in space. Muscle expression of the IGF-I transgene has demonstrated about a 20% increase in hind limb muscle mass over control nontransgenic litter mates. These recent experiments, also establish the utility of hind-limb suspension in mice as a workable model to study atrophy in weight bearing muscles. Thus, transgenic mice will be used in hind-limb suspension models to determine the role of GH/IGF-I on maintenance of muscle mass and whether concentric exercises might act in synergy with hormone treatment. As a means to engineer and ensure long-term protein production that would be workable in humans, gene therapy technology will be used by to monitor muscle mass preservation during hind-limb suspension, after direct intramuscular injection of a genetically engineered muscle-specific vector expressing GHRH. Effects of this gene-based therapy will be assessed in both fast twitch (medial gastrocnemius) and slow twitch muscle (soleus). End-points include muscle size, ultrastructure, fiber type, and contractile function, in normal animals, hind limb suspension, and reambutation.

Schwartz, Robert J.↗

DOE OSTI2022

Artificial intelligence based analysis of nanoindentation load–displacement data using a genetic algorithm

In this work, we developed an automated tool, Nanoindentation Neo package for the analysis of nanoindentation load–displacement curves using a Genetic Algorithm (GA) applied to the Oliver-Pharr method (Oliver et al.,1992). For some materials, such as polycrystalline isotropic graphites, Least Squares Fitting (LSF) of the unload curve can produce unrealistic fit parameters. These graphites exhibit sharply peaked unloading curves not easily fit using the LSF, which tends to overestimate the indenter tip geometry parameter. To tackle this problem, we extended our general materials characterization tool Neo for EXAFS analysis (Terry et al., 2021) to fit nanoindentation data. Nanoindentation Neo automatically processes and analyzes nanoindentation data with minimal user input while producing meaningful fit parameters. GA, a robust metaheuristic method, begins with a population of temporary solutions using model parameters called chromosomes; from these we evaluate a fitness value for each solution, and select the best solutions to mix with random solutions producing the next generation. A mutation operator then modifies existing solutions by random perturbations, and the optimal solution is selected. We tested the GA method using Silica and Al reference standards. We fit samples of graphite and a high entropy alloy (HEA) consisting of BCC and FCC phases.

42 ENGINEERING↗

DOE OSTI2021

Systematic discovery of pseudomonad genetic factors involved in sensitivity to tailocins

Tailocins are bactericidal protein complexes produced by a wide variety of bacteria that kill closely related strains and may play a role in microbial community structure. Thanks to their high specificity, tailocins have been proposed as precision antibacterial agents for therapeutic applications. Compared to tailed phages, with whom they share an evolutionary and morphological relationship, bacterially produced tailocins kill their host upon production but producing strains display resistance to self-intoxication. Though lipopolysaccharide (LPS) has been shown to act as a receptor for tailocins, the breadth of factors involved in tailocin sensitivity, and the mechanisms behind resistance to self-intoxication, remain unclear. Here, we employed genome-wide screens in four non-model pseudomonads to identify mutants with altered fitness in the presence of tailocins produced by closely related pseudomonads. Our mutant screens identified O-antigen composition and display as most important in defining sensitivity to our tailocins. In addition, the screens suggest LPS thinning as a mechanism by which resistant strains can become more sensitive to tailocins. Furthermore, we validate many of these novel findings, and extend these observations of tailocin sensitivity to 130 genome-sequenced pseudomonads. This work offers insights into tailocin–bacteria interactions, informing the potential use of tailocins in microbiome manipulation and antibacterial therapy.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2003

Tactical Maneuvering Using Immunized Sequence Selection

This paper describes a tactical maneuvering system that uses an artificial immune system based approach for selecting maneuver sequences. This approach combines the problem solving abilities of genetic algorithms with the memory retention characteristics of an immune system. Of significant importance here is the fact that the tactical maneuvering system can make time-critical decisions to accomplish near-term objectives within a dynamic environment. These objectives can be received from a human operator, autonomous executive, or various flight planning specialists. Simulation tests were performed using a high performance military aircraft model. Results demonstrate the potential of using immunized sequence selection in order to accomplish tactical maneuvering objectives ranging from flying to a location while avoiding unforeseen obstacles, to performing relative positioning in support of air combat maneuvering.

Kaneshige, John↗

DOE OSTI2022

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Abstract Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Genotoxicity evaluation using primary hepatocytes isolated from rhesus macaque ( Macaca mulatta )

Non-human primates (NHPs) have played a vital role in fundamental, pre-clinical, and translational studies because of their high physiological and genetic similarity to humans. We report a method to isolate primary hepatocytes from the livers of rhesus macaques (Macaca mulatta) after in situ whole liver perfusion. Isolated primary macaque hepatocytes (PMHs) were treated with various compounds known to have different pathways of genotoxicity/carcinogenicity and the resulting DNA damage was evaluated using the high-throughput CometChip assay. The comet data were quantified using benchmark dose (BMD) modeling and the BMD50 values for treatments of PMHs were compared with those generated from primary human hepatocytes (PHHs) in our previous study. The results showed that despite varying CYP450 enzyme activities, PMHs had the same sensitivity and specificity as PHHs in detecting four indirect-acting (i.e., requiring metabolic activation) and seven direct-acting genotoxicants/carcinogens, as well as five non-carcinogens that are negative or equivocal for genotoxicity in vivo. The BMD50 estimates and their confidence intervals revealed species differences for DNA damage potency, especially for direct-acting compounds. The present study provides a practical method for maximizing the use of animal tissues by isolating primary hepatocytes from NHPs. Our data support the use of PMHs as a reliable surrogate of PHHs for evaluating the genotoxic hazards of chemical substances for humans.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Gene loss during a transition to multicellularity

Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited novel interactions, suggesting that gene loss could play a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest gene loss, in addition to gene gain and co-option, may be important for the evolution developmental complexity.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1999

Molecular Mechanisms Regulating Muscle Fiber Composition Under Microgravity

The overall goal of this project is to reveal the molecular mechanisms underlying the selective and debilitating atrophy of specific skeletal muscle fiber types that accompanies sustained conditions of microgravity. Since little is currently known about the regulation of fiber-specific gene expression programs in mammalian muscle, elucidation of the basic mechanisms of fiber diversification is a necessary prerequisite to the generation of therapeutic strategies for attenuation of muscle atrophy on earth or in space. Vertebrate skeletal muscle development involves the fusion of undifferentiated mononucleated myoblasts to form multinucleated myofibers, with a concomitant activation of muscle-specific genes encoding proteins that form the force-generating contractile apparatus. The regulatory circuitry controlling skeletal muscle gene expression has been well studied in a number of vertebrate animal systems. The goal of this project has been to achieve a similar level of understanding of the mechanisms underlying the further specification of muscles into different fiber types, and the role played by innervation and physical activity in the maintenance and adaptation of different fiber phenotypes into adulthood. Our recent research on the genetic basis of fiber specificity has focused on the emergence of mature fiber types and have implicated a group of transcriptional regulatory proteins, known as E proteins, in the control of fiber specificity. The restriction of E proteins to selected muscle fiber types is an attractive hypothetical mechanism for the generation of muscle fiber-specific patterns of gene expression. To date our results support a model wherein different E proteins are selectively expressed in muscle cells to determine fiber-restricted gene expression. These studies are a first step to define the molecular mechanisms responsible for the shifts in fiber type under conditions of microgravity, and to determine the potential importance of E proteins as upstream targets for the effects of weightlessness. In the past year we have determined that the expression of E Proteins is restricted to specific fiber types by post-transcriptional mechanisms. By far, the most prevalent mechanism of cellular control for achieving post-transcriptional regulation of gene expression is selective proteolysis -through the ubiquitin -proteasome pathway. Steady-state levels of HEB message are similar in all fast and slow skeletal muscle fiber types, yet the protein is restricted to Type IIX fibers. HEB appears to be a nodal point for regulating fiber-specific transcription, as expression of the transcription factor is regulated at the post-transcriptional level. It is not clear at present whether the regulation is at the level of protein synthesis or degradation. We are now poised to evaluate the biological role of ubiquitination in fiber specific-gene expression by controlling the post-transcriptional expression of E Proteins. The use of metabolic labelling and pharmacological inhibitors of the ubiquitin pathway will be used to identify the mode of regulation of the Type IIX expression pattern. The potential role of specific kinases in effecting the restriction of HEB expression will be examined by using both inhibitors and activators. The results of these studies will provide the necessary information to evaluate the biological role of E proteins in controlling fiber type transitions, and in potentially attenuating the atrophic effects of microgravity conditions. We have also recently shown that ectopic expression of the HEB protein transactivates the Type IIX-specific skeletal a-actin reporter. The 218 bp skeletal a-actin promoter drives transgene expression solely in mature Type IIX fibers. A mouse also carrying the transgene MLCI/HEB (which ectopically expresses the E Protein HEB in Type IIB fibers) forces expression of the skeletal a-actin reporter gene in Type IIB fibers. We can now dissect the composition of this fiber-specific cis-element. The skeletal a-actin promoter is quite compact and has been extensively characterized in vitro for activity and binding factors. The single E box may act as a binding target of myogenic factor/HEB heterodimer to allow for IIX expression. The HEB transcription factor may recognize either the precise flanking sequences of the E Box, or perhaps interacting with other proteins bound nearby, and activating expression in Type IIX fibers. This E box will be both ablated, and alternatively, as ablation may well destroy any muscle-specific transcriptional activity, flanking sequences substituted with those surrounding the E box (El) of the myogenin promoter. Modification of fiber-specific transgene expression will be tested in transgenic mice. The results of these studies will provide basic information on the regulatory circuitry underlying fiber specificity, and will form the basis for building appropriate transgenic regulatory cassettes to effect fiber transitions in subsequent experimental manipulations on unweighted muscles.

Rosenthal, Nadia A.↗

DOE OSTI2020

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Microbial Metabolite 3-Indolepropionic Acid Mediates Immunosuppression

The microbial-derived metabolite, 3-indolepropionic acid (3-IPA), has been intensely studied since its origins were discovered in 2009; however, 3-IPA’s role in immunosuppression has had limited attention. Untargeted metabolomic analyses of T-cell exhaustion and immunosuppression, represented by dysfunctional under-responsive CD8+ T cells, reveal a potential role of 3-IPA in these responses. T-cell exhaustion was examined via infection of two genetically related mouse strains, DBA/1J and DBA/2J, with lymphocytic choriomeningitis virus (LCMV) Clone 13 (Cl13). The different mouse strains produced disparate outcomes driven by their T-cell responses. Infected DBA/2J presented with exhausted T cells and persistent infection, and DBA/1J mice died one week after infection from cytotoxic T lymphocytes (CTLs)-mediated pulmonary failure. Metabolomics revealed over 70 metabolites were altered between the DBA/1J and DBA/2J models over the course of the infection, most of them in mice with a fatal outcome. Cognitive-driven prioritization combined with statistical significance and fold change were used to prioritize the metabolites. 3-IPA, a tryptophan-derived metabolite, was identified as a high-priority candidate for testing. To test its activity 3-IPA was added to the drinking water of the mouse models during LCMV Cl13 infection, with the results showing that 3-IPA allowed the mice to survive longer. This negative immune-modulation effect might be of interest for the modulation of CTL responses in events such as autoimmune diseases, type I diabetes or even COVID-19. Moreover, 3-IPA’s bacterial origin raises the possibility of targeting the microbiome to enhance CTL responses in diseases such as cancer and chronic infection.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1993

Partial melting of ordinary chondrites: Lost City (H) and St. Severin (LL)

Eucrites and diogenites are examples of asteroidal basalts and orthopyroxenites, respectively. As they are found intermingled in howardites, which are inferred to be regolith breccias, eucrites and diogenites are thought to be genetically related. But the details of this relationship and of their individual origins remain controversial. Work by Jurewicz et al. showed that 1170-1180 C partial melts of the (anhydrous) Murchison (CM) chondrite have major element compositions extremely similar to primitive eucrites, such as Sioux County. However, the MnO contents of these melts were about half that of Sioux County, a problem for the simple partial melting model. In addition, partial melting of Murchison could not produce diogenites, because residual pyroxenes in the Murchison experiments were too Fe- and Ca-rich and were minor phases at all but the lowest temperatures. A parent magma for diogenites needs an expanded low-calcium pyroxene field. In their partial melting study of an L6 chondrite, Kushiro and Mysen found that ordinary chondrites did have an expanded low-Ca pyroxene field over that of CV chondrites (i.e., Allende), probably because ordinary chondrites have lower Mg/Si ratios. This study expands that of both Kushiro and Mysen and Jurewicz et al. to the Lost City (H) and St. Severin (LL) chondrites at temperatures ranging from 1170 to 1325 C, at an fO2 of one log unit below the iron-wuestite buffer (IW-1).

Jurewicz, Amy J. G.↗

NASA NTRS1998

The molecular basis of ethylene signalling in Arabidopsis

The simple gas ethylene profoundly influences plants at nearly every stage of growth and development. In the past ten years, the use of a genetic approach, based on the triple response phenotype, has been a powerful tool for investigating the molecular events that underlie these effects. Several fundamental elements of the pathway have been described: a receptor with homology to bacterial two-component histidine kinases (ETR1), elements of a MAP kinase cascade (CTR1) and a putative transcription factor (EIN3). Taken together, these elements can be assembled into a simple, linear model for ethylene signalling that accounts for most of the well-characterized ethylene mediated responses.

Non-NASA Center↗

DOE OSTI2024

Molecular Dynamical and Quantum Mechanical Exploration of the Site-Specific Dynamics of Cy3 Dimers Internally Linked to dsDNA

Performing spectroscopic measurements on biomolecules labeled with fluorescent probes is a powerful approach to locating the molecular behavior and dynamics of large systems at specific sites within their local environments. The indocarbocyanine dye Cy3 has emerged as one of the most commonly used chromophores. The incorporation of Cy3 dimers into DNA enhances experimental resolution owing to the spectral characteristics influenced by the geometric orientation of excitonically coupled monomeric units. Various theoretical models and simulations have been utilized to aid in the interpretation of the experimental spectra. In this study, we employ all-atom molecular dynamics simulations to study the structural dynamics of Cy3 dimers internally linked to the dsDNA backbone. We used quantum mechanical calculations to derive insights from both the linear absorption spectra and the circular dichroism data. Furthermore, we explore potential limitations within a commonly used force field for cyanine dyes. The molecular dynamics simulations suggest the presence of four possible Cy3 dimeric populations. The spectral simulations on the four populations show one of them to agree better with the experimental signatures, suggesting it to be the dominant population. Furthermore, the relative orientation of Cy3 in this population compares very well with previous predictions from the Holstein–Frenkel Hamiltonian model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1998

The effect of track structure on the induction of chromosomal aberrations in murine cells

PURPOSE: To measure chromosome aberrations in C3H 10T1/2 mouse fibroblasts using FISH painting at the first mitosis following exposure to 30 keV/microm hydrogen or neon ions. MATERIALS AND METHODS: Cells in plateau-phase were irradiated with 0.86 MeV protons at the TTT-3 Tandem accelerator in Naples (Italy), or with 400 MeV/n Ne ions at the HIMAC accelerator in Chiba (Japan). Colcemid-blocked cells were harvested at the first mitosis following exposure, and chromosome spreads were hybridized in situ with a fluorescein-labelled composite mouse DNA probe specific for chromosomes 2 and 8. RESULTS: Protons were more efficient than neon ions at the same LET in the induction of chromosome interchanges and breaks. Yields of complex exchanges were similar for both particles at the same dose, but protons produced mostly insertions, while with Ne exposure non-reciprocal exchanges were the most frequent complex-type exchange. CONCLUSIONS: Charged particles with the same LET produce different yields of chromosome aberrations, and some observed differences can be explained based on the available track-structure models.

NASA Center JSC↗

DOE OSTI2020

Coexpressed subunits of dual genetic origin define a conserved supercomplex mediating essential protein import into chloroplasts

Significance Chloroplasts are of vital importance in photosynthetic eukaryotic organisms. Like mitochondria, they contain their own genomes. Nevertheless, most chloroplast proteins are encoded by nuclear genes, translated in the cytosol, and must cross the chloroplast envelope membranes to reach their proper destinations inside the organelle. Despite its fundamental role, our knowledge of the machinery catalyzing chloroplast protein import is incomplete and controversial. Here, we address the evolutionary conservation, composition, and function of the chloroplast protein import machinery using the green alga Chlamydomonas reinhardtii as a model system. Our findings help clarify the current debate regarding the composition of the chloroplast protein import machinery, provide evidence for cross-compartmental coordination of its biogenesis, and open promising avenues for its structural characterization.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

State-dependent motion of a genetically encoded fluorescent biosensor

Genetically encoded biosensors can measure biochemical properties such as small-molecule concentrations with single-cell resolution, even in vivo. Despite their utility, these sensors are “black boxes”: Very little is known about the structures of their low- and high-fluorescence states or what features are required to transition between them. We used LiLac, a lactate biosensor with a quantitative fluorescence-lifetime readout, as a model system to address these questions. X-ray crystal structures and engineered high-affinity metal bridges demonstrate that LiLac exhibits a large interdomain twist motion that pulls the fluorescent protein away from a “sealed,” high-lifetime state in the absence of lactate to a “cracked,” low-lifetime state in its presence. Understanding the structures and dynamics of LiLac will help to think about and engineer other fluorescent biosensors.

Rosen, Paul C. (ORCID:000000017414454X)↗

DOE OSTI2021

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism

Abstract Background Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737. Results We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10 −3 ), and combined dataset (p = 1.1 × 10 −4 ). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10 −35 , loss-of-function p = 2.26 × 10 −13 ) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10 −6 , OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia. Conclusions In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Systematic identification of cis -regulatory variants that cause gene expression differences in a yeast cross

Sequence variation in regulatory DNA alters gene expression and shapes genetically complex traits. However, the identification of individual, causal regulatory variants is challenging. Here, we used a massively parallel reporter assay to measure the cis-regulatory consequences of 5832 natural DNA variants in the promoters of 2503 genes in the yeast Saccharomyces cerevisiae. We identified 451 causal variants, which underlie genetic loci known to affect gene expression. Several promoters harbored multiple causal variants. In five promoters, pairs of variants showed non-additive, epistatic interactions. Causal variants were enriched at conserved nucleotides, tended to have low derived allele frequency, and were depleted from promoters of essential genes, which is consistent with the action of negative selection. Causal variants were also enriched for alterations in transcription factor binding sites. Models integrating these features provided modest, but statistically significant, ability to predict causal variants. This work revealed a complex molecular basis for cis-acting regulatory variation.

60 APPLIED LIFE SCIENCES↗