Genome-wide analysis of gene expression profiles during the kernel development of maize (Zea mays L.)

Recent developments in multi-omics and breeding strategies for abiotic stress tolerance in maize (Zea mays L.)

High-throughput sequencing technologies (HSTs) have revolutionized crop breeding. The advent of these technologies has enabled the identification of beneficial quantitative trait loci (QTL), genes, and alleles for crop improvement. Climate change have made a significant effect on the global maize yield. To date, the well-known omic approaches such as genomics, transcriptomics, proteomics, and metabolomics are being incorporated in maize breeding studies. These approaches have identified novel biological markers that are being utilized for maize improvement against various abiotic stresses. This review discusses the current information on the morpho-physiological and molecular mechanism of abiotic stress tolerance in maize. The utilization of omics approaches to improve abiotic stress tolerance in maize is highlighted. As compared to single approach, the integration of multi-omics offers a great potential in addressing the challenges of abiotic stresses of maize productivity.

Comparative Transcriptome Analysis Reveals Mechanisms of Folate Accumulation in Maize Grains

Previously, the complexity of folate accumulation in the early stages of maize kernel development has been reported, but the mechanisms of folate accumulation are unclear. Two maize inbred lines, DAN3130 and JI63, with different patterns of folate accumulation and different total folate contents in mature kernels were used to investigate the transcriptional regulation of folate metabolism during late stages of kernel formation by comparative transcriptome analysis. The folate accumulation during DAP 24 to mature kernels could be controlled by circumjacent pathways of folate biosynthesis, such as pyruvate metabolism, glutamate metabolism, and serine/glycine metabolism. In addition, the folate variation between these two inbred lines was related to those genes among folate metabolism, such as genes in the pteridine branch, para-aminobenzoate branch, serine/tetrahydrofolate (THF)/5-methyltetrahydrofolate cycle, and the conversion of THF monoglutamate to THF polyglutamate. The findings provided insight into folate accumulation mechanisms during maize kernel formation to promote folate biofortification.

Gene Expression in Parthenogenic Maize Proembryos

Angiosperm plants reproduce both sexually and asexually (by apomixis). In apomictic plants, the embryo and endosperm develop without fertilization. Modern maize seems to have a broken apomixis-triggering mechanism, which still works in Tripsacum and in Tripsacum-maize hybrids. For the first time, maize lines characterized by pronounced and inheritable high-frequency maternal parthenogenesis were generated 40 years ago, but there are no data on gene expression in parthenogenic maize proembryos. Here we examined for the first time gene expression in parthenogenic proembryos isolated from unpollinated embryo sacs (ESs) of a parthenogenic maize line (AT-4). The DNA-methylation genes (dmt103, dmt105) and the genes coding for the chromatin-modifying enzymes (chr106, hdt104, hon101) were expressed much higher in parthenogenic proembryos than in unpollinated ESs. The expression of the fertilization-independent endosperm (fie1) genes was found for the first time in parthenogenic proembryos and unpollinated ESs. In parthenogenic proembryos, the Zm_fie2 gene was expressed up to two times higher than it was expressed in unpollinated ESs.

A proteomic analysis of grain yield-related traits in wheat

Grain yield, which is mainly contributed by tillering capacity as well as kernel number and weight, is the most important trait to plant breeders and agronomists. Label-free quantitative proteomics was used to analyse yield-contributing organs in wheat. These were leaf sample, tiller initiation, spike initiation, ovary and three successive kernel development stages at 5, 10 and 15 days after anthesis (DAA). We identified 3182 proteins across all samples. The largest number was obtained for spike initiation (1673), while the smallest was kernel sample at 15 DAA (709). Of the 3182 proteins, 296 of them were common to all seven organs. Organ-specific proteins ranged from 148 in ovary to 561 in spike initiation. When relative protein abundances were compared to that of leaf sample, 347 and 519 proteins were identified as differentially abundant in tiller initiation and spike initiation, respectively. When compared with ovary, 81, 35 and 96 proteins were identified as differentially abundant in kernels sampled at 5, 10 and 15 DAA, respectively. Our study indicated that two Argonaute proteins were solely expressed in spike initiation. Of the four expansin proteins detected, three of them were mainly expressed during the first 10 days of kernel development after anthesis. We also detected cell wall invertases and sucrose and starch synthases mainly during the kernel development period. The manipulation of these proteins could lead to increases in tillers, kernels per spike or final grain weight, and is worth exploring in future studies.

SNP-based mixed model association of growth- and yield-related traits in popcorn

The identification of the genes responsible for complex traits is highly promising to accelerate crop breeding, but such information is still limited for popcorn. Thus, in the present study, a mixed linear model-based association analysis (MLMA) was applied for six important popcorn traits: plant and ear height, 100-grain weight, popping expansion, grain yield and expanded popcorn volume per hectare. To this end, 196 plants of the open-pollinated popcorn population UENF-14 were sampled, selfed (S1), and then genotyped with a panel of 10,507 single nucleotide polymorphisms (SNPs) markers distributed throughout the genome. The six traits were studied under two environments [Campos dos Goytacazes-RJ (ENV1) and Itaocara-RJ (ENV2)] in an incomplete block design. Based on the phenotypic data of the S1 progenies and on the genetic characteristics of the parents, the MLMA was performed. Thereafter, genes annotated in the MaizeGDB platform were screened for potential linkage disequilibrium with the SNPs associated to the six evaluated traits. Overall, seven and eight genes were identified as associated with the traits in ENV1 and ENV2, respectively, and proteins encoded by these genes were evaluated for their function. The results obtained here contribute to increase knowledge on the genetic architecture of the six evaluated traits and might be used for marker-assisted selection in breeding programs.

