The wheat and barley vernalization gene VRN3 is an orthologue of FT

HFR1, a bHLH Transcriptional Regulator from Arabidopsis thaliana, Improves Grain Yield, Shade and Osmotic Stress Tolerances in Common Wheat

Common wheat, Triticum aestivum, is the most widely grown staple crop worldwide. To catch up with the increasing global population and cope with the changing climate, it is valuable to breed wheat cultivars that are tolerant to abiotic or shade stresses for density farming. Arabidopsis LONG HYPOCOTYL IN FAR-RED 1 (AtHFR1), a photomorphogenesis-promoting factor, is involved in multiple light-related signaling pathways and inhibits seedling etiolation and shade avoidance. We report that overexpression of AtHFR1 in wheat inhibits etiolation phenotypes under various light and shade conditions, leading to shortened plant height and increased spike number relative to non-transgenic plants in the field. Ectopic expression of AtHFR1 in wheat increases the transcript levels of TaCAB and TaCHS as observed previously in Arabidopsis, indicating that the AtHFR1 transgene can activate the light signal transduction pathway in wheat. AtHFR1 transgenic seedlings significantly exhibit tolerance to osmotic stress during seed germination compared to non-transgenic wheat. The AtHFR1 transgene represses transcription of TaFT1, TaCO1, and TaCO2, delaying development of the shoot apex and heading in wheat. Furthermore, the AtHFR1 transgene in wheat inhibits transcript levels of PHYTOCHROME-INTERACTING FACTOR 3-LIKEs (TaPIL13, TaPIL15-1B, and TaPIL15-1D), downregulating the target gene STAYGREEN (TaSGR), and thus delaying dark-induced leaf senescence. In the field, grain yields of three AtHFR1 transgenic lines were 18.2–48.1% higher than those of non-transgenic wheat. In summary, genetic modification of light signaling pathways using a photomorphogenesis-promoting factor has positive effects on grain yield due to changes in plant architecture and resource allocation and enhances tolerances to osmotic stress and shade avoidance response.

DNA hypermethylation promotes the flowering of orchardgrass during vernalization

Vernalization, influenced by environmental factors, is an essential process associated with the productivity of temperate crops, during which epigenetic regulation of gene expression plays an important role. Although DNA methylation is one of the major epigenetic mechanisms associated with the control of gene expression, global changes in DNA methylation in the regulation of gene expression during vernalization-induced flowering of temperate plants remain largely undetermined. To characterize vernalization-associated DNA methylation dynamics, we performed whole-genome bisulfite-treated sequencing and transcriptome sequencing in orchardgrass (Dactylis glomerata) during vernalization. The results revealed that increased levels of genome DNA methylation during the early vernalization of orchardgrass were associated with transcriptional changes in DNA methyltransferase and demethylase genes. Upregulated expression of vernalization-related genes during early vernalization was attributable to an increase in mCHH in the promoter regions of these genes. Application of an exogenous DNA methylation accelerator or overexpression of orchardgrass NUCLEAR POLY(A) POLYMERASE (DgPAPS4) promoted earlier flowering, indicating that DNA hypermethylation plays an important role in vernalization-induced flowering. Collectively, our findings revealed that vernalization-induced hypermethylation is responsible for floral primordium initiation and development. These observations provide a theoretical foundation for further studies on the molecular mechanisms underlying the control of vernalization in temperate grasses.

Overexpression of TaLBD16-4D alters plant architecture and heading date in transgenic wheat

Lateral organ boundaries domain (LBD) proteins, a class of plant-specific transcription factors with a special domain of lateral organ boundaries (LOB), play essential roles in plant growth and development. However, there is little known about the functions of these genes in wheat to date. Our previous study demonstrated that TaLBD16-4D is conducive to increasing lateral root number in wheat. In the present work, we further examined important agronomical traits of the aerial part of transgenic wheat overexpressing TaLBD16-4D. Interestingly, it was revealed that overexpressing TaLBD16-4D could lead to early heading and multiple alterations of plant architecture, including decreased plant height, increased flag leaf size and stem diameter, reduced spike length and tillering number, improved spike density and grain width, and decreased grain length. Moreover, auxin-responsive experiments demonstrated that the expression of TaLBD16-4D in wild-type (WT) wheat plants showed a significant upregulation through 2,4-D treatment. TaLBD16-4D-overexpression lines displayed a hyposensitivity to 2,4-D treatment and reduced shoot gravitropic response. The expressions of a set of auxin-responsive genes were markedly different between WT and transgenic plants. In addition, overexpressing TaLBD16-4D affected the transcript levels of flowering-related genes (TaGI, TaCO1, TaHd1, TaVRN1, TaVRN2, and TaFT1). Notably, the expression of TaGI, TaCO1, TaHd1, TaVRN1, and TaFT1 displayed significant upregulation under IAA treatment. Collectively, our observations indicated that overexpressing TaLBD16-4D could affect aerial architecture and heading time possibly though participating in the auxin pathway.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Flowering time runs hot and cold

Evidence suggests that anthropogenically-mediated global warming results in accelerated flowering for many plant populations. However, the fact that some plants are late flowering or unaffected by warming, underscores the complex relationship between phase change, temperature, and phylogeny. In this review, we present an emerging picture of how plants sense temperature changes, and then discuss the independent recruitment of ancient flowering pathway genes for the evolution of ambient, low, and high temperature-regulated reproductive development. As well as revealing areas of research required for a better understanding of how past thermal climates have shaped global patterns of plasticity in plant phase change, we consider the implications for these phenological thermal responses in light of climate change.

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.

Expression Patterns of Key Genes in the Photoperiod and Vernalization Flowering Pathways in Lilium longiflorum with Different Bulb Sizes

Lilium longiflorum is a wild Lilium, and its flowering transition requires a long period of cold exposure to meet the demand of vernalization. The responses of different sized bulbs to cold exposure and photoperiod are different, and the floral transition pathways of small and large bulbs are different. In this study, small and large bulbs were placed in cold storage for different weeks and then cultured at a constant ambient temperature of 25 °C under long day (LD) and short day (SD) conditions. Then, the flowering characteristics and expression patterns of key genes related to the vernalization and photoperiod pathways in different groups were calculated and analyzed. The results showed that the floral transition of Lilium longiflorum was influenced by both vernalization and photoperiod, that vernalization and LD conditions can significantly improve the flowering rate of Lilium longiflorum, and that the time from planting to visible flowering buds’ appearance was decreased. The flowering time and rate of large bulbs were greatly influenced by cold exposure, and the vernalization pathway acted more actively at the floral transition stage. The floral transition of small bulbs was affected more by the photoperiod pathway. Moreover, it was speculated that cold exposure may promote greater sensitivity of the small bulbs to LD conditions. In addition, the expression of LlVRN1, LlFKF1, LlGI, LlCO5, LlCO7, LlCO16, LlFT1, LlFT3 and LlSOC1 was high during the process of floral transition, and LlCO13, LlCO14 and LlCO15 were highly expressed in the vegetative stage. The expression of LlCO13 and LlCO14 was different under different lighting conditions, and the flowering induction function of LlCO9 and LlFT3 was related to vernalization. Moreover, LlFKF1, LlGI, LlCO5, LlCO16, LlSOC1 and LlFT2 were involved in the entire growth process of plants, while LlCO6, LlCO16 and LlFT1 are involved in the differentiation and formation of small bulblets of plants after the inflorescence stage, and this process is also closely related to LD conditions. This study has great significance for understanding the molecular mechanisms of the vernalization and photoperiod flowering pathways of Lilium longiflorum.

