High throughput SNP discovery and genotyping in hexaploid wheat

Development of a next generation SNP genotyping array for wheat

High-throughput genotyping arrays have provided a cost-effective, reliable and interoperable system for genotyping hexaploid wheat and its relatives. Existing, highly cited arrays including our 35K Wheat Breeder's array and the Illumina 90K array were designed based on a limited amount of varietal sequence diversity and with imperfect knowledge of SNP positions. Recent progress in wheat sequencing has given us access to a vast pool of SNP diversity, whilst technological improvements have allowed us to fit significantly more probes onto a 384-well format Axiom array than previously possible. Here we describe a novel Axiom genotyping array, the 'Triticum aestivum Next Generation' array (TaNG), largely derived from whole genome skim sequencing of 204 elite wheat lines and 111 wheat landraces taken from the Watkins 'Core Collection'. We used a novel haplotype optimization approach to select SNPs with the highest combined varietal discrimination and a design iteration step to test and replace SNPs which failed to convert to reliable markers. The final design with 43 372 SNPs contains a combination of haplotype-optimized novel SNPs and legacy cross-platform markers. We show that this design has an improved distribution of SNPs compared to previous arrays and can be used to generate genetic maps with a significantly higher number of distinct bins than our previous array. We also demonstrate the improved performance of TaNGv1.1 for Genome-wide association studies (GWAS) and its utility for Copy Number Variation (CNV) analysis. The array is commercially available with supporting marker annotations and initial genotyping results freely available.

GeCKO: user-friendly workflows for genotyping complex genomes using target enrichment capture. A use case on the large tetraploid durum wheat genome

BACKGROUND

Genotyping of individuals plays a pivotal role in various biological analyses, with technology choice influenced by multiple factors including genomic constraints, number of targeted loci and individuals, cost considerations, and the ease of sample preparation and data processing. Target enrichment capture of specific polymorphic regions has emerged as a flexible and cost-effective genomic reduction method for genotyping, especially adapted to the case of very large genomes. However, this approach necessitates complex bioinformatics treatment to extract genotyping data from raw reads. Existing workflows predominantly cater to phylogenetic inference, leaving a gap in user-friendly tools for genotyping analysis based on capture methods. In response to these challenges, we have developed GeCKO (Genotyping Complexity Knocked-Out). To assess the effectiveness of combining target enrichment capture with GeCKO, we conducted a case study on durum wheat domestication history, involving sequencing, processing, and analyzing variants in four relevant durum wheat groups.

RESULTS

GeCKO encompasses four distinct workflows, each designed for specific steps of genomic data processing: (i) read demultiplexing and trimming for data cleaning, (ii) read mapping to align sequences to a reference genome, (iii) variant calling to identify genetic variants, and (iv) variant filtering. Each workflow in GeCKO can be easily configured and is executable across diverse computational environments. The workflows generate comprehensive HTML reports including key summary statistics and illustrative graphs, ensuring traceable, reproducible results and facilitating straightforward quality assessment. A specific innovation within GeCKO is its 'targeted remapping' feature, specifically designed for efficient treatment of targeted enrichment capture data. This process consists of extracting reads mapped to the targeted regions, constructing a smaller sub-reference genome, and remapping the reads to this sub-reference, thereby enhancing the efficiency of subsequent steps.

CONCLUSIONS

The case study results showed the expected intra-group diversity and inter-group differentiation levels, confirming the method's effectiveness for genotyping and analyzing genetic diversity in species with complex genomes. GeCKO streamlined the data processing, significantly improving computational performance and efficiency. The targeted remapping enabled straightforward SNP calling in durum wheat, a task otherwise complicated by the species' large genome size. This illustrates its potential applications in various biological research contexts.

Wheat dwarfing reshapes plant and fungal development in arbuscular mycorrhizal symbiosis

The introduction of Reduced height (Rht) dwarfing genes into elite wheat varieties has contributed to enhanced yield gain in high input agrosystems by preventing lodging. Yet, how modern selection for dwarfing has affected symbiosis remains poorly documented. In this study, we evaluated the response of both the plant and the arbuscular mycorrhizal fungus to plant genetic variation at a major Quantitative Trait Locus called QTL 4B2, known to harbor a Rht dwarfing gene, when forming the symbiosis. We used twelve inbred genotypes derived from a diversity base broadened durum wheat Evolutionary Pre-breeding Population and genotyped with a high-throughput Single Nucleotide Polymorphism (SNP) genotyping array. In a microcosm setup segregating roots and the extra-radical mycelium, each wheat genotype was grown with or without the presence of Rhizophagus irregularis. To characterize arbuscular mycorrhizal symbiosis, we assessed hyphal density, root colonization, spore production, and plant biomass. Additionally, we split the variation of these variables due either to genotypes or to the Rht dwarfing genes alone. The fungus exhibited greater development in the roots of Dwarf plants compared to non-Dwarf plants, showing increases of 27%, 37% and 51% in root colonization, arbuscules, and vesicles, respectively. In addition, the biomass of the extra-radical fungal structures increased by around 31% in Dwarf plants. The biomass of plant roots decreased by about 43% in mycorrhizal Dwarf plants. Interestingly, extraradical hyphal production was found to be partly genetically determined with no significant effect of Rht, as for plant biomasses. In contrast, variations in root colonization, arbuscules and extraradical spore production were explained by Rht dwarfing genes. Finally, when mycorrhizal, Dwarf plants had significantly lower total P content, pointing towards a less beneficial symbiosis for the plant and increased profit for the fungus. These results highlight the effect of Rht dwarfing genes on both root and fungal development. This calls for further research into the molecular mechanisms governing these effects, as well as changes in plant physiology, and their implications for fostering arbuscular mycorrhizal symbiosis in sustainable agrosystems.

Reviewing the essential roles of remote phenotyping, GWAS and explainable AI in practical marker-assisted selection for drought-tolerant winter wheat breeding

Marker-assisted selection (MAS) plays a crucial role in crop breeding improving the speed and precision of conventional breeding programmes by quickly and reliably identifying and selecting plants with desired traits. However, the efficacy of MAS depends on several prerequisites, with precise phenotyping being a key aspect of any plant breeding programme. Recent advancements in high-throughput remote phenotyping, facilitated by unmanned aerial vehicles coupled to machine learning, offer a non-destructive and efficient alternative to traditional, time-consuming, and labour-intensive methods. Furthermore, MAS relies on knowledge of marker-trait associations, commonly obtained through genome-wide association studies (GWAS), to understand complex traits such as drought tolerance, including yield components and phenology. However, GWAS has limitations that artificial intelligence (AI) has been shown to partially overcome. Additionally, AI and its explainable variants, which ensure transparency and interpretability, are increasingly being used as recognised problem-solving tools throughout the breeding process. Given these rapid technological advancements, this review provides an overview of state-of-the-art methods and processes underlying each MAS, from phenotyping, genotyping and association analyses to the integration of explainable AI along the entire workflow. In this context, we specifically address the challenges and importance of breeding winter wheat for greater drought tolerance with stable yields, as regional droughts during critical developmental stages pose a threat to winter wheat production. Finally, we explore the transition from scientific progress to practical implementation and discuss ways to bridge the gap between cutting-edge developments and breeders, expediting MAS-based winter wheat breeding for drought tolerance.

A diverse panel of 755 bread wheat accessions harbors untapped genetic diversity in landraces and reveals novel genetic regions conferring powdery mildew resistance

KEY MESSAGE

