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

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.

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.

Evolution of wheat blast resistance gene Rmg8 accompanied by differentiation of variants recognizing the powdery mildew fungus

Wheat blast, a devastating disease having spread recently from South America to Asia and Africa, is caused by Pyricularia oryzae (synonym of Magnaporthe oryzae) pathotype Triticum, which first emerged in Brazil in 1985. Rmg8 and Rmg7, genes for resistance to wheat blast found in common wheat and tetraploid wheat, respectively, recognize the same avirulence gene, AVR-Rmg8. Here we show that an ancestral resistance gene, which had obtained an ability to recognize AVR-Rmg8 before the differentiation of Triticum and Aegilops, has expanded its target pathogens. Molecular cloning revealed that Rmg7 was an allele of Pm4, a gene for resistance to wheat powdery mildew on 2AL, whereas Rmg8 was its homoeologue on 2BL ineffective against wheat powdery mildew. Rmg8 variants with the ability to recognize AVR-Rmg8 were distributed not only in Triticum spp. but also in Aegilops speltoides, Aegilops umbellulata and Aegilops comosa. This result suggests that the origin of resistance gene(s) recognizing AVR-Rmg8 dates back to the time before differentiation of A, B, S, U and M genomes, that is, ~5 Myr before the emergence of its current target, the wheat blast fungus. Phylogenetic analyses suggested that, in the evolutionary process thereafter, some of their variants gained the ability to recognize the wheat powdery mildew fungus and evolved into genes controlling dual resistance to wheat powdery mildew and wheat blast.

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.

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.

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.

Positional cloning and characterization reveal the role of TaSRN-3D and TaBSR1 in the regulation of seminal root number in wheat

Seminal roots play a critical role in water and nutrient absorption, particularly in the early developmental stages of wheat. However, the genes responsible for controlling SRN in wheat remain largely unknown. Genetic mapping and functional analyses identified a candidate gene (TraesCS3D01G137200, TaSRN-3D) encoding a Ser/Thr kinase glycogen synthase kinase 3 (STKc_GSK3) that regulated SRN in wheat. Additionally, experiments involving hormone treatment, nitrate absorption and protein interaction were conducted to explore the regulatory mechanism of TaSRN-3D. Results showed that the TaSRN-3D4332 allele inhibited seminal roots initiation and development, while loss-of-function mutants showed significantly higher seminal root number (SRN). Exogenous application of epi-brassinolide could increase the SRN in a HS2-allelic background. Furthermore, chlorate sensitivity and 15N uptake assays revealed that a higher number of seminal roots promoted nitrate accumulation. TaBSR1 (BIN2-related SRN Regulator 1, orthologous to OsGRF4/GL2 in rice) acts as an interactor of TaSRN-3D and promotes TaBSR1 degradation to reduce SRN. This study provides valuable insights into understanding the genetic basis and regulatory network of SRN in wheat, highlighting their roles as potential targets for root-based improvement in wheat breeding.

Genomic insights into glume pubescence in durum wheat: GWAS and haplotype analysis implicates TdELD1-1A as a candidate gene

Glume pubescence is an important morphological trait for the characterization of wheat cultivars. It shows tolerance to biotic and abiotic stresses to some extent. Hg1 (formerly named Hg) locus on chromosome 1AS controls glume pubescence in wheat. Its genetic analysis, fine-mapping and candidate gene analysis have been widely studied recently, however, the cloning of Hg1 has not yet been reported. Here, we conducted a GWAS between a dense panel of 171,103 SNPs and glume pubescence (Gp) in a durum wheat population of 145 lines, and further analyzed the candidate genes of Hg1 combined with the gene expression, functional annotation, and haplotype analysis. As a results, TRITD0Uv1G104670 (TdELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was detected as the most promising candidate gene of Hg1 for glume pubescence in durum wheat. Our findings not only contribute to a deeper understanding of its cloning and functional validation but also underscore the significance of accurate genome sequences and annotations. Additionally, our study highlights the relevance of unanchored sequences in chrUn and the application of bioinformatics analysis for gene discovery in durum wheat.

Characterization of Quantitative Trait Loci for Leaf Rust Resistance in the Uzbekistani Wheat Landrace Teremai Bugdai

Leaf rust, caused by Puccinia triticina, is a major cause of wheat yield losses globally, and novel leaf rust resistance genes are needed to enhance wheat leaf rust resistance. Teremai Bugdai is a landrace from Uzebekistan that is highly resistant to many races of P. triticina in the United States. To unravel leaf rust resistance loci in Teremai Bugdai, a recombinant inbred line (RIL) population of Teremai Bugdai × TAM 110 was evaluated for response to P. triticina race Pt54-1 (TNBGJ) and genotyped using single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS). Quantitative trait loci (QTL) analysis using 5,130 high-quality GBS-SNPs revealed three QTLs, QLr-Stars-2DS, QLr-Stars-6BL, and QLr.Stars-7BL, for leaf rust resistance in two experiments. QLr-Stars-2DS, which is either a new Lr2 allele or a new resistance locus, was delimited to an ∼19.47-Mb interval between 46.4 and 65.9 Mb on 2DS and explained 31.3 and 33.2% of the phenotypic variance in the two experiments. QLr-Stars-6BL was mapped in an ∼84.0-kb interval between 719.48 and 719.56 Mb on 6BL, accounting for 33 to 36.8% of the phenotypic variance in two experiments. QLr.Stars-7BL was placed in a 350-kb interval between 762.41 and 762.76 Mb on 7BL and explained 4.4 to 5.3% of the phenotypic variance. Nine GBS-SNPs flanking these QTLs were converted to kompetitive allele specific PCR (KASP) markers, and these markers can be used to facilitate their introgression into locally adapted wheat lines.

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.

