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 (GWAS) to be associated with plant height. Triple allele mutant plants produced by CRISPR/Cas9 show increased plant height, particularly at 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 with 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 new regulator for plant height and heading date, thereby providing an important target for wheat yield improvement.

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.

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.

TaSPL14-7A is a conserved regulator controlling plant architecture and yield traits in common wheat (Triticum aestivum L.).

Plant architecture is a crucial influencing factor of wheat yield and adaptation. In this study, we cloned and characterized TaSPL14, a homologous gene of the rice ideal plant architecture gene OsSPL14 in wheat. TaSPL14 homoeologs (TaSPL14-7A, TaSPL14-7B and TaSPL14-7D) exhibited similar expression patterns, and they were all preferentially expressed in stems at the elongation stage and in young spikes. Moreover, the expression level of TaSPL14-7A was higher than that of TaSPL14-7B and TaSPL14-7D. Overexpression of TaSPL14-7A in wheat resulted in significant changes in plant architecture and yield traits, including decreased tiller number and increased kernel size and weight. Three TaSPL14-7A haplotypes were identified in Chinese wheat core collection, and haplotype-based association analysis showed that TaSPL14-7A-Hap1/2 were significantly correlated with fewer tillers, larger kernels and higher kernel weights in modern cultivars. The haplotype effect resulted from a difference in TaSPL14-7A expression levels among genotypes, with TaSPL14-7A-Hap1/2 leading to higher expression levels than TaSPL14-7A-Hap3. As favorable haplotypes, TaSPL14-7A-Hap1/2 underwent positive selection during global wheat breeding over the last century. Together, the findings of our study provide insight into the function and genetic effects of TaSPL14 and provide a useful molecular marker for wheat breeding.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat.

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

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.

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.

Transcriptome Analysis of Developing Wheat Grains at Rapid Expanding Phase Reveals Dynamic Gene Expression Patterns.

Grain development, as a vital process in the crop's life cycle, is crucial for determining crop quality and yield. The wheat grain expanding phase is the early process involving the rapid morphological changes and initiation of grain filling. However, little is known about the molecular basis of grain development at this stage. Here, we provide a time-series transcriptome profile of developing wheat grain at 0, 2, 4, 6, 8, and 10 days after pollination of the wheat landrace Chinese Spring. A total of 26,892 differentially expressed genes, including 1468 transcription factors, were found between adjacent time points. Co-expression cluster analysis and Gene Ontology enrichment revealed dynamic expressions of cell division and starch biosynthesis related structural genes and transcription factors. Moreover, diverse, differential and drastically varied expression trends of the key genes related to hormone metabolism were identified. Furthermore, ~30% of triads showed unbalanced expression patterns enriching for genes in multiple pivotal metabolic pathways. Hormone metabolism related genes, such as YUC10 (YUCCA flavin-containing monooxygenase 10), AOS2 (allene oxide synthase 2), CYP90D2 (cytochrome P450 90D2), and CKX1 (cytokinin dehydrogenase 1), were dominantly contributed by A or D homoeologs of the triads. Our study provided a systematic picture of transcriptional regulation of wheat grains at the early grain expanding phase which should deepen our understanding of wheat grain development and help in wheat yield improvement.

Matrilineal empowers wheat pollen with haploid induction potency by triggering postmitosis reactive oxygen species activity.

Reactive oxygen species (ROS) play important roles during anther and pollen development. DNA damage may cause chromosome fragmentation that is considered to underlie chromosome elimination for haploid induction by matrilineal pollen, a key step in MATRILINEAL-based double haploid breeding technology. But when and how DNA damage occurs is unknown. We performed comparative studies of wheat pollens from the wild-type and the CRISPR/Cas9 edited matrilineal mutant (mMTL). Chemical assays detected a second wave of ROS in mMTL pollen at the three-nuclei-stage and subsequently, along with reduced antioxidant enzyme activities. RNA-seq analysis revealed disturbed expression of genes for fatty acid biosynthesis and ROS homoeostasis. Gas chromatography-mass spectrometry measurement identified abnormal fatty acid metabolism that may contribute to defective mMTL pollen walls as observed using electron microscopy, consistent with the function of MTL as a phospholipase. Moreover, DNA damage was identified using TdT-mediated dUTP nick-end labelling and quantified using comet assays. Velocity patterns showed that ROS increments preceded that of DNA damage over the course of pollen maturation. Our work hypothesises that mMTL-triggered later-stage-specific ROS causes DNA damage that may contribute to chromosome fragmentation and hence chromosome elimination during haploid induction. These findings may provide more ways to accelerate double haploid-based plant breeding.

