Natural variations in the promoter of Awn Length Inhibitor 1 (ALI-1) are associated with awn elongation and grain length in common wheat.

Wheat awn plays a vital role in photosynthesis, grain production, and drought tolerance. However, the systematic identification or cloning of genes controlling wheat awn development is seldom reported. Here, we conducted a genome-wide association study (GWAS) with 364 wheat accessions and identified 26 loci involved in awn length development, including previously characterized B1, B2, Hd, and several rice homologs. The dominant awn suppressor B1 was fine mapped to a 125-kb physical interval, and a C2 H2 zinc finger protein Awn Length Inhibitor 1 (ALI-1) was confirmed to be the underlying gene of the B1 locus through the functional complimentary test with native awnless allele. ALI-1 expresses predominantly in the developing spike of awnless individuals, transcriptionally suppressing downstream genes. ALI-1 reduces cytokinin content and simultaneously restrains cytokinin signal transduction, leading to a stagnation of cell proliferation and reduction of cell numbers during awn development. Polymorphisms of four single nucleotide polymorphisms (SNPs) located in ALI-1 promoter region are diagnostic for the B1/b1 genotypes, and these SNPs are associated with awn length (AL), grain length (GL) and thousand-grain weight (TGW). More importantly, ali-1 was observed to increase grain length in wheat, which is a valuable attribute of awn on grain weight, aside from photosynthesis. Therefore, ALI-1 pleiotropically regulates awn and grain development, providing an alternative for grain yield improvement and addressing future climate changes.

The MYB family transcription factor TuODORANT1 from Triticum urartu and the homolog TaODORANT1 from Triticum aestivum inhibit seed storage protein synthesis in wheat.

Seed storage proteins (SSPs) are determinants of wheat end-product quality. SSP synthesis is mainly regulated at the transcriptional level. Few transcriptional regulators of SSP synthesis have been identified in wheat and this study aims to identify novel SSP gene regulators. Here, the R2R3 MYB transcription factor TuODORANT1 from Triticum urartu was found to be preferentially expressed in the developing endosperm during grain filling. In common wheat (Triticum aestivum) overexpressing TuODORANT1, the transcription levels of all the SSP genes tested by RNA-Seq analysis were reduced by 49.71% throughout grain filling, which contributed to 13.38%-35.60% declines in the total SSP levels of mature grains. In in vitro assays, TuODORANT1 inhibited both the promoter activities and the transcription of SSP genes by 1- to 13-fold. The electrophoretic mobility shift assay (EMSA) and ChIP-qPCR analysis demonstrated that TuODORANT1 bound to the cis-elements 5'-T/CAACCA-3' and 5'-T/CAACT/AG-3' in SSP gene promoters both in vitro and in vivo. Similarly, the homolog TaODORANT1 in common wheat hindered both the promoter activities and the transcription of SSP genes by 1- to 112-fold in vitro. Knockdown of TaODORANT1 in common wheat led to 14.73%-232.78% increases in the transcription of the tested SSP genes, which contributed to 11.43%-19.35% elevation in the total SSP levels. Our data show that both TuODORANT1 and TaODORANT1 are repressors of SSP synthesis.

A novel NAC family transcription factor SPR suppresses seed storage protein synthesis in wheat.

The synthesis of seed storage protein (SSP) is mainly regulated at the transcriptional level. However, few transcriptional regulators of SSP synthesis have been characterized in common wheat (Triticum aestivum) owing to the complex genome. As the A genome donor of common wheat, Triticum urartu could be an elite model in wheat research considering its simple genome. Here, a novel NAC family transcription factor TuSPR from T. urartu was found preferentially expressed in developing endosperm during grain-filling stages. In common wheat transgenically overexpressing TuSPR, the content of total SSPs was reduced by c. 15.97% attributed to the transcription declines of SSP genes. Both in vitro and in vivo assays showed that TuSPR bound to the cis-element 5'-CANNTG-3' distributed in SSP gene promoters and suppressed the transcription. The homolog in common wheat TaSPR shared a conserved function with TuSPR on SSP synthesis suppression. The knock-down of TaSPR in common wheat resulted in 7.07%-20.34% increases in the total SSPs. Both TuSPR and TaSPR could be superior targets in genetic engineering to manipulate SSP content in wheat, and this work undoubtedly expands our knowledge of SSP gene regulation.

