Understanding the evolution of miRNA biogenesis machinery in plants with special focus on rice.

miRNA biogenesis process is an intricate and complex event consisting of many proteins working in a highly coordinated fashion. Most of these proteins have been studied in Arabidopsis; however, their orthologs and functions have not been explored in other plant species. In the present study, we have manually curated all the experimentally verified information present in the literature regarding these proteins and found a total of 98 genes involved in miRNA biogenesis in Arabidopsis. The conservation pattern of these proteins was identified in other plant species ranging from dicots to lower organisms, and we found that a major proportion of proteins involved in the pri-miRNA processing are conserved. However, nearly 20% of the genes, mostly involved in either transcription or functioning of the miRNAs, were absent in the lower organisms. Further, we manually curated a regulatory network of the core components of the biogenesis process and found that nearly half (46%) of the proteins interact with them, indicating that the processing step is perhaps the most under surveillance/regulation. We have subsequently attempted to characterize the orthologs identified in Oryza sativa, on the basis of transcriptome and epigenetic modifications under field drought conditions in order to assess the impact of drought on the process. We found several participating genes to be differentially expressed and/or epigenetically methylated under drought, although the core components like DCL1, SE, and HYL1 remain unaffected by the stress itself. The study enhances our present understanding of the biogenesis process and its regulation.

Delineating the tissue-mediated drought stress governed tuning of conserved miR408 and its targets in rice.

