Deep profiling of potential substrate atlas of porcine epidemic diarrhea virus 3C-like protease.

Coronavirus (CoV) 3C-like protease (3CLpro) is essential for viral replication and is involved in immune escape by proteolyzing host proteins. Deep profiling the 3CLpro substrates in the host proteome extends our understanding of viral pathogenesis and facilitates antiviral drug discovery. Here, 3CLpro from porcine epidemic diarrhea virus (PEDV), an enteropathogenic CoV, was used as a model which to identify the potential 3CLpro cleavage motifs in all porcine proteins. We characterized the selectivity of PEDV 3CLpro at sites P5-P4'. We then compiled the 3CLpro substrate preferences into a position-specific scoring matrix and developed a 3CLpro profiling strategy to delineate the protein substrate landscape of CoV 3CLpro. We identified 1,398 potential targets in the porcine proteome containing at least one putative cleavage site and experimentally validated the reliability of the substrate degradome. The PEDV 3CLpro-targeted pathways are involved in mRNA processing, translation, and key effectors of autophagy and the immune system. We also demonstrated that PEDV 3CLpro suppresses the type 1 interferon (IFN-I) cascade via the proteolysis of multiple signaling adaptors in the retinoic acid-inducible gene I (RIG-I) signaling pathway. Our composite method is reproducible and accurate, with an unprecedented depth of coverage for substrate motifs. The 3CLpro substrate degradome establishes a comprehensive substrate atlas that will accelerate the investigation of CoV pathogenicity and the development of anti-CoV drugs.IMPORTANCECoronaviruses (CoVs) are major pathogens that infect humans and animals. The 3C-like protease (3CLpro) encoded by CoV not only cleaves the CoV polyproteins but also degrades host proteins and is considered an attractive target for the development of anti-CoV drugs. However, the comprehensive characterization of an atlas of CoV 3CLpro substrates is a long-standing challenge. Using porcine epidemic diarrhea virus (PEDV) 3CLpro as a model, we developed a method that accurately predicts the substrates of 3CLpro and comprehensively maps the substrate degradome of PEDV 3CLpro. Interestingly, we found that 3CLpro may simultaneously degrade multiple molecules responsible for a specific function. For instance, it cleaves at least four adaptors in the RIG-I signaling pathway to suppress type 1 interferon production. These findings highlight the complexity of the 3CLpro substrate degradome and provide new insights to facilitate the development of anti-CoV drugs.

In Vivo Cleavage Map Illuminates the Central Role of RNase E in Coding and Non-coding RNA Pathways.

Understanding RNA processing and turnover requires knowledge of cleavages by major endoribonucleases within a living cell. We have employed TIER-seq (transiently inactivating an endoribonuclease followed by RNA-seq) to profile cleavage products of the essential endoribonuclease RNase E in Salmonella enterica. A dominating cleavage signature is the location of a uridine two nucleotides downstream in a single-stranded segment, which we rationalize structurally as a key recognition determinant that may favor RNase E catalysis. Our results suggest a prominent biogenesis pathway for bacterial regulatory small RNAs whereby RNase E acts together with the RNA chaperone Hfq to liberate stable 3' fragments from various precursor RNAs. Recapitulating this process in vitro, Hfq guides RNase E cleavage of a representative small-RNA precursor for interaction with a mRNA target. In vivo, the processing is required for target regulation. Our findings reveal a general maturation mechanism for a major class of post-transcriptional regulators.

Degradomic Identification of Membrane Type 1-Matrix Metalloproteinase as an ADAMTS9 and ADAMTS20 Substrate.

