Allelic variation at the VRN-1 promoter region in polyploid wheat

Gene-based model to predict heading date in wheat based on allelic characterization and environmental drivers

While numerous wheat phenology prediction models are available, most of them are constrained to using variety-dependent coefficients. The overarching objective of this study was to calibrate a gene-based model to predict wheat heading date that allows breeders to select specific gene combinations that would head within the optimal window for a given environment independently of varietal genetic background. A dataset with a total of 49 Argentine wheat cultivars and two recombinant inbred lines was chosen to cover a wide range of allelic combinations for major vernalization, photoperiod, and earliness per se genes. The model was validated using independent data from an Argentine wheat trial network that includes sites from a wide latitudinal range. Ultimately, using this gene-based model, simulations were made to identify optimal gene combinations (ideotypes) × site combinations in contrasting locations. The selected model accurately predicted heading date with an overall median error of 4.6 d. This gene-based crop model for wheat phenology allowed the identification of groups of gene combinations predicted to produce heads within a low-risk window and can be adapted to predict other phenological stages based on accessible climatic information and publicly available molecular markers, facilitating its adoption in wheat-growing regions worldwide.

Allelic Variations in Phenology Genes in Club Wheat (Triticum compactum) and Their Association with Heading Date

The allelic diversity within genes controlling the vernalization requirement (VRN1) and photoperiod response (PPD1) determines the ability of wheat to adapt to a wide range of environmental conditions and influences grain yield. In this study, allelic variations at the VRN-A1, VRN-B1, VRN-D1 and PPD-D1 genes were studied for 89 accessions of Triticum compactum from different eco-geographical regions of the world. The collection was evaluated for heading date in both field and greenhouse experiments under a long photoperiod and without vernalization. Based on heading date characteristics, 52 (58.4%) of the genotypes had a spring growth habit, and all of them carried at least one dominant VRN1 allele, while 37 (41.6%) accessions had a winter growth habit and carried the triple recessive allele combination. The photoperiod-sensitive Ppd-D1b allele was detected in 85 (95.5%) accessions and the insensitive Ppd-D1a allele in four (4.5%) accessions. A total of 10 phenology gene profiles (haplotypes) were observed at four major genes in the T. compactum germplasm collection. The LSD test revealed significant differences in the mean heading date among the different spring phenology gene profiles, both in greenhouse and field conditions. In addition, 21 microsatellite markers (simple sequence repeats, SSRs) were used to assess the genetic diversity in the collection. The 21 SSR markers amplified a total of 183 alleles across all the genotypes, with a mean of 3.2 alleles per locus. The polymorphic information content ranged from 0.49 to 0.94, with a mean of 0.84. The results of this study may be useful for both T. compactum and common wheat breeding programs as a source of agronomic traits.

Enhancing the breeding gene pool of wheat using accessions in gene banks as demonstrated by the Watkins collection

KEY MESSAGE

This study demonstrates an approach to identify accessions carrying desirable alleles affecting traits of breeding interest, paving the way for the development of environmentally resilient cultivars. Recent advances in genome technology have opened opportunities to explore old germplasm collections for accessions possessing traits valuable in cultivar development. This study analyzed genotypic and phenotypic data from the pre-Green Revolution Watkins collection, to find loci linked to climate adaptation and rust resistance (leaf, stem, and stripe) using genome scans and single-nucleotide polymorphism (SNP)-level fixation index (FST) analysis. Key signature regions were revealed on chromosomes 2A and 3B. Consistent with the past findings, a highly differentiated region was detected on 6B in the cultivar subset of the collection, highlighting this region a breeding target. Chromosome-based scans located a region previously associated with RVA peak viscosity and breakdown, spanning 12.5 Mb and 13 Mb on 7D, harboring 13 genes, including the most pleiotropic TRAESCS7D02G026700 (ET2/AP2/ERF). SNP-FST analysis noticeably differentiated several loci between spring and winter wheat, including 5A:582,550,290 linked to the Vrn-A1 gene. Likewise, SNP-level FST analysis between rust resistant and susceptible accessions detected loci associated with resistance such as 1B:670,137,479-670,543,997 linked to Lr46/Yr29, and 2B:763,926,560 in a region harboring the Yr7/Yr5/YrSP gene cluster. Significant phenotypic variation was observed between the resistant and susceptible alleles at the high-FST loci indicating the role of these loci in disease resistance. At all the top ten high-FST loci, eight and thirteen accessions carried the resistance alleles for leaf and stripe rust, respectively, suggesting these accessions as potential pre-breeding candidates. This study highlights how historic germplasm collections, combined with their genetic and environmental data, can enhance breeding diversity and identify pre-breeding materials for developing resilient cultivars.

Effect of low temperature acclimation on developmental regulation of redox responses and phytohormones metabolism in lines of crosses between spring and winter wheat

Low temperature (LT) acclimation in winter wheat (Triticum aestivum L.) was related to developmental regulation of transcriptome and metabolome for balancing growth and responses. In this study, six wheat lines from the F8 generation, derived from crosses between spring wheat (Pishtaz) and winter wheat (Claire) with distinct growth habits (based on the Vrn-1 loci) were planted under field conditions. The final leaf number (FLN) and double ridge (DR) formation showed that genotypes without vernalization requirement, including Pishtaz parent, and lines 8041 and 8044 transitioned rapidly into the reproductive stage. They also had lower LT tolerance, antioxidants activity and abscisic acid (ABA) content among genotypes. In these genotypes, cytokinin (CK) and gibberellin (GA3) contents and expression levels of gibberellin 20 oxidase (GA20ox) and gibberellin 3 oxidase (GA3ox) genes, were more active than other genotypes. Facultative lines 8020 and 8025 had higher antioxidants activity and lower hydrogen peroxide (H2O2) and malondialdehyde (MDA) contents compared to spring types. Winter genotypes, including Claire parent, lines 8011 and 8015 had a strong vernalization requirement resulted in prolonged vegetative phase, accompanied by increased LT tolerance, antioxidants activity and expression of ABA biosynthetic genes, confirming that the duration of the vegetative phase plays a key role in determining wheat's winter survival capacity. Higher LT tolerance was effectively related to retarded reproductive phase, minimized redox damages through co-regulating phytohormone-metabolites under developmental periods in winter wheat.

Winter Wheat Vernalization Alleles and Freezing Tolerance at the Seedling and Jointing Stages

This study explores the relationship between allelic variation of the vernalization genes (VRN) and the freezing tolerance at the seedling and jointing stages of winter wheat growth. It provides a basis for molecular marker development for freezing tolerance breeding of winter wheat. A total of 435 wheat accessions were used to identify and evaluate the freezing tolerance at the seedling stage using field tests, while 192 wheat accessions were used to evaluate the freezing tolerance at the jointing stage in climate chamber tests. The VRN genes of the wheat accessions were detected using allele-specific markers of the VRN-A1, VRN-B1, VRN-D1 and VRN-B3 loci, and the relationship between VRN genotype and freezing tolerance at the two developmental stages was tested. There were significant differences in freezing tolerance between the wheat accessions. Assessing the freezing tolerance of 52 wheat accessions at both the seedling and jointing stages revealed no significant correlation between tolerance at these two stages. The genotypic analysis found that Vrn-D1 was the most frequent dominant allele in winter wheat, while no accession contained the dominant alleles Vrn-A1 and Vrn-B3. Notably, freezing tolerance showed stage-specific genetic regulation; seedling-stage freezing tolerance strongly correlated with vernalization gene allelic combinations (p < 0.05), whereas jointing-stage freezing tolerance exhibited no such association. The presence of all recessive alleles vrn-A1, vrn-B1, vrn-D1 and vrn-B3 was required for strong seedling-stage freezing tolerance. The VRN-D1 marker was effective for screening freezing tolerance materials under the premise that vrn-A1 and vrn-B1 alleles are recessive at winter wheat seedling stage.

Effects of five allelic variants of the wheat vernalization gene VRN-B1 on heading date and vernalization requirements