Pleiotropic Effect of the compactum Gene and Its Combined Effects with Other Loci for Spike and Grain-Related Traits in Wheat

Club wheat (Triticum aestivum ssp. compactum) with a distinctly compact spike morphology was conditioned by the dominant compactum (C) locus on chromosome 2D and resulted in a redistribution of spike yield components. The disclosure of the genetic basis of club wheat was a prerequisite for the development of widely adapted, agronomically competitive club wheat cultivars. In this study, we used a recombinant inbred line population derived from a cross between club wheat Hiller and modern cultivar Yangmai 158 to construct a genetic linkage map and identify quantitative trait loci associated with 15 morphological traits. The club allele acted in a semi-dominant manner and the C gene was mapped to 370.12–406.29 Mb physical region on the long arm of 2D. Apart from compact spikes, C exhibited a pleiotropic effect on ten other agronomic traits, including plant height, three spike-related traits and six grain-related traits. The compact spike phenotype was correlated with decreased grain size and weight, but with an increase in floret fertility and grain number. These pleiotropic effects make club wheat have compatible spike weight with a normal spike from common wheat. The genetic effects of various gene combinations of C with four yield-related genes, including Ppd-D1, Vrn-D3, Rht-B1b and Rht8, were evaluated. C had no epistatic interaction with any of these genes, indicating that their combinations would have an additive effect on other agronomically important traits. Our research provided a theoretical foundation for the potentially effective deployment of C gene into modern breeding varieties in combination with other favorable alleles.

Chromosome-scale assembly and annotation of the perennial ryegrass genome

BACKGROUND

The availability of chromosome-scale genome assemblies is fundamentally important to advance genetics and breeding in crops, as well as for evolutionary and comparative genomics. The improvement of long-read sequencing technologies and the advent of optical mapping and chromosome conformation capture technologies in the last few years, significantly promoted the development of chromosome-scale genome assemblies of model plants and crop species. In grasses, chromosome-scale genome assemblies recently became available for cultivated and wild species of the Triticeae subfamily. Development of state-of-the-art genomic resources in species of the Poeae subfamily, which includes important crops like fescues and ryegrasses, is lagging behind the progress in the cereal species.

RESULTS

Here, we report a new chromosome-scale genome sequence assembly for perennial ryegrass, obtained by combining PacBio long-read sequencing, Illumina short-read polishing, BioNano optical mapping and Hi-C scaffolding. More than 90% of the total genome size of perennial ryegrass (approximately 2.55 Gb) is covered by seven pseudo-chromosomes that show high levels of collinearity to the orthologous chromosomes of Triticeae species. The transposon fraction of perennial ryegrass was found to be relatively low, approximately 35% of the total genome content, which is less than half of the genome repeat content of cultivated cereal species. We predicted 54,629 high-confidence gene models, 10,287 long non-coding RNAs and a total of 8,393 short non-coding RNAs in the perennial ryegrass genome.

CONCLUSIONS

The new reference genome sequence and annotation presented here are valuable resources for comparative genomic studies in grasses, as well as for breeding applications and will expedite the development of productive varieties in perennial ryegrass and related species.

Dissecting the Genetic Architecture of Phenology Affecting Adaptation of Spring Bread Wheat Genotypes to the Major Wheat-Producing Zones in India

Spring bread wheat adaptation to diverse environments is supported by various traits such as phenology and plant architecture. A large-scale genome-wide association study (GWAS) was designed to investigate and dissect the genetic architecture of phenology affecting adaptation. It used 48 datasets from 4,680 spring wheat lines. For 8 years (2014-2021), these lines were evaluated for days to heading (DH) and maturity (DM) at three sites: Jabalpur, Ludhiana, and Samastipur (Pusa), which represent the three major Indian wheat-producing zones: the Central Zone (CZ), North-Western Plain Zone (NWPZ), and North-Eastern Plain Zone (NEPZ), respectively. Ludhiana had the highest mean DH of 103.8 days and DM of 148.6 days, whereas Jabalpur had the lowest mean DH of 77.7 days and DM of 121.6 days. We identified 119 markers significantly associated with DH and DM on chromosomes 5B (76), 2B (18), 7D (10), 4D (8), 5A (1), 6B (4), 7B (1), and 3D (1). Our results clearly indicated the importance of the photoperiod-associated gene (Ppd-B1) for adaptation to the NWPZ and the Vrn-B1 gene for adaptation to the NEPZ and CZ. A maximum variation of 21.1 and 14% was explained by markers 2B_56134146 and 5B_574145576 linked to the Ppd-B1 and Vrn-B1 genes, respectively, indicating their significant role in regulating DH and DM. The results provide important insights into the genomic regions associated with the two phenological traits that influence adaptation to the major wheat-producing zones in India.

The genetic architecture of flowering time changes in pea from wild to crop

Change in phenology has been an important component in crop evolution, and selection for earlier flowering through a reduction in environmental sensitivity has helped broaden adaptation in many species. Natural variation for flowering in domesticated pea (Pisum sativum L.) has been noted and studied for decades, but there has been no clear account of change relative to its wild progenitor. Here we examined the genetic control of differences in flowering time between wild P. sativum ssp. humile and a typical late-flowering photoperiodic P. s. sativum accession in a recombinant inbred population under long and short photoperiods. Our results confirm the importance of the major photoperiod sensitivity locus Hr/PsELF3a and identify two other loci on chromosomes 1 (DTF1) and 3 (DTF3) that contribute to earlier flowering in the domesticated line under both photoperiods. The domesticated allele at a fourth locus on chromosome 6 (DTF6) delays flowering under long days only. Map positions, inheritance patterns, and expression analyses in near-isogenic comparisons imply that DTF1, DTF3, and DTF6 represent gain-of-function alleles of the florigen/antiflorigen genes FTa3, FTa1, and TFL1c/LF, respectively. This echoes similar variation in chickpea and lentil, and suggests a conserved route to reduced photoperiod sensitivity and early phenology in temperate pulses.

Major niche transitions in Pooideae correlate with variation in photoperiodic flowering and evolution of CCT domain genes

The external cues that trigger timely flowering vary greatly across tropical and temperate plant taxa, the latter relying on predictable seasonal fluctuations in temperature and photoperiod. In the grass family (Poaceae) for example, species of the subfamily Pooideae have become specialists of the northern temperate hemisphere, generating the hypothesis that their progenitor evolved a flowering response to long days from a short-day or day-neutral ancestor. Sampling across the Pooideae, we found support for this hypothesis, and identified several secondary shifts to day-neutral flowering and one to short-day flowering in a tropical highland clade. To explain the proximate mechanisms for the secondary transition back to short-day-regulated flowering, we investigated the expression of CCT domain genes, some of which are known to repress flowering in cereal grasses under specific photoperiods. We found a shift in CONSTANS 1 and CONSTANS 9 expression that coincides with the derived short-day photoperiodism of our exemplar species Nassella pubiflora. This sets up the testable hypothesis that trans- or cis-regulatory elements of these CCT domain genes were the targets of selection for major niche shifts in Pooideae grasses.

Association of the Recessive Allele vrn-D1 With Winter Frost Tolerance in Bread Wheat