Wheat breeding history reveals synergistic selection of pleiotropic genomic sites for plant architecture and grain yield.

Diversity surveys of crop germplasm are important for gaining insights into the genomic basis for plant architecture and grain yield improvement, which is still poorly understood in wheat. In this study, we exome sequenced 287 wheat accessions that were collected in the past 100 years. Population genetics analysis identified that 6.7% of the wheat genome falls within the selective sweeps between landraces and cultivars, which harbors the genes known for yield improvement. These regions were asymmetrically distributed on the A and B subgenomes with regulatory genes being favorably selected. Genome-wide association study (GWAS) identified genomic loci associated with traits for yield potential, and two underlying genes, TaARF12 encoding an auxin response factor and TaDEP1 encoding the G-protein γ-subunit, were located and characterized to pleiotropically regulate both plant height and grain weight. Elite single-nucleotide haplotypes with increased allele frequency in cultivars relative to the landraces were identified and found to have accumulated over the course of breeding. Interestingly, we found that TaARF12 and TaDEP1 function in epistasis with the classical plant height Rht-1 locus, leading to propose a "Green Revolution"-based working model for historical wheat breeding. Collectively, our study identifies selection signatures that fine-tune the gibberellin pathway during modern wheat breeding and provides a wealth of genomic diversity resources for the wheat research community.

The wheat AGL6-like MADS-box gene is a master regulator for floral organ identity and a target for spikelet meristem development manipulation.

The AGAMOUS-LIKE6 (AGL6)-like genes are ancient MADS-box genes and are functionally studied in a few model plants. The knowledge of these genes in wheat remains limited. Here, by studying a 'double homoeolog mutant' of the AGL6 gene in tetraploid wheat, we showed that AGL6 was required for the development of all four whorls of floral organs with dosage-dependent effect on floret fertility. Yeast two-hybrid analyses detected interactions of AGL6 with all classes of MADS-box proteins in the ABCDE model for floral organ development. AGL6 was found to interact with several additional proteins, including the G protein ß and γ (DEP1) subunits. Analysis of the DEP1-B mutant showed a significant reduction in spikelet number per spike in tetraploid wheat, while overexpression of AGL6 in common wheat increased the spikelet number per spike and hence the grain number per spike. RNA-seq analysis identified the regulation of several meristem activity genes by AGL6, such as FUL2 and TaMADS55. Our work therefore extensively updated the wheat ABCDE model and proposed an alternative approach to improve wheat grain yield by manipulating the AGL6 gene.

TaIAA21 represses TaARF25-mediated expression of TaERFs required for grain size and weight development in wheat.

Auxin signaling is essential for the development of grain size and grain weight, two important components for crop yield. However, no auxin/indole acetic acid repressor (Aux/IAA) has been functionally characterized to be involved in the development of wheat (Triticum aestivum L.) grains to date. Here, we identified a wheat Aux/IAA gene, TaIAA21, and studied its regulatory pathway. We found that TaIAA21 mutation significantly increased grain length, grain width, and grain weight. Cross-sections of mutant grains revealed elongated outer pericarp cells compared to those of the wild type, where the expression of TaIAA21 was detected by in situ hybridization. Screening of auxin response factor (ARF) genes highly expressed in early developing grains revealed that TaARF25 interacts with TaIAA21. In contrast, mutation of the tetraploid wheat (Triticum turgidum) ARF25 gene significantly reduced grain size and weight. RNA sequencing analysis revealed upregulation of several ethylene response factor genes (ERFs) in taiaa21 mutants which carried auxin response cis-elements in their promoter. One of them, ERF3, was upregulated in the taiaa21 mutant and downregulated in the ttarf25 mutant. Transactivation assays showed that ARF25 promotes ERF3 transcription, while mutation of TtERF3 resulted in reduced grain size and weight. Analysis of natural variations identified three TaIAA21-A haplotypes with increased allele frequencies in cultivars relative to landraces, a signature of breeding selection. Our work demonstrates that TaIAA21 works as a negative regulator of grain size and weight development via the ARF25-ERFs module and is useful for yield improvement in wheat.