Winter wheat must undergo vernalization to flower, while spring wheat does not require vernalization. The requirement for vernalization in wheat is primarily controlled by vernalization genes. VRN-1 are the most important vernalization genes. The recessive vrn-1 alleles have a strict vernalization requirement, while dominant mutations in Vrn-1 eliminate or reduce this requirement. In this study, the near-isogenic lines for several VRN-B1 allelic variants (Vrn-B1a, Vrn-B1b, Vrn-B1c, Vrn-B1 d and vrn-B1) were generated in two winter wheat backgrounds. Under field conditions, the four dominant Vrn-B1 allelic variants (Vrn-B1a, Vrn-B1b, Vrn-B1c, and Vrn-B1 d) resulted in an advancement in the heading date by 3-5 days. Using an artificially controlled gradient vernalization treatment (4-5 ℃, ranging from 0 to 45 days with 5-day intervals), the vernalization requirements of VRN-B1 allelic variants were analyzed. The relative effects on vernalization requirements were found to be vrn-B1 > Vrn-B1a = Vrn-B1 d > Vrn-B1b = Vrn-B1c (opposite to the heading date). Gene expression analysis indicates that the earlier heading associated with the dominant Vrn-B1 allelic variants is linked to their open expression under non-vernalization conditions. There may be an expression threshold at the VRN-B1 locus that eliminates the vernalization requirement, and this threshold should be lower than the vrn-B1 levels observed under saturated vernalization conditions. Furthermore, once this hypothesized threshold is reached, there appears to be no dosage effect on VRN-B1 expression. These results deepen our understanding of wheat vernalization genes and provide a theoretical basis for utilizing these genes in breeding programs aimed at improving wheat adaptability.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-025-01565-1.

Molecular genetic regulation of the vegetative-generative transition in wheat from an environmental perspective

The key to the wide geographical distribution of wheat is its high adaptability. One of the most commonly used methods for studying adaptation is investigation of the transition between the vegetative-generative phase and the subsequent intensive stem elongation process. These processes are determined largely by changes in ambient temperature, the diurnal and annual periodicity of daylength, and the composition of the light spectrum. Many genes are involved in the perception of external environmental signals, forming a complex network of interconnections that are then integrated by a few integrator genes. This hierarchical cascade system ensures the precise occurrence of the developmental stages that enable maximum productivity. This review presents the interrelationship of molecular-genetic pathways (Earliness per se, circadian/photoperiod length, vernalization - cold requirement, phytohormonal - gibberellic acid, light perception, ambient temperature perception and ageing - miRNA) responsible for environmental adaptation in wheat. Detailed molecular genetic mapping of wheat adaptability will allow breeders to incorporate new alleles that will create varieties best adapted to local environmental conditions.

The E3 ligase TaE3V-B1 ubiquitinates proteins encoded by the vernalization gene TaVRN1 and regulates developmental processes in wheat

In wheat (Triticum aestivum), early maturity is desired to avoid the hot and dry summer season, especially in view of climate change. Here, we report that TaE3V1, a C3H2C3 RING-type E3 ligase that interacts with TaVRN1, is associated with early development. Aside from its RING domain, TaE3V1 does not harbor any domains that are conserved in other RING-type or other E3 ligase proteins. TaE3V-B1b, encoded by the functional TaE3V1 allele, interacts with and ubiquitinates TaVRN1. In contrast, TaE3V-B1a, encoded by a natural nonfunctional TaE3V1 allele, neither interacts with TaVRN1 nor has E3 ligase activity. TaE3V-B1b activity decreases with plant age under warmer temperatures, but not under the low temperatures required for vernalization. We employed a gene editing method to simultaneously inactivate the 3 homoeologous TaE3V1 genes to validate their functions. Overall, our results suggest that the naturally mutated and edited TaE3V1 alleles can accelerate wheat development and aid adaptation to warming climates.

Genome-wide characterization and evolution analysis of miniature inverted-repeat transposable elements in Barley (Hordeum vulgare)

Miniature inverted-repeat transposable elements (MITEs) constitute a class of class II transposable elements (TEs) that are abundant in plant genomes, playing a crucial role in their evolution and diversity. Barley (Hordeum vulgare), the fourth-most important cereal crop globally, is widely used for brewing, animal feed, and human consumption. However, despite their significance, the mechanisms underlying the insertion or amplification of MITEs and their contributions to barley genome evolution and diversity remain poorly understood. Through our comprehensive analysis, we identified 32,258 full-length MITEs belonging to 2,992 distinct families, accounting for approximately 0.17% of the barley genome. These MITE families can be grouped into four well-known superfamilies (Tc1/Mariner-like, PIF/Harbinger-like, hAT-like, and Mutator-like) and one unidentified superfamily. Notably, we observed two major expansion events in the barley MITE population, occurring approximately 12-13 million years ago (Mya) and 2-3 Mya. Our investigation revealed a strong preference of MITEs for gene-related regions, particularly in promoters, suggesting their potential involvement in regulating host gene expression. Additionally, we discovered that 7.73% miRNAs are derived from MITEs, thereby influencing the origin of certain miRNAs and potentially exerting a significant impact on post-transcriptional gene expression control. Evolutionary analysis demonstrated that MITEs exhibit lower conservation compared to genes, consistent with their dynamic mobility. We also identified a series of MITE insertions or deletions associated with domestication, highlighting these regions as promising targets for crop improvement strategies. These findings significantly advance our understanding of the fundamental characteristics and evolutionary patterns of MITEs in the barley genome. Moreover, they contribute to our knowledge of gene regulatory networks and provide valuable insights for crop improvement endeavors.

Identification of Quantitative Trait Loci Associated with Plant Adaptation Traits Using Nested Association Mapping Population

This study evaluated 290 recombinant inbred lines (RILs) of the nested association mapping (NAM) population from the UK. The population derived from 24 families, where a common parent was "Paragon," one of the UK's spring wheat cultivar standards. All genotypes were tested in two regions of Kazakhstan at the Kazakh Research Institute of Agriculture and Plant Industry (KRIAPI, Almaty region, Southeast Kazakhstan, 2019-2022 years) and Alexandr Barayev Scientific-Production Center for Grain Farming (SPCGF, Shortandy, Akmola region, Northern Kazakhstan, 2019-2022 years). The studied traits consisted of plant adaptation-related traits, including heading date (HD, days), seed maturation date (SMD, days), plant height (PH, cm), and peduncle length (PL, cm). In addition, the yield per m2 was analyzed in both regions. Based on a field evaluation of the population in northern and southeastern Kazakhstan and using 10,448 polymorphic SNP (single-nucleotide polymorphism) markers, the genome-wide association study (GWAS) allowed for detecting 74 QTLs in four studied agronomic traits (HD, SMD, PH, and PL). The literature survey suggested that 16 of the 74 QTLs identified in our study had also been detected in previous QTL mapping studies and GWASs for all studied traits. The results will be used for further studies related to the adaptation and productivity of wheat in breeding projects for higher grain productivity.

Genetic basis of an elite wheat cultivar Guinong 29 with harmonious improvement between multiple diseases resistance and other comprehensive traits

Fungal diseases, such as powdery mildew and rusts, significantly affect the quality and yield of wheat. Pyramiding diverse types of resistance genes into cultivars represents the preferred strategy to combat these diseases. Moreover, achieving collaborative improvement between diseases resistance, abiotic stress, quality, and agronomic and yield traits is difficult in genetic breeding. In this study, the wheat cultivar, Guinong 29 (GN29), showed high resistance to powdery mildew and stripe rust at both seedling and adult plant stages, and was susceptible to leaf rust at the seedling stage but slow resistance at the adult-plant stage. Meanwhile, it has elite agronomic and yield traits, indicating promising coordination ability among multiple diseases resistance and other key breeding traits. To determine the genetic basis of these elite traits, GN29 was tested with 113 molecular markers for 98 genes associated with diseases resistance, stress tolerance, quality, and adaptability. The results indicated that two powdery mildew resistance (Pm) genes, Pm2 and Pm21, confirmed the outstanding resistance to powdery mildew through genetic analysis, marker detection, genomic in situ hybridization (GISH), non-denaturing fluorescence in situ hybridization (ND-FISH), and homology-based cloning; the stripe rust resistance (Yr) gene Yr26 and leaf rust resistance (Lr) genes Lr1 and Lr46 conferred the stripe rust and slow leaf rust resistance in GN29, respectively. Meanwhile, GN29 carries dwarfing genes Rht-B1b and Rht-D1a, vernalization genes vrn-A1, vrn-B1, vrn-D1, and vrn-B3, which were consistent with the phenotypic traits in dwarf characteristic and semi-winter property; carries genes Dreb1 and Ta-CRT for stress tolerance to drought, salinity, low temperature, and abscisic acid (ABA), suggesting that GN29 may also have elite stress-tolerance ability; and carries two low-molecular-weight glutenin subunit genes Glu-B3b and Glu-B3bef which contributed to high baking quality. This study not only elucidated the genetic basis of the elite traits in GN29 but also verified the capability for harmonious improvement in both multiple diseases resistance and other comprehensive traits, offering valuable information for breeding breakthrough-resistant cultivars.

