Akhilesh K. Tyagi

F-Box Proteins in Rice. Genome-Wide Analysis, Classification, Temporal and Spatial Gene Expression during Panicle and Seed Development, and Regulation by Light and Abiotic Stress

F-box proteins constitute a large family in eukaryotes and are characterized by a conserved F-box motif (approximately 40 amino acids). As components of the Skp1p-cullin-F-box complex, F-box proteins are critical for the controlled degradation of cellular proteins. We have identified 687 potential F-box proteins in rice (Oryza sativa), the model monocotyledonous plant, by a reiterative database search. Computational analysis revealed the presence of several other functional domains, including leucine-rich repeats, kelch repeats, F-box associated domain, domain of unknown function, and tubby domain in F-box proteins. Based upon their domain composition, they have been classified into 10 subfamilies. Several putative novel conserved motifs have been identified in F-box proteins, which do not contain any other known functional domain. An analysis of a complete set of F-box proteins in rice is presented, including classification, chromosomal location, conserved motifs, and phylogenetic relationship. It appears that the expansion of F-box family in rice, in large part, might have occurred due to localized gene duplications. Furthermore, comprehensive digital expression analysis of F-box protein-encoding genes has been complemented with microarray analysis. The results reveal specific and/or overlapping expression of rice F-box protein-encoding genes during floral transition as well as panicle and seed development. At least 43 F-box protein-encoding genes have been found to be differentially expressed in rice seedlings subjected to different abiotic stress conditions. The expression of several F-box protein-encoding genes is also influenced by light. The structure and function of F-box proteins in plants is discussed in light of these results and the published information. These data will be useful for prioritization of F-box proteins for functional validation in rice.

De novo assembly of chickpea transcriptome using short reads for gene discovery and marker identification.

Chickpea ranks third among the food legume crops production in the world. However, the genomic resources available for chickpea are still very limited. In the present study, the transcriptome of chickpea was sequenced with short reads on Illumina Genome Analyzer platform. We have assessed the effect of sequence quality, various assembly parameters and assembly programs on the final assembly output. We assembled ∼107million high-quality trimmed reads using Velvet followed by Oases with optimal parameters into a non-redundant set of 53 409 transcripts (≥100 bp), representing about 28 Mb of unique transcriptome sequence. The average length of transcripts was 523 bp and N50 length of 900 bp with coverage of 25.7 rpkm (reads per kilobase per million). At the protein level, a total of 45 636 (85.5%) chickpea transcripts showed significant similarity with unigenes/predicted proteins from other legumes or sequenced plant genomes. Functional categorization revealed the conservation of genes involved in various biological processes in chickpea. In addition, we identified simple sequence repeat motifs in transcripts. The chickpea transcripts set generated here provides a resource for gene discovery and development of functional molecular markers. In addition, the strategy for de novo assembly of transcriptome data presented here will be helpful in other similar transcriptome studies.

MADS-box gene family in rice: genome-wide identification, organization and expression profiling during reproductive development and stress

BackgroundMADS-box transcription factors, besides being involved in floral organ specification, have also been implicated in several aspects of plant growth and development. In recent years, there have been reports on genomic localization, protein motif structure, phylogenetic relationships, gene structure and expression of the entire MADS-box family in the model plant system, Arabidopsis. Though there have been some studies in rice as well, an analysis of the complete MADS-box family along with a comprehensive expression profiling was still awaited after the completion of rice genome sequencing. Furthermore, owing to the role of MADS-box family in flower development, an analysis involving structure, expression and functional aspects of MADS-box genes in rice and Arabidopsis was required to understand the role of this gene family in reproductive development.ResultsA genome-wide molecular characterization and microarray-based expression profiling of the genes encoding MADS-box transcription factor family in rice is presented. Using a thorough annotation exercise, 75 MADS-box genes have been identified in rice and categorized into MIKCc, MIKC*, Mα, Mβ and Mγ groups based on phylogeny. Chromosomal localization of these genes reveals that 16 MADS-box genes, mostly MIKCc-type, are located within the duplicated segments of the rice genome, whereas most of the M-type genes, 20 in all, seem to have resulted from tandem duplications. Nine members belonging to the Mβ group, which was considered absent in monocots, have also been identified. The expression profiles of all the MADS-box genes have been analyzed under 11 temporal stages of panicle and seed development, three abiotic stress conditions, along with three stages of vegetative development. Transcripts for 31 genes accumulate preferentially in the reproductive phase, of which, 12 genes are specifically expressed in seeds, and six genes show expression specific to panicle development. Differential expression of seven genes under stress conditions is also evident. An attempt has been made to gain insight into plausible functions of rice MADS-box genes by collating the expression data of functionally validated genes in rice and Arabidopsis.ConclusionOnly a limited number of MADS genes have been functionally validated in rice. A comprehensive annotation and transcriptome profiling undertaken in this investigation adds to our understanding of the involvement of MADS-box family genes during reproductive development and stress in rice and also provides the basis for selection of candidate genes for functional validation studies.

