Three Decades of Rust Surveys in the United States Reveal Drastic Virulence Changes in Oat Crown Rust

To better understand how the pathogenicity of the oat crown rust pathogen Puccinia coronata f. sp. avenae (Pca) has changed in the United States, 30 years of United States Department of Agriculture (USDA) survey isolates (n = 5,456) tested on 30 to 40 differential lines were analyzed for overall and Pc-resistance-gene-specific virulence trends and correlations. Pca is incredibly pathologically diverse, with 88% of races represented by a single isolate. There are a slightly higher proportion of unique races from the Northern region of the United States, and for one fourth of the years, Northern region isolates were significantly more virulent than Southern isolates, which supports the idea that sexual recombination in this region is mediated by the alternate host as a major factor in creating new races. However, there is also support for regular isolate movement between North and South regions as isolates in the United States are steadily accumulating virulences at a rate of 0.35 virulences per year. Virulence significantly increased for 23 and decreased for four of the 40 differential lines. In the past few years, virulence has reached 90% or greater for 16 differential lines. There were also strong correlations in virulence for certain Pc genes that are likely identical, allelic, or target the same or closely linked pathogen effectors (e.g., Pc39, Pc55, and Pc71), and the results were largely in concordance with recent genome-wide association study (GWAS) effector studies using USDA isolate subsets. Understanding changes in Pca pathogenicity is essential for the responsible deployment and management of Pc resistance genes for sustainable and profitable oat production.

Genome-Enabled Analysis of Population Dynamics and Virulence-Associated Loci in the Oat Crown Rust Fungus Puccinia coronata f. sp. avenae

Puccinia coronata f. sp. avenae (Pca) is an important fungal pathogen causing crown rust that impacts oat production worldwide. Genetic resistance for crop protection against Pca is often overcome by the rapid virulence evolution of the pathogen. This study investigated the factors shaping adaptive evolution of Pca using pathogen populations from distinct geographic regions within the United States and South Africa. Phenotypic and genome-wide sequencing data of these diverse Pca collections, including 217 isolates, uncovered phylogenetic relationships and established distinct genetic composition between populations from northern and southern regions from the United States and South Africa. The population dynamics of Pca involve a bidirectional movement of inoculum between northern and southern regions of the United States and contributions from clonality and sexuality. The population from South Africa is solely clonal. A genome-wide association study (GWAS) employing a haplotype-resolved Pca reference genome was used to define 11 virulence-associated loci corresponding to 25 oat differential lines. These regions were screened to determine candidate Avr effector genes. Overall, the GWAS results allowed us to identify the underlying genetic factors controlling pathogen recognition in an oat differential set used in the United States to assign pathogen races (pathotypes). Key GWAS findings support complex genetic interactions in several oat lines, suggesting allelism among resistance genes or redundancy of genes included in the differential set, multiple resistance genes recognizing genetically linked Avr effector genes, or potentially epistatic relationships. A careful evaluation of the composition of the oat differential set accompanied by the development or implementation of molecular markers is recommended. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Mapping crown rust resistance in the oat diploid accession PI 258731 (Avena strigosa)

Oat crown rust, caused by Puccinia coronata Corda f. sp. avenae Eriks. (Pca), is a major biotic impediment to global oat production. Crown rust resistance has been described in oat diploid species A. strigosa accession PI 258731 and resistance from this accession has been successfully introgressed into hexaploid A. sativa germplasm. The current study focuses on 1) mapping the location of QTL containing resistance and evaluating the number of quantitative trait loci (QTL) conditioning resistance in PI 258731; 2) understanding the relationship between the original genomic location in A. strigosa and the location of the introgression in the A. sativa genome; 3) identifying molecular markers tightly linked with PI 258731 resistance loci that could be used for marker assisted selection and detection of this resistance in diverse A. strigosa accessions. To achieve this, A. strigosa accessions, PI 258731 and PI 573582 were crossed to produce 168 F5:6 recombinant inbred lines (RILs) through single seed descent. Parents and RILs were genotyped with the 6K Illumina SNP array which generated 168 segregating SNPs. Seedling reactions to two isolates of Pca (races TTTG, QTRG) were conditioned by two genes (0.6 cM apart) in this population. Linkage mapping placed these two resistant loci to 7.7 (QTRG) to 8 (TTTG) cM region on LG7. Field reaction data was used for QTL analysis and the results of interval mapping (MIM) revealed a major QTL (QPc.FD-AS-AA4) for field resistance. SNP marker assays were developed and tested in 125 diverse A. strigosa accessions that were rated for crown rust resistance in Baton Rouge, LA and Gainesville, FL and as seedlings against races TTTG and QTRG. Our data proposed SNP marker GMI_ES17_c6425_188 as a candidate for use in marker-assisted selection, in addition to the marker GMI_ES02_c37788_255 suggested by Rine's group, which provides an additional tool in facilitating the utilization of this gene in oat breeding programs.

Virulence patterns of oat crown rust in Australia - season 2022

Puccinia coronata f. sp. avenae (Pca) is an important foliar pathogen of oat which causes crown rust disease. The virulence profile of 48 Pca isolates derived from different locations in Australia was characterised using a collection of oat lines often utilised in rust surveys in the USA and Australia. This analysis indicates that Pca populations in Eastern Australia are broadly virulent, in contrast to the population in Western Australia (WA). Several oat lines/Pc genes are effective against all rust samples collected from WA, suggesting they may provide useful resistance in this region if deployed in combination. We identified 19 lines from the USA oat differential set that display disease resistance to Pca in WA, some in agreement with previous rust survey reports. We adopted the 10-letter nomenclature system to define oat crown rust races in Australia and compare the frequency of those virulence traits to published data from the USA. Based on this nomenclature, 42 unique races were detected among the 48 isolates, reflecting the high diversity of virulence phenotypes for Pca in Australia. Nevertheless, the Pca population in the USA is substantially more broadly virulent than that of Australia. Close examination of resistance profiles for the oat differential set lines after infection with Pca supports hypotheses of allelism or redundancy among Pc genes or the presence of several resistance genes in some oat differential lines. These findings illustrate the need to deconvolute the oat differential set using molecular tools.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

Genetic Analysis and Physical Mapping of Oat Adult Plant Resistance Loci Against Puccinia coronata f. sp. avenae

