Search NASA⌕ Search

SEARCH · Search NASA

Results for “Genome-wide association studies”

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

Climate adaptation in Populus trichocarpa : key adaptive loci identified for stomata and leaf traits

We investigated adaptive genetic variation in Populus trichocarpa, a potential biofuel feedstock crop, to better understand how physiological traits may influence tolerance to water limitation. Our study focused on leaf and stomatal traits, given their roles in plant–water relations and adaptation. Using a diversity panel of over 1300 genotypes, we measured 14 leaf and stomatal traits under control (well-watered) and drought (water-limited) conditions. We conducted genome-wide association studies (GWAS), climate association analyses, and transcriptome (RNA-seq) profiling to identify genetic loci associated with phenotypic variation and adaptation. Stomatal traits, including size and density, were correlated with the climate of origin, with genotypes from more arid regions tending to have smaller but denser stomata. GWAS identified multiple loci associated with trait variation, including a major-effect region on chromosome 10 linked to stomatal size and abaxial contact angle. This locus overlapped with a tandem array of 3-ketoacyl-CoA synthase (KCS) genes and showed strong allele–climate and gene expression associations. Our findings reveal genetic and phenotypic variation consistent with local adaptation and suggest that future climates may favor alleles associated with smaller stomata, particularly under increasing aridity. This work provides insights into climate adaptation and breeding strategies for resilience in perennial crops.

Populus trichocarpa↗
DOE OSTI2025

Identification of a QTL region for tomato brown rugose fruit virus resistance in Solanum pimpinellifolium

Abstract Tomato (Solanum lycopersicumL.), one of the most widely grown vegetables in the world, has been seriously impacted in the past decade by the emerging tomato brown rugose fruit virus (ToBRFV). ToBRFV is a seed-borne tobamovirus, with ability to overcome the commonly usedTm-2 2 resistance gene in tomato. The objective of this study was to conduct quantitative trait locus (QTL) mapping and identify single-nucleotide polymorphism (SNP) markers associated with ToBRFV resistance in tomato. Two F 2 populations were used for QTL mapping: One derived from a cross betweenS. pimpinellifoliumUSVL333 (PI 390718) × USVL332 (PI 390717) and another from ‘Moneymaker’ × USVL332 (PI 390717), with population sizes of 195 and 79 plants, respectively. The resistance trait was derived from theS. pimpinellifoliumaccession USVL332 (PI 390717). A major QTL for ToBRFV resistance was identified on chromosome 11 (SL4.0ch11), with the peak located at approximately 46.84 Mbp. This QTL spans a 22-kb interval between 46,825,788 bp and 46,847,421 bp, as determined through both genome-wide association study (GWAS) and QTL linkage mapping. Three SNP markers, SL4.0ch11_46825788, SL4.0ch11_46847421, and SL4.0ch11_46850215, demonstrated the most significant association with high LOD values (LOD = 13 in the Blink model) in GWAS analysis. In this genomic region, two disease resistance gene analogs, Solyc11g062150 (TIR-NBS-LRR resistance protein, Toll-Interleukin receptor) and Solyc11g062180 (disease resistance protein, leucine-rich repeat), were identified, which may serve as candidates for ToBRFV resistance. The QTL identified in this study could be valuable for plant breeders in facilitating tomato breeding with ToBRFV resistance.

Agriculture↗
DOE OSTI2026

Populus_trichocarpa_Breeding_Population_SNPs

