Anil Rai

Synergistic regulatory networks mediated by microRNAs and transcription factors under drought, heat and salt stresses in Oryza Sativa spp.

BACKGROUND Transcription factors (TFs) and microRNAs (miRNAs) are primary gene regulators within the cell. Regulatory mechanisms of these two main regulators are of great interest to biologists and may provide insights into the abiotic and biotic stresses. However, the interaction between miRNAs and TFs in a gene regulatory network (GRN) still remains uncovered. Previous research has been mostly directed at inferring either miRNA or TF regulatory networks from data. However, networks involving a single type of regulator may not fully reveal the complex gene regulatory mechanisms, therefore study of interplay among these two regulators in gene regulation is important towards explaining the mechanism of different abiotic stresses. RESULT Oligonucleotide microarrays containing 51,279 transcripts were used to identify total 133 salt responsive target genes regulated by 11 TFs that are also differentially regulated by miRNA under salinity, heat and drought stresses in Oryza sativa. TFs-target interactions which are most enriched in their downstream regulation were also identified. Many genes whose encoded proteins are implicated in response to light and radiation stimulus, hormone stimuli, oxidative stress, copper ion binding and electron transport were found to be enriched. However the majority were novel for the combined abiotic stress, which indicates that there are a great number of genes induced after the exposure these abiotic stresses and regulated by miRNA. CONCLUSION Analysis of the expression profile data of Oryza provides clues regarding some putative cellular and molecular processes that are undertaken in response to these stresses. The study also identified a large number of candidate functional genes that appear to be constitutively involved in salt, drought and heat stresses tolerance. Further examination of these genes may enable the molecular basis of abiotic stress tolerance in Oryza, to be elucidated.

First whole genome based microsatellite DNA marker database of tomato for mapping and variety identification

