Wheat genomics comes of age

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

A conditional mutation in a wheat (Triticum aestivum L.) gene regulating root morphology

KEY MESSAGE

Characterisation and genetic mapping of a key gene defining root morphology in bread wheat. Root morphology is central to plants for the efficient uptake up of soil water and mineral nutrients. Here we describe a conditional mutant of hexaploid wheat (Triticum aestivum L.) that when grown in soil with high Ca2+ develops a larger rhizosheath accompanied with shorter roots than the wild type. In wheat, rhizosheath size is a reliable surrogate for root hair length and this was verified in the mutant which possessed longer root hairs than the wild type when grown in high Ca2+ soil. We named the mutant Stumpy and showed it to be due to a single semi-dominant mutation. The short root phenotype at high Ca2+ was due to reduced cellular elongation which might also explain the long root hair phenotype. Analysis of root cell walls showed that the polysaccharide composition of Stumpy roots is remodelled when grown at non-permissive (high) Ca2+ concentrations. The mutation mapped to chromosome 7B and sequencing of the 7B chromosomes in both wild type and Stumpy identified a candidate gene underlying the Stumpy mutation. As part of the process to determine whether the candidate gene was causative, we identified wheat lines in a Cadenza TILLING population with large rhizosheaths but accompanied with normal root length. This finding illustrates the potential of manipulating the gene to disconnect root length from root hair length as a means of developing wheat lines with improved efficiency of nutrient and water uptake. The Stumpy mutant will be valuable for understanding the mechanisms that regulate root morphology in wheat.

Exome-based new allele-specific PCR markers and transferability for sodicity tolerance in bread wheat (Triticum aestivum L.)

