A haplotype-led approach to increase the precision of wheat breeding

Wheat2035: Integrating pan-omics and advanced biotechnology for future wheat design

Wheat (Triticum aestivum) production is vital for global food security, providing energy and protein to millions of people worldwide. Recent advancements in wheat research have led to significant increases in production, fueled by technological and scientific innovation. Here, we summarize the major advancements in wheat research, particularly the integration of biotechnologies and a deeper understanding of wheat biology. The shift from multi-omics to pan-omics approaches in wheat research has greatly enhanced our understanding of the complex genome, genomic variations, and regulatory networks to decode complex traits. We also outline key scientific questions, potential research directions, and technological strategies for improving wheat over the next decade. Since global wheat production is expected to increase by 60% in 2050, continued innovation and collaboration are crucial. Integrating biotechnologies and a deeper understanding of wheat biology will be essential for addressing future challenges in wheat production, ensuring sustainable practices and improved productivity.

Haplotype-resolved genome and population genomics provide insights into dioscin biosynthesis and evolutionary history of the medicinal species Dioscorea nipponica

Dioscorea nipponica, a perennial herb widely distributed in the Sino-Japanese Floristic Region, is renowned for its medicinal properties, particularly its ability to produce dioscin. Here, we present a haplotype-resolved genome assembly of the diploid D. nipponica, comprising 511.41 Mb for Haplotype A and 498.29 Mb for Haplotype B, each organized into 10 chromosomes. The two haplotypes exhibited high similarity, with only 2.75% of the allelic genes exhibiting specific expression. Key genes in the dioscin biosynthesis pathway were identified, and expression analysis revealed that the majority (16/21) of genes involved in the first two stages were highly expressed in rhizomes. Notably, significant expansion of the CYP90, CYP94, and UGT73 gene families was observed in dioscin-producing species, highlighting their critical roles in dioscin biosynthesis. Additionally, genome size estimation and population genomic analyses of diverse D. nipponica accessions identified four principal clades in D. nipponica, corresponding to diploid, di-tetra-octoploid, tetraploid, and hexaploid accessions from various geographic regions, with clade A (diploids) further divided into five subclades. Demographic analysis of the diploid D. nipponica revealed a prolonged decline in effective population size from the Pleistocene to the Last Glacial Maximum, with population splits occurring during the mid-to-late Pleistocene. Selective sweep analysis identified key genes, including HD-Zip I, ADH1, SMT1, and CYPs that may contribute to adaptation to high-latitude environments and the geographical variations in dioscin content. Overall, this study enhances our understanding of the genomic architecture, biosynthetic pathways, and evolutionary dynamics of D. nipponica, providing valuable insights into its medicinal potential and evolution.

Exploring intra- and intergenomic variation in haplotype-resolved pangenomes

With advances in long-read sequencing and assembly techniques, haplotype-resolved (phased) genome assemblies are becoming more common, also in the field of plant genomics. Computational tools to effectively explore these phased genomes, particularly for polyploid genomes, are currently limited. Here we describe a new strategy adopting a pangenome approach. To analyse both intra- and intergenomic variation in phased genome assemblies, we have made the software package PanTools ploidy-aware by updating the pangenome graph representation and adding several novel functionalities to assess synteny and gene retention, profile repeats and calculate synonymous and nonsynonymous mutation rates. Using PanTools, we constructed and analysed a pangenome comprising of one diploid and four tetraploid potato cultivars, and a pangenome of five diploid apple species. Both pangenomes show high intra- and intergenomic allelic diversity in terms of gene absence/presence, SNPs, indels and larger structural variants. Our findings show that the new functionalities and visualizations are useful to discover introgressions and detect likely misassemblies in phased genomes. PanTools is available at https://git.wur.nl/bioinformatics/pantools.

Drought-tolerant wheat for enhancing global food security

Wheat is among the most produced grain crops of the world and alone provides a fifth of the world's calories and protein. Wheat has played a key role in food security since the crop served as a Neolithic founder crop for the establishment of world agriculture. Projections showing a decline in global wheat yields in changing climates imply that food security targets could be jeopardized. Increased frequency and intensity of drought occurrence is evident in major wheat-producing regions worldwide, and notably, the wheat-producing area under drought is projected to swell globally by 60% by the end of the 21st century. Wheat yields are significantly reduced due to changes in plant morphological, physiological, biochemical, and molecular activities in response to drought stress. Advances in wheat genetics, multi-omics technologies and plant phenotyping have enhanced the understanding of crop responses to drought conditions. Research has elucidated key genomic regions, candidate genes, signalling molecules and associated networks that orchestrate tolerance mechanisms under drought stress. Robust and low-cost selection tools are now available in wheat for screening genetic variations for drought tolerance traits. New breeding techniques and selection tools open a unique opportunity to tailor future wheat crop with optimal trait combinations that help withstand extreme drought. Adoption of the new wheat varieties will increase crop diversity in rain-fed agriculture and ensure sustainable improvements in crop yields to safeguard the world's food security in drier environments.

Wheat genomics: genomes, pangenomes, and beyond

There is an urgent need to improve wheat for upcoming challenges, including biotic and abiotic stresses. Sustainable wheat improvement requires the introduction of new genes and alleles in high-yielding wheat cultivars. Using new approaches, tools, and technologies to identify and introduce new genes in wheat cultivars is critical. High-quality genomes, transcriptomes, and pangenomes provide essential resources and tools to examine wheat closely to identify and manipulate new and targeted genes and alleles. Wheat genomics has improved excellently in the past 5 years, generating multiple genomes, pangenomes, and transcriptomes. Leveraging these resources allows us to accelerate our crop improvement pipelines. This review summarizes the progress made in wheat genomics and trait discovery in the past 5 years.

Exploring Novel Genomic Loci and Candidate Genes Associated with Plant Height in Bulgarian Bread Wheat via Multi-Model GWAS

