Mapping causal mutations by exome sequencing in a wheat TILLING population: a tall mutant case study

Targeting Induced Local Lesions in Genomes (TILLING): advances and opportunities for fast tracking crop breeding

The intensification of food production via conventional crop breeding alone is inadequate to cater for global hunger. The development of precise and expeditious high throughput reverse genetics approaches has hugely benefited modern plant breeding programs. Targeting Induced Local Lesions in Genomes (TILLING) is one such reverse genetics approach which employs chemical/physical mutagenesis to create new genetic sources and identifies superior/novel alleles. Owing to technical limitations and sectional applicability of the original TILLING protocol, it has been timely modified. Successions include: EcoTILLING, Double stranded EcoTILLING (DEcoTILLING), Self-EcoTILLING, Individualized TILLING (iTILLING), Deletion-TILLING (De-TILLING), PolyTILLING, and VeggieTILLING. This has widened its application to a variety of crops and needs. They can characterize mutations in coding as well as non-coding regions and can overcome complexities associated with the large genomes. Combining next generation sequencing tools with the existing TILLING protocols has enabled screening of huge germplasm collections and mutant populations for the target genes. In silico TILLING platforms have transformed TILLING into an exciting breeding approach. The present review outlines these multifarious TILLING modifications for precise mutation detection and their application in advance breeding programmes together with relevant case studies. Appropriate use of these protocols will open up new avenues for crop improvement in the twenty first century.

Advancements in genetic techniques and functional genomics for enhancing crop traits and agricultural sustainability

Genetic variability is essential for the development of new crop varieties with economically beneficial traits. The traits can be inherited from wild relatives or induced through mutagenesis. Novel genetic elements can then be identified and new gene functions can be predicted. In this study, forward and reverse genetics approaches were described, in addition to their applications in modern crop improvement programs and functional genomics. By using heritable phenotypes and linked genetic markers, forward genetics searches for genes by using traditional genetic mapping and allele frequency estimation. Despite recent advances in sequencing technology, omics and computation, genetic redundancy remains a major challenge in forward genetics. By analyzing close-related genes, we will be able to dissect their functional redundancy and predict possible traits and gene activity patterns. In addition to these predictions, sophisticated reverse gene editing tools can be used to verify them, including TILLING, targeted insertional mutagenesis, gene silencing, gene targeting and genome editing. By using gene knock-down, knock-up and knock-out strategies, these tools are able to detect genetic changes in cells. In addition, epigenome analysis and editing enable the development of novel traits in existing crop cultivars without affecting their genetic makeup by increasing epiallelic variants. Our understanding of gene functions and molecular dynamics of various biological phenomena has been revised by all of these findings. The study also identifies novel genetic targets in crop species to improve yields and stress tolerances through conventional and non-conventional methods. In this article, genetic techniques and functional genomics are specifically discussed and assessed for their potential in crop improvement.

Dissection of a novel major stable QTL on chromosome 7D for grain hardness and its breeding value estimation in bread wheat

Grain hardness (Gh) is important for wheat processing and end-product quality. Puroindolines polymorphism explains over 60% of Gh variation and the novel genetic factors remain to be exploited. In this study, a total of 153 quantitative trait loci (QTLs), clustered into 12 genomic intervals (C1-C12), for 13 quality-related traits were identified using a recombinant inbred line population derived from the cross of Zhongkemai138 (ZKM138) and Chuanmai44 (CM44). Among them, C7 (harboring eight QTLs for different quality-related traits) and C8 (mainly harboring QGh.cib-5D.1 for Gh) were attributed to the famous genes, Rht-D1 and Pina, respectively, indicating that the correlation of involved traits was supported by the pleotropic or linked genes. Notably, a novel major stable QTL for Gh was detected in the C12, QGh.cib-7D, with ZKM138-derived allele increasing grain hardness, which was simultaneously mapped by the BSE-Seq method. The geographic pattern and transmissibility of this locus revealed that the increasing-Gh allele is highly frequently present in 85.79% of 373 worldwide wheat varieties and presented 99.31% transmissibility in 144 ZKM138-derivatives, indicating the non-negative effect on yield performance and that its indirect passive selection has happened during the actual breeding process. Thus, the contribution of this new Gh-related locus was highlighted in consideration of improving the efficiency and accuracy of the soft/hard material selection in the molecular marker-assisted process. Further, TraesCS7D02G099400, TraesCS7D02G098000, and TraesCS7D02G099500 were initially deduced to be the most potential candidate genes of QGh.cib-7D. Collectively, this study provided valuable information of elucidating the genetic architecture of Gh for wheat quality improvement.

Genetic analyses and prediction for lodging‑related traits in a diverse Iranian hexaploid wheat collection

Lodging is one of the most important limiting environmental factors for achieving the maximum yield and quality of grains in cereals, including wheat. However, little is known about the genetic foundation underlying lodging resistance (LR) in wheat. In this study, 208 landraces and 90 cultivars were phenotyped in two cropping seasons (2018-2019 and 2019-2020) for 19 LR-related traits. A genome-wide association study (GWAS) and genomics prediction were carried out to dissect the genomic regions of LR. The number of significant marker pairs (MPs) was highest for genome B in both landraces (427,017) and cultivars (37,359). The strongest linkage disequilibrium (LD) between marker pairs was found on chromosome 4A (0.318). For stem lodging-related traits, 465, 497, and 478 marker-trait associations (MTAs) and 45 candidate genes were identified in year 1, year 2, and pooled. Gene ontology exhibited genomic region on Chr. 2B, 6B, and 7B control lodging. Most of these genes have key roles in defense response, calcium ion transmembrane transport, carbohydrate metabolic process, nitrogen compound metabolic process, and some genes harbor unknown functions that, all together may respond to lodging as a complex network. The module associated with starch and sucrose biosynthesis was highlighted. Regarding genomic prediction, the GBLUP model performed better than BRR and RRBLUP. This suggests that GBLUP would be a good tool for wheat genome selection. As a result of these findings, it has been possible to identify pivotal QTLs and genes that could be used to improve stem lodging resistance in Triticum aestivum L.

Identification and high-resolution mapping of a novel tiller number gene (tin6) by combining forward genetics screen and MutMap approach in bread wheat

