Dynamic transcriptome landscape of maize embryo and endosperm development

Opaque-2 Regulates a Complex Gene Network Associated with Cell Differentiation and Storage Functions of Maize Endosperm

Development of the cereal endosperm involves cell differentiation processes that enable nutrient uptake from the maternal plant, accumulation of storage products, and their utilization during germination. However, little is known about the regulatory mechanisms that link cell differentiation processes with those controlling storage product synthesis and deposition, including the activation of zein genes by the maize (Zea mays) bZIP transcription factor Opaque-2 (O2). Here, we mapped in vivo binding sites of O2 in B73 endosperm and compared the results with genes differentially expressed in B73 and B73o2 We identified 186 putative direct O2 targets and 1677 indirect targets, encoding a broad set of gene functionalities. Examination of the temporal expression patterns of O2 targets revealed at least two distinct modes of O2-mediated gene activation. Two O2-activated genes, bZIP17 and NAKED ENDOSPERM2 (NKD2), encode transcription factors, which can in turn coactivate other O2 network genes with O2. NKD2 (with its paralog NKD1) was previously shown to be involved in regulation of aleurone development. Collectively, our results provide insights into the complexity of the O2-regulated network and its role in regulation of endosperm cell differentiation and function.

Maize Oxalyl-CoA Decarboxylase1 Degrades Oxalate and Affects the Seed Metabolome and Nutritional Quality

The organic acid oxalate occurs in microbes, animals, and plants; however, excessive oxalate accumulation in vivo is toxic to cell growth and decreases the nutritional quality of certain vegetables. However, the enzymes and functions required for oxalate degradation in plants remain largely unknown. Here, we report the cloning of a maize (Zea mays) opaque endosperm mutant that encodes oxalyl-CoA decarboxylase1 (EC4.1.1.8; OCD1). Ocd1 is generally expressed and is specifically induced by oxalate. The ocd1 mutant seeds contain a significantly higher level of oxalate than the wild type, indicating that the ocd1 mutants have a defect in oxalate catabolism. The maize classic mutant opaque7 (o7) was initially cloned for its high lysine trait, although the gene function was not understood until its homolog in Arabidopsis thaliana was found to encode an oxalyl-CoA synthetase (EC 6.2.1.8), which ligates oxalate and CoA to form oxalyl-CoA. Our enzymatic analysis showed that ZmOCD1 catalyzes oxalyl-CoA, the product of O7, into formyl-CoA and CO₂ for degradation. Mutations in ocd1 caused dramatic alterations in the metabolome in the endosperm. Our findings demonstrate that ZmOCD1 acts downstream of O7 in oxalate degradation and affects endosperm development, the metabolome, and nutritional quality in maize seeds.

AgriSeqDB: an online RNA-Seq database for functional studies of agriculturally relevant plant species

BACKGROUND

The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops.

DESCRIPTION

AgriSeqDB ( https://expression.latrobe.edu.au/agriseqdb ) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications.

CONCLUSION

Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Maize opaque mutants are no longer so opaque

The endosperm of angiosperms is a zygotic seed organ that stores nutrient reserves to support embryogenesis and seed germination. Cereal endosperm is also a major source of human calories and an industrial feedstock. Maize opaque endosperm mutants commonly exhibit opaque, floury kernels, along with other abnormal seed and/or non-seed phenotypes. The opaque endosperm phenotype is sometimes accompanied by a soft kernel texture and increased nutritional quality, including a higher lysine content, which are valuable agronomic traits that have drawn attention of maize breeders. Recently, an increasing number of genes that underlie opaque mutants have been cloned, and their characterization has begun to shed light on the molecular basis of the opaque endosperm phenotype. These mutants are categorized by disruption of genes encoding zein or non-zein proteins localized to protein bodies, enzymes involved in endosperm metabolic processes, or transcriptional regulatory proteins associated with endosperm storage programs.

Oil biosynthesis and transcriptome profiles in developing endosperm and oil characteristic analyses in Paeonia ostii var. lishizhenii

Paeonia ostii var. lishizhenii, a well-known medicinal and horticultural plant, is indigenous to China. Recent studies have shown that its seed has a high oil content, and it was approved as a novel resource of edible oil with a high level of α-linolenic acid by the Chinese Government. This study measured the seed oil contents and fatty acid components of P. ostii var. lishizhenii and six other peonies, P. suffruticosa, P. ludlowii, P. decomposita, P. rockii, and P. lactiflora Pall. 'Heze' and 'Gansu'. The results show that P. ostii var. lishizhenii exhibits the average oil characteristics of tested peonies, with an oil content of 21.3%, α-linolenic acid 43.8%, and unsaturated fatty acids around 92.1%. Hygiene indicators for the seven peony seed oils met the Chinese national food standards. P. ostii var. lishizhenii seeds were used to analyze transcriptome gene regulation networks on endosperm development and oil biosynthesis. In total, 124,117 transcripts were obtained from six endosperm developing stages (S0-S5). The significant changes in differential expression genes (DEGs) clarify three peony endosperm developmental phases: the endosperm cell mitotic phase (S0-S1), the TAG biosynthesis phase (S1-S4), and the mature phase (S5). The DEGs in plant hormone signal transduction, DNA replication, cell division, differentiation, transcription factors, and seed dormancy pathways regulate the endosperm development process. Another 199 functional DEGs participate in glycolysis, pentose phosphate pathway, citrate cycle, FA biosynthesis, TAG assembly, and other pathways. A key transcription factor (WRI1) and some important target genes (ACCase, FATA, LPCAT, FADs, and DGAT etc.) were found in the comprehensive genetic networks of oil biosynthesis.

Insights into the correlation between Physiological changes in and seed development of tartary buckwheat (Fagopyrum tataricum Gaertn.)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum Gaertn.) is a widely cultivated medicinal and edible crop with excellent economic and nutritional value. The development of tartary buckwheat seeds is a very complex process involving many expression-dependent physiological changes and regulation of a large number of genes and phytohormones. In recent years, the gene regulatory network governing the physiological changes occurring during seed development have received little attention.

RESULTS

Here, we characterized the seed development of tartary buckwheat using light and electron microscopy and measured phytohormone and nutrient accumulation by using high performance liquid chromatography (HPLC) and by profiling the expression of key genes using RNA sequencing with the support of the tartary buckwheat genome. We first divided the development of tartary buckwheat seed into five stages that include complex changes in development, morphology, physiology and phytohormone levels. At the same time, the contents of phytohormones (gibberellin, indole-3-acetic acid, abscisic acid, and zeatin) and nutrients (rutin, starch, total proteins and soluble sugars) at five stages were determined, and their accumulation patterns in the development of tartary buckwheat seeds were analyzed. Second, gene expression patterns of tartary buckwheat samples were compared during three seed developmental stages (13, 19, and 25 days postanthesis, DPA), and 9 765 differentially expressed genes (DEGs) were identified. We analyzed the overlapping DEGs in different sample combinations and measured 665 DEGs in the three samples. Furthermore, expression patterns of DEGs related to phytohormones, flavonoids, starch, and storage proteins were analyzed. Third, we noted the correlation between the trait (physiological changes, nutrient changes) and metabolites during seed development, and discussed the key genes that might be involved in the synthesis and degradation of each of them.

CONCLUSION

We provided abundant genomic resources for tartary buckwheat and Polygonaceae communities and revealed novel molecular insights into the correlations between the physiological changes and seed development of tartary buckwheat.

Evolutionary, structural and expression analysis of core genes involved in starch synthesis

Starch is the main storage carbohydrate in plants and an important natural resource for food, feed and industrial raw materials. However, the details regarding the pathway for starch biosynthesis and the diversity of biosynthetic enzymes involved in this process are poorly understood. This study uses a comprehensive phylogenetic analysis of 74 sequenced plant genomes to revisit the evolutionary history of the genes encoding ADP-glucose pyrophosphorylase (AGPase), starch synthase (SS), starch branching enzyme (SBE) and starch de-branching enzyme (DBE). Additionally, the protein structures and expression patterns of these four core genes in starch biosynthesis were studied to determine their functional differences. The results showed that AGPase, SS, SBE and DBE have undergone complicated evolutionary processes in plants and that gene/genome duplications are responsible for the observed differences in isoform numbers. A structure analysis of these proteins suggested that the deletion/mutation of amino acids in some active sites resulted in not only structural variation but also sub-functionalization or neo-functionalization. Expression profiling indicated that AGPase-, SS-, SBE- and DBE-encoding genes exhibit spatio-temporally divergent expression patterns related to the composition of functional complexes in starch biosynthesis. This study provides a comprehensive atlas of the starch biosynthetic pathway, and these data should support future studies aimed at increasing understanding of starch biosynthesis and the functional evolutionary divergence of AGPase, SS, SBE, and DBE in plants.