Identification and validation of two major QTL for grain number per spike on chromosomes 2B and 2D in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Major QTL for grain number per spike were identified on chromosomes 2B and 2D. Haplotypes and candidate genes of QGns.cib-2B.1 were analyzed. Grain number per spike (GNS) is one of the main components of wheat yield. Genetic dissection of their regulatory factors is essential to improve the yield potential. In present study, a recombinant inbred line population comprising 180 lines developed from the cross between a high GNS line W7268 and a cultivar Chuanyu12 was employed to identify quantitative trait loci (QTL) associated with GNS across six environments. Two major QTL, QGns.cib-2B.1 and QGns.cib-2D.1, were detected in at least four environments with the phenotypic variations of 12.99-27.07% and 8.50-13.79%, respectively. And significant interactions were observed between the two major QTL. In addition, QGns.cib-2B.1 is a QTL cluster for GNS, grain number per spikelet and fertile tiller number, and they were validated in different genetic backgrounds using Kompetitive Allele Specific PCR (KASP) markers. QGns.cib-2B.1 showed pleotropic effects on other yield-related traits including plant height, spike length, and spikelet number per spike, but did not significantly affect thousand grain weight which suggested that it might be potentially applicable in breeding program. Comparison analysis suggested that QGns.cib-2B.1 might be a novel QTL. Furthermore, haplotype analysis of QGns.cib-2B.1 indicated that it is a hot spot of artificial selection during wheat improvement. Based on the expression patterns, gene annotation, orthologs analysis and sequence variations, the candidate genes of QGns.cib-2B.1 were predicted. Collectively, the major QTL and KASP markers reported here provided a wealth of information for the genetic basis of GNS and grain yield improvement.

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.

Analysis of the Effects of the Vrn-1 and Ppd-1 Alleles on Adaptive and Agronomic Traits in Common Wheat (Triticum aestivum L.)

Wheat heading time is primarily governed by two loci: VRN-1 (response to vernalization) and PPD-1 (response to photoperiod). Five sets of near-isogenic lines (NILs) were studied with the aim of investigating the effect of the aforementioned genes on wheat vegetative period duration and 14 yield-related traits. Every NIL was sown in the hydroponic greenhouse of the Institute of Cytology and Genetics, SB RAS. To assess their allelic composition at the VRN-1 and PPD-1 loci, molecular markers were used. It was shown that HT in plants with the Vrn-A1vrn-B1vrn-D1 genotype was reduced by 29 and 21 days (p < 0.001) in comparison to HT in plants with the vrn-A1Vrn-B1vrn-D1 and the vrn-A1vrn-B1Vrn-D1 genotypes, respectively. In our study, we noticed a decrease in spike length as well as spikelet number per spike parameter for some NIL carriers of the Vrn-A1a allele in comparison to carriers of the Vrn-B1 allele. PCA revealed three first principal components (PC), together explaining more than 70% of the data variance. Among the studied genetic traits, the Vrn-A1a and Ppd-D1a alleles showed significant correlations with PCs. Regarding genetic components, significant correlations were calculated between PC3 and Ppd-B1a (-0.26, p < 0.05) and Vrn-B1 (0.57, p < 0.05) alleles. Thus, the presence of the Vrn-A1a allele affects heading time, while Ppd-D1a is associated with plant height reduction.

Identification of a Rye Spring Mutant Derived from a Winter Rye Variety by High-Altitude Environment Screening Using RNA Sequencing Technology

Wintergrazer-70 and Ganyin No1 are high-yield forage varieties suitable for cultivation in high-altitude areas of Tibet (4300 m above sea level). Ganyin No1 was developed from Wintergrazer-70, with the latter serving as its parent variety. Ganyin No1 was identified as a spring variety, and subsequent RNA sequencing was conducted. RNA sequencing analysis identified 4 differentially expressed genes related to vernalization and 28 genes related to photoperiod regulation. The Sc7296g5-i1G3 gene is related to the flowering inhibition of rye, which may be related to the phenotypic difference in the Ganyin No1 variety in winter and spring. This finding provides valuable insights for future research on Ganyin No1, especially in addressing feed shortages in Tibet during winter and spring.

Systematic identification of wheat spike developmental regulators by integrated multi-omics, transcriptional network, GWAS, and genetic analyses

