A population of deletion mutants and an integrated mapping and exome-seq pipeline for gene discovery in maize

To better understand maize endosperm filling and maturation, we used gamma irradiation of the B73 maize reference line to generate mutants with opaque endosperm and reduced kernel fill phenotypes, and created a population of 1788 lines including 47 Mo17 × F2s showing stable, segregating and viable kernel phenotypes. For molecular characterization of the mutants, we developed a novel functional genomics platform that combined Bulked Segregant RNA and Exome sequencing (BSREx-seq) to map causative mutations and identify candidate genes within mapping intervals. To exemplify the utility of the mutants and provide proof-ofconcept for the bioinformatics platform, we present detailed characterization of line 937; an opaque mutant harboring a 6203-bp in-frame deletion covering six exons within the Opaque-1 gene. In addition, we describe mutant line 146 which contains a 4.8 Kb intra-gene deletion within the Sugary-1 gene and line 916 in which an 8.6 Kb deletion knocks out a Cyclin A2 gene. The publically available algorithm developed in this work improves the identification of causative deletions and its corresponding gaps within mapping peaks. This study demonstrates the utility of gamma irradiation for forward genetics in large non-dense genomes such as maize since deletions often affect single genes. Furthermore, we show how this classical mutagenesis method becomes applicable for functional genomics when combined with state-of-the-art genomics tools

A population of deletion mutants and an integrated mapping and exome-seq pipeline for gene discovery in maize

To better understand maize endosperm filling and maturation, we used gamma irradiation of the B73 maize reference line to generate mutants with opaque endosperm and reduced kernel fill phenotypes, and created a population of 1788 lines including 47 Mo17 × F2s showing stable, segregating and viable kernel phenotypes. For molecular characterization of the mutants, we developed a novel functional genomics platform that combined Bulked Segregant RNA and Exome sequencing (BSREx-seq) to map causative mutations and identify candidate genes within mapping intervals. To exemplify the utility of the mutants and provide proof-ofconcept for the bioinformatics platform, we present detailed characterization of line 937; an opaque mutant harboring a 6203-bp in-frame deletion covering six exons within the Opaque-1 gene. In addition, we describe mutant line 146 which contains a 4.8 Kb intra-gene deletion within the Sugary-1 gene and line 916 in which an 8.6 Kb deletion knocks out a Cyclin A2 gene. The publically available algorithm developed in this work improves the identification of causative deletions and its corresponding gaps within mapping peaks. This study demonstrates the utility of gamma irradiation for forward genetics in large non-dense genomes such as maize since deletions often affect single genes. Furthermore, we show how this classical mutagenesis method becomes applicable for functional genomics when combined with state-of-the-art genomics tools

Identification of seedling resistance in wild oat relatives against oat crown rust

Faculty advisor: Melania FigueroaThis research was supported by the Undergraduate Research Opportunities Program (UROP)

Endosperm Tolerance of Paternal Aneuploidy Allows Radiation Hybrid Mapping of the Wheat D-Genome and a Measure of γ Ray-Induced Chromosome Breaks

Physical mapping and genome sequencing are underway for the ≈17 Gb wheat genome. Physical mapping methods independent of meiotic recombination, such as radiation hybrid (RH) mapping, will aid precise anchoring of BAC contigs in the large regions of suppressed recombination in Triticeae genomes. Reports of endosperm development following pollination with irradiated pollen at dosages that cause embryo abortion prompted us to investigate endosperm as a potential source of RH mapping germplasm. Here, we report a novel approach to construct RH based physical maps of all seven D-genome chromosomes of the hexaploid wheat 'Chinese Spring', simultaneously. An 81-member subset of endosperm samples derived from 20-Gy irradiated pollen was genotyped for deletions, and 737 markers were mapped on seven D-genome chromosomes. Analysis of well-defined regions of six chromosomes suggested a map resolution of ~830 kb could be achieved; this estimate was validated with assays of markers from a sequenced contig. We estimate that the panel contains ~6,000 deletion bins for D-genome chromosomes and will require ~18,000 markers for high resolution mapping. Map-based deletion estimates revealed a majority of 1-20 Mb interstitial deletions suggesting mutagenic repair of double-strand breaks in pollen provides a useful resource for RH mapping and map based cloning studie

Physical mapping resources for large plant genomes: radiation hybrids for wheat D-genome progenitor Aegilops tauschii