Alleles on locus chromosome 4B from different parents confer tiller number and the yield-associated traits in wheat

Pleiotropy is frequently detected in agronomic traits of wheat (Triticum aestivum). A locus on chromosome 4B, QTn/Ptn/Sl/Sns/Al/Tgw/Gl/Gw.caas-4B, proved to show pleiotropic effects on tiller, spike, and grain traits using a recombinant inbred line (RIL) population of Qingxinmai × 041133. The allele from Qingxinmai increased tiller numbers, and the allele from line 041133 produced better performances of spike traits and grain traits. Another 52 QTL for the eight traits investigated were detected on 18 chromosomes, except for chromosomes 5D, 6D, and 7B. Several genes in the genomic interval of the locus on chromosome 4B were differentially expressed in crown and inflorescence samples between Qingxinmai and line 041133. The development of the KASP marker specific for the locus on chromosome 4B is useful for molecular marker-assisted selection in wheat breeding.

Characterization of a new greenbug resistance gene Gb9 in a synthetic hexaploid wheat

Greenbug [Schizaphis graminum (Rondani)] is a serious insect pest that not only damages cereal crops, but also transmits several destructive viruses. The emergence of new greenbug biotypes in the field makes it urgent to identify novel greenbug resistance genes in wheat. CWI 76364 (PI 703397), a synthetic hexaploid wheat (SHW) line, exhibits greenbug resistance. Evaluation of an F2:3 population from cross OK 14319 × CWI 76364 indicated that a dominant gene, designated Gb9, conditions greenbug resistance in CWI 76364. Selective genotyping of a subset of F2 plants with contrasting phenotypes by genotyping-by-sequencing identified 25 SNPs closely linked to Gb9 on chromosome arm 7DL. Ten of these SNPs were converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers for genotyping the entire F2 population. Genetic analysis delimited Gb9 to a 0.6-Mb interval flanked by KASP markers located at 599,835,668 bp (Stars-KASP872) and 600,471,081 bp (Stars-KASP881) on 7DL. Gb9 was 0.5 cM distal to Stars-KASP872 and 0.5 cM proximal to Stars-KASP881. Allelism tests indicated that Gb9 is a new greenbug resistance gene which confers resistance to greenbug biotypes C, E, H, I, and TX1. TX1 is one of the most widely virulent biotypes and has overcome most known wheat greenbug resistance genes. The introgression of Gb9 into locally adapted wheat cultivars is of economic importance, and the KASP markers developed in this study can be used to tag Gb9 in cultivar development.

Identification of candidate genes for Fusarium head blight resistance from QTLs using RIL population in wheat

Fusarium head blight (FHB) stands out as one of the most devastating wheat diseases and leads to significantly grain yield losses and quality reductions in epidemic years. Exploring quantitative trait loci (QTL) for FHB resistance is a critical step for developing new FHB-resistant varieties. We previously constructed a genetic map of unigenes (UG-Map) according to the physical positions using a set of recombinant-inbred lines (RILs) derived from the cross of 'TN18 × LM6' (TL-RILs). Here, the number of diseased spikelets (NDS) and relative disease index (RDI) for FHB resistance were investigated under four environments using TL-RILs, which were distributed across 13 chromosomes. A number of 36 candidate genes for NDS and RDI from of 19 stable QTLs were identified. The average number of candidate genes per QTL was 1.89, with 14 (73.7%), two (10.5%), and three (15.8%) QTLs including one, two, and 3-10 candidate genes, respectively. Among the 24 candidate genes annotated in the reference genome RefSeq v1.1, the homologous genes of seven candidate genes, including TraesCS4B02G227300 for QNds/Rdi-4BL-4553, TraesCS5B02G303200, TraesCS5B02G303300, TraesCS5B02G303700, TraesCS5B02G303800 and TraesCS5B02G304000 for QNds/Rdi-5BL-9509, and TraesCS7A02G568400 for QNds/Rdi-7AL-14499, were previously reported to be related to FHB resistance in wheat, barely or Brachypodium distachyon. These genes should be closely associated with FHB resistance in wheat. In addition, the homologous genes of five genes, including TraesCS1A02G037600LC for QNds-1AS-2225, TraesCS1D02G017800 and TraesCS1D02G017900 for QNds-1DS-527, TraesCS1D02G018000 for QRdi-1DS-575, and TraesCS4B02G227400 for QNds/Rdi-4BL-4553, were involved in plant defense responses against pathogens. These genes should be likely associated with FHB resistance in wheat.