The spike architecture of wheat plays a crucial role in determining grain number, making it a key trait for optimization in wheat breeding programs. In this study, we used a multi-omic approach to analyze the transcriptome and epigenome profiles of the young spike at eight developmental stages, revealing coordinated changes in chromatin accessibility and H3K27me3 abundance during the flowering transition. We constructed a core transcriptional regulatory network (TRN) that drives wheat spike formation and experimentally validated a multi-layer regulatory module involving TaSPL15, TaAGLG1, and TaFUL2. By integrating the TRN with genome-wide association studies, we identified 227 transcription factors, including 42 with known functions and 185 with unknown functions. Further investigation of 61 novel transcription factors using multiple homozygous mutant lines revealed 36 transcription factors that regulate spike architecture or flowering time, such as TaMYC2-A1, TaMYB30-A1, and TaWRKY37-A1. Of particular interest, TaMYB30-A1, downstream of and repressed by WFZP, was found to regulate fertile spikelet number. Notably, the excellent haplotype of TaMYB30-A1, which contains a C allele at the WFZP binding site, was enriched during wheat breeding improvement in China, leading to improved agronomic traits. Finally, we constructed a free and open access Wheat Spike Multi-Omic Database (http://39.98.48.156:8800/#/). Our study identifies novel and high-confidence regulators and offers an effective strategy for dissecting the genetic basis of wheat spike development, with practical value for wheat breeding.

The TaSOC1-TaVRN1 module integrates photoperiod and vernalization signals to regulate wheat flowering

Wheat needs different durations of vernalization, which accelerates flowering by exposure to cold temperature, to ensure reproductive development at the optimum time, as that is critical for adaptability and high yield. TaVRN1 is the central flowering regulator in the vernalization pathway and encodes a MADS-box transcription factor (TF) that usually works by forming hetero- or homo-dimers. We previously identified that TaVRN1 bound to an MADS-box TF TaSOC1 whose orthologues are flowering activators in other plants. The specific function of TaSOC1 and the biological implication of its interaction with TaVRN1 remained unknown. Here, we demonstrated that TaSOC1 was a flowering repressor in the vernalization and photoperiod pathways by overexpression and knockout assays. We confirmed the physical interaction between TaSOC1 and TaVRN1 in wheat protoplasts and in planta, and further validated their genetic interplay. A Flowering Promoting Factor 1-like gene TaFPF1-2B was identified as a common downstream target of TaSOC1 and TaVRN1 through transcriptome and chromatin immunoprecipitation analyses. TaSOC1 competed with TaVRT2, another MADS-box flowering regulator, to bind to TaVRN1; their coding genes synergistically control TaFPF1-2B expression and flowering initiation in response to photoperiod and low temperature. We identified major haplotypes of TaSOC1 and found that TaSOC1-Hap1 conferred earlier flowering than TaSOC1-Hap2 and had been subjected to positive selection in wheat breeding. We also revealed that wheat SOC1 family members were important domestication loci and expanded by tandem and segmental duplication events. These findings offer new insights into the regulatory mechanism underlying flowering control along with useful genetic resources for wheat improvement.

Allelic Variations in Vernalization (Vrn) Genes in Triticum spp

Rapid climate changes, with higher warming rates during winter and spring seasons, dramatically affect the vernalization requirements, one of the most critical processes for the induction of wheat reproductive growth, with severe consequences on flowering time, grain filling, and grain yield. Specifically, the Vrn genes play a major role in the transition from vegetative to reproductive growth in wheat. Recent advances in wheat genomics have significantly improved the understanding of the molecular mechanisms of Vrn genes (Vrn-1, Vrn-2, Vrn-3, and Vrn-4), unveiling a diverse array of natural allelic variations. In this review, we have examined the current knowledge of Vrn genes from a functional and structural point of view, considering the studies conducted on Vrn alleles at different ploidy levels (diploid, tetraploid, and hexaploid). The molecular characterization of Vrn-1 alleles has been a focal point, revealing a diverse array of allelic forms with implications for flowering time. We have highlighted the structural complexity of the different allelic forms and the problems linked to the different nomenclature of some Vrn alleles. Addressing these issues will be crucial for harmonizing research efforts and enhancing our understanding of Vrn gene function and evolution. The increasing availability of genome and transcriptome sequences, along with the improvements in bioinformatics and computational biology, offers a versatile range of possibilities for enriching genomic regions surrounding the target sites of Vrn genes, paving the way for innovative approaches to manipulate flowering time and improve wheat productivity.

Genetic dissection of major QTL for grain number per spike on chromosomes 5A and 6A in bread wheat (Triticum aestivum L.)

Grain number per spike (GNS) is a crucial component of grain yield and plays a significant role in improving wheat yield. To identify quantitative trait loci (QTL) associated with GNS, a recombinant inbred line (RIL) population derived from the cross of Zhongkemai 13F10 and Chuanmai 42 was employed to conduct QTL mapping across eight environments. Based on the bulked segregant exome sequencing (BSE-Seq), genomic regions associated with GNS were detected on chromosomes 5A and 6A. According to the constructed genetic maps, two major QTL QGns.cib-5A (LOD = 4.35-8.16, PVE = 8.46-14.43%) and QGns.cib-6A (LOD = 3.82-30.80, PVE = 5.44-12.38%) were detected in five and four environments, respectively. QGns.cib-6A is a QTL cluster for other seven yield-related traits. QGns.cib-5A and QGns.cib-6A were further validated using linked Kompetitive Allele Specific PCR (KASP) markers in different genetic backgrounds. QGns.cib-5A exhibited pleiotropic effects on productive tiller number (PTN), spike length (SL), fertile spikelet number per spike (FSN), and ratio of grain length to grain width (GL/GW) but did not significantly affect thousand grain weight (TGW). Haplotype analysis revealed that QGns.cib-5A and QGns.cib-6A were the targets of artificial selection during wheat improvement. Candidate genes for QGns.cib-5A and QGns.cib-6A were predicted by analyzing gene annotation, spatiotemporal expression patterns, and orthologous and sequence differences. These findings will be valuable for fine mapping and map-based cloning of genes underlying QGns.cib-5A and QGns.cib-6A.

A molecular mechanism for embryonic resetting of winter memory and restoration of winter annual growth habit in wheat

The staple food crop winter bread wheat (Triticum aestivum) acquires competence to flower in late spring after experiencing prolonged cold in temperate winter seasons, through the physiological process of vernalization. Prolonged cold exposure results in transcriptional repression of the floral repressor VERNALIZATION 2 (TaVRN2) and activates the expression of the potent floral promoter VERNALIZATION 1 (TaVRN1). Cold-induced TaVRN1 activation and TaVRN2 repression are maintained in post-cold vegetative growth and development, leading to an epigenetic 'memory of winter cold', enabling spring flowering. When and how the cold memory is reset in wheat is essentially unknown. Here we report that the cold-induced TaVRN1 activation is inherited by early embryos, but reset in subsequent embryo development, whereas TaVRN2 remains silenced through seed development, but is reactivated rapidly by light during seed germination. We further found that a chromatin reader mediates embryonic resetting of TaVRN1 and that chromatin modifications play an important role in the regulation of TaVRN1 expression and thus the floral transition, in response to developmental state and environmental cues. The findings define a two-step molecular mechanism for re-establishing vernalization requirement in common wheat, ensuring that each generation must experience winter cold to acquire competence to flower in spring.

Evaluation of the Allelic Variations in Vernalisation (VRN1) and Photoperiod (PPD1) Genes and Genetic Diversity in a Spanish Spelt Wheat Collection

Allelic variation within genes controlling the vernalisation requirement (VRN1) and photoperiod response (PPD1) determines the adaptation of wheat to different environmental growing conditions as well as influences other traits related to grain yield. This study aimed to screen a Spanish spelt wheat collection using gene-specific molecular markers for VRN-A1, VRN-B1, VRN-D1, and PPD-D1 loci and to phenotype for heading date (HD) in both field and greenhouse experiments under a long photoperiod and without vernalisation. Fifty-five spelt genotypes (91.7%) exhibited a spring growth habit, and all of them carried at least one dominant VRN1 allele, whereas five (8.3%) genotypes had a winter growth habit, and they carried the triple recessive allele combination. The Vrn-D1s was the most frequent allele in the studied set of spelt accessions, and it was found in combination with both the dominant Vrn-A1b and/or Vrn-B1a alleles in 88.3% of the spelt accessions tested. All spelt accessions carried the photoperiod-sensitive Ppd-D1b allele, which may explain the late heading of spelt germplasm compared to the commercial spring bread wheat Setenil used as a control. The least significant difference test showed significant differences between allelic combinations, the earliest accessions being those carrying two or three dominant alleles, followed by the one-gene combinations. In addition, the genetic diversity was evaluated through capillary electrophoresis using 15 wheat simple sequence repeat (SSR) markers. Most markers had high levels of polymorphism, producing 95 different alleles which ranged between 53 and 279 bp in size. Based on the polymorphic information content values obtained (from 0.51 to 0.97), 12 out of the 15 SSRs were catalogued as informative markers (values > 0.5). According to the dendrogram generated, the spelt accessions clustered as a separate group from the commercial bread wheat Setenil. Knowledge of VRN1 and PPD1 alleles, heading time, and genetic variability using SSR markers is valuable for spelt wheat breeding programs.

Deciphering spike architecture formation towards yield improvement in wheat

Wheat is the most widely grown crop globally, providing 20% of the daily consumed calories and protein content around the world. With the growing global population and frequent occurrence of extreme weather caused by climate change, ensuring adequate wheat production is essential for food security. The architecture of the inflorescence plays a crucial role in determining the grain number and size, which is a key trait for improving yield. Recent advances in wheat genomics and gene cloning techniques have improved our understanding of wheat spike development and its applications in breeding practices. Here, we summarize the genetic regulation network governing wheat spike formation, the strategies used for identifying and studying the key factors affecting spike architecture, and the progress made in breeding applications. Additionally, we highlight future directions that will aid in the regulatory mechanistic study of wheat spike determination and targeted breeding for grain yield improvement.

Analysis of the Structural Organization and Expression of the Vrn-D1 Gene Controlling Growth Habit (Spring vs. Winter) in Aegilops tauschii Coss

The duration of the vegetative period is an important agronomic characteristic of cereal crops. It is mainly influenced by the Vrn (response to vernalization) and Ppd (response to photoperiod) genes. In this work, we searched for alleles of several known genes of these two systems of response to external conditions in 15 accessions of Aegilops tauschii Coss. (syn. Ae. squarrosa L.), with the aim of studying the impact these alleles have on the vegetative period duration and growth habit. As a result, three allelic variants have been found for the Vrn-D1 gene: (i) one intact (winter type), (ii) one with a 5437 bp deletion in the first intron and (iii) one previously undescribed allele with a 3273 bp deletion in the first intron. It has been shown that the spring growth habit of Ae. tauschii can be developed due to the presence of a new allele of the Vrn-D1 gene. Significant differences in expression levels between the new allelic variant of the Vrn-D1 gene and the intact allele vrn-D1 were confirmed by qPCR. The new allele can be introgressed into common wheat to enhance the biodiversity of the spring growth habit and vegetative period duration of plants.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Deciphering genetic basis of developmental and agronomic traits by integrating high-throughput optical phenotyping and genome-wide association studies in wheat

Dissecting the genetic basis of complex traits such as dynamic growth and yield potential is a major challenge in crops. Monitoring the growth throughout growing season in a large wheat population to uncover the temporal genetic controls for plant growth and yield-related traits has so far not been explored. In this study, a diverse wheat panel composed of 288 lines was monitored by a non-invasive and high-throughput phenotyping platform to collect growth traits from seedling to grain filling stage and their relationship with yield-related traits was further explored. Whole genome re-sequencing of the panel provided 12.64 million markers for a high-resolution genome-wide association analysis using 190 image-based traits and 17 agronomic traits. A total of 8327 marker-trait associations were detected and clustered into 1605 quantitative trait loci (QTLs) including a number of known genes or QTLs. We identified 277 pleiotropic QTLs controlling multiple traits at different growth stages which revealed temporal dynamics of QTLs action on plant development and yield production in wheat. A candidate gene related to plant growth that was detected by image traits was further validated. Particularly, our study demonstrated that the yield-related traits are largely predictable using models developed based on i-traits and provide possibility for high-throughput early selection, thus to accelerate breeding process. Our study explored the genetic architecture of growth and yield-related traits by combining high-throughput phenotyping and genotyping, which further unravelled the complex and stage-specific contributions of genetic loci to optimize growth and yield in wheat.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

GWAS for Early-Establishment QTLs and Their Linkage to Major Phenology-Affecting Genes (Vrn, Ppd, and Eps) in Bread Wheat

Farmers in northern and central Indian regions prefer to plant wheat early in the season to take advantage of the remaining soil moisture. By planting crops before the start of the season, it is possible to extend the time frame for spring wheat. The early-wheat-establishment experiment began in the 2017 growing season at the Borlaug Institute for South Asia (BISA) in Ludhiana, India, and, after three years of intensive study, numerous agronomic, physiological, and yield data points were gathered. This study aimed to identify wheat lines suitable for early establishment through an analysis of the agro-morphological traits and the genetic mapping of associated genes or quantitative trait loci (QTLs). Advancing the planting schedule by two-three weeks proved to be advantageous in terms of providing a longer duration for crop growth and reducing the need for irrigation. This is attributed to the presence of residual soil moisture resulting from the monsoon season. Early sowing facilitated the selection of genotypes able to withstand early elevated temperatures and a prolonged phenological period. The ideotype, which includes increased photo-growing degree days for booting and heading, as well as a longer grain-filling period, is better suited to early planting than timely planting. Senescence was delayed in combination with a slower rate of canopy temperature rise, which was an excellent trait for early-adapted ideotypes. Thus, a novel approach to wheat breeding would include a screening of genotypes for early planting and an ideotype design with consistent and appropriate features. A genome-wide association study (GWAS) revealed multiple QTLs linked to early adaptation in terms of the yield and its contributing traits. Among them, 44 novel QTLs were also found along with known loci. Furthermore, the study discovered that the phenology regulatory genes, such as Vrn and Ppd, are in the same genomic region, thereby contributing to early adaptation.

Fine mapping and genetic analysis identified a C2H2-type zinc finger as a candidate gene for heading date regulation in wheat

KEY MESSAGE

A minor-effect QTL, Qhd.2AS, that affects heading date in wheat was mapped to a genomic interval of 1.70-Mb on 2AS, and gene analysis indicated that the C2H2-type zinc finger protein gene TraesCS2A02G181200 is the best candidate for Qhd.2AS. Heading date (HD) is a complex quantitative trait that determines the regional adaptability of cereal crops, and identifying the underlying genetic elements with minor effects on HD is important for improving wheat production in diverse environments. In this study, a minor QTL for HD that we named Qhd.2AS was detected on the short arm of chromosome 2A by Bulked Segregant Analysis and validated in a recombinant inbred population. Using a segregating population of 4894 individuals, Qhd.2AS was further delimited to an interval of 0.41 cM, corresponding to a genomic region spanning 1.70 Mb (from 138.87 to 140.57 Mb) that contains 16 high-confidence genes based on IWGSC RefSeq v1.0. Analyses of sequence variations and gene transcription indicated that TraesCS2A02G181200, which encodes a C2H2-type zinc finger protein, is the best candidate gene for Qhd.2AS that influences HD. Screening a TILLING mutant library identified two mutants with premature stop codons in TraesCS2A02G181200, both of which exhibited a delay in HD of 2-4 days. Additionally, variations in its putative regulatory sites were widely present in natural accession, and we also identified the allele which was positively selected during wheat breeding. Epistatic analyses indicated that Qhd.2AS-mediated HD variation is independent of VRN-B1 and environmental factors. Phenotypic investigation of homozygous recombinant inbred lines (RILs) and F2:3 families showed that Qhd.2AS has no negative effect on yield-related traits. These results provide important cues for refining HD and therefore improving yield in wheat breeding programs and will deepen our understanding of the genetic regulation of HD in cereal plants.

Long-Term In Situ Conservation Drove Microevolution of Solina d'Abruzzo Wheat on Adaptive, Agronomic and Qualitative Traits

Solina is an example of a bread wheat landrace that has been conserved in situ for centuries in Central Italy. A core collection of Solina lines sampled in areas at different altitudes and climatic conditions was obtained and genotyped. A clustering analysis based on a wide SNP dataset generated from DArTseq analysis outlined the existence of two main groups, which, after F_st analysis, showed polymorphism in genes associated with vernalization and photoperiod response. Starting from the hypothesis that the different pedoclimatic environments in which Solina lines were conserved may have shaped the population, some phenotypic characteristics were studied in the Solina core collection. Growth habit, low-temperature resistance, allelic variations at major loci involved in vernalization response, and sensitivity to photoperiod were evaluated, together with seed morphologies, grain colour, and hardness. The two Solina groups showed different responses to low temperatures and to photoperiod-specific allelic variations as well as the different morphology and technological characteristics of the grain. In conclusion, the long-term in situ conservation of Solina in environments sited at different altitudes has had an impact on the evolution of this landrace which, despite its high genetic diversity, remains clearly identifiable and distinct so as to be included in conservation varieties.

Identification of QTL for reducing loss of grain yield under salt stress conditions in bi-parental populations derived from wheat landrace Hongmangmai

A novel QTL (QSt.nftec-2BL) was mapped to a 0.7 cM interval on chromosome 2B. Plants carrying QSt.nftec-2BL produced higher grain yields by up to 21.4% than otherwise in salinized fields. Wheat yield has been limited by soil salinity in many wheat-growing areas globally. The wheat landrace Hongmangmai (HMM) possesses salt tolerance as it produced higher grain yields than other tested wheat varieties including Early Premium (EP) under salt stresses. To detect QTL underlying this tolerance, wheat cross EP × HMM was chosen to serve as mapping population that was homozygous at Ppd (photoperiod response gene), Rht (reduced plant height gene) and Vrn (vernalization gene); thus, interference with QTL detection by these loci could be minimized. QTL mapping was conducted firstly using 102 recombinant inbred lines (RILs) that were selected from the EP × HMM population (827 RILs) for similarity in grain yield under non-saline condition. Under salt stresses, however, the 102 RILs varied significantly in grain yield. These RILs were genotyped using a 90 K SNP (single nucleotide polymorphism) array; consequently, a QTL (QSt.nftec-2BL) was detected on chromosome 2B. Then, using 827 RILs and new simple sequence repeat (SSR) markers developed according to the reference sequence IWGSC RefSeq v1.0, location of QSt.nftec-2BL was refined to a 0.7 cM (6.9 Mb) interval flanked by SSR markers 2B-557.23 and 2B-564.09. Selection for QSt.nftec-2BL was performed based on the flanking markers using two bi-parental wheat populations. Trials for validating effectiveness of the selection were conducted in salinized fields in two geographical areas and two crop seasons, demonstrating that wheat plants with the salt-tolerant allele in homozygous status at QSt.nftec-2BL produced higher grain yields by up to 21.4% than otherwise.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Wild emmer wheat, the progenitor of modern bread wheat, exhibits great diversity in the VERNALIZATION1 gene

Wild emmer wheat is an excellent reservoir of genetic variability that can be utilized to improve cultivated wheat to address the challenges of the expanding world population and climate change. Bearing this in mind, we have collected a panel of 263 wild emmer wheat (WEW) genotypes across the Fertile Crescent. The genotypes were grown in different locations and phenotyped for heading date. Genome-wide association mapping (GWAS) was carried out, and 16 SNPs were associated with the heading date. As the flowering time is controlled by photoperiod and vernalization, we sequenced the VRN1 gene, the most important of the vernalization response genes, to discover new alleles. Unlike most earlier attempts, which characterized known VRN1 alleles according to a partial promoter or intron sequences, we obtained full-length sequences of VRN-A1 and VRN-B1 genes in a panel of 95 wild emmer wheat from the Fertile Crescent and uncovered a significant sequence variation. Phylogenetic analysis of VRN-A1 and VRN-B1 haplotypes revealed their evolutionary relationships and geographic distribution in the Fertile Crescent region. The newly described alleles represent an attractive resource for durum and bread wheat improvement programs.

Contemplation on wheat vernalization

Vernalization is a period of low non-freezing temperatures, which provides the competence to flower. This mechanism ensures that plants sown before winter develop reproductive organs in more favourable conditions during spring. Such an evolutionary mechanism has evolved in both monocot and eudicot plants. Studies in monocots, represented by temperate cereals like wheat and barley, have identified and proposed the VERNALIZATION1 (VRN1) gene as a key player in the vernalization response. VRN1 belongs to MADS-box transcription factors and is expressed in the leaves and the apical meristem, where it subsequently promotes flowering. Despite substantial research advancement in the last two decades, there are still gaps in our understanding of the vernalization mechanism. Here we summarise the present knowledge of wheat vernalization. We discuss VRN1 allelic variation, review vernalization models, talk VRN1 copy number variation and devernalization phenomenon. Finally, we suggest possible future directions of the vernalization research in wheat.

Allelic variations of Vrn-1 and Ppd-1 genes in Japanese wheat varieties reveal the genotype-environment interaction for heading time

The timing of heading is largely affected by environmental conditions. In wheat, Vrn-1 and Ppd-1 have been identified as the major genes involved in vernalization requirement and photoperiod sensitivity, respectively. To compare the effects of Vrn-1 and Ppd-1 alleles on heading time under different environments, we genotyped Vrn-1 and Ppd-1 homoeologues and measured the heading time at Morioka, Tsukuba and Chikugo in Japan for two growing seasons. A total of 128 Japanese and six foreign varieties, classified into four populations based on the 519 genome-wide SNPs, were used for analysis. Varieties with the spring alleles (Vrn-D1a or Vrn-D1b) at the Vrn-D1 locus and insensitive allele (Hapl-I) at the Ppd-D1 locus were found in earlier heading varieties. The effects of Vrn-D1 and Ppd-D1 on heading time were stronger than those of the other Vrn-1 and Ppd-1 homoeologues. Analysis of variance revealed that heading time was significantly affected by the genotype-environment interactions. Some Vrn-1 and Ppd-1 alleles conferred earlier or later heading in specific environments, indicating that the effect of both alleles on the timing of heading depends on the environment. Information on Vrn-1 and Ppd-1 alleles, together with heading time in various environments, provide useful information for wheat breeding.

Allelic Variations in Phenology Genes of Eastern U.S. Soft Winter and Korean Winter Wheat and Their Associations with Heading Date

Wheat heading time is genetically controlled by phenology genes including vernalization (Vrn), photoperiod (Ppd) and earliness per se (Eps) genes. Characterization of the existing genetic variation in the phenology genes of wheat would provide breeding programs with valuable genetic resources necessary for the development of wheat varieties well-adapted to the local environment and early-maturing traits suitable for double-cropping system. One hundred forty-nine eastern U.S. soft winter (ESW) and 32 Korean winter (KW) wheat genotypes were characterized using molecular markers for Vrn, Ppd, Eps and reduced-height (Rht) genes, and phenotyped for heading date (HD) in the eastern U.S. region. The Ppd-D1 and Rht-D1 genes exhibited the highest genetic diversity in ESW and KW wheat, respectively. The genetic variations for HD of ESW wheat were largely contributed by Ppd-B1, Ppd-D1 and Vrn-D3 genes. The Rht-D1 gene largely contributed to the genetic variation for HD of KW wheat. KW wheat headed on average 14 days earlier than ESW wheat in each crop year, largely due to the presence of the one-copy vrn-A1 allele in the former. The development of early-maturing ESW wheat varieties could be achieved by selecting for the one-copy vrn-A1 and vrn-D3a alleles in combination with Ppd-B1a and Ppd-D1a photoperiod insensitive alleles.

Advancing understanding of oat phenology for crop adaptation

Oat (Avena sativa) is an annual cereal grown for forage, fodder and grain. Seasonal flowering behaviour, or phenology, is a key contributor to the success of oat as a crop. As a species, oat is a vernalization-responsive long-day plant that flowers after winter as days lengthen in spring. Variation in both vernalization and daylength requirements broadens adaptation of oat and has been used to breed modern cultivars with seasonal flowering behaviours suited to different regions, sowing dates and farming practices. This review examines the importance of variation in oat phenology for crop adaptation. Strategies to advance understanding of the genetic basis of oat phenology are then outlined. These include the potential to transfer knowledge from related temperate cereals, particularly wheat (Triticum aestivum) and barley (Hordeum vulgare), to provide insights into the potential molecular basis of variation in oat phenology. Approaches that use emerging genomic resources to directly investigate the molecular basis of oat phenology are also described, including application of high-resolution genome-wide diversity surveys to map genes linked to variation in flowering behaviour. The need to resolve the contribution of individual phenology genes to crop performance by developing oat genetic resources, such as near-isogenic lines, is emphasised. Finally, ways that deeper knowledge of oat phenology can be applied to breed improved varieties and to inform on-farm decision-making are outlined.

Genetic incorporation of genes for the optimal plant architecture in common wheat

Wheat grain yield is affected by plant height, which is the total length of spike, the uppermost internode, and other elongated internodes. In this study, a population of recombinant inbred lines generated from a cross between two advanced winter wheat breeding lines were phenotyped over four locations/years and genotyped by using markers of genotyping-by-sequencing (GBS) and Diversity Array Technology (DArT) for mapping of genes for three traits, including spike length, the uppermost internode length, and plant height. Five genomic regions or quantitative trait loci (QTLs) were associated with candidate genes for these traits. A major QTL was associated with Q5A, and two novel haplotypes of Q5A were identified, one for a single nucleotide polymorphism (SNP) at position -2,149 in promoter region and the other for copy number variation. Compared with one copy Q5A on chromosome 5A in Chinese Spring, the novel haplotype of Q5A with two copies Q5A was found to generate spikes that are extremely compacted. A major QTL was associated with allelic variation in the recessive vrn-A1 alleles involving in protein sequences, and this QTL was associated with increased uppermost internode length but not with plant height. A major QTL for plant height was associated with Rht-B1b on chromosome 4B, but its effects could be compromised by two new minor QTLs on chromosome 7. Collectively, the favorable alleles from the four loci can be used to establish the optimal plant height in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01336-2.

Genetic networks underlying salinity tolerance in wheat uncovered with genome-wide analyses and selective sweeps

A genetic framework underpinning salinity tolerance at reproductive stage was revealed by genome-wide SNP markers and major adaptability genes in synthetic-derived wheats, and trait-associated loci were used to predict phenotypes. Using wild relatives of crops to identify genes related to improved productivity and resilience to climate extremes is a prioritized area of crop genetic improvement. High salinity is a widespread crop production constraint, and development of salt-tolerant cultivars is a sustainable solution. We evaluated a panel of 294 wheat accessions comprising synthetic-derived wheat lines (SYN-DERs) and modern bread wheat advanced lines under control and high salinity conditions at two locations. The GWAS analysis revealed a quantitative genetic framework of more than 200 loci with minor effect underlying salinity tolerance at reproductive stage. The significant trait-associated SNPs were used to predict phenotypes using a GBLUP model, and the prediction accuracy (r²) ranged between 0.57 and 0.74. The r² values for flag leaf weight, days to flowering, biomass, and number of spikes per plant were all above 0.70, validating the phenotypic effects of the loci discovered in this study. Furthermore, the germplasm sets were compared to identify selection sweeps associated with salt tolerance loci in SYN-DERs. Six loci associated with salinity tolerance were found to be differentially selected in the SYN-DERs (12.4 Mb on chromosome (chr)1B, 7.1 Mb on chr2A, 11.2 Mb on chr2D, 200 Mb on chr3D, 600 Mb on chr6B, and 700.9 Mb on chr7B). A total of 228 reported markers and genes, including 17 well-characterized genes, were uncovered using GWAS and EigenGWAS. A linkage disequilibrium (LD) block on chr5A, including the Vrn-A1 gene at 575 Mb and its homeologs on chr5D, were strongly associated with multiple yield-related traits and flowering time under salinity stress conditions. The diversity panel was screened with more than 68 kompetitive allele-specific PCR (KASP) markers of functional genes in wheat, and the pleiotropic effects of superior alleles of Rht-1, TaGASR-A1, and TaCwi-A1 were revealed under salinity stress. To effectively utilize the extensive genetic information obtained from the GWAS analysis, a genetic interaction network was constructed to reveal correlations among the investigated traits. The genetic network data combined with GWAS, selective sweeps, and the functional gene survey provided a quantitative genetic framework for identifying differentially retained loci associated with salinity tolerance in wheat.

Rapid identification and deployment of major genes for flowering time and awn traits in common wheat

Molecular markers are developed to accelerate deployment of genes for desirable traits segregated in a bi-parental population of recombinant inbred lines (RILs) or doubled haplotype (DH) lines for mapping. However, it would be the most effective if such markers for multiple traits could be identified in an F₂ population. In this study, single nucleotide polymorphisms (SNP) chips were used to identify major genes for heading date and awn in an F₂ population without developing RILs or DH lines. The population was generated from a cross between a locally adapted spring wheat cultivar "Ningmaizi119" and a winter wheat cultivar "Tabasco" with a diverse genetic background. It was found that the dominant Vrn-D1 allele could make Ningmaizi119 flowered a few months earlier than Tabasco in the greenhouse and without vernalization. The observed effects of the allele were validated in F₃ populations. It was also found that the dominant Ali-A1 allele for awnless trait in Tabasco or the recessive ali-A1 allele for awn trait in Ningmaizi119 was segregated in the F₂ population. The allelic variation in the ALI-A1 gene relies not only on the DNA polymorphisms in the promoter but also on gene copy number, with one copy ali-A1 in Ningmaizi119 but two copies Ali-A1 in Tabasco based on RT-PCR results. According to wheat genome sequences, cultivar "Mattis" has two copies Ali-A1 and cultivar "Spelta" has four copies Ali-A in a chromosome that was uncharacterized (ChrUN), in addition to one copy on chromosome 5A. This study rapidly characterized the effects of the dominant Vrn-D1 allele and identified the haplotype of Ali-A1 in gene copy number in the F₂ segregation population of common wheat will accelerate their deployment in cycling lines in breeding.

Identification of QTLs for wheat heading time across multiple-environments

KEY MESSAGE

The genetic response to changing climatic factors selects consistent across the tested environments and location-specific thermo-sensitive and photoperiod susceptible alleles in lower and higher altitudes, respectively, for starting flowering in winter wheat. Wheat breeders select heading date to match the most favorable conditions for their target environments and this is favored by the extensive genetic variation for this trait that has the potential to be further explored. In this study, we used a germplasm with broad geographic distribution and tested it in multi-location field trials across Germany over three years. The genotypic response to the variation in the climatic parameters depending on location and year uncovered the effect of photoperiod and spring temperatures in accelerating heading date in higher and lower latitudes, respectively. Spring temperature dominates other factors in inducing heading, whereas the higher amount of solar radiation delays it. A genome-wide scan of marker-trait associations with heading date detected two QTL: an adapted allele at locus TaHd102 on chromosome 5A that has a consistent effect on HD in German cultivars in multiple environments and a non-adapted allele at locus TaHd044 on chromosome 3A that accelerates flowering by 5.6 days. TaHd102 and TaHd044 explain 13.8% and 33% of the genetic variance, respectively. The interplay of the climatic variables led to the detection of environment specific association responding to temperature in lower latitudes and photoperiod in higher ones. Another locus TaHd098 on chromosome 5A showed epistatic interactions with 15 known regulators of flowering time when non-adapted cultivars from outside Germany were included in the analysis.

Long-Amplicon Single-Molecule Sequencing Reveals Novel, Trait-Associated Variants of VERNALIZATION1 Homoeologs in Hexaploid Wheat

The gene VERNALIZATION1 (VRN1) is a key controller of vernalization requirement in wheat. The genome of hexaploid wheat (Triticum aestivum) harbors three homoeologous VRN1 loci on chromosomes 5A, 5B, and 5D. Structural sequence variants including small and large deletions and insertions and single nucleotide polymorphisms (SNPs) in the three homoeologous VRN1 genes not only play an important role in the control of vernalization requirement, but also have been reported to be associated with other yield related traits of wheat. Here we used single-molecule sequencing of barcoded long-amplicons to assay the full-length sequences (∼13 kbp plus 700 bp from the promoter sequence) of the three homoeologous VRN1 genes in a panel of 192 predominantly European winter wheat cultivars. Long read sequences revealed previously undetected duplications, insertions and single-nucleotide polymorphisms in the three homoeologous VRN1 genes. All the polymorphisms were confirmed by Sanger sequencing. Sequence analysis showed the predominance of the winter alleles vrn-A1, vrn-B1, and vrn-D1 across the investigated cultivars. Associations of SNPs and structural variations within the three VRN1 genes with 20 economically relevant traits including yield, nodal root-angle index and quality related traits were evaluated at the levels of alleles, haplotypes, and copy number variants. Cultivars carrying structural variants within VRN1 genes showed lower grain yield, protein yield and biomass compared to those with intact genes. Cultivars carrying a single vrn-A1 copy and a unique haplotype with a high number of SNPs were found to have elevated grain yield, kernels per spike and kernels per m2 along with lower grain sedimentation values. In addition, we detected a novel SNP polymorphism within the G-quadruplex region of the promoter of vrn-A1 that was associated with deeper roots in winter wheat. Our findings show that multiplex, single-molecule long-amplicon sequencing is a useful tool for detecting variants in target genes within large plant populations, and can be used to simultaneously assay sequence variants among target multiple gene homoeologs in polyploid crops. Numerous novel VRN1 haplotypes and alleles were identified that showed significantly associations to economically important traits. These polymorphisms were converted into PCR or KASP assays for use in marker-assisted breeding.

Plant Development of Early-Maturing Spring Wheat (Triticum aestivum L.) under Inoculation with Bacillus sp. V2026

The effect of a plant growth-promoting bacterium (PGPB) Bacillus sp. V2026, a producer of indolyl-3-acetic acid (IAA) and gibberellic acid (GA), on the ontogenesis and productivity of four genotypes of early-maturing spring wheat was studied under controlled conditions. The inoculation of wheat plants with Bacillus sp. V2026 increased the levels of endogenous IAA and GA in wheat of all genotypes and the level of trans-Zeatin in Sonora 64 and Leningradskaya rannyaya cvs but decreased it in AFI177 and AFI91 ultra-early lines. Interactions between the factors “genotype” and “inoculation” were significant for IAA, GA, and trans-Zeatin concentrations in wheat shoots and roots. The inoculation increased the levels of chlorophylls and carotenoids and reduced lipid peroxidation in leaves of all genotypes. The inoculation resulted in a significant increase in grain yield (by 33–62%), a reduction in the time for passing the stages of ontogenesis (by 2–3 days), and an increase in the content of macro- and microelements and protein in the grain. Early-maturing wheat genotypes showed a different response to inoculation with the bacterium Bacillus sp. V2026. Cv. Leningradskaya rannyaya was most responsive to inoculation with Bacillus sp. V2026.

Population structure, allelic variation at Rht-B1 and Ppd-A1 loci and its effects on agronomic traits in Argentinian durum wheat

Exploring the genetic variability in yield and yield-related traits is essential to continue improving genetic gains. Fifty-nine Argentinian durum wheat cultivars were analyzed for important agronomic traits in three field experiments. The collection was genotyped with 3565 genome-wide SNPs and functional markers in order to determine the allelic variation at Rht-B1 and Ppd-A1 genes. Population structure analyses revealed the presence of three main groups, composed by old, modern and genotypes with European or CIMMYT ancestry. The photoperiod sensitivity Ppd-A1b allele showed higher frequency (75%) than the insensitivity one Ppd-A1a (GS105). The semi-dwarfism Rht-B1b and the Ppd-A1a (GS105) alleles were associated with increases in harvest index and decreases in plant height, grain protein content and earlier heading date, although only the varieties carrying the Rht-B1 variants showed differences in grain yield. Out of the two main yield components, grain number per plant was affected by allelic variants at Rht-B1 and Ppd-A1 loci, while no differences were observed in thousand kernel weight. The increases in grain number per spike associated with Rht-B1b were attributed to a higher grain number per spikelet, whereas Ppd-A1a (GS105) was associated with higher grain number per spikelet, but also with lower spikelets per spike.

Developing core marker sets for effective genomic-assisted selection in wheat and barley breeding programs

Wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) are widely cultivated temperate crops. In breeding programs with these crops in Japan, effective genomic-assisted selection was performed by selecting core marker sets from thousands of genome-wide amplicon sequencing markers. The core sets consist of 768 and 960 markers for barley and wheat, respectively. These markers are distributed evenly across the genomes and effectively detect widely distributed polymorphisms in the chromosomes. The core set utility was assessed using 1,032 barley and 1,798 wheat accessions across the country. Minor allele frequency and chromosomal distributions showed that the core sets could effectively capture polymorphisms across the entire genome, indicating that the core sets are applicable to highly-related advanced breeding materials. Using the core sets, we also assessed the trait value predictability. As observed via fivefold cross-validation, the prediction accuracies of six barley traits ranged from 0.56-0.74 and 0.62 on average, and the corresponding values for eight wheat traits ranged from 0.44-0.83 and 0.65 on average. These data indicate that the established core marker sets enable breeding processes to be accelerated in a cost-effective manner and provide a strong foundation for further research on genomic selection in crops.