These data are from the manuscript “Application of Genomic Prediction in a Populus trichocarpa Breeding Program”, by Brian J. Stanton, David Macaya-Sanz, Chanaka Roshan Abeyratne, David Kainer, Kathy Haiby, Austin Himes, Carlos Gantz, Gerald A. Tuskan, and Stephen P. DiFazio. The data are based on genome resequencing to approximately 10X depth on two collections of Populus trichocarpa trees from Oregon, Washington, California, and British Columbia. The first collection consists of 293 genets collected by Poplar Innovations LLC for a breeding program. The second collection consists of 961 trees collected for the purpose of genome-wide association studies. These genets were sequenced using short, paired-end Illumina sequence reads (Chhetri et al. 2019). Reads were aligned to the P. trichocarpa ′Stettler-14′ reference (Hofmeister et al. 2020), with minor modifications to correct mis-assemblies (Zhou et al. 2020), and variants were called as per methods described in (Abeyratne et al. 2023). Identified variants were filtered using GATK’s VariantFiltration tool (DePristo et al. 2011), with filter expression flag set to “AF < 0.01 || AF > 0.99 || QD < 10.0 || ExcessHet > 20.0 || FS > 10.0 || MQ < 58.0”. SNPs with severe departures from Hardy−Weinberg expectations (exact-test p< 0.01) were also removed using vcftools --hwe flag (Danecek et al. 2011), resulting in 15,627,211 bi-allelic SNPs. The data included here consist of 141,903 high quality bi-allelic genome-wide SNPs obtained by further filtering the original SNP dataset using vcftools with flags --maf 0.05, --max-maf 0.95, --max-missing 0.95, --min-meanDP 10.75, --max-meanDP 43.00, --thin 2000. Collectively, these filtering parameters removed SNPs with 1) a minor allele frequency ≤ 0.05; 2) proportion of missing data for individual loci exceeding 5%; 3) sequencing depth more than 2X mean-depth or less than 0.5X mean-depth; or 4) a distance of

09 BIOMASS FUELS↗
DOE OSTI2024

An orphan gene BOOSTER enhances photosynthetic efficiency and plant productivity

Organelle-to-nucleus DNA transfer is an ongoing process playing an important role in the evolution of eukaryotic life. Here, genome-wide association studies (GWAS) of non-photochemical quenching parameters in 743 Populus trichocarpa accessions identified a nuclear-encoded genomic region associated with variation in photosynthesis under fluctuating light. The identified gene, BOOSTER (BSTR), comprises three exons, two with apparent endophytic origin and the third containing a large fragment of plastid-encoded Rubisco large subunit. Higher expression of BSTR facilitated anterograde signaling between nucleus and plastid, which corresponded to enhanced expression of Rubisco, increased photosynthesis, and up to 35% greater plant height and 88% biomass in poplar accessions under field conditions. Overexpression of BSTR in Populus tremula × P. alba achieved up to a 200% in plant height. Similarly, Arabidopsis plants heterologously expressing BSTR gained up to 200% in biomass and up to 50% increase in seed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations

Abstract Motivation Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data. Results We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization. Availability and implementation It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Du, Lin (ORCID:0009000289869812)↗
DOE OSTI2026

A single genomic region controls primocane fruiting in tetraploid blackberry

The fresh-market blackberry ( Rubus subgenus Rubus ) industry has expanded dramatically in the past 2 decades, driven in part by improved cultivars. Introgression of the primocane-fruiting (PF; annual flowering) trait into elite germplasm has enabled dual cropping in a single year, season extension, and cultivation in tropical and subtropical regions. Despite its economic performance, the genetic basis of PF is not well understood. It has been proposed that the PF trait is controlled by a major recessive locus, but its genomic location is unclear. Here, a genome-wide association study (GWAS) of 365 tetraploid blackberry genotypes identified a single genomic region on chromosome Ra03 (∼33 Mb) strongly associated with PF. Genetic linkage analysis in a biparental population confirmed that the same interval (32–35 Mb) was linked to the PF phenotype. Ten putative candidate genes were identified in this region. Allele mining using whole-genome resequencing of 17 genotypes highlighted 2 high-priority candidates: a CCCH-type zinc finger gene and an ubiquitin-specific protease gene. Use of an improved Rubus argutus “Hillquist” genome annotation (v1.2) enabled refined variant interpretation, including identification of regulatory 3′ UTR polymorphisms in the zinc finger homolog. Two diagnostic KASP markers (PF1 and PF2), designed from the most significant GWAS SNPs, predicted the PF phenotype with over 96% accuracy in a validation panel of 494 tetraploid blackberries from multiple breeding programs. Together, these results provide the first high-resolution mapping of the PF locus in blackberry, identify candidate genes for flowering regulation in Rubus , and deliver diagnostic markers that can be immediately deployed in breeding programs.

