Shifting the limits in wheat research and breeding using a fully annotated reference genome

DArTseq genotyping facilitates the transfer of "exotic" chromatin from a Secale cereale × S. strictum hybrid into wheat

Cultivated and wild species of the genus rye (Secale) are important but underexploited gene sources for increasing the genetic diversity of bread wheat. Gene transfer is possible via bridge genetic materials derived from intergeneric hybrids. During this process, it is essential to precisely identify the rye chromatin in the wheat genetic background. In the present study, backcross generation BC2F8 from a cross between Triticum aestivum (Mv9kr1) and S. cereanum ('Kriszta,' a cultivar from the artificial hybrid of S. cereale and S. strictum) was screened using in-situ hybridization (GISH and FISH) and analyzed by DArTseq genotyping in order to select potentially agronomically useful genotypes for prebreeding purposes. Of the 329,267 high-quality short sequence reads generated, 27,822 SilicoDArT and 8,842 SNP markers specific to S. cereanum 1R-7R chromosomes were identified. Heatmaps of the marker densities along the 'Lo7' rye reference pseudomolecules revealed subtle differences between the FISH- and DArTseq-based results. This study demonstrates that the "exotic" rye chromatin of S. cereanum introgressed into wheat can be reliably identified by high-throughput DArTseq genotyping. The Mv9kr1-'Kriszta' addition and translocation lines presented here may serve as valuable prebreeding genetic materials for the development of stress-tolerant or disease-resistant wheat varieties.

Expression patterns of candidate genes for the Lr46/Yr29 "slow rust" locus in common wheat (Triticum aestivum L.) and associated miRNAs inform of the gene conferring the Puccinia triticina resistance trait

Leaf rust caused by Puccinia triticina (Pt) is one of the most impactful diseases causing substantial losses in common wheat (Triticum aestivum L.) crops. In adult plants resistant to Pt, a horizontal adult plant resistance (APR) is observed: APR protects the plant against multiple pathogen races and is distinguished by durable persistence under production conditions. The Lr46/Yr29 locus was mapped to chromosome 1B of common wheat genome, but the identity of the underlying gene has not been demonstrated although several candidate genes have been proposed. This study aimed to analyze the expression of nine candidate genes located at the Lr46/Yr29 locus and their four complementary miRNAs (tae-miR5384-3p, tae-miR9780, tae-miR9775, and tae-miR164), in response to Pt infection. The plant materials tested included five reference cultivars in which the molecular marker csLV46G22 associated with the Lr46/Yr29-based Pt resistance was identified, as well as one susceptible control cultivar. Biotic stress was induced in adult plants by inoculation with fungal spores under controlled conditions. Plant material was sampled before and at 6, 12, 24, 48 hours post inoculation (hpi). Differences in expression of candidate genes at the Lr46/Yr29 locus were analyzed by qRT-PCR and showed that the expression of the genes varied at the analyzed time points. The highest expression of Lr46/Yr29 candidate genes (Lr46-Glu1, Lr46-Glu2, Lr46-Glu3, Lr46-RLK1, Lr46-RLK2, Lr46-RLK3, Lr46-RLK4, Lr46-Snex, and Lr46-WRKY) occurred at 12 and 24 hpi and such expression profiles were obtained only for one candidate gene among the nine genes analyzed (Lr46-Glu2), indicating that it may be a contributing factor in the resistance response to Pt infection.

Identification and characterization of sulphotransferase (SOT) genes for tolerance against drought and heat in wheat and six related species

BACKGROUND

Sulphotransferase (SOT) enzyme (encoded by a conserved family of SOT genes) is involved in sulphonation of a variety of compounds, through transfer of a sulphuryl moiety from 3'phosphoadenosine- 5'phosphosulphate (PAPS) to a variety of secondary metabolites. The PAPS itself is derived from 3'adenosine-5'phosphosulphate (APS) that is formed after uptake of sulphate ions from the soil. The process provides tolerance against abiotic stresses like drought and heat in plants. Therefore, a knowledge of SOT genes in any crop may help in designing molecular breeding methods for improvement of tolerance for drought and heat.

METHODS

Sequences of rice SOT genes and SOT domain (PF00685) of corresponding proteins were both used for identification of SOT genes in wheat and six related species (T. urartu, Ae. tauschii, T. turgidum, Z. mays, B. distachyon and Hordeum vulgare), although detailed analysis was conducted only in wheat. The wheat genes were mapped on individual chromosomes and also subjected to synteny and collinearity analysis. The proteins encoded by these genes were examined for the presence of a complete SOT domain using 'Conserved Domain Database' (CDD) search tool at NCBI.

RESULTS

In wheat, 107 TaSOT genes, ranging in length from 969 bp to 7636 bp, were identified and mapped onto individual chromosomes. SSRs (simple sequence repeats), microRNAs, long non-coding RNAs (lncRNAs) and their target sites were also identified in wheat SOT genes. SOT proteins were also studied in detail. An expression assay of TaSOT genes via wheat RNA-seq data suggested engagement of these genes in growth, development and responses to various hormones and biotic/abiotic stresses.

CONCLUSIONS

The results of the present study should help in further functional characterization of SOT genes in wheat and other related crops.

Wheat Cybrid Plants, OryzaWheat, Regenerated from Wheat-Rice Hybrid Zygotes via in Vitro Fertilization System Possess Wheat-Rice Hybrid Mitochondria

Hybridization generates biodiversity, and wide hybridization plays a pivotal role in enhancing and broadening the useful attributes of crops. The hybridization barrier between wheat and rice, the two most important cereals, was recently overcome by in vitro production of allopolyploid wheat-rice hybrid zygotes, which can develop and grow into mature plants. In the study, genomic sequences and compositions of the possible hybrid plants were investigated through short- and long-read sequencing analyses and fluorescence in situ hybridization (FISH)-based visualization. The possible hybrid possessed whole wheat nuclear and cytoplasmic DNAs and rice mitochondrial (mt) DNA, along with variable retention rates of rice mtDNA ranging from 11% to 47%. The rice mtDNA retained in the wheat cybrid, termed Oryzawheat, can be transmitted across generations. In addition to mitochondrial hybridization, translocation of rice chromosome 1 into wheat chromosome 6A was detected in a F1 hybrid individual. OryzaWheat can provide a new horizon for utilizing inter-subfamily genetic resources among wheat and rice belonging to different subfamilies, Pooideae and Ehrhartoideae, respectively.

Dissecting the genetic basis of resistance to Soil-borne cereal mosaic virus (SBCMV) in durum wheat by bi-parental mapping and GWAS

Soil-borne cereal mosaic virus (SBCMV), the causative agent of wheat mosaic, is a Furovirus challenging wheat production all over Europe. Differently from bread wheat, durum wheat shows greater susceptibility and stronger yield penalties, so identification and genetic characterization of resistance sources are major targets for durum genetics and breeding. The Sbm1 locus providing high level of resistance to SBCMV was mapped in bread wheat to the 5DL chromosome arm (Bass in Genome 49:1140-1148, 2006). This excluded the direct use of Sbm1 for durum wheat improvement. Only one major QTL has been mapped in durum wheat, namely QSbm.ubo-2B, on the 2BS chromosome region coincident with Sbm2, already known in bread wheat as reported (Bayles in HGCA Project Report, 2007). Therefore, QSbm.ubo-2B = Sbm2 is considered a pillar for growing durum in SBCMV-affected areas. Herein, we report the fine mapping of Sbm2 based on bi-parental mapping and GWAS, using the Infinium 90 K SNP array and high-throughput KASP®. Fine mapping pointed out a critical haploblock of 3.2 Mb defined by concatenated SNPs successfully converted to high-throughput KASP® markers coded as KUBO. The combination of KUBO-27, wPt-2106-ASO/HRM, KUBO-29, and KUBO-1 allows unequivocal tracing of the Sbm2-resistant haplotype. The interval harbors 52 high- and 41 low-confidence genes, encoding 17 cytochrome p450, three receptor kinases, two defensins, and three NBS-LRR genes. These results pave the way for Sbm2 positional cloning. Importantly, the development of Sbm2 haplotype tagging KASP® provides a valuable case study for improving efficacy of the European variety testing system and, ultimately, the decision-making process related to varietal characterization and choice.

