Genome-Wide Identification and Analysis of Expression Profiles of Maize Mitogen-Activated Protein Kinase Kinase Kinase

Mitogen-activated protein kinase (MAPK) cascades are highly conserved signal transduction model in animals, yeast and plants. Plant MAPK cascades have been implicated in development and stress responses. Although MAPKKKs have been investigated in several plant species including Arabidopsis and rice, no systematic analysis has been conducted in maize. In this study, we performed a bioinformatics analysis of the entire maize genome and identified 74 MAPKKK genes. Phylogenetic analyses of MAPKKKs from maize, rice and Arabidopsis have classified them into three subgroups, which included Raf, ZIK and MEKK. Evolutionary relationships within subfamilies were also supported by exon-intron organizations and the conserved protein motifs. Further expression analysis of the MAPKKKs in microarray databases revealed that MAPKKKs were involved in important signaling pathways in maize different organs and developmental stages. Our genomics analysis of maize MAPKKK genes provides important information for evolutionary and functional characterization of this family in maize.     

Oxysterol binding proteins (OSBPs) and their encoding genes in Arabidopsis and rice

Cell wall deposition, biosynthesis of steroid hormones, and maintenance of membrane composition and integrity, are some of the crucial functions of sterols in plants. Followed by their synthesis in the endoplasmic reticulum, the sterols accumulate in the plasma membrane. The concept of sterol trafficking in plant cell is not well understood. The oxysterol binding proteins are implicated in sterol transport in non-plant systems. In the study, the oxysterol binding proteins in Arabidopsis and rice are described and classified. The Arabidopsis genome encodes 12 oxysterol binding proteins-related proteins (ORPs) as compared to 6 oxysterol binding proteins (OSBPs/ORPs) in rice. The protein alignment studies reveal that amino acid sequences for oxysterol binding proteins are relatively well conserved in Arabidopsis and rice. The rice OSBPs are classified based on their phylogenetic relationship with Arabidopsis ORPs. The sequence LOGO built on LOC_Os03g16690 indicated presence of fingerprint region of amino acids “EQVSHHPP” for Arabidopsis and rice OSBPs/ORPs. The organization of pleckstrin homology domain is identified in several OSBPs/ORPs in Arabidopsis and rice. The Arabidopsis oligonucleotide array data is explored to understand the expression patterns of ORPs under 17 different experimental conditions. The analysis showed the expression of ORPs in Arabidopsis is necessarily under the control of biotic stress, chemical, elicitor, hormone, light intensity, abiotic stress, and temperature conditions. The linear mean signal values for Arabidopsis ORPs revealed their relative expression patterns in different developmental stages. The genes for ORP3C and ORP3B are highly expressed in all developmental stages that were analyzed. The present study thus indicates crucial functional role of the individual members of this gene family in different environmental stress conditions.     

Phosphorylation of Arabidopsis thaliana MEKK1 via Ca2+ signaling as a part of the cold stress response

The Arabidopsis mitogen activated protein kinase kinase kinase (MEKK1) plays an important role in stress signaling. However, little is known about the upstream pathways of MEKK1. This report describes the regulation of MEKK1 activity during cold signaling. Immunoprecipitated MEKK1 from cold-treated Arabidopsis seedlings showed elevated kinase activity towards mitogen activated protein kinase kinase2 (MKK2), one of the candidate MEKK1 substrates. To clarify how MEKK1 becomes active in response to cold stress signaling, MEKK1 phosphorylation was monitored by an enzyme extracted from the seedlings grown under cold stress with or without EGTA. MEKK1 was phosphorylated after cold stress, but EGTA inhibited the phosphorylation. MKK2 was also phosphorylated by the same extract, but only when EGTA was absent. These results suggested that Ca²⁺ signaling occurred upstream of the MEKK1–MKK2 pathway. Full-length MEKK1 showed almost no activity but MEKK1 without the N-terminal region (MEKK1 KD) that retained the kinase domain had a strong ability to phosphorylate MKK2, demonstrating the inhibitory role of the N-terminal region of MEKK1. In addition, MEKK1 was phosphorylated by calcium/calmodulin-regulated receptor-like kinase (CRLK1), which suggested that CRLK1 is one of candidates located upstream of MEKK1.     

Comparison of Peptide Array Substrate Phosphorylation of c-Raf and Mitogen Activated Protein Kinase Kinase Kinase 8

