Unveiling the Micronome of Cassava (Manihot esculenta Crantz)

In silico identification of papaya genome-encoded microRNAs to target begomovirus genes in papaya leaf curl disease

Papaya leaf curl disease (PaLCuD) is widespread and classified in the genus begomovirus (Geminiviridae), disseminated by the vector whitefly Bemisia tabaci. RNA interference (RNAi)-based antiviral innate immunity stands as a pivotal defense mechanism and biological process in limiting viral genomes to manage plant diseases. The current study aims to identify and analyze Carica Papaya locus-derived capa-microRNAs with predicted potential for targeting divergent begomovirus species-encoded mRNAs using a 'four integrative in silico algorithms' approach. This research aims to experimentally activate the RNAi catalytic pathway using in silico-predicted endogenous capa-miRNAs and create papaya varieties capable of assessing potential resistance against begomovirus species and monitoring antiviral capabilities. This study identified 48 predicted papaya locus-derived candidates from 23 miRNA families, which were further investigated for targeting begomovirus genes. Premised all the four algorithms combined, capa-miR5021 was the most anticipated miRNA followed by capa-miR482, capa-miR5658, capa-miR530b, capa-miR3441.2, and capa-miR414 'effective' papaya locus-derived candidate capa-miRNA and respected putative binding sites for targets at the consensus nucleotide position. It was predicted to bind and target mostly to AC1 gene of the complementary strand and the AV1 gene of the virion strand of different begomovirus isolates, which were associated with replication-associated protein and encapsidation, respectively, during PaLCuD. These miRNAs were also found targeting betaC1 gene of betasatellite which were associated with retardation in leaf growth and developmental abnormalities with severe symptoms during begomovirus infection. To validate target prediction accuracy, we created an integrated Circos plot for comprehensive visualization of host-virus interaction. In silico-predicted papaya genome-wide miRNA-mediated begomovirus target gene regulatory network corroborated interactions that permit in vivo analysis, which could provide biological material and valuable evidence, leading to the development of begomovirus-resistant papaya plants. The integrative nature of our research positions it at the forefront of efforts to ensure the sustainable cultivation of papaya, particularly in the face of evolving pathogenic threats. As we move forward, the knowledge gained from this study provides a solid foundation for continued exploration and innovation in the field of papaya virology, and to the best of our knowledge, this study represents a groundbreaking endeavor, undertaken for the first time in the context of PaLCuD research.

Genome-wide identification and expression analysis of Hsp70 family genes in Cassava (Manihot esculenta Crantz)

Hsp70 proteins function as molecular chaperones, regulating various cellular processes in plants. In this study, a genome-wide analysis led to the identification of 22 Hsp70 (MeHsp70) genes in cassava. Phylogenetic relationship studies with other Malpighiales genomes (Populus trichocarpa, Ricinus communis and Salix purpurea) classified MeHsp70 proteins into eight groups (Ia, Ib, Ic, Id, Ie, If, IIa and IIb). Promoter analysis of MeHsp70 genes revealed the presence of tissue-specific, light, biotic and abiotic stress-responsive cis-regulatory elements showing their functional importance in cassava. Meta-analysis of publically available RNA-seq transcriptome datasets showed constitutive, tissue-specific, biotic and abiotic stress-specific expression patterns among MeHsp70s in cassava. Among 22 Hsp70, six MeHsp70s viz., MecHsp70-3, MecHsp70-6, MeBiP-1, MeBiP-2, MeBiP-3 and MecpHsp70-2 displayed constitutive expression, while three MecHsp70s were induced under both drought and cold stress conditions. Five MeHsp70s, MecHsp70-7, MecHsp70-11, MecHsp70-12, MecHsp70-13, and MecHsp70-14 were induced under drought stress conditions. We predicted that 19 MeHsp70 genes are under the regulation of 24 miRNAs. This comprehensive genome-wide analysis of the Hsp70 gene family in cassava provided valuable insights into their functional roles and identified various potential Hsp70 genes associated with stress tolerance and adaptation to environmental stimuli.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03760-3.

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

Genome-Wide Screening and Characterization of Non-Coding RNAs in Coffea canephora