Winter frost has been considered the primary limiting factor in wheat production. Shimai 12 is an elite wheat cultivar grown in central and southern Hebei province of China, but sensitive to winter frost. In this study, the winter frost tolerant cultivar Lunxuan 103 was bred by introducing the recessive allele vrn-D1 from winter wheat Shijiazhuang 8 (frost tolerance) into Shimai 12 using marker-assisted selection (MAS). Different from Shimai 12, Lunxuan 103 exhibited a winter growth habit with strong winter frost tolerance. In the Shimai 12 × Shijiazhuang 8 population, the winter progenies (vrn-D1vrn-D1) had significantly lower winter-killed seedling/tiller rates than spring progenies (Vrn-D1aVrn-D1a), and the consistent result was observed in an association population. Winter frost damage caused a significant decrease in grain yield and spike number/m2 in Shimai 12, but not in Lunxuan 103 and Shijiazhuang 8. The time-course expression analysis showed that the transcript accumulation levels of the cold-responsive genes were higher in Lunxuan 103 and Shijiazhuang 8 than in Shimai 12. Lunxuan 103 possessed the same alleles as its parents in the loci for plant height, vernalization, and photoperiod, except for the vernalization gene Vrn-D1. An analysis of genomic composition showed that the two parents contributed similar proportions of genetic compositions to Lunxuan 103. This study provides an example of the improvement of winter frost tolerance by introducing the recessive vernalization gene in bread wheat.

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.

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.

Genome analysis in Avena sativa reveals hidden breeding barriers and opportunities for oat improvement

Oat (Avena sativa L.) is an important and nutritious cereal crop, and there is a growing need to identify genes that contribute to improved oat varieties. Here we utilize a newly sequenced and annotated oat reference genome to locate and characterize quantitative trait loci (QTLs) affecting agronomic and grain-quality traits in five oat populations. We find strong and significant associations between the positions of candidate genes and QTL that affect heading date, as well as those that influence the concentrations of oil and β-glucan in the grain. We examine genome-wide recombination profiles to confirm the presence of a large, unbalanced translocation from chromosome 1 C to 1 A, and a possible inversion on chromosome 7D. Such chromosome rearrangements appear to be common in oat, where they cause pseudo-linkage and recombination suppression, affecting the segregation, localization, and deployment of QTLs in breeding programs.

Tinker et al. identified the position and effects of major QTLs relative to a new fully annotated reference genome in five recombinant inbred line populations representing nine diverse oat (Avena sativa) varieties. They also characterized two major chromosome rearrangements that may impact breeding targets affected by QTL that are located in these regions.

A single nucleotide deletion in the third exon of FT-D1 increases the spikelet number and delays heading date in wheat (Triticum aestivum L.)

The spikelet number and heading date are two crucial and correlated traits for yield in wheat. Here, a quantitative trait locus (QTL) analysis was conducted in F₈ recombinant inbred lines (RILs) derived from crossing two common wheats with different spikelet numbers. A total of 15 stable QTL influencing total spikelet number (TSN) and heading date (HD) were detected. Notably, FT-D1, a well-known flowering time gene in wheat, was located within the finely mapped interval of a major QTL on 7DS (QTsn/Hd.cau-7D). A causal indel of one G in the third exon of FT-D1 was significantly associated with total spikelet number and heading date. Consistently, CRISPR/Cas9 mutant lines with homozygous mutations in FT-D1 displayed an increase in total spikelet number and heading date when compared with wild type. Moreover, one simple and robust marker developed according to the polymorphic site of FT-D1 revealed that this one G indel had been preferentially selected to adapt to different environments. Collectively, these data provide further insights into the genetic basis of spikelet number and heading date, and the diagnostic marker of FT-D1 will be useful for marker-assisted pyramiding in wheat breeding.

Comparative analysis of buds transcriptome and identification of two florigen gene AkFTs in Amorphophallus konjac

Leaves and flowers of Amorphophallus konjac do not develop simultaneously thus unique features can be elucidated through study of flowering transformation in A. konjac. In this study, transcriptome libraries of A. konjac leaf buds (LB) and flower buds (FB) were constructed followed by high-throughput sequencing. A total of 68,906 unigenes with an average length of 920 bp were obtained after library assembly. Out of these genes, 24,622 unigenes had annotation information. A total of 6859 differentially expressed genes (DEGs) were identified through differential expression analysis using LB as control. Notably, 2415 DEGs were upregulated whereas 4444 DEGs were downregulated in the two transcriptomes. Go and KEGG analysis showed that the DEGs belonged to 44 functional categories and were implicated in 98 metabolic pathways and 38 DEGs involved in plant hormone signal transduction. Several genes were mined that may be involved in A. konjac flower bud differentiation and flower organ development. Eight DEGs were selected for verification of RNA-seq results using qRT-PCR analysis. Two FLOWERING LOCUS T (FT) genes named AkFT1 and AkFT2 were identified though homologous analysis may be the florigen gene implicated in modulation of A. konjac flowering. These genes were significantly upregulated in flower buds compared with the expression levels on leaf buds. Overexpression of AkFT genes though heterologous expression in Arabidopsis showed that the transgenics flowered at a very early stage relative to wild type plants. These findings indicate that AkFT1 and AkFT2 function as regulation genes in A. konjac flowering development and the two genes may present similar functions during flowering transition.

Identification and validation of a QTL for spikelet number on chromosome arm 6BL of common wheat (Triticum aestivum L.)

To provide food security for a growing world population, it will be necessary to increase yields of staple crops such as wheat (Triticum aestivum L.). Yield is a complex, polygenic trait influenced by grain weight and number, which are negatively correlated with one another. Spikelet number is an important determinant of grain number, but allelic variants impacting its expression are often associated with heading date, constraining their use in wheat germplasm that must be adapted for specific environments. Identification and characterization of genetic variants affecting spikelet number will increase selection efficiency through their deployment in breeding programs. In this study, a quantitative trait locus (QTL) on chromosome arm 6BL for spikelet number was identified and validated using an association mapping panel, a recombinant inbred line population, and seven derived heterogeneous inbred families. The superior allele, QSn.csu-6Bb, was associated with an increase of 0.248 to 0.808 spikelets per spike across multiple environments that varied for mean spikelet number. Despite epistatic interactions between QSn.csu-6B and three other loci (WAPO-A1, VRN-D3, and PPD-B1), genotypes with a greater number of superior alleles at these loci consistently exhibit higher spikelet number. The frequency of superior alleles at these loci varies among winter wheat varieties adapted to different latitudes of the US Great Plains, revealing opportunities for breeders to select for increased spikelet number using simple molecular markers. This work lays the foundation for the positional cloning of the genetic variant underlying the QSn.csu-6B QTL to strengthen our understanding of spikelet number determination in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01288-7.

MiR172-APETALA2-like genes integrate vernalization and plant age to control flowering time in wheat