De novo transcriptome assembly of the Chinese pearl barley, adlay, by full-length isoform and short-read RNA sequencing

Adlay (Coix lacryma-jobi) is a tropical grass that has long been used in traditional Chinese medicine and is known for its nutritional benefits. Recent studies have shown that vitamin E compounds in adlay protect against chronic diseases such as cancer and heart disease. However, the molecular basis of adlay's health benefits remains unknown. Here, we generated adlay gene sets by de novo transcriptome assembly using long-read isoform sequencing (Iso-Seq) and short-read RNA-Sequencing (RNA-Seq). The gene sets obtained from Iso-seq and RNA-seq contained 31,177 genes and 57,901 genes, respectively. We confirmed the validity of the assembled gene sets by experimentally analyzing the levels of prolamin and vitamin E biosynthesis-associated proteins in adlay plant tissues and seeds. We compared the screened adlay genes with known gene families from closely related plant species, such as rice, sorghum and maize. We also identified tissue-specific genes from the adlay leaf, root, and young and mature seed, and experimentally validated the differential expression of 12 randomly-selected genes. Our study of the adlay transcriptome will provide a valuable resource for genetic studies that can enhance adlay breeding programs in the future.

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.

Genome-wide trait-trait dynamics correlation study dissects the gene regulation pattern in maize kernels

BACKGROUND

Dissecting the genetic basis and regulatory mechanisms for the biosynthesis and accumulation of nutrients in maize could lead to the improved nutritional quality of this crop. Gene expression is regulated at the genomic, transcriptional, and post-transcriptional levels, all of which can produce diversity among traits. However, the expression of most genes connected with a particular trait usually does not have a direct association with the variation of that trait. In addition, expression profiles of genes involved in a single pathway may vary as the intrinsic cellular state changes. To work around these issues, we utilized a statistical method, liquid association (LA) to investigate the complex pattern of gene regulation in maize kernels.

RESULTS

We applied LA to the expression profiles of 28,769 genes to dissect dynamic trait-trait correlation patterns in maize kernels. Among the 1000 LA pairs (LAPs) with the largest LA scores, 686 LAPs were identified conditional correlation. We also identified 830 and 215 LA-scouting leaders based on the positive and negative LA scores, which were significantly enriched for some biological processes and molecular functions. Our analysis of the dynamic co-expression patterns in the carotene biosynthetic pathway clearly indicated the important role of lcyE, CYP97A, ZEP1, and VDE in this pathway, which may change the direction of carotene biosynthesis by controlling the influx and efflux of the substrate. The dynamic trait-trait correlation patterns between gene expression and oil concentration in the fatty acid metabolic pathway and its complex regulatory network were also assessed. 23 of 26 oil-associated genes were correlated with oil concentration conditioning on 580 LA-scoutinggenes, and 5% of these LA-scouting genes were annotated as enzymes in the oil metabolic pathway.

CONCLUSIONS

By focusing on the carotenoid and oil biosynthetic pathways in maize, we showed that a genome-wide LA analysis provides a novel and effective way to detect transcriptional regulatory relationships. This method will help us understand the biological role of maize kernel genes and will benefit maize breeding programs.

Signaling in Early Maize Kernel Development

Developing the next plant generation within the seed requires the coordination of complex programs driving pattern formation, growth, and differentiation of the three main seed compartments: the embryo (future plant), the endosperm (storage compartment), representing the two filial tissues, and the surrounding maternal tissues. This review focuses on the signaling pathways and molecular players involved in early maize kernel development. In the 2 weeks following pollination, functional tissues are shaped from single cells, readying the kernel for filling with storage compounds. Although the overall picture of the signaling pathways regulating embryo and endosperm development remains fragmentary, several types of molecular actors, such as hormones, sugars, or peptides, have been shown to be involved in particular aspects of these developmental processes. These molecular actors are likely to be components of signaling pathways that lead to transcriptional programming mediated by transcriptional factors. Through the integrated action of these components, multiple types of information received by cells or tissues lead to the correct differentiation and patterning of kernel compartments. In this review, recent advances regarding the four types of molecular actors (hormones, sugars, peptides/receptors, and transcription factors) involved in early maize development are presented.

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.

Genome-Wide Identification and Characterization of microRNAs in Developing Grains of Zea mays L

The development and maturation of maize kernel involves meticulous and fine gene regulation at transcriptional and post-transcriptional levels, and miRNAs play important roles during this process. Although a number of miRNAs have been identified in maize seed, the ones involved in the early development of grains and in different lines of maize have not been well studied. Here, we profiled four small RNA libraries, each constructed from groups of immature grains of Zea mays inbred line Chang 7–2 collected 4–6, 7–9, 12–14, and 18–23 days after pollination (DAP). A total of 40 known (containing 111 unique miRNAs) and 162 novel (containing 196 unique miRNA candidates) miRNA families were identified. For conserved and novel miRNAs with over 100 total reads, 44% had higher accumulation before the 9^th DAP, especially miR166 family members. 42% of miRNAs had highest accumulation during 12–14 DAP (which is the transition stage from embryogenesis to nutrient storage). Only 14% of miRNAs had higher expression 18–23 DAP. Prediction of potential targets of all miRNAs showed that 165 miRNA families had 377 target genes. For miR164 and miR166, we showed that the transcriptional levels of their target genes were significantly decreased when co-expressed with their cognate miRNA precursors in vivo. Further analysis shows miR159, miR164, miR166, miR171, miR390, miR399, and miR529 families have putative roles in the embryogenesis of maize grain development by participating in transcriptional regulation and morphogenesis, while miR167 and miR528 families participate in metabolism process and stress response during nutrient storage. Our study is the first to present an integrated dynamic expression pattern of miRNAs during maize kernel formation and maturation.