The secreted metalloproteases ADAMTS9 and ADAMTS20 are implicated in extracellular matrix proteolysis and primary cilium biogenesis. Here, we show that clonal gene-edited RPE-1 cells in which ADAMTS9 was inactivated, and which constitutively lack ADAMTS20 expression, have morphologic characteristics distinct from parental RPE-1 cells. To investigate underlying proteolytic mechanisms, a quantitative terminomics method, terminal amine isotopic labeling of substrates was used to compare the parental and gene-edited RPE-1 cells and their medium to identify ADAMTS9 substrates. Among differentially abundant neo-amino (N) terminal peptides arising from secreted and transmembrane proteins, a peptide with lower abundance in the medium of gene-edited cells suggested cleavage at the Tyr314-Gly315 bond in the ectodomain of the transmembrane metalloprotease membrane type 1-matrix metalloproteinase (MT1-MMP), whose mRNA was also reduced in gene-edited cells. This cleavage, occurring in the MT1-MMP hinge, that is, between the catalytic and hemopexin domains, was orthogonally validated both by lack of an MT1-MMP catalytic domain fragment in the medium of gene-edited cells and restoration of its release from the cell surface by reexpression of ADAMTS9 and ADAMTS20 and was dependent on hinge O-glycosylation. A C-terminally semitryptic MT1-MMP peptide with greater abundance in WT RPE-1 medium identified a second ADAMTS9 cleavage site in the MT1-MMP hemopexin domain. Consistent with greater retention of MT1-MMP on the surface of gene-edited cells, pro-MMP2 activation, which requires cell surface MT1-MMP, was increased. MT1-MMP knockdown in gene-edited ADAMTS9/20-deficient cells restored focal adhesions but not ciliogenesis. The findings expand the web of interacting proteases at the cell surface, suggest a role for ADAMTS9 and ADAMTS20 in regulating cell surface activity of MT1-MMP, and indicate that MT1-MMP shedding does not underlie their observed requirement in ciliogenesis.

Small RNA and Degradome Deep Sequencing Reveal Regulatory Roles of MicroRNAs in Response to Sugarcane Mosaic Virus Infection on Two Contrasting Sugarcane Cultivars.

MicroRNAs (miRNAs) play an essential regulatory role in plant-virus interaction. However, few studies have focused on the roles of miRNAs and their targets after sugarcane mosaic virus (SCMV) infection in sugarcane. To address this issue, we conducted small RNA (sRNA) and degradome sequencing on two contrasting sugarcanes (SCMV-resistant 'Fuoguo1' [FG1] and susceptible 'Badila') infected by SCMV at five time points. A total of 1,578 miRNAs were profiled from 30 sRNA libraries, comprising 660 known miRNAs and 380 novel miRNAs. Differential expression analysis of miRNAs revealed that most were highly expressed during the SCMV exponential phase in Badila at 18 h postinfection, with expression profiles positively correlated with virus replication dynamics as observed through clustering. Analysis of degradome data indicated a higher number of differential miRNA targets in Badila compared to FG1 at 18 h postinfection. Gene ontology (GO) enrichment analysis significantly enriched the stimulus-response pathway, suggesting negative regulatory roles to SCMV resistance. Specifically, miR160 upregulated expression patterns and validated in Badila through quantitative real-time PCR (qRT-PCR) in the early stages of SCMV multiplication. Our research provides new insights into the dynamic response of plant miRNA and virus replication and contributes valuable information on the intricate interplay between miRNAs and SCMV infection dynamics. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Yeast Dxo1 is required for 25S rRNA maturation and acts as a transcriptome-wide distributive exonuclease.

The Dxo1/Rai1/DXO family of decapping and exonuclease enzymes can catalyze the in vitro removal of chemically diverse 5' ends from RNA. Specifically, these enzymes act poorly on RNAs with a canonical 7mGpppN cap, but instead prefer RNAs with a triphosphate, monophosphate, hydroxyl, or nonconventional cap. In each case, these enzymes generate an RNA with a 5' monophosphate, which is then thought to be further degraded by Rat1/Xrn1 5' exoribonucleases. For most Dxo1/Rai1/DXO family members, it is not known which of these activities is most important in vivo. Here we describe the in vivo function of the poorly characterized cytoplasmic family member, yeast Dxo1. Using RNA-seq of 5' monophosphate ends, we show that Dxo1 can act as a distributive exonuclease, removing a few nucleotides from endonuclease or decapping products. We also show that Dxo1 is required for the final 5' end processing of 25S rRNA, and that this is the primary role of Dxo1. While Dxo1/Rai1/DXO members were expected to act upstream of Rat1/Xrn1, this order is reversed in 25S rRNA processing, with Dxo1 acting downstream from Rat1. Such a hand-off from a processive to a distributive exonuclease may be a general phenomenon in the precise maturation of RNA ends.

First Insight into the Degradome of Aspergillus ochraceus: Novel Secreted Peptidases and Their Inhibitors.