Plants possess regulatory mechanisms that allow them to flower under conditions that maximize reproductive success. Selection of natural variants affecting those mechanisms has been critical in agriculture to modulate the flowering response of crops to specific environments and to increase yield. In the temperate cereals, wheat and barley, the photoperiod and vernalization pathways explain most of the natural variation in flowering time. However, other pathways also participate in fine-tuning the flowering response. In this work, we integrate the conserved microRNA miR172 and its targets APETALA2-like (AP2L) genes into the temperate grass flowering network involving VERNALIZATION 1 (VRN1), VRN2 and FLOWERING LOCUS T 1 (FT1 = VRN3) genes. Using mutants, transgenics and different growing conditions, we show that miR172 promotes flowering in wheat, while its target genes AP2L1 (TaTOE1) and AP2L5 (Q) act as flowering repressors. Moreover, we reveal that the miR172-AP2L pathway regulates FT1 expression in the leaves, and that this regulation is independent of VRN2 and VRN1. In addition, we show that the miR172-AP2L module and flowering are both controlled by plant age through miR156 in spring cultivars. However, in winter cultivars, flowering and the regulation of AP2L1 expression are decoupled from miR156 downregulation with age, and induction of VRN1 by vernalization is required to repress AP2L1 in the leaves and promote flowering. Interestingly, the levels of miR172 and both AP2L genes modulate the flowering response to different vernalization treatments in winter cultivars. In summary, our results show that conserved and grass specific gene networks interact to modulate the flowering response, and that natural or induced mutations in AP2L genes are useful tools for fine-tuning wheat flowering time in a changing environment.

The wild grass Brachypodium distachyon as a developmental model system

The arrival of cheap and high-throughput sequencing paired with efficient gene editing technologies allows us to use non-traditional model systems and mechanistically approach biological phenomena beyond what was conceivable just a decade ago. Venturing into different model systems enables us to explore for example clade-specific environmental responses to changing climates or the genetics and development of clade-specific organs, tissues and cell types. We-both early career researchers working with the wild grass model Brachypodium distachyon-want to use this review to (1) highlight why we think B. distachyon is a fantastic grass developmental model system, (2) summarize the tools and resources that have enabled discoveries made in B. distachyon, and (3) discuss a handful of developmental biology vignettes made possible by using B. distachyon as a model system. Finally, we want to conclude by (4) relating our personal stories with this emerging model system and (5) share what we think is important to consider before starting work with an emerging model system.

Hybrids Provide More Options for Fine-Tuning Flowering Time Responses of Winter Barley

Crop adaptation requires matching resource availability to plant development. Tight coordination of the plant cycle with prevailing environmental conditions is crucial to maximizing yield. It is expected that winters in temperate areas will become warmer, so the vernalization requirements of current cultivars can be desynchronized with the environment's vernalizing potential. Therefore, current phenological ideotypes may not be optimum for future climatic conditions. Major genes conferring vernalization sensitivity and phenological responses in barley (Hordeum vulgare L.) are known, but some allelic combinations remain insufficiently evaluated. Furthermore, there is a lack of knowledge about flowering time in a hybrid context. To honor the promise of increased yield potentials, hybrid barley phenology must be studied, and the knowledge deployed in new cultivars. A set of three male and two female barley lines, as well as their six F₁ hybrids, were studied in growth chambers, subjected to three vernalization treatments: complete (8 weeks), moderate (4 weeks), and low (2 weeks). Development was recorded up to flowering, and expression of major genes was assayed at key stages. We observed a gradation in responses to vernalization, mostly additive, concentrated in the phase until the initiation of stem elongation, and proportional to the allele constitution and dosage present in VRN-H1. These responses were further modulated by the presence of PPD-H2. The duration of the late reproductive phase presented more dominance toward earliness and was affected by the rich variety of alleles at VRN-H3. Our results provide further opportunities for fine-tuning total and phasal growth duration in hybrid barley, beyond what is currently feasible in inbred cultivars.

MIG-seq is an effective method for high-throughput genotyping in wheat (Triticum spp.)

MIG-seq (Multiplexed inter-simple sequence repeats genotyping by sequencing) has been developed as a low cost genotyping technology, although the number of polymorphisms obtained is assumed to be minimal, resulting in the low application of this technique to analyses of agricultural plants. We applied MIG-seq to 12 plant species that include various crops and investigated the relationship between genome size and the number of bases that can be stably sequenced. The genome size and the number of loci, which can be sequenced by MIG-seq, are positively correlated. This is due to the linkage between genome size and the number of simple sequence repeats (SSRs) through the genome. The applicability of MIG-seq to population structure analysis, linkage mapping, and quantitative trait loci (QTL) analysis in wheat, which has a relatively large genome, was further evaluated. The results of population structure analysis for tetraploid wheat showed the differences among collection sites and subspecies, which agreed with previous findings. Additionally, in wheat biparental mapping populations, over 3,000 SNPs/indels with low deficiency were detected using MIG-seq, and the QTL analysis was able to detect recognized flowering-related genes. These results revealed the effectiveness of MIG-seq for genomic analysis of agricultural plants with large genomes, including wheat.

Photoperiod Control of Plant Growth: Flowering Time Genes Beyond Flowering

Fluctuations in environmental conditions greatly influence life on earth. Plants, as sessile organisms, have developed molecular mechanisms to adapt their development to changes in daylength, or photoperiod. One of the first plant features that comes to mind as affected by the duration of the day is flowering time; we all bring up a clear image of spring blossom. However, for many plants flowering happens at other times of the year, and many other developmental aspects are also affected by changes in daylength, which range from hypocotyl elongation in Arabidopsis thaliana to tuberization in potato or autumn growth cessation in trees. Strikingly, many of the processes affected by photoperiod employ similar gene networks to respond to changes in the length of light/dark cycles. In this review, we have focused on developmental processes affected by photoperiod that share similar genes and gene regulatory networks.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

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.

Responses of Barley to High Ambient Temperature Are Modulated by Vernalization

Ambient temperatures are increasing due to climate change. Cereal crops development and production will be affected consequently. Flowering time is a key factor for adaptation of small grain cereals and, therefore, exploring developmental responses of barley to rising temperatures is required. In this work, we studied phasic growth, and inflorescence traits related to yield, in eight near isogenic lines of barley (Hordeum vulgare L.) differing at the VRN-H1, VRN-H2 and PPD-H1 genes, representing different growth habits. The lines were grown in contrasting vernalization treatments, under two temperature regimes (18 and 25°C), in long days. Lines with recessive ppd-H1 presented delayed development compared to lines with the sensitive PPD-H1 allele, across the two growth phases considered. High temperature delayed flowering in all unvernalized plants, and in vernalized spring barleys carrying the insensitive ppd-H1 allele, whilst it accelerated flowering in spring barleys with the sensitive PPD-H1 allele. This finding evidenced an interaction between PPD-H1, temperature and vernalization. At the high temperature, PPD-H1 lines in spring backgrounds (VRN-H1-7) yielded more, whereas lines with ppd-H1 were best in vrn-H1 background. Our study revealed new information that will support breeding high-yielding cultivars with specific combinations of major adaptation genes tailored to future climatic conditions.

Population genomic analysis of Aegilops tauschii identifies targets for bread wheat improvement

Aegilops tauschii, the diploid wild progenitor of the D subgenome of bread wheat, is a reservoir of genetic diversity for improving bread wheat performance and environmental resilience. Here we sequenced 242 Ae. tauschii accessions and compared them to the wheat D subgenome to characterize genomic diversity. We found that a rare lineage of Ae. tauschii geographically restricted to present-day Georgia contributed to the wheat D subgenome in the independent hybridizations that gave rise to modern bread wheat. Through k-mer-based association mapping, we identified discrete genomic regions with candidate genes for disease and pest resistance and demonstrated their functional transfer into wheat by transgenesis and wide crossing, including the generation of a library of hexaploids incorporating diverse Ae. tauschii genomes. Exploiting the genomic diversity of the Ae. tauschii ancestral diploid genome permits rapid trait discovery and functional genetic validation in a hexaploid background amenable to breeding.

