Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing.

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na+ -K+ /H+ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

Modulation of miRNA expression in natural populations of A. thaliana along a wide altitudinal gradient of Indian Himalayas.

Plant populations growing along an altitudinal gradient are exposed to different environmental conditions. They are excellent resources to study regulatory mechanisms adopted by plants to respond to different environmental stresses. Regulation by miRNA is one of such strategies. Here, we report how different miRNAs are preferentially expressed in the three natural populations of A. thaliana originating from a wide altitudinal range. The expression level of miRNAs was mostly governed by temperature and radiation. Majority of the identified miRNAs expressed commonly in the three populations. However, 30 miRNAs expressed significantly at different level between the low and the high altitude populations. Most of these miRNAs regulate the genes associated with different developmental processes, abiotic stresses including UV, cold, secondary metabolites, etc. Further, the expression of miR397 and miR858 involved in lignin biosynthesis and regulation of secondary metabolites respectively, may be regulated by light intensity. A few miRNAs expressed at increasing level with the increase in the altitude of the site indicating environment driven tight regulation of these miRNAs. Further, several novel miRNAs and isomiR diversity specific to the Himalayas are reported which might have an adaptive advantage. To the best of our knowledge, this is the first report on miRNA expression from natural plant populations.

Global gene expression and pigment analysis of two contrasting flower color cultivars of Canna.

Development of flower color in plants is a complex process. Among others, it is an important trait for ornamental flowering plants. Canna is a flowering ornamental plant of family Cannaceae. To understand the molecular mechanism of flower color development in Canna, RNA sequencing from flower tissues of two contrasting flower color cultivars, Red President (RP) and Tropical Sunrise (TS) was performed. More than 27.0 million and 19.0 million clean reads were obtained from RP and TS, respectively. The combined clean reads were assembled into 147,295 unigenes. The Canna unigenes showed maximum homology with Populus trichocarpa (26.79%). A total of 2702 unigenes expressed differentially between the two cultivars of which 1972 were up-regulated and 730 were down-regulated in RP. Phenylpropanoid and flavonoid biosynthetic processes were the significant processes in RP. Expression of a vast number of transcription factors including MYB, bHLH, ARF, and WRKY were higher in RP than TS. The expression analysis of RNA sequencing data was validated by qRT-PCR analysis. Further, concentration of measured anthocyanidins and flavonols were very low or absent in TS, corroborating largely with our transcriptome data. These findings may help in understanding flower color development in Canna and in future crop breeding program.

High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition.

BACKGROUND: Population differentiation and their adaptation to a particular environment depend on their ability to respond to a new environment. This, in turn is governed to an extent, by the degree of phenotypic plasticity exhibited by the populations. The populations of same species inhabiting different climatic conditions may differ in their phenotypic plasticity. Himalayan populations of Arabidopsis thaliana originating from a steep altitude are exposed to different climatic conditions ranging from sub-tropical to temperate. Thus they might have experienced different selection pressures during evolution and may respond differently under common environmental condition. RESULTS: Phenotypic plasticity and differentiation of natural populations of A. thaliana grown under common garden and controlled conditions were determined. A total of seventeen morphological traits, their plasticity, association between traits and environment were performed using 45 accessions from three populations. Plants from different altitudes differed in phenotypes, their selection and fitness under two conditions. Under both the conditions lower altitude population was characterized by higher leaf count and larger silique than higher and middle altitude population. Flowering time of high altitude population increased while that of low and medium altitude decreased under controlled condition compared to open field. An increase in seed weight and germination was observed for all the population under open field than controlled. Rosette area was under divergent selection in both the condition. The heritability of lower altitude population was the highest under both the conditions, where as it was the least for higher altitude population further indicating that the high altitude populations are more responsive towards phenotypic changes under new environmental conditions. Ninety-nine percent of variability in traits and their plasticity co-varied with the altitude of their origin. The population of high altitude was more plastic and differentiated as compared to the lower altitude one. CONCLUSIONS: Arabidopsis thaliana population native to different altitudes of the west Himalaya responds differently when grown under common environments. The success of high altitude population is more in common garden than the controlled conditions. The significant variability in phenotype and its association with altitude of origin predicts for non-random genetic differentiation among the populations.