GCTTCA as a novel motif for regulating mesocarp-specific expression of the oil palm (Elaeis guineensis Jacq.) stearoyl-ACP desaturase gene

TAAAAT and a novel motif, GCTTCA found in the oil palm stearoyl-ACP desaturase (SAD1) promoter are involved in regulating mesocarp-specific expression. Two key fatty acid biosynthetic genes, stearoyl-ACP desaturase (SAD1), and acyl-carrier protein (ACP3) in Elaeis guineensis (oil palm) showed high level of expression during the period of oil synthesis in the mesocarp [12-19 weeks after anthesis (w.a.a.)] and kernel (12-15 w.a.a.). Both genes are expressed in spear leaves at much lower levels and the expression increased by 1.5-fold to 2.5-fold following treatments with ethylene and abscisic acid (ABA). Both SAD1 and ACP3 promoters contain phytohormone-responsive, light-responsive, abiotic factors/wounding-responsive, endosperm specificity and fruit maturation/ripening regulatory motifs. The activities of the full length and six 5' deletion fragments of the SAD1 promoter were analyzed in transiently transformed oil palm tissues by quantitative β-glucuronidase (GUS) fluorometric assay. The highest SAD1 promoter activity was observed in the mesocarp followed by kernel and the least in the leaves. GUS activity in the D3 deletion construct (- 486 to + 108) was the highest, while the D2 (- 535 to + 108) gave the lowest suggesting the presence of negative cis-acting regulatory element(s) in the deleted - 535 to - 486 (49 bp). It was found that the 49-bp region binds to the nuclear protein extract from mesocarp but not from leaves in electrophoretic mobility shift assay (EMSA). Further fine-tuned analysis of this 49-bp region using truncated DNA led to the identification of GCTTCA as a novel motif in the SAD1 promoter. Interestingly, another known fruit ripening-related motif, LECPLEACS2 (TAAAAT) was found to be required for effective binding of the novel motif to the mesocarp nuclear protein extract.

Sucrose metabolism in developing oil-rich tubers of Cyperus esculentus: comparative transcriptome analysis

BACKGROUND

Cyperus esculentus is unique in that it can accumulate significant amounts of oil, starch and sugar as major storage reserves in tubers with high tuber yield and therefore considered as a novel model to study carbon allocation into different storage reserves in underground sink tissues such as tubers and roots. Sucrose (Suc) plays a central role in control of carbon flux toward biosynthesis of different storage reserves; however, it remains unclear for the molecular mechanism underlying Suc metabolism in underground oil-rich storage tissues. In the present study, a comprehensive transcriptome analysis of C. esculentus oil tuber compared to other plant oil- or carbohydrate-rich storage tissues was made for the expression patterns of genes related to the Suc metabolism.

RESULTS

The results revealed some species-specific features of gene transcripts in oil tuber of C. esculentus, indicating that: (i) the expressions of genes responsible for Suc metabolism are developmentally regulated and displayed a pattern dissimilar to other plant storage tissues; (ii) both of Suc breakdown and biosynthesis processes might be the major pathways associated with Suc metabolism; (iii) it was probably that Suc degradation could be primarily through the action of Suc synthase (SUS) other than invertase (INV) during tuber development. The orthologs of SUS1, SUS3 and SUS4 are the main SUS isoforms catalyzing Suc breakdown while the vacuolar INV (VIN) is the leading determinant controlling sugar composition; (iv) cytosolic hexose phosphorylation possibly relies more on fructose as substrate and uridine diphosphate glucose pyrophosphorylase (UGP) plays an important role in this pathway; (v) it is Suc-phosphate synthase (SPS) B- and C-family members rather than SPS A that are the principal contributors to SPS enzymes and play crucial roles in Suc biosynthesis pathway.

CONCLUSIONS

We have successfully identified the Suc metabolic pathways in C. esculentus tubers, highlighting several conserved and distinct expressions that might contribute to sugar accumulation in this unique underground storage tissue. The specific and differential expression genes revealed in this study might indicate the special molecular mechanism and transcriptional regulation of Suc metabolism occurred in oil tubers of C. esculentus.

The HTPmod Shiny application enables modeling and visualization of large-scale biological data

The wave of high-throughput technologies in genomics and phenomics are enabling data to be generated on an unprecedented scale and at a reasonable cost. Exploring the large-scale data sets generated by these technologies to derive biological insights requires efficient bioinformatic tools. Here we introduce an interactive, open-source web application (HTPmod) for high-throughput biological data modeling and visualization. HTPmod is implemented with the Shiny framework by integrating the computational power and professional visualization of R and including various machine-learning approaches. We demonstrate that HTPmod can be used for modeling and visualizing large-scale, high-dimensional data sets (such as multiple omics data) under a broad context. By reinvestigating example data sets from recent studies, we find not only that HTPmod can reproduce results from the original studies in a straightforward fashion and within a reasonable time, but also that novel insights may be gained from fast reinvestigation of existing data by HTPmod.

Dijun Chen et al. present HTPmod, a Shiny web application for modeling and visualization of large-scale genomic and phenomic datasets. The authors show that HTPmod can quickly reproduce analyses of high-throughput biological datasets and produce publication-quality figures.

Distinct tissue-specific transcriptional regulation revealed by gene regulatory networks in maize

BACKGROUND

Transcription factors (TFs) are proteins that can bind to DNA sequences and regulate gene expression. Many TFs are master regulators in cells that contribute to tissue-specific and cell-type-specific gene expression patterns in eukaryotes. Maize has been a model organism for over one hundred years, but little is known about its tissue-specific gene regulation through TFs. In this study, we used a network approach to elucidate gene regulatory networks (GRNs) in four tissues (leaf, root, SAM and seed) in maize. We utilized GENIE3, a machine-learning algorithm combined with large quantity of RNA-Seq expression data to construct four tissue-specific GRNs. Unlike some other techniques, this approach is not limited by high-quality Position Weighed Matrix (PWM), and can therefore predict GRNs for over 2000 TFs in maize.

RESULTS

Although many TFs were expressed across multiple tissues, a multi-tiered analysis predicted tissue-specific regulatory functions for many transcription factors. Some well-studied TFs emerged within the four tissue-specific GRNs, and the GRN predictions matched expectations based upon published results for many of these examples. Our GRNs were also validated by ChIP-Seq datasets (KN1, FEA4 and O2). Key TFs were identified for each tissue and matched expectations for key regulators in each tissue, including GO enrichment and identity with known regulatory factors for that tissue. We also found functional modules in each network by clustering analysis with the MCL algorithm.

CONCLUSIONS