Major Genomic Regions for Wheat Grain Weight as Revealed by QTL Linkage Mapping and Meta-Analysis

Grain weight is a key determinant for grain yield potential in wheat, which is highly governed by a type of quantitative genetic basis. The identification of major quantitative trait locus (QTL) and functional genes are urgently required for molecular improvements in wheat grain yield. In this study, major genomic regions and putative candidate genes for thousand grain weight (TGW) were revealed by integrative approaches with QTL linkage mapping, meta-analysis and transcriptome evaluation. Forty-five TGW QTLs were detected using a set of recombinant inbred lines, explaining 1.76-12.87% of the phenotypic variation. Of these, ten stable QTLs were identified across more than four environments. Meta-QTL (MQTL) analysis were performed on 394 initial TGW QTLs available from previous studies and the present study, where 274 loci were finally refined into 67 MQTLs. The average confidence interval of these MQTLs was 3.73-fold less than that of initial QTLs. A total of 134 putative candidate genes were mined within MQTL regions by combined analysis of transcriptomic and omics data. Some key putative candidate genes similar to those reported early for grain development and grain weight formation were further discussed. This finding will provide a better understanding of the genetic determinants of TGW and will be useful for marker-assisted selection of high yield in wheat breeding.

Allelic Variation Analysis at the Vernalization Response and Photoperiod Genes in Russian Wheat Varieties Identified Two Novel Alleles of Vrn-B3

Heading time is an important agronomic trait affecting the adaptability and productivity of common wheat. In this study, 95 common wheat varieties from Russia and the late-maturing breeding line ‘Velut’ were tested for allelic diversity of genes having the strongest effect on heading. In this research, allelic variation at the Ppd-D1, Vrn-A1, Vrn-B1, Vrn-D1, and Vrn-B3 loci was tested. The Vrn-B1 and Vrn-B3 loci provided the largest contribution to genetic diversity. We found two novel allelic variants of the Vrn-B3 gene in the studied varieties. Ten varieties carried a 160 bp insertion in the promoter region, and the breeding line ‘Velut’ carried a 1617 bp insertion. These alleles were designated Vrn-B3e and Vrn-B3d, respectively. The analysis of the sequences showed the recent insertion of a retrotransposon homologous to the LTR retrotransposon (RLX_Hvul_Dacia_ RND-1) in the Vrn-B3d allele. Plants with the Vrn-B3e and the ‘Velut’ line with the Vrn-B3d allele headed later than the plants with the wild-type allele; among these plants, ‘Velut’ is the latest maturing wheat variety. Analysis of the gene expression of two groups of lines differing by the Vrn-B3 alleles (Vrn-B3d or vrn-B3) from the F₂ population with ‘Velut’ as a parental line did not reveal a significant difference in the expression level between the groups. Additional research is required to study the reasons for the late maturation of the ‘Velut’ line. However, the studied wheat varieties could be used as a potential source of natural variation in genes controlling heading times.

Juvenile Heat Tolerance in Wheat for Attaining Higher Grain Yield by Shifting to Early Sowing in October in South Asia

Farmers in northwestern and central India have been exploring to sow their wheat much earlier (October) than normal (November) to sustain productivity by escaping terminal heat stress and to utilize the available soil moisture after the harvesting of rice crop. However, current popular varieties are poorly adapted to early sowing due to the exposure of juvenile plants to the warmer temperatures in the month of October and early November. Therefore, a study was undertaken to identify wheat genotypes suited to October sowing under warmer temperatures in India. A diverse collection of 3322 bread wheat varieties and elite lines was prepared in CIMMYT, Mexico, and planted in the 3rd week of October during the crop season 2012-2013 in six locations (Ludhiana, Karnal, New Delhi, Indore, Pune and Dharwad) spread over northwestern plains zone (NWPZ) and central and Peninsular zone (CZ and PZ; designated as CPZ) of India. Agronomic traits data from the seedling stage to maturity were recorded. Results indicated substantial diversity for yield and yield-associated traits, with some lines showing indications of higher yields under October sowing. Based on agronomic performance and disease resistance, the top 48 lines (and two local checks) were identified and planted in the next crop season (2013-2014) in a replicated trial in all six locations under October sowing (third week). High yielding lines that could tolerate higher temperature in October sowing were identified for both zones; however, performance for grain yield was more promising in the NWPZ. Hence, a new trial of 30 lines was planted only in NWPZ under October sowing. Lines showing significantly superior yield over the best check and the most popular cultivars in the zone were identified. The study suggested that agronomically superior wheat varieties with early heat tolerance can be obtained that can provide yield up to 8 t/ha by planting in the third to fourth week of October.

The interphase period "germination-heading" of 8x and 6x triticale with different dominant Vrn genes

The existing spring forms of wheat-rye amphiploids are characterized by late maturity due to the long duration of the interphase period “germination–heading”. The manifestation of this trait is inf luenced by Vrn-1 genes. Their dominant alleles also determine the spring type of development. The results of studying the interphase period “germination– heading” of spring octaploid and hexaploid forms of triticale created for use in research and breeding programs under the conditions of forest-steppe of Western Siberia are given in this article. The interphase period of the primary forms 8xVrnA1, 8xVrnB1 and 8xVrnD1 obtained by artif icial doubling of the chromosome number of the wheat-rye hybrids made by pollination of three lines of the soft wheat ‘Triple Dirk’ – donors of different dominant Vrn-1 genes – by a winter rye variety ‘Korotkostebel’naya 69’ was determined under the f ield conditions in the nursery of octaploid (8x) triticale. In the nursery of hexaploid triticale, this trait was studied in the populations of hybrids obtained by hybridization of these three primary forms of octaploid triticale with the hexaploid winter triticale variety ‘Sears 57’. In the offspring of crossing 8хVrnD1 × ‘Sears 57’, spring genotypes of 6x triticale bearing Vrn-D1 were selected. This fact was determined by PСR. It means that the genetic material from the chromosome of the f ifth homeologous group of the D genome of the bread wheat is included in the plant genotypes. This genome is absent in the winter 6x triticale ‘Sears 57’. The grain content of spikes of the created hexaploid forms of triticale is superiour to that of the maternal octaploid triticale forms. It was shown that plants of the hybrid populations 8xVrnA1 × ‘Sears 57’ and 8xVrnD1 × ‘Sears 57’ carrying the dominant alleles Vrn-A1a and Vrn-D1a, respectively, have a shorter duration of the “germination–heading” interphase period than the initial parental forms of primary 8x triticale. The short interphase period of “germination–heading” of the 6x triticale is a valuable breading trait for the creation of early maturing and productive genotypes of triticale.

The Relationships between Plant Developmental Traits and Winter Field Survival in Rye (Secale cereale L.)

