Genetic analysis for detection of genes associated to drought tolerance in rice accessions belonging to north east India.

INTRODUCTION: The North East (NE) India is rich in biodiversity and also considered as the secondary centre for origin of rice. The NE rice accessions was characterized previously using genetic markers and morphological traits. Simultaneously, genome-wide association studies (GWAS) reveal significant marker-trait associations for the drought tolerance traits. METHODS AND RESULTS: The genetic diversity and population structure of 296 NE rice accessions were studied using 96,712 single nucleotide polymorphism (SNP) markers distributed across 12 chromosomes. The accessions were clustered into two major sub-groups (SG). A total of 91 accessions were assembled as SG1 and 114 accessions as SG2, while the remaining 91 were admixture genotypes. A total of 200 genotypes belonging to different groups were phenotyped for yield component traits under drought and control conditions. The GWAS was performed to identify significant marker-trait associations (MTAs). Consequently, 47 MTAs were detected under drought, exhibiting 0.02-9.95% of phenotypic variance (P.V.). Whereas 58 MTAs were discovered under control conditions, showing a 0.01-9.74% contribution to the phenotype. Through in-silico mining of QTLs, 2999 genes were identified. Among these; only 22 genes were directly associated with stress response. CONCLUSION: These QTLs/genes may be deployed for marker-assisted pyramiding to improve drought tolerance in popular drought susceptible rice varieties.

Identification and characterization of genes for drought tolerance in upland rice cultivar 'Banglami' of North East India.

INTRODUCTION: Rice is a major crop in Assam, North East (NE) India. The rice accessions belonging to NE India possess unique traits of breeder's interest, i.e., tolerant to biotic and abiotic stresses. In the present research programme, the stress responsive genes were identified within the QTLs associated with drought tolerance. The differential expression profiling of genes were performed under drought stress and control conditions. Thus, the 'candidate genes' associated with drought tolerance were recognised and may be deployed in a breeding programme. METHODS AND RESULTS: A drought-tolerant traditional rice cultivar, Banglami, was crossed with a high-yielding, drought-susceptible variety, Ranjit. The mapping population (F4) was raised through the single seed descent (SSD) method and used in QTL analysis. Under drought stress, a total of 4752 genes were identified through in-silico mining of QTLs. Among these, only 21 genes primarily associated with the stress response. The maximum of four stress-responsive genes were located within the QTLs, qNOG12.1 and qGY1.1. However, under control conditions, 2088 genes were identified, out of which, only 15 were categorised as the major stress responsive genes. The functional characterization of genes recognized 24 different types of proteins. Among these, peroxidase and heat shock proteins (Hsp) are the principal proteins encoded during stress. In addition to that, OsbZIP23, inorganic pyrophosphatase, universal stress protein, serine threonine kinase, NADPH oxidoreductase, and proteins belonging to the ABC1 family were also produced during stress condition. The differential expression profiling showed a profound expression pattern of three candidate genes under drought stress condition, i.e., OsI_32199 (Ascorbate peroxidase), OsI_37694 (Universal stress protein) and OsI_32167 (Heat shock protein 81 - 1). CONCLUSION: The novel candidate genes identified for drought tolerance, may be used in the breeding programme for the development of 'climate smart rice varieties'.

In silico characterization and differential expression analysis of 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR) of Centella asiatica.

The 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR; EC1.1.1.267), an NADPH-dependent reductase, plays a pivotal role in the methylerythritol 4-phosphate pathway (MEP), in the conversion of 1-deoxy-d-xylulose-5-phosphate (DXP) into MEP. Photochemical profiles, as well as pharmaceutical activities of Centella asiatica (L.), one of the most valuable medicinal plants, divulge the presence of secondary metabolites called Centellosides. Despite well-studied pharmaceutical activities, not much is known about the genes responsible for the synthesis of these compounds. In the present study, the full-length DXR gene sequence (JQ965955) of Centella submitted in NCBI was characterized using various bioinformatics tools and tissue specific differential expression studies were also carried out. The full-length CDNA of CaDXR contains an open reading frame (ORF) of 1425 bp which encodes a peptide of 474 amino acids. The molecular weight of this protein was found to be 51.5 kDa with isoelectric point of 6.33. The protein contains three conserved domain, namely NADPH (GSTGSIGT and LAAGSNV), substrate binding (LPADSEHSAI and NKGLEVIEAHY) and Cys-Ser-(Ala/Met/Val/Thr) cleavage-site domains. Phylogenetic studies of CaDXR sequence show close homology with DXR sequence of Angelica sinensis and Daucus carota subsp sativus as they all belong to Apiaceae family. In silico analysis predicted that CaDXR protein contains 21 α-helix and 11 ß-sheets and further DXR protein model was validated by Ramachandran plot analysis. The results of molecular dynamics (MD) simulations unveil dynamic stability of the proposed model and docking studies suggest that the NDP cofactor tightly binds in the active site of the protein with a strong network of hydrogen and hydrophobic interactions. The expression studies by semi-RT followed by qRT-PCR suggests that CaDXR is differentially expressed in different tissues (with maximal expression in the node and lowest in the roots). Thus, characterization and structure-function analysis of DXR gene in Centella facilitate us to understand not only the functions of DXR gene but also regulatory mechanisms involved in the MEP pathway in C. asiatica plant at the molecular level. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-02723-w.