Regulation of root growth and elongation in wheat

Currently, the control of rhizosphere selection on farms has been applied to achieve enhancements in phenotype, extending from improvements in single root characteristics to the dynamic nature of entire crop systems. Several specific signals, regulatory elements, and mechanisms that regulate the initiation, morphogenesis, and growth of new lateral or adventitious root species have been identified, but much more work remains. Today, phenotyping technology drives the development of root traits. Available models for simulation can support all phenotyping decisions (root trait improvement). The detection and use of markers for quantitative trait loci (QTLs) are effective for enhancing selection efficiency and increasing reproductive genetic gains. Furthermore, QTLs may help wheat breeders select the appropriate roots for efficient nutrient acquisition. Single-nucleotide polymorphisms (SNPs) or alignment of sequences can only be helpful when they are associated with phenotypic variation for root development and elongation. Here, we focus on major root development processes and detail important new insights recently generated regarding the wheat genome. The first part of this review paper discusses the root morphology, apical meristem, transcriptional control, auxin distribution, phenotyping of the root system, and simulation models. In the second part, the molecular genetics of the wheat root system, SNPs, TFs, and QTLs related to root development as well as genome editing (GE) techniques for the improvement of root traits in wheat are discussed. Finally, we address the effect of omics strategies on root biomass production and summarize existing knowledge of the main molecular mechanisms involved in wheat root development and elongation.

Future of durum wheat research and breeding: Insights from early career researchers

Durum wheat (Triticum turgidum ssp. durum) is globally cultivated for pasta, couscous, and bulgur production. With the changing climate and growing world population, the need to significantly increase durum production to meet the anticipated demand is paramount. This review summarizes recent advancements in durum research, encompassing the exploitation of existing and novel genetic diversity, exploration of potential new diversity sources, breeding for climate-resilient varieties, enhancements in production and management practices, and the utilization of modern technologies in breeding and cultivar development. In comparison to bread wheat (T. aestivum), the durum wheat community and production area are considerably smaller, often comprising many small-family farmers, notably in African and Asian countries. Public breeding programs such as the International Maize and Wheat Improvement Center (CIMMYT) and the International Center for Agricultural Research in the Dry Areas (ICARDA) play a pivotal role in providing new and adapted cultivars for these small-scale growers. We spotlight the contributions of these and others in this review. Additionally, we offer our recommendations on key areas for the durum research community to explore in addressing the challenges posed by climate change while striving to enhance durum production and sustainability. As part of the Wheat Initiative, the Expert Working Group on Durum Wheat Genomics and Breeding recognizes the significance of collaborative efforts in advancing toward a shared objective. We hope the insights presented in this review stimulate future research and deliberations on the trajectory for durum wheat genomics and breeding.

Cloning a novel reduced-height (Rht) gene TaOSCA1.4 from a QTL in wheat

Reducing plant height (PH) is one of the core contents of the "Green Revolution", which began in the 1960s in wheat. A number of 27 reduced-height (Rht) genes have been identified and a great number of quantitative trait loci (QTLs) for PH have been mapped on all 21 chromosomes. Nonetheless, only several genes regulated PH have been cloned. In this study, we found the interval of QTL QPh-1B included an EST-SSR marker swes1079. According to the sequence of swes1079, we cloned the TaOSCA1.4 gene. We developed a CAPS marker to analyze the variation across a natural population. The result showed that the PH was significantly different between the two haplotypes of TaOSCA1.4-1B under most of the 12 environments and the average values of irrigation and rainfed conditions. This result further demonstrated that TaOSCA1.4 was associated with PH. Then, we validated the TaOSCA1.4 via RNAi technology. The average PHs of the wild-type (WT), RNAi lines 1 (Ri-1) and 2 (Ri-2) were 94.6, 83.6 and 79.2 cm, respectively, with significant differences between the WT and Ri-1 and Ri-2. This result indicated that the TaOSCA1.4 gene controls PH. TaOSCA1.4 is a constitutively expressed gene and its protein localizes to the cell membrane. TaOSCA1.4 gene is a member of the OSCA gene family, which regulates intracellular Ca2+ concentration. We hypothesized that knock down mutants of TaOSCA1.4 gene reduced regulatory ability of Ca2+, thus reducing the PH. Furthermore, the cell lengths of the knock down mutants are not significantly different than that of WT. We speculate that TaOSCA1.4 gene is not directly associated with gibberellin (GA), which should be a novel mechanism for a wheat Rht gene.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Genomic insights into the origin and evolution of spelt (Triticum spelta L.) as a valuable gene pool for modern wheat breeding