Overwintering cereals accumulate low temperature tolerance (LTT) during cold acclimation in the autumn. Simultaneously, the plants adjust to the colder season by making developmental changes at the shoot apical meristem. These processes lead to higher winter hardiness in winter rye varieties (Secale cereale L.) adapted to Northern latitudes as compared to other cereal crops. To dissect the winter-hardiness trait in rye, a panel of 96 genotypes of different origins and growth habits was assessed for winter field survival (WFS), LTT, and six developmental traits. Best Linear Unbiased Estimates for WFS determined from five field trials correlated strongly with LTT (r = 0.90, p < 0.001); thus, cold acclimation efficiency was the major contributor to WFS. WFS also correlated strongly (p < 0.001) with final leaf number (r = 0.80), prostrate growth habit (r = 0.61), plant height (r = 0.34), but showed weaker associations with top internode length (r = 0.30, p < 0.01) and days to anthesis (r = 0.25, p < 0.05). The heritability estimates (h²) for WFS-associated traits ranged from 0.45 (prostrate growth habit) to 0.81 (final leaf number) and were overall higher than for WFS (h² = 0.48). All developmental traits associated with WFS and LTT are postulated to be regulated by phytohormone levels at shoot apical meristem.

In-Depth Sequence Analysis of Bread Wheat VRN1 Genes

The VERNALIZATION1 (VRN1) gene encodes a MADS-box transcription factor and plays an important role in the cold-induced transition from the vegetative to reproductive stage. Allelic variability of VRN1 homoeologs has been associated with large differences in flowering time. The aim of this study was to investigate the genetic variability of VRN1 homoeologs (VRN-A1, VRN-B1 and VRN-D1). We performed an in-depth sequence analysis of VRN1 homoeologs in a panel of 105 winter and spring varieties of hexaploid wheat. We describe the novel allele Vrn-B1f with an 836 bp insertion within intron 1 and show its specific expression pattern associated with reduced heading time. We further provide the complete sequence of the Vrn-A1b allele, revealing a 177 bp insertion in intron 1, which is transcribed into an alternative splice variant. Copy number variation (CNV) analysis of VRN1 homoeologs showed that VRN-B1 and VRN-D1 are present in only one copy. The copy number of recessive vrn-A1 ranged from one to four, while that of dominant Vrn-A1 was one or two. Different numbers of Vrn-A1a copies in the spring cultivars Branisovicka IX/49 and Bastion did not significantly affect heading time. We also report on the deletion of secondary structures (G-quadruplex) in promoter sequences of cultivars with more vrn-A1 copies.

Temporal expression study of miRNAs in the crown tissues of winter wheat grown under natural growth conditions

BACKGROUND

Winter wheat requires prolonged exposure to low temperature to initiate flowering (vernalization). Shoot apical meristem of the crown is the site of cold perception, which produces leaf primordia during vegetative growth before developing into floral primordia at the initiation of the reproductive phase. Although many essential genes for winter wheat cold acclimation and floral initiation have been revealed, the importance of microRNA (miRNA) meditated post-transcriptional regulation in crowns is not well understood. To understand the potential roles of miRNAs in crown tissues, we performed a temporal expression study of miRNAs in crown tissues at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage of winter wheat grown under natural growth conditions.

RESULTS

In total, 348 miRNAs belonging to 298 miRNA families, were identified in wheat crown tissues. Among them, 92 differentially expressed miRNAs (DEMs) were found to be significantly regulated from the three-leaf stage to the jointing stage. Most of these DEMs were highly expressed at the three-leaf stage and winter dormancy stage, and then declined in later stages. Six DEMs, including miR156a-5p were markedly induced during the winter dormancy stage. Eleven DEMs, including miR159a.1, miR390a-5p, miR393-5p, miR160a-5p, and miR1436, were highly expressed at the green-up stage. Twelve DEMs, such as miR172a-5p, miR394a, miR319b-3p, and miR9676-5p were highly induced at the jointing stage. Moreover, 14 novel target genes of nine wheat or Pooideae-specific miRNAs were verified using RLM-5' RACE assay. Notably, six mTERFs and two Rf1 genes, which are associated with mitochondrial gene expression, were confirmed as targets of three wheat-specific miRNAs.

CONCLUSIONS

The present study not only confirmed the known miRNAs associated with phase transition and floral development, but also identified a number of wheat or Pooideae-specific miRNAs critical for winter wheat cold acclimation and floral development. Most importantly, this study provided experimental evidence that miRNA could regulate mitochondrial gene expression by targeting mTERF and Rf1 genes. Our study provides valuable information for further exploration of the mechanism of miRNA mediated post-transcriptional regulation during winter wheat vernalization and inflorescent initiation.

Yield-Related QTL Clusters and the Potential Candidate Genes in Two Wheat DH Populations

In the present study, four large-scale field trials using two doubled haploid wheat populations were conducted in different environments for two years. Grain protein content (GPC) and 21 other yield-related traits were investigated. A total of 227 QTL were mapped on 18 chromosomes, which formed 35 QTL clusters. The potential candidate genes underlying the QTL clusters were suggested. Furthermore, adding to the significant correlations between yield and its related traits, correlation variations were clearly shown within the QTL clusters. The QTL clusters with consistently positive correlations were suggested to be directly utilized in wheat breeding, including 1B.2, 2A.2, 2B (4.9–16.5 Mb), 2B.3, 3B (68.9–214.5 Mb), 4A.2, 4B.2, 4D, 5A.1, 5A.2, 5B.1, and 5D. The QTL clusters with negative alignments between traits may also have potential value for yield or GPC improvement in specific environments, including 1A.1, 2B.1, 1B.3, 5A.3, 5B.2 (612.1–613.6 Mb), 7A.1, 7A.2, 7B.1, and 7B.2. One GPC QTL (5B.2: 671.3–672.9 Mb) contributed by cultivar Spitfire was positively associated with nitrogen use efficiency or grain protein yield and is highly recommended for breeding use. Another GPC QTL without negatively pleiotropic effects on 2A (50.0–56.3 Mb), 2D, 4D, and 6B is suggested for quality wheat breeding.

Identification and validation of two major QTLs for spike compactness and length in bread wheat (Triticum aestivum L.) showing pleiotropic effects on yield-related traits

Two major and stable QTLs for spike compactness and length were detected and validated in multiple genetic backgrounds and environments, and their pleiotropic effects on yield-related traits were analyzed. Spike compactness (SC) and length (SL) are greatly associated with wheat (Triticum aestivum L.) grain yield. To detect quantitative trait loci (QTL) associated with SC and SL, two biparental populations derived from crosses of Chuanmai42/Kechengmai1 and Chuanmai42/Chuannong16 were employed to perform QTL mapping in five environments. A total of 34 QTLs were identified, in which six major QTLs were repeatedly detected in more than four environments and the best linear unbiased prediction datasets, explaining 7.13-33.6% of phenotypic variation. These major QTLs were co-located in two genomic regions on chromosome 5A and 6A, namely QSc/Sl.cib-5A and QSc/Sl.cib-6A, respectively. By developing kompetitive allele-specific PCR (KASP) markers that linked to them, the two loci were validated in different genetic backgrounds, and their interactions were also analyzed. Comparison analysis showed that QSc/Sl.cib-5A was not Vrn-A1 and Q, and QSc/Sl.cib-6A was likely a new locus for SC and SL. Both QSc/Sl.cib-5A and QSc/Sl.cib-6A had pleiotropic effects on other yield-related traits including plant height, thousand grain weight and grain length. Therefore, the two loci combined with the developed KASP markers might be potentially applicable in wheat breeding. Furthermore, based on the spatiotemporal expression patterns, gene annotation, orthologous search and sequence differences, TraesCS5A01G301400 and TraesCS6A01G090300 were considered as potential candidates for QSc/Sl.cib-5A and QSc/Sl.cib-6A, respectively. These results provided valuable information for fine mapping and cloning of the two loci in the future.