Asymmetric transcriptomic signatures between the cob and florets in the maize ear under optimal- and low-nitrogen conditions at silking, and functional characterization of amino acid transporters ZmAAP4 and ZmVAAT3

Coordinated functioning of the cob and florets of the maize ear confers grain yield. The cob is critical for carbon partitioning and assimilated nitrogen (N) supply for grain development. However, molecular recognition of the cob and peripheral florets, characterization of genes mediating translocation of N assimilates, and responses of these two tissues to low N (LN) remain elusive. Transcriptional profiling of the ear of a maize hybrid at silking in the field revealed 1864 differentially expressed genes between the cob and florets, with 1314 genes up-regulated in the cob and 550 genes up-regulated in florets. The cob was characterized by striking enrichment of genes that are involved in carbon/N transport and metabolism, consistent with the physiological role of the cob in carbon/N storage and transfer during ear development. The florets were characterized by enrichment of hormone signalling components and development related genes. We next examined the response of the cob and florets to LN stress. LN caused differential expression of 588 genes in the cob and only 195 genes in the florets, indicating that the cob dominated the response of the ear to LN at the transcriptional level. LN caused comprehensive alterations such as carbon/N metabolism or partitioning, hormone signalling and protein phosphorylation in terms of gene expression in the cob and/or florets. Fourteen genes responsive specifically to LN provided potential molecular markers for N-efficient maize breeding. We further functionally characterized two newly identified broad-spectrum amino acid transporters, ZmAAP4 and ZmVAAT3, that showed distinct expression patterns in the cob and florets and potentially important roles in amino-N mobilization in the ear. While both proteins could transport various amino acids into yeast or Arabidopsis cells, ZmAAP4 appeared to have higher efficiencies than ZmVAAT3 in transporting seven out of 22 examined amino acids.

Genome-wide identification and transcriptional analysis of folate metabolism-related genes in maize kernels

BACKGROUND

Maize is a major staple food crop globally and contains various concentrations of vitamins. Folates are essential water-soluble B-vitamins that play an important role as one-carbon (C1) donors and acceptors in organisms. To gain an understanding of folate metabolism in maize, we performed an intensive in silico analysis to screen for genes involved in folate metabolism using publicly available databases, followed by examination of the transcript expression patterns and profiling of the folate derivatives in the kernels of two maize inbred lines.

RESULTS

A total of 36 candidate genes corresponding to 16 folate metabolism-related enzymes were identified. The maize genome contains all the enzymes required for folate and C1 metabolism, characterized by highly conserved functional domains across all the other species investigated. Phylogenetic analysis revealed that these enzymes in maize are conserved throughout evolution and have a high level of similarity with those in sorghum and millet. The LC-MS analyses of two maize inbred lines demonstrated that 5-methyltetrahydrofolate was the major form of folate derivative in young seeds, while 5-formyltetrahydrofolate in mature seeds. Most of the genes involved in folate and C1 metabolism exhibited similar transcriptional expression patterns between these two maize lines, with the highest transcript abundance detected on day after pollination (DAP) 6 and the decreased transcript abundance on DAP 12 and 18. Compared with the seeds on DAP 30, 5-methyltetrahydrofolate was decreased and 5-formyltetrahydrofolate was increased sharply in the mature dry seeds.

CONCLUSIONS

The enzymes involved in folate and C1 metabolism are conserved between maize and other plant species. Folate and C1 metabolism is active in young developing maize seeds at transcriptional levels.

Genome-wide analysis reveals phytohormone action during cassava storage root initiation

Development of storage roots is a process associated with a phase change from cell division and elongation to radial growth and accumulation of massive amounts of reserve substances such as starch. Knowledge of the regulation of cassava storage root formation has accumulated over time; however, gene regulation during the initiation and early stage of storage root development is still poorly understood. In this study, transcription profiling of fibrous, intermediate and storage roots at eight weeks old were investigated using a 60-mer-oligo microarray. Transcription and gene expression were found to be the key regulating processes during the transition stage from fibrous to intermediate roots, while homeostasis and signal transduction influenced regulation from intermediate roots to storage roots. Clustering analysis of significant genes and transcription factors (TF) indicated that a number of phytohormone-related TF were differentially expressed; therefore, phytohormone-related genes were assembled into a network of correlative nodes. We propose a model showing the relationship between KNOX1 and phytohormones during storage root initiation. Exogeneous treatment of phytohormones N (6) -benzylaminopurine and 1-Naphthaleneacetic acid were used to induce the storage root initiation stage and to investigate expression patterns of the genes involved in storage root initiation. The results support the hypothesis that phytohormones are acting in concert to regulate the onset of cassava storage root development. Moreover, MeAGL20 is a factor that might play an important role at the onset of storage root initiation when the root tip becomes swollen.

Role of B3 domain transcription factors of the AFL family in maize kernel filling