RNAi technology for plant protection and its application in wheat.

The RNAi technology takes advantage of the intrinsic RNA interference (RNAi) mechanism that exists in nearly all eukaryotes in which target mRNAs are degraded or functionally suppressed. Significant progress has been made in recent years where RNAi technology is applied to several crops and economic plants for protection against diseases like fungi, pests, and nematode. RNAi technology is also applied in controlling pathogen damages in wheat, one of the most important crops in the world. In this review, we first give a brief introduction of the RNAi technology and the underneath mechanism. We then review the recent progress of its utilization in crops, particular wheat. Finally, we discuss the existing challenges and prospect future development of this technology in crop protection.

Comprehensive analysis of Q gene near-isogenic lines reveals key molecular pathways for wheat domestication and improvement.

The wheat AP2-like transcription factor gene Q has played a major role in domestication by conferring the free-threshing character and pleiotropically affecting numerous other traits. However, little information is known regarding the molecular mechanisms associated with the regulation of these traits by Q, especially for the structural determination of threshability. Here, transcriptome analysis of immature spike tissues in three lines nearly isogenic for Q revealed over 3000 differentially expressed genes (DEGs) involved in a number of pathways. Using phenotypic, microscopic, transcriptomic, and tissue-specific gene expression analyses, we demonstrated that Q governs threshability through extensive modification of wheat glumes including their structure, cell wall thickness, and chemical composition. Critical DEGs and pathways involved in secondary cell wall synthesis and regulation of the chemical composition of glumes were identified. We also showed that the mutation giving rise to the Q allele synchronized the expression of genes for micro-sporogenesis that affected pollen fertility, and may determine the final grain number for wheat spikes. Transcriptome dissection of genes and genetic pathways regulated by Q should further our understanding of wheat domestication and improvement.

The genome-wide transcriptional consequences of the nullisomic-tetrasomic stocks for homoeologous group 7 in bread wheat.

BACKGROUND: Hexaploid bread wheat (Triticum aestivum L) arose by two polyploidisation events from three diploid species with homoeologous genomes. Nullisomic-tetrasomic (nulli-tetra or NT) lines are aneuploid wheat plants lacking two and adding two of six homoeologous chromosomes. These plants can grow normally, but with significantly morphological variations because the adding two chromosomes or the remaining four chromosomes compensate for those absent. Despite these interesting phenomena, detailed molecular mechanisms underlying dosage deletion and compensation in these useful genetic materials have not been determined. RESULTS: By sequencing the transcriptomes of leaves in two-week-old seedlings, we showed that the profiles of differentially expressed genes between NT stocks for homoeologous group 7 and the parent hexaploid Chinese Spring (CS) occurred throughout the whole genome with a subgenome and chromosome preference. The deletion effect of nulli-chromosomes was compensated partly by the tetra-chromosomes via the dose level of expressed genes, according to the types of homoeologous genes. The functions of differentially regulated genes primarily focused on carbon metabolic process, photosynthesis process, hormone metabolism, and responding to stimulus, and etc., which might be related to the defective phenotypes that included reductions in plant height, flag leaf length, spikelet number, and kernels per spike. CONCLUSIONS: The perturbation of the expression levels of transcriptional genes among the NT stocks for homoeologous group 7 demonstrated the gene dosage effect of the subgenome at the genome-wide level. The gene dosage deletion and compensation can be used as a model to elucidate the functions of the subgenomes in modern polyploid plants.

DNA methylation dynamics during the interaction of wheat progenitor Aegilops tauschii with the obligate biotrophic fungus Blumeria graminis f. sp. tritici.