Wheat (Triticum aestivum) is one of the most important food crops worldwide, providing up to 20% of the caloric intake per day. Developing high-yielding wheat cultivars with tolerance against abiotic and biotic stresses is important to keep up with the increasing human population. Tiller number is one of the major yield-related traits, directly affecting the number of grains produced per plant; however, only a small number of QTL and underlining genes have been identified for this important factor. Identification of novel genetic variation underlying contrasting traits and their precise genetic mapping in wheat is considered difficult due to the complexity and size of the genome; however, advancements in genomic resources have made efficient gene localization more possible. In this study, we report the characterization of a novel tillering number gene using a mutant identified in the forward genetic screen of an ethyl methane sulfonate (EMS)-treated population of cv. "Jagger." By crossing the low tillering mutant with the Jagger wild-type plant, we generated an F2 population and used the MutMap approach to identify a novel physical interval on 11 Mb on chromosome 2DS. Using an F₂ population of 442 gametes and polymorphic SNP markers, we were able to delineate the tin6 locus to a 2.1 Mb region containing 22 candidate genes.

Genetic approaches to exploit landraces for improvement of Triticum turgidum ssp. durum in the age of climate change

Addressing the challenges of climate change and durum wheat production is becoming an important driver for food and nutrition security in the Mediterranean area, where are located the major producing countries (Italy, Spain, France, Greece, Morocco, Algeria, Tunisia, Turkey, and Syria). One of the emergent strategies, to cope with durum wheat adaptation, is the exploration and exploitation of the existing genetic variability in landrace populations. In this context, this review aims to highlight the important role of durum wheat landraces as a useful genetic resource to improve the sustainability of Mediterranean agroecosystems, with a focus on adaptation to environmental stresses. We described the most recent molecular techniques and statistical approaches suitable for the identification of beneficial genes/alleles related to the most important traits in landraces and the development of molecular markers for marker-assisted selection. Finally, we outline the state of the art about landraces genetic diversity and signature of selection, already identified from these accessions, for adaptability to the environment.

Mapping and characterization of a novel adult-plant leaf rust resistance gene LrYang16G216 via bulked segregant analysis and conventional linkage method

A novel adult-plant leaf rust resistance gene LrYang16G216 on wheat chromosome 6BL was identified and mapped to a 0.59 cM genetic interval by BSA and conventional linkage method. Leaf rust (Puccinia triticina) is one of the most devastating fungal diseases of wheat (Triticum aestivum L.). Discovery and identification of new resistance genes is essential to develop disease-resistant cultivars. An advanced breeding line Yang16G216 was previously identified to confer adult-plant resistance (APR) to leaf rust. In this research, a recombinant inbred line (RIL) population was constructed from the cross between Yang16G216 and a highly susceptible line Yang16M6393, and genotyped with exome capture sequencing and 55 K SNP array. Through bulked segregant analysis (BSA) and genetic linkage mapping, a stable APR gene, designated as LrYang16G216, was detected and mapped to the distal region of chromosome arm 6BL with a genetic interval of 2.8 cM. For further verification, another RIL population derived from the cross between Yang16G216 and a susceptible wheat variety Yangmai 29 was analyzed using the enriched markers in the target interval, and LrYang16G216 was further narrowed to a 0.59 cM genetic interval flanked by the KASP markers Ax109403980 and Ax95083494, corresponding to the physical position 712.34-713.94 Mb in the Chinese Spring reference genome, in which twenty-six disease resistance-related genes were annotated. Based on leaf rust resistance spectrum, mapping data and physical location, LrYang16G216 was identified to be a novel and effective APR gene. The LrYang16G216 with linked markers will be useful for marker-assisted selection in wheat resistance breeding.

Genome-wide association mapping and genomic prediction of agronomical traits and breeding values in Iranian wheat under rain-fed and well-watered conditions

BACKGROUND

The markers detected by genome-wide association study (GWAS) make it possible to dissect genetic structure and diversity at many loci. This can enable a wheat breeder to reveal and used genomic loci controlling drought tolerance. This study was focused on determining the population structure of Iranian 208 wheat landraces and 90 cultivars via genotyping-by-sequencing (GBS) and also on detecting marker-trait associations (MTAs) by GWAS and genomic prediction (GS) of wheat agronomic traits for drought-tolerance breeding. GWASs were conducted using both the original phenotypes (pGWAS) and estimated breeding values (eGWAS). The bayesian ridge regression (BRR), genomic best linear unbiased prediction (gBLUP), and ridge regression-best linear unbiased prediction (rrBLUP) approaches were used to estimate breeding values and estimate prediction accuracies in genomic selection.

RESULTS

Population structure analysis using 2,174,975 SNPs revealed four genetically distinct sub-populations from wheat accessions. D-Genome harbored the lowest number of significant marker pairs and the highest linkage disequilibrium (LD), reflecting different evolutionary histories of wheat genomes. From pGWAS, BRR, gBLUP, and rrBLUP, 284, 363, 359 and 295 significant MTAs were found under normal and 195, 365, 362 and 302 under stress conditions, respectively. The gBLUP with the most similarity (80.98 and 71.28% in well-watered and rain-fed environments, correspondingly) with the pGWAS method in the terms of discovered significant SNPs, suggesting the potential of gBLUP in uncovering SNPs. Results from gene ontology revealed that 29 and 30 SNPs in the imputed dataset were located in protein-coding regions for well-watered and rain-fed conditions, respectively. gBLUP model revealed genetic effects better than other models, suggesting a suitable tool for genome selection in wheat.

CONCLUSION

We illustrate that Iranian landraces of bread wheat contain novel alleles that are adaptive to drought stress environments. gBLUP model can be helpful for fine mapping and cloning of the relevant QTLs and genes, and for carrying out trait introgression and marker-assisted selection in both normal and drought environments in wheat collections.

Fine Mapping of Stripe-Rust-Resistance Gene YrJ22 in Common Wheat by BSR-Seq and MutMap-Based Sequencing

Identification and accurate mapping of new resistance genes are essential for gene pyramiding in wheat breeding. The YrJ22 gene is a dominant stripe-rust-resistance gene located at the distal end of chromosome 2AL, which was identified in a leading Chinese-wheat variety, Jimai 22, showing high resistance to CYR32, a prevalent race of Puccinia striiformis tritici (Pst) in China. In the current study, 15 F₁ and 2273 F₂ plants derived from the cross of Jimai 22/Avocet S were used for the fine-mapping of YrJ22. The RNA-Seq of resistant and susceptible bulks of F₂ plants (designated BSR-Seq) identified 10 single-nucleotide polymorphisms (SNP) in a 12.09 Mb physical interval on chromosome 2AL. A total of 1022 EMS-induced M₃ lines of Jimai 22 were screened, to identify susceptible mutants for MutMap analysis. Four CAPS markers were developed from SNPs identified using BSR-Seq and MutMap. A linkage map for YrJ22 was constructed with 11 CAPS/STS and three SSR markers. YrJ22 was located at a 0.9 cM genetic interval flanked by markers H736 and H400, corresponding to a 340.46 kb physical region (768.7-769.0 Mb), including 13 high-confidence genes based on the Chinese Spring reference genome. TraesCS2A01G573200 is a potential candidate-gene, according to linkage and quantitative real-time PCR (qPCR) analyses. The CAPS marker H732 designed from an SNP in TraesCS2A01G573200 co-segregated with YrJ22. These results provide a useful stripe-rust-resistance gene and molecular markers for marker-assisted selection in wheat breeding and for further cloning of the gene.

