Eric W. Jackson

Model SNP development for complex genomes based on hexaploid oat using high-throughput 454 sequencing technology

BackgroundGenetic markers are pivotal to modern genomics research; however, discovery and genotyping of molecular markers in oat has been hindered by the size and complexity of the genome, and by a scarcity of sequence data. The purpose of this study was to generate oat expressed sequence tag (EST) information, develop a bioinformatics pipeline for SNP discovery, and establish a method for rapid, cost-effective, and straightforward genotyping of SNP markers in complex polyploid genomes such as oat.ResultsBased on cDNA libraries of four cultivated oat genotypes, approximately 127,000 contigs were assembled from approximately one million Roche 454 sequence reads. Contigs were filtered through a novel bioinformatics pipeline to eliminate ambiguous polymorphism caused by subgenome homology, and 96 in silico SNPs were selected from 9,448 candidate loci for validation using high-resolution melting (HRM) analysis. Of these, 52 (54%) were polymorphic between parents of the Ogle1040 × TAM O-301 (OT) mapping population, with 48 segregating as single Mendelian loci, and 44 being placed on the existing OT linkage map. Ogle and TAM amplicons from 12 primers were sequenced for SNP validation, revealing complex polymorphism in seven amplicons but general sequence conservation within SNP loci. Whole-amplicon interrogation with HRM revealed insertions, deletions, and heterozygotes in secondary oat germplasm pools, generating multiple alleles at some primer targets. To validate marker utility, 36 SNP assays were used to evaluate the genetic diversity of 34 diverse oat genotypes. Dendrogram clusters corresponded generally to known genome composition and genetic ancestry.ConclusionsThe high-throughput SNP discovery pipeline presented here is a rapid and effective method for identification of polymorphic SNP alleles in the oat genome. The current-generation HRM system is a simple and highly-informative platform for SNP genotyping. These techniques provide a model for SNP discovery and genotyping in other species with complex and poorly-characterized genomes.

Using Genotyping-By-Sequencing (GBS) for Genomic Discovery in Cultivated Oat

Advances in next-generation sequencing offer high-throughput and cost-effective genotyping alternatives, including genotyping-by-sequencing (GBS). Results have shown that this methodology is efficient for genotyping a variety of species, including those with complex genomes. To assess the utility of GBS in cultivated hexaploid oat (Avena sativa L.), seven bi-parental mapping populations and diverse inbred lines from breeding programs around the world were studied. We examined technical factors that influence GBS SNP calls, established a workflow that combines two bioinformatics pipelines for GBS SNP calling, and provided a nomenclature for oat GBS loci. The high-throughput GBS system enabled us to place 45,117 loci on an oat consensus map, thus establishing a positional reference for further genomic studies. Using the diversity lines, we estimated that a minimum density of one marker per 2 to 2.8 cM would be required for genome-wide association studies (GWAS), and GBS markers met this density requirement in most chromosome regions. We also demonstrated the utility of GBS in additional diagnostic applications related to oat breeding. We conclude that GBS is a powerful and useful approach, which will have many additional applications in oat breeding and genomic studies.

SNP Discovery and Chromosome Anchoring Provide the First Physically-Anchored Hexaploid Oat Map and Reveal Synteny with Model Species

A physically anchored consensus map is foundational to modern genomics research; however, construction of such a map in oat (Avena sativa L., 2n = 6x = 42) has been hindered by the size and complexity of the genome, the scarcity of robust molecular markers, and the lack of aneuploid stocks. Resources developed in this study include a modified SNP discovery method for complex genomes, a diverse set of oat SNP markers, and a novel chromosome-deficient SNP anchoring strategy. These resources were applied to build the first complete, physically-anchored consensus map of hexaploid oat. Approximately 11,000 high-confidence in silico SNPs were discovered based on nine million inter-varietal sequence reads of genomic and cDNA origin. GoldenGate genotyping of 3,072 SNP assays yielded 1,311 robust markers, of which 985 were mapped in 390 recombinant-inbred lines from six bi-parental mapping populations ranging in size from 49 to 97 progeny. The consensus map included 985 SNPs and 68 previously-published markers, resolving 21 linkage groups with a total map distance of 1,838.8 cM. Consensus linkage groups were assigned to 21 chromosomes using SNP deletion analysis of chromosome-deficient monosomic hybrid stocks. Alignments with sequenced genomes of rice and Brachypodium provide evidence for extensive conservation of genomic regions, and renewed encouragement for orthology-based genomic discovery in this important hexaploid species. These results also provide a framework for high-resolution genetic analysis in oat, and a model for marker development and map construction in other species with complex genomes and limited resources.

