Identification and validation of novel low-tiller number QTL in common wheat

Genetic Mapping by 55K Single-Nucleotide Polymorphism Array Reveals Candidate Genes for Tillering Trait in Wheat Mutant dmc

BACKGROUND

The tiller number is a key agronomic trait for increasing the yield potential of wheat (Triticum aestivum L.). A number of quantitative trait loci (QTLs) and key genes controlling tillering have been identified, but the regulatory mechanisms remain unclear.

METHODS

In this study, we utilized the dwarf-monoculm mutant (dmc) obtained from the ethyl methane sulfonate (EMS)-treated wheat cultivar Guomai 301. The F2 populations were constructed using the dmc mutant crossed to multiple tiller parents. The F2 populations were surveyed for tillering traits at the critical fertility stage for genetic analyses. The extreme-tillering-phenotype plants from the F2 population were used to construct mixing pools that were analyzed by a wheat 55K SNP array. The tillering genes of dmc were mapped using the wheat 55K SNP array combined with transcriptomic data.

RESULTS

The results showed that the genetic phenotype of dmc is controlled by two dominant genes. The tillering genes of dmc were mapped on the 60-100 Mb region of chromosome 5B and the 135-160 Mb region of chromosome 7A. A total of sixteen candidate genes associated with the tillering trait of dmc were identified. Two candidate genes, TraesCS5B02G058800 and TraesCS7A02G184200, were predicted to be involved in indole acetic acid (IAA) response and transport, which were considered as potential regulatory genes.

CONCLUSIONS

This study elucidated the genetic basis of the dmc mutant and provided two valuable reference genes for studying the development and regulatory mechanisms of wheat tillering.

Identification and validation of novel plant compactness QTL in common wheat

BACKGROUND

Plant compactness (PC) is a crucial agronomic trait that affects plant density in wheat, which in turn influences biomass and grain yield potential. The canopy of high-yielding wheat varieties should exhibit appropriate aboveground plant architecture. In this study, three recombinant inbred line (RIL) populations were generated to identify and validate quantitative trait loci (QTL) related to plant compactness. The erect and compact tillering genotype SN05525 was used as a common parent. A total of 193 F8 RILs from the cross SN05525/SN22 were genotyped using the high-density Illumina iSelect 90 K single nucleotide polymorphism (SNP) assay.

RESULTS

A linkage map with 7180 SNP loci was constructed, revealing six QTL on chromosomes 3B (Qpc.sdau-3B.1, Qpc.sdau-3B.2),5B (Qpc.sdau-5B), 5D (Qpc.sdau-5D), 2 A (Qpc.sdau-2 A), and 7 A (Qpc.sdau-7 A) that control compact tillers. Qpc.sdau-3B.1, Qpc.sdau-5D, and Qpc.sdau-5B accounted for up to 16.70%, 16.89%, and 14.56% of the phenotypic variance, respectively. Kompetitive allele-specific PCR (KASP) markers were developed for these QTL, and their effects were validated in two additional RIL populations, SN05525/Luyuan 502 and SN05525/Xinong 511. Significant effects of Qpc.sdau-3B.1 and Qpc.sdau-5D on compactness were observed in the validation populations.

CONCLUSION

Three major QTL loci closely related to plant compactness in wheat were successfully identified, and their effects were validated in two additional RIL populations across multi environments. Plant architecture plays a crucial role in enhancing yield and economic value in wheat. In the process of molecular marker-assisted selection breeding, the closely linked KASP markers could potentially be utilized in molecular marker-assisted selection for adjusting plant compactness and for further characterization of the underlying gene(s).

Chromosome genomics facilitates the marker development and selection of wheat-Aegilops biuncialis addition, substitution and translocation lines

The annual goatgrass, Aegilops biuncialis is a rich source of genes with considerable agronomic value. This genetic potential can be exploited for wheat improvement through interspecific hybridization to increase stress resistance, grain quality and adaptability. However, the low throughput of cytogenetic selection hampers the development of alien introgressions. Using the sequence of flow-sorted chromosomes of diploid progenitors, the present study enabled the development of chromosome-specific markers. In total, 482 PCR markers were validated on wheat (Mv9kr1) and Ae. biuncialis (MvGB642) crossing partners, and 126 on wheat-Aegilops additions. Thirty-two markers specific for U- or M-chromosomes were used in combination with GISH and FISH for the screening of 44 Mv9kr1 × Ae. biuncialis BC3F3 genotypes. The predominance of chromosomes 4M and 5M, as well as the presence of chromosomal aberrations, may indicate that these chromosomes have a gametocidal effect. A new wheat-Ae. biuncialis disomic 4U addition, 4M(4D) and 5M(5D) substitutions, as well as several introgression lines were selected. Spike morphology and fertility indicated that the Aegilops 4M or 5M compensated well for the loss of 4D and 5D, respectively. The new cytogenetic stocks represent valuable genetic resources for the introgression of key genes alleles into wheat.

Genome-wide identification and characterization of circRNAs in wheat tiller

KEY MESSAGE

This study found that the intergenic circRNAs of wheat were more abundant than those of other plants. More importantly, a circRNA-mediated network associated with tillering was constructed for the first time. Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs with covalently closed circular structures, which play an important role in transcriptional and post-transcriptional regulation. Tiller is an important agronomic trait that determines plant morphological architecture and affects spike number in wheat. However, no studies on the characteristics and functions of circRNAs involved in the regulation of wheat tiller. Here, we performed a genome-wide identification of circRNAs using ribosomal-depleted RNA-seq from wheat tiller of two pairs near-isogenic lines. A total of 686 circRNAs were identified and distributed on 21 chromosomes of wheat, of which 537 were novel circRNAs. Unlike other plants, the majority of these circRNAs (61.8%) were derived from intergenic regions. One circRNA-mediated network associated with tillering was constructed through weighted gene co-expression network analysis, including 323 circRNAs, 117 miRNAs, and 968 mRNAs. GO and pathway enrichment analysis of mRNAs suggested that these circRNAs are involved in cell cycle, ncRNA export from nucleus, developmental process, plant hormone signal transduction, MAPK signaling pathway, RNA degradation. Of these circRNAs, ten circRNAs are associated with known tillering/branching genes in rice or Arabidopsis thaliana, including OsCesA7, EBR1, DTE1, CRD1, LPA1, PAY1, LRK1, OsNR2, OsCCA1, OsBZR1. In summary, we present the first study of the identification and characterization of circRNAs in wheat tiller, and the results suggest these circRNAs associated with tillering could play an important role in wheat tiller formation and development.