Coffea canephora grains are highly traded commodities worldwide. Non-coding RNAs (ncRNAs) are transcriptional products involved in genome regulation, environmental responses, and plant development. There is not an extensive genome-wide analysis that uncovers the ncRNA portion of the C. canephora genome. This study aimed to provide a curated characterization of six ncRNA classes in the Coffea canephora genome. For this purpose, we employed a combination of similarity-based and structural-based computational approaches with stringent curation. Candidate ncRNA loci had expression evidence analyzed using sRNA-seq libraries. We identified 7455 ncRNA loci (6976 with transcriptional evidence) in the C. canephora genome. This comprised of total 115 snRNAs, 1031 snoRNAs, 92 miRNA precursors, 602 tRNAs, 72 rRNAs, and 5064 lncRNAs. For miRNAs, we identified 159 putative high-confidence targets. This study was the most extensive genomic catalog of curated ncRNAs in the Coffea genus. This data might help elaborating more robust hypotheses in future comparative genomic studies as well as gene regulation and genome dynamics, helping to understand the molecular basis of domestication, environmental adaptation, resistance to pests and diseases, and coffee productivity.

RNA silencing of South African cassava mosaic virus in transgenic cassava expressing AC1/AC4 hp- RNA induces tolerance

•

A hairpin-RNA construct targeting South African cassava mosaic virus AC1/AC4 overlapping genome region confers tolerance in cassava.

•

Viral load and symptom severity were significantly lower in transgenic lines transformed with AC1/AC4 construct compared to wildtype cv.60444.

•

SACMV transgenic lines infected with SACMV displayed a recovery phenotype similar to cassava mosaic disease-tolerant cassava landrace TME3.

Cassava mosaic disease (CMD), caused by geminiviruses, is a major hurdle to cassava production. Due to the heterozygous nature of cassava, breeding for virus resistance is difficult, but cassava has been shown to be a good candidate for genetic engineering using RNA interference (RNAi). T This study reports on the ability of a transgene-derived RNA hairpin, homologous to an overlapping region of the SACMV replication associated protein and putative virus suppressor of silencing protein (AC1/AC4), to confer tolerance in the CMD-susceptible model cassava cultivar 60444. Three of the fourteen transgenic lines expressing SACMV AC1/AC4 hairpin-derived siRNAs showed decreased symptoms and viral loads compared to untransformed control plants. Expression of SACMV AC1/AC4 homologous siRNAs showed that this tolerance is most likely associated with post-transcriptional gene silencing of the virus. This is the first report of targeting the overlapping AC1 and AC4 genes of SACMV conferring CMD tolerance in cassava.

Genome-Wide Identification of Putative MicroRNAs in Cassava (Manihot esculenta Crantz) and Their Functional Landscape in Cellular Regulation

MicroRNAs are small noncoding RNAs, involved in the regulation of many cellular processes in plants. Hundreds of miRNAs have been identified in cassava by various techniques, yet these identifications were constrained by a lack of miRNA templates and the narrow range of conditions in transcriptome study. In this research, we conducted genome-wide analysis identification, whereby miRNAs from cassava genome were thoroughly screened using bioinformatics approach independent of predefined templates and studied conditions. Our work provided a catalog of putative mature miRNAs and explored the landscape of miRNAome in cassava. These putative miRNAs were validated using statistical analysis as well as available cassava expression data. We showed that the crowded locations of cassava miRNAs are consistent with other plants and animals and hypothesized to have the same evolutionary origin. At least 10 conserved miRNAs were identified in cassava based on the comparative study of miRNA conservation. Finally, investigation of miRNAs and target gene relationships enabled us to envisage the complexities of cellular regulatory systems modulated at posttranscriptional level.

Identifying High Confidence microRNAs in the Developing Seeds of Jatropha curcas

MicroRNAs (miRNAs) are endogenously short noncoding regulatory RNAs implicated in plant development and physiology. Nine small RNA (sRNA) libraries from three typical seed developmental stages (young, intermediate, and mature) were generated by deep sequencing to identify the miRNAs of J. curcas, a potential oilseed crop for the production of renewable oil. Strict criteria were adopted to identify 93 high confidence miRNAs including 48 conserved miRNAs and 45 novel miRNAs. Target genes of these miRNAs were involved in a broad range of physiological functions, including gene expression regulation, primary & secondary metabolism, growth & development, signal transduction, and stress response. About one third (29 out of 93) miRNAs showed significant changes in expression levels during the seed developmental process, indicating that the miRNAs might regulate its targets by their changes of transcription levels in seed development. However, most miRNAs were found differentially expressed in the late stage of seed development, suggesting that miRNAs play more important roles in the stage when seed accumulating organic matters and suffering dehydration stress. This study presents the first large scale identification of high confidence miRNAs in the developing seeds of J. curcas.

Plant IsomiR Atlas: Large Scale Detection, Profiling, and Target Repertoire of IsomiRs in Plants