High light intensity plays a major role in emergence of population level variation in Arabidopsis thaliana along an altitudinal gradient.

Environmental conditions play an important role in the emergence of genetic variations in natural populations. We identified genome-wide patterns of nucleotide variations in the coding regions of natural Arabidopsis thaliana populations. These populations originated from 700 m to 3400 m a.m.s.l. in the Western Himalaya. Using a pooled RNA-Seq approach, we identified the local and global level population-specific SNPs. The biological functions of the SNP-containing genes were primarily related to the high light intensity prevalent at high-altitude regions. The novel SNPs identified in these genes might have arisen de novo in these populations. In another approach, the FSTs of SNP-containing genes were correlated with the corresponding climatic factors. 'Radiation in the growing season' was the only environmental factor found to be strongly correlated with the gene-level FSTs. In both the approaches, the high light intensity was identified as the primary abiotic stress associated with the variations in these populations. The differential gene expression analysis between field and controlled condition grown plants also showed high light intensity as the primary abiotic stress, particularly for the high altitude populations. Our results provide a genome-wide perspective of nucleotide variations in populations along altitudinal gradient and their putative role in emergence of these variations.

Identification and Expression Analyses of miRNAs from Two Contrasting Flower Color Cultivars of Canna by Deep Sequencing.

miRNAs are endogenous small RNA (sRNA) that play critical roles in plant development processes. Canna is an ornamental plant belonging to family Cannaceae. Here, we report for the first time the identification and differential expression of miRNAs in two contrasting flower color cultivars of Canna, Tropical sunrise and Red president. A total of 313 known miRNAs belonging to 78 miRNA families were identified from both the cultivars. Thirty one miRNAs (17 miRNA families) were specific to Tropical sunrise and 43 miRNAs (10 miRNA families) were specific to Red president. Thirty two and 18 putative new miRNAs were identified from Tropical sunrise and Red president, respectively. One hundred and nine miRNAs were differentially expressed in the two cultivars targeting 1343 genes. Among these, 16 miRNAs families targeting 60 genes were involved in flower development related traits and five miRNA families targeting five genes were involved in phenyl propanoid and pigment metabolic processes. We further validated the expression analysis of a few miRNA and their target genes by qRT-PCR. Transcription factors were the major miRNA targets identified. Target validation of a few randomly selected miRNAs by RLM-RACE was performed but was successful with only miR162. These findings will help in understanding flower development processes, particularly the color development in Canna.

Genetic diversity and population structure of Arabidopsis thaliana along an altitudinal gradient.

The natural genetic variation within a plant species is primarily a consequence of its phylogeography and evolutionary history. This variation largely determines its present-day population structure. Arabidopsis thaliana, as a model plant, has been studied in great detail including its probable origin, local as well as global genetic diversity pattern, population structure, adaptation, etc. However, no such studies have so far been reported from the Indian Himalayan region. Here, we describe a comprehensive study on the genetic diversity and population structure of A. thaliana from an altitudinal range of 700-3400 m above mean sea level the highest altitudinal range reported so far. We also compare these populations with previously reported worldwide populations. A total of 48 accessions representing six populations were analysed using 19 microsatellites and 11 chloroplast markers. Genetic diversity analysis indicated populations to be highly diverse and comparable with worldwide populations. STRUCTURE, principal coordinate and isolation by distance (IBD) analyses showed that genetic variation in different populations is structured at geographical and altitudinal level. Further analyses indicate that these populations are genetically distinct from the rest of the world populations. Different parameters of the demographic expansion model support a rapid expansion. Based on mismatch distribution, the initial time of expansion of west Himalayan populations was found to be about 130 000 years. Bayesian analysis of divergence time indicated that these populations have a long evolutionary history in this region. Based on the results of genetic diversity parameters, demographic expansion and divergence time estimation, it appears that west Himalayan populations may be the source of the west-east expansion model.

