Search NASA⌕ Search

SEARCH · Search NASA

Results for “SNPs”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integrating functional scoring and regulatory data to predict the effect of non-coding SNPs in a complex neurological disease

Abstract Most SNPs associated with complex diseases seem to lie in non-coding regions of the genome; however, their contribution to gene expression and disease phenotype remains poorly understood. Here, we established a workflow to provide assistance in prioritising the functional relevance of non-coding SNPs of candidate genes as susceptibility loci in polygenic neurological disorders. To illustrate the applicability of our workflow, we considered the multifactorial disorder migraine as a model to follow our step-by-step approach. We annotated the overlap of selected SNPs with regulatory elements and assessed their potential impact on gene expression based on publicly available prediction algorithms and functional genomics information. Some migraine risk loci have been hypothesised to reside in non-coding regions and to be implicated in the neurotransmission pathway. In this study, we used a set of 22 non-coding SNPs from neurotransmission and synaptic machinery-related genes previously suggested to be involved in migraine susceptibility based on our candidate gene association studies. After prioritising these SNPs, we focused on non-reported ones that demonstrated high regulatory potential: (1) VAMP2_rs1150 (3′ UTR) was predicted as a target of hsa-mir-5010-3p miRNA, possibly disrupting its own gene expression; (2) STX1A_rs6951030 (proximal enhancer) may affect the binding affinity of zinc-finger transcription factors (namely ZNF423) and disturb TBL2 gene expression; and (3) SNAP25_rs2327264 (distal enhancer) expected to be in a binding site of ONECUT2 transcription factor. This study demonstrated the applicability of our practical workflow to facilitate the prioritisation of potentially relevant non-coding SNPs and predict their functional impact in multifactorial neurological diseases.

Felício, Daniela↗
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 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↗
NASA NTRSDate not supplied

Benchmarking Computational Tools for Calling SNPs and Indels in Complex Microbial Populations

The NASA BioNutrients missions seek to understand the suitability of microorganisms for bioproduction during space flight. One topic of interest is the stability of microbial genomes during long-term ambient storage and subsequent rehydration and growth. To address these questions, samples from 8 species were flown to ISS for 5 years of desiccated storage at ambient temperature (Stasis Packs) and 2 species were packaged along with powdered media inside a bioreactor system to allow hydration and growth in microgravity (Production Packs). For both systems, Whole Genome Sequencing (WGS) of the DNA extracted from the returned samples and paired ground controls will be conducted to identify changes in genome stability due to time, storage conditions and growth in space. Across the technical replicates, ground controls, 10 timepoints, and multiple experimental conditions, ~300 samples have been selected for initial analysis with WGS sequencing to 100x coverage. A flexible and resource efficient mutation calling pipeline is needed to process this large dataset and allow for comparisons between species. Many bioinformatics tools for calling Indels and Single Nucleotide Variants (SNVs) are designed for use with pure isolates, where true variations from the reference genome are expected to dominate the reads aligning to the location of mutation. In contrast, DNA from the Stasis Pack (SP) samples was collected directly after recovery from desiccated storage and the Production Pack (PP) samples were collected after fermentation. In this context, reads with mutations are expected to be less frequent than reads that align with the reference genome, as each sample will include multiple lines of cells. Thus, BioNutrients samples are expected to be similar to samples from cancer cell or “pooled” sequencing approaches. In preparation for the analysis of the BioNutrients samples, we have tested three mutation calling tools (GATK for Microbes, BreSeq and DiscoSNP) designed for complex samples. A challenge of validating mutation identification pipelines is a lack of “Ground Truth” datasets, especially for complex samples. To compare these three tools, we sought to identify mutations in pre-existing WGS data collected from populations of Chlamydomonas reinhardtii that were exposed to UV mutagenesis and growth in LEO as part of the Space Algae-1 mission. Here we present a summary of these tools against the analysis originally conducted using the CRISP tool. Critical metrics are compared such as runtime, the number of SNPs, the number and size of Indels, and patterns of transversion and transitions identified by each tool are reported. By sharing these benchmarking results collected in support of the BioNutrients mission, we aim to guide others seeking to identify SNVs in similarly complex microbial samples.