A bread wheat panel reveals rich genetic diversity in Turkish, Pakistani and Iranian landraces and novel resistance loci to diverse powdery mildew isolates via subsetting approaches in association studies. Wheat breeding for disease resistance relies on the availability and use of diverse genetic resources. More than 800,000 wheat accessions are globally conserved in gene banks, but they are mostly uncharacterized for the presence of resistance genes and their potential for agriculture. Based on the selective reduction of previously assembled collections for allele mining for disease resistance, we assembled a trait-customized panel of 755 geographically diverse bread wheat accessions with a focus on landraces, called the LandracePLUS panel. Population structure analysis of this panel based on the TaBW35K SNP array revealed an increased genetic diversity compared to 632 landraces genotyped in an earlier study and 17 high-quality sequenced wheat accessions. The additional genetic diversity found here mostly originated from Turkish, Iranian and Pakistani landraces. We characterized the LandracePLUS panel for resistance to ten diverse isolates of the fungal pathogen powdery mildew. Performing genome-wide association studies and dividing the panel further by a targeted subsetting approach for accessions of distinct geographical origin, we detected several known and already cloned genes, including the Pm2a gene. In addition, we identified 22 putatively novel powdery mildew resistance loci that represent useful sources for resistance breeding and for research on the mildew-wheat pathosystem. Our study shows the value of assembling trait-customized collections and utilizing a diverse range of pathogen races to detect novel loci. It further highlights the importance of integrating landraces of different geographical origins into future diversity studies.

Different wheat loci are associated to heritable free asparagine content in grain grown under different water and nitrogen availability

KEY MESSAGE

Different wheat QTLs were associated to the free asparagine content of grain grown in four different conditions. Environmental effects are a key factor when selecting for low acrylamide-forming potential. The amount of free asparagine in grain of a wheat genotype determines its potential to form harmful acrylamide in derivative food products. Here, we explored the variation in the free asparagine, aspartate, glutamine and glutamate contents of 485 accessions reflecting wheat worldwide diversity to define the genetic architecture governing the accumulation of these amino acids in grain. Accessions were grown under high and low nitrogen availability and in water-deficient and well-watered conditions, and plant and grain phenotypes were measured. Free amino acid contents of grain varied from 0.01 to 1.02 mg g-1 among genotypes in a highly heritable way that did not correlate strongly with grain yield, protein content, specific weight, thousand-kernel weight or heading date. Mean free asparagine content was 4% higher under high nitrogen and 3% higher in water-deficient conditions. After genotyping the accessions, single-locus and multi-locus genome-wide association study models were used to identify several QTLs for free asparagine content located on nine chromosomes. Each QTL was associated with a single amino acid and growing environment, and none of the QTLs colocalised with genes known to be involved in the corresponding amino acid metabolism. This suggests that free asparagine content is controlled by several loci with minor effects interacting with the environment. We conclude that breeding for reduced asparagine content is feasible, but should be firmly based on multi-environment field trials.

KEY MESSAGE

Different wheat QTLs were associated to the free asparagine content of grain grown in four different conditions. Environmental effects are a key factor when selecting for low acrylamide-forming potential.

TaRECQ4 contributes to maintain both homologous and homoeologous recombination during wheat meiosis

Introduction

Meiotic recombination (or crossover, CO) is essential for gamete fertility as well as for alleles and genes reshuffling that is at the heart of plant breeding. However, CO remains a limited event, which strongly hampers the rapid production of original and improved cultivars. RecQ4 is a gene encoding a helicase protein that, when mutated, contributes to improve recombination rate in all species where it has been evaluated so far.

Methods

In this study, we developed wheat (Triticum aestivum L.) triple mutant (TM) for the three homoeologous copies of TaRecQ4 as well as mutants for two copies and heterozygous for the last one (Htz-A, Htz-B, Htz-D).

Results

Phenotypic observation revealed a significant reduction of fertility and pollen viability in TM and Htz-B plants compared to wild type plants suggesting major defects during meiosis. Cytogenetic analyses of these plants showed that complete absence of TaRecQ4 as observed in TM plants, leads to chromosome fragmentation during the pachytene stage, resulting in problems in the segregation of chromosomes during meiosis. Htz-A and Htz-D mutants had an almost normal meiotic progression indicating that both TaRecQ4-A and TaRecQ4-D copies are functional and that there is no dosage effect for TaRecQ4 in bread wheat. On the contrary, the TaRecQ4-B copy seems knocked-out, probably because of a SNP leading to a Threonine>Alanine change at position 539 (T539A) of the protein, that occurs in the crucial helicase ATP bind/DEAD/ResIII domain which unwinds nucleic acids. Occurrence of numerous multivalents in TM plants suggests that TaRecQ4 could also play a role in the control of homoeologous recombination.

Discussion

These findings provide a foundation for further molecular investigations into wheat meiosis regulation to fully understand the underlying mechanisms of how TaRecQ4 affects chiasma formation, as well as to identify ways to mitigate these defects and enhance both homologous and homoeologous recombination efficiency in wheat.

Plant trait relationships are maintained within a major crop species: lack of artificial selection signal and potential for improved agronomic performance

The exploration of phenotypic spaces of large sets of plant species has considerably increased our understanding of diversification processes in the plant kingdom. Nevertheless, such advances have predominantly relied on interspecific comparisons that hold several limitations. Here, we grew in the field a unique set of 179 inbred lines of durum wheat, Triticum turgidum spp. durum, characterized by variable degrees of artificial selection. We measured aboveground and belowground traits as well as agronomic traits to explore the functional and agronomic trait spaces and to investigate trait-to-agronomic performance relationships. We showed that the wheat functional trait space shared commonalities with global cross-species spaces previously described, with two main axes of variation: a root foraging axis and a slow-fast trade-off axis. Moreover, we detected a clear signature of artificial selection on the variation of agronomic traits, unlike functional traits. Interestingly, we identified alternative phenotypic combinations that can optimize crop performance. Our work brings insightful knowledge about the structure of phenotypic spaces of domesticated plants and the maintenance of phenotypic trade-offs in response to artificial selection, with implications for trade-off-free and multi-criteria selection in plant breeding.

Genetic diversity, linkage disequilibrium, and population structure of tetraploid wheat landraces originating from Europe and Asia

BACKGROUND

Durum wheat is one of the most important crops, especially in the Mediterranean region. Insight into the genetic diversity of germplasm can improve the breeding program management in various traits. This study was done using single nucleotide polymorphisms (SNP) markers to characterize the genetic distinctiveness and differentiation of tetraploid wheat landraces collected from nine European and Asian countries. A sum of 23,334 polymorphic SNPs was detected in 126 tetraploid wheat landraces in relation to the reference genome.

RESULTS

The number of identified SNPs was 11,613 and 11,721 in A and B genomes, respectively. The highest and lowest diversity was on 6B and 6 A chromosomes, respectively. Structure analysis classified the landraces into two distinct subpopulations (K = 2). Evaluating the principal coordinate analysis (PCoA) and weighted pair-group method using arithmetic averages (WPGMA) clustering results demonstrated that landraces (99.2%) are categorized into one of the two chief subpopulations. Therefore, the grouping pattern did not clearly show the presence of a clear pattern of relationships between genetic diversity and their geographical derivation. Part of this result could be due to the historical exchange between different germplasms. Although the result did not separate landraces based on their region of origin, the landraces collected from Iran were classified into the same group and cluster. Analysis of molecular variance (AMOVA) also confirmed the results of population structure. Finally, Durum wheat landraces in some countries, including Turkey, Russia, Ukraine, and Afghanistan, were highly diverse, while others, including Iran and China, were low-diversity.

CONCLUSION

The recent study concluded that the 126 tetraploid wheat genotypes and their GBS-SNP markers are very appropriate for quantitative trait loci (QTLs) mapping and genome-wide association studies (GWAS). The core collection comprises two distinct subpopulations. Subpopulation II genotypes are the most diverse genotypes, and if they possess desired traits, they may be used in future breeding programs. The degree of diversity in the landraces of countries can provide the ground for the improvement of new cultivars with international cooperation. linkage disequilibrium (LD) hotspot distribution across the genome was investigated, which provides useful information about the genomic regions that contain intriguing genes.

Comparison of genomic-enabled cross selection criteria for the improvement of inbred line breeding populations

