Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture.

KEY MESSAGE: A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Function analysis of transcription factor OSR1 regulating osmotic stress resistance in maize.

BACKGROUND: Maize is a major cereal crop world widely, however, the yield of maize is frequently limited by dehydration and even death of plants, which resulted from osmotic stress such as drought and salinity. Dissection of molecular mechanisms controlling stress tolerance will enable plant scientists and breeders to increase crops yield by manipulating key regulatory components. METHODS: The candidate OSR1 gene was identified by map-based cloning. The expression level of OSR1 was verified by qRT-PCR and digital PCR in WT and osr1 mutant. Electrophoretic mobility shift assay, transactivation activity assay, subcellular localization, transcriptome analysis and physiological characters measurements were conducted to analyze the function of OSR1 in osmotic stress resistance in maize. RESULTS: The osr1 mutant was significantly less sensitive to osmotic stress than the WT plants and displayed stronger water-holding capacity, and the OSR1 homologous mutant in Arabidopsis showed a phenotype similar with maize osr1 mutant. Differentially expressed genes (DEGs) were identified between WT and osr1 under osmotic stress by transcriptome analysis, the expression levels of many genes, such as LEA, auxin-related factors, PPR family members, and TPR family members, changed notably, which may primarily involve in osmotic stress or promote root development. CONCLUSIONS: OSR1 may serve as a negative regulatory factor in response to osmotic stress in maize. The present study sheds new light on the molecular mechanisms of osmotic stress in maize.

A platform for whole-genome speed introgression from Aegilops tauschii to wheat for breeding future crops.

Breeding new and sustainable crop cultivars of high yields and desirable traits has been a major challenge for ensuring food security for the growing global human population. For polyploid crops such as wheat, introducing genetic variation from wild relatives of its subgenomes is a key strategy to improve the quality of their breeding pools. Over the past decades, considerable progress has been made in speed breeding, genome sequencing, high-throughput phenotyping and genomics-assisted breeding, which now allows us to realize whole-genome introgression from wild relatives to modern crops. Here, we present a standardized protocol to rapidly introgress the entire genome of Aegilops tauschii, the progenitor of the D subgenome of bread wheat, into elite wheat backgrounds. This protocol integrates multiple modern high-throughput technologies and includes three major phases: development of synthetic octaploid wheat, generation of hexaploid A. tauschii-wheat introgression lines (A-WIs) and homozygosis of the generated A-WIs. Our approach readily generates stable introgression lines in 2 y, thus greatly accelerating the generation of A-WIs and the introduction of desirable genes from A. tauschii to wheat cultivars. These A-WIs are valuable for wheat-breeding programs and functional gene discovery. The current protocol can be easily modified and used for introgressing the genomes of wild relatives to other polyploid crops.

A maize epimerase modulates cell wall synthesis and glycosylation during stomatal morphogenesis.

The unique dumbbell-shape of grass guard cells (GCs) is controlled by their cell walls which enable their rapid responses to the environment. The molecular mechanisms regulating the synthesis and assembly of GC walls are as yet unknown. Here we have identified BZU3, a maize gene encoding UDP-glucose 4-epimerase that regulates the supply of UDP-glucose during GC wall synthesis. The BZU3 mutation leads to significant decreases in cellular UDP-glucose levels. Immunofluorescence intensities reporting levels of cellulose and mixed-linkage glucans are reduced in the GCs, resulting in impaired local wall thickening. BZU3 also catalyzes the epimerization of UDP-N-acetylgalactosamine to UDP-N-acetylglucosamine, and the BZU3 mutation affects N-glycosylation of proteins that may be involved in cell wall synthesis and signaling. Our results suggest that the spatiotemporal modulation of BZU3 plays a dual role in controlling cell wall synthesis and glycosylation via controlling UDP-glucose/N-acetylglucosamine homeostasis during stomatal morphogenesis. These findings provide insights into the mechanisms controlling formation of the unique morphology of grass stomata.

Characterization of Glutenin Genes in Bread Wheat by Third-Generation RNA Sequencing and the Development of a Glu-1Dx5 Marker Specific for the Extra Cysteine Residue.