Biology↗
DOE OSTI2022

Evaluating the Impact of gRNA SNPs in CasRx Activity for Reducing Viral RNA in HCoV-OC43

Viruses within a given family often share common essential genes that are highly conserved due to their critical role for the virus’s replication and survival. In this work, we developed a proof-of-concept for a pan-coronavirus CRISPR effector system by designing CRISPR targets that are cross-reactive among essential genes of different human coronaviruses (HCoV). We designed CRISPR targets for both the RNA-dependent RNA polymerase (RdRp) gene as well as the nucleocapsid (N) gene in coronaviruses. Using sequencing alignment, we determined the most highly conserved regions of these genes to design guide RNA (gRNA) sequences. In regions that were not completely homologous among HCoV species, we introduced mismatches into the gRNA sequence and tested the efficacy of CasRx, a Cas13d type CRISPR effector, using reverse transcription quantitative polymerase chain reaction (RT-qPCR) in HCoV-OC43. We evaluated the effect that mismatches in gRNA sequences has on the cleavage activity of CasRx and found that this CRISPR effector can tolerate up to three mismatches while still maintaining its nuclease activity in HCoV-OC43 viral RNA. This work highlights the need to evaluate off-target effects of CasRx with gRNAs containing up to three mismatches in order to design safe and effective CRISPR experiments.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

One-Carbon Metabolism and SANS 2022 Update

Spaceflight Associated Neuro-ocular Syndrome, or SANS, affects a subset of astronauts (1, 2), and biochemical evidence has documented differences in those astronauts (3). Specifically, they had higher circulating concentrations of metabolites of the one-carbon metabolic pathway (1C), including homocysteine, and these concentrations were higher before flight (3). After ruling out many potential confounding factors in these otherwise healthy individuals (e.g., sex, kidney function, vitamin status, coffee consumption), a study of genetics was warranted. In an initial pilot effort, we documented a genetic predisposition to develop ophthalmic changes after long-duration space flight (4). That is, from a limited study of 5 single-nucleotide polymorphisms (SNPs), we found that the G allele for the MTRR A66G SNP was associated with a greater risk of choroidal folds and cotton-wool spots after flight, and the C allele for SHMT1 C1420T was protective against optic disc edema (4). These data provide a potential pathway for understanding why some individuals develop SANS, while others do not. The initial pilot study of 5 SNPs yielded striking findings, but the 1C pathway is far more complex. An effort was undertaken to examine more than 500 1C SNPs to see if a broader examination could help illuminate this association. That work is ongoing. The astronaut findings led us to advocate for the inclusion of 1C pathway genetic and biochemistry testing on other SANS-related projects, noting that genetics might help identify responders, non-responders, or outliers. The first such effort yielded evidence of an association of specific forms of the MTRR and SHMT-1 SNPs and vitamin B12 status with end-tidal CO2 after acute carbon dioxide exposure (5). The second such effort led to the identification that individuals exposed to strict head-down tilt and CO2 for 30-d who developed optic disc edema also had risk alleles for the two SNPs described above (6). Additionally, we identified a clinical population with many characteristics either attributed or purported to be involved in the ocular changes seen in affected astronauts: women with polycystic ovary syndrome (PCOS). PCOS is a condition of androgen excess and anovulatory menstrual cycles. The shared characteristics and clinical findings between SANS and PCOS generally include higher circulating homocysteine concentrations, increased retinal nerve fiber layer thickness, increased androgen concentrations (or responses), and altered carbohydrate metabolism. To our knowledge, no study has examined whether women with PCOS have asymptomatic ophthalmic anomalies observed in astronauts with SANS. While researchers have evaluated the one-carbon metabolism pathway polymorphisms of PCOS patients, and initial studies show an association with certain one-carbon polymorphisms, none have looked at the set of SNPs identified in our studies that are associated with ophthalmic changes in astronauts. Accordingly, we designed a study to evaluate the association of one-carbon pathway SNPs and ophthalmic findings in patients with PCOS and/or IIH compared to controls. Subjects provided blood samples for vitamin and one carbon biochemistry analyses, an extensive analysis of >500 SNPs associated with one carbon metabolism and had eye examinations and ocular imaging. Data analysis are underway. The data collected to date have shown associations between one carbon pathway biochemistry and genetics and incidence of SANS. The mechanisms for SANS has yet to be identified, although many hypotheses exist. Based on our data, we have developed (8, 9) and expanded (6) a multi-hit hypothesis for how these seemingly disparate findings could be linked. While intriguing, the hypothesis represents the starting point for further research. We aim to clarify the relationship between B-vitamin status and genetics with regard to the risk of SANS. Ultimately, understanding the mechanism(s) behind this will provide a means to predict, prevent, or treat these ophthalmologic pathologies in astronauts, and terrestrial populations.