GWAS↗
DOE OSTI2026

Transcriptomic and functional analyses uncover a conserved effector driving genotype-dependent virulence in the Sphaerulina musiva-Populus trichocarpa interaction

The introduction of invasive microbes compromises the structure, biodiversity, and function of naïve ecosystems. Sphaerulina musiva, a hemibiotrophic pathogen that causes leaf spot and stem cankers in Populus species, exemplifies an invasive fungal pathogen spread by human activities. However, the genetic mechanisms of pathogenicity and virulence are poorly understood, impeding mitigation strategies. We utilized RNA sequencing to identify fungal effectors linked to stem canker formation, informing the development of future strategies for effective disease management. Our analysis revealed 70 genes differentially expressed at 2 weeks and 110 genes at 3 weeks between inoculated trees and controls. Notably, the gene with the highest expression at 2 weeks and the second highest at 3 weeks was homologous to Extracellular protein 2 (Ecp2). Complementary genome-wide association studies linked sequence polymorphisms in this locus to phenotypic variation in disease severity. Infiltration of S. musiva Ecp2 into Populus trichocarpa leaves induced necrosis in susceptible genotypes. Gene disruption using a CRISPR-Cas9 RNP system resulted in a genotype-dependent reduction of stem canker and disease severity. Tracing the evolutionary history of this effector across the fungal kingdom, we uncovered clade-specific gene-family expansions and orthologs in new species. These findings raise questions about the function and adaptive significance of these gene families in fungal lifestyles. Our study provides the first tractable target for breeding resistant poplar genotypes, addressing the challenges of managing S. musiva and uncovering mechanisms that drive its virulence, and provides deeper insights into the evolutionary dynamics of a conserved small-secreted protein with a diversity of functions.

Sondreli, Kelsey L [Oregon State University]↗
DOE OSTI2024

Chromosome-level genome assemblies and genetic maps reveal heterochiasmy and macrosynteny in endangered Atlantic Acropora

Abstract Background Over their evolutionary history, corals have adapted to sea level rise and increasing ocean temperatures, however, it is unclear how quickly they may respond to rapid change. Genome structure and genetic diversity contained within may highlight their adaptive potential. Results We present chromosome-scale genome assemblies and linkage maps of the critically endangered Atlantic acroporids,Acropora palmataandA. cervicornis. Both assemblies and linkage maps were resolved into 14 chromosomes with their gene content and colinearity. Repeats and chromosome arrangements were largely preserved between the species. The family Acroporidae and the genusAcroporaexhibited many phylogenetically significant gene family expansions. Macrosynteny decreased with phylogenetic distance. Nevertheless, scleractinians shared six of the 21 cnidarian ancestral linkage groups as well as numerous fission and fusion events compared to other distantly related cnidarians. Genetic linkage maps were constructed from oneA. palmatafamily and 16A. cervicornisfamilies using a genotyping array. The consensus maps span 1,013.42 cM and 927.36 cM forA. palmataandA. cervicornis, respectively. Both species exhibited high genome-wide recombination rates (3.04 to 3.53 cM/Mb) and pronounced sex-based differences, known as heterochiasmy, with 2 to 2.5X higher recombination rates estimated in the female maps. Conclusions Together, the chromosome-scale assemblies and genetic maps we present here are the first detailed look at the genomic landscapes of the critically endangered Atlantic acroporids. These data sets revealed that adaptive capacity of Atlantic acroporids is not limited by their recombination rates. The sister species maintain macrosynteny with few genes with high sequence divergence that may act as reproductive barriers between them. In the AtlanticAcropora, hybridization between the two sister species yields an F1 hybrid with limited fertility despite the high levels of macrosynteny and gene colinearity of their genomes. Together, these resources now enable genome-wide association studies and discovery of quantitative trait loci, two tools that can aid in the conservation of these species.

Biotechnology & Applied Microbiology↗
DOE OSTI2025

