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 37 records · Page 2

Genome-wide approaches delineate the additive, epistatic, and pleiotropic nature of variants controlling fatty acid composition in peanut ( Arachis hypogaea L.)

Abstract The fatty acid composition of seed oil is a major determinant of the flavor, shelf-life, and nutritional quality of peanuts. Major QTLs controlling high oil content, high oleic content, and low linoleic content have been characterized in several seed oil crop species. Here, we employ genome-wide association approaches on a recently genotyped collection of 787 plant introduction accessions in the USDA peanut core collection, plus selected improved cultivars, to discover markers associated with the natural variation in fatty acid composition, and to explain the genetic control of fatty acid composition in seed oils. Overall, 251 single nucleotide polymorphisms (SNPs) had significant trait associations with the measured fatty acid components. Twelve SNPs were associated with two or three different traits. Of these loci with apparent pleiotropic effects, 10 were associated with both oleic (C18:1) and linoleic acid (C18:2) content at different positions in the genome. In all 10 cases, the favorable allele had an opposite effect—increasing and lowering the concentration, respectively, of oleic and linoleic acid. The other traits with pleiotropic variant control were palmitic (C16:0), behenic (C22:0), lignoceric (C24:0), gadoleic (C20:1), total saturated, and total unsaturated fatty acid content. One hundred (100) of the significantly associated SNPs were located within 1000 kbp of 55 genes with fatty acid biosynthesis functional annotations. These genes encoded, among others: ACCase carboxyl transferase subunits, and several fatty acid synthase II enzymes. With the exception of gadoleic (C20:1) and lignoceric (C24:0) acid content, which occur at relatively low abundance in cultivated peanuts, all traits had significant SNP interactions exceeding a stringent Bonferroni threshold (α = 1%). We detected 7682 pairwise SNP interactions affecting the relative abundance of fatty acid components in the seed oil. Of these, 627 SNP pairs had at least one SNP within 1000 kbp of a gene with fatty acid biosynthesis functional annotation. We evaluated 168 candidate genes underlying these SNP interactions. Functional enrichment and protein-to-protein interactions supported significant interactions (P-value < 1.0E−16) among the genes evaluated. These results show the complex nature of the biology and genes underlying the variation in seed oil fatty acid composition and contribute to an improved genotype-to-phenotype map for fatty acid variation in peanut seed oil.

Otyama, Paul I.↗

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

59 BASIC BIOLOGICAL SCIENCES↗

MINE: maximally informative next experiment—toward a new GWAS experimental design and methodology

Abstract The computational methodology of Genome Wide Association Studies (GWAS) currently has several limitations: (i) the number of observations (rows) on a quantitative trait tends to be smaller than the number of single nucleotide polymorphisms (SNPs) (columns) in the design matrix; (ii) each SNP is usually modeled separately, failing to acknowledge interaction between each other (ie epistasis); (iii) there is implicit linkage disequilibrium (LD) between neighboring SNPs due to their linkage. To overcome these issues, we developed a tool that uses ensemble methods to fit mixed linear models to GWAS data, and these ensemble methods include the development of a new experimental design approach in GWAS, which uses the resultant models and data to select the next informative experiment over time. This new adaptive and staged approach for GWAS experimental design was developed and tested in a 3 yr adaptive model-guided discovery experiment against a fixed classical design. In Sorghum bicolor a total of 79, 86, and 78 accessions were tested in years 1, 2, and 3, respectively out of 343 accessions available in the Bioenergy Association Panel (BAP) each identified for 232,303 SNPs, 1 every 2–3 kb in the genomes. We demonstrated the feasibility of MINE enacted with 8 people in the field per year over 3 yr vs in 1 large classical design enacted with 20 people in 1 yr. The MINE results for chromosomal regions identified controlling dry weight were confirmed against results from previous sorghum GWAS experiments and 1 large classical design for the BAP panel.

Genetics & Heredity↗