Utility of Triti-Map for bulk-segregated mapping of causal genes and regulatory elements in Triticeae

Triticeae species, including wheat, barley, and rye, are critical for global food security. Mapping agronomically important genes is crucial for elucidating molecular mechanisms and improving crops. However, Triticeae includes many wild relatives with desirable agronomic traits, and frequent introgressions occurred during Triticeae evolution and domestication. Thus, Triticeae genomes are generally large and complex, making the localization of genes or functional elements that control agronomic traits challenging. Here, we developed Triti-Map, which contains a suite of user-friendly computational packages specifically designed and optimized to overcome the obstacles of gene mapping in Triticeae, as well as a web interface integrating multi-omics data from Triticeae for the efficient mining of genes or functional elements that control particular traits. The Triti-Map pipeline accepts both DNA and RNA bulk-segregated sequencing data as well as traditional QTL data as inputs for locating genes and elucidating their functions. We illustrate the usage of Triti-Map with a combination of bulk-segregated ChIP-seq data to detect a wheat disease-resistance gene with its promoter sequence that is absent from the reference genome and clarify its evolutionary process. We hope that Triti-Map will facilitate gene isolation and accelerate Triticeae breeding.

Genetic Analysis of Adult Plant Resistance to Stripe Rust in Common Wheat Cultivar "Pascal"

Wheat stripe rust is an important foliar disease that affects the wheat yield globally. Breeding for resistant wheat varieties is one of the most economically and environmentally effective ways to control this disease. The common wheat (Triticum aestivum L.) cultivar "Pascal" exhibited susceptibility to stripe rust at the seedling stage but it showed high resistance to stripe rust at the adult plant stage over 20 years in Gansu, a hotspot of the disease in northwestern China. To understand the genetic mechanism of stripe rust resistance in this cultivar, a 55K SNP array was used to analyze the two parents and the 220 recombinant inbred lines (RILs) derived from the cross of "Huixianhong" × "Pascal." We detected three new stripe rust adult plant resistance (APR) quantitative trait locus (QTL) contributed by Pascal, viz. QYr.gaas-1AL, QYr.gaas-3DL, and QYr.gaas-5AS, using the inclusive composite interval mapping method. They were flanked by SNP markers AX-111218361-AX-110577861, AX-111460455-AX-108798599, and AX-111523523-AX-110028503, respectively, and explained the phenotypic variation ranging from 11.0 to 23.1%. Bulked segregant exome capture sequencing (BSE-Seq) was used for fine mapping of QYr.gaas-1AL and selection of candidate genes, TraesCS1A02G313700, TraesCS1A02G313800, and TraesCS1A02G314900 for QYr.gaas-1AL. KASP markers BSE-1A-12 and HXPA-3D for QYr.gaas-1AL and QYr.gaas-3DL were developed for breeders to develop durable stripe rust-resistant wheat varieties.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches

Background

Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking.

Results

We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs.

Conclusion

We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

Reverse genetic approaches for breeding nutrient-rich and climate-resilient cereal and food legume crops

In the last decade, advancements in genomics tools and techniques have led to the discovery of many genes. Most of these genes still need to be characterized for their associated function and therefore, such genes remain underutilized for breeding the next generation of improved crop varieties. The recent developments in different reverse genetic approaches have made it possible to identify the function of genes controlling nutritional, biochemical, and metabolic traits imparting drought, heat, cold, salinity tolerance as well as diseases and insect-pests. This article focuses on reviewing the current status and prospects of using reverse genetic approaches to breed nutrient-rich and climate resilient cereal and food legume crops.

A Multi-Omics Approach for Rapid Identification of Large Genomic Lesions at the Wheat Dense Spike (wds) Locus

Optimal spike architecture provides a favorable structure for grain development and yield improvement. However, the number of genes cloned to underlie wheat spike architecture is extremely limited. Here, we obtained a wheat dense spike mutant (wds) induced by 60Co treatment of a common wheat landrace Huangfangzhu that exhibited significantly reduced spike and grain lengths. The shortened spike length was caused by longitudinal reduction in number and length of rachis cells. We adopted a multi-omics approach to identify the genomic locus underlying the wds mutant. We performed Exome Capture Sequencing (ECS) and identified two large deletion segments, named 6BL.1 at 334.8∼424.3 Mb and 6BL.2, 579.4∼717.8 Mb in the wds mutant. RNA-seq analysis confirmed that genes located in these regions lost their RNA expression. We then found that the 6BL.2 locus was overlapping with a known spike length QTL, qSL6B.2. Totally, 499 genes were located within the deleted region and two of them were found to be positively correlated with long spike accessions but not the ones with short spike. One of them, TraesCS6B01G334600, a well-matched homolog of the rice OsBUL1 gene that works in the Brassinosteroids (BR) pathway, was identified to be involved in cell size and number regulation. Further transcriptome analysis of young spikes showed that hormone-related genes were enriched among differentially expressed genes, supporting TraesCS6B01G334600 as a candidate gene. Our work provides a strategy to rapid locate genetic loci with large genomic lesions in wheat and useful resources for future wheat study.

Effects of the different waxy proteins on starch biosynthesis, starch physicochemical properties and Chinese noodle quality in wheat

Noodles are an important food in Asia. Wheat starch is the most important component in Chinese noodles. Loss of the waxy genes leads to lower activity of starch synthesis enzymes and decreased amylose content that further affects starch properties and noodle quality. To study the effects of different waxy (Wx) protein subunits on starch biosynthesis and processing quality, the high-yielding wheat cultivar Jimai 22 was treated with the mutagen ethyl methane sulfonate (EMS) to produce a population of Wx lines and chosen 7 Wx protein combinations. The amylose content increased but swelling power decreased as the number of Wx proteins increased. Both GBSS activity and gene expression were the lowest for the waxy mutant, followed by the mutants with 1 Wx protein. The combinations of these mutant alleles lead to reductions in both RNA expression and protein levels. Noodles made from materials with 2 Wx protein subunits had the highest score, which agreed with peak viscosity. The influence of the Wx-B1 protein on amylose synthesis and noodle quality was the highest, whereas the influence of Wx-A1 protein was the lowest. Mutants with lower amylose content caused by the absence of 1 subunit, especially the Wx-B1 subunit, had superior noodle quality. Additionally, the identified mutant lines can be used as intermediate materials to improve wheat quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01292-x.

