Comparison of orthologous loci from small grass genomes Brachypodium and rice: implications for wheat genomics and grass genome annotation

Plant Metabolomics: An Overview of the Role of Primary and Secondary Metabolites against Different Environmental Stress Factors

Several environmental stresses, including biotic and abiotic factors, adversely affect the growth and development of crops, thereby lowering their yield. However, abiotic factors, e.g., drought, salinity, cold, heat, ultraviolet radiations (UVr), reactive oxygen species (ROS), trace metals (TM), and soil pH, are extremely destructive and decrease crop yield worldwide. It is expected that more than 50% of crop production losses are due to abiotic stresses. Moreover, these factors are responsible for physiological and biochemical changes in plants. The response of different plant species to such stresses is a complex phenomenon with individual features for several species. In addition, it has been shown that abiotic factors stimulate multi-gene responses by making modifications in the accumulation of the primary and secondary metabolites. Metabolomics is a promising way to interpret biotic and abiotic stress tolerance in plants. The study of metabolic profiling revealed different types of metabolites, e.g., amino acids, carbohydrates, phenols, polyamines, terpenes, etc, which are accumulated in plants. Among all, primary metabolites, such as amino acids, carbohydrates, lipids polyamines, and glycine betaine, are considered the major contributing factors that work as osmolytes and osmoprotectants for plants from various environmental stress factors. In contrast, plant-derived secondary metabolites, e.g., phenolics, terpenoids, and nitrogen-containing compounds (alkaloids), have no direct role in the growth and development of plants. Nevertheless, such metabolites could play a significant role as a defense by protecting plants from biotic factors such as herbivores, insects, and pathogens. In addition, they can enhance the resistance against abiotic factors. Therefore, metabolomics practices are becoming essential and influential in plants by identifying different phytochemicals that are part of the acclimation responses to various stimuli. Hence, an accurate metabolome analysis is important to understand the basics of stress physiology and biochemistry. This review provides insight into the current information related to the impact of biotic and abiotic factors on variations of various sets of metabolite levels and explores how primary and secondary metabolites help plants in response to these stresses.

Systematic annotation of conservation states provides insights into regulatory regions in rice

Plant genomes contain a large fraction of noncoding sequences. The discovery and annotation of conserved noncoding sequences (CNSs) in plants is an ongoing challenge. Here we report the application of comparative genomics to systematically identify CNSs in 50 well-annotated Gramineae genomes using rice (Oryza sativa) as the reference. We conduct multiple-way whole-genome alignments to the rice genome. The rice genome is annotated as 20 conservation states (CSs) at single-nucleotide resolution using a multivariate hidden Markov model (ConsHMM) based on the multiple-genome alignments. Different states show distinct enrichments for various genomic features, and the conservation scores of CSs are highly correlated with the level of associated chromatin accessibility. We find that at least 33.5% of the rice genome is highly under selection, with more than 70% of the sequence lying outside of coding regions. A catalog of 855,366 regulatory CNSs is generated, and they significantly overlapped with putative active regulatory elements such as promoters, enhancers, and transcription factor binding sites. Collectively, our study provides a resource for elucidating functional noncoding regions of the rice genome and an evolutionary aspect of regulatory sequences in higher plants.

Cold-inducible promoter-driven knockdown of Brachypodium antifreeze proteins confers freezing and phytopathogen susceptibility

The model forage crop, Brachypodium distachyon, has a cluster of ice recrystallization inhibition (BdIRI) genes, which encode antifreeze proteins that function by adsorbing to ice crystals and inhibiting their growth. The genes were targeted for knockdown using a cold-induced promoter from rice (prOsMYB1R35) to drive miRNA. The transgenic lines showed no apparent pleiotropic developmental defects but had reduced antifreeze activity as assessed by assays for ice-recrystallization inhibition, thermal hysteresis, electrolyte leakage, and leaf infrared thermography. Strikingly, the number of cold-acclimated transgenic plants that survived freezing at -8°C was reduced by half or killed entirely, depending on the line, compared with cold-acclimated wild type plants. In addition, more leaf damage was apparent at subzero temperatures in knockdowns after infection with an ice nucleating pathogen, Pseudomonas syringae. Although antifreeze proteins have been studied for almost 60 years, this is the first unequivocal demonstration of their function by knockdown in any organism, and their dual contribution to freeze protection as well as pathogen susceptibility, independent of obvious developmental defects. These proteins are thus of potential interest in a wide range of biotechnological applications from cryopreservation, to frozen product additives, to the engineering of transgenic crops with enhanced pathogen and freezing tolerance.

The origin and evolution of a two-component system of paralogous genes encoding the centromeric histone CENH3 in cereals

BACKGROUND

The cereal family Poaceae is one of the largest and most diverse angiosperm families. The central component of centromere specification and function is the centromere-specific histone H3 (CENH3). Some cereal species (maize, rice) have one copy of the gene encoding this protein, while some (wheat, barley, rye) have two. We applied a homology-based approach to sequenced cereal genomes, in order to finally trace the mutual evolution of the structure of the CENH3 genes and the nearby regions in various tribes.

RESULTS

We have established that the syntenic group or the CENH3 locus with the CENH3 gene and the boundaries defined by the CDPK2 and bZIP genes first appeared around 50 Mya in a common ancestor of the subfamilies Bambusoideae, Oryzoideae and Pooideae. This locus came to Pooideae with one copy of CENH3 in the most ancient tribes Nardeae and Meliceae. The βCENH3 gene as a part of the locus appeared in the tribes Stipeae and Brachypodieae around 35-40 Mya. The duplication was accompanied by changes in the exon-intron structure. Purifying selection acts mostly on αCENH3s, while βCENH3s form more heterogeneous structures, in which clade-specific amino acid motifs are present. In barley species, the βCENH3 gene assumed an inverted orientation relative to αCENH3 and the CDPK2 gene was substituted with LHCB-l. As the evolution and domestication of plant species went on, the locus was growing in size due to an increasing distance between αCENH3 and βCENH3 because of a massive insertion of the main LTR-containing retrotransposon superfamilies, gypsy and copia, without any evolutionary preference on either of them. A comparison of the molecular structure of the locus in the A, B and D subgenomes of the hexaploid wheat T. aestivum showed that invasion by mobile elements and concomitant rearrangements took place in an independent way even in evolutionarily close species.

CONCLUSIONS

The CENH3 duplication in cereals was accompanied by changes in the exon-intron structure of the βCENH3 paralog. The observed general tendency towards the expansion of the CENH3 locus reveals an amazing diversity of ways in which different species implement the scenario described in this paper.

GWAS revealed effect of genotype × environment interactions for grain yield of Nebraska winter wheat

Background

Improving grain yield in cereals especially in wheat is a main objective for plant breeders. One of the main constrains for improving this trait is the G × E interaction (GEI) which affects the performance of wheat genotypes in different environments. Selecting high yielding genotypes that can be used for a target set of environments is needed. Phenotypic selection can be misleading due to the environmental conditions. Incorporating information from phenotypic and genomic analyses can be useful in selecting the higher yielding genotypes for a group of environments.

Results

A set of 270 F_3:6 wheat genotypes in the Nebraska winter wheat breeding program was tested for grain yield in nine environments. High genetic variation for grain yield was found among the genotypes. G × E interaction was also highly significant. The highest yielding genotype differed in each environment. The correlation for grain yield among the nine environments was low (0 to 0.43). Genome-wide association study revealed 70 marker traits association (MTAs) associated with increased grain yield. The analysis of linkage disequilibrium revealed 16 genomic regions with a highly significant linkage disequilibrium (LD). The candidate parents’ genotypes for improving grain yield in a group of environments were selected based on three criteria; number of alleles associated with increased grain yield in each selected genotype, genetic distance among the selected genotypes, and number of different alleles between each two selected parents.

Conclusion

Although G × E interaction was present, the advances in DNA technology provided very useful tools and analyzes. Such features helped to genetically select the highest yielding genotypes that can be used to cross grain production in a group of environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07308-0.

Genome-wide association study of agronomic traits in bread wheat reveals novel putative alleles for future breeding programs