A SNP-based GWAS and functional haplotype-based GWAS of flag leaf-related traits and their influence on the yield of bread wheat (Triticum aestivum L.).

KEY MESSAGE: The genetic architecture of five flag leaf morphology traits was dissected by the functional haplotype-based GWAS and a standard SNP-based GWAS in a diverse population consisting of 197 varieties. Flag leaf morphology (FLM) is a critical factor affecting plant architecture and grain yield in wheat. The genetic architecture of FLM traits has been extensively studied with QTL mapping in bi-parental populations, while few studies exploited genome-wide association studies (GWAS) in diverse populations. In this study, a panel of 197 elite and historical varieties from China was evaluated for five FLM traits including the length (FLL), width (FLW), ratio (FLR), area (FLA) and angle (FLANG) as well as yield in nine environments. Based on the phenotypic correlation between yield and FLL (-0.43), FLA (- 0.32) and FLW (0.11), an empirical FLM index combining the three FLM traits proved to be a good predictor for yield. Two GWAS approaches were applied to dissect the genetic architecture of five FLM traits with a Wheat660K SNP array. The functional haplotype-based GWAS revealed 6, 5 and 7 QTL for FLANG, FLL and FLR, respectively, whereas two QTL for FLW and one for FLR were identified by the standard SNP-based GWAS. Due to co-localization, there were 18 independent QTL and 10 of them were close to known ones. One co-localized QTL on chromosome 5A was associated with FLL, FLANG and FLR. Moreover, both GWAS approaches identified a novel QTL for FLR on chromosome 6B which was not reported in previous studies. This study provides new insights into the relationship between FLM and yield and broadens our understanding of the genetic architecture of FLM traits in wheat.

Identification of quantitative trait loci for Fusarium head blight (FHB) resistance in the cross between wheat landrace N553 and elite cultivar Yangmai 13.

Fusarium head blight (FHB) of wheat poses a serious threat to food security in the Yellow-Huai River Valley Winter Wheat Region (YHW) of China. Discovery of new resistant quantitative trait loci (QTLs) or genes and application of them to highly susceptible varieties in the YHW are of great significance for ensuring the grain yield. Here, 160 recombinant inbred lines (RILs) from the cross between N553 (resistant) and Yangmai 13 (moderately susceptible) were used to evaluate FHB resistance by point inoculation, spray inoculation, and natural infection. A high-density genetic map was constructed by using a 15K SNP array and 128 polymorphism SSR markers. A total of 1452 polymorphic markers were identified, which formed 21 linkage groups and covered a total of 3555.1 cM in length. Two and four QTLs respectively related to type I and type II resistance were detected, among which QFhb-hnau.3BS.1 and QFhb-hnau.2DL were stably identified in most environments in Yangzhou and Zhengzhou, whereas QFhbn-hnau.5AL was only identified under natural infection in Jianyang. Based on the physical position (IWGSC RefSeq v1.0), QFhb-hnau.3BS.1 from the landrace N553 is likely to be Fhb1, while QFhb-hnau.2DL from Yangmai 13 may be a novel QTL. Significantly higher FHB resistance was observed in the lines with both QFhb-hnau.3BS.1 and QFhb-hnau.2DL, indicating that these two QTLs have apparent additive effects, and the RILs harboring both the two QTLs may have great application potential for the improvement of FHB resistance in wheat breeding. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01220-5.

Genetic Mapping of ms1s, a Recessive Gene for Male Sterility in Common Wheat.

