Heavy metal pollution is a serious global environmental problem as it adversely affects plant growth and genetic variation. It also alters the composition and activity of soil microbial communities. The objectives of this study were to determine the soil microbial diversity, bermudagrass genetic variation in Cd contaminated or uncontaminated soils from Hunan province of China, and to evaluate Cd-tolerance of bermudagrass at different soils. The Biolog method, hydroponic experiments and simple sequence repeat markers were used to assess the functional diversity of microorganisms, Cd-tolerance and the genetic diversity of bermudagrass, respectively. Four of the sampling sites were heavily contaminated with heavy metals. The total bioactivity, richness, and microbial diversity decreased with increasing concentration of heavy metal. The hydroponic experiment revealed that bermudagrass populations collected from polluted sites have evolved, encompassing the feature of a higher resistance to Cd toxicity. Higher genetic diversity was observed to be more in contaminated populations than in uncontaminated populations. Heavy metal pollution can result in adverse effects on plant growth, soil microbial diversity and activity, and apparently has a stronger impact on the genetic structure. The results of this study provide new insights and a background to produce a genetic description of populations in a species that is suitable for use in phytoremediation practices.

Drought is the most serious abiotic stress limiting rice production, and deep root is the key contributor to drought avoidance. However, the genetic mechanism regulating the development of deep roots is largely unknown. In this study, the transcriptomes of 74 root samples from 37 rice varieties, representing the extreme genotypes of shallow or deep rooting, were surveyed by RNA-seq. The 13,242 differentially expressed genes (DEGs) between deep rooting and shallow rooting varieties (H vs. L) were enriched in the pathway of genetic information processing and metabolism, while the 1,052 DEGs between the deep roots and shallow roots from each of the plants (D vs. S) were significantly enriched in metabolic pathways especially energy metabolism. Ten quantitative trait transcripts (QTTs) were identified and some were involved in energy metabolism. Forty-nine candidate DEGs were confirmed by qRT-PCR and microarray. Through weighted gene co-expression network analysis (WGCNA), we found 18 hub genes. Surprisingly, all these hub genes expressed higher in deep roots than in shallow roots, furthermore half of them functioned in energy metabolism. We also estimated that the ATP production in the deep roots was faster than shallow roots. Our results provided a lot of reliable candidate genes to improve deep rooting, and firstly highlight the importance of energy metabolism to the development of deep roots.

Tea leaf color is not only important from an aesthetics standpoint but is also related to tea quality. To investigate the molecular mechanisms that determine tea leaf color, we examined Camellia sinensis cv. 'Anjin Baicha' (an albino tea cultivar) by tandem mass tag isobaric labeling to generate a high-resolution proteome and acetyl-proteome atlas of three leaf developmental stages. We identified a total of 7,637 proteins and quantified 6,256; of these, 3,232 were classified as differentially accumulated proteins (DAPs). We also identified 3,161 lysine acetylation sites in 1,752 proteins and quantified 2,869 in 1,612 proteins. The acetylation levels at 468 sites were significantly altered across the three developmental stages during periodic albinism; the corresponding proteins were associated with a variety of biological processes. Interestingly, a large number of DAPs and acetylated proteins with increased/decreased acetylation were related to photosynthesis and secondary metabolite biosynthetic pathways, suggesting that the accumulation or acetylation level of these proteins regulates periodic albinism in 'Anjin Baicha.' Additionally, overlap between succinylome and acetylome among three 'Anjin Baicha' developmental stages were found. These data provide important insight into the mechanisms of leaf coloration in the tea plant. The mass spectrometry data have been deposited to Proteome X change via the PRIDE partner repository with the data set identifier PXD008134.