Aspergillus fungi constitute a pivotal element within ecosystems, serving as both contributors of biologically active compounds and harboring the potential to cause various diseases across living organisms. The organism's proteolytic enzyme complex, termed the degradome, acts as an intermediary in its dynamic interaction with the surrounding environment. Using techniques such as genome and transcriptome sequencing, alongside protein prediction methodologies, we identified putative extracellular peptidases within Aspergillus ochraceus VKM-F4104D. Following manual annotation procedures, a total of 11 aspartic, 2 cysteine, 2 glutamic, 21 serine, 1 threonine, and 21 metallopeptidases were attributed to the extracellular degradome of A. ochraceus VKM-F4104D. Among them are enzymes with promising applications in biotechnology, potential targets and agents for antifungal therapy, and microbial antagonism factors. Thus, additional functionalities of the extracellular degradome, extending beyond mere protein substrate digestion for nutritional purposes, were demonstrated.

TRUEE; a bioinformatic pipeline to define the functional microRNA targetome of Arabidopsis.

Central to plant microRNA (miRNA) biology is the identification of functional miRNA-target interactions (MTIs). However, the complementarity basis of bioinformatic target prediction results in mostly false positives, and the degree of complementarity does not equate with regulation. Here, we develop a bioinformatic workflow named TRUEE (Targets Ranked Using Experimental Evidence) that ranks MTIs on the extent to which they are subjected to miRNA-mediated cleavage. It sorts predicted targets into high (HE) and low evidence (LE) groupings based on the frequency and strength of miRNA-guided cleavage degradome signals across multiple degradome experiments. From this, each target is assigned a numerical value, termed a Category Score, ranking the extent to which it is subjected to miRNA-mediated cleavage. As a proof-of-concept, the 428 Arabidopsis miRNAs annotated in miRBase were processed through the TRUEE pipeline to determine the miRNA 'targetome'. The majority of high-ranking Category Score targets corresponded to highly conserved MTIs, validating the workflow. Very few Arabidopsis-specific, Brassicaceae-specific, or Conserved-passenger miRNAs had HE targets with high Category Scores. In total, only several hundred MTIs were found to have Category Scores characteristic of currently known physiologically significance MTIs. Although non-exhaustive, clearly the number of functional MTIs is much narrower than many studies claim. Therefore, using TRUEE to numerically rank targets directly on experimental evidence has given insights into the scope of the functional miRNA targetome of Arabidopsis.

Integration of miRNA dynamics and drought tolerant QTLs in rice reveals the role of miR2919 in drought stress response.

To combat drought stress in rice, a major threat to global food security, three major quantitative trait loci for 'yield under drought stress' (qDTYs) were successfully exploited in the last decade. However, their molecular basis still remains unknown. To understand the role of secondary regulation by miRNA in drought stress response and their relation, if any, with the three qDTYs, the miRNA dynamics under drought stress was studied at booting stage in two drought tolerant (Sahbaghi Dhan and Vandana) and one drought sensitive (IR 20) cultivars. In total, 53 known and 40 novel differentially expressed (DE) miRNAs were identified. The primary drought responsive miRNAs were Osa-MIR2919, Osa-MIR3979, Osa-MIR159f, Osa-MIR156k, Osa-MIR528, Osa-MIR530, Osa-MIR2091, Osa-MIR531a, Osa-MIR531b as well as three novel ones. Sixty-one target genes that corresponded to 11 known and 4 novel DE miRNAs were found to be co-localized with the three qDTYs, out of the 1746 target genes identified. We could validate miRNA-mRNA expression under drought for nine known and three novel miRNAs in eight different rice genotypes showing varying degree of tolerance. From our study, Osa-MIR2919, Osa-MIR3979, Osa-MIR528, Osa-MIR2091-5p and Chr01_11911S14Astr and their target genes LOC_Os01g72000, LOC_Os01g66890, LOC_Os01g57990, LOC_Os01g56780, LOC_Os01g72834, LOC_Os01g61880 and LOC_Os01g72780 were identified as the most promising candidates for drought tolerance at booting stage. Of these, Osa-MIR2919 with 19 target genes in the qDTYs is being reported for the first time. It acts as a negative regulator of drought stress tolerance by modulating the cytokinin and brassinosteroid signalling pathway.

Target Landscape of Conserved Plant MicroRNAs and the Complexities of Their Ancient MicroRNA-Binding Sites.