Genomic Survey and Gene Expression Analysis of the Basic Leucine Zipper Transcription Factor Family in Rice

The basic leucine (Leu) zipper (bZIP) proteins compose a family of transcriptional regulators present exclusively in eukaryotes. The bZIP proteins characteristically harbor a bZIP domain composed of two structural features: a DNA-binding basic region and the Leu zipper dimerization region. They have been shown to regulate diverse plant-specific phenomena, including seed maturation and germination, floral induction and development, and photomorphogenesis, and are also involved in stress and hormone signaling. We have identified 89 bZIP transcription factor-encoding genes in the rice (Oryza sativa) genome. Their chromosomal distribution and sequence analyses suggest that the bZIP transcription factor family has evolved via gene duplication. The phylogenetic relationship among rice bZIP domains as well as with bZIP domains from other plant bZIP factors suggests that homologous bZIP domains exist in plants. Similar intron/exon structural patterns were observed in the basic and hinge regions of their bZIP domains. Detailed sequence analysis has been done to identify additional conserved motifs outside the bZIP domain and to predict their DNA-binding site specificity as well as dimerization properties, which has helped classify them into different groups and subfamilies, respectively. Expression of bZIP transcription factor-encoding genes has been analyzed by full-length cDNA and expressed sequence tag-based expression profiling. This expression profiling was complemented by microarray analysis. The results indicate specific or coexpression patterns of rice bZIP transcription factors starting from floral transition to various stages of panicle and seed development. bZIP transcription factor-encoding genes in rice also displayed differential expression patterns in rice seedlings in response to abiotic stress and light irradiation. An effort has been made to link the structure and expression pattern of bZIP transcription factor-encoding genes in rice to their function, based on the information obtained from our analyses and earlier known results. This information will be important for functional characterization of bZIP transcription factors in rice.

A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.).

Cicer arietinum L. (chickpea) is the third most important food legume crop. We have generated the draft sequence of a desi-type chickpea genome using next-generation sequencing platforms, bacterial artificial chromosome end sequences and a genetic map. The 520-Mb assembly covers 70% of the predicted 740-Mb genome length, and more than 80% of the gene space. Genome analysis predicts the presence of 27,571 genes and 210 Mb as repeat elements. The gene expression analysis performed using 274 million RNA-Seq reads identified several tissue-specific and stress-responsive genes. Although segmental duplicated blocks are observed, the chickpea genome does not exhibit any indication of recent whole-genome duplication. Nucleotide diversity analysis provides an assessment of a narrow genetic base within the chickpea cultivars. We have developed a resource for genetic markers by comparing the genome sequences of one wild and three cultivated chickpea genotypes. The draft genome sequence is expected to facilitate genetic enhancement and breeding to develop improved chickpea varieties.

Molecular characterization and differential expression of cytokinin-responsive type-A response regulators in rice (Oryza sativa)