Genotype selection identified elite lines through quantitative trait loci mapping of agronomically important traits in wheat

Wheat is one of the most important staple foods in the world. Genetic characterization of wheat agronomically important traits is crucial for yield improvement through molecular breeding. In this study, a recombinant inbred line (RIL) population was developed by crossing a local adapted high yield variety Jimai 22 (JM22) with an external variety Cunmai no.1 (CM1). A high-density genetic map containing 7,359 single nucleotide polymorphism (SNP) markers was constructed. Quantitative trait loci (QTL) mapping identified 61 QTL for eight yield-related traits under six environments (years). Among them, 17 QTL affecting spike number per plant, grain number per spike and thousand grain weight showed high predictability for theoretical yield per plant (TYP), of which, 12 QTL alleles positively contributed to TYP. Nine promising candidate genes for seven of the 12 QTL were identified including three known wheat genes and six rice orthologs. Four elite lines with TYP increased by 5.6%-15.2% were identified through genotype selection which carried 7-9 favorable alleles from JM22 and 2-3 favorable alleles from CM1 of the 12 QTL. Moreover, the linked SNPs of the 12 QTL were converted to high-throughput kompetitive allele-specific PCR (KASP) markers and validated in the population. The mapped QTL, identified promising candidate genes, developed elite lines and KASP markers are highly valuable in future genotype selection to improve wheat yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01496-3.

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.

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

Rmg10, a Novel Wheat Blast Resistance Gene Derived from Aegilops tauschii

Wheat blast, caused by Pyricularia oryzae (syn. Magnaporthe oryzae) pathotype Triticum (MoT), is a devastating disease that can result in up to 100% yield loss in affected fields. To find new resistance genes against wheat blast, we screened 199 accessions of Aegilops tauschii, the D genome progenitor of common wheat (Triticum aestivum), by seedling inoculation assays with Brazilian MoT isolate Br48 and found 14 resistant accessions. A synthetic hexaploid wheat line (Ldn/KU-2097) derived from a cross between the T. turgidum 'Langdon' (Ldn) and resistant A. tauschii accession KU-2097 exhibited resistance in seedlings and spikes against Br48. In an F2 population derived from 'Chinese Spring' × Ldn/KU-2097, resistant and susceptible individuals segregated in a 3:1 ratio, suggesting that the resistance from KU-2097 is controlled by a single dominant gene. We designated this gene Rmg10. Genetic mapping using an F2:3 population from the same cross mapped the RMG10 locus to the short arm of chromosome 2D. Rmg10 was ineffective against Bangladesh isolates but effective against Brazilian isolates. Field tests in Bolivia showed increased spike resistance in a synthetic octaploid wheat line produced from a cross between common wheat cultivar 'Gladius' and KU-2097. These results suggest that Rmg10 would be beneficial in farmers' fields in South America.

Elucidation of the genetic architecture of water absorption capacity in hard winter wheat through genome wide association study

Water absorption capacity (WAC) influences various aspects of bread making, such as loaf volume, bread yield, and shelf life. Despite its importance in the baking process and end-product quality, its genetic determinants are less explored. To address this limitation, a genome-wide association study was conducted on 337 hard wheat (Triticum aestivum L.) genotypes evaluated over 5 years in multi-environmental trials. Phenotyping was done using the solvent retention capacity (SRC) test with water (SRC-water), sucrose (SRC-sucrose), lactic acid (SRC-lactic acid), and sodium carbonate (SRC-carbonate) as solvents. Individuals were genotyped using genotyping-by-sequencing to detect single nucleotide polymorphisms across the wheat genome. To detect the genomic regions that underline the SRCs and gluten performance index (GPI), a genome-wide association study was performed using six multi-locus models using the mrMLM package in R. Adjusted means for SRC-water ranged from 54.1% to 66.5%, while SRC-carbonate exhibited a narrow range from 84.9% to 93.9%. Moderate to high genomic heritability values were observed for SRCs and GPI, ranging from h2 = 0.61 to 0.88. The genome-wide association study identified a total of 42 quantitative trait nucleotides (QTNs), of which five explained over 10% of the phenotypic variation (R2 ≥ 10%). Most of the QTNs were detected on chromosomes 1A, 1B, 3B, and 5B. Few QTNs, such as S1A_5190318, S1B_3282665, S4D_472908721, and S7A_37433960, were located near gliadin, glutenin starch synthesis, and galactosyltransferase genes. Overall, these results show WAC to be under polygenic genetic control, with genes involved in the synthesis of key flour components influencing overall water absorption.

Ploidy variation induces butterfly effect on chromatin topology in wheat

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

Karyotypic and phenotypic condensation in allotetraploid wheats accompanied with reproductive strategy transformation: from natural evolution to domestication

MAIN CONCLUSION

Allotetraploid wheat reflects evolutionary divergence and domestication convergence in the karyotypic and phenotypic evolution, accompanied with the transformation from r- strategy to K- strategy in reproductive fitness. Allotetraploid wheat, the progenitor of hexaploidy bread wheat, has undergone 300,000 years of natural evolution and 10,000 years of domestication. The variations in karyotype and phenotype as well as fertility fitness have not been systematically linked. Here, by combining fluorescent in situ hybridization with the quantification of phenotypic and reproductive traits, we compared the karyotype, vegetative growth phenotype and reproductive fitness among synthesized, wild and domesticated accessions of allotetraploid wheat. We detected that the wild accessions showed dramatically high frequencies of homologous recombination and copy number variations of simple sequence repeats (SSR) comparing with synthetic and domesticated accessions. The phenotypic traits reflected significant differences among the populations shaped by distinct evolutionary processes. The diversity observed in wild accessions was significantly greater than that in domesticated ones, particularly in traits associated with vegetative growth and spike morphology. We found that the active pollen of domesticated accessions exhibited greater potential of germination, despite a lower rate of active pollen compared with the wild accessions, indicating a transformation in reproductive fitness strategy for pollen development in domesticated accessions compared to the wild accessions, from r-strategy to K-strategy. Our results demonstrate the condensation of karyotype and phenotype from natural wild accessions to domesticated accessions in allotetraploid wheats. Ecological strategy transformation should be seriously considered from evolution to domestication in polyploid plants, especially crops, which may provide a perspective on the adaptive evolution of polyploid plants.

Identification of a recessive gene RgM4G52 conferring red glume, stem, and rachis in a Triticum boeoticum mutant