S M Smith↗
DOE OSTI2025

Identification of functional non-coding variants associated with orofacial cleft

Oral facial cleft (OFC) comprises cleft lip with or without cleft palate (CL/P) or cleft palate only. Genome wide association studies (GWAS) of isolated OFC have identified common single nucleotide polymorphisms (SNPs) in many genomic loci where the presumed effector gene (for example, IRF6 in the 1q32 locus) is expressed in embryonic oral epithelium. To identify candidates for functional SNPs at eight such loci we conduct a massively parallel reporter assay in a fetal oral epithelial cell line, revealing SNPs with allele-specific effects on enhancer activity. We filter these SNPs against chromatin-mark evidence of enhancers and test a subset in traditional reporter assays, which support the candidacy of SNPs at loci containing FOXE1, IRF6, MAFB, TFAP2A, and TP63. For two SNPs near IRF6 and one near FOXE1, we engineer the genome of induced pluripotent stem cells, differentiate the cells into embryonic oral epithelium, and discover allele-specific effects on the levels of effector gene expression, and, in two cases, the binding affinity of transcription factors FOXE1 or ETS2. Conditional analyses of GWAS data suggest the two functional SNPs near IRF6 account for the majority of risk for CL/P at this locus. This study connects genetic variation associated with OFC to mechanisms of pathogenesis.

Kumari, Priyanka↗
DOE OSTI2022

Linkage map construction using limited parental genotypic information

Abstract Genetic linkage maps based on single nucleotide polymorphisms (SNPs) represent an essential tool for a variety of genomic analyses. Today, next-generation sequencing (NGS) enables rapid genotyping of different mapping populations based on thousands of SNPs and the construction of highly saturated linkage maps. Nevertheless, missing data in the genotyping of the parental lines creates a bottleneck that determines the number of SNPs that can be used for the linkage map. As a proof of concept, a highly saturated genetic linkage map was constructed using the imputed genotypic data of a recombinant inbred line (RIL) population and the limited genotypic information of its parental lines. Two ABH genotype files were created from a pseudo-parental genotypic data set that includes all the SNPs present in the RIL population. In the first ABH file pseudo-parental 1 was considered parental A, while in the second pseudo-parental 1 was considered parental B. These two duplicate ABH genotype files were merged by chromosome and subjected to linkage map analysis. Since the ABH data were duplicated, two mirrored linkage groups were generated per chromosome. The correct linkage map was identified and selected based on the partial genotypic data of the parental lines. This strategy was effective for constructing a highly saturated linkage map of 33,421 SNPs based on the genotyping of 205 RILs and a limited number of 100 SNPs present in the parental lines. This strategy enables the use of all the NGS SNP data obtained from a low-coverage sequencing experiment in the mapping population.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Fast and accurate metagenotyping of the human gut microbiome with GT-Pro

Single nucleotide polymorphisms (SNPs) in metagenomics are used to quantify population structure, track strains and identify genetic determinants of microbial phenotypes. However, existing alignment-based approaches for metagenomic SNP detection require high-performance computing and enough read coverage to distinguish SNPs from sequencing errors. To address these issues, we developed the GenoTyper for Prokaryotes (GT-Pro), a suite of methods to catalog SNPs from genomes and use unique k-mers to rapidly genotype these SNPs from metagenomes. Compared to methods that use read alignment, GT-Pro is more accurate and two orders of magnitude faster. Here, using high-quality genomes, we constructed a catalog of 104 million SNPs in 909 human gut species and used unique k-mers targeting this catalog to characterize the global population structure of gut microbes from 7,459 samples. GT-Pro enables fast and memory-efficient metagenotyping of millions of SNPs on a personal computer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Particle dynamics of nanoplastics suspended in water with soil microparticles: insights from small angle neutron scattering (SANS) and ultra-SANS