Whole-exome sequencing of selected bread wheat recombinant inbred lines as a useful resource for allele mining and bulked segregant analysis

Although wheat (Triticum aestivum L.) is the main staple crop in the world and a major source of carbohydrates and proteins, functional genomics and allele mining are still big challenges. Given the advances in next-generation sequencing (NGS) technologies, the identification of causal variants associated with a target phenotype has become feasible. For these reasons, here, by combining sequence capture and target-enrichment methods with high-throughput NGS re-sequencing, we were able to scan at exome-wide level 46 randomly selected bread wheat individuals from a recombinant inbred line population and to identify and classify a large number of single nucleotide polymorphisms (SNPs). For technical validation of results, eight randomly selected SNPs were converted into Kompetitive Allele-Specific PCR (KASP) markers. This resource was established as an accessible and reusable molecular toolkit for allele data mining. The dataset we are making available could be exploited for novel studies on bread wheat genetics and as a foundation for starting breeding programs aimed at improving different key agronomic traits.

Whole Exome Sequencing-Based Identification of a Novel Gene Involved in Root Hair Development in Barley (Hordeum vulgare L.)

Root hairs play a crucial role in anchoring plants in soil, interaction with microorganisms and nutrient uptake from the rhizosphere. In contrast to Arabidopsis, there is a limited knowledge of root hair morphogenesis in monocots, including barley (Hordeum vulgare L.). We have isolated barley mutant rhp1.e with an abnormal root hair phenotype after chemical mutagenesis of spring cultivar ‘Sebastian’. The development of root hairs was initiated in the mutant but inhibited at the very early stage of tip growth. The length of root hairs reached only 3% of the length of parent cultivar. Using a whole exome sequencing (WES) approach, we identified G1674A mutation in the HORVU1Hr1G077230 gene, located on chromosome 1HL and encoding a cellulose synthase-like C1 protein (HvCSLC1) that might be involved in the xyloglucan (XyG) synthesis in root hairs. The identified mutation led to the retention of the second intron and premature termination of the HvCSLC1 protein. The mutation co-segregated with the abnormal root hair phenotype in the F₂ progeny of rhp1.e mutant and its wild-type parent. Additionally, different substitutions in HORVU1Hr1G077230 were found in four other allelic mutants with the same root hair phenotype. Here, we discuss the putative role of HvCSLC1 protein in root hair tube elongation in barley.

Delaying or delivering: identification of novel NAM-1 alleles that delay senescence to extend wheat grain fill duration

Senescence is a complex trait under genetic and environmental control, in which resources are remobilized from vegetative tissue into grain. Delayed senescence, or 'staygreen' traits, can confer stress tolerance, with extended photosynthetic activity hypothetically sustaining grain filling. The genetics of senescence regulation are largely unknown, with senescence variation often correlated with phenological traits. Here, we confirm staygreen phenotypes of two Triticum aestivum cv. Paragon ethyl methane sulfonate mutants previously identified during a forward genetic screen and selected for their agronomic performance, similar phenology, and differential senescence phenotypes. Grain filling experiments confirmed a positive relationship between onset of senescence and grain fill duration, reporting an associated ~14% increase in final dry grain weight for one mutant (P<0.05). Recombinant inbred line (RIL) populations segregating for the timing of senescence were developed for trait mapping purposes and phenotyped over multiple years under field conditions. Quantification and comparison of senescence metrics aided RIL selection, facilitating exome capture-enabled bulk segregant analysis (BSA). Using BSA we mapped our two staygreen traits to two independent, dominant, loci of 4.8 and 16.7 Mb in size encompassing 56 and 142 genes, respectively. Combining association analysis with variant effect prediction, we identified single nucleotide polymorphisms encoding self-validating mutations located in NAM-1 homoeologues, which we propose as gene candidates.

Mapping and Characterization of a Wheat Stem Rust Resistance Gene in Durum Wheat "Kronos"

Wheat stem (or black) rust is one of the most devastating fungal diseases, threatening global wheat production. Identification, mapping, and deployment of effective resistance genes are critical to addressing this challenge. In this study, we mapped and characterized one stem rust resistance (Sr) gene from the tetraploid durum wheat variety Kronos (temporary designation SrKN). This gene was mapped on the long arm of chromosome 2B and confers resistance to multiple virulent Pgt races, such as TRTTF and BCCBC. Using a large mapping population (3,366 gametes), we mapped SrKN within a 0.29 cM region flanked by the sequenced-based markers pku4856F2R2 and pku4917F3R3, which corresponds to 5.6- and 7.2-Mb regions in the Svevo and Chinese Spring reference genomes, respectively. Both regions include a cluster of nucleotide binding leucine-repeat (NLR) genes that likely includes the candidate gene. An allelism test failed to detect recombination between SrKN and the previously mapped Sr9e gene. This result, together with the similar seedling resistance responses and resistance profiles, suggested that SrKN and Sr9e may represent the same gene. We introgressed SrKN into common wheat and developed completely linked markers to accelerate its deployment in the wheat breeding programs. SrKN can be a valuable component of transgenic cassettes or gene pyramids that includes multiple resistance genes to control this devastating disease.

The Triple Jags of Dietary Fibers in Cereals: How Biotechnology Is for High

Cereals represent an important source of beneficial compounds for human health, such as macro- and micronutrients, vitamins, and bioactive molecules. Generally, the consumption of whole-grain products is associated with significant health benefits, due to the elevated amount of dietary fiber (DF). However, the consumption of whole-grain foods is still modest compared to more refined products. In this sense, it is worth focusing on the increase of DF fractions inside the inner compartment of the seed, the endosperm, which represents the main part of the derived flour. The main components of the grain fiber are arabinoxylan (AX), β-glucan (βG), and resistant starch (RS). These three components are differently distributed in grains, however, all of them are represented in the endosperm. AX and βG, classified as non-starch polysaccharides (NSP), are in cell walls, whereas, RS is in the endosperm, being a starch fraction. As the chemical structure of DFs influences their digestibility, the identification of key actors involved in their metabolism can pave the way to improve their function in human health. Here, we reviewed the main achievements of plant biotechnologies in DFs manipulation in cereals, highlighting new genetic targets to be exploited, and main issues to face to increase the potential of cereals in fighting malnutrition.