Kinases are pivotal regulators of cellular physiology. The human genome contains more than 500 putative kinases, which exert their action via the phosphorylation of specific substrates. The determinants of this specificity are still only partly understood and as a consequence it is difficult to predict kinase substrate preferences from the primary structure, hampering the understanding of kinase function in physiology and prompting the development of technologies that allow easy assessment of kinase substrate consensus sequences. Hence, we decided to explore the usefulness of phosphorylation of peptide arrays comprising of 1176 different peptide substrates with recombinant kinases for determining kinase substrate preferences, based on the contribution of individual amino acids to total array phosphorylation. Employing this technology, we were able to determine the consensus peptide sequences for substrates of both c-Raf and Mitogen Activated Protein Kinase Kinase Kinase 8, two highly homologous kinases with distinct signalling roles in cellular physiology. The results show that although consensus sequences for these two kinases identified through our analysis share important chemical similarities, there is still some sequence specificity that could explain the different biological action of the two enzymes. Thus peptide arrays are a useful instrument for deducing substrate consensus sequences and highly homologous kinases can differ in their requirement for phosphorylation events.     

The rice Raf-like MAPKKK OsILA1 confers broad-spectrum resistance to bacterial blight by suppressing the OsMAPKK4–OsMAPK6 cascade

Mitogen-activated protein kinase kinase kinase (MAPKKK) are the first components of MAPK cascades, which play pivotal roles in signaling during plant development and physiological processes. The genome of rice encodes 75 MAPKKKs, of which 43 are Raf-like MAPKKKs. The functions and action modes of most of the Raf-like MAPKKKs, whether they function as bona fide MAPKKKs and which are their downstream MAPKKs, are largely unknown. Here, we identified the osmapkkk43 mutant, which conferred broad-spectrum resistance to Xanthomonas oryzae pv. oryzae (Xoo), the destructive bacterial pathogen of rice. Oryza sativa (Os)MAPKKK43 encoding a Raf-like MAPKKK was previously known as Increased Leaf Angle 1 (OsILA1). Genetic analysis indicated that OsILA1 functioned as a negative regulator and acted upstream of the OsMAPKK4–OsMAPK6 cascade in rice–Xoo interactions. Unlike classical MAPKKKs, OsILA1 mainly phosphorylated the threonine 34 site at the N-terminal domain of OsMAPKK4, which possibly influenced the stability of OsMAPKK4. The N-terminal domain of OsILA1 is required for its homodimer formation and its full phosphorylation capacity. Taken together, our findings reveal that OsILA1 acts as a negative regulator of the OsMAPKK4–OsMAPK6 cascade and is involved in rice–Xoo interactions.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Genome-Wide Analysis of Two-Component Systems and Prediction of Stress-Responsive Two-Component System Members in Soybean

In plants, the two-component systems (TCSs) play important roles in regulating diverse biological processes, including responses to environmental stress stimuli. Within the soybean genome, the TCSs consist of at least 21 histidine kinases, 13 authentic and pseudo-phosphotransfers and 18 type-A, 15 type-B, 3 type-C and 11 pseudo-response regulator proteins. Structural and phylogenetic analyses of soybean TCS members with their Arabidopsis and rice counterparts revealed similar architecture of their TCSs. We identified a large number of closely homologous soybean TCS genes, which likely resulted from genome duplication. Additionally, we analysed tissue-specific expression profiles of those TCS genes, whose data are available from public resources. To predict the putative regulatory functions of soybean TCS members, with special emphasis on stress-responsive functions, we performed comparative analyses from all the TCS members of soybean, Arabidopsis and rice and coupled these data with annotations of known abiotic stress-responsive cis-elements in the promoter region of each soybean TCS gene. Our study provides insights into the architecture and a solid foundation for further functional characterization of soybean TCS elements. In addition, we provide a new resource for studying the conservation and divergence among the TCSs within plant species and/or between plants and other organisms.     

A root-specific wall-associated kinase gene, HvWAK1, regulates root growth and is highly divergent in barley and other cereals

Wall-associated receptor-like kinases (WAKs) are important candidates for directly linking the extracellular matrix with intracellular compartments and are involved in developmental processes and stress response. WAK gene family has been identified in plants such as Arabidopsis and rice. Here, we present a detailed analysis of the WAK1 gene from barley cv. Golden Promise, mapped to chromosome 5H. Three BAC clones corresponding to the WAK fragment were sequenced and the full-length WAK1 gene was characterized. The gene has three exons and two short introns with a coding region of 2,178 bp encoding a protein of 725 amino acids. A regulatory region was analyzed in ?1,000 bp sequence upstream to start codon. Using conserved domains database and SMART, various conserved domains such as GUB WAK Bind, epidermal growth factor CA, and protein kinase C as well as other regions like signal peptides, active sites, and transmembrane domains were identified. The gene organization of HvWAK1 was compared with wheat (TaWAK1) and Arabidopsis (AtWAK1), suggesting that the WAK1 gene organization has remained highly conserved. Nonetheless, WAK1 was found to be highly divergent when compared with sequences available from barley cv. Haruna Nijo (50 %), rice (46 %), wheat (21 %), Arabidopsis (25 %), and maize (19 %). This divergence may have facilitated a better adaptation to surrounding environments due to its role in communication between the extracellular matrix, cell, and outer environment. Semiquantitative RT-PCR-based expression analysis indicates HvWAK1 expression is specific to roots. Significant differences in root growth between GP wild type and GP-Ds mutant seedlings were observed under control and salt stress conditions.     