Tiller Number1 encodes an ankyrin repeat protein that controls tillering in bread wheat

Wheat (Triticum aestivum L.) is a major staple food for more than one-third of the world's population. Tiller number is an important agronomic trait in wheat, but only few related genes have been cloned. Here, we isolate a wheat mutant, tiller number1 (tn1), with much fewer tillers. We clone the TN1 gene via map-based cloning: TN1 encodes an ankyrin repeat protein with a transmembrane domain (ANK-TM). We show that a single amino acid substitution in the third conserved ankyrin repeat domain causes the decreased tiller number of tn1 mutant plants. Resequencing and haplotype analysis indicate that TN1 is conserved in wheat landraces and modern cultivars. Further, we reveal that the expression level of the abscisic acid (ABA) biosynthetic gene TaNCED3 and ABA content are significantly increased in the shoot base and tiller bud of the tn1 mutants; TN1 but not tn1 could inhibit the binding of TaPYL to TaPP2C via direct interaction with TaPYL. Taken together, we clone a key wheat tiller number regulatory gene TN1, which promotes tiller bud outgrowth probably through inhibiting ABA biosynthesis and signaling.

Precision phenotyping across the life cycle to validate and decipher drought-adaptive QTLs of wild emmer wheat (Triticum turgidum ssp. dicoccoides) introduced into elite wheat varieties

Drought events or the combination of drought and heat conditions are expected to become more frequent due to global warming, and wheat yields may fall below their long-term average. One way to increase climate-resilience of modern high-yielding varieties is by their genetic improvement with beneficial alleles from crop wild relatives. In the present study, the effect of two beneficial QTLs introgressed from wild emmer wheat and incorporated in the three wheat varieties BarNir, Zahir and Uzan was studied under well-watered conditions and under drought stress using non-destructive High-throughput Phenotyping (HTP) throughout the life cycle in a single pot-experiment. Plants were daily imaged with RGB top and side view cameras and watered automatically. Further, at two time points, the quantum yield of photosystem II was measured with a top view FluorCam. The QTL carrying near isogenic lines (NILs) were compared with their corresponding parents by t-test for all non-invasively obtained traits and for the manually determined agronomic and yield parameters. Data quality of phenotypic traits (repeatability) in the controlled HTP experiment was above 85% throughout the life cycle and at maturity. Drought stress had a strong effect on growth in all wheat genotypes causing biomass reduction from 2% up to 70% at early and late points in the drought period, respectively. At maturity, the drought caused 47-55% decreases in yield-related traits grain weight, straw weight and total biomass and reduced TKW by 10%, while water use efficiency (WUE) increased under drought by 29%. The yield-enhancing effect of the introgressed QTLs under drought conditions that were previously demonstrated under field/screenhouse conditions in Israel, could be mostly confirmed in a greenhouse pot experiment using HTP. Daily precision phenotyping enabled to decipher the mode of action of the QTLs in the different genetic backgrounds throughout the entire wheat life cycle. Daily phenotyping allowed a precise determination of the timing and size of the QTLs effect (s) and further yielded information about which image-derived traits are informative at which developmental stage of wheat during the entire life cycle. Maximum height and estimated biovolume were reached about a week after heading, so experiments that only aim at exploring these traits would not need a longer observation period. To obtain information on different onset and progress of senescence, the CVa curves represented best the ongoing senescence of plants. The QTL on 7A in the BarNir background was found to improve yield under drought by increased biomass growth, a higher photosynthetic performance, a higher WUE and a "stay green effect."

Elucidation of gene action and combining ability for productive tillering in spring barley

The purpose of the present study is to identify breeding and genetic peculiarities for productive tillering in spring barley genotypes of different origin, purposes of usage and botanical affiliation, as well as to identify effective genetic sources to further improving of the trait. There were created two complete (6 × 6) diallel crossing schemes. Into the Scheme I elite Ukrainian (MIP Tytul and Avhur) and Western European (Datcha, Quench, Gladys, and Beatrix) malting spring barley varieties were involved. Scheme II included awnless covered barley varieties Kozyr and Vitrazh bred at the Plant Production Institute named after V. Y. Yuriev of NAAS of Ukraine, naked barley varieties Condor and CDC Rattan from Canada, as well as awned feed barley variety MIP Myroslav created at MIW and malting barley variety Sebastian from Denmark. For more reliable and informative characterization of barley varieties and their progeny for productive tillering in terms of inheritance, parameters of genetic variation and general combining ability (GCA) statistical analyses of experimental data from different (2019 and 2020) growing seasons were conducted. Accordingly to the indicator of phenotypic dominance all possible modes of inheritance were detected, except for negative dominance in the Scheme I in 2020. The degree of phenotypic dominance significantly varied depending on both varieties involved in crossing schemes and conditions of the years of trials. There was overdominance in loci in both schemes in both years. The other parameters of genetic variation showed significant differences in gene action for productive tillering between crossing Schemes. In Scheme I in both years the dominance was mainly unidirectional and due to dominant effects. In the Scheme II in both years there was multidirectional dominance. In Scheme I compliance with the additive-dominant system was revealed in 2019, but in 2020 there was a strong epistasis. In Scheme II in both years non-allelic interaction was identified. In general, the mode of gene action showed a very complex gene action for productive tillering in barley and a significant role of non-genetic factors in phenotypic manifestation of the trait. Despite this, the level of heritability in the narrow sense in both Schemes pointed to the possibility of the successful selection of individuals with genetically determined increased productive tillering in the splitting generations. In Scheme I the final selection for productive tillering will be more effective in later generations, when dominant alleles become homozygous. In Scheme II it is theoretically possible to select plants with high productive tillering on both recessive and dominant basis. In both schemes the non-allelic interaction should be taken into consideration. Spring barley varieties Beatrix, Datcha, MIP Myroslav and Kozyr can be used as effective genetic sources for involvement in crossings aimed at improving the productive tillering. The results of present study contribute to further development of studies devoted to evaluation of gene action for yield-related traits in spring barley, as well as identification of new genetic sources for plant improvement.