Emmer Wheat Eco-Geographic and Genomic Congruence Shapes Phenotypic Performance under Mediterranean Climate

Emmer wheat (Triticum turgidum ssp. dicoccum) is one of the world’s oldest domesticated crops, and it harbors a potentially rich reservoir of agronomic and nutritional quality trait variations. The growing global demand for plant-based health-food niche markets has promoted new commercial interest in ancient grains, including Emmer wheat. Although T. dicoccum can also perform well under harsh environments, its cultivation along the Mediterranean agro-ecosystems is sparse. Here, we analyze a unique tetraploid wheat collection (n = 121) representing a wide geographic range of Emmer accessions, using 9897 DArTseq markers and on-field phenotypic characterization to quantify the extent of diversity among populations and the interactions between eco-geographic, genetic, and phenotypic attributes. Population genomic inferences based on the DArTseq data indicated that the collection could be split into four distinguished clusters in accordance with their eco-geographic origin although significant phenotypic variation was observed within clusters. Superior early vegetative vigor, shorter plant height, and early phenology were observed among emmer wheat accessions from Ethiopia compared to accessions from northern regions. This adaptive advantage highlights the potential of emmer wheat as an exotic germplasm for wheat improvement through breeding. The direct integration of such germplasm into conventional or organic farming agro-systems under the Mediterranean basin climate is also discussed.