BACKGROUND: Development of a high quality reference sequence is a daunting task in crops like wheat with large (~17Gb), highly repetitive (>80%) and polyploid genome. To achieve complete sequence assembly of such genomes, development of a high quality physical map is a necessary first step. However, due to the lack of recombination in certain regions of the chromosomes, genetic mapping, which uses recombination frequency to map marker loci, alone is not sufficient to develop high quality marker scaffolds for a sequence ready physical map. Radiation hybrid (RH) mapping, which uses radiation induced chromosomal breaks, has proven to be a successful approach for developing marker scaffolds for sequence assembly in animal systems. Here, the development and characterization of a RH panel for the mapping of D-genome of wheat progenitor Aegilops tauschii is reported. RESULTS: Radiation dosages of 350 and 450 Gy were optimized for seed irradiation of a synthetic hexaploid (AABBDD) wheat with the D-genome of Ae. tauschii accession AL8/78. The surviving plants after irradiation were crossed to durum wheat (AABB), to produce pentaploid RH(1)s (AABBD), which allows the simultaneous mapping of the whole D-genome. A panel of 1,510 RH(1) plants was obtained, of which 592 plants were generated from the mature RH(1) seeds, and 918 plants were rescued through embryo culture due to poor germination (<3%) of mature RH(1) seeds. This panel showed a homogenous marker loss (2.1%) after screening with SSR markers uniformly covering all the D-genome chromosomes. Different marker systems mostly detected different lines with deletions. Using markers covering known distances, the mapping resolution of this RH panel was estimated to be <140kb. Analysis of only 16 RH lines carrying deletions on chromosome 2D resulted in a physical map with cM/cR ratio of 1:5.2 and 15 distinct bins. Additionally, with this small set of lines, almost all the tested ESTs could be mapped. A set of 399 most informative RH lines with an average deletion frequency of ~10% were identified for developing high density marker scaffolds of the D-genome. CONCLUSIONS: The RH panel reported here is the first developed for any wild ancestor of a major cultivated plant species. The results provided insight into various aspects of RH mapping in plants, including the genetically effective cell number for wheat (for the first time) and the potential implementation of this technique in other plant species. This RH panel will be an invaluable resource for mapping gene based markers, developing a complete marker scaffold for the whole genome sequence assembly, fine mapping of markers and functional characterization of genes and gene networks present on the D-genome

Genome-Wide Association Mapping for Yield and Related Traits Under Drought Stressed and Non-stressed Environments in Wheat

Understanding the genetics of drought tolerance in hard red spring wheat (HRSW) in northern USA is a prerequisite for developing drought-tolerant cultivars for this region. An association mapping (AM) study for drought tolerance in spring wheat in northern USA was undertaken using 361 wheat genotypes and Infinium 90K single-nucleotide polymorphism (SNP) assay. The genotypes were evaluated in nine different locations of North Dakota (ND) for plant height (PH), days to heading (DH), yield (YLD), test weight (TW), and thousand kernel weight (TKW) under rain-fed conditions. Rainfall data and soil type of the locations were used to assess drought conditions. A mixed linear model (MLM), which accounts for population structure and kinship (PC+K), was used for marker–trait association. A total of 69 consistent QTL involved with drought tolerance-related traits were identified, with p ≤ 0.001. Chromosomes 1A, 3A, 3B, 4B, 4D, 5B, 6A, and 6B were identified to harbor major QTL for drought tolerance. Six potential novel QTL were identified on chromosomes 3D, 4A, 5B, 7A, and 7B. The novel QTL were identified for DH, PH, and TKW. The findings of this study can be used in marker-assisted selection (MAS) for drought-tolerance breeding in spring wheat

De novo assembly and phasing of dikaryotic genomes from two isolates of puccinia coronata f. Sp. avenae, the causal agent of oat crown rust

Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenae.This work was funded by the USDA-ARS and the University of Minnesota Standard Cooperative Agreement (grant 3002-11031-00053115 shared by S.F.K. and M.F.), the University of Minnesota Experimental Station USDA-NIFA Hatch/Figueroa project MIN22-058, and an Organization for Economic Cooperation and Development Fellowship to M.F. M.E.M. was partially supported by a USDA-NIFA Postdoctoral Fellowship Award (2017-67012-26117). J.S. was supported by an OCE Postdoctoral Fellowship. R.F.P. receives funding from the Australian Grains Research Development Corporation (grant US00067). J.M.P. was supported by the Northern Research Station of the USDA Forest Service