BACKGROUND

Identification of loci for agronomic traits and characterization of their genetic architecture are crucial in marker-assisted selection (MAS). Genome-wide association studies (GWAS) have increasingly been used as potent tools in identifying marker-trait associations (MTAs). The introduction of new adaptive alleles in the diverse genetic backgrounds may help to improve grain yield of old or newly developed varieties of wheat to balance supply and demand throughout the world. Landraces collected from different climate zones can be an invaluable resource for such adaptive alleles.

RESULTS

GWAS was performed using a collection of 298 Iranian bread wheat varieties and landraces to explore the genetic basis of agronomic traits during 2016-2018 cropping seasons under normal (well-watered) and stressed (rain-fed) conditions. A high-quality genotyping by sequencing (GBS) dataset was obtained using either all original single nucleotide polymorphism (SNP, 10938 SNPs) or with additional imputation (46,862 SNPs) based on W7984 reference genome. The results confirm that the B genome carries the highest number of significant marker pairs in both varieties (49,880, 27.37%) and landraces (55,086, 28.99%). The strongest linkage disequilibrium (LD) between pairs of markers was observed on chromosome 2D (0.296). LD decay was lower in the D genome, compared to the A and B genomes. Association mapping under two tested environments yielded a total of 313 and 394 significant (-log10 P >3) MTAs for the original and imputed SNP data sets, respectively. Gene ontology results showed that 27 and 27.5% of MTAs of SNPs in the original set were located in protein-coding regions for well-watered and rain-fed conditions, respectively. While, for the imputed data set 22.6 and 16.6% of MTAs represented in protein-coding genes for the well-watered and rain-fed conditions, respectively.

CONCLUSIONS

Our finding suggests that Iranian bread wheat landraces harbor valuable alleles that are adaptive under drought stress conditions. MTAs located within coding genes can be utilized in genome-based breeding of new wheat varieties. Although imputation of missing data increased the number of MTAs, the fraction of these MTAs located in coding genes were decreased across the different sub-genomes.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Comparative transcriptome analysis provides key insights into seedling development in switchgrass (Panicum virgatum L.)

BACKGROUND

Switchgrass (Panicum virgatum L.), a warm-season perennial C4 plant, can be used as a forage plant, a soil and water conservation plant, a windbreak plant, and as a good source of biofuels and alternative energy with low planting costs. However, switchgrass exhibits low rates of seedling development compared to other crops, which means it is typically out-competed by weeds. There is a large variation in seedling development rates among different plantlets in switchgrass, which limits its usefulness for large-scale cultivation. Little is currently known about the molecular reasons for slow seedling growth.

RESULTS

Characterization of the seedling development process via growth indices indicated a relatively stagnant growth stage in switchgrass. A total of 678 differentially expressed genes (DEGs) were identified from the comparison of transcriptomes from slowly developed (sd) and rapidly developed (rd) switchgrass seedlings. Gene ontology and pathway enrichment analysis showed that DEGs were enriched in diterpenoid biosynthesis, thiamine metabolism, and circadian rhythm. Transcription factor enrichment and expression analyses showed MYB-related, bHLH and NAC family genes were essential for seedling growth. The transcriptome results were consistent with those of quantitative real-time polymerase chain reaction. Then, the expression profiles of maize and switchgrass were compared during seedling leaf development. A total of 128 DEGs that play key roles in seedling growth were aligned to maize genes. Transcriptional information and physiological indices suggested that several genes involved in the circadian rhythm, thiamine metabolism, energy metabolism, gibberellic acid biosynthesis, and signal transduction played important roles in seedling development.

CONCLUSIONS

The seedling development process of switchgrass was characterized, and the molecular differences between slowly developed and rapidly developed seedlings were discussed. This study provides new insights into the reasons for slow seedling development in switchgrass and will be useful for the genetic improvement of switchgrass and other crops.

The highly divergent Jekyll genes, required for sexual reproduction, are lineage specific for the related grass tribes Triticeae and Bromeae

Phylogenetically related groups of species contain lineage-specific genes that exhibit no sequence similarity to any genes outside the lineage. We describe here that the Jekyll gene, required for sexual reproduction, exists in two much diverged allelic variants, Jek1 and Jek3. Despite low similarity, the Jek1 and Jek3 proteins share identical signal peptides, conserved cysteine positions and direct repeats. The Jek1/Jek3 sequences are located at the same chromosomal locus and inherited in a monogenic Mendelian fashion. Jek3 has a similar expression as Jek1 and complements the Jek1 function in Jek1-deficient plants. Jek1 and Jek3 allelic variants were almost equally distributed in a collection of 485 wild and domesticated barley accessions. All domesticated barleys harboring the Jek1 allele belong to single haplotype J1-H1 indicating a genetic bottleneck during domestication. Domesticated barleys harboring the Jek3 allele consisted of three haplotypes. Jekyll-like sequences were found only in species of the closely related tribes Bromeae and Triticeae but not in other Poaceae. Non-invasive magnetic resonance imaging revealed intrinsic grain structure in Triticeae and Bromeae, associated with the Jekyll function. The emergence of Jekyll suggests its role in the separation of the Bromeae and Triticeae lineages within the Poaceae and identifies the Jekyll genes as lineage-specific.

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.

Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family

Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.

Components of Brachypodium distachyon resistance to nonadapted wheat stripe rust pathogens are simply inherited

Phytopathogens have a limited range of host plant species that they can successfully parasitise ie. that they are adapted for. Infection of plants by nonadapted pathogens often results in an active resistance response that is relatively poorly characterised because phenotypic variation in this response often does not exist within a plant species, or is too subtle for genetic dissection. In addition, complex polygenic inheritance often underlies these resistance phenotypes and mutagenesis often does not impact upon this resistance, presumably due to genetic or mechanistic redundancy. Here it is demonstrated that phenotypic differences in the resistance response of Brachypodium distachyon to the nonadapted wheat stripe rust pathogen Puccinia striiformis f. sp. tritici (Pst) are genetically tractable and simply inherited. Two dominant loci were identified on B. distachyon chromosome 4 that each reduce attempted Pst colonisation compared with sib and parent lines without these loci. One locus (Yrr1) is effective against diverse Australian Pst isolates and present in two B. distachyon mapping families as a conserved region that was reduced to 5 candidate genes by fine mapping. A second locus, Yrr2, shows Pst race-specificity and encodes a disease resistance gene family typically associated with host plant resistance. These data indicate that some components of resistance to nonadapted pathogens are genetically tractable in some instances and may mechanistically overlap with host plant resistance to avirulent adapted pathogens.

Silencing of a Unique Integrated Domain Nucleotide-Binding Leucine-Rich Repeat Gene in Wheat Abolishes Diuraphis noxia Resistance

Plants respond in a similar manner to aphid feeding as to pathogen attack. Diuraphis noxia is a specialist aphid, feeding only on selected grasses that include wheat, barley, and oats. The wheat-Diuraphis noxia interaction is characterized by responses very similar to those seen in wheat-pathogen interactions with none of the underlying resistance pathways and genes characterized yet. From wheat harboring the Dn1 resistance gene, we have identified a nucleotide-binding leucine-rich repeat (NLR) gene containing two integrated domains (IDs). These are three C-terminus ankyrin repeat domains and an N-terminus WRKY domain. The NLR core of the gene can be traced through speciation events within the grass family, with a recent WRKY domain integration that is Triticum-specific. Virus-induced gene silencing of the gene in a resistant wheat line resulted in the abolishment of the resistance response and induced a highly susceptible phenotype. Silenced plants supported a higher number of aphids, similar to the susceptible near-isogenic line (NIL), and the intrinsic rate of increase of the aphids matched that of aphids feeding on the susceptible NIL. The presence of the gene is necessary for Dn1 resistance and we have named the gene Associated with Dn resistance 1 (Adnr1) to reflect this function.

Application of Brachypodium genotypes to the analysis of type II resistance to Fusarium head blight (FHB)

