An Expanded Maize Gene Expression Atlas based on RNA Sequencing and its Use to Explore Root Development

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

Genome-wide identification, expression profiling, and network analysis of AT-hook gene family in maize

AT-hook motif nuclear localized (AHL) genes have diverse but poorly understood biological functions. We identified and analyzed 37 AHL genes in maize. We also discovered four and one additional AHLs in rice and sorghum, respectively, besides those reported earlier. The maize AHLs were classified into two clades (A and B) and three distinct types (I, II, and III) as also reported in Arabidopsis. Phylogenetic and ortholog analyses showed that, while the evolutionary classification was conserved in plants, expansion of the AHL gene family in maize was accompanied with new biological functions. Gene structure analysis showed that, while all but one Type-I AHLs lacked an intron, origin of Type-II and Type-III AHLs was associated with the gain of introns suggesting evolutionarily distinct temporal and spatial expression patterns and, likely, neofunctionalization. Gene duplication analysis revealed that AHLs in maize expanded via dispersive duplication further supporting their functional diversity. To discern these functions, we analyzed 71 transcriptomes from diverse tissues and developmental stages of maize and classified AHLs into eight groups with distinct temporal/spatial expression profiles. Coexpression analysis implicated 5 AHLs and 33 novel genes in networks specific to endosperm, seed, root, leaf, and reproductive tissues indicating their role in the development of these organs. Major processes coregulated by AHLs include pollen development, drought response, senescence, and wound response. We also identified interactions of AHL proteins in coregulating important processes including stress response. These novel insights into the role of AHLs in plant development provide a platform for functional analyses in maize and related grasses.

A maize polygalacturonase functions as a suppressor of programmed cell death in plants

BACKGROUND

The hypersensitive defense response (HR) in plants is a fast, localized necrotic response around the point of pathogen ingress. HR is usually triggered by a pathogen recognition event mediated by a nucleotide-binding site, leucine-rich repeat (NLR) protein. The autoactive maize NLR gene Rp1-D21 confers a spontaneous HR response in the absence of pathogen recognition. Previous work identified a set of loci associated with variation in the strength of Rp1-D21-induced HR. A polygalacturonase gene homolog, here termed ZmPGH1, was identified as a possible causal gene at one of these loci on chromosome 7.

RESULTS

Expression of ZmPGH1 inhibited the HR-inducing activity of both Rp1-D21 and that of another autoactive NLR, RPM1(D505V), in a Nicotiana benthamiana transient expression assay system. Overexpression of ZmPGH1 in a transposon insertion line of maize was associated with suppression of chemically-induced programmed cell death and with suppression of HR induced by Rp1-D21 in maize plants grown in the field.

CONCLUSIONS

ZmPGH1 functions as a suppressor of programmed cell death induced by at least two autoactive NLR proteins and by two chemical inducers. These findings deepen our understanding of the control of the HR in plants.

Discovering and Constructing ceRNA-miRNA-Target Gene Regulatory Networks during Anther Development in Maize

The “competing endogenous RNA (ceRNA) hypothesis” has recently been proposed for a new type of gene regulatory model in many organisms. Anther development is a crucial biological process in plant reproduction, and its gene regulatory network (GRN) has been gradually revealed during the past two decades. However, it is still unknown whether ceRNAs contribute to anther development and sexual reproduction in plants. We performed RNA and small RNA sequencing of anther tissues sampled at three developmental stages in two maize lines. A total of 28,233 stably transcribed loci, 61 known and 51 potentially novel microRNAs (miRNAs) were identified from the transcriptomes. Predicted ceRNAs and target genes were found to conserve in sequences of recognition sites where their corresponding miRNAs bound. We then reconstructed 79 ceRNA-miRNA-target gene regulatory networks consisting of 51 known miRNAs, 28 potentially novel miRNAs, 619 ceRNA-miRNA pairs, and 869 miRNA-target gene pairs. More than half of the regulation pairs showed significant negative correlations at transcriptional levels. Several well-studied miRNA-target gene pairs associated with plant flower development were located in some networks, including miR156-SPL, miR159-MYB, miR160-ARF, miR164-NAC, miR172-AP2, and miR319-TCP pairs. Six target genes in the networks were found to be orthologs of functionally confirmed genes participating in anther development in plants. Our results provide an insight that the ceRNA-miRNA-target gene regulatory networks likely contribute to anther development in maize. Further functional studies on a number of ceRNAs, miRNAs, and target genes will facilitate our deep understanding on mechanisms of anther development and sexual plants reproduction.

Chromatin interaction maps reveal genetic regulation for quantitative traits in maize

Chromatin loops connect regulatory elements to their target genes. They serve as bridges between transcriptional regulation and phenotypic variation in mammals. However, spatial organization of regulatory elements and its impact on gene expression in plants remain unclear. Here, we characterize epigenetic features of active promoter proximal regions and candidate distal regulatory elements to construct high-resolution chromatin interaction maps for maize via long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET). The maps indicate that chromatin loops are formed between regulatory elements, and that gene pairs between promoter proximal regions tend to be co-expressed. The maps also demonstrated the topological basis of quantitative trait loci which influence gene expression and phenotype. Many promoter proximal regions are involved in chromatin loops with distal regulatory elements, which regulate important agronomic traits. Collectively, these maps provide a high-resolution view of 3D maize genome architecture, and its role in gene expression and phenotypic variation.

Spatial organization of regulatory elements and its impact on gene expression in plants remain unclear. Here, the authors construct maize chromatin interaction maps using chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and show their associations with gene expression and agronomic traits.

A maize cytochrome b-c1 complex subunit protein ZmQCR7 controls variation in the hypersensitive response

The gene GRMZM2G318346 which encodes a cytochrome b-c1 complex subunit 7 is associated with variation in strength of the hypersensitive response in maize. We previously identified a QTL at 3,545,354 bp (B73 reference genome V2) on maize chromosome 5 associated with variation in the hypersensitive response (HR) conferred by the autoactive R-gene Rp1-D21 (Olukolu et al. in PLoS Genet 10:e1004562 2014). In this study, we show that a gene at this locus, GRMZM2G318346 which encodes a cytochrome b-c1 complex subunit seven (ZmQCR7), an important part of the mitochondrial electron transport chain, can suppress HR mediated by Rp1-D21 in a transient expression assay. ZmQCR7 alleles from two maize lines, W22 and B73 differ for the encoded proteins at just two sites, amino acid 27 (threonine and alanine in B73 and W22, respectively) and amino acid 109 (asparagine and serine), however, the B73 allele is much more effective at suppressing HR. We show that variation at amino acid 27 controlled this variation in HR-suppressing effects. We furthermore demonstrate that the B73 allele of ZmQCR7 can suppress HR induced by RPM1(D505 V), another autoactive R-gene, and that Arabidopsis homologs of ZmQCR7 can also suppress NLR-induced HR. The implications of these findings are discussed.

Evolutionarily informed deep learning methods for predicting relative transcript abundance from DNA sequence