Development of a D genome specific marker resource for diploid and hexaploid wheat

Background: Mapping and map-based cloning of genes that control agriculturally and economically important traits remain great challenges for plants with complex highly repetitive genomes such as those within the grass tribe, Triticeae. Mapping limitations in the Triticeae are primarily due to low frequencies of polymorphic gene markers and poor genetic recombination in certain genetic regions. Although the abundance of repetitive sequence may pose common problems in genome analysis and sequence assembly of large and complex genomes, they provide repeat junction markers with random and unbiased distribution throughout chromosomes. Hence, development of a high-throughput mapping technology that combine both gene-based and repeat junction-based markers is needed to generate maps that have better coverage of the entire genome. Results: In this study, the available genomics resource of the diploid Aegilop tauschii, the D genome donor of bread wheat, were used to develop genome specific markers that can be applied for mapping in modern hexaploid wheat. A NimbleGen array containing both gene-based and repeat junction probe sequences derived from Ae. tauschii was developed and used to map the Chinese Spring nullisomic-tetrasomic lines and deletion bin lines of the D genome chromosomes. Based on these mapping data, we have now anchored 5,171 repeat junction probes and 10,892 gene probes, corresponding to 5,070 gene markers, to the delineated deletion bins of the D genome. The order of the gene-based markers within the deletion bins of the Chinese Spring can be inferred based on their positions on the Ae. tauschii genetic map. Analysis of the probe sequences against the Chinese Spring chromosome sequence assembly database facilitated mapping of the NimbleGen probes to the sequence contigs and allowed assignment or ordering of these sequence contigs within the deletion bins. The accumulated length of anchored sequence contigs is about 155 Mb, representing ~ 3.2 % of the D genome. A specific database was developed to allow user to search or BLAST against the probe sequence information and to directly download PCR primers for mapping specific genetic loci. Conclusions: In bread wheat, aneuploid stocks have been extensively used to assign markers linked with genes/traits to chromosomes, chromosome arms, and their specific bins. Through this study, we added thousands of markers to the existing wheat chromosome bin map, representing a significant step forward in providing a resource to navigate the wheat genome. The database website (http://probes.pw.usda.gov/ATRJM/) provides easy access and efficient utilization of the data. The resources developed herein can aid map-based cloning of traits of interest and the sequencing of the D genome of hexaploid wheat.Keywords: Recombination, Repeat junction markers, Genetic map, NimbleGen array, Molecular markers, Wheat deletion bin

DNA repair and crossing over favor similar chromosome regions as discovered in radiation hybrid of Triticum

Background: The uneven distribution of recombination across the length of chromosomes results in inaccurate estimates of genetic to physical distances. In wheat (Triticum aestivum L.) chromosome 3B, it has been estimated that 90% of the cross over events occur in distal sub-telomeric regions representing 40% of the chromosome. Radiation hybrid (RH) mapping which does not rely on recombination is a strategy to map genomes and has been widely employed in animal species and more recently in some plants. RH maps have been proposed to provide i) higher and ii) more uniform resolution than genetic maps, and iii) to be independent of the distribution patterns observed for meiotic recombination. An in vivo RH panel was generated for mapping chromosome 3B of wheat in an attempt to provide a complete scaffold for this ~1 Gb segment of the genome and compare the resolution to previous genetic maps. Results: A high density RH map with 541 marker loci anchored to chromosome 3B spanning a total distance of 1871.9 cR was generated. Detailed comparisons with a genetic map of similar quality confirmed that i) the overall resolution of the RH map was 10.5 fold higher and ii) six fold more uniform. A significant interaction (r = 0.879 at p = 0.01) was observed between the DNA repair mechanism and the distribution of crossing-over events. This observation could be explained by accepting the possibility that the DNA repair mechanism in somatic cells is affected by the chromatin state in a way similar to the effect that chromatin state has on recombination frequencies in gametic cells. Conclusions: The RH data presented here support for the first time in vivo the hypothesis of non-casual interaction between recombination hot-spots and DNA repair. Further, two major hypotheses are presented on how chromatin compactness could affect the DNA repair mechanism. Since the initial RH application 37 years ago, we were able to show for the first time that the iii) third hypothesis of RH mapping might not be entirely correct.Keywords: Physical mapping, Deletion mutant, Non homologous end joining, Wheat chromosome 3B, Chromatin\, , Radiation hybrid, Gamma radiatio