In the context of crop breeding, plant height (PH) plays a pivotal role in determining straw and grain yield. Although extensive research has explored the genetic control of PH in wheat, there remains an opportunity for further advancements by integrating genomics with growth-related phenomics. Our study utilizes the latest genome-wide association scan (GWAS) techniques to unravel the genetic basis of temporal variation in PH across 179 Bulgarian bread wheat accessions, including landraces, tall historical, and semi-dwarf modern varieties. A GWAS was performed with phenotypic data from three growing seasons, the calculated best linear unbiased estimators, and the leveraging genotypic information from the 25K Infinium iSelect array, using three statistical methods (MLM, FarmCPU, and BLINK). Twenty-five quantitative trait loci (QTL) associated with PH were identified across fourteen chromosomes, encompassing 21 environmentally stable quantitative trait nucleotides (QTNs), and four haplotype blocks. Certain loci (17) on chromosomes 1A, 1B, 1D, 2A, 2D, 3A, 3B, 4A, 5B, 5D, and 6A remain unlinked to any known Rht (Reduced height) genes, QTL, or GWAS loci associated with PH, and represent novel regions of potential breeding significance. Notably, these loci exhibit varying effects on PH, contribute significantly to natural variance, and are expressed during seedling to reproductive stages. The haplotype block on chromosome 6A contains five QTN loci associated with reduced height and two loci promoting height. This configuration suggests a substantial impact on natural variation and holds promise for accurate marker-assisted selection. The potentially novel genomic regions harbor putative candidate gene coding for glutamine synthetase, gibberellin 2-oxidase, auxin response factor, ethylene-responsive transcription factor, and nitric oxide synthase; cell cycle-related genes, encoding cyclin, regulator of chromosome condensation (RCC1) protein, katanin p60 ATPase-containing subunit, and expansins; genes implicated in stem mechanical strength and defense mechanisms, as well as gene regulators such as transcription factors and protein kinases. These findings enrich the pool of semi-dwarfing gene resources, providing the potential to further optimize PH, improve lodging resistance, and achieve higher grain yields in bread wheat.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

Accelerating wheat improvement through trait characterization: advances and perspectives

Wheat (Triticum spp.) is a primary dietary staple food for humanity. Many wheat genetic resources with variable genomes have a record of domestication history and are widespread throughout the world. To develop elite wheat varieties, agronomical and stress-responsive trait characterization is foremost for evaluating existing germplasm to promote breeding. However, genomic complexity is one of the primary impediments to trait mining and characterization. Multiple reference genomes and cutting-edge technologies like haplotype mapping, genomic selection, precise gene editing tools, high-throughput phenotyping platforms, high-efficiency genetic transformation systems, and speed-breeding facilities are transforming wheat functional genomics research to understand the genomic diversity of polyploidy. This review focuses on the research achievements in wheat genomics, the available omics approaches, and bioinformatic resources developed in the past decades. Advances in genomics and system biology approaches are highlighted to circumvent bottlenecks in genomic and phenotypic selection, as well as gene transfer. In addition, we propose conducting precise functional genomic studies and developing sustainable breeding strategies for wheat. These developments in understanding wheat traits have speed up the creation of high-yielding, stress-resistant, and nutritionally enhanced wheat varieties, which will help in addressing global food security and agricultural sustainability in the era of climate change.

Harnessing landrace diversity empowers wheat breeding

Harnessing genetic diversity in major staple crops through the development of new breeding capabilities is essential to ensure food security1. Here we examined the genetic and phenotypic diversity of the A. E. Watkins landrace collection2 of bread wheat (Triticum aestivum), a major global cereal, by whole-genome re-sequencing of 827 Watkins landraces and 208 modern cultivars and in-depth field evaluation spanning a decade. We found that modern cultivars are derived from two of the seven ancestral groups of wheat and maintain very long-range haplotype integrity. The remaining five groups represent untapped genetic sources, providing access to landrace-specific alleles and haplotypes for breeding. Linkage disequilibrium-based haplotypes and association genetics analyses link Watkins genomes to the thousands of identified high-resolution quantitative trait loci and significant marker-trait associations. Using these structured germplasm, genotyping and informatics resources, we revealed many Watkins-unique beneficial haplotypes that can confer superior traits in modern wheat. Furthermore, we assessed the phenotypic effects of 44,338 Watkins-unique haplotypes, introgressed from 143 prioritized quantitative trait loci in the context of modern cultivars, bridging the gap between landrace diversity and current breeding. This study establishes a framework for systematically utilizing genetic diversity in crop improvement to achieve sustainable food security.

Haplotype-based association mapping of genomic regions associated with Zymoseptoria tritici resistance using 217 diverse wheat genotypes

BACKGROUND

Septoria tritici blotch (STB) is considered to be one of the most destructive foliar wheat diseases and is caused by Zymoseptoria tritici. The yield losses are severe and in Northwestern Europe can reach up to 50%. The efficacy of fungicides is diminishing due to changes in the genetic structure of the pathogen. Therefore, resistance breeding is the most effective strategy of disease management. Recently, genome-wide association studies (GWAS) have become more popular due to their robustness in dissecting complex traits, including STB resistance in wheat. This was made possible by the use of large mapping populations and new sequencing technologies. High-resolution mapping benefits from historical recombination and greater allele numbers in GWAS.

RESULTS

In our study, 217 wheat genotypes of diverse origin were phenotyped against five Z. tritici isolates (IPO323, IPO88004, IPO92004, IPO86036 and St1-03) and genotyped on the DArTseq platform. In polytunnel tests two disease parameters were evaluated: the percentage of leaf area covered by necrotic lesions (NEC) and the percentage of leaf area covered by lesions bearing pycnidia (PYC). The disease escape parameters heading date (Hd) and plant height (Ht) were also measured. Pearson's correlation showed a positive effect between disease parameters, providing additional information. The Structure analysis indicated four subpopulations which included from 28 (subpopulation 2) to 79 genotypes (subpopulation 3). All of the subpopulations showed a relatively high degree of admixture, which ranged from 60% of genotypes with less than 80% of proportions of the genome attributed to assigned subpopulation for group 2 to 85% for group 4. Haplotype-based GWAS analysis allowed us to identify 27 haploblocks (HBs) significantly associated with analysed traits with a p-value above the genome-wide significance threshold (5%, which was -log10(p) > 3.64) and spread across the wheat genome. The explained phenotypic variation of identified significant HBs ranged from 0.2% to 21.5%. The results of the analysis showed that four haplotypes (HTs) associated with disease parameters cause a reduction in the level of leaf coverage by necrosis and pycnidia, namely: Chr3A_HB98_HT2, Chr5B_HB47_HT1, Chr7B_HB36_HT1 and Chr5D_HB10_HT3.

CONCLUSIONS

GWAS analysis enabled us to identify four significant chromosomal regions associated with a reduction in STB disease parameters. The list of valuable HBs and wheat varieties possessing them provides promising material for further molecular analysis of resistance loci and development of breeding programmes.

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

