Transcriptome analysis of near-isogenic lines provides molecular insights into starch biosynthesis in maize kernel

Nitrogen deficiency tolerance conferred by introgression of a QTL derived from wild emmer into bread wheat

KEY MESSAGE

Genetic dissection of a QTL from wild emmer wheat, QGpc.huj.uh-5B.2, introgressed into bread wheat, identified candidate genes associated with tolerance to nitrogen deficiency, and potentially useful for improving nitrogen-use efficiency. Nitrogen (N) is an important macronutrient critical to wheat growth and development; its deficiency is one of the main factors causing reductions in grain yield and quality. N availability is significantly affected by drought or flooding, that are dependent on additional factors including soil type or duration and severity of stress. In a previous study, we identified a high grain protein content QTL (QGpc.huj.uh-5B.2) derived from the 5B chromosome of wild emmer wheat, that showed a higher proportion of explained variation under water-stress conditions. We hypothesized that this QTL is associated with tolerance to N deficiency as a possible mechanism underlying the higher effect under stress. To validate this hypothesis, we introgressed the QTL into the elite bread wheat var. Ruta, and showed that under N-deficient field conditions the introgression IL99 had a 33% increase in GPC (p < 0.05) compared to the recipient parent. Furthermore, evaluation of IL99 response to severe N deficiency (10% N) for 14 days, applied using a semi-hydroponic system under controlled conditions, confirmed its tolerance to N deficiency. Fine-mapping of the QTL resulted in 26 homozygous near-isogenic lines (BC4F5) segregating to N-deficiency tolerance. The QTL was delimited from - 28.28 to - 1.29 Mb and included 13 candidate genes, most associated with N-stress response, N transport, and abiotic stress responses. These genes may improve N-use efficiency under severely N-deficient environments. Our study demonstrates the importance of WEW as a source of novel candidate genes for sustainable improvement in tolerance to N deficiency in wheat.

Molecular mechanisms underlying low temperature inhibition of grain filling in maize (Zea mays L.): coordination of growth and cold responses

Low temperature (LT) greatly restricts grain filling in maize (Zea mays L.), but the relevant molecular mechanisms are not fully understood. To better understand the effect of LT on grain development, 17 hybrids were subjected to LT stress in field trials over 3 years, and two hybrids of them with contrasting LT responses were exposed to 30/20°C and 20/10°C for 7 days during grain filling in a greenhouse. At LT, thousand-kernel weight declined, especially in LT-sensitive hybrid FM985, while grain-filling rate was on average about 48% higher in LT-tolerant hybrid DK159 than FM985. LT reduced starch synthesis in kernel mainly by suppression of transcript levels and enzyme activities for sucrose synthase and hexokinase. Brassinolide (BR) was abundant in DK159 kernel, and genes involved in BR and cytokinin signals were inducible by stress. LT downregulated the genes in light-harvesting complex and photosystem I/II subunits, accompanied by reduced photosynthetic rate and Fv/Fm in ear leaf. The LT-tolerant hybrid could maintain a high soluble sugar content and fast interconversion between sucrose and hexose in the stem internode and cob, improving assimilate allocation to kernel at LT stress and paving the way for simultaneous growth and LT stress responses.

ZmDRR206 Regulates Nutrient Accumulation in Endosperm through Its Role in Cell Wall Biogenesis during Maize Kernel Development