The utilization of heterosis is an important way to improve wheat yield, and the production of wheat hybrid seeds mainly relies on male-sterile lines. Male sterility in line 15 Fan 03 derived from a cross of 72,180 and Xiaoyan 6 is controlled by a single recessive gene. The gene was mapped to the distal region of chromosome 4BS in a genetic interval of 1.4 cM and physical distance of 6.57 Mb between SSR markers Ms4BS42 and Ms4BS199 using an F2 population with 1205 individuals. Sterile individuals had a deletion of 4.57 Mb in the region presumed to carry the Ms1 locus. The allele for sterility was therefore named ms1s. Three CAPS markers were developed and verified from the region upstream of the deleted fragment and can be used for ms1s marker-assisted selection in wheat hybrid breeding. This work will enrich the utilization of male sterility genetic resources.

A wheat dominant dwarfing line with Rht12, which reduces stem cell length and affects gibberellic acid synthesis, is a 5AL terminal deletion line.

Dwarfing and semi-dwarfing are important agronomic traits that have great potential for the improvement of wheat yields. Rht12, a dominant gibberellic acid (GA)-responsive dwarfing gene from the gamma-ray-induced wheat mutant Karcagi 522M7K, is located in the long arm of chromosome 5A, which is closely linked with the locus Xwmc410. Rht12 is likely an ideal gene for GA biosynthesis and deactivation research in common wheat. However, information on the Rht12 locus and sequence is lacking. In this study, Rht12 significantly shortened stem cell length and decreased GA biosynthetic components. Using bulked segregant RNA-Seq, wheat 660k single nucleotide polymorphism chip detection, and newly developed simple sequence repeat markers, Rht12 was mapped to a 11.21-Mb region at the terminal end of chromosome 5AL, and was found to be closely linked with the Xw5ac207SSR marker with a 10.73-Mb fragment deletion in all of the homologous dwarfing plants. Transcriptome analyses of the remaining 483-kb region showed significantly higher expression of the TraesCS5A01G543100 gene encoding the GA metabolic enzyme GA 2-ß-dioxygenase in dwarfing plants than in high stalk plants, suggesting that Rht12 reduces plant height by activating TaGA2ox-A14. Taken together, our findings will promote cloning and functional studies of Rht12 in common wheat.

TubZIP28, a novel bZIP family transcription factor from Triticum urartu, and TabZIP28, its homologue from Triticum aestivum, enhance starch synthesis in wheat.

Starch in wheat grain provides humans with carbohydrates and influences the quality of wheaten food. However, no transcriptional regulator of starch synthesis has been identified first in common wheat (Triticum aestivum) due to the complex genome. Here, a novel basic leucine zipper (bZIP) family transcription factor TubZIP28 was found to be preferentially expressed in the endosperm throughout grain-filling stages in Triticum urartu, the A genome donor of common wheat. When TubZIP28 was overexpressed in common wheat, the total starch content increased by c. 4%, which contributed to c. 5% increase in the thousand kernel weight. The grain weight per plant of overexpression wheat was also elevated by c. 9%. Both in vitro and in vivo assays showed that TubZIP28 bound to the promoter of cytosolic AGPase and enhanced both the transcription and activity of the latter. Knockout of the homologue TabZIP28 in common wheat resulted in declines of both the transcription and activity of cytosolic AGPase in developing endosperms and c. 4% reduction of the total starch in mature grains. To the best of our knowledge, TubZIP28 and TabZIP28 are transcriptional activators of starch synthesis first identified in wheat, and they could be superior targets to improve the starch content and yield potential of wheat.

TaCKX gene family, at large, is associated with thousand-grain weight and plant height in common wheat.