Identification of One Major QTL and a Novel Gene OsIAA17q5 Associated with Tiller Number in Rice Using QTL Analysis

Rice tillers are one of the most important traits for the yield and development of rice, although little is known about its mode of inheritance. Tiller numbers were recorded every 7 days a total of nine times, starting 30 days after transplantation. Quantitative trait locus (QTL) based analysis on a set of double haploid population derivatives of a cross between the Cheongcheong and Nagdong varieties identified a major effect of locus RM18130-RM3381 on chromosome 5, which was expressed in eight different growth stages. Within the target region RM18130-RM3381 (physical distance: 2.08 Mb), 61 candidate genes were screened by annotation. Among the candidate genes, Os05g0230700 (named OsIAA17q5), which belongs to the family of auxin-responsive genes, was selected as a target. Auxin promotes cell division and meristem maintenance and is an effective plant regulator which influences plant growth and development by altering the expression of various genes. OsIAA17q5 is expected to control the number of tillers. The present study provides further understanding of the basic genetic mechanisms that selectively express the control of tiller numbers in different growth stages, as well as provides valuable information for future research aimed at cloning the target gene. These results may contribute to developing a comprehensive understanding of the basic genetic processes regulating the developmental behavior of tiller numbers in rice.

Fine mapping of the tiller inhibition gene TIN4 contributing to ideal plant architecture in common wheat

A tiller inhibition gene, TIN4, was mapped to an approximately 311 kb genomic interval on chromosome arm 2DL of wheat. The tiller is one of the key components of plant morphological architecture and a central agronomic trait affecting spike number in wheat. Low tiller number has been proposed as a major component of crop ideotypes for high yield potential. In this study, we characterized the development of tillering in near-isogenic lines (NIL7A and NIL7B), indicating that the TIN4 gene inhibited the growth of tillering buds and negatively regulated tiller number. Low-tillering was controlled by a single gene (TIN4) located on chromosome 2DL by genetic analysis and bulked segregant RNA-seq analysis. A total of 17 new polymorphic markers were developed in this study, and 61 recombinants were identified in the secondary F2 population containing 4,266 individuals. TIN4 was finally mapped on a 0.35 cM interval, co-segregated with molecular marker M380, within a 311 kb genomic interval of the wheat cultivar Chinese Spring reference genome sequence that contained twelve predicted genes. Yield experiments showed that the yield of low-tillering lines was higher than that of high-tillering lines at a higher density. Overall, this study provides a foundation for the construction of a low-tillering ideotype for improving wheat yield and further cloning TIN4 by map-based cloning approach.

Identification and Validation of Quantitative Trait Loci Mapping for Spike-Layer Uniformity in Wheat

Spike-layer uniformity (SLU), the consistency of the spike distribution in the vertical space, is an important trait. It directly affects the yield potential and appearance. Revealing the genetic basis of SLU will provide new insights into wheat improvement. To map the SLU-related quantitative trait loci (QTL), 300 recombinant inbred lines (RILs) that were derived from a cross between H461 and Chinese Spring were used in this study. The RILs and parents were tested in fields from two continuous years from two different pilots. Phenotypic analysis showed that H461 was more consistent in the vertical spatial distribution of the spike layer than in Chinese Spring. Based on inclusive composite interval mapping, four QTL were identified for SLU. There were two major QTL on chromosomes 2BL and 2DL and two minor QTL on chromosomes 1BS and 2BL that were identified. The additive effects of QSlu.sicau-1B, Qslu.sicau-2B-2, and QSlu.sicau-2D were all from the parent, H461. The major QTL, QSlu.sicau-2B-2 and QSlu.sicau-2D, were detected in each of the conducted trials. Based on the best linear unbiased prediction values, the two loci explained 23.97% and 15.98% of the phenotypic variation, respectively. Compared with previous studies, the two major loci were potentially novel and the two minor loci were overlapped. Based on the kompetitive allele-specific PCR (KASP) marker, the genetic effects for QSlu.sicau-2B-2 were validated in an additional RIL population. The genetic effects ranged from 26.65% to 32.56%, with an average value of 30.40%. In addition, QSlu.sicau-2B-2 showed a significant (p < 0.01) and positive influence on the spike length, spikelet number, and thousand kernel weight. The identified QTL and the developed KASP marker will be helpful for fine-mapping these loci, finally contributing to wheat breeding programs in a marker-assisted selection way.

Genome wide association study of plant height and tiller number in hulless barley

Hulless barley (Hordeum vulgare L. var. nudum), also called naked barley, is a unique variety of cultivated barley. The genome-wide specific length amplified fragment sequencing (SLAF-seq) method is a rapid deep sequencing technology that is used for the selection and identification of genetic loci or markers. In this study, we collected 300 hulless barley accessions and used the SLAF-seq method to identify candidate genes involved in plant height (PH) and tiller number (TN). We obtained a total of 1407 M paired-end reads, and 228,227 SLAF tags were developed. After filtering using an integrity threshold of >0.8 and a minor allele frequency of >0.05, 14,504,892 single-nucleotide polymorphisms (SNP) loci were screened out. The remaining SNPs were used for the construction of a neighbour-joining phylogenetic tree, and the three subcluster members showed no obvious differentiation among regional varieties. We used a genome wide association study approach to identify 1006 and 113 SNPs associated with TN and PH, respectively. Based on best linear unbiased predictors (BLUP), 41 and 29 SNPs associated with TN and PH, respectively. Thus, several of genes, including Hd3a and CKX5, may be useful candidates for the future genetic breeding of hulless barley. Taken together, our results provide insight into the molecular mechanisms controlling barley architecture, which is important for breeding and yield.

QTL Analysis of Five Morpho-Physiological Traits in Bread Wheat Using Two Mapping Populations Derived from Common Parents