High-molecular-weight glutenin subunits (HMW-GS) and low-molecular-weight glutenin subunits (LMW-GS) in a mature grain play important roles in the formation of a glutenin macropolymer and gluten quality. To characterize the expressed glutenin genes of the bread wheat variety Xinmai 26 during seed development, a total of 18 full-length transcripts were obtained by the newly emerged third-generation RNA sequencing of the PacBio Sequel II platform, including 5 transcripts of HMW-GS genes and 13 transcripts of LMW-GS genes (8 intact genes and 5 pseudogenes). Combined with the patterns of sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS), allelic types of the obtained glutenin genes were, respectively, determined, wherein molecular characterization deduced by transcript1528 (1Dx5) and transcript907 (Glu-A3c) indicated their great influence on dough quality. In addition, a specific functional marker dCAPS5 was developed for the single-nucleotide substitution at position 353 of the 1Dx5 subunit, which was further intensively compared with the other proposed markers to efficiently utilize the 1Dx5 subunit with the extra cysteine residue. This study provides an efficient method to accurately identify and utilize glutenin genes in bread wheat, which is helpful in understanding the contributions of glutenin genes to wheat quality.

Fine mapping of Pm58 from Aegilops tauschii conferring powdery mildew resistance.

KEY MESSAGE: The powdery mildew resistance gene Pm58 was traced to a 141.3-kb interval with the co-segregating marker Xkasp68500 in wheat breeding. Pm58 is a powdery mildew resistance gene identified in Aegilops tauschii accession TA1662 and effective in a common wheat background. To finely map Pm58, an F2 population of 676 plants derived from the cross T093 × TA1662 was used for recombinant screening. We obtained 13 recombinants that occurred between the flanking markers Xhnu670 and Xhnu186. Genotyping and phenotyping these recombinant F2:3 families delimited Pm58 to a 0.22-cM interval (Xsts20220-Xkasp61553) on chromosome arm 2DS. The region carrying the Pm58 locus was approximately 141.3-kb, which contained eight annotated genes according to the reference genome sequence of Ae. tauschii AL8/78. Haplotype analysis of 178 Ae. tauschii accessions using the candidate gene-specific markers identified a disease resistance gene AET2Gv20068500 as a candidate for Pm58. Comparative mapping of the Pm58-containing interval revealed two presence/absence variations (PAVs) between AL8/78 and common wheat Chinese Spring. PAV-1 resides in the 3'-end of AET2Gv20068500. The majority of 158 common wheat cultivars (84.8%) displayed the absence of a 14.1-kb fragment in the PAV-1 region, which was confirmed by aligning the targeted genome sequences of the other sequenced Ae. tauschii accessions and common wheat cultivars. A co-segregating marker Xkasp68500 developed from AET2Gv20068500 can distinguish TA1662 from all randomly selected common wheat cultivars and will be instrumental for tracking Pm58 in breeding programs.

The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata.

The unique morphology of grass stomata enables rapid responses to environmental changes. Deciphering the basis for these responses is critical for improving food security. We have developed a planta platform of single-nucleus RNA-sequencing by combined fluorescence-activated nuclei flow sorting, and used it to identify cell types in mature and developing stomata from 33,098 nuclei of the maize epidermis-enriched tissues. Guard cells (GCs) and subsidiary cells (SCs) displayed differential expression of genes, besides those encoding transporters, involved in the abscisic acid, CO2, Ca2+, starch metabolism, and blue light signaling pathways, implicating coordinated signal integration in speedy stomatal responses, and of genes affecting cell wall plasticity, implying a more sophisticated relationship between GCs and SCs in stomatal development and dumbbell-shaped guard cell formation. The trajectory of stomatal development identified in young tissues, and by comparison to the bulk RNA-seq data of the MUTE defective mutant in stomatal development, confirmed known features, and shed light on key participants in stomatal development. Our study provides a valuable, comprehensive, and fundamental foundation for further insights into grass stomatal function.

New insights into the dispersion history and adaptive evolution of taxon Aegilops tauschii in China.

Aegilops tauschii, the wild progenitor of wheat D-genome and a valuable germplasm for wheat improvement, has a wide natural distribution from eastern Turkey to China. However, the phylogenetic relationship and dispersion history of Ae. tauschii in China has not been scientifically clarified. In this study, we genotyped 208 accessions (with 104 in China) using ddRAD sequencing and 55K SNP array, and classified the population into six sublineages. Three possible spreading routes or events were identified, resulting in specific distribution patterns, with four sublineages found in Xinjiang, one in Qinghai, two in Shaanxi and one in Henan. We also established the correlation of SNP-based, karyotype-based and spike-morphology-based techniques to demonstrate the internal classification of Ae. tauschii, and developed consensus dataset with 1245 putative accessions by merging data previously published. Our analysis suggested that eight inter-lineage accessions could be assigned to the putative Lineage 3 and these accessions would help to conserve the genetic diversity of the species. By developing the consensus phylogenetic relationships of Ae. tauschii, our work validated the hypothesis on the dispersal history of Ae. tauschii in China, and contributed to the efficient and comprehensive germplasm-mining of the species.