In plants, microRNA (miRNA)-target interactions (MTIs) require high complementarity, a feature from which bioinformatic programs have predicted numerous and diverse targets for any given miRNA, promoting the idea of complex miRNA networks. Opposing this is a hypothesis of constrained miRNA specificity, in which functional MTIs are restricted to the few targets whose required expression output is compatible with the expression of the miRNA. To explore these opposing views, the bioinformatic pipeline Targets Ranked Using Experimental Evidence was applied to strongly conserved miRNAs to identity their high-evidence (HE) targets across species. For each miRNA family, HE targets predominantly consisted of homologs from one conserved target gene family (primary family). These primary families corresponded to the known canonical miRNA-target families, validating the approach. Very few additional HE target families were identified (secondary family), and if so, they were likely functionally related to the primary family. Many primary target families contained highly conserved nucleotide sequences flanking their miRNA-binding sites that were enriched in HE homologs across species. A number of these flanking sequences are predicted to form conserved RNA secondary structures that preferentially base pair with the miRNA-binding site, implying that these sites are highly structured. Our findings support a target landscape view that is dominated by the conserved primary target families, with a minority of either secondary target families or non-conserved targets. This is consistent with the constrained hypothesis of functional miRNA specificity, which potentially in part is being facilitated by features beyond complementarity.

Cytological observation and RNA-seq analysis reveal novel miRNAs high expression associated with the pollen fertility of neo-tetraploid rice.

BACKGROUND: Neo-tetraploid rice lines exhibit high fertility and strong heterosis and harbor novel specific alleles, which are useful germplasm for polyploid rice breeding. However, the mechanism of the fertility associated with miRNAs remains unknown. In this study, a neo-tetraploid rice line, termed Huaduo21 (H21), was used. Cytological observation and RNA-sequencing were employed to identify the fertility-related miRNAs in neo-tetraploid rice. RESULTS: H21 showed high pollen fertility (88.08%), a lower percentage of the pollen mother cell (PMC) abnormalities, and lower abnormalities during double fertilization and embryogenesis compared with autotetraploid rice. A total of 166 non-additive miRNAs and 3108 non-additive genes were detected between H21 and its parents. GO and KEGG analysis of non-additive genes revealed significant enrichments in the DNA replication, Chromosome and associated proteins, and Replication and repair pathways. Comprehensive multi-omics analysis identified 32 pairs of miRNA/target that were associated with the fertility in H21. Of these, osa-miR408-3p and osa-miR528-5p displayed high expression patterns, targeted the phytocyanin genes, and were associated with high pollen fertility. Suppression of osa-miR528-5p in Huaduo1 resulted in a low seed set and a decrease in the number of grains. Moreover, transgenic analysis implied that osa-MIR397b-p3, osa-miR5492, and osa-MIR5495-p5 might participate in the fertility of H21. CONCLUSION: Taken together, the regulation network of fertility-related miRNAs-targets pairs might contribute to the high seed setting in neo-tetraploid rice. These findings enhance our understanding of the regulatory mechanisms of pollen fertility associated with miRNAs in neo-tetraploid rice.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress.

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

Small RNA and Degradome Sequencing Reveal Roles of miRNAs in the Petal Color Fading of Malus Crabapple.

The Malus crabapple is an important woody ornamental plant. The fading of petals during its development significantly affects their ornamental value. Petal color is related to anthocyanin content and miRNAs play an important role in the post-transcriptional regulation of anthocyanin synthesis. However, the mechanisms underlying miRNA regulation of petal fading have rarely been studied. Transcriptome and small RNA sequencing of petals from the blooming phases of Malus. 'Indian Summer' varieties S1 (small bud), S2 (initial-flowering), and S3 (late-flowering) allowed us to identify 230 known miRNAs and 17 novel miRNAs, including 52 differentially expressed miRNAs which targeted 494 genes and formed 823 miRNA-target pairs. Based on the target gene annotation results, miRNA-target pairs were screened that may be involved in the fading process of Malus crabapple petals through three different pathways: anthocyanin synthesis, transport, and degradation, involving mcr-miR858-MYB1\MYB5 and mcr-miR396-McCHI inhibiting anthocyanin synthesis; mcr-miR167, mcr-miR390, mcr-miR535, and mcr-miR858 inhibiting anthocyanin transport from the cytoplasm to the vacuole by targeting ABC transporter genes (ABCB, ABCC, ABCD, and ABCG); and mcr-miR398 targeting the superoxide dismutase genes (CZSOD2 and CCS) to accelerate anthocyanin degradation. These findings offer a novel approach to understanding the mechanism of petal fading and serve as a reference for other plants with floral fading.