Traits such as plant height (PH), juvenile growth habit (GH), heading date (HD), and tiller number are important for both increasing yield potential and improving crop adaptation to climate change. In the present study, these traits were investigated by using the same bi-parental population at early (F₂ and F₂-derived F₃ families) and late (F₆ and F₇, recombinant inbred lines, RILs) generations to detect quantitative trait loci (QTLs) and search for candidate genes. A total of 176 and 178 lines were genotyped by the wheat Illumina 25K Infinium SNP array. The two genetic maps spanned 2486.97 cM and 3732.84 cM in length, for the F₂ and RILs, respectively. QTLs explaining the highest phenotypic variation were found on chromosomes 2B, 2D, 5A, and 7D for HD and GH, whereas those for PH were found on chromosomes 4B and 4D. Several QTL detected in the early generations (i.e., PH and tiller number) were not detected in the late generations as they were due to dominance effects. Some of the identified QTLs co-mapped to well-known adaptive genes (i.e., Ppd-1, Vrn-1, and Rht-1). Other putative candidate genes were identified for each trait, of which PINE1 and PIF4 may be considered new for GH and TTN in wheat. The use of a large F₂ mapping population combined with NGS-based genotyping techniques could improve map resolution and allow closer QTL tagging.

Identification and validation of genetic loci for tiller angle in bread wheat

KEY MESSAGE

Two major QTL for tiller angle were identified on chromosomes 1AL and 5DL, and TaTAC-D1 might be the candidate gene for QTA.caas-5DL. An ideal plant architecture is important for achieving high grain yield in crops. Tiller angle (TA) is an important factor influencing yield. In the present study, 266 recombinant inbred lines (RILs) derived from a cross between Zhongmai 871 (ZM871) and its sister line Zhongmai 895 (ZM895) was used to map TA by extreme pool-genotyping and inclusive composite interval mapping (ICIM). Two quantitative trait loci (QTL) on chromosomes 1AL and 5DL were identified with reduced tiller angle alleles contributed by ZM895. QTA.caas-1AL was detected in six environments, explaining 5.4-11.2% of the phenotypic variances. The major stable QTL, QTA.caas-5DL, was identified in all eight environments, accounting for 13.8-24.8% of the phenotypic variances. The two QTL were further validated using BC1F4 populations derived from backcrosses ZM871/ZM895//ZM871 (121 lines) and ZM871/ZM895//ZM895 (175 lines). Gene TraesCS5D02G322600, located in the 5DL QTL and designated TaTAC-D1, had a SNP in the third exon with 'A' and 'G' in ZM871 and ZM895, respectively, resulting in a Thr169Ala amino acid change. A KASP marker based on this SNP was validated in two sets of germplasm, providing further evidence for the significant effects of TaTAC-D1 on TA. Thus extreme pool-genotyping can be employed to detect QTL for plant architecture traits and KASP markers tightly linked with the QTL can be used in wheat breeding programs targeting improved plant architecture.

Identification and Validation of a Novel Major Quantitative Trait Locus for Plant Height in Common Wheat (Triticum aestivum L.)

Plant height (PH) plays a pivotal role in plant morphological architecture and is associated with yield potential in wheat. For the quantitative trait locus (QTL) analysis, a recombinant inbred line population was developed between varieties differing significantly in PH. Two major QTL were identified on chromosomes 4B (QPh.sicau-4B) and 6D (QPh.sicau-6D) in multiple environments, which were then validated in two different backgrounds by using closely linked markers. QPh.sicau-4B explained 10.1-21.3% of the phenotypic variance, and the location corresponded to the dwarfing gene Rht-B1. QPh.sicau-6D might be a novel QTL for PH, explaining 6.6-13.6% of the phenotypic variance and affecting spike length, thousand-kernel weight, and spikelet compactness. Three candidate genes associated with plant growth and development were identified in the physical interval of QPh.sicau-6D. Collectively, we identified a novel stable and major PH QTL, QPh.sicau-6D, which could aid in the development of closely linked markers for marker-assisted breeding and cloning genes underlying this QTL.

Genome-Wide Association Study Reveals Genomic Regions Associated With Ten Agronomical Traits in Wheat Under Late-Sown Conditions

Poor understanding of the genetic and molecular basis of heat tolerance component traits is a major bottleneck in designing heat tolerant wheat cultivars. The impact of terminal heat stress is generally reported in the case of late sown wheat. In this study, our aim was to identify genomic regions for various agronomic traits under late sown conditions by using genome-wide association approach. An association mapping panel of 205 wheat accessions was evaluated under late sown conditions at three different locations in India. Genotyping of the association panel revealed 15,886 SNPs, out of which 11,911 SNPs with exact physical locations on the wheat reference genome were used in association analysis. A total of 69 QTLs (10 significantly associated and 59 suggestive) were identified for ten different traits including productive tiller number (17), grain yield (14), plant height (12), grain filling rate (6), grain filling duration (5), days to physiological maturity (4), grain number (3), thousand grain weight (3), harvest index (3), and biomass (2). Out of these associated QTLs, 17 were novel for traits, namely PTL (3), GY (2), GFR (6), HI (3) and GNM (3). Moreover, five consistent QTLs across environments were identified for GY (4) and TGW (1). Also, 11 multi-trait SNPs and three hot spot regions on Chr1Ds, Chr2BS, Chr2DS harboring many QTLs for many traits were identified. In addition, identification of heat tolerant germplasm lines based on favorable alleles HD2888, IC611071, IC611273, IC75240, IC321906, IC416188, and J31-170 would facilitate their targeted introgression into popular wheat cultivars. The significantly associated QTLs identified in the present study can be further validated to identify robust markers for utilization in marker-assisted selection (MAS) for development of heat tolerant wheat cultivars.

Detection of genomic regions associated with tiller number in Iranian bread wheat under different water regimes using genome-wide association study

Two of the important traits for wheat yield are tiller and fertile tiller number, both of which have been thought to increase cereal yield in favorable and unfavorable environments. A total of 6,349 single nucleotide polymorphism (SNP) markers from the 15 K wheat Infinium array were employed for genome-wide association study (GWAS) of tillering number traits, generating a physical distance of 14,041.6 Mb based on the IWGSC wheat genome sequence. GWAS analysis using Fixed and random model Circulating Probability Unification (FarmCPU) identified a total of 47 significant marker-trait associations (MTAs) for total tiller number (TTN) and fertile tiller number (FTN) in Iranian bread wheat under different water regimes. After applying a 5% false discovery rate (FDR) threshold, a total of 13 and 11 MTAs distributed on 10 chromosomes were found to be significantly associated with TTN and FTN, respectively. Linked single nucleotide polymorphisms for IWB39005 (2A) and IWB44377 (7A) were highly significantly associated (FDR < 0.01) with TTN and FTN traits. Moreover, to validate GWAS results, meta-analysis was performed and 30 meta-QTL regions were identified on 11 chromosomes. The integration of GWAS and meta-QTLs revealed that tillering trait in wheat is a complex trait which is conditioned by the combined effects of minor changes in multiple genes. The information provided by this study can enrich the currently available candidate genes and genetic resources pools, offering evidence for subsequent analysis of genetic adaptation of wheat to different climatic conditions of Iran and other countries.