Introgressing the Aegilops tauschii genome into wheat as a basis for cereal improvement.

Increasing crop production is necessary to feed the world's expanding population, and crop breeders often utilize genetic variations to improve crop yield and quality. However, the narrow diversity of the wheat D genome seriously restricts its selective breeding. A practical solution is to exploit the genomic variations of Aegilops tauschii via introgression. Here, we established a rapid introgression platform for transferring the overall genetic variations of A. tauschii to elite wheats, thereby enriching the wheat germplasm pool. To accelerate the process, we assembled four new reference genomes, resequenced 278 accessions of A. tauschii and constructed the variation landscape of this wheat progenitor species. Genome comparisons highlighted diverse functional genes or novel haplotypes with potential applications in wheat improvement. We constructed the core germplasm of A. tauschii, including 85 accessions covering more than 99% of the species' overall genetic variations. This was crossed with elite wheat cultivars to generate an A. tauschii-wheat synthetic octoploid wheat (A-WSOW) pool. Laboratory and field analysis with two examples of the introgression lines confirmed its great potential for wheat breeding. Our high-quality reference genomes, genomic variation landscape of A. tauschii and the A-WSOW pool provide valuable resources to facilitate gene discovery and breeding in wheat.

Characterization of PmHHXM, a New Broad-Spectrum Powdery Mildew Resistance Gene in Chinese Wheat Landrace Honghuaxiaomai.

Powdery mildew, caused by fungal pathogen Blumeria graminis f. sp. tritici, is an agronomically important and widespread wheat disease causing severe yield losses. Deployment of broad-spectrum disease resistance genes is the preferred strategy to prevent this pathogen. Chinese wheat landrace Honghuaxiaomai (HHXM) was resistant to all 23 tested B. graminis f. sp. tritici isolates at the seedling stage. The F1, F2, and F2:3 progenies derived from the cross HHXM × Yangmai 158 were used in this study, and genetic analysis revealed that a single dominant gene, designated PmHHXM, conferred resistance to B. graminis f. sp. tritici isolate E09. Bulked segregant analysis and molecular mapping initially located PmHHXM to the distal region of chromosome 4AL. To fine map PmHHXM, we identified two critical recombinants from 592 F2 plants and delimited PmHHXM to a 0.18-cM Xkasp475200 to Xhnu552 interval covering 1.77 Mb, in which a number of disease resistance-related gene clusters were annotated. Comparative mapping of this interval revealed a perturbed synteny among Triticeae species. This study reports the new powdery mildew resistance gene PmHHXM, which seems different from three known quantitative trait loci/genes identified on chromosome 4AL and has significant values for further genetic improvement. Analysis of the polymorphisms of 13 cosegregating markers between HHXM and 170 modern wheat cultivars indicates that Xhnu227 and Xsts478700 developed here are ideal for marker-assisted introgression of this locus in wheat breeding.

Comparative Transcriptome Analysis of Two Aegilops tauschii with Contrasting Drought Tolerance by RNA-Seq.

As the diploid progenitor of common wheat, Aegilops tauschii is considered to be a valuable resistance source to various biotic and abiotic stresses. However, little has been reported concerning the molecular mechanism of drought tolerance in Ae. tauschii. In this work, the drought tolerance of 155 Ae. tauschii accessions was firstly screened on the basis of their coleoptile lengths under simulated drought stress. Subsequently, two accessions (XJ002 and XJ098) with contrasting coleoptile lengths were selected and intensively analyzed on rate of water loss (RWL) as well as physiological characters, confirming the difference in drought tolerance at the seedling stage. Further, RNA-seq was utilized for global transcriptome profiling of the two accessions seedling leaves under drought stress conditions. A total of 6969 differentially expressed genes (DEGs) associated with drought tolerance were identified, and their functional annotations demonstrated that the stress response was mediated by pathways involving alpha-linolenic acid metabolism, starch and sucrose metabolism, peroxisome, mitogen-activated protein kinase (MAPK) signaling, carbon fixation in photosynthetic organisms, and glycerophospholipid metabolism. In addition, DEGs with obvious differences between the two accessions were intensively analyzed, indicating that the expression level of DEGs was basically in alignment with the physiological changes of Ae. tauschii under drought stress. The results not only shed fundamental light on the regulatory process of drought tolerance in Ae. tauschii, but also provide a new gene resource for improving the drought tolerance of common wheat.

