Comparative cellular, physiological and transcriptome analyses reveal the potential easy dehulling mechanism of rice-tartary buckwheat (Fagopyrum Tararicum).

BACKGROUND: Tartary buckwheat has gained popularity in the food marketplace due to its abundant nutrients and high bioactive flavonoid content. However, its difficult dehulling process has severely restricted its food processing industry development. Rice-tartary buckwheat, a rare local variety, is very easily dehulled, but the cellular, physiological and molecular mechanisms responsible for this easy dehulling remains largely unclear. RESULTS: In this study, we integrated analyses of the comparative cellular, physiological, transcriptome, and gene coexpression network to insight into the reason that rice-tartary buckwheat is easy to dehull. Compared to normal tartary buckwheat, rice-tartary buckwheat has significantly brittler and thinner hull, and thinner cell wall in hull sclerenchyma cells. Furthermore, the cellulose, hemicellulose, and lignin contents of rice-tartary buckwheat hull were significantly lower than those in all or part of the tested normal tartary buckwheat cultivars, respectively, and the significant difference in cellulose and hemicellulose contents between rice-tartary buckwheat and normal tartary buckwheat began at 10 days after pollination (DAP). Comparative transcriptome analysis identified a total of 9250 differentially expressed genes (DEGs) between the rice- and normal-tartary buckwheat hulls at four different development stages. Weighted gene coexpression network analysis (WGCNA) of all DEGs identified a key module associated with the formation of the hull difference between rice- and normal-tartary buckwheat. In this specific module, many secondary cell wall (SCW) biosynthesis regulatory and structural genes, which involved in cellulose and hemicellulose biosynthesis, were identified as hub genes and displayed coexpression. These identified hub genes of SCW biosynthesis were significantly lower expression in rice-tartary buckwheat hull than in normal tartary buckwheat at the early hull development stages. Among them, the expression of 17 SCW biosynthesis relative-hub genes were further verified by quantitative real-time polymerase chain reaction (qRT-PCR). CONCLUSIONS: Our results showed that the lower expression of SCW biosynthesis regulatory and structural genes in rice-tartary buckwheat hull in the early development stages contributes to its easy dehulling by reducing the content of cell wall chemical components, which further effects the cell wall thickness of hull sclerenchyma cells, and hull thickness and mechanical strength.

Transcriptome Analysis Reveals Key Seed-Development Genes in Common Buckwheat (Fagopyrum esculentum).

Seed development is an essential and complex process, which is involved in seed size change and various nutrients accumulation, and determines crop yield and quality. Common buckwheat (Fagopyrum esculentum Moench) is a widely cultivated minor crop with excellent economic and nutritional value in temperate zones. However, little is known about the molecular mechanisms of seed development in common buckwheat (Fagopyrum esculentum). In this study, we performed RNA-Seq to investigate the transcriptional dynamics and identify the key genes involved in common buckwheat seed development at three different developmental stages. A total of 4619 differentially expressed genes (DEGs) were identified. Based on the results of Gene Ontology (GO) and KEGG analysis of DEGs, many key genes involved in the seed development, including the Ca2+ signal transduction pathway, the hormone signal transduction pathways, transcription factors (TFs), and starch biosynthesis-related genes, were identified. More importantly, 18 DEGs were identified as the key candidate genes for seed size through homologous query using the known seed size-related genes from different seed plants. Furthermore, 15 DEGs from these identified as the key genes of seed development were selected to confirm the validity of the data by using quantitative real-time PCR (qRT-PCR), and the results show high consistency with the RNA-Seq results. Taken together, our results revealed the underlying molecular mechanisms of common buckwheat seed development and could provide valuable information for further studies, especially for common buckwheat seed improvement.

Distribution and haplotype diversity of WKS resistance genes in wild emmer wheat natural populations.