In contrast to wild species, drought-tolerance in crops requires a fully functional metabolism during drought (particularly photosynthetic processes). However, the link between drought-tolerance, photosynthetic regulation during drought, and the associated transcript and metabolic foundation, remains largely unknown. For this study, we used two rice cultivars with contrasting drought-tolerance (the drought-intolerant cultivar IRAT109 and the drought-tolerant cultivar IAC1246) to explore transcript and metabolic responses to long-term drought. The drought-tolerant cultivar represented higher osmotic adjustment and antioxidant capacity, as well as higher relative photosynthesis rate under a progressive drought stress occurred in a modified field with shallow soil-layers. A total of 4059 and 2677 differentially expressed genes (DEGs) were identified in IRAT109 and IAC1246 between the drought and well-watered conditions, respectively. A total of 69 and 47 differential metabolites (DMs) were identified between the two treatments in IRAT109 and IAC1246, respectively. Compared to IRAT109, the DEGs of IAC1246 displayed enhanced regulatory amplitude during drought. We found significant correlations between DEGs and the osmolality and total antioxidant capacity (AOC) of both cultivars. During the early stages of drought, we detected up-regulation of DEGs in IAC1246 related to photosynthesis, in accordance with its higher relative photosynthesis rate. The contents of six differential metabolites were correlated with the osmotic potential and AOC. Moreover, they were differently regulated between the two cultivars. Particularly, up-regulations of 4-hydroxycinnamic acid and ferulic acid were consistent with the performance of photosynthesis-related DEGs at the early stages of drought in IAC1246. Therefore, 4-hydroxycinnamic acid and ferulic acid were considered as key metabolites for rice drought-tolerance. DEGs involved in pathways of these metabolites are expected to be good candidate genes to improve drought-tolerance. In conclusion, well-maintained photosynthesis under drought should contribute to improved drought-tolerance in rice. Metabolites play vital roles in protecting photosynthesis under dehydration via osmotic adjustments and/or antioxidant mechanisms. A metabolite-based method was thus an effective way to explore drought candidate genes. Metabolic accompanied by transcript responses to drought stress should be further studied to find more useful metabolites, pathways, and genes.

Melatonin (N-acetyl-5-methoxytryptamine) plays critical roles in plant growth and development and during the response to multiple abiotic stresses. However, the roles of melatonin in plant response to K+ deficiency remain largely unknown. In the present study, we observed that the endogenous melatonin contents in bermudagrass were remarkably increased by low K+ (LK) treatment, suggesting that melatonin was involved in bermudagrass response to LK stress. Further phenotype analysis revealed that exogenous melatonin application conferred Bermudagrass enhanced tolerance to LK stress. Interestingly, exogenous melatonin application also promoted bermudagrass growth and development at normal condition. Furthermore, the K+ contents measurement revealed that melatonin-treated plants accumulated more K+ in both shoot (under both control and LK condition) and root tissues (under LK condition) compared with those of melatonin non-treated plants. Expression analysis indicated that the transcripts of K+ transport genes were significantly induced by exogenous melatonin treatment in bermudagrass under both control and LK stress conditions, especially under a combined treatment of LK stress and melatonin, which may increase accumulation of K+ content profoundly under LK stress and thereby contributed to the LK-tolerant phenotype. In addition, we investigated the role of melatonin in the regulation of photosystem II (PSII) activities under LK stress. The chlorophyll fluorescence transient (OJIP) curves were obviously higher in plants grown in LK with melatonin (LK+Mel) than those of plants grown in LK medium without melatonin application for 1 or 2 weeks, suggesting that melatonin plays important roles in PSII against LK stress. After a combined treatment of LK stress and melatonin, the values for performance indexes (PIABS, PITotal, and PICS), flux ratios (φP0, ΨE0, and φE0) and specific energy fluxes (ETO/RC) were significantly improved compared with those of LK stress alone, suggesting that melatonin plays positive roles in protecting PSII activity under LK stress. Collectively, this study reveals an important role of melatonin in regulating bermudagrass response to LK stress.

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.