Six quantitative trait loci (QTLs) for adult plant resistance against oat crown rust (Puccinia coronata f. sp. avenae) were identified from mapping three recombinant inbred populations. Using genotyping-by-sequencing with markers called against the OT3098 v1 reference genome, the QTLs were mapped on six different chromosomes: Chr1D, Chr4D, Chr5A, Chr5D, Chr7A, and Chr7C. Composite interval mapping with marker cofactor selection showed that the phenotypic variance explained by all identified QTLs for coefficient of infection range from 12.2 to 46.9%, whereas heritability estimates ranged from 0.11 to 0.38. The significant regions were narrowed down to intervals of 3.9 to 25 cM, equivalent to physical distances of 11 to 133 Mb. At least two flanking single-nucleotide polymorphism markers were identified within 10 cM of each QTL that could be used in marker-assisted introgression, pyramiding, and selection. The additive effects of the QTLs in each population were determined using single-nucleotide polymorphism haplotype data, which showed a significantly lower coefficient of infection in lines homozygous for the resistant alleles. Analysis of pairwise linkage disequilibrium also revealed high correlation of markers and presence of linkage blocks in the significant regions. To further facilitate marker-assisted breeding, polymerase chain reaction allelic competitive extension (PACE) markers for the adult plant resistance loci were developed. Putative candidate genes were also identified in each of the significant regions, which include resistance gene analogs that encode for kinases, ligases, and predicted receptors of avirulence proteins from pathogens.

A reference-anchored oat linkage map reveals quantitative trait loci conferring adult plant resistance to crown rust (Puccinia coronata f. sp. avenae)

We mapped three adult plant resistance (APR) loci on oat chromosomes 4D and 6C and developed flanking KASP/PACE markers for marker-assisted selection and gene pyramiding. Using sequence orthology search and the available oat genomic and transcriptomic data, we surveyed these genomic regions for genes that may control disease resistance. Sources of durable disease resistance are needed to minimize yield losses in cultivated oat caused by crown rust (Puccinia coronata f. sp. avenae). In this study, we developed five oat recombinant inbred line mapping populations to identify sources of adult plant resistance from crosses between five APR donors and Otana, a susceptible variety. The preliminary bulk segregant mapping based on allele frequencies showed two regions in linkage group Mrg21 (Chr4D) that are associated with the APR phenotype in all five populations. Six markers from these regions in Chr4D were converted to high-throughput allele specific PCR assays and were used to genotype all individuals in each population. Simple interval mapping showed two peaks in Chr4D, named QPc.APR-4D.1 and QPc.APR-4D.2, which were detected in the OtanaA/CI4706-2 and OtanaA/CI9416-2 and in the Otana/PI189733, OtanaD/PI260616, and OtanaA/CI8000-4 populations, respectively. These results were validated by mapping two entire populations, Otana/PI189733 and OtanaA/CI9416, genotyped using Illumina HiSeq, in which polymorphisms were called against the OT3098 oat reference genome. Composite interval mapping results confirmed the presence of the two quantitative trait loci (QTL) located on oat chromosome 4D and an additional QTL with a smaller effect located on chromosome 6C. This mapping approach also narrowed down the physical intervals to between 5 and 19 Mb, and indicated that QPc.APR-4D.1, QPc.APR-4D.2, and QPc.APR-6C explained 43.4%, 38.5%, and 21.5% of the phenotypic variation, respectively. In a survey of the gene content of each QTL, several clusters of disease resistance genes that may contribute to APR were found. The allele specific PCR markers developed for these QTL regions would be beneficial for marker-assisted breeding, gene pyramiding, and future cloning of resistance genes from oat.

Archetypes of inflorescence: genome-wide association networks of panicle morphometric, growth, and disease variables in a multiparent oat population

There is limited information regarding the morphometric relationships of panicle traits in oat (Avena sativa) and their contribution to phenology and growth, physiology, and pathology traits important for yield. To model panicle growth and development and identify genomic regions associated with corresponding traits, 10 diverse spring oat mapping populations (n = 2,993) were evaluated in the field and 9 genotyped via genotyping-by-sequencing. Representative panicles from all progeny individuals, parents, and check lines were scanned, and images were analyzed using manual and automated techniques, resulting in over 60 unique panicle, rachis, and spikelet variables. Spatial modeling and days to heading were used to account for environmental and phenological variances, respectively. Panicle variables were intercorrelated, providing reproducible archetypal and growth models. Notably, adult plant resistance for oat crown rust was most prominent for taller, stiff stalked plants having a more open panicle structure. Within and among family variance for panicle traits reflected the moderate-to-high heritability and mutual genome-wide associations (hotspots) with numerous high-effect loci. Candidate genes and potential breeding applications are discussed. This work adds to the growing genetic resources for oat and provides a unique perspective on the genetic basis of panicle architecture in cereal crops.

A chromosome-level, fully phased genome assembly of the oat crown rust fungus Puccinia coronata f. sp. avenae: a resource to enable comparative genomics in the cereal rusts

Advances in sequencing technologies as well as development of algorithms and workflows have made it possible to generate fully phased genome references for organisms with nonhaploid genomes such as dikaryotic rust fungi. To enable discovery of pathogen effectors and further our understanding of virulence evolution, we generated a chromosome-scale assembly for each of the 2 nuclear genomes of the oat crown rust pathogen, Puccinia coronata f. sp. avenae (Pca). This resource complements 2 previously released partially phased genome references of Pca, which display virulence traits absent in the isolate of historic race 203 (isolate Pca203) which was selected for this genome project. A fully phased, chromosome-level reference for Pca203 was generated using PacBio reads and Hi-C data and a recently developed pipeline named NuclearPhaser for phase assignment of contigs and phase switch correction. With 18 chromosomes in each haplotype and a total size of 208.10 Mbp, Pca203 has the same number of chromosomes as other cereal rust fungi such as Puccinia graminis f. sp. tritici and Puccinia triticina, the causal agents of wheat stem rust and wheat leaf rust, respectively. The Pca203 reference marks the third fully phased chromosome-level assembly of a cereal rust to date. Here, we demonstrate that the chromosomes of these 3 Puccinia species are syntenous and that chromosomal size variations are primarily due to differences in repeat element content.

Identification of Candidate Susceptibility Genes to Puccinia graminis f. sp. tritici in Wheat

Wheat stem rust disease caused by Puccinia graminis f. sp. tritici (Pgt) is a global threat to wheat production. Fast evolving populations of Pgt limit the efficacy of plant genetic resistance and constrain disease management strategies. Understanding molecular mechanisms that lead to rust infection and disease susceptibility could deliver novel strategies to deploy crop resistance through genetic loss of disease susceptibility. We used comparative transcriptome-based and orthology-guided approaches to characterize gene expression changes associated with Pgt infection in susceptible and resistant Triticum aestivum genotypes as well as the non-host Brachypodium distachyon. We targeted our analysis to genes with differential expression in T. aestivum and genes suppressed or not affected in B. distachyon and report several processes potentially linked to susceptibility to Pgt, such as cell death suppression and impairment of photosynthesis. We complemented our approach with a gene co-expression network analysis to identify wheat targets to deliver resistance to Pgt through removal or modification of putative susceptibility genes.