A crucial step in inbred plant breeding is the choice of mating design to derive high-performing inbred varieties while also maintaining a competitive breeding population to secure sufficient genetic gain in future generations. In practice, the mating design usually relies on crosses involving the best parental inbred lines to ensure high mean progeny performance. This excludes crosses involving lower performing but more complementary parents in terms of favorable alleles. We predicted the ability of crosses to produce putative outstanding progenies (high mean and high variance progeny distribution) using genomic prediction models. This study compared the benefits and drawbacks of 7 genomic cross selection criteria (CSC) in terms of genetic gain for 1 trait and genetic diversity in the next generation. Six CSC were already published, and we propose an improved CSC that can estimate the proportion of progeny above a threshold defined for the whole mating plan. We simulated mating designs optimized using different CSC. The 835 elite parents came from a real breeding program and were evaluated between 2000 and 2016. We applied constraints on parental contributions and genetic similarities between selected parents according to usual breeder practices. Our results showed that CSC based on progeny variance estimation increased the genetic value of superior progenies by up to 5% in the next generation compared to CSC based on the progeny mean estimation (i.e. parental genetic values) alone. It also increased the genetic gain (up to 4%) and/or maintained more genetic diversity at QTLs (up to 4% more genic variance when the marker effects were perfectly estimated).

Genetic regions determine tolerance to nitrogen deficiency in European elite bread wheats grown under contrasting nitrogen stress scenarios

KEY MESSAGE

Clustering 24 environments in four contrasting nitrogen stress scenarios enabled the detection of genetic regions determining tolerance to nitrogen deficiency in European elite bread wheats. Increasing the nitrogen use efficiency of wheat varieties is an important goal for breeding. However, most genetic studies of wheat grown at different nitrogen levels in the field report significant interactions with the genotype. The chromosomal regions possibly involved in these interactions are largely unknown. The objective of this study was to quantify the response of elite bread wheat cultivars to different nitrogen field stress scenarios and identify genomic regions involved in this response. For this purpose, 212 elite bread wheat varieties were grown in a multi-environment trial at different nitrogen levels. Genomic regions associated with grain yield, protein concentration and grain protein deviation responses to nitrogen deficiency were identified. Environments were clustered according to adjusted means for grain yield, yield components and grain protein concentration. Four nitrogen availability scenarios were identified: optimal condition, moderate early deficiency, severe late deficiency, and severe continuous deficiency. A large range of tolerance to nitrogen deficiency was observed among varieties, which were ranked differently in different nitrogen deficiency scenarios. The well-known negative correlation between grain yield and grain protein concentration also existed between their respective tolerance indices. Interestingly, the tolerance indices for grain yield and grain protein deviation were either null or weakly positive meaning that breeding for the two traits should be less difficult than expected. Twenty-two QTL regions were identified for the tolerance indices. By selecting associated markers, these regions may be selected separately or combined to improve the tolerance to N deficiency within a breeding programme.

The genetic identity of neighboring plants in intraspecific mixtures modulates disease susceptibility of both wheat and rice

Mixing crop cultivars has long been considered as a way to control epidemics at the field level and is experiencing a revival of interest in agriculture. Yet, the ability of mixing to control pests is highly variable and often unpredictable in the field. Beyond classical diversity effects such as dispersal barrier generated by genotypic diversity, several understudied processes are involved. Among them is the recently discovered neighbor-modulated susceptibility (NMS), which depicts the phenomenon that susceptibility in a given plant is affected by the presence of another healthy neighboring plant. Despite the putative tremendous importance of NMS for crop science, its occurrence and quantitative contribution to modulating susceptibility in cultivated species remains unknown. Here, in both rice and wheat inoculated in greenhouse conditions with foliar fungal pathogens considered as major threats, using more than 200 pairs of intraspecific genotype mixtures, we experimentally demonstrate the occurrence of NMS in 11% of the mixtures grown in experimental conditions that precluded any epidemics. Thus, the susceptibility of these 2 major crops results from indirect effects originating from neighboring plants. Quite remarkably, the levels of susceptibility modulated by plant-plant interactions can reach those conferred by intrinsic basal immunity. These findings open new avenues to develop more sustainable agricultural practices by engineering less susceptible crop mixtures thanks to emergent but now predictable properties of mixtures.

Multivariate Genome-Wide Association Study of Concentrations of Seven Elements in Seeds Reveals Four New Loci in Russian Wheat Lines

Wheat is a cereal grain that plays an important role in the world's food industry. The identification of the loci that change the concentration of elements in wheat seeds is an important challenge nowadays especially for genomic selection and breeding of novel varieties. In this study, we performed a multivariate genome-wide association study (GWAS) of the seven traits-concentrations of Zn, Mg, Mn, Ca, Cu, Fe, and K in grain-of the Russian collection of common wheat Triticum aestivum (N = 149 measured in two years in two different fields). We replicated one known locus associated with the concentration of Zn (IAAV1375). We identified four novel loci-BS00022069_51 (associated with concentrations of Ca and K), RFL_Contig6053_3082 (associated with concentrations of Fe and Mn), Kukri_rep_c70864_329 (associated with concentrations of all elements), and IAAV8416 (associated with concentrations of Fe and Mn)-three of them were located near the genes TraesCS6A02G375400, TraesCS7A02G094800, and TraesCS5B02G325400. Our result adds novel information on the loci involved in wheat grain element contents and may be further used in genomic selection.

Symbiotic Variations among Wheat Genotypes and Detection of Quantitative Trait Loci for Molecular Interaction with Auxin-Producing Azospirillum PGPR

Crop varieties differ in their ability to interact with Plant Growth-Promoting Rhizobacteria (PGPR), but the genetic basis for these differences is unknown. This issue was addressed with the PGPR Azospirillum baldaniorum Sp245, using 187 wheat accessions. We screened the accessions based on the seedling colonization by the PGPR and the expression of the phenylpyruvate decarboxylase gene ppdC (for synthesis of the auxin indole-3-acetic acid), using gusA fusions. Then, the effects of the PGPR on the selected accessions stimulating Sp245 (or not) were compared in soil under stress. Finally, a genome-wide association approach was implemented to identify the quantitative trait loci (QTL) associated with PGPR interaction. Overall, the ancient genotypes were more effective than the modern genotypes for Azospirillum root colonization and ppdC expression. In non-sterile soil, A. baldaniorum Sp245 improved wheat performance for three of the four PGPR-stimulating genotypes and none of the four non-PGPR-stimulating genotypes. The genome-wide association did not identify any region for root colonization but revealed 22 regions spread on 11 wheat chromosomes for ppdC expression and/or ppdC induction rate. This is the first QTL study focusing on molecular interaction with PGPR bacteria. The molecular markers identified provide the possibility to improve the capacity of modern wheat genotypes to interact with Sp245, as well as, potentially, other Azospirillum strains.

A framework for improving wheat spike development and yield based on the master regulatory TOR and SnRK gene systems

The low rates of yield gain in wheat breeding programs create an ominous situation for the world. Amongst the reasons for this low rate are issues manifested in spike development that result in too few spikelets, fertile florets, and therefore grains being produced. Phases in spike development are particularly sensitive to stresses of various kinds and origins, and these are partly responsible for the deficiencies in grain production and slow rates of gain in yield. The diversity of developmental processes, stresses, and the large numbers of genes involved make it particularly difficult to prioritize approaches in breeding programs without an overarching, mechanistic framework. Such a framework, introduced here, is provided around the master regulator target of rapamycin and sucrose non-fermenting-1 (SNF1)-related protein kinase complexes and their control by trehalose-6-phosphate and other molecules. Being master regulators of the balance between growth and growth inhibition under stress, these provide genetic targets for creating breakthroughs in yield enhancement. Examples of potential targets and experimental approaches are described.

Wheat Omics: Advancements and Opportunities

Plant omics, which includes genomics, transcriptomics, metabolomics and proteomics, has played a remarkable role in the discovery of new genes and biomolecules that can be deployed for crop improvement. In wheat, great insights have been gleaned from the utilization of diverse omics approaches for both qualitative and quantitative traits. Especially, a combination of omics approaches has led to significant advances in gene discovery and pathway investigations and in deciphering the essential components of stress responses and yields. Recently, a Wheat Omics database has been developed for wheat which could be used by scientists for further accelerating functional genomics studies. In this review, we have discussed various omics technologies and platforms that have been used in wheat to enhance the understanding of the stress biology of the crop and the molecular mechanisms underlying stress tolerance.

Diversity matters in wheat mixtures: A genomic survey of the impact of genetic diversity on the performance of 12 way durum wheat mixtures grown in two contrasted and controlled environments