microRNAs (miRNAs) play an important role as key regulators controlling the post-transcriptional events in plants across development, abiotic and biotic stress, tissue polarity and also in defining the evolutionary basis of the origin of the post-transcriptional machinery. Identifying patterns of regulated and co-regulated small RNAs, in particular miRNAs and their sequence variants with the availability of next generation sequencing approaches has widely demonstrated the role of miRNAs and their temporal regulation in maintaining plant development and their response to stress conditions. Although the role of canonical miRNAs has been widely explored and functional diversity is revealed, those works for isomiRs are still limited and urgent to be carried out across plants. This relative lack of information with respect to isomiRs might be attributed to the non-availability of large-scale detection of isomiRs across wide plant species. In the present research, we addressed this by developing Plant isomiR Atlas, which provides large-scale detection of isomiRs across 23 plant species utilizing 677 smallRNAs datasets and reveals a total of 98,374 templated and non-templated isomiRs from 6,167 precursors. Plant isomiR Atlas provides several visualization features such as species specific isomiRs, isomiRs and canonical miRNAs overlap, terminal modification classifications, target identification using psRNATarget and TargetFinder and also canonical miRNAs:target interactions. Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate to understand the functional role of isomiRs. Plant isomiR Atlas is available at www.mcr.org.in/isomir. One Sentence Summary  Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate the understanding and diversity of functional targets of plants isomiRs.

MicroRNA-like RNAs from the same miRNA precursors play a role in cassava chilling responses

MicroRNAs (miRNAs) are known to play important roles in various cellular processes and stress responses. MiRNAs can be identified by analyzing reads from high-throughput deep sequencing. The reads realigned to miRNA precursors besides canonical miRNAs were initially considered as sequencing noise and ignored from further analysis. Here we reported a small-RNA species of phased and half-phased miRNA-like RNAs different from canonical miRNAs from cassava miRNA precursors detected under four distinct chilling conditions. They can form abundant multiple small RNAs arranged along precursors in a tandem and phased or half-phased fashion. Some of these miRNA-like RNAs were experimentally confirmed by re-amplification and re-sequencing, and have a similar qRT-PCR detection ratio as their cognate canonical miRNAs. The target genes of those phased and half-phased miRNA-like RNAs function in process of cell growth metabolism and play roles in protein kinase. Half-phased miR171d.3 was confirmed to have cleavage activities on its target gene P-glycoprotein 11, a broad substrate efflux pump across cellular membranes, which is thought to provide protection for tropical cassava during sharp temperature decease. Our results showed that the RNAs from miRNA precursors are miRNA-like small RNAs that are viable negative gene regulators and may have potential functions in cassava chilling responses.

Deciphering microRNAs and Their Associated Hairpin Precursors in a Non-Model Plant, Abelmoschus esculentus

MicroRNAs (miRNAs) are crucial regulatory RNAs, originated from hairpin precursors. For the past decade, researchers have been focusing extensively on miRNA profiles in various plants. However, there have been few studies on the global profiling of precursor miRNAs (pre-miRNAs), even in model plants. Here, for the first time in a non-model plant—Abelmoschus esculentus with negligible genome information—we are reporting the global profiling to characterize the miRNAs and their associated pre-miRNAs by applying a next generation sequencing approach. Preliminarily, we performed small RNA (sRNA) sequencing with five biological replicates of leaf samples to attain 207,285,863 reads; data analysis using miRPlant revealed 128 known and 845 novel miRNA candidates. With the objective of seizing their associated hairpin precursors, we accomplished pre-miRNA sequencing to attain 83,269,844 reads. The paired end reads are merged and adaptor trimmed, and the resulting 40–241 nt (nucleotide) sequences were picked out for analysis by using perl scripts from the miRGrep tool and an in-house built shell script for Minimum Fold Energy Index (MFEI) calculation. Applying the stringent criteria of the Dicer cleavage pattern and the perfect stem loop structure, precursors for 57 known miRNAs of 15 families and 18 novel miRNAs were revealed. Quantitative Real Time (qRT) PCR was performed to determine the expression of selected miRNAs.