Increased virulence of Puccinia coronata f. sp.avenae populations through allele frequency changes at multiple putative Avr loci

Pathogen populations are expected to evolve virulence traits in response to resistance deployed in agricultural settings. However, few temporal datasets have been available to characterize this process at the population level. Here, we examined two temporally separated populations of Puccinia coronata f. sp. avenae (Pca), which causes crown rust disease in oat (Avena sativa) sampled from 1990 to 2015. We show that a substantial increase in virulence occurred from 1990 to 2015 and this was associated with a genetic differentiation between populations detected by genome-wide sequencing. We found strong evidence for genetic recombination in these populations, showing the importance of the alternate host in generating genotypic variation through sexual reproduction. However, asexual expansion of some clonal lineages was also observed within years. Genome-wide association analysis identified seven Avr loci associated with virulence towards fifteen Pc resistance genes in oat and suggests that some groups of Pc genes recognize the same pathogen effectors. The temporal shift in virulence patterns in the Pca populations between 1990 and 2015 is associated with changes in allele frequency in these genomic regions. Nucleotide diversity patterns at a single Avr locus corresponding to Pc38, Pc39, Pc55, Pc63, Pc70, and Pc71 showed evidence of a selective sweep associated with the shift to virulence towards these resistance genes in all 2015 collected isolates.

Epigenetic regulation of gene expression improves Fusarium head blight resistance in durum wheat

Eight advanced durum-breeding lines were treated with 5-methyl-azacytidine to test the feasibility of generating sources of Fusarium head blight (FHB) resistance. Of the 800 treated seeds, 415 germinated and were advanced up to four (M4) generations by selfing. Thirty-two of the resulting 415 M4 lines were selected following preliminary screening and were further tested for FHB resistance for three years at two field locations, and in the greenhouse. Five of the 32 M4 lines showed less than 30% disease severity, as compared to the parental lines and susceptible checks. Fusarium-damaged kernels and deoxynivalenol analyses supported the findings of the field and greenhouse disease assessments. Two of the most resistant M4 lines were crossed to a susceptible parent, advanced to third generation (BC₁:F₃) and were tested for stability and inheritance of the resistance. About, one third of the BC₁:F₃ lines showed FHB resistance similar to their M4 parents. The overall methylation levels (%) were compared using FASTmC method, which did not show a significant difference between M4 and parental lines. However, transcriptome analysis of one M4 line revealed significant number of differentially expressed genes related to biosynthesis of secondary metabolites, MAPK signaling, photosynthesis, starch and sucrose metabolism, plant hormone signal transduction and plant-pathogen interaction pathways, which may have helped in improved FHB resistance.

Detection of Race-Specific Resistance Against Puccinia coronata f. sp. avenae in Brachypodium Species

Oat crown rust caused by Puccinia coronata f. sp. avenae is the most destructive foliar disease of cultivated oat. Characterization of genetic factors controlling resistance responses to Puccinia coronata f. sp. avenae in nonhost species could provide new resources for developing disease protection strategies in oat. We examined symptom development and fungal colonization levels of a collection of Brachypodium distachyon and B. hybridum accessions infected with three North American P. coronata f. sp. avenae isolates. Our results demonstrated that colonization phenotypes are dependent on both host and pathogen genotypes, indicating a role for race-specific responses in these interactions. These responses were independent of the accumulation of reactive oxygen species. Expression analysis of several defense-related genes suggested that salicylic acid and ethylene-mediated signaling but not jasmonic acid are components of resistance reaction to P. coronata f. sp. avenae. Our findings provide the basis to conduct a genetic inheritance study to examine whether effector-triggered immunity contributes to nonhost resistance to P. coronata f. sp. avenae in Brachypodium spp.

Comparative analysis of transcriptome in two wheat genotypes with contrasting levels of drought tolerance

Drought tolerance is a complex trait that is governed by multiple genes. The study presents differential transcriptome analysis between drought-tolerant (Triticum aestivum Cv. C306) and drought-sensitive (Triticum aestivum Cv. WL711) genotypes, using Affymetrix GeneChip® Wheat Genome Array. Both genotypes exhibited diverse global transcriptional responses under control and drought conditions. Pathway analysis suggested significant induction or repression of genes involved in secondary metabolism, nucleic acid synthesis, protein synthesis, and transport in C306, as compared to WL711. Significant up- and downregulation of transcripts for enzymes, hormone metabolism, and stress response pathways were observed in C306 under drought. The elevated expression of plasma membrane intrinsic protein 1 and downregulation of late embryogenesis abundant in the leaf tissues could play an important role in delayed wilting in C306. The other regulatory genes such as MT, FT, AP2, SKP1, ABA2, ARF6, WRKY6, AOS, and LOX2 are involved in defense response in C306 genotype. Additionally, transcripts with unknown functions were identified as differentially expressed, which could participate in drought responses.

A Novel Approach to Assess Salt Stress Tolerance in Wheat Using Hyperspectral Imaging

Salinity stress has significant adverse effects on crop productivity and yield. The primary goal of this study was to quantitatively rank salt tolerance in wheat using hyperspectral imaging. Four wheat lines were assayed in a hydroponic system with control and salt treatments (0 and 200 mM NaCl). Hyperspectral images were captured one day after salt application when there were no visual symptoms. Subsequent to necessary preprocessing tasks, two endmembers, each representing one of the treatment, were identified in each image using successive volume maximization. To simplify image analysis and interpretation, similarity of all pixels to the salt endmember was calculated by a technique proposed in this study, referred to as vector-wise similarity measurement. Using this approach allowed high-dimensional hyperspectral images to be reduced to one-dimensional gray-scale images while retaining all relevant information. Two methods were then utilized to analyze the gray-scale images: minimum difference of pair assignments and Bayesian method. The rankings of both methods were similar and consistent with the expected ranking obtained by conventional phenotyping experiments and historical evidence of salt tolerance. This research highlights the application of machine learning in hyperspectral image analysis for phenotyping of plants in a quantitative, interpretable, and non-invasive manner.

High-resolution crossover mapping reveals similarities and differences of male and female recombination in maize

Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male–female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.

Sex-specific meiotic crossover (CO) landscapes have been identified in multiple species. Here, the authors show that male and female meioses in maize have similar CO landscapes, and differences between COs in the two sexes only exists in their location relative to transcription start sites and some chromatin marks.

Puccinia coronata f. sp. avenae: a threat to global oat production