Deep learning methodologies have revolutionized prediction in many fields and show potential to do the same in molecular biology and genetics. However, applying these methods in their current forms ignores evolutionary dependencies within biological systems and can result in false positives and spurious conclusions. We developed two approaches that account for evolutionary relatedness in machine learning models: (i) gene-family-guided splitting and (ii) ortholog contrasts. The first approach accounts for evolution by constraining model training and testing sets to include different gene families. The second approach uses evolutionarily informed comparisons between orthologous genes to both control for and leverage evolutionary divergence during the training process. The two approaches were explored and validated within the context of mRNA expression level prediction and have the area under the ROC curve (auROC) values ranging from 0.75 to 0.94. Model weight inspections showed biologically interpretable patterns, resulting in the hypothesis that the 3' UTR is more important for fine-tuning mRNA abundance levels while the 5' UTR is more important for large-scale changes.

Dynamic Patterns of Gene Expression Additivity and Regulatory Variation throughout Maize Development

Gene expression variation is a key component underlying phenotypic variation and heterosis. Transcriptome profiling was performed on 23 different tissues or developmental stages of two maize inbreds, B73 and Mo17, as well as their F₁ hybrid. The obtained large-scale datasets provided opportunities to monitor the developmental dynamics of differential expression, additivity for gene expression, and regulatory variation. The transcriptome can be divided into ∼30 000 genes that are expressed in at least one tissue of one inbred and an additional ∼10 000 ″silent" genes that are not expressed in any tissue of any genotype, 90% of which are non-syntenic relative to other grasses. Many (∼74%) of the expressed genes exhibit differential expression in at least one tissue. However, the majority of genes with differential expression do not exhibit consistent differential expression in different tissues. These genes often exhibit tissue-specific differential expression with equivalent expression in other tissues, and in many cases they switch the directionality of differential expression in different tissues. This suggests widespread variation for tissue-specific regulation of gene expression between the two maize inbreds B73 and Mo17. Nearly 5000 genes are expressed in only one parent in at least one tissue (single parent expression) and 97% of these genes are expressed at mid-parent levels or higher in the hybrid, providing extensive opportunities for hybrid complementation in heterosis. In general, additive expression patterns are much more common than non-additive patterns, and this trend is more pronounced for genes with strong differential expression or single parent expression. There is relatively little evidence for non-additive expression patterns that are maintained in multiple tissues. The analysis of allele-specific expression allowed classification of cis- and trans-regulatory variation. Genes with cis-regulatory variation often exhibit additive expression and tend to have more consistent regulatory variation throughout development. In contrast, genes with trans-regulatory variation are enriched for non-additive patterns and often show tissue-specific differential expression. Taken together, this study provides a deeper understanding of regulatory variation and the degree of additive gene expression throughout maize development. The dynamic nature of differential expression, additivity, and regulatory variation imply abundant variability for tissue-specific regulatory mechanisms and suggest that connections between transcriptome and phenome will require expression data from multiple tissues.

Genome-Wide Association and Genomic Prediction Models of Tocochromanols in Fresh Sweet Corn Kernels

Sweet corn ( L.), a highly consumed fresh vegetable in the United States, varies for tocochromanol (tocopherol and tocotrienol) levels but makes only a limited contribution to daily intake of vitamin E and antioxidants. We performed a genome-wide association study of six tocochromanol compounds and 14 derivative traits across a sweet corn inbred line association panel to identify genes associated with natural variation for tocochromanols and vitamin E in fresh kernels. Concordant with prior studies in mature maize kernels, an association was detected between γ-tocopherol methyltransferase (vte4) and α-tocopherol content, along with () and () for tocotrienol variation. Additionally, two kernel starch synthesis genes, () and (), were associated with tocotrienols, with the strongest evidence for in combination with fixed, strong and alleles, accounting for the greater amount of tocotrienols in and lines. In prediction models with genome-wide markers, predictive abilities were higher for tocotrienols than tocopherols, and these models were superior to those that used marker sets targeting a priori genes involved in the biosynthesis and/or genetic control of tocochromanols. Through this quantitative genetic analysis, we have established a key step for increasing tocochromanols in fresh kernels of sweet corn for human health and nutrition.

RicetissueTFDB: A Genome-Wide Identification of Tissue-Specific Transcription Factors in Rice

Transcription factors (TFs) regulate plant gene expression in different tissues. To investigate TF genes in rice ( L.), a genome-wide TF identification was conducted with the japonica rice genome. This study identified 3078 putative TFs in 59 families. The TF number of the top 10 TF families accounted for 58% of the 3078 rice TFs. The three largest TF families were the myeloblastosis (MYB) superfamily, basic helix-loop-helix (bHLH), and far-red-impaired response (FAR1), which contained 413, 228, and 210 TF members, respectively. The expression profiles of the 3078 TF genes were surveyed with the RNA sequencing (RNA-seq) data of 13 rice tissue types. Based on these expression profiles, we validated 1087 TFs expressed in 13 rice tissue types, which accounted for 35.32% of the 3078 putative TFs. We further analyzed the tissue-specific TFs in rice. In total, 28, 14, 14, 10, 9, 5, 5, 4, 3, 3, 2, 11, and 1 tissue-specific TF sequences were identified in the dry seed, pistil, spikelet, aleurone, anther, ovules, embryo 25 d after pollination (DAP), seed 5 DAP, root, leaf, seed 10 DAP, shoot, and endosperm 25 DAP, respectively. Moreover, we constructed RicetissueTFDB (), a comprehensive and public rice TF database that integrates tissue expression characters, genomic location, and Gene Ontology (GO) terms for each TF. The RicetissueTFDB database will facilitate the identification of target TFs and the functional studies about rice TFs.

An updated gene atlas for maize reveals organ-specific and stress-induced genes

Maize (Zea mays L.), a model species for genetic studies, is one of the two most important crop species worldwide. The genome sequence of the reference genotype, B73, representative of the stiff stalk heterotic group was recently updated (AGPv4) using long-read sequencing and optical mapping technology. To facilitate the use of AGPv4 and to enable functional genomic studies and association of genotype with phenotype, we determined expression abundances for replicated mRNA-sequencing datasets from 79 tissues and five abiotic/biotic stress treatments revealing 36 207 expressed genes. Characterization of the B73 transcriptome across six organs revealed 4154 organ-specific and 7704 differentially expressed (DE) genes following stress treatment. Gene co-expression network analyses revealed 12 modules associated with distinct biological processes containing 13 590 genes providing a resource for further association of gene function based on co-expression patterns. Presence-absence variants (PAVs) previously identified using whole genome resequencing data from 61 additional inbred lines were enriched in organ-specific and stress-induced DE genes suggesting that PAVs may function in phenological variation and adaptation to environment. Relative to core genes conserved across the 62 profiled inbreds, PAVs have lower expression abundances which are correlated with their frequency of dispersion across inbreds and on average have significantly fewer co-expression network connections suggesting that a subset of PAVs may be on an evolutionary path to pseudogenization. To facilitate use by the community, we developed the Maize Genomics Resource website (maize.plantbiology.msu.edu) for viewing and data-mining these resources and deployed two new views on the maize electronic Fluorescent Pictograph Browser (bar.utoronto.ca/efp_maize).

Cell Type-Specific Transcriptomics of Lateral Root Formation and Plasticity