Dirigent proteins (DIRs) contribute to plant fitness by dynamically reorganizing the cell wall and/or by generating defense compounds during plant growth, development, and interactions with environmental stresses. ZmDRR206 is a maize DIR, it plays a role in maintaining cell wall integrity during seedling growth and defense response in maize, but its role in regulating maize kernel development is unclear. Association analysis of candidate genes indicated that the natural variations of ZmDRR206 were significantly associated with maize hundred-kernel weight (HKW). ZmDRR206 plays a dominant role in storage nutrient accumulation in endosperm during maize kernel development, ZmDRR206 overexpression resulted in small and shrunken maize kernel with significantly reduced starch content and significantly decreased HKW. Cytological characterization of the developing maize kernels revealed that ZmDRR206 overexpression induced dysfunctional basal endosperm transfer layer (BETL) cells, which were shorter with less wall ingrowth, and defense response was constitutively activated in developing maize kernel at 15 and 18 DAP by ZmDRR206 overexpression. The BETL-development-related genes and auxin signal-related genes were down-regulated, while cell wall biogenesis-related genes were up-regulated in developing BETL of the ZmDRR206-overexpressing kernel. Moreover, the developing ZmDRR206-overexpressing kernel had significantly reduced contents of the cell wall components such as cellulose and acid soluble lignin. These results suggest that ZmDRR206 may play a regulatory role in coordinating cell development, storage nutrient metabolism, and stress responses during maize kernel development through its role in cell wall biogenesis and defense response, and provides new insights into understanding the mechanisms of kernel development in maize.

Transcriptome analysis reveals new insights in the starch biosynthesis of non-waxy and waxy broomcorn millet (Panicum miliaceum L.)

Broomcorn millet is a popular cereal with health benefits, and its grains are rich in starch. However, the differences in the pathway and key genes involved in starch biosynthesis of waxy and non-waxy broomcorn millet grain remain unclear. Therefore, the grain and starch physicochemical index and transcriptomic analyses of two genotypes of broomcorn millet were conducted at 3, 6, 9, 12, 15, 18, and 21 days after pollination. The phenotypic and physiological results indicated that the starch synthetic process of non-waxy and waxy broomcorn millet was significantly different. The amylose, amylopectin, and total starch contents of non-waxy broomcorn millet were 1.99, 4.74, and 6.73 mg/grain, while those of waxy broomcorn millet were 0.34, 5.94, and 6.28 mg/grain, respectively. The transcriptomic analysis revealed that 106 differentially expressed genes were identified, which were mainly enriched in the "amino sugar and nucleotide sugar metabolism", "pyruvate metabolism", "galactose metabolism", and "starch and sucrose metabolism" pathways. The WGCNA suggested that a total of 31 hub genes were correlated with starch biosynthesis. These findings provide a new approach to studying the starch synthesis in broomcorn millet.

Genetic dissection of QTLs for starch content in four maize DH populations

Starch is the principal carbohydrate source in maize kernels. Understanding the genetic basis of starch content (SC) benefits greatly in improving maize yield and optimizing end-use quality. Here, four double haploid (DH) populations were generated and were used to identify quantitative trait loci (QTLs) associated with SC. The phenotype of SC exhibited continuous and approximate normal distribution in each population. A total of 13 QTLs for SC in maize kernels was detected in a range of 3.65-16.18% of phenotypic variation explained (PVE). Among those, only some partly overlapped with QTLs previously known to be related to SC. Meanwhile, 12 genes involved in starch synthesis and metabolism located within QTLs were identified in this study. These QTLs will lay the foundation to explore candidate genes regulating SC in maize kernel and facilitate the application of molecular marker-assisted selection for a breeding program to cultivate maize varieties with a deal of grain quality.

Transcriptome analysis of near-isogenic lines for glume hairiness of wheat

Hairiness, which is a phenotypic trait common among land plants, primarily affects the stem, leaf, and floral organs. Plant hairiness is associated with complex functions. For example, glume hairiness in wheat is related to the resistance to biotic and abiotic stresses, and may also influence human health. In the present study, two pairs of near-isogenic lines (NILs) for glume hairiness, which were derived from a cross between a Tibetan semi-wild wheat accession (Triticum aestivum ssp. tibetanum Q1028) and a common wheat cultivar (T. aestivum 'Zhengmai 9023'), underwent a glume transcriptome analysis. We detected 27,935 novel genes, of which 18,027 were annotated. Additionally, 488 and 600 differentially expressed genes (DEGs) were detected in NIL1 and NIL2, respectively, with 37 DEGs detected in both NIL pairs. Moreover, 987 and 1584 single nucleotide polymorphisms (SNPs) were detected in NIL1 and NIL2, respectively, with 39 SNPs detected in both NIL pairs, of which most were located in the Hairy glume (Hg) gene region on chromosome arm 1AS. The annotation of the DEGs with gene ontology terms revealed that genes associated with hairiness in Arabidopsis and rice were similarly enriched. The possible functions of these genes related to glume hairiness were examined. The study results provide useful information for identifying candidate genes and the fine-mapping of Hg in the wheat genome.

