Transcriptome Analysis Reveals Important Candidate Genes Related to Nutrient Reservoir, Carbohydrate Metabolism, and Defence Proteins during Grain Development of Hexaploid Bread Wheat and Its Diploid Progenitors

Starch and storage protein dynamics in the developing and matured grains of durum wheat and diploid progenitor species

Durum wheat, less immunogenically intolerant than bread wheat, originates from diploid progenitors known for nutritional quality and stress tolerance. Present study involves the analysis of major grain parameters, viz. size, weight, sugar, starch, and protein content of Triticum durum (AABB genome) and its diploid progenitors, Triticum monococcum (AA genome) and Aegilops speltoides (BB genome). Samples were collected during 2-5 weeks after anthesis (WAA), and at maturity. The investigation revealed that T. durum displayed the maximum grain size and weight. Expression analysis of Grain Weight 2 (GW2) and Glutamine Synthase (GS2), negative and positive regulators of grain weight and size, respectively, revealed higher GW2 expression in Ae. speltoides and higher GS2 expression in T. durum. Further we explored total starch, sugar and protein content, observing higher levels of starch and sugar in durum wheat while AA genome species exhibited higher protein content dominated by the fractions of albumin/globulin. HPLC profiling revealed unique sub-fractions in all three genome species. Additionally, a comparative transcriptome analysis also corroborated with the starch and protein content in the grains. This study provides valuable insights into the genetic and biochemical distinctions among durum wheat and its diploid progenitors, offering a foundation for their nutritional composition.

Genetic dissection of protein and starch during wheat grain development using QTL mapping and GWAS

Protein, starch, and their components are important for wheat grain yield and end-products, which are affected by wheat grain development. Therefore, QTL mapping and a genome-wide association study (GWAS) of grain protein content (GPC), glutenin macropolymer content (GMP), amylopectin content (GApC), and amylose content (GAsC) were performed on wheat grain development at 7, 14, 21, and 28 days after anthesis (DAA) in two environments using a recombinant inbred line (RIL) population of 256 stable lines and a panel of 205 wheat accessions. A total of 29 unconditional QTLs, 13 conditional QTLs, 99 unconditional marker-trait associations (MTAs), and 14 conditional MTAs significantly associated (p < 10^-4) with four quality traits were found to be distributed on 15 chromosomes, with the phenotypic variation explained (PVE) ranging from 5.35% to 39.86%. Among these genomic variations, three major QTLs [QGPC3B, QGPC2A, and QGPC(S3|S2)3B] and SNP clusters on the 3A and 6B chromosomes were detected for GPC, and the SNP TA005876-0602 was stably expressed during the three periods in the natural population. The QGMP3B locus was detected five times in three developmental stages in two environments with 5.89%-33.62% PVE, and SNP clusters for GMP content were found on the 3A and 3B chromosomes. For GApC, the QGApC3B.1 locus had the highest PVE of 25.69%, and SNP clusters were found on chromosomes 4A, 4B, 5B, 6B, and 7B. Four major QTLs of GAsC were detected at 21 and 28 DAA. Most interestingly, both QTL mapping and GWAS analysis indicated that four chromosomes (3B, 4A, 6B, and 7A) were mainly involved in the development of protein, GMP, amylopectin, and amylose synthesis. Of these, the wPt-5870-wPt-3620 marker interval on chromosome 3B seemed to be most important because it played an important role in the synthesis of GMP and amylopectin before 7 DAA, in the synthesis of protein and GMP from 14 to 21 DAA, and in the development of GApC and GAsC from 21 to 28 DAA. Using the annotation information of IWGSC Chinese Spring RefSeq v1.1 genome assembly, we predicted 28 and 69 candidate genes for major loci from QTL mapping and GWAS, respectively. Most of them have multiple effects on protein and starch synthesis during grain development. These results provide new insights and information for the potential regulatory network between grain protein and starch synthesis.

TaMADS29 interacts with TaNF-YB1 to synergistically regulate early grain development in bread wheat