Climate Change-The Rise of Climate-Resilient Crops

Climate change disrupts food production in many regions of the world. The accompanying extreme weather events, such as droughts, floods, heat waves, and cold snaps, pose threats to crops. The concentration of carbon dioxide also increases in the atmosphere. The United Nations is implementing the climate-smart agriculture initiative to ensure food security. An element of this project involves the breeding of climate-resilient crops or plant cultivars with enhanced resistance to unfavorable environmental conditions. Modern agriculture, which is currently homogeneous, needs to diversify the species and cultivars of cultivated plants. Plant breeding programs should extensively incorporate new molecular technologies, supported by the development of field phenotyping techniques. Breeders should closely cooperate with scientists from various fields of science.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

GreenHill: a de novo chromosome-level scaffolding and phasing tool using Hi-C

Chromosome-level haplotype-resolved genome assembly is an important resource in molecular biology. However, current de novo haplotype assemblers require parental data or reference genomes and often fail to provide chromosome-level results. We present GreenHill, a novel scaffolding and phasing tool that considers various assemblers' contigs as input to reconstruct chromosome-level haplotypes using Hi-C without parental or reference data. Its unique functions include new error correction based on Hi-C contacts and the simultaneous use of Hi-C and long reads. Benchmarks reveal that GreenHill outperforms other approaches in contiguity and phasing accuracy, and the majority of chromosome arms are entirely phased.

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

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

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

Wheat EARLY FLOWERING 3 affects heading date without disrupting circadian oscillations

Plant breeders have indirectly selected for variation at circadian-associated loci in many of the world's major crops, when breeding to increase yield and improve crop performance. Using an eight-parent Multiparent Advanced Generation Inter-Cross (MAGIC) population, we investigated how variation in circadian clock-associated genes contributes to the regulation of heading date in UK and European winter wheat (Triticum aestivum) varieties. We identified homoeologues of EARLY FLOWERING 3 (ELF3) as candidates for the Earliness per se (Eps) D1 and B1 loci under field conditions. We then confirmed a single-nucleotide polymorphism within the coding region of TaELF3-B1 as a candidate polymorphism underlying the Eps-B1 locus. We found that a reported deletion at the Eps-D1 locus encompassing TaELF3-D1 is, instead, an allele that lies within an introgression region containing an inversion relative to the Chinese Spring D genome. Using Triticum turgidum cv. Kronos carrying loss-of-function alleles of TtELF3, we showed that ELF3 regulates heading, with loss of a single ELF3 homoeologue sufficient to alter heading date. These studies demonstrated that ELF3 forms part of the circadian oscillator; however, the loss of all homoeologues was required to affect circadian rhythms. Similarly, loss of functional LUX ARRHYTHMO (LUX) in T. aestivum, an orthologue of a protein partner of Arabidopsis (Arabidopsis thaliana) ELF3, severely disrupted circadian rhythms. ELF3 and LUX transcripts are not co-expressed at dusk, suggesting that the structure of the wheat circadian oscillator might differ from that of Arabidopsis. Our demonstration that alterations to ELF3 homoeologues can affect heading date separately from effects on the circadian oscillator suggests a role for ELF3 in cereal photoperiodic responses that could be selected for without pleiotropic deleterious alterations to circadian rhythms.

A 'wiring diagram' for sink strength traits impacting wheat yield potential

Identifying traits for improving sink strength is a bottleneck to increasing wheat yield. The interacting processes determining sink strength and yield potential are reviewed and visualized in a set of 'wiring diagrams', covering critical phases of development (and summarizing known underlying genetics). Using this framework, we reviewed and assembled the main traits determining sink strength and identified research gaps and potential hypotheses to be tested for achieving gains in sink strength. In pre-anthesis, grain number could be increased through: (i) enhanced spike growth associated with optimized floret development and/or a reduction in specific stem-internode lengths and (ii) improved fruiting efficiency through an accelerated rate of floret development, improved partitioning between spikes, or optimized spike cytokinin levels. In post-anthesis, grain, sink strength could be augmented through manipulation of grain size potential via ovary size and/or endosperm cell division and expansion. Prospects for improving spike vascular architecture to support all rapidly growing florets, enabling the improved flow of assimilate, are also discussed. Finally, we considered the prospects for enhancing grain weight realization in relation to genetic variation in stay-green traits as well as stem carbohydrate remobilization. The wiring diagrams provide a potential workspace for breeders and crop scientists to achieve yield gains in wheat and other field crops.

An autoactive NB-LRR gene causes Rht13 dwarfism in wheat

Semidwarfing genes have greatly increased wheat yields globally, yet the widely used gibberellin (GA)-insensitive genes Rht-B1b and Rht-D1b have disadvantages for seedling emergence. Use of the GA-sensitive semidwarfing gene Rht13 avoids this pleiotropic effect. Here, we show that Rht13 encodes a nucleotide-binding site/leucine-rich repeat (NB-LRR) gene. A point mutation in the semidwarf Rht-B13b allele autoactivates the NB-LRR gene and causes a height reduction comparable with Rht-B1b and Rht-D1b in diverse genetic backgrounds. The autoactive Rht-B13b allele leads to transcriptional up-regulation of pathogenesis-related genes including class III peroxidases associated with cell wall remodeling. Rht13 represents a new class of reduced height (Rht) gene, unlike other Rht genes, which encode components of the GA signaling or metabolic pathways. This discovery opens avenues to use autoactive NB-LRR genes as semidwarfing genes in a range of crop species, and to apply Rht13 in wheat breeding programs using a perfect genetic marker.

Deciphering diversity at er loci for diversification of powdery mildew resistance in pea

Agricultural biotechnology aims to scrutinize the field crops which feed half of the world's population by improving their agronomic traits using various biotechnological tools. Pea- an important cash crop, rich in nutrients, but frequently infected with powdery mildew (fungal disease caused by Erysiphe pisi) that destroys the whole crop and causes economic loss for growers. We, therefore, targeted this research to find the pathogen-resistant pea lines and further decipher the diversity at er locus among resistant pea lines. Screening for resistant pea lines was done with Erysiphe pisi isolates (Genebank submission: KX455922.1) under the net house and greenhouse conditions. Molecular studies revealed that the Erysiphe resistant (er1) gene was present in 40 lines out of selected 50 pea lines and the mutational character was conferred up to 36 genotypes with 11 haplotype groups. The haplotype (gene) diversity (Hd) was found to be 0.5571 ± 0.099 SD and the nucleotide diversity (Pi) was 0.0160 ± 0.0042 SD Majority of resistant lines (67%) occurred in Hap-1, other remaining haplotypes (Hap 2-10) having 33% resistant lines, each showing characteristic nucleotide substitutions with respect to reference PsMLO1 gene; genotypes from these divergent haplotypes can be used in pea resistance breeding to avoid genetic homogeneity and genetic vulnerability.