The resistance to Fusarium head blight (FHB) in wheat is mainly via the restrain of fungal expansion through spike rachis (type II resistance). In order to unravel the resistance mechanisms, Brachypodium distachyon 21 (Bd21), a monocotyledonous model plant, was previously proved to interact with F. graminearum, while the disease development in spike still needs to be explored in detail. Herein, it is found that the fungal spores mainly germinate on pistil of Bd21, then the hyphae rapidly extend to the bottom of floret and enter spike rachis, similar with the infection progress in wheat. However, structural difference of spike rachis was found between Brachypodium and wheat. It was found that the spread of the fungus through the rachis node of inoculated spikelets is an important index for the evaluation of type II FHB resistance in Brachypodium under optimal conditions at 28 °C and 50%-70% humidity. To verify the feasibility of this strategy, the transcription factor TaTGA2 was overexpressed in Bd21, and transgenic plants were found to show improved resistance to F. graminearum in both spikes and detached leaves, which was further supported by the increased disease severity when silencing TaTGA2 in the wheat cultivar "Sumai 3" or in tilling "Kronos" mutants. Except for Bd21, another 49 Brachypodium germplasms were further screened for FHB resistance, and three moderately susceptible germplasms, namely, PI 317418, W6-39284, and PI 254868, feasible for transformation, were determined to be better hosts than Bd21 when evaluating heterologous genes that positively regulate FHB resistance. The present study also observed variations in the levels of FHB resistance between coleoptiles and spikes or transgenic plants and natural germplasms.

Identification and characterization of wheat stem rust resistance gene Sr21 effective against the Ug99 race group at high temperature

Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is a devastating foliar disease. The Ug99 race group has combined virulence to most stem rust (Sr) resistance genes deployed in wheat and is a threat to global wheat production. Here we identified a coiled-coil, nucleotide-binding leucine-rich repeat protein (NLR) completely linked to the Ug99 resistance gene Sr21 from Triticum monococcum. Loss-of-function mutations and transgenic complementation confirmed that this gene is Sr21. Sr21 transcripts were significantly higher at high temperatures, and this was associated with significant upregulation of pathogenesis related (PR) genes and increased levels of resistance at those temperatures. Introgression of Sr21 into hexaploid wheat resulted in lower levels of resistance than in diploid wheat, but transgenic hexaploid wheat lines with high levels of Sr21 expression showed high levels of resistance. Sr21 can be a valuable component of transgenic cassettes or gene pyramids combining multiple resistance genes against Ug99.

Wheat stem rust is a devastating disease that is threatening global wheat production. The emergence of new virulent races of this pathogen in Africa, including the Ug99 race group, has prompted global efforts to find effective resistance genes. We report here the identification of stem rust resistance gene Sr21 that is effective against the Ug99 race group. We developed a diagnostic marker to accelerate its deployment in wheat breeding programs and demonstrated that the introduction of two Sr21 copies in transgenic wheat results in high levels of resistance. An unusual characteristic of Sr21 is its increased resistance to stem rust at high temperatures. We show here that this is associated with the ability of Sr21 to coordinate the upregulation of multiple pathogenesis related genes at high temperatures. These genes slow down the growth of the pathogen and result in the characteristic Sr21 intermediate resistance reaction at high temperatures. A better understanding of this temperature dependent resistance mechanism will be useful for controlling the rust pathogens in our changing environments.

Identification and characterization of Sr13, a tetraploid wheat gene that confers resistance to the Ug99 stem rust race group

The Puccinia graminis f. sp. tritici (Pgt) Ug99 race group is virulent to most stem rust resistance genes currently deployed in wheat and poses a threat to global wheat production. The durum wheat (Triticum turgidum ssp. durum) gene Sr13 confers resistance to Ug99 and other virulent races, and is more effective at high temperatures. Using map-based cloning, we delimited a candidate region including two linked genes encoding coiled-coil nucleotide-binding leucine-rich repeat proteins designated CNL3 and CNL13. Three independent truncation mutations identified in each of these genes demonstrated that only CNL13 was required for Ug99 resistance. Transformation of an 8-kb genomic sequence including CNL13 into the susceptible wheat variety Fielder was sufficient to confer resistance to Ug99, confirming that CNL13 is Sr13CNL13 transcripts were slightly down-regulated 2-6 days after Pgt inoculation and were not affected by temperature. By contrast, six pathogenesis-related (PR) genes were up-regulated at high temperatures only when both Sr13 and Pgt were present, suggesting that they may contribute to the high temperature resistance mechanism. We identified three Sr13-resistant haplotypes, which were present in one-third of cultivated emmer and durum wheats but absent in most tested common wheats (Triticum aestivum). These results suggest that Sr13 can be used to improve Ug99 resistance in a large proportion of modern wheat cultivars. To accelerate its deployment, we developed a diagnostic marker for Sr13 The identification of Sr13 expands the number of Pgt-resistance genes that can be incorporated into multigene transgenic cassettes to control this devastating disease.

A Highly Efficient and Reproducible Fusarium spp. Inoculation Method for Brachypodium distachyon

Fusarium spp. are devastating fungal pathogens which cause significant losses in many cereal crops like wheat, maize, and barley. Genetic improvement of disease resistance requires an improved understanding of defense-associated processes operating in the host in response to an attack by Fusarium spp. Brachypodium distachyon is emerging as a model where host-cereal-infecting pathogen interactions can be studied conveniently. However, this requires developing an efficient infection assay that facilitates quick screening of germplasm (e.g., mutant lines). Here, we provide an efficient and reproducible Fusarium infection assay for Brachypodium. We believe this method will help further develop Brachypodium as a model for genetic improvement of disease resistance in cereals against Fusarium pathogens.

Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number

In crop plants, a high-density genetic linkage map is essential for both genetic and genomic researches. The complexity and the large size of wheat genome have hampered the acquisition of a high-resolution genetic map. In this study, we report a high-density genetic map based on an individual mapping population using the Affymetrix Wheat660K single-nucleotide polymorphism (SNP) array as a probe in hexaploid wheat. The resultant genetic map consisted of 119 566 loci spanning 4424.4 cM, and 119 001 of those loci were SNP markers. This genetic map showed good collinearity with the 90 K and 820 K consensus genetic maps and was also in accordance with the recently released wheat whole genome assembly. The high-density wheat genetic map will provide a major resource for future genetic and genomic research in wheat. Moreover, a comparative genomics analysis among gramineous plant genomes was conducted based on the high-density wheat genetic map, providing an overview of the structural relationships among theses gramineous plant genomes. A major stable quantitative trait locus (QTL) for kernel number per spike was characterized, providing a solid foundation for the future high-resolution mapping and map-based cloning of the targeted QTL.

Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number

In crop plants, a high-density genetic linkage map is essential for both genetic and genomic researches. The complexity and the large size of wheat genome have hampered the acquisition of a high-resolution genetic map. In this study, we report a high-density genetic map based on an individual mapping population using the Affymetrix Wheat660K single-nucleotide polymorphism (SNP) array as a probe in hexaploid wheat. The resultant genetic map consisted of 119 566 loci spanning 4424.4 cM, and 119 001 of those loci were SNP markers. This genetic map showed good collinearity with the 90 K and 820 K consensus genetic maps and was also in accordance with the recently released wheat whole genome assembly. The high-density wheat genetic map will provide a major resource for future genetic and genomic research in wheat. Moreover, a comparative genomics analysis among gramineous plant genomes was conducted based on the high-density wheat genetic map, providing an overview of the structural relationships among theses gramineous plant genomes. A major stable quantitative trait locus (QTL) for kernel number per spike was characterized, providing a solid foundation for the future high-resolution mapping and map-based cloning of the targeted QTL.

Brachypodium: A Potential Model Host for Fungal Pathogens of Turfgrasses