Grain development is a crucial determinant of yield and quality in bread wheat (Triticum aestivum L.). However, the regulatory mechanisms underlying wheat grain development remain elusive. Here we report how TaMADS29 interacts with TaNF-YB1 to synergistically regulate early grain development in bread wheat. The tamads29 mutants generated by CRISPR/Cas9 exhibited severe grain filling deficiency, coupled with excessive accumulation of reactive oxygen species (ROS) and abnormal programmed cell death that occurred in early developing grains, while overexpression of TaMADS29 increased grain width and 1,000-kernel weight. Further analysis revealed that TaMADS29 interacted directly with TaNF-YB1; null mutation in TaNF-YB1 caused grain developmental deficiency similar to tamads29 mutants. The regulatory complex composed of TaMADS29 and TaNF-YB1 exercises its possible function that inhibits the excessive accumulation of ROS by regulating the genes involved in chloroplast development and photosynthesis in early developing wheat grains and prevents nucellar projection degradation and endosperm cell death, facilitating transportation of nutrients into the endosperm and wholly filling of developing grains. Collectively, our work not only discloses the molecular mechanism of MADS-box and NF-Y TFs in facilitating bread wheat grain development, but also indicates that caryopsis chloroplast might be a central regulator of grain development rather than merely a photosynthesis organelle. More importantly, our work offers an innovative way to breed high-yield wheat cultivars by controlling the ROS level in developing grains.

Transcriptome and Proteome Analysis Revealed the Influence of High-Molecular-Weight Glutenin Subunits (HMW-GSs) Deficiency on Expression of Storage Substances and the Potential Regulatory Mechanism of HMW-GSs

The processing quality of wheat is affected by seed storage substances, such as protein and starch. High-molecular-weight glutenin subunits (HMW-GSs) are the major components of wheat seed storage proteins (SSPs); they are also key determinators of wheat end-use quality. However, the effects of HMW-GSs absence on the expression of other storage substances and the regulation mechanism of HMW-GSs are still limited. Previously, a wheat transgenic line LH-11 with complete deletions of HMW-GSs was obtained through introducing an exogenous gene Glu-1E^bx to the wild-type cultivar Bobwhite by transgenic approach. In this study, comparative seed transcriptomics and proteomics of transgenic and non-transgenic lines at different seed developmental stages were carried out to explore the changes in genes and proteins and the underlying regulatory mechanism. Results revealed that a number of genes, including genes related to SSPs, carbohydrates metabolism, amino acids metabolism, transcription, translation, and protein process were differentially enriched. Seed storage proteins displayed differential expression patterns between the transgenic and non-transgenic line, a major rise in the expression levels of gliadins were observed at 21 and 28 days post anthesis (DPA) in the transgenic line. Changes in expressions of low-molecular-weight glutenins (LMW-GSs), avenin-like proteins (ALPs), lipid transfer proteins (LTPs), and protease inhibitors (PIs) were also observed. In addition, genes related to carbohydrate metabolism were differentially expressed, which probably leads to a difference in starch component and deposition. A list of gene categories participating in the accumulation of SSPs was proposed according to the transcriptome and proteome data. Six genes from the MYB and eight genes from the NAC transcription families are likely important regulators of HMW-GSs accumulation. This study will provide data support for understanding the regulatory network of wheat storage substances. The screened candidate genes can lay a foundation for further research on the regulation mechanism of HMW-GSs.

Understanding plant-microbe interaction of rice and soybean with two contrasting diazotrophic bacteria through comparative transcriptome analysis