Small-angle neutron scattering (SANS) and ultra-SANS (USANS) were employed to understand the aggregation behavior and observe the size reduction for nanoplastics (NPs) formed from a biodegradable mulch film, and microparticles of vermiculite (V), an artificial soil, suspended in water in the presence of low convective shear (ex situ stirring) prior to measurements. Neutron contrast matching was employed to minimize the signal of V (by 100-fold) and thereby isolate the signal due to NPs in the neutron beam, as the contrast match point (CMP) for V (67 vol% deuteration of water) differed from that of NPs by more than 20%. The original NPs' size distribution was bimodal: <200 nm and 500–1200 nm, referred to as small and large NPs, i.e., SNPs and LNPs, respectively. In the absence of V, SNPs formed homoaggregates at higher concentrations that decreased with stirring time, while the size of LNPs remained unchanged. The presence of V at 2-fold lower concentration than NPs did not change the size of SNPs but reduced the size of LNPs by nearly 2-fold as stirring time increased. Because the size of SNPs and LNPs did not differ substantially between CMP and 100% D 2 O solvents, it is evident that SNPs and LNPs are mainly composed of NPs and not V. In conclusion, the results suggest that LNPs are susceptible to size reduction through collisions with soil microparticles via convection, yielding SNPs near soil–water interfaces within vadose zones.

Nanoplastics↗
DOE OSTI2022

Genome‐wide association studies reveal genetic basis of ionomic variation in cassava

SUMMARY As one of the most important food crops, cassava ( Manihot esculenta ) is the main dietary source of micronutrients for about 1 billion people. However, the ionomic variation in cassava and the underlying genetic mechanisms remain unclear so far. Herein, genome‐wide association studies were performed to reveal the specific single nucleotide polymorphisms (SNPs) that affect the ionomic variation in cassava. We identified 164 SNPs with P‐ values lower than the threshold located in 88 loci associated with divergent ionomic variations. Among them, 13 SNPs are related to both calcium (Ca) and magnesium (Mg), and many loci for different ionomic traits seem to be clustered on specific chromosome regions. Moreover, we identified the peak SNPs in the promoter regions of Sc10g003170 (encoding methionyl‐tRNA synthetase [MetRS]) and Sc18g015190 (encoding the transcriptional regulatory protein AlgP) for nitrogen (N) and phosphorus (P) accumulation, respectively. Notably, these two SNPs (chr10_32807962 and chr18_31343738) were directly correlated with the transcript levels of Sc10g003170 ( MetRS ) and Sc18g015190 ( AlgP ), which positively modulated N accumulation and P concentration in cassava, respectively. Taken together, this study provides important insight into the genetic basis of cassava natural ionomic variation, which will promote genetic breeding to improve nutrient use and accumulation of elements in cassava.

Yin, Hongyan↗
DOE OSTI2023

Genome‐wide association and genomic prediction for yield and component traits of Miscanthus sacchariflorus

Abstract Accelerating biomass improvement is a major goal of Miscanthus breeding. The development and implementation of genomic‐enabled breeding tools, like marker‐assisted selection (MAS) and genomic selection, has the potential to improve the efficiency of Miscanthus breeding. The present study conducted genome‐wide association (GWA) and genomic prediction of biomass yield and 14 yield‐components traits in Miscanthus sacchariflorus . We evaluated a diversity panel with 590 accessions of M. sacchariflorus grown across 4 years in one subtropical and three temperate locations and genotyped with 268,109 single‐nucleotide polymorphisms (SNPs). The GWA study identified a total of 835 significant SNPs and 674 candidate genes across all traits and locations. Of the significant SNPs identified, 280 were localized in mapped quantitative trait loci intervals and proximal to SNPs identified for similar traits in previously reported Miscanthus studies, providing additional support for the importance of these genomic regions for biomass yield. Our study gave insights into the genetic basis for yield‐component traits in M. sacchariflorus that may facilitate marker‐assisted breeding for biomass yield. Genomic prediction accuracy for the yield‐related traits ranged from 0.15 to 0.52 across all locations and genetic groups. Prediction accuracies within the six genetic groupings of M. sacchariflorus were limited due to low sample sizes. Nevertheless, the Korea/NE China/Russia ( N = 237) genetic group had the highest prediction accuracy of all genetic groups (ranging 0.26–0.71), suggesting that with adequate sample sizes, there is strong potential for genomic selection within the genetic groupings of M. sacchariflorus . This study indicated that MAS and genomic prediction will likely be beneficial for conducting population‐improvement of M. sacchariflorus .