The miRNA-mRNA Regulatory Modules of Pinus massoniana Lamb. in Response to Drought Stress.

Masson pine (Pinus massoniana Lamb.) is a major fast-growing woody tree species and pioneer species for afforestation in barren sites in southern China. However, the regulatory mechanism of gene expression in P. massoniana under drought remains unclear. To uncover candidate microRNAs, their expression profiles, and microRNA-mRNA interactions, small RNA-seq was used to investigate the transcriptome from seedling roots under drought and rewatering in P. massoniana. A total of 421 plant microRNAs were identified. Pairwise differential expression analysis between treatment and control groups unveiled 134, 156, and 96 differential expressed microRNAs at three stages. These constitute 248 unique microRNAs, which were subsequently categorized into six clusters based on their expression profiles. Degradome sequencing revealed that these 248 differentially expressed microRNAs targeted 2069 genes. Gene Ontology enrichment analysis suggested that these target genes were related to translational and posttranslational regulation, cell wall modification, and reactive oxygen species scavenging. miRNAs such as miR482, miR398, miR11571, miR396, miR166, miRN88, and miRN74, along with their target genes annotated as F-box/kelch-repeat protein, 60S ribosomal protein, copper-zinc superoxide dismutase, luminal-binding protein, S-adenosylmethionine synthase, and Early Responsive to Dehydration Stress may play critical roles in drought response. This study provides insights into microRNA responsive to drought and rewatering in Masson pine and advances the understanding of drought tolerance mechanisms in Pinus.

Integrative Transcriptome, miRNAs, Degradome, and Phytohormone Analysis of Brassica rapa L. in Response to Plasmodiophora brassicae.

Clubroot is an infectious root disease caused by Plasmodiophora brassicae in Brassica crops, which can cause immeasurable losses. We analyzed integrative transcriptome, small RNAs, degradome, and phytohormone comprehensively to explore the infection mechanism of P. brassicae. In this study, root samples of Brassica rapa resistant line material BrT24 (R-line) and susceptible line material Y510-9 (S-line) were collected at four different time points for cytological, transcriptome, miRNA, and degradome analyses. We found the critical period of disease resistance and infection were at 0-3 DAI (days after inoculation) and 9-20 DAI, respectively. Based on our finding, we further analyzed the data of 9 DAI vs. 20 DAI of S-line and predicted the key genes ARF8, NAC1, NAC4, TCP10, SPL14, REV, and AtHB, which were related to clubroot disease development and regulating disease resistance mechanisms. These genes are mainly related to auxin, cytokinin, jasmonic acid, and ethylene cycles. We proposed a regulatory model of plant hormones under the mRNA-miRNA regulation in the critical period of P. brassicae infection by using the present data of the integrative transcriptome, small RNAs, degradome, and phytohormone with our previously published results. Our integrative analysis provided new insights into the regulation relationship of miRNAs and plant hormones during the process of disease infection with P. brassicae.

Integration of mRNA and miRNA Analysis Reveals the Post-Transcriptional Regulation of Salt Stress Response in Hemerocallis fulva.

MicroRNAs (miRNAs) belong to non-coding small RNAs which have been shown to take a regulatory function at the posttranscriptional level in plant growth development and response to abiotic stress. Hemerocallis fulva is an herbaceous perennial plant with fleshy roots, wide distribution, and strong adaptability. However, salt stress is one of the most serious abiotic stresses to limit the growth and production of Hemerocallis fulva. To identify the miRNAs and their targets involved in the salt stress resistance, the salt-tolerant H. fulva with and without NaCl treatment were used as materials, and the expression differences of miRNAs-mRNAs related to salt-tolerance were explored and the cleavage sites between miRNAs and targets were also identified by using degradome sequencing technology. In this study, twenty and three significantly differential expression miRNAs (p-value < 0.05) were identified in the roots and leaves of H. fulva separately. Additionally, 12,691 and 1538 differentially expressed genes (DEGs) were also obtained, respectively, in roots and leaves. Moreover, 222 target genes of 61 family miRNAs were validated by degradome sequencing. Among the DE miRNAs, 29 pairs of miRNA targets displayed negatively correlated expression profiles. The qRT-PCR results also showed that the trends of miRNA and DEG expression were consistent with those of RNA-seq. A gene ontology (GO) enrichment analysis of these targets revealed that the calcium ion pathway, oxidative defense response, microtubule cytoskeleton organization, and DNA binding transcription factor responded to NaCl stress. Five miRNAs, miR156, miR160, miR393, miR166, and miR396, and several hub genes, squamosa promoter-binding-like protein (SPL), auxin response factor 12 (ARF), transport inhibitor response 1-like protein (TIR1), calmodulin-like proteins (CML), and growth-regulating factor 4 (GRF4), might play central roles in the regulation of NaCl-responsive genes. These results indicate that non-coding small RNAs and their target genes that are related to phytohormone signaling, Ca2+ signaling, and oxidative defense signaling pathways are involved in H. fulva's response to NaCl stress.