Understanding the beneficial plant-microbe interactions is becoming extremely critical for deploying microbes imparting plant fitness and achieving sustainability in agriculture. Diazotrophic bacteria have the unique ability to survive without external sources of nitrogen and simultaneously promote host plant growth, but the mechanisms of endophytic interaction in cereals and legumes have not been studied extensively. We have studied the early interaction of two diazotrophic bacteria, Gluconacetobacter diazotrophicus (GAB) and Bradyrhizobium japonicum (BRH), in 15-day-old seedlings of rice and soybean up to 120 h after inoculation (hai) under low-nitrogen medium. Root colonization of GAB in rice was higher than that of BRH, and BRH colonization was higher in soybean roots as observed from the scanning electron microscopy at 120 hai. Peroxidase enzyme was significantly higher at 24 hai but thereafter was reduced sharply in soybean and gradually in rice. The roots of rice and soybean inoculated with GAB and BRH harvested from five time points were pooled, and transcriptome analysis was executed along with control. Two pathways, "Plant pathogen interaction" and "MAPK signaling," were specific to Rice-Gluconacetobacter (RG), whereas the pathways related to nitrogen metabolism and plant hormone signaling were specific to Rice-Bradyrhizobium (RB) in rice. Comparative transcriptome analysis of the root tissues revealed that several plant-diazotroph-specific differentially expressed genes (DEGs) and metabolic pathways of plant-diazotroph-specific transcripts, viz., chitinase, brassinosteroid, auxin, Myeloblastosis (MYB), nodulin, and nitrate transporter (NRT), were common in all plant-diazotroph combinations; three transcripts, viz., nitrate transport accessory protein (NAR), thaumatin, and thionin, were exclusive in rice and another three transcripts, viz., NAC (NAM: no apical meristem, ATAF: Arabidopsis thaliana activating factor, and CUC: cup-shaped cotyledon), ABA (abscisic acid), and ammonium transporter, were exclusive in soybean. Differential expression of these transcripts and reduction in pathogenesis-related (PR) protein expression show the early interaction. Based on the interaction, it can be inferred that the compatibility of rice and soybean is more with GAB and BRH, respectively. We propose that rice is unable to identify the diazotroph as a beneficial microorganism or a pathogen from an early response. So, it expressed the hypersensitivity-related transcripts along with PR proteins. The molecular mechanism of diazotrophic associations of GAB and BRH with rice vis-à-vis soybean will shed light on the basic understanding of host responses to beneficial microorganisms.

Heat Stress Tolerance 2 confers basal heat stress tolerance in allohexaploid wheat (Triticum aestivum L.)

Heat stress substantially reduces the yield potential of wheat (Triticum aestivum L.), one of the most widely cultivated staple crops, and greatly threatens global food security in the context of global warming. However, few studies have explored the heat stress tolerance (HST)-related genetic resources in wheat. Here, we identified and fine-mapped a wheat HST locus, TaHST2, which is indispensable for HST in both the vegetative and reproductive stages of the wheat life cycle. The studied pair of near isogenic lines (NILs) exhibited diverse morphologies under heat stress, based on which we mapped TaHST2 to a 485 kb interval on chromosome arm 4DS. Under heat stress, TaHST2 confers a superior conversion rate from soluble sugars to starch in wheat grains, resulting in faster grain filling and a higher yield potential. A further exploration of genetic resources indicated that TaHST2 underwent strong artificial selection during wheat domestication, suggesting it is an essential locus for basal HST in wheat. Our findings provide deeper insights into the genetic basis of wheat HST and might be useful for global efforts to breed heat-stress-tolerant cultivars.

Grain Transcriptome Dynamics Induced by Heat in Commercial and Traditional Bread Wheat Genotypes

High temperature (HT) events have negative impact on wheat grains yield and quality. Transcriptome profiles of wheat developing grains of commercial genotypes (Antequera and Bancal) and landraces (Ardito and Magueija) submitted to heatwave-like treatments during grain filling were evaluated. Landraces showed significantly more differentially expressed genes (DEGs) and presented more similar responses than commercial genotypes. DEGs were more associated with transcription and RNA and protein synthesis in Antequera and with metabolism alterations in Bancal and landraces. Landraces upregulated genes encoding proteins already described as HT responsive, like heat shock proteins and cupins. Apart from the genes encoding HSP, two other genes were upregulated in all genotypes, one encoding for Adenylate kinase, essential for the cellular homeostasis, and the other for ferritin, recently related with increased tolerance to several abiotic stress in Arabidopsis. Moreover, a NAC transcription factor involved in plant development, known to be a negative regulator of starch synthesis and grain yield, was found to be upregulated in both commercial varieties and downregulated in Magueija landrace. The detected diversity of molecular processes involved in heat response of commercial and traditional genotypes contribute to understand the importance of genetic diversity and relevant pathways to cope with these extreme events.