Brachypodium distachyon is a C3 grass that is an attractive model host system for studying pathogenicity of major turfgrass pathogens due to its genetic similarity to many cool-season turfgrasses. Infection assays with two or more isolates of the casual agents of dollar spot, brown patch, and Microdochium patch resulted in compatible interactions with B. distachyon inbred line Bd21-3. The symptoms produced by these pathogens on Bd21-3 closely resembled those observed on the natural turfgrass host (creeping bentgrass), demonstrating that B. distachyon is susceptible to the fungal pathogens that cause dollar spot, brown patch, and Microdochium patch on turfgrasses. The interaction between Sclerotinia homoeocarpa isolates and Brachypodium ecotypes was also investigated. Interestingly, differential responses of these ecotypes to S. homoeocarpa isolates was found, particularly when comparing B. distachyon to B. hybridum ecotypes. Taken together, these findings demonstrate that B. distachyon can be used as a model host system for these turfgrass diseases and leveraged for studies of molecular mechanisms contributing to host resistance.

Mapping-by-sequencing in complex polyploid genomes using genic sequence capture: a case study to map yellow rust resistance in hexaploid wheat

Previously we extended the utility of mapping-by-sequencing by combining it with sequence capture and mapping sequence data to pseudo-chromosomes that were organized using wheat-Brachypodium synteny. This, with a bespoke haplotyping algorithm, enabled us to map the flowering time locus in the diploid wheat Triticum monococcum L. identifying a set of deleted genes (Gardiner et al., 2014). Here, we develop this combination of gene enrichment and sliding window mapping-by-synteny analysis to map the Yr6 locus for yellow stripe rust resistance in hexaploid wheat. A 110 MB NimbleGen capture probe set was used to enrich and sequence a doubled haploid mapping population of hexaploid wheat derived from an Avalon and Cadenza cross. The Yr6 locus was identified by mapping to the POPSEQ chromosomal pseudomolecules using a bespoke pipeline and algorithm (Chapman et al., 2015). Furthermore the same locus was identified using newly developed pseudo-chromosome sequences as a mapping reference that are based on the genic sequence used for sequence enrichment. The pseudo-chromosomes allow us to demonstrate the application of mapping-by-sequencing to even poorly defined polyploidy genomes where chromosomes are incomplete and sub-genome assemblies are collapsed. This analysis uniquely enabled us to: compare wheat genome annotations; identify the Yr6 locus - defining a smaller genic region than was previously possible; associate the interval with one wheat sub-genome and increase the density of SNP markers associated. Finally, we built the pipeline in iPlant, making it a user-friendly community resource for phenotype mapping.

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.

Recent trends and perspectives of molecular markers against fungal diseases in wheat

Wheat accounts for 19% of the total production of major cereal crops in the world. In view of ever increasing population and demand for global food production, there is an imperative need of 40-60% increase in wheat production to meet the requirement of developing world in coming 40 years. However, both biotic and abiotic stresses are major hurdles for attaining the goal. Among the most important diseases in wheat, fungal diseases pose serious threat for widening the gap between actual and attainable yield. Fungal disease management, mainly, depends on the pathogen detection, genetic and pathological variability in population, development of resistant cultivars and deployment of effective resistant genes in different epidemiological regions. Wheat protection and breeding of resistant cultivars using conventional methods are time-consuming, intricate and slow processes. Molecular markers offer an excellent alternative in development of improved disease resistant cultivars that would lead to increase in crop yield. They are employed for tagging the important disease resistance genes and provide valuable assistance in increasing selection efficiency for valuable traits via marker assisted selection (MAS). Plant breeding strategies with known molecular markers for resistance and functional genomics enable a breeder for developing resistant cultivars of wheat against different fungal diseases.

Pushing the boundaries of resistance: insights from Brachypodium-rust interactions

The implications of global population growth urge transformation of current food and bioenergy production systems to sustainability. Members of the family Poaceae are of particular importance both in food security and for their applications as biofuel substrates. For centuries, rust fungi have threatened the production of valuable crops such as wheat, barley, oat, and other small grains; similarly, biofuel crops can also be susceptible to these pathogens. Emerging rust pathogenic races with increased virulence and recurrent rust epidemics around the world point out the vulnerability of monocultures. Basic research in plant immunity, especially in model plants, can make contributions to understanding plant resistance mechanisms and improve disease management strategies. The development of the grass Brachypodium distachyon as a genetically tractable model for monocots, especially temperate cereals and grasses, offers the possibility to overcome the experimental challenges presented by the genetic and genomic complexities of economically valuable crop plants. The numerous resources and tools available in Brachypodium have opened new doors to investigate the underlying molecular and genetic bases of plant-microbe interactions in grasses and evidence demonstrating the applicability and advantages of working with B. distachyon is increasing. Importantly, several interactions between B. distachyon and devastating plant pathogens, such rust fungi, have been examined in the context of non-host resistance. Here, we discuss the use of B. distachyon in these various pathosystems. Exploiting B. distachyon to understand the mechanisms underpinning disease resistance to non-adapted rust fungi may provide effective and durable approaches to fend off these pathogens. The close phylogenetic relationship among Brachypodium spp. and grasses with industrial and agronomic value support harnessing this model plant to improve cropping systems and encourage its use in translational research.

Fine mapping and identification of a candidate gene for the barley Un8 true loose smut resistance gene

The candidate gene for the barley Un8 true loose smut resistance gene encodes a deduced protein containing two tandem protein kinase domains. In North America, durable resistance against all known isolates of barley true loose smut, caused by the basidiomycete pathogen Ustilago nuda (Jens.) Rostr. (U. nuda), is under the control of the Un8 resistance gene. Previous genetic studies mapped Un8 to the long arm of chromosome 5 (1HL). Here, a population of 4625 lines segregating for Un8 was used to delimit the Un8 gene to a 0.108 cM interval on chromosome arm 1HL, and assign it to fingerprinted contig 546 of the barley physical map. The minimal tilling path was identified for the Un8 locus using two flanking markers and consisted of two overlapping bacterial artificial chromosomes. One gene located close to a marker co-segregating with Un8 showed high sequence identity to a disease resistance gene containing two kinase domains. Sequence of the candidate gene from the parents of the segregating population, and in an additional 19 barley lines representing a broader spectrum of diversity, showed there was no intron in alleles present in either resistant or susceptible lines, and fifteen amino acid variations unique to the deduced protein sequence in resistant lines differentiated it from the deduced protein sequences in susceptible lines. Some of these variations were present within putative functional domains which may cause a loss of function in the deduced protein sequences within susceptible lines.

Exploring the utility of Brachypodium distachyon as a model pathosystem for the wheat pathogen Zymoseptoria tritici

BACKGROUND

Zymoseptoria tritici, the causative organism of Septoria tritici blotch disease is a prevalent biotic stressor of wheat production, exerting substantial economic constraints on farmers, requiring intensive chemical control to protect yields. A hemibiotrophic pathogen with a long asymptomless phase of up to 11 days post inoculation (dpi) before a rapid switch to necrotrophy; a deficit exists in our understanding of the events occurring within the host during the two phases of infection. Brachypodium distachyon has demonstrated its potential as a model species for the investigation of fungal disease resistance in cereal and grass species. The aim of this study was to assess the physical interaction between Z. tritici (strain IPO323) and B. distachyon and examine its potential as a model pathosystem for Z. tritici.

RESULTS

Septoria tritici blotch symptoms developed on the wheat cultivar Riband from 12 dpi with pycnidial formation abundant by 20 dpi. Symptoms on B. distachyon ecotype Bd21-1 were visible from 1 dpi: characteristic pale, water soaked lesions which developed into blotch-like lesions by 4 dpi. These lesions then became necrotic with chlorotic regions expanding up to 7 dpi. Sporulation on B. distachyon tissues was not observed and no evidence of fungal penetration could be obtained, indicating that Z. tritici was unable to complete its life cycle within B. distachyon ecotypes. However, observation of host responses to the Z. tritici strain IPO323 in five B. distachyon ecotypes revealed a variation in resistance responses, ranging from immunity to a chlorotic/necrotic phenotype.

CONCLUSIONS

The observed interactions suggest that B. distachyon is an incompatible host for Z. tritici infection, with STB symptom development on B. distachyon comparable to that observed during the early infection stages on the natural host, wheat. However first visible symptoms occurred more rapidly on B. distachyon; from 1 dpi in comparison to 12 dpi in wheat. Consequently, we propose that the interaction between B. distachyon and Z. tritici as observed in this study could serve as a suitable model pathosystem with which to investigate mechanisms underpinning an incompatible host response to Z. tritici.