An Insight Into Structure, Function, and Expression Analysis of 3-Hydroxy-3-Methylglutaryl-CoA Reductase of Cymbopogon winterianus.

Citronella (Cymbopogon winterianus) is one of the richest sources of high-value isoprenoid aromatic compounds used as flavour, fragrance, and therapeutic elements. These isoprenoid compounds are synthesized by 2 independent pathways: mevalonate pathway and 2-C-methyl-d-erythritol-4-phosphate pathway. Evidence suggests that 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) is a rate-controlling enzyme for the synthesis of variety of isoprenoids. This study reports the isolation, characterization, and tissue-specific expression analysis of HMGR from citronella. The modelled HMGR is a class I type of HMGR enzyme with 3-domain architecture. The active site comprises a cofactor (nicotinamide adenine dinucleotide phosphate) and the substrate-binding motifs. The real-time and quantitative reverse transcription-polymerase chain reaction results revealed equal expression level in both leaf sheath and root tissue. The results from our study shall be a valuable resource for future molecular intervention to alter the metabolic flux towards improvement of key active ingredient in this important medicinal plant.

Molecular Modeling and Dynamics Simulation Analysis of KATNAL1 for Identification of Novel Inhibitor of Sperm Maturation.

BACKGROUND: Hormone based birth control often causes various side effects. A recent study revealed that temporary infertility without changing hormone levels can be attained by inhibiting Katanin p60 ATPase-containing subunit A-like 1 protein (KATNAL1) which is critical for sperm maturation in the testes. OBJECTIVE: This study aimed at attaining the most energetically stable three dimensional (3D) structure of KATNAL1 protein using comparative modeling followed by screening of a ligand library of known natural spermicidal compounds for their binding affinity with KATNAL1. This in turn may inhibit the development of mature sperm in the seminiferous epithelium. METHOD: A series of computational techniques were used for building the 3D structure of KATNAL1 which was further optimized by molecular dynamics (MD) simulation. For revealing the ATP binding mode of KATNAL1, docking study was carried out using the optimized model obtained from the MD simulation. The docking study was also employed to test the binding efficiency of the ligand library. RESULTS: Molecular docking study confirmed the ATP binding of KATNAL1 with various hydrophobic and hydrogen bond interactions. Binding efficiency of the ligand library suggested that calotropin, a cardenolide of Calotropis procera showed the highest binding efficiency against the target protein without toxicity. MD simulation of the docked complex validated the results of the docking study. CONCLUSION: This study revealed the ATP binding mode of KATNAL1 and identified calotropin as a potential lead molecule against it showing high binding efficiency with good bioavailability and no mutagenicity. Further in vitro and in vivo bioassay of calotropin could facilitate the development of novel non-hormonal male-specific contraceptive in near future.

Genome wide transcriptome profiling reveals differential gene expression in secondary metabolite pathway of Cymbopogon winterianus.

Advances in transcriptome sequencing provide fast, cost-effective and reliable approach to generate large expression datasets especially suitable for non-model species to identify putative genes, key pathway and regulatory mechanism. Citronella (Cymbopogon winterianus) is an aromatic medicinal grass used for anti-tumoral, antibacterial, anti-fungal, antiviral, detoxifying and natural insect repellent properties. Despite of having number of utilities, the genes involved in terpenes biosynthetic pathway is not yet clearly elucidated. The present study is a pioneering attempt to generate an exhaustive molecular information of secondary metabolite pathway and to increase genomic resources in Citronella. Using high-throughput RNA-Seq technology, root and leaf transcriptome was analysed at an unprecedented depth (11.7 Gb). Targeted searches identified majority of the genes associated with metabolic pathway and other natural product pathway viz. antibiotics synthesis along with many novel genes. Terpenoid biosynthesis genes comparative expression results were validated for 15 unigenes by RT-PCR and qRT-PCR. Thus the coverage of these transcriptome is comprehensive enough to discover all known genes of major metabolic pathways. This transcriptome dataset can serve as important public information for gene expression, genomics and function genomics studies in Citronella and shall act as a benchmark for future improvement of the crop.