Integrated transcriptome, small RNA, and degradome analysis reveals the complex network regulating starch biosynthesis in maize

Background

Starch biosynthesis in endosperm is a key process influencing grain yield and quality in maize. Although a number of starch biosynthetic genes have been well characterized, the mechanisms by which the expression of these genes is regulated, especially in regard to microRNAs (miRNAs), remain largely unclear.

Results

Sequence data for small RNAs, degradome, and transcriptome of maize endosperm at 15 and 25 d after pollination (DAP) from inbred lines Mo17 and Ji419, which exhibit distinct starch content and starch granule structure, revealed the mediation of starch biosynthetic pathways by miRNAs. Transcriptome analysis of these two lines indicated that 33 of 40 starch biosynthetic genes were differentially expressed, of which 12 were up-regulated in Ji419 at 15 DAP, one was up-regulated in Ji419 at 25 DAP, 14 were up-regulated in Ji419 at both 15 and 25 DAP, one was down-regulated in Ji419 at 15 DAP, two were down-regulated in Ji419 at 25 DAP, and three were up-regulated in Ji419 at 15 DAP and down-regulated in Ji419 at 25 DAP, compared with Mo17. Through combined analyses of small RNA and degradome sequences, 22 differentially expressed miRNAs were identified, including 14 known and eight previously unknown miRNAs that could target 35 genes. Furthermore, a complex co-expression regulatory network was constructed, in which 19 miRNAs could modulate starch biosynthesis in endosperm by tuning the expression of 19 target genes. Moreover, the potential operation of four miRNA-mediated pathways involving transcription factors, miR169a-NF-YA1-GBSSI/SSIIIa and miR169o-GATA9-SSIIIa/SBEIIb, was validated via analyses of expression pattern, transient transformation assays, and transactivation assays.

Conclusion

Our results suggest that miRNAs play a critical role in starch biosynthesis in endosperm, and that miRNA-mediated networks could modulate starch biosynthesis in this tissue. These results have provided important insights into the molecular mechanism of starch biosynthesis in developing maize endosperm.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5945-1) contains supplementary material, which is available to authorized users.

Maize transcriptomic repertoires respond to gibberellin stimulation

Phytohormone gibberellin (GA) serves as hub modulator of diverse biological events. Understanding the transcriptomic features of GA-mediated processes has scientific significance. The transcriptomic landscapes of cereal crops upon GA stimulation remains largely unknown. Herein, to reveal the transcriptomic changes in cereal crop maize under GA treatment, we first selected normal height and GA-sensitive maize dwarf plants from advanced backcross population for GA treatment. RNA-seq analysis discovered multiple protein-coding transcripts that were differentially expressed in GA-treated samples compared to distilled water-treated ones. Some differentially expressed transcripts, namely GA-responsive transcripts in this study, encoded the components of GA pathway, including CPS, KS, and KO enzymes for GA biosynthesis, GA2ox enzymes for GA degradation, DELLA repressors and GID1 receptor for GA signaling. A total of 214 shared GA-responsive transcripts were identified both in GA₃-treated normal height and GA-sensitive dwarf samples. Shared GA-responsive transcripts were involved in GA signaling, auxin biosynthesis, ethylene response, the composition and structure of cell wall, chlorophyll biogenesis, and sugar homeostasis. In addition, the convergence and divergence in expression of shared GA-responsive transcripts were observed in GA₃-treated normal height and GA-sensitive dwarf plants. Interaction network modeling indicated that some shared GA-responsive transcripts tended to be co-regulated, which increases the complexity of GA-triggered regulation at transcriptomic layer. Results presented here will extend our knowledge of GA-mediated regulatory cascade, and enhance our ability to apply hormone GA knowledge in agricultural practice.