Anthocyanins are plant secondary metabolites belonging to the polyphenol class of natural water-soluble phytopigments. The accumulation of anthocyanins in different plant tissues can improve plant survival under adverse conditions. In addition, plants with the resulting colorful morphology can be utilized as landscape plants. Triticum boeoticum (syn. Triticum monococcum ssp. aegilopoides, 2n=2x=14, AbAb) serves as a valuable genetic resource for the improvement of its close relative common wheat in terms of enhancing resilience to various biotic and abiotic stresses. In our previous study, the EMS-mutagenized mutant Z2921 with a red glume, stem, and rachis was generated from T. boeoticum G52, which has a green glume, stem, and rachis. In this study, the F1, F2, and F2:3 generations of a cross between mutant-type Z2921 and wild-type G52 were developed. A single recessive gene, tentatively designated RgM4G52, was identified in Z2921 via genetic analysis. Using bulked segregant exome capture sequencing (BSE-Seq) analysis, RgM4G52 was mapped to chromosome 6AL and was flanked by the markers KASP-58 and KASP-26 within a 3.40-cM genetic interval corresponding to 1.71-Mb and 1.61-Mb physical regions in the Chinese Spring (IWGSC RefSeq v1.1) and Triticum boeoticum (TA299) reference genomes, respectively, in which seven and four genes related to anthocyanin synthesis development were annotated. Unlike previously reported color morphology-related genes, RgM4G52 is a recessive gene that can simultaneously control the color of glumes, stems, and rachis in wild einkorn. In addition, a synthetic Triticum dicoccum-T. boeoticum amphiploid Syn-ABAb-34, derived from the colchicine treatment of F1 hybrids between tetraploid wheat PI 352367 (T. dicoccum, AABB) and Z2921, expressed the red stems of Z2921. The flanking markers of RgM4G52 developed in this study could be useful for developing additional common wheat lines with red stems, laying the foundation for marker-assisted breeding and the fine mapping of RgM4G52.

Preferentially expressed endosperm genes reveal unique activities in wheat endosperm during grain filling

BACKGROUND

Bread wheat (Triticum aestivum L.) endosperm contains starch and proteins, which determine the final yield, quality, and nutritional value of wheat grain. The preferentially expressed endosperm genes can precisely provide targets in the endosperm for improving wheat grain quality and nutrition using modern bioengineering technologies. However, the genes specifically expressed in developing endosperms remain largely unknown.

RESULTS

In this study, 315 preferentially expressed endosperm genes (PEEGs) in the spring wheat landrace, Chinese Spring, were screened using data obtained from an open bioinformatics database, which reveals a unique grain reserve deposition process and special signal transduction in a developing wheat endosperm. Furthermore, transcription and accumulation of storage proteins in the wheat cultivar, XC26 were evaluated. The results revealed that 315 PEEG plays a critical role in storage protein fragment deposition and is a potential candidate for modifying grain quality and nutrition.

CONCLUSION

These results provide new insights into endosperm development and candidate genes and promoters for improving wheat grain quality through genetic engineering and plant breeding techniques.

Differential brainstem connectivity according to sex and menopausal status in healthy men and women

Background Brainstem nuclei play a critical role in both ascending monoaminergic modulation of cortical function and arousal, and in descending bulbospinal pain modulation. Even though sex-related differences in the function of both systems have been reported in animal models, a complete understanding of sex differences, as well as menopausal effects, in brainstem connectivity in humans is lacking. This study evaluated resting-state connectivity of the dorsal raphe nucleus (DRN), right and left locus coeruleus complex (LCC), and periaqueductal gray (PAG) according to sex and menopausal status in healthy individuals. In addition, relationships between systemic estrogen levels and brainstem-network connectivity were examined in a subset of participants. Methods Resting-state fMRI was performed in 50 healthy men (age, 31.2 ± 8.0 years), 53 healthy premenopausal women (age, 24.7 ± 7.3 years; 22 in the follicular phase, 31 in the luteal phase), and 20 postmenopausal women (age, 54.6 ± 7.2 years). Permutation Analysis of Linear Models (5000 permutations) was used to evaluate differences in brainstem-network connectivity according to sex and menopausal status, controlling for age. In 10 men and 17 women (9 premenopausal; 8 postmenopausal), estrogen and estrogen metabolite levels in plasma and stool were determined by liquid chromatography-mass spectrometry/mass spectrometry. Relationships between estrogen levels and brainstem-network connectivity were evaluated by partial least squares analysis. Results Left LCC-executive control network (ECN) connectivity showed an overall sex difference (p = 0.02), with higher connectivity in women than in men; however, this was mainly due to differences between men and pre-menopausal women (p = 0.008). Additional sex differences were dependent on menopausal status: PAG-default mode network (DMN) connectivity was higher in postmenopausal women than in men (p = 0.04), and PAG-sensorimotor network (SMN) connectivity was higher in premenopausal women than in men (p = 0.03) and postmenopausal women (p = 0.007). Notably, higher free 2-hydroxyestrone levels in stool were associated with higher PAG-SMN and PAG-DMN connectivity in premenopausal women (p < 0.01). Conclusions Healthy women show higher brainstem-network connectivity involved in cognitive control, sensorimotor function, and self-relevant processes than men, dependent on their menopausal status. Further, 2-hydroxyestrone, implicated in pain, may modulate PAG connectivity in premenopausal women. These findings may relate to differential vulnerabilities to chronic stress-sensitive disorders at different life stages.

β-diketone accumulation in response to drought stress is weakened in modern bread wheat varieties (Triticum aestivum L.)

Cuticular waxes coating leaf surfaces can help plants tolerate drought events by reducing non-stomatal water loss. Despite their role in drought tolerance, little is known about how cuticular wax composition has changed during breeding in Canadian bread wheat (Triticum aestivum L.) varieties. To fill in this gap, flag leaves of the Canadian Heritage Bread Wheat Panel, which include 30 varieties released between 1842 and 2018, were surveyed to determine if and how cuticular wax composition in wheat has changed at two breeding ecozones over this period. Following this, a subset of varieties was subjected to drought conditions to compare their responses. As expected, modern varieties outperformed old varieties with a significantly larger head length and reaching maturity earlier. Yet, when challenged with drought, old varieties were able to significantly increase the accumulation of β-diketones to a higher extent than modern varieties. Furthermore, RNAseq was performed on the flag leaf of four modern varieties to identify potential markers that could be used for selection of higher accumulation of cuticular waxes. This analysis revealed that the W1 locus is a good candidate for selecting higher accumulation of β-diketones. These findings indicate that the variation in cuticular waxes upon drought could be further incorporated in breeding of future bread wheat varieties.

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.

Isocitrate lyase promotes Puccinia striiformis f. sp. tritici susceptibility in wheat (Triticum aestivum) by suppressing accumulation of glyoxylate cycle intermediates

Plant fungal parasites manipulate host metabolism to support their own survival. Among the many central metabolic pathways altered during infection, the glyoxylate cycle is frequently upregulated in both fungi and their host plants. Here, we examined the response of the glyoxylate cycle in bread wheat (Triticum aestivum) to infection by the obligate biotrophic fungal pathogen Puccinia striiformis f. sp. tritici (Pst). Gene expression analysis revealed that wheat genes encoding the two unique enzymes of the glyoxylate cycle, isocitrate lyase (TaICL) and malate synthase, diverged in their expression between susceptible and resistant Pst interactions. Focusing on TaICL, we determined that the TaICL B homoeolog is specifically upregulated during early stages of a successful Pst infection. Furthermore, disruption of the B homoeolog alone was sufficient to significantly perturb Pst disease progression. Indeed, Pst infection of the TaICL-B disruption mutant (TaICL-BY400*) was inhibited early during initial penetration, with the TaICL-BY400* line also accumulating high levels of malic acid, citric acid, and aconitic acid. Exogenous application of malic acid or aconitic acid also suppressed Pst infection, with trans-aconitic acid treatment having the most pronounced effect by decreasing fungal biomass 15-fold. Thus, enhanced TaICL-B expression during Pst infection may lower accumulation of malic acid and aconitic acid to promote Pst proliferation. As exogenous application of aconitic acid and malic acid has previously been shown to inhibit other critical pests and pathogens, we propose TaICL as a potential target for disruption in resistance breeding that could have wide-reaching protective benefits for wheat and beyond.

Genome-wide analysis of salicylic acid and jasmonic acid signalling marker gene families in wheat