Targeted exome-based genotype by sequencing (t-GBS), a sequencing technology that tags SNPs and haplotypes in gene-rich regions was used in previous genome-wide association studies (GWAS) for sodicity tolerance in bread wheat. Thirty-nine novel SNPs including 18 haplotypes for yield and yield-components were identified. The present study aimed at developing SNP-derived markers by precisely locating new SNPs on ~180 bp allelic sequence of t-GBS, marker validation, and SNP functional characterization based on its exonic location. We identified unknown locations of significant SNPs/haplotypes by aligning allelic sequences on to IWGSC RefSeqv1.0 on respective chromosomes. Eighteen out of the target 39 SNP locations fulfilled the criteria for producing PCR markers, among which only eight produced polymorphic signals. These eight markers associated with yield, plants m-2, heads m-2, and harvest index, including a pleiotropic marker for yield, harvest index, and grains/head were validated for its amplification efficiency and phenotypic effects in focused identification germplasm strategy (FIGS) wheat set and a doubled haploid (DH) population (Scepter/IG107116). The phenotypic variation explained by these markers are in the range of 4.1-37.6 in the FIGS population. High throughput PCR-based genotyping using new markers and association with phenotypes in FIGS wheat set and DH population validated the effect of functional SNP on closely associated genes-calcineurin B-like- and dirigent protein, basic helix-loop-helix (BHLH-), plant homeodomain (PHD-) and helix-turn-helix myeloblastosis (HTH myb) type -transcription factor. Further, genome-wide SNP annotation using SnpEff tool confirmed that these SNPs are in gene regulatory regions (upstream, 3'-UTR, and intron) modifying gene expression and protein-coding. This integrated approach of marker design for t-GBS alleles, SNP functional annotation, and high-throughput genotyping of functional SNP offers translation solutions across crops and complex traits in crop improvement programs.

Optimizing genomic selection of agricultural traits using K-wheat core collection

The agricultural traits that constitute basic plant breeding information are usually quantitative or complex in nature. This quantitative and complex combination of traits complicates the process of selection in breeding. This study examined the potential of genome-wide association studies (GWAS) and genomewide selection (GS) for breeding ten agricultural traits by using genome-wide SNPs. As a first step, a trait-associated candidate marker was identified by GWAS using a genetically diverse 567 Korean (K)-wheat core collection. The accessions were genotyped using an Axiom^® 35K wheat DNA chip, and ten agricultural traits were determined (awn color, awn length, culm color, culm length, ear color, ear length, days to heading, days to maturity, leaf length, and leaf width). It is essential to sustain global wheat production by utilizing accessions in wheat breeding. Among the traits associated with awn color and ear color that showed a high positive correlation, a SNP located on chr1B was significantly associated with both traits. Next, GS evaluated the prediction accuracy using six predictive models (G-BLUP, LASSO, BayseA, reproducing kernel Hilbert space, support vector machine (SVM), and random forest) and various training populations (TPs). With the exception of the SVM, all statistical models demonstrated a prediction accuracy of 0.4 or better. For the optimization of the TP, the number of TPs was randomly selected (10%, 30%, 50% and 70%) or divided into three subgroups (CC-sub 1, CC-sub 2 and CC-sub 3) based on the subpopulation structure. Based on subgroup-based TPs, better prediction accuracy was found for awn color, culm color, culm length, ear color, ear length, and leaf width. A variety of Korean wheat cultivars were used for validation to evaluate the prediction ability of populations. Seven out of ten cultivars showed phenotype-consistent results based on genomics-evaluated breeding values (GEBVs) calculated by the reproducing kernel Hilbert space (RKHS) predictive model. Our research provides a basis for improving complex traits in wheat breeding programs through genomics assisted breeding. The results of our research can be used as a basis for improving wheat breeding programs by using genomics-assisted breeding.

Current Progress and Future Prospect of Wheat Genetics Research towards an Enhanced Nitrogen Use Efficiency

To improve the yield and quality of wheat is of great importance for food security worldwide. One of the most effective and significant approaches to achieve this goal is to enhance the nitrogen use efficiency (NUE) in wheat. In this review, a comprehensive understanding of the factors involved in the process of the wheat nitrogen uptake, assimilation and remobilization of nitrogen in wheat were introduced. An appropriate definition of NUE is vital prior to its precise evaluation for the following gene identification and breeding process. Apart from grain yield (GY) and grain protein content (GPC), the commonly recognized major indicators of NUE, grain protein deviation (GPD) could also be considered as a potential trait for NUE evaluation. As a complex quantitative trait, NUE is affected by transporter proteins, kinases, transcription factors (TFs) and micro RNAs (miRNAs), which participate in the nitrogen uptake process, as well as key enzymes, circadian regulators, cross-talks between carbon metabolism, which are associated with nitrogen assimilation and remobilization. A series of quantitative genetic loci (QTLs) and linking markers were compiled in the hope to help discover more efficient and useful genetic resources for breeding program. For future NUE improvement, an exploration for other criteria during selection process that incorporates morphological, physiological and biochemical traits is needed. Applying new technologies from phenomics will allow high-throughput NUE phenotyping and accelerate the breeding process. A combination of multi-omics techniques and the previously verified QTLs and molecular markers will facilitate the NUE QTL-mapping and novel gene identification.

A modified Agrobacterium-mediated transformation for two oomycete pathogens

Oomycetes are a group of filamentous microorganisms that include some of the biggest threats to food security and natural ecosystems. However, much of the molecular basis of the pathogenesis and the development in these organisms remains to be learned, largely due to shortage of efficient genetic manipulation methods. In this study, we developed modified transformation methods for two important oomycete species, Phytophthora infestans and Plasmopara viticola, that bring destructive damage in agricultural production. As part of the study, we established an improved Agrobacterium-mediated transformation (AMT) method by prokaryotic expression in Agrobacterium tumefaciens of AtVIP1 (VirE2-interacting protein 1), an Arabidopsis bZIP gene required for AMT but absent in oomycetes genomes. Using the new method, we achieved an increment in transformation efficiency in two P. infestans strains. We further obtained a positive GFP transformant of P. viticola using the modified AMT method. By combining this method with the CRISPR/Cas12a genome editing system, we successfully performed targeted mutagenesis and generated loss-of-function mutations in two P. infestans genes. We edited a MADS-box transcription factor-encoding gene and found that a homozygous mutation in MADS-box results in poor sporulation and significantly reduced virulence. Meanwhile, a single-copy avirulence effector-encoding gene Avr8 in P. infestans was targeted and the edited transformants were virulent on potato carrying the cognate resistance gene R8, suggesting that loss of Avr8 led to successful evasion of the host immune response by the pathogen. In summary, this study reports on a modified genetic transformation and genome editing system, providing a potential tool for accelerating molecular genetic studies not only in oomycetes, but also other microorganisms.

WheatCENet: A Database for Comparative Co-expression Networks Analysis of Allohexaploid Wheat and Its Progenitors

Genetic and epigenetic changes after polyploidization events could result in variable gene expression and modified regulatory networks. Here, using large-scale transcriptome data, we constructed co-expression networks for diploid, tetraploid, and hexaploid wheat species, and built a platform for comparing co-expression networks of allohexaploid wheat and its progenitors, named WheatCENet. WheatCENet is a platform for searching and comparing specific functional co-expression networks, as well as identifying the related functions of the genes clustered therein. Functional annotations like pathways, gene families, protein-protein interactions, microRNAs (miRNAs), and several lines of epigenome data are integrated into this platform, and Gene Ontology (GO) annotation, gene set enrichment analysis (GSEA), motif identification, and other useful tools are also included. Using WheatCENet, we found that the network of WHEAT ABERRANT PANICLE ORGANIZATION 1 (WAPO1) has more co-expressed genes related to spike development in hexaploid wheat than its progenitors. We also found a novel motif of CCWWWWWWGG (CArG) specifically in the promoter region of WAPO-A1, suggesting that neofunctionalization of the WAPO-A1 gene affects spikelet development in hexaploid wheat. WheatCENet is useful for investigating co-expression networks and conducting other analyses, and thus facilitates comparative and functional genomic studies in wheat. WheatCENet is freely available at http://bioinformatics.cpolar.cn/WheatCENet and http://bioinformatics.cau.edu.cn/WheatCENet.

Transcriptome Dynamic Analysis Reveals New Candidate Genes Associated with Resistance to Fusarium Head Blight in Two Chinese Contrasting Wheat Genotypes

In recent years, Fusarium head blight (FHB) has developed into a global disease that seriously affects the yield and quality of wheat. Effective measures to solve this problem include exploring disease-resistant genes and breeding disease-resistant varieties. In this study, we conducted a comparative transcriptome analysis to identify the important genes that are differentially expressed in FHB medium-resistant (Nankang 1) and FHB medium-susceptible (Shannong 102) wheat varieties for various periods after Fusarium graminearum infection using RNA-seq technology. In total, 96,628 differentially expressed genes (DEGs) were identified, 42,767 from Shannong 102 and 53,861 from Nankang 1 (FDR < 0.05 and |log2FC| > 1). Of these, 5754 and 6841 genes were found to be shared among the three time points in Shannong 102 and Nankang 1, respectively. After inoculation for 48 h, the number of upregulated genes in Nankang 1 was significantly lower than that of Shannong 102, but at 96 h, the number of DEGs in Nankang 1 was higher than that in Shannong 102. This indicated that Shannong 102 and Nankang 1 had different defensive responses to F. graminearum in the early stages of infection. By comparing the DEGs, there were 2282 genes shared at the three time points between the two strains. GO and KEGG analyses of these DEGs showed that the following pathways were associated with disease resistance genes: response to stimulus pathway in GO, glutathione metabolism, phenylpropanoid biosynthesis, plant hormone signal transduction, and plant-pathogen interaction in KEGG. Among them, 16 upregulated genes were identified in the plant-pathogen interaction pathway. There were five upregulated genes, TraesCS5A02G439700, TraesCS5B02G442900, TraesCS5B02G443300, TraesCS5B02G443400, and TraesCS5D02G446900, with significantly higher expression levels in Nankang 1 than in Shannong 102, and these genes may have an important role in regulating the resistance of Nankang 1 to F. graminearum infection. The PR proteins they encode are PR protein 1-9, PR protein 1-6, PR protein 1-7, PR protein 1-7, and PR protein 1-like. In addition, the number of DEGs in Nankang 1 was higher than that in Shannong 102 on almost all chromosomes, except chromosomes 1A and 3D, but especially on chromosomes 6B, 4B, 3B, and 5A. These results indicate that gene expression and the genetic background must be considered for FHB resistance in wheat breeding.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

The quest for optimal plant architecture

Changes in plant architecture can improve cereal crop yield.

Progenitor species hold untapped diversity for potential climate-responsive traits for use in wheat breeding and crop improvement

Climate change will have numerous impacts on crop production worldwide necessitating a broadening of the germplasm base required to source and incorporate novel traits. Major variation exists in crop progenitor species for seasonal adaptation, photosynthetic characteristics, and root system architecture. Wheat is crucial for securing future food and nutrition security and its evolutionary history and progenitor diversity offer opportunities to mine favourable functional variation in the primary gene pool. Here we provide a review of the status of characterisation of wheat progenitor variation and the potential to use this knowledge to inform the use of variation in other cereal crops. Although significant knowledge of progenitor variation has been generated, we make recommendations for further work required to systematically characterise underlying genetics and physiological mechanisms and propose steps for effective use in breeding. This will enable targeted exploitation of useful variation, supported by the growing portfolio of genomics and accelerated breeding approaches. The knowledge and approaches generated are also likely to be useful across wider crop improvement.

Elucidating the Drought Responsiveness in Wheat Genotypes

Drought is one of the major abiotic constraints on wheat yields and also for sustainability of production levels around the world. In the near future, the occurrence likelihood of droughts is predicted to become more common, due to changing climatic conditions, thereby posing a serious threat to the food security system. Heterogeneity, in its time of occurrence and severity levels, is likely to further augment the complexity of drought conditions. Although wheat crop growth has progressively risen to good levels, as evident by notable increases in both area and production, the expected wheat demand for the ever-growing population is quite high. Besides crop yield volatility in the era of climate change and dwindling resources, “trait-based” breeding programs are required, so as to develop high yielding, climate resilient and stable genotypes, at a faster pace. For this to happen, a broad genetic base and wider adaptability to suit varied agro-ecologies would provide enough scope for their quicker spread. The current review places emphasis on making distinct categories of the wheat cultivars/advanced breeding lines, as tolerant, moderately tolerant or susceptible to drought stresses, duly supported by an extensive up-to-date literature base and will be useful for wheat researchers, in order to choose the best potential donors as parents, coupled with the associated traits for the development of drought-tolerant wheat varieties, and also to facilitate molecular studies.

The Antarctic Moss Pohlia nutans Genome Provides Insights Into the Evolution of Bryophytes and the Adaptation to Extreme Terrestrial Habitats

The Antarctic continent has extreme natural environment and fragile ecosystem. Mosses are one of the dominant floras in the Antarctic continent. However, their genomic features and adaptation processes to extreme environments remain poorly understood. Here, we assembled the high-quality genome sequence of the Antarctic moss (Pohlia nutans) with 698.20 Mb and 22 chromosomes. We found that the high proportion of repeat sequences and a recent whole-genome duplication (WGD) contribute to the large size genome of P. nutans when compared to other bryophytes. The genome of P. nutans harbors the signatures of massive segmental gene duplications and large expansions of gene families, likely facilitating neofunctionalization. Genomic characteristics that may support the Antarctic lifestyle of this moss comprise expanded gene families involved in phenylpropanoid biosynthesis, unsaturated fatty acid biosynthesis, and plant hormone signal transduction. Additional contributions include the significant expansion and upregulation of several genes encoding DNA photolyase, antioxidant enzymes, flavonoid biosynthesis enzymes, possibly reflecting diverse adaptive strategies. Notably, integrated multi-omic analyses elucidate flavonoid biosynthesis may function as the reactive oxygen species detoxification under UV-B radiation. Our studies provide insight into the unique features of the Antarctic moss genome and their molecular responses to extreme terrestrial environments.

Resistance of the Wheat Cultivar ‘Renan’ to Septoria Leaf Blotch Explained by a Combination of Strain Specific and Strain Non-Specific QTL Mapped on an Ultra-Dense Genetic Map

Quantitative resistance is considered more durable than qualitative resistance as it does not involve major resistance genes that can be easily overcome by pathogen populations, but rather a combination of genes with a lower individual effect. This durability means that quantitative resistance could be an interesting tool for breeding crops that would not systematically require phytosanitary products. Quantitative resistance has yet to reveal all of its intricacies. Here, we delve into the case of the wheat/Septoria tritici blotch (STB) pathosystem. Using a population resulting from a cross between French cultivar Renan, generally resistant to STB, and Chinese Spring, a cultivar susceptible to the disease, we built an ultra-dense genetic map that carries 148,820 single nucleotide polymorphism (SNP) markers. Phenotyping the interaction was done with two different Zymoseptoria tritici strains with contrasted pathogenicities on Renan. A linkage analysis led to the detection of three quantitative trait loci (QTL) related to resistance in Renan. These QTL, on chromosomes 7B, 1D, and 5D, present with an interesting diversity as that on 7B was detected with both fungal strains, while those on 1D and 5D were strain-specific. The resistance on 7B was located in the region of Stb8 and the resistance on 1D colocalized with Stb19. However, the resistance on 5D was new, so further designated Stb20q. Several wall-associated kinases (WAK), nucleotide-binding and leucine-rich repeats (NB-LRR) type, and kinase domain carrying genes were present in the QTL regions, and some of them were expressed during the infection. These results advocate for a role of Stb genes in quantitative resistance and for resistance in the wheat/STB pathosystem being as a whole quantitative and polygenic.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Histone acetyltransferase TaHAG1 acts as a crucial regulator to strengthen salt tolerance of hexaploid wheat

Polyploidy occurs prevalently and plays an important role during plant speciation and evolution. This phenomenon suggests polyploidy could develop novel features that enable them to adapt wider range of environmental conditions compared with diploid progenitors. Bread wheat (Triticum aestivum L., BBAADD) is a typical allohexaploid species and generally exhibits greater salt tolerance than its tetraploid wheat progenitor (BBAA). However, little is known about the underlying molecular basis and the regulatory pathway of this trait. Here, we show that the histone acetyltransferase TaHAG1 acts as a crucial regulator to strengthen salt tolerance of hexaploid wheat. Salinity-induced TaHAG1 expression was associated with tolerance variation in polyploidy wheat. Overexpression, silencing, and CRISPR-mediated knockout of TaHAG1 validated the role of TaHAG1 in salinity tolerance of wheat. TaHAG1 contributed to salt tolerance by modulating reactive oxygen species (ROS) production and signal specificity. Moreover, TaHAG1 directly targeted a subset of genes that are responsible for hydrogen peroxide production, and enrichment of TaHAG1 triggered increased H3 acetylation and transcriptional upregulation of these loci under salt stress. In addition, we found the salinity-induced TaHAG1-mediated ROS production pathway is involved in salt tolerance difference of wheat accessions with varying ploidy. Our findings provide insight into the molecular mechanism of how an epigenetic regulatory factor facilitates adaptability of polyploidy wheat and highlights this epigenetic modulator as a strategy for salt tolerance breeding in bread wheat.

Myb10-D confers PHS-3D resistance to pre-harvest sprouting by regulating NCED in ABA biosynthesis pathway of wheat

Pre-harvest sprouting (PHS), the germination of grain before harvest, is a serious problem resulting in wheat yield and quality losses. Here, we mapped the PHS resistance gene PHS-3D from synthetic hexaploid wheat to a 2.4 Mb presence-absence variation (PAV) region and found that its resistance effect was attributed to the pleiotropic Myb10-D by integrated omics and functional analyses. Three haplotypes were detected in this PAV region among 262 worldwide wheat lines and 16 Aegilops tauschii, and the germination percentages of wheat lines containing Myb10-D was approximately 40% lower than that of the other lines. Transcriptome and metabolome profiling indicated that Myb10-D affected the transcription of genes in both the flavonoid and abscisic acid (ABA) biosynthesis pathways, which resulted in increases in flavonoids and ABA in transgenic wheat lines. Myb10-D activates 9-cis-epoxycarotenoid dioxygenase (NCED) by biding the secondary wall MYB-responsive element (SMRE) to promote ABA biosynthesis in early wheat seed development stages. We revealed that the newly discovered function of Myb10-D confers PHS resistance by enhancing ABA biosynthesis to delay germination in wheat. The PAV harboring Myb10-D associated with grain color and PHS will be useful for understanding and selecting white grained PHS resistant wheat cultivars.

Ploidy influences wheat mesophyll cell geometry, packing and leaf function

Leaf function is influenced by leaf structure, which is itself related not only to the spatial arrangement of constituent mesophyll cells, but also their size and shape. In this study, we used confocal microscopy to image leaves of Triticum genotypes varying in ploidy level to extract 3D information on individual mesophyll cell size and geometry. Combined with X-ray Computed Tomography and gas exchange analysis, the effect of changes in wheat mesophyll cell geometry upon leaf structure and function were investigated. Mesophyll cell size and shape were found to have changed during the course of wheat evolution. An unexpected linear relationship between mesophyll cell surface area and volume was discovered, suggesting anisotropic scaling of mesophyll cell geometry with increasing ploidy. Altered mesophyll cell size and shape were demonstrated to be associated with changes in mesophyll tissue architecture. Under experimental growth conditions, CO2 assimilation did not vary with ploidy, but stomatal conductance was lower in hexaploid plants, conferring a greater instantaneous water-use efficiency. We propose that as wheat mesophyll cells have become larger with increased ploidy, this has been accompanied by changes in cell geometry and packing which limit water loss while maintaining carbon assimilation.

Differences between diploid donors are the main contributing factor for subgenome asymmetry measured in either gene ratio or relative diversity in allopolyploids

Subgenome asymmetry (SA) has routinely been attributed to different responses between the subgenomes of a polyploid to various stimuli during evolution. Here, we compared subgenome differences in gene ratio and relative diversity between artificial and natural genotypes of several allopolyploid species. Surprisingly, consistent differences were not detected between these two types of polyploid genotypes, although they differ in times exposed to evolutionary selection. The estimated ratio of shared genes between a subgenome and its diploid donor was invariably higher for the artificial allopolyploid genotypes than those for the natural genotypes, which is expected as it is now well-known that many genes in a species are not shared among all individuals. As the exact diploid parent for a given subgenome is unknown, the estimated ratios of shared genes for the natural genotypes would also include difference among individual genotypes of the diploid donor species. Further, we detected the presence of SA in genotypes before the completion of the polyploidization events as well as in those which were not formed via polyploidization. These results indicate that SA may, to a large degree, reflect differences between its diploid donors or that changes occurred during polyploid evolution are defined by their donor genomes.

Genome-Wide and Exome-Capturing Sequencing of a Gamma-Ray-Induced Mutant Reveals Biased Variations in Common Wheat

Induced mutagenesis is a powerful approach for the creation of novel germplasm and the improvement of agronomic traits. The evaluation of mutagenic effects and functional variations in crops is needed for breeding mutant strains. To investigate the mutagenic effects of gamma-ray irradiation in wheat, this study characterized genomic variations of wheat early heading mutant (eh1) as compared to wild-type (WT) Zhongyuan 9 (ZY9). Whole-genome resequencing of eh1 and ZY9 produced 737.7 Gb sequencing data and identified a total of 23,537,117 homozygous single nucleotide polymorphism (SNP) and 1,608,468 Indel. Analysis of SNP distribution across the chromosome suggests that mutation hotspots existed in certain chromosomal regions. Among the three subgenomes, the variation frequency in subgenome D was significantly lower than in subgenomes A and B. A total of 27.8 Gb data were obtained by exome-capturing sequencing, while 217,948 SNP and 13,554 Indel were identified. Variation annotation in the gene-coding sequences demonstrated that 5.0% of the SNP and 5.3% of the Indel were functionally important. Characterization of exomic variations in 12 additional gamma-ray-induced mutant lines further provided additional insights into the mutagenic effects of this approach. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) analysis suggested that genes with functional variations were enriched in several metabolic pathways, including plant-pathogen interactions and ADP binding. Kompetitive allele-specific PCR (KASP) genotyping with selected SNP within functional genes indicated that 85.7% of the SNPs were polymorphic between the eh1 and wild type. This study provides a basic understanding of the mechanism behind gamma-ray irradiation in hexaploid wheat.