KEY MESSAGE: We used SMRT sequencing and explored the haplotypes of TaCKX genes, linked with thousand-grain weight and plant height, and developed the functionally validated markers, which can be used in the marker-assisted breeding program. Cytokinin oxidase/dehydrogenase (CKX) enzymes catalyze the permanent degradation of cytokinins. Identification of the TaCKX alleles associated with yield traits and the development of functional markers is the first step in using these alleles in marker-assisted breeding program. To identify the alleles, we sequenced the genome fragments, containing TaCKX genes from 48 wheat genotypes, by PacBio® sequencing. Six out of 22 TaCKX genes were found polymorphic, forming 14 distinct haplotypes. Functional markers were developed and validated for all the polymorphic TaCKX genes. Four specific haplotypes, i.e., TaCKX2A_2, TaCKX4A_2, TaCKX5A_3, and TaCKX9A_2, were found significantly associated with high thousand-grain weight (TGW) and short plant height (PH) in Chinese wheat micro-core collection (MCC) and GWAS open population (GWAS-OP), whereas TaCKX1B_2 in GWAS-OP and TaCKX11A_3 in MCC were significantly associated with high TGW and short PH. The mean values of TGW and PH for cumulative favorable haplotypes from chromosome 3A, i.e., TaCKX2A_2, TaCKX4A_2, and TaCKX5A_3, were significantly higher as compared to the cumulative unfavored haplotypes, and the change was additive in manner. Frequency distribution analysis revealed that since the 1960s, the frequency of the favorable haplotypes and TGW has gradually increased in Chinese wheat cultivars. Expression profiling in the seed tissue excised at 2, 4, 6, and 8 days after anthesis depicted that the favorable haplotypes are significantly less expressive as compared to the unfavored haplotypes. We conclude that the functional markers developed in this study can be used to select the favorable haplotypes of TaCKX genes in wheat marker-assisted breeding programs.

Interactive effects of multiple vernalization (Vrn-1)- and photoperiod (Ppd-1)-related genes on the growth habit of bread wheat and their association with heading and flowering time.

BACKGROUND: The precise identification of Winterness/Springness (growth habit) for bread wheat, which is determined by genes involved in vernalization and photoperiod, will contribute to the effective utilization of bread wheat varieties. Here, 198 varieties from the Yellow and Huai wheat production region (YHW) in China were collected to identify their vernalization (Vrn-1) and photoperiod (Ppd-1) gene composition via a series of functional markers and their association with vernalization and photoperiod requirements at three locations during two years of experiments. The growth habits were measured during the spring sowing season. RESULTS: The results showed that the semi-winter varieties (grades1-4) were most prevalent in the population. The relative effects of single Vrn alleles on the growth period, such as heading date (HD) and/or flowering date (FD), were as follows: Vrn-B1b > Vrn-B1a > Vrn-D1b > Vrn-D1a > vrn-D1 = vrn-B1. The interactive effects of Vrn-B1 and Vrn-D1 on HD and FD were identical to those of Vrn-B1b. Approximately 35.3% of the cultivars carried Ppd-B1a (photoperiod-insensitive) and exhibited the earliest HD and FD. The Ppd-D1a-insensitive allele (Hapl II) was carried by just 0.5% of the varieties; however, the other two sensitive alleles were present at a higher frequency, and their effects were slightly weaker than those of Ppd-B1a. In addition, strong interactive effects between Ppd-B1 and Ppd-D1 were detected. In terms of mean values among various genotypes, the effects followed the order of Vrn-1 > Ppd-1. CONCLUSIONS: According to the results of ANOVA and least significant range (LSR) tests, we can conclude that Vrn-1 rather than Ppd-1 played a major role in controlling vernalization and photoperiod responses in this region. This research will be helpful for precisely characterizing and evaluating the HD, FD and even growth habit of varieties in the YHW at molecular levels.

Mechanisms, origin and heredity of Glu-1Ay silencing in wheat evolution and domestication.