Spelt (Triticum aestivum ssp. spelta) is an important wheat subspecies mainly cultivated in Europe before the 20th century that has contributed to modern wheat breeding as a valuable genetic resource. However, relatively little is known about the origins and maintenance of spelt populations. Here, using resequencing data from 416 worldwide wheat accessions, including representative spelt wheat, we demonstrate that European spelt emerged when primitive hexaploid wheat spread to the west and hybridized with pre-settled domesticated emmer, the putative maternal donor. Genomic introgression regions from domesticated emmer confer spelt's primitive morphological characters used for species taxonomy, such as tenacious glumes and later flowering. We propose a haplotype-based "spelt index" to identify spelt-type wheat varieties and to quantify utilization of the spelt gene pool in modern wheat cultivars. This study reveals the genetic basis for the establishment of the spelt wheat subspecies in a specific ecological niche and the vital role of the spelt gene pool as a unique germplasm resource in modern wheat breeding.

Natural variant of Rht27, a dwarfing gene, enhances yield potential in wheat

KEY MESSAGE

Discovery of Rht27, a dwarf gene in wheat, showed potential in enhancing grain yield by reducing plant height. Plant height plays a crucial role in crop architecture and grain yield, and semi-dwarf Reduced Height (Rht) alleles contribute to lodging resistance and were important in "Green Revolution." However, the use of these alleles is associated with some negative side effects in some environments, such as reduced coleoptile length, low nitrogen use efficiency, and reduced yield. Therefore, novel dwarf gene resources are needed to pave an alternative route to overcome these side effects. In this study, a super-dwarf mutant rht27 was obtained by the mutagenesis of G1812 (Triticum urartu, the progenitor of the A sub-genome of common wheat). Genetic analysis revealed that the dwarf phenotype was regulated by a single recessive genetic factor. The candidate region for Rht27 was narrowed to a 1.55 Mb region on chromosome 3, within which we found two potential candidate genes that showed polymorphisms between the mutant and non-mutagenized G1812. Furthermore, the natural variants and elite haplotypes of the two candidates were investigated in a natural population of common wheat. The results showed that the natural variants affect grain yield components, and the dwarf haplotypes show the potential in improving agronomic traits and grain yield. Although the mutation in Rht27 results in severe dwarf phenotype in T. urartu, the natural variants in common wheat showed desirable phenotype, which suggests that Rht27 has the potential to improve wheat yield by utilizing its weak allelic mutation or fine-tuning its expression level.

Mapping QTLs for adult-plant resistance to powdery mildew and stripe rust using a recombinant inbred line population derived from cross Qingxinmai × 041133

A recombinant inbred line (RIL) population derived from wheat landrace Qingxinmai and breeding line 041133 exhibited segregation in resistance to powdery mildew and stripe rust in five and three field tests, respectively. A 16K genotyping by target sequencing (GBTS) single-nucleotide polymorphism (SNP) array-based genetic linkage map was used to dissect the quantitative trait loci (QTLs) for disease resistance. Four and seven QTLs were identified for adult-plant resistance (APR) against powdery mildew and stripe rust. QPm.caas-1B and QPm.caas-5A on chromosomes 1B and 5A were responsible for the APR against powdery mildew in line 041133. QYr.caas-1B, QYr.caas-3B, QYr.caas-4B, QYr.caas-6B.1, QYr.caas-6B.2, and QYr.caas-7B detected on the five B-genome chromosomes of line 041133 conferred its APR to stripe rust. QPm.caas-1B and QYr.caas.1B were co-localized with the pleiotropic locus Lr46/Yr29/Sr58/Pm39/Ltn2. A Kompetitive Allele Specific Polymorphic (KASP) marker KASP_1B_668028290 was developed to trace QPm/Yr.caas.1B. Four lines pyramiding six major disease resistance loci, PmQ, Yr041133, QPm/Yr.caas-1B, QPm.caas-2B.1, QYr.caas-3B, and QPm.caas-6B, were developed. They displayed effective resistance against both powdery mildew and stripe rust at the seedling and adult-plant stages.

An Accurate Representation of the Number of bZIP Transcription Factors in the Triticum aestivum (Wheat) Genome and the Regulation of Functional Genes during Salt Stress

