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

Genome-wide association study of long COVID

Abstract Infections can lead to persistent symptoms and diseases such as shingles after varicella zoster or rheumatic fever after streptococcal infections. Similarly, severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) infection can result in long coronavirus disease (COVID), typically manifesting as fatigue, pulmonary symptoms and cognitive dysfunction. The biological mechanisms behind long COVID remain unclear. We performed a genome-wide association study for long COVID including up to 6,450 long COVID cases and 1,093,995 population controls from 24 studies across 16 countries. We discovered an association ofFOXP4with long COVID, independent of its previously identified association with severe COVID-19. The signal was replicated in 9,500 long COVID cases and 798,835 population controls. Given the transcription factor FOXP4’s role in lung physiology and pathology, our findings highlight the importance of lung function in the pathophysiology of long COVID.

Genetics & Heredity↗
DOE OSTI2025

Combining genome-wide association studies and expression quantitative trait nucleotide mapping with molecular and genetic validations to identify transcriptional networks regulating drought tolerance in Populus

Objectives: (i). To deploy a large-scale experimental drought trial for up to 1000 unique genotypes of Populus equipping the sites with controlled irrigation and drought treatments that are fully automated and monitored. FULLY COMPLETED (ii) To test the hypothesis that a suite of traits identified for drought tolerance in P. nigra can be measured in drought and control treatments in the wide germplasm collection of P. trichocarpa. FULLY COMPLETED (iii) To use established and novel GWAS model approaches to identify gene loci linked to drought tolerance traits on interest in P. trichocarpa. FULLY COMPLETED (iv) To undertake comparative analysis of GWAS results for drought tolerance traits in P. nigra and P. trichocarpa. PARTIALLY COMPLETED – remains active (v) Using RNAseq in P. trichocarpa, in droughted and control treatments to identify cis- and trans-regulated eQTN. FULLY COMPLETED (vi) Validate up to 50 cis-QTNs, from network hubs using transient protoplast assays. FULLY COMPLETED (vii) To establish Agrobacterium-based gene editing protocols in Populus. FULLY COMPLETED (viii) To utilize early leads from previous research to investigate at least 6 candidate genes for drought tolerance in Populus. FULLY COMPLETED (ix) To validate up to 20 candidate genes for drought tolerance in P. trichocarpa refined from the long-list tested in the transient assays for cis-acting hub gene targets. PARTIALLY COMPLETED- remains active.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

SAIGE-GPU: accelerating genome- and phenome-wide association studies using GPUs

Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.