Identification and validation of stable quantitative trait loci for grain filling rate in common wheat (Triticum aestivum L.)

We identified and validated two stable grain filling rate (GFR) quantitative trait loci (QTL) in wheat that positively influenced several yield-related traits. Among them, QGfr.sicau-7D.1 was a novel GFR QTL. The grain filling rate (GFR) plays a crucial role in determining grain yield. To advance the current understanding of the genetic characteristics underlying the GFR in common wheat, three recombinant inbred line populations were used to map and validate GFR quantitative trait loci (QTL). Using a high-density genetic linkage map, 10 GFR QTL were detected. They were located on chromosomes 2D, 4A, 4B, 5B, 6D, 7A and 7D, explained 4.99-12.62% of the phenotypic variation. Two of them, QGfr.sicau-6D and QGfr.sicau-7D.1, were detected in all four environments tested and their genetic effect was validated by closely linked kompetitive allele specific PCR (KASP) markers in different genetic backgrounds. The effects of these two GFR QTL on other yield-related traits were also estimated. QGfr.sicau-6D had a significant positive influence (p < 0.01) on thousand kernel weight, kernel width, kernel volume, and kernel surface area. QGfr.sicau-7D.1 had a significant positive influence (p < 0.01) on thousand kernel weight and kernel length. Furthermore, QGfr.sicau-7D.1 was a completely novel QTL for GFR; several genes associated with grain growth and development were predicted in its physical interval. These results will facilitate molecular marker-assisted selection of wheat with high-confidence QTL for GFR and fine mapping of genes associated with GFR, thereby contributing to yield improvement.

Characterization of a major quantitative trait locus on the short arm of chromosome 4B for spike number per unit area in common wheat (Triticum aestivum L.)

An InDel marker closely linked with a major and stable quantitative trait locus (QTL) on chromosome 4BS, QSnpa.cau-4B, controlling spike number per unit area will benefit wheat yield improvement. Spike number per unit area (SNPA) is an essential yield-related trait, and analyzing its genetic basis is important for cultivar improvement in wheat (Triticum aestivum L.). In this study, we used the F₂ population derived from a cross between two wheat accessions displaying significant differences in SNPA to perform quantitative trait locus (QTL) analysis. Through bulked segregant analysis, a major and stable QTL that explained 18.11-82.11% of the phenotypic variation was identified on chromosome 4BS. The QTL interval was validated using F_4:5 and F_6:7 families and narrowed it to a 24.91-38.36 Mb region of chromosome 4BS according to the 'Chinese Spring' reference genome sequence. In this region, variations in 16 genes caused amino acid changes and three genes were present in only one parent. Among these, we annotated a gene orthologous to TB1 in maize (Zea mays), namely TraesCS4B01G042700, which carried a 44-bp deletion in its promoter in the higher-SNPA parent. An InDel marker based on the insertion/deletion polymorphism was designed and used to diagnose the allelic distribution within a natural population. The frequency of the 44-bp deletion allele associated with higher SNPA was relatively low (13.24%), implying that this favorable allele has not been widely utilized and could be valuable for wheat yield improvement. In summary, we identified a major and stable QTL for SNPA and developed a diagnostic marker for the more-spiked trait, which will be beneficial for molecular-assisted breeding in wheat.

A Genome-Wide Association Study of Coleoptile Length in Different Chinese Wheat Landraces

From the perspective of wheat yield improvement, the coleoptile is vital for successful crop establishment, and long coleoptile lengths (CLs) are preferred in wheat-growing regions where deep planting is practiced. To determine the genetic basis of CL, we performed a genome-wide association study on a set of 707 Chinese wheat landraces using 18,594 single-nucleotide polymorphisms and 38,678 diversity array technology sequencing markers. We accordingly detected a total of 29 significant markers [-log₁₀ ^(P) > 4.76] distributed on chromosomes 2B, 2D, 3A, 4A, 5A, 6A, 6B, 6D, and 7B. Based on linkage disequilibrium decay distance, we identified a total of 17 quantitative trait loci associated with CL, among which QCl.sicau-6B.2, located at 508.17-509.26 Mb on chromosome 6B, was recognized as a novel major locus. We subsequently developed a high-resolution melt marker for QCl.sicau-6B.2, which was validated in an F _{2 : 3} population. Our findings provide important insights into the genetic mechanisms underlying coleoptile growth and could be applied to marker-assisted wheat selection.

Identification of lncRNAs involved in wheat tillering development in two pairs of near-isogenic lines

Emerging evidence demonstrates that lncRNAs participate in various developmental processes in plants via post-transcription regulation. However, few lncRNAs have been identified as regulators of tiller development in wheat (Triticum aestivum L.). In this study, high-throughput ribosomal depleted RNA sequencing was performed on the tillering nodes of two pairs of near-isogenic lines that differed only in the tillering trait. We identified 5399 lncRNA transcripts using bioinformational analyses. KEGG pathway analysis revealed 74 common differentially expressed lncRNAs substantially enriched in photosynthesis-related, phenylpropanoid biosynthesis, phosphatidylinositol signaling, brassinosteroid biosynthesis, zeatin biosynthesis, and carotenoid biosynthesis pathways. Detailed functional annotations of target genes were used to identify 27 tillering-associated lncRNAs. Among these, 10 were in photosynthesis-related pathways; 15 were in secondary metabolite pathways; and 8 were in plant hormone pathways, with 6 enriched in two kinds of pathways. These findings contribute to identifying tillering-associated lncRNAs in wheat and enable further investigation into the functions and roles of key candidate lncRNAs, and more experimental evidence was also needed if breeders wanted to utilize these candidate lncRNAs in wheat crop yield improvement in the future.