Climate change is dramatically increasing the overall area of saline soils around the world, which is increasing by approximately two million hectares each year. Soil salinity decreases crop yields and, thereby, makes farming less profitable, potentially causing increased poverty and hunger in many areas. A solution to this problem is increasing the salt tolerance of crop plants. Transcription factors (TFs) within crop plants represent a key to understanding salt tolerance, as these proteins play important roles in the regulation of functional genes linked to salt stress. The basic leucine zipper (bZIP) TF has a well-documented role in the regulation of salt tolerance. To better understand how bZIP TFs are linked to salt tolerance, we performed a genome-wide analysis in wheat using the Chinese spring wheat genome, which has been assembled by the International Wheat Genome Sequencing Consortium. We identified 89 additional bZIP gene sequences, which brings the total of bZIP gene sequences in wheat to 237. The majority of these 237 sequences included a single bZIP protein domain; however, different combinations of five other domains also exist. The bZIP proteins are divided into ten subfamily groups. Using an in silico analysis, we identified five bZIP genes (ABF2, ABF4, ABI5, EMBP1, and VIP1) that were involved in regulating salt stress. By scrutinizing the binding properties to the 2000 bp upstream region, we identified putative functional genes under the regulation of these TFs. Expression analyses of plant tissue that had been treated with or without 100 mM NaCl revealed variable patterns between the TFs and functional genes. For example, an increased expression of ABF4 was correlated with an increased expression of the corresponding functional genes in both root and shoot tissues, whereas VIP1 downregulation in root tissues strongly decreased the expression of two functional genes. Identifying strategies to sustain the expression of the functional genes described in this study could enhance wheat's salt tolerance.

Identification and characterization of the CRK gene family in the wheat genome and analysis of their expression profile in response to high temperature-induced male sterility

Cysteine-rich receptor-like kinases (CRKs) play many important roles during plant development, including defense responses under both biotic and abiotic stress, reactive oxygen species (ROS) homeostasis, callose deposition and programmed cell death (PCD). However, there are few studies on the involvement of the CRK family in male sterility due to heat stress in wheat (Triticum aestivum L.). In this study, a genome-wide characterization of the CRK family was performed to investigate the structural and functional attributes of the wheat CRKs in anther sterility caused by heat stress. A total of 95 CRK genes were unevenly distributed on 18 chromosomes, with the most genes distributed on chromosome 2B. Paralogous homologous genes with Ka/Ks ratios less than 1 may have undergone strong purifying selection during evolution and are more functionally conserved. The collinearity analysis results of CRK genes showed that wheat and Arabidopsis (A. thaliana), foxtail millet, Brachypodium distachyon (B. distachyon), and rice have three, 12, 15, and 11 pairs of orthologous genes, respectively. In addition, the results of the network interactions of genes and miRNAs showed that five miRNAs were in the hub of the interactions map, namely tae-miR9657b-5p, tae-miR9780, tae-miR9676-5p, tae-miR164, and tae-miR531. Furthermore, qRT-PCR validation of the six TaCRK genes showed that they play key roles in the development of the mononuclear stage anthers, as all six genes were expressed at highly significant levels in heat-stressed male sterile mononuclear stage anthers compared to normal anthers. We hypothesized that the TaCRK gene is significant in the process of high-temperature-induced sterility in wheat based on the combination of anther phenotypes, paraffin sections, and qRT-PCR data. These results improve our understanding of their relationship.

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture

KEY MESSAGE

A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Mapping and characterization of rust resistance genes Lr53 and Yr35 introgressed from Aegilops species

KEY MESSAGE

The rust resistance genes Lr53 and Yr35 were introgressed into bread wheat from Aegilops longissima or Aegilops sharonensis or their S-genome containing species and mapped to the telomeric region of chromosome arm 6BS. Wheat leaf and stripe rusts are damaging fungal diseases of wheat worldwide. Breeding for resistance is a sustainable approach to control these two foliar diseases. In this study, we used SNP analysis, sequence comparisons, and cytogenetic assays to determine that the chromosomal segment carrying Lr53 and Yr35 was originated from Ae.longissima or Ae. sharonensis or their derived species. In seedling tests, Lr53 conferred strong resistance against all five Chinese Pt races tested, and Yr35 showed effectiveness against Pst race CYR34 but susceptibility to race CYR32. Using a large population (3892 recombinant gametes) derived from plants homozygous for the ph1b mutation obtained from the cross 98M71 × CSph1b, both Lr53 and Yr35 were successfully mapped to a 6.03-Mb telomeric region of chromosome arm 6BS in the Chinese Spring reference genome v1.1. Co-segregation between Lr53 and Yr35 was observed within this large mapping population. Within the candidate region, several nucleotide-binding leucine-rich repeat genes and protein kinases were identified as candidate genes. Marker pku6B3127 was completely linked to both genes and accurately predicted the absence or presence of alien segment harboring Lr53 and Yr35 in 87 tetraploid and 149 hexaploid wheat genotypes tested. We developed a line with a smaller alien segment (< 6.03 Mb) to reduce any potential linkage drag and demonstrated that it conferred resistance levels similar to those of the original donor parent 98M71. The newly developed introgression line and closely linked PCR markers will accelerate the deployment of Lr53 and Yr35 in wheat breeding programs.

Pre-breeding of spontaneous Robertsonian translocations for density planting architecture by transferring Agropyron cristatum chromosome 1P into wheat

KEY MESSAGE