Lateral roots are a major determinant of root architecture and are instrumental for the efficient uptake of water and nutrients. Lateral roots consist of multiple cell types each expressing a unique transcriptome at a given developmental stage. Therefore, transcriptome analyses of complete lateral roots provide only average gene expression levels integrated over all cell types. Such analyses have the risk to mask genes, pathways and networks specifically expressed in a particular cell type during lateral root formation. Cell type-specific transcriptomics paves the way for a holistic understanding of the programming and re-programming of cells such as pericycle cells, involved in lateral root initiation. Recent discoveries have advanced the molecular understanding of the intrinsic genetic control of lateral root initiation and elongation. Moreover, the impact of nitrate availability on the transcriptional regulation of lateral root formation in Arabidopsis and cereals has been studied. In this review, we will focus on the systemic dissection of lateral root formation and its interaction with environmental nitrate through cell type-specific transcriptome analyses. These novel discoveries provide a better mechanistic understanding of postembryonic lateral root development in plants.

Identification and characterization of genes encoding the nuclear envelope LINC complex in the monocot species Zea mays

The linker of nucleoskeleton to cytoskeleton (LINC) complex is an essential multi-protein structure spanning the nuclear envelope. It connects the cytoplasm to the nucleoplasm, functions to maintain nuclear shape and architecture and regulates chromosome dynamics during cell division. Knowledge of LINC complex composition and function in the plant kingdom is primarily limited to Arabidopsis, but critically missing from the evolutionarily distant monocots, which include grasses, the most important agronomic crops worldwide. To fill this knowledge gap, we identified and characterized 22 maize genes, including a new grass-specific KASH gene family. By using bioinformatic, biochemical and cell biological approaches, we provide evidence that representative KASH candidates localize to the nuclear periphery and interact with Zea mays (Zm)SUN2 in vivo FRAP experiments using domain deletion constructs verified that this SUN-KASH interaction was dependent on the SUN but not the coiled-coil domain of ZmSUN2. A summary working model is proposed for the entire maize LINC complex encoded by conserved and divergent gene families. These findings expand our knowledge of the plant nuclear envelope in a model grass species, with implications for both basic and applied cellular research.This article has an associated First Person interview with the first author of the paper.

Genome-wide association analysis of stalk biomass and anatomical traits in maize

BACKGROUND

Maize stover is an important source of crop residues and a promising sustainable energy source in the United States. Stalk is the main component of stover, representing about half of stover dry weight. Characterization of genetic determinants of stalk traits provide a foundation to optimize maize stover as a biofuel feedstock. We investigated maize natural genetic variation in genome-wide association studies (GWAS) to detect candidate genes associated with traits related to stalk biomass (stalk diameter and plant height) and stalk anatomy (rind thickness, vascular bundle density and area).

RESULTS

Using a panel of 942 diverse inbred lines, 899,784 RNA-Seq derived single nucleotide polymorphism (SNP) markers were identified. Stalk traits were measured on 800 members of the panel in replicated field trials across years. GWAS revealed 16 candidate genes associated with four stalk traits. Most of the detected candidate genes were involved in fundamental cellular functions, such as regulation of gene expression and cell cycle progression. Two of the regulatory genes (Zmm22 and an ortholog of Fpa) that were associated with plant height were previously shown to be involved in regulating the vegetative to floral transition. The association of Zmm22 with plant height was confirmed using a transgenic approach. Transgenic lines with increased expression of Zmm22 showed a significant decrease in plant height as well as tassel branch number, indicating a pleiotropic effect of Zmm22.

CONCLUSION

Substantial heritable variation was observed in the association panel for stalk traits, indicating a large potential for improving useful stalk traits in breeding programs. Genome-wide association analyses detected several candidate genes associated with multiple traits, suggesting common regulatory elements underlie various stalk traits. Results of this study provide insights into the genetic control of maize stalk anatomy and biomass.

ZmRAP2.7, an AP2 Transcription Factor, Is Involved in Maize Brace Roots Development

In maize, shoot-borne roots dominate the whole root system and play essential roles in water and nutrient acquisition and lodging tolerance. Shoot-borne roots initiate at shoot nodes, including crown roots from the belowground nodes and brace roots from aboveground nodes. In contrast to crown roots, few genes for brace roots development have been identified. Here, we characterized a maize AP2/ERF transcription factor, ZmRAP2.7, to be involved in brace roots development. ZmRAP2.7 expressed in all types of roots, and the encoded protein localized in the nucleus with transcriptional activation activity. A maize transposon insert mutant RAP2.7-Mu defective in ZmRAP2.7 expression revealed a decreased number of brace roots but not crown roots. Maize Corngrass1 mutant, which showed an elevated expression of ZmRAP2.7, however, revealed an increased number of brace roots. The ZmRAP2.7-based association analysis in a maize panel further identified a SNP marker at the fifth exon of gene to be associated with number of brace roots. These results uncovered a function of ZmRAP2.7 in brace roots development and provided the valuable gene and allele for genetic improvement of maize root systems.

Genome-wide association study of maize plant architecture using F1 populations

Genome-wide association study of maize plant architecture using F₁ populations can better dissect various genetic effects that can provide precise guidance for genetic improvement in maize breeding. Maize grain yield has increased at least eightfold during the past decades. Plant architecture, including plant height, leaf angle, leaf length, and leaf width, has been changed significantly to adapt to higher planting density. Although the genetic architecture of these traits has been dissected using different populations, the genetic basis remains unclear in the F₁ population. In this work, we perform a genome-wide association study of the four traits using 573 F₁ hybrids with a mixed linear model approach and QTXNetwork mapping software. A total of 36 highly significant associated quantitative trait SNPs were identified for these traits, which explained 51.86-79.92% of the phenotypic variation and were contributed mainly by additive, dominance, and environment-specific effects. Heritability as a result of environmental interaction was more important for leaf angle and leaf length, while major effects (a, aa, and d) were more important for leaf width and plant height. The potential breeding values of the superior lines and superior hybrids were also predicted, and these values can be applied in maize breeding by direct selection of superior genotypes for the associated quantitative trait SNPs. A total of 108 candidate genes were identified for the four traits, and further analysis was performed to screen the potential genes involved in the development of maize plant architecture. Our results provide new insights into the genetic architecture of the four traits, and will be helpful in marker-assisted breeding for maize plant architecture.

Isolation and identification of a vegetative organ-specific promoter from maize

To avoid the unregulated overexpression of the exogenous genes, specific or inducible expression is necessary for some exogenous genes in transgenic plants. But little is known about organ- or tissue-specific promoters in maize. In the present study, the expression of a maize pentatricopeptide repeat (PPR) protein encoding gene, GRMZM2G129783, was analyzed by RNA-sequencing data and confirmed by quantitative real time PCR. The results showed that the PPR _{GRMZM2G129783} gene specifically expressed in vegetative organs. Consequently, a 1830 bp sequence upstream of the start codon of the promoter for GRMZM2G129783 gene was isolated from maize genome (P ₁₈₃₀ ). To validate whether the promoter possesses the vegetative organ-specificity, the full-length and three 5'-end deletion fragments of P ₁₈₃₀ of different length (1387, 437, and 146 bp) were fused with glucuronidase (GUS) gene to generate promoter::GUS constructs and transformed into tobacco. The transient expression and fluorometric GUS assay in transgenic tobacco showed that all promoter could drive the expression of the GUS gene, the - 437 to - 146 bp region possessed some crucial elements for root-specific expression, and the shortest and optimal sequence to maintain transcription activity was 146 and 437 bp in length, respectively. These results indicate that the promoter of the PPR _{GRMZM2G129783} gene is a vegetative organ-specific promoter and will be useful in transgenic modification of commercial crops for moderate specific expression after further evaluation in monocotyledons.