Combined Omics Approaches Reveal Distinct Mechanisms of Resistance and/or Susceptibility in Sugar Beet Double Haploid Genotypes at Early Stages of Beet Curly Top Virus Infection.

Sugar beet is susceptible to Beet curly top virus (BCTV), which significantly reduces yield and sugar production in the semi-arid growing regions worldwide. Sources of genetic resistance to BCTV is limited and control depends upon insecticide seed treatments with neonicotinoids. Through double haploid production and genetic selection, BCTV resistant breeding lines have been developed. Using BCTV resistant (R) [KDH13; Line 13 and KDH4-9; Line 4] and susceptible (S) [KDH19-17; Line 19] lines, beet leafhopper mediated natural infection, mRNA/sRNA sequencing, and metabolite analyses, potential mechanisms of resistance against the virus and vector were identified. At early infection stages (2- and 6-days post inoculation), examples of differentially expressed genes highly up-regulated in the 'R' lines (vs. 'S') included EL10Ac5g10437 (inhibitor of trypsin and hageman factor), EL10Ac6g14635 (jasmonate-induced protein), EL10Ac3g06016 (ribosome related), EL10Ac2g02812 (probable prolyl 4-hydroxylase 10), etc. Pathway enrichment analysis showed differentially expressed genes were predominantly involved with peroxisome, amino acids metabolism, fatty acid degradation, amino/nucleotide sugar metabolism, etc. Metabolite analysis revealed significantly higher amounts of specific isoflavonoid O-glycosides, flavonoid 8-C glycosides, triterpenoid, and iridoid-O-glycosides in the leaves of the 'R' lines (vs. 'S'). These data suggest that a combination of transcriptional regulation and production of putative antiviral metabolites might contribute to BCTV resistance. In addition, genome divergence among BCTV strains differentially affects the production of small non-coding RNAs (sncRNAs) and small peptides which may potentially affect pathogenicity and disease symptom development.

Multi-omics analysis to visualize the dynamic roles of defense genes in the response of tea plants to gray blight.

Gray blight (GB) is one of the most destructive diseases of tea plants, causing considerable damage and productivity losses; however, the dynamic roles of defense genes during pathogen infection remain largely unclear. To explore the numerous molecular interactions associated with GB stress in tea plants, we employed transcriptome, sRNAome and degradome sequencing from 1 to 13 days post-inoculation (dpi) at 3-day intervals. The transcriptomics results showed that differentially expressed genes (DEGs) related to flavonoid synthesis, such as chalcone synthase (CHS) and phenylalanine ammonia-lyase (PAL), were particularly induced at 4 dpi. Consistent with this, the contents of catechins (especially gallocatechin), which are the dominant flavonoids in tea plants, also increased in the leaves of tea plants infected with GB. Combined analysis of the sRNAome and degradome revealed that microRNAs could mediate tea plant immunity by regulating DEG expression at the post-transcriptional level. Co-expression network analysis demonstrated that miR530b-ethylene responsive factor 96 (ERF96) and miRn211-thaumatin-like protein (TLP) play crucial roles in the response to GB. Accordingly, gene-specific antisense oligonucleotide assays suggested that suppressing ERF96 decreased the levels of reactive oxygen species (ROS), whereas suppressing TLP increased the levels of ROS. Furthermore, ERF96 was induced, but TLP was suppressed, in susceptible tea cultivars. Our results collectively demonstrate that ERF96 is a negative regulator and TLP is a positive regulator in the response of tea plants to GB. Taken together, our comprehensive integrated analysis reveals a dynamic regulatory network linked to GB stress in tea plants and provides candidate genes for improvement of tea plants.