In ecology, an increase in genetic diversity within a community in natural ecosystems increases its productivity, while in evolutionary biology, kinship selection predicts that relatedness on social traits improves fitness. Varietal mixtures, where different genotypes are grown together, show contrasting results, especially for grain yield where both positive and negative effects of mixtures have been reported. To understand the effect of diversity on field performance, we grew 96 independent mixtures each composed with 12 durum wheat (Triticum turgidum ssp. durum Thell.) inbred lines, under two contrasting environmental conditions for water availability. Using dense genotyping, we imputed allelic frequencies and a genetic diversity index on more than 96000 loci for each mixture. We then analyzed the effect of genetic diversity on agronomic performance using a genome-wide approach. We explored the stress gradient hypothesis, which proposes that the greater the unfavourable conditions, the more beneficial the effect of diversity on mixture performance. We found that diversity on average had a negative effect on yield and its components while it was beneficial on grain weight. There was little support for the stress gradient theory. We discuss how to use genomic data to improve the assembly of varietal mixtures.

Multiomics reveals an essential role of long-distance translocation in regulating plant cadmium resistance and grain accumulation in allohexaploid wheat (Triticum aestivum)

Cadmium (Cd) is a highly toxic heavy metal that readily enters cereals, such as wheat, via the roots and is translocated to the shoots and grains, thereby posing high risks to human health. However, the vast and complex genome of allohexaploid wheat makes it challenging to understand Cd resistance and accumulation. In this study, a Cd-resistant cultivar of wheat, 'ZM1860', and a Cd-sensitive cultivar, 'ZM32', selected from a panel of 442 accessions, exhibited significantly different plant resistance and grain accumulation. We performed an integrated comparative analysis of the morpho-physiological traits, ionomic and phytohormone profiles, genomic variations, transcriptomic landscapes, and gene functionality in order to identify the mechanisms underlying these differences. Under Cd toxicity, 'ZM1860' outperformed 'ZM32', which showed more severe leaf chlorosis, poorer root architecture, higher accumulation of reactive oxygen species, and disordered phytohormone homeostasis. Ionomics showed that 'ZM32' had a higher root-to-shoot translocation coefficient of Cd and accumulated more Cd in the grains than 'ZM1860'. Whole-genome re-sequencing (WGS) and transcriptome sequencing identified numerous DNA variants and differentially expressed genes involved in abiotic stress responses and ion transport between the two genotypes. Combined ionomics, transcriptomics, and functional gene analysis identified the plasma membrane-localized heavy metal ATPase TaHMA2b-7A as a crucial Cd exporter regulating long-distance Cd translocation in wheat. WGS- and PCR-based analysis of sequence polymorphisms revealed a 25-bp InDel site in the promoter region of TaHMA2b-7A, and this was probably responsible for the differential expression. Our multiomics approach thus enabled the identification of a core transporter involved in long-distance Cd translocation in wheat, and it may provide an elite genetic resource for improving plant Cd resistance and reducing grain Cd accumulation in wheat and other cereal crops.

Bread wheat (T. aestivum) variability: Phenotypic and genotypic data from 75 varieties

Most bread wheat is consumed after processing, which mainly depends on the quantity and quality of protein in the grain. Storage protein content and composition particularly influence the end use quality of milled grain products. Storage proteins are components of the gluten network that confer dough viscoelasticity, an essential property for processing. To explore grain storage protein diversity, 75 bread wheat accessions were grown with two replicates each at two locations. Grains were harvested at maturity and samples were phenotyped for each site and each replicate plant. Grain hardness, thousand-kernel weight and grain nitrogen content were measured. The protein composition of flour from each replicate was characterised by reverse phase-high performance liquid chromatography (RP-HPLC). The molecular distribution of flour polymers was determined by asymmetric flow field-flow fractionation (AF4) and dough technological properties were assessed using a Glutomatic system and a Chopin alveograph. In addition, the 75 accessions were genotyped by the BreedWheat 35k genotyping array (Axiom TaBW35K) containing 34,746 single nucleotide polymorphism markers (SNPs). The dataset produced by this work includes six files with raw data, two files with protocols and figures. Data show the genotypic and phenotypic variabilities of the material used and can be used to explore genetic and environmental effects on traits involved in grain protein quality. This dataset is associated to the research article "Differences in bread protein digestibility traced to wheat cultivar traits" [1].

Discovery of DNA polymorphisms via genome-resequencing and development of molecular markers between two barley cultivars

Genome resequencing uncovers genome-wide DNA polymorphisms that are useful for the development of high-density InDel markers between two barley cultivars. Discovering genomic variations and developing genetic markers are crucial for genetics studies and molecular breeding in cereal crops. Although InDels (insertions and deletions) have become popular because of their abundance and ease of detection, discovery of genome-wide DNA polymorphisms and development of InDel markers in barley have lagged behind other cereal crops such as rice, maize and wheat. In this study, we re-sequenced two barley cultivars, Golden Promise (GP, a classic British spring barley variety) and Hua30 (a Chinese spring barley variety), and mapped clean reads to the reference Morex genome, and identified in total 13,933,145 single nucleotide polymorphisms (SNPs) and 1,240,456 InDels for GP with Morex, 11,297,100 SNPs and 781,687 InDels for Hua30 with Morex, and 13,742,399 SNPs and 1,191,597 InDels for GP with Hua30. We further characterized distinct types, chromosomal distribution patterns, genome location, functional effect, and other features of these DNA polymorphisms. Additionally, we revealed the functional relevance of these identified SNPs/InDels regarding different flowering times between Hua30 and GP within 17 flowering time genes. Furthermore, we developed a series of InDel markers and validated them experimentally in 43 barley core accessions, respectively. Finally, we rebuilt population structure and phylogenetic tree of these 43 barley core accessions. Collectively, all of these genetic resources will facilitate not only the basic research but also applied research in barley.

CGIAR Barley Breeding Toolbox: A diversity panel to facilitate breeding and genomic research in the developing world

Breeding programs in developing countries still cannot afford the new genotyping technologies, hindering their research. We aimed to assemble an Association Mapping panel to serve as CGIAR Barley Breeding Toolbox (CBBT), especially for the Developing World. The germplasm had to be representative of the one grown in the Developing World; with high genetic variability and be of public domain. For it, we genotyped with the Infinium iSelect 50K chip, a Global Barley Panel (GBP) of 530 genotypes representing a wide range of row-types, end-uses, growth habits, geographical origins and environments. 40,342 markers were polymorphic with an average polymorphism information content of 0.35 and 66% of them exceeding 0.25. The analysis of the population structure identified 8 subpopulations mostly linked to geographical origin, four of them with significant ICARDA origin. The 16 allele combinations at 4 major flowering genes (HvVRN-H3, HvPPD-H1, HvVRN-H1 and HvCEN) explained 11.07% genetic variation and were linked to the geographic origins of the lines. ICARDA material showed the widest diversity as revealed by the highest number of polymorphic loci (99.76% of all polymorphic SNPs in GBP), number of private alleles and the fact that ICARDA lines were present in all 8 subpopulations and carried all 16 allelic combinations. Due to their genetic diversity and their representativity of the germplasm adapted to the Developing World, ICARDA-derived lines and cultivated landraces were pre-selected to form the CBBT. Using the Mean of Transformed Kinships method, we assembled a panel capturing most of the allelic diversity in the GBP. The CBBT (N=250) preserves good balance between row-types and good representation of both phenology allelic combinations and subpopulations of the GBP. The CBBT and its genotypic data is available to researchers worldwide as a collaborative tool to underpin the genetic mechanisms of traits of interest for barley cultivation.

A high-throughput skim-sequencing approach for genotyping, dosage estimation and identifying translocations