Jasmonic acid (JA) and salicylic acid (SA) phytohormone pathways are important regulators of stress tolerance. Knowledge regarding the diversity, phylogeny and functionality of wheat genes involved in JA and SA response is limited. Using Arabidopsis, rice and wheat genomic and wheat disease transcriptomic data, we deduced the size, phylogenetic diversity and pathogen-responsiveness of seven hormone-responsive gene families, and thus selected 14 candidates as potential hormone responsive gene markers. Gene-specific expression studies assessed the impact of exogenous JA and SA on their transcriptional activation in leaves of two distinct wheat cultivars. RNAseq data were interrogated to assess their disease responsiveness and tissue-specific expression. This study elucidated the number, phylogeny and pathogen-responsiveness of wheat genes from seven families, including 12 TaAOS, 6 TaJAMyb, 256 TaWRKY group III, 85 TaPR1, 205 TaPR2, 76 TaPR3 and 124 TaPR5. This included the first description of the wheat AOS, JAMyb, PR2, PR3 and PR5 gene families. Gene expression studies delineated TaAOS1-5B and TaJAMyb-4A as JA-responsive in leaves, but not significantly responsive to SA treatment, while TaWRKY45-B was a SA- but not a JA-responsive marker. Other candidate genes were either unresponsive or non-specific to SA or JA. Our findings highlight that all seven gene families are greatly expanded in wheat as compared to other plants (up to 7.6-fold expansion), and demonstrate disparity in the response to biotic stress between some homoeologous and paralogous sequences within these families. The SA- and JA-responsive marker genes identified herein will prove useful tools to monitor these signalling pathways in wheat.

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.

Association mapping of tan spot and septoria nodorum blotch resistance in cultivated emmer wheat

KEY MESSAGE

A total of 65 SNPs associated with resistance to tan spot and septoria nodorum blotch were identified in a panel of 180 cultivated emmer accessions through association mapping Tan spot and septoria nodorum blotch (SNB) are foliar diseases caused by the respective fungal pathogens Pyrenophora tritici-repentis and Parastagonospora nodorum that affect global wheat production. To find new sources of resistance, we evaluated a panel of 180 cultivated emmer wheat (Triticum turgidum ssp. dicoccum) accessions for reactions to four P. tritici-repentis isolates Pti2, 86-124, 331-9 and DW5, two P. nodorum isolate, Sn4 and Sn2000, and four necrotrophic effectors (NEs) produced by the pathogens. About 8-36% of the accessions exhibited resistance to the four P. tritici-repentis isolates, with five accessions demonstrating resistance to all isolates. For SNB, 64% accessions showed resistance to Sn4, 43% to Sn2000 and 36% to both isolates, with Spain (11% accessions) as the most common origin of resistance. To understand the genetic basis of resistance, association mapping was performed using SNP (single nucleotide polymorphism) markers generated by genotype-by-sequencing and the 9 K SNP Infinium array. A total of 46 SNPs were significantly associated with tan spot and 19 SNPs with SNB resistance or susceptibility. Six trait loci on chromosome arms 1BL, 3BL, 4AL (2), 6BL and 7AL conferred resistance to two or more isolates. Known NE sensitivity genes for disease development were undetected except Snn5 for Sn2000, suggesting novel genetic factors are controlling host-pathogen interaction in cultivated emmer. The emmer accessions with the highest levels of resistance to the six pathogen isolates (e.g., CItr 14133-1, PI 94634-1 and PI 377672) could serve as donors for tan spot and SNB resistance in wheat breeding programs.

Three novel QTLs for FHB resistance identified and mapped in spring wheat PI672538 by bulked segregant analysis of the recombinant inbred line

Introduction

Fusarium head blight (FHB) has a large influence on both the yield and quality of wheat grain worldwide. Host resistance is the most effective method for controlling FHB, but unfortunately, very few genetic resources on FHB resistance are available; therefore, identifying novel resistance genes or quantitative trait loci (QTLs) is valuable.

Methods

Here, a recombinant inbred line (RIL) population containing 451 lines derived from the cross L661/PI672538 was sown in four different environments (2019CZa, 2019CZb, 2021QL and 2021WJ).

Results

Five QTLs, consisting of two previously reported QTLs (FhbL693a and FhbL693b) and three new QTLs (FhbL693c, FhbL693d and FhbL693e), were identified. Further investigation revealed that FhbL693b, FhbL693c and FhbL693d could be detected in all four environments, and FhbL693a and FhbL693e were detected only in 2019CZb and 2021WJ, respectively. Among the QTLs, the greatest contribution (10.5%) to the phenotypic variation effect (PVE) was FhbL693d in 2021WJ, while the smallest (1.2%) was FhbL693e and FhbL693a in 2019CZb. The selection of 5Dindel-4 for FhbL693d, 4Aindel-7 for FhbL693c and 3Bindel-24 for FhbL693b decreased the number of damaged spikelets by 2.1, and a new line resistant to FHB named H140-2 was developed by marker-assisted selection (MAS).

Discussion

These results could help to further improve FHB resistance in the future.

Chromosome-scale assembly of the wild cereal relative Elymus sibiricus

Elymus species, belonging to Triticeae tribe, is a tertiary gene pool for improvement of major cereal crops. Elymus sibiricus, a tetraploid with StH genome, is a typical species in the genus Elymus, which is widely utilized as a high-quality perennial forage grass in template regions. In this study, we report the construction of a chromosome-scale reference assembly of E. sibiricus line Gaomu No. 1 based on PacBio HiFi reads and chromosome conformation capture. Subgenome St and H were well phased by assisting with kmer and subgenome-specific repetitive sequence. The total assembly size was 6.929 Gb with a contig N50 of 49.518 Mb. In total, 89,800 protein-coding genes were predicted. The repetitive sequences accounted for 82.49% of the genome in E. sibiricus. Comparative genome analysis confirmed a major species-specific 4H/6H reciprocal translocation in E. sibiricus. The E. sibiricus assembly will be much helpful to exploit genetic resource of StH species in genus Elymus, and provides an important tool for E. sibiricus domestication.

The addition of Psathyrostachys Huashanica Keng 6Ns large segment chromosomes has positive impact on stripe rust resistance and plant spikelet number of common wheat

BACKGROUND

Developing novel germplasm by using wheat wild related species is an effective way to rebuild the wheat resource bank. The Psathyrostachys huashanica Keng (P. huashanica, 2n = 2x = 14, NsNs) is regarded as a superior species to improve wheat breeding because of its multi-resistance, early maturation and numerous tiller traits. Introducing genetic components of P. huashanica into the common wheat background is the most important step in achieving the effective use. Therefore, the cytogenetic characterization and influence of the introgressed P. huashanica large segment chromosomes in the wheat background is necessary to be explored.

RESULTS

In this study, we characterized a novel derived line, named D88-2a, a progeny of the former characterized wheat-P. huashanica partial amphiploid line H8911 (2n = 7x = 49, AABBDDNs). Cytological identification showed that the chromosomal composition of D88-2a was 2n = 44 = 22II, indicating the addition of exogenous chromosomes. Genomic in situ hybridization demonstrated that the supernumerary chromosomes were a pair of homologues from the P. huashanica and could be stably inherited in the common wheat background. Molecular markers and 15 K SNP array indicated that the additional chromosomes were derived from the sixth homoeologous group (i.e., 6Ns) of P. huashanica. Based on the distribution of the heterozygous single-nucleotide polymorphism sites and fluorescence in situ hybridization karyotype of each chromosome, this pair of additional chromosomes was confirmed as P. huashanica 6Ns large segment chromosomes, which contained the entire short arm and the proximal centromere portion of the long arm. In terms of the agronomic traits, the addition line D88-2a exhibited enhanced stripe rust resistance, improved spike characteristics and increased protein content than its wheat parent line 7182.