Global transcriptome analysis uncovers the gene co-expression regulation network and key genes involved in grain development of wheat (Triticum aestivum L.)

Wheat grain development is a robust biological process that largely determines grain quality and yield. In this study, we investigated the grain transcriptome of winter wheat cv. Xiaoyan-6 at four developmental stages (5, 10, 15, and 20 days post-anthesis), using high-throughput RNA sequencing (RNA-Seq). We identified 427 grain-specific transcription factors (TFs) and 1653 differentially expressed TFs during grain development as well as a grain co-expression regulation network (GrainNet) of the TFs and their predicted co-expressed genes. Our study identified ten putative key TFs and the predicted regulatory genes of these TFs in wheat grain development of Xiaoyan-6. The analysis was given a firm basis through the study of additional wheat tissues, including root, stem, leaf, flag leaf, grain, spikes (from wheat plants at booting or heading stages) to provide a dataset of 92,478 high-confidence protein-coding genes that were mostly evenly distributed among subgenomes, but unevenly distributed across each of the chromosomes or each of the seven homeologous groups. Within this larger framework of the transcriptomes, we identified 4659 grain-specific genes (SEGs) and 26,500 differentially expressed genes (DEGs) throughout grain development stages tested. The SEGs identified mainly associate with regulation and signaling-related biological processes, while the DEGs mainly involve in cellular component organization or biogenesis and nutrient reservoir activity during grain development of Xiaoyan-6. This study establishes new targets for modifying genes related to grain development and yield, to fine-tune expression in different varieties and environments.

Candidate loci for the kernel row number in maize revealed by a combination of transcriptome analysis and regional association mapping

BACKGROUND

The kernel row number (KRN) of an ear is an important trait related to yield and domestication in maize. Exploring the underlying genetic mechanisms of KRN has great research significance and application potential.

RESULTS

In the present study, N531 with a KRN of 18-22 and SLN with a KRN of 4-6 were used as the recurrent parent and the donor parent, respectively, to develop two introgression lines (ILs), IL_A and IL_B, both of which have common negative-effect alleles from SLN on chromosomes 1, 5 and 10 and significantly reduced inflorescence meristem (IM) diameter and KRN compared with those of N531. We used RNA-Seq to investigate the transcriptome profiles of 5-mm immature ears of N531, IL_A and IL_B. We identified a total of 2872 differentially expressed genes (DEGs) between N531 and IL_A, 2428 DEGs between N531 and IL_B and 1811 DEGs between IL_A and IL_B. A total of 1252 DEGs were detected as overlapping DEGs, and 89 DEGs were located on the common introgression fragments. Furthermore, three DEGs (Zm00001d013277, Zm00001d015310 and Zm00001d015377) containing three SNPs associated with KRN were identified using regional association mapping.

CONCLUSIONS

These results will facilitate our understanding of ear development and provide important candidate genes for further study on KRN.

Seeds as oil factories

Studying seed oil metabolism. The seeds of higher plants represent valuable factories capable of converting photosynthetically derived sugars into a variety of storage compounds, including oils. Oils are the most energy-dense plant reserves and fatty acids composing these oils represent an excellent nutritional source. They supply humans with much of the calories and essential fatty acids required in their diet. These oils are then increasingly being utilized as renewable alternatives to petroleum for the chemical industry and for biofuels. Plant oils therefore represent a highly valuable agricultural commodity, the demand for which is increasing rapidly. Knowledge regarding seed oil production is extensively exploited in the frame of breeding programs and approaches of metabolic engineering for oilseed crop improvement. Complementary aspects of this research include (1) the study of carbon metabolism responsible for the conversion of photosynthetically derived sugars into precursors for fatty acid biosynthesis, (2) the identification and characterization of the enzymatic actors allowing the production of the wide set of fatty acid structures found in seed oils, and (3) the investigation of the complex biosynthetic pathways leading to the production of storage lipids (waxes, triacylglycerols). In this review, we outline the most recent developments in our understanding of the underlying biochemical and molecular mechanisms of seed oil production, focusing on fatty acids and oils that can have a significant impact on the emerging bioeconomy.