Identification of a major and stable QTL on chromosome 5A confers spike length in wheat (Triticum aestivum L.)

Spike length (SL) is the key determinant of plant architecture and yield potential. In this study, 193 recombinant inbred lines (RILs) derived from a cross between 13F10 and Chuanmai 42 (CM42) were evaluated for spike length in six environments. Sixty RILs consisting of 30 high and 30 low SLs were genotyped using the bulked segregant analysis exome sequencing (BSE-Seq) analysis for preliminary quantitative trait locus (QTL) mapping. A 6.69 Mb (518.43-525.12 Mb) region on chromosome 5AL was found to have a significant effect on the SL trait. Fifteen competitive allele-specific PCR (KASP) markers were successfully converted from the single nucleotide polymorphisms (SNPs) in the SL target region. Combined with four novel simple sequence repeat (SSR) markers, a genetic linkage map spanning 21.159 cM was constructed. The mapping result confirmed the identity of a major and stable QTL named QSl.cib-5A in the targeted region that explained 7.88-26.60% of the phenotypic variation in SL. QSl.cib-5A was narrowed to a region of 4.84 cM interval corresponding to a 4.67 Mb (516.60-521.27 Mb) physical region in the Chinese Spring RefSeq v2.0 containing 17 high-confidence genes with 25 transcripts. In addition, this QTL exhibited pleiotropic effects on spikelet density (SD), with the phenotypic variances proportion ranging from 11.34 to 19.92%. This study provides a foundational step for cloning the QSl.cib-5A, which is involved in the regulation of spike morphology in common wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01249-6.

Random mutagenesis in vegetatively propagated crops: opportunities, challenges and genome editing prospects

In order to meet the growing human food and nutrition demand a perpetual process of crop improvement is idealized. It has seen changing trends and varying concepts throughout human history; from simple selection to complex gene-editing. Among these techniques, random mutagenesis has been shown to be a promising technology to achieve desirable genetic gain with less time and minimal efforts. Over the decade, several hundred varieties have been released through random mutagenesis, but the production is falling behind the demand. Several food crops like banana, potato, cassava, sweet potato, apple, citrus, and others are vegetatively propagated. Since such crops are not propagated through seed, genetic improvement through classical breeding is impractical for them. Besides, in the case of polyploids, accomplishment of allelic homozygosity requires a considerable land area, extensive fieldwork with huge manpower, and hefty funding for an extended period of time. Apart from induction, mapping of induced genes to facilitate the knowledge of biological processes has been performed only in a few selected facultative vegetative crops like banana and cassava which can form a segregating population. During the last few decades, there has been a shift in the techniques used for crop improvement. With the introduction of the robust technologies like meganucleases, zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats (CRISPR) more and more crops are being subjected to gene editing. However, more work needs to be done in case of vegetatively propagated crops.

The Sequencing-Based Mapping Method for Effectively Cloning Plant Mutated Genes

A forward genetic approach is a powerful tool for identifying the genes underlying the phenotypes of interest. However, the conventional map-based cloning method is lengthy, requires a large mapping population and confirmation of many candidate genes in a broad genetic region to clone the causal variant. The whole-genome sequencing method clones the variants with a certain failure probability for multiple reasons, especially for heterozygotes, and could not be used to clone the mutation of epigenetic modifications. Here, we applied the highly complementary characteristics of these two methods and developed a sequencing-based mapping method (SBM) for identifying the location of plant variants effectively with a small population and low cost, which is very user-friendly for most popular laboratories. This method used the whole-genome sequencing data of two pooled populations to screen out enough markers. These markers were used to identify and narrow the candidate region by analyzing the marker-indexes and recombinants. Finally, the possible mutational sites were identified using the whole-genome sequencing data and verified in individual mutants. To elaborate the new method, we displayed the cloned processes in one Arabidopsis heterozygous mutant and two rice homozygous mutants. Thus, the sequencing-based mapping method could clone effectively different types of plant mutations and was a powerful tool for studying the functions of plant genes in the species with known genomic sequences.

A TILLING by sequencing approach to identify induced mutations in sunflower genes

The Targeting Induced Local Lesions in Genomes (TILLING) technology is a reverse genetic strategy broadly applicable to every kind of genome and represents an attractive tool for functional genomic and agronomic applications. It consists of chemical random mutagenesis followed by high-throughput screening of point mutations in targeted genomic regions. Although multiple methods for mutation discovery in amplicons have been described, next-generation sequencing (NGS) is the tool of choice for mutation detection because it quickly allows for the analysis of a large number of amplicons. The aim of the present work was to screen a previously generated sunflower TILLING population and identify alterations in genes involved in several important and complex physiological processes. Twenty-one candidate sunflower genes were chosen as targets for the screening. The TILLING by sequencing strategy allowed us to identify multiple mutations in selected genes and we subsequently validated 16 mutations in 11 different genes through Sanger sequencing. In addition to addressing challenges posed by outcrossing, our detection and validation of mutations in multiple regulatory loci highlights the importance of this sunflower population as a genetic resource.

Technologies enabling rapid crop improvements for sustainable agriculture: example pennycress (Thlaspi arvense L.)

Growing concerns over food insecurity and ecosystems health related to population growth and climate change have challenged scientists to develop new crops, employing revolutionary technologies in combination with traditional methods. In this review, we discuss the domestication of the oilseed-producing cover crop pennycress, which along with the development of other new crops and improvements to farming practices can provide sustainable solutions to address malnutrition and environmental impacts of production agriculture. We highlight some of the new technologies such as bioinformatics-enabled next-generation sequencing and CRISPR genome editing in combination with traditional mutation breeding that has accelerated pennycress development as a new crop and a potential model system. Furthermore, we provide a brief overview of the technologies that can be integrated for improving pennycress and other crops and the status of pennycress development using these technologies.

VAP-RELATED SUPPRESSORS OF TOO MANY MOUTHS (VST) family proteins are regulators of root system architecture

Root system architecture (RSA) is a key factor in the efficiency of nutrient capture and water uptake in plants. Understanding the genetic control of RSA will be useful in minimizing fertilizer and water usage in agricultural cropping systems. Using a hydroponic screen and a gel-based imaging system, we identified a rice (Oryza sativa) gene, VAP-RELATED SUPPRESSOR OF TOO MANY MOUTHS1 (OsVST1), which plays a key role in controlling RSA. This gene encodes a homolog of the VAP-RELATED SUPPRESSORS OF TOO MANY MOUTHS (VST) proteins in Arabidopsis (Arabidopsis thaliana), which promote signaling in stomata by mediating plasma membrane-endoplasmic reticulum contacts. OsVST1 mutants have shorter primary roots, decreased root meristem size, and a more compact RSA. We show that the Arabidopsis VST triple mutants have similar phenotypes, with reduced primary root growth and smaller root meristems. Expression of OsVST1 largely complements the short root length and reduced plant height in the Arabidopsis triple mutant, supporting conservation of function between rice and Arabidopsis VST proteins. In a field trial, mutations in OsVST1 did not adversely affect grain yield, suggesting that modulation of this gene could be used as a way to optimize RSA without an inherent yield penalty.