Functional mutation allele mining of plant architecture and yield-related agronomic traits and characterization of their effects in wheat

BACKGROUND

Wheat mutant resources with phenotypic variation have been developed in recent years. These mutants might carry favorable mutation alleles, which have the potential to be utilized in the breeding process. Plant architecture and yield-related features are important agronomic traits for wheat breeders and mining favorable alleles of these traits will improve wheat characteristics.

RESULTS

Here we used 190 wheat phenotypic mutants as material and by analyzing their SNP variation and phenotypic data, mutation alleles for plant architecture and yield-related traits were identified, and the genetic effects of these alleles were evaluated. In total, 32 mutation alleles, including three pleiotropic alleles, significantly associated with agronomic traits were identified from the 190 wheat mutant lines. The SNPs were distributed on 12 chromosomes and were associated with plant height (PH), tiller number, flag leaf angle (FLA), thousand grain weight (TGW), and other yield-related traits. Further phenotypic analysis of multiple lines carrying the same mutant allele was performed to determine the effect of the allele on the traits of interest. PH-associated SNPs on chromosomes 2BL, 3BS, 3DL, and 5DL might show additive effects, reducing PH by 10.0 cm to 31.3 cm compared with wild type, which means that these alleles may be favorable for wheat improvement. Only unfavorable mutation alleles that reduced TGW and tiller number were identified. A region on chromosome 5DL with mutation alleles for PH and TGW contained several long ncRNAs, and their sequences shared more than 90% identity with cytokinin oxidase/dehydrogenase genes. Some of the mutation alleles we mined were colocalized with previously reported QTLs or genes while others were novel; these novel alleles could also result in phenotypic variation.

CONCLUSION

Our results demonstrate that favorable mutation alleles are present in mutant resources, and the region between 409.5 to 419.8 Mb on chromosome 5DL affects wheat plant height and thousand grain weight.

Agronomic and Physiological Traits, and Associated Quantitative Trait Loci (QTL) Affecting Yield Response in Wheat (Triticum aestivum L.): A Review

Enhanced grain yield has been achieved in bread wheat (Triticum aestivum L.) through development and cultivation of superior genotypes incorporating yield-related agronomic and physiological traits derived from genetically diverse and complementary genetic pool. Despite significant breeding progress, yield levels in wheat have remained relatively low and stagnant under marginal growing environments. There is a need for genetic improvement of wheat using yield-promoting morpho-physiological attributes and desired genotypes under the target production environments to meet the demand for food and feed. This review presents breeding progress in wheat for yield gains using agronomic and physiological traits. Further, the paper discusses globally available wheat genetic resources to identify and select promising genotypes possessing useful agronomic and physiological traits to enhance water, nutrient-, and radiation-use efficiency to improve grain yield potential and tolerance to abiotic stresses (i.e. elevated CO2, high temperature, and drought stresses). Finally, the paper highlights quantitative trait loci (QTL) linked to agronomic and physiological traits to aid breeding of high-performing wheat genotypes.

Comprehensive transcriptomics, proteomics, and metabolomics analyses of the mechanisms regulating tiller production in low-tillering wheat

Tiller development in low-tillering wheat is related to several differentially expressed genes, proteins, and metabolites, as determined by an integrated omics approach combining transcriptome analysis, iTRAQ, and HPLC-MS on multiple NILs. Tillering is an important aspect of plant morphology that affects spike number, thereby contributing to the final crop yield. However, the mechanisms inhibiting tiller production in low-tillering wheat are poorly characterized. To investigate this aspect of wheat biology, two pairs of near-isogenic lines were developed, and an integrated omics approach combining transcriptome analysis, isobaric tags for relative and absolute quantification, and high-performance liquid chromatography-mass spectrometry were used to compare the free-tillering and low-tillering caused by an allele at Qltn.sicau-2D in wheat samples. Overall, 474 genes, 166 proteins, and 28 metabolites were identified as tillering-associated differentially expressed genes, proteins, and metabolites (DEGs, DEPs, and DEMs, respectively). Functional analysis indicated that the abundance of DEGs/DEPs/DEMs was related to lignin and cellulose metabolism, cell division, cell cycle processes, and glycerophospholipid metabolism; three transcription factor families, GRAS, GRF, and REV, might be related to the decrease in tillering in low-tillering wheat. These findings contribute to improve our understanding of the mechanisms responsible for the inhibition of tiller development in low-tillering wheat cultivars.

Optimizing plant density and nitrogen application to manipulate tiller growth and increase grain yield and nitrogen-use efficiency in winter wheat

The growth of wheat tillers and plant nitrogen-use efficiency (NUE) will gradually deteriorate in response to high plant density and over-application of N. Therefore, in this study, a 2-year field study was conducted with three levels of plant densities (75 ×10⁴plants ha⁻¹, D1; 300 ×10⁴plants ha⁻¹, D2; 525 ×10⁴plants ha⁻¹, D3) and three levels of N application rates (120 kg N ha⁻¹, N1; 240 kg N ha⁻¹, N2; 360 kg N ha⁻¹, N3) to determine how to optimize plant density and N application to regulate tiller growth and to assess the contribution of such measures to enhancing grain yield (GY) and NUE. The results indicated that an increase in plant density significantly increased the number of superior tillers and the number of spikes per m²(SN), resulting in a higher GY and higher partial factor productivity of applied N (PFP_N). However, there was no significant difference in GY and PFP_N between plant densities D2 and D3. Increasing the N application rate significantly increased the vascular bundle number (NVB) and area (AVB), however, excess N application (N3) did not significantly improve these parameters. N application significantly increased GY, whereas there was a significant decrease in PFP_N in response to an increase in N application rate. The two years results suggested that increasing the plant density (from 75 ×10⁴plants ha⁻¹to 336 ×10⁴plants ha⁻¹) in conjunction with the application of 290 kg N ha⁻¹N will maximize GY, and also increase PFP_N(39.7 kg kg⁻¹), compared with the application of 360 kg N ha⁻¹N. Therefore, an appropriate combination of increased planting density with reduced N application could regulate tiller number and favor the superior tiller group, to produce wheat populations with enhanced yield and NUE.