By combining publicly available genome-wide expression data and network analysis, we can uncover GRNs at tissue-level resolution in maize. Since ChIP-Seq and PWMs are still limited in several model organisms, our study provides a uniform platform that can be adapted to any species with genome-wide expression data to construct GRNs. We also present a publicly available database, maize tissue-specific GRN (mGRN, https://www.bio.fsu.edu/mcginnislab/mgrn/ ), for easy querying. All source code and data are available at Github ( https://github.com/timedreamer/maize_tissue-specific_GRN ).

Proliferation of Regulatory DNA Elements Derived from Transposable Elements in the Maize Genome

Genomic regions free of nucleosomes, which are hypersensitive to DNase I digestion, are known as DNase I hypersensitive sites (DHSs) and frequently contain cis-regulatory DNA elements. To investigate their prevalence and characteristics in maize (Zea mays), we developed high-resolution genome-wide DHS maps using a modified DNase-seq technique. Maize DHSs exhibit depletion of nucleosomes and low levels of DNA methylation and are enriched with conserved noncoding sequences (CNSs). We developed a protoplast-based transient transformation assay to assess the potential gene expression enhancer and/or promoter functions associated with DHSs, which showed that more than 80% of DHSs overlapping with CNSs showed an enhancer function. Strikingly, nearly 25% of maize DHSs were derived from transposable elements (TEs), including both class I and class II transposons. Interestingly, TE-derived DHSs (teDHSs) homologous to retrotransposons were enriched with sequences related to the intrinsic cis-regulatory elements within the long terminal repeats of retrotransposons. We demonstrate that more than 80% of teDHSs can drive transcription of a reporter gene in protoplast assays. These results reveal the widespread occurrence of TE-derived cis-regulatory sequences and suggest that teDHSs play a major role in transcriptional regulation in maize.

Evolution of intron-poor clades and expression patterns of the glycosyltransferase family 47

A large-scale bioinformatics analysis revealed the origin and evolution of GT47 gene family, and identified two clades of intron-poor genes with putative functions in drought stress responses and seed development in maize. Glycosyltransferase family 47 (GT47) genes encode β-galactosyltransferases and β-glucuronyltransferases that synthesize pectin, xyloglucans and xylan, which are important components of the plant cell wall. In this study, we performed a systematic and large-scale bioinformatics analysis of GT47 gene family using 352 GT47 proteins from 15 species ranging from cyanobacteria to seed plants. The analysis results showed that GT47 family may originate in cyanobacteria and expand along the evolutionary trajectory to moss. Further analysis of 47 GT47 genes in maize revealed that they can divide into five clades with diverse exon-intron structures. Among these five clades, two were mainly composed with intron-poor genes, which may originate in the moss. Gene duplication analysis revealed that the expansion of GT47 gene family in maize was significantly driven from tandem duplication events and segmental duplication events. Significantly, almost all duplicated genes are intron-poor genes. Expression analysis indicated that several intron-poor GT47 genes may be involved in the drought stress response and seed development in maize. This work provides insight into the origin and evolutionary process, expansion mechanisms and expression patterns of GT47 genes, thus facilitating their functional investigations in the future.

Sequential gene activation and gene imprinting during early embryo development in maize

Gene imprinting is a widely observed epigenetic phenomenon in maize endosperm; however, whether it also occurs in the maize embryo remains controversial. Here, we used high-throughput RNA sequencing on laser capture microdissected and manually dissected maize embryos from reciprocal crosses between inbred lines B73 and Mo17 at six time points (3-13 days after pollination, DAP) to analyze allelic gene expression patterns. Co-expression analysis revealed sequential gene activation during maize embryo development. Gene imprinting was observed in maize embryos, and a greater number of imprinted genes were identified at early embryo stages. Sixty-four strongly imprinted genes were identified (at the threshold of 9:1) on manually dissected embryos 5-13 DAP (more imprinted genes at 5 DAP). Forty-one strongly imprinted genes were identified from laser capture microdissected embryos at 3 and 5 DAP (more imprinted genes at 3 DAP). Furthermore, of the 56 genes that were completely imprinted (at the threshold of 99:1), 36 were not previously identified as imprinted genes in endosperm or embryos. In situ hybridization demonstrated that most of the imprinted genes were expressed abundantly in maize embryonic tissue. Our results shed lights on early maize embryo development and provide evidence to support that gene imprinting occurs in maize embryos.

The genetic architecture of amylose biosynthesis in maize kernel

Starch is the most abundant storage carbohydrate in maize kernel. The content of amylose and amylopectin confers unique properties in food processing and industrial application. Thus, the resurgent interest has been switched to the study of individual amylose or amylopectin rather than total starch, whereas the enzymatic machinery for amylose synthesis remains elusive. We took advantage of the phenotype of amylose content and the genotype of 9,007,194 single nucleotide polymorphisms from 464 inbred maize lines. The genome-wide association study identified 27 associated loci involving 39 candidate genes that were linked to amylose content including transcription factors, glycosyltransferases, glycosidases, as well as hydrolases. Except the waxy gene that encodes the granule-bound starch synthase, the remaining candidate genes were located in the upstream pathway of amylose synthesis, while the downstream members were already known from prior studies. The linked candidate genes could be transferred to manipulate amylose content and thus add value to maize kernel in the breeding programme.

Defective Kernel 39 encodes a PPR protein required for seed development in maize

RNA editing is a posttranscriptional process that is important in mitochondria and plastids of higher plants. All RNA editing-specific trans-factors reported so far belong to PLS-class of pentatricopeptide repeat (PPR) proteins. Here, we report the map-based cloning and molecular characterization of a defective kernel mutant dek39 in maize. Loss of Dek39 function leads to delayed embryogenesis and endosperm development, reduced kernel size, and seedling lethality. Dek39 encodes an E sub-class PPR protein that targets to both mitochondria and chloroplasts, and is involved in RNA editing in mitochondrial NADH dehydrogenase3 (nad3) at nad3-247 and nad3-275. C-to-U editing of nad3-275 is not conserved and even lost in Arabidopsis, consistent with the idea that no close DEK39 homologs are present in Arabidopsis. However, the amino acids generated by editing nad3-247 and nad3-275 are highly conserved in many other plant species, and the reductions of editing at these two sites decrease the activity of mitochondria NADH dehydrogenase complex I, indicating that the alteration of amino acid sequence is necessary for Nad3 function. Our results indicate that Dek39 encodes an E sub-class PPR protein that is involved in RNA editing of multiple sites and is necessary for seed development of maize.

Functions of maize genes encoding pyruvate phosphate dikinase in developing endosperm

Maize opaque2 (o2) mutations are beneficial for endosperm nutritional quality but cause negative pleiotropic effects for reasons that are not fully understood. Direct targets of the bZIP transcriptional regulator encoded by o2 include pdk1 and pdk2 that specify pyruvate phosphate dikinase (PPDK). This enzyme reversibly converts AMP, pyrophosphate, and phosphoenolpyruvate to ATP, orthophosphate, and pyruvate and provides diverse functions in plants. This study addressed PPDK function in maize starchy endosperm where it is highly abundant during grain fill. pdk1 and pdk2 were inactivated individually by transposon insertions, and both genes were simultaneously targeted by endosperm-specific RNAi. pdk2 accounts for the large majority of endosperm PPDK, whereas pdk1 specifies the abundant mesophyll form. The pdk1- mutation is seedling-lethal, indicating that C4 photosynthesis is essential in maize. RNAi expression in transgenic endosperm eliminated detectable PPDK protein and enzyme activity. Transgenic kernels weighed the same on average as nontransgenic siblings, with normal endosperm starch and total N contents, indicating that PPDK is not required for net storage compound synthesis. An opaque phenotype resulted from complete PPDK knockout, including loss of vitreous endosperm character similar to the phenotype conditioned by o2-. Concentrations of multiple glycolytic intermediates were elevated in transgenic endosperm, energy charge was altered, and starch granules were more numerous but smaller on average than normal. The data indicate that PPDK modulates endosperm metabolism, potentially through reversible adjustments to energy charge, and reveal that o2- mutations can affect the opaque phenotype through regulation of PPDK in addition to their previously demonstrated effects on storage protein gene expression.

High-Throughput Sequencing of Small RNA Transcriptomes in Maize Kernel Identifies miRNAs Involved in Embryo and Endosperm Development

Maize kernel development is a complex biological process that involves the temporal and spatial expression of many genes and fine gene regulation at a transcriptional and post-transcriptional level, and microRNAs (miRNAs) play vital roles during this process. To gain insight into miRNA-mediated regulation of maize kernel development, a deep-sequencing technique was used to investigate the dynamic expression of miRNAs in the embryo and endosperm at three developmental stages in B73. By miRNA transcriptomic analysis, we characterized 132 known miRNAs and six novel miRNAs in developing maize kernel, among which, 15 and 14 miRNAs were commonly differentially expressed between the embryo and endosperm at 9 days after pollination (DAP), 15 DAP and 20 DAP respectively. Conserved miRNA families such as miR159, miR160, miR166, miR390, miR319, miR528 and miR529 were highly expressed in developing embryos; miR164, miR171, miR393 and miR2118 were highly expressed in developing endosperm. Genes targeted by those highly expressed miRNAs were found to be largely related to a regulation category, including the transcription, macromolecule biosynthetic and metabolic process in the embryo as well as the vitamin biosynthetic and metabolic process in the endosperm. Quantitative reverse transcription-PCR (qRT-PCR) analysis showed that these miRNAs displayed a negative correlation with the levels of their corresponding target genes. Importantly, our findings revealed that members of the miR169 family were highly and dynamically expressed in the developing kernel, which will help to exploit new players functioning in maize kernel development.

Highly interwoven communities of a gene regulatory network unveil topologically important genes for maize seed development

The complex interactions between transcription factors (TFs) and their target genes in a spatially and temporally specific manner are crucial to all cellular processes. Reconstruction of gene regulatory networks (GRNs) from gene expression profiles can help to decipher TF-gene regulations in a variety of contexts; however, the inevitable prediction errors of GRNs hinder optimal data mining of RNA-Seq transcriptome profiles. Here we perform an integrative study of Zea mays (maize) seed development in order to identify key genes in a complex developmental process. First, we reverse engineered a GRN from 78 maize seed transcriptome profiles. Then, we studied collective gene interaction patterns and uncovered highly interwoven network communities as the building blocks of the GRN. One community, composed of mostly unknown genes interacting with opaque2, brittle endosperm1 and shrunken2, contributes to seed phenotypes. Another community, composed mostly of genes expressed in the basal endosperm transfer layer, is responsible for nutrient transport. We further integrated our inferred GRN with gene expression patterns in different seed compartments and at various developmental stages and pathways. The integration facilitated a biological interpretation of the GRN. Our yeast one-hybrid assays verified six out of eight TF-promoter bindings in the reconstructed GRN. This study identified topologically important genes in interwoven network communities that may be crucial to maize seed development.

Transcriptome dynamics revealed by a gene expression atlas of the early Arabidopsis embryo

During early plant embryogenesis, precursors for all major tissues and stem cells are formed. While several components of the regulatory framework are known, how cell fates are instructed by genome-wide transcriptional activity remains unanswered-in part because of difficulties in capturing transcriptome changes at cellular resolution. Here, we have adapted a two-component transgenic labelling system to purify cell-type-specific nuclear RNA and generate a transcriptome atlas of early Arabidopsis embryo development, with a focus on root stem cell niche formation. We validated the dataset through gene expression analysis, and show that gene activity shifts in a spatio-temporal manner, probably signifying transcriptional reprogramming, to induce developmental processes reflecting cell states and state transitions. This atlas provides the most comprehensive tissue- and cell-specific description of genome-wide gene activity in the early plant embryo, and serves as a valuable resource for understanding the genetic control of early plant development.

Transcriptome dynamics revealed by a gene expression atlas of the early Arabidopsis embryo

During early plant embryogenesis, precursors for all major tissues and stem cells are formed. While several components of the regulatory framework are known, how cell fates are instructed by genome-wide transcriptional activity remains unanswered-in part because of difficulties in capturing transcriptome changes at cellular resolution. Here, we have adapted a two-component transgenic labelling system to purify cell-type-specific nuclear RNA and generate a transcriptome atlas of early Arabidopsis embryo development, with a focus on root stem cell niche formation. We validated the dataset through gene expression analysis, and show that gene activity shifts in a spatio-temporal manner, probably signifying transcriptional reprogramming, to induce developmental processes reflecting cell states and state transitions. This atlas provides the most comprehensive tissue- and cell-specific description of genome-wide gene activity in the early plant embryo, and serves as a valuable resource for understanding the genetic control of early plant development.

Transcriptomic Analysis of Long Non-Coding RNAs and Coding Genes Uncovers a Complex Regulatory Network That Is Involved in Maize Seed Development

Long non-coding RNAs (lncRNAs) have been reported to be involved in the development of maize plant. However, few focused on seed development of maize. Here, we identified 753 lncRNA candidates in maize genome from six seed samples. Similar to the mRNAs, lncRNAs showed tissue developmental stage specific and differential expression, indicating their putative role in seed development. Increasing evidence shows that crosstalk among RNAs mediated by shared microRNAs (miRNAs) represents a novel layer of gene regulation, which plays important roles in plant development. Functional roles and regulatory mechanisms of lncRNAs as competing endogenous RNAs (ceRNA) in plants, particularly in maize seed development, are unclear. We combined analyses of consistently altered 17 lncRNAs, 840 mRNAs and known miRNA to genome-wide investigate potential lncRNA-mediated ceRNA based on “ceRNA hypothesis”. The results uncovered seven novel lncRNAs as potential functional ceRNAs. Functional analyses based on their competitive coding-gene partners by Gene Ontology (GO) and KEGG biological pathway demonstrated that combined effects of multiple ceRNAs can have major impacts on general developmental and metabolic processes in maize seed. These findings provided a useful platform for uncovering novel mechanisms of maize seed development and may provide opportunities for the functional characterization of individual lncRNA in future studies.

The Maize Imprinted Gene Floury3 Encodes a PLATZ Protein Required for tRNA and 5S rRNA Transcription through Interaction with RNA Polymerase III

Maize (Zea mays) floury3 (fl3) is a classic semidominant negative mutant that exhibits severe defects in the endosperm but fl3 plants otherwise appear normal. We cloned the fl3 gene and determined that it encodes a PLATZ (plant AT-rich sequence and zinc binding) protein. The mutation in fl3 resulted in an Asn-to-His replacement in the conserved PLATZ domain, creating a dominant allele. Fl3 is specifically expressed in starchy endosperm cells and regulated by genomic imprinting, which leads to the suppressed expression of fl3 when transmitted through the male, perhaps as a consequence the semidominant behavior. Yeast two-hybrid screening and bimolecular luciferase complementation experiments revealed that FL3 interacts with the RNA polymerase III subunit 53 (RPC53) and transcription factor class C 1 (TFC1), two critical factors of the RNA polymerase III (RNAPIII) transcription complex. In the fl3 endosperm, the levels of many tRNAs and 5S rRNA that are transcribed by RNAPIII are significantly reduced, suggesting that the incorrectly folded fl3 protein may impair the function of RNAPIII. The transcriptome is dramatically altered in fl3 mutants, in which the downregulated genes are primarily enriched in pathways related to translation, ribosome, misfolded protein responses, and nutrient reservoir activity. Collectively, these changes may lead to defects in endosperm development and storage reserve filling in fl3 seeds.

Seed maturation associated transcriptional programs and regulatory networks underlying genotypic difference in seed dormancy and size/weight in wheat (Triticum aestivum L.)

BACKGROUND

Maturation forms one of the critical seed developmental phases and it is characterized mainly by programmed cell death, dormancy and desiccation, however, the transcriptional programs and regulatory networks underlying acquisition of dormancy and deposition of storage reserves during the maturation phase of seed development are poorly understood in wheat. The present study performed comparative spatiotemporal transcriptomic analysis of seed maturation in two wheat genotypes with contrasting seed weight/size and dormancy phenotype.

RESULTS

The embryo and endosperm tissues of maturing seeds appeared to exhibit genotype-specific temporal shifts in gene expression profile that might contribute to the seed phenotypic variations. Functional annotations of gene clusters suggest that the two tissues exhibit distinct but genotypically overlapping molecular functions. Motif enrichment predicts genotypically distinct abscisic acid (ABA) and gibberellin (GA) regulated transcriptional networks contribute to the contrasting seed weight/size and dormancy phenotypes between the two genotypes. While other ABA responsive element (ABRE) motifs are enriched in both genotypes, the prevalence of G-box-like motif specifically in tissues of the dormant genotype suggests distinct ABA mediated transcriptional mechanisms control the establishment of dormancy during seed maturation. In agreement with this, the bZIP transcription factors that co-express with ABRE enriched embryonic genes differ with genotype. The enrichment of SITEIIATCYTC motif specifically in embryo clusters of maturing seeds irrespective of genotype predicts a tissue specific role for the respective TCP transcription factors with no or minimal contribution to the variations in seed dormancy.

CONCLUSION

The results of this study advance our understanding of the seed maturation associated molecular mechanisms underlying variation in dormancy and weight/size in wheat seeds, which is a critical step towards the designing of molecular strategies for enhancing seed yield and quality.

iTRAQ-Based Proteomics Analysis and Network Integration for Kernel Tissue Development in Maize

Grain weight is one of the most important yield components and a developmentally complex structure comprised of two major compartments (endosperm and pericarp) in maize (Zea mays L.), however, very little is known concerning the coordinated accumulation of the numerous proteins involved. Herein, we used isobaric tags for relative and absolute quantitation (iTRAQ)-based comparative proteomic method to analyze the characteristics of dynamic proteomics for endosperm and pericarp during grain development. Totally, 9539 proteins were identified for both components at four development stages, among which 1401 proteins were non-redundant, 232 proteins were specific in pericarp and 153 proteins were specific in endosperm. A functional annotation of the identified proteins revealed the importance of metabolic and cellular processes, and binding and catalytic activities for the tissue development. Three and 76 proteins involved in 49 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were integrated for the specific endosperm and pericarp proteins, respectively, reflecting their complex metabolic interactions. In addition, four proteins with important functions and different expression levels were chosen for gene cloning and expression analysis. Different concordance between mRNA level and the protein abundance was observed across different proteins, stages, and tissues as in previous research. These results could provide useful message for understanding the developmental mechanisms in grain development in maize.

Maize network analysis revealed gene modules involved in development, nutrients utilization, metabolism, and stress response

BACKGROUND

The advent of big data in biology offers opportunities while poses challenges to derive biological insights. For maize, a large amount of publicly available transcriptome datasets have been generated but a comprehensive analysis is lacking.

RESULTS

We constructed a maize gene co-expression network based on the graphical Gaussian model, using massive RNA-seq data. The network, containing 20,269 genes, assembles into 964 gene modules that function in a variety of plant processes, such as cell organization, the development of inflorescences, ligules and kernels, the uptake and utilization of nutrients (e.g. nitrogen and phosphate), the metabolism of benzoxazionids, oxylipins, flavonoids, and wax, and the response to stresses. Among them, the inflorescences development module is enriched with domestication genes (like ra1, ba1, gt1, tb1, tga1) that control plant architecture and kernel structure, while multiple other modules relate to diverse agronomic traits. Contained within these modules are transcription factors acting as known or potential expression regulators for the genes within the same modules, suggesting them as candidate regulators for related biological processes. A comparison with an established Arabidopsis network revealed conserved gene association patterns for specific modules involved in cell organization, nutrients uptake & utilization, and metabolism. The analysis also identified significant divergences between the two species for modules that orchestrate developmental pathways.

CONCLUSIONS

This network sheds light on how gene modules are organized between different species in the context of evolutionary divergence and highlights modules whose structure and gene content can provide important resources for maize gene functional studies with application potential.

Genome-wide Nucleosome Occupancy and Organization Modulates the Plasticity of Gene Transcriptional Status in Maize

Nucleosomes are fundamental units of chromatin that play critical roles in gene regulation by modulating DNA accessibility. However, their roles in regulating tissue-specific gene transcription are poorly understood. Here, we present genome-wide nucleosome maps of maize shoot and endosperm generated by sequencing the micrococcal nuclease digested nucleosomal DNA. The changes of gene transcriptional status between shoot and endosperm were accompanied by preferential nucleosome loss from the promoters and shifts in the first nucleosome downstream of the transcriptional start sites (+1 nucleosome) and upstream of transcriptional termination sites (-1 nucleosome). Intrinsically DNA-encoded nucleosome organization was largely associated with the capacity of a gene to alter its transcriptional status among different tissues. Compared with tissue-specific genes, constitutively expressed genes showed more pronounced 5' and 3' nucleosome-depleted regions as well as further +1 nucleosome to transcriptional start sites and -1 nucleosome to transcriptional termination sites. Moreover, nucleosome organization was more highly correlated with the plasticity of gene transcriptional status than with its expression level when examined using in vivo and predicted nucleosome data. In addition, the translational efficiencies of tissue-specific genes appeared to be greater than those of constitutively expressed genes. Taken together, our results indicate that intrinsically DNA-encoded nucleosome organization is important, beyond its role in regulating gene expression levels, in determining the plasticity of gene transcriptional status.

Spatio-temporal analysis of coding and long noncoding transcripts during maize endosperm development

The maize endosperm consists of three major compartmentalized cell types: the starchy endosperm (SE), the basal endosperm transfer cell layer (BETL), and the aleurone cell layer (AL). Differential genetic programs are activated in each cell type to construct functionally and structurally distinct cells. To compare gene expression patterns involved in maize endosperm cell differentiation, we isolated transcripts from cryo-dissected endosperm specimens enriched with BETL, AL, or SE at 8, 12, and 16 days after pollination (DAP). We performed transcriptome profiling of coding and long noncoding transcripts in the three cell types during differentiation and identified clusters of the transcripts exhibiting spatio-temporal specificities. Our analysis uncovered that the BETL at 12 DAP undergoes the most dynamic transcriptional regulation for both coding and long noncoding transcripts. In addition, our transcriptome analysis revealed spatio-temporal regulatory networks of transcription factors, imprinted genes, and loci marked with histone H3 trimethylated at lysine 27. Our study suggests that various regulatory mechanisms contribute to the genetic networks specific to the functions and structures of the cell types of the endosperm.

Systems biology of seeds: deciphering the molecular mechanisms of seed storage, dormancy and onset of germination

Seeds are heterogeneous storage reserves with wide array of storage compounds that include various soluble carbohydrates, starch polymer, storage proteins and lipids. These stored reserves comprise 70% of the world's caloric intake in the form of food and animal feed produced through sustainable agriculture, which contributes to food and nutritional security. Seed systems biology remains an enigmatic subject in understanding seed storage processes, maturation and pre-germinative metabolism. The reviews and research articles covered in this special issue of Plant Cell Reports highlight recent advances made in the area of seed biology that cover various systems biology applications such as gene regulatory networks, metabolomics, epigenetics and the role of micro-RNA in seed development.

Dek35 Encodes a PPR Protein that Affects cis-Splicing of Mitochondrial nad4 Intron 1 and Seed Development in Maize

In higher plants, the splicing of organelle-encoded mRNA involves a complex collaboration with nuclear-encoded proteins. Pentatricopeptide repeat (PPR) proteins have been implicated in these RNA-protein interactions. In this study, we performed the cloning and functional characterization of maize Defective kernel 35 (Dek35). The dek35-ref mutant is a lethal-seed mutant with developmental deficiency. Dek35 was cloned through Mutator tag isolation and further confirmed by four additional independent mutant alleles. Dek35 encodes an P-type PPR protein that targets the mitochondria. The dek35 mutation causes significant reduction in the accumulation of DEK35 proteins and reduced splicing efficiency of mitochondrial nad4 intron 1. Analysis of mitochondrial complex in dek35 immature seeds indicated severe deficiency in the complex I assembly and NADH dehydrogenase activity. Transcriptome analysis of dek35 endosperm revealed enhanced expression of genes involved in the alternative respiratory pathway and extensive differentially expressed genes related to mitochondrial function and activity. Collectively, these results indicate that Dek35 encodes an PPR protein that affects the cis-splicing of mitochondrial nad4 intron 1 and is required for mitochondrial function and seed development.

Distant eQTLs and Non-coding Sequences Play Critical Roles in Regulating Gene Expression and Quantitative Trait Variation in Maize

A detailed understanding of genetic architecture of mRNA expression by millions of genetic variants is important for studying quantitative trait variation. In this study, we identified 1.25M SNPs with a minor allele frequency greater than 0.05 by combining reduced genome sequencing (GBS), high-density array technologies (600K), and previous deep RNA-sequencing data from 368 diverse inbred lines of maize. The balanced allelic frequencies and distributions in a relatively large and diverse natural panel helped to identify expression quantitative trait loci (eQTLs) associated with more than 18 000 genes (63.4% of tested genes). We found that distant eQTLs were more frequent (∼75% of all eQTLs) across the whole genome. Thirteen novel associated loci affecting maize kernel oil concentration were identified using the new dataset, among which one intergenic locus affected the kernel oil variation by controlling expression of three other known oil-related genes. Altogether, this study provides resources for expanding our understanding of cellular regulatory mechanisms of transcriptome variation and the landscape of functional variants within the maize genome, thereby enhancing the understanding of quantitative variations.

Aberrant splicing in maize rough endosperm3 reveals a conserved role for U12 splicing in eukaryotic multicellular development

RNA splicing of U12-type introns functions in human cell differentiation, but it is not known whether this class of introns has a similar role in plants. The maize ROUGH ENDOSPERM3 (RGH3) protein is orthologous to the human splicing factor, ZRSR2. ZRSR2 mutations are associated with myelodysplastic syndrome (MDS) and cause U12 splicing defects. Maize rgh3 mutants have aberrant endosperm cell differentiation and proliferation. We found that most U12-type introns are retained or misspliced in rgh3 Genes affected in rgh3 and ZRSR2 mutants identify cell cycle and protein glycosylation as common pathways disrupted. Transcripts with retained U12-type introns can be found in polysomes, suggesting that splicing efficiency can alter protein isoforms. The rgh3 mutant protein disrupts colocalization with a known ZRSR2-interacting protein, U2AF2. These results indicate conserved function for RGH3/ZRSR2 in U12 splicing and a deeply conserved role for the minor spliceosome to promote cell differentiation from stem cells to terminal fates.

zmsbt1 and zmsbt2, two new subtilisin-like serine proteases genes expressed in early maize kernel development

Two subtilisin-like proteases show highly specific and complementary expression patterns in developing grains. These genes label the complete surface of the filial-maternal interface, suggesting a role in filial epithelial differentiation. The cereal endosperm is the most important source of nutrition and raw materials for mankind, as well as the storage compartment enabling initial growth of the germinating plantlets. The development of the different cell types in this tissue is regulated environmentally, genetically and epigenetically, resulting in the formation of top-bottom, adaxial-abaxial and surface-central axes. However, the mechanisms governing the interactions among the different inputs are mostly unknown. We have screened a kernel cDNA library for tissue-specific transcripts as initial step to identify genes relevant in cell differentiation. We report here on the isolation of two maize subtilisin-related genes that show grain-specific, surficial expression. zmsbt1 (Zea mays Subtilisin1) is expressed at the developing aleurone in a time-regulated manner, while zmsbt2 concentrates at the pedicel in front of the endosperm basal transfer layer. We have shown that their presence, early in the maize caryopsis development, is dependent on proper initial tissue determination, and have isolated their promoters to produce transgenic reporter lines that assist in the study of their regulation.

A Comprehensive Analysis of Alternative Splicing in Paleopolyploid Maize

Identifying and characterizing alternative splicing (AS) enables our understanding of the biological role of transcript isoform diversity. This study describes the use of publicly available RNA-Seq data to identify and characterize the global diversity of AS isoforms in maize using the inbred lines B73 and Mo17, and a related species, sorghum. Identification and characterization of AS within maize tissues revealed that genes expressed in seed exhibit the largest differential AS relative to other tissues examined. Additionally, differences in AS between the two genotypes B73 and Mo17 are greatest within genes expressed in seed. We demonstrate that changes in the level of alternatively spliced transcripts (intron retention and exon skipping) do not solely reflect differences in total transcript abundance, and we present evidence that intron retention may act to fine-tune gene expression across seed development stages. Furthermore, we have identified temperature sensitive AS in maize and demonstrate that drought-induced changes in AS involve distinct sets of genes in reproductive and vegetative tissues. Examining our identified AS isoforms within B73 × Mo17 recombinant inbred lines (RILs) identified splicing QTL (sQTL). The 43.3% of cis-sQTL regulated junctions are actually identified as alternatively spliced junctions in our analysis, while 10 Mb windows on each side of 48.2% of trans-sQTLs overlap with splicing related genes. Using sorghum as an out-group enabled direct examination of loss or conservation of AS between homeologous genes representing the two subgenomes of maize. We identify several instances where AS isoforms that are conserved between one maize homeolog and its sorghum ortholog are absent from the second maize homeolog, suggesting that these AS isoforms may have been lost after the maize whole genome duplication event. This comprehensive analysis provides new insights into the complexity of AS in maize.

Natural variations in stearoyl-acp desaturase genes affect the conversion of stearic to oleic acid in maize kernel

We identified 11 SAD genes, and mined their natural variations associated with the conservation of stearic to oleic acid, especially ZmSAD1 supported by both the QTL and an expression QTL. Maize oil is generally regarded as a healthy vegetable oil owing to its low abundance of saturated fatty acids. Stearoyl-ACP desaturase (SAD) is a key rate-limiting enzyme for the conservation of stearic (C18:0) to oleic (C18:1) acid. Here, 11 maize SAD genes were identified to have more divergent functions than Arabidopsis SAD genes. The genomic regional associations in a maize panel including 508 inbred lines identified 6 SAD genes significantly associated (P < 0.01) with the C18:0/C18:1 ratio or the level of C18:0 or C18:1, one gene of which co-localized with a quantitative trait locus (QTL) and 5 of which co-localized with an expression QTL. ZmSAD1, supported by both the QTL and an expression QTL, had the largest effect on C18:0/C18:1. One nonsynonymous single-nucleotide polymorphism in exon 3 and one 5-bp insertion/deletion in the 3' untranslated region were further shown to contribute to the natural variation in C18:0/C18:1 according to ZmSAD1-based association mapping. Finally, selection tests of ZmSAD1 in teosinte, regular maize, and high-oil maize indicated that ZmSAD1 was not a selection target during the process of maize domestication and high-oil maize development. These results will guide the manipulation of the ratio between saturated and unsaturated fatty acids in maize.

Natural variations in stearoyl-acp desaturase genes affect the conversion of stearic to oleic acid in maize kernel

We identified 11 SAD genes, and mined their natural variations associated with the conservation of stearic to oleic acid, especially ZmSAD1 supported by both the QTL and an expression QTL. Maize oil is generally regarded as a healthy vegetable oil owing to its low abundance of saturated fatty acids. Stearoyl-ACP desaturase (SAD) is a key rate-limiting enzyme for the conservation of stearic (C18:0) to oleic (C18:1) acid. Here, 11 maize SAD genes were identified to have more divergent functions than Arabidopsis SAD genes. The genomic regional associations in a maize panel including 508 inbred lines identified 6 SAD genes significantly associated (P < 0.01) with the C18:0/C18:1 ratio or the level of C18:0 or C18:1, one gene of which co-localized with a quantitative trait locus (QTL) and 5 of which co-localized with an expression QTL. ZmSAD1, supported by both the QTL and an expression QTL, had the largest effect on C18:0/C18:1. One nonsynonymous single-nucleotide polymorphism in exon 3 and one 5-bp insertion/deletion in the 3' untranslated region were further shown to contribute to the natural variation in C18:0/C18:1 according to ZmSAD1-based association mapping. Finally, selection tests of ZmSAD1 in teosinte, regular maize, and high-oil maize indicated that ZmSAD1 was not a selection target during the process of maize domestication and high-oil maize development. These results will guide the manipulation of the ratio between saturated and unsaturated fatty acids in maize.

Genome-wide association study (GWAS) reveals the genetic architecture of four husk traits in maize

Background

Maize (Zea mays) husk referring to the leafy outer enclosing the ear, plays an important role in grain production by directly contributing photosynthate and protecting ear from pathogen infection. Although the physiological functions related to husk have been extensively studied, little is known about its morphological variation and genetic basis in natural population.

Results

Here we utilized a maize association panel including 508 inbred lines with tropical, subtropical and temperate backgrounds to decipher the genetic architecture attributed to four husk traits, i.e. number of layers, length, width and thickness. Evaluating the phenotypic diversity at two different environments showed that four traits exhibit broadly natural variations and moderate levels of heritability with 0.64, 0.74, 0.49 and 0.75 for number, length, width and thickness, respectively. Diversity analysis indicated that different traits have dissimilar responses to subpopulation effects. A series of significantly positive or negative correlations between husk phenotypes and other agronomic traits were identified, indicating that husk growth is coordinated with other developmental processes. Combining husk traits with about half of a million of single nucleotide polymorphisms (SNPs) via genome-wide association study revealed a total of 9 variants significantly associated with traits at P < 1.04 × 10^-5, which are implicated in multiple functional categories, such as cellular trafficking, transcriptional regulation and metabolism.

Conclusions

These results provide instrumental information for understanding the genetic basis of husk development, and further studies on identified candidate genes facilitate to illuminate molecular pathways regulating maize husk growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3229-6) contains supplementary material, which is available to authorized users.