KEY MESSAGE: Allotetraploidization drives Glu-1Ay silencing in polyploid wheat. The high-molecular-weight glutenin subunit gene, Glu-1Ay, is always silenced in common wheat via elusive mechanisms. To investigate its silencing and heredity during wheat polyploidization and domestication, the Glu-1Ay gene was characterized in 1246 accessions containing diploid and polyploid wheat worldwide. Eight expressed Glu-1Ay alleles (in 71.81% accessions) and five silenced alleles with a premature termination codon (PTC) were identified in Triticum urartu; 4 expressed alleles (in 41.21% accessions), 13 alleles with PTCs and 1 allele with a WIS 2-1A retrotransposon were present in wild tetraploid wheat; and only silenced alleles with PTC or WIS 2-1A were in cultivated tetra- and hexaploid wheat. Both the PTC number and position in T. urartu Glu-1Ay alleles (one in the N-terminal region) differed from its progeny wild tetraploid wheat (1-5 PTCs mainly in the repetitive domain). The WIS 2-1A insertion occurred ~ 0.13 million years ago in wild tetraploid wheat, much later than the allotetraploidization event. The Glu-1Ay alleles with PTCs or WIS 2-1A that arose in wild tetraploid wheat were fully succeeded to cultivated tetraploid and hexaploid wheat. In addition, the Glu-1Ay gene in wild einkorn inherited to cultivated einkorn. Our data demonstrated that the silencing of Glu-1Ay in tetraploid and hexaploid wheat was attributed to the new PTCs and WIS 2-1A insertion in wild tetraploid wheat, and most silenced alleles were delivered to the cultivated tetraploid and hexaploid wheat, providing a clear evolutionary history of the Glu-1Ay gene in the wheat polyploidization and domestication processes.

Genome-Wide Identification and Analysis of HAK/KUP/KT Potassium Transporters Gene Family in Wheat (Triticum aestivum L.).

In plants, the HAK (high-affinity K⁺)/KUP (K⁺ uptake)/KT (K⁺ transporter) family represents a large group of potassium transporters that play important roles in plant growth and environmental adaptation. Although HAK/KUP/KT genes have been extensively investigated in many plant species, they remain uncharacterized in wheat, especially those involved in the response to environmental stresses. In this study, 56 wheat HAK/KUP/KT (hereafter called TaHAKs) genes were identified by a genome-wide search using recently released wheat genomic data. Phylogenetic analysis grouped these genes into four clusters (Ι, II, III, IV), containing 22, 19, 7 and 8 genes, respectively. Chromosomal distribution, gene structure, and conserved motif analyses of the 56 TaHAK genes were subsequently performed. In silico RNA-seq data analysis revealed that TaHAKs from clusters II and III are constitutively expressed in various wheat tissues, while most genes from clusters I and IV have very low expression levels in the examined tissues at different developmental stages. qRT-PCR analysis showed that expression levels of TaHAK genes in wheat seedlings were significantly up- or downregulated when seedlings were exposed to K⁺ deficiency, high salinity, or dehydration. Furthermore, we functionally characterized TaHAK1b-2BL and showed that it facilitates K⁺ transport in yeast. Collectively, these results provide valuable information for further functional studies of TaHAKs, and contribute to a better understanding of the molecular basis of wheat development and stress tolerance.

Genetic diversity, population structure and marker-trait associations for agronomic and grain traits in wild diploid wheat Triticum urartu.

BACKGROUND: Wild diploid wheat, Triticum urartu (T. urartu) is the progenitor of bread wheat, and understanding its genetic diversity and genome function will provide considerable reference for dissecting genomic information of common wheat. RESULTS: In this study, we investigated the morphological and genetic diversity and population structure of 238 T. urartu accessions collected from different geographic regions. This collection had 19.37 alleles per SSR locus and its polymorphic information content (PIC) value was 0.76, and the PIC and Nei's gene diversity (GD) of high-molecular-weight glutenin subunits (HMW-GSs) were 0.86 and 0.88, respectively. UPGMA clustering analysis indicated that the 238 T. urartu accessions could be classified into two subpopulations, of which Cluster I contained accessions from Eastern Mediterranean coast and those from Mesopotamia and Transcaucasia belonged to Cluster II. The wide range of genetic diversity along with the manageable number of accessions makes it one of the best collections for mining valuable genes based on marker-trait association. Significant associations were observed between simple sequence repeats (SSR) or HMW-GSs and six morphological traits: heading date (HD), plant height (PH), spike length (SPL), spikelet number per spike (SPLN), tiller angle (TA) and grain length (GL). CONCLUSIONS: Our data demonstrated that SSRs and HMW-GSs were useful markers for identification of beneficial genes controlling important traits in T. urartu, and subsequently for their conservation and future utilization, which may be useful for genetic improvement of the cultivated hexaploid wheat.