CONCLUSIONS

The new wheat germplasm D88-2a is a novel cytogenetically stable wheat-P. huashanica 6Ns large segment addition line, and the introgressed large segment alien chromosome has positive impact on plant spikelet number and stripe rust resistance. Thus, this germplasm can be used for genetic improvement of cultivated wheat and the study of functional alien chromosome segment.

Unveiling the Mechanism of Phenamacril Resistance in F. graminearum: Computational and Experimental Insights into the C423A Mutation in FgMyoI

Phenamacril (PHA) is a highly selective fungicide for controlling fusarium head blight (FHB) mainly caused by F. graminearum and F. asiaticum. However, the C423A mutation in myosin I of F. graminearum (FgMyoI) leads to natural resistance to PHA. Here, based on the computational approaches and biochemical validation, we elucidate the atomic-level mechanism behind the natural resistance of F. graminearum to the fungicide PHA due to the C423A mutation in FgMyoI. The mutation leads to a rearrangement of pocket residues, resulting in increased size and flexibility of the binding pocket, which impairs the stable binding of PHA. MST experiments confirm that the mutant protein FgMyoIC423A exhibits significantly reduced affinity for PHA compared to wild-type FgMyoI and the nonresistant C423K mutant. This decreased binding affinity likely underlies the development of PHA resistance in F. graminearum. Conversely, the nonresistant C423K mutant retains sensitivity to PHA due to the introduction of a strong hydrogen bond donor, which facilitates stable binding of PHA in the pocket. These findings shed light on the molecular basis of PHA resistance and provide new directions for the creation of new myosin inhibitors.

The 4T and 7T introgressions from Amblyopyrum muticum and the 5Au introgression from Triticum urartu increases grain zinc and iron concentrations in Malawian wheat backgrounds

Micronutrient deficiencies (MNDs) particularly zinc (Zn) and iron (Fe) remain widespread in sub-Saharan Africa (SSA) due to low dietary intake. Wheat is an important source of energy globally, although cultivated wheat is inherently low in grain micronutrient concentrations. Malawian wheat/Am. muticum and Malawian wheat/T. urartu BC1F3 introgression lines, developed by crossing three Malawian wheat varieties (Kenya nyati, Nduna and Kadzibonga) with DH-348 (wheat/Am. muticum) and DH-254 (wheat/T. urartu), were phenotyped for grain Zn and Fe, and associated agronomic traits in Zn-deficient soils, in Malawi. 98% (47) of the BC1F3 introgression lines showed higher Zn above the checks Paragon, Chinese Spring, Kadzibonga, Kenya Nyati and Nduna. 23% (11) of the introgression lines showed a combination of high yields and an increase in grain Zn by 16-30 mg kg -1 above Nduna and Kadzibonga, and 11-25 mg kg -1 above Kenya nyati, Paragon and Chinese Spring. Among the 23%, 64% (7) also showed 8-12 mg kg -1 improvement in grain Fe compared to Nduna and Kenya nyati. Grain Zn concentrations showed a significant positive correlation with grain Fe, whilst grain Zn and Fe negatively and significantly correlated with TKW and grain yield. This work will contribute to the efforts of increasing mineral nutrient density in wheat, specifically targeting countries in the SSA.

TaMYB44-5A reduces drought tolerance by repressing transcription of TaRD22-3A in the abscisic acid signaling pathway

MAIN CONCLUSION

TaMYB44-5A identified as a transcription factor negatively regulates drought tolerance in transgenic Arabidopsis. Drought can severely reduce yields throughout the wheat-growing season. Many studies have shown that R2R3-MYB transcription factors are involved in drought stress responses. In this study, the R2R3-MYB transcription factor MYB44-5A was identified in wheat (Triticum aestivum L.) and functionally analyzed. Three homologs of TaMYB44 were isolated, all of which localized to the nucleus. Overexpression of TaMYB44-5A reduced drought tolerance in Arabidopsis thaliana. Further analysis showed that TaMYB44-5A reduced the sensitivity of transgenic Arabidopsis to ABA. Genetic and transcriptional regulation analyses demonstrated that the expression levels of drought- and ABA-responsive genes were downregulated by TaMYB44-5A, and TaMYB44-5A directly bound to the MYB-binding site on the promoter to repress the transcription level of TaRD22-3A. Our results provide insights into a novel molecular pathway in which the R2R3-MYB transcription factor negatively regulates ABA signaling in response to drought stress.

Genome-wide atlas of rust resistance loci in wheat

Rust diseases, including leaf rust, stripe/yellow rust, and stem rust, significantly impact wheat (Triticum aestivum L.) yields, causing substantial economic losses every year. Breeding and deployment of cultivars with genetic resistance is the most effective and sustainable approach to control these diseases. The genetic toolkit for wheat breeders to select for rust resistance has rapidly expanded with a multitude of genetic loci identified using the latest advances in genomics, mapping and cloning strategies. The goal of this review was to establish a wheat genome atlas that provides a comprehensive summary of reported loci associated with rust resistance. Our atlas provides a summary of mapped quantitative trait loci (QTL) and characterised genes for the three rusts from 170 publications over the past two decades. A total of 920 QTL or resistance genes were positioned across the 21 chromosomes of wheat based on the latest wheat reference genome (IWGSC RefSeq v2.1). Interestingly, 26 genomic regions contained multiple rust loci suggesting they could have pleiotropic effects on two or more rust diseases. We discuss a range of strategies to exploit this wealth of genetic information to efficiently utilise sources of resistance, including genomic information to stack desirable and multiple QTL to develop wheat cultivars with enhanced resistance to rust disease.

High-Copy Transposons from a Pathogen Give Rise to a Conserved sRNA Family with a Novel Host Immunity Target

Small RNAs (sRNAs) are involved in gene silencing in multiple ways, including through cross-kingdom transfers from parasites to their hosts. Little is known about the evolutionary mechanisms enabling eukaryotic microbes to evolve functional mimics of host small regulatory RNAs. Here, we describe the identification and functional characterization of SINE_sRNA1, an sRNA family derived from highly abundant short interspersed nuclear element (SINE) retrotransposons in the genome of the wheat powdery mildew pathogen. SINE_sRNA1 is encoded by a sequence motif that is conserved in multiple SINE families and corresponds to a functional plant microRNA (miRNA) mimic targeting Tae_AP1, a wheat gene encoding an aspartic protease only found in monocots. Tae_AP1 has a novel function enhancing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), thereby contributing to the cross activation of plant defenses. We conclude that SINE_sRNA1 and Tae_AP1 are functional innovations, suggesting the contribution of transposons to the evolutionary arms race between a parasite and its host. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Beyond pathways: Accelerated flavonoids candidate identification and novel exploration of enzymatic properties using combined mapping populations of wheat

Although forward-genetics-metabolomics methods such as mGWAS and mQTL have proven effective in providing myriad loci affecting metabolite contents, they are somehow constrained by their respective constitutional flaws such as the hidden population structure for GWAS and insufficient recombinant rate for QTL. Here, the combination of mGWAS and mQTL was performed, conveying an improved statistical power to investigate the flavonoid pathways in common wheat. A total of 941 and 289 loci were, respectively, generated from mGWAS and mQTL, within which 13 of them were co-mapped using both approaches. Subsequently, the mGWAS or mQTL outputs alone and their combination were, respectively, utilized to delineate the metabolic routes. Using this approach, we identified two MYB transcription factor encoding genes and five structural genes, and the flavonoid pathway in wheat was accordingly updated. Moreover, we have discovered some rare-activity-exhibiting flavonoid glycosyl- and methyl-transferases, which may possess unique biological significance, and harnessing these novel catalytic capabilities provides potentially new breeding directions. Collectively, we propose our survey illustrates that the forward-genetics-metabolomics approaches including multiple populations with high density markers could be more frequently applied for delineating metabolic pathways in common wheat, which will ultimately contribute to metabolomics-assisted wheat crop improvement.