Identification of novel genomic regions associated with resistance to Pyrenophora tritici-repentis races 1 and 5 in spring wheat landraces using association analysis

Tan spot, caused by Pyrenophora tritici-repentis, is a major foliar disease of wheat worldwide. Host plant resistance is the best strategy to manage this disease. Traditionally, bi-parental mapping populations have been used to identify and map quantitative trait loci (QTL) affecting tan spot resistance in wheat. The association mapping (AM) could be an alternative approach to identify QTL based on linkage disequilibrium (LD) within a diverse germplasm set. In this study, we assessed resistance to P. tritici-repentis races 1 and 5 in 567 spring wheat landraces from the USDA-ARS National Small Grains Collection (NSGC). Using 832 diversity array technology (DArT) markers, QTL for resistance to P. tritici-repentis races 1 and 5 were identified. A linear model with principal components suggests that at least seven and three DArT markers were significantly associated with resistance to P. tritici-repentis races 1 and 5, respectively. The DArT markers associated with resistance to race 1 were detected on chromosomes 1D, 2A, 2B, 2D, 4A, 5B, and 7D and explained 1.3–3.1% of the phenotypic variance, while markers associated with resistance to race 5 were distributed on 2D, 6A and 7D, and explained 2.2–5.9% of the phenotypic variance. Some of the genomic regions identified in this study correspond to previously identified loci responsible for resistance to P. tritici-repentis, offering validation for our AM approach. Other regions identified were novel and could possess genes useful for resistance breeding. Some DArT markers associated with resistance to race 1 also were localized in the same regions of wheat chromosomes where QTL for resistance to yellow rust, leaf rust and powdery mildew, have been mapped previously. This study demonstrates that AM can be a useful approach to identify and map novel genomic regions involved in resistance to P. tritici-repentis.

Association Mapping of Quantitative Trait Loci in Spring Wheat Landraces Conferring Resistance to Bacterial Leaf Streak and Spot Blotch

Bacterial leaf streak (BLS), caused by Xanthomonas translucens pv. undulosa (Smith et al.) Bragard et al., and spot blotch (SB), caused by Cochliobolus sativus (S. Ito & Kurib.) Drechs. ex Dastur, are two emerging diseases of wheat (Triticum aestivum L.). To achieve sustainable disease management strategies and reduce yield losses, identifying new genes that confer quantitative resistance would benefit resistance breeding efforts. The main objective of this study was to use association mapping (AM) with 832 polymorphic Diversity Arrays Technology (DArT) markers to identify genomic regions associated with resistance to BLS and SB in 566 spring wheat landraces. From data analysis of this diverse panel of wheat accessions, we discovered five novel genomic regions significantly associated with resistance to BLS on chromosomes 1A, 4A, 4B, 6B, and 7D. Similarly, four genomic regions were found to be associated with resistance to SB on chromosomes 1A, 3B, 7B, and 7D. A high degree of linkage disequilibrium (LD) decayed over short genetic distance in the set of wheat accessions studied, and some of these genomic regions appear to be involved in multiple disease resistance (MDR). These results suggest that the AM approach provides a platform for discovery of resistance conditioned by multiple genes with quantitative effects, which could be validated and deployed in wheat breeding programs.

Characterization and Mapping of Oat Crown Rust Resistance Genes Using Three Assessment Methods

ABSTRACT Resistance is the primary means of control for crown rust of oat (Avena sativa L.), caused by Puccinia coronata f. sp. avenae, and better knowledge of the genetics of resistance will enhance resistance breeding. Disease data were generated in the field and greenhouse for parents and recombinant inbred lines of the Ogle/TAM O-301 (OT) oat mapping population using (i) a new quantitative assay that employs quantitative real-time polymerase chain reaction (q-PCR) to estimate fungal growth in the host, (ii) digital image analysis, and (iii) visual ratings. The objectives of this study were to evaluate each assessment methods ability to map a major gene from cv. Ogle and potential quantitative trait loci (QTL) contributed by Ogle and TAM O-301. All three assessment methods identified the major gene in Ogle, which was mapped to linkage group OT6. The resolution produced by q-PCR, however, enabled more precise mapping of the major gene. Quantitative analysis indicated that 64% of the phenotypic variation was accounted for using q-PCR, whereas 41 and 52% were accounted for using visual and digital assessments, respectively. Data generated by q-PCR permitted identification of QTL on linkage groups OT32, accounting for 6% of the phenotypic variation, and OT2, accounting for 4% of the variation. QTL on both OT32 and OT2 were conferred by TAM O-301, one of which (OT2) was indiscernible using data from the visual and digital assessments. The new method of precisely phenotyping crown rust resistance provided a more accurate and thorough means of dissecting resistance in the OT mapping population. Similar methods could be developed and applied to other important cereal rust diseases.