Proteases Underground: Analysis of the Maize Root Apoplast Identifies Organ Specific Papain-Like Cysteine Protease Activity

Plant proteases are key regulators of plant cell processes such as seed development, immune responses, senescence and programmed cell death (PCD). Apoplastic papain-like cysteine proteases (PL) are hubs in plant-microbe interactions and play an important role during abiotic stresses. The apoplast is a crucial interface for the interaction between plant and microbes. So far, apoplastic maize PL and their function have been mostly described for aerial parts. In this study, we focused on apoplastic PLCPs in the roots of maize plants. We have analyzed the phylogeny of maize PLCPs and investigated their protein abundance after salicylic acid (SA) treatment. Using activity-based protein profiling (ABPP) we have identified a novel root-specific PLCP belonging to the RD21-like subfamily, as well as three SA activated PLCPs. The root specific PLCP CP1C shares sequence and structural similarities to known CP1-like proteases. Biochemical analysis of recombinant CP1C revealed different substrate specificities and inhibitor affinities compared to the related proteases. This study characterized a root-specific PLCP and identifies differences between the SA-dependent activation of PLCPs in roots and leaves.

Genome-wide analysis of the basic Helix-Loop-Helix (bHLH) transcription factor family in maize

BACKGROUND

In plants, the basic helix-loop-helix (bHLH) transcription factors play key roles in diverse biological processes. Genome-wide comprehensive and systematic analyses of bHLH proteins have been well conducted in Arabidopsis, rice, tomato and other plant species. However, only few of bHLH family genes have been functional characterized in maize.

RESULTS

In this study, our genome-wide analysis identified 208 putative bHLH family proteins (ZmbHLH proteins) in maize (Zea mays). We classified these proteins into 18 subfamilies by comparing the ZmbHLHs with Arabidopsis thaliana bHLH proteins. Phylogenetic analysis, conserved protein motifs, and exon-intron patterns further supported the evolutionary relationships among these bHLH proteins. Genome distribution analysis found that the 208 ZmbHLH loci were located non-randomly on the ten maize chromosomes. Further, analysis of conserved cis-elements in the promoter regions, protein interaction networks, and expression patterns in roots, leaves, and seeds across developmental stages, suggested that bHLH family proteins in maize are probably involved in multiple physiological processes in plant growth and development.

CONCLUSION

We performed a genome-wide, systematic analysis of bHLH proteins in maize. This comprehensive analysis provides a useful resource that enables further investigation of the physiological roles and molecular functions of the ZmbHLH transcription factors.

The RNA-Seq-based high resolution gene expression atlas of chickpea (Cicer arietinum L.) reveals dynamic spatio-temporal changes associated with growth and development

Chickpea is one of the world's largest cultivated food legumes and is an excellent source of high-quality protein to the human diet. Plant growth and development are controlled by programmed expression of a suite of genes at the given time, stage, and tissue. Understanding how the underlying genome sequence translates into specific plant phenotypes at key developmental stages, information on gene expression patterns is crucial. Here, we present a comprehensive Cicer arietinum Gene Expression Atlas (CaGEA) across different plant developmental stages and organs covering the entire life cycle of chickpea. One of the widely used drought tolerant cultivars, ICC 4958 has been used to generate RNA-Seq data from 27 samples at 5 major developmental stages of the plant. A total of 816 million raw reads were generated and of these, 794 million filtered reads after quality control (QC) were subjected to downstream analysis. A total of 15,947 unique number of differentially expressed genes across different pairwise tissue combinations were identified. Significant differences in gene expression patterns contributing in the process of flowering, nodulation, and seed and root development were inferred in this study. Furthermore, differentially expressed candidate genes from "QTL-hotspot" region associated with drought stress response in chickpea were validated.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Comparative analysis of the male inflorescence transcriptome profiles of an ms22 mutant of maize

In modern agricultural production, maize is the most successful crop utilizing heterosis. 712C-ms22 is an important male sterile material in maize. In this study, we performed transcriptome sequencing analysis of the V10 stage of male inflorescence. Through this analysis, 27.63 million raw reads were obtained, and trimming of the raw data revealed 26.63 million clean reads, with an average match rate of 94.64%. Using Tophat software, we matched these clean reads to the maize reference genome. The abundance of 39,622 genes was measured, and 35,399 genes remained after filtering out the non-expressed genes across all the samples. These genes were classified into 19 categories by clusters of orthologous groups of protein annotation. Transcriptome sequencing analysis of the male sterile and fertile 712C-ms22 maize revealed some key DEGs that may be related to metabolic pathways. qRT-PCR analysis validated the gene expression patterns identified by RNA-seq. This analysis revealed some of the essential genes responsible for pollen development and for pollen tube elongation. Our findings provide useful markers of male sterility and new insights into the global mechanisms mediating male sterility in maize.

Expansion and diversification of the gibberellin receptor GIBBERELLIN INSENSITIVE DWARF1 (GID1) family in land plants

Here we uncover the major evolutionary events shaping the evolution of the GID1 family of gibberellin receptors in land plants at the sequence, structure and gene expression levels. Gibberellic acid (gibberellin, GA) controls key developmental processes in the life cycle of land plants. By interacting with the GIBBERELLIN INSENSITIVE DWARF1 (GID1) receptor, GA regulates the expression of a wide range of genes through different pathways. Here we report the systematic identification and classification of GID1s in 54 plants genomes, encompassing from bryophytes and lycophytes, to several monocots and eudicots. We investigated the evolutionary relationship of GID1s using a comparative genomics framework and found strong support for a previously proposed phylogenetic classification of this family in land plants. We identified lineage-specific expansions of particular subfamilies (i.e. GID1ac and GID1b) in different eudicot lineages (e.g. GID1b in legumes). Further, we found both, shared and divergent structural features between GID1ac and GID1b subgroups in eudicots that provide mechanistic insights on their functions. Gene expression data from several species show that at least one GID1 gene is expressed in every sampled tissue, with a strong bias of GID1b expression towards underground tissues and dry legume seeds (which typically have low GA levels). Taken together, our results indicate that GID1ac retained canonical GA signaling roles, whereas GID1b specialized in conditions of low GA concentrations. We propose that this functional specialization occurred initially at the gene expression level and was later fine-tuned by mutations that conferred greater GA affinity to GID1b, including a Phe residue in the GA-binding pocket. Finally, we discuss the importance of our findings to understand the diversification of GA perception mechanisms in land plants.

Comparative Genomics, Whole-Genome Re-sequencing and Expression Profile Analysis of Nucleobase:Cation Symporter 2 (NCS2) Genes in Maize