The success of particular Acinetobacter baumannii clones: accumulating resistance and virulence inside a sugary shield

Abstract Background In Portugal, carbapenem-resistant Acinetobacter baumannii (CRAB) has been associated with ST98, ST103 and ST208 (Oxford Scheme, Oxf) and a clone has usually been associated with a particular period of time. These clonal shifts were primarily explained by an increased antimicrobial resistance profile. Here we explore genomic and biochemical differences among these and more recent clones, which could further explain the diversity and evolution of this species. Methods A total of 116 CRAB isolates (2010–15), together with representatives of a previously described CRAB collection (4 isolates, 2001–06) were characterized by attenuated total reflection Fourier transform infrared spectroscopy (FTIR-ATR) and MLST. Representatives of different FTIR-ATR/MLST clusters were selected for WGS (n = 13), which allowed the in silico extraction of resistance and virulence genes, capsule locus and SNP analysis. Results A. baumannii clonal shifts of OXA-58-producing ST103Oxf (2001–04), OXA-40-producing ST98Oxf (2002–06), OXA-23-producing ST208Oxf (2006–10) and OXA-23-producing ST218Oxf (2010–15) were accompanied by an increase in AMR genes and virulence factors. FTIR-ATR clustering was congruent with sugar composition predicted from the capsular locus: a fucosamine cluster comprising ST98Oxf, ST103Oxf and a single ST218Oxf isolate; a pseudaminic acid cluster of ST208Oxf and ST1557Oxf isolates; and legionaminic acid, resembling the sialic acid from mammalian cells, in a cluster comprising ST218Oxf isolates. The whole-genome phylogenetic tree was congruent with MLST, with isolates presenting 5–28 938 SNPs. ST208Oxf and ST218Oxf presented ∼1900 SNPs while ST103Oxf and ST1557Oxf showed a greater number of SNPs (∼28 000). Conclusions Clonal shifts of CRAB were promoted, in our country, by consecutive virulence and AMR gene pool enlargement, together with features increasing pathogen–host adaptation. Worldwide dominance of ST218Oxf is supported by the combination of high AMR and virulence levels.

Silva, Liliana↗

Unveiling shared genetic regulators of plant architectural and biomass yield traits in the Sorghum Association Panel

Abstract Sorghum is emerging as an ideal genetic model for designing high-biomass bioenergy crops. Biomass yield, a complex trait influenced by various plant architectural characteristics, is typically regulated by numerous genes. This study aimed to dissect the genetic regulators underlying 14 plant architectural traits and 10 biomass yield traits in the Sorghum Association Panel across two growing seasons. We identified 321 associated loci through genome-wide association studies (GWAS), involving 234 264 single nucleotide polymorphisms (SNPs). These loci include genes with known associations to biomass traits, such as maturity, dwarfing (Dw), and leafbladeless1, as well as several uncharacterized loci not previously linked to these traits. We also identified 22 pleiotropic loci associated with variation in multiple phenotypes. Three of these loci, located on chromosomes 3 (S03_15463061), 6 (S06_42790178; Dw2), and 9 (S09_57005346; Dw1), exerted significant and consistent effects on multiple traits across both growing seasons. Additionally, we identified three genomic hotspots on chromosomes 6, 7, and 9, each containing multiple SNPs associated with variation in plant architecture and biomass yield traits. Chromosome-wise correlation analyses revealed multiple blocks of positively associated SNPs located near or within the same genomic regions. Finally, genome-wide correlation-based network analysis showed that loci associated with flowering, plant height, leaf traits, plant density, and tiller number per plant were highly interconnected with other genetic loci influencing plant architectural and biomass yield traits. The pyramiding of favorable alleles related to these traits holds promise for enhancing the future development of bioenergy sorghum crops.

Singh, Anuradha (ORCID:0000000197149095)↗

Genomic description of human clinical Aspergillus fumigatus isolates, California, 2020

