Transcriptional Dynamics of Grain Development in Barley (Hordeum vulgare L.)

The transcriptome landscape of developing barley seeds

Cereal grains are an important source of food and feed. To provide comprehensive spatiotemporal information about biological processes in developing seeds of cultivated barley (Hordeum vulgare L. subsp. vulgare), we performed a transcriptomic study of the embryo, endosperm, and seed maternal tissues collected from grains 4-32 days after pollination. Weighted gene co-expression network and motif enrichment analyses identified specific groups of genes and transcription factors (TFs) potentially regulating barley seed tissue development. We defined a set of tissue-specific marker genes and families of TFs for functional studies of the pathways controlling barley grain development. Assessing selected groups of chromatin regulators revealed that epigenetic processes are highly dynamic and likely play a major role during barley endosperm development. The repressive H3K27me3 modification is globally reduced in endosperm tissues and at specific genes related to development and storage compounds. Altogether, this atlas uncovers the complexity of developmentally regulated gene expression in developing barley grains.

Transcriptome Landscape Analyses of the Regulatory Network for Zygotic Embryo Development in Paeonia ostii

Paeonia ostii is a worldwide ornamental flower and an emerging oil crop. Zyotic embryogenesis is a critical process during seed development, and it can provide a basis for improving the efficiency of somatic embryogenesis (SE). In this study, transcriptome sequencing of embryo development was performed to investigate gene expression profiling in P. ostii and identified Differentially expressed genes (DEGs) related to transcription factors, plant hormones, and antioxidant enzymes. The results indicated that IAA (Indole-3-acetic acid), GA (Gibberellin), BR (Brassinosteroid) and ETH (Ethylene) were beneficial to early embryonic morphogenesis, while CTK (Cytokinin) and ABA (Abscisic Acid) promoted embryo morphogenesis and maturation. The antioxidant enzymes' activity was the highest in early embryos and an important participant in embryo formation. The high expression of the genes encoding fatty acid desaturase was beneficial to fast oil accumulation. Representative DEGs were selected and validated using qRT-PCR. Protein-protein interaction network (PPI) was predicted, and six central node proteins, including AUX1, PIN1, ARF6, LAX3, ABCB19, PIF3, and PIF4, were screened. Our results provided new insights into the formation of embryo development and even somatic embryo development in tree peonies.

Identification of the genes associated with β-glucan synthesis and accumulation during grain development in barley

The presence of β-glucan in barley grains is one of its important quality traits. Lower β-glucan content is required for the barley used in beer and feed production, while higher β-glucan content is beneficial for food barley. Although intensive research has been carried out on the genotypic and environmental differences in β-glucan content in barley grains, little information is available on the molecular mechanisms underlying their genotypic differences and genetic regulation of β-glucan synthesis and accumulation. In this study, RNA sequencing analysis was conducted to compare the transcriptome profiles of two barley genotypes (BCS192 and BCS297) that greatly differ in grain β-glucan content, in order to identify the key genes responsible for β-glucan synthesis and accumulation during grain development. The results showed that carbohydrate metabolic processes and starch and sucrose metabolism play significant roles in β-glucan synthesis. The identified differently expressed genes (DEGs), which are closely associated with grain β-glucan content, are mainly involved in hydrolase activity and glucan metabolic processes. In addition, β-glucan accumulation in barley grains is predominantly regulated by photosynthesis and carbon metabolism. The DEGs identified in this study and their functions may provide new insights into the molecular mechanisms of β-glucan synthesis and genotypic differences in barley grains.

Integrative metabolomic and transcriptomic analyses reveal the mechanisms of Tibetan hulless barley grain coloration