Nucleobase:cation symporter 2 (NCS2) proteins are important for the transport of free nucleobases, participating in diverse plant growth and developmental processes, as well as response to abiotic stress. To date, a comprehensive analysis of the NCS2 gene family has not been performed in maize. In this study, we conducted a comparative genomics analysis of NCS2 genes in 28 plant species, ranging from aquatic algae to land plants, concentrating mainly on maize. Gene duplication events contributed to the expansion of NCS2 genes from lower aquatic plants to higher angiosperms, and whole-genome/segmental and single-gene duplication events were responsible for the expansion of the maize NCS2 gene family. Phylogenetic construction showed three NCS2 subfamilies, I, II, and III. According to homology-based relationships, members of subfamily I are NCS2/AzgA-like genes, whereas those in subfamilies II and III are NCS2/NATs. Moreover, subfamily I exhibited ancient origins. A motif compositional analysis showed that one symbolic motif (motif 4) of the NCS2/NAT genes was absent in subfamily I. In maize, three NCS2/AzgA-like and 21 NCS2/NAT genes were identified, and purifying selection influenced the duplication of maize NCS2 genes. Additionally, a population genetic analysis of NCS2 genes revealed that ZmNCS2-21 showed the greatest diversity between the 78 inbred and 22 wild surveyed maize populations. An expression profile analysis using transcriptome data and quantitative real-time PCR revealed that NCS2 genes in maize are involved in diverse developmental processes and responses to abiotic stresses, including abscisic acid, salt (NaCl), polyethylene glycol, and low (4°C) and high (42°C) temperatures. ZmNCS2 genes with relatively close relationships had similar expression patterns, strongly indicating functional redundancy. Finally, ZmNCS2-16 and ZmNCS2-23 localize in the plasma membrane, which confirmed their predicted membrane structures. These results provide a foundation for future studies regarding the functions of ZmNCS2 proteins, particularly those with potentially important roles in plant responses to abiotic stresses.

Genome wide association study for gray leaf spot resistance in tropical maize core

Gray leaf spot is a maize foliar disease with worldwide distribution and can drastically reduce the production in susceptible genotypes. Published works indicate that resistance to gray leaf spot is a complex trait controlled by multiple genes, with additive effect and influenced by environment. The aim of this study was to identify genomic regions, including putative genes, associated with resistance to gray leaf spot under natural conditions of disease occurrence. A genome wide association study was conducted with 355,972 single nucleotide polymorphism markers on a phenotypic data composed by 157 tropical maize inbred lines, evaluated at Maringá -Brazil. Seven single nucleotide polymorphisms were significantly associated with gray leaf spot, some of which were localized to previously reported quantitative trait loci regions. Three gene models linked to the associated single nucleotide polymorphism were expressed at flowering time and tissue related with gray leaf spot infection, explaining a considerable proportion of the phenotypic variance, ranging from 0.34 to 0.38. The gene model GRMZM2G073465 (bin 10.07) encodes a cysteine protease3 protein, gene model GRMZM2G007188 (bin 1.02) expresses a rybosylation factor-like protein and the gene model GRMZM2G476902 (bin 4.08) encodes an armadillo repeat protein. These three proteins are related with plant defense pathway. Once these genes are validated in next studies, they will be useful for marker-assisted selection and can help improve the understanding of maize resistance to gray leaf spot.

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 ).

A calcium-dependent protein kinase, ZmCPK32, specifically expressed in maize pollen to regulate pollen tube growth

Calcium-dependent protein kinases (CPKs) play an essential role in the regulation of pollen tube growth. Although CPK genes have been identified in maize, and some have been functionally characterized, the molecular function of ZmCPKs associated with pollen tube development remains less well studied. Here, we report that a pollen-specific CPK, ZmCPK32, is involved in the regulation of pollen germination and tube extension. ZmCPK32 exhibited CPK activity and was localized on the plasma membrane and punctate internal membrane compartments via N-terminal acylation. In situ hybridization and real-time PCR revealed that ZmCPK32 transcripts accumulated in pollen and expression was dramatically upregulated during shedding. To elucidate the function of this gene, we transiently expressed a ZmCPK32-GFP fusion protein in tobacco pollen using microparticle bombardment. ZmCPK32 accumulation inhibited pollen germination and reduced pollen tube growth, but this effect was abolished when the kinase-inactive variant was expressed, indicating that kinase activity is critical for its regulatory function. In addition, the plasma membrane localization of ZmCPK32 is essential for regulating polar growth, as pollen expressing the cytosol-localized kinase displayed reduced tube length but germinated well. Moreover, the constitutively active form of ZmCPK32 enhanced the reduction in the germination rate, indicating that the specific activation of ZmCPK32 via calcium ions at the cortical growth point is essential for regulating appropriate germination. The results suggest that ZmCPK32 is functionally associated with pollen tube growth, and could represent a potential target for breeding male-sterile maize.

Drought-tolerant and drought-sensitive genotypes of maize (Zea mays L.) differ in contents of endogenous brassinosteroids and their drought-induced changes

The contents of endogenous brassinosteroids (BRs) together with various aspects of plant morphology, water management, photosynthesis and protection against cell damage were assessed in two maize genotypes that differed in their drought sensitivity. The presence of 28-norbrassinolide in rather high quantities (1-2 pg mg-1 fresh mass) in the leaves of monocot plants is reported for the first time. The intraspecific variability in the presence/content of the individual BRs in drought-stressed plants is also described for the first time. The drought-resistant genotype was characterised by a significantly higher content of total endogenous BRs (particularly typhasterol and 28-norbrassinolide) compared with the drought-sensitive genotype. On the other hand, the drought-sensitive genotype showed higher levels of 28-norcastasterone. Both genotypes also differed in the drought-induced reduction/elevation of the levels of 28-norbrassinolide, 28-norcastasterone, 28-homocastasterone and 28-homodolichosterone. The differences observed between both genotypes in the endogenous BR content are probably correlated with their different degrees of drought sensitivity, which was demonstrated at various levels of plant morphology, physiology and biochemistry.

A comparative transcriptional landscape of maize and sorghum obtained by single-molecule sequencing

Maize and sorghum are both important crops with similar overall plant architectures, but they have key differences, especially in regard to their inflorescences. To better understand these two organisms at the molecular level, we compared expression profiles of both protein-coding and noncoding transcripts in 11 matched tissues using single-molecule, long-read, deep RNA sequencing. This comparative analysis revealed large numbers of novel isoforms in both species. Evolutionarily young genes were likely to be generated in reproductive tissues and usually had fewer isoforms than old genes. We also observed similarities and differences in alternative splicing patterns and activities, both among tissues and between species. The maize subgenomes exhibited no bias in isoform generation; however, genes in the B genome were more highly expressed in pollen tissue, whereas genes in the A genome were more highly expressed in endosperm. We also identified a number of splicing events conserved between maize and sorghum. In addition, we generated comprehensive and high-resolution maps of poly(A) sites, revealing similarities and differences in mRNA cleavage between the two species. Overall, our results reveal considerable splicing and expression diversity between sorghum and maize, well beyond what was reported in previous studies, likely reflecting the differences in architecture between these two species.

Genome-wide association studies of doubled haploid exotic introgression lines for root system architecture traits in maize (Zea mays L.)