KEY MESSAGE: The wheat stripe rust resistance gene Yr36 ( WKS1 ) with a unique kinase-START domain architecture is highly conserved in wild emmer wheat natural populations. Wild emmer wheat (Triticum dicoccoides) populations have developed various resistance strategies against the stripe rust pathogen Puccinia striiformis f. sp. tritici (Pst). The wild emmer gene, Yr36 (WKS1), which confers partial resistance to a broad spectrum of Pst races, is composed of a kinase and a START lipid-binding domain, a unique gene architecture found only in the Triticeae tribe. The analysis of 435 wild emmer accessions from a broad range of natural habitats revealed that WKS1 and its paralogue WKS2 are present only in the southern distribution range of wild emmer in the Fertile Crescent, supporting the idea that wheat domestication occurred in the northern populations. An analysis of full-length WKS1 sequence from 54 accessions identified 15 different haplotypes and very low-nucleotide diversity (π = 0.00019). The high level of WKS1 sequence conservation among wild emmer populations is in contrast to the high level of diversity previously observed in NB-LRR genes (e.g., Lr10 and Pm3). This phenomenon may reflect the different resistance mechanisms and different evolutionary pathways that shaped these genes, and may shed light on the evolution of genes that confer partial resistance to stripe rust. Only five WKS1 coding sequence haplotypes were revealed among all tested accessions, encoding four different putative WKS1 proteins (designated P0, P1, P2, and P3). Infection tests showed that P0, P1, and P3 haplotypes display a resistance response, while P2 displayed a susceptible response. These results show that the WKS1 proteins (P0, P1, and P3) can be useful to improve wheat resistance to stripe rust.

Transcriptome analysis revealed gene regulatory network involved in PEG-induced drought stress in Tartary buckwheat (Fagopyrum Tararicum).

Tartary buckwheat is a nutritious pseudo-cereal crop that is resistant to abiotic stresses, such as drought. However, the buckwheat's mechanisms for responding to drought stress remains unknown. We investigated the changes in physiology and gene expression under drought stress, which was simulated by treatment with polyethylene glycol (PEG). Five physiological indexes, namely MDA content, H2O2 content, CAT activity, SOD activity, and POD activity, were measured over time after 20% PEG treatment. All indexes showed dramatic changes in response to drought stress. A total of 1,190 differentially expressed genes (DEGs) were identified using RNA-seq and the most predominant were related to a number of stress-response genes and late embryogenesis abundant (LEA) proteins. DEGs were gathered into six clusters and were found to be involved in the ABA biosynthesis and signal pathway based on hierarchical clustering and GO and KEGG pathway enrichment. Transcription factors, such as NAC and bZIP, also took part in the response to drought stress. We determined an ABA-dependent and ABA-independent pathway in the regulation of drought stress in Tartary buckwheat. To the best of our knowledge, this is the first transcriptome analysis of drought stress in Tartary buckwheat, and our results provide a comprehensive gene regulatory network of this crop in response to drought stress.

Molecular characterization of two silenced y-type genes for Glu-B1 in Triticum aestivum ssp. yunnanese and ssp. tibetanum.

The high molecular weight glutenin subunits (HMW-GSs) are a major class of common wheat storage proteins. The bread-making quality of common wheat flour is influenced by the composition of HMW-GSs. In the present study, two unexpressed 1By genes from Triticum aesitvum L.ssp.yunnanese AS332 and T. aesitvum ssp.tibetanum AS908 were respectively cloned and characterized. The results indicated that both of the silenced 1By genes in AS332 and AS908 were 1By9. In contrast to previously reported mechanisms for silenced genes 1Ax and 1Ay, which was due to the insertion of transposon elements or the presence of premature stop codon via base substitution of C-->T transition in trinucleotides CAA or CAG, the silence of 1By9 genes was caused by premature stop codons via the deletion of base A in trinucleotide CAA, which lead to frameshift mutation and indirectly produced several premature stop codons (TAG) downstream of the coding sequence.

Global transcriptome analysis and identification of genes involved in nutrients accumulation during seed development of rice tartary buckwheat (Fagopyrum Tararicum).

Tartary buckwheat seeds are rich in various nutrients, such as storage proteins, starch, and flavonoids. To get a good knowledge of the transcriptome dynamics and gene regulatory mechanism during the process of seed development and nutrients accumulation, we performed a comprehensive global transcriptome analysis using rice tartary buckwheat seeds at different development stages, namely pre-filling stage, filling stage, and mature stage. 24 819 expressed genes, including 108 specifically expressed genes, and 11 676 differentially expressed genes (DEGs) were identified. qRT-PCR analysis was performed on 34 DEGs to validate the transcriptome data, and a good consistence was obtained. Based on their expression patterns, the identified DEGs were classified to eight clusters, and the enriched GO items in each cluster were analyzed. In addition, 633 DEGs related to plant hormones were identified. Furthermore, genes in the biosynthesis pathway of nutrients accumulation were analyzed, including 10, 20, and 23 DEGs corresponding to the biosynthesis of seed storage proteins, flavonoids, and starch, respectively. This is the first transcriptome analysis during seed development of tartary buckwheat. It would provide us a comprehensive understanding of the complex transcriptome dynamics during seed development and gene regulatory mechanism of nutrients accumulation.