In the dicot Arabidopsis thaliana, the B3 transcription factors, ABA-INSENSITIVE 3 (ABI3), FUSCA 3 (FUS3) and LEAFY COTYLEDON 2 (LEC2) are key regulators of seed maturation. This raises the question of the role of ABI3/FUS3/LEC2 (AFL) proteins in cereals, where not only the embryo but also the persistent endosperm accumulates reserve substances. Among the five ZmAFL genes identified in the maize genome, ZmAFL2 and ZmAFL3/ZmVp1 closely resemble FUS3 and ABI3, respectively, in terms of their sequences, domain structure and gene activity profiles. Of the three genes that fall into the LEC2 phylogenetic sub-clade, ZmAFL5 and ZmAFL6 have constitutive gene activity, whereas ZmAFL4, like LEC2, has preferential gene activity in pollen and seed, although its seed gene activity is restricted to the endosperm during reserve accumulation. Knock down of ZmAFL4 gene activity perturbs carbon metabolism and reduces starch content in the developing endosperm at 20 DAP. ZmAFL4 and ZmAFL3/ZmVp1 trans-activate a maize oleosin promoter in a heterologous moss system. In conclusion our results suggest, based on gene activity profiles, that the functions of FUS3 and ABI3 could be conserved between dicot and monocot species. In contrast, LEC2 function may have partially diverged in cereals where our findings provide first evidence of the specialization of ZmAFL4 for roles in the endosperm.

UNCLES: method for the identification of genes differentially consistently co-expressed in a specific subset of datasets

Background

Collective analysis of the increasingly emerging gene expression datasets are required. The recently proposed binarisation of consensus partition matrices (Bi-CoPaM) method can combine clustering results from multiple datasets to identify the subsets of genes which are consistently co-expressed in all of the provided datasets in a tuneable manner. However, results validation and parameter setting are issues that complicate the design of such methods. Moreover, although it is a common practice to test methods by application to synthetic datasets, the mathematical models used to synthesise such datasets are usually based on approximations which may not always be sufficiently representative of real datasets.

Results

Here, we propose an unsupervised method for the unification of clustering results from multiple datasets using external specifications (UNCLES). This method has the ability to identify the subsets of genes consistently co-expressed in a subset of datasets while being poorly co-expressed in another subset of datasets, and to identify the subsets of genes consistently co-expressed in all given datasets. We also propose the M-N scatter plots validation technique and adopt it to set the parameters of UNCLES, such as the number of clusters, automatically. Additionally, we propose an approach for the synthesis of gene expression datasets using real data profiles in a way which combines the ground-truth-knowledge of synthetic data and the realistic expression values of real data, and therefore overcomes the problem of faithfulness of synthetic expression data modelling. By application to those datasets, we validate UNCLES while comparing it with other conventional clustering methods, and of particular relevance, biclustering methods. We further validate UNCLES by application to a set of 14 real genome-wide yeast datasets as it produces focused clusters that conform well to known biological facts. Furthermore, in-silico-based hypotheses regarding the function of a few previously unknown genes in those focused clusters are drawn.

Conclusions

The UNCLES method, the M-N scatter plots technique, and the expression data synthesis approach will have wide application for the comprehensive analysis of genomic and other sources of multiple complex biological datasets. Moreover, the derived in-silico-based biological hypotheses represent subjects for future functional studies.

Electronic supplementary material

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

The genetic architecture of maize (Zea mays L.) kernel weight determination

Individual kernel weight is an important trait for maize yield determination. We have identified genomic regions controlling this trait by using the B73xMo17 population; however, the effect of genetic background on control of this complex trait and its physiological components is not yet known. The objective of this study was to understand how genetic background affected our previous results. Two nested stable recombinant inbred line populations (N209xMo17 and R18xMo17) were designed for this purpose. A total of 408 recombinant inbred lines were genotyped and phenotyped at two environments for kernel weight and five other traits related to kernel growth and development. All traits showed very high and significant (P < 0.001) phenotypic variability and medium-to-high heritability (0.60-0.90). When N209xMo17 and R18xMo17 were analyzed separately, a total of 23 environmentally stable quantitative trait loci (QTL) and five epistatic interactions were detected for N209xMo17. For R18xMo17, 59 environmentally stable QTL and 17 epistatic interactions were detected. A joint analysis detected 14 stable QTL regardless of the genetic background. Between 57 and 83% of detected QTL were population specific, denoting medium-to-high genetic background effects. This percentage was dependent on the trait. A meta-analysis including our previous B73xMo17 results identified five relevant genomic regions deserving further characterization. In summary, our grain filling traits were dominated by small additive QTL with several epistatic and few environmental interactions and medium-to-high genetic background effects. This study demonstrates that the number of detected QTL and additive effects for different physiologically related grain filling traits need to be understood relative to the specific germplasm.

Dynamic transcriptome landscape of maize embryo and endosperm development

Maize (Zea mays) is an excellent cereal model for research on seed development because of its relatively large size for both embryo and endosperm. Despite the importance of seed in agriculture, the genome-wide transcriptome pattern throughout seed development has not been well characterized. Using high-throughput RNA sequencing, we developed a spatiotemporal transcriptome atlas of B73 maize seed development based on 53 samples from fertilization to maturity for embryo, endosperm, and whole seed tissues. A total of 26,105 genes were found to be involved in programming seed development, including 1,614 transcription factors. Global comparisons of gene expression highlighted the fundamental transcriptomic reprogramming and the phases of development. Coexpression analysis provided further insight into the dynamic reprogramming of the transcriptome by revealing functional transitions during maturation. Combined with the published nonseed high-throughput RNA sequencing data, we identified 91 transcription factors and 1,167 other seed-specific genes, which should help elucidate key mechanisms and regulatory networks that underlie seed development. In addition, correlation of gene expression with the pattern of DNA methylation revealed that hypomethylation of the gene body region should be an important factor for the expressional activation of seed-specific genes, especially for extremely highly expressed genes such as zeins. This study provides a valuable resource for understanding the genetic control of seed development of monocotyledon plants.