The development of next-generation sequencing (NGS) enabled a shift from array-based genotyping to directly sequencing genomic libraries for high-throughput genotyping. Even though whole-genome sequencing was initially too costly for routine analysis in large populations such as breeding or genetic studies, continued advancements in genome sequencing and bioinformatics have provided the opportunity to capitalize on whole-genome information. As new sequencing platforms can routinely provide high-quality sequencing data for sufficient genome coverage to genotype various breeding populations, a limitation comes in the time and cost of library construction when multiplexing a large number of samples. Here we describe a high-throughput whole-genome skim-sequencing (skim-seq) approach that can be utilized for a broad range of genotyping and genomic characterization. Using optimized low-volume Illumina Nextera chemistry, we developed a skim-seq method and combined up to 960 samples in one multiplex library using dual index barcoding. With the dual-index barcoding, the number of samples for multiplexing can be adjusted depending on the amount of data required, and could be extended to 3,072 samples or more. Panels of doubled haploid wheat lines (Triticum aestivum, CDC Stanley x CDC Landmark), wheat-barley (T. aestivum x Hordeum vulgare) and wheat-wheatgrass (Triticum durum x Thinopyrum intermedium) introgression lines as well as known monosomic wheat stocks were genotyped using the skim-seq approach. Bioinformatics pipelines were developed for various applications where sequencing coverage ranged from 1 × down to 0.01 × per sample. Using reference genomes, we detected chromosome dosage, identified aneuploidy, and karyotyped introgression lines from the skim-seq data. Leveraging the recent advancements in genome sequencing, skim-seq provides an effective and low-cost tool for routine genotyping and genetic analysis, which can track and identify introgressions and genomic regions of interest in genetics research and applied breeding programs.

Fine mapping of a stripe rust resistance gene YrZM175 in bread wheat

A stripe rust resistance gene YrZM175 in Chinese wheat cultivar Zhongmai 175 was mapped to a genomic interval of 636.4 kb on chromosome arm 2AL, and a candidate gene was predicted. Stripe rust, caused by Puccinia striiformis f. sp. tritici (PST), is a worldwide wheat disease that causes large losses in production. Fine mapping and cloning of resistance genes are important for accurate marker-assisted breeding. Here, we report the fine mapping and candidate gene analysis of stripe rust resistance gene YrZM175 in a Chinese wheat cultivar Zhongmai 175. Fifteen F₁, 7,325 F₂ plants and 117 F_2:3 lines derived from cross Avocet S/Zhongmai 175 were inoculated with PST race CYR32 at the seedling stage in a greenhouse, and F_2:3 lines were also evaluated for stripe rust reaction in the field using mixed PST races. Bulked segregant RNA-seq (BSR-seq) analyses revealed 13 SNPs in the region 762.50-768.52 Mb on chromosome arm 2AL. By genome mining, we identified SNPs and InDels between the parents and contrasting bulks and mapped YrZM175 to a 0.72-cM, 636.4-kb interval spanned by YrZM175-InD1 and YrZM175-InD2 (763,452,916-764,089,317 bp) including two putative disease resistance genes based on IWGSC RefSeq v1.0. Collinearity analysis indicated similar target genomic intervals in Chinese Spring, Aegilops tauschii (2D: 647.7-650.5 Mb), Triticum urartu (2A: 750.7-752.3 Mb), Triticum dicoccoides (2A: 771.0-774.5 Mb), Triticum turgidum (2B: 784.7-788.2 Mb), and Triticum aestivum cv. Aikang 58 (2A: 776.3-778.9 Mb) and Jagger (2A: 789.3-791.7 Mb). Through collinearity analysis, sequence alignments of resistant and susceptible parents and gene expression level analysis, we predicted TRITD2Bv1G264480 from Triticum turgidum to be a candidate gene for map-based cloning of YrZM175. A gene-specific marker for TRITD2Bv1G264480 co-segregated with the resistance gene. Molecular marker analysis and stripe rust response data revealed that YrZM175 was different from genes Yr1, Yr17, Yr32, and YrJ22 located on chromosome 2A. Fine mapping of YrZM175 lays a solid foundation for functional gene analysis and marker-assisted selection for improved stripe rust resistance in wheat.

Fine mapping of a stripe rust resistance gene YrZM175 in bread wheat

A stripe rust resistance gene YrZM175 in Chinese wheat cultivar Zhongmai 175 was mapped to a genomic interval of 636.4 kb on chromosome arm 2AL, and a candidate gene was predicted. Stripe rust, caused by Puccinia striiformis f. sp. tritici (PST), is a worldwide wheat disease that causes large losses in production. Fine mapping and cloning of resistance genes are important for accurate marker-assisted breeding. Here, we report the fine mapping and candidate gene analysis of stripe rust resistance gene YrZM175 in a Chinese wheat cultivar Zhongmai 175. Fifteen F₁, 7,325 F₂ plants and 117 F_2:3 lines derived from cross Avocet S/Zhongmai 175 were inoculated with PST race CYR32 at the seedling stage in a greenhouse, and F_2:3 lines were also evaluated for stripe rust reaction in the field using mixed PST races. Bulked segregant RNA-seq (BSR-seq) analyses revealed 13 SNPs in the region 762.50-768.52 Mb on chromosome arm 2AL. By genome mining, we identified SNPs and InDels between the parents and contrasting bulks and mapped YrZM175 to a 0.72-cM, 636.4-kb interval spanned by YrZM175-InD1 and YrZM175-InD2 (763,452,916-764,089,317 bp) including two putative disease resistance genes based on IWGSC RefSeq v1.0. Collinearity analysis indicated similar target genomic intervals in Chinese Spring, Aegilops tauschii (2D: 647.7-650.5 Mb), Triticum urartu (2A: 750.7-752.3 Mb), Triticum dicoccoides (2A: 771.0-774.5 Mb), Triticum turgidum (2B: 784.7-788.2 Mb), and Triticum aestivum cv. Aikang 58 (2A: 776.3-778.9 Mb) and Jagger (2A: 789.3-791.7 Mb). Through collinearity analysis, sequence alignments of resistant and susceptible parents and gene expression level analysis, we predicted TRITD2Bv1G264480 from Triticum turgidum to be a candidate gene for map-based cloning of YrZM175. A gene-specific marker for TRITD2Bv1G264480 co-segregated with the resistance gene. Molecular marker analysis and stripe rust response data revealed that YrZM175 was different from genes Yr1, Yr17, Yr32, and YrJ22 located on chromosome 2A. Fine mapping of YrZM175 lays a solid foundation for functional gene analysis and marker-assisted selection for improved stripe rust resistance in wheat.

Phenomic selection in wheat breeding: prediction of the genotype-by-environment interaction in multi-environment breeding trials

Phenomic prediction of wheat grain yield and heading date in different multi-environmental trial scenarios is accurate. Modelling the genotype-by-environment interaction effect using phenomic data is a potentially low-cost complement to genomic prediction. The performance of wheat cultivars in multi-environmental trials (MET) is difficult to predict because of the genotype-by-environment interactions (G × E). Phenomic selection is supposed to be efficient for modelling the G × E effect because it accounts for non-additive effects. Here, phenomic data are near-infrared (NIR) spectra obtained from plant material. While phenomic selection has recently been shown to accurately predict wheat grain yield in single environments, its accuracy needs to be investigated for MET. We used four datasets from two winter wheat breeding programs to test and compare the predictive abilities of phenomic and genomic models for grain yield and heading date in different MET scenarios. We also compared different methods to model the G × E using different covariance matrices based on spectra. On average, phenomic and genomic prediction abilities are similar in all different MET scenarios. Better predictive abilities were obtained when G × E effects were modelled with NIR spectra than without them, and it was better to use all the spectra of all genotypes in all environments for modelling the G × E. To facilitate the implementation of phenomic prediction, we tested MET designs where the NIR spectra were measured only on the genotype-environment combinations phenotyped for the target trait. Missing spectra were predicted with a weighted multivariate ridge regression. Intermediate predictive abilities for grain yield were obtained in a sparse testing scenario and for new genotypes, which shows that phenomic selection is an efficient and practicable prediction method for dealing with G × E.

Exploring the Genetic Diversity and Population Structure of Wheat Landrace Population Conserved at ICARDA Genebank