We developed T1AL·1PS and T1AS·1PL Robertsonian translocations by breakage-fusion mechanism based on wheat-A. cristatum 1P(1A) substitution line with smaller leaf area, shorter plant height, and other excellent agronomic traits Agropyron cristatum, a wild relative of wheat, is a valuable germplasm resource for improving wheat genetic diversity and yield. Our previous study confirmed that the A. cristatum chromosome 1P carries alien genes that reduce plant height and leaf size in wheat. Here, we developed T1AL·1PS and T1AS·1PL Robertsonian translocations (RobTs) by breakage-fusion mechanism based on wheat-A. cristatum 1P (1A) substitution line II-3-1c. Combining molecular markers and cytological analysis, we identified 16 spontaneous RobTs from 911 F2 individuals derived from the cross of Jimai22 and II-3-1c. Fluorescence in situ hybridization (FISH) was applied to detect the fusion structures of the centromeres in wheat and A. cristatum chromosomes. Resequencing results indicated that the chromosomal junction point was located at the physical position of Triticum aestivum chromosome 1A (212.5 Mb) and A. cristatum chromosome 1P (230 Mb). Genomic in situ hybridization (GISH) in pollen mother cells showed that the produced translocation lines could form stable ring bivalent. Introducing chromosome 1PS translocation fragment into wheat significantly increased the number of fertile tillers, grain number per spike, and grain weight and reduced the flag leaf area. However, introducing chromosome 1PL translocation fragment into wheat significantly reduced flag leaf area and plant height with a negative effect on yield components. The pre-breeding of two spontaneous RobTs T1AL·1PS and T1AS·1PL was important for wheat architecture improvement.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Diploid wild oat Avena longiglumis has nutritional and adaptive traits which are valuable for common oat (A. sativa) breeding. The combination of Illumina, Nanopore and Hi-C data allowed us to assemble a high-quality chromosome-level genome of A. longiglumis (ALO), evidenced by contig N50 of 12.68 Mb with 99% BUSCO completeness for the assembly size of 3,960.97 Mb. A total of 40,845 protein-coding genes were annotated. The assembled genome was composed of 87.04% repetitive DNA sequences. Dotplots of the genome assembly (PI657387) with two published ALO genomes were compared to indicate the conservation of gene order and equal expansion of all syntenic blocks among three genome assemblies. Two recent whole-genome duplication events were characterized in genomes of diploid Avena species. These findings provide new knowledge for the genomic features of A. longiglumis, give information about the species diversity, and will accelerate the functional genomics and breeding studies in oat and related cereal crops.

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 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.

The potassium transporter TaNHX2 interacts with TaGAD1 to promote drought tolerance via modulating stomatal aperture in wheat

Drought is a major global challenge in agriculture that decreases crop production. γ-Aminobutyric acid (GABA) interfaces with drought stress in plants; however, a mechanistic understanding of the interaction between GABA accumulation and drought response remains to be established. Here we showed the potassium/proton exchanger TaNHX2 functions as a positive regulator in drought resistance in wheat by mediating cross-talk between the stomatal aperture and GABA accumulation. TaNHX2 interacted with glutamate decarboxylase TaGAD1, a key enzyme that synthesizes GABA from glutamate. Furthermore, TaNHX2 targeted the C-terminal auto-inhibitory domain of TaGAD1, enhanced its activity, and promoted GABA accumulation under drought stress. Consistent with this, the tanhx2 and tagad1 mutants showed reduced drought tolerance, and transgenic wheat with enhanced TaNHX2 expression had a yield advantage under water deficit without growth penalty. These results shed light on the plant stomatal movement mechanism under drought stress and the TaNHX2-TaGAD1 module may be harnessed for amelioration of negative environmental effects in wheat as well as other crops.

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

FIGL1 prevents aberrant chromosome associations and fragmentation and limits crossovers in polyploid wheat meiosis

Meiotic crossovers (COs) generate genetic diversity and are crucial for viable gamete production. Plant COs are typically limited to 1-3 per chromosome pair, constraining the development of improved varieties, which in wheat is exacerbated by an extreme distal localisation bias. Advances in wheat genomics and related technologies provide new opportunities to investigate, and possibly modify, recombination in this important crop species. Here, we investigate the disruption of FIGL1 in tetraploid and hexaploid wheat as a potential strategy for modifying CO frequency/position. We analysed figl1 mutants and virus-induced gene silencing lines cytogenetically. Genetic mapping was performed in the hexaploid. FIGL1 prevents abnormal meiotic chromosome associations/fragmentation in both ploidies. It suppresses class II COs in the tetraploid such that CO/chiasma frequency increased 2.1-fold in a figl1 msh5 quadruple mutant compared with a msh5 double mutant. It does not appear to affect class I COs based on HEI10 foci counts in a hexaploid figl1 triple mutant. Genetic mapping in the triple mutant suggested no significant overall increase in total recombination across examined intervals but revealed large increases in specific individual intervals. Notably, the tetraploid figl1 double mutant was sterile but the hexaploid triple mutant was moderately fertile, indicating potential utility for wheat breeding.