RNA sequencing reveals the complex regulatory network in the maize kernel

RNA sequencing can simultaneously identify exonic polymorphisms and quantitate gene expression. Here we report RNA sequencing of developing maize kernels from 368 inbred lines producing 25.8 billion reads and 3.6 million single-nucleotide polymorphisms. Both the MaizeSNP50 BeadChip and the Sequenom MassArray iPLEX platforms confirm a subset of high-quality SNPs. Of these SNPs, we have mapped 931,484 to gene regions with a mean density of 40.3 SNPs per gene. The genome-wide association study identifies 16,408 expression quantitative trait loci. A two-step approach defines 95.1% of the eQTLs to a 10-kb region, and 67.7% of them include a single gene. The establishment of relationships between eQTLs and their targets reveals a large-scale gene regulatory network, which include the regulation of 31 zein and 16 key kernel genes. These results contribute to our understanding of kernel development and to the improvement of maize yield and nutritional quality.

Genome-wide analysis of coordinated transcript abundance during seed development in different Brassica rapa morphotypes

BACKGROUND

Brassica seeds are important as basic units of plant growth and sources of vegetable oil. Seed development is regulated by many dynamic metabolic processes controlled by complex networks of spatially and temporally expressed genes. We conducted a global microarray gene co-expression analysis by measuring transcript abundance of developing seeds from two diverse B. rapa morphotypes: a pak choi (leafy-type) and a yellow sarson (oil-type), and two of their doubled haploid (DH) progenies, (1) to study the timing of metabolic processes in developing seeds, (2) to explore the major transcriptional differences in developing seeds of the two morphotypes, and (3) to identify the optimum stage for a genetical genomics study in B. rapa seed.

RESULTS

Seed developmental stages were similar in developing seeds of pak choi and yellow sarson of B. rapa; however, the colour of embryo and seed coat differed among these two morphotypes. In this study, most transcriptional changes occurred between 25 and 35 DAP, which shows that the timing of seed developmental processes in B. rapa is at later developmental stages than in the related species B. napus. Using a Weighted Gene Co-expression Network Analysis (WGCNA), we identified 47 "gene modules", of which 27 showed a significant association with temporal and/or genotypic variation. An additional hierarchical cluster analysis identified broad spectra of gene expression patterns during seed development. The predominant variation in gene expression was according to developmental stages rather than morphotype differences. Since lipids are the major storage compounds of Brassica seeds, we investigated in more detail the regulation of lipid metabolism. Four co-regulated gene clusters were identified with 17 putative cis-regulatory elements predicted in their 1000 bp upstream region, either specific or common to different lipid metabolic pathways.

CONCLUSIONS

This is the first study of genome-wide profiling of transcript abundance during seed development in B. rapa. The identification of key physiological events, major expression patterns, and putative cis-regulatory elements provides useful information to construct gene regulatory networks in B. rapa developing seeds and provides a starting point for a genetical genomics study of seed quality traits.

Kinetics and functional diversity among the five members of the NADP-malic enzyme family from Zea mays, a C4 species

NADP-malic enzyme (NADP-ME) is involved in different metabolic pathways in several organisms due to the relevant physiological functions of the substrates and products of its reaction. In plants, it is one of the most important proteins that were recruited to fulfil key roles in C4 photosynthesis. Recent advances in genomics allowed the characterization of the complete set of NADP-ME genes from some C3 species, as Arabidopsis thaliana and Oryza sativa; however, the characterization of the complete NADP-ME family from a C4 species has not been performed yet. In this study, while taking advantage of the complete Zea mays genome sequence recently released, the characterization of the whole NADP-ME family is presented. The maize NADP-ME family is composed of five genes, two encoding plastidic NADP-MEs (ZmC4- and ZmnonC4-NADP-ME), and three cytosolic enzymes (Zmcyt1-, Zmcyt2-, and Zmcyt3-NADP-ME). The results presented clearly show that each maize NADP-ME displays particular organ distribution, response to stress stimuli, and differential biochemical properties. Phylogenetic footprinting studies performed with the NADP-MEs from several grasses, indicate that four members of the maize NADP-ME family share conserved transcription factor binding motifs with their orthologs, indicating conserved physiological functions for these genes in monocots. Based on the results obtained in this study, and considering the biochemical plasticity shown by the NADP-ME, it is discussed the relevance of the presence of a multigene family, in which each member encodes an isoform with particular biochemical properties, in the evolution of the C4 NADP-ME, improved to fulfil the requirements for an efficient C4 mechanism.

Seed-development programs: a systems biology-based comparison between dicots and monocots