Proteomic analysis of maize grain development using iTRAQ reveals temporal programs of diverse metabolic processes

BACKGROUND

Grain development in maize is an essential process in the plant's life cycle and is vital for use of the plant as a crop for animals and humans. However, little is known regarding the protein regulatory networks that control grain development. Here, isobaric tag for relative and absolute quantification (iTRAQ) technology was used to analyze temporal changes in protein expression during maize grain development.

RESULTS

Maize grain proteins and changes in protein expression at eight developmental stages from 3 to 50 d after pollination (DAP) were performed using iTRAQ-based proteomics. Overall, 4751 proteins were identified; 2639 of these were quantified and 1235 showed at least 1.5-fold changes in expression levels at different developmental stages and were identified as differentially expressed proteins (DEPs). The DEPs were involved in different cellular and metabolic processes with a preferential distribution to protein synthesis/destination and metabolism categories. A K-means clustering analysis revealed coordinated protein expression associated with different functional categories/subcategories at different development stages.

CONCLUSIONS

Our results revealed developing maize grain display different proteomic characteristics at distinct stages, such as numerous DEPs for cell growth/division were highly expressed during early stages, whereas those for starch biosynthesis and defense/stress accumulated in middle and late stages, respectively. We also observed coordinated expression of multiple proteins of the antioxidant system, which are essential for the maintenance of reactive oxygen species (ROS) homeostasis during grain development. Particularly, some DEPs, such as zinc metallothionein class II, pyruvate orthophosphate dikinase (PPDK) and 14-3-3 proteins, undergo major changes in expression at specific developmental stages, suggesting their roles in maize grain development. These results provide a valuable resource for analyzing protein function on a global scale and also provide new insights into the potential protein regulatory networks that control grain yield and quality.

