Search NASA⌕ Search

SEARCH · Search NASA

Results for “Single nucleotide polymorphisms”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Development and validation of a quantitative PCR assay method of assessing relative resistance of oat ( Avena sativa ) to crown rust ( Puccinia coronata f. sp. avenae )

Evaluation of oat crown rust resistance is usually based on visual assessment of disease severity or infection types. Visual assessment is subjective, prone to rater bias and requires expert knowledge. PCR‐based quantitative assays can overcome challenges associated with visual assessment. New TaqMan primers and probes were designed from Puccinia coronata f. sp. avenae (Pca) sequences. The primer–probe sets were specific to Pca, amplified using as little as 0.5 pg fungal DNA (fDNA) and allowed for scaling to variation in sample total DNA quantity. The quantitative PCR (qPCR) assay was validated using oat recombinant inbred lines (RILs) from the Provena × 94197A1‐9‐2‐2‐2‐5 cross evaluated under a controlled environment. For comparison with fDNA load, inoculation with the Pca race LCBB provided segregation data on the hypersensitive response, while Pca race LSLG provided data on segregation for reduced pustule number. fDNA content was positively correlated with both pustule number and infection type (IT). Composite interval mapping identified two quantitative trait loci (QTLs) on oat linkage groups Mrg12 and Mrg20 using visual and qPCR assessments (pustule number, IT and fDNA). In this study a qPCR assay method that can be used to assess the relative resistance of oat to crown rust was refined and validated, and single nucleotide polymorphisms (SNPs) closely linked with two QTLs derived from the crown rust resistant line 94197A1‐9‐2‐2‐2‐5 were identified.

Admassu‐Yimer, B.↗
DOE OSTI2021

The soybean β‐expansin gene GmINS1 contributes to nodule development in response to phosphate starvation

Abstract Legume biological nitrogen fixation (BNF) is the most important N source in agricultural ecosystems. Nodule organogenesis from the primordia to the development of mature nodules with the ability to fix N 2 largely determines BNF capacity. However, nodule growth is often limited by low phosphorus (P) availability, while the mechanisms underlying nodule development responses to P deficiency remain largely unknown. In this study, we found that nodule enlargement is severely inhibited by P deficiency, as reflected by the smaller individual nodule size from a soybean core collection in the field. Wide‐ranging natural diversity in nodule size was further identified in soybeans reared in low P soils, with the FC‐1 genotype outperforming FC‐2 in assessments of nodulation under low P conditions. Among β‐expansin members, GmINS1 expression is most abundantly enhanced by P deficiency in FC‐1 nodules, and its transcript level is further displayed to be tightly associated with nodule enlargement. Four single nucleotide polymorphisms discovered in the GmINS1 promoter distinguished the FC‐1 and FC‐2 genotypes and accounted for the differential expression levels of GmINS1 responses to P deficiency. GmINS1 overexpression led to increases in nodule size, infected cell abundance, and N 2 fixation capacity and subsequently promoted increases in N and P content, soybean biomass, and yield. Our findings provide a candidate gene for optimizing BNF capacity responses to low P stress in soybean molecular breeding programs.