Molecular cloning, characterization and expression analysis of 3-hydroxy-3-methylglutaryl coenzyme A reductase gene from Centella asiatica L.

3-Hydroxy-3-methylglutaryl-CoA reductases (HMGR) plays an important role in catalyzing the first committed step of isoprenoid biosynthesis in the mevelonic (MVA) pathway (catalyzes the conversion of HMG-CoA to MVA) in plants. The present manuscript reports the full length cDNA cloning of HMGR (CaHMGR, GenBank accession number: KJ939450.2) and its characterization from Centella asiatica. Sequence analysis indicated that the cDNA was of 1965 bp, which had an open reading frame of 1617 bp and encoded a protein containing 539 amino-acids with a mol wt of 57.9 kDa. A BLASTp search against non-redundant (nr) protein sequence showed that C. asiatica HMGR (CaHMGR) has 65-81% identity with HMGRs from different plant species and multi-alignment comparison analysis showed the presence of two motif each corresponding to HMG-CoA-binding and NADP(H)-binding. The Conserved Domain Database analysis predicted that CaHMGR belongs to Class I hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase. Three-dimensional modeling confirmed the novelty of CaHMGR with a spatial structure similar to Homo sapiens (PDB id: 1IDQ8_A). Tissue Expression analysis indicates that CaHMGR is ubiquitous albeit differentially expressed among different tissues analysed, Strong expression was recorded in the nodes and leaves and low in the roots. The present investigation confirmed that nodes are vital to terpenoid synthesis in C. asiatica. Thus, the cloning of full length CDS, characterization and structure-function analysis of HMGR gene in Centella facilitate to understand the HMGR's functions and regulatory mechanisms involved in mevalonate pathway in C. asiatica at genetic level.

Structural comparison, substrate specificity, and inhibitor binding of AGPase small subunit from monocot and dicot: present insight and future potential.

ADP-glucose pyrophosphorylase (AGPase) is the first rate limiting enzyme of starch biosynthesis pathway and has been exploited as the target for greater starch yield in several plants. The structure-function analysis and substrate binding specificity of AGPase have provided enormous potential for understanding the role of specific amino acid or motifs responsible for allosteric regulation and catalytic mechanisms, which facilitate the engineering of AGPases. We report the three-dimensional structure, substrate, and inhibitor binding specificity of AGPase small subunit from different monocot and dicot crop plants. Both monocot and dicot subunits were found to exploit similar interactions with the substrate and inhibitor molecule as in the case of their closest homologue potato tuber AGPase small subunit. Comparative sequence and structural analysis followed by molecular docking and electrostatic surface potential analysis reveal that rearrangements of secondary structure elements, substrate, and inhibitor binding residues are strongly conserved and follow common folding pattern and orientation within monocot and dicot displaying a similar mode of allosteric regulation and catalytic mechanism. The results from this study along with site-directed mutagenesis complemented by molecular dynamics simulation will shed more light on increasing the starch content of crop plants to ensure the food security worldwide.

Computational identification and characterization of conserved miRNAs and their target genes in garlic (Allium sativum L.) expressed sequence tags.

The endogenous small non-coding functional microRNAs (miRNAs) are short in size, range from ~21 to 24 nucleotides in length, play a pivotal role in gene expression in plants and animals by silencing genes either by destructing or blocking of translation of homologous mRNA. Although various high-throughput, time consuming and expensive techniques like forward genetics and direct cloning are employed to detect miRNAs in plants but comparative genomics complemented with novel bioinformatic tools pave the way for efficient and cost-effective identification of miRNAs through homologous sequence search with previously known miRNAs. In this study, an attempt was made to identify and characterize conserved miRNAs in garlic expressed sequence tags (ESTs) through computational means. For identification of novel miRNAs in garlic, a total 3227 known mature miRNAs of plant kingdom Viridiplantae were searched for homology against 21,637 EST sequences resulting in identification of 6 potential miRNA candidates belonging to 6 different miRNA families. The psRNATarget server predicted 33 potential target genes and their probable functions for the six identified miRNA families in garlic. Most of the garlic miRNA target genes seem to encode transcription factors as well as genes involved in stress response, metabolism, plant growth and development. The results from the present study will shed more light on the understanding of molecular mechanisms of miRNA in garlic which may aid in the development of novel and precise techniques to understand some post-transcriptional gene silencing mechanism in response to stress tolerance.