An ethyl methanesulfonate-induced neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice

Spontaneous mutants are mainly obtained from tissue culture or natural occurrences in plants. The traditional strategy for identifying spontaneously mutated genes is to continuously backcross these mutants to another variety and develop a near-isogenic F₂ population for map-based cloning or bulked segregant analysis. However, this strategy is time-consuming. Here, we have developed a new method to efficiently accelerate the identification process. The chemical mutagen ethyl methanesulfonate was first used to treat the wild type of the spontaneous mutants to induce thousands of neutral mutations. An induced individual without any statistically significant phenotypic changes which was compared with the wild type was chosen as the neutral mutant. The spontaneous mutant was then crossed with the neutral mutant to develop a pseudo-near-isogenic F₂ population in which only the induced neutral mutations and the causal mutation were segregated in the genome. This population ensures that the variation of the mutated trait is controlled only by the spontaneously mutated gene. Finally, after sequencing the neutral mutant and the mutant-type DNA pool of the F₂ population the spontaneous mutation will be identified quickly by bioinformatics analysis. Using this method, two spontaneously mutated genes were identified successfully. Therefore, the neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice, and its value in other plants is discussed.

Next Generation Sequencing Based Forward Genetic Approaches for Identification and Mapping of Causal Mutations in Crop Plants: A Comprehensive Review

The recent advancements in forward genetics have expanded the applications of mutation techniques in advanced genetics and genomics, ahead of direct use in breeding programs. The advent of next-generation sequencing (NGS) has enabled easy identification and mapping of causal mutations within a short period and at relatively low cost. Identifying the genetic mutations and genes that underlie phenotypic changes is essential for understanding a wide variety of biological functions. To accelerate the mutation mapping for crop improvement, several high-throughput and novel NGS based forward genetic approaches have been developed and applied in various crops. These techniques are highly efficient in crop plants, as it is relatively easy to grow and screen thousands of individuals. These approaches have improved the resolution in quantitative trait loci (QTL) position/point mutations and assisted in determining the functional causative variations in genes. To be successful in the interpretation of NGS data, bioinformatics computational methods are critical elements in delivering accurate assembly, alignment, and variant detection. Numerous bioinformatics tools/pipelines have been developed for such analysis. This article intends to review the recent advances in NGS based forward genetic approaches to identify and map the causal mutations in the crop genomes. The article also highlights the available bioinformatics tools/pipelines for reducing the complexity of NGS data and delivering the concluding outcomes.

Genetic architecture underpinning yield component traits in wheat

Genetic atlas, reliable QTL and candidate genes of yield component traits in wheat were figured out, laying concrete foundations for map-based gene cloning and dissection of regulatory mechanisms underlying yield. Mining genetic loci for yield is challenging due to the polygenic nature, large influence of environment and complex relationship among yield component traits (YCT). Many genetic loci related to wheat yield have been identified, but its genetic architecture and key genetic loci for selection are largely unknown. Wheat yield potential can be determined by three YCT, thousand kernel weight, kernel number per spike and spike number. Here, we summarized the genetic loci underpinning YCT from QTL mapping, association analysis and homology-based gene cloning. The major loci determining yield-associated agronomic traits, such as flowering time and plant height, were also included in comparative analyses with those for YCT. We integrated yield-related genetic loci onto chromosomes based on their physical locations. To identify the major stable loci for YCT, 58 QTL-rich clusters (QRC) were defined based on their distribution on chromosomes. Candidate genes in each QRC were predicted according to gene annotation of the wheat reference genome and previous information on validation of those genes in other species. Finally, a technological route was proposed to take full advantage of the resultant resources for gene cloning, molecular marker-assisted breeding and dissection of molecular regulatory mechanisms underlying wheat yield.

Longin CF, Ganal M, Himmelbach A, Reif JC, Schulthess AW. Exome association analysis sheds light onto leaf rust (Puccinia triticina) resistance genes currently used in wheat breeding (Triticum aestivum L.)

Resistance breeding is crucial for a sustainable control of leaf rust (Puccinia triticina) in wheat (Triticum aestivum L.) while directly targeting functional variants is the Holy Grail for efficient marker-assisted selection and map-based cloning. We assessed the limits and prospects of exome association analysis for severity of leaf rust in a large hybrid wheat population of 1574 single-crosses plus their 133 parents. After imputation and quality control, exome sequencing revealed 202 875 single-nucleotide polymorphisms (SNPs) covering 19.7% of the high-confidence annotated gene space. We performed intensive data mining and found significant associations for 2171 SNPs corresponding to 50 different loci. Some of these associations mapped in the proximity of the already known resistance genes Lr21, Lr34-B, Lr1 and Lr10, while other associated genomic regions, such as those on chromosomes 1A and 3D, harboured several annotated genes putatively involved in resistance. Validation with an independent population helped to narrow down the list of putative resistance genes that should be targeted by fine-mapping. We expect that the proposed strategy of intensive data mining coupled with validation will significantly influence research in plant genetics and breeding.

Exome sequencing of bulked segregants identified a novel TaMKK3-A allele linked to the wheat ERA8 ABA-hypersensitive germination phenotype

Using bulked segregant analysis of exome sequence, we fine-mapped the ABA-hypersensitive mutant ERA8 in a wheat backcross population to the TaMKK3-A locus of chromosome 4A. Preharvest sprouting (PHS) is the germination of mature grain on the mother plant when it rains before harvest. The ENHANCED RESPONSE TO ABA8 (ERA8) mutant increases seed dormancy and, consequently, PHS tolerance in soft white wheat 'Zak.' ERA8 was mapped to chromosome 4A in a Zak/'ZakERA8' backcross population using bulked segregant analysis of exome sequenced DNA (BSA-exome-seq). ERA8 was fine-mapped relative to mutagen-induced SNPs to a 4.6 Mb region containing 70 genes. In the backcross population, the ERA8 ABA-hypersensitive phenotype was strongly linked to a missense mutation in TaMKK3-A-G1093A (LOD 16.5), a gene associated with natural PHS tolerance in barley and wheat. The map position of ERA8 was confirmed in an 'Otis'/ZakERA8 but not in a 'Louise'/ZakERA8 mapping population. This is likely because Otis carries the same natural PHS susceptible MKK3-A-A660^S allele as Zak, whereas Louise carries the PHS-tolerant MKK3-A-C660^R allele. Thus, the variation for grain dormancy and PHS tolerance in the Louise/ZakERA8 population likely resulted from segregation of other loci rather than segregation for PHS tolerance at the MKK3 locus. This inadvertent complementation test suggests that the MKK3-A-G1093A mutation causes the ERA8 phenotype. Moreover, MKK3 was a known ABA signaling gene in the 70-gene 4.6 Mb ERA8 interval. None of these 70 genes showed the differential regulation in wild-type Zak versus ERA8 expected of a promoter mutation. Thus, the working model is that the ERA8 phenotype results from the MKK3-A-G1093A mutation.

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.