An assessment of heavy ion irradiation mutagenesis for reverse genetics in wheat (Triticum aestivum L.)

Reverse genetic techniques harnessing mutational approaches are powerful tools that can provide substantial insight into gene function in plants. However, as compared to diploid species, reverse genetic analyses in polyploid plants such as bread wheat can present substantial challenges associated with high levels of sequence and functional similarity amongst homoeologous loci. We previously developed a high-throughput method to identify deletions of genes within a physically mutagenized wheat population. Here we describe our efforts to combine multiple homoeologous deletions of three candidate disease susceptibility genes (TaWRKY11, TaPFT1 and TaPLDß1). We were able to produce lines featuring homozygous deletions at two of the three homoeoloci for all genes, but this was dependent on the individual mutants used in crossing. Intriguingly, despite extensive efforts, viable lines possessing homozygous deletions at all three homoeoloci could not be produced for any of the candidate genes. To investigate deletion size as a possible reason for this phenomenon, we developed an amplicon sequencing approach based on synteny to Brachypodium distachyon to assess the size of the deletions removing one candidate gene (TaPFT1) in our mutants. These analyses revealed that genomic deletions removing the locus are relatively large, resulting in the loss of multiple additional genes. The implications of this work for the use of heavy ion mutagenesis for reverse genetic analyses in wheat are discussed.

Genetic and physical mapping of powdery mildew resistance gene MlHLT in Chinese wheat landrace Hulutou

A powdery mildew resistance gene MlHLT derived from a Chinese wheat landrace maps within a 3.6 centimorgan (cM) genetic interval spanning a 13.4 megabase (Mb) physical genomic region on chromosome 1DS. Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt) is a devastating disease that can cause severe yield losses. Chinese wheat landrace Hulutou confers nearly immune resistance against prevailing Bgt isolate E09 in Beijing. Genetic analysis indicate that the powdery mildew resistance of Hulutou is controlled by a single dominant gene, provisionally designated MlHLT. Bulked segregant analysis(BSA) and simple sequence repeat (SSR) mapping showed that MlHLT is located on chromosome arm 1DS between markers Xgwm337 and Xcfd83/Xcfd72. By applying comparative genomics analysis, collinearity genomic regions of the MlHLT locus on Aegilops tauschii chromosome 1DS were identified in Brachypodium distachyon chromosome 2, rice chromosome 5 and sorghum chromosome 9, respectively. Three new polymorphic markers were developed using the draft genome sequences and the extended single nucleotide polymorphism (SNP) marker sequences of Ae. tauschii accession AL8/78, as well as the Triticum aestivum cv. Chinese Spring 454 contig sequences and the International Wheat Genome Sequencing Consortium (IWGSC) survey sequences. MlHLT mapped into a 3.6 cM genetic interval spanning 13.4 Mb physical genomic region containing seven contigs (ctg220, ctg4623, ctg1063, ctg5929, ctg3163, ctg699 and ctg1065) on 1DS that has synteny with a 369.8 kb genomic region in Brachypodium, a 380.8 kb genomic region in rice and a 298.4 kb genomic region in sorghum. The genetic and physical maps of MlHLT provide framework for map-based cloning and marker-assisted selection (MAS) of the powdery mildew resistance gene MlHLT in Hulutou.

Genome-wide association for grain yield under rainfed conditions in historical wheat cultivars from Pakistan

Genome-wide association studies (GWAS) were undertaken to identify SNP markers associated with yield and yield-related traits in 123 Pakistani historical wheat cultivars evaluated during 2011-2014 seasons under rainfed field conditions. The population was genotyped by using high-density Illumina iSelect 90K single nucleotide polymorphism (SNP) assay, and finally 14,960 high quality SNPs were used in GWAS. Population structure examined using 1000 unlinked markers identified seven subpopulations (K = 7) that were representative of different breeding programs in Pakistan, in addition to local landraces. Forty four stable marker-trait associations (MTAs) with -log p > 4 were identified for nine yield-related traits. Nine multi-trait MTAs were found on chromosomes 1AL, 1BS, 2AL, 2BS, 2BL, 4BL, 5BL, 6AL, and 6BL, and those on 5BL and 6AL were stable across two seasons. Gene annotation and syntey identified that 14 trait-associated SNPs were linked to genes having significant importance in plant development. Favorable alleles for days to heading (DH), plant height (PH), thousand grain weight (TGW), and grain yield (GY) showed minor additive effects and their frequencies were slightly higher in cultivars released after 2000. However, no selection pressure on any favorable allele was identified. These genomic regions identified have historically contributed to achieve yield gains from 2.63 million tons in 1947 to 25.7 million tons in 2015. Future breeding strategies can be devised to initiate marker assisted breeding to accumulate these favorable alleles of SNPs associated with yield-related traits to increase grain yield. Additionally, in silico identification of 454-contigs corresponding to MTAs will facilitate fine mapping and subsequent cloning of candidate genes and functional marker development.

Nonhost resistance to rust pathogens - a continuation of continua

The rust fungi (order: Pucciniales) are a group of widely distributed fungal plant pathogens, which can infect representatives of all vascular plant groups. Rust diseases significantly impact several crop species and considerable research focuses on understanding the basis of host specificity and nonhost resistance. Like many pathogens, rust fungi vary considerably in the number of hosts they can infect, such as wheat leaf rust (Puccinia triticina), which can only infect species in the genera Triticum and Aegilops, whereas Asian soybean rust (Phakopsora pachyrhizi) is known to infect over 95 species from over 42 genera. A greater understanding of the genetic basis determining host range has the potential to identify sources of durable resistance for agronomically important crops. Delimiting the boundary between host and nonhost has been complicated by the quantitative nature of phenotypes in the transition between these two states. Plant-pathogen interactions in this intermediate state are characterized either by (1) the majority of accessions of a species being resistant to the rust or (2) the rust only being able to partially complete key components of its life cycle. This leads to a continuum of disease phenotypes in the interaction with different plant species, observed as a range from compatibility (host) to complete immunity within a species (nonhost). In this review we will highlight how the quantitative nature of disease resistance in these intermediate interactions is caused by a continuum of defense barriers, which a pathogen needs to overcome for successfully establishing itself in the host. To illustrate continua as this underlying principle, we will discuss the advances that have been made in studying nonhost resistance towards rust pathogens, particularly cereal rust pathogens.

Using genic sequence capture in combination with a syntenic pseudo genome to map a deletion mutant in a wheat species

Mapping-by-sequencing analyses have largely required a complete reference sequence and employed whole genome re-sequencing. In species such as wheat, no finished genome reference sequence is available. Additionally, because of its large genome size (17 Gb), re-sequencing at sufficient depth of coverage is not practical. Here, we extend the utility of mapping by sequencing, developing a bespoke pipeline and algorithm to map an early-flowering locus in einkorn wheat (Triticum monococcum L.) that is closely related to the bread wheat genome A progenitor. We have developed a genomic enrichment approach using the gene-rich regions of hexaploid bread wheat to design a 110-Mbp NimbleGen SeqCap EZ in solution capture probe set, representing the majority of genes in wheat. Here, we use the capture probe set to enrich and sequence an F2 mapping population of the mutant. The mutant locus was identified in T. monococcum, which lacks a complete genome reference sequence, by mapping the enriched data set onto pseudo-chromosomes derived from the capture probe target sequence, with a long-range order of genes based on synteny of wheat with Brachypodium distachyon. Using this approach we are able to map the region and identify a set of deleted genes within the interval.

Identification of candidate genes, regions and markers for pre-harvest sprouting resistance in wheat (Triticum aestivum L.)

BACKGROUND

