Grain yield improvement by genome editing of TaARF12 that decoupled peduncle and rachis development trajectories via differential regulation of gibberellin signalling in wheat.

Plant breeding is constrained by trade-offs among different agronomic traits by the pleiotropic nature of many genes. Genes that contribute to two or more favourable traits with no penalty on yield are rarely reported, especially in wheat. Here, we describe the editing of a wheat auxin response factor TaARF12 by using CRISPR/Cas9 that rendered shorter plant height with larger spikes. Changes in plant architecture enhanced grain number per spike up to 14.7% with significantly higher thousand-grain weight and up to 11.1% of yield increase under field trials. Weighted Gene Co-Expression Network Analysis (WGCNA) of spatial-temporal transcriptome profiles revealed two hub genes: RhtL1, a DELLA domain-free Rht-1 paralog, which was up-regulated in peduncle, and TaNGR5, an organ size regulator that was up-regulated in rachis, in taarf12 plants. The up-regulation of RhtL1 in peduncle suggested the repression of GA signalling, whereas up-regulation of TaNGR5 in spike may promote GA response, a working model supported by differential expression patterns of GA biogenesis genes in the two tissues. Thus, TaARF12 complemented plant height reduction with larger spikes that gave higher grain yield. Manipulation of TaARF12 may represent a new strategy in trait pyramiding for yield improvement in wheat.

InDels Identification and Association Analysis with Spike and Awn Length in Chinese Wheat Mini-Core Collection.

Diversity surveys of germplasm are important for gaining insight into the genomic basis for crop improvement; especially InDels, which are poorly understood in hexaploid common wheat. Here, we describe a map of 89,923 InDels from exome sequencing of 262 accessions of a Chinese wheat mini-core collection. Population structure analysis, principal component analysis and selective sweep analysis between landraces and cultivars were performed. Further genome-wide association study (GWAS) identified five QTL (Quantitative Trait Loci) that were associated with spike length, two of them, on chromosomes 2B and 6A, were detected in 10 phenotypic data sets. Assisted with RNA-seq data, we identified 14 and 21 genes, respectively that expressed in spike and rachis within the two QTL regions that can be further investigated for candidate genes discovery. Moreover, InDels were found to be associated with awn length on chromosomes 5A, 6B and 4A, which overlapped with previously reported genetic loci B1 (Tipped 1), B2 (Tipped 2) and Hd (Hooded). One of the genes TaAGL6 that was previously shown to affect floral organ development was found at the B2 locus to affect awn length development. Our study shows that trait-associated InDels may contribute to wheat improvement and may be valuable molecular markers for future wheat breeding.

Genome-wide analysis of wheat DNA-binding with one finger (Dof) transcription factor genes: evolutionary characteristics and diverse abiotic stress responses.

BACKGROUND: DNA binding with one finger (Dof) transcription factors play important roles in plant growth and abiotic stress responses. Although genome-wide identification and analysis of the DOF transcription factor family has been reported in other species, no relevant studies have emerged in wheat. The aim of this study was to investigate the evolutionary and functional characteristics associated with plant growth and abiotic stress responses by genome-wide analysis of the wheat Dof transcription factor gene family. RESULTS: Using the recently released wheat genome database (IWGSC RefSeq v1.0), we identified 96 wheat Dof gene family members, which were phylogenetically clustered into five distinct subfamilies. Gene duplication analysis revealed a broad and heterogeneous distribution of TaDofs on the chromosome groups 1 to 7, and obvious tandem duplication genes were present on chromosomes 2 and 3.Members of the same gene subfamily had similar exon-intron structures, while members of different subfamilies had obvious differences. Functional divergence analysis indicated that type-II functional divergence played a major role in the differentiation of the TaDof gene family. Positive selection analysis revealed that the Dof gene family experienced different degrees of positive selection pressure during the process of evolution, and five significant positive selection sites (30A, 31 T, 33A, 102G and 104S) were identified. Additionally, nine groups of coevolving amino acid sites, which may play a key role in maintaining the structural and functional stability of Dof proteins, were identified. The results from the RNA-seq data and qRT-PCR analysis revealed that TaDof genes exhibited obvious expression preference or specificity in different organs and developmental stages, as well as in diverse abiotic stress responses. Most TaDof genes were significantly upregulated by heat, PEG and heavy metal stresses. CONCLUSIONS: The genome-wide analysis and identification of wheat DOF transcription factor family and the discovery of important amino acid sites are expected to provide new insights into the structure, evolution and function of the plant Dof gene family.