Genetic dissection of a major QTL for kernel weight spanning the Rht-B1 locus in bread wheat

Genetic dissection uncovered a major QTL QTKW.caas-4BS corresponding with a 483 kb deletion that included genes ZnF, EamA and Rht-B1. This deletion was associated with increased grain weight and semi-dwarf phenotype. Previous studies identified quantitative trait loci (QTL) for thousand kernel weight (TKW) in the region spanning the Rht-B1 locus in wheat (Triticum aestivum L.). We recently mapped a major QTL QTKW.caas-4BS for TKW spanning the Rht-B1 locus in a recombinant inbred line (RIL) population derived from Doumai/Shi 4185 using the wheat 90K array. The allele from Doumai at QTKW.caas-4BS significantly increased TKW and kernel number per spike, and conferred semi-dwarf trait, which was beneficial to improve grain yield without a penalty to lodging. To further dissect QTKW.caas-4BS, we firstly re-investigated the genotypes and phenotypes of the RILs and confirmed the QTL using cleaved amplified polymorphic sequence (CAPS) markers developed from flanking SNP markers IWA102 and IWB54814. The target sequences of the CAPS markers were used as queries to BLAST the wheat reference genome RefSeq v1.0 and hit an approximate 10.4 Mb genomic region. Based on genomic mining and SNP loci from the wheat 660K SNP array in the above genomic region, we developed eight new markers and narrowed QTKW.caas-4BS to a genetic interval of 1.5 cM. A 483 kb deletion in Doumai corresponded with QTKW.caas-4BS genetically, including three genes ZnF, EamA and Rht-B1. The other 15 genes with either differential expressions and/or sequence variations between parents were also potential candidate genes for QTKW.caas-4BS. The findings not only provide a toolkit for marker-assisted selection of QTKW.caas-4BS but also defined candidate genes for further functional analysis.

Genetic Screening for Mutants with Altered Seminal Root Numbers in Hexaploid Wheat Using a High-Throughput Root Phenotyping Platform

Roots are the main channel for water and nutrient uptake in plants. Optimization of root architecture provides a viable strategy to improve nutrient and water uptake efficiency and maintain crop productivity under water-limiting and nutrient-poor conditions. We know little, however, about the genetic control of root development in wheat, a crop supplying 20% of global calorie and protein intake. To improve our understanding of the genetic control of seminal root development in wheat, we conducted a high-throughput screen for variation in seminal root number using an exome-sequenced mutant population derived from the hexaploid wheat cultivar Cadenza. The screen identified seven independent mutants with homozygous and stably altered seminal root number phenotypes. One mutant, Cadenza0900, displays a recessive extra seminal root number phenotype, while six mutants (Cadenza0062, Cadenza0369, Cadenza0393, Cadenza0465, Cadenza0818 and Cadenza1273) show lower seminal root number phenotypes most likely originating from defects in the formation and activation of seminal root primordia. Segregation analysis in F2 populations suggest that the phenotype of Cadenza0900 is controlled by multiple loci whereas the Cadenza0062 phenotype fits a 3:1 mutant:wild-type segregation ratio characteristic of dominant single gene action. This work highlights the potential to use the sequenced wheat mutant population as a forward genetic resource to uncover novel variation in agronomic traits, such as seminal root architecture.

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.

Development of model web-server for crop variety identification using throughput SNP genotyping data