Seeds develop differently in dicots and monocots, especially with respect to the major storage organs. High-resolution transcriptome data have provided the first insights into the molecular networks and pathway interactions that function during the development of individual seed compartments. Here, we review mainly recent data obtained by systems biology-based approaches, which have allowed researchers to construct and model complex metabolic networks and fluxes and identify key limiting steps in seed development. Comparative coexpression network analyses define evolutionarily conservative (FUS3/ABI3/LEC1) and divergent (LEC2) networks in dicots and monocots. Finally, we discuss the determination of seed size--an important yield-related characteristic--as mediated by a number of processes (maternal and epigenetic factors, fine-tuned regulation of cell death in distinct seed compartments, and endosperm growth) and underlying genes defined through mutant analyses. Altogether, systems approaches can make important contributions toward a more complete and holistic knowledge of seed biology and thus support strategies for knowledge-based molecular breeding.

Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Carefully analyzed expression profiles can serve as a valuable reference for deciphering gene functions. We exploited the potential of whole genome microarrays to measure the spatial and temporal expression profiles of rice genes in 19 stages of vegetative and reproductive development. We could verify expression of 22,980 genes in at least one of the tissues. Differential expression analysis with respect to five vegetative tissues and preceding stages of development revealed reproductive stage-preferential/-specific genes. By using subtractive logic, we identified 354 and 456 genes expressing specifically during panicle and seed development, respectively. The metabolic/hormonal pathways and transcription factor families playing key role in reproductive development were elucidated after overlaying the expression data on the public databases and manually curated list of transcription factors, respectively. During floral meristem differentiation (P1) and male meiosis (P3), the genes involved in jasmonic acid and phenylpropanoid biosynthesis were significantly upregulated. P6 stage of panicle, containing mature gametophytes, exhibited enrichment of transcripts involved in homogalacturonon degradation. Genes regulating auxin biosynthesis were induced during early seed development. We validated the stage-specificity of regulatory regions of three panicle-specific genes, OsAGO3, OsSub42, and RTS, and an early seed-specific gene, XYH, in transgenic rice. The data generated here provides a snapshot of the underlying complexity of the gene networks regulating rice reproductive development.

Controlling lipid accumulation in cereal grains

Plant oils have so far been mostly directed toward food and feed production. Nowadays however, these oils are more and more used as competitive alternatives to mineral hydrocarbon-based products. This increasing demand for vegetable oils has led to a renewed interest in elucidating the metabolism of storage lipids and its regulation in various plant systems. Cereal grains store carbon in the form of starch in a large endosperm and as oil in an embryo of limited size. Complementary studies on kernel development and metabolism have paved the way for breeding or engineering new varieties with higher grain oil content. This could be achieved either by increasing the relative proportion of the oil-rich embryo within the grain, or by enhancing oil synthesis and accumulation in embryonic structures. For instance, diacylglycerol acyltransferase (DGAT) that catalyses the ultimate reaction in the biosynthesis of triacylglycerol appears to be a promising target for increasing oil content in maize embryos. Similarly, over-expression of the maize transcriptional regulators ZmLEAFY COTYLEDON1 and ZmWRINKLED1 efficiently stimulates oil accumulation in the kernels of transgenic lines. Redirecting carbon from starch to oil in the endosperm, though not yet realized, is discussed.

Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Carefully analyzed expression profiles can serve as a valuable reference for deciphering gene functions. We exploited the potential of whole genome microarrays to measure the spatial and temporal expression profiles of rice genes in 19 stages of vegetative and reproductive development. We could verify expression of 22,980 genes in at least one of the tissues. Differential expression analysis with respect to five vegetative tissues and preceding stages of development revealed reproductive stage-preferential/-specific genes. By using subtractive logic, we identified 354 and 456 genes expressing specifically during panicle and seed development, respectively. The metabolic/hormonal pathways and transcription factor families playing key role in reproductive development were elucidated after overlaying the expression data on the public databases and manually curated list of transcription factors, respectively. During floral meristem differentiation (P1) and male meiosis (P3), the genes involved in jasmonic acid and phenylpropanoid biosynthesis were significantly upregulated. P6 stage of panicle, containing mature gametophytes, exhibited enrichment of transcripts involved in homogalacturonon degradation. Genes regulating auxin biosynthesis were induced during early seed development. We validated the stage-specificity of regulatory regions of three panicle-specific genes, OsAGO3, OsSub42, and RTS, and an early seed-specific gene, XYH, in transgenic rice. The data generated here provides a snapshot of the underlying complexity of the gene networks regulating rice reproductive development.

Transcriptional profile analysis of E3 ligase and hormone-related genes expressed during wheat grain development

BACKGROUND

Wheat grains are an important source of food, stock feed and raw materials for industry, but current production levels cannot meet world needs. Elucidation of the molecular mechanisms underlying wheat grain development will contribute valuable information to improving wheat cultivation. One of the most important mechanisms implicated in plant developmental processes is the ubiquitin-proteasome system (UPS). Among the different roles of the UPS, it is clear that it is essential to hormone signaling. In particular, E3 ubiquitin ligases of the UPS have been shown to play critical roles in hormone perception and signal transduction.

RESULTS

A NimbleGen microarray containing 39,179 UniGenes was used to study the kinetics of gene expression during wheat grain development from the early stages of cell division to the mid-grain filling stage. By comparing 11 consecutive time-points, 9284 differentially expressed genes were identified and annotated during this study. A comparison of the temporal profiles of these genes revealed dynamic transcript accumulation profiles with major reprogramming events that occurred during the time intervals of 80-120 and 220-240°Cdays. The list of the genes expressed differentially during these transitions were identified and annotated. Emphasis was placed on E3 ligase and hormone-related genes. In total, 173 E3 ligase coding genes and 126 hormone-related genes were differentially expressed during the cell division and grain filling stages, with each family displaying a different expression profile.

CONCLUSIONS

The differential expression of genes involved in the UPS and plant hormone pathways suggests that phytohormones and UPS crosstalk might play a critical role in the wheat grain developmental process. Some E3 ligase and hormone-related genes seem to be up- or down-regulated during the early and late stages of the grain development.

Genome-wide analysis of the complex transcriptional networks of rice developing seeds