Degradome comparison between wild and cultivated rice identifies differential targeting by miRNAs.

BACKGROUND: Small non-coding (s)RNAs are involved in the negative regulation of gene expression, playing critical roles in genome integrity, development and metabolic pathways. Targeting of RNAs by ribonucleoprotein complexes of sRNAs bound to Argonaute (AGO) proteins results in cleaved RNAs having precise and predictable 5` ends. While tools to study sliced bits of RNAs to confirm the efficiency of sRNA-mediated regulation are available, they are sub-optimal. In this study, we provide an improvised version of a tool with better efficiency to accurately validate sRNA targets. RESULTS: Here, we improvised the CleaveLand tool to identify additional micro (mi)RNA targets that belong to the same family and also other targets within a specified free energy cut-off. These additional targets were otherwise excluded during the default run. We employed these tools to understand the sRNA targeting efficiency in wild and cultivated rice, sequenced degradome from two rice lines, O. nivara and O. sativa indica Pusa Basmati-1 and analyzed variations in sRNA targeting. Our results indicate the existence of multiple miRNA-mediated targeting differences between domesticated and wild species. For example, Os5NG4 was targeted only in wild rice that might be responsible for the poor secondary wall formation when compared to cultivated rice. We also identified differential mRNA targets of secondary sRNAs that were generated after miRNA-mediated cleavage of primary targets. CONCLUSIONS: We identified many differentially targeted mRNAs between wild and domesticated rice lines. In addition to providing a step-wise guide to generate and analyze degradome datasets, we showed how domestication altered sRNA-mediated cascade silencing during the evolution of indica rice.

Degradome sequencing reveals an integrative miRNA-mediated gene interaction network regulating rice seed vigor.

BACKGROUND: It is well known that seed vigor is essential for agricultural production and rice (Oryza sativa L.) is one of the most important crops in the world. Though we previously reported that miR164c regulates rice seed vigor, but whether and how other miRNAs cooperate with miR164c to regulate seed vigor is still unknown. RESULTS: Based on degradome data of six RNA samples isolated from seeds of the wild-type (WT) indica rice cultivar 'Kasalath' as well as two modified lines in 'Kasalath' background (miR164c-silenced line [MIM164c] and miR164c overexpression line [OE164c]), which were subjected to either no aging treatment or an 8-day artificial aging treatment, 1247 different target transcripts potentially cleaved by 421 miRNAs were identified. The miRNA target genes were functionally annotated via GO and KEGG enrichment analyses. By STRING database assay, a miRNA-mediated gene interaction network regulating seed vigor in rice was revealed, which comprised at least four interconnected pathways: the miR5075-mediated oxidoreductase related pathway, the plant hormone related pathway, the miR164e related pathway, and the previously reported RPS27AA related pathway. Knockout and overexpression of the target gene Os02g0817500 of miR5075 decreased and enhanced seed vigor, respectively. By Y2H assay, the proteins encoded by five seed vigor-related genes, Os08g0295100, Os07g0633100, REFA1, OsPER1 and OsGAPC3, were identified to interact with Os02g0817500. CONCLUSIONS: miRNAs cooperate to regulate seed vigor in rice via an integrative gene interaction network comprising miRNA target genes and other functional genes. The result provided a basis for fully understanding the molecular mechanisms of seed vigor regulation.

Cultivar-specific miRNA-mediated RNA silencing in grapes.

MAIN CONCLUSION: In-depth comparative degradome analysis of two domesticated grape cultivars with diverse secondary metabolite accumulation reveals differential miRNA-mediated targeting. Small (s)RNAs such as micro(mi)RNAs and secondary small interfering (si) often work as negative switches of gene expression. In plants, it is well known that miRNAs target and cleave mRNAs that have high sequence complementarity. However, it is not known if there are variations in miRNA-mediated targeting between subspecies and cultivars that have been subjected to vast genetic modifications through breeding and other selections. Here, we have used PAREsnip2 tool for analysis of degradome datasets derived from two contrasting domesticated grape cultivars having varied fruit color, habit and leaf shape. We identified several interesting variations in sRNA targeting using degradome and 5'RACE analysis between two contrasting grape cultivars that was further correlated using RNA-seq analysis. Several of the differences we identified are associated with secondary metabolic pathways. We propose possible means by which sRNAs might contribute to diversity in secondary metabolites and other development pathways between two domesticated cultivars of grapes.