Yang, Zhaojun↗
DOE OSTI2015

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high-coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single-nucleotide polymorphisms (SNPs), of which 1 590 653 high-confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype-restricted, non-synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near-reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up-CNVs), 112 630 genes present in fewer copies than in the reference genome (down-CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland-specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.

09 BIOMASS FUELS↗
DOE OSTI2019

Construction and comparison of three reference-quality genome assemblies for soybean

Here, we report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

SNPeffect: identifying functional roles of SNPs using metabolic networks

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Alternative splicing of OsGS1;1 affects nitrogen‐use efficiency, grain development, and amylose content in rice

SUMMARY Excessive nitrogen fertilizer application is harmful to the environment and reduces the quality of cereal crops. Maintaining crop yields under low nitrogen (LN) conditions and improving quality are important goals for cereal crop breeding. Although the effects of nitrogen assimilation on crop nitrogen‐use efficiency (NUE) have been intensively studied, natural variations of the key assimilation genes underlying grain development and quality are largely unclear. Here, we identified an NUE‐associated gene, OsGS1;1 , encoding glutamine synthase, through genome‐wide association analysis, followed by validation experiments and functional analysis. Fifteen single‐nucleotide polymorphisms in the OsGS1;1 region led to alternative splicing that generated two functional transcripts: OsGS1;1a and OsGS1;1b . The elite haplotype of OsGS1;1 showed high OsGS1;1b activity, which improved NUE, affected grain development, and reduced amylose content. The results show that OsGS1;1 , which is induced under LN conditions, affects grain formation by regulating sugar metabolism and may provide a new avenue for the breeding of high‐yield and high‐quality rice ( Oryza sativa ).

Liu, Xiaolan↗
DOE OSTI2022

Sorghum Association Panel whole‐genome sequencing establishes cornerstone resource for dissecting genomic diversity

SUMMARY Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions. We report the whole‐genome sequencing (WGS) of 400 sorghum ( Sorghum bicolor (L.) Moench) accessions from the Sorghum Association Panel (SAP) at an average coverage of 38× (25–72×), enabling the development of a high‐density genomic marker set of 43 983 694 variants including single‐nucleotide polymorphisms (approximately 38 million), insertions/deletions (indels) (approximately 5 million), and copy number variants (CNVs) (approximately 170 000). We observe slightly more deletions among indels and a much higher prevalence of deletions among CNVs compared to insertions. This new marker set enabled the identification of several novel putative genomic associations for plant height and tannin content, which were not identified when using previous lower‐density marker sets. WGS identified and scored variants in 5‐kb bins where available genotyping‐by‐sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six F st peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has served and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

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

SU-D-204-06: Integration of Machine Learning and Bioinformatics Methods to Analyze Genome-Wide Association Study Data for Rectal Bleeding and Erectile Dysfunction Following Radiotherapy in Prostate Cancer

Purpose: We investigated whether integration of machine learning and bioinformatics techniques on genome-wide association study (GWAS) data can improve the performance of predictive models in predicting the risk of developing radiation-induced late rectal bleeding and erectile dysfunction in prostate cancer patients. Methods: We analyzed a GWAS dataset generated from 385 prostate cancer patients treated with radiotherapy. Using genotype information from these patients, we designed a machine learning-based predictive model of late radiation-induced toxicities: rectal bleeding and erectile dysfunction. The model building process was performed using 2/3 of samples (training) and the predictive model was tested with 1/3 of samples (validation). To identify important single nucleotide polymorphisms (SNPs), we computed the SNP importance score, resulting from our random forest regression model. We performed gene ontology (GO) enrichment analysis for nearby genes of the important SNPs. Results: After univariate analysis on the training dataset, we filtered out many SNPs with p>0.001, resulting in 749 and 367 SNPs that were used in the model building process for rectal bleeding and erectile dysfunction, respectively. On the validation dataset, our random forest regression model achieved the area under the curve (AUC)=0.70 and 0.62 for rectal bleeding and erectile dysfunction, respectively. We performed GO enrichment analysis for the top 25%, 50%, 75%, and 100% SNPs out of the select SNPs in the univariate analysis. When we used the top 50% SNPs, more plausible biological processes were obtained for both toxicities. An additional test with the top 50% SNPs improved predictive power with AUC=0.71 and 0.65 for rectal bleeding and erectile dysfunction. A better performance was achieved with AUC=0.67 when age and androgen deprivation therapy were added to the model for erectile dysfunction. Conclusion: Our approach that combines machine learning and bioinformatics techniques enabled designing better models and identifying more plausible biological processes associated with the outcomes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-207B-02: A Bayesian Network Approach for Joint Prediction of Tumor Control and Radiation Pneumonitis (RP) in Non-Small-Cell Lung Cancer (NSCLC)

Purpose: NSCLC radiotherapy treatment is a trade-off between controlling the tumor while limiting radiation-induced toxicities. Here we identify hierarchical biophysical relationships that could simultaneously influence both local control (LC) and RP by using an integrated Bayesian Networks (BN) approach. Methods: We studied 79 NSCLC patients treated on prospective protocol with 56 cases of LC and 21 events of RP. Beyond dosimetric information, each patient had 193 features including 12 clinical factors, 60 circulating blood cytokines before and during radiotherapy, 62 microRNAs, and 59 single-nucleotide polymorphisms (SNPs). The most relevant biophysical predictors for both LC and RP were identified using a Markov blanket local discovery algorithm and the corresponding BN was constructed using a score-learning algorithm. The area under the free-response receiver operating characteristics (AU-FROC) was used for performance evaluation. Cross-validation was employed to guard against overfitting pitfalls. Results: A BN revealing the biophysical interrelationships jointly in terms of LC and RP was developed and evaluated. The integrated BN included two SNPs, one microRNA, one clinical factor, three pre-treatment cytokines, relative changes of two cytokines between pre and during-treatment, and gEUDs of the GTV (a=-20) and lung (a=1). On cross-validation, the AUC prediction of independent LC was 0.85 (95% CI: 0.75–0.95) and RP was 0.83 (0.73–0.92). The AU-FROC of the integrated BN to predict both LC/RP was 0.81 (0.71–0.90) based on 2000 stratified bootstrap, indicating minimal loss in joint prediction power. Conclusions: We developed a new approach for multiple outcome utility application in radiotherapy based on integrated BN techniques. The BN developed from large-scale retrospective data is able to simultaneously predict LC and RP in NSCLC treatments based on individual patient characteristics. The joint prediction is only slightly compromised compared to independent predictions. Our approach shows promise for use in clinical decision support system for personalized radiotherapy subject to multiple endpoints. These studies were supported by a grant from the NCI/NIH P01-CA59827.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Listeria monocytogenes in Stone Fruits Linked to a Multistate Outbreak: Enumeration of Cells and Whole-Genome Sequencing

In 2014, the identification of stone fruits contaminated with led to the subsequent identification of a multistate outbreak. Simultaneous detection and enumeration of L. monocytogenes were performed on 105 fruits, each weighing 127 to 145 g, collected from 7 contaminated lots. The results showed that 53.3% of the fruits yielded L. monocytogenes (lower limit of detection, 5 CFU/fruit), and the levels ranged from 5 to 2,850 CFU/fruit, with a geometric mean of 11.3 CFU/fruit (0.1 CFU/g of fruit). Two serotypes, IVb-v1 and 1/2b, were identified by a combination of PCR- and antiserum-based serotyping among isolates from fruits and their packing environment; certain fruits contained a mixture of both serotypes. Single nucleotide polymorphism (SNP)-based whole-genome sequencing (WGS) analysis clustered isolates from two case-patients with the serotype IVb-v1 isolates and distinguished outbreak-associated isolates from pulsed-field gel electrophoresis (PFGE)-matched, but epidemiologically unrelated, clinical isolates. The outbreak-associated isolates differed by up to 42 SNPs. All but one serotype 1/2b isolate formed another WGS cluster and differed by up to 17 SNPs. Fully closed genomes of isolates from the stone fruits were used as references to maximize the resolution and to increase our confidence in prophage analysis. Putative prophages were conserved among isolates of each WGS cluster. All serotype IVb-v1 isolates belonged to singleton sequence type 382 (ST382); all but one serotype 1/2b isolate belonged to clonal complex 5.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Laboratory evolution in Novosphingobium aromaticivorans enables rapid catabolism of a model lignin-derived aromatic dimer

Lignin contains a variety of interunit linkages, leading to a range of potential decomposition products that can be used as carbon and energy sources by microbes. β-O-4 linkages are the most common in native lignin, and associated catabolic pathways have been well characterized. However, the fate of the mono-aromatic intermediates that result from β-O-4 dimer cleavage has not been fully elucidated. Here, we used experimental evolution to identify mutant strains of Novosphingobium aromaticivorans with improved catabolism of a model aromatic dimer containing a β-O-4 linkage, guaiacylglycerol-β-guaiacyl ether (GGE). We identified several parallel causal mutations, including a single nucleotide polymorphism in the promoter of an uncharacterized gene that roughly doubled the growth yield with GGE. We characterized the associated enzyme and demonstrated that it oxidizes an intermediate in GGE catabolism, β-hydroxypropiovanillone, to vanilloyl acetaldehyde. Identification of this enzyme and its key role in GGE catabolism furthers our understanding of catabolic pathways for lignin-derived aromatic compounds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

The Source and Evolutionary History of a Microbial Contaminant Identified Through Soil Metagenomic Analysis

ABSTRACT In this study, strain-resolved metagenomics was used to solve a mystery. A 6.4-Mbp complete closed genome was recovered from a soil metagenome and found to be astonishingly similar to that of Delftia acidovorans SPH-1, which was isolated in Germany a decade ago. It was suspected that this organism was not native to the soil sample because it lacked the diversity that is characteristic of other soil organisms; this suspicion was confirmed when PCR testing failed to detect the bacterium in the original soil samples. D. acidovorans was also identified in 16 previously published metagenomes from multiple environments, but detailed-scale single nucleotide polymorphism analysis grouped these into five distinct clades. All of the strains indicated as contaminants fell into one clade. Fragment length anomalies were identified in paired reads mapping to the contaminant clade genotypes only. This finding was used to establish that the DNA was present in specific size selection reagents used during sequencing. Ultimately, the source of the contaminant was identified as bacterial biofilms growing in tubing. On the basis of direct measurement of the rate of fixation of mutations across the period of time in which contamination was occurring, we estimated the time of separation of the contaminant strain from the genomically sequenced ancestral population within a factor of 2. This research serves as a case study of high-resolution microbial forensics and strain tracking accomplished through metagenomics-based comparative genomics. The specific case reported here is unusual in that the study was conducted in the background of a soil metagenome and the conclusions were confirmed by independent methods. IMPORTANCE It is often important to determine the source of a microbial strain. Examples include tracking a bacterium linked to a disease epidemic, contaminating the food supply, or used in bioterrorism. Strain identification and tracking are generally approached by using cultivation-based or relatively nonspecific gene fingerprinting methods. Genomic methods have the ability to distinguish strains, but this approach typically has been restricted to isolates or relatively low-complexity communities. We demonstrate that strain-resolved metagenomics can be applied to extremely complex soil samples. We genotypically defined a soil-associated bacterium and identified it as a contaminant. By linking together snapshots of the bacterial genome over time, it was possible to estimate how long the contaminant had been diverging from a likely source population. The results are congruent with the derivation of the bacterium from a strain isolated in Germany and sequenced a decade ago and highlight the utility of metagenomics in strain tracking.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Application of quasimetagenomics methods to define microbial diversity and subtype Listeria monocytogenes in dairy and seafood production facilities

Microorganisms frequently colonize surfaces within food production facilities. Detection of Listeria monocytogenes in this setting relies on culture-dependent methods, but the complex dynamics of bacterial interactions within these environments and their impact on pathogen detection remain largely unexplored. To address this challenge, we applied both 16S rRNA and shotgun quasimetagenomic (enriched microbiome) sequencing of swab culture enrichments from five seafood and seven dairy production environments. Utilizing 16S rRNA amplicon sequencing, we observed variability between 355 samples taken from these 12 production facilities and a distinctive microbiome for each environment. With shotgun quasimetagenomic sequencing, we were able to assemble L. monocytogenes metagenome-assembled genomes (MAGs) from 28 of the 32 culture-positive samples. We compared these MAGs to their corresponding whole-genome sequencing assemblies, which resulted in two polyphyletic clades consisting of L. monocytogenes lineages I and II with 13,195 and 25,556 single-nucleotide polymorphism sites, respectively. The remaining four MAGs did not produce sufficient genome coverage. To understand and establish limits for pathogen detection and subtyping using shotgun quasimetagenomics, these same data sets were downsampled in slilico to produce a titration series of abundances of L. monocytogenes and analyzed. Pathogen detection was achieved for all downsampled data sets, even those with only 3× genome coverage. This study contributes to the understanding of microbial diversity within food production environments and presents insights into the level of genome coverage needed in a metagenome sequencing data set to detect, subtype, and source track a foodborne pathogen.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Recent Advances in Understanding the Role of Genomic and Epigenomic Factors in Noncommunicable Diseases

The NGS based whole genome and exome sequencing endeavors have identified a plethora of mutations in the coding regions of the human genome that are held responsible as the causative factors for various human noncommunicable diseases including cancer. These mutations encompass various nonsynonymous nucleotide substitutions causing missense, nonsense, and frameshift changes in the protein-coding genes. For these types of mutations it is easier to establish a causal link between these mutations and disease phenotypes as they tend to change the amino acids in the proteins resulting in their loss-of-function (of tumor suppressor genes) or gain-of-function (of oncogenes) which might be directly responsible for the disease such as cancer. Apart from these disease-causing mutations, vast majority (more than 80%) of single-nucleotide polymorphisms were distributed throughout the human genome, mostly in the noncoding regions such as introns and intergenic regions. Large body of evidence suggested that variants in the noncoding regions can be considered as genetic predisposition to many noncommunicable diseases including cancer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Exploring the Genetic Basis of Wild Boar ( Sus scrofa ) and Its Connection to Classical Swine Fever Spread

Classical swine fever (CSF) is the one of the most devastating contagious diseases in domestic swine and wild boar/pigs (Sus scrofa). Population genetics is often used to estimate animal dispersal and can also help evaluate host population connectivity, which is crucial for understanding pathogen dispersal. We surveyed genetic population structure of boars using MIG-seq analysis to clarify the geographic barriers that influence boar dispersal in north-central Japan and to demonstrate the relationship between the spread of CSF infection among boars and their population structure. We obtained 382 single-nucleotide polymorphisms from 348 wild boar samples, and the results of STRUCTURE analysis indicated that the highest ΔK value was at K = 2, followed by K = 4. Based on these results, it is evident that the Abukuma river, a major river within north-central Japan, does not act as a barrier to the gene flow of boars, but rather that human infrastructure hinders their dispersal. Further, according to the time series change in the capture site of CSF-infected wild boar and the sum of the probability of belonging to each of the four clades in individual CSF-infected wild boar, our results indicated that the genetic structure of boar populations was correlated with the outbreak pathway of CSF across our study region. Our study suggests that predictions of disease spread, especially for widely distributed host species, is challenging because of the risk of cryptic breaks and changes in wide range connectivity; however, understanding the genetic population structure of wild boar can be a useful tool for predicting the spread of CSF. We concluded that genetic analysis of host population structure may have the possibility to improve predictions of the future dynamics of disease spread.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

High-performance epistasis detection in quantitative trait GWAS

epiSNP is a program for identifying pairwise single nucleotide polymorphism (SNP) interactions (epistasis) in quantitative-trait genome-wide association studies (GWAS). A parallel MPI version (EPISNPmpi) was created in 2008 to address this computationally expensive analysis on large data sets with many quantitative traits and SNP markers. However, the falling cost of genotyping has led to an explosion of large-scale GWAS data sets that challenge EPISNPmpi’s ability to compute results in a reasonable amount of time. Therefore, we optimized epiSNP for modern multi-core and highly parallel many-core processors to efficiently handle these large data sets. This paper describes the serial optimizations, dynamic load balancing using MPI-3 RMA operations, and shared-memory parallelization with OpenMP to further enhance load balancing and allow execution on the Intel Xeon Phi coprocessor (MIC). For a large GWAS data set, our optimizations provided a 38.43× speedup over EPISNPmpi on 126 nodes using 2 MICs on TACC’s Stampede Supercomputer. We also describe a Coarray Fortran (CAF) version that demonstrates the suitability of PGAS languages for problems with this computational pattern. We show that the Coarray version performs competitively with the MPI version on the NERSC Edison Cray XC30 supercomputer. Finally, the performance benefits of hyper-threading for this application on Edison (average 1.35× speedup) are demonstrated.

97 MATHEMATICS AND COMPUTING↗