Accounting for heading date gene effects allows detection of small-effect QTL associated with resistance to Septoria nodorum blotch in wheat

In humid and temperate areas, Septoria nodorum blotch (SNB) is a major fungal disease of common wheat (Triticum aestivum L.) in which grain yield is reduced when the pathogen, Parastagonospora nodorum, infects leaves and glumes during grain filling. Foliar SNB susceptibility may be associated with sensitivity to P. nodorum necrotrophic effectors (NEs). Both foliar and glume susceptibility are quantitative, and the underlying genetics are not understood in detail. We genetically mapped resistance quantitative trait loci (QTL) to leaf and glume blotch using a double haploid (DH) population derived from the cross between the moderately susceptible cultivar AGS2033 and the resistant breeding line GA03185-12LE29. The population was evaluated for SNB resistance in the field in four successive years (2018–2021). We identified major heading date (HD) and plant height (PH) variants on chromosomes 2A and 2D, co-located with SNB escape mechanisms. Five QTL with small effects associated with adult plant resistance to SNB leaf and glume blotch were detected on 1A, 1B, and 6B linkage groups. These QTL explained a relatively small proportion of the total phenotypic variation, ranging from 5.6 to 11.8%. The small-effect QTL detected in this study did not overlap with QTL associated with morphological and developmental traits, and thus are sources of resistance to SNB.