Understanding the Effects of Growing Seasons, Genotypes, and Their Interactions on the Anthesis Date of Wheat Sown in North China

Simple Summary

Wheat anthesis date is an important turning point for wheat from vegetative growth stage to reproductive growth stage, which is crucial for wheat to adapt to environment and increase grain yield. In this study, a panel of adaptive wheat varieties including historical varieties from the 1940s and current varieties was used to understand the contribution of growing season, genotypes and their interaction effects to anthesis date. Based on our results, we can conclude that growing seasons contributed tremendously to the anthesis date of wheat. In future wheat breeding, more consideration should be given to growing seasons, and the gene combination with the strongest adaptability to the growing seasons should be selected.

Abstract

Quantitative studies on the effects of growing season, genotype (including photoperiod genes and vernalization genes), and their interaction (GGI) on the anthesis date of winter wheat (Triticum aestivum L.) are helpful to provide a scientific reference for selecting or developing adaptive varieties in target environments. In this study, we collected 100 winter wheat varieties with ecological adaptability in North China and identified the anthesis date under field conditions for three consecutive years from 2016 to 2019 with mapped photoperiod and vernalization alleles. Our results showed that the number of the photoperiod-insensitive Ppd-D1a allele increased with variety replacement, while the haplotype Ppd-A1b + Ppd-D1b + vrn-D1 (A4B2) decreased from the 1940s to 2000s. The anthesis date of A4B2 was significantly delayed due to the photoperiod-insensitive alleles Ppd-A1b and Ppd-D1b. The additive main effect and multiplicative interaction (AMMI) model and GGI biplot analysis were used for data analysis. A large portion of the total variation was explained by growing seasons (66.3%), while genotypes and GGIs explained 21.9% and 10.1% of the anthesis dates, respectively. The varieties from the 1940s and before had a great influence on the anthesis date, suggesting these germplasms tend to avoid premature anthesis and could facilitate the development of phenological resilient varieties.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Large-scale integration of meta-QTL and genome-wide association study discovers the genomic regions and candidate genes for yield and yield-related traits in bread wheat

Based on the large-scale integration of meta-QTL and Genome-Wide  Association Study, 76 high-confidence MQTL regions and 237 candidate genes that affected wheat yield and yield-related traits were discovered. Improving yield and yield-related traits are key goals in wheat breeding program. The integration of accumulated wheat genetic resources provides an opportunity to uncover important genomic regions and candidate genes that affect wheat yield. Here, a comprehensive meta-QTL analysis was conducted on 2230 QTL of yield-related traits obtained from 119 QTL studies. These QTL were refined into 145 meta-QTL (MQTL), and 89 MQTL were verified by GWAS with different natural populations. The average confidence interval (CI) of these MQTL was 2.92 times less than that of the initial QTL. Furthermore, 76 core MQTL regions with a physical distance less than 25 Mb were detected. Based on the homology analysis and expression patterns, 237 candidate genes in the MQTL involved in photoperiod response, grain development, multiple plant growth regulator pathways, carbon and nitrogen metabolism and spike and flower organ development were determined. A novel candidate gene TaKAO-4A was confirmed to be significantly associated with grain size, and a CAPS marker was developed based on its dominant haplotype. In summary, this study clarified a method based on the integration of meta-QTL, GWAS and homology comparison to reveal the genomic regions and candidate genes that affect important yield-related traits in wheat. This work will help to lay a foundation for the identification, transfer and aggregation of these important QTL or candidate genes in wheat high-yield breeding.

Molecular responses to chilling in a warming climate and their impacts on plant reproductive development and yield

Responses to prolonged winter chilling are universal in temperate plants which use seasonal temperature cues in the seed, vegetative and reproductive phases to align development with the earth's orbit. Climate change is driving a decline in reliable winter chill and affecting the sub-tropical extent of cultivation for temperate over-wintering crops. Here we explore molecular aspects of plant responses to winter chill including seasonal bud break and flowering, and how variation in the intensity of winter chilling or de-vernalisation can lead to effects on post-chilling plant development, including that of structures necessary for crop yields.

Single-Molecule Real-Time and Illumina-Based RNA Sequencing Data Identified Vernalization-Responsive Candidate Genes in Faba Bean (Vicia faba L.)

Faba bean (Vicia faba L.) is one of the most widely grown cool season legume crops in the world. Winter faba bean normally has a vernalization requirement, which promotes an earlier flowering and pod setting than unvernalized plants. However, the molecular mechanisms of vernalization in faba bean are largely unknown. Discovering vernalization-related candidate genes is of great importance for faba bean breeding. In this study, the whole transcriptome of faba bean buds was profiled by using next-generation sequencing (NGS) and single-molecule, real-time (SMRT) full-length transcriptome sequencing technology. A total of 29,203 high-quality non-redundant transcripts, 21,098 complete coding sequences (CDS), 1,045 long non-coding RNAs (lncRNAs), and 12,939 simple sequence repeats (SSRs) were identified. Furthermore, 4,044 differentially expressed genes (DEGs) were identified through pairwise comparisons. By Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, these differentially expressed transcripts were found to be enriched in binding and transcription factor activity, electron carrier activity, rhythmic process, and receptor activity. Finally, 50 putative vernalization-related genes that played important roles in the vernalization of faba bean were identified; we also found that the levels of vernalization-responsive transcripts showed significantly higher expression levels in cold-treated buds. The expression of VfSOC1, one of the candidate genes, was sensitive to vernalization. Ectopic expression of VfSOC1 in Arabidopsis brought earlier flowering. In conclusion, the abundant vernalization-related transcripts identified in this study will provide a basis for future researches on the vernalization and faba bean breeding and established a reference full-length transcriptome for future studies on faba bean.

TaVrt2, an SVP-like gene, cooperates with TaVrn1 to regulate vernalization-induced flowering in wheat

TaVrn1, encoding a MADS-box transcription factor (TF), is the central regulator of wheat vernalization-induced flowering. Considering that the MADS-box TF usually works by forming hetero- or homodimers, we conducted yeast-two-hybrid screening and identified an SVP-like MADS-box protein TaVrt2 interacting with TaVrn1. However, the specific function of TaVrt2 and the biological implication of its interaction with TaVrn1 remained unknown. We validated the function of TaVrt2 and TaVrn1 by wheat transgenic experiments and their interaction through multiple protein-binding assays. Population genetic analysis also was used to display their interplay. Transcriptomic sequencing and chromatin immunoprecipitation assays were performed to identify their common targets. TaVrt2 and TaVrn1 are flowering promoters in the vernalization pathway and interact physically in vitro, in planta and in wheat cells. Additionally, TaVrt2 and TaVrn1 were significantly induced in leaves by vernalization, suggesting their spatio-temporal interaction during vernalization. Genetic analysis indicated that TaVrt2 and TaVrn1 had significant epistatic effects on flowering time. Furthermore, native TaVrn1 was up-regulated significantly in TaVrn1-OE (overexpression) and TaVrt2-OE lines. Moreover, TaVrt2 could bind with TaVrn1 promoter directly. A TaVrt2-mediated positive feedback loop of TaVrn1 during vernalization was proposed, providing additional understanding on the regulatory mechanism underlying vernalization-induced flowering.