A SNP Genotyping Array for Hexaploid Oat

Recognizing a need in cultivated hexaploid oat (Avena sativa L.) for a reliable set of reference single nucleotide polymorphisms (SNPs), we have developed a 6000 (6K) BeadChip design containing 257 Infinium I and 5486 Infinium II designs corresponding to 5743 SNPs. Of those, 4975 SNPs yielded successful assays after array manufacturing. These SNPs were discovered based on a variety of bioinformatics pipelines in complementary DNA (cDNA) and genomic DNA originating from 20 or more diverse oat cultivars. The array was validated in 1100 samples from six recombinant inbred line (RIL) mapping populations and sets of diverse oat cultivars and breeding lines, and provided approximately 3500 discernible Mendelian polymorphisms. Here, we present an annotation of these SNPs, including methods of discovery, gene identification and orthology, population‐genetic characteristics, and tentative positions on an oat consensus map. We also evaluate a new cluster‐based method of calling SNPs. The SNP design sequences are made publicly available, and the full SNP genotyping platform is available for commercial purchase from an independent third party.

A consensus map in cultivated hexaploid oat reveals conserved grass synteny with substantial subgenome rearrangement

We constructed a hexaploid oat consensus map from 12 populations representing 19 parents. The map represents the most common physical chromosome arrangements in oat. Deviations from the consensus map may indicate physical rearrangements. Large chromosomal translocations vary among different varieties. There is regional synteny with rice but considerable subgenome rearrangement.

Association Mapping of Quantitative Resistance to Phaeosphaeria nodorum in Spring Wheat Landraces from the USDA National Small Grains Collection

Stagonospora nodorum blotch (SNB), caused by Phaeosphaeria nodorum, is a destructive disease of wheat (Triticum aestivum) found throughout the United States. Host resistance is the only economically feasible option for managing the disease; however, few SNB-resistant wheat cultivars are known to exist. In this study, we report findings from an association mapping (AM) of resistance to P. nodorum in 567 spring wheat landraces of diverse geographic origin. The accessions were evaluated for seedling resistance to P. nodorum in a greenhouse. Phenotypic data and 625 polymorphic diversity array technology (DArT) markers have been used for linkage disequilibrium (LD) and association analyses. The results showed that seven DArT markers on five chromosomes (2D, 3B, 5B, 6A, and 7A) were significantly associated with resistance to P. nodorum. Genetic regions on 2D, 3B, and 5B correspond to previously mapped quantitative trait loci (QTL) conferring resistance to P. nodorum whereas the remaining QTL appeared to be novel. These results demonstrate that the use of AM is an effective method for identifying new genomic regions associated with resistance to P. nodorum in spring wheat landraces. Additionally, the novel resistance found in this study could be useful in wheat breeding aimed at controlling SNB.

Qualitative and quantitative trait loci conditioning resistance to Puccinia coronata pathotypes NQMG and LGCG in the oat (Avena sativa L.) cultivars Ogle and TAM O-301

Mapping disease resistance loci relies on the type and precision of phenotypic measurements. For crown rust of oat, disease severity is commonly assessed based on visual ratings of infection types (IT) and/or diseased leaf area (DLA) of infected plants in the greenhouse or field. These data can be affected by several variables including; (i) non-uniform disease development in the field; (ii) atypical symptom development in the greenhouse; (iii) the presence of multiple pathogenic races or pathotypes in the field, and (iv) rating bias. To overcome these limitations, we mapped crown rust resistance to single isolates in the Ogle/TAM O-301 (OT) recombinant inbred line (RIL) population using detailed measurements of IT, uredinia length (UL) and relative fungal DNA (FDNA) estimates determined by q-PCR. Measurements were taken on OT parents and recombinant inbred lines (RIL) inoculated with Puccinia coronata pathotypes NQMG and LGCG in separate greenhouse and field tests. Qualitative mapping identified an allele conferred by TAM O-301 on linkage group (LG) OT-11, which produced a bleached fleck phenotype to both NQMG and LGCG. Quantitative mapping identified two major quantitative trait loci (QTL) originating from TAM O-301 on LGs OT-11 and OT-32 which reduced UL and FDNA of both isolates in all experiments. Additionally, minor QTLs that reduced UL and FDNA were detected on LGs OT-15 and OT-8, originating from TAM O-301, and on LG OT-27, originating from Ogle. Detailed assessments of the OT population using two pathotypes in both the greenhouse and field provided comprehensive information to effectively map the genes responsible for crown rust resistance in Ogle and TAM O-301 to NQMG and LGCG.