A wiring diagram to integrate physiological traits of wheat yield potential

As crop yields are pushed closer to biophysical limits, achieving yield gains becomes increasingly challenging and will require more insight into deterministic pathways to yields. Here, we propose a wiring diagram as a platform to illustrate the interrelationships of the physiological traits that impact wheat yield potential and to serve as a decision support tool for crop scientists. The wiring diagram is based on the premise that crop yield is a function of photosynthesis (source), the investment of assimilates into reproductive organs (sinks) and the underlying processes that enable expression of both. By illustrating these linkages as coded wires, the wiring diagram can show connections among traits that may not have been apparent, and can inform new research hypotheses and guide crosses designed to accumulate beneficial traits and alleles in breeding. The wiring diagram can also serve to create an ever-richer common point of reference for refining crop models in the future.

Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat

Allopolyploidy greatly expands the range of possible regulatory interactions among functionally redundant homoeologous genes. However, connection between the emerging regulatory complexity and expression and phenotypic diversity in polyploid crops remains elusive. Here, we use diverse wheat accessions to map expression quantitative trait loci (eQTL) and evaluate their effects on the population-scale variation in homoeolog expression dosage. The relative contribution of cis- and trans-eQTL to homoeolog expression variation is strongly affected by both selection and demographic events. Though trans-acting effects play major role in expression regulation, the expression dosage of homoeologs is largely influenced by cis-acting variants, which appear to be subjected to selection. The frequency and expression of homoeologous gene alleles showing strong expression dosage bias are predictive of variation in yield-related traits, and have likely been impacted by breeding for increased productivity. Our study highlights the importance of genomic variants affecting homoeolog expression dosage in shaping agronomic phenotypes and points at their potential utility for improving yield in polyploid crops.