Major flowering time genes of barley: allelic diversity, effects, and comparison with wheat

This review summarizes the allelic series, effects, interactions between genes and with the environment, for the major flowering time genes that drive phenological adaptation of barley. The optimization of phenology is a major goal of plant breeding addressing the production of high-yielding varieties adapted to changing climatic conditions. Flowering time in cereals is regulated by genetic networks that respond predominately to day length and temperature. Allelic diversity at these genes is at the basis of barley wide adaptation. Detailed knowledge of their effects, and genetic and environmental interactions will facilitate plant breeders manipulating flowering time in cereal germplasm enhancement, by exploiting appropriate gene combinations. This review describes a catalogue of alleles found in QTL studies by barley geneticists, corresponding to the genetic diversity at major flowering time genes, the main drivers of barley phenological adaptation: VRN-H1 (HvBM5A), VRN-H2 (HvZCCTa-c), VRN-H3 (HvFT1), PPD-H1 (HvPRR37), PPD-H2 (HvFT3), and eam6/eps2 (HvCEN). For each gene, allelic series, size and direction of QTL effects, interactions between genes and with the environment are presented. Pleiotropic effects on agronomically important traits such as grain yield are also discussed. The review includes brief comments on additional genes with large effects on phenology that became relevant in modern barley breeding. The parallelisms between flowering time allelic variation between the two most cultivated Triticeae species (barley and wheat) are also outlined. This work is mostly based on previously published data, although we added some new data and hypothesis supported by a number of studies. This review shows the wide variety of allelic effects that provide enormous plasticity in barley flowering behavior, which opens new avenues to breeders for fine-tuning phenology of the barley crop.

QTL mapping for grain yield and three yield components in a population derived from two high-yielding spring wheat cultivars

Four genomic regions on chromosomes 4A, 6A, 7B, and 7D were discovered, each with multiple tightly linked QTL (QTL clusters) associated with two to three yield components. The 7D QTL cluster was associated with grain yield, fertile spikelet number per spike, thousand kernel weight, and heading date. It was located in the flanking region of FT-D1, a homolog gene of Arabidopsis FLOWERING LOCUS T, a major gene that regulates wheat flowering. Genetic manipulation of yield components is an important approach to increase grain yield in wheat (Triticum aestivum). The present study used a mapping population comprised of 181 doubled haploid lines derived from two high-yielding spring wheat cultivars, UI Platinum and LCS Star. The two cultivars and the derived population were assessed for six traits in eight field trials primarily in Idaho in the USA. The six traits were grain yield, fertile spikelet number per spike, productive tiller number per unit area, thousand kernel weight, heading date, and plant height. Quantitative Trait Locus (QTL) analysis of the six traits was conducted using 14,236 single-nucleotide polymorphism (SNP) markers generated from the wheat 90 K SNP and the exome and promoter capture arrays. Of the 19 QTL detected, 14 were clustered in four chromosomal regions on 4A, 6A, 7B and 7D. Each of the four QTL clusters was associated with multiple yield component traits, and these traits were often negatively correlated with one another. As a result, additional QTL dissection studies are needed to optimize trade-offs among yield component traits for specific production environments. Kompetitive allele-specific PCR markers for the four QTL clusters were developed and assessed in an elite spring wheat panel of 170 lines, and eight of the 14 QTL were validated. The two parents contain complementary alleles for the four QTL clusters, suggesting the possibility of improving grain yield via genetic recombination of yield component loci.

Reproductive developmental transcriptome analysis of Tripidium ravennae (Poaceae)

BACKGROUND

Tripidium ravennae is a cold-hardy, diploid species in the sugarcane complex (Poaceae subtribe Saccharinae) with considerable potential as a genetic resource for developing improved bioenergy and ornamental grasses. An improved understanding of the genetic regulation of reproductive processes (e.g., floral induction, inflorescence development, and seed development) will enable future applications of precision breeding and gene editing of floral and seed development. In particular, the ability to silence reproductive processes would allow for developing seedless forms of valuable but potentially invasive plants. The objective of this research was to characterize the gene expression environment of reproductive development in T. ravennae.

RESULTS

During the early phases of inflorescence development, multiple key canonical floral integrators and pathways were identified. Annotations of type II subfamily of MADS-box transcription factors, in particular, were over-represented in the GO enrichment analyses and tests for differential expression (FDR p-value < 0.05). The differential expression of floral integrators observed in the early phases of inflorescence development diminished prior to inflorescence determinacy regulation. Differential expression analysis did not identify many unique genes at mid-inflorescence development stages, though typical biological processes involved in plant growth and development expressed abundantly. The increase in inflorescence determinacy regulatory elements and putative homeotic floral development unigenes at mid-inflorescence development coincided with the expression of multiple meiosis annotations and multicellular organism developmental processes. Analysis of seed development identified multiple unigenes involved in oxidative-reductive processes.

CONCLUSION

Reproduction in grasses is a dynamic system involving the sequential coordination of complex gene regulatory networks and developmental processes. This research identified differentially expressed transcripts associated with floral induction, inflorescence development, and seed development in T. ravennae. These results provide insights into the molecular regulation of reproductive development and provide a foundation for future investigations and analyses, including genome annotation, functional genomics characterization, gene family evolutionary studies, comparative genomics, and precision breeding.

Transcriptome profiling of developing leaf and shoot apices to reveal the molecular mechanism and co-expression genes responsible for the wheat heading date

Background

Wheat is one of the most widely planted crops worldwide. The heading date is important for wheat environmental adaptability, as it not only controls flowering time but also determines the yield component in terms of grain number per spike.

Results

In this research, homozygous genotypes with early and late heading dates derived from backcrossed progeny were selected to conduct RNA-Seq analysis at the double ridge stage (W2.0) and androgynous primordium differentiation stage (W3.5) of the leaf and apical meristem, respectively. In total, 18,352 differentially expressed genes (DEGs) were identified, many of which are strongly associated with wheat heading date genes. Gene Ontology (GO) enrichment analysis revealed that carbohydrate metabolism, trehalose metabolic process, photosynthesis, and light reaction are closely related to the flowering time regulation pathway. Based on MapMan metabolic analysis, the DEGs are mainly involved in the light reaction, hormone signaling, lipid metabolism, secondary metabolism, and nucleotide synthesis. In addition, 1,225 DEGs were annotated to 45 transcription factor gene families, including LFY, SBP, and MADS-box transcription factors closely related to flowering time. Weighted gene co-expression network analysis (WGCNA) showed that 16, 336, 446, and 124 DEGs have biological connections with Vrn1-5 A, Vrn3-7B, Ppd-1D, and WSOC1, respectively. Furthermore, TraesCS2D02G181400 encodes a MADS-MIKC transcription factor and is co-expressed with Vrn1-5 A, which indicates that this gene may be related to flowering time.

Conclusions

RNA-Seq analysis provided transcriptome data for the wheat heading date at key flower development stages of double ridge (W2.0) and androgynous primordium differentiation (W3.5). Based on the DEGs identified, co-expression networks of key flowering time genes in Vrn1-5 A, Vrn3-7B, WSOC1, and Ppd-1D were established. Moreover, we discovered a potential candidate flowering time gene, TraesCS2D02G181400. Taken together, these results serve as a foundation for further study on the regulatory mechanism of the wheat heading date.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07797-7.