Expression of rabies glycoprotein and ricin toxin B chain (RGP-RTB) fusion protein in tomato hairy roots: a step towards oral vaccination for rabies.

Transgenic hairy roots of Solanum lycopersicum were engineered to express a recombinant protein containing a fusion of rabies glycoprotein and ricin toxin B chain (rgp-rtxB) antigen under the control of constitutive CaMV35S promoter. Asialofetuin-mediated direct ELISA of transgenic hairy root extracts was performed using polyclonal anti-rabies antibodies (Ab1) and epitope-specific peptidal anti-RGP (Ab2) antibodies which confirmed the expression of functionally viable RGP-RTB fusion protein. Direct ELISA based on asialofetuin-binding activity was used to screen crude protein extracts from five transgenic hairy root lines. Expressions of RGP-RTB fusion protein in different tomato hairy root lines varied between 1.4 and 8 µg in per gram of tissue. Immunoblotting assay of RGP-RTB fusion protein from these lines showed a protein band on monomeric size of ~84 kDa after denaturation. Tomato hairy root line H03 showed highest level of RGP-RTB protein expression (1.14 %) and was used further in bench-top bioreactor for the optimization of scale-up process to produce large quantity of recombinant protein. Partially purified RGP-RTB fusion protein was able to induce the immune response in BALB/c mice after intra-mucosal immunization. In the present investigation, we have not only successfully scaled up the hairy root culture but also established the utility of this system to produce vaccine antigen which subsequently will reduce the total production cost for implementing rabies vaccination programs in developing nations. This study in a way aims to provide consolidated base for low-cost preparation of improved oral vaccine against rabies.

The internal transcribed spacer (ITS) region and trnH-psbA [corrected] are suitable candidate loci for DNA barcoding of tropical tree species of India.

BACKGROUND: DNA barcoding as a tool for species identification has been successful in animals and other organisms, including certain groups of plants. The exploration of this new tool for species identification, particularly in tree species, is very scanty from biodiversity-rich countries like India. rbcL and matK are standard barcode loci while ITS, and trnH-psbA are considered as supplementary loci for plants. METHODOLOGY AND PRINCIPAL FINDINGS: Plant barcode loci, namely, rbcL, matK, ITS, trnH-psbA, and the recently proposed ITS2, were tested for their efficacy as barcode loci using 300 accessions of tropical tree species. We tested these loci for PCR, sequencing success, and species discrimination ability using three methods. rbcL was the best locus as far as PCR and sequencing success rate were concerned, but not for the species discrimination ability of tropical tree species. ITS and trnH-psbA were the second best loci in PCR and sequencing success, respectively. The species discrimination ability of ITS ranged from 24.4 percent to 74.3 percent and that of trnH-psbA was 25.6 percent to 67.7 percent, depending upon the data set and the method used. matK provided the least PCR success, followed by ITS2 (59. 0%). Species resolution by ITS2 and rbcL ranged from 9.0 percent to 48.7 percent and 13.2 percent to 43.6 percent, respectively. Further, we observed that the NCBI nucleotide database is poorly represented by the sequences of barcode loci studied here for tree species. CONCLUSION: Although a conservative approach of a success rate of 60-70 percent by both ITS and trnH-psbA may not be considered as highly successful but would certainly help in large-scale biodiversity inventorization, particularly for tropical tree species, considering the standard success rate of plant DNA barcode program reported so far. The recommended matK and rbcL primers combination may not work in tropical tree species as barcode markers.

Universal plant DNA barcode loci may not work in complex groups: a case study with Indian berberis species.