Wheat NAC transcription factor NAC5-1 is a positive regulator of senescence

Wheat (Triticum aestivum L.) is an important source of both calories and protein in global diets, but there is a trade-off between grain yield and protein content. The timing of leaf senescence could mediate this trade-off as it is associated with both declines in photosynthesis and nitrogen remobilization from leaves to grain. NAC transcription factors play key roles in regulating senescence timing. In rice, OsNAC5 expression is correlated with increased protein content and upregulated in senescing leaves, but the role of the wheat ortholog in senescence had not been characterized. We verified that NAC5-1 is the ortholog of OsNAC5 and that it is expressed in senescing flag leaves in wheat. To characterize NAC5-1, we combined missense mutations in NAC5-A1 and NAC5-B1 from a TILLING mutant population and overexpressed NAC5-A1 in wheat. Mutation in NAC5-1 was associated with delayed onset of flag leaf senescence, while overexpression of NAC5-A1 was associated with slightly earlier onset of leaf senescence. DAP-seq was performed to locate transcription factor binding sites of NAC5-1. Analysis of DAP-seq and comparison with other studies identified putative downstream target genes of NAC5-1 which could be associated with senescence. This work showed that NAC5-1 is a positive transcriptional regulator of leaf senescence in wheat. Further research is needed to test the effect of NAC5-1 on yield and protein content in field trials, to assess the potential to exploit this senescence regulator to develop high-yielding wheat while maintaining grain protein content.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

A k-mer-based pangenome approach for cataloging seed-storage-protein genes in wheat to facilitate genotype-to-phenotype prediction and improvement of end-use quality

Wheat is a staple food for more than 35% of the world's population, with wheat flour used to make hundreds of baked goods. Superior end-use quality is a major breeding target; however, improving it is especially time-consuming and expensive. Furthermore, genes encoding seed-storage proteins (SSPs) form multi-gene families and are repetitive, with gaps commonplace in several genome assemblies. To overcome these barriers and efficiently identify superior wheat SSP alleles, we developed "PanSK" (Pan-SSP k-mer) for genotype-to-phenotype prediction based on an SSP-based pangenome resource. PanSK uses 29-mer sequences that represent each SSP gene at the pangenomic level to reveal untapped diversity across landraces and modern cultivars. Genome-wide association studies with k-mers identified 23 SSP genes associated with end-use quality that represent novel targets for improvement. We evaluated the effect of rye secalin genes on end-use quality and found that removal of ω-secalins from 1BL/1RS wheat translocation lines is associated with enhanced end-use quality. Finally, using machine-learning-based prediction inspired by PanSK, we predicted the quality phenotypes with high accuracy from genotypes alone. This study provides an effective approach for genome design based on SSP genes, enabling the breeding of wheat varieties with superior processing capabilities and improved end-use quality.

A novel regulator of wheat tillering LT1 identified by using an upgraded BSA method, uni-BSA

Branching/tillering is a critical process for plant architecture and grain yield. However, Branching is intricately controlled by both endogenous and environmental factors. The underlying mechanisms of tillering in wheat remain poorly understood. In this study, we identified Less Tiller 1 (LT1) as a novel regulator of wheat tillering using an enhanced bulked segregant analysis (BSA) method, uni-BSA. This method effectively reduces alignment noise caused by the high repetitive sequence content in the wheat genome. Loss-of-function of LT1 results in fewer tillers due to defects in axillary meristem initiation and bud outgrowth. We mapped LT1 to a 6 Mb region on the chromosome 2D short arm and validated a nucleotide-binding (NB) domain encoding gene as LT1 using CRISPR/Cas9. Furthermore, the lower sucrose concentration in the shoot bases of lt1 might result in inadequate bud outgrowth due to disturbances in the sucrose biosynthesis pathways. Co-expression analysis suggests that LT1 controls tillering by regulating TaROX/TaLAX1, the ortholog of the Arabidopsis tiller regulator REGULATOR OF AXILLARY MERISTEM FORMATION (ROX) or the rice axillary meristem regulator LAX PANICLE1 (LAX1). This study not only offers a novel genetic resource for cultivating optimal plant architecture but also underscores the importance of our innovative BSA method. This uni-BSA method enables the swift and precise identification of pivotal genes associated with significant agronomic traits, thereby hastening gene cloning and crop breeding processes in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01484-7.

Chromosome-level assembly of the synthetic hexaploid wheat-derived cultivar Chuanmai 104

Synthetic hexaploid wheats (SHWs) are effective genetic resources for transferring agronomically important genes from wild relatives to common wheat (Triticum aestivum L.). Dozens of reference-quality pseudomolecule assemblies of hexaploid wheat have been generated, but none is reported for SHW-derived cultivars. Here, we generated a chromosome-scale assembly for the SHW-derived cultivar 'Chuanmai 104' based on PacBio HiFi reads and chromosome conformation capture sequencing. The total assembly size was 14.81 Gb with a contig N50 length of 58.25 Mb. A BUSCO analysis yielded a completeness score of 99.30%. In total, repetitive elements comprised 81.36% of the genome and 122,554 high-confidence protein-coding gene models were predicted. In summary, the first chromosome-level assembly for a SHW-derived cultivar presents a promising outlook for the study and utilization of SHWs in wheat improvement, which is essential to meet the global food demand.

Genome-wide association study and linkage mapping reveal TaqW-6B associated with water-extractable arabinoxylan content in wheat grain

KEY MESSAGE

A novel QTL, TaqW-6B of water-extractable arabinoxylan content in the wheat grain on chromosome 6BL was identified and fine mapped in a narrow region 3.8 Mb. Water-extractable arabinoxylan (WE-AX), an important component of hemicellulose, is associated with various abundant health benefits. In this study, QTLs for WE-AX content were detected in two populations: (1) a recombinant inbred line (RIL) population with 164 lines derived from a cross between Avocet and Chilero (AC population) genotyped with diversity array technology (DArT), and (2) a natural population of 243 varieties (CH population) genotyped with the Axiom wheat 660 K single-nucleotide polymorphism (SNP) array. A stable QTL Qwe-ax.haust-6B, explaining 8.51-15.59% of the phenotypic variance, was mapped in the physical interval 459.38-572.09 Mb on the long arm of chromosome 6B in the AC population, tightly linked with DArT markers 3,944,740 and 4,991,038 under three experimental conditions. The Qwe-ax.haust-6B was further narrowed down to be delimited in the physical interval 516.47-571.58 Mb on chromosome 6BL, explaining 5.86-16.27% of the phenotypic variance in the CH population. Furthermore, we developed high-throughput kompetitive allele-specific PCR (KASP) markers to reconstruct the genetic linkage map in the AC population, and Qwe-ax.haust-6B was fine mapped into a narrow region named TaqW-6B, which was compressed between KASP-6B-3 and KASP-6B-6 at a physical distance of 3.8 Mb. In the meanwhile, the markers were also validated in a natural population of 160 wheat lines (NP population). Consequently, this study is of great importance to provide the theoretical basis for cloning the key gene and developing functional markers for molecular breeding.

Genome-wide association study reveals 18 QTL for major agronomic traits in a Nordic-Baltic spring wheat germplasm