Complexity of rice Hsp100 gene family: lessons from rice genome sequence data

Elucidation of genome sequence provides an excellent platform to understand detailed complexity of the various gene families. Hsp100 is an important family of chaperones in diverse living systems. There are eight putative gene loci encoding for Hsp100 proteins in Arabidopsis genome. In rice, two full-length Hsp100 cDNAs have been isolated and sequenced so far. Analysis of rice genomic sequence by in silico approach showed that two isolated rice Hsp100 cDNAs correspond to Os05g44340 and Os02g32520 genes in the rice genome database. There appears to be three additional proteins (encoded by Os03g31300, Os04g32560 and Os04g33210 gene loci) that are variably homologous to Os05g44340 and Os02g32520 throughout the entire amino acid sequence. The above five rice Hsp100 genes show significant similarities in the signature sequences known to be conserved among Hsp100 proteins. While Os05g44340 encodes cytoplasmic Hsp100 protein, those encoded by the other four genes are predicted to have chloroplast transit peptides.     

Small RNA in rice genome

Rice has many characteristics of a model plant. The recent completion of the draft of the rice genome represents an important advance in our knowledge of plant biology and also has an important contribution to the understanding of general genomic evolution. Besides the rice genome finishing map, the next urgent step for rice researchers is to annotate the genes and noncoding functional sequences. The recent work shows that noncoding RNAs (ncRNAs) play significant roles in biological systems. We have explored all the known small RNAs (a kind of ncRNA) within rice genome and other six species sequences, including Arabidopsis, maize, yeast, worm, mouse and pig. As a result we find 160 out of 552 small RNAs (sRNAs) in database have homologs in 108 rice scaffolds, and almost all of them (99.41%) locate in intron regions of rice by gene predication. 19 sRNAs only appear in rice. More importantly, we find two special LJ14 sRNAs: one is located in a set of sRNA ZMU14SNR9(s) which only appears in three plants, 86% sequences of them can be compared as the same sequence in rice, Arabidopsis and maize; the other conserved sRNA XLHS7CU14 has a segment which appears in almost all these species from plants to animals. All these results indicate that sRNA do not have evident borderline between plants and animals.     

Genome-Wide Analysis of Mitogen-Activated Protein Kinase Gene Family in Maize

Mitogen-activated protein kinase (MAPK) cascades are universal signal transduction modules in plants. As the last component of the MAPK cascade (MAPKKK–MAPKK–MAPK), MAPK plays important roles in linking upstream kinases and downstream substrates. The MAPK proteins belong to a complex gene family in plants, with 20 MAPK genes in the Arabidopsis genome, 17 in the rice genome, and 21 in the poplar genome. Although the maize genome sequencing has been completed, no comprehensive study has been reported thus far for the MAPK gene family in maize. In this study, we identified 19 MAPK genes in maize. These ZmMPK genes belong to four groups (A–D) found in other plants. The phylogeny, chromosomal location, gene structure, and the functional relevancy of ZmMPK genes were analyzed. Moreover, we discuss the evolutionary divergence of MAPK genes in maize. Furthermore, we analyzed the expression profiles of ZmMPKs using the public microarray data and performed expression analyses in maize seedlings and adult plants. The data obtained from our study contribute to a better understanding of the complexity of MAPKs in plants and provide a useful reference for further functional analysis of MAPK genes in maize.     

A segment of the apospory-specific genomic region is highly microsyntenic not only between the apomicts Pennisetum squamulatum and buffelgrass, but also with a rice chromosome 11 centromeric-proximal genomic region