Discovery and Validation of Grain Shape Loci in U.S. Rice Germplasm Through Haplotype Characterization

Rice grain shape is a major determinant of rice market value and the end-use. We mapped quantitative trait loci (QTL) for grain shape traits in a bi-parental recombinant inbred line population (Trenasse/Jupiter) and discovered two major grain length QTLs—qGL3.1 and qGL7.1. Previously, a major grain shape gene GS3 was reported in the qGL3.1 region and grain length gene GL7 was reported to be encompassing qGL7.1 locus. The re-sequencing SNP data on the International Rice Research Institute (IRRI) 3K Rice Genome Project (RGP) panel were obtained from the IRRI SNP-Seek database for both genes and haplotype diversity was characterized for each gene in this diverse panel. United States rice germplasm was not well represented in the IRRI 3K RGP database. Therefore, a minimum SNP set was identified for each gene that could differentiate all the characterized haplotypes. These haplotypes in the 3K RGP panel were screened across 323 elite U.S. genotypes using the minimum SNP set. The screening of haplotypes and phenotype association confirmed the role of GS3 under qGL3.1. However, screening of the GL7 haplotypes in the U.S. germplasm panel showed that GL7 did not play a role in qGL7.1, and in addition, GL7.1 did not segregate in the Trenasse/Jupiter RIL population. This concluded that qGL7.1 is a novel QTL discovered on chr7 for grain shape in the Trenasse/Jupiter RIL population. A high-throughput KASP-based SNP marker for each locus (GS3 and qGL7.1) was identified and validated in elite U.S. rice germplasm to be used in an applied rice breeding program.

TaNAC020 homoeologous genes are associated with higher thousand kernel weight and kernel length in Chinese wheat

NAC proteins constitute one of the largest plant-specific transcription factor (TF) families and play significant roles in plant growth and development. In the present study, three TaNAC020 homoeologous genes located on chromosomes 7A, 7B, and 7D were isolated from wheat (Triticum aestivum L.). TaNAC020s were predominantly expressed in developing grains. The developed transgenic rice lines for TaNAC020-B showed higher starch density and lower amylose contents than those of the wild type (WT). Sequence polymorphism studies showed seven and eight SNPs in TaNAC020-A/B, making three and two haplotypes, respectively. No sequence polymorphism was identified in TaNAC020-D. Association analysis revealed that HAP-2 of TaNAC020-A and TaNAC020-B was the favored haplotype for higher thousand kernel weight and length. Geographic distribution and allelic frequency showed that our favored haplotype experienced strong selection in China, and likewise, diversity increased in TaNAC020s during wheat polyploidization. The results obtained in this study demonstrate that TaNAC020s positively influence starch synthesis and accumulation and are one of the key regulators of the kernel (seed) size and kernel number and have the potential for utilization in wheat breeding to improve grain yield. Molecular markers developed in this study stand instrumental in marker-assisted selection for genetic improvement and germplasm enhancement in wheat.

Dispersed emergence and protracted domestication of polyploid wheat uncovered by mosaic ancestral haploblock inference

Major crops are all survivors of domestication bottlenecks. Studies have focused on the genetic loci related to the domestication syndrome, while the contribution of ancient haplotypes remains largely unknown. Here, an ancestral genomic haploblock dissection method is developed and applied to a resequencing dataset of 386 tetraploid/hexaploid wheat accessions, generating a pan-ancestry haploblock map. Together with cytoplastic evidences, we reveal that domesticated polyploid wheat emerged from the admixture of six founder wild emmer lineages, which contributed the foundation of ancestral mosaics. The key domestication-related loci, originated over a wide geographical range, were gradually pyramided through a protracted process. Diverse stable-inheritance ancestral haplotype groups of the chromosome central zone are identified, revealing the expanding routes of wheat and the trends of modern wheat breeding. Finally, an evolution model of polyploid wheat is proposed, highlighting the key role of wild-to-crop and interploidy introgression, that increased genomic diversity following bottlenecks introduced by domestication and polyploidization.

Characterization of Expression and Epigenetic Features of Core Genes in Common Wheat

The availability of multiple wheat genome sequences enables us to identify core genes and characterize their genetic and epigenetic features, thereby advancing our understanding of their biological implications within individual plant species. It is, however, largely understudied in wheat. To this end, we reanalyzed genome sequences from 16 different wheat varieties and identified 62,299 core genes. We found that core and non-core genes have different roles in subgenome differentiation. Meanwhile, according to their expression profiles, these core genes can be classified into genes related to tissue development and stress responses, including 3376 genes highly expressed in both spikelets and at high temperatures. After associating with six histone marks and open chromatin, we found that these core genes can be divided into eight sub-clusters with distinct epigenomic features. Furthermore, we found that ca. 51% of the expressed transcription factors (TFs) were marked with both H3K27me3 and H3K4me3, indicative of the bivalency feature, which can be involved in tissue development through the TF-centered regulatory network. Thus, our study provides a valuable resource for the functional characterization of core genes in stress responses and tissue development in wheat.

Trade-offs in the genetic control of functional and nutritional quality traits in UK winter wheat

A complex network of trade-offs exists between wheat quality and nutritional traits. We investigated the correlated relationships among several milling and baking traits as well as mineral density in refined white and whole grain flour. Our aim was to determine their pleiotropic genetic control in a multi-parent population over two trial years with direct application to practical breeding. Co-location of major quantitative trait loci (QTL) and principal component based multi-trait QTL mapping increased the power to detect QTL and revealed pleiotropic effects explaining many complementary and antagonistic trait relationships. High molecular weight glutenin subunit genes explained much of the heritable variation in important dough rheology traits, although additional QTL were detected. Several QTL, including one linked to the TaGW2 gene, controlled grain size and increased flour extraction rate. The semi-dwarf Rht-D1b allele had a positive effect on Hagberg falling number, but reduced grain size, specific weight, grain protein content and flour water absorption. Mineral nutrient concentrations were lower in Rht-D1b lines for many elements, in wholemeal and white flour, but potassium concentration was higher in Rht-D1b lines. The presence of awns increased calcium content without decreasing extraction rate, despite the negative correlation between these traits. QTL were also found that affect the relative concentrations of key mineral nutrients compared to phosphorus which may help increase bioavailability without associated anti-nutritional effects of phytic acid. Taken together these results demonstrate the potential for marker-based selection to optimise trait trade-offs and enhance wheat nutritional value by considering pleiotropic genetic effects across multiple traits.