Development of an integrated linkage map of einkorn wheat and its application for QTL mapping and genome sequence anchoring.

KEY MESSAGE: An integrated genetic map was constructed for einkorn wheat A genome and provided valuable information for QTL mapping and genome sequence anchoring. Wheat is one of the most widely grown food grain crops in the world. The construction of a genetic map is a key step to organize biologically or agronomically important traits along the chromosomes. In the present study, an integrated linkage map of einkorn wheat was developed using 109 recombinant inbred lines (RILs) derived from an inter sub-specific cross, KT1-1 (T. monococcum ssp. boeoticum) × KT3-5 (T. monococcum ssp. monococcum). The map contains 926 molecular markers assigned to seven linkage groups, and covers 1,377 cM with an average marker interval of 1.5 cM. A quantitative trait locus (QTL) analysis of five agronomic traits identified 16 stable QTL on all seven chromosomes, except 6A. The total phenotypic variance explained by these stable QTL using multiple regressions varied across environments from 8.8 to 87.1 % for days to heading, 24.4-63.0 % for spike length, 48.2-79.6 % for spikelet number per spike, 13.1-48.1 % for plant architecture, and 12.2-26.5 % for plant height, revealing that much of the RIL phenotypic variation had been genetically dissected. Co-localizations of closely linked QTL for different traits were frequently observed, especially on 3A and 7A. The QTL on 3A, 5A and 7A were closely associated with Eps-A m 3, Vrn1 and Vrn3 loci, respectively. Furthermore, this genetic map facilitated the anchoring of 237 T. urartu scaffolds onto seven chromosomes with a physical length of 26.15 Mb. This map and the QTL data provide valuable genetic information to dissect important agronomic and developmental traits in diploid wheat and contribute to the genetic ordering of the genome assembly.

Composition, variation, expression and evolution of low-molecular-weight glutenin subunit genes in Triticum urartu.

BACKGROUND: Wheat (AABBDD, 2n = 6x = 42) is a major dietary component for many populations across the world. Bread-making quality of wheat is mainly determined by glutenin subunits, but it remains challenging to elucidate the composition and variation of low-molecular-weight glutenin subunits (LMW-GS) genes, the major components for glutenin subunits in hexaploid wheat. This problem, however, can be greatly simplified by characterizing the LMW-GS genes in Triticum urartu, the A-genome donor of hexaploid wheat. In the present study, we exploited the high-throughput molecular marker system, gene cloning, proteomic methods and molecular evolutionary genetic analysis to reveal the composition, variation, expression and evolution of LMW-GS genes in a T. urartu population from the Fertile Crescent region. RESULTS: Eight LMW-GS genes, including four m-type, one s-type and three i-type, were characterized in the T. urartu population. Six or seven genes, the highest number at the Glu-A3 locus, were detected in each accession. Three i-type genes, each containing more than six allelic variants, were tightly linked because of their co-segregation in every accession. Only 2-3 allelic variants were detected for each m- and s-type gene. The m-type gene, TuA3-385, for which homologs were previously characterized only at Glu-D3 locus in common wheat and Aegilops tauschii, was detected at Glu-A3 locus in T. urartu. TuA3-460 was the first s-type gene identified at Glu-A3 locus. Proteomic analysis showed 1-4 genes, mainly i-type, expressed in individual accessions. About 62% accessions had three active i-type genes, rather than one or two in common wheat. Southeastern Turkey might be the center of origin and diversity for T. urartu due to its abundance of LMW-GS genes/genotypes. Phylogenetic reconstruction demonstrated that the characterized T. urartu might be the direct donor of the Glu-A3 locus in common wheat varieties. CONCLUSIONS: Compared with the Glu-A3 locus in common wheat, a large number of highly diverse LMW-GS genes and active genes were characterized in T. urartu, demonstrating that this progenitor might provide valuable genetic resources for LMW-GS genes to improve the quality of common wheat. The phylogenetic analysis provided molecular evidence and confirmed that T. urartu was the A-genome donor of hexaploid wheat.