An Efficient Brome mosaic virus-Based Gene Silencing Protocol for Hexaploid Wheat (Triticum aestivum L.)

Virus-induced gene silencing (VIGS) is a rapid and powerful method to evaluate gene function, especially for species like hexaploid wheat that have large, redundant genomes and are difficult and time-consuming to transform. The Brome mosaic virus (BMV)-based VIGS vector is widely used in monocotyledonous species but not wheat. Here we report the establishment of a simple and effective VIGS procedure in bread wheat using BMVCP5, the most recently improved BMV silencing vector, and wheat genes PHYTOENE DESATURASE (TaPDS) and PHOSPHATE2 (TaPHO2) as targets. Time-course experiments revealed that smaller inserts (~100 nucleotides, nt) were more stable in BMVCP5 and conferred higher silencing efficiency and longer silencing duration, compared with larger inserts. When using a 100-nt insert and a novel coleoptile inoculation method, BMVCP5 induced extensive silencing of TaPDS transcript and a visible bleaching phenotype in the 2nd to 5th systemically-infected leaves from nine to at least 28 days post inoculation (dpi). For TaPHO2, the ability of BMVCP5 to simultaneously silence all three homoeologs was demonstrated. To investigate the feasibility of BMV VIGS in wheat roots, ectopically expressed enhanced GREEN FLUORESCENT PROTEIN (eGFP) in a transgenic wheat line was targeted for silencing. Silencing of eGFP fluorescence was observed in both the maturation and elongation zones of roots. BMVCP5 mediated significant silencing of eGFP and TaPHO2 mRNA expression in roots at 14 and 21 dpi, and TaPHO2 silencing led to the doubling of inorganic phosphate concentration in the 2nd through 4th systemic leaves. All 54 wheat cultivars screened were susceptible to BMV infection. BMVCP5-mediated TaPDS silencing resulted in the expected bleaching phenotype in all eight cultivars examined, and decreased TaPDS transcript was detected in all three cultivars examined. This BMVCP5 VIGS technology may serve as a rapid and effective functional genomics tool for high-throughput gene function studies in aerial and root tissues and in many wheat cultivars.