BACKGROUND: The concept of DNA barcoding for species identification has gained considerable momentum in animals because of fairly successful species identification using cytochrome oxidase I (COI). In plants, matK and rbcL have been proposed as standard barcodes. However, barcoding in complex genera is a challenging task. METHODOLOGY AND PRINCIPAL FINDINGS: We investigated the species discriminatory power of four reportedly most promising plant DNA barcoding loci (one from nuclear genome--ITS, and three from plastid genome--trnH-psbA, rbcL and matK) in species of Indian Berberis L. (Berberidaceae) and two other genera, Ficus L. (Moraceae) and Gossypium L. (Malvaceae). Berberis species were delineated using morphological characters. These characters resulted in a well resolved species tree. Applying both nucleotide distance and nucleotide character-based approaches, we found that none of the loci, either singly or in combinations, could discriminate the species of Berberis. ITS resolved all the tested species of Ficus and Gossypium and trnH-psbA resolved 82% of the tested species in Ficus. The highly regarded matK and rbcL could not resolve all the species. Finally, we employed amplified fragment length polymorphism test in species of Berberis to determine their relationships. Using ten primer pair combinations in AFLP, the data demonstrated incomplete species resolution. Further, AFLP analysis showed that there was a tendency of the Berberis accessions to cluster according to their geographic origin rather than species affiliation. CONCLUSIONS/SIGNIFICANCE: We reconfirm the earlier reports that the concept of universal barcode in plants may not work in a number of genera. Our results also suggest that the matK and rbcL, recommended as universal barcode loci for plants, may not work in all the genera of land plants. Morphological, geographical and molecular data analyses of Indian species of Berberis suggest probable reticulate evolution and thus barcode markers may not work in this case.

Oligonucleotide frequencies of barcoding loci can discriminate species across kingdoms.

BACKGROUND: DNA barcoding refers to the use of short DNA sequences for rapid identification of species. Genetic distance or character attributes of a particular barcode locus discriminate the species. We report an efficient approach to analyze short sequence data for discrimination between species. METHODOLOGY AND PRINCIPAL FINDINGS: A new approach, Oligonucleotide Frequency Range (OFR) of barcode loci for species discrimination is proposed. OFR of the loci that discriminates between species was characteristic of a species, i.e., the maxima and minima within a species did not overlap with that of other species. We compared the species resolution ability of different barcode loci using p-distance, Euclidean distance of oligonucleotide frequencies, nucleotide-character based approach and OFR method. The species resolution by OFR was either higher or comparable to the other methods. A short fragment of 126 bp of internal transcribed spacer region in ribosomal RNA gene was sufficient to discriminate a majority of the species using OFR. CONCLUSIONS/SIGNIFICANCE: Oligonucleotide frequency range of a barcode locus can discriminate between species. Ability to discriminate species using very short DNA fragments may have wider applications in forensic and conservation studies.

Characterization of a Taura syndrome virus isolate originating from the 2004 Texas epizootic in cultured shrimp.

A comprehensive investigation of the Taura syndrome virus (TSV) isolate that caused epizootics in shrimp farms in Texas in 2004 (Texas isolate) revealed that this virus was more virulent in laboratory bioassays than the TSV reference isolate, Hawaii 1994, causing severe symptom development and rapid mortality. Histopathology of moribund animals demonstrated epithelial necrosis within the stomach, appendages, general body cuticle and gills, and the surviving animals demonstrated moderate to numerous lymphoid organ spheroids. Purified virions showed icosahedral morphology, with a diameter of 31 nm. Comparative genome analysis showed that the Texas isolate is more closely related to TSV isolates from Thailand and China than to the Hawaii isolate. The predicted tertiary structures of the inhibition of apoptosis protein (IAP) and protease domains of the Texas isolate are very similar to those of the Hawaii isolate. However, the RNA-dependent RNA polymerase (RdRp) of the Texas isolate has significant structural differences from the Hawaii isolate due to point mutation(s) in the RdRp gene. Changes in the RdRp tertiary structure might contribute to the replication fidelity, virulence and ecological adaptability of the Texas isolate.

Rabies glycoprotein fused with B subunit of cholera toxin expressed in tobacco plants folds into biologically active pentameric protein.