BackgroundThe cultivated tomato is second most consumed vegetable of the world and is an important part of a diverse and balanced diet as a rich source of vitamins, minerals, phenolic antioxidants and antioxidant lycopene having anti-cancer properties. To reap benefit of genomics of the domestic tomato (Solanum lycopersicum L.) unravelled by Tomato Genome Consortium (The Tomato Genome Consortium, 2012), the bulk mining of its markers in totality is imperative and critically required. The solgenomics has limited number of microsatellite DNA markers (2867) pertaining to solanaceae family. As these markers are of linkage map having relative distance, the choice of selected markers based on absolute distance as of physical map is missing. Only limited microsatellite markers with limitations are reported for variety identification thus there is a need for more markers supplementing DUS test and also for traceability of product in global market.DescriptionWe present here the first whole genome based microsatellite DNA marker database of tomato, TomSatDB (Tomato MicroSatellite Database) with more than 1.4 million markers mined in-silico, using MIcroSAtellite (MISA) tool. To cater the customized needs of wet lab, features with a novelty of an automated primer designing tool is added. TomSatDB (http://cabindb.iasri.res.in/tomsatdb), a user-friendly and freely accessible tool offers chromosome wise as well as location wise search of primers. It is an online relational database based on “three-tier architecture” that catalogues information of microsatellites in MySQL and user-friendly interface developed using PHP (Hypertext Pre Processor).ConclusionBesides abiotic stress, tomato is known to have biotic stress due to its susceptibility over 200 diseases caused by pathogenic fungi, bacteria, viruses and nematodes. These markers are expected to pave the way of germplasm management over abiotic and biotic stress as well as improvement through molecular breeding, leading to increased tomato productivity in India as well as other parts of the world. In era of IPR the new variety can be identified based on allelic variation among varieties supplementing DUS test and product traceability.

How far in-silico computing meets real experiments. A study on the structure and dynamics of spin labeled vinculin tail protein by molecular dynamics simulations and EPR spectroscopy

BackgroundInvestigation of conformational changes in a protein is a prerequisite to understand its biological function. To explore these conformational changes in proteins we developed a strategy with the combination of molecular dynamics (MD) simulations and electron paramagnetic resonance (EPR) spectroscopy. The major goal of this work is to investigate how far computer simulations can meet the experiments.MethodsVinculin tail protein is chosen as a model system as conformational changes within the vinculin protein are believed to be important for its biological function at the sites of cell adhesion. MD simulations were performed on vinculin tail protein both in water and in vacuo environments. EPR experimental data is compared with those of the simulated data for corresponding spin label positions.ResultsThe calculated EPR spectra from MD simulations trajectories of selected spin labelled positions are comparable to experimental EPR spectra. The results show that the information contained in the spin label mobility provides a powerful means of mapping protein folds and their conformational changes.ConclusionsThe results suggest the localization of dynamic and flexible regions of the vinculin tail protein. This study shows MD simulations can be used as a complementary tool to interpret experimental EPR data.

PIPEMicroDB: microsatellite database and primer generation tool for pigeonpea genome

Molecular markers play a significant role for crop improvement in desirable characteristics, such as high yield, resistance to disease and others that will benefit the crop in long term. Pigeonpea (Cajanus cajan L.) is the recently sequenced legume by global consortium led by ICRISAT (Hyderabad, India) and been analysed for gene prediction, synteny maps, markers, etc. We present PIgeonPEa Microsatellite DataBase (PIPEMicroDB) with an automated primer designing tool for pigeonpea genome, based on chromosome wise as well as location wise search of primers. Total of 123 387 Short Tandem Repeats (STRs) were extracted from pigeonpea genome, available in public domain using MIcroSAtellite tool (MISA). The database is an online relational database based on ‘three-tier architecture’ that catalogues information of microsatellites in MySQL and user-friendly interface is developed using PHP. Search for STRs may be customized by limiting their location on chromosome as well as number of markers in that range. This is a novel approach and is not been implemented in any of the existing marker database. This database has been further appended with Primer3 for primer designing of selected markers with left and right flankings of size up to 500 bp. This will enable researchers to select markers of choice at desired interval over the chromosome. Furthermore, one can use individual STRs of a targeted region over chromosome to narrow down location of gene of interest or linked Quantitative Trait Loci (QTLs). Although it is an in silico approach, markers’ search based on characteristics and location of STRs is expected to be beneficial for researchers. Database URL: http://cabindb.iasri.res.in/pigeonpea/

Time series analysis of groundwater levels and projection of future trend

A study was under taken for identifying the trends in pre and post-monsoon groundwater levels using Mann- Kendall test and Sen’s slope estimator, and for time series modelling of groundwater levels for forecasting the pre and post-monsoon water levels in Karnal district of Haryana. Results showed that the groundwater levels had significantly declined during 1974 to 2010. Average rates of water level decline were 0.228 and 0.267 m/yr during pre and postmonsoon seasons, respectively. There was rapid decline in water level between 2001 and 2010. The ARIMA (0, 1, 2) was identified as the appropriate model for time series modelling and forecasting. Results showed that the pre and postmonsoon groundwater level in 2050 would decline by 12.97 m and 12.00 m over the observed water level in 2010, and reach to a level of 29.95 m and 28.14 m below ground surface. The average rate of decline of pre and post-monsoon groundwater level in the district during this period would be 0.32 and 0.30 m/yr, respectively.

In silico mining of putative microsatellite markers from whole genome sequence of water buffalo (Bubalus bubalis) and development of first BuffSatDB

BackgroundThough India has sequenced water buffalo genome but its draft assembly is based on cattle genome BTau 4.0, thus de novo chromosome wise assembly is a major pending issue for global community. The existing radiation hybrid of buffalo and these reported STR can be used further in final gap plugging and “finishing” expected in de novo genome assembly. QTL and gene mapping needs mining of putative STR from buffalo genome at equal interval on each and every chromosome. Such markers have potential role in improvement of desirable characteristics, such as high milk yields, resistance to diseases, high growth rate. The STR mining from whole genome and development of user friendly database is yet to be done to reap the benefit of whole genome sequence.DescriptionBy in silico microsatellite mining of whole genome, we have developed first STR database of water buffalo, BuffSatDb (Buffalo MicroSatellite Database (http://cabindb.iasri.res.in/buffsatdb/) which is a web based relational database of 910529 microsatellite markers, developed using PHP and MySQL database. Microsatellite markers have been generated using MIcroSAtellite tool. It is simple and systematic web based search for customised retrieval of chromosome wise and genome-wide microsatellites. Search has been enabled based on chromosomes, motif type (mono-hexa), repeat motif and repeat kind (simple and composite). The search may be customised by limiting location of STR on chromosome as well as number of markers in that range. This is a novel approach and not been implemented in any of the existing marker database. This database has been further appended with Primer3 for primer designing of the selected markers enabling researcher to select markers of choice at desired interval over the chromosome. The unique add-on of degenerate bases further helps in resolving presence of degenerate bases in current buffalo assembly.ConclusionBeing first buffalo STR database in the world , this would not only pave the way in resolving current assembly problem but shall be of immense use for global community in QTL/gene mapping critically required to increase knowledge in the endeavour to increase buffalo productivity, especially for third world country where rural economy is significantly dependent on buffalo productivity.

Identification of Putative RuBisCo Activase (TaRca1)—The Catalytic Chaperone Regulating Carbon Assimilatory Pathway in Wheat (Triticum aestivum) under the Heat Stress

RuBisCo activase (Rca) is a catalytic chaperone involved in modulating the activity of RuBisCo (key enzyme of photosynthetic pathway). Here, we identified eight novel transcripts from wheat through data mining predicted to be Rca and cloned a transcript of 1.4 kb from cv. HD2985, named as TaRca1 (GenBank acc. no. KC776912). Single copy number of TaRca1 was observed in wheat genome. Expression analysis in diverse wheat genotypes (HD2985, Halna, PBW621, and HD2329) showed very high relative expression of TaRca1 in Halna under control and HS-treated, as compared to other cultivars at different stages of growth. TaRca1 protein was predicted to be chloroplast-localized with numerous potential phosphorylation sites. Northern blot analysis showed maximum accumulation of TaRca1 transcript in thermotolerant cv. during mealy-ripe stage, as compared to thermosusceptible. Decrease in the photosynthetic parameters was observed in all the cultivars, except PBW621 in response to HS. We observed significant increase in the Rca activity in all the cultivars under HS at different stages of growth. HS causes decrease in the RuBisCo activity; maximum reduction was observed during pollination stage in thermosusceptible cvs. as validated through immunoblotting. We observed uniform carbon distribution in different tissues of thermotolerant cvs., as compared to thermosusceptible. Similarly, tolerance level of leaf was observed maximum in Halna having high Rca activity under HS. A positive correlation was observed between the transcript and activity of TaRca1 in HS-treated Halna. Similarly, TaRca1 enzyme showed positive correlation with the activity of RuBisCo. There is, however, need to manipulate the thermal stability of TaRca1 enzyme through protein engineering for sustaining the photosynthetic rate under HS—a novel approach toward development of “climate-smart” crop.

Design and development of portal for biological database in agriculture

The application of novel and modern techniques in genetic engineering and genomics has resulted in information explosion in genomics. Three major genome databases under International Nucleotide Sequence Database collaboration NCBI, DDBJ and EMBL have been providing a convenient platform for submission of sequences which they share among themselves. Many institutes in India under Indian Council of Agricultural Research have scientists working on biotechnology and bioinformatics research. The various studies conducted by them, generate massive data related to biological information of plants, animals, insects, microbes and fisheries. These scientists are dependent on NCBI, EMBL, DDBJ and other portals for their sequence submissions, analysis and other data mining tasks. Due to various limitations imposed on these sites and the poor connectivity problem prevents them to conduct their studies on these open domain databases. The valued information generated by them needs to be shared by the scientific communities to eliminate the duplication of efforts and expedite their knowledge extended towards new findings. A secured common submission portal system with user-friendly interfaces, integrated help and error checking facilities has been developed in such a way that the database at the backend consists of a union of the items available on the above mentioned databases. Standard database management concepts have been employed for their systematic storage management. Extensive hardware resources in the form of high performance computing facility are being installed for deployment of this portal. Availability http://cabindb.iasri.res.in:8080/sequence_portal/

Synonymous codon usage in Thermosynechococcus elongatus (cyanobacteria) identifies the factors shaping codon usage variation.

Analysis of synonymous codon usage pattern in the genome of a thermophilic cyanobacterium, Thermosynechococcus elongatus BP-1 using multivariate statistical analysis revealed a single major explanatory axis accounting for codon usage variation in the organism. This axis is correlated with the GC content at third base of synonymous codons (GC3s) in correspondence analysis taking T. elongatus genes. A negative correlation was observed between effective number of codons i.e. Nc and GC3s. Results suggested a mutational bias as the major factor in shaping codon usage in this cyanobacterium. In comparison to the lowly expressed genes, highly expressed genes of this organism possess significantly higher proportion of pyrimidine-ending codons suggesting that besides, mutational bias, translational selection also influenced codon usage variation in T. elongatus. Correspondence analysis of relative synonymous codon usage (RSCU) with A, T, G, C at third positions (A3s, T3s, G3s, C3s, respectively) also supported this fact and expression levels of genes and gene length also influenced codon usage. A role of translational accuracy was identified in dictating the codon usage variation of this genome. Results indicated that although mutational bias is the major factor in shaping codon usage in T. elongatus, factors like translational selection, translational accuracy and gene expression level also influenced codon usage variation.

Computational analysis of miRNA-target community network reveals cross talk among different metabolisms

To date, only a few conserved miRNAs have been predicted in hexaploid (AABBDD) bread wheat and till now community behavior among miRNA is still in dark. Analysis of publically available 1287279 ESTs from NCBI resulted 262 putative pre-miRNAs and 39 novel mature miRNAs. A total 22,468 targets were identified on 21 chromosomes. MiRNA target community was identified for genomes with different levels of cross talks. Gene ontology of these community targets suggests their differential involvement in different metabolisms along with common and stringent involvement in nitrogen metabolism.