Puccinia coronata f. sp. avenae (Pca) causes crown rust disease in cultivated and wild oat (Avena spp.). The significant yield losses inflicted by this pathogen make crown rust the most devastating disease in the oat industry. Pca is a basidiomycete fungus with an obligate biotrophic lifestyle, and is classified as a typical macrocyclic and heteroecious fungus. The asexual phase in the life cycle of Pca occurs in oat, whereas the sexual phase takes place primarily in Rhamnus species as the alternative host. Epidemics of crown rust happens in areas with warm temperatures (20-25 °C) and high humidity. Infection by the pathogen leads to plant lodging and shrivelled grain of poor quality. Disease symptoms: Infection of susceptible oat varieties gives rise to orange-yellow round to oblong uredinia (pustules) containing newly formed urediniospores. Pustules vary in size and can be larger than 5 mm in length. Infection occurs primarily on the surfaces of leaves, although occasional symptoms develop in the oat leaf sheaths and/or floral structures, such as awns. Symptoms in resistant oat varieties vary from flecks to small pustules, typically accompanied by chlorotic halos and/or necrosis. The pycnial and aecial stages are mostly present in the leaves of Rhamnus species, but occasionally symptoms can also be observed in petioles, young stems and floral structures. Aecial structures display a characteristic hypertrophy and can differ in size, occasionally reaching more than 5 mm in diameter. Taxonomy: Pca belongs to the kingdom Fungi, phylum Basidiomycota, class Pucciniomycetes, order Pucciniales and family Pucciniaceae. Host range: Puccinia coronata sensu lato can infect 290 species of grass hosts. Pca is prevalent in all oat-growing regions and, compared with other cereal rusts, displays a broad telial host range. The most common grass hosts of Pca include cultivated hexaploid oat (Avena sativa) and wild relatives, such as bluejoint grass, perennial ryegrass and fescue. Alternative hosts include several species of Rhamnus, with R. cathartica (common buckthorn) as the most important alternative host in Europe and North America.

CONTROL

Most crown rust management strategies involve the use of rust-resistant crop varieties and the application of fungicides. The attainment of the durability of resistance against Pca is difficult as it is a highly variable pathogen with a great propensity to overcome the genetic resistance of varieties. Thus, adult plant resistance is often exploited in oat breeding programmes to develop new crown rust-resistant varieties. Useful website: https://www.ars.usda.gov/midwest-area/st-paul-mn/cereal-disease-lab/docs/cereal-rusts/race-surveys/.

Identification, introgression, and molecular marker genetic analysis and selection of a highly effective novel oat crown rust resistance from diploid oat, Avena strigosa

Oat crown rust is one of the most damaging diseases of oat. We identified a new source of resistance and developed KASP and TaqMan markers for selection in breeding programs. A new highly effective resistance to oat crown rust (Puccinia coronata f. sp. avenae) was identified in the diploid oat Avena strigosa PI 258731 and introgressed into hexaploid cultivated oat. Young plants with this resistance show moderate susceptibility, whereas older plant tissues and adult plants are resistant with no virulent isolates encountered in over 8 years of testing. Resistance was incorporated into hexaploid oat by embryo rescue, colchicine chromosome doubling followed by backcrosses with a hexaploid parent, and selection for stable transmission of resistance. To mitigate flag leaf and panicle chlorosis/necrosis associated with the resistance, crosses were made with derived resistant lines to breeding lines of divergent parentage followed by selection. Subsequently, two F₂ sister lines, termed MNBT1020-1 and MNBT1021-1, were identified in which the chlorosis/necrosis was reduced. These two lines performed well in replicated multi-location state trials in 2015 and 2016 out-yielding all cultivar entries. Segregating F_2:3 plants resulting from crosses of MNBT lines to susceptible parents were genotyped with the oat 6K SNP array, and SNP loci with close linkage to the resistance were identified. KASP assays generated from linked SNPs showed accurate discrimination of the resistance in derivatives of the resistant MNBT lines crossed to susceptible breeding lines. A TaqMan marker was developed and correctly identified homozygous resistance in over 95% of 379 F₄ plants when rust was scored in F_4:5 plants in the field. Thus, a novel highly effective resistance and associated molecular markers are available for use in breeding, genetic analysis, and functional studies.

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. avenaeIMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae, which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae, resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae.

Optimization and Comparative Analysis of Plant Organellar DNA Enrichment Methods Suitable for Next-generation Sequencing

Plant organellar genomes contain large, repetitive elements that may undergo pairing or recombination to form complex structures and/or sub-genomic fragments. Organellar genomes also exist in admixtures within a given cell or tissue type (heteroplasmy), and an abundance of subtypes may change throughout development or when under stress (sub-stoichiometric shifting). Next-generation sequencing (NGS) technologies are required to obtain deeper understanding of organellar genome structure and function. Traditional sequencing studies use several methods to obtain organellar DNA: (1) If a large amount of starting tissue is used, it is homogenized and subjected to differential centrifugation and/or gradient purification. (2) If a smaller amount of tissue is used (i.e., if seeds, material, or space is limited), the same process is performed as in (1), followed by whole-genome amplification to obtain sufficient DNA. (3) Bioinformatics analysis can be used to sequence the total genomic DNA and to parse out organellar reads. All these methods have inherent challenges and tradeoffs. In (1), it may be difficult to obtain such a large amount of starting tissue; in (2), whole-genome amplification could introduce a sequencing bias; and in (3), homology between nuclear and organellar genomes could interfere with assembly and analysis. In plants with large nuclear genomes, it is advantageous to enrich for organellar DNA to reduce sequencing costs and sequence complexity for bioinformatics analyses. Here, we compare a traditional differential centrifugation method with a fourth method, an adapted CpG-methyl pulldown approach, to separate the total genomic DNA into nuclear and organellar fractions. Both methods yield sufficient DNA for NGS, DNA that is highly enriched for organellar sequences, albeit at different ratios in mitochondria and chloroplasts. We present the optimization of these methods for wheat leaf tissue and discuss major advantages and disadvantages of each approach in the context of sample input, protocol ease, and downstream application.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

The role of mitochondria in plant development and stress tolerance