09 BIOMASS FUELS↗
DOE OSTI2023

Data for Genome-Wide Association and Genomic Prediction for Yield and Component Traits of Miscanthus sacchariflorus

Accelerating biomass improvement is a major goal of miscanthus breeding. The development and implementation of genomic-enabled breeding tools, like marker-assisted selection (MAS) and genomic selection, has the potential to improve the efficiency of miscanthus breeding. The present study conducted genome-wide association (GWA) and genomic prediction of biomass yield and 14 yield-components traits in Miscanthus sacchariflorus . We evaluated a diversity panel with 590 accessions of M. sacchariflorus grown across four years in one subtropical and three temperate locations and genotyped with 268,109 single-nucleotide polymorphisms (SNPs). The GWA study identified a total of 835 significant SNPs and 674 candidate genes across all traits and locations. Of the significant SNPs identified, 280 were localized in mapped quantitative trait loci intervals and proximal to SNPs identified for similar traits in previously reported miscanthus studies, providing additional support for the importance of these genomic regions for biomass yield. Our study gave insights into the genetic basis for yield-component traits in M. sacchariflorus that may facilitate marker-assisted breeding for biomass yield. Genomic prediction accuracy for the yield-related traits ranged from 0.15 to 0.52 across all locations and genetic groups. Prediction accuracies within the six genetic groupings of M. sacchariflorus were limited due to low sample sizes. Nevertheless, the Korea/NE China/Russia (N = 237) genetic group had the highest prediction accuracy of all genetic groups (ranging 0.26–0.71), suggesting that with adequate sample sizes, there is strong potential for genomic selection within the genetic groupings of M. sacchariflorus . This study indicated that MAS and genomic prediction will likely be beneficial for conducting population-improvement of M. sacchariflorus .

Biomass Analytics↗
DOE OSTI2020

Novel Mode of Molybdate Inhibition of Desulfovibrio vulgaris Hildenborough

Sulfate-reducing microorganisms (SRM) are found in multiple environments and play a major role in global carbon and sulfur cycling. Because of their growth capabilities and association with metal corrosion, controlling the growth of SRM has become of increased interest. One such mechanism of control has been the use of molybdate (MoO 4 2− ), which is thought to be a specific inhibitor of SRM. The way in which molybdate inhibits the growth of SRM has been enigmatic. It has been reported that molybdate is involved in a futile energy cycle with the sulfate-activating enzyme, sulfate adenylyl transferase (Sat), which results in loss of cellular ATP. However, we show here that a deletion of this enzyme in the model SRM, Desulfovibrio vulgaris Hildenborough, remained sensitive to molybdate. We performed several subcultures of the ∆ sat strain in the presence of increasing concentrations of molybdate and obtained a culture with increased resistance to the inhibitor (up to 3 mM). The culture was re-sequenced and three single nucleotide polymorphisms (SNPs) were identified that were not present in the parental strain. Two of the SNPs seemed unlikely candidates for molybdate resistance due to a lack of conservation of the mutated residues in homologous genes of closely related strains. The remaining SNP was located in DVU2210, a protein containing two domains: a YcaO-like domain and a tetratricopeptide-repeat domain. The SNP resulted in a change of a serine residue to arginine in the ATP-hydrolyzing motif of the YcaO-like domain. Deletion mutants of each of the three genes apparently enriched with SNPs in the presence of inhibitory molybdate and combinations of these genes were generated in the Δ sat and wild-type strains. Strains lacking both sat and DVU2210 became more resistant to molybdate. Deletions of the other two genes in which SNPs were observed did not result in increased resistance to molybdate. YcaO-like proteins are distributed across the bacterial and archaeal domains, though the function of these proteins is largely unknown. The role of this protein in D. vulgaris is unknown. Due to the distribution of YcaO-like proteins in prokaryotes, the veracity of molybdate as a specific SRM inhibitor should be reconsidered.