Engineering drought tolerance in rice needs to focus on regulators that enhance tolerance while boosting plant growth and vigor. The present study delineated the concealed function and tissue-mediated interplay of the miR408/target module in imparting drought stress tolerance in rice. The plant miR408 family comprises three dominant mature forms (21 nt), including a distinct monocot variant (F-7 with 5' C) and is divided into six groups. miR408 majorly cleaves genes belonging to the blue copper protein in addition to several other species-specific targets in plants. Comparative sequence analysis in 4726 rice accessions identified 22 sequence variants (SNP and InDELs) in its promoter (15) and pre-miR408 region. Haplotype analysis of the sequence variants indicated eight haplotypes (three: Japonica-specific and five: Indica-specific) of the miR408 promoter. In drought-tolerant Nagina 22, miR408 follows flag leaf preferential expression. Under drought conditions, its levels are upregulated in flag leaf and roots which seems to be regulated by a differential fraction of methylated cytosines (mCs) in the precursor region. The active pool of miR408 regulated targets under control and drought conditions is impacted by the tissue type. Comparative expression analysis of the miR408/target module under different sets of conditions features 83 targets exhibiting antagonistic expression in rice, out of which 12 genes, including four PLANTACYANINS (OsUCL6, 7, 9 and 30), PIRIN, OsLPR1, OsCHUP1, OsDOF12, OsBGLU1, glycine-rich cell wall gene, OsDUT, and OsERF7, are among the high confidence targets. Further, overexpression of MIR408 in drought-sensitive rice cultivar (PB1) leads to the massive enhancement of vegetative growth in rice with improved ETR and Y(II) and enhanced dehydration stress tolerance. The above results suggest that miR408 is likely to act as a positive regulator of growth and vigor, as well as dehydration stress, making it a potential candidate for engineering drought tolerance in rice.

Variety-specific transcript accumulation during reproductive stage in drought-stressed rice.

The divergence of natural stress tolerance mechanisms between species is an intriguing phenomenon. To study it in rice, a comparative transcriptome analysis was carried out in 'heading' stage tissue (flag leaf, panicles and roots) of Nagina 22 (N22; drought-tolerant) and IR64 (drought-sensitive) plants subjected to field drought. Interestingly, N22 showed almost double the number of differentially expressed genes (DEGs) than IR64. Many DEGs colocalized within drought-related QTLs responsible for grain yield and drought tolerance and also associated with drought tolerance and critical drought-related plant traits such as leaf rolling, trehalose content, sucrose and cellulose content. Besides, co-expression analysis of the DEGs revealed several 'hub' genes known to actively regulate drought stress response. Strikingly, 1366 DEGs, including 21 'hub' genes, showed a distinct opposite regulation in the two rice varieties under similar drought conditions. Annotation of these variety-specific DEGs (VS-DEGs) revealed that they are distributed in various biological pathways. Furthermore, 103 VS-DEGs were found to physically interact with over 1300 genes, including 32 that physically interact with other VS-DEGs as well. The promoter region of these genes has sequence variations among the two rice varieties, which might be in part responsible for their unique expression pattern.

Emerging Molecular Strategies for Improving Rice Drought Tolerance.

Rice occupies a pre-eminent position as a food crop in the world. Its production, how- ever, entails up to 3000 liters of water per kilogram of grain produced. Such high demand makes rice prone to drought easily. Sustainable rice cultivation with limited water resources requires the deployment of a suitable strategy for better water use efficiency and improved drought tolerance. Several drought-related genes have been evaluated in rice for their mode of action in conferring drought tolerance. Manipulation of components of abscisic acid signal transduction, stomatal density, deposition of cuticular wax, and protein modification pathways are emerging as priority targets. Gene reprogramming by microRNAs is also being explored to achieve drought tolerance. Genetically dissected Quantitative Trait Loci (QTLs) and their constituent genes are being deployed to develop drought-tolerant rice varieties. Progressive research and challenges include a better understanding of crucial components of drought response and search for new targets and the deployment of improved varieties in the field.

Genome-wide analysis of H3K4me3 and H3K27me3 modifications due to Lr28 for leaf rust resistance in bread wheat (Triticum aestivum).

KEY MESSAGE: Present study revealed a complex relationship among histone H3 methylation (examined using H3K4/K27me3 marks), cytosine DNA methylation and differential gene expression during Lr28 mediated leaf rust resistance in wheat. During the present study, genome-wide histone modifications were examined in a pair of near isogenic lines (NILs) (with and without Lr28 in the background of cv. HD2329). The two histone marks used included H3K4me3 (an activation mark) and H3K27me3 (a repression mark). The results were compared with levels of expression (using RNA-seq) and DNA methylation (MeDIP) data obtained using the same pair of NILs. Some of the salient features of the present study include the following: (i) large scale differential binding sites (DBS) were available for only H3K4me3 in the susceptible cultivar, but for both H3K4me3 and H3K27me3 in its resistant NIL; (ii) DBSs for H3K27me3 mark were more abundant (> 80%) in intergenic regions, whereas DBSs for H3K4me3 were distributed in all genomic regions including exons, introns, intergenic, TTS (transcription termination sites) and promoters; (iii) fourteen (14) genes associated with DBSs showed co-localization for both the marks; (iv) only a small fraction (7% for H3K4me3 and 12% for H3K27me3) of genes associated with DBSs matched with the levels of gene expression inferred from RNA-seq data; (v) validation studies using qRT-PCR were conducted on 26 selected representative genes; results for only 11 genes could be validated. The proteins encoded by important genes involved in promoting infection included domains generally carried by R gene proteins such as Mlo like protein, protein kinases and purple acid phosphatase. Similarly, proteins encoded by genes involved in resistance included those carrying domains for lectin kinase, R gene, aspartyl protease, etc. Overall, the results suggest a very complex network of downstream genes that are expressed during compatible and incompatible interactions; some of the genes identified during the present study may be used in future validation studies involving RNAi/overexpression approaches.

Investigation into the miRNA/5' isomiRNAs function and drought-mediated miRNA processing in rice.

MicroRNAs lie at the core of biological regulatory networks in plants. The recent discovery of isomiRs that are length variants of the annotated mature miRNAs has further unveiled the complexity of miRNome. Delineation of their functional relevance is critical to understand the complete functional spectrum of the miRNome. To apprehend the role of 5' isomiRs in rice, we performed a comprehensive analysis of the annotated miRNA pool using 8 deep-sequencing datasets from flag leaf and spikelet tissues from two cultivars of rice viz. N22 and PB1 grown under control and drought conditions. The products of the 5' start site variability termed as "5' isomiRs" were found to be widespread in all the datasets. It was possible to identify several 5' isomiRs that were highly distinct and abundant and supported by more than 90% of the tags that map in the region. Majority of miRNA/5' isomiR pair share similar tissue and drought-mediated expression dynamics. Analysis of the degradome data identified targets for several of these 5' isomiRs, thereby confirming their biological activity. Since the isomiRs are length variants at the 5' end, the target sites were found to be accordingly shifted as compared to the target site of the annotated miRNA. Further we also observed that drought affects the processing accuracy of several miRNAs across all tissues of both the cultivars leading to differential accumulation of 5' isomiR/miRNA pair.

H3K4/K9 acetylation and Lr28-mediated expression of six leaf rust responsive genes in wheat (Triticum aestivum).

Development of leaf rust-resistant cultivars is a priority during wheat breeding, since leaf rust causes major losses in yield. Resistance against leaf rust due to Lr genes is partly controlled by epigenetic modifications including histone acetylation that is known to respond to biotic/abiotic stresses. In the present study, enrichment of H3K4ac and H3K9ac in promoters of six defense responsive genes (N-acetyltransferase, WRKY 40, WRKY 70, ASR1, Peroxidase 12 and Sarcosine oxidase) was compared with their expression in a pair of near-isogenic lines (NILs) for the gene Lr28 following inoculation with leaf rust pathotype '77-5'; ChIP-qPCR was used for this purpose. The proximal and distal promoters of these genes contained a number of motifs that are known to respond to biotic stresses. The enrichment of two acetylation marks changed with passage of time; changes in expression of two of the six genes (N-acetyltransferase and peroxidase12), largely matched with changes in H3K4/H3K9 acetylation patterns of the two promoter regions. For example, enrichment of both the marks matched with higher expression of N-acetyltransferase gene in susceptible NIL and the deacetylation (H3K4ac) largely matched with reduced gene expression in resistant NIL. In peroxidase12, enrichment of H3K4ac and H3K9ac largely matched with higher expression in both the NILs. In the remaining four genes, changes in H3 acetylation did not always match with gene expression levels. This indicated complexity in the regulation of the expression of these remaining four genes, which may be controlled by other epigenetic/genetic regulatory mechanisms that need further analysis.

Evaluating five different loci (rbcL, rpoB, rpoC1, matK, and ITS) for DNA barcoding of Indian orchids.

Orchidaceae, one of the largest families of angiosperms, is represented in India by 1600 species distributed in diverse habitats. Orchids are in high demand owing to their beautiful flowers and therapeutic properties. Overexploitation and habitat destruction have made many orchid species endangered. In the absence of effective identification methods, illicit trade of orchids continues unabated. Considering DNA barcoding as a potential identification tool, species discrimination capability of five loci, ITS, matK, rbcL, rpoB, and rpoC1, was tested in 393 accessions of 94 Indian orchid species belonging to 47 genera, including one listed in Appendix I of CITES and 26 medicinal species. ITS provided the highest species discrimination rate of 94.9%. While, among the chloroplast loci, matK provided the highest species discrimination rate of 85.7%. None of the tested loci individually discriminated 100% of the species. Therefore, multi-locus combinations of up to five loci were tested for their species resolution capability. Among two-locus combinations, the maximum species resolution (86.7%) was provided by ITS+matK. ITS and matK sequences of the medicinal orchids were species specific, thus providing unique molecular identification tags for their identification and detection. These observations emphasize the need for the inclusion of ITS in the core barcode for plants, whenever required and available.

Manually curated database of rice proteins.

'Manually Curated Database of Rice Proteins' (MCDRP) available at http://www.genomeindia.org/biocuration is a unique curated database based on published experimental data. Semantic integration of scientific data is essential to gain a higher level of understanding of biological systems. Since the majority of scientific data is available as published literature, text mining is an essential step before the data can be integrated and made available for computer-based search in various databases. However, text mining is a tedious exercise and thus, there is a large gap in the data available in curated databases and published literature. Moreover, data in an experiment can be perceived from several perspectives, which may not reflect in the text-based curation. In order to address such issues, we have demonstrated the feasibility of digitizing the experimental data itself by creating a database on rice proteins based on in-house developed data curation models. Using these models data of individual experiments have been digitized with the help of universal ontologies. Currently, the database has data for over 1800 rice proteins curated from >4000 different experiments of over 400 research articles. Since every aspect of the experiment such as gene name, plant type, tissue and developmental stage has been digitized, experimental data can be rapidly accessed and integrated.

Unique miRNome during anthesis in drought-tolerant indica rice var. Nagina 22.

MAIN CONCLUSION: Drought-tolerant rice variety, Nagina 22 (N22), has a unique spikelet miRNome during anthesis stage drought as well as transition from heading to anthesis. Molecular characterization of genetic diversity of rice is essential to understand the evolution and molecular basis of various agronomically important traits such as drought tolerance. miRNAs play an important role in regulating plant development as well as stress response such as drought. In this study, we characterized the yet unexplored dynamics of the spikelet miRNA population during developmental transition from 'heading' to 'anthesis' as well as anthesis stage drought stress in a drought-tolerant indica rice variety, N22. A significant proportion of miRNA population (~20 %) in N22 spikelets is modulated during transition from heading to anthesis indicating a unique miRNome at anthesis, a developmental stage highly sensitive to stress (drought/heat). Based on the analysis of degradome data, majority of differentially regulated miRNAs appear to regulate transcription factors, some of which are implicated in regulation of development and fertilization. Similarly, drought during anthesis leads to a global change in miRNA expression pattern including those which regulate ROS homeostasis. It was possible to identify several miRNAs that were not reported to be drought responsive in earlier studies. Interestingly, a significant proportion of the drought-regulated miRNAs co-localize within QTLs related to drought tolerance and associated traits. Comparison of the expression profiles between N22 and Pusa Basmati 1 (drought sensitive) identified miRNAs with variety-specific expression patterns during phase transition (miR164, miR396, miR812, and miR1881) as well as drought stress (miR1881) indicating an evolution of a distinct and variety-specific regulatory mechanism. The promoters of these miRNAs contain LREs (light-responsive elements) and are induced by dark treatment. It was also possible to identify 4 novel miRNAs including an intronic miRNA that was conserved in both rice varieties.

Massive gene acquisitions in Mycobacterium indicus pranii provide a perspective on mycobacterial evolution.

Understanding the evolutionary and genomic mechanisms responsible for turning the soil-derived saprophytic mycobacteria into lethal intracellular pathogens is a critical step towards the development of strategies for the control of mycobacterial diseases. In this context, Mycobacterium indicus pranii (MIP) is of specific interest because of its unique immunological and evolutionary significance. Evolutionarily, it is the progenitor of opportunistic pathogens belonging to M. avium complex and is endowed with features that place it between saprophytic and pathogenic species. Herein, we have sequenced the complete MIP genome to understand its unique life style, basis of immunomodulation and habitat diversification in mycobacteria. As a case of massive gene acquisitions, 50.5% of MIP open reading frames (ORFs) are laterally acquired. We show, for the first time for Mycobacterium, that MIP genome has mosaic architecture. These gene acquisitions have led to the enrichment of selected gene families critical to MIP physiology. Comparative genomic analysis indicates a higher antigenic potential of MIP imparting it a unique ability for immunomodulation. Besides, it also suggests an important role of genomic fluidity in habitat diversification within mycobacteria and provides a unique view of evolutionary divergence and putative bottlenecks that might have eventually led to intracellular survival and pathogenic attributes in mycobacteria.

Understanding the relationship between cerebellum and the frontal-cortex region of C9orf72-related amyotrophic lateral sclerosis: A comparative analysis of genetic features.

BACKGROUND: Amyotrophic lateral sclerosis (ALS) is a relentlessly progressive and fatal neurodegenerative diseases for which at present no cure is available. Despite the extensive research the progress from diagnosis to prognosis in ALS and frontotemporal dementia (FTD) has been slow which represents suboptimal understanding of disease pathophysiological processes. In recent studies, several genes have been associated with the ALS and FTD diseases such as SOD1, TDP43, and TBK1, whereas the hexanucleotide GGGGCC repeat expansion (HRE) in C9orf72 gene is a most frequent cause of ALS and FTD, that has changed the understanding of these diseases. METHODS: The goal of this study was to identify and spatially determine differential gene expression signature differences between cerebellum and frontal cortex in C9orf72-associated ALS (C9-ALS), to study the network properties of these differentially expressed genes, and to identify miRNAs targeting the common differentially expressed genes in both the tissues. This study thus highlights underlying differential cell susceptibilities to the disease mechanisms in C9-ALS and suggesting therapeutic target selection in C9-ALS. RESULTS: In this manuscript, we have identified that the genes involved in neuron development, protein localization and transcription are mostly enriched in cerebellum of C9-ALS patients, while the UPR-related genes are enriched in the frontal cortex. Of note, UPR pathway genes were mostly dysregulated both in the C9-ALS cerebellum and frontal cortex. Overall, the data presented here show that defects in normal RNA processing and the UPR pathway are the pathological hallmarks of C9-ALS. Interestingly, the cerebellum showed more strong transcriptome changes than the frontal cortex. CONCLUSION: Interestingly, the cerebellum region showed more significant transcriptomic changes as compared to the frontal cortex region suggesting its active participation in the disease process. This nuanced understanding may offer valuable insights for the development of targeted therapeutic strategies aimed at mitigating disease progression in C9-ALS.

The SARS-CoV-2 targeted human RNA binding proteins network biology to investigate COVID-19 associated manifestations.

The global coronavirus disease 2019 (COVID-19) pandemic caused by the SARS-CoV-2 virus has had unprecedented social and economic ramifications. Identifying targets for drug repurposing could be an effective means to present new and fast treatments. Furthermore, the risk of morbidity and mortality from COVID-19 goes up when there are coexisting medical conditions, however, the underlying mechanisms remain unclear. In the current study, we have adopted a network-based systems biology approach to investigate the RNA binding proteins (RBPs)-based molecular interplay between COVID-19, various human cancers, and neurological disorders. The network based on RBPs commonly involved in the three disease conditions consisted of nine RBPs connecting 10 different cancer types, 22 brain disorders, and COVID-19 infection, ultimately hinting at the comorbidities and complexity of COVID-19. Further, we underscored five miRNAs with reported antiviral properties that target all of the nine shared RBPs and are thus therapeutically valuable. As a strategy to improve the clinical conditions in comorbidities associated with COVID-19, we propose perturbing the shared RBPs by drug repurposing. The network-based analysis presented hereby contributes to a better knowledge of the molecular underpinnings of the comorbidities associated with COVID-19.

Multi-stage transcriptome profiling of the neglected food-borne echinostome Artyfechinostomum sufrartyfex reveal potential diagnostic and drug targets.

Lack of effective surveillance and control methods for neglected helminth diseases particularly in context of rural areas in India is a serious concern in terms of public health. With regard to the emerging food-borne echinostomid Artyfechinostomum sufrartyfex infection in the country, the current study is an in silico attempt to screen for plausible diagnostic and drug targets against the trematode. Transcriptome of adult, encysted and excysted metacercaria stages of the parasite was generated using Illumina sequencing platform. A de-novo assembly strategy utilizing transcriptome data generated from the three lifecycle stages was followed to generate the representative transcripts. Longest open reading frames identified for the transcripts were further conceptually translated into their respective protein sequences. Detailed analysis of this dataset through various bioinformatics pipelines and tools eventually identified 14 credible diagnostic and 10 drug targets along with their FDA-approved and ZINC molecules. Some of the important diagnostic candidates include thioredoxin peroxidase, haemoglobinase, cathepsin L, cathepsin L-like and B-like cysteine proteases. Among the drug targets, uncharacterized sodium dependent transporter and bifunctional protein Aas were identified as top targets exhibiting significant interaction with Rifamycin and ZINC02820058 molecule, respectively. Further, B-cell epitope analysis of the diagnostic targets revealed unique epitopes for 10 of them thus indicating their potential role in specific diagnosis of the parasite. The diagnostic candidates along with a number of lesser known drug targets and their ligand molecules identified in this study provides a reasonable basis for evaluation and development of future intervention strategies against A. sufrartyfex.

Whole genome sequence of the rifamycin B-producing strain Amycolatopsis mediterranei S699.

Amycolatopsis mediterranei S699 is an actinomycete that produces an important antibiotic, rifamycin B. Semisynthetic derivatives of rifamycin B are used for the treatment of tuberculosis, leprosy, and AIDS-related mycobacterial infections. Here, we report the complete genome sequence (10.2 Mb) of A. mediterranei S699, with 9,575 predicted coding sequences.

Investigations on Regulation of MicroRNAs in Rice Reveal [Ca2+]cyt Signal Transduction Regulated MicroRNAs.

MicroRNAs (miRNAs) are critical components of the multidimensional regulatory networks in eukaryotic systems. Given their diverse spectrum of function, it is apparent that the transcription, processing, and activity of the miRNAs themselves, is very dynamically regulated. One of the most important and universally implicated signaling molecules is [Ca2+]cyt. It is known to regulate a plethora of developmental and metabolic processes in both plants and animals; however, its impact on the regulation of miRNA expression is relatively less explored. The current study employed a combination of internal and external calcium channel inhibitors to establishing that [Ca2+]cyt signatures actively regulate miRNA expression in rice. Involvement of [Ca2+]cyt in the regulation of miRNA expression was further confirmed by treatment with calcimycin, the calcium ionophore. Modulation of the cytosolic calcium levels was also found to regulate the drought-responsive expression as well as ABA-mediated response of miRNA genes in rice seedlings. The study further establishes the role of calmodulins and Calmodulin-binding Transcription Activators (CAMTAs) as important components of the signal transduction schema that regulates miRNA expression. Yeast one-hybrid assay established that OsCAMTA4 & 6 are involved in the transcriptional regulation of miR156a and miR167h. Thus, the study was able to establish that [Ca2+]cyt is actively involved in regulating the expression of miRNA genes both under control and stress conditions.

The water-deficit stress- and red-rot-related genes in sugarcane.

Sugarcane is an important international commodity as a valuable agricultural crop especially in developing countries. Sequencing was carried out to generate >35,000 expressed sequence tags (ESTs) from healthy as well as red-rot-infected tissue of Indian subtropical variety of sugarcane. Subsequent clustering with existing sugarcane ESTs in public databases identified 4,087 clusters, including 85 clusters that preferentially express upon Colletotrichum falcatum (red-rot) infection, which were previously unreported. Real-time reverse transcription-PCR profiling of selected EST clusters identified several sugarcane clusters that show differential expression in response to biotic and abiotic stress conditions. Twenty-five stress-related clusters showed >2-fold relative expression during water-deficit stress in sugarcane. Similarly, EST clusters could be identified, which exhibit association with red-rot disease when assessed in red-rot-susceptible and red-rot-resistant varieties of sugarcane. Such EST clusters are good candidates for in-depth analysis to elucidate stress-responsive pathways in sugarcane and facilitate genetic manipulation to tailor this crop for tolerance to various stresses.

Development of EST-based new SSR markers in seabuckthorn.

EST-based SSR markers were developed by screening a collection of 1584 clustered ESTs of seabuckthorn (Hippophae rhamnoides). PCR primers were designed for the amplification of 30 microsatellite loci. Two to five allelic bands were displayed by nine primer pairs in H. rhamnoides genotypes and by eleven primer pairs in H. salicifolia genotypes. None of the thirty primer pairs detected polymorphism in H. tibetana genotypes. Considering the high polymorphism detected in the tested genotypes and their direct origin from the genic regions, these EST-SSR markers hold immense promise in seabuckthorn genome analysis, molecular breeding and population genetics.

A study of transcriptome in leaf rust infected bread wheat involving seedling resistance gene Lr28.

Leaf rust disease causes severe yield losses in wheat throughout the world. During the present study, high-throughput RNA-Seq analysis was used to gain insights into the role of Lr28 gene in imparting seedling leaf rust resistance in wheat. Differential expression analysis was conducted using a pair of near-isogenic lines (NILs) (HD 2329 and HD 2329+Lr28) at early (0h before inoculation (hbi), 24 and 48h after inoculation (hai)) and late stages (72, 96 and 168 hai) after inoculation with a virulent pathotype of pathogen Puccinia triticina. Expression of a large number of genes was found to be affected due to the presence/absence of Lr28. Gene ontology analysis of the differentially expressed transcripts suggested enrichment of transcripts involved in carbohydrate and amino acid metabolism, oxidative stress and hormone metabolism, in resistant and/or susceptible NILs. Genes encoding receptor like kinases (RLKs) (including ATP binding; serine threonine kinases) and other kinases were the most abundant class of genes, whose expression was affected. Genes involved in reactive oxygen species (ROS) homeostasis and several genes encoding transcription factors (TFs) (most abundant being WRKY TFs) were also identified along with some ncRNAs and histone variants. Quantitative real-time PCR was also used for validation of 39 representative selected genes. In the long term, the present study should prove useful in developing leaf rust resistant wheat cultivars through molecular breeding.