QTL Mapping of Adult Plant Resistance to Powdery Mildew in Chinese Wheat Landrace Baidatou

Wheat powdery mildew, caused by the biotrophic fungus Blumeria graminis f. sp. tritici (Bgt), is one of the most devastating diseases affecting wheat throughout the world. Breeding and growing resistant wheat cultivars is one of the most economic and effective methods to control the disease, and as such, identifying and mapping the new and effective resistance genes is critical. Baidatou, a Chinese wheat landrace, shows excellent field resistance to powdery mildew. To identify the resistance gene(s) in Baidatou, 170 F7:8 recombinant inbred lines (RILs) derived from the cross Mingxian 169/Baidatou were evaluated for powdery mildew response at the adult-plant stage in the experimental fields in Yangling (YL) of Shaanxi Province and Tianshui (TS) in Gansu Province in 2019, 2020, and 2021. The relative area under disease progress curve (rAUDPC) of Mingxian 169/Baidatou F7:8 RILs indicated that the resistance of Baidatou to powdery mildew was controlled by quantitative trait loci (QTLs). Based on bulk segregation analysis combined with the 660K single nucleotide polymorphism (SNP) array and genotyping by target sequencing (16K SNP) of the entire RIL population, two QTLs, QPmbdt.nwafu-2AS and QPmbdt.nwafu-3AS, were identified, and these accounted for up to 44.5% of the phenotypic variation. One of the QTLs was located on the 3.32 cM genetic interval on wheat chromosome 2AS between the kompetitive allele-specific PCR markers AX-111012288 and AX_174233809, and another was located on the 9.6 cM genetic interval on chromosome 3AS between the SNP markers 3A_684044820 and 3A_686681822. These markers could be useful for successful breeding of powdery mildew resistance in wheat.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Population epigenetics: DNA methylation in the plant omics era

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

High-throughput field phenotyping reveals that selection in breeding has affected the phenology and temperature response of wheat in the stem elongation phase

Crop growth and phenology are driven by seasonal changes in environmental variables, with temperature as one important factor. However, knowledge about genotype-specific temperature response and its influence on phenology is limited. Such information is fundamental to improve crop models and adapt selection strategies. We measured the increase in height of 352 European winter wheat varieties in 4 years to quantify phenology, and fitted an asymptotic temperature response model. The model used hourly fluctuations in temperature to parameterize the base temperature (Tmin), the temperature optimum (rmax), and the steepness (lrc) of growth responses. Our results show that higher Tmin and lrc relate to an earlier start and end of stem elongation. A higher rmax relates to an increased final height. Both final height and rmax decreased for varieties originating from the continental east of Europe towards the maritime west. A genome-wide association study (GWAS) indicated a quantitative inheritance and a large degree of independence among loci. Nevertheless, genomic prediction accuracies (GBLUPs) for Tmin and lrc were low (r≤0.32) compared with other traits (r≥0.59). As well as known, major genes related to vernalization, photoperiod, or dwarfing, the GWAS indicated additional, as yet unknown loci that dominate the temperature response.

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.

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.

Genome-Wide Association Study to Identify Marker-Trait Associations for Seed Color in Colored Wheat (Triticum aestivum L.)

This study conducted phenotypic evaluations on a wheat F3 population derived from 155 F2 plants. Traits related to seed color, including chlorophyll a, chlorophyll b, carotenoid, anthocyanin, L*, a*, and b*, were assessed, revealing highly significant correlations among various traits. Genotyping using 81,587 SNP markers resulted in 3969 high-quality markers, revealing a genome-wide distribution with varying densities across chromosomes. A genome-wide association study using fixed and random model circulating probability unification (FarmCPU) and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK) identified 11 significant marker-trait associations (MTAs) associated with L*, a*, and b*, and chromosomal distribution patterns revealed predominant locations on chromosomes 2A, 2B, and 4B. A comprehensive annotation uncovered 69 genes within the genomic vicinity of each MTA, providing potential functional insights. Gene expression analysis during seed development identified greater than 2-fold increases or decreases in expression in colored wheat for 16 of 69 genes. Among these, eight genes, including transcription factors and genes related to flavonoid and ubiquitination pathways, exhibited distinct expression patterns during seed development, providing further approaches for exploring seed coloration. This comprehensive exploration expands our understanding of the genetic basis of seed color and paves the way for informed discussions on the molecular intricacies contributing to this phenotypic trait.

A new model construction based on the knowledge graph for mining elite polyphenotype genes in crops