Abstract Aspergillus fumigatus, an environmental mold, causes life-threatening infections. Studies on the phylogenetic structure of human clinical A. fumigatus isolates are limited. Here, we performed whole genome sequencing of 24 A. fumigatus isolates collected from 18 patients in U.S. healthcare facilities in California. Single-nucleotide polymorphism (SNP) differences between patient isolates ranged from 187 to 70 829 SNPs. For five patients with multiple isolates, we calculated the within-host diversities. Three patients had a within-host diversity that ranged from 4 to 10 SNPs and two patients ranged from 2 to 16 977 SNPs. Findings revealed highly diverse A. fumigatus strains among patients and two patterns of diversity for isolates that come from the same patient, low and extremely high diversity.

Misas, Elizabeth (ORCID:0000000162437716)↗

Genetic Identification of Lamprey Genera and Anadromous Ecotypes in Watersheds of the Northeastern Pacific Ocean

ABSTRACT Nonparasitic, nonmigratory Western Brook Lamprey (WBL; Lampetra ayresii ), and parasitic, anadromous Western River Lamprey (WRL; L. ayresii ) are sympatric lampreys that likely represent different life history variations of a single species. Novel genetic tools are critical for differentiating WBL and WRL, whose larvae preclude morphological identification (ID) and will enable comprehensive assessment of imperiled native lampreys of the Northeastern Pacific (including WBL, WRL, and Pacific Lamprey, Entosphenus tridentatus ). We developed 47 candidate single nucleotide polymorphism (SNP) markers using whole genome resequencing of WBL ( N = 24) and WRL ( N = 15) from Ksi Ts'oohl Ts'ap Creek (Nass River, British Columbia, Canada) which are likely ecotypes distinguished by few divergent SNPs across multiple chromosomes. We used five novel candidate SNPs to perform genetic ID of WBL and WRL ecotypes in collections of mixed native lampreys from lower Columbia River tributaries ( N = 1474), Ksi Ts'oohl Ts'ap Creek ( N = 352), and ocean phase WRL from the Georgia Basin (Salish Sea, British Columbia, Canada; N = 91). Two previously published SNPs were used to ID genera, Entosphenus versus Lampetra . Morphological ID utilized photographs collected from a subset of genotyped lampreys, and high concordance was demonstrated between ID methods for genera (99%) and Lampetra ecotypes (> 98%). We characterized spatial and temporal composition of lamprey genera and ecotypes surveyed across NE Pacific tributaries under the expectation these compositions would be similar across nearby sites and across years at the same site. Proportions of lamprey genera were highly variable within regions and across years; however, Lampetra ecotypic proportions were spatially and temporally stable. WRL were rare in lower Columbia tributaries (~1% average rate among Lampetra ) and common further north (> 40% of Lampetra ). Genetic ID methods are powerful monitoring tools that create the novel ability to ascertain genera and ecotypes regardless of life stage, while increasing the efficiency of surveys by eliminating time‐intensive morphological data collection.