Genetic Mapping and Validation of Loci for Kernel-Related Traits in Wheat (Triticum aestivum L.)

Kernel size (KS) and kernel weight play a key role in wheat yield. Phenotypic data from six environments and a Wheat55K single-nucleotide polymorphism array-based constructed genetic linkage map from a recombinant inbred line population derived from the cross between the wheat line 20828 and the line SY95-71 were used to identify quantitative trait locus (QTL) for kernel length (KL), kernel width (KW), kernel thickness (KT), thousand-kernel weight (TKW), kernel length-width ratio (LWR), KS, and factor form density (FFD). The results showed that 65 QTLs associated with kernel traits were detected, of which the major QTLs QKL.sicau-2SY-1B, QKW.sicau-2SY-6D, QKT.sicau-2SY-2D, and QTKW.sicau-2SY-2D, QLWR.sicau-2SY-6D, QKS.sicau-2SY-1B/2D/6D, and QFFD.sicau-2SY-2D controlling KL, KW, KT, TKW, LWR, KS, and FFD, and identified in multiple environments, respectively. They were located on chromosomes 1BL, 2DL, and 6DS and formed three QTL clusters. Comparison of genetic and physical interval suggested that only QKL.sicau-2SY-1B located on chromosome 1BL was likely a novel QTL. A Kompetitive Allele Specific Polymerase chain reaction (KASP) marker, KASP-AX-109379070, closely linked to this novel QTL was developed and used to successfully confirm its effect in two different genetic populations and three variety panels consisting of 272 Chinese wheat landraces, 300 Chinese wheat cultivars most from the Yellow and Huai River Valley wheat region, and 165 Sichuan wheat cultivars. The relationships between kernel traits and other agronomic traits were detected and discussed. A few predicted genes involved in regulation of kernel growth and development were identified in the intervals of these identified major QTL. Taken together, these stable and major QTLs provide valuable information for understanding the genetic composition of kernel yield and provide the basis for molecular marker-assisted breeding.

A Collinearity-Incorporating Homology Inference Strategy for Connecting Emerging Assemblies in the Triticeae Tribe as a Pilot Practice in the Plant Pangenomic Era