RatXcan: A framework for cross-species integration of genome-wide association and gene expression data

Genome-wide association studies (GWAS) have implicated specific alleles and genes as risk factors for numerous complex traits. However, translating GWAS results into biologically and therapeutically meaningful discoveries remains extremely challenging. Most GWAS results identify noncoding regions of the genome, suggesting that differences in gene regulation are the major driver of trait variability. To better integrate GWAS results with gene regulatory polymorphisms, we previously developed PrediXcan (also known as “transcriptome-wide association studies” orTWAS), which maps SNPs to predicted gene expression using GWAS data. In this study, we developed RatXcan, a framework that extends this methodology to outbred heterogeneous stock (HS) rats. RatXcan accounts for the close familial relationships among HS rats by modeling the relatedness with a random effect that encodes the genetic relatedness. RatXcan also corrects for polygenic-driven inflation because of the equivalence between a relatedness random effect and the infinitesimal polygenic model. To develop RatXcan, we trained transcript predictors for 8,934 genes using reference genotype and expression data from five rat brain regions. We found that the cis genetic architecture of gene expression in both rats and humans was sparse and similar across brain tissues. We tested the association between predicted expression in rats and two example traits (body length and BMI) using phenotype and genotype data from 5,401 densely genotyped HS rats and identified a significant enrichment between the genes associated with rat and human body length and BMI. Thus, RatXcan represents a valuable tool for identifying the relationship between gene expression and phenotypes across species and paves the way to explore shared biological mechanisms of complex traits.

Genetics & Heredity↗
NASA NTRSDate not supplied

BLOOD-BASED MULTI-SCALE MODEL FOR CANCER RISK FROM GCR IN GENETICALLY DIVERSE POPULATIONS