Spring wheat (Triticum aestivum L.) remains an important alternative to winter wheat cultivation at Northern latitudes due to high risk of overwintering or delayed sowing of winter wheat. We studied nine major agronomic traits in a set of 299 spring wheat genotypes in trials across 12-year-site combinations in Lithuania, Latvia, Estonia, and Norway for three consecutive years. The dataset analyzed here consisted of previously published phenotypic data collected in 2021 and 2022, supplemented with additional phenotypic data from the 2023 field season collected in this study. We combined these phenotypic datasets with previously published genotypic data generated using a 25K single nucleotide polymorphism (SNP) array that yielded 18,467 markers with a minor allele frequency above 0.05. Analysis of these datasets via genome-wide association study revealed 18 consistent quantitative trait loci (QTL) replicated in two or more trials that explained more than 5% of phenotypic variance for plant height, grain protein content, thousand kernel weight, or heading date. The most consistent markers across the tested environments were detected for plant height, thousand kernel weight, and days to heading in eight, five, and six trials, respectively. No beneficial effect of the semi-dwarfing alleles Rht-B1b and Rht-D1b on grain yield performance was observed across the 12 tested trials. Moreover, the cultivars carrying these alleles were low yielding in general. Based on principal component analysis, wheat genotypes developed in the Northern European region clustered separately from those developed at the southern latitudes, and markers associated with the clustering were identified. Important phenotypic traits, such as grain yield, days to heading, grain protein content, and thousand kernel weight were associated with this clustering of the genotype sets. Interestingly, despite being adapted to the Nordic environment, genotypes in the Northern set demonstrated lower grain yield performance across all tested environments. The results indicate that spring wheat germplasm harbors valuable QTL/alleles, and the identified trait-marker associations might be useful in improving Nordic-Baltic spring wheat germplasm under global warming conditions.

Characterising Biological and Physiological Drought Signals in Diverse Parents of a Wheat Mapping Population

Water deficit affects the growth as well as physiological and biochemical processes in plants. The aim of this study was to determine differences in physiological and biochemical responses to drought stress in two wheat cultivars-Chinese Spring (CS) and SQ1 (which are parents of a mapping population of doubled haploid lines)-and to relate these responses to final yield and agronomic traits. Drought stress was induced by withholding water for 14 days, after which plants were re-watered and maintained until harvest. Instantaneous gas exchange parameters were evaluated on the 3rd, 5th, 10th, and 14th days of seedling growth under drought. After 14 days, water content and levels of chlorophyll a+b, carotenoids, malondialdehyde, soluble carbohydrates, phenolics, salicylic acid, abscisic acid (ABA), and polyamines were measured. At final maturity, yield components (grain number and weight), biomass, straw weight, and harvest index were evaluated. Physiological and biochemical parameters of CS responded more than those of SQ1 to the 14-day drought, reflected in a greater reduction in final biomass and yield in CS. Marked biochemical differences between responses of CS and SQ1 to the drought were found for soluble carbohydrates and polyamines. These would be good candidates for testing in the mapping population for the coincidence of the genetic control of these traits and final biomass and yield.

Diversity in bread and durum wheat stigma morphology and linkage of increased stigma length to dwarfing gene Rht14

KEY MESSAGE

The dwarfing allele Rht14 of durum wheat associates with greater stigma length, an important trait for hybrid breeding, whilst major dwarfing alleles Rht-B1b and Rht-D1b showed little to no effect. Although much understudied in wheat, the stigma is a crucial component for attaining grain set, the fundamental basis for yield, particularly in hybrid production systems where successful grain set relies on wind-driven pollen dispersal by the male parent and effective pollen capture by the female parent. Females with long stigma that exsert early are thought to be advantageous. Using glasshouse-grown lines, we examined variation in Total Stigma Length (TSL) across diverse panels comprising 27 durum and 116 bread wheat genotypes. Contrasting genotypes were selected for population development and genetic analysis. Quantitative trait loci (QTL) analysis was performed on a durum F2 population and a bread wheat recombinant inbred line (RIL) population. Contrasting with studies of anther length, we found no large effect on TSL of the GA-insensitive semi-dwarfing genes Rht-B1 and Rht-D1 in either durum or bread wheat. However, in durum cultivar Italo, we identified a region on chromosome 6A which is robustly associated with larger TSL and contains the Rht14 allele for reduced plant height, a trait that is favourable for female line development in hybrid systems. This dual effect locus explained 25.2 and 19.2% of TSL phenotypic variation in experiments across two growing seasons, with preliminary results suggesting this locus may increase TSL when transferred to bread wheat. In a bread wheat, RIL population minor QTL on 1A and 2A was indicated, but the strongest association was with Ppd-B1. Methods developed here, and the identification of a TSL-enhancing locus provides advances and further opportunities in the study of wheat stigma.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

Seed development-related genes contribute to high yield heterosis in integrated utilization of elite autotetraploid and neo-tetraploid rice

Introduction

Autotetraploid rice holds high resistance to abiotic stress and substantial promise for yield increase, but it could not be commercially used because of low fertility. Thus, our team developed neo-tetraploid rice with high fertility and hybrid vigor when crossed with indica autotetraploid rice. Despite these advances, the molecular mechanisms underlying this heterosis remain poorly understood.

Methods

An elite indica autotetraploid rice line (HD11) was used to cross with neo-tetraploid rice, and 34 hybrids were obtained to evaluate agronomic traits related to yield. WE-CLSM, RNA-seq, and CRISPR/Cas9 were employed to observe endosperm structure and identify candidate genes from two represent hybrids.

Results and discussion

These hybrids showed high seed setting and an approximately 55% increase in 1000-grain weight, some of which achieved grain yields comparable to those of the diploid rice variety. The endosperm observations indicated that the starch grains in the hybrids were more compact than those in paternal lines. A total of 119 seed heterosis related genes (SHRGs) with different expressions were identified, which might contribute to high 1000-grain weight heterosis in neo-tetraploid hybrids. Among them, 12 genes had been found to regulate grain weight formation, including OsFl3, ONAC023, OsNAC024, ONAC025, ONAC026, RAG2, FLO4, FLO11, OsISA1, OsNF-YB1, NF-YC12, and OsYUC9. Haplotype analyses of these 12 genes revealed the various effects on grain weight among different haplotypes. The hybrids could polymerize more dominant haplotypes of above grain weight regulators than any homozygous cultivar. Moreover, two SHRGs (OsFl3 and SHRG2) mutants displayed a significant reduction in 1000-grain weight and an increase in grain chalkiness, indicating that OsFl3 and SHRG2 positively regulate grain weight. Our research has identified a valuable indica autotetraploid germplasm for generating strong yield heterosis in combination with neo-tetraploid lines and gaining molecular insights into the regulatory processes of heterosis in tetraploid rice.

Molecular characterization of QTL for grain zinc and iron concentrations in wheat landrace Chinese Spring

KEY MESSAGE

Three stable QTL for grain zinc concentration were identified in wheat landrace Chinese Spring. Favorable alleles were more frequent in landraces than in modern wheat cultivars. Wheat is a major source of dietary energy for the growing world population. Developing cultivars with enriched zinc and iron can potentially alleviate human micronutrient deficiency. In this study, a recombinant inbred line (RIL) population with 245 lines derived from cross Zhou 8425B/Chinese Spring was used to detect quantitative trait loci (QTL) for grain zinc concentration (GZnC) and grain iron concentration (GFeC) across four environments. Three stable QTL for GZnC with all favorable alleles from Chinese Spring were identified on chromosomes 3BL, 5AL, and 5BL. These QTL explaining maxima of 8.7%, 5.8%, and 7.1% of phenotypic variances were validated in 125 resequenced wheat accessions encompassing both landraces and modern cultivars using six kompetitive allele specific PCR (KASP) assays. The frequencies of favorable alleles for QGZnCzc.caas-3BL, QGZnCzc.caas-5AL and QGZnCzc.caas-5BL were higher in landraces (90.4%, 68.0%, and 100.0%, respectively) compared to modern cultivars (45.9%, 35.4%, and 40.9%), suggesting they were not selected in breeding programs. Candidate gene association studies on GZnC in the cultivar panel further delimited the QTL into 8.5 Mb, 4.1 Mb, and 47.8 Mb regions containing 46, 4, and 199 candidate genes, respectively. The 5BL QTL located in a region where recombination was suppressed. Two stable and three less stable QTL for GFeC with favorable alleles also from Chinese Spring were identified on chromosomes 4BS (Rht-B1a), 4DS (Rht-D1a), 1DS, 3AS, and 6DS. This study sheds light on the genetic basis of GZnC and GFeC in Chinese Spring and provides useful molecular markers for wheat biofortification.