Unravelling consensus genomic regions associated with quality traits in wheat using meta-analysis of quantitative trait loci

Meta-analysis in wheat for three major quality traits identified 110 meta-QTL (MQTL) with reduced confidence interval (CI). Five GWAS validated MQTL (viz., 1A.1, 1B.2, 3B.4, 5B.2, and 6B.2), each involving more than 20 initial QTL and reduced CI (95%) (< 2 cM), were selected for quality breeding programmes. Functional characterization including candidate gene mining and expression analysis discovered 44 high confidence candidate genes associated with quality traits. A meta-analysis of quantitative trait loci (QTL) associated with dough rheology properties, nutritional traits, and processing quality traits was conducted in wheat. For this purpose, as many as 2458 QTL were collected from 50 interval mapping studies published during 2013-2020. Of the total QTL, 1126 QTL were projected onto the consensus map saturated with 249,603 markers which led to the identification of 110 meta-QTL (MQTL). These MQTL exhibited an 18.84-fold reduction in the average CI compared to the average CI of the initial QTL (ranging from 14.87 to 95.55 cM with an average of 40.35 cM). Of the 110, 108 MQTL were physically anchored to the wheat reference genome, including 51 MQTL verified with marker-trait associations (MTAs) reported from earlier genome-wide association studies. Candidate gene (CG) mining allowed the identification of 2533 unique gene models from the MQTL regions. In-silico expression analysis discovered 439 differentially expressed gene models with > 2 transcripts per million expressions in grains and related tissues, which also included 44 high-confidence CGs involved in the various cellular and biochemical processes related to quality traits. Nine functionally characterized wheat genes associated with grain protein content, high-molecular-weight glutenin, and starch synthase enzymes were also found to be co-localized with some of the MQTL. Synteny analysis between wheat and rice MQTL regions identified 23 wheat MQTL syntenic to 16 rice MQTL associated with quality traits. Furthermore, 64 wheat orthologues of 30 known rice genes were detected in 44 MQTL regions. Markers flanking the MQTL identified in the present study can be used for marker-assisted breeding and as fixed effects in the genomic selection models for improving the prediction accuracy during quality breeding. Wheat orthologues of rice genes and other CGs available from MQTLs can be promising targets for further functional validation and to better understand the molecular mechanism underlying the quality traits in wheat.

Genome-wide identification of seed storage protein gene regulators in wheat through coexpression analysis

Only a few transcriptional regulators of seed storage protein (SSP) genes have been identified in common wheat (Triticum aestivum L.). Coexpression analysis could be an efficient approach to characterize novel transcriptional regulators at the genome-scale considering the correlated expression between transcriptional regulators and target genes. As the A genome donor of common wheat, Triticum urartu is more suitable for coexpression analysis than common wheat considering the diploid genome and single gene copy. In this work, the transcriptome dynamics in endosperm of T. urartu throughout grain filling were revealed by RNA-Seq analysis. In the coexpression analysis, a total of 71 transcription factors (TFs) from 23 families were found to be coexpressed with SSP genes. Among these TFs, TuNAC77 enhanced the transcription of SSP genes by binding to cis-elements distributed in promoters. The homolog of TuNAC77 in common wheat, TaNAC77, shared an identical function, and the total SSPs were reduced by about 24% in common wheat when TaNAC77 was knocked down. This is the first genome-wide identification of transcriptional regulators of SSP genes in wheat, and the newly characterized transcriptional regulators will undoubtedly expand our knowledge of the transcriptional regulation of SSP synthesis.