Evaluation of assembly strategies using RNA-seq data associated with grain development of wheat (Triticum aestivum L.).

Wheat (Triticum aestivum L.) is one of the most important crops cultivated worldwide. Identifying the complete transcriptome of wheat grain could serve as foundation for further study of wheat seed development. However, the relatively large size and the polyploid complexity of the genome have been substantial barriers to molecular genetics and transcriptome analysis of wheat. Alternatively, RNA sequencing has provided some useful information about wheat genes. However, because of the large number of short reads generated by RNA sequencing, factors that are crucial to transcriptome assembly, including software, candidate parameters and assembly strategies, need to be optimized and evaluated for wheat data. In the present study, four cDNA libraries associated with wheat grain development were constructed and sequenced. A total of 14.17 Gb of high-quality reads were obtained and used to assess different assembly strategies. The most successful approach was to filter the reads with Q30 prior to de novo assembly using Trinity, merge the assembled contigs with genes available in wheat cDNA reference data sets, and combine the resulting assembly with an assembly from a reference-based strategy. Using this approach, a relatively accurate and nearly complete transcriptome associated with wheat grain development was obtained, suggesting that this is an effective strategy for generation of a high-quality transcriptome from RNA sequencing data.

Production of aneuhaploid and euhaploid sporocytes by meiotic restitution in fertile hybrids between durum wheat Langdon chromosome substitution lines and Aegilops tauschii.

Fertile F(1) hybrids were obtained between durum wheat (Triticum durum Desf.) Langdon (LDN) and its 10 disomic substitution (LDN DS) lines with Aegilops tauschii accession AS60 without embryo rescue. Selfed seedset rates for hybrids of LDN with AS60 were 36.87% and 49.45% in 2005 and 2006, respectively. Similar or higher selfed seedset rates were observed in the hybrids of 1D (1A), 1D (1B), 3D (3A), 4D (4B), 7D (7A), and 2D (2B) with AS60, while lower in hybrids of 3D (3B) + 3BL, 5D (5A) + 5AL, 5D (5B) + 5B and 6D (6B) + 6BS with AS60 compared with the hybrids of LDN with AS60. Observation of male gametogenesis showed that meiotic restitution, both first-division restitution (FDR) and single-division meiosis (SDM) resulted in the formation of functional unreduced gametes, which in turn produced seeds. Both euhaploid and aneuhaploid gametes were produced in F(1) hybrids. This suggested a strategy to simultaneously transfer and locate major genes from the ancestral species T. turgidum or Ae. tauschii. Moreover, there was no significant difference in the aneuhaploid rates between the F(1) hybrids of LDN and LDN DS lines with AS60, suggesting that meiotic pairing between the two D chromosomes in the hybrids of LDN DS lines with AS60 did not promote the formation of aneuhaploid gametes.

Comparison of newly synthetic hexaploid wheat with its donors on SSR products.

Microsatellites or SSRs as powerful genetic markers have widely been used in genetics and evolutionary biology in common wheat. Because of the high polymorphism, newly synthesized hexaploid wheat has been used in the construction of genetic segregation population for SSR markers. However, data on the evolution of microsatellites during the polyploidization event of hexaploid wheat are limited. In this study, 66 pairs of specific to A/B genome SSR patterns among newly synthesized hexaploid wheat, the donor tetraploid wheat and Aegilops tauschii were compared. The results indicated that most SSR markers were conserved during the polyploidization events of newly synthetic hexaploid wheat, from Triticum turgidum and Ae. tauschii. Over 70% A/B genome specific SSR markers could amplify the SSR sequences from the D genome of Ae. tauschii. Most amplified fragments from Ae. tauschii were detected in synthetic hexaploid at corresponding positions with the same sizes and patterns as in its parental Ae. tauschii. This suggested that these SSR markers, specific for A/B genome in common wheat, could amplify SSR products of D genome besides A/B genome in the newly synthesized hexaploid wheat, that is, these SSR primers specific for A/B genome in common wheat were nonspecific for the A/B genome in the synthetic hexaploid wheat. In addition, one amplified Ae. tauschii product was not detected in the newly synthetic hexaploid wheat. An extra-amplified product was found in the newly synthetic hexaploid wheat. These results suggested that caution should be taken when using SSR marker to genotype newly synthetic hexaploid wheat.