Characterization of putative calcium-dependent protein kinase-1 (TaCPK-1) gene: hubs in signalling and tolerance network of wheat under terminal heat

Calcium-dependent protein kinase (CDPK) is member of one of the most important signalling cascades operating inside the plant system due to its peculiar role as thermo-sensor. Here, we identified 28 full length putative CDPKs from wheat designated as TaCDPK (1-28). Based on digital gene expression, we cloned full length TaCPK-1 gene of 1691 nucleotides with open reading frame (ORF) of 548 amino acids (accession number OP125853). The expression of TaCPK-1 was observed maximum (3.1-fold) in leaf of wheat cv. HD2985 (thermotolerant) under T2 (38 ± 3 °C, 2 h), as compared to control. A positive correlation was observed between the expression of TaCPK-1 and other stress-associated genes (MAPK6, CDPK4, HSFA6e, HSF3, HSP17, HSP70, SOD and CAT) involved in thermotolerance. Global protein kinase assay showed maximum activity in leaves, as compared to root, stem and spike under heat stress. Immunoblot analysis showed abundance of CDPK protein in wheat cv. HD2985 (thermotolerant) in response to T2 (38 ± 3 °C, 2 h), as compared to HD2329 (thermosusceptible). Calcium ion (Ca2+), being inducer of CDPK, showed strong Ca-signature in the leaf tissue (Ca-622 ppm) of thermotolerant wheat cv. under heat stress, whereas it was minimum (Ca-201 ppm) in spike tissue. We observed significant variations in the ionome of wheat under HS. To conclude, TaCPK-1 plays important role in triggering signaling network and in modulation of HS-tolerance in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03989-6.

Identification of the causal mutation in early heading mutant of bread wheat (Triticum aestivum L.) using MutMap approach

In bread wheat (Triticum aestivum L.), fine-tuning the heading time is essential to maximize grain yield. Photoperiod-1 (Ppd-1) and VERNALIZATION 1 (Vrn-1) are major genes affecting photoperiod sensitivity and vernalization requirements, respectively. These genes have predominantly governed heading timing. However, Ppd-1 and Vrn-1 significantly impact heading dates, necessitating another gene that can slightly modify heading dates for fine-tuning. In this study, we developed an early heading mutant from the ethyl methanesulfonate-mutagenized population of the Japanese winter wheat cultivar "Kitahonami." MutMap analysis identified a nonsense mutation in the clock component gene Wheat PHYTOCLOCK 1/LUX ARRHYTHMO (WPCL-D1) as the probable SNP responsible for the early heading mutant on chromosome 3D. Segregation analysis using F2 and F3 populations confirmed that plants carrying the wpcl-D1 allele headed significantly earlier than those with the functional WPCL-D1. The early heading mutant exhibited increased expression levels of Ppd-1 and circadian clock genes, such as WPCL1 and LATE ELONGATED HYPOCOTYL (LHY). Notably, the transcript accumulation levels of Ppd-A1 and Ppd-D1 were influenced by the copy number of the functional WPCL1 gene. These results suggest that a loss-of-function mutation in WPCL-D1 is the causal mutation for the early heading phenotype. Adjusting the functional copy number of WPCL1 will be beneficial in fine-tuning of heading dates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01478-5.

First investigation into the genetic control of meiosis in sugarcane

The sugarcane (Saccharum spp.) genome is one of the most complex of all. Modern varieties are highly polyploid and aneuploid as a result of hybridization between Saccharum officinarum and S. spontaneum. Little research has been done on meiotic control in polyploid species, with the exception of the wheat Ph1 locus harboring the ZIP4 gene (TaZIP4-B2) which promotes pairing between homologous chromosomes while suppressing crossover between homeologs. In sugarcane, despite its interspecific origin, bivalent association is favored, and multivalents, if any, are resolved at the end of prophase I. Thus, our aim herein was to investigate the purported genetic control of meiosis in the parental species and in sugarcane itself. We investigated the ZIP4 gene and immunolocalized meiotic proteins, namely synaptonemal complex proteins Zyp1 and Asy1. The sugarcane ZIP4 gene is located on chromosome 2 and expressed more abundantly in flowers, a similar profile to that found for TaZIP4-B2. ZIP4 expression is higher in S. spontaneum a neoautopolyploid, with lower expression in S. officinarum, a stable octoploid species. The sugarcane Zip4 protein contains a TPR domain, essential for scaffolding. Its 3D structure was also predicted, and it was found to be very similar to that of TaZIP4-B2, reflecting their functional relatedness. Immunolocalization of the Asy1 and Zyp1 proteins revealed that S. officinarum completes synapsis. However, in S. spontaneum and SP80-3280 (a modern variety), no nuclei with complete synapsis were observed. Importantly, our results have implications for sugarcane cytogenetics, genetic mapping, and genomics.

Variation in TaSPL6-D confers salinity tolerance in bread wheat by activating TaHKT1;5-D while preserving yield-related traits

Na+ exclusion from above-ground tissues via the Na+-selective transporter HKT1;5 is a major salt-tolerance mechanism in crops. Using the expression genome-wide association study and yeast-one-hybrid screening, we identified TaSPL6-D, a transcriptional suppressor of TaHKT1;5-D in bread wheat. SPL6 also targeted HKT1;5 in rice and Brachypodium. A 47-bp insertion in the first exon of TaSPL6-D resulted in a truncated peptide, TaSPL6-DIn, disrupting TaHKT1;5-D repression exhibited by TaSPL6-DDel. Overexpressing TaSPL6-DDel, but not TaSPL6-DIn, led to inhibited TaHKT1;5-D expression and increased salt sensitivity. Knockout of TaSPL6-DDel in two wheat genotypes enhanced salinity tolerance, which was attenuated by a further TaHKT1;5-D knockdown. Spike development was preserved in Taspl6-dd mutants but not in Taspl6-aabbdd mutants. TaSPL6-DIn was mainly present in landraces, and molecular-assisted introduction of TaSPL6-DIn from a landrace into a leading wheat cultivar successfully improved yield on saline soils. The SPL6-HKT1;5 module offers a target for the molecular breeding of salt-tolerant crops.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

The wheat powdery mildew resistance gene Pm4 also confers resistance to wheat blast

Wheat blast, caused by the fungus Magnaporthe oryzae, threatens global cereal production since its emergence in Brazil in 1985 and recently spread to Bangladesh and Zambia. Here we demonstrate that the AVR-Rmg8 effector, common in wheat-infecting isolates, is recognized by the gene Pm4, previously shown to confer resistance to specific races of Blumeria graminis f. sp. tritici, the cause of powdery mildew of wheat. We show that Pm4 alleles differ in their recognition of different AVR-Rmg8 alleles, and some confer resistance only in seedling leaves but not spikes, making it important to select for those alleles that function in both tissues. This study has identified a gene recognizing an important virulence factor present in wheat blast isolates in Bangladesh and Zambia and represents an important first step towards developing durably resistant wheat cultivars for these regions.