The pentameric B subunit of cholera toxin (CtxB) is an efficient mucosal adjuvant for vaccines. We report the expression of a chimeric protein comprising the synthetic cholera toxin B subunit fused at its C-terminal with rabies surface glycoprotein (G protein) in tobacco plants. The approximately 80.3 kDa fusion polypeptide expressed at 0.4% of the total soluble protein in leaves of the selected transgenic lines. The fusion protein formed a approximately 403 kDa pentameric protein which was functionally active in binding to GM1 receptor. The plant-made protein had a higher affinity for GM1 receptor than the native bacterial CtxB. The pentameric fusion protein was recognized by the anti-cholera toxin as well as anti-rabies antibodies. Its immuno-protective ability against rabies remains to be examined.

High level expression of a functionally active cholera toxin B: rabies glycoprotein fusion protein in tobacco seeds.

A synthetic DNA construct containing cholera toxin B subunit, genetically fused to the surface glycoprotein of rabies virus was expressed in tobacco plants from a seed specific (legumin) promoter. Seed specific expression was monitored by real-time PCR, GM1-ELISA and Western blot analyses. The fusion protein accumulated in tobacco seeds at up to 1.22% of the total seed protein. It was functionally active in binding to the GM1-ganglioside receptors, suggesting its assembly into pentamers in seeds of the transgenic plants. Immunoblot analysis confirmed that the approximately 80.6 kDa monomeric fusion polypeptide was expressed in tobacco seeds and accumulated as an approximately 403 kDa pentamer. Evaluation of its immunoprotective ability against rabies and cholera is to be examined.

Molecular cloning and functional identification of a ribosome inactivating/antiviral protein from leaves of post-flowering stage of Celosia cristata and its expression in E. coli.

A full-length cDNA clone, encoding a ribosome inactivating/antiviral protein (RIP/AVP) was isolated from the cDNA library of post-flowering stage of Celosia cristata leaves. The full-length cDNA consisted of 1015 nucleotides, with an open reading frame encoding 283 amino acids. The deduced amino acid sequence had a putative active site domain conserved in other ribosome inactivating/antiviral proteins (RIPs/AVPs). The coding region of the cDNA was amplified by polymerase chain reaction (PCR), cloned and expressed in Escherichia coli as recombinant protein of 72 kDa. The expressed fusion product was confirmed by Western analysis and purification by affinity chromatography. Both the recombinant protein (reCCP-27) and purified expressed protein (eCCP-27) inhibited translation in rabbit reticulocytes showing IC50 values at 95 ng and 45 ng, respectively. The native purified nCCP-27 has IC50 at 25 ng. The purified product also showed N-glycosidase activity towards tobacco ribosomes and antiviral activity towards tobacco mosaic virus (TMV) and sunnhemp rosette virus (SRV).

Purification, characterization and cloning of antiviral/ribosome inactivating protein from Amaranthus tricolor leaves.

An antiviral protein (AVP), imparting high level of resistance against sunnhemp rosette virus (SRV) was purified from the dried leaves of Amaranthus tricolor. The purified protein (AAP-27) exhibited approximately 98% inhibition of local lesion formation at a concentration range of approximately 30 microg ml(-1). The protein was found to be highly basic glycoprotein monomer (pI approximately 9.8) of Mr 27 kDa, with neutral sugar content of 4%. The purified protein exhibited N-glycosidase and RNase activities. We have also isolated full-length cDNA clone, encoding this protein designated as A. tricolor antiviral protein-1 (AAP-1). Two primers, one designed on the basis of N-terminal sequence of the purified protein and the other from the conserved active peptides of other AVPs/RIPs were used for PCR amplification of double stranded cDNA, isolated from the leaves of A. tricolor. The amplified fragment was used as a probe for library screening. The isolated full-length cDNA consisted of 1058 nucleotides with an open reading frame encoding a polypeptide of 297 amino acids. The deduced amino acid sequence of AAP-1 has a putative active domain conserved in other AVPs/RIPs and shows varying homology to the RIPs from other plant species.