Silver, G. S. [Columbia River Inter‐Tribal Fish Co↗

Enhanced Resistance Pines for Improved Renewable Biofuel and Chemical Production (Technical Report)

We completed phenotyping constitutive and inducible oleoresin flow across two seasons, constitutive resin canal number and density and wood terpene content in our ADEPT2 and CCLONES populations. We completed genetic association between 19 oleoresin phenotypes and a total of 523,192 SNP markers from ADEPT2 and 13,883 SNP markers in CCLONES using four mixed linear models. A total of 293 significant SNPs (FDR = 0.20) were identified. We used the MENTOR tool to mine mechanistic connections from a multiplex network constructed from poplar multi-omic data to construct a conceptual model for a subset of these significant SNPs. Our model contains 6 transcriptional regulators in addition to 3 monoterpene synthases. To generate more lines of evidence for these significant SNPs, we completed a time course RNAseq experiment after inducing vascular zone cells to differentiate into new resin canals with a methyl jasmonate treatment, a single nuclei RNAseq that identified differentiating resin canal epithelial cells and are completing analysis for a QTL study in a hybrid pine population. The time course identified 4634 significantly down and 1890 significantly up regulated transcripts after treatment with methyl jasmonate, an inducer of new resin canal formation in the vascular cambial meristem. To analyze this large set of differentially regulated genes, we created a predictive expression network and analyzed it with random walk restart using 6 seed genes coding for transcription factors regulating xylem differentiation in poplar. Of the top ranked 200 transcripts, 119 transcripts were significant differentially expressed supporting these transcripts as potential candidates regulating resin canal formation. Analysis of single nuclei sequencing of shoot tips that contain differentiating resin canals, identified 10 clusters. One cluster was highly enriched in transcripts coding for 9 of the enzymes in the MEP pathway 3 prenyl synthetases, and 3 monoterpene synthases strongly suggesting that this cluster represents resin canal epithelial cells. We are mining the additional transcripts to create a trajectory analysis. In summary, we have identified > 10 novel genes that are strongly supported candidates for further analysis in breeding lines and for genetic engineering over- and under- expressing lines to increase wood terpene content to improve resistance to insect and fungal pathogens while simultaneously increasing terpene supplies for renewable chemicals and biofuels.

59 BASIC BIOLOGICAL SCIENCES↗

Comparing Deep Learning Approaches for Understanding Genotype × Phenotype Interactions in Biomass Sorghum

We explore the use of deep convolutional neural networks (CNNs) trained on overhead imagery of biomass sorghum to ascertain the relationship between single nucleotide polymorphisms (SNPs), or groups of related SNPs, and the phenotypes they control. We consider both CNNs trained explicitly on the classification task of predicting whether an image shows a plant with a reference or alternate version of various SNPs as well as CNNs trained to create data-driven features based on learning features so that images from the same plot are more similar than images from different plots, and then using the features this network learns for genetic marker classification. We characterize how efficient both approaches are at predicting the presence or absence of a genetic markers, and visualize what parts of the images are most important for those predictions. We find that the data-driven approaches give somewhat higher prediction performance, but have visualizations that are harder to interpret; and we give suggestions of potential future machine learning research and discuss the possibilities of using this approach to uncover unknown genotype × phenotype relationships.

Zhang, Zeyu↗

Identification of Novel Genomic Associations and Gene Candidates for Grain Starch Content in Sorghum

Starch accumulated in the endosperm of cereal grains as reserve energy for germination serves as a staple in human and animal nutrition. Unraveling genetic control for starch metabolism is important for breeding grains with high starch content. In this study, we used a sorghum association panel with 389 individuals and 141,557 single nucleotide polymorphisms (SNPs) to fit linear mixed models (LMM) for identifying genomic regions and potential candidate genes associated with starch content. Three associated genomic regions, one in chromosome (chr) 1 and two novel associations in chr-8, were identified using combination of LMM and Bayesian sparse LMM. All significant SNPs were located within protein coding genes, with SNPs ∼ 52 Mb of chr-8 encoding a Casperian strip membrane protein (CASP)-like protein (Sobic.008G111500) and a heat shock protein (HSP) 90 (Sobic.008G111600) that were highly expressed in reproductive tissues including within the embryo and endosperm. The HSP90 is a potential hub gene with gene network of 75 high-confidence first interactors that is enriched for five biochemical pathways including protein processing. The first interactors of HSP90 also showed high transcript abundance in reproductive tissues. The candidates of this study are likely involved in intricate metabolic pathways and represent candidate gene targets for source-sink activities and drought and heat stress tolerance during grain filling.

genomics↗

Whole-Genome SNP Analysis Identifies Putative Mycobacterium bovis Transmission Clusters in Livestock and Wildlife in Catalonia, Spain

The high-resolution WGS analyses of MTBC strains have provided useful insight for determining sources of infection for animal tuberculosis. In Spain, tuberculosis in livestock is caused by Mycobacterium bovis and Mycobacterium caprae, where wildlife reservoirs play an important role. We analyzed a set of 125 M. bovis isolates obtained from livestock and wildlife from Catalonia to investigate strain diversity and identify possible sources and/or causes of infection. Whole-genome SNP profiles were used for phylogenetic reconstruction and pairwise SNP distance analysis. Additionally, SNPs were investigated to identify virulence and antimicrobial resistance factors to investigate clade-specific associations. Putative transmission clusters (≤12 SNPs) were identified, and associated epidemiological metadata were used to determine possible explanatory factors for transmission. M. bovis distribution was heterogeneous, with 7 major clades and 21 putative transmission clusters. In order of importance, the explanatory factors associated were proximity and neighborhood, residual infection, livestock-wildlife interaction, shared pasture, and movement. Genes related to lipid transport and metabolism showed the highest number of SNPs. All isolates were pyrazinamide resistant, and five were additionally resistant to isoniazid, but no clade-specific associations could be determined. Our findings highlight the importance of high-resolution molecular surveillance to monitor bovine tuberculosis dynamics in a low-prevalence setting.

59 BASIC BIOLOGICAL SCIENCES↗

Development and Assessment of SNP Genotyping Arrays for Citrus and Its Close Relatives

Rapid advancements in technologies provide various tools to analyze fruit crop genomes to better understand genetic diversity and relationships and aid in breeding. Genome-wide single nucleotide polymorphism (SNP) genotyping arrays offer highly multiplexed assays at a relatively low cost per data point. We report the development and validation of 1.4M SNP Axiom® Citrus HD Genotyping Array (Citrus 15AX 1 and Citrus 15AX 2) and 58K SNP Axiom® Citrus Genotyping Arrays for Citrus and close relatives. SNPs represented were chosen from a citrus variant discovery panel consisting of 41 diverse whole-genome re-sequenced accessions of Citrus and close relatives, including eight progenitor citrus species. SNPs chosen mainly target putative genic regions of the genome and are accurately called in both Citrus and its closely related genera while providing good coverage of the nuclear and chloroplast genomes. Reproducibility of the arrays was nearly 100%, with a large majority of the SNPs classified as the most stringent class of markers, “PolyHighResolution” (PHR) polymorphisms. Concordance between SNP calls in sequence data and array data average 98%. Phylogenies generated with array data were similar to those with comparable sequence data and little affected by 3 to 5% genotyping error. Both arrays are publicly available.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

59 BASIC BIOLOGICAL SCIENCES↗

Using low-coverage whole genome sequencing (genome skimming) to delineate three introgressed species of buffalofish ( Ictiobus )

Consumption of buffalofish has been sporadically associated with Haff disease-like illnesses involving sudden onset muscle pain and weakness due to skeletal muscle rhabdomyolysis, but determination of precisely which species are associated with these illnesses has been impeded by a lack of species-specific DNA-based markers. Here, three closely related species of buffalofish native to the Mississippi River Basin (Ictiobus bubalus, Ictiobus cyprinellus and Ictiobus niger) that have previously proven genetically indistinguishable using both mitochondrial and nuclear single-locus sequencing were reliably discriminated using low-coverage whole genome sequencing (‘genome skimming’). Using 44 specimens representing the three species collected from the mid/upper (Missouri) and lower (Louisiana) regions of the species’ native ranges, the SISRS (Site Identification from Short Read Sequences) bioinformatics pipeline was adapted to (1) identify over 620Mbp of putatively homologous nuclear sequence data and (2) isolate over 140,000 single-nucleotide polymorphisms (SNPs) that supported accurate species delimitation, all without the use of a reference genome or annotation data. These sites were used to classify Ictiobus spp. samples with genome-skim data, along with a larger set (n = 67) where ultraconserved elements (UCEs) were sequenced. Analyses of whole mitochondrial data revealed more limited signal. Nearly all samples matched their purported species based on morphologic identification, but two Missouri samples morphologically identified as I. niger grouped with samples of I. bubalus, albeit with significant enrichment of I. niger SNPs. To our knowledge this is the first report of a DNA-based tool to reliably discriminate these three morphologically distinct species.

59 BASIC BIOLOGICAL SCIENCES↗

A view of the pan‐genome of domesticated Cowpea ( Vigna unguiculata [L.] Walp.)

Abstract Cowpea, Vigna unguiculata L . Walp., is a diploid warm‐season legume of critical importance as both food and fodder in sub‐Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K‐499‐35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

59 BASIC BIOLOGICAL SCIENCES↗

The genetic basis for panicle trait variation in switchgrass ( Panicum virgatum )

Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

59 BASIC BIOLOGICAL SCIENCES↗

A variant selection framework for genome graphs

Abstract Motivation Variation graph representations are projected to either replace or supplement conventional single genome references due to their ability to capture population genetic diversity and reduce reference bias. Vast catalogues of genetic variants for many species now exist, and it is natural to ask which among these are crucial to circumvent reference bias during read mapping. Results In this work, we propose a novel mathematical framework for variant selection, by casting it in terms of minimizing variation graph size subject to preserving paths of length α with at most δ differences. This framework leads to a rich set of problems based on the types of variants [e.g. single nucleotide polymorphisms (SNPs), indels or structural variants (SVs)], and whether the goal is to minimize the number of positions at which variants are listed or to minimize the total number of variants listed. We classify the computational complexity of these problems and provide efficient algorithms along with their software implementation when feasible. We empirically evaluate the magnitude of graph reduction achieved in human chromosome variation graphs using multiple α and δ parameter values corresponding to short and long-read resequencing characteristics. When our algorithm is run with parameter settings amenable to long-read mapping (α = 10 kbp, δ = 1000), 99.99% SNPs and 73% SVs can be safely excluded from human chromosome 1 variation graph. The graph size reduction can benefit downstream pan-genome analysis. Availability and implementation https://github.com/AT-CG/VF. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

Population Genomics of Pseudocercospora griseola Reveals New Groups in the Middle American Clade and the Presence of the Endophytic Bacterium Achromobacter xylosoxidans

Angular leaf spot (ALS), caused by Pseudocercospora griseola is an important disease of common beans. P. griseola, is highly variable and has co-evolved with its host. In this study, 48 isolates of P. griseola from Puerto Rico, Guatemala, Honduras and Tanzania were sequenced (3RADseq), resulting in the de novo assembly of 42,214 contigs. Phylogenomic, population genetic structure and principal component analyses using 1,260 SNPs divided these isolates into two populations, Andean and Middle American, while the Middle American population was further divided into three sub-populations. There were moderate to high levels of differentiation between P. griseola populations, with pairwise Fst values ranging from 0.11 to 0.95. The Andean population was composed of isolates from Tanzania, and was separated from the Middle American population (Fst = 0.95). The Middle American population was separated into 3 subpopulations including isolates from: 1. Guatemala and Honduras, 2. Tanzania, and 3. Puerto Rico. Pathogenicity testing of 27 isolates from Puerto Rico, using 12 common bean differential lines, identified ten races, but these races were not associated with SNPs found in virulence genes. DNA of an endophytic bacterium (Achromobacter xylosoxidans) was found in seven mildly virulent isolates suggesting a possible role of the bacterium in the observed virulence patterns. To understand the evolution and diversity of P. griseola, further study of the virulence genes and the interactions among the endophytic bacterium, the fungus, and the host plant is required. Such information is critical to inform breeding strategies for the development of resistant germplasm and cultivars.

Serrato-Diaz, Luz M. [U.S. Department of Agricultu↗