Cereal grains accumulate anthocyanin during developmental process. The anthocyanin content increases at grain filling stages to develop grain coloration in cereals. However, anthocyanin biosynthesis responsible for grain coloring and its regulatory mechanisms controlled by structural and functional genes remain unclear. Therefore, this study aimed to explore the global map of metabolic changes linked to grain coloration of Tibetan hulless barley (qingke) using an integrative metabolome and transcriptome approach. Grains from three colored qingke cultivars at different developmental stages were considered for molecular and metabolic investigations. A total of 120 differentially accumulated metabolites (DAMs) and 8,327 differentially expressed genes (DEGs) were filtered. DEGs were mainly enriched in the phenylpropanoid and flavonoid pathways. The transcript levels of anthocyanin biosynthesis genes (PAL, C4H, 4CL, CHS, FLS, F3H, F3'H, DFR, ANS, GT, OMT, and MAT) significantly upregulate in colored qingke compared to the non-colored variety. During grain development and maturation, the strong correlation of HvMYC2 expression with anthocyanin contents and anthocyanin biosynthesis genes suggested it as a critical gene in anthocyanin accumulation. Further results confirmed that HvMYC2 could be activated by HvMYB and be a positive regulator of UV-B and cold tolerance in qingke. In addition, verification based on enzymatic assays indicated that six key modifier enzymes could catalyze glycosylation, malonylation, and methylation of anthocyanins, thereby dissecting the major anthocyanin modification pathway in colored qingke. Overall, our study provides global insight into anthocyanin accumulation and the mechanism underlying grain coloration in qingke.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Genome-Wide Identification and Transcriptional Expression Profiles of PP2C in the Barley (Hordeum vulgare L.) Pan-Genome

The gene family protein phosphatase 2C (PP2C) is related to developmental processes and stress responses in plants. Barley (Hordeum vulgare L.) is a popular cereal crop that is primarily utilized for human consumption and nutrition. However, there is little knowledge regarding the PP2C gene family in barley. In this study, a total of 1635 PP2C genes were identified in 20 barley pan-genome accessions. Then, chromosome localization, physical and chemical feature predictions and subcellular localization were systematically analyzed. One wild barley accession (B1K-04-12) and one cultivated barley (Morex) were chosen as representatives to further analyze and compare the differences in HvPP2Cs between wild and cultivated barley. Phylogenetic analysis showed that these HvPP2Cs were divided into 12 subgroups. Additionally, gene structure, conserved domain and motif, gene duplication event detection, interaction networks and gene expression profiles were analyzed in accessions Morex and B1K-04-12. In addition, qRT-PCR experiments in Morex indicated that seven HvMorexPP2C genes were involved in the response to aluminum and low pH stresses. Finally, a series of positively selected homologous genes were identified between wild accession B1K-04-12 and another 14 cultivated materials, indicating that these genes are important during barley domestication. This work provides a global overview of the putative physiological and biological functions of PP2C genes in barley. We provide a broad framework for understanding the domestication- and evolutionary-induced changes in PP2C genes between wild and cultivated barley.

Genome-wide investigation and functional analysis of RNA editing sites in wheat

Wheat is an important cereal and half of the world population consumed it. Wheat faces environmental stresses and different techniques (CRISPR, gene silencing, GWAS, etc.) were used to enhance its production but RNA editing (RESs) is not fully explored in wheat. RNA editing has a special role in controlling environmental stresses. The genome-wide identification and functional characterization of RESs in different types of wheat genotypes was done. We employed six wheat genotypes by RNA-seq analyses to achieve RESs. The findings revealed that RNA editing events occurred on all chromosomes equally. RNA editing sites were distributed randomly and 10–12 types of RESs were detected in wheat genotypes. Higher number of RESs were detected in drought-tolerant genotypes. A-to-I RNA editing (2952, 2977, 1916, 2576, 3422, and 3459) sites were also identified in six wheat genotypes. Most of the genes were found to be engaged in molecular processes after a Gene Ontology analysis. PPR (pentatricopeptide repeat), OZ1 (organelle zinc-finger), and MORF/RIP gene expression levels in wheat were also examined. Normal growth conditions diverge gene expression of these three different gene families, implying that normal growth conditions for various genotypes can modify RNA editing events and have an impact on gene expression levels. While the expression of PPR genes was not change. We used Variant Effect Predictor (VEP) to annotate RNA editing sites, and Local White had the highest RESs in the CDS region of the protein. These findings will be useful for prediction of RESs in other crops and will be helpful in drought tolerance development in wheat.

Genome-wide identification of seed storage protein gene regulators in wheat through coexpression analysis

Only a few transcriptional regulators of seed storage protein (SSP) genes have been identified in common wheat (Triticum aestivum L.). Coexpression analysis could be an efficient approach to characterize novel transcriptional regulators at the genome-scale considering the correlated expression between transcriptional regulators and target genes. As the A genome donor of common wheat, Triticum urartu is more suitable for coexpression analysis than common wheat considering the diploid genome and single gene copy. In this work, the transcriptome dynamics in endosperm of T. urartu throughout grain filling were revealed by RNA-Seq analysis. In the coexpression analysis, a total of 71 transcription factors (TFs) from 23 families were found to be coexpressed with SSP genes. Among these TFs, TuNAC77 enhanced the transcription of SSP genes by binding to cis-elements distributed in promoters. The homolog of TuNAC77 in common wheat, TaNAC77, shared an identical function, and the total SSPs were reduced by about 24% in common wheat when TaNAC77 was knocked down. This is the first genome-wide identification of transcriptional regulators of SSP genes in wheat, and the newly characterized transcriptional regulators will undoubtedly expand our knowledge of the transcriptional regulation of SSP synthesis.

Short Time-Series Expression Transcriptome Data Reveal the Gene Expression Patterns of Dairy Cow Mammary Gland as Milk Yield Decreased Process

The existing research on dairy cow mammary gland genes is extensive, but there have been few reports about dynamic changes in dairy cow mammary gland genes as milk yield decrease. For the first time, transcriptome analysis based on short time-series expression miner (STEM) and histological observations were performed using the Holstein dairy cow mammary gland to explore gene expression patterns in this process of decrease (at peak, mid-, and late lactation). Histological observations suggested that the number of mammary acinous cells at peak/mid-lactation was significantly higher than that at mid-/late lactation, and the lipid droplets area secreted by dairy cows was almost unaltered across the three stages of lactation (p > 0.05). Totals of 882 and 1439 genes were differentially expressed at mid- and late lactation, respectively, compared to peak lactation. Function analysis showed that differentially expressed genes (DEGs) were mainly related to apoptosis and energy metabolism (fold change ≥ 2 or fold change ≤ 0.5, p-value ≤ 0.05). Transcriptome analysis based on STEM identified 16 profiles of differential gene expression patterns, including 5 significant profiles (false discovery rate, FDR ≤ 0.05). Function analysis revealed DEGs involved in milk fat synthesis were downregulated in Profile 0 and DEGs in Profile 12 associated with protein synthesis. These findings provide a foundation for future studies on the molecular mechanisms underlying mammary gland development in dairy cows.

Genome-Wide Identification of Barley Long Noncoding RNAs and Analysis of Their Regulatory Interactions during Shoot and Grain Development

Long noncoding RNAs (lncRNAs) are a class of RNA molecules with gene regulatory functions in plant development and the stress response. Although the number of lncRNAs identified in plants is rapidly increasing, very little is known about their role in barley development. In this study, we performed global identification of barley lncRNAs based on 53 RNAseq libraries derived from nine different barley tissues and organs. In total, 17,250 lncRNAs derived from 10,883 loci were identified, including 8954 novel lncRNAs. Differential expression of lncRNAs was observed in the developing shoot apices and grains, the two organs that have a direct influence on the final yield. The regulatory interaction of differentially expressed lncRNAs with the potential target genes was evaluated. We identified 176 cis-acting lncRNAs in shoot apices and 424 in grains, while the number of trans-acting lncRNAs in these organs was 1736 and 540, respectively. The potential target protein-coding genes were identified, and their biological function was annotated using MapMan ontology. This is the first insight into the roles of lncRNAs in barley development on the genome-wide scale, and our results provide a solid background for future functional studies.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

Quantitative Proteome Profiling Provides Insight into the Proteins Associated with β-Glucan Accumulation in Hull-less Barley Grains

The hull-less barley (Qingke) is widely planted as a staple food crop in the Tibetan area, China, and the grains contains high content of β-glucan (BG). The mechanisms of BG synthesis and accumulation in qingke has not been studied at the protein level. This study characterized the proteins associated with BG synthesis and accumulation during qingke seed development. The proteome profiles of qingke seeds taken at 20, 30, and 40 days after flowering were compared using the TMT-based quantitative proteomics. A total of 4283 proteins were identified, with 759 being differentially expressed (DEPs) throughout seed development. Comparisons of protein expression pattern, functions, and pathway enrichment tests highlight cell wall modification, carbon and energy metabolism, polysaccharide metabolism, post-transcriptional modifications, and vesicular transport as critical biological processes related to qingke BG accumulation. Furthermore, induction of starch synthase, starch branching enzyme, pectin acetyl esterases, beta-glucosidases, beta-amylases, 1,4-beta-xylan, xyloglucan, α-amylase inhibitors, and glycosyltransferases underpinned BG synthesis. The results also indicated that the proteins involved in glycolytic, gluconeogenesis, and glyoxylate bypass pathways provided energy and reducing power for BG storage. Parallel reaction monitoring (PRM) and quantitative real-time PCR (qPCR) analyses confirmed the expression profile of the proteins obtained by TMT-based proteomics. The current results provided an insight into the mechanisms of BG synthesis and accumulation during qingke seed development.

The Impact of Zinc Oxide Nanoparticles on Cytotoxicity, Genotoxicity, and miRNA Expression in Barley (Hordeum vulgare L.) Seedlings

Zinc oxide nanoparticles are one of the most commonly engineered nanomaterials and necessarily enter the environment because of the large quantities produced and their widespread application. Understanding the impacts of nanoparticles on plant growth and development is crucial for the assessment of probable environmental risks to food safety and human health, because plants are a fundamental living component of the ecosystem and the most important source in the human food chain. The objective of this study was to examine the impact of different concentrations of zinc oxide nanoparticles on barley Hordeum vulgare L. seed germination, seedling morphology, root cell viability, stress level, genotoxicity, and expression of miRNAs. The results demonstrate that zinc oxide nanoparticles enhance barley seed germination, shoot/root elongation, and H₂O₂ stress level and decrease root cell viability and genomic template stability and up- and downregulated miRNAs in barley seedlings.

Genome Wide Analysis of the Transcriptional Profiles in Different Regions of the Developing Rice Grains

BACKGROUND

Rice is an important food source for humans worldwide. Because of its nutritional and agricultural significance, a number of studies addressed various aspects of rice grain development and grain filling. Nevertheless, the molecular processes underlying grain filling and development, and in particular the contributions of different grain tissues to these processes, are not understood.

MAIN TEXT

Using RNA-sequencing, we profiled gene expression activity in grain tissues comprised of cross cells (CC), the nucellar epidermis (NE), ovular vascular trace (OVT), endosperm (EN) and the aleurone layer (AL). These tissues were dissected using laser capture microdissection (LCM) at three distinct grain development stages. The mRNA expression datasets offer comprehensive and new insights into the gene expression patterns in different rice grain tissues and their contributions to grain development. Comparative analysis of the different tissues revealed their similar and/or unique functions, as well as the spatio-temporal regulation of common and tissue-specific genes. The expression patterns of genes encoding hormones and transporters indicate an important role of the OVT tissue in metabolite transport during grain development. Gene co-expression network prediction on OVT-specific genes identified several distinct and common development-specific transcription factors. Further analysis of enriched DNA sequence motifs proximal to OVT-specific genes revealed known and novel DNA sequence motifs relevant to rice grain development.

CONCLUSION

Together, the dataset of gene expression in rice grain tissues is a novel and useful resource for further work to dissect the molecular and metabolic processes during rice grain development.

DOF, MYB and TCP transcription factors: Their possible roles on barley germination and seedling establishment

Seed germination is a multi-staged complex process during seed plant life cycle, and it is tightly regulated through a coordinated expression of diverse genes in diverse tissues. As regulatory molecules of gene expression, determination of transcription factors is crucial to understanding molecular basis and regulatory network of germination process and seedling establishment. However, limited data on the contributions of these transcription factors to the germination of crop barley (Hordeum vulgare L.) are available. Here, we investigated the expression profiles of selected transcription factors from different families (DOF, MYB and TCP) with qRT-PCR analysis in various tissues including coleoptiles, leaves and roots following the germination. Analysis of MYB and DOF gene expression profiles indicated that there were differing expressions in different aged tissues, HvMYB5 and HvDOF2 being the most outstanding one in the oldest tissue, 15-day-old root. On the other hand, investigated TCP genes were lowly expressed compared to selected MYB and DOF genes, except HvTCP3, where the highest expression was observed in 15-day-old root tissue. The obtained expression profiles illustrate the importance of potential regulatory roles of transcription factors in early developmental stages of barley germination and seedling establishment.

Heritable temporal gene expression patterns correlate with metabolomic seed content in developing hexaploid oat seed

Oat ranks sixth in world cereal production and has a higher content of health-promoting compounds compared with other cereals. However, there is neither a robust oat reference genome nor transcriptome. Using deeply sequenced full-length mRNA libraries of oat cultivar Ogle-C, a de novo high-quality and comprehensive oat seed transcriptome was assembled. With this reference transcriptome and QuantSeq 3' mRNA sequencing, gene expression was quantified during seed development from 22 diverse lines across six time points. Transcript expression showed higher correlations between adjacent time points. Based on differentially expressed genes, we identified 22 major temporal co-expression (TCoE) patterns of gene expression and revealed enriched gene ontology biological processes. Within each TCoE set, highly correlated transcripts, putatively commonly affected by genetic background, were clustered and termed genetic co-expression (GCoE) sets. Seventeen of the 22 TCoE sets had GCoE sets with median heritabilities higher than 0.50, and these heritability estimates were much higher than that estimated from permutation analysis, with no divergence observed in cluster sizes between permutation and non-permutation analyses. Linear regression between 634 metabolites from mature seeds and the PC1 score of each of the GCoE sets showed significantly lower p-values than permutation analysis. Temporal expression patterns of oat avenanthramides and lipid biosynthetic genes were concordant with previous studies of avenanthramide biosynthetic enzyme activity and lipid accumulation. This study expands our understanding of physiological processes that occur during oat seed maturation and provides plant breeders the means to change oat seed composition through targeted manipulation of key pathways.

Harnessing Transcription Factors as Potential Tools to Enhance Grain Size Under Stressful Abiotic Conditions in Cereal Crops

Predicted climate change is widely cited to significantly reduce yields of the major cereal crop species in a period where demand is rapidly rising due to a growing global population. This requires exhaustive research to develop genetic resources in order to address the expected production deficiencies which will largely be driven by abiotic stress. Modification of multiple genes is an approach that can address the predicted challenges; however, it is time-consuming and costly to modify multiple genes simultaneously. Transcription factors represent a group of proteins regulating multiple genes simultaneously and are therefore promising targets to concurrently improve multiple traits concurrently, such as abiotic stress tolerance and grain size (a contributor to yield). Many studies have identified the complex role that transcription factors of multiple families have contributed toward abiotic stress tolerance or grain size, although research addressing both simultaneously is in its infancy despite its potential significance for cereal crop improvement. Here we discuss the potential role that transcription factors may contribute toward improving cereal crop productivity under adverse environmental conditions and offer research objectives that need to be addressed before the modification of transcription factors becomes routinely used to positively manipulate multiple target traits.

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 Ca²⁺ 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.