Functional Roles of microRNAs in Agronomically Important Plants-Potential as Targets for Crop Improvement and Protection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have recently emerged as important regulators of gene expression, mainly through cleavage and/or translation inhibition of the target mRNAs during or after transcription. miRNAs play important roles by regulating a multitude of biological processes in plants which include maintenance of genome integrity, development, metabolism, and adaptive responses toward environmental stresses. The increasing population of the world and their food demands requires focused efforts for the improvement of crop plants to ensure sustainable food production. Manipulation of mRNA transcript abundance via miRNA control provides a unique strategy for modulating differential plant gene expression and miRNAs are thus emerging as the next generation targets for genetic engineering for improvement of the agronomic properties of crops. However, a deeper understanding of its potential and the mechanisms involved will facilitate the design of suitable strategies to obtain the desirable traits with minimum trade-offs in the modified crops. In this regard, this review highlights the diverse roles of conserved and newly identified miRNAs in various food and industrial crops and recent advances made in the uses of miRNAs to improve plants of agronomically importance so as to significantly enhance crop yields and increase tolerance to various environmental stress agents of biotic-or abiotic origin.

isomiR2Function: An Integrated Workflow for Identifying MicroRNA Variants in Plants

In plants, post transcriptional regulation by non-coding RNAs (ncRNAs), in particular miRNAs (19-24 nt) has been involved in modulating the transcriptional landscape in developmental, biotic and abiotic interactions. In past few years, considerable focus has been leveraged on delineating and deciphering the role of miRNAs and their canonical isomiRs in plants. However, proper classification and accurate prediction of plant isomiRs taking into account the relative features by which we define isomiRs, such as templated or non-templated is still lacking. In the present research, we present isomiR2Function, a standalone easily deployable tool that allows for the robust and high-throughput discovery of templated and non-templated isomiRs. Additionally, isomiR2Function allows for identification of differentially expressed isomiRs and in parallel target prediction based on both transcripts or PARE-Seq either using Targetfinder or Cleaveland. isomiR2Function allows for the functional enrichment of the detected targets using TopGO package. Benchmarking of isomiR2Function revealed highly accurate prediction and classification of isomiRs as compared to the previously developed isomiR prediction tools. Additionally, the downstream implementation of additional features allows isomiR2Function to be classified as a single standalone tool for isomiR profiling from discovery to functional roles. All in all, isomiR2Function allows the streamline processing of the miRNA-seq for the identification and characterization of isomiRs with minimal efforts. isomiR2Function can be accessed through: https://github.com/347033139/isomiR2Function.

Small RNA and methylation responses in susceptible and tolerant landraces of cassava infected with South African cassava mosaic virus

Endogenous small RNAs (sRNAs) associated with gene regulatory mechanisms respond to virus infection, and virus-derived small RNAs (vsRNAs) have been implicated in recovery or symptom remission in some geminivirus-host interactions. Transcriptional gene silencing (TGS) (24 nt vsRNAs) and post transcriptional gene silencing (PTGS) (21-23 nt vsRNAs) have been associated with geminivirus intergenic (IR) and coding regions, respectively. In this Illumina deep sequencing study, we compared for the first time, the small RNA response to South African cassava mosaic virus (SACMV) of cassava landrace TME3 which shows a recovery and tolerant phenotype, and T200, a highly susceptible landrace. Interestingly, different patterns in the percentage of SACMV-induced normalized total endogenous sRNA reads were observed between T200 and TME3. Notably in virus-infected T200 there was an increase in 21 nt sRNAs during the early pre-symptomatic response (12dpi) compared to mock, while in TME3, the 22 nt sRNA size class was predominant at 32days post infection with SACMV. While vsRNAs of 21-24 nt size classes mapped to the entire SACMV DNA-A and DNA-B genome components in T200 and TME3, vsRNA population counts were lower at 32 (symptomatic stage) and 67 dpi (recovery stage) in tolerant TME3 compared with T200 (non-recovery). It is suggested that the high accumulation of primary vsRNAs, which correlated with high virus titers and severe symptoms in susceptible T200, may be due to failure to target SACMV-derived mRNA. Likewise, in contrast, in TME3 low vsRNA counts may represent efficient PTGS of viral mRNA, leading to a depletion/sequestration of vsRNA populations, supporting a role for PTGS in tolerance/recovery in TME3. Notably, in TME3 at recovery (67 dpi) the percentage (expressed as a percentage of total vsRNA counts) of redundant and non-redundant (unique) 24 nt vsRNAs increased dramatically. Since methylation of the SACMV genome was not detected by bisulfite sequencing, and vsRNA counts targeting the intergenic region (where the promoters reside) were very low in both the tolerant or susceptible landraces, we could not provide conclusive evidence that 24 nt vsRNA-mediated RNA directed genome methylation plays a central role in disease phenotype in these landraces, notwithstanding recognition for a possible role in histone modification in TME3. This work represents an important step toward understanding variable roles of sRNAs in different cassava genotype-geminivirus interactions.