Plant Hexokinases are Multifaceted Proteins

Sugars are the main carbon and energy source in cells, but they can also act as signaling molecules that affect the whole plant life cycle. Certain tissues can produce sugars and supply them to others, and this plant tissue heterogeneity makes sugar signaling a highly complex process that requires elements capable of perceiving changes in sugar concentrations among different tissues, cell compartments and developmental stages. In plants, the regulatory effects of glucose (Glc) have been the most studied to date. The first Glc sensor identified in plants was hexokinase (HXK), which is currently recognized as a dual-function protein. In addition to its catalytic activity, this enzyme can also repress the expression of some photosynthetic genes in response to high internal Glc concentrations. Additionally, the catalytic activity of HXKs has a profound impact on cell metabolism and other sugar signaling pathways that depend on phosphorylated hexoses and intermediate glycolytic products. HXKs are the only proteins that are able to phosphorylate Glc in plants, since no evidence has been provided to date concerning the existence of a glucokinase. Moreover, the intracellular localization of HXKs seems to be crucial to their activity and sensor functions. Recently, two new and surprising functions have been described for HXKs. In this review, we discuss the versatility of HXKs in regard to their catalytic and glucose sensor activities, intracellular location, protein-protein and hormone interactions, as well as how these HXK characteristics influence plant growth and development, in an effort to understand this enzyme's role in improving plant productivity.

Inheritance and quantitative trail loci mapping of adventitious root numbers in cucumber seedlings under waterlogging conditions

The hypocotyl-derived adventitious root (AR) is an important morphological acclimation to waterlogging stress; however, its genetic basis has not been adequately understood. In the present study, a mixed major gene plus polygene inheritance model was used to analyze AR numbers (ARN) 7 days after waterlogging treatment in six generations (P₁, P₂, F₁, B₁, B_2, and F₂), using cucumber waterlogging tolerant line Zaoer-N and sensitive Pepino as parents. The results showed that the genetic model D-4, mixed one negative dominance major gene and additive-dominance polygenes, is the best-fitting genetic model for waterlogging-triggered ARN phenotype. A genetic linkage map spanning 550.8 cM and consisting of 149 simple sequence repeat (SSR) markers segregating into seven linkage groups was constructed. Three QTLs (ARN3.1, ARN5.1, and ARN6.1) distributed on chromosomes 3, 5, and 6 were identified by composite interval mapping. The major-effect QTL, ARN6.1, located between SSR12898 and SSR04751, was the only locus detected in three seasons, with least likelihood (LOD) scores of 8.8, 10.4, and 9.5 and account for 17.6, 24, and 19.8% of the phenotypic variance, respectively. Using five additional single nucleotide polymorphism (SNP) makers, the ARN6.1 was narrowed down to a 0.79 Mb interval franked by SSR12898 and SNP25558853. Illumina RNA-sequencing data generated on hypocotyls of two parents 48 h after waterlogging treatment revealed 15 genes in the 0.79 Mb interval were differentially expressed, including Csa6G503880 encoding a salicylic acid methyl transferase-like protein, Csa6G504590 encoding a cytochrome P450 monooxygenase, and Csa6G505230 encoding a heavy metal-associated protein. Our findings shed light on the genetic architecture underlying adventitious rooting during waterlogging stress in cucumber, and provide a list of potential gene targets for further elucidating waterlogging tolerance in plants.