Root growth and the architecture of the root system in Arabidopsis are largely determined by root meristematic activity. Legume roots show strong developmental plasticity in response to both abiotic and biotic stimuli, including symbiotic rhizobia. However, a global analysis of gene regulation in the root meristem of soybean plants is lacking. In this study, we performed a global analysis of the small RNA transcriptome of root tips from soybean seedlings grown under normal and salt stress conditions. In total, 71 miRNA candidates, including known and novel variants of 59 miRNA families, were identified. We found 66 salt-responsive miRNAs in the soybean root meristem; among them, 22 are novel miRNAs. Interestingly, we found auxin-responsive cis-elements in the promoters of many salt-responsive miRNAs, implying that these miRNAs may be regulated by auxin and auxin signaling plays a key role in regulating the plasticity of the miRNAome and root development in soybean. A functional analysis of miR399, a salt-responsive miRNA in the root meristem, indicates the crucial role of this miRNA in modulating soybean root developmental plasticity. Our data provide novel insight into the miRNAome-mediated regulatory mechanism in soybean root growth under salt stress.

Panicle degradation causes severe yield reduction in rice. There are two main types of panicle degradation: apical spikelet abortion and basal degeneration. In this study, we isolated and characterized the apical panicle abortion mutant apical spikelet abortion (asa), which exhibits degeneration and defects in the apical spikelets. This mutant had a pleiotropic phenotype, characterized by reduced plant height, increased tiller number, and decreased pollen fertility. Map-based cloning revealed that OsASA encodes a boric acid channel protein that showed the highest expression in the inflorescence, peduncle, and anther. RNA-seq analysis of the asa mutant vs wild-type (WT) plants revealed that biological processes related to reactive oxygen species (ROS) homeostasis and salicylic acid (SA) metabolism were significantly affected. Furthermore, the asa mutants had an increased SA level and H2O2 accumulation in the young panicles compared to the WT plants. Moreover, the SA level and the expression of OsPAL3, OsPAL4, and OsPAL6 genes (related to SA biosynthesis) were significantly increased under boron-deficient conditions in the asa mutant and in OsASA-knockout plants. Collectively, these results suggest that the boron distribution maintained by OsASA is required for normal panicle development in a process that involves modulating ROS homeostasis and SA biosynthesis.

Tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most important economic crops with multiple mutants. Recently, we found a special tea germplasm that has an aberrant tissue on its branches. To figure out whether this aberrant tissue is associated with floral bud (FB) or dormant bud (DB), we performed tissue section, transcriptome sequencing, and metabolomic analysis of these tissues. Longitudinal sections indicated the aberrant tissue internal structure was more like a special bud (SB), but was similar to that of DB. Transcriptome data analysis showed that the number of heterozygous and homozygous SNPs was significantly different in the aberrant tissue compared with FB and DB. Further, by aligning the unmapped sequences of the aberrant tissue to the Non-Redundant Protein Sequences (NR) database, we observed that 36.13% of unmapped sequences were insect sequences, which suggested that the aberrant tissue might be a variation of dormant bud tissue influenced by the interaction of tea plants and insects or pathogens. Metabolomic analysis showed that the differentially expressed metabolites (DEMs) between the aberrant tissue and DB were significantly enriched in the metabolic pathways of biosynthesis of plant hormones and biosynthesis of phenylpropanoids. Subsequently, we analyzed the differentially expressed genes (DEGs) in the above mentioned two tissues, and the results indicated that photosynthetic capacity in the aberrant tissue was reduced, whereas the ethylene, salicylic acid and jasmonic acid signaling pathways were activated. We speculated that exogenous infection induced programmed cell death (PCD) and increased the lignin content in dormant buds of tea plants, leading to the formation of this aberrant tissue. This study advanced our understanding of the interaction between plants and insects or pathogens, providing important clues about biotic stress factors and key genes that lead to mutations and formation of the aberrant tissue.

There are numerous non-volatile metabolites in the fresh shoots of tea plants. However, we know little about the complex relationship between the content of these metabolites and their gene expression levels. In investigating this, this study involved non-volatile metabolites from 68 accessions of tea plants that were detected and identified using untargeted metabolomics. The tea accessions were divided into three groups from the results of a principal component analysis based on the relative content of the metabolites. There were differences in variability between the primary and secondary metabolites. Furthermore, correlations among genes, gene metabolites, and metabolites were conducted based on Pearson's correlation coefficient (PCC) values. This study offered several significant insights into the co-current network of genes and metabolites in the global genetic background. Thus, the study is useful for providing insights into the regulatory relationship of the genetic basis for predominant metabolites in fresh tea shoots.