Mitochondrial Function and Maize Kernel Development Requires Dek2, a Pentatricopeptide Repeat Protein Involved in nad1 mRNA Splicing

In flowering plants, many respiration-related proteins are encoded by the mitochondrial genome and the splicing of mitochondrion-encoded messenger RNA (mRNA) involves a complex collaboration with nuclear-encoded proteins. Pentatricopeptide repeat (PPR) proteins have been implicated in these RNA-protein interactions. Maize defective kernel 2 (dek2) is a classic mutant with small kernels and delayed development. Through positional cloning and allelic confirmation, we found Dek2 encodes a novel P-type PPR protein that targets mitochondria. Mitochondrial transcript analysis indicated that dek2 mutation causes reduced splicing efficiency of mitochondrial nad1 intron 1. Mitochondrial complex analysis in dek2 immature kernels showed severe deficiency of complex I assembly. Dramatically up-regulated expression of alternative oxidases (AOXs), transcriptome data, and TEM analysis results revealed that proper splicing of nad1 is critical for mitochondrial functions and inner cristaes morphology. This study indicated that Dek2 is a new PPR protein that affects the splicing of mitochondrial nad1 intron 1 and is required for mitochondrial function and kernel development.

Deciphering Transcriptional Programming during Pod and Seed Development Using RNA-Seq in Pigeonpea (Cajanus cajan)