Fine Mapping a Broad-Spectrum Powdery Mildew Resistance Gene in Chinese Landrace Datoumai, PmDTM, and Its Relationship with Pm24.

Powdery mildew, caused by the biotrophic fungal pathogen Blumeria graminis f. sp. tritici (Bgt), is a globally important wheat disease causing severe yield losses, and deployment of resistant varieties is the preferred choice for managing this disease. Chinese wheat landrace Datoumai was resistant to 22 of 23 Bgt isolates at the seedling stage. Genetic analysis based on the inoculation of Bgt isolate E09 on the F1, F2, and F2:3 populations derived from the cross Datoumai × Huixianhong revealed that the powdery mildew resistance of Datoumai is controlled by a single dominant gene, temporarily designated as PmDTM. Bulked segregant analysis and simple sequence repeat mapping with 200 F2 plants showed that PmDTM was located in the same genetic region as Pm24 on chromosome 1DS. To fine map PmDTM, 12 critical recombinants were identified from 1,192 F2 plants and delimited PmDTM to a 0.5-cM Xhnu58800 to Xhnu59000 interval covering 180.5 Kb (38,728,125 to 38,908,656 bp) on chromosome 1DS, and only one highly confident gene, TraesCS1D02G058900, was annotated within this region. TraesCS1D02G058900 encodes a receptor-like serine/threonine-protein kinase (STK), and a 6-bp deletion in exon 5 may confer the resistance to powdery mildew. Allele frequency analysis indicated that the STK allele with 6-bp deletion was only present in three landraces (Datoumai, Chiyacao [Pm24], and Hulutou) and was absent in all of the 353 Chinese modern cultivars and 147 foreign cultivars. These results demonstrate that PmDTM is mapped to the same locus as Pm24 and can be widely used to enhance powdery mildew resistance in wheat growing regions worldwide.

Transcriptome analysis reveals key genes involved in the regulation of nicotine biosynthesis at early time points after topping in tobacco (Nicotiana tabacum L.).

BACKGROUND: Nicotiana tabacum is an important economic crop. Topping, a common agricultural practice employed with flue-cured tobacco, is designed to increase leaf nicotine contents by increasing nicotine biosynthesis in roots. Many genes are found to be differentially expressed in response to topping, particularly genes involved in nicotine biosynthesis, but comprehensive analyses of early transcriptional responses induced by topping are not yet available. To develop a detailed understanding of the mechanisms regulating nicotine biosynthesis after topping, we have sequenced the transcriptomes of Nicotiana tabacum roots at seven time points following topping. RESULTS: Differential expression analysis revealed that 4830 genes responded to topping across all time points. Amongst these, nine gene families involved in nicotine biosynthesis and two gene families involved in nicotine transport showed significant changes during the immediate 24 h period following topping. No obvious preference to the parental species was detected in the differentially expressed genes (DEGs). Significant changes in transcript levels of nine genes involved in nicotine biosynthesis and phytohormone signal transduction were validated by qRT-PCR assays. 549 genes encoding transcription factors (TFs), found to exhibit significant changes in gene expression after topping, formed 15 clusters based on similarities of their transcript level time-course profiles. 336 DEGs involved in phytohormone signal transduction, including genes functionally related to the phytohormones jasmonic acid, abscisic acid, auxin, ethylene, and gibberellin, were identified at the earliest time point after topping. CONCLUSIONS: Our research provides the first detailed analysis of the early transcriptional responses to topping in N. tabacum, and identifies excellent candidates for further detailed studies concerning the regulation of nicotine biosynthesis in tobacco roots.

The genome of Populus alba x Populus tremula var. glandulosa clone 84K.

Poplar 84K (Populus alba x P. tremula var. glandulosa) is a fast-growing poplar hybrid. Originated in South Korea, this hybrid has been extensively cultivated in northern China. Due to the economic and ecological importance of this hybrid and high transformability, we now report the de novo sequencing and assembly of a male individual of poplar 84K using PacBio and Hi-C technologies. The final reference nuclear genome (747.5 Mb) has a contig N50 size of 1.99 Mb and a scaffold N50 size of 19.6 Mb. Complete chloroplast and mitochondrial genomes were also assembled from the sequencing data. Based on similarities to the genomes of P. alba var. pyramidalis and P. tremula, we were able to identify two subgenomes, representing 356 Mb from P. alba (subgenome A) and 354 Mb from P. tremula var. glandulosa (subgenome G). The phased assembly allowed us to detect the transcriptional bias between the two subgenomes, and we found that the subgenome from P. tremula displayed dominant expression in both 84K and another widely used hybrid, P. tremula x P. alba. This high-quality poplar 84K genome will be a valuable resource for poplar breeding and for molecular biology studies.