Novel natural allelic variations at the Rht-1 loci in wheat.

Plant height is an important agronomic trait. Dramatic increase in wheat yield during the "green revolution" is mainly due to the widespread utilization of the Reduced height (Rht)-1 gene. We analyzed the natural allelic variations of three homoeologous loci Rht-A1, Rht-B1, and Rht-D1 in Chinese wheat (Triticum aestivum L.) micro-core collections and the Rht-B1/D1 genotypes in over 1,500 bred cultivars and germplasms using a modified EcoTILLING. We identified six new Rht-A1 allelic variations (Rht-A1b-g), eight new Rht-B1 allelic variations (Rht-B1h-o), and six new Rht-D1 allelic variations (Rht-D1e-j). These allelic variations contain single nucleotide polymorphisms (SNPs) or small insertions and deletions in the coding or uncoding regions, involving two frame-shift mutations and 15 missenses. Of which, Rht-D1e and Rht-D1h resulted in the loss of interactions of GID1-DELLA-GID2, Rht-B1i could increase plant height. We found that the Rht-B1h contains the same SNPs and 197 bp fragment insertion as reported in Rht-B1c. Further detection of Rht-B1h in Tibet wheat germplasms and wheat relatives indicated that Rht-B1c may originate from Rht-B1h. These results suggest rich genetic diversity at the Rht-1 loci and provide new resources for wheat breeding.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement.

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Genome-Wide Identification and Functional Characterization of the Chloride Channel TaCLC Gene Family in Wheat (Triticum aestivum L.).

In plants, chloride channels (CLC) are involved in a series of specific functions, such as regulation of nutrient transport and stress tolerance. Members of the wheat Triticum aestivum L. CLC (TaCLC) gene family have been proposed to encode anion channels/transporters that may be related to nitrogen transportation. To better understand their roles, TaCLC family was screened and 23 TaCLC gene sequences were identified using a Hidden Markov Model in conjunction with wheat genome database. Gene structure, chromosome location, conserved motif, and expression pattern of the resulting family members were then analyzed. Phylogenetic analysis showed that the TaCLC family can be divided into two subclasses (I and II) and seven clusters (-a, -c1, -c2, -e, -f1, -f2, and -g2). Using a wheat RNA-seq database, the expression pattern of TaCLC family members was determined to be an inducible expression type. In addition, seven genes from seven different clusters were selected for quantitative real-time PCR (qRT-PCR) analysis under low nitrogen stress or salt stress conditions, respectively. The results indicated that the gene expression levels of this family were up-regulated under low nitrogen stress and salt stress, except the genes of TaCLC-c2 cluster which were from subfamily -c. The yeast complementary experiments illustrated that TaCLC-a-6AS-1, TaCLC-c1-3AS, and TaCLC-e-3AL all had anion transport functions for NO3 - or Cl-, and compensated the hypersensitivity of yeast GEF1 mutant strain YJR040w (Δgef1) in restoring anion-sensitive phenotype. This study establishes a theoretical foundation for further functional characterization of TaCLC genes and provides an initial reference for better understanding nitrate nitrogen transportation in wheat.

Genome-wide association, RNA-seq and iTRAQ analyses identify candidate genes controlling radicle length of wheat.