Seed development is an important event in plant life cycle that has interested humankind since ages, especially in crops of economic importance. Pigeonpea is an important grain legume of the semi-arid tropics, used mainly for its protein rich seeds. In order to understand the transcriptional programming during the pod and seed development, RNA-seq data was generated from embryo sac from the day of anthesis (0 DAA), seed and pod wall (5, 10, 20 and 30 DAA) of pigeonpea variety "Asha" (ICPL 87119) using Illumina HiSeq 2500. About 684 million sequencing reads have been generated from nine samples, which resulted in the identification of 27,441 expressed genes after sequence analysis. These genes have been studied for their differentially expression, co-expression, temporal and spatial gene expression. We have also used the RNA-seq data to identify important seed-specific transcription factors, biological processes and associated pathways during seed development process in pigeonpea. The comprehensive gene expression study from flowering to mature pod development in pigeonpea would be crucial in identifying candidate genes involved in seed traits directly or indirectly related to yield and quality. The dataset will serve as an important resource for gene discovery and deciphering the molecular mechanisms underlying various seed related traits.

Transcriptome Dynamics during Maize Endosperm Development

The endosperm is a major organ of the seed that plays vital roles in determining seed weight and quality. However, genome-wide transcriptome patterns throughout maize endosperm development have not been comprehensively investigated to date. Accordingly, we performed a high-throughput RNA sequencing (RNA-seq) analysis of the maize endosperm transcriptome at 5, 10, 15 and 20 days after pollination (DAP). We found that more than 11,000 protein-coding genes underwent alternative splicing (AS) events during the four developmental stages studied. These genes were mainly involved in intracellular protein transport, signal transmission, cellular carbohydrate metabolism, cellular lipid metabolism, lipid biosynthesis, protein modification, histone modification, cellular amino acid metabolism, and DNA repair. Additionally, 7,633 genes, including 473 transcription factors (TFs), were differentially expressed among the four developmental stages. The differentially expressed TFs were from 50 families, including the bZIP, WRKY, GeBP and ARF families. Further analysis of the stage-specific TFs showed that binding, nucleus and ligand-dependent nuclear receptor activities might be important at 5 DAP, that immune responses, signalling, binding and lumen development are involved at 10 DAP, that protein metabolic processes and the cytoplasm might be important at 15 DAP, and that the responses to various stimuli are different at 20 DAP compared with the other developmental stages. This RNA-seq analysis provides novel, comprehensive insights into the transcriptome dynamics during early endosperm development in maize.