Root system architecture (RSA) is becoming recognized as important for water and nutrient acquisition in plants. This study focuses on finding single nucleotide polymorphisms (SNPs) associated with seedling RSA traits from 300 doubled haploid (DH) lines derived from crosses between Germplasm Enhancement of Maize (GEM) accessions and inbred lines PHB47 and PHZ51. These DH lines were genotyped using 62,077 SNP markers, while root and shoot phenotype data were collected from 14-day old seedlings. Genome-wide association studies (GWAS) were conducted using three models to offset false positives/negatives. Multiple SNPs associated with seedling root traits were detected, some of which were within or linked to gene models that showed expression in seedling roots. Significant trait associations involving the SNP S5_152926936 on Chromosome 5 were detected in all three models, particularly the trait network area. The SNP is within the gene model GRMZM2G021110, which is expressed in roots at seedling stage. SNPs that were significantly associated with seedling root traits, and closely linked to gene models that encode proteins associated with root development were also detected. This study shows that the GEM-DH panel may be a source of allelic diversity for genes controlling seedling root development.

A plastidial pantoate transporter with a potential role in pantothenate synthesis

The pantothenate (vitamin B5) synthesis pathway in plants is not fully defined because the subcellular site of its ketopantoate → pantoate reduction step is unclear. However, the pathway is known to be split between cytosol, mitochondria, and potentially plastids, and inferred to involve mitochondrial or plastidial transport of ketopantoate or pantoate. No proteins that mediate these transport steps have been identified. Comparative genomic and transcriptomic analyses identified Arabidopsis thaliana BASS1 (At1g78560) and its maize (Zea mays) ortholog as candidates for such a transport role. BASS1 proteins belong to the bile acid : sodium symporter family and share similarity with the Salmonella enterica PanS pantoate/ketopantoate transporter and with predicted bacterial transporters whose genes cluster on the chromosome with pantothenate synthesis genes. Furthermore, Arabidopsis BASS1 is co-expressed with genes related to metabolism of coenzyme A, the cofactor derived from pantothenate. Expression of Arabidopsis or maize BASS1 promoted the growth of a S. enterica panB panS mutant strain when pantoate, but not ketopantoate, was supplied, and increased the rate of [3H]pantoate uptake. Subcellular localization of green fluorescent protein fusions in Nicotiana tabacum BY-2 cells demonstrated that Arabidopsis BASS1 is targeted solely to the plastid inner envelope. Two independent Arabidopsis BASS1 knockout mutants accumulated pantoate ∼10-fold in leaves and had smaller seeds. Taken together, these data indicate that BASS1 is a physiologically significant plastidial pantoate transporter and that the pantoate reduction step in pantothenate biosynthesis could be at least partly localized in plastids.

Maize leaf PPDK regulatory protein isoform-2 is specific to bundle sheath chloroplasts and paradoxically lacks a Pi-dependent PPDK activation activity

In C4 plants, the pyruvate phosphate dikinase regulatory protein (PDRP) regulates the C4 pathway enzyme pyruvate phosphate dikinase (PPDK) in response to changes in incident light intensity. In maize (Zea mays) leaves, two distinct isoforms of PDRP are expressed, ZmPDRP1 and ZmPDRP2. The properties and C4 function of the ZmPDRP1 isoform are well understood. However, the PDRP2 isoform has only recently been identified and its properties and function(s) in maize leaves are unknown. We therefore initiated an investigation into the maize PDRP2 isoform by performing a side by side comparison of its enzyme properties and cell-specific distribution with PDRP1. In terms of enzyme functionality, PDRP2 was found to possess the same protein kinase-specific activity as PDRP1. However, the PDRP2 isoform was found to lack the phosphotransferase activity of the bifunctional PDRP1 isoform except when PDRP2 in the assays is elevated 5- to 10-fold. A primarily immuno-based approach was used to show that PDRP1 is strictly expressed in mesophyll cells and PDRP2 is strictly expressed in bundle sheath strand cells (BSCs). Additionally, using in situ immunolocalization, we establish a regulatory target for PDRP2 by showing a significant presence of C4 PPDK in BSC chloroplasts. However, a metabolic role for PPDK in this compartment is obscure, assuming PPDK accumulating in this compartment would be irreversibly inactivated each dark cycle by a monofunctional PDRP2.

Sorghum root-system classification in contrasting P environments reveals three main rooting types and root-architecture-related marker-trait associations

BACKGROUND AND AIMS

Roots facilitate acquisition of macro- and micronutrients, which are crucial for plant productivity and anchorage in the soil. Phosphorus (P) is rapidly immobilized in the soil and hardly available for plants. Adaptation to P scarcity relies on changes in root morphology towards rooting systems well suited for topsoil foraging. Root-system architecture (RSA) defines the spatial organization of the network comprising primary, lateral and stem-derived roots and is important for adaptation to stress conditions. RSA phenotyping is a challenging task and essential for understanding root development.

METHODS

In this study, 19 traits describing RSA were analysed in a diversity panel comprising 194 sorghum genotypes, fingerprinted with a 90-k single-nucleotide polymorphism (SNP) array and grown under low and high P availability.

KEY RESULTS

Multivariate analysis was conducted and revealed three different RSA types: (1) a small root system; (2) a compact and bushy rooting type; and (3) an exploratory root system, which might benefit plant growth and development if water, nitrogen (N) or P availability is limited. While several genotypes displayed similar rooting types in different environments, others responded to P scarcity positively by developing more exploratory root systems, or negatively with root growth suppression. Genome-wide association studies revealed significant quantitative trait loci (P < 2.9 × 10-6) on chromosomes SBI-02, SBI-03, SBI-05 and SBI-09. Co-localization of significant and suggestive (P < 5.7 × 10-5) associations for several traits indicated hotspots controlling root-system development on chromosomes SBI-02 and SBI-03.

CONCLUSIONS

Sorghum genotypes with a compact, bushy and shallow root system provide potential adaptation to P scarcity in the field by allowing thorough topsoil foraging, while genotypes with an exploratory root system may be advantageous if N or water is the limiting factor, although such genotypes showed highest P uptake levels under the artificial conditions of the present study.

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.

Control of Maize Vegetative and Reproductive Development, Fertility, and rRNAs Silencing by HISTONE DEACETYLASE 108

Histone deacetylases (HDACs) catalyze the removal of acetyl groups from acetylated histone tails that consequently interact more closely with DNA, leading to chromatin state refractory to transcription. Zea mays HDA108 belongs to the Rpd3/HDA1 HDAC family and is ubiquitously expressed during development. The newly isolated hda108/hda108 insertional mutant exhibited many developmental defects: significant reduction in plant height, alterations of shoot and leaf development, and alterations of inflorescence patterning and fertility. Western blot analyses and immunolocalization experiments revealed an evident increase in histone acetylation, accompanied by a marked reduction in H3K9 dimethylation, in mutant nuclei. The DNA methylation status, in the CHG sequence context, and the transcript level of ribosomal sequences were also affected in hda108 mutants, while enrichment in H3 and H4 acetylation characterizes both repetitive and nonrepetitive transcriptional up-regulated loci. RNA-Seq of both young leaf and anthers indicated that transcription factor expression is highly affected and that the pollen developmental program is disrupted in hda108 mutants. Crosses between hda108/hda108 and epiregulator mutants did not produce any double mutant progeny indicating possible genetic interactions of HDA108 with distinct epigenetic pathways. Our findings indicate that HDA108 is directly involved in regulation of maize development, fertility, and epigenetic regulation of genome activity.