Identification and characterization of QTL for spike morphological traits, plant height and heading date derived from the D genome of natural and resynthetic allohexaploid wheat

QHd.cau-7D.1 for heading date was delimited into the physical interval of approximately 17.38 Mb harboring three CONSTANS-like zinc finger genes. Spike morphological traits, plant height and heading date play important roles in yield improvement of wheat. To reveal the genetic factors that controlling spike morphological traits, plant height and heading date on the D genome, we conducted analysis of quantitative traits locus (QTL) using 198 F_7:8 recombinant inbred lines (RILs) derived from a cross between the common wheat TAA10 and resynthesized allohexaploid wheat XX329 with similar AABB genomes. A total of 23 environmentally stable QTL on the D sub-genome for spike length (SL), fertile spikelet number per spike (FSN), sterile spikelet number per spike (SSN), total spikelet number per spike (TSN), spike compactness (SC), plant height (PHT) and heading date (HD) were detected, among which eight appeared to be novel QTL. Furthermore, QHd.cau-7D.1 and QPht.cau-7D.2 shared identical confidence interval and were delimited into the physical interval of approximately 17.38 Mb with 145 annotated genes, including three CONSTANS-like zinc finger genes (TraesCS7D02G209000, TraesCS7D02G213000 and TraesCS7D02G220300). This study will help elucidate the molecular mechanism of the seven traits (SL, FSN, SSN, TSN, SC, PHT and HD) and provide a potentially valuable resource for genetic improvement.

Time-course transcriptome profiling revealed the specific expression patterns of MADS-box genes associated with the distinct developmental processes between winter and spring wheat

The shoot apex is a region where new cells are produced and elongate. The developmental state of the wheat shoot apex under low temperature affects its cold resistance. In this study, the morphology of shoot apex before overwintering was characterized for 24 wheat line with different winter and spring characteristics. Our research showed that the shoot apex of autumn-sown spring wheat lines reached the temperature sensitive double-ridge stage before overwintering, whereas shoot apex of winter wheat lines are found in temperature-insensitive vegetative or elongation stages. In order to explore how gene expression is associated with shoot apex differentiation in winter and spring wheat, we used strand-specific RNA sequencing to profile the gene expression patterns at four time-points between 14 after germination and 45 days after germination in the winter wheat cultivar Dongnongdongmai No. 1 (DM1) and in the spring wheat cultivar China Spring (CS). We identified 11,848 differentially expressed genes between the two cultivars. Most up-regulated genes in CS were involved in energy metabolism and transport during the seedling stage, whereas up-regulated genes in DM1 were involved in protein and DNA synthesis. MADS-box genes affect plant growth and development. In this study, MADS-boxes with differential expression between CS and DM1 were screened and evolutionary tree analysis was conducted. During all sampling periods, CS highly expressed MADS-box genes that induce flowering promotion genes such as VRN1, VRT and AG, while lowly expressed MADS-box genes that induce flowering-inhibiting homologous genes such as SVP. TaVRN1 composition in DM1 and CS was vrn-A1, vrn-B1, and Vrn-D1b. Analysis of the sequence of TaVRN1 (TraesCS5A01G391700) from DM1 and CS revealed 5 SNP differences in the promoter regions and 3 SNP deletions in the intron regions. The expression levels of cold resistant genes in DM1 were significantly higher than those in CS at seedling stage (neither DM1 nor CS experienced cold in this study), including CBF, cold induced protein,acid desaturase and proline rich proteins. Additionally, the expression levels of auxin-related genes were significantly higher in CS than those in DM1 at 45 days after germination. Our study identified candidate genes associated with the process of differentiation of the shoot apex in winter and spring wheat at the seedling stage and also raised an internal stress tolerance model for winter wheat to endogenously anticipate the coming stressful conditions in winter.