OBJECTIVES AND METHODS This project addresses the challenge of understanding and predicting individual radiation sensitivity by integrating genetics, demographics and biomarker characteristics across species (mice and humans). We hypothesize that ex vivo DNA repair response to GCR components is a central determinant of cancer risk from space radiation and can serve as a biomarker of radiation risk in combination with genetics. Automated image quantification of 53BP1+ radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET in non-immortalized primary skin fibroblasts derived from 76 mice across 15 strains (5 inbred reference strains and 10 collaborative-cross strains) exposed to X rays (0.1, 1 and 4 Gy), 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100μm2), as well as in peripheral blood mononuclear cells (PBMCs) from 768 healthy donors (matched ethnicity, 50/50 male/female, 18-70 years old) exposed to gamma rays (0.1 and 1 Gy), 350 MeV/n 28Si, 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100μm2). QUANTIFICATION OF 53BP1+ FOCI IN VITRO AND ASSOCIATIONS TO IN VIVO RADIATION SUSCEPTIBILITY IN 15 MOUSE STRAINS We reported in vitro repair kinetic and repairable fractions of RIF for the 15 mouse strains and introduced a mathematical model for RIF as a function of time, dose and LET. We noted that the metabolic activity of cells modulates the RIF response, and we introduced the open access tool terRIFic (Tool for Enhanced Results of RIF In Cells, https://radbiolab.shinyapps.io/terrific/) to correct for such bias using confluence level. Notably, at 4h post-irradiation, RIF/Gy decreased with dose or LET: as the dose or LET increases, so does the proximity of DNA double-strand-breaks (DSB) and our data suggest that proximal DSBs are brought together inside isolated RIF for repair. The RIF/Gy trend was inverted at 24h, suggesting RIF with high DSB content are more difficult to repair. We showed that in vitro metrics correlate with in vivo measurements in the same 15 mouse strains, such as survival levels of immune cells or spontaneous cancer incidence, suggesting a relationship between the efficiency of DSB repair and cancer risk or radiation toxicity. In addition to the efficiency of repair and persistent RIF, the amount of spontaneous foci before irradiation was also found to be strain dependent. Finally, we performed genome-wide association study in the same 15 mouse strains using all RIF phenotypes measured in vitro, identifying genes of interest and validating RIF as an ideal biomarker for individual radiation sensitivity. BASELINE 53BP1+ FOCI PREDICTS INDIVIDUAL HUMAN RESPONSE TO GCR COMPONENTS Based on the analysis of radiation responses of 576 donor PBMCs (using quantification of 53BP1+ foci, oxidative stress and cell death), we observed a wide variability of subject- and LET-dependent radiation responses, with radiation-induced DNA repair foci increasing with LET, though oxidative stress being notably reduced by high-LET irradiation, potentially due to a switch between hydrogen peroxide and oxygen radical-based mechanisms. We identified a relationship between few spontaneous DNA foci at baseline and increased DNA repair after irradiation, accompanied by an alteration in immunoregulatory cytokine secretion, which might be adapted as biomarkers to predict ionizing radiation sensitivity. Among demographic variables, only latent cytomegalovirus infection and age were predictive of high baseline foci formation. Finally, we have performed low-throughput whole genome sequencing of all samples and are currently in the process of identifying the genes and pathways associated with low and high-LET ionizing radiation sensitivity in humans.

53BP1↗
NASA NTRSDate not supplied

An Open-Science Approach to Address Individual Response to Simulated GCR In Genetically Diverse Populations of Mice and Humans

This project addresses the challenge of understanding and predicting individual radiation sensitivity by integrating genetics, demographics and biomarker characteristics across species (mice and humans). We hypothesize that ex vivo DNA repair response to GCR components is a central determinant of cancer risk from space radiation and can serve as a biomarker of radiation risk in combination with genetics. Automated image quantification of 53BP1+ radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET in non-immortalized primary skin fibroblasts derived from 76 mice across 15 strains (5 inbred reference strains and 10 collaborative-cross strains) exposed to X rays (0.1, 1 and 4 Gy), 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100sq. μm), as well as in peripheral blood mononuclear cells (PBMCs) from 768 healthy donors (matched ethnicity, 50/50 male/female, 18-70 years old) exposed to gamma rays (0.1 and 1 Gy), 350 MeV/n 28Si, 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100sq. μm). A genome-wide association study (GWAS) was performed on the mouse strains between DNA damage responses to space radiation and single nucleotide polymorphisms (SNPs). We found SNPs, which were significantly associated to the RIF phenotype, mapped to genes and pathways that are functionally linked to health hazards for deep space exploration (e.g. carcinogenesis, nervous system damage and immune dysfunction). Some of these SNPs were located within protein coding regions, potentially interfering with protein functions and providing promising genetic targets for countermeasures. We also found correlations between both spontaneous and radiation-induced DNA damage and SNPs mapped to pathways associated with cellular metabolism. GWAS is undergoing for the human data. All data have been made available via the NASA Space Biology Open-Science database (genelab.nasa.gov) and we will discuss how various genomic and transcriptomic datasets can be accessed for modeling and integrated using machine learning methods for discovering new radiation biology.

Sylvain V Costes↗
NASA NTRS2018

Variability in Galactic Cosmic Radiation- Induced DNA Damage Response in Inbred Mice Is Modulated by Genetics