The gibberellin-responsive dwarfing gene Rht12 can significantly reduce plant height without changing seedling vigor and substantially increase ear fertility in bread wheat (Triticum aestivum. L). However, Rht12 delays heading date and anthesis date, hindering the use of Rht12 in wheat improvement. To promote early flowering of the Rht12 dwarf plants, the photoperiod-insensitive allele Ppd-D1a was introduced through a cross between Jinmai47 (Ppd-D1a) and Karcagi (Rht12). The results showed that Ppd-D1a can rescue the delaying effect of Rht12 on flowering time and promote earlier flowering by 9.0 days (163.2°Cd) in the Rht12 dwarf plants by shortening the late reproduction phase. Plant height was reduced by Rht12 (43.2%) and Ppd-D1a (10.9%), achieving dwarf plants with higher lodging resistance. Ear fertility, like the grain number per spike, was significantly increased by Rht12 (21.3%), while it was reduced by Ppd-D1a (6.5%). However, thousand kernel weight was significantly reduced by Rht12 (12.9%) but significantly increased by Ppd-D1a (16.9%). Finally, plant yield was increased by 16.4 and 8.2%, and harvest index was increased by 24.9 and 15.4% in the Rht12 dwarf lines and tall lines with Ppd-D1a, respectively. Clearly, there was an additive interaction between Rht12 and Ppd-D1 and the introduction of Ppd-D1a advanced the flowering time and improved the yield traits of Rht12 dwarf plants, suggesting that the combination of Rht12 and Ppd-D1a would be conducive to the successful use of Rht12 in wheat breeding programs.

As a typical warm-season grass, Bermudagrass [Cynodon dactylon (L).Pers.] is widely applied in turf systems and animal husbandry. However, cold temperature is a key factor limiting resource utilization for Bermudagrass. Therefore, it is relevant to study the mechanisms by which Burmudagrass responds to cold. Melatonin is a crucial animal and plant hormone that is responsible for plant abiotic stress responses. The objective of this study was to investigate the role of melatonin in cold stress response of Bermudagrass. Wild Bermudagrass pre-treated with 100 μM melatonin was subjected to different cold stress treatments (-5°C for 8 h with or without cold acclimation). The results showed lower malondialdehyde (MDA) and electrolyte leakage (EL) values, higher levels of chlorophyll, and greater superoxide dismutase and peroxidase activities after melatonin treatment than those in non-melatonin treatment under cold stress. Analysis of chlorophyll a revealed that the chlorophyll fluorescence transient (OJIP) curves were higher after treatment with melatonin than that of non-melatonin treated plants under cold stress. The values of photosynthetic fluorescence parameters increased after treatment with melatonin under cold stress. The analysis of metabolism showed alterations in 46 metabolites in cold-stressed plants after melatonin treatment. Among the measured metabolites, five sugars (arabinose, mannose, glucopyranose, maltose, and turanose) and one organic acid (propanoic acid) were significantly increased. However, valine and threonic acid contents were reduced in melatonin-treated plants. In summary, melatonin maintained cell membrane stability, increased antioxidant enzymes activities, improved the process of photosystem II, and induced alterations in Bermudagrass metabolism under cold stress.

As a representative warm-season grass, Bermudagrass [Cynodon dactylon (L). Pers.] is widely used in turf systems. However, low temperature remarkably limits its growth and distribution. ABA is a crucial phytohormone that has been reported to regulate much important physiological and biochemical processes in plants under abiotic stress. Therefore, the objective of this study was to figure out the effects of ABA on the cold-sensitive (S) and cold-resistant (R) Bermudagrass genotypes response to cold stress. In this study, the plants were treated with 100 μM ABA solution and exposed to 4°C temperature. After 7 days of cold treatment, the electrolyte leakage (EL), malonaldehyde (MDA) and H2O2 content were significantly increased in both genotypes compared with control condition, and these values were higher in R genotype than those of S genotype, respectively. By contrast, exogenous ABA application decreased the electrolyte leakage (EL), MDA and H2O2 content in both genotypes compared with those plants without ABA treatment under cold treatment condition. In addition, exogenous ABA application increased the levels of chlorophyll a fluorescence transient curve for both genotypes, and it was higher in R genotype than that of S genotype. Analysis of photosynthetic fluorescence parameters revealed that ABA treatment improved the performance of photosystem II under cold condition, particularly for the R genotype. Moreover, cold stress significantly increased δ13C values for both genotypes, while it was alleviated by exogenous ABA. Additionally, exogenous ABA application altered the expression of ABA- or cold related genes, including ABF1, CBF1, and LEA. In summary, exogenous ABA application enhanced cold resistance of both genotypes by maintaining cell membrane stability, improving the process of photosystem II, increasing carbon isotopic fractionation under cold stress, and more prominently in R genotype compared with S genotype.