BZU2/ZmMUTE controls symmetrical division of guard mother cell and specifies neighbor cell fate in maize.

Intercellular communication in adjacent cell layers determines cell fate and polarity, thus orchestrating tissue specification and differentiation. Here we use the maize stomatal apparatus as a model to investigate cell fate determination. Mutations in ZmBZU2 (bizui2, bzu2) confer a complete absence of subsidiary cells (SCs) and normal guard cells (GCs), leading to failure of formation of mature stomatal complexes. Nuclear polarization and actin accumulation at the interface between subsidiary mother cells (SMCs) and guard mother cells (GMCs), an essential pre-requisite for asymmetric cell division, did not occur in Zmbzu2 mutants. ZmBZU2 encodes a basic helix-loop-helix (bHLH) transcription factor, which is an ortholog of AtMUTE in Arabidopsis (BZU2/ZmMUTE). We found that a number of genes implicated in stomatal development are transcriptionally regulated by BZU2/ZmMUTE. In particular, BZU2/ZmMUTE directly binds to the promoters of PAN1 and PAN2, two early regulators of protodermal cell fate and SMC polarization, consistent with the low levels of transcription of these genes observed in bzu2-1 mutants. BZU2/ZmMUTE has the cell-to-cell mobility characteristic similar to that of BdMUTE in Brachypodium distachyon. Unexpectedly, BZU2/ZmMUTE is expressed in GMC from the asymmetric division stage to the GMC division stage, and especially in the SMC establishment stage. Taken together, these data imply that BZU2/ZmMUTE is required for early events in SMC polarization and differentiation as well as for the last symmetrical division of GMCs to produce the two GCs, and is a master determinant of the cell fate of its neighbors through cell-to-cell communication.

Transcriptome-wide analysis of pseudouridylation of mRNA and non-coding RNAs in Arabidopsis.

Pseudouridine (Ψ) is widely distributed in mRNA and various non-coding RNAs in yeast and mammals, and the specificity of its distribution has been determined. However, knowledge about Ψs in the RNAs of plants, particularly in mRNA, is lacking. In this study, we performed genome-wide pseudouridine-sequencing in Arabidopsis and for the first time identified hundreds of Ψ sites in mRNA and multiple Ψ sites in non-coding RNAs. Many predicted and novel Ψ sites in rRNA and tRNA were detected. mRNA was extensively pseudouridylated, but with Ψs being under-represented in 3'-untranslated regions and enriched at position 1 of triple codons. The phenylalanine codon UUC was the most frequently pseudouridylated site. Some Ψs present in chloroplast 23S, 16S, and 4.5S rRNAs in wild-type Col-0 were absent in plants with a mutation of SVR1 (Suppressor of variegation 1), a chloroplast pseudouridine synthase gene. Many plastid ribosomal proteins and photosynthesis-related proteins were significantly reduced in svr1 relative to the wild-type, indicating the roles of SVR1 in chloroplast protein biosynthesis in Arabidopsis. Our results provide new insights into the occurrence of pseudouridine in Arabidopsis RNAs and the biological functions of SVR1, and will pave the way for further exploiting the mechanisms underlying Ψ modifications in controlling gene expression and protein biosynthesis in plants.

Molecular Characterization and Variation of the Celiac Disease Epitope Domains among α-Gliadin Genes in Aegilops tauschii.

To explore the distribution and quantity of toxic epitopes in α-gliadins from Aegilops tauschii, a total of 133 complete α-gliadin coding sequences were obtained, including 69 pseudogenes with at least one premature stop codon and 64 genes with complete open reading frames (ORFs). Plenty of deletions and single amino acid substitutions were found in the 4 celiac disease (CD) toxic epitope domains through multiple alignments, in which the sequence of DQ2.5-glia-α2 demonstrated the most significant changes. Interestingly, 7 of the 59 α-gliadins were free of any kind of intact CD toxic epitopes, providing potential gene resources for low CD toxicity breeding of common wheat. Analysis of the neighbor-joining tree demonstrates that 2 of the totally 7 α-gliadins cluster within the homologues of Triticum (A genome), and the other 5 group with those of Aegilops Sitopsis (B genome). This result implies that the 7 α-gliadin genes may be originated from the ancestor species of Ae. tauschii, evolved by the homoploid hybrid of Triticum and Aegilops Sitopsis. The remaining 52 α-gliadins form a separate clade from other homologues of A and B genomes, suggesting a recent rapid gene expansion by gene duplication associated with the species adaptation.