In radiation biology, the ability to predict cancer risk associated with exposure to low doses of high-LET (Linear Energy Transfer) ionizing radiation remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual sensitivity. We have hypothesized that DNA repair capacity is the primary factor differentiating peoples radiation sensitivity. We previously showed in immortalized human cell lines that characterizing the dose and time dependence of p53-binding protein 1 (53BP1) foci formation in the nucleus following X-rays exposure is sufficient to predict DNA repair response to any other LET in the same cell line. We now tested this hypothesis across a population of mice with different genetic background. Fibroblast cells were extracted and cultivated from 76 individual mice from 15 different strains and exposed to HZE (high (H) atomic number (Z) and energy (E) galactic cosmic ray particles) particles and X-rays. Individual radiation sensitivities were investigated by high throughput measurement of DNA repair kinetics that evaluated 53bp1 foci numbers as a surrogate for DNA double-strand breaks at various times post-irradiation. Instead of just counting foci which can be hard to distinguish for high-LET or high doses, we also took into account the track structure of high-LET particles to compute the remaining number of unrepaired tracks as a function of time post-irradiation. As expected, the percentage of unrepaired track over a 48 hours follow-up period increased with LET. In addition, repair rate was modulated by genetics, with animals from the same strain showing small variance while large rate differences were observed between strains. Radiation strain sensitivity ranking was estimated based on repair rates from exposure to each LET evaluated in this work, and ranking for high-LET correlated better with ranking from high dose of X-ray, not low dose. At the in-vivo level, drops in T-cells and B-cells number measured 24 hours after 0.1 Gy (Gray) X-ray exposure, correlated with slower DNA repair kinetic in fibroblast cells of the same strains of mice. At the genomic level, mouse genome wide association (GWA) analysis identified seven significant genetic loci on chromosomes 2, 3, 7, 10, 11, 13 and 19 with different significance depending on the LET. Interestingly, for the two highest LET, a common locus on Chromosome 10 was identified with high enrichment for DNA repair associated genes.Overall, this work suggests that repair kinetics of primary skin fibroblasts is a good surrogate marker for in-vivo radiation sensitivities in other tissues and that this response is modulated by genetics. Our study also confirms that DNA repair kinetics following high doses of X-ray can be used to predict radiation sensitivity to high-LET.

High-LET Radiation↗
NASA NTRS2019

DNA Damage Response to Low and High-LET in a Large Cohort of Mice and Humans and Latest Advancement in NASA Space Omics

This presentation will first focus on a thorough evaluation of the DNA damage response to both low and high-LET in a cohort of 76 mice primary skin fibroblast derived from 15 different strains or in human blood mononuclear cells derived from 550 healthy donors. In both the human and mice work, we have hypothesized that DNA repair capacity can be used as a marker to evaluate and differentiate individual radiation sensitivity. More specifically, this work is based on the concept that the combined time-dose dependence of radiation-induced foci (RIF) of p53-binding protein 1 (53BP1) following low-LET exposure contains sufficient information to infer sensitivity to any other LET. This work is one of the most extensive studies on the kinetics and possible genetic underpinnings of radiation-induced DNA damage and repair. Results on humans are still preliminary as we are still in the process of collecting and isolating primary blood mononuclear cells from 500 to 800 healthy subjects of European descent, 18-75 years of age, 50/50 male/female distribution. We have analyzed 53BP1+ RIF formation as well as oxidative stress and cell death in primary cells from 192 subjects in response to the same HZE particles as used in mice: 600 MeV/n Fe, 350 MeV/n Ar and 350 MeV/n Si, 1.1 and 3 particles/100m2, 4 and 24 hours after irradiation. The second part of the talk will focus on describing GeneLab: The NASA Systems Biology Platform for Space Omics Repository, Analysis and Visualization. NASA GeneLab is an open-access repository for omics datasets generated by biological experiments conducted in space or experiments relevant to spaceflight (e.g. simulated cosmic radiation, simulated microgravity, bed rest studies). Started as a repository designed to archive precious omics from space experiments, GeneLab has expanded its scope to maximize the intelligibility of the raw data (e.g. RNAseq, microarray, WGBS, metagenome), particularly for users with limited bioinformatics knowledge. As such GeneLab is now providing processed data derived from the raw data covering a large spectrum of omics (genome, epigenome, transcriptome, epitranscriptome, proteome, metabolome), to help users explore important questions: Which genes or proteins are expressed differently in space for various living organisms? What are the consequences arising from these changes? What specifics DNA mutations or epigenetic changes happen in space? What species or genetic features lead to better adaption to such a unique environment? In this presentation, we will report on the current and future objectives for GeneLab, and review recent published studies relating molecular changes observed in various animal models and tissue with microgravity, radiation, circadian rhythm, hydration and carbon dioxide conditions.

DNA repair kinetics↗
NASA NTRS2019

Predicting Cancer Risk from Ionizing Radiation