Comprehensive Comparative Analysis of the JAZ Gene Family in Common Wheat (Triticum aestivum) and Its D-Subgenome Donor Aegilops tauschii.

JASMONATE-ZIM DOMAIN (JAZ) repressor proteins work as co-receptors in the jasmonic acid (JA) signalling pathway and are essential for plant development and environmental adaptation. Despite wheat being one of the main staple food crops, until recently, comprehensive analysis of its JAZ gene family has been limited due to the lack of complete and high-quality reference genomes. Here, using the latest reference genome, we identified 17 JAZ genes in the wheat D-genome donor Aegilops tauschii. Then, 54 TaJAZs were identified in common wheat. A systematic examination of the gene structures, conserved protein domains, and phylogenetic relationships of this gene family was performed. Five new JAZ genes were identified as being derived from tandem duplication after wheat divergence from other species. We integrated RNA-seq data and yield QTL information and found that tandemly duplicated TaJAZ genes were prone to association with spike-related traits. Moreover, 12 TaJAZ genes were located within breeding selection sweeps, including 9 tandemly duplicated ones. Haplotype variation analysis of selected JAZ genes showed significant association of TaJAZ7A and TaJAZ13A with thousand-grain weight. Our work provides a clearer picture of wheat JAZ gene evolution and puts forward the possibility of using these genes for wheat yield improvement.

Genome-editing of a circadian clock gene TaPRR95 facilitates wheat peduncle growth and heading date.

Plant height and heading date are important agronomic traits in wheat (Triticum aestivum L.) that affect final grain yield. In wheat, knowledge of pseudo-response regulator (PRR) genes on agronomic traits is limited. Here, we identify a wheat TaPRR95 gene by genome-wide association study to be associated with plant height. Triple allele mutant plants produced by CRISPR/Cas9 show increased plant height, particularly the peduncle, with an earlier heading date. The longer peduncle is mainly caused by the increased cell elongation at its upper section, whilst the early heading date is accompanied by elevated expression of flowering genes, such as TaFT and TaCO1. A peduncle-specific transcriptome analysis reveals up-regulated photosynthesis genes and down-regulated IAA/Aux genes for auxin signaling in prr95aabbdd plants that may act as a regulatory mechanism to promote robust plant growth. A haplotype analysis identifies a TaPRR95-B haplotype (Hap2) to be closely associated with reduced plant height and increased thousand-grain weight. Moreover, the Hap2 frequency is higher in cultivars than that in landraces, suggesting the artificial selection on the allele during wheat breeding. These findings suggest that TaPRR95 is a regulator for plant height and heading date, thereby providing an important target for wheat yield improvement.

A Multi-Omics Approach for Rapid Identification of Large Genomic Lesions at the Wheat Dense Spike (wds) Locus.

Optimal spike architecture provides a favorable structure for grain development and yield improvement. However, the number of genes cloned to underlie wheat spike architecture is extremely limited. Here, we obtained a wheat dense spike mutant (wds) induced by 60Co treatment of a common wheat landrace Huangfangzhu that exhibited significantly reduced spike and grain lengths. The shortened spike length was caused by longitudinal reduction in number and length of rachis cells. We adopted a multi-omics approach to identify the genomic locus underlying the wds mutant. We performed Exome Capture Sequencing (ECS) and identified two large deletion segments, named 6BL.1 at 334.8∼424.3 Mb and 6BL.2, 579.4∼717.8 Mb in the wds mutant. RNA-seq analysis confirmed that genes located in these regions lost their RNA expression. We then found that the 6BL.2 locus was overlapping with a known spike length QTL, qSL6B.2. Totally, 499 genes were located within the deleted region and two of them were found to be positively correlated with long spike accessions but not the ones with short spike. One of them, TraesCS6B01G334600, a well-matched homolog of the rice OsBUL1 gene that works in the Brassinosteroids (BR) pathway, was identified to be involved in cell size and number regulation. Further transcriptome analysis of young spikes showed that hormone-related genes were enriched among differentially expressed genes, supporting TraesCS6B01G334600 as a candidate gene. Our work provides a strategy to rapid locate genetic loci with large genomic lesions in wheat and useful resources for future wheat study.