DNA methylation is dynamically involved in plant immunity, but little information is known about its roles in plant interactions with biotrophic fungi, especially in temperate grasses such as wheat (Triticum aestivum). Using wheat diploid progenitor Aegilops tauschii accession AL8/78, the genome of which has been sequenced, we assessed the extent of DNA methylation in response to infection with Blumeria graminis f. sp. tritici (Bgt), which causes powdery mildew. Upon Bgt infection, ARGONAUTE4a (AGO4a) was significantly downregulated in A. tauschii, which was accompanied by a substantial reduction in AGO4a-sorted 24-nt siRNA levels, especially for genes near transposable elements (TAGs). Bisulfite sequencing revealed abundant differentially methylated regions (DMRs) with CHH hypomethylation. TAGs bearing CHH-hypomethylated DMRs were enriched for 'response to stress' functions, including receptor kinase, peroxidase, and pathogenesis-related genes. Virus-induced gene silencing (VIGS) of a DOMAINS REARRANGED METHYLASE 2 (DRM2) homolog enhanced plant resistance to Bgt. The effect of CHH hypomethylation was exemplified by the upregulation of a pathogenesis-related ß-1,3-glucanse gene implicated in Bgt defense. These findings support the idea that dynamic DNA methylation represents a regulatory layer in the complex mechanism of plant immunity, which could be exploited to improve disease resistance in common wheat.

Genome-Wide Identification and Expression Profiling of the TCP Family Genes in Spike and Grain Development of Wheat (Triticum aestivum L.).

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

Transcriptome Profiling of Wheat Inflorescence Development from Spikelet Initiation to Floral Patterning Identified Stage-Specific Regulatory Genes.

Early reproductive development in cereals is crucial for final grain number per spike and hence the yield potential of the crop. To date, however, no systematic analyses of gene expression profiles during this important process have been conducted for common wheat (Triticum aestivum). Here, we studied the transcriptome profiles at four stages of early wheat reproductive development, from spikelet initiation to floral organ differentiation. K-means clustering and stage-specific transcript identification detected dynamically expressed homeologs of important transcription regulators in spikelet and floral meristems that may be involved in spikelet initiation, floret meristem specification, and floral organ patterning, as inferred from their homologs in model plants. Small RNA transcriptome sequencing discovered key microRNAs that were differentially expressed during wheat inflorescence development alongside their target genes, suggesting that miRNA-mediated regulatory mechanisms for floral development may be conserved in cereals and Arabidopsis. Our analysis was further substantiated by the functional characterization of the ARGONAUTE1d (AGO1d) gene, which was initially expressed in stamen primordia and later in the tapetum during anther maturation. In agreement with its stage-specific expression pattern, the loss of function of the predominantly expressed B homeolog of AGO1d in a tetraploid durum wheat mutant resulted in smaller anthers with more infertile pollens than the wild type and a reduced grain number per spike. Together, our work provides a first glimpse of the gene regulatory networks in wheat inflorescence development that may be pivotal for floral and grain development, highlighting potential targets for genetic manipulation to improve future wheat yields.

Regulation of FT splicing by an endogenous cue in temperate grasses.

Appropriate flowering timing is crucial for plant reproductive success. The florigen, FLOWERING LOCUS T (FT), interacts with 14-3-3 proteins and the bZIP transcription factor FD, functioning at core nodes in multiple flowering pathways. There are two FT homologues, FT1 and FT2, in Brachypodium distachyon. Here we show that FT2 undergoes age-dependent alternative splicing (AS), resulting in two splice variants (FT2α and FT2ß). The FT2ß-encoded protein cannot interact with FD or 14-3-3s but is able to form heterodimers with FT2α and FT1, thereby interfering with the florigen-mediated assembly of the flowering initiation complex. Notably, transgenic plants overproducing FT2ß exhibit delayed flowering, while transgenic plants in which FT2ß is silenced by an artificial microRNA display accelerated flowering, demonstrating a dominant-negative role of FT2ß in flowering induction. Furthermore, we show that the AS splicing of FT2 is conserved in important cereal crops, such as barley and wheat. Collectively, these findings reveal a novel posttranscriptional mode of FT regulation in temperate grasses.