Background

The development of rice (Oryza sativa) seed is closely associated with assimilates storage and plant yield, and is fine controlled by complex regulatory networks. Exhaustive transcriptome analysis of developing rice embryo and endosperm will help to characterize the genes possibly involved in the regulation of seed development and provide clues of yield and quality improvement.

Principal Findings

Our analysis showed that genes involved in metabolism regulation, hormone response and cellular organization processes are predominantly expressed during rice development. Interestingly, 191 transcription factor (TF)-encoding genes are predominantly expressed in seed and 59 TFs are regulated during seed development, some of which are homologs of seed-specific TFs or regulators of Arabidopsis seed development. Gene co-expression network analysis showed these TFs associated with multiple cellular and metabolism pathways, indicating a complex regulation of rice seed development. Further, by employing a cold-resistant cultivar Hanfeng (HF), genome-wide analyses of seed transcriptome at normal and low temperature reveal that rice seed is sensitive to low temperature at early stage and many genes associated with seed development are down-regulated by low temperature, indicating that the delayed development of rice seed by low temperature is mainly caused by the inhibition of the development-related genes. The transcriptional response of seed and seedling to low temperature is different, and the differential expressions of genes in signaling and metabolism pathways may contribute to the chilling tolerance of HF during seed development.

Conclusions

These results provide informative clues and will significantly improve the understanding of rice seed development regulation and the mechanism of cold response in rice seed.

Genome-wide transcriptome analysis of gametophyte development in Physcomitrella patens

BACKGROUND

Regulation of gene expression plays a pivotal role in controlling the development of multicellular plants. To explore the molecular mechanism of plant developmental-stage transition and cell-fate determination, a genome-wide analysis was undertaken of sequential developmental time-points and individual tissue types in the model moss Physcomitrella patens because of the short life cycle and relative structural simplicity of this plant.

RESULTS

Gene expression was analyzed by digital gene expression tag profiling of samples taken from P. patens protonema at 3, 14 and 24 days, and from leafy shoot tissues at 30 days, after protoplast isolation, and from 14-day-old caulonemal and chloronemal tissues. In total, 4333 genes were identified as differentially displayed. Among these genes, 4129 were developmental-stage specific and 423 were preferentially expressed in either chloronemal or caulonemal tissues. Most of the differentially displayed genes were assigned to functions in organic substance and energy metabolism or macromolecule biosynthetic and catabolic processes based on gene ontology descriptions. In addition, some regulatory genes identified as candidates might be involved in controlling the developmental-stage transition and cell differentiation, namely MYB-like, HB-8, AL3, zinc finger family proteins, bHLH superfamily, GATA superfamily, GATA and bZIP transcription factors, protein kinases, genes related to protein/amino acid methylation, and auxin, ethylene, and cytokinin signaling pathways.

CONCLUSIONS

These genes that show highly dynamic changes in expression during development in P. patens are potential targets for further functional characterization and evolutionary developmental biology studies.

Expression profile of maize (Zea mays L.) embryonic axes during germination: translational regulation of ribosomal protein mRNAs

Seed germination is a critical developmental period for plant propagation. Information regarding gene expression within this important period is relevant for understanding the main biochemical processes required for successful germination, particularly in maize, one of the most important cereals in the world. The present research focuses on the global microarray analysis of differential gene expression between quiescent and germinated maize embryo stages. This analysis revealed that a large number of mRNAs stored in the quiescent embryonic axes (QEAs) were differentially regulated during germination in the 24 h germinated embryonic axes (GEAs). These genes belong to 14 different functional categories and most of them correspond to metabolic processes, followed by transport, transcription and translation. Interestingly, the expression of mRNAs encoding ribosomal proteins [(r)-proteins], required for new ribosome formation during this fast-growing period, remains mostly unchanged throughout the germination process, suggesting that these genes are not regulated at the transcriptional level during this developmental period. To investigate this issue further, comparative microarray analyses on polysomal mRNAs from growth-stimulated and non-stimulated GEAs were performed. The results revealed that (r)-protein mRNAs accumulate to high levels in polysomes of the growth-stimulated tissues, indicating a translational control mechanism to account for the rapid (r)-protein synthesis observed within this period. Bioinformatic analysis of (r)-protein mRNAs showed that 5' TOP (tract of pyrimidines)-like sequences are present only in the 5'-untranslated region set of up-regulated (r)-protein mRNAs. This overall approach to the germination process allows an in-depth view of molecular changes, enabling a broader understanding of the regulatory mechanisms that occur during this process.

Identification of differentially expressed genes at two key endosperm development stages using two maize inbreds with large and small grain and integration with detected QTL for grain weight

Maize endosperm accounts for more than 80% of the grain weight. Cell division and grain filling are the two key stages for endosperm development. Previous studies showed that gene expression during differential stages in endosperm development is greatly different. However, information on systematic identification and characterization of the differentially expressed genes between the two stages are limited. In this study, suppression subtractive hybridization (SSH) was used to generate four subtracted cDNA libraries for the two stages using two maize inbreds with large and small grain. Totally, 4,784 differentially expressed sequence tags (ESTs) were sequenced and 902 were non-redundant, which consisted of 344 unique ESTs. Among them 192 had high sequence similarity to the GenBank entries and represent diverse of functional categories, such as metabolism, cell growth/division, transcription, signal transduction, protein destination/storage, protein synthesis and others. The expression patterns of 75.7% SSH-derived cDNAs were confirmed by reverse Northern blot and semi-quantitative reverse transcription polymerase chain reaction, and exhibited the similar results (75.0%). Genes differentially expressed between two key stages for the two inbreds were involved in diverse physiological process pathway, which might be responsible for the formation of grain weight. 43.8% (70 of the 160 unique ESTs) of the identified ESTs were assigned to 39 chromosome bins distributed over all ten maize chromosomes. Eleven ESTs were found to co-localize with previous detected QTLs for grain weight, which might be considered as the candidate genes of grain weight for further study.