Landraces are considered a valuable source of potential genetic diversity that could be used in the selection process in any plant breeding program. Here, we assembled a population of 600 bread wheat landraces collected from eight different countries, conserved at the ICARDA's genebank, and evaluated the genetic diversity and the population structure of the landraces using single nucleotide polymorphism (SNP) markers. A total of 11,830 high-quality SNPs distributed across the genomes A (40.5%), B (45.9%), and D (13.6%) were used for the final analysis. The population structure analysis was evaluated using the model-based method (STRUCTURE) and distance-based methods [discriminant analysis of principal components (DAPC) and principal component analysis (PCA)]. The STRUCTURE method grouped the landraces into two major clusters, with the landraces from Syria and Turkey forming two clusters with high proportions of admixture, whereas the DAPC and PCA analysis grouped the population into three subpopulations mostly according to the geographical information of the landraces, i.e., Syria, Iran, and Turkey with admixture. The analysis of molecular variance revealed that the majority of the variation was due to genetic differences within the populations as compared with between subpopulations, and it was the same for both the cluster-based and distance-based methods. Genetic distance analysis was also studied to estimate the differences between the landraces from different countries, and it was observed that the maximum genetic distance (0.389) was between the landraces from Spain and Palestine, whereas the minimum genetic distance (0.013) was observed between the landraces from Syria and Turkey. It was concluded from the study that the model-based methods (DAPC and PCA) could dissect the population structure more precisely when compared with the STRUCTURE method. The population structure and genetic diversity analysis of the bread wheat landraces presented here highlight the complex genetic architecture of the landraces native to the Fertile Crescent region. The results of this study provide useful information for the genetic improvement of hexaploid wheat and facilitate the use of landraces in wheat breeding programs.

QTL Mapping of Adult Plant Resistance to Stripe Rust in a Doubled Haploid Wheat Population