The limited role of differential fractionation in genome content variation and function in maize (Zea mays L.) inbred lines

Maize is a diverse paleotetraploid species with considerable presence/absence variation and copy number variation. One mechanism through which presence/absence variation can arise is differential fractionation. Fractionation refers to the loss of duplicate gene pairs from one of the maize subgenomes during diploidization. Differential fractionation refers to non-shared gene loss events between individuals following a whole-genome duplication event. We investigated the prevalence of presence/absence variation resulting from differential fractionation in the syntenic portion of the genome using two whole-genome de novo assemblies of the inbred lines B73 and PH207. Between these two genomes, syntenic genes were highly conserved with less than 1% of syntenic genes being subject to differential fractionation. The few variably fractionated syntenic genes that were identified are unlikely to contribute to functional phenotypic variation, as there is a significant depletion of these genes in annotated gene sets. In further comparisons of 60 diverse inbred lines, non-syntenic genes were six times more likely to be variable than syntenic genes, suggesting that comparisons among additional genome assemblies are not likely to result in the discovery of large-scale presence/absence variation among syntenic genes.

Genetic Control of Root System Development in Maize

The maize root system comprises structurally and functionally different root types. Mutant analyses have revealed that root-type-specific genetic regulators intrinsically determine the maize root system architecture. Molecular cloning of these genes has demonstrated that key elements of auxin signal transduction, such as LOB domain (LBD) and Aux/IAA proteins, are instrumental for seminal, shoot-borne, and lateral root initiation. Moreover, genetic analyses have demonstrated that genes related to exocytotic vesicle docking, cell wall loosening, and cellulose synthesis and organization control root hair elongation. The identification of upstream regulators, protein interaction partners, and downstream targets of these genes together with cell-type-specific transcriptome analyses have provided novel insights into the regulatory networks controlling root development and architecture in 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.

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.

Characterization of the maize lipoxygenase gene family in relation to aflatoxin accumulation resistance

Maize (Zea mays L.) is a globally important staple food crop prone to contamination by aflatoxin, a carcinogenic secondary metabolite produced by the fungus Aspergillus flavus. An efficient approach to reduce accumulation of aflatoxin is the development of germplasm resistant to colonization and toxin production by A. flavus. Lipoxygenases (LOXs) are a group of non-heme iron containing dioxygenase enzymes that catalyze oxygenation of polyunsaturated fatty acids (PUFAs). LOX derived oxylipins play critical roles in plant defense against pathogens including A. flavus. The objectives of this study were to summarize sequence diversity and expression patterns for all LOX genes in the maize genome, and map their effect on aflatoxin accumulation via linkage and association mapping. In total, 13 LOX genes were identified, characterized, and mapped. The sequence of one gene, ZmLOX10, is reported from 5 inbred lines. Genes ZmLOX1/2, 5, 8, 9, 10 and 12 (GRMZM2G156861, or V4 numbers ZM00001D042541 and Zm00001D042540, GRMZM2G102760, GRMZM2G104843, GRMZM2G017616, GRMZM2G015419, and GRMZM2G106748, respectively) fell under previously published QTL in one or more mapping populations and are linked to a measurable reduction of aflatoxin in maize grains. Association mapping results found 28 of the 726 SNPs tested were associated with reduced aflatoxin levels at p ≤ 9.71 x 10-4 according to association statistics. These fell within or near nine of the ZmLOX genes. This work confirms the importance of some lipoxygenases for resistance to aflatoxin accumulation and may be used to direct future genetic selection in maize.

Genome-wide identification, phylogeny and expression analyses of SCARECROW-LIKE(SCL) genes in millet (Setaria italica)

As a member of the GRAS gene family, SCARECROW-LIKE (SCL) genes encode transcriptional regulators that are involved in plant information transmission and signal transduction. In this study, 44 SCL genes including two SCARECROW genes in millet were identified to be distributed on eight chromosomes, except chromosome 6. All the millet genes contain motifs 6-8, indicating that these motifs are conserved during the evolution. SCL genes of millet were divided into eight groups based on the phylogenetic relationship and classification of Arabidopsis SCL genes. Several putative millet orthologous genes in Arabidopsis, maize and rice were identified. High throughput RNA sequencing revealed that the expressions of millet SCL genes in root, stem, leaf, spica, and along leaf gradient varied greatly. Analyses combining the gene expression patterns, gene structures, motif compositions, promoter cis-elements identification, alternative splicing of transcripts and phylogenetic relationship of SCL genes indicate that the these genes may play diverse functions. Functionally characterized SCL genes in maize, rice and Arabidopsis would provide us some clues for future characterization of their homologues in millet. To the best of our knowledge, this is the first study of millet SCL genes at the genome wide level. Our work provides a useful platform for functional analysis of SCL genes in millet, a model crop for C₄ photosynthesis and bioenergy studies.

Integration of quantitated expression estimates from polyA-selected and rRNA-depleted RNA-seq libraries

BACKGROUND

The availability of fast alignment-free algorithms has greatly reduced the computational burden of RNA-seq processing, especially for relatively poorly assembled genomes. Using these approaches, previous RNA-seq datasets could potentially be processed and integrated with newly sequenced libraries. Confounding factors in such integration include sequencing depth and methods of RNA extraction and selection. Different selection methods (typically, either polyA-selection or rRNA-depletion) omit different RNAs, resulting in different fractions of the transcriptome being sequenced. In particular, rRNA-depleted libraries sample a broader fraction of the transcriptome than polyA-selected libraries. This study aimed to develop a systematic means of accounting for library type that allows data from these two methods to be compared.

RESULTS

The method was developed by comparing two RNA-seq datasets from ovine macrophages, identical except for RNA selection method. Gene-level expression estimates were obtained using a two-part process centred on the high-speed transcript quantification tool Kallisto. Firstly, a set of reference transcripts was defined that constitute a standardised RNA space, with expression from both datasets quantified against it. Secondly, a simple ratio-based correction was applied to the rRNA-depleted estimates. The outcome is an almost perfect correlation between gene expression estimates, independent of library type and across the full range of levels of expression.

CONCLUSION

A combination of reference transcriptome filtering and a ratio-based correction can create equivalent expression profiles from both polyA-selected and rRNA-depleted libraries. This approach will allow meta-analysis and integration of existing RNA-seq data into transcriptional atlas projects.

E+ subgroup PPR protein defective kernel 36 is required for multiple mitochondrial transcripts editing and seed development in maize and Arabidopsis