The ability to predict cancer risk associated with exposure to low doses of high-LET ionizing radiation (IR) remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual variance in IR sensitivity. We have hypothesized that DNA repair capacity can be used as a marker to evaluate and differentiate individual radiation sensitivity. More specifically, this work is based on the concept that the combined time-dose dependence of radiation-induced foci (RIF) of p53-binding protein 1 (53BP1) following low-LET exposure contains sufficient information to infer sensitivity to any other LET. Our hypothesis was tested in 15 different mouse strains as well as in primary human immune cells. We first approached individual ionizing radiation sensitivity in a mouse model by culturing primary skin fibroblasts extracted from 76 mice of 15 different genetic backgrounds and exposing them to HZE particles and X-rays. This work is one of the most extensive studies on the kinetics and possible genetic underpinnings of radiation-induced DNA damage and repair. Our results is in agreement with a DNA repair model we previously postulated, where nearby DNA double strand breaks (DSB) in the nucleus are brought together for more efficient repair, leading to RIF clustering. Such mechanism was evidenced by a specific dose and LET dependence of RIF numbers. Briefly, RIF quantification after low-LET X-ray exposure showed an asymptotic saturation for doses between 1 Gy and 4 Gy 4 hours post-irradiation across all 15 strains. The clustering of DSB across all strains also led to more RIF/Gy for lower LET (X-ray and 350 MeV/n Ar) than for higher LET (600 MeV/n Fe) 4 hours post-exposure. Considering the fact that the number of DSB/Gy should be independent of LET, our data suggest there are more DSB in individual RIF as the LET increases. RIF numbers for 24 and 48 hours post-exposure led to the inverse trend, with more remaining RIF/Gy for higher LET (by 600 MeV/n Fe). This result suggests cells have more difficulty resolving RIF from higher LET as they the number DSB/RIF increases. Note that for most conditions, the variance of RIF/Gy was small within individual animals of the same strain and large between strains, suggesting a strong genetics component. Furthermore, we present our preliminary data from an ongoing study on human genetic associations with IR sensitivity. To address the human variability in responses to HZE particle irradiation in a maximally comprehensive manner, we are in the process of collecting and isolating primary blood mononuclear cells from 768 healthy subjects of European descent, 18-75 years of age, 50/50 male/female distribution. We have analyzed 53BP1+ RIF formation as well as oxidative stress and cell death in primary cells from 192 subjects in response to the same HZE particles as used in mice: 600 MeV/n Fe, 350 MeV/n Ar and 350 MeV/n Si, 1.1 and 3 particles/100m2, 4 and 24 hours after irradiation. We will next complete the quantification of HZE particle-induced DNA and cellular damage in the remaining subjects and compare it to their responses to low-LET irradiation. Finally, we will perform GWAS analysis to identify human genomic associations with IR sensitivity and potential targets for biomarker development.

High-LET radiation↗
DOE OSTI2025

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm’s behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl .

Gorstein, Evan↗
DOE OSTI2026

Assessing the Application of a Genomic Network Analysis in Population Ecology: Inferring Patterns of Dispersal and Geographic Structure in the Emerging Pathogen, Coccidioides

A challenge in population ecology studies is identifying how to best group individuals into populations, especially when individual origin is unknown. Machine learning has improved upon traditional methods of identifying population structure and is more efficient at handling large, complex datasets. We demonstrate the applicability of a machine learning method to identify hierarchical population structure in an emerging pathogen, Coccidioides spp., the causative agent of Valley fever. We compared the network clusters to structure identified by traditional tools as a validation of the network performance. We used publicly available whole-genome data for 48 C. immitis and 102 C. posadasii, resulting in 168,211 genome-wide SNPs among the two species. The network analysis grouped samples into populations comparable to the literature for these species but also identified fine-scale geographic structure and travel-associated cases not reported thus far. Exploring different resolutions in the network made it easy to identify unique genotypes specific to California and possibly Nevada, as well as Phoenix- and Tucson-acquired infections in non-endemic areas, regardless of reported travel history. The present study provides a promising example of how a ML-based network analysis can improve our ability to understand pathogen ecology, group cases into populations and infer travel-associated infections.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Editorial: Transcriptional and epigenetic landscapes of abiotic stress response in plants