Pre-harvest sprouting (PHS) of wheat grain leads to a reduction in grain yield and quality. The availability of markers for marker-assisted selection (MAS) of PHS resistance will serve to enhance breeding selection and advancement of lines for cultivar development. The aim of this study was to identify candidate regions and develop molecular markers for PHS resistance in wheat. This was achieved via high density mapping of single nucleotide polymorphism (SNP) markers from an Illumina 90 K Infinium Custom Beadchip in a doubled haploid (DH) population derived from a RL4452/'AC Domain' cross and subsequent detection of quantitative trait loci (QTL) for PHS related traits (falling number [FN], germination index [GI] and sprouting index [SI]). SNP marker sequences flanking QTL were used to locate colinear regions in Brachypodium and rice, and identify genic markers associated with PHS resistance that can be utilized for MAS in wheat.

RESULTS

A linkage map spanning 2569.4 cM was constructed with a total of 12,201 SNP, simple sequence repeat (SSR), diversity arrays technology (DArT) and expressed sequence tag (EST) markers. QTL analyses using Multiple Interval Mapping (MIM) identified four QTL for PHS resistance traits on chromosomes 3B, 4A, 7B and 7D. Sequences of SNPs flanking these QTL were subject to a BLASTN search on the International Wheat Genome Sequencing Consortium (IWGSC) database (http://wheat-urgi.versailles.inra.fr/Seq-Repository). Best survey sequence hits were subject to a BLASTN search on Gramene (www.gramene.org) against both Brachypodium and rice databases, and candidate genes and regions for PHS resistance were identified. A total of 18 SNP flanking sequences on chromosomes 3B, 4A, 7B and 7D were converted to KASP markers and validated with matching genotype calls of Infinium SNP data.

CONCLUSIONS

Our study identified candidate genes involved in abscissic acid (ABA) and gibberellin (GA) metabolism, and flowering time in four genomic regions of Brachypodium and rice respectively, in addition to 18 KASP markers for PHS resistance in wheat. These markers can be deployed in future genetic studies of PHS resistance and might also be useful in the evaluation of PHS in germplasm and breeding material.

Dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium distachyon: comparison with common wheat and Aegilops peregrina

BACKGROUND

Thorough understanding of seed starch biosynthesis and accumulation mechanisms is of great importance for agriculture and crop improvement strategies. We conducted the first comprehensive study of the dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium distachyon and compared the findings with those reported for common wheat (Chinese Spring, CS) and Aegilops peregrina.

RESULTS

Only B-granules were identified in Brachypodium Bd21, and the shape variation and development of starch granules were similar in the B-granules of CS and Bd21. Phylogenetic analysis showed that most of the Bd21 starch synthesis-related genes were more similar to those in wheat than in rice. Early expression of key genes in Bd21 starch biosynthesis mediate starch synthesis in the pericarp; intermediate-stage expression increases the number and size of starch granules. In contrast, these enzymes in CS and Ae. peregrina were mostly expressed at intermediate stages, driving production of new B-granules and increasing the granule size, respectively. Immunogold labeling showed that granule-bound starch synthase (GBSSI; related to amylose synthesis) was mainly present in starch granules: at lower levels in the B-granules of Bd21 than in CS. Furthermore, GBSSI was phosphorylated at threonine 183 and tyrosine 185 in the starch synthase catalytic domain in CS and Ae. peregrina, but neither site was phosphorylated in Bd21, suggesting GBSSI phosphorylation could improve amylose biosynthesis.

CONCLUSIONS

Bd21 contains only B-granules, and the expression of key genes in the three studied genera is consistent with the dynamic development of starch granules. GBSSI is present in greater amounts in the B-granules of CS than in Bd21; two phosphorylation sites (Thr183 and Tyr185) were found in Triticum and Aegilops; these sites were not phosphorylated in Bd21. GBSSI phosphorylation may reflect its importance in amylose synthesis.

Fine physical and genetic mapping of powdery mildew resistance gene MlIW172 originating from wild emmer (Triticum dicoccoides)

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is one of the most important wheat diseases in the world. In this study, a single dominant powdery mildew resistance gene MlIW172 was identified in the IW172 wild emmer accession and mapped to the distal region of chromosome arm 7AL (bin7AL-16-0.86-0.90) via molecular marker analysis. MlIW172 was closely linked with the RFLP probe Xpsr680-derived STS marker Xmag2185 and the EST markers BE405531 and BE637476. This suggested that MlIW172 might be allelic to the Pm1 locus or a new locus closely linked to Pm1. By screening genomic BAC library of durum wheat cv. Langdon and 7AL-specific BAC library of hexaploid wheat cv. Chinese Spring, and after analyzing genome scaffolds of Triticum urartu containing the marker sequences, additional markers were developed to construct a fine genetic linkage map on the MlIW172 locus region and to delineate the resistance gene within a 0.48 cM interval. Comparative genetics analyses using ESTs and RFLP probe sequences flanking the MlIW172 region against other grass species revealed a general co-linearity in this region with the orthologous genomic regions of rice chromosome 6, Brachypodium chromosome 1, and sorghum chromosome 10. However, orthologous resistance gene-like RGA sequences were only present in wheat and Brachypodium. The BAC contigs and sequence scaffolds that we have developed provide a framework for the physical mapping and map-based cloning of MlIW172.

Identification of a robust molecular marker for the detection of the stem rust resistance gene Sr45 in common wheat

Fine mapping of the Ug99 effective stem rust resistance gene Sr45 introgressed into common wheat from the D -genome goatgrass Aegilops tauschii. Stem rust resistance gene Sr45, discovered in Aegilops tauschii, the progenitor of the D -genome of wheat, is effective against commercially important Puccinia graminis f. sp. tritici races prevalent in Australia, South Africa and the Ug99 race group. A synthetic hexaploid wheat (RL5406) generated by crossing Ae. tauschii accession RL5289 (carrying Sr45 and the leaf rust resistance gene Lr21) with a tetraploid experimental line 'TetraCanthatch' was previously used as the source in the transfer of these rust resistance genes to other hexaploid cultivars. Previous genetic studies on hexaploid wheats mapped Sr45 on the short arm of chromosome 1D with the following gene order: centromere-Sr45-Sr33-Lr21-telomere. To identify closely linked markers, we fine mapped the Sr45 region in a large mapping population generated by crossing CS1D5406 (disomic substitution line with chromosome 1D of RL5406 substituted for Chinese Spring 1D) with Chinese Spring. Closely linked markers based on 1DS-specific microsatellites, expressed sequence tags and AFLP were useful in the delineation of the Sr45 region. Sequences from an AFLP marker amplified a fragment that was linked with Sr45 at a distance of 0.39 cM. The fragment was located in a bacterial artificial chromosome clone of contig (ctg)2981 of the Ae. tauschii accession AL8/78 physical map. A PCR marker derived from clone MI221O11 of ctg2981 amplified 1DS-specific sequence that harboured an 18-bp indel polymorphism that specifically tagged the Sr45 carrying haplotype. This new Sr45 marker can be combined with a previously reported marker for Lr21, which will facilitate selecting Sr45 and Lr21 in breeding populations.

Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe

Using Roche/454 technology, we sequenced the chloroplast genomes of 12 Triticeae species, including bread wheat, barley and rye, as well as the diploid progenitors and relatives of bread wheat Triticum urartu, Aegilops speltoides and Ae. tauschii. Two wild tetraploid taxa, Ae. cylindrica and Ae. geniculata, were also included. Additionally, we incorporated wild Einkorn wheat Triticum boeoticum and its domesticated form T. monococcum and two Hordeum spontaneum (wild barley) genotypes. Chloroplast genomes were used for overall sequence comparison, phylogenetic analysis and dating of divergence times. We estimate that barley diverged from rye and wheat approximately 8–9 million years ago (MYA). The genome donors of hexaploid wheat diverged between 2.1–2.9 MYA, while rye diverged from Triticum aestivum approximately 3–4 MYA, more recently than previously estimated. Interestingly, the A genome taxa T. boeoticum and T. urartu were estimated to have diverged approximately 570,000 years ago. As these two have a reproductive barrier, the divergence time estimate also provides an upper limit for the time required for the formation of a species boundary between the two. Furthermore, we conclusively show that the chloroplast genome of hexaploid wheat was contributed by the B genome donor and that this unknown species diverged from Ae. speltoides about 980,000 years ago. Additionally, sequence alignments identified a translocation of a chloroplast segment to the nuclear genome which is specific to the rye/wheat lineage. We propose the presented phylogeny and divergence time estimates as a reference framework for future studies on Triticeae.

TaWRKY68 responses to biotic stresses are revealed by the orthologous genes from major cereals

WRKY transcription factors have been extensively characterized in the past 20 years, but in wheat, studies on WRKY genes and their function are lagging behind many other species. To explore the function of wheat WRKY genes, we identified a TaWRKY68 gene from a common wheat cultivar. It encodes a protein comprising 313 amino acids which harbors 19 conserved motifs or active sites. Gene expression patterns were determined by analyzing microarray data of TaWRKY68 in wheat and of orthologous genes from maize, rice and barley using Genevestigator. TaWRKY68 orthologs were identified and clustered using DELTA-BLAST and COBALT programs available at NCBI. The results showed that these genes, which are expressed in all tissues tested, had relatively higher levels in the roots and were up-regulated in response to biotic stresses. Bioinformatics results were confirmed by RT-PCR experiments using wheat plants infected by Agrobacterium tumefaciens and Blumeria graminis, or treated with Deoxynivalenol, a Fusarium graminearum-induced mycotoxin in wheat or barley. In summary, TaWRKY68 functions differ during plant developmental stages and might be representing a hub gene function in wheat responses to various biotic stresses. It was also found that including data from major cereal genes in the bioinformatics analysis gave more accurate and comprehensive predictions of wheat gene functions.

Comparative high-resolution mapping of the wax inhibitors Iw1 and Iw2 in hexaploid wheat

The wax (glaucousness) on wheat leaves and stems is mainly controlled by two sets of genes: glaucousness loci (W1 and W2) and non-glaucousness loci (Iw1 and Iw2). The non-glaucousness (Iw) loci act as inhibitors of the glaucousness loci (W). High-resolution comparative genetic linkage maps of the wax inhibitors Iw1 originating from Triticum dicoccoides, and Iw2 from Aegilops tauschii were developed by comparative genomics analyses of Brachypodium, sorghum and rice genomic sequences corresponding to the syntenic regions of the Iw loci in wheat. Eleven Iw1 and eight Iw2 linked EST markers were developed and mapped to linkage maps on the distal regions of chromosomes 2BS and 2DS, respectively. The Iw1 locus mapped within a 0.96 cM interval flanked by the BE498358 and CA499581 EST markers that are collinear with 122 kb, 202 kb, and 466 kb genomic regions in the Brachypodium 5S chromosome, the sorghum 6S chromosome and the rice 4S chromosome, respectively. The Iw2 locus was located in a 4.1 to 5.4-cM interval in chromosome 2DS that is flanked by the CJ886319 and CJ519831 EST markers, and this region is collinear with a 2.3 cM region spanning the Iw1 locus on chromosome 2BS. Both Iw1 and Iw2 co-segregated with the BF474014 and CJ876545 EST markers, indicating they are most likely orthologs on 2BS and 2DS. These high-resolution maps can serve as a framework for chromosome landing, physical mapping and map-based cloning of the wax inhibitors in wheat.

The physical map of wheat chromosome 1BS provides insights into its gene space organization and evolution

Background

The wheat genome sequence is an essential tool for advanced genomic research and improvements. The generation of a high-quality wheat genome sequence is challenging due to its complex 17 Gb polyploid genome. To overcome these difficulties, sequencing through the construction of BAC-based physical maps of individual chromosomes is employed by the wheat genomics community. Here, we present the construction of the first comprehensive physical map of chromosome 1BS, and illustrate its unique gene space organization and evolution.

Results

Fingerprinted BAC clones were assembled into 57 long scaffolds, anchored and ordered with 2,438 markers, covering 83% of chromosome 1BS. The BAC-based chromosome 1BS physical map and gene order of the orthologous regions of model grass species were consistent, providing strong support for the reliability of the chromosome 1BS assembly. The gene space for chromosome 1BS spans the entire length of the chromosome arm, with 76% of the genes organized in small gene islands, accompanied by a two-fold increase in gene density from the centromere to the telomere.

Conclusions

This study provides new evidence on common and chromosome-specific features in the organization and evolution of the wheat genome, including a non-uniform distribution of gene density along the centromere-telomere axis, abundance of non-syntenic genes, the degree of colinearity with other grass genomes and a non-uniform size expansion along the centromere-telomere axis compared with other model cereal genomes. The high-quality physical map constructed in this study provides a solid basis for the assembly of a reference sequence of chromosome 1BS and for breeding applications.

Molecular characterisation and evolution of HMW glutenin subunit genes in Brachypodium distachyon L

Brachypodium distachyon, a small wild grass within the Pooideae family, is a new model organism for exploring the functional genomics of cereal crops. It was shown to have close relationships to wheat, barley and rice. Here, we describe the molecular characterisation and evolutionary relationships of high molecular weight glutenin subunits (HMW-GS) genes from B. distachyon. Sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE), high performance capillary electrophoresis (HPCE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses demonstrated that there was no HMW-GS expression in the Brachypodium grains due to the silencing of their encoding genes. Through allele-specific polymerase chain reaction (AS-PCR) amplification and cloning, a total of 13 HMW-GS encoding genes from diploid, tetraploid and hexaploid Brachypodium species were obtained, and all of them had typical structural features of y-type HMW-GS genes from common wheat and related species, particularly more similar to the 1Dy12 gene. However, the presence of an in-frame premature stop codon (TAG) at position 1521 in the coding region resulted in the conversion of all the genes to pseudogenes. Further, quantitative real-time PCR (qRT-PCR) analysis revealed that HMW-GS genes in B. distachyon displayed a similar trend, but with a low transcriptional expression profile during grain development due to the occurrence of the stop codon. Phylogenetic analysis showed that the highly conserved Glu-1-2 loci were presented in B. distachyon, which displayed close phylogenetic evolutionary relationships with Triticum and related species.

Description of durum wheat linkage map and comparative sequence analysis of wheat mapped DArT markers with rice and Brachypodium genomes

Background

The importance of wheat to the world economy, together with progresses in high-throughput next-generation DNA sequencing, have accelerated initiatives of genetic research for wheat improvement. The availability of high density linkage maps is crucial to identify genotype-phenotype associations, but also for anchoring BAC contigs to genetic maps, a strategy followed for sequencing the wheat genome.

Results

Here we report a genetic linkage map in a durum wheat segregating population and the study of mapped DArT markers. The linkage map consists of 126 gSSR, 31 EST-SSR and 351 DArT markers distributed in 24 linkage groups for a total length of 1,272 cM. Through bioinformatic approaches we have analysed 327 DArT clones to reveal their redundancy, syntenic and functional aspects. The DNA sequences of 174 DArT markers were assembled into a non-redundant set of 60 marker clusters. This explained the generation of clusters in very small chromosome regions across genomes. Of these DArT markers, 61 showed highly significant (Expectation < E-10) BLAST similarity to gene sequences in public databases of model species such as Brachypodium and rice. Based on sequence alignments, the analysis revealed a mosaic gene conservation, with 54 and 72 genes present in rice and Brachypodium species, respectively.

Conclusions

In the present manuscript we provide a detailed DArT markers characterization and the basis for future efforts in durum wheat map comparing.

Next generation characterisation of cereal genomes for marker discovery

Cereal crops form the bulk of the world’s food sources, and thus their importance cannot be understated. Crop breeding programs increasingly rely on high-resolution molecular genetic markers to accelerate the breeding process. The development of these markers is hampered by the complexity of some of the major cereal crop genomes, as well as the time and cost required. In this review, we address current and future methods available for the characterisation of cereal genomes, with an emphasis on faster and more cost effective approaches for genome sequencing and the development of markers for trait association and marker assisted selection (MAS) in crop breeding programs.

Dynamic evolution of rht-1 homologous regions in grass genomes

Hexaploid bread wheat contains A, B, and D three subgenomes with its well-characterized ancestral genomes existed at diploid and tetraploid levels, making the wheat act as a good model species for studying evolutionary genomic dynamics. Here, we performed intra- and inter-species comparative analyses of wheat and related grass genomes to examine the dynamics of homologous regions surrounding Rht-1, a well-known "green revolution" gene. Our results showed that the divergence of the two A genomes in the Rht-1 region from the diploid and tetraploid species is greater than that from the tetraploid and hexaploid wheat. The divergence of D genome between diploid and hexaploid is lower than those of A genome, suggesting that D genome diverged latter than others. The divergence among the A, B and D subgenomes was larger than that among different ploidy levels for each subgenome which mainly resulted from genomic structural variation of insertions and, perhaps deletions, of the repetitive sequences. Meanwhile, the repetitive sequences caused genome expansion further after the divergence of the three subgenomes. However, several conserved non-coding sequences were identified to be shared among the three subgenomes of wheat, suggesting that they may have played an important role to maintain the homolog of three subgenomes. This is a pilot study on evolutionary dynamics across the wheat ploids, subgenomes and differently related grasses. Our results gained new insights into evolutionary dynamics of Rht-1 region at sequence level as well as the evolution of wheat during the plolyploidization process.

Brachypodium sylvaticum, a model for perennial grasses: transformation and inbred line development

Perennial species offer significant advantages as crops including reduced soil erosion, lower energy inputs after the first year, deeper root systems that access more soil moisture, and decreased fertilizer inputs due to the remobilization of nutrients at the end of the growing season. These advantages are particularly relevant for emerging biomass crops and it is projected that perennial grasses will be among the most important dedicated biomass crops. The advantages offered by perennial crops could also prove favorable for incorporation into annual grain crops like wheat, rice, sorghum and barley, especially under the dryer and more variable climate conditions projected for many grain-producing regions. Thus, it would be useful to have a perennial model system to test biotechnological approaches to crop improvement and for fundamental research. The perennial grass Brachypodiumsylvaticum is a candidate for such a model because it is diploid, has a small genome, is self-fertile, has a modest stature, and short generation time. Its close relationship to the annual model Brachypodiumdistachyon will facilitate comparative studies and allow researchers to leverage the resources developed for B. distachyon. Here we report on the development of two keystone resources that are essential for a model plant: high-efficiency transformation and inbred lines. Using Agrobacterium tumefaciens-mediated transformation we achieved an average transformation efficiency of 67%. We also surveyed the genetic diversity of 19 accessions from the National Plant Germplasm System using SSR markers and created 15 inbred lines.

Infection of Brachypodium distachyon with selected grass rust pathogens

The model temperate grass Brachypodium distachyon is considered a nonhost for wheat rust diseases caused by Puccinia graminis f. sp. tritici, P. triticina, and P. striiformis. Up to 140 Brachypodium accessions were infected with these three rust species, in addition to P. graminis ff. spp. avena and phalaridis. Related B. distachyon lines showed similar cytological nonhost resistance (NHR) phenotypes, and an inverse relationship between P. graminis f. sp. tritici and P. striiformis growth was observed in many lines, with accessions that allowed the most growth of P. graminis f. sp. tritici showing the least P. striiformis development and vice versa. Callose deposition patterns during infection by all three rust species showed similarity to the wheat basal defense response while cell death that resulted in autofluorescence did not appear to be a major component of the defense response. Infection of B. distachyon with P. graminis f. sp. avena and P. graminis f. sp. phalaridis produced much greater colonization, indicating that P. graminis rusts with Poeae hosts show greater ability to infect B. distachyon than those with Triticeae hosts. P. striiformis infection of progeny from two B. distachyon families demonstrated that these NHR phenotypes are highly heritable and appear to be under relatively simple genetic control, making this species a powerful tool for elucidating the molecular basis of NHR to cereal rust pathogens.

Transcriptome analysis in sheepgrass (Leymus chinensis): a dominant perennial grass of the Eurasian Steppe

BACKGROUND

Sheepgrass [Leymus chinensis (Trin.) Tzvel.] is an important perennial forage grass across the Eurasian Steppe and is known for its adaptability to various environmental conditions. However, insufficient data resources in public databases for sheepgrass limited our understanding of the mechanism of environmental adaptations, gene discovery and molecular marker development.

RESULTS

The transcriptome of sheepgrass was sequenced using Roche 454 pyrosequencing technology. We assembled 952,328 high-quality reads into 87,214 unigenes, including 32,416 contigs and 54,798 singletons. There were 15,450 contigs over 500 bp in length. BLAST searches of our database against Swiss-Prot and NCBI non-redundant protein sequences (nr) databases resulted in the annotation of 54,584 (62.6%) of the unigenes. Gene Ontology (GO) analysis assigned 89,129 GO term annotations for 17,463 unigenes. We identified 11,675 core Poaceae-specific and 12,811 putative sheepgrass-specific unigenes by BLAST searches against all plant genome and transcriptome databases. A total of 2,979 specific freezing-responsive unigenes were found from this RNAseq dataset. We identified 3,818 EST-SSRs in 3,597 unigenes, and some SSRs contained unigenes that were also candidates for freezing-response genes. Characterizations of nucleotide repeats and dominant motifs of SSRs in sheepgrass were also performed. Similarity and phylogenetic analysis indicated that sheepgrass is closely related to barley and wheat.

CONCLUSIONS

This research has greatly enriched sheepgrass transcriptome resources. The identified stress-related genes will help us to decipher the genetic basis of the environmental and ecological adaptations of this species and will be used to improve wheat and barley crops through hybridization or genetic transformation. The EST-SSRs reported here will be a valuable resource for future gene-phenotype studies and for the molecular breeding of sheepgrass and other Poaceae species.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

An atlas of over 90,000 conserved noncoding sequences provides insight into crucifer regulatory regions

Despite the central importance of noncoding DNA to gene regulation and evolution, understanding of the extent of selection on plant noncoding DNA remains limited compared to that of other organisms. Here we report sequencing of genomes from three Brassicaceae species (Leavenworthia alabamica, Sisymbrium irio and Aethionema arabicum) and their joint analysis with six previously sequenced crucifer genomes. Conservation across orthologous bases suggests that at least 17% of the Arabidopsis thaliana genome is under selection, with nearly one-quarter of the sequence under selection lying outside of coding regions. Much of this sequence can be localized to approximately 90,000 conserved noncoding sequences (CNSs) that show evidence of transcriptional and post-transcriptional regulation. Population genomics analyses of two crucifer species, A. thaliana and Capsella grandiflora, confirm that most of the identified CNSs are evolving under medium to strong purifying selection. Overall, these CNSs highlight both similarities and several key differences between the regulatory DNA of plants and other species.

Towards positional isolation of three quantitative trait loci conferring resistance to powdery mildew in two Spanish barley landraces

Three quantitative trait loci (QTL) conferring broad spectrum resistance to powdery mildew, caused by the fungus Blumeria graminis f. sp. hordei, were previously identified on chromosomes 7HS, 7HL and 6HL in the Spanish barley landrace-derived lines SBCC097 and SBCC145. In the present work, a genome-wide putative linear gene index of barley (Genome Zipper) and the first draft of the physical, genetic and functional sequence of the barley genome were used to go one step further in the shortening and explicit demarcation on the barley genome of these regions conferring resistance to powdery mildew as well as in the identification of candidate genes. First, a comparative analysis of the target regions to the barley Genome Zippers of chromosomes 7H and 6H allowed the development of 25 new gene-based molecular markers, which slightly better delimit the QTL intervals. These new markers provided the framework for anchoring of genetic and physical maps, figuring out the outline of the barley genome at the target regions in SBCC097 and SBCC145. The outermost flanking markers of QTLs on 7HS, 7HL and 6HL defined a physical area of 4 Mb, 3.7 Mb and 3.2 Mb, respectively. In total, 21, 10 and 16 genes on 7HS, 7HL and 6HL, respectively, could be interpreted as potential candidates to explain the resistance to powdery mildew, as they encode proteins of related functions with respect to the known pathogen defense-related processes. The majority of these were annotated as belonging to the NBS-LRR class or protein kinase family.