QTL identification and KASP marker development for productive tiller and fertile spikelet numbers in two high-yielding hard white spring wheat cultivars

Selecting high-yielding wheat cultivars with more productive tillers per unit area (PTN) combined with more fertile spikelets per spike (fSNS) is difficult. QTL mapping of these traits may aid understanding of this bottleneck and accelerate precision breeding for high yield via marker-assisted selection. PTN and fSNS were assessed in four to five trials from 2015 to 2017 in a doubled haploid population derived from two high-yielding cultivars "UI Platinum" and "SY Capstone." Two QTL for PTN (QPTN.uia-4A and QPTN.uia-6A) and four QTL for fSNS (QfSNS.uia-4A, QfSNS.uia-5A, QfSNS.uia-6A, and QfSNS.uia-7A) were identified. The effects of the QTL were primarily additive and, therefore, pyramiding of multiple QTL may increase PTN and fSNS. However, the two QTL for PTN were positioned in the flanking regions for the two QTL for fSNS on chromosomes 4A and 6A, respectively, suggesting either possible pleiotropic effect of the same QTL or tightly linked QTL and explaining the difficulty of selecting both high PTN and fSNS in phenotypic selection. Kompetitive allele-specific PCR (KASP) markers for all identified QTL were developed and validated in a recombinant inbred line (RIL) population derived from the same two cultivars. In addition, KASP markers for three of the QTL (QPTN.uia-6A, QfSNS.uia-6A, and QfSNS.uia-7A) were further validated in a diverse spring wheat panel, indicating their usefulness under different genetic backgrounds. These KASP markers could be used by wheat breeders to select high PTN and fSNS.

Quantitative Changes in the Transcription of Phytohormone-Related Genes: Some Transcription Factors Are Major Causes of the Wheat Mutant dmc Not Tillering

Tiller number is an important agronomic trait for grain yield of wheat (Triticum aestivum L.). A dwarf-monoculm wheat mutant (dmc) was obtained from cultivar Guomai 301 (wild type, WT). Here, we explored the molecular basis for the restrained tiller development of the mutant dmc. Two bulked samples of the mutant dmc (T1, T2 and T3) and WT (T4, T5 and T6) with three biological replicates were comparatively analyzed at the transcriptional level by bulked RNA sequencing (RNA-Seq). In total, 68.8 Gb data and 463 million reads were generated, 80% of which were mapped to the wheat reference genome of Chinese Spring. A total of 4904 differentially expressed genes (DEGs) were identified between the mutant dmc and WT. DEGs and their related major biological functions were characterized based on GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) categories. These results were confirmed by quantitatively analyzing the expression profiles of twelve selected DEGs via real-time qRT-PCR. The down-regulated gene expressions related to phytohormone syntheses of auxin, zeatin, cytokinin and some transcription factor (TF) families of TALE, and WOX might be the major causes of the mutant dmc, not tillering. Our work provides a foundation for subsequent tiller development research in the future.

Global QTL Analysis Identifies Genomic Regions on Chromosomes 4A and 4B Harboring Stable Loci for Yield-Related Traits Across Different Environments in Wheat (Triticum aestivum L.)

Major advances in wheat production are needed to address global food insecurity under future climate conditions, such as high temperatures. The grain yield of bread wheat (Triticum aestivum L.) is a quantitatively inherited complex trait that is strongly influenced by interacting genetic and environmental factors. Here, we conducted global QTL analysis for five yield-related traits, including spike yield, yield components and plant height (PH), in the Nongda3338/Jingdong6 doubled haploid (DH) population using a high-density SNP and SSR-based genetic map. A total of 12 major genomic regions with stable QTL controlling yield-related traits were detected on chromosomes 1B, 2A, 2B, 2D, 3A, 4A, 4B, 4D, 5A, 6A, and 7A across 12 different field trials with timely sown (normal) and late sown (heat stress) conditions. Co-location of yield components revealed significant tradeoffs between thousand grain weight (TGW) and grain number per spike (GNS) on chromosome 4A. Dissection of a "QTL-hotspot" region for grain weight on chromosome 4B was helpful in marker-assisted selection (MAS) breeding. Moreover, this study identified a novel QTL for heat susceptibility index of thousand grain weight (HSITGW) on chromosome 4BL that explains approximately 10% of phenotypic variation. QPh.cau-4B.2, QPh.cau-4D.1 and QPh.cau-2D.3 were coincident with the dwarfing genes Rht1, Rht2, and Rht8, and haplotype analysis revealed their pleiotropic architecture with yield components. Overall, our findings will be useful for elucidating the genetic architecture of yield-related traits and developing new wheat varieties with high and stable yield.

Utilization of a Wheat55K SNP Array for Mapping of Major QTL for Temporal Expression of the Tiller Number

Maximum tiller number and productive tiller number are important traits for wheat grain yield, but research involving the temporal expression of tiller number at different quantitative trait loci (QTL) levels is limited. In the present study, a population set of 371 recombined inbred lines derived from a cross between Chuan-Nong18 and T1208 was used to construct a high-density genetic map using a Wheat55K SNP Array and to perform dynamic QTL analysis of the tiller number at four growth stages. A high-density genetic map containing 11,583 SNP markers and 59 SSR markers that spanned 4,513.95 cM and was distributed across 21 wheat chromosomes was constructed. A total of 28 single environmental QTL were identified in the recombined inbred lines population, and among these, seven QTL were stable and used for multi-environmental and dynamic analysis. These QTL were mapped to chromosomes 2D, 4A, 4D, 5A, 5D, and 7D, respectively. Each QTL explained 1.63-21.22% of the observed phenotypic variation, with an additive effect from -20.51 to 11.59. Dynamic analysis showed that cqTN-2D.2 can be detected at four growth stages of tillering, explaining 4.92-17.16% of the observed phenotypic variations and spanning 13.71 Mb (AX-109283238-AX-110544009: 82189047-95895626) according to the physical location of the flanking markers. The effects of the stable QTL were validated in the recombined inbred lines population, and the beneficial alleles could be utilized in future marker-assisted selection. Several candidate genes for MTN and PTN were predicted. The results provide a better understanding of the QTL selectively expressing the control of tiller number and will facilitate future map-based cloning. 9.17% SNP markers showed best hits to the Chinese Spring contigs. It was indicated that Wheat55K Array was efficient and valid to construct a high-density wheat genetic map.