59 BASIC BIOLOGICAL SCIENCES↗
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↗
DOE OSTI2022

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress ( Thlaspi arvense )

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A sorghum practical haplotype graph facilitates genome‐wide imputation and cost‐effective genomic prediction

Abstract Successful management and utilization of increasingly large genomic datasets is essential for breeding programs to accelerate cultivar development. To help with this, we developed a Sorghum bicolor Practical Haplotype Graph (PHG) pangenome database that stores haplotypes and variant information. We developed two PHGs in sorghum that were used to identify genome‐wide variants for 24 founders of the Chibas sorghum breeding program from 0.01x sequence coverage. The PHG called single nucleotide polymorphisms (SNPs) with 5.9% error at 0.01x coverage—only 3% higher than PHG error when calling SNPs from 8x coverage sequence. Additionally, 207 progenies from the Chibas genomic selection (GS) training population were sequenced and processed through the PHG. Missing genotypes were imputed from PHG parental haplotypes and used for genomic prediction. Mean prediction accuracies with PHG SNP calls range from .57–.73 and are similar to prediction accuracies obtained with genotyping‐by‐sequencing or targeted amplicon sequencing (rhAmpSeq) markers. This study demonstrates the use of a sorghum PHG to impute SNPs from low‐coverage sequence data and shows that the PHG can unify genotype calls across multiple sequencing platforms. By reducing input sequence requirements, the PHG can decrease the cost of genotyping, make GS more feasible, and facilitate larger breeding populations. Our results demonstrate that the PHG is a useful research and breeding tool that maintains variant information from a diverse group of taxa, stores sequence data in a condensed but readily accessible format, unifies genotypes across genotyping platforms, and provides a cost‐effective option for genomic selection.

Jensen, Sarah E.↗
DOE OSTI2024

Genome‐wide association studies on resistance to powdery mildew in cultivated emmer wheat

Abstract Powdery mildew, caused by the fungal pathogenBlumeria graminis(DC.) E. O. Speer f. sp.triticiEm. Marchal (Bgt), is a constant threat to global wheat (Triticum aestivumL.) production. Although ∼100 powdery mildew (Pm) resistance genes and alleles have been identified in wheat and its relatives, more is needed to minimizeBgt’s fast evolving virulence. In tetraploid wheat (Triticum turgidumL.), wild emmer wheat [T. turgidumssp.dicoccoides(Körn. ex Asch. & Graebn.) Thell.] accessions from Israel have contributed manyPmresistance genes. However, the diverse genetic reservoirs of cultivated emmer wheat [T. turgidumssp.dicoccum(Schrank ex Schübl.) Thell.] have not been fully exploited. In the present study, we evaluated a diverse panel of 174 cultivated emmer accessions for their reaction toBgtisolateOKS(14)‐B‐3‐1and found that 66% of accessions, particularly those of Ethiopian (30.5%) and Indian (6.3%) origins, exhibited high resistance. To determine the genetic basis ofBgtresistance in the panel, genome‐wide association studies were performed using 46,383 single nucleotide polymorphisms (SNPs) from genotype‐by‐sequencing and 4331 SNPs from the 9K SNP Infinium array. Twenty‐five significant SNP markers were identified to be associated withBgtresistance, of which 21 SNPs are likely novel loci, whereas four possibly represent emmer derivedPm4a,Pm5a,PmG16, andPm64. Most novel loci exhibited minor effects, whereas three novel loci on chromosome arms 2AS, 3BS, and 5AL had major effect on the phenotypic variance. This study demonstrates cultivated emmer as a rich source of powdery mildew resistance, and the resistant accessions and novel loci found herein can be utilized in wheat breeding programs to enhanceBgtresistance in wheat.

Genetics & Heredity↗