Maize endosperm-specific transcription factors O2 and PBF network the regulation of protein and starch synthesis

The maize endosperm-specific transcription factors opaque2 (O2) and prolamine-box binding factor (PBF) regulate storage protein zein genes. We show that they also control starch synthesis. The starch content in the PbfRNAi and o2 mutants was reduced by ∼5% and 11%, respectively, compared with normal genotypes. In the double-mutant PbfRNAi;o2, starch was decreased by 25%. Transcriptome analysis reveals that >1,000 genes were affected in each of the two mutants and in the double mutant; these genes were mainly enriched in sugar and protein metabolism. Pyruvate orthophosphate dikinase 1 and 2 (PPDKs) and starch synthase III (SSIII) are critical components in the starch biosynthetic enzyme complex. The expression of PPDK1, PPDK2, and SSIII and their protein levels are further reduced in the double mutants as compared with the single mutants. When the promoters of these genes were analyzed, we found a prolamine box and an O2 box that can be additively transactivated by PBF and O2. Starch synthase IIa (SSIIa, encoding another starch synthase for amylopectin) and starch branching enzyme 1 (SBEI, encoding one of the two main starch branching enzymes) are not directly regulated by PBF and O2, but their protein levels are significantly decreased in the o2 mutant and are further decreased in the double mutant, indicating that o2 and PbfRNAi may affect the levels of some other transcription factor(s) or mRNA regulatory factor(s) that in turn would affect the transcript and protein levels of SSIIa and SBEI These findings show that three important traits-nutritional quality, calories, and yield-are linked through the same transcription factors.

The green ash transcriptome and identification of genes responding to abiotic and biotic stresses

Background

To develop a set of transcriptome sequences to support research on environmental stress responses in green ash (Fraxinus pennsylvanica), we undertook deep RNA sequencing of green ash tissues under various stress treatments. The treatments, including emerald ash borer (EAB) feeding, heat, drought, cold and ozone, were selected to mimic the increasing threats of climate change and invasive pests faced by green ash across its native habitat.

Results