Bacterial artificial chromosome (BAC) clones from apomicts Pennisetum squamulatum and buffelgrass (Cenchrus ciliaris), isolated with the apospory-specific genomic region (ASGR) marker ugt197, were assembled into contigs that were extended by chromosome walking. Gene-like sequences from contigs were identified by shotgun sequencing and BLAST searches, and used to isolate orthologous rice contigs. Additional gene-like sequences in the apomicts' contigs were identified by bioinformatics using fully sequenced BACs from orthologous rice contigs as templates, as well as by interspecies, whole-contig cross-hybridizations. Hierarchical contig orthology was rapidly assessed by constructing detailed long-range contig molecular maps showing the distribution of gene-like sequences and markers, and searching for microsyntenic patterns of sequence identity and spatial distribution within and across species contigs. We found microsynteny between P. squamulatum and buffelgrass contigs. Importantly, this approach also enabled us to isolate from within the rice (Oryza sativa) genome contig Rice A, which shows the highest microsynteny and is most orthologous to the ugt197-containing C1C buffelgrass contig. Contig Rice A belongs to the rice genome database contig 77 (according to the current September 12, 2003, rice fingerprint contig build) that maps proximal to the chromosome 11 centromere, a feature that interestingly correlates with the mapping of ASGR-linked BACs proximal to the centromere or centromere-like sequences. Thus, relatedness between these two orthologous contigs is supported both by their molecular microstructure and by their centromeric-proximal location. Our discoveries promote the use of a microsynteny-based positional-cloning approach using the rice genome as a template to aid in constructing the ASGR toward the isolation of genes underlying apospory.     

MEKK2 Kinase Association with 14-3-3 Protein Regulates Activation of c-Jun N-terminal Kinase

MEKK2 (MAP/ERK kinase kinase-2) is a serine/threonine kinase that belongs to the MEKK/STE11 family of MAP kinase kinase kinases (MAP(3)Ks). MEKK2 integrates stress and mitogenic signals to the activation of NF-κB, JNK1/2, p38, and ERK5 pathways. We have found that MEKK2 is regulated through a phosphorylation-dependent association with 14-3-3, a group of adapters that modulate dimerization and association between proteins. We found that MEKK2 was phosphorylated at Thr-283, which resulted in decreased activation loop phosphorylation at Ser-519 and consequently reduced activity. Mechanistically, we found that MEKK2 associated with inactive MEKK2 in the absence of 14-3-3 binding, which led to trans-autophosphorylation of Ser-519. Enforced binding with 14-3-3 reduced Ser-519 trans-autophosphorylation. Expression of T283A MEKK2 within a MEKK2^−/− background enhanced stress-activated c-Jun N-terminal kinase activity while elevating IL-6 expression, but also reduced ERK activation with a corresponding reduced proliferation rate. These results indicate that Thr-283 phosphorylation is an important regulatory mechanism for MEKK2 activation.     

Genome-Wide Analysis of Repetitive Elements in Papaya

Papaya (Carica papaya L.) is an important fruit crop cultivated in tropical and subtropical regions worldwide. A first draft of its genome sequence has been recently released. Together with Arabidopsis, rice, poplar, grapevine and other genomes in the pipeline, it represents a good opportunity to gain insight into the organization of plant genomes. Here we report a detailed analysis of repetitive elements in the papaya genome, including transposable elements (TEs), tandemly-arrayed sequences, and high copy number genes. These repetitive sequences account for ～56% of the papaya genome with TEs being the most abundant at 52%, tandem repeats at 1.3% and high copy number genes at 3%. Most common types of TEs are represented in the papaya genome with retrotransposons being the dominant class, accounting for 40% of the genome. The most prevalent retrotransposons are Ty3-gypsy (27.8%) and Ty1-copia (5.5%). Among the tandem repeats, microsatellites are the most abundant in number, but represent only 0.19% of the genome. Minisatellites and satellites are less abundant, but represent 0.68% and 0.43% of the genome, respectively, due to greater repeat length. Despite an overall smaller gene repertoire in papaya than many other angiosperms, a significant fraction of genes (>2%) are present in large gene families with copy number greater than 20. This repeat database clarified a major part of the papaya genome organization and partly explained the lower gene repertoire in papaya than in Arabidopsis.     

From genomics to functional markers in the era of next-generation sequencing

The availability of complete genome sequences, along with other genomic resources for Arabidopsis, rice, pigeon pea, soybean and other crops, has revolutionized our understanding of the genetic make-up of plants. Next-generation DNA sequencing (NGS) has facilitated single nucleotide polymorphism discovery in plants. Functionally-characterized sequences can be identified and functional markers (FMs) for important traits can be developed at an ever-increasing ease. FMs are derived from sequence polymorphisms found in allelic variants of a functional gene. Linkage disequilibrium-based association mapping and homologous recombinants have been developed for identification of “perfect” markers for their use in crop improvement practices. Compared with many other molecular markers, FMs derived from the functionally characterized sequence genes using NGS techniques and their use provide opportunities to develop high-yielding plant genotypes resistant to various stresses at a fast pace.