Eukaryotic cells require orchestrated communication between nuclear and organellar genomes, perturbations in which are linked to stress response and disease in both animals and plants. In addition to mitochondria, which are found across eukaryotes, plant cells contain a second organelle, the plastid. Signaling both among the organelles (cytoplasmic) and between the cytoplasm and the nucleus (i.e. nuclear-cytoplasmic interactions (NCI)) is essential for proper cellular function. A deeper understanding of NCI and its impact on development, stress response, and long-term health is needed in both animal and plant systems. Here we focus on the role of plant mitochondria in development and stress response. We compare and contrast features of plant and animal mitochondrial genomes (mtDNA), particularly highlighting the large and highly dynamic nature of plant mtDNA. Plant-based tools are powerful, yet underutilized, resources for enhancing our fundamental understanding of NCI. These tools also have great potential for improving crop production. Across taxa, mitochondria are most abundant in cells that have high energy or nutrient demands as well as at key developmental time points. Although plant mitochondria act as integrators of signals involved in both development and stress response pathways, little is known about plant mtDNA diversity and its impact on these processes. In humans, there are strong correlations between particular mitotypes (and mtDNA mutations) and developmental differences (or disease). We propose that future work in plants should focus on defining mitotypes more carefully and investigating their functional implications as well as improving techniques to facilitate this research.

Fast-forward genetics by radiation hybrids to saturate the locus regulating nuclear-cytoplasmic compatibility in Triticum

The nuclear-encoded species cytoplasm specific (scs) genes control nuclear-cytoplasmic compatibility in wheat (genus Triticum). Alloplasmic cells, which have nucleus and cytoplasm derived from different species, produce vigorous and vital organisms only when the correct version of scs is present in their nucleus. In this study, bulks of in vivo radiation hybrids segregating for the scs phenotype have been genotyped by sequencing with over 1.9 million markers. The high marker saturation obtained for a critical region of chromosome 1D allowed identification of 3318 reads that mapped in close proximity of the scs. A novel in silico approach was deployed to extend these short reads to sequences of up to 70 Kb in length and identify candidate open reading frames (ORFs). Markers were developed to anchor the short contigs containing ORFs to a radiation hybrid map of 650 individuals with resolution of 288 Kb. The region containing the scs locus was narrowed to a single Bacterial Artificial Chromosome (BAC) contig of Aegilops tauschii. Its sequencing and assembly by nano-mapping allowed rapid identification of a rhomboid gene as the only ORF existing within the refined scs locus. Resequencing of this gene from multiple germplasm sources identified a single nucleotide mutation, which gives rise to a functional amino acid change. Gene expression characterization revealed that an active copy of this rhomboid exists on all homoeologous chromosomes of wheat, and depending on the specific cytoplasm each copy is preferentially expressed. Therefore, a new methodology was applied to unique genetic stocks to rapidly identify a strong candidate gene for the control of nuclear-cytoplasmic compatibility in wheat.

Dissection of Genetic Factors underlying Wheat Kernel Shape and Size in an Elite × Nonadapted Cross using a High Density SNP Linkage Map

Wheat kernel shape and size has been under selection since early domestication. Kernel morphology is a major consideration in wheat breeding, as it impacts grain yield and quality. A population of 160 recombinant inbred lines (RIL), developed using an elite (ND 705) and a nonadapted genotype (PI 414566), was extensively phenotyped in replicated field trials and genotyped using Infinium iSelect 90K assay to gain insight into the genetic architecture of kernel shape and size. A high density genetic map consisting of 10,172 single nucleotide polymorphism (SNP) markers, with an average marker density of 0.39 cM/marker, identified a total of 29 genomic regions associated with six grain shape and size traits; ∼80% of these regions were associated with multiple traits. The analyses showed that kernel length (KL) and width (KW) are genetically independent, while a large number (∼59%) of the quantitative trait loci (QTL) for kernel shape traits were in common with genomic regions associated with kernel size traits. The most significant QTL was identified on chromosome 4B, and could be an ortholog of major rice grain size and shape gene or . Major and stable loci also were identified on the homeologous regions of Group 5 chromosomes, and in the regions of (6A) and (7A) genes. Both parental genotypes contributed equivalent positive QTL alleles, suggesting that the nonadapted germplasm has a great potential for enhancing the gene pool for grain shape and size. This study provides new knowledge on the genetic dissection of kernel morphology, with a much higher resolution, which may aid further improvement in wheat yield and quality using genomic tools.

Novel nuclear-cytoplasmic interaction in wheat (Triticum aestivum) induces vigorous plants

Interspecific hybridization can be considered an accelerator of evolution, otherwise a slow process, solely dependent on mutation and recombination. Upon interspecific hybridization, several novel interactions between nuclear and cytoplasmic genomes emerge which provide additional sources of diversity. The magnitude and essence of intergenomic interactions between nuclear and cytoplasmic genomes remain unknown due to the direction of many crosses. This study was conducted to address the role of nuclear-cytoplasmic interactions as a source of variation upon hybridization. Wheat (Triticum aestivum) alloplasmic lines carrying the cytoplasm of Aegilops mutica along with an integrated approach utilizing comparative quantitative trait locus (QTL) and epigenome analysis were used to dissect this interaction. The results indicate that cytoplasmic genomes can modify the magnitude of QTL controlling certain physiological traits such as dry matter weight. Furthermore, methylation profiling analysis detected eight polymorphic regions affected by the cytoplasm type. In general, these results indicate that novel nuclear-cytoplasmic interactions can potentially trigger an epigenetic modification cascade in nuclear genes which eventually change the genetic network controlling physiological traits. These modified genetic networks can serve as new sources of variation to accelerate the evolutionary process. Furthermore, this variation can synthetically be produced by breeders in their programs to develop epigenomic-segregating lines.

Dissecting Plant Chromosomes by the Use of Ionizing Radiation

Radiation treatment of genomes is used to generate chromosome breaks for numerous applications. This protocol describes the preparation of seeds and the determination of the optimal level of irradiation dosage for the creation of a radiation hybrid (RH) population. These RH lines can be used to generate high-resolution physical maps for the assembly of sequenced genomes as well as the fine mapping of genes. This procedure can also be used for mutation breeding and forward/reverse genetics.

Toward a better understanding of the genomic region harboring Fusarium head blight resistance QTL Qfhs.ndsu-3AS in durum wheat