Crop varieties or genotypes of a given species are pivotal for agricultural production and ownership, management and improvement of their germplasm is a great challenge. Its morphological identification requires time, cost and descriptors are often compromised statistically due to phenotypic plasticity. Development of DNA based signature of varieties can overcome these limitations. There is a global need to implement world trade organization (WTO) and intellectual property rights (IPR) guidelines of Plant Breeders Rights (PBR) where DUS (distinctness, uniformity and stability) testing can be supplemented by DNA profile. Universalization and minimization of SNP number without compromising identification accuracy is the major challenge in development of varietal profile by rapid genotype assay. Besides this, there is no server-based approach reducing computational skill with global accessibility of referral phenotypic and genotypic data. We report world’s first model web server for crop variety identification using >350 Indian wheat varieties and Axiom 35 K SNP chip data. Standard filtering and linkage disequilibrium approach were used to develop varietal signature in Linux using HTML, Java, PHP and MySQL with provision of QR code generator to facilitate bar-coding. Phylogenetic tree constructed by selected SNPs confirms six major trait based clusters of varieties and their pedigree. Our user friendly server based tool, VISTa (Variety Identification System of Triticum aestivum) (http://webtom.cabgrid.res.in/vista) can be used in DUS testing having dispute resolution of sovereignty and access benefit sharing (ABS) issues. This model approach can be used in other crops with pan-global level management of crop germplasm in endeavour of crop productivity.

Applying the latest advances in genomics and phenomics for trait discovery in polyploid wheat

Improving traits in wheat has historically been challenging due to its large and polyploid genome, limited genetic diversity and in-field phenotyping constraints. However, within recent years many of these barriers have been lowered. The availability of a chromosome-level assembly of the wheat genome now facilitates a step-change in wheat genetics and provides a common platform for resources, including variation data, gene expression data and genetic markers. The development of sequenced mutant populations and gene-editing techniques now enables the rapid assessment of gene function in wheat directly. The ability to alter gene function in a targeted manner will unmask the effects of homoeolog redundancy and allow the hidden potential of this polyploid genome to be discovered. New techniques to identify and exploit the genetic diversity within wheat wild relatives now enable wheat breeders to take advantage of these additional sources of variation to address challenges facing food production. Finally, advances in phenomics have unlocked rapid screening of populations for many traits of interest both in greenhouses and in the field. Looking forwards, integrating diverse data types, including genomic, epigenetic and phenomics data, will take advantage of big data approaches including machine learning to understand trait biology in wheat in unprecedented detail.

Phenotypic and transcriptomic characterization of a wheat tall mutant carrying an induced mutation in the C-terminal PFYRE motif of RHT-B1b

BACKGROUND

As central regulators of the gibberellic acid (GA) signaling pathway in plants, DELLA proteins function as growth repressors and affect diverse biological processes. The wheat RHT-B1b and RHT-D1b semi-dwarfing alleles, which encode GA-insensitive DELLA proteins, have been widely adopted in modern wheat varieties to improve lodging tolerance and harvest index. However, the molecular mechanisms by which DELLA modulates these responses in wheat remain largely unknown.

RESULTS

We identified a tall tetraploid wheat mutant line carrying an induced missense mutation (E529K) in the PFYRE motif of RHT-B1b that partially suppressed the semi-dwarf phenotype. The height-increasing effect of RHT-B1bE529K relative to RHT-B1b (19 cm or 21% increase) was significantly smaller than the effect of RHT-B1a (33 cm or 34% increase) relative to RHT-B1b in the same field experiment. The RHT-B1bE529K mutation was also associated with length increases in coleoptiles, seedling shoots, and stem internodes relative to the RHT-B1b allele. We detected no significant differences in germination rate, seedling root length, tiller number, flag leaf size, spike length, or yield components. Using RNA-seq, we compared gene expression profiles of plants encoding RHT-B1b and RHT-B1bE529K in coleoptile, first leaf, and elongating peduncles. We detected limited overlap among tissues of the genes differentially regulated by the two genotypes, and more genes upregulated (77%) than downregulated (23%) in RHT-B1bE529K relative to RHT-B1b. These results suggest that the wheat DELLA protein affects the transcriptome in a tissue-specific manner and that the mutation mainly eliminates or reduces repression functions of the RHT-B1 protein. Our study identified distinct sets of potential DELLA direct or indirect target genes involved in cell wall and carbohydrate metabolisms, cell cycle/division, and hormone pathways.

CONCLUSIONS

We identified the hypomorphic RHT-B1bE529K allele that confers an intermediate plant height and coleoptile elongation. This allele can be useful in rain-fed wheat breeding programs where the strong reduction in height and biomass associated with RHT-B1b has detrimental effects. Transcriptomic characterization of different tissues from the plants encoding RHT-B1bE529K and RHT-B1b provided valuable information for identifying DELLA downstream GA response genes in wheat.

Identification and characterization of Rht25, a locus on chromosome arm 6AS affecting wheat plant height, heading time, and spike development

This study identified Rht25, a new plant height locus on wheat chromosome arm 6AS, and characterized its pleiotropic effects on important agronomic traits. Understanding genes regulating wheat plant height is important to optimize harvest index and maximize grain yield. In modern wheat varieties grown under high-input conditions, the gibberellin-insensitive semi-dwarfing alleles Rht-B1b and Rht-D1b have been used extensively to confer lodging tolerance and improve harvest index. However, negative pleiotropic effects of these alleles (e.g., poor seedling emergence and reduced biomass) can cause yield losses in hot and dry environments. As part of current efforts to diversify the dwarfing alleles used in wheat breeding, we identified a quantitative trait locus (QHt.ucw-6AS) affecting plant height in the proximal region of chromosome arm 6AS (< 0.4 cM from the centromere). Using a large segregating population (~ 2800 gametes) and extensive progeny tests (70-93 plants per recombinant family), we mapped QHt.ucw-6AS as a Mendelian locus to a 0.2 cM interval (144.0-148.3 Mb, IWGSC Ref Seq v1.0) and show that it is different from Rht18. QHt.ucw-6AS is officially designated as Rht25, with Rht25a representing the height-increasing allele and Rht25b the dwarfing allele. The average dwarfing effect of Rht25b was found to be approximately half of the effect observed for Rht-B1b and Rht-D1b, and the effect is greater in the presence of the height-increasing Rht-B1a and Rht-D1a alleles than in the presence of the dwarfing alleles. Rht25b is gibberellin-sensitive and shows significant pleiotropic effects on coleoptile length, heading date, spike length, spikelet number, spikelet density, and grain weight. Rht25 represents a new alternative dwarfing locus that should be evaluated for its potential to improve wheat yield in different environments.

Identification of induced mutations in hexaploid wheat genome using exome capture assay

Wheat is a staple food crop of many countries. Improving resilience to biotic and abiotic stresses remain key breeding targets. Among these, rust diseases are the most detrimental in terms of depressing wheat production. In the present study, chemical mutagenesis was used to induce mutations in the wheat variety NN-Gandum-1. This cultivar is moderately resistant to leaf and yellow rust. The aim of mutagenesis was to improve resistance to the disease as well as to study function of genes conferring resistance to the disease. In the present investigation, a 0.8% EMS dose was found optimum for supporting 45-55% germination of NN-Gandum-1. A total of 3,634 M2 fertile plants were produced from each of the M1 plant. Out of these, 33 (0.91%) and 20 plants (0.55%) showed absolute resistance to leaf and yellow rust, respectively. While 126 (3.46%) and 127 plants (3.49%) exhibited high susceptibility to the leaf and yellow rust, respectively. In the M4 generation, a total of 11 M4 lines (nine absolute resistant and two highly susceptible) and one wild type were selected for NGS-based exome capture assay. A total of 104,779 SNPs were identified that were randomly distributed throughout the wheat sub genomes (A, B and D). Induced mutations in intronic sequences predominated. The highest total number of SNPs detected in this assay were mapped to chr.2B (14,273 SNPs), which contains the highest number of targeted base pairs in the assay. The average mutation density across all regions interrogated was estimated to be one mutation per 20.91 Mb. The highest mutation frequency was found in chr.2D (1/11.7 kb) and the lowest in chr.7D (1/353.4 kb). Out of the detected mutations, 101 SNPs were filtered using analysis criteria aimed to enrich for mutations that may affect gene function. Out of these, one putative SNP detected in Lr21 were selected for further analysis. The SNP identified in chimeric allele (Lr21) of a resistant mutant (N1-252) was located in a NBS domain of chr.1BS at 3.4 Mb position. Through computational analysis, it was demonstrated that this identified SNP causes a substitution of glutamic acid with alanine, resulting in a predicted altered protein structure. This mutation, therefore, is a candidate for contributing to the resistance phenotype in the mutant line. Based on this work, we conclude that the wheat mutant resource developed is useful as a source of novel genetic variation for forward-genetic screens and also as a useful tool for gaining insights into the important biological circuits of different traits of complex genomes like wheat.