Structural analysis and molecular dynamics simulations of novel δ-endotoxin Cry1Id from Bacillus thuringiensis to pave the way for development of novel fusion proteins against insect pests of crops.

The theoretical three-dimensional structure of a novel δ-endotoxin Cry1Id (81 kDa) belonging to Cry1I class, toxic to many of the lepidopteran pests has been investigated through comparative modeling. Molecular dynamics (MD) simulations was carried out to characterize its structural and dynamical features at 10 ns in explicit solvent using the GROMACS version 4.5.4. Finally the simulated model was validated by the SAVES, WHAT IF, MetaMQAP, ProQ, ModFOLD and MolProbity servers. Despite low sequence identity with its structural homologs, Cry1Id not only resembles the previously reported Cry structures but also shares the common five conserved blocks of amino acid residues. Although the domain II of Cry1Id superpose well with its closest structural homolog Cry8Ea1, variation of amino acids and length in the apical loop2 of domain II was observed. In this work, we have hypothesized that the variations in apical loop2 might be the sole factor for providing variable surface accessibility to Cry1Id protein that could be important in receptor recognition. MD simulation showed the proposed endotoxin retains its stable conformation in aqueous solution. The result from this study is expected to aid in the development hybrid Cry proteins and new potent fusion proteins with novel specificities against different insect pests for improved pest management of crop plants.

Construction of cDNA library and preliminary analysis of expressed sequence tags from tea plant [Camellia sinensis (L) O. Kuntze].

Tea is the most popular non-alcoholic and healthy beverage across the world. The understanding of the genetic organization and molecular biology of tea plant, which is very poorly understood at present, is required for quantum increase in productivity and efficient use of germplasm for either cultivation or breeding program. Single-pass sequencing of randomly selected cDNA clones is the most widely accepted technique for gene identification and cloning. In the present study, a good quality cDNA library was constructed and preliminary analysis of ESTs was carried out. The titers of unamplified and amplified libraries were 1.4 × 10(6)pfu/ml and 5.27 × 10(8)pfu/ml respectively. A total of 210 cDNA clones from the constructed cDNA library were sequenced and analyzed. A total of 84 high quality Expressed Sequence Tags (ESTs) were generated, among which 71 ESTs had significant homology with sequences in NCBI non-redundant protein database by BLAST X analysis. About 80% ESTs had poly (A) tail at 3' end indicating that the cDNAs were full length. The database-matched ESTs were classified into putative cellular roles, viz. energy-related category (corresponding to 20% of total BLAST X matched ESTs), Transcription (14.2%), protein synthesis (14.2%) cell growth and division (8.6%), cell structure (5.7%), signal transduction (5.7%), transporters (2.9%), disease and defenses (2.9%), secondary metabolism (2.9%) and gene regulation (2.9%). This study provides an overview of the mRNA expression profile and first hand information of gene sequence expressed in tender leaves and apical buds of tea plant.

Isolation and characterization of cDNA clones encoding ADP-glucose pyrophosphorylase (AGPase) large and small subunits from chickpea (Cicer arietinum L.).

Four cDNA clones encoding two large subunits and two small subunits of the starch regulatory enzyme ADP-glucose pyrophosphorylase (AGPase) were isolated from a chickpea (Cicer arietinum L.) stem cDNA library. DNA sequence and Southern blot analyses of these clones, designated CagpL1, CagpL2 (large subunits) and CagpS1 and CagpS2 (small subunits), revealed that these isoforms represented different AGPase large and small subunits. RNA expression analysis indicated that CagpL1 was expressed strongly in leaves with reduced expression in the stem. No detectable expression was observed in seeds and roots. CagpL2 was expressed moderately in seeds followed by weak expression in leaves, stems and roots. Similar analysis showed that CagpS1 and CagpS2 displayed a spatial expression pattern similar to that observed for CagpL2 with the exception that CagpS1 showed a much higher expression in seeds than CagpS2. The spatial expression patterns of these different AGPase subunit sequences indicate that different AGPase isoforms are used to control starch biosynthesis in different organs during chickpea development.