New molecular markers were developed and mapped to the FHB resistance QTL region in high resolution. Micro-collinearity of the QTL region with rice and Brachypodium was revealed for a better understanding of the genomic region. The wild emmer wheat (Triticum dicoccoides)-derived Fusarium head blight (FHB) resistance quantitative trait locus (QTL) Qfhs.ndsu-3AS previously mapped to the short arm of chromosome 3A (3AS) in a population of recombinant inbred chromosome lines (RICLs). This study aimed to attain a better understanding of the genomic region harboring Qfhs.ndsu-3AS and to improve the utility of the QTL in wheat breeding. Micro-collinearity of the QTL region with rice chromosome 1 and Brachypodium chromosome 2 was identified and used for marker development in saturation mapping. A total of 42 new EST-derived sequence tagged site (STS) and simple sequence repeat (SSR) markers were developed and mapped to the QTL and nearby regions on 3AS. Further comparative analysis revealed a complex collinearity of the 3AS genomic region with their collinear counterparts of rice and Brachypodium. Fine mapping of the QTL region resolved five co-segregating markers (Xwgc1186/Xwgc716/Xwgc1143/Xwgc501/Xwgc1204) into three distinct loci proximal to Xgwm2, a marker previously reported to be closely linked to the QTL. Four other markers (Xwgc1226, Xwgc510, Xwgc1296, and Xwgc1301) mapped farther proximal to the above markers in the QTL region with a higher resolution. Five homozygous recombinants with shortened T. dicoccoides chromosomal segments in the QTL region were recovered by molecular marker analysis and evaluated for FHB resistance. Qfhs.ndsu-3AS was positioned to a 5.2 cM interval flanked by the marker Xwgc501 and Xwgc510. The recombinants containing Qfhs.ndsu-3AS and new markers defining the QTL will facilitate utilization of this resistance source in wheat breeding.

Gene Evolutionary Trajectories and GC Patterns Driven by Recombination in Zea mays

Recombination occurring during meiosis is critical for creating genetic variation and plays an essential role in plant evolution. In addition to creating novel gene combinations, recombination can affect genome structure through altering GC patterns. In maize (Zea mays) and other grasses, another intriguing GC pattern exists. Maize genes show a bimodal GC content distribution that has been attributed to nucleotide bias in the third, or wobble, position of the codon. Recombination may be an underlying driving force given that recombination sites are often associated with high GC content. Here we explore the relationship between recombination and genomic GC patterns by comparing GC gene content at each of the three codon positions (GC₁, GC₂, and GC₃, collectively termed GC_x) to instances of a variable GC-rich motif that underlies double strand break (DSB) hotspots and to meiocyte-specific gene expression. Surprisingly, GC_x bimodality in maize cannot be fully explained by the codon wobble hypothesis. High GC_x genes show a strong overlap with the DSB hotspot motif, possibly providing a mechanism for the high evolutionary rates seen in these genes. On the other hand, genes that are turned on in meiosis (early prophase I) are biased against both high GC_x genes and genes with the DSB hotspot motif, possibly allowing important meiotic genes to avoid DSBs. Our data suggests a strong link between the GC-rich motif underlying DSB hotspots and high GC_x genes.

Radiation hybrid maps of the D-genome of Aegilops tauschii and their application in sequence assembly of large and complex plant genomes

BACKGROUND

The large and complex genome of bread wheat (Triticum aestivum L., ~17 Gb) requires high resolution genome maps with saturated marker scaffolds to anchor and orient BAC contigs/ sequence scaffolds for whole genome assembly. Radiation hybrid (RH) mapping has proven to be an excellent tool for the development of such maps for it offers much higher and more uniform marker resolution across the length of the chromosome compared to genetic mapping and does not require marker polymorphism per se, as it is based on presence (retention) vs. absence (deletion) marker assay.

METHODS

In this study, a 178 line RH panel was genotyped with SSRs and DArT markers to develop the first high resolution RH maps of the entire D-genome of Ae. tauschii accession AL8/78. To confirm map order accuracy, the AL8/78-RH maps were compared with:1) a DArT consensus genetic map constructed using more than 100 bi-parental populations, 2) a RH map of the D-genome of reference hexaploid wheat 'Chinese Spring', and 3) two SNP-based genetic maps, one with anchored D-genome BAC contigs and another with anchored D-genome sequence scaffolds. Using marker sequences, the RH maps were also anchored with a BAC contig based physical map and draft sequence of the D-genome of Ae. tauschii.

RESULTS

A total of 609 markers were mapped to 503 unique positions on the seven D-genome chromosomes, with a total map length of 14,706.7 cR. The average distance between any two marker loci was 29.2 cR which corresponds to 2.1 cM or 9.8 Mb. The average mapping resolution across the D-genome was estimated to be 0.34 Mb (Mb/cR) or 0.07 cM (cM/cR). The RH maps showed almost perfect agreement with several published maps with regard to chromosome assignments of markers. The mean rank correlations between the position of markers on AL8/78 maps and the four published maps, ranged from 0.75 to 0.92, suggesting a good agreement in marker order. With 609 mapped markers, a total of 2481 deletions for the whole D-genome were detected with an average deletion size of 42.0 Mb. A total of 520 markers were anchored to 216 Ae. tauschii sequence scaffolds, 116 of which were not anchored earlier to the D-genome.

CONCLUSION

This study reports the development of first high resolution RH maps for the D-genome of Ae. tauschii accession AL8/78, which were then used for the anchoring of unassigned sequence scaffolds. This study demonstrates how RH mapping, which offered high and uniform resolution across the length of the chromosome, can facilitate the complete sequence assembly of the large and complex plant genomes.

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.

Radiation Hybrid Map of Barley Chromosome 3H

Assembly of the barley (Hordeum vulgare L.) genome is complicated by its large size (5.1 Gb) and proportion of repetitive elements (84%). This process is facilitated by high resolution maps for aligning bacterial artificial chromosome (BAC) contigs along chromosomes. Available genetic maps, however, do not provide accurate information on the physical position of a large portion of the genome located in recombination-poor regions. Radiation hybrid (RH) mapping is an alternative approach, which is based on radiation-induced deletions along the length of chromosomes. In this study, the first RH map for barley chromosome 3H was developed. In total, 373 in vivo RH lines were generated by irradiating wheat (Triticum aestivum L.)-barley chromosome 3H addition lines and crossing them to a normal wheat cultivar. Each RH informative line (containing deletions) had, on average, three deletions. The induced deletion size varied from 36.58 Kb to 576.00 Mb, with an average length of 52.42 Mb. This initial chromosome 3H radiation hybrid (3H-RH) map had a 9.53× higher resolution than an analogous genetic map, reaching a maximum of >262.40× resolution in regions around the centromere. The final RH map was 3066.1 cR in length, with a 0.76 Mb resolution. It was estimated that the map resolution can be improved to an average of 30.34 Kb by saturating the 3H-RH map with molecular markers. The generated RH panel enabled alignment of BAC and sequenced contigs as small as 1.50 Kb in size. The high resolution and the coverage of poor-recombination regions make RH maps an ideal resource for barley genome assembly, as well as other genetic studies.

Reliable Radiation Hybrid Maps: An Efficient Scalable Clustering-Based Approach