Mitochondria are semi-autonomous organelles that are the powerhouse of the cells. Plant mitochondrial RNA editing guided by pentatricopeptide repeat (PPR) proteins is essential for energy production. We identify a maize defective kernel mutant dek36, which produces small and collapsed kernels, leading to embryos and/or seedlings lethality. Seed filling in dek36 is drastically impaired, in line with the defects observed in the organization of endosperm transfer tissue. Positional cloning reveals that DEK36, encoding a mitochondria-targeted E+ subgroup PPR protein, is required for mitochondrial RNA editing at atp4-59, nad7-383 and ccmF_N -302, thus resulting in decreased activities of mitochondrial complex I, complex III and complex IV in dek36. Loss-of-function of its Arabidopsis ortholog At DEK36 causes arrested embryo and endosperm development, leading to embryo lethality. At_dek36 also has RNA editing defects in atp4, nad7, ccmF_N1 and ccmF_N2 , but at the nonconserved sites. Importantly, efficiency of all editing sites in ccmF_N1 , ccmF_N2 and rps12 is severely decreased in At_dek36, probably caused by the impairment of their RNA stabilization. These results suggest that the DEK36 orthologue pair are essential for embryo and endosperm development in both maize and Arabidopsis, but through divergent function in regulating RNA metabolism of their mitochondrial targets.

Structural Characterization of Maize SIRK1 Kinase Domain Reveals an Unusual Architecture of the Activation Segment

Kinases are primary regulators of plant metabolism and excellent targets for plant breeding. However, most kinases, including the abundant receptor-like kinases (RLK), have no assigned role. SIRK1 is a leucine-rich repeat receptor-like kinase (LRR-RLK), the largest family of RLK. In Arabidopsis thaliana, SIRK1 (AtSIRK1) is phosphorylated after sucrose is resupplied to sucrose-starved seedlings and it modulates the sugar response by phosphorylating several substrates. In maize, the ZmSIRK1 expression is altered in response to drought stress. In neither Arabidopsis nor in maize has the function of SIRK1 been completely elucidated. As a first step toward the biochemical characterization of ZmSIRK1, we obtained its recombinant kinase domain, demonstrated that it binds AMP-PNP, a non-hydrolysable ATP-analog, and solved the structure of ZmSIRK1- AMP-PNP co-crystal. The ZmSIRK1 crystal structure revealed a unique conformation for the activation segment. In an attempt to find inhibitors for ZmSIRK1, we screened a focused small molecule library and identified six compounds that stabilized ZmSIRK1 against thermal melt. ITC analysis confirmed that three of these compounds bound to ZmSIRK1 with low micromolar affinity. Solving the 3D structure of ZmSIRK1-AMP-PNP co-crystal provided information on the molecular mechanism of ZmSIRK1 activity. Furthermore, the identification of small molecules that bind this kinase can serve as initial backbone for development of new potent and selective ZmSIRK1 antagonists.

Cyclopiazonic Acid Is a Pathogenicity Factor for Aspergillus flavus and a Promising Target for Screening Germplasm for Ear Rot Resistance

Aspergillus flavus, an opportunistic pathogen, contaminates maize and other key crops with carcinogenic aflatoxins (AFs). Besides AFs, A. flavus makes many more secondary metabolites (SMs) whose toxicity in insects or vertebrates has been studied. However, the role of SMs in the invasion of plant hosts by A. flavus remains to be investigated. Cyclopiazonic acid (CPA), a neurotoxic SM made by A. flavus, is a nanomolar inhibitor of endoplasmic reticulum calcium ATPases (ECAs) and a potent inducer of cell death in plants. We hypothesized that CPA, by virtue of its cytotoxicity, may serve as a key pathogenicity factor that kills plant cells and supports the saprophytic life style of the fungus while compromising the host defense response. This proposal was tested by two complementary approaches. A comparison of CPA levels among A. flavus isolates indicated that CPA may be a determinant of niche adaptation, i.e., isolates that colonize maize make more CPA than those restricted only to the soil. Further, mutants in the CPA biosynthetic pathway are less virulent in causing ear rot than their wild-type parent in field inoculation assays. Additionally, genes encoding ECAs are expressed in developing maize seeds and are induced by A. flavus infection. Building on these results, we developed a seedling assay in which maize roots were exposed to CPA, and cell death was measured as Evans Blue uptake. Among >40 maize inbreds screened for CPA tolerance, inbreds with proven susceptibility to ear rot were also highly CPA sensitive. The publicly available data on resistance to silk colonization or AF contamination for many of the lines was also broadly correlated with their CPA sensitivity. In summary, our studies show that i) CPA serves as a key pathogenicity factor that enables the saprophytic life style of A. flavus and ii) maize inbreds are diverse in their tolerance to CPA. Taking advantage of this natural variation, we are currently pursuing both genome-wide and candidate gene approaches to identify novel components of maize resistance to Aspergillus ear rot.

Sparse panicle1 is required for inflorescence development in Setaria viridis and maize

Setaria viridis is a rapid-life-cycle model panicoid grass. To identify genes that may contribute to inflorescence architecture and thus have the potential to influence grain yield in related crops such as maize, we conducted an N-nitroso-N-methylurea (NMU) mutagenesis of S. viridis and screened for visible inflorescence mutant phenotypes. Of the approximately 2,700 M₂ families screened, we identified four recessive sparse panicle mutants (spp1-spp4) characterized by reduced and uneven branching of the inflorescence. To identify the gene underlying the sparse panicle1 (spp1) phenotype, we performed bulked segregant analysis and deep sequencing to fine map it to an approximately 1 Mb interval. Within this interval, we identified disruptive mutations in two genes. Complementation tests between spp1 and spp3 revealed they were allelic, and deep sequencing of spp3 identified an independent disruptive mutation in SvAUX1 (AUXIN1), one of the two genes in the ∼1 Mb interval and the only gene disruption shared between spp1 and spp3. SvAUX1 was found to affect both inflorescence development and root gravitropism in S. viridis. A search for orthologous mutant alleles in maize confirmed a very similar role of ZmAUX1 in maize, which highlights the utility of S. viridis in accelerating functional genomic studies in maize.

MATRILINEAL, a sperm-specific phospholipase, triggers maize haploid induction

Sexual reproduction in flowering plants involves double fertilization, the union of two sperm from pollen with two sex cells in the female embryo sac. Modern plant breeders increasingly seek to circumvent this process to produce doubled haploid individuals, which derive from the chromosome-doubled cells of the haploid gametophyte. Doubled haploid production fixes recombinant haploid genomes in inbred lines, shaving years off the breeding process. Costly, genotype-dependent tissue culture methods are used in many crops, while seed-based in vivo doubled haploid systems are rare in nature and difficult to manage in breeding programmes. The multi-billion-dollar maize hybrid seed business, however, is supported by industrial doubled haploid pipelines using intraspecific crosses to in vivo haploid inducer males derived from Stock 6, first reported in 1959 (ref. 5), followed by colchicine treatment. Despite decades of use, the mode of action remains controversial. Here we establish, through fine mapping, genome sequencing, genetic complementation, and gene editing, that haploid induction in maize (Zea mays) is triggered by a frame-shift mutation in MATRILINEAL (MTL), a pollen-specific phospholipase, and that novel edits in MTL lead to a 6.7% haploid induction rate (the percentage of haploid progeny versus total progeny). Wild-type MTL protein localizes exclusively to sperm cytoplasm, and pollen RNA-sequence profiling identifies a suite of pollen-specific genes overexpressed during haploid induction, some of which may mediate the formation of haploid seed. These findings highlight the importance of male gamete cytoplasmic components to reproductive success and male genome transmittance. Given the conservation of MTL in the cereals, this discovery may enable development of in vivo haploid induction systems to accelerate breeding in crop plants.

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.