QTL identification and KASP marker development for productive tiller and fertile spikelet numbers in two high-yielding hard white spring wheat cultivars

Selecting high-yielding wheat cultivars with more productive tillers per unit area (PTN) combined with more fertile spikelets per spike (fSNS) is difficult. QTL mapping of these traits may aid understanding of this bottleneck and accelerate precision breeding for high yield via marker-assisted selection. PTN and fSNS were assessed in four to five trials from 2015 to 2017 in a doubled haploid population derived from two high-yielding cultivars "UI Platinum" and "SY Capstone." Two QTL for PTN (QPTN.uia-4A and QPTN.uia-6A) and four QTL for fSNS (QfSNS.uia-4A, QfSNS.uia-5A, QfSNS.uia-6A, and QfSNS.uia-7A) were identified. The effects of the QTL were primarily additive and, therefore, pyramiding of multiple QTL may increase PTN and fSNS. However, the two QTL for PTN were positioned in the flanking regions for the two QTL for fSNS on chromosomes 4A and 6A, respectively, suggesting either possible pleiotropic effect of the same QTL or tightly linked QTL and explaining the difficulty of selecting both high PTN and fSNS in phenotypic selection. Kompetitive allele-specific PCR (KASP) markers for all identified QTL were developed and validated in a recombinant inbred line (RIL) population derived from the same two cultivars. In addition, KASP markers for three of the QTL (QPTN.uia-6A, QfSNS.uia-6A, and QfSNS.uia-7A) were further validated in a diverse spring wheat panel, indicating their usefulness under different genetic backgrounds. These KASP markers could be used by wheat breeders to select high PTN and fSNS.

Exogenous Hormonal Application Regulates the Occurrence of Wheat Tillers by Changing Endogenous Hormones

Plant hormones play important roles in regulating the occurrence of crop tillers. However, little is known about the relationships and the underlying mechanisms between endogenous hormones and the occurrence of wheat tillers induced by exogenous hormones. In this study, two winter wheat cultivars, Xinong 979 and Xiaoyan 22, were used to investigate the effects of the exogenous application of indole-3-acetic acid (IAA) and zeatin (Z) on the occurrence of wheat tillers and investigate underlying mechanisms regulating the occurrence of tillers. The results showed that the application of IAA inhibited the occurrence of tillers, and external Z application promoted the occurrence rate of tillers under low nitrogen conditions. Further analysis of the results showed that exogenous IAA completely inhibited the growth of tiller buds, while exogenous Z significantly promoted the growth rate of tiller buds in low nitrogen conditions. Endogenous hormones exhibit important functions in regulating the growth of tiller buds, which contents were affected by exogenous hormones. Furthermore, according to the principal component analysis and correlation analysis, the growth of tiller buds was significantly positively correlated with the content of endogenous Z, whereas it was significantly negatively correlated with the ratios of endogenous IAA to endogenous Z (IAA:Z) and endogenous abscisic acid (ABA) to endogenous Z (ABA:Z). Moreover, no significant correlation was observed between the growth of the tiller buds and the endogenous IAA, endogenous gibberellins (GAs), and endogenous ABA content. These results suggested that Z played key roles in regulating the tiller occurrence, and exogenous hormones regulated the growth of wheat tiller buds via affecting the Z contents, thus regulating the occurrence of wheat tiller.

Molecular mapping and genetic analysis of a QTL controlling spike formation rate and tiller number in wheat

Spike formation rate (SR), which is based on maximum tiller number per unit area and spike number per unit area, is an important yield-related trait in wheat. Increasing the spike formation rate reduces growth competition and wastage of photosynthate from ineffective tillers. Unfortunately, research studies involving quantitative trait locus (QTL) mapping for wheat spike formation rate are limited. In the present study, a set of 371 recombinant inbreed line (RIL) population, which were derived from 1BL/LRS wheat-rye translocation lines CN18 and T1208, was analysed by simple sequence repeat (SSR) markers. Genetic analysis showed that a stable and major QTL (QSR.sicau-4D) for spike formation rate was localized to chromosome 4D and explained 18.24% and 24.48% of the observed phenotypic variance in 2015 and 2016, respectively. This QTL was closely linked to SSR marker Xcfd23, and the genetic distance between the flank markers was 3.28cM. Furthermore, QSR.sicau-4D might be a novel pleiotropic QTL, which also controlled maximum tiller number per unit area (QMTN.sicau-4D) and tiller number during pre-winter per unit area (QTNW.sicau-4D). The marker Xcfd23 associated with SR may be utilized in marker-assisted selection in wheat breeding.

Characterization of a common wheat (Triticum aestivum L.) high-tillering dwarf mutant

A novel high-tillering dwarf mutant in common wheat Wangshuibai was characterized and mapped to facilitate breeding for plant height and tiller and the future cloning of the causal gene. Tiller number and plant height are two major agronomic traits in cereal crops affecting plant architecture and grain yield. NAUH167, a mutant of common wheat landrace Wangshuibai induced by ethylmethyl sulfide (EMS) treatment, exhibits higher tiller number and reduced plant height. Microscope observation showed that the dwarf phenotype was attributed to the decrease in the number of cells and their length. The same as the wild type, the mutant was sensitive to exogenous gibberellins. Genetic analysis showed that the high-tillering number and dwarf phenotype were related and controlled by a partial recessive gene. Using a RIL_2:6 population derived from the cross NAUH167/Sumai3, a molecular marker-based genetic map was constructed. The map consisted of 283 loci, spanning a total length of 1007.98 cM with an average markers interval of 3.56 cM. By composite interval mapping, a stable major QTL designated QHt.nau-2D controlling both traits, was mapped to the short arm of chromosome 2D flanked by markers Xcfd11 and Xgpw361. To further map the QHt.nau-2D loci, another population consisted of 180 F₂ progeny from a cross 2011I-78/NAUH167 was constructed. Finally, QHt.nau-2D was located within a genetic region of 0.8 cM between markers QHT239 and QHT187 covering a predicted physical distance of 6.77 Mb. This research laid the foundation for map-based cloning of QHt.nau-2D and would facilitate the characterization of plant height and tiller number in wheat.