The process of mapping markers from radiation hybrid mapping (RHM) experiments is equivalent to the traveling salesman problem and, thereby, has combinatorial complexity. As an additional problem, experiments typically result in some unreliable markers that reduce the overall quality of the map. We propose a clustering approach for addressing both problems efficiently by eliminating unreliable markers without the need for mapping the complete set of markers. Traditional approaches for eliminating markers use resampling of the full data set, which has an even higher computational complexity than the original mapping problem. In contrast, the proposed approach uses a divide-and-conquer strategy to construct framework maps based on clusters that exclude unreliable markers. Clusters are ordered using parallel processing and are then combined to form the complete map. We present three algorithms that explore the trade-off between the number of markers included in the map and placement accuracy. Using an RHM data set of the human genome, we compare the framework maps from our proposed approaches with published physical maps and with the results of using the Carthagene tool. Overall, our approaches have a very low computational complexity and produce solid framework maps with good chromosome coverage and high agreement with the physical map marker order.

Analysis of ATP6 sequence diversity in the Triticum-Aegilops species group reveals the crucial role of rearrangement in mitochondrial genome evolution

Mutation and chromosomal rearrangements are the two main forces of increasing genetic diversity for natural selection to act upon, and ultimately drive the evolutionary process. Although genome evolution is a function of both forces, simultaneously, the ratio of each can be varied among different genomes and genomic regions. It is believed that in plant mitochondrial genome, rearrangements play a more important role than point mutations, but relatively few studies have directly addressed this phenomenon. To address this issue, we isolated and sequenced the ATP6-1 and ATP6-2 genes from 46 different euplasmic and alloplasmic wheat lines. Four different ATP6-1 orthologs were detected, two of them reported for the first time. Expression analysis revealed that all four orthologs are transcriptionally active. Results also indicated that both point mutation and genomic rearrangement are involved in the evolution of ATP6. However, rearrangement is the predominant force that triggers drastic variation. Data also indicated that speciation of domesticated wheat cultivars were simultaneous with the duplication of this gene. These results directly support the notion that rearrangement plays a significant role in driving plant mitochondrial genome evolution.

Mitochondrial dynamics and the cell cycle

Nuclear-mitochondrial (NM) communication impacts many aspects of plant development including vigor, sterility, and viability. Dynamic changes in mitochondrial number, shape, size, and cellular location takes place during the cell cycle possibly impacting the process itself and leading to distribution of this organelle into daughter cells. The genes that underlie these changes are beginning to be identified in model plants such as Arabidopsis. In animals disruption of the drp1 gene, a homolog to the plant drp3A and drp3B, delays mitochondrial division. This mutation results in increased aneuploidy due to chromosome mis-segregation. It remains to be discovered if a similar outcome is observed in plants. Alloplasmic lines provide an opportunity to understand the communication between the cytoplasmic organelles and the nucleus. Examples of studies in these lines, especially from the extensive collection in wheat, point to the role of mitochondria in chromosome movement, pollen fertility and other aspects of development.

Accelerated evolution of the mitochondrial genome in an alloplasmic line of durum wheat

BACKGROUND

Wheat is an excellent plant species for nuclear mitochondrial interaction studies due to availability of large collection of alloplasmic lines. These lines exhibit different vegetative and physiological properties than their parents. To investigate the level of sequence changes introduced into the mitochondrial genome under the alloplasmic condition, three mitochondrial genomes of the Triticum-Aegilops species were sequenced: 1) durum alloplasmic line with the Ae. longissima cytoplasm that carries the T. turgidum nucleus designated as (lo) durum, 2) the cytoplasmic donor line, and 3) the nuclear donor line.

RESULTS

The mitochondrial genome of the T. turgidum was 451,678 bp in length with high structural and nucleotide identity to the previously characterized T. aestivum genome. The assembled mitochondrial genome of the (lo) durum and the Ae. longissima were 431,959 bp and 399,005 bp in size, respectively. The high sequence coverage for all three genomes allowed analysis of heteroplasmy within each genome. The mitochondrial genome structure in the alloplasmic line was genetically distant from both maternal and paternal genomes. The alloplasmic durum and the Ae. longissima carry the same versions of atp6, nad6, rps19-p, cob and cox2 exon 2 which are different from the T. turgidum parent. Evidence of paternal leakage was also observed by analyzing nad9 and orf359 among all three lines. Nucleotide search identified a number of open reading frames, of which 27 were specific to the (lo) durum line.

CONCLUSIONS

Several heteroplasmic regions were observed within genes and intergenic regions of the mitochondrial genomes of all three lines. The number of rearrangements and nucleotide changes in the mitochondrial genome of the alloplasmic line that have occurred in less than half a century was significant considering the high sequence conservation between the T. turgidum and the T. aestivum that diverged from each other 10,000 years ago. We showed that the changes in genes were not limited to paternal leakage but were sufficiently significant to suggest that other mechanisms, such as recombination and mutation, were responsible. The newly formed ORFs, differences in gene sequences and copy numbers, heteroplasmy, and substoichiometric changes show the potential of the alloplasmic condition to accelerate evolution towards forming new mitochondrial genomes.

Deletion mutagenesis identifies a haploinsufficient role for γ-zein in opaque2 endosperm modification

Quality Protein Maize (QPM) is a hard kernel variant of the high-lysine mutant opaque2. Using γ-irradiation, we created opaque QPM variants to identify opaque2 modifier genes and to investigate deletion mutagenesis combined with Illumina sequencing as a maize (Zea mays) functional genomics tool. A K0326Y QPM deletion mutant was null for the 27- and 50-kD γ-zeins and abolished vitreous endosperm formation. Illumina exon and RNA sequencing revealed a 1.2-megabase pair deletion encompassing the 27- and 50-kD γ-zein genes on chromosome 7 and a deletion of at least 232 kb on chromosome 9. Protein body number was reduced by over 90%, while protein body size is similar to the wild type. Kernels hemizygous for the γ-zein deletion had intermediate 27- and 50-kD γ-zein levels and were semivitreous, indicating haploinsufficiency of these gene products in opaque2 endosperm modification. The γ-zein deletion further increased lysine in QPM in its homozygous and hemizygous states. This work identifies 27-kD γ-zein as an opaque2 modifier gene within the largest QPM quantitative trait locus and may suggest the 50-kD γ-zein also contributes to this quantitative trait locus. It further demonstrates that genome-wide deletions in nonreference maize lines can be identified through a combination of assembly of Illumina reads against the B73 genome and integration of RNA sequencing data.

Radiation hybrid QTL mapping of Tdes2 involved in the first meiotic division of wheat