The radicle, present in the embryo of a seed, is the first root to emerge at germination, and its rapid growth is essential for establishment and survival of the seedling. However, there are few studies on the critical mechanisms underlying radicle and then radicle length in wheat seedlings, despite its importance as a food crop throughout the world. In the present study, 196 wheat accessions from the Huanghuai Wheat Region were screened to measure radicle length under 4 hydroponic culture environments over 3 years. Different expression genes and proteins (DEGs/DEPs) between accessions with extremely long [Yunong 949 (WRL1), Zhongyu 9,302 (WRL2)] and short roots [Yunong 201 (WRS1), Beijing 841 (WRS2)] were identified in 12 sets of root tissue samples by RNA-seq and iTRAQ (Isobaric tags for relative and absolute quantification). Phenotypic results showed that the elongation zone was significantly longer in root accessions with long roots compared to the short-rooted accessions. A genome-wide association study (GWAS) identified four stable chromosomal regions significantly associated with radicle length, among which 1A, 4A, and 7A chromosomes regions explained 7.17% to12.93% of the phenotypic variation. The omics studies identified the expression patterns of 24 DEGs/DEPs changed at both the transcriptional and protein levels. These DEGs/DEPs were mainly involved in carbon fixation in photosynthetic organisms, photosynthesis and phenylpropanoid biosynthesis pathways. TraesCS1A02G104100 and TraesCS2B02G519100 were involved in the biosynthesis of tricin-lignins in cell walls and may affect the extension of cell walls in the radicle elongation zone. A combination of GWAS and RNA-seq analyses revealed 19 DEGs with expression changes in the four accessions, among which, TraesCS1A02G422700 (a cysteine-rich receptor-like protein kinase 6, CRK6) also showed upregulation in the comparison group by RNA-seq, iTRAQ, and qRT-PCR. BSMV-mediated gene silencing also showed that TaCRK6 improves root development in wheat. Our data suggest that TaCRK6 is a candidate gene regulating radicle length in wheat.

Comparative analysis of mitochondrial genomes between a wheat K-type cytoplasmic male sterility (CMS) line and its maintainer line.

BACKGROUND: Plant mitochondria, semiautonomous organelles that function as manufacturers of cellular ATP, have their own genome that has a slow rate of evolution and rapid rearrangement. Cytoplasmic male sterility (CMS), a common phenotype in higher plants, is closely associated with rearrangements in mitochondrial DNA (mtDNA), and is widely used to produce F1 hybrid seeds in a variety of valuable crop species. Novel chimeric genes deduced from mtDNA rearrangements causing CMS have been identified in several plants, such as rice, sunflower, pepper, and rapeseed, but there are very few reports about mtDNA rearrangements in wheat. In the present work, we describe the mitochondrial genome of a wheat K-type CMS line and compare it with its maintainer line. RESULTS: The complete mtDNA sequence of a wheat K-type (with cytoplasm of Aegilops kotschyi) CMS line, Ks3, was assembled into a master circle (MC) molecule of 647,559 bp and found to harbor 34 known protein-coding genes, three rRNAs (18 S, 26 S, and 5 S rRNAs), and 16 different tRNAs. Compared to our previously published sequence of a K-type maintainer line, Km3, we detected Ks3-specific mtDNA (> 100 bp, 11.38%) and repeats (> 100 bp, 29 units) as well as genes that are unique to each line: rpl5 was missing in Ks3 and trnH was absent from Km3. We also defined 32 single nucleotide polymorphisms (SNPs) in 13 protein-coding, albeit functionally irrelevant, genes, and predicted 22 unique ORFs in Ks3, representing potential candidates for K-type CMS. All these sequence variations are candidates for involvement in CMS. A comparative analysis of the mtDNA of several angiosperms, including those from Ks3, Km3, rice, maize, Arabidopsis thaliana, and rapeseed, showed that non-coding sequences of higher plants had mostly divergent multiple reorganizations during the mtDNA evolution of higher plants. CONCLUSION: The complete mitochondrial genome of the wheat K-type CMS line Ks3 is very different from that of its maintainer line Km3, especially in non-coding sequences. Sequence rearrangement has produced novel chimeric ORFs, which may be candidate genes for CMS. Comparative analysis of several angiosperm mtDNAs indicated that non-coding sequences are the most frequently reorganized during mtDNA evolution in higher plants.