Plant genome sequencing has dramatically increased, and some species even have multiple high-quality reference versions. Demands for clade-specific homology inference and analysis have increased in the pangenomic era. Here we present a novel method, GeneTribe (https://chenym1.github.io/genetribe/), for homology inference among genetically similar genomes that incorporates gene collinearity and shows better performance than traditional sequence-similarity-based methods in terms of accuracy and scalability. The Triticeae tribe is a typical allopolyploid-rich clade with complex species relationships that includes many important crops, such as wheat, barley, and rye. We built Triticeae-GeneTribe (http://wheat.cau.edu.cn/TGT/), a homology database, by integrating 12 Triticeae genomes and 3 outgroup model genomes and implemented versatile analysis and visualization functions. With macrocollinearity analysis, we were able to construct a refined model illustrating the structural rearrangements of the 4A-5A-7B chromosomes in wheat as two major translocation events. With collinearity analysis at both the macro- and microscale, we illustrated the complex evolutionary history of homologs of the wheat vernalization gene Vrn2, which evolved as a combined result of genome translocation, duplication, and polyploidization and gene loss events. Our work provides a useful practice for connecting emerging genome assemblies, with awareness of the extensive polyploidy in plants, and will help researchers efficiently exploit genome sequence resources.

Multiple wheat genomes reveal global variation in modern breeding

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome¹, and the lack of genome-assembly data for multiple wheat lines^2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses^4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm1⁶, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Comparison of multiple genome assemblies from wheat reveals extensive diversity that results from the complex breeding history of wheat and provides a basis for further potential improvements to this important food crop.

Population structure of Nepali spring wheat (Triticum aestivum L.) germplasm

BACKGROUND

Appropriate information about genetic diversity and population structure of germplasm improves the efficiency of plant breeding. The low productivity of Nepali bread wheat (Triticum aestivum L.) is a major concern particularly since Nepal is ranked the 4th most vulnerable nation globally to climate change. The genetic diversity and population structure of Nepali spring wheat have not been reported. This study aims to improve the exploitation of more diverse and under-utilized genetic resources to contribute to current and future breeding efforts for global food security.

RESULTS

We used genotyping-by-sequencing (GBS) to characterize a panel of 318 spring wheat accessions from Nepal including 166 landraces, 115 CIMMYT advanced lines, and 34 Nepali released varieties. We identified 95 K high-quality SNPs. The greatest genetic diversity was observed among the landraces, followed by CIMMYT lines, and released varieties. Though we expected only 3 groupings corresponding to these 3 seed origins, the population structure revealed two large, distinct subpopulations along with two smaller and scattered subpopulations in between, with significant admixture. This result was confirmed by principal component analysis (PCA) and UPGMA distance-based clustering. The pattern of LD decay differed between subpopulations, ranging from 60 to 150 Kb. We discuss the possibility that germplasm explorations during the 1970s-1990s may have mistakenly collected exotic germplasm instead of local landraces and/or collected materials that had already cross-hybridized since exotic germplasm was introduced starting in the 1950s.

CONCLUSION

We suggest that only a subset of wheat "landraces" in Nepal are authentic which this study has identified. Targeting these authentic landraces may accelerate local breeding programs to improve the food security of this climate-vulnerable nation. Overall, this study provides a novel understanding of the genetic diversity of wheat in Nepal and this may contribute to global wheat breeding initiatives.

High-Resolution Genome-wide Association Study Identifies Genomic Regions and Candidate Genes for Important Agronomic Traits in Wheat

Wheat (Triticum aestivum) is a major staple food crop worldwide. Genetic dissection of important agronomic traits is essential for continuous improvement of wheat yield to meet the demand of the world's growing population. We conducted a large-scale genome-wide association study (GWAS) using a panel of 768 wheat cultivars that were genotyped with 327 609 single-nucleotide polymorphisms generated by genotyping-by-sequencing and detected 395 quantitative trait loci (QTLs) for 12 traits under 7 environments. Among them, 273 QTLs were delimited to ≤1.0-Mb intervals and 7 of them are either known genes (Rht-D, Vrn-B1, and Vrn-D1) that have been cloned or known QTLs (TaGA2ox8, APO1, TaSus1-7B, and Rht12) that were previously mapped. Eight putative candidate genes were identified for three QTLs that enhance spike seed setting and grain size using gene expression data and were validated in three bi-parental populations. Protein sequence analysis identified 33 putative wheat orthologs that have high identity with rice genes in QTLs affecting similar traits. Large r² values for additive effects observed among the QTLs for most traits indicated that the phenotypes of these identified QTLs were highly predictable. Results from this study demonstrated that significantly increasing GWAS population size and marker density greatly improves detection and identification of candidate genes underlying a QTL, solidifying the foundation for large-scale QTL fine mapping, candidate gene validation, and developing functional markers for genomics-based breeding in wheat.

Recent Progress in Germplasm Evaluation and Gene Mapping to Enable Breeding of Drought-Tolerant Wheat

There is a need to increase wheat productivity to meet the food demands of the ever-growing human population. However, accelerated development of high yielding varieties is hindered by drought, which is worsening due to climate change. In this context, germplasm diversity is central to the development of drought-tolerant wheat. Extensive collections of these genetic resources are conserved in national and international genebanks. In addition to phenotypic assessments, the use of advanced molecular techniques (e.g., genotype by sequencing) to identify quantitative trait loci (QTLs) for drought tolerance related traits is useful for genome- and marker-assisted selection based approaches. Therefore, to assist wheat breeders at a critical time, we searched the recent peer-reviewed literature (2011-current), first, to identify wheat germplasm observed to be useful genetic sources for drought tolerance, and second, to report QTLs associated with drought tolerance. Though many breeders limit the parents used in breeding programs to a familiar core collection, the results of this review show that larger germplasm collections have been sources of useful genes for drought tolerance in wheat. The review also demonstrates that QTLs for drought tolerance in wheat are associated with diverse physio-morphological traits, at different growth stages. Here, we also briefly discuss the potential of genome engineering/editing to improve drought tolerance in wheat. The use of CRISPR-Cas9 and other gene-editing technologies can be used to fine-tune the expression of genes controlling drought adaptive traits, while high throughput phenotyping (HTP) techniques can potentially accelerate the selection process. These efforts are empowered by wheat researcher consortia.

Molecular and genetic bases of heat stress responses in crop plants and breeding for increased resilience and productivity

To ensure the food security of future generations and to address the challenge of the 'no hunger zone' proposed by the FAO (Food and Agriculture Organization), crop production must be doubled by 2050, but environmental stresses are counteracting this goal. Heat stress in particular is affecting agricultural crops more frequently and more severely. Since the discovery of the physiological, molecular, and genetic bases of heat stress responses, cultivated plants have become the subject of intense research on how they may avoid or tolerate heat stress by either using natural genetic variation or creating new variation with DNA technologies, mutational breeding, or genome editing. This review reports current understanding of the genetic and molecular bases of heat stress in crops together with recent approaches to creating heat-tolerant varieties. Research is close to a breakthrough of global relevance, breeding plants fitter to face the biggest challenge of our time.

The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat

The rapid development and application of molecular marker assays have facilitated genomic selection and genome-wide linkage and association studies in wheat breeding. Although PCR-based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker-assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single-nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome-specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost-effective and may be the best choice for targeted genotyping and marker-assisted selection in wheat genetic improvement.

Genomic and functional genomics analyses of gluten proteins and prospect for simultaneous improvement of end-use and health-related traits in wheat

KEY MESSAGE

Recent genomic and functional genomics analyses have substantially improved the understanding on gluten proteins, which are important determinants of wheat grain quality traits. The new insights obtained and the availability of precise, versatile and high-throughput genome editing technologies will accelerate simultaneous improvement of wheat end-use and health-related traits. Being a major staple food crop in the world, wheat provides an indispensable source of dietary energy and nutrients to the human population. As worldwide population grows and living standards rise in both developed and developing countries, the demand for wheat with high quality attributes increases globally. However, efficient breeding of high-quality wheat depends on critically the knowledge on gluten proteins, which mainly include several families of prolamin proteins specifically accumulated in the endospermic tissues of grains. Although gluten proteins have been studied for many decades, efficient manipulation of these proteins for simultaneous enhancement of end-use and health-related traits has been difficult because of high complexities in their expression, function and genetic variation. However, recent genomic and functional genomics analyses have substantially improved the understanding on gluten proteins. Therefore, the main objective of this review is to summarize the genomic and functional genomics information obtained in the last 10 years on gluten protein chromosome loci and genes and the cis- and trans-factors regulating their expression in the grains, as well as the efforts in elucidating the involvement of gluten proteins in several wheat sensitivities affecting genetically susceptible human individuals. The new insights gathered, plus the availability of precise, versatile and high-throughput genome editing technologies, promise to speed up the concurrent improvement of wheat end-use and health-related traits and the development of high-quality cultivars for different consumption needs.

Where the Wild Things Are: Transposable Elements as Drivers of Structural and Functional Variations in the Wheat Genome

Transposable elements (TEs) are major contributors to genome plasticity and thus are likely to have a dramatic impact on genetic diversity and speciation. Recent technological developments facilitated the sequencing and assembly of the wheat genome, opening the gate for whole genome analysis of TEs in wheat, which occupy over 80% of the genome. Questions that have been long unanswered regarding TE dynamics throughout the evolution of wheat, are now being addressed more easily, while new questions are rising. In this review, we discuss recent advances in the field of TE dynamics in wheat and possible future directions.

Biolistic Delivery of Programmable Nuclease (CRISPR/Cas9) in Bread Wheat

The discovery of site-specific programmable nucleases has led to a major breakthrough in the area of genome editing. In the past few years, CRISPR/Cas system has been utilized for genome editing of a large number of crops including cereals like wheat, rice, maize, and barley. In terms of consumption, wheat is second only to rice as the most important crop of the world. In the present chapter, we describe biolistic delivery method of ribonucleoprotein (RNP) complexes of programmable nuclease (CRISPR/Cas9) for targeted genome editing and selection-free screening of transformants in wheat. The method not only overcomes the problem of random integration into the genome but also reduces the off-targets. Besides the step-by-step protocol, plausible challenges and ways to overcome them are also discussed. By using the described method of biolistic delivery of CRISPR/Cas9 in plant systems, genome-edited plants can be identified within 11 weeks.

Combining a New Exome Capture Panel With an Effective varBScore Algorithm Accelerates BSA-Based Gene Cloning in Wheat

The discovery of functional genes underlying agronomic traits is of great importance for wheat improvement. Here we designed a new wheat exome capture probe panel based on IWGSC RefSeq v1.0 genome sequence information and developed an effective algorithm, varBScore, that can sufficiently reduce the background noise in gene mapping and identification. An effective method, termed bulked segregant exome capture sequencing (BSE-Seq) for identifying causal mutations or candidate genes was established by combining the use of a newly designed wheat exome capture panel, sequencing of bulked segregant pools from segregating populations, and the robust algorithm varBScore. We evaluated the effectiveness of varBScore on SNP calling using the published dataset for mapping and cloning the yellow rust resistance gene Yr7 in wheat. Furthermore, using BSE-Seq, we rapidly identified a wheat yellow leaf mutant gene, ygl1, in an ethyl methanesulfonate (EMS) mutant population and found that a single mutation of G to A at 921 position in the wild type YGL1 gene encoding magnesium-chelatase subunit chlI caused the leaf yellowing phenotype. We further showed that mutation of YGL1 through CRISPR/Cas9 gene editing led to a yellow phenotype on the leaves of transgenic wheat, indicating that ygl1 is the correct causal gene responsible for the mutant phenotype. In summary, our approach is highly efficient for discovering causal mutations and gene cloning in wheat.

The power of model-to-crop translation illustrated by reducing seed loss from pod shatter in oilseed rape

Elucidation of key regulators in Arabidopsis fruit patterning has facilitated knowledge-translation into crop species to address yield loss caused by premature seed dispersal (pod shatter). In the 1980s, plant scientists descended on a small weed Arabidopsis thaliana (thale cress) and developed it into a powerful model system to study plant biology. The massive advances in genetics and genomics since then have allowed us to obtain incredibly detailed knowledge on specific biological processes of Arabidopsis growth and development, its genome sequence and the function of many of the individual genes. This wealth of information provides immense potential for translation into crops to improve their performance and address issues of global importance such as food security. Here, we describe how fundamental insight into the genetic mechanism by which seed dispersal occurs in members of the Brassicaceae family can be exploited to reduce seed loss in oilseed rape (Brassica napus). We demonstrate that by exploiting data on gene function in model species, it is possible to adjust the pod-opening process in oilseed rape, thereby significantly increasing yield. Specifically, we identified mutations in multiple paralogues of the INDEHISCENT and GA4 genes in B. napus and have overcome genetic redundancy by combining mutant alleles. Finally, we present novel software for the analysis of pod shatter data that is applicable to any crop for which seed dispersal is a serious problem. These findings highlight the tremendous potential of fundamental research in guiding strategies for crop improvement.

Identification of quantitative trait loci for kernel traits in a wheat cultivar Chuannong16

Background

Kernel length (KL), kernel width (KW) and thousand-kernel weight (TKW) are key agronomic traits in wheat breeding. Chuannong16 (‘CN16’) is a commercial cultivar with significantly longer kernels than the line ‘20828’. To identify and characterize potential alleles from CN16 controlling KL, the previously developed recombinant inbred line (RIL) population derived from the cross ‘20828’ × ‘CN16’ and the genetic map constructed by the Wheat55K SNP array and SSR markers were used to perform quantitative trait locus/loci (QTL) analyses for kernel traits.

Results

A total of 11 putative QTL associated with kernel traits were identified and they were located on chromosomes 1A (2 QTL), 2B (2 QTL), 2D (3 QTL), 3D, 4A, 6A, and 7A, respectively. Among them, three major QTL, QKL.sicau-2D, QKW.sicau-2D and QTKW.sicau-2D, controlling KL, KW and TKW, respectively, were detected in three different environments. Respectively, they explained 10.88–18.85%, 17.21–21.49% and 10.01–23.20% of the phenotypic variance. Further, they were genetically mapped in the same interval on chromosome 2DS. A previously developed kompetitive allele-specific PCR (KASP) marker KASP-AX-94721936 was integrated in the genetic map and QTL re-mapping finally located the three major QTL in a 1- cM region flanked by AX-111096297 and KASP-AX-94721936. Another two co-located QTL intervals for KL and TKW were also identified. A few predicted genes involved in regulation of kernel growth and development were identified in the intervals of these identified QTL. Significant relationships between kernel traits and spikelet number per spike and anthesis date were detected and discussed.

Conclusions

Three major and stably expressed QTL associated with KL, KW, and TKW were identified. A KASP marker tightly linked to these three major QTL was integrated. These findings provide information for subsequent fine mapping and cloning the three co-localized major QTL for kernel traits.

Use of near-isogenic lines to precisely map and validate a major QTL for grain weight on chromosome 4AL in bread wheat (Triticum aestivum L.)

This study precisely mapped and validated a major quantitative trait locus (QTL) on chromosome 4AL for thousand-grain weight in wheat using multiple near-isogenic lines. Thousand-grain weight (TGW) is an essential yield component. Following the previous identification of a major QTL for TGW within the interval of 15.7 cM (92.7-108.4 cM) on chromosome 4AL using the Nongda3338 (ND3338)/Jingdong6 (JD6) doubled haploid population, the aim of this study was to perform more precise mapping and validate the genetic effect of the QTL. Multiple near-isogenic lines (NILs) were developed using ND3338 as the recurrent parent through marker-assisted selection. Based on five independent BC₃F_3:4 segregating populations derived from BC₃F₃ plants with different heterozygous segments for the target QTL site and the results of genotyping analysis performed using the Wheat660 K SNP array, it was possible to delimit the QTL region to a physical interval of approximately 6.5 Mb (677.11-683.61 Mb, IWGSC Ref Seq v1.0). Field trials across multiple environments showed that NILs^JD6 had a consistent effect on increasing the TGW by 5.16-27.48% and decreasing the grain number per spike (GNS) by 3.98-32.91% compared to the corresponding NILs^ND3338, which exhibited locus-specific TGW-GNS trade-offs. Moreover, by using RNA sequencing (RNA-Seq) of whole grains at 10 days after pollination stage of multiple NILs, we found that differentially expressed genes between the NIL pairs were significantly enriched for cell cycle and the replication of chromosome-related genes, hence affecting cell division and cell proliferation. Overall, our results provide a basis for map-based cloning of the major QTL and determining the mechanisms underlying TGW-GNS trade-offs in wheat, which would help to fine-tune these two components and maximize the grain yield for breeders.

Identification of a Dominant Chlorosis Phenotype Through a Forward Screen of the Triticum turgidum cv. Kronos TILLING Population

Durum wheat (Triticum turgidum) derives from a hybridization event approximately 400,000 years ago which led to the creation of an allotetraploid genome. The evolutionary recent origin of durum wheat means that its genome has not yet been fully diploidised. As a result, many of the genes present in the durum genome act in a redundant fashion, where loss-of-function mutations must be present in both gene copies to observe a phenotypic effect. Here, we use a novel set of induced variation within the cv. Kronos TILLING population to identify a locus controlling a dominant, environmentally dependent chlorosis phenotype. We carried out a forward screen of the sequenced cv. Kronos TILLING lines for senescence phenotypes and identified a line with a dominant early senescence and chlorosis phenotype. Mutant plants contained less chlorophyll throughout their development and displayed premature flag leaf senescence. A segregating population was classified into discrete phenotypic groups and subjected to bulked-segregant analysis using exome capture followed by next-generation sequencing. This allowed the identification of a single region on chromosome 3A, Yellow Early Senescence 1 (YES-1), which was associated with the mutant phenotype. While this phenotype was consistent across 4 years of field trials in the United Kingdom, the mutant phenotype was not observed when grown in Davis, CA (United States). To obtain further SNPs for fine-mapping, we isolated chromosome 3A using flow sorting and sequenced the entire chromosome. By mapping these reads against both the cv. Chinese Spring reference sequence and the cv. Kronos assembly, we could identify high-quality, novel EMS-induced SNPs in non-coding regions within YES-1 that were previously missed in the exome capture data. This allowed us to fine-map YES-1 to 4.3 Mb, containing 59 genes. Our study shows that populations containing induced variation can be sources of novel dominant variation in polyploid crop species, highlighting their importance in future genetic screens. We also demonstrate the value of using cultivar-specific genome assemblies alongside the gold-standard reference genomes particularly when working with non-coding regions of the genome. Further fine-mapping of the YES-1 locus will be pursued to identify the causal SNP underpinning this dominant, environmentally dependent phenotype.

A large nested association mapping population for breeding and quantitative trait locus mapping in Ethiopian durum wheat

The Ethiopian plateau hosts thousands of durum wheat (Triticum turgidum subsp. durum) farmer varieties (FV) with high adaptability and breeding potential. To harness their unique allelic diversity, we produced a large nested association mapping (NAM) population intercrossing fifty Ethiopian FVs with an international elite durum wheat variety (Asassa). The Ethiopian NAM population (EtNAM) is composed of fifty interconnected bi-parental families, totalling 6280 recombinant inbred lines (RILs) that represent both a powerful quantitative trait loci (QTL) mapping tool, and a large pre-breeding panel. Here, we discuss the molecular and phenotypic diversity of the EtNAM founder lines, then we use an array featuring 13 000 single nucleotide polymorphisms (SNPs) to characterize a subset of 1200 EtNAM RILs from 12 families. Finally, we test the usefulness of the population by mapping phenology traits and plant height using a genome wide association (GWA) approach. EtNAM RILs showed high allelic variation and a genetic makeup combining genetic diversity from Ethiopian FVs with the international durum wheat allele pool. EtNAM SNP data were projected on the fully sequenced AB genome of wild emmer wheat, and were used to estimate pairwise linkage disequilibrium (LD) measures that reported an LD decay distance of 7.4 Mb on average, and balanced founder contributions across EtNAM families. GWA analyses identified 11 genomic loci individually affecting up to 3 days in flowering time and more than 1.6 cm in height. We argue that the EtNAM is a powerful tool to support the production of new durum wheat varieties targeting local and global agriculture.

A major QTL co-localized on chromosome 6BL and its epistatic interaction for enhanced wheat stripe rust resistance

Co-localization of a major QTL for wheat stripe rust resistance to a 3.9-cM interval on chromosome 6BL across both populations and another QTL on chromosome 2B with epistatic interaction. Cultivars with diverse resistance are the optimal strategy to minimize yield losses caused by wheat stripe rust (Puccinia striiformis f. sp. tritici). Two wheat populations involving resistant wheat lines P10078 and Snb"S" from CIMMYT were evaluated for stripe rust response in multiple environments. Pool analysis by Wheat660K SNP array showed that the overlapping interval on chromosome 6B likely harbored a major QTL between two populations. Then, linkage maps were constructed using KASP markers, and a co-localized locus with large effect on chromosome 6BL was detected using QTL analysis in both populations. The coincident QTL, named QYr.nwafu-6BL.2, explained 59.7% of the phenotypic maximum variation in the Mingxian 169 × P10078 and 52.5% in the Zhengmai 9023 × Snb"S" populations, respectively. This co-localization interval spanning 3.9 cM corresponds to ~ 30.5-Mb genomic region of the newest common wheat reference genome (IWGSC RefSeq v.1.0). In addition, another QTL was also detected on chromosome 2B in Zhengmai 9023 × Snb"S" population and it can accelerate expression of QYr.nwafu-6BL.2 to enhance resistance with epistatic interaction. Allowing for Pst response, marker genotypes, pedigree analysis and relative genetic distance, QYr.nwafu-6BL.2 is likely to be a distinct adult plant resistance QTL. Haplotype analysis of QYr.nwafu-6BL.2 revealed specific SNPs or alleles in the target region from a diversity panel of 176 unrelated wheat accessions. This QTL region provides opportunity for further map-based cloning, and haplotypes analysis enables pyramiding favorable alleles into commercial cultivars by marker-assisted selection.

A reductionist approach to dissecting grain weight and yield in wheat

Grain yield is a highly polygenic trait that is influenced by the environment and integrates events throughout the life cycle of a plant. In wheat, the major grain yield components often present compensatory effects among them, which alongside the polyploid nature of wheat, makes their genetic and physiological study challenging. We propose a reductionist and systematic approach as an initial step to understand the gene networks regulating each individual yield component. Here, we focus on grain weight and discuss the importance of examining individual sub-components, not only to help in their genetic dissection, but also to inform our mechanistic understanding of how they interrelate. This knowledge should allow the development of novel combinations, across homoeologs and between complementary modes of action, thereby advancing towards a more integrated strategy for yield improvement. We argue that this will break barriers in terms of phenotypic variation, enhance our understanding of the physiology of yield, and potentially deliver improved on-farm yield.

From markers to genome-based breeding in wheat

Recent technological advances in wheat genomics provide new opportunities to uncover genetic variation in traits of breeding interest and enable genome-based breeding to deliver wheat cultivars for the projected food requirements for 2050. There has been tremendous progress in development of whole-genome sequencing resources in wheat and its progenitor species during the last 5 years. High-throughput genotyping is now possible in wheat not only for routine gene introgression but also for high-density genome-wide genotyping. This is a major transition phase to enable genome-based breeding to achieve progressive genetic gains to parallel to projected wheat production demands. These advances have intrigued wheat researchers to practice less pursued analytical approaches which were not practiced due to the short history of genome sequence availability. Such approaches have been successful in gene discovery and breeding applications in other crops and animals for which genome sequences have been available for much longer. These strategies include, (i) environmental genome-wide association studies in wheat genetic resources stored in genbanks to identify genes for local adaptation by using agroclimatic traits as phenotypes, (ii) haplotype-based analyses to improve the statistical power and resolution of genomic selection and gene mapping experiments, (iii) new breeding strategies for genome-based prediction of heterosis patterns in wheat, and (iv) ultimate use of genomics information to develop more efficient and robust genome-wide genotyping platforms to precisely predict higher yield potential and stability with greater precision. Genome-based breeding has potential to achieve the ultimate objective of ensuring sustainable wheat production through developing high yielding, climate-resilient wheat cultivars with high nutritional quality.

Developing a High-Throughput SNP-Based Marker System to Facilitate the Introgression of Traits From Aegilops Species Into Bread Wheat (Triticum aestivum)

The genus Aegilops contains a diverse collection of wild species exhibiting variation in geographical distribution, ecological adaptation, ploidy and genome organization. Aegilops is the most closely related genus to Triticum which includes cultivated wheat, a globally important crop that has a limited gene pool for modern breeding. Aegilops species are a potential future resource for wheat breeding for traits, such as adaptation to different ecological conditions and pest and disease resistance. This study describes the development and application of the first high-throughput genotyping platform specifically designed for screening wheat relative species. The platform was used to screen multiple accessions representing all species in the genus Aegilops. Firstly, the data was demonstrated to be useful for screening diversity and examining relationships within and between Aegilops species. Secondly, markers able to characterize and track introgressions from Aegilops species in hexaploid wheat were identified and validated using two different approaches.

An efficient and reproducible Agrobacterium-mediated transformation method for hexaploid wheat (Triticum aestivum L.)

BACKGROUND

Despite wheat being a worldwide staple, it is still considered the most difficult to transform out of the main cereal crops. Therefore, for the wheat research community, a freely available and effective wheat transformation system is still greatly needed.

RESULTS

We have developed and optimised a reproducible Agrobacterium-mediated transformation system for the spring wheat cv 'Fielder' that yields transformation efficiencies of up to 25%. We report on some of the important factors that influence transformation efficiencies. In particular, these include donor plant health, stage of the donor material, pre-treatment by centrifugation, vector type and selection cassette. Transgene copy number data for independent plants regenerated from the same original immature embryo suggests that multiple transgenic events arise from single immature embryos, therefore, actual efficiencies might be even higher than those reported.

CONCLUSION

We reported here a high-throughput, highly efficient and repeatable transformation system for wheat and this system has been used successfully to introduce genes of interest, for RNAi, over-expression and for CRISPR-Cas9 based genome editing.

Molecular characterization of the TaWTG1 in bread wheat (Triticum aestivum L.)

OsWTG1 (LOC_Os08g42540.1) functions as an important factor determining grain size and shape in rice. Our understanding on its ortholog in wheat, TaWTG1, is limited. Here, we identified and mapped TaWTG1 in wheat, characterized its gene and protein structures, predicted transcription factor binding sites of its promoter, and the expression patterns was also analysed bases on real-time quantitative PCR and public available microarray data. The WTG1 orthologs in barley (HvWTG1), rice (OsWTG1), Aegilops tauschii (AtWTG1), Triticum urartu (TuWTG1), Triticum turgidum (TtWTG1) and Brachypodium distachyon (BdWTG1) were also identified for comparative analyses. TaWTG1 was mapped onto the short arms of group 7 chromosomes (7AS, 7BS, and 7DS). Multiple alignments indicated that WTG1 possesses eight exons and seven introns in all of the orthologs, except for the orthologs on 7A of wild emmer and on 7D of A. tauschii (seven exons and six introns). An exon-intron junction composed of intron 2 to intron 3 and exon 2 to exon 4 was highly conserved. The protein of WTG1 exists a conserved domain (Peptidase_C65). WTG1 was mainly expressed in wheat roots, spikes and grains, in barley caryopsis and roots, and in rice anthers. Drought and heat stresses significantly regulated the expression of TaWTG1 in wheat. In barley, WTG1 was significantly down-regulated under Fusarium at late stage. In addition, significant correlations between the expression patterns of predicted transcription factors and WTG1 were also detected. Overall, the results presented here broaden our knowledge on WTG1 and will be helpful for its manipulation aiming at dissecting its function in plants.

Genome-wide association study of heading and flowering dates and construction of its prediction equation in Chinese common wheat

Heading date is one of the most important traits in wheat breeding as it affects adaptation and yield potential. A genome-wide association study (GWAS) using the 90 K iSelect SNP genotyping assay indicated that a total of 306 loci were significantly associated with heading and flowering dates in 13 environments in Chinese common wheat from the Yellow and Huai wheat region. Of these, 105 loci were significantly correlated with both heading and flowering dates and were found in clusters on chromosomes 2, 5, 6, and 7. Based on differences in distribution of the vernalization and photoperiod genes among chromosomes, arms, or block regions, 13 novel, environmentally stable genetic loci were associated with heading and flowering dates, including RAC875_c41145_189 on 1DS, RAC875_c50422_299 on 2BL, and RAC875_c48703_148 on 2DS, that accounted for more than 20% phenotypic variance explained (PVE) of the heading/flowering date in at least four environments. GWAS and t test of a combination of SNPs and vernalization and photoperiod alleles indicated that the Vrn-B1, Vrn-D1, and Ppd-D1 genes significantly affect heading and flowering dates in Chinese common wheat. Based on the association of heading and flowering dates with the vernalization and photoperiod alleles at seven loci and three significant SNPs, optimal linear regression equations were established, which show that of the seven loci, the Ppd-D1 gene plays the most important role in modulating heading and flowering dates in Chinese wheat, followed by Vrn-B1 and Vrn-D1. Additionally, three novel genetic loci (RAC875_c41145_189, Excalibur_c60164_137, and RAC875_c50422_299) also show important effect on heading and flowering dates. Therefore, Ppd-D1, Vrn-B1, Vrn-D1, and the novel genetic loci should be further investigated in terms of improving heading and flowering dates in Chinese wheat. Further quantitative analysis of an F₁₀ recombinant inbred lines population identified a major QTL that controls heading and flowering dates within the Ppd-D1 locus with PVEs of 28.4% and 34.0%, respectively; this QTL was also significantly associated with spike length, peduncle length, fertile spikelets number, cold resistance, and tiller number.

Gene editing and mutagenesis reveal inter-cultivar differences and additivity in the contribution of TaGW2 homoeologues to grain size and weight in wheat

KEY MESSAGE

CRISPR-Cas9-based genome editing and EMS mutagenesis revealed inter-cultivar differences and additivity in the contribution of TaGW2 homoeologues to grain size and weight in wheat. The TaGW2 gene homoeologues have been reported to be negative regulators of grain size (GS) and thousand grain weight (TGW) in wheat. However, the contribution of each homoeologue to trait variation among different wheat cultivars is not well documented. We used the CRISPR-Cas9 system and TILLING to mutagenize each homoeologous gene copy in cultivars Bobwhite and Paragon, respectively. Plants carrying single-copy nonsense mutations in different genomes showed different levels of GS/TGW increase, with TGW increasing by an average of 5.5% (edited lines) and 5.3% (TILLING mutants). In any combination, the double homoeologue mutants showed higher phenotypic effects than the respective single-genome mutants. The double mutants had on average 12.1% (edited) and 10.5% (TILLING) higher TGW with respect to wild-type lines. The highest increase in GS and TGW was shown for triple mutants of both cultivars, with increases in 16.3% (edited) and 20.7% (TILLING) in TGW. The additive effects of the TaGW2 homoeologues were also demonstrated by the negative correlation between the functional gene copy number and GS/TGW in Bobwhite mutants and an F2 population. The highest single-genome increases in GS and TGW in Paragon and Bobwhite were obtained by mutations in the B and D genomes, respectively. These inter-cultivar differences in the phenotypic effects between the TaGW2 gene homoeologues coincide with inter-cultivar differences in the homoeologue expression levels. These results indicate that GS/TGW variation in wheat can be modulated by the dosage of homoeologous genes with inter-cultivar differences in the magnitude of the individual homoeologue effects.

Conserved and differential transcriptional responses of peroxisome associated pathways to drought, dehydration and ABA

Plant peroxisomes are important components of cellular antioxidant networks, dealing with ROS generated by multiple metabolic pathways. Peroxisomes respond to environmental and cellular conditions by changing their size, number, and proteomic content. To investigate the role of peroxisomes in response to drought, dehydration and ABA treatment we took an evolutionary and comparative genomics approach. Colonisation of land required evolution of dehydration tolerance in the absence of subsequent anatomical adaptations. Therefore, the model bryophyte Physcomitrella patens, the model dicot Arabidopsis thaliana and wheat (Tricitcum aestivum), a globally important cereal crop were compared. Three sets of genes namely 'PTS1 genes' (a proxy for genes encoding peroxisome targeted proteins), PEX genes (involved in peroxisome biogenesis) and genes involved in plant antioxidant networks were identified in all 3 species and their expression compared under drought (dehydration) and ABA treatment. Genes encoding enzymes of β-oxidation and gluconeogenesis, antioxidant enzymes including catalase and glutathione reductase and PEX3 and PEX11 isoforms showed conserved up-regulation, and peroxisome proliferation was induced by ABA in moss. Interestingly, expression of some of these genes differed between drought sensitive and resistant genotypes of wheat in line with measured photosynthetic and biochemical differences. These results point to an underappreciated role for peroxisomes in drought response.

Comparative genome-wide mapping versus extreme pool-genotyping and development of diagnostic SNP markers linked to QTL for adult plant resistance to stripe rust in common wheat

A major stripe rust resistance QTL on chromosome 4BL was localized to a 4.5-Mb interval using comparative QTL mapping methods and validated in 276 wheat genotypes by haplotype analysis. CYMMIT-derived wheat line P10103 was previously identified to have adult plant resistance (APR) to stripe rust in the greenhouse and field. The conventional approach for QTL mapping in common wheat is laborious. Here, we performed QTL detection of APR using a combination of genome-wide scanning and extreme pool-genotyping. SNP-based genetic maps were constructed using the Wheat55 K SNP array to genotype a recombinant inbred line (RIL) population derived from the cross Mingxian 169 × P10103. Five stable QTL were detected across multiple environments. A fter comparing SNP profiles from contrasting, extreme DNA pools of RILs six putative QTL were located to approximate chromosome positions. A major QTL on chromosome 4B was identified in F_2:4 contrasting pools from cross Zhengmai 9023 × P10103. A consensus QTL (LOD = 26-40, PVE = 42-55%), named QYr.nwafu-4BL, was defined and localized to a 4.5-Mb interval flanked by SNP markers AX-110963704 and AX-110519862 in chromosome arm 4BL. Based on stripe rust response, marker genotypes, pedigree analysis and mapping data, QYr.nwafu-4BL is likely to be a new APR QTL. The applicability of the SNP-based markers flanking QYr.nwafu-4BL was validated on a diversity panel of 276 wheat lines. The additional minor QTL on chromosomes 4A, 5A, 5B and 6A enhanced the level of resistance conferred by QYr.nwafu-4BL. Marker-assisted pyramiding of QYr.nwafu-4BL and other favorable minor QTL in new wheat cultivars should improve the level of APR to stripe rust.