Since the dawn of wheat cytogenetics, chromosome 3B has been known to harbor a gene(s) that, when removed, causes chromosome desynapsis and gametic sterility. The lack of natural genetic diversity for this gene(s) has prevented any attempt to fine map and further characterize it. Here, gamma radiation treatment was used to create artificial diversity for this locus. A total of 696 radiation hybrid lines were genotyped with a custom mini array of 140 DArT markers, selected to evenly span the whole 3B chromosome. The resulting map spanned 2,852 centi Ray with a calculated resolution of 0.384 Mb. Phenotyping for the occurrence of meiotic desynapsis was conducted by measuring the level of gametic sterility as seeds produced per spikelet and pollen viability at booting. Composite interval mapping revealed a single QTL with LOD of 16.2 and r (2) of 25.6 % between markers wmc326 and wPt-8983 on the long arm of chromosome 3B. By independent analysis, the location of the QTL was confirmed to be within the deletion bin 3BL7-0.63-1.00 and to correspond to a single gene located ~1.4 Mb away from wPt-8983. The meiotic behavior of lines lacking this gene was characterized cytogenetically to reveal striking similarities with mutants for the dy locus, located on the syntenic chromosome 3 of maize. This represents the first example to date of employing radiation hybrids for QTL analysis. The success achieved by this approach provides an ideal starting point for the final cloning of this interesting gene involved in meiosis of cereals.

A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus

The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."

A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus

The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."

Wheat Zapper: a flexible online tool for colinearity studies in grass genomes

In the course of evolution, the genomes of grasses have maintained an observable degree of gene order conservation. The information available for already sequenced genomes can be used to predict the gene order of nonsequenced species by means of comparative colinearity studies. The "Wheat Zapper" application presented here performs on-demand colinearity analysis between wheat, rice, Sorghum, and Brachypodium in a simple, time efficient, and flexible manner. This application was specifically designed to provide plant scientists with a set of tools, comprising not only synteny inference, but also automated primer design, intron/exon boundaries prediction, visual representation using the graphic tool Circos 0.53, and the possibility of downloading FASTA sequences for downstream applications. Quality of the "Wheat Zapper" prediction was confirmed against the genome of maize, with good correlation (r > 0.83) observed between the gene order predicted on the basis of synteny and their actual position on the genome. Further, the accuracy of "Wheat Zapper" was calculated at 0.65 considering the "Genome Zipper" application as the "gold" standard. The differences between these two tools are amply discussed, making the point that "Wheat Zapper" is an accurate and reliable on-demand tool that is sure to benefit the cereal scientific community. The Wheat Zapper is available at http://wge.ndsu.nodak.edu/wheatzapper/ .

PIECE: a database for plant gene structure comparison and evolution

Gene families often show degrees of differences in terms of exon–intron structures depending on their distinct evolutionary histories. Comparative analysis of gene structures is important for understanding their evolutionary and functional relationships within plant species. Here, we present a comparative genomics database named PIECE (http://wheat.pw.usda.gov/piece) for Plant Intron and Exon Comparison and Evolution studies. The database contains all the annotated genes extracted from 25 sequenced plant genomes. These genes were classified based on Pfam motifs. Phylogenetic trees were pre-constructed for each gene category. PIECE provides a user-friendly interface for different types of searches and a graphical viewer for displaying a gene structure pattern diagram linked to the resulting bootstrapped dendrogram for each gene family. The gene structure evolution of orthologous gene groups was determined using the GLOOME, Exalign and GECA software programs that can be accessed within the database. PIECE also provides a web server version of the software, GSDraw, for drawing schematic diagrams of gene structures. PIECE is a powerful tool for comparing gene sequences and provides valuable insights into the evolution of gene structure in plant genomes.

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 studies.

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 RH1s (AABBD), which allows the simultaneous mapping of the whole D-genome. A panel of 1,510 RH1 plants was obtained, of which 592 plants were generated from the mature RH1 seeds, and 918 plants were rescued through embryo culture due to poor germination (<3%) of mature RH1 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.

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 RH1s (AABBD), which allows the simultaneous mapping of the whole D-genome. A panel of 1,510 RH1 plants was obtained, of which 592 plants were generated from the mature RH1 seeds, and 918 plants were rescued through embryo culture due to poor germination (<3%) of mature RH1 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.

A molecular linkage map of cultivated oat

A molecular linkage map of cultivated oat composed of 561 loci has been developed using 71 recombinant inbred lines from a cross between Avena byzantina cv. Kanota and A. sativa cv. Ogle. The loci are mainly restriction fragment length polymorphisms detected by oat cDNA clones from leaf, endosperm, and root tissue, as well as by barley leaf cDNA clones. The loci form 38 linkage groups ranging in size from 0.0 to 122.1 cM (mean, 39 cM) and consist of 2-51 loci each (mean, 14). Twenty-nine loci remain unlinked. The current map size is 1482 cM and the total size, on the basis of the number of unlinked loci, is estimated to be 2932.0 cM. This indicates that this map covers at least 50% of the cultivated oat genome. Comparisons with an A-genome diploid oat map and between linkage groups exhibiting homoeology to each other indicate that several major chromosomal rearrangements exist in cultivated oat. This map provides a tool for marker-assisted selection, quantitative trait loci analyses, and studies of genome organization in oat.

Population structure and linkage disequilibrium in Lupinus albus L. germplasm and its implication for association mapping

White lupin (Lupinus albus L.) has been around since 300 B.C. and is recognized for its ability to grow on poor soils and application as green manure in addition to seed harvest. The seed has very high levels of protein (33-47 %) and oil (6-13 %). It also has many secondary metabolites that are potentially of nutraceutical value to animals and humans. Despite such a great potential, lupins role in modern agriculture began only in the twentieth century. Although a large collection of Lupinus germplasm accessions is available worldwide, rarely have they been genetically characterized. Additionally, scarce genomic resources in terms of recombinant populations and genome information have been generated for L. albus. With the advancement in association mapping methods, the natural populations have the potential to replace the recombinant populations in gene mapping and marker-trait associations. Therefore, we studied the genetic similarity, population structure and marker-trait association in a USDA germplasm collection for their current and future application in this crop improvement. A total of 122 PI (Plant Inventory) lines were screened with 18 AFLP primer pairs that generated 2,277 fragments. A subset of 892 polymorphic markers with MAF >0.05 (minor allele frequency) were used for association mapping. The cluster analysis failed to group accessions on the basis of their passport information, and a weak structure and low linkage disequilibrium (LD) were observed indicating the usefulness of the collection for association mapping. Moreover, we were also able to identify two markers (a p value of 1.53 × 10(-4) and 2.3 × 10(-4)) that explained 22.69 and 20.5 % of seed weight variation determined using R (LR) (2) . The implications of lack of geographic clustering, population structure, low LD and the ability of AFLP to map seed weight trait using association mapping and the usefulness of the PI collections in breeding programs are discussed.

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.

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.