Chromosome-Level Haplotype Assembly for Equus asinu

Background: Haplotype provides significant insights into understanding genomes at both individual and population levels. However, research on many non-model organisms is still based on independent genetic variations due to the lack of haplotype.

Results: We conducted haplotype assembling for Equus asinu, a non-model organism that plays a vital role in human civilization. We described the hybrid single individual assembled haplotype of the Dezhou donkey based on the high-depth sequencing data from single-molecule real-time sequencing (×30), Illumina short-read sequencing (×211), and high-throughput chromosome conformation capture (×56). We assembled a near-complete haplotype for the high-depth sequenced Dezhou donkey individual and a phased cohort for the resequencing data of the donkey population.

Conclusion: Here, we described the complete chromosome-scale haplotype of the Dezhou donkey with more than a 99.7% phase rate. We further phased a cohort of 156 donkeys to form a donkey haplotype dataset with more than 39 million genetic variations.

High-resolution detection of quantitative trait loci for seven important yield-related traits in wheat (Triticum aestivum L.) using a high-density SLAF-seq genetic map

BACKGROUND

Yield-related traits including thousand grain weight (TGW), grain number per spike (GNS), grain width (GW), grain length (GL), plant height (PH), spike length (SL), and spikelet number per spike (SNS) are greatly associated with grain yield of wheat (Triticum aestivum L.). To detect quantitative trait loci (QTL) associated with them, 193 recombinant inbred lines derived from two elite winter wheat varieties Chuanmai42 and Chuanmai39 were employed to perform QTL mapping in six/eight environments.

RESULTS

A total of 30 QTLs on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 4A, 5A, 5B, 6A, 6D, 7A, 7B and 7D were identified. Among them, six major QTLs QTgw.cib-6A.1, QTgw.cib-6A.2, QGw.cib-6A, QGl.cib-3A, QGl.cib-6A, and QSl.cib-2D explaining 5.96-23.75% of the phenotypic variance were detected in multi-environments and showed strong and stable effects on corresponding traits. Three QTL clusters on chromosomes 2D and 6A containing 10 QTLs were also detected, which showed significant pleiotropic effects on multiple traits. Additionally, three Kompetitive Allele Specific PCR (KASP) markers linked with five of these major QTLs were developed. Candidate genes of QTgw.cib-6A.1/QGl.cib-6A and QGl.cib-3A were analyzed based on the spatiotemporal expression patterns, gene annotation, and orthologous search.

CONCLUSIONS

Six major QTLs for TGW, GL, GW and SL were detected. Three KASP markers linked with five of these major QTLs were developed. These QTLs and KASP markers will be useful for elucidating the genetic architecture of grain yield and developing new wheat varieties with high and stable yield in wheat.

Long-read and chromosome-scale assembly of the hexaploid wheat genome achieves high resolution for research and breeding

BACKGROUND

The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost.

RESULTS

Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan.

CONCLUSIONS

We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

Genome sequences of three Aegilops species of the section Sitopsis reveal phylogenetic relationships and provide resources for wheat improvement

Aegilops is a close relative of wheat (Triticum spp.), and Aegilops species in the section Sitopsis represent a rich reservoir of genetic diversity for the improvement of wheat. To understand their diversity and advance their utilization, we produced whole-genome assemblies of Aegilops longissima and Aegilops speltoides. Whole-genome comparative analysis, along with the recently sequenced Aegilops sharonensis genome, showed that the Ae. longissima and Ae. sharonensis genomes are highly similar and are most closely related to the wheat D subgenome. By contrast, the Ae. speltoides genome is more closely related to the B subgenome. Haplotype block analysis supported the idea that Ae. speltoides genome is closest to the wheat B subgenome, and highlighted variable and similar genomic regions between the three Aegilops species and wheat. Genome-wide analysis of nucleotide-binding leucine-rich repeat (NLR) genes revealed species-specific and lineage-specific NLR genes and variants, demonstrating the potential of Aegilops genomes for wheat improvement.

ggComp enables dissection of germplasm resources and construction of a multiscale germplasm network in wheat