BackgroundThe response regulators represent the elements of bacterial two-component system and have been characterized from dicot plants like Arabidopsis but little information is available on the monocots, including the cereal crops. The aim of this study was to characterize type-A response regulator genes from rice, and to investigate their expression in various organs as well as in response to different hormones, including cytokinin, and environmental stimuli.ResultsBy analysis of the whole genome sequence of rice, we have identified ten genes encoding type-A response regulators based upon their high sequence identity within the receiver domain. The exon-intron organization, intron-phasing as well as chromosomal location of all the RT-PCR amplified rice (Oryza sativa) response regulator (OsRR) genes have been analyzed. The transcripts of OsRR genes could be detected by real-time PCR in all organs of the light- and dark-grown rice seedlings/plants, although there were quantitative differences. The steady-state transcript levels of most of the OsRR genes increased rapidly (within 15 min) on exogenous cytokinin application even in the presence of cycloheximide. Moreover, the expression of the OsRR6 gene was enhanced in rice seedlings exposed to salinity, dehydration and low temperature stress.ConclusionTen type-A response regulator genes identified in rice, the model monocot plant, show overlapping/differential expression patterns in various organs and in response to light. The induction of OsRR genes by cytokinin even in the absence of de novo protein synthesis qualifies them to be primary cytokinin response genes. The induction of OsRR6 in response to different environmental stimuli indicates its role in cross-talk between abiotic stress and cytokinin signaling. These results provide a foundation for further investigations on specific as well as overlapping cellular functions of type-A response regulators in rice.

Gene Discovery and Tissue-Specific Transcriptome Analysis in Chickpea with Massively Parallel Pyrosequencing and Web Resource Development

Chickpea (Cicer arietinum) is an important food legume crop but lags in the availability of genomic resources. In this study, we have generated about 2 million high-quality sequences of average length of 372 bp using pyrosequencing technology. The optimization of de novo assembly clearly indicated that hybrid assembly of long-read and short-read primary assemblies gave better results. The hybrid assembly generated a set of 34,760 transcripts with an average length of 1,020 bp representing about 4.8% (35.5 Mb) of the total chickpea genome. We identified more than 4,000 simple sequence repeats, which can be developed as functional molecular markers in chickpea. Putative function and Gene Ontology terms were assigned to at least 73.2% and 71.0% of chickpea transcripts, respectively. We have also identified several chickpea transcripts that showed tissue-specific expression and validated the results using real-time polymerase chain reaction analysis. Based on sequence comparison with other species within the plant kingdom, we identified two sets of lineage-specific genes, including those conserved in the Fabaceae family (legume specific) and those lacking significant similarity with any non chickpea species (chickpea specific). Finally, we have developed a Web resource, Chickpea Transcriptome Database, which provides public access to the data and results reported in this study. The strategy for optimization of de novo assembly presented here may further facilitate the transcriptome sequencing and characterization in other organisms. Most importantly, the data and results reported in this study will help to accelerate research in various areas of genomics and implementing breeding programs in chickpea.

Structure and expression analysis of early auxin-responsive Aux/IAA gene family in rice (Oryza sativa)

Auxin exerts pleiotropic effects on plant growth and development by regulating the expression of early auxin-responsive genes of auxin/indoleacetic acid (Aux/IAA), small auxin-up RNA, and GH3 classes. These genes have been studied extensively in dicots like soybean and Arabidopsis. We had earlier characterized a cDNA of the first monocot member of Aux/IAA family from rice. The achievement of the large scale rice genome sequencing combined with the availability of full-length cDNA sequences from Knowledge-based Oryza Molecular Biological Encyclopedia provided us the opportunity to draw up the first comprehensive list of Aux/IAA genes in a monocot. By screening the available databases, we have identified 31 Aux/IAA genes having high sequence identity within the conserved domains I, II, III, and IV. The genomic organization as well as chromosomal location of all the Oryza sativa indoleacetic acid (OsIAA) genes is reported. The rice Aux/IAA proteins can be classified in two groups (A and B) on the basis of their phylogenetic relationship with Arabidopsis Aux/IAA proteins. An evolutionary pattern of the rice Aux/IAA genes has been discussed by analyzing their structure (exon/intron organization) and duplications. Interestingly, the duplication of rice Aux/IAA genes was found to be associated with chromosomal block duplication events in rice. The in-silico analysis has been complemented with real-time polymerase chain reaction analysis to quantify transcript levels of all Aux/IAA family members. OsIAA genes showed differential and overlapping organ-specific expression patterns in light- and dark-grown seedlings/plants. Although auxin enhanced the transcript abundance of most of the OsIAA genes, the effect was more pronounced on OsIAA9, 14, 19, 20, 24, and 31. These results provide a foundation for future studies on elucidating the precise role of rice Aux/IAA genes in early steps of auxin signal transduction.