We report the generation and assembly of RNA sequences from 55 green ash samples into 107,611 putative unique transcripts (PUTs). 52,899 open reading frames were identified. Functional annotation of the PUTs by comparison to the Uniprot protein database identified matches for 63 % of transcripts and for 98 % of transcripts with ORFs. Further functional annotation identified conserved protein domains and assigned gene ontology terms to the PUTs. Examination of transcript expression across different RNA libraries revealed that expression patterns clustered based on tissues regardless of stress treatment. The transcripts from stress treatments were further examined to identify differential expression. Tens to hundreds of differentially expressed PUTs were identified for each stress treatment. A set of 109 PUTs were found to be consistently up or down regulated across three or more different stress treatments, representing basal stress response candidate genes in green ash. In addition, 1956 simple sequence repeats were identified in the PUTs, of which we identified 465 high quality DNA markers and designed flanking PCR primers.

Conclusions

North American native ash trees have suffered extensive mortality due to EAB infestation, creating a need to breed or select for resistant green ash genotypes. Stress from climate change is an additional concern for longevity of native ash populations. The use of genomics could accelerate management efforts. The green ash transcriptome we have developed provides important sequence information, genetic markers and stress-response candidate genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3052-0) contains supplementary material, which is available to authorized users.

Deciphering and modifying LAFL transcriptional regulatory network in seed for improving yield and quality of storage compounds

Increasing yield and quality of seed storage compounds in a sustainable way is a key challenge for our societies. Genome-wide analyses conducted in both monocot and dicot angiosperms emphasized drastic transcriptional switches that occur during seed development. In Arabidopsis thaliana, a reference species, genetic and molecular analyses have demonstrated the key role of LAFL (LEC1, ABI3, FUS3, and LEC2) transcription factors (TFs), in controlling gene expression programs essential to accomplish seed maturation and the accumulation of storage compounds. Here, we summarize recent progress obtained in the characterization of these LAFL proteins, their regulation, partners and target genes. Moreover, we illustrate how these evolutionary conserved TFs can be used to engineer new crops with altered seed compositions and point out the current limitations. Last, we discuss about the interest of investigating further the environmental and epigenetic regulation of this network for the coming years.

Identification of line-specific strategies for improving carotenoid production in synthetic maize through data-driven mathematical modeling

Plant synthetic biology is still in its infancy. However, synthetic biology approaches have been used to manipulate and improve the nutritional and health value of staple food crops such as rice, potato and maize. With current technologies, production yields of the synthetic nutrients are a result of trial and error, and systematic rational strategies to optimize those yields are still lacking. Here, we present a workflow that combines gene expression and quantitative metabolomics with mathematical modeling to identify strategies for increasing production yields of nutritionally important carotenoids in the seed endosperm synthesized through alternative biosynthetic pathways in synthetic lines of white maize, which is normally devoid of carotenoids. Quantitative metabolomics and gene expression data are used to create and fit parameters of mathematical models that are specific to four independent maize lines. Sensitivity analysis and simulation of each model is used to predict which gene activities should be further engineered in order to increase production yields for carotenoid accumulation in each line. Some of these predictions (e.g. increasing Zmlycb/Gllycb will increase accumulated β-carotenes) are valid across the four maize lines and consistent with experimental observations in other systems. Other predictions are line specific. The workflow is adaptable to any other biological system for which appropriate quantitative information is available. Furthermore, we validate some of the predictions using experimental data from additional synthetic maize lines for which no models were developed.

Genetic Diversity and Molecular Evolution of a Violaxanthin De-epoxidase Gene in Maize

Violaxanthin de-epoxidase (VDE) has a critical role in the carotenoid biosynthesis pathway, which is involved in protecting the photosynthesis apparatus from damage caused by excessive light. Here, a VDE gene in maize, ZmVDE1, was cloned and shown to have functional domains in common with the gramineous VDE protein. Candidate gene association analysis indicated that no polymorphic sites in ZmVDE1 were significant association with any of the examined carotenoid-related traits at P = 0.05 in an association panel containing 155 maize inbred lines. Nucleotide diversity analysis of VDE1 in maize and teosinte indicated that its exon had less genetic variation, consistent with the conserved function of VDE1 in plants. In addition, dramatically reduced nucleotide diversity, fewer haplotypes and a significantly negative parameter deviation for Tajima’s D test of ZmVDE1 in maize and teosinte suggested that a potential selective force had acted across the ZmVDE1 locus. We further identified a 4.2 Mb selective sweep with low recombination surrounding the ZmVDE1 locus that resulted in severely reduced nucleotide diversity on chromosome 2. Collectively, natural selection and the conserved domains of ZmVDE1 might show an important role in the xanthophyll cycle of the carotenoid biosynthesis pathway.

MODEM: multi-omics data envelopment and mining in maize

MODEM is a comprehensive database of maize multidimensional omics data, including genomic, transcriptomic, metabolic and phenotypic information from the cellular to individual plant level. This initial release contains approximately 1.06 M high quality SNPs for 508 diverse inbred lines obtained by combining variations from RNA sequencing on whole kernels (15 days after pollination) of 368 lines and a 50 K array for all 508 individuals. As all of these data were derived from the same diverse panel of lines, the database also allows various types of genetic mapping (including characterization of phenotypic QTLs, pQTLs; expression QTLs, eQTLs and metabolic QTLs, mQTLs). MODEM is thus designed to promote a better understanding of maize genetic architecture and deep functional annotation of the complex maize genome (and potentially those of other crop plants) and to explore the genotype-phenotype relationships and regulation of maize kernel development at multiple scales, which is also comprehensive for developing novel methods. MODEM is additionally designed to link with other databases to make full use of current resources, and it provides visualization tools for easy browsing. All of the original data and the related mapping results are freely available for easy query and download. This platform also provides helpful tools for general analyses and will be continually updated with additional materials, features and public data related to maize genetics or regulation as they become available.Database URL: (http://modem.hzau.edu.cn).

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

Starch is the major component in maize kernels, providing a stable carbohydrate source for humans and livestock as well as raw material for the biofuel industry. Increasing maize kernel starch content will help meet industry demands and has the potential to increase overall yields. We developed a pair of maize near-isogenic lines (NILs) with different alleles for a starch quantitative trait locus on chromosome 3 (qHS3), resulting in different kernel starch content. To investigate the candidate genes for qHS3 and elucidate their effects on starch metabolism, RNA-Seq was performed for the developing kernels of the NILs at 14 and 21 d after pollination (DAP). Analysis of genomic and transcriptomic data identified 76 genes with nonsynonymous single nucleotide polymorphisms and 384 differentially expressed genes (DEGs) in the introgressed fragment, including a hexokinase gene, ZmHXK3a, which catalyzes the conversion of glucose to glucose-6-phosphate and may play a key role in starch metabolism. The expression pattern of all DEGs in starch metabolism shows that altered expression of the candidate genes for qHS3 promoted starch synthesis, with positive consequences for kernel starch content. These results expand the current understanding of starch biosynthesis and accumulation in maize kernels and provide potential candidate genes to increase starch content.

Divergent Transactivation of Maize Storage Protein Zein Genes by the Transcription Factors Opaque2 and OHPs

Maize transcription factors (TFs) opaque2 (O2) and the O2 heterodimerizing proteins (OHP1 and OHP2) originated from an ancient segmental duplication. The 22-kDa (z1C) and 19-kDa (z1A, z1B, and z1D) α-zeins are the most abundant storage proteins in maize endosperm. O2 is known to regulate α-zein gene expression, but its target motifs in the 19-kDa α-zein gene promoters have not been identified. The mechanisms underlying the regulation of α-zein genes by these TFs are also not well understood. In this study, we found that the O2 binding motifs in the α-zein gene promoters are quite flexible, with ACGT being present in the z1C and z1A promoters and a variant, ACAT, being present in the z1B and z1D promoters. OHPs recognized and transactivated all of the α-zein promoters, although to much lower levels than did O2. In the presence of O2, the suppression of OHPs did not cause a significant reduction in the transcription of α-zein genes, but in the absence of O2, OHPs were critical for the expression of residual levels of α-zeins. These findings demonstrated that O2 is the primary TF and that OHPs function as minor TFs in this process. This relationship is the converse of that involved in 27-kDa γ-zein gene regulation, indicating that the specificities of O2 and the OHPs for regulating zein genes diverged after gene duplication. The prolamine-box binding factor by itself has limited transactivation activity, but it promotes the binding of O2 to O2 motifs, resulting in the synergistic transactivation of α-zein genes.