The type III secretion system (T3SS) in many Gram-negative bacterial pathogens is regarded as the most critical virulence determinant and an attractive target for novel anti-virulence drugs. In this study, we constructed a T3SS secretion reporter containing the β-lactamase gene fused with a signal peptide sequence of the T3SS effector gene, and established a high-throughput screening system for T3SS inhibitors in the plant pathogenic bacterium Acidovorax citrulli. From a library of 12,000 chemical compounds, we identified a series of benzyloxy carbonimidoyl dicyanide (BCD) derivatives that effectively blocked T3SS-dependent β-lactamase secretion. Substitution of halogens or nitro groups at the para-position on the benzene ring contributed to an increased inhibitory activity. One representative compound, BCD03 (3,4-dichloro-benzyloxy carbonimidoyl dicyanide), dramatically reduced pathogenicity of A. citrulli on melon seedlings, and attenuated hypersensitive responses in the non-host Nicotiana tabacum caused by pathogenic bacteria A. citrulli, Xanthomonas oryzae pv. oryzae and Pseudomonas syringae pv. tomato at sub-MIC concentrations. Western blotting assay further confirmed that BCD03 inhibited effector secretion from the above bacteria via T3SS in the liquid medium. Taken together, our data suggest that BCD derivatives act as novel inhibitors of T3SS in multiple plant bacterial pathogens.

Stomatal density is an important factor that determines the efficiency of plant gas exchange and water transpiration. Through forward genetics, we screened a mutant rice stomata developmental defect 1 (rsd1-1) with decreased stomatal density and clustered stomata in rice (Oryza sativa). After the first asymmetric division, some of the larger sister cells undergo an extra asymmetric division to produce a small cell neighboring guard mother cell. Some of these small cells develop into stomata, which leads to stomatal clustering, and the rest arrested or developed into pavement cell. After map-based cloning, we found the protein encoded by this gene containing DUF630 and DUF632 domains. Evolutionary analysis showed that the DUF630/632 gene family differentiated earlier in land plants. It was found that the deletion of RSD1 would lead to the disorder of gene expression regarding stomatal development, especially the expression of stomatal density and distribution 1 (OsSDD1). Through the construction of OsSDD1 deletion mutants by CRISPR-Cas9, we found that, similar to rsd1 mutants, the ossdd1 mutants have clustered stomata and extra small cells adjacent to the stomata. OsSDD1 and RSD1 are both required for inhibiting ectopic asymmetric cell divisions (ACDs) and clustered stomata. By dehydration stress assay, the decreased stomatal density of rsd1 mutants enhanced their dehydration avoidance. This study characterized the functions of RSD1 and OsSDD1 in rice stomatal development. Our findings will be helpful in developing drought-resistant crops through controlling the stomatal density.

Continuous identification and application of novel resistance genes against stripe rust are of great importance for wheat breeding. Wild emmer wheat, Triticum dicoccoides, has adapted to a broad range of environments and is a valuable genetic resource that harbors important beneficial traits, including resistance to stripe rust caused by Puccinia striiformis f. sp. tritici (Pst). However, there has been a lack of systematic exploration of genes against Pst races in wild emmer wheat.