Genome-wide identification, organization and phylogenetic analysis of Dicer-like, Argonaute and RNA-dependent RNA Polymerase gene families and their expression analysis during reproductive development and stress in rice

BackgroundImportant developmental processes in both plants and animals are partly regulated by genes whose expression is modulated at the post-transcriptional level by processes such as RNA interference (RNAi). Dicers, Argonautes and RNA-dependent RNA polymerases (RDR) form the core components that facilitate gene silencing and have been implicated in the initiation and maintenance of the trigger RNA molecules, central to process of RNAi. Investigations in eukaryotes have revealed that these proteins are encoded by variable number of genes with plants showing relatively higher number in each gene family. To date, no systematic expression profiling of these genes in any of the organisms has been reported.ResultsIn this study, we provide a complete analysis of rice Dicer-like, Argonaute and RDR gene families including gene structure, genomic localization and phylogenetic relatedness among gene family members. We also present microarray-based expression profiling of these genes during 14 stages of reproductive and 5 stages of vegetative development and in response to cold, salt and dehydration stress. We have identified 8 Dicer-like (OsDCLs), 19 Argonaute (OsAGOs) and 5 RNA-dependent RNA polymerase (OsRDRs) genes in rice. Based on phylogeny, each of these genes families have been categorized into four subgroups. Although most of the genes express both in vegetative and reproductive organs, 2 OsDCLs, 14 OsAGOs and 3 OsRDRs were found to express specifically/preferentially during stages of reproductive development. Of these, 2 OsAGOs exhibited preferential up-regulation in seeds. One of the Argonautes (OsAGO2) also showed specific up-regulation in response to cold, salt and dehydration stress.ConclusionThis investigation has identified 23 rice genes belonging to DCL, Argonaute and RDR gene families that could potentially be involved in reproductive development-specific gene regulatory mechanisms. These data provide an insight into probable domains of activity of these genes and a basis for further, more detailed investigations aimed at understanding the contribution of individual components of RNA silencing machinery during reproductive phase of plant development.

Expression analysis of calcium-dependent protein kinase gene family during reproductive development and abiotic stress conditions in rice (Oryza sativa L. ssp. indica)

Calcium-dependent protein kinases (CDPKs) are important sensors of Ca+2 flux in plants, which control plant development and responses by regulating downstream components of calcium signaling pathways. Availability of the whole genome sequence and microarray platform allows investigation of genome-wide organization and expression profile of CDPK genes in rice with a view to ultimately define their function in plant systems. Genome-wide analysis led to identification of 31 CDPK genes in rice after a thorough annotation exercise based upon HMM profiles. Twenty-nine already identified CDPK genes were verified and two new members were added to the CDPK gene family of rice. Relative expression of all these genes has been analyzed by using Affymetrix rice genome array™ during three vegetative stages, six stages of panicle (P1–P6) and five stages of seed (S1–S5) development along with three abiotic stress conditions, viz. cold, salt and desiccation, given to seedling. Thirty-one CDPK genes were found to express in at least one of the experimental stages studied. Of these, transcripts for twenty three genes accumulated differentially during reproductive developmental stages; nine of them were preferentially up-regulated only in panicle, five were up-regulated in stages of panicles as well as seed development, whereas, expression of one gene was found to be specific to the S1 stage of seed development. Eight genes were found to be down-regulated during the panicle and seed developmental stages. Six CDPK genes were found to be induced while the expression of one gene was down-regulated under stress conditions. The differential expression of CDPK genes during reproductive development and stress is suggestive of their involvement in the underlying signal transduction pathways. Furthermore, up-regulation of common genes both during reproductive development as well as stress responses is indicative of common element between reproduction and stress.