Stripe rust caused by Puccinia striiformis Westend. f. sp. tritici. is a major bread wheat disease worldwide with yield losses of up to 100% under severe disease pressure. The deployment of resistant cultivars with adult plant resistance to the disease provides a long-term solution to stripe rust of wheat. An advanced line from the International Winter Wheat Improvement Program (IWWIP) 130675 (Avd/Vee#1//1-27-6275/Cf 1770/3/MV171-C-17466) showed a high level of adult plant resistance to stripe rust in the field. To identify the adult plant resistance genes in this elite line, a mapping population of 190 doubled haploid (DH) lines was developed from a cross between line 130675 and the universal stripe rust-susceptible variety Avocet S. The DH population was evaluated at precision wheat stripe rust phenotyping platform, in Izmir during 2019, 2020, and 2021 cropping seasons under artificial inoculations. Composite interval mapping (CIM) identified two stable QTLs QYr.rcrrc-3B.1, and QYr.rcrrc-3B.2, which were detected in multiple years. In addition to these two QTLs, five more QTLs, QYr.rcrrc-1B, QYr.rcrrc-2A, QYr.rcrrc-3A, QYr.rcrrc-5A, and QYr.rcrrc-7D, were identified, which were specific to the cropping year (environment). All QTLs were derived from the resistant parent, except QYr.rcrrc-3A. The significant QTLs explained 3.4-20.6% of the phenotypic variance. SNP markers flanking the QTL regions can be amenable to marker-assisted selection. The best DH lines with high yield, end-use quality, and stripe rust resistance can be used for further selection for improved germplasm. SNP markers flanking the QTL regions can aid in identifying such lines.

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.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

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

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

Genome-wide Association Study for Starch Pasting Properties in Chinese Spring Wheat

In order to understand the genetic basis of starch pasting viscosity characteristics of Chinese spring wheat, we assessed the genetic variation of RVA parameters determined by the Rapid Visco Analyser in a panel of 192 Chinese spring wheat accessions grown in Er’shi, Shihezi and Zhaosu during 2012 and 2013 cropping seasons. A genome-wide association study with 47,362 single nucleotide polymorphism (SNP) markers was conducted to detect marker-trait associations using mixed linear model. Phenotypic variations of RVA parameters ranged from 1.6 to 30.7% and broad-sense heritabilities ranged from 0.62 to 0.91. Forty-one SNP markers at 25 loci were significantly associated with seven RVA traits in at least two environments; among these, 20 SNPs were located in coding sequences (CDS) of 18 annotation genes, which can lead to discovering novel genes underpinning starch gelatinization in spring wheat. Haplotype analysis revealed one block for breakdown (BD) on chromosome 3B and two blocks for pasting temperature (T) on chromosome 7B. Cultivars with superior haplotypes at these loci showed better starch pasting viscosity than the average of all cultivars surveyed. The identified loci and associated markers provide valuable sources for future functional characterization and genetic improvement of starch quality in wheat.

Genetic Analysis of Platform-Phenotyped Root System Architecture of Bread and Durum Wheat in Relation to Agronomic Traits

Roots are essential for water and nutrient uptake but are rarely the direct target of breeding efforts. To characterize the genetic variability of wheat root architecture, the root and shoot traits of 200 durum and 715 bread wheat varieties were measured at a young stage on a high-throughput phenotyping platform. Heritability of platform traits ranged from 0.40 for root biomass in durum wheat to 0.82 for the number of tillers. Field phenotyping data for yield components and SNP genotyping were already available for all the genotypes. Taking differences in earliness into account, several significant correlations between root traits and field agronomic performances were found, suggesting that plants investing more resources in roots in some stressed environments favored water and nutrient uptake, with improved wheat yield. We identified 100 quantitative trait locus (QTLs) of root traits in the bread wheat panels and 34 in the durum wheat panel. Most colocalized with QTLs of traits measured in field conditions, including yield components and earliness for bread wheat, but only in a few environments. Stress and climatic indicators explained the differential effect of some platform QTLs on yield, which was positive, null, or negative depending on the environmental conditions. Modern breeding has led to deeper rooting but fewer seminal roots in bread wheat. The number of tillers has been increased in bread wheat, but decreased in durum wheat, and while the root-shoot ratio for bread wheat has remained stable, for durum wheat it has been increased. Breeding for root traits or designing ideotypes might help to maintain current yield while adapting to specific drought scenarios.

From cultivar mixtures to allelic mixtures: opposite effects of allelic richness between genotypes and genotype richness in wheat

Agroecosystem diversification through increased crop genetic diversity could provide multiple services such as improved disease control or increased productivity. However, we still poorly understand how genetic diversity affects agronomic performance. We grew 179 inbred lines of durum wheat in pure stands and in 202 binary mixtures in field conditions. We then tested the effect of allelic richness between genotypes and genotype richness on grain yield and Septoria tritici blotch disease. Allelic richness was tested at 19K single nucleotide polymorphisms distributed along the durum wheat genome. Both genotype richness and allelic richness could be equal to 1 or 2. Mixtures were overall more productive and less diseased than their pure stand components. Yet, we identified one locus at which allelic richness between genotypes was associated with increased disease severity and decreased grain yield. The effect of allelic richness at this locus was stronger than the effect of genotype richness on grain yield (-7.6% vs +5.7%). Our results suggest that positive effects of crop diversity can be reversed by unfavourable allelic associations. This highlights the need to integrate genomic data into crop diversification strategies. More generally, investigating plant-plant interactions at the genomic level is promising to better understand biodiversity-ecosystem functioning relationships.

Identification of QTLs affecting post-anthesis heat stress responses in European bread wheat

The response of a large panel of European elite wheat varieties to post-anthesis heat stress is influenced by 17 QTL linked to grain weight or the stay-green phenotype. Heat stress is a critical abiotic stress for winter bread wheat (Triticum aestivum L.) especially at the flowering and grain filling stages, limiting its growth and productivity in Europe and elsewhere. The breeding of new high-yield and stress-tolerant wheat varieties requires improved understanding of the physiological and genetic bases of heat tolerance. To identify genomic areas associated with plant and grain characteristics under heat stress, a panel of elite European wheat varieties (N = 199) was evaluated under controlled conditions in 2016 and 2017. A split-plot design was used to test the effects of high temperature for ten days after flowering. Flowering time, leaf chlorophyll content, the number of productive spikes, grain number, grain weight and grain size were measured, and the senescence process was modeled. Using genotyping data from a 280 K SNP chip, a genome-wide association study was carried out to test the main effect of each SNP and the effect of SNP × treatment interaction. Genotype × treatment interactions were mainly observed for grain traits measured on the main shoots and tillers. We identified 10 QTLs associated with the main effect of at least one trait and seven QTLs associated with the response to post-anthesis heat stress. Of these, two main QTLs associated with the heat tolerance of thousand-kernel weight were identified on chromosomes 4B and 6B. These QTLs will be useful for breeders to improve grain yield in environments where terminal heat stress is likely to occur.

Phenomic selection in wheat breeding: identification and optimisation of factors influencing prediction accuracy and comparison to genomic selection

Phenomic selection is a promising alternative or complement to genomic selection in wheat breeding. Models combining spectra from different environments maximise the predictive ability of grain yield and heading date of wheat breeding lines. Phenomic selection (PS) is a recent breeding approach similar to genomic selection (GS) except that genotyping is replaced by near-infrared (NIR) spectroscopy. PS can potentially account for non-additive effects and has the major advantage of being low cost and high throughput. Factors influencing GS predictive abilities have been intensively studied, but little is known about PS. We tested and compared the abilities of PS and GS to predict grain yield and heading date from several datasets of bread wheat lines corresponding to the first or second years of trial evaluation from two breeding companies and one research institute in France. We evaluated several factors affecting PS predictive abilities including the possibility of combining spectra collected in different environments. A simple H-BLUP model predicted both traits with prediction ability from 0.26 to 0.62 and with an efficient computation time. Our results showed that the environments in which lines are grown had a crucial impact on predictive ability based on the spectra acquired and was specific to the trait considered. Models combining NIR spectra from different environments were the best PS models and were at least as accurate as GS in most of the datasets. Furthermore, a GH-BLUP model combining genotyping and NIR spectra was the best model of all (prediction ability from 0.31 to 0.73). We demonstrated also that as for GS, the size and the composition of the training set have a crucial impact on predictive ability. PS could therefore replace or complement GS for efficient wheat breeding programs.

Whole genome re-sequencing reveals the genetic diversity and evolutionary patterns of Eucommia ulmoides

Eucommia ulmoides (E. ulmoides) is a deciduous perennial tree belonging to the order Garryales, and is known as "living fossil" plant, along with ginkgo (Ginkgo biloba), metaspaca (Metasequoia glyptostroboides) and dove tree (Davidia involucrata Baill). However, the genetic diversity and population structure of E. ulmoides are still  ambiguous nowdays. In this study, we re-sequenced the genomes of 12 E. ulmoides accessions from different major climatic geography regions in China to elucidate the genetic diversity, population structure and evolutionary pattern. By integration of phylogenetic analysis, principal component analysis and population structure analysis based on a number of high-quality SNPs, a total of 12 E. ulmoides accessions were clustered into four different groups. This result is consistent with their geographical location except for group samples from Shanghai and Hunan province. E. ulmoides accessions from Hunan province exhibited a closer genetic relationship with E. ulmoides accessions from Shanghai in China compared with other regions, which is also supported by the result of population structure analyses. Genetic diversity analysis further revealed that E. ulmoides samples in Shanghai and Hunan province were with higher genetic diversity than those in other regions in this study. In addition, we treated the E. ulmoides materials from Shanghai and Hunan province as group A, and the other materials from other places as group B, and then analyzed the evolutionary pattern of E. ulmoides. The result showed the significant differentiation (Fst = 0.1545) between group A and group B. Some candidate highly divergent genome regions were identified in group A by selective sweep analyses, and the function analysis of candidate genes in these regions showed that biological regulation processes could be correlated with the Eu-rubber biosynthesis. Notably, nine genes were identified from selective sweep regions. They were involved in the Eu-rubber biosynthesis and expressed in rubber containing tissues. The genetic diversity research and evolution model of E. ulmoides were preliminarily explored in this study, which laid the foundation for the protection of germplasm resources and the development and utilization of multipurpose germplasm resources in the future.

A genome‐wide analysis suggests pleiotropic effects of Green Revolution genes on shade avoidance in wheat

A classic example of phenotypic plasticity in plants is the suit of phenotypic responses induced by a change in the ratio of red to far-red light (R∶FR) as a result of shading, also known as the shade avoidance syndrome (SAS). While the adaptive consequences of this syndrome have been extensively discussed in natural ecosystems, how SAS varies within crop populations and how SAS evolved during crop domestication and breeding remain poorly known. In this study, we grew a panel of 180 durum wheat (Triticum turgidum ssp. durum) genotypes spanning diversity from wild, early domesticated, and elite genetic compartments under two light treatments: low R:FR light (shaded treatment) and high R:FR light (unshaded treatment). We first quantified the genetic variability of SAS, here measured as a change in plant height at the seedling stage. We then dissected the genetic basis of this variation through genome-wide association mapping. Genotypes grown in shaded conditions were taller than those grown under unshaded conditions. Interaction between light quality and genotype did not affect plant height. We found six QTLs affecting plant height. Three significantly interacted with light quality among which the well-known Rht1 gene introgressed in elite germplasm during the Green Revolution. Interestingly at three loci, short genotypes systematically expressed reduced SAS, suggesting a positive genetic correlation between plant height and plant height plasticity. Overall, our study sheds light on the evolutionary history of crops and illustrates the relevance of genetic approaches to tackle agricultural challenges.

Breeding for Economically and Environmentally Sustainable Wheat Varieties: An Integrated Approach from Genomics to Selection

There is currently a strong societal demand for sustainability, quality, and safety in bread wheat production. To address these challenges, new and innovative knowledge, resources, tools, and methods to facilitate breeding are needed. This starts with the development of high throughput genomic tools including single nucleotide polymorphism (SNP) arrays, high density molecular marker maps, and full genome sequences. Such powerful tools are essential to perform genome-wide association studies (GWAS), to implement genomic and phenomic selection, and to characterize the worldwide diversity. This is also useful to breeders to broaden the genetic basis of elite varieties through the introduction of novel sources of genetic diversity. Improvement in varieties particularly relies on the detection of genomic regions involved in agronomical traits including tolerance to biotic (diseases and pests) and abiotic (drought, nutrient deficiency, high temperature) stresses. When enough resolution is achieved, this can result in the identification of candidate genes that could further be characterized to identify relevant alleles. Breeding must also now be approached through in silico modeling to simulate plant development, investigate genotype × environment interactions, and introduce marker-trait linkage information in the models to better implement genomic selection. Breeders must be aware of new developments and the information must be made available to the world wheat community to develop new high-yielding varieties that can meet the challenge of higher wheat production in a sustainable and fluctuating agricultural context. In this review, we compiled all knowledge and tools produced during the BREEDWHEAT project to show how they may contribute to face this challenge in the coming years.

Generation of Deletion Lines in Allohexaploid Bread Wheat

Positional cloning in bread wheat (Triticum aestivum L.) remains a daunting task because of its large genome, high density of repeats, low recombination rate especially in pericentromeric regions and its allopolyploidy. One way to face this challenge is to decrease the size of the interval bearing the gene of interest both genetically and physically, in order to reduce significantly the number of potential candidate genes. In this chapter, we describe a technical approach to produce chromosome-specific deletion lines to locate precisely genes of interest onto wheat chromosomes, a step forward to their cloning.

Analyzing the Economic Effectiveness of Genomic Selection Relative to Conventional Breeding Approaches

Comparing the economic efficiency of alternative strategies for breeding requires to compare the genetic gain obtained with breeding schemes that represent the same total investment. In this chapter, we present a generic method to assess this economic efficiency for alternative breeding schemes. After presenting the baseline framework and the necessity of comparing breeding schemes with equivalent total investment, we propose one illustrative example on wheat breeding. In this application, we compare the use of conventional breeding and genomic selection. With this example, we explain the requirements and the different steps to implement this method. At last, we discuss several extensions of the baseline model.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Scripting Analyses of Genomes in Ensembl Plants

Ensembl Plants ( http://plants.ensembl.org ) offers genome-scale information for plants, with four releases per year. As of release 47 (April 2020) it features 79 species and includes genome sequence, gene models, and functional annotation. Comparative analyses help reconstruct the evolutionary history of gene families, genomes, and components of polyploid genomes. Some species have gene expression baseline reports or variation across genotypes. While the data can be accessed through the Ensembl genome browser, here we review specifically how our plant genomes can be interrogated programmatically and the data downloaded in bulk. These access routes are generally consistent across Ensembl for other non-plant species, including plant pathogens, pests, and pollinators.

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

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

Novel Design of Imputation-Enabled SNP Arrays for Breeding and Research Applications Supporting Multi-Species Hybridization

Array-based single nucleotide polymorphism (SNP) genotyping platforms have low genotype error and missing data rates compared to genotyping-by-sequencing technologies. However, design decisions used to create array-based SNP genotyping assays for both research and breeding applications are critical to their success. We describe a novel approach applicable to any animal or plant species for the design of cost-effective imputation-enabled SNP genotyping arrays with broad utility and demonstrate its application through the development of the Illumina Infinium Wheat Barley 40K SNP array Version 1.0. We show that the approach delivers high quality and high resolution data for wheat and barley, including when samples are jointly hybridised. The new array aims to maximally capture haplotypic diversity in globally diverse wheat and barley germplasm while minimizing ascertainment bias. Comprising mostly biallelic markers that were designed to be species-specific and single-copy, the array permits highly accurate imputation in diverse germplasm to improve the statistical power of genome-wide association studies (GWAS) and genomic selection. The SNP content captures tetraploid wheat (A- and B-genome) and Aegilops tauschii Coss. (D-genome) diversity and delineates synthetic and tetraploid wheat from other wheat, as well as tetraploid species and subgroups. The content includes SNP tagging key trait loci in wheat and barley, as well as direct connections to other genotyping platforms and legacy datasets. The utility of the array is enhanced through the web-based tool, Pretzel (https://plantinformatics.io/) which enables the content of the array to be visualized and interrogated interactively in the context of numerous genetic and genomic resources to be connected more seamlessly to research and breeding. The array is available for use by the international wheat and barley community.

Meta-QTLs and candidate genes for stripe rust resistance in wheat

In bread wheat, meta-QTL analysis was conducted using 353 QTLs that were available from earlier studies. When projected onto a dense consensus map comprising 76,753 markers, only 184 QTLs with the required information, could be utilized leading to identification of 61 MQTLs spread over 18 of the 21 chromosomes (barring 5D, 6D and 7D). The range for mean R2 (PVE %) was 1.9% to 48.1%, and that of CI was 0.02 to 11.47 cM; these CIs also carried 37 Yr genes. Using these MQTLs, 385 candidate genes (CGs) were also identified. Out of these CGs, 241 encoded known R proteins and 120 showed differential expression due to stripe rust infection at the seedling stage; the remaining 24 CGs were common in the sense that they encoded R proteins as well as showed differential expression. The proteins encoded by CGs carried the following widely known domains: NBS-LRR domain, WRKY domains, ankyrin repeat domains, sugar transport domains, etc. Thirteen breeders' MQTLs (PVE > 20%) including four pairs of closely linked MQTLs are recommended for use in wheat molecular breeding, for future studies to understand the molecular mechanism of stripe rust resistance and for gene cloning.

Dissecting Bread Wheat Heterosis through the Integration of Agronomic and Physiological Traits

Simple Summary

To meet the challenge of feeding almost 10 billion people by 2050, wheat yield has to double by 2050. However, over the past 20 years, yield increase has slowed down and even stagnated in the main producing countries. Similar to what has been observed in maize, hybrids have been suggested as a solution to overcome yield stagnation in wheat. However, wheat heterosis, i.e., the fact that a progeny surpasses the performances of its parents, is still limited and poorly understood. To better characterize this phenomenon, we developed and phenotyped for physiological and agronomic traits 91 hybrids and their nineteen female and sixteen male parents. We showed that hybrids had a longer grain filling phase that led to bigger grains and an increased thousand kernel weight. This resulted in a better yield for 86% of hybrids compared to the average yield of their parents. In addition, hybrids appeared to be less affected by the negative correlation between protein content and yield compared to pure lines. These results shed light on the physiological bases underlying yield heterosis in wheat, paving new ways to breed for better wheat hybrids that can help to meet agriculture’s challenges.

Abstract

To meet the challenge of feeding almost 10 billion people by 2050, wheat yield has to double by 2050. However, over the past 20 years, yield increase has slowed down and even stagnated in the main producing countries. Following the example of maize, hybrids have been suggested as a solution to overcome yield stagnation in wheat. However, wheat heterosis is still limited and poorly understood. Gaining a better understanding of hybrid vigor holds the key to breed for better varieties. To this aim, we have developed and phenotyped for physiological and agronomic traits an incomplete factorial design consisting of 91 hybrids and their nineteen female and sixteen male parents. Monitoring the plant development with normalized difference vegetation index revealed that 89% of the hybrids including the five higher yielding hybrids had a longer grain filling phase with a delayed senescence that results in larger grain size. This average increase of 7.7% in thousand kernel weight translated to a positive mid-parent heterosis for grain yield for 86% of hybrids. In addition, hybrids displayed a positive grain protein deviation leading to a +4.7% heterosis in protein yield. These results shed light on the physiological bases underlying yield heterosis in wheat, paving new ways to breed for better wheat hybrids.

Dissection of the impact of prioritized QTL-linked and -unlinked SNP markers on the accuracy of genomic selection1

Background

Use of genomic information has resulted in an undeniable improvement in prediction accuracies and an increase in genetic gain in animal and plant genetic selection programs in spite of oversimplified assumptions about the true biological processes. Even for complex traits, a large portion of markers do not segregate with or effectively track genomic regions contributing to trait variation; yet it is not clear how genomic prediction accuracies are impacted by such potentially nonrelevant markers. In this study, a simulation was carried out to evaluate genomic predictions in the presence of markers unlinked with trait-relevant QTL. Further, we compared the ability of the population statistic F_ST and absolute estimated marker effect as preselection statistics to discriminate between linked and unlinked markers and the corresponding impact on accuracy.

Results

We found that the accuracy of genomic predictions decreased as the proportion of unlinked markers used to calculate the genomic relationships increased. Using all, only linked, and only unlinked marker sets yielded prediction accuracies of 0.62, 0.89, and 0.22, respectively. Furthermore, it was found that prediction accuracies are severely impacted by unlinked markers with large spurious associations. F_ST-preselected marker sets of 10 k and larger yielded accuracies 8.97 to 17.91% higher than those achieved using preselection by absolute estimated marker effects, despite selecting 5.1 to 37.7% more unlinked markers and explaining 2.4 to 5.0% less of the genetic variance. This was attributed to false positives selected by absolute estimated marker effects having a larger spurious association with the trait of interest and more negative impact on predictions. The Pearson correlation between F_ST scores and absolute estimated marker effects was 0.77 and 0.27 among only linked and only unlinked markers, respectively. The sensitivity of F_ST scores to detect truly linked markers is comparable to absolute estimated marker effects but the consistency between the two statistics regarding false positives is weak.

Conclusion

Identification and exclusion of markers that have little to no relevance to the trait of interest may significantly increase genomic prediction accuracies. The population statistic F_ST presents an efficient and effective tool for preselection of trait-relevant markers.

Evolution of Recombination Landscapes in Diverging Populations of Bread Wheat

Reciprocal exchanges of DNA (crossovers) that occur during meiosis are mandatory to ensure the production of fertile gametes in sexually reproducing species. They also contribute to shuffle parental alleles into new combinations thereby fueling genetic variation and evolution. However, due to biological constraints, the recombination landscape is highly heterogeneous along the genome which limits the range of allelic combinations and the adaptability of populations. An approach to better understand the constraints on the recombination process is to study how it evolved in the past. In this work, we tackled this question by constructing recombination profiles in four diverging bread wheat (Triticum aestivum L.) populations established from 371 landraces genotyped at 200,062 SNPs. We used linkage disequilibrium (LD) patterns to estimate in each population the past distribution of recombination along the genome and characterize its fine-scale heterogeneity. At the megabase scale, recombination rates derived from LD patterns were consistent with family-based estimates obtained from a population of 406 recombinant inbred lines. Among the four populations, recombination landscapes were positively correlated between each other and shared a statistically significant proportion of highly recombinant intervals. However, this comparison also highlighted that the similarity in recombination landscapes between populations was significantly decreasing with their genetic differentiation in most regions of the genome. This observation was found to be robust to SNPs ascertainment and demography and suggests a relatively rapid evolution of factors determining the fine-scale localization of recombination in bread wheat.