Rodriguez, Alex [Argonne National Laboratory (ANL)↗
DOE OSTI2025

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

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

Science & Technology - Other Topics↗
DOE OSTI2026

Genetic variations and their interaction with thirdhand smoke exposure on anxiety and memory in Collaborative Cross mice

Thirdhand smoke (THS) is linked to adverse health effects, but the effect of genetic variations on behavioral outcomes is poorly understood. To investigate this, we assessed anxiety- and memory-related behaviors in 820 mice from 21 strains of the genetically diverse Collaborative Cross (CC) mouse that were exposed to THS from 4 through 10 weeks of age. Anxiety was evaluated with a light/dark box assay with a previously established risk score system. Females were generally more sensitive: THS reduced anxiety risk in strains CC013, CC019, and CC051, but increased risk in CC036 and CC061, while males showed no significant effects. Memory was tested using passive avoidance: impairments were observed in both sexes in CC016 and CC019, with sex-dependent effects in CC002 and CC051. A genome-wide association study identified 2,347 SNPs associated with anxiety and 1,568 SNPs with memory, with 32 and 85 SNPs, respectively, interacting with THS exposure. Enrichment analyses revealed distinct biological processes underlying susceptibility, including axonogenesis, synapse organization, cognition, and learning and memory. KEGG pathway analysis identified distinct genetic pathways, including GTPase binding and GTPase regulatory activity, that act as critical molecular switches in the brain that regulate synaptic plasticity, dendritic spine structure, and neuronal signaling, directly influencing anxiety-like behaviors and memory formation. These findings show that THS exposure affects neurobehavioral outcomes in a sex- and genotype-dependent manner, highlighting critical gene-environment interactions and providing a foundation for mechanistic insights into THS neurotoxicity

Anxiety↗
DOE OSTI2026

Constitutive and inducible oleoresin defenses share genetic architectures and mechanisms in Pinus taeda

The oleoresin defense system of loblolly pine (Pinus taeda) protects trees from insects and pathogens and is an important source of renewable biofuels and chemicals, but the genetic basis of oleoresin production is poorly understood. We characterized the genetic architecture of oleoresin flow, resin canal number, stem wood terpene content, and monoterpene composition in two clonal populations of P. taeda. We used quantitative genetic analyses, genome-wide association studies (GWASs), multiplex network learning, and gene expression profiling to elucidate shared gene networks underlying defense traits and to identify high-quality candidates for breeding and engineering loblolly pine. Genetic analyses revealed polygenic inheritance and trait-to-trait correlations provide strong evidence for shared genes regulating constitutive and induced oleoresin flow. We identified 236 single nucleotide polymorphisms associated with oleoresin flow, resin canal number, and terpene composition and highlight candidate genes likely involved in terpene biosynthesis, cambial meristem reprogramming, and pathogen perception and immune signaling. Fourteen GWAS candidates were methyl jasmonate-responsive in tissues where resin canals initiate and terpene production occurs. Integrating quantitative genetics, GWAS, gene expression, and multiplex network analyses enabled the prioritization of high-quality candidate genes. This work advances the development of more resilient loblolly pine optimized for ecological performance, renewable chemical, and biofuel production.

genome-wide association study↗
DOE OSTI2026

Unraveling plant phenotype to genotype associations with daily hyperspectral traits in Populus trichocarpa

Hyperspectral remote sensing is a powerful, high-throughput phenotyping tool that quantifies physiologically and structurally relevant wavelengths across diverse genotypes and over varying temporal scales. In this study, we combined tower-based continuous hyperspectral sensing with genome-wide association studies to analyze 1423 wavebands (400-900 nm) and derivative vegetation indices across 505 genotypes and the genetic architecture of hyperspectral phenotypes over time in Populus trichocarpa Torr. & Gray grown under field conditions. Wavelengths related to chlorophyll and carotenoid absorption spectra exhibited the strongest genetic variation resulting in 98 significant SNP associations. Notably, we found substantial overlap in genetic association between the blue and red spectral regions, indicative of carotenoids and chlorophyll, respectively, and identified more than 10 candidate genes associated with chloroplast function, underpinning photosynthetic activity. Furthermore, fluctuations in associations for vegetative indices, such as the chlorophyll:carotenoid index (CCI), across the growing season reveal a temporally dynamic genetic architecture of physiological traits associated with fall senescence of this temperate tree species. Finally, we also observed correlations (spearman rho = 0.3, p < 1x10 −8 ) between individual wavebands or vegetative indices and growth rate, assessed as the relative change of tree height over the growing season. The growth rate prediction was substantially improved by a regularization multivariate model (spearman rho>0.5, p < 1x10 −16 ), reinforcing the value of hyperspectral measurements for predicting traits linked to tree productivity. These findings highlight the potential of high-throughput, rapid, hyperspectral genome wide association studies GWAS to uncover physiologically meaningful genetic variation and offer promising insights for future acceleration for plant breeding.

09 BIOMASS FUELS↗
DOE OSTI2025

Transcripts and genomic intervals associated with variation in metabolite abundance in maize leaves under field conditions

Abstract Plants exhibit extensive environment-dependent intraspecific metabolic variation, which likely plays a role in determining variation in whole plant phenotypes. However, much of the work seeking to use natural variation to link genes and transcript’s impacts on plant metabolism has employed data from controlled environments. Here, we generated and analyzed data on the variation in the abundance of 26 metabolites across 660 maize inbred lines under field conditions. We employ these data and previously published transcript and whole plant phenotype data reported for the same field experiment to identify both genomic intervals (through genome-wide association studies (GWAS)) and transcripts (using both transcriptome-wide association studies (TWAS) and an explainable artificial intelligence (AI) approach based on random forest (RF)) associated with variation in metabolite abundance. Both genome-wide association and random forest-based methods identified substantial numbers of significant associations including genes with plausible links to the metabolites they are associated with. In contrast, the transcriptome-wide association identified only six significant associations. In three cases, genetic markers associated with metabolic variation in our study colocalized with markers linked to variation in non-metabolic traits scored in the same experiment. We speculate that the poor performance of transcriptome-wide association studies in identifying transcript-metabolite associations may reflect a high prevalence of non-linear interactions between transcripts and metabolites and/or a bias towards rare transcripts playing a large role in determining intraspecific metabolic variation.

Mathivanan, Ramesh Kanna↗
DOE OSTI2025

Elucidating the drought-responsive changes in Poplar cuticular waxes: A GWAS analysis of genes involved in fatty acid biosynthesis

Drought and episodic drought events are major impending impacts of climate change, limiting the productivity of plants and especially trees due to their inherent high transpiration rates. One common mechanism used by plants to cope with drought stress is to change the composition of their leaf cuticular waxes. Cuticular waxes are essential for controlling non-stomatal water loss and are typically composed of a homologous series of very-long-chain fatty acid-derived compounds, as well as flavonoids, tocopherols, triterpenoids, and phytosterols. In this study, we compared the cuticular waxes of 339 natural accessions of Populus trichocarpa (black cottonwood) grown under control and drought conditions in a common garden. A Genome-Wide Association Study (GWAS) was then used to identify candidate genes associated with cuticular wax biosynthesis and/or its regulation. Although no major differences were observed in total wax load when subject to drought conditions, the amounts of the individual wax constituents were indeed responsive to drought. Specifically, changes in alkenes, alcohols, esters, and aldehydes were evident, and suggest that they contribute to the drought response/tolerance in poplar. GWAS uncovered several genes linked to fatty acid biosynthesis, including CER1, CER3, CER4, FATB, FAB1, FAR3, FAR4, KCS, and a homolog of SOH1, as well as other candidate genes that may be involved in coordinating the drought responses in poplar trees. Our findings provide new evidence that genotype-specific shifts in wax composition, rather than total wax accumulation, contribute to drought adaptation in poplar. Additionally, we show that genetic variation in key wax biosynthetic genes drives cuticular wax plasticity in P. trichocarpa under drought, identifying putative molecular targets for improving stress resilience in trees. This study expands our understanding of the adaptative mechanisms of the cuticle and their potential for enhancing drought tolerance in poplar species.

Alkanes↗
DOE OSTI2025

A multifunctional sesquiterpene synthase integrates with cytochrome P450s to reinforce the terpenoid defense network in maize

Terpenoids, the largest and most structurally diverse class of plant natural products, play essential roles in maize defense and ecological interactions. In this study, we identified and functionally characterized a sesquiterpenoid-based defense pathway in maize centered on α-santalenoic acid, a pathogen-inducible sesquiterpenoid antibiotic. Using a combination of metabolite-based genome-wide association studies (mGWAS), linkage mapping, and heterologous expression assays, we identified ZmTPS9 as a multiproduct terpene synthase that primarily produces α-santalene and β-bisabolene. Sequence analysis and site-directed mutagenesis revealed that threonine at position 413 is critical for enzyme activity, with its deletion resulting in a complete loss of enzyme activity. The sesquiterpene hydrocarbons produced by ZmTPS9 are further oxidized by three cytochrome P450 monooxygenases, ZmCYP71Z16, ZmCYP71Z18, and ZmCYP71Z19, to yield antimicrobial metabolites including α-santalenoic acid, zealexin D1 (ZD1), and zealexin D2 (ZD2). Together, these findings demonstrate a convergent biosynthetic strategy in maize, where multiproduct terpene synthases and promiscuous P450s collaboratively generate a flexible and robust terpenoid defense network.

a-santalenoic acid↗
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 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 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↗
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↗