An expression analysis of 57 transcription factors derived from ESTs of developing seeds in Maize (Zea mays)

Maize seeds are an important source of food, animal feed, and industrial raw materials. To understand global gene expression and regulation during maize seed development, a normalized cDNA library, covering most of the developmental stages of maize seeds, was constructed. Sequencing analysis of 10,848 randomly selected clones identified 6,630 unique ESTs. Among them, 57 putative transcription factors (TFs) were identified. The TFs belong to seven different super-families, specifically 17 Zinc-finger, 13 bZIP, 8 bHLH, 6 MADS, 7 MYB, 3 Homedomain, and 3 AP2/EREBP. The spatial and temporal expression of the TFs was analyzed by semi-quantitative RT-PCR with representative tissue types and seeds at different developmental stages, revealing their diverse expression patterns and expression levels. One-third (19) of the maize TFs was found their putative orthologs in Arabidopsis. Similar expression patterns were observed in both maize and Arabidopsis for the majority of orthologous pairs (15 out of 19), suggesting their conserved functions during seed development. In conclusion, the systematic analysis of maize seed TFs has provided valuable insight into transcriptional regulation during maize seed development.

Dissecting grain yield pathways and their interactions with grain dry matter content by a two-step correlation approach with maize seedling transcriptome

BACKGROUND

The importance of maize for human and animal nutrition, but also as a source for bio-energy is rapidly increasing. Maize yield is a quantitative trait controlled by many genes with small effects, spread throughout the genome. The precise location of the genes and the identity of the gene networks underlying maize grain yield is unknown. The objective of our study was to contribute to the knowledge of these genes and gene networks by transcription profiling with microarrays.

RESULTS

We assessed the grain yield and grain dry matter content (an indicator for early maturity) of 98 maize hybrids in multi-environment field trials. The gene expression in seedlings of the parental inbred lines, which have four different genetic backgrounds, was assessed with genome-scale oligonucleotide arrays. We identified genes associated with grain yield and grain dry matter content using a newly developed two-step correlation approach and found overlapping gene networks for both traits. The underlying metabolic pathways and biological processes were elucidated. Genes involved in sucrose degradation and glycolysis, as well as genes involved in cell expansion and endocycle were found to be associated with grain yield.

CONCLUSIONS

Our results indicate that the capability of providing energy and substrates, as well as expanding the cell at the seedling stage, highly influences the grain yield of hybrids. Knowledge of these genes underlying grain yield in maize can contribute to the development of new high yielding varieties.

Flux of transcript patterns during soybean seed development

BACKGROUND

To understand gene expression networks leading to functional properties of the soybean seed, we have undertaken a detailed examination of soybean seed development during the stages of major accumulation of oils, proteins, and starches, as well as the desiccating and mature stages, using microarrays consisting of up to 27,000 soybean cDNAs. A subset of these genes on a highly-repetitive 70-mer oligonucleotide microarray was also used to support the results.

RESULTS

It was discovered that genes related to cell growth and maintenance processes, as well as energy processes like photosynthesis, decreased in expression levels as the cotyledons approached the mature, dry stage. Genes involved with some storage proteins had their highest expression levels at the stage of highest fresh weight. However, genes encoding many transcription factors and DNA binding proteins showed higher expression levels in the desiccating and dry seeds than in most of the green stages.

CONCLUSIONS

Data on 27,000 cDNAs have been obtained over five stages of soybean development, including the stages of major accumulation of agronomically-important products, using two different types of microarrays. Of particular interest are the genes found to peak in expression at the desiccating and dry seed stages, such as those annotated as transcription factors, which may indicate the preparation of pathways that will be needed later in the early stages of imbibition and germination.

Epistatic interactions between Opaque2 transcriptional activator and its target gene CyPPDK1 control kernel trait variation in maize

Association genetics is a powerful method to track gene polymorphisms responsible for phenotypic variation, since it takes advantage of existing collections and historical recombination to study the correlation between large genetic diversity and phenotypic variation. We used a collection of 375 maize (Zea mays ssp. mays) inbred lines representative of tropical, American, and European diversity, previously characterized for genome-wide neutral markers and population structure, to investigate the roles of two functionally related candidate genes, Opaque2 and CyPPDK1, on kernel quality traits. Opaque2 encodes a basic leucine zipper transcriptional activator specifically expressed during endosperm development that controls the transcription of many target genes, including CyPPDK1, which encodes a cytosolic pyruvate orthophosphate dikinase. Using statistical models that correct for population structure and individual kinship, Opaque2 polymorphism was found to be strongly associated with variation of the essential amino acid lysine. This effect could be due to the direct role of Opaque2 on either zein transcription, zeins being major storage proteins devoid of lysine, or lysine degradation through the activation of lysine ketoglutarate reductase. Moreover, we found that a polymorphism in the Opaque2 coding sequence and several polymorphisms in the CyPPDK1 promoter nonadditively interact to modify both lysine content and the protein-versus-starch balance, thus revealing the role in quantitative variation in plants of epistatic interactions between a transcriptional activator and one of its target genes.