Identifying polyphenotype genes that simultaneously regulate important agronomic traits (e.g., plant height, yield, and disease resistance) is critical for developing novel high-quality crop varieties. Predicting the associations between genes and traits requires the organization and analysis of multi-dimensional scientific data. The existing methods for establishing the relationships between genomic data and phenotypic data can only elucidate the associations between genes and individual traits. However, there are relatively few methods for detecting elite polyphenotype genes. In this study, a knowledge graph for traits regulating-genes was constructed by collecting data from the PubMed database and eight other databases related to the staple food crops rice, maize, and wheat as well as the model plant Arabidopsis thaliana. On the basis of the knowledge graph, a model for predicting traits regulating-genes was constructed by combining the data attributes of the gene nodes and the topological relationship attributes of the gene nodes. Additionally, a scoring method for predicting the genes regulating specific traits was developed to screen for elite polyphenotype genes. A total of 125,591 nodes and 547,224 semantic relationships were included in the knowledge graph. The accuracy of the knowledge graph-based model for predicting traits regulating-genes was 0.89, the precision rate was 0.91, the recall rate was 0.96, and the F1 value was 0.94. Moreover, 4,447 polyphenotype genes for 31 trait combinations were identified, among which the rice polyphenotype gene IPA1 and the A. thaliana polyphenotype gene CUC2 were verified via a literature search. Furthermore, the wheat gene TraesCS5A02G275900 was revealed as a potential polyphenotype gene that will need to be further characterized. Meanwhile, the result of venn diagram analysis between the polyphenotype gene datasets (consists of genes that are predicted by our model) and the transcriptome gene datasets (consists of genes that were differential expression in response to disease, drought or salt) showed approximately 70% and 54% polyphenotype genes were identified in the transcriptome datasets of Arabidopsis and rice, respectively. The application of the model driven by knowledge graph for predicting traits regulating-genes represents a novel method for detecting elite polyphenotype genes.

The genetics of Ug99 stem rust resistance in spring wheat variety 'Linkert'

Wheat stem rust caused by Puccinia graminis f. sp. tritici (Pgt) threatens wheat production worldwide. The objective of this study was to characterize wheat stem rust resistance in 'Linkert', a variety with adult plant resistance effective to emerging wheat stem rust pathogen strain Ug99. Two doubled haploid (DH) populations and one recombinant inbred line (RIL) population were developed with 'Linkert' as a stem rust resistant parent. Hard red spring wheat variety 'Forefront' and genetic stock 'LMPG' were used as stem rust susceptible parents of the DH populations. Breeding line 'MN07098-6' was used as a susceptible parent of the RIL population. Both DH and RIL populations with their parents were evaluated both at the seedling stage and in the field against Pgt races. Genotyping data of the DH populations were generated using the wheat iSelect 90k SNP assay. The RIL population was genotyped by genotyping-by-sequencing. We found QTL consistently associated with wheat stem rust resistance on chromosome 2BS for the Linkert/Forefront DH population and the Linkert/MN07098-6 RIL population both in Ethiopia and Kenya. Additional reliable QTL were detected on chromosomes 5BL (125.91 cM) and 4AL (Sr7a) for the Linkert/LMPG population in Ethiopia and Kenya. Different QTL identified in the populations reflect the importance of examining the genetics of resistance in populations derived from adapted germplasm (Forefront and MN07098-6) in addition to a genetic stock (LMPG). The associated markers in this study could be used to track and select for the identified QTL in wheat breeding programs.

A new winter wheat genetic resource harbors untapped diversity from synthetic hexaploid wheat

KEY MESSAGE

The NIAB_WW_SHW_NAM population, a large nested association mapping panel, is a useful resource for mapping QTL from synthetic hexaploid wheat that can improve modern elite wheat cultivars. The allelic richness harbored in progenitors of hexaploid bread wheat (Triticum aestivum L.) is a useful resource for addressing the genetic diversity bottleneck in modern cultivars. Synthetic hexaploid wheat (SHW) is created through resynthesis of the hybridisation events between the tetraploid (Triticum turgidum subsp. durum Desf.) and diploid (Aegilops tauschii Coss.) bread wheat progenitors. We developed a large and diverse winter wheat nested association mapping (NAM) population (termed the NIAB_WW_SHW_NAM) consisting of 3241 genotypes derived from 54 nested back-cross 1 (BC1) populations, each formed via back-crossing a different primary SHW into the UK winter wheat cultivar 'Robigus'. The primary SHW lines were created using 15 T. durum donors and 47 Ae. tauschii accessions that spanned the lineages and geographical range of the species. Primary SHW parents were typically earlier flowering, taller and showed better resistance to yellow rust infection (Yr) than 'Robigus'. The NIAB_WW_SHW_NAM population was genotyped using a single nucleotide polymorphism (SNP) array and 27 quantitative trait loci (QTLs) were detected for flowering time, plant height and Yr resistance. Across multiple field trials, a QTL for Yr resistance was found on chromosome 4D that corresponded to the Yr28 resistance gene previously reported in other SHW lines. These results demonstrate the value of the NIAB_WW_SHW_NAM population for genetic mapping and provide the first evidence of Yr28 working in current UK environments and genetic backgrounds. These examples, coupled with the evidence of commercial wheat breeders selecting promising genotypes, highlight the potential value of the NIAB_WW_SHW_NAM to variety improvement.