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.

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.

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.

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.

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.

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.

The abundance of homoeologue transcripts is disrupted by hybridization and is partially restored by genome doubling in synthetic hexaploid wheat.

BACKGROUND: The formation of an allopolyploid is a two step process, comprising an initial wide hybridization event, which is later followed by a whole genome doubling. Both processes can affect the transcription of homoeologues. Here, RNA-Seq was used to obtain the genome-wide leaf transcriptome of two independent Triticum turgidum × Aegilops tauschii allotriploids (F1), along with their spontaneous allohexaploids (S1) and their parental lines. The resulting sequence data were then used to characterize variation in homoeologue transcript abundance. RESULTS: The hybridization event strongly down-regulated D-subgenome homoeologues, but this effect was in many cases reversed by whole genome doubling. The suppression of D-subgenome homoeologue transcription resulted in a marked frequency of parental transcription level dominance, especially with respect to genes encoding proteins involved in photosynthesis. Singletons (genes where no homoeologues were present) were frequently transcribed at both the allotriploid and allohexaploid plants. CONCLUSIONS: The implication is that whole genome doubling helps to overcome the phenotypic weakness of the allotriploid, restoring a more favourable gene dosage in genes experiencing transcription level dominance in hexaploid wheat.

A MADS-box gene NtSVP regulates pedicel elongation by directly suppressing a KNAT1-like KNOX gene NtBPL in tobacco (Nicotiana tabacum L.).

Optimal inflorescence architecture is important for plant reproductive success by affecting the ultimate number of flowers that set fruits and for plant competitiveness when interacting with biotic or abiotic conditions. The pedicel is one of the key contributors to inflorescence architecture diversity. To date, knowledge about the molecular mechanisms of pedicel development is derived from Arabidopsis. Not much is known regarding other plants. Here, an SVP family MADS-box gene, NtSVP, in tobacco (Nicotiana tabacum) that is required for pedicel elongation was identified. It is shown that knockdown of NtSVP by RNA interference (RNAi) caused elongated pedicels, while overexpression resulted in compact inflorescences with much shortened pedicels. Moreover, an Arabidopsis BREVIPEDECELLUS/KNAT1 homologue NtBP-Like (NtBPL) was significantly up-regulated in NtSVP-RNAi plants. Disruption of NtBPL decreased pedicel lengths and shortened cortex cells. Consistent with the presence of a CArG-box at the NtBPL promoter, the direct binding of NtSVP to the NtBPL promoter was demonstrated by yeast one-hybrid assay, electrophoretic mobility shift assay, and dual-luciferase assay, in which NtSVP may act as a repressor of NtBPL. Microarray analysis showed that down-regulation of NtBPL resulted in differential expression of genes associated with a number of hormone biogenesis and signalling genes such as those for auxin and gibberellin. These findings together suggest the function of a MADS-box transcription factor in plant pedicel development, probably via negative regulation of a BP-like class I KNOX gene. The present work thus postulates the conservation and divergence of the molecular regulatory pathways underlying the development of plant inflorescence architecture.

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.

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.

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.

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.

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.

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.

Genome-wide quantitative trait locus mapping identifies multiple major loci for brittle rachis and threshability in Tibetan semi-wild wheat (Triticum aestivum ssp. tibetanum Shao).

Tibetan semi-wild wheat (Triticum aestivum ssp. tibetanum Shao) is a semi-wild hexaploid wheat resource that is only naturally distributed in the Qinghai-Tibet Plateau. Brittle rachis and hard threshing are two important characters of Tibetan semi-wild wheat. A whole-genome linkage map of T. aestivum ssp. tibetanum was constructed using a recombinant inbred line population (Q1028×ZM9023) with 186 lines, 564 diversity array technology markers, and 117 simple sequence repeat markers. Phenotypic data on brittle rachis and threshability, as two quantitative traits, were evaluated on the basis of the number of average spike rachis fragments per spike and percent threshability in 2012 and 2013, respectively. Quantitative trait locus (QTL) mapping performed using inclusive composite interval mapping analysis clearly identified four QTLs for brittle rachis and three QTLs for threshability. However, three loci on 2DS, 2DL, and 5AL showed pleiotropism for brittle rachis and threshability; they respectively explained 5.3%, 18.6%, and 18.6% of phenotypic variation for brittle rachis and 17.4%, 13.2%, and 35.2% of phenotypic variation for threshability. A locus on 3DS showed an independent effect on brittle rachis, which explained 38.7% of the phenotypic variation. The loci on 2DS and 3DS probably represented the effect of Tg and Br1, respectively. The locus on 5AL was in very close proximity to the Q gene, but was different from the predicted q in Tibetan semi-wild wheat. To our knowledge, the locus on 2DL has never been reported in common wheat but was prominent in T. aestivum ssp. tibetanum accession Q1028. It remarkably interacted with the locus on 5AL to affect brittle rachis. Several major loci for brittle rachis and threshability were identified in Tibetan semi-wild wheat, improving the understanding of these two characters and suggesting the occurrence of special evolution in Tibetan semi-wild wheat.