In nature, plants constantly face various biotic and abiotic stresses that impact their growth, development, and productivity. Among these, abiotic stresses often have a more severe impact than biotic stresses. For instance, drought has been reported to cause greater yield losses than the combined impact of all plant pathogens (Gupta et al., 2020). Abiotic stresses are the immediate outcome of climate change, and the magnitude of these stresses has gradually increased every year with the rise in global temperatures. Thus, it has become imperative to study the impact of these stresses on plants and how plants respond to them at different levels to show resilient traits. This includes analysing the plants at morpho-physiological, biochemical, and molecular levels. Researchers often compare stressed plants to control (non-stressed) plants or evaluate contrasting genotypes, such as tolerant and sensitive lines, to elucidate the mechanisms underlying stress responses. While these studies have provided some insights, a comprehensive understanding of the intricate mechanisms governing plant responses to abiotic stress remains largely unknown. Recent advances in next-generation tools and technologies have enabled researchers to dissect the molecular basis of plant stress responses at genomic, transcriptomic, proteomic, metabolomic, epigenetic and epigenomic levels. Among these, knowledge of the transcriptional/epigenomic landscape of the trait-associated variations is limited. Given the importance of transcriptional changes and histone modifications in abiotic stress responses, this Research Topic was edited to collage the knowledge available on transcriptional and epigenetic landscapes of abiotic stress response in plants. The Research Topic features eight original research articles and one review, covering various aspects of transcriptome and epigenetic reprogramming in plants during abiotic stresses. Four of the research articles employ transcriptomics integrated with other omics approaches to explore transcriptome reprogramming, candidate gene identification, and the role of long non-coding RNA during different stresses. Two articles focus on the functional characterization of specific candidate genes involved in stress response, while another provides a genome-wide analysis of a stress-responsive gene family. Additionally, one study investigates genome-wide histone modifications, specifically H3K4me3 and H3K27me3, in response to abiotic stresses.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Prenatal 2g Exposure Alters Placental Expression of Stress-Related Genes

Extra-terrestrial colonization is of growing interest to space agencies and private entities, emphasizing the importance of research on reproduction and development in the absence of Earth's 1G. Maternal stressors can modify offspring development, exerting significant lifespan and crossgenerational changes through prenatal programming. The space environment is stressful, therefore exposure to altered gravity during pregnancy may impact later life outcomes in offspring. In ground-based studies, we exposed pregnant rats to continuous +G (above Earth gravity), and observed overweight and elevated anxiety in adult male (but not female) offspring, common phenotypes associated with prenatal maternal stress. Here we hypothesize that exposure to increased gravity during pregnancy elicits changes in the expression of stress-related genes in placenta that may mediate emergence of later life outcomes. While the placenta transports maternal factors to the fetus and produces endogenous fetal hormones, stress-induced changes at the placental-uterine interface may also alter communication between mother and fetus, facilitating prenatal transmission of unfavorable later life outcomes and cross-generational epigenetic alterations. Maternal stress elevates maternal glucocorticoids however placental 11b-hydroxysteroid dehydrogenase type 2 (HSD11B2) buffers fetal exposure by converting cortisol/corticosterone into inactive metabolites. Maternal stress during pregnancy down-regulates this enzyme and can induce epigenetic changes in placental and fetal tissues accounting for heightened adult HPA reactivity. Past studies have shown a placenta-specific increase in DNA methyltransferase (DNMT3a) mRNA in stressed mothers, an effect with implications for genome-wide epigenetic changes that may account for diverse phenotypic outcomes following maternal stress. Here we exposed groups of pregnant rats to one of five gravity loads (1, 1.5, 1.75 and 2G) and analyzed placental samples during late gestation. We predicted a systematic dose-response relationship between gravity load and the expression of the HSD11B2 and DNMT3 genes, thereby linking maternal exposure to altered gravity during pregancy with maternal stress.

Kalotia, Simranjit K.↗