Laccases are widely distributed in plant kingdom catalyzing the polymerization of lignin monolignols. Rosmarinic acid (RA) has a lignin monolignol-like structure and is converted into salvianolic acid B (SAB), which is a representatively effective hydrophilic compound of a well-known medicinal plant Salvia miltiorrhiza and also the final compound of phenolic acids metabolism pathway in the plant. But the roles of laccases in the biosynthesis of SAB are poorly understood. This work systematically characterizes S. miltiorrhiza laccase (SmLAC) gene family and identifies the SAB-specific candidates. Totally, 29 laccase candidates (SmLAC1-SmLAC29) are found to contain three signature Cu-oxidase domains. They present relatively low sequence identity and diverse intron-exon patterns. The phylogenetic clustering of laccases from S. miltiorrhiza and other ten plants indicates that the 29 SmLACs can be divided into seven groups, revealing potential distinct functions. Existence of diverse cis regulatory elements in the SmLACs promoters suggests putative interactions with transcription factors. Seven SmLACs are found to be potential targets of miR397. Putative glycosylation sites and phosphorylation sites are identified in SmLAC amino acid sequences. Moreover, the expression profile of SmLACs in different organs and tissues deciphers that 5 SmLACs (SmLAC7/8/20/27/28) are expressed preferentially in roots, adding the evidence that they may be involved in the phenylpropanoid metabolic pathway. Besides, silencing of SmLAC7, SmLAC20 and SmLAC28, and overexpression of SmLAC7 and SmLAC20 in the hairy roots of S. miltiorrhiza result in diversification of SAB, signifying that SmLAC7 and SmLAC20 take roles in SAB biosynthesis. The results of this study lay a foundation for further elucidation of laccase functions in S. miltiorrhiza, and add to the knowledge for SAB biosynthesis in S. miltiorrhiza.

The TCP family genes are plant-specific transcription factors and play important roles in plant development. TCPs have been evolutionarily and functionally studied in several plants. Although common wheat (Triticum aestivum L.) is a major staple crop worldwide, no systematic analysis of TCPs in this important crop has been conducted. Here, we performed a genome-wide survey in wheat and found 66 TCP genes that belonged to 22 homoeologous groups. We then mapped these genes on wheat chromosomes and found that several TCP genes were duplicated in wheat including the ortholog of the maize TEOSINTE BRANCHED 1. Expression study using both RT-PCR and in situ hybridization assay showed that most wheat TCP genes were expressed throughout development of young spike and immature seed. Cis-acting element survey along promoter regions suggests that subfunctionalization may have occurred for homoeologous genes. Moreover, protein-protein interaction experiments of three TCP proteins showed that they can form either homodimers or heterodimers. Finally, we characterized two TaTCP9 mutants from tetraploid wheat. Each of these two mutant lines contained a premature stop codon in the A subgenome homoeolog that was dominantly expressed over the B subgenome homoeolog. We observed that mutation caused increased spike and grain lengths. Together, our analysis of the wheat TCP gene family provides a start point for further functional study of these important transcription factors in wheat.

In wheat breeding, improved quality traits, including grain quality and dough rheological properties, have long been a critical goal. To understand the genetic basis of key quality traits of wheat, two single-locus and five multi-locus GWAS models were performed for six grain quality traits and three dough rheological properties based on 19, 254 SNPs in 267 bread wheat accessions. As a result, 299 quantitative trait nucleotides (QTNs) within 105 regions were identified to be associated with these quality traits in four environments. Of which, 40 core QTN regions were stably detected in at least three environments, 19 of which were novel. Compared with the previous studies, these novel QTN regions explained smaller phenotypic variation, which verified the advantages of the multi-locus GWAS models in detecting important small effect QTNs associated with complex traits. After characterization of the function and expression in-depth, 67 core candidate genes involved in protein/sugar synthesis, histone modification and the regulation of transcription factor were observed to be associated with the formation of grain quality, which showed that multi-level regulations influenced wheat grain quality. Finally, a preliminary network of gene regulation that may affect wheat quality formation was inferred. This study verified the power and reliability of multi-locus GWAS methods in wheat quality trait research, and increased the understanding of wheat quality formation mechanisms. The detected QTN regions and candidate genes in this study could be further used for gene cloning and marker-assisted selection in high-quality breeding of bread wheat.