Accurate germplasm characterization is a vital step for accelerating crop genetic improvement, which remains largely infeasible for crops such as bread wheat (Triticum aestivum L.), which has a complex genome that undergoes frequent introgression and contains many structural variations. Here, we propose a genomic strategy called ggComp, which integrates resequencing data with copy number variations and stratified single-nucleotide polymorphism densities to enable unsupervised identification of pairwise germplasm resource-based Identity-By-Descent (gIBD) blocks. The reliability of ggComp was verified in wheat cultivar Nongda5181 by dissecting parental-descent patterns represented by inherited genomic blocks. With gIBD blocks identified among 212 wheat accessions, we constructed a multi-scale genomic-based germplasm network. At the whole-genome level, the network helps to clarify pedigree relationship, demonstrate genetic flow, and identify key founder lines. At the chromosome level, we were able to trace the utilization of 1RS introgression in modern wheat breeding by hitchhiked segments. At the single block scale, the dissected germplasm-based haplotypes nicely matched with previously identified alleles of "Green Revolution" genes and can guide allele mining and dissect the trajectory of beneficial alleles in wheat breeding. Our work presents a model-based framework for precisely evaluating germplasm resources with genomic data. A database, WheatCompDB (http://wheat.cau.edu.cn/WheatCompDB/), is available for researchers to exploit the identified gIBDs with a multi-scale network.

Genomic prediction of morphometric and colorimetric traits in Solanaceous fruits

Selection of high-performance lines with respect to traits of interest is a key step in plant breeding. Genomic prediction allows to determine the genomic estimated breeding values of unseen lines for trait of interest using genetic markers, e.g. single-nucleotide polymorphisms (SNPs), and machine learning approaches, which can therefore shorten breeding cycles, referring to genomic selection (GS). Here, we applied GS approaches in two populations of Solanaceous crops, i.e. tomato and pepper, to predict morphometric and colorimetric traits. The traits were measured by using scoring-based conventional descriptors (CDs) as well as by Tomato Analyzer (TA) tool using the longitudinally and latitudinally cut fruit images. The GS performance was assessed in cross-validations of classification-based and regression-based machine learning models for CD and TA traits, respectively. The results showed the usage of TA traits and tag SNPs provide a powerful combination to predict morphology and color-related traits of Solanaceous fruits. The highest predictability of 0.89 was achieved for fruit width in pepper, with an average predictability of 0.69 over all traits. The multi-trait GS models are of slightly better predictability than single-trait models for some colorimetric traits in pepper. While model validation performs poorly on wild tomato accessions, the usage as many as one accession per wild species in the training set can increase the transferability of models to unseen populations for some traits (e.g. fruit shape for which predictability in unseen scenario increased from zero to 0.6). Overall, GS approaches can assist the selection of high-performance Solanaceous fruits in crop breeding.

Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome

Common wheat (Triticum aestivum, BBAADD) is a major staple food crop worldwide. The diploid progenitors of the A and D subgenomes have been unequivocally identified; that of B, however, remains ambiguous and controversial but is suspected to be related to species of Aegilops, section Sitopsis. Here, we report the assembly of chromosome-level genome sequences of all five Sitopsis species, namely Aegilops bicornis, Ae. longissima, Ae. searsii, Ae. sharonensis, and Ae. speltoides, as well as the partial assembly of the Amblyopyrum muticum (synonym Aegilops mutica) genome for phylogenetic analysis. Our results reveal that the donor of the common wheat B subgenome is a distinct, and most probably extinct, diploid species that diverged from an ancestral progenitor of the B lineage to which the still extant Ae. speltoides and Am. muticum belong. In addition, we identified interspecific genetic introgressions throughout the evolution of the Triticum/Aegilops species complex. The five Sitopsis species have various assembled genome sizes (4.11-5.89 Gb) with high proportions of repetitive sequences (85.99%-89.81%); nonetheless, they retain high collinearity with other genomes or subgenomes of species in the Triticum/Aegilops complex. Differences in genome size were primarily due to independent post-speciation amplification of transposons. We also identified a set of Sitopsis genes pertinent to important agronomic traits that can be harnessed for wheat breeding. These newly assembled genome resources provide a new roadmap for evolutionary and genetic studies of the Triticum/Aegilops complex, as well as for wheat improvement.

Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines

Many wild-relative species are being used in prebreeding programs to increase the genetic diversity of wheat (Triticum aestivum L.). Genotyping tools such as single nucleotide polymorphism (SNP)-based arrays and molecular markers have been widely used to characterize wheat-wild relative introgression lines. However, due to the polyploid nature of the recipient wheat genome, it is difficult to develop SNP-based Kompetitive allele-specific polymerase chain reaction (KASP) markers that are codominant to track the introgressions from the wild species. Previous attempts to develop KASP markers have involved both exome- and polymerase chain reaction (PCR)-amplicon-based sequencing of the wild species. But chromosome-specific KASP assays have been hindered by homoeologous SNPs within the wheat genome. This study involved whole genome sequencing of the diploid wheat wild relative Amblyopyrum muticum (Boiss.) Eig and development of a de novo SNP discovery pipeline that generated ∼38,000 SNPs in unique wheat genome sequences. New assays were designed to increase the density of Am. muticum polymorphic KASP markers. With a goal of one marker per 60 Mbp, 335 new KASP assays were validated as diagnostic for Am. muticum in a wheat background. Together with assays validated in previous studies, 498 well distributed chromosome-specific markers were used to recharacterize previously genotyped wheat-Am. muticum doubled haploid (DH) introgression lines. The chromosome-specific nature of the KASP markers allowed clarification of which wheat chromosomes were involved with recombination events or substituted with Am. muticum chromosomes and the higher density of markers allowed detection of new small introgressions in these DH lines.

The 3366 chickpea genomes for research and breeding

Genome sequences provide an unprecedented resource to rapidly develop modern crops. A recent paper by Varshney et al. provides genome variation maps of 3366 chickpea accessions. Here, we highlight how this breakthrough research can fundamentally change breeding practices of chickpea and potentially other crops.

Wheat Breeding, Transcription Factories, and Genetic Interactions: New Perspectives

Epistatic interactions and negative heterosis have been shown to be associated with interchromosomal interactions in wheat. Physical gene-gene interactions between co-regulated genes clustered in "transcription factories" have been documented, and a genome-wide atlas of functionally paired, interacting regulatory elements and genes of wheat recently produced. Integration of these new studies on gene and regulatory element interactions, co-regulation of gene expression in "transcription factories," and epigenetics generates new perspectives for wheat breeding and trait enhancement.

Trend, population structure, and trait mapping from 15 years of national varietal trials of UK winter wheat

There are now a rich variety of genomic and genotypic resources available to wheat researchers and breeders. However, the generation of high-quality and field-relevant phenotyping data which is required to capture the complexities of gene × environment interactions remains a major bottleneck. Historical datasets from national variety performance trials (NVPT) provide sufficient dimensions, in terms of numbers of years and locations, to examine phenotypic trends and study gene × environment interactions. Using NVPT for winter wheat varieties grown in the United Kingdom between 2002 and 2017, we examined temporal trends for eight traits related to yield, adaptation, and grain quality performance. We show a non-stationary linear trend for yield, grain protein content, Hagberg Falling Number (HFN), and days to ripening. Our data also show high environmental stability for yield, grain protein content, and specific weight in UK winter wheat varieties and high environmental sensitivity for HFN. We also show that UK varieties released within this period cluster into four main population groups. Using the historical NVPT data in a genome-wide association analysis, we uncovered a significant marker-trait association peak on wheat chromosome 6A spanning the NAM-A1 gene that have been previously associated with early senescence. Together, our results show the value of utilizing the data routinely collected during national variety evaluation process for examining breeding progress and the genetic architecture of important traits.

Different Shades of Kale-Approaches to Analyze Kale Variety Interrelations

Brassica oleracea is a vegetable crop with an amazing morphological diversity. Among the various crops derived from B. oleracea, kale has been in the spotlight globally due to its various health-benefitting compounds and many different varieties. Knowledge of the existing genetic diversity is essential for the improved breeding of kale. Here, we analyze the interrelationships, population structures, and genetic diversity of 72 kale and cabbage varieties by extending our previous diversity analysis and evaluating the use of summed potential lengths of shared haplotypes (SPLoSH) as a new method for such analyses. To this end, we made use of the high-density Brassica 60K SNP array, analyzed SNPs included in an available Brassica genetic map, and used these resources to generate and evaluate the information from SPLoSH data. With our results we could consistently differentiate four groups of kale across all analyses: the curly kale varieties, Italian, American, and Russian varieties, as well as wild and cultivated types. The best results were achieved by using SPLoSH information, thus validating the use of this information in improving analyses of interrelations in kale. In conclusion, our definition of kale includes the curly varieties as the kales in a strict sense, regardless of their origin. These results contribute to a better understanding of the huge diversity of kale and its interrelations.

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.

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

Genome-wide association study uncovers major genetic loci associated with seed flooding tolerance in soybean

BACKGROUND

Seed flooding stress is one of the threatening environmental stressors that adversely limits soybean at the germination stage across the globe. The knowledge on the genetic basis underlying seed-flooding tolerance is limited. Therefore, we performed a genome-wide association study (GWAS) using 34,718 single nucleotide polymorphism (SNPs) in a panel of 243 worldwide soybean collections to identify genetic loci linked to soybean seed flooding tolerance at the germination stage.

RESULTS

In the present study, GWAS was performed with two contrasting models, Mixed Linear Model (MLM) and Multi-Locus Random-SNP-Effect Mixed Linear Model (mrMLM) to identify significant SNPs associated with electrical conductivity (EC), germination rate (GR), shoot length (ShL), and root length (RL) traits at germination stage in soybean. With MLM, a total of 20, 40, 4, and 9 SNPs associated with EC, GR, ShL and RL, respectively, whereas in the same order mrMLM detected 27, 17, 13, and 18 SNPs. Among these SNPs, two major SNPs, Gm_08_11971416, and Gm_08_46239716 were found to be consistently connected with seed-flooding tolerance related traits, namely EC and GR across two environments. We also detected two SNPs, Gm_05_1000479 and Gm_01_53535790 linked to ShL and RL, respectively. Based on Gene Ontology enrichment analysis, gene functional annotations, and protein-protein interaction network analysis, we predicted eight candidate genes and three hub genes within the regions of the four SNPs with Cis-elements in promoter regions which may be involved in seed-flooding tolerance in soybeans and these warrant further screening and functional validation.

CONCLUSIONS

Our findings demonstrate that GWAS based on high-density SNP markers is an efficient approach to dissect the genetic basis of complex traits and identify candidate genes in soybean. The trait associated SNPs could be used for genetic improvement in soybean breeding programs. The candidate genes could help researchers better understand the molecular mechanisms underlying seed-flooding stress tolerance in soybean.

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

Fast-forward breeding for a food-secure world

Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.

Crossover-active regions of the wheat genome are distinguished by DMC1, the chromosome axis, H3K27me3, and signatures of adaptation

The hexaploid bread wheat genome comprises over 16 gigabases of sequence across 21 chromosomes. Meiotic crossovers are highly polarized along the chromosomes, with elevation in the gene-dense distal regions and suppression in the Gypsy retrotransposon-dense centromere-proximal regions. We profiled the genomic landscapes of the meiotic recombinase DMC1 and the chromosome axis protein ASY1 in wheat and investigated their relationships with crossovers, chromatin state, and genetic diversity. DMC1 and ASY1 chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed strong co-enrichment in the distal, crossover-active regions of the wheat chromosomes. Distal ChIP-seq enrichment is consistent with spatiotemporally biased cytological immunolocalization of DMC1 and ASY1 close to the telomeres during meiotic prophase I. DMC1 and ASY1 ChIP-seq peaks show significant overlap with genes and transposable elements in the Mariner and Mutator superfamilies. However, DMC1 and ASY1 ChIP-seq peaks were detected along the length of each chromosome, including in low-crossover regions. At the fine scale, crossover elevation at DMC1 and ASY1 peaks and genes correlates with enrichment of the Polycomb histone modification H3K27me3. This indicates a role for facultative heterochromatin, coincident with high DMC1 and ASY1, in promoting crossovers in wheat and is reflected in distalized H3K27me3 enrichment observed via ChIP-seq and immunocytology. Genes with elevated crossover rates and high DMC1 and ASY1 ChIP-seq signals are overrepresented for defense-response and immunity annotations, have higher sequence polymorphism, and exhibit signatures of selection. Our findings are consistent with meiotic recombination promoting genetic diversity, shaping host–pathogen co-evolution, and accelerating adaptation by increasing the efficiency of selection.

Data-driven decentralized breeding increases prediction accuracy in a challenging crop production environment

Crop breeding must embrace the broad diversity of smallholder agricultural systems to ensure food security to the hundreds of millions of people living in challenging production environments. This need can be addressed by combining genomics, farmers' knowledge, and environmental analysis into a data-driven decentralized approach (3D-breeding). We tested this idea as a proof-of-concept by comparing a durum wheat (Triticum durum Desf.) decentralized trial distributed as incomplete blocks in 1,165 farmer-managed fields across the Ethiopian highlands with a benchmark representing genomic prediction applied to conventional breeding. We found that 3D-breeding could double the prediction accuracy of the benchmark. 3D-breeding could identify genotypes with enhanced local adaptation providing superior productive performance across seasons. We propose this decentralized approach to leverage the diversity in farmer fields and complement conventional plant breeding to enhance local adaptation in challenging crop production environments.

Positional-based cloning 'fail-safe' approach is overpowered by wheat chromosome structural variation

Positional-based cloning is a foundational method for understanding the genes and gene networks that control valuable agronomic traits such as grain yield components. In this study, we sought to positionally clone the causal genetic variant of a 1000-grain weight (TGW) quantitative trait loci (QTL) on wheat (Triticum aestivum L.) chromosome arm 5AL. We developed heterogenous inbred families (HIFs) (>5,000 plants) for enhanced genotypic resolution and fine-mapped the QTL to a 10-Mbp region. The transcriptome of developing grains from positive and negative control HIF haplotypes revealed presence-absence chromosome arm 5AS structural variation and unexpectedly no differential expression of genes within the chromosome arm 5AL candidate region. Evaluation of genomic, transcriptomic, and phenotypic data, and predicted function of genes, identified that the 5AL QTL was the result of strong linkage disequilibrium (LD) with chromosome arm 5AS presence or absence (HIF r² = 0.91). Structural variation is common in wheat, and our results highlight that the redundant polyploid genome's masking of such variation is a significant barrier to positional cloning. We propose recommendations for more efficient and robust detection of structural variation, including transitioning from a single nucleotide polymorphism (SNP) to a haplotype-based approach to identify positional cloning targets. We also present nine candidate genes for grain yield components based on chromosome arm 5AS presence or absence, which may unveil hidden variation of homoeolog dosage-dependent genes across the group five chromosome short arms. Taken together, our discovery demonstrates the phenotypic resiliency of polyploid genomic structural variation and highlights a considerable challenge to routine positional cloning in wheat.

Utilization of a Wheat50K SNP Microarray-Derived High-Density Genetic Map for QTL Mapping of Plant Height and Grain Traits in Wheat

Plant height is significantly correlated with grain traits, which is a component of wheat yield. The purpose of this study is to investigate the main quantitative trait loci (QTLs) that control plant height and grain-related traits in multiple environments. In this study, we constructed a high-density genetic linkage map using the Wheat50K SNP Array to map QTLs for these traits in 198 recombinant inbred lines (RILs). The two ends of the chromosome were identified as recombination-rich areas in all chromosomes except chromosome 1B. Both the genetic map and the physical map showed a significant correlation, with a correlation coefficient between 0.63 and 0.99. However, there was almost no recombination between 1RS and 1BS. In terms of plant height, 1RS contributed to the reduction of plant height by 3.43 cm. In terms of grain length, 1RS contributed to the elongation of grain by 0.11 mm. A total of 43 QTLs were identified, including eight QTLs for plant height (PH), 11 QTLs for thousand grain weight (TGW), 15 QTLs for grain length (GL), and nine QTLs for grain width (GW), which explained 1.36–33.08% of the phenotypic variation. Seven were environment-stable QTLs, including two loci (Qph.nwafu-4B and Qph.nwafu-4D) that determined plant height. The explanation rates of phenotypic variation were 7.39–12.26% and 20.11–27.08%, respectively. One QTL, Qtgw.nwafu-4B, which influenced TGW, showed an explanation rate of 3.43–6.85% for phenotypic variation. Two co-segregating KASP markers were developed, and the physical locations corresponding to KASP_AX-109316968 and KASP_AX-109519968 were 25.888344 MB and 25.847691 MB, respectively. Qph.nwafu-4B, controlling plant height, and Qtgw.nwafu-4B, controlling TGW, had an obvious linkage relationship, with a distance of 7–8 cM. Breeding is based on molecular markers that control plant height and thousand-grain weight by selecting strains with low plant height and large grain weight. Another QTL, Qgw.nwafu-4D, which determined grain width, had an explanation rate of 3.43–6.85%. Three loci that affected grain length were Qgl.nwafu-5A, Qgl.nwafu-5D.2, and Qgl.nwafu-6B, illustrating the explanation rates of phenotypic variation as 6.72–9.59%, 5.62–7.75%, and 6.68–10.73%, respectively. Two QTL clusters were identified on chromosomes 4B and 4D.

Designing Future Crops: Genomics-Assisted Breeding Comes of Age

Over the past decade, genomics-assisted breeding (GAB) has been instrumental in harnessing the potential of modern genome resources and characterizing and exploiting allelic variation for germplasm enhancement and cultivar development. Sustaining GAB in the future (GAB 2.0) will rely upon a suite of new approaches that fast-track targeted manipulation of allelic variation for creating novel diversity and facilitate their rapid and efficient incorporation in crop improvement programs. Genomic breeding strategies that optimize crop genomes with accumulation of beneficial alleles and purging of deleterious alleles will be indispensable for designing future crops. In coming decades, GAB 2.0 is expected to play a crucial role in breeding more climate-smart crop cultivars with higher nutritional value in a cost-effective and timely manner.

Crop breeding for a changing climate: integrating phenomics and genomics with bioinformatics

Safeguarding crop yields in a changing climate requires bioinformatics advances in harnessing data from vast phenomics and genomics datasets to translate research findings into climate smart crops in the field. Climate change and an additional 3 billion mouths to feed by 2050 raise serious concerns over global food security. Crop breeding and land management strategies will need to evolve to maximize the utilization of finite resources in coming years. High-throughput phenotyping and genomics technologies are providing researchers with the information required to guide and inform the breeding of climate smart crops adapted to the environment. Bioinformatics has a fundamental role to play in integrating and exploiting this fast accumulating wealth of data, through association studies to detect genomic targets underlying key adaptive climate-resilient traits. These data provide tools for breeders to tailor crops to their environment and can be introduced using advanced selection or genome editing methods. To effectively translate research into the field, genomic and phenomic information will need to be integrated into comprehensive clade-specific databases and platforms alongside accessible tools that can be used by breeders to inform the selection of climate adaptive traits. Here we discuss the role of bioinformatics in extracting, analysing, integrating and managing genomic and phenomic data to improve climate resilience in crops, including current, emerging and potential approaches, applications and bottlenecks in the research and breeding pipeline.

Limited haplotype diversity underlies polygenic trait architecture across 70 years of wheat breeding

Background

Selection has dramatically shaped genetic and phenotypic variation in bread wheat. We can assess the genomic basis of historical phenotypic changes, and the potential for future improvement, using experimental populations that attempt to undo selection through the randomizing effects of recombination.

Results

We bred the NIAB Diverse MAGIC multi-parent population comprising over 500 recombinant inbred lines, descended from sixteen historical UK bread wheat varieties released between 1935 and 2004. We sequence the founders’ genes and promoters by capture, and the MAGIC population by low-coverage whole-genome sequencing. We impute 1.1 M high-quality SNPs that are over 99% concordant with array genotypes. Imputation accuracy only marginally improves when including the founders’ genomes as a haplotype reference panel. Despite capturing 73% of global wheat genetic polymorphism, 83% of genes cluster into no more than three haplotypes. We phenotype 47 agronomic traits over 2 years and map 136 genome-wide significant associations, concentrated at 42 genetic loci with large and often pleiotropic effects. Around half of these overlap known quantitative trait loci. Most traits exhibit extensive polygenicity, as revealed by multi-locus shrinkage modelling.

Conclusions

Our results are consistent with a gene pool of low haplotypic diversity, containing few novel loci of large effect. Most past, and projected future, phenotypic changes arising from existing variation involve fine-scale shuffling of a few haplotypes to recombine dozens of polygenic alleles of small effect. Moreover, extensive pleiotropy means selection on one trait will have unintended consequences, exemplified by the negative trade-off between yield and protein content, unless selection and recombination can break unfavorable trait-trait associations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02354-7.