Analysis of phytoplasma-responsive sRNAs provide insight into the pathogenic mechanisms of mulberry yellow dwarf disease

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Molecular perspectives on age-related resistance of plants to (viral) pathogens

Age-related resistance to microbe invasion is a commonly accepted concept in plant pathology. However, the impact of such age-dependent interactive phenomena is perhaps not yet sufficiently recognized by the broader plant science community. Toward cataloging an understanding of underlying mechanisms, this review explores recent molecular studies and their relevance to the concept. Examples describe differences in genetic background, transcriptomics, hormonal balances, protein-mediated events, and the contribution by short RNA-controlled gene silencing events. Throughout, recent findings with viral systems are highlighted as an illustration of the complexity of the interactions. It will become apparent that instead of uncovering a unifying explanation, we unveiled only trends. Nevertheless, with a degree of confidence, we propose that the process of plant age-related defenses is actively regulated at multiple levels. The overarching goal of this control for plants is to avoid a constitutive waste of resources, especially at crucial metabolically draining early developmental stages.

Small RNA Profiling of Aster Yellows Phytoplasma-Infected Catharanthus roseus Plants Showing Different Symptoms

Micropropagated Catharantus roseus plants infected with 'Candidatus Phytoplasma asteris' showed virescence symptoms, witches' broom symptoms, or became asymptomatic after their planting in pots. Nine plants were grouped into three categories according to these symptoms, which were then employed for investigation. The phytoplasma concentration, as determined by qPCR, correlated well with the severity of symptoms. To reveal the changes in the small RNA profiles in these plants, small RNA high-throughput sequencing (HTS) was carried out. The bioinformatics comparison of the micro (mi) RNA and small interfering (si) RNA profiles of the symptomatic and asymptomatic plants showed changes, which could be correlated to some of the observed symptoms. These results complement previous studies on phytoplasmas and serve as a starting point for small RNA-omic studies in phytoplasma research.

The miR393-Target Module Regulates Plant Development and Responses to Biotic and Abiotic Stresses

MicroRNAs (miRNAs), a class of endogenous small RNAs, are broadly involved in plant development, morphogenesis and responses to various environmental stresses, through manipulating the cleavage, translational expression, or DNA methylation of target mRNAs. miR393 is a conserved miRNA family present in many plants, which mainly targets genes encoding the transport inhibitor response1 (TIR1)/auxin signaling F-box (AFB) auxin receptors, and thus greatly affects the auxin signal perception, Aux/IAA degradation, and related gene expression. This review introduces the advances made on the miR393/target module regulating plant development and the plant’s responses to biotic and abiotic stresses. This module is valuable for genetic manipulation of optimized conditions for crop growth and development and would also be helpful in improving crop yield through molecular breeding.

Preparation of Wood-Based Panel Composites with Poplar Veneer as the Surface Layer Modified by In-Situ Polymerization of Active Monomers

Wood-based panels covered by melamine-impregnated paper are widely used in floors and furniture, due to its good surface texture, hardness, wear resistance, and waterproof function. However, there are still some problems, such as formaldehyde release from the impregnated resin, non-wood touch, and complex preparation processes. Therefore, this study designed glycidyl methacrylate (GMA) and ethyleneglycol dimethacrylate (EGDMA), combined with maleic anhydride (MAN) as a reactive catalyst, to build an active monomers system. It was first impregnated into poplar veneers, and then in-situ polymerized within the veneer using a hot pressing process, which realized the gluing of the veneer onto the wood-based panel substrate, synchronously. Such treatment aims to obtain wood-based panel composites decorated by the modified veneer, with real solid wood touch feeling, satisfied surface properties, and environment friendly glue bonding. The results indicated that the optimized reaction ratio of the active monomers (GMA:EGDMA) was 2:1 (molar ratio), and the maleic anhydride addition accounted for 6 wt.% of the total monomers. Under the optimized hot pressing condition, the modified veneer closely bonded to the wood-based panel substrate without obvious interfacial gaps. The hardness, abrasion resistance, modulus of rupture, and water resistance of the composites were significantly improved. Such results indicate that the treatment realized the perfect merging of solid wood touch feeling, environment friendly feature, and excellent properties of the composite. It was highly expected to replace the traditional melamine-impregnated paper to decorate wood-based panels, and could be potentially applied as surface decorating materials in wide areas of desktop, floor, cupboard, wardrobe, and so on.

Cell Wall Bulking by Maleic Anhydride for Wood Durability Improvement

Wood is susceptible to swelling deformation and decay fungi due to moisture adsorption that originates from the dynamic nanopores of the cell wall and the abundant hydroxyl groups in wood components. This study employed as a modifier maleic anhydride (MAn), with the help of acetone as solvent, to diffuse into the wood cell wall, bulk nanopores, and further chemically bond to the hydroxyl groups of wood components, reducing the numbers of free hydroxyl groups and weakening the diffusion of water molecules into the wood cell wall. The derived MAn-bulked wood, compared to the control wood, presented a reduction in water absorptivity (RWA) of ~23% as well as an anti-swelling efficiency (ASE) of ~39% after immersion in water for 228 h, and showed an improvement in decay resistance of 81.42% against white-rot fungus and 69.79% against brown-rot fungus, respectively. The method of combined cell wall bulking and hydroxyl group bonding could effectively improve the dimensional stability and decay resistance with lower doses of modifier, providing a new strategy for wood durability improvement.

Time to Fight: Molecular Mechanisms of Age-Related Resistance

Plant age is a crucial factor in determining the outcome of a host-pathogen interaction. In successive developmental stages throughout their life cycles, plants face dynamic changes in biotic and abiotic conditions that create distinct ecological niches for host-pathogen interactions. As an adaptive strategy, plants have evolved intrinsic regulatory networks that integrate developmental signals with those from pathogen perception and defense activation. As a result, amplitude and timing of defense responses are optimized, so as to balance the cost and benefit of immunity activation. A general term "age-related resistance" refers to a gain of disease resistance against a certain pathogen when plants reach a relatively mature stage. Age-related resistance is a common observation on fruits, vegetables, and row crops for their resistance against viruses, bacteria, fungi, oomycetes pathogens, and insects. This review focuses on the recent advances in understanding the molecular mechanisms of how plants coordinate developmental timing and immune response.

Plant Hormones in Phytoplasma Infected Plants

Phytoplasmas are bacterial plant pathogens that need a plant host and an insect vector for their spread and survival. In plants, the physiological responses that phytoplasmas trigger result in symptom development through effects on hormonal, nutritional, and stress signaling pathways, and the interactions between these. In this review, recent advances on the involvement of plant hormones together with their known and deduced roles in plants infected with phytoplasmas are discussed. Several studies have directly, or in many cases indirectly, addressed plant hormone systems in phytoplasma-infected plants. These have provided accumulating evidence that phytoplasmas extensively affect plant hormone pathways. Phytoplasmas thus, with disturbing complex plant hormone networks, suppress plant immunity and modify plant structure, while optimizing their nutrient acquisition and facilitating their colonization of the plants, and their dissemination among plants by their insect vectors.

Alterations of plant architecture and phase transition by the phytoplasma virulence factor SAP11

As key mediators linking developmental processes with plant immunity, TCP (TEOSINTE-BRANCHED, CYCLOIDEA, PROLIFERATION FACTOR 1 and 2) transcription factors have been increasingly shown to be targets of pathogenic effectors. We report here that TB/CYC (TEOSINTE-BRANCHED/CYCLOIDEA)-TCPs are destabilized by phytoplasma SAP11 effectors, leading to the proliferation of axillary meristems. Although a high degree of sequence diversity was observed among putative SAP11 effectors identified from evolutionarily distinct clusters of phytoplasmas, these effectors acquired fundamental activity in destabilizing TB/CYC-TCPs. In addition, we demonstrate that miR156/SPLs and miR172/AP2 modules, which represent key regulatory hubs involved in plant phase transition, were modulated by Aster Yellows phytoplasma strain Witches' Broom (AY-WB) protein SAP11. A late-flowering phenotype with significant changes in the expression of flowering-related genes was observed in transgenic Arabidopsis plants expressing SAP11AYWB. These morphological and molecular alterations were correlated with the ability of SAP11 effectors to destabilize CIN (CINCINNATA)-TCPs. Although not all putative SAP11 effectors display broad-spectrum activities in modulating morphological and physiological changes in host plants, they serve as core virulence factors responsible for the witches' broom symptom caused by phytoplasmas.

Integrated Phloem Sap mRNA and Protein Expression Analysis Reveals Phytoplasma-infection Responses in Mulberry

To gain insight into the response of mulberry to phytoplasma-infection, the expression profiles of mRNAs and proteins in mulberry phloem sap were examined. A total of 955 unigenes and 136 proteins were found to be differentially expressed between the healthy and infected phloem sap. These differentially expressed mRNAs and proteins are involved in signaling, hormone metabolism, stress responses, etc. Interestingly, we found that both the mRNA and protein levels of the major latex protein-like 329 (MuMLPL329) gene were increased in the infected phloem saps. Expression of the MuMLPL329 gene was induced by pathogen inoculation and was responsive to jasmonic acid. Ectopic expression of MuMLPL329 in Arabidopsis enhances transgenic plant resistance to Botrytis cinerea, Pseudomonas syringae pv tomato DC3000 (Pst. DC3000) and phytoplasma. Further analysis revealed that MuMLPL329 can enhance the expression of some defense genes and might be involved in altering flavonoid content resulting in increased resistance of plants to pathogen infection. Finally, the roles of the differentially expressed mRNAs and proteins and the potential molecular mechanisms of their changes were discussed. It was likely that the phytoplasma-responsive mRNAs and proteins in the phloem saps were involved in multiple pathways of mulberry responses to phytoplasma-infection, and their changes may be partially responsible for some symptoms in the phytoplasma infected plants.

ceRNA Cross-Talk in Paulownia Witches' Broom Disease

Long noncoding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA) are important in the regulation of life activities. However, their function is unclear in Paulownia fortunei. To identify lncRNAs, circRNAs, and miRNA, and investigate their roles in the infection progress of Paulownia witches’ broom (PaWB) disease, we performed RNA sequencing of healthy and infected P. fortunei. A total of 3126 lncRNAs, 1634 circRNAs, and 550 miRNAs were identified. Among them, 229 lncRNAs, 65 circRNAs, and 65 miRNAs were differentially expressed in a significant manner. We constructed a competing endogenous RNA (ceRNA) network, which contains 5 miRNAs, 4 circRNAs, 5 lncRNAs, and 15 mRNAs, all of which were differentially expressed between healthy and infected P. fortunei. This study provides the first catalog of candidate ceRNAs in Paulownia and gives a revealing insight into the molecular mechanism responsible for PaWB.

miRVIT: A Novel miRNA Database and Its Application to Uncover Vitis Responses to Flavescence dorée Infection

Micro(mi)RNAs play crucial roles in plant developmental processes and in defense responses to biotic and abiotic stresses. In the last years, many works on small RNAs in grapevine (Vitis spp.) were published, and several conserved and putative novel grapevine-specific miRNAs were identified. In order to reorganize the high quantity of available data, we produced "miRVIT," the first database of all novel grapevine miRNA candidates characterized so far, and still not deposited in miRBase. To this aim, each miRNA accession was renamed, repositioned in the last version of the grapevine genome, and compared with all the novel and conserved miRNAs detected in grapevine. Conserved and novel miRNAs cataloged in miRVIT were then used for analyzing Vitis vinifera plants infected by Flavescence dorée (FD), one of the most severe phytoplasma diseases affecting grapevine. The analysis of small RNAs from healthy, recovered (plants showing spontaneous and stable remission of symptoms), and FD-infected "Barbera" grapevines showed that FD altered the expression profiles of several miRNAs, including those involved in cell development and photosynthesis, jasmonate signaling, and disease resistance response. The application of miRVIT in a biological context confirmed the effectiveness of the followed approach, especially for the identification of novel miRNA candidates in grapevine. miRVIT database is available at http://mirvit.ipsp.cnr.it. Highlights: The application of the newly produced database of grapevine novel miRNAs to the analysis of plants infected by Flavescence dorée reveals key roles of miRNAs in photosynthesis and jasmonate signaling.

Plant Responses to Pathogen Attack: Small RNAs in Focus

Small RNAs (sRNA) are a significant group of gene expression regulators for multiple biological processes in eukaryotes. In plants, many sRNA silencing pathways produce extensive array of sRNAs with specialized roles. The evidence on record advocates for the functions of sRNAs during plant microbe interactions. Host sRNAs are reckoned as mandatory elements of plant defense. sRNAs involved in plant defense processes via different pathways include both short interfering RNA (siRNA) and microRNA (miRNA) that actively regulate immunity in response to pathogenic attack via tackling pathogen-associated molecular patterns (PAMPs) and other effectors. In response to pathogen attack, plants protect themselves with the help of sRNA-dependent immune systems. That sRNA-mediated plant defense responses play a role during infections is an established fact. However, the regulations of several sRNAs still need extensive research. In this review, we discussed the topical advancements and findings relevant to pathogen attack and plant defense mediated by sRNAs. We attempted to point out diverse sRNAs as key defenders in plant systems. It is hoped that sRNAs would be exploited as a mainstream player to achieve food security by tackling different plant diseases.

MiRNA-seq-based profiles of miRNAs in mulberry phloem sap provide insight into the pathogenic mechanisms of mulberry yellow dwarf disease

A wide range of miRNAs have been identified as phloem-mobile molecules that play important roles in coordinating plant development and physiology. Phytoplasmas are associated with hundreds of plant diseases, and the pathogenesis involved in the interactions between phytoplasmas and plants is still poorly understood. To analyse the molecular mechanisms of phytoplasma pathogenicity, the miRNAs profiles in mulberry phloem saps were examined in response to phytoplasma infection. A total of 86 conserved miRNAs and 19 novel miRNAs were identified, and 30 conserved miRNAs and 13 novel miRNAs were differentially expressed upon infection with phytoplasmas. The target genes of the differentially expressed miRNAs are involved in diverse signalling pathways showing the complex interactions between mulberry and phytoplasma. Interestingly, we found that mul-miR482a-5p was up-regulated in the infected phloem saps, and grafting experiments showed that it can be transported from scions to rootstock. Based on the results, the complexity and roles of the miRNAs in phloem sap and the potential molecular mechanisms of their changes were discussed. It is likely that the phytoplasma-responsive miRNAs in the phloem sap modulate multiple pathways and work cooperatively in response to phytoplasma infection, and their expression changes may be responsible for some symptoms in the infected plants.

Plant Small Non-coding RNAs and Their Roles in Biotic Stresses

Non-coding RNAs (ncRNAs) have emerged as critical components of gene regulatory networks across a plethora of plant species. In particular, the 20-30 nucleotide small ncRNAs (sRNAs) play important roles in mediating both developmental processes and responses to biotic stresses. Based on variation in their biogenesis pathways, a number of different sRNA classes have been identified, and their specific functions have begun to be characterized. Here, we review the current knowledge of the biogenesis of the primary sRNA classes, microRNA (miRNA) and small nuclear RNA (snRNA), and their respective secondary classes, and discuss the roles of sRNAs in plant-pathogen interactions. sRNA mobility between species is also discussed along with potential applications of sRNA-plant-pathogen interactions in crop improvement technologies.

Effective enhancement of resistance to Phytophthora infestans by overexpression of miR172a and b in Solanum lycopersicum

Overexpression of miR172a and b in tomato ( Solanum lycopersicum ) Zaofen No. 2 increased resistance to Phytophthora infestans infection by suppressing of an AP2/ERF transcription factor. The miR172 family has been shown to participate in the growth phase transition, flowering time control, abiotic and biotic stresses by regulating the expression of a small group of AP2/ERF transcription factors. In this study, the precursors of miR172a and b were cloned from tomato, Solanum pimpinellifolium L3708. We used the degradome sequencing to determine the cleavage site of miR172 to a member of the AP2/ERF transcription factor family (Solyc11g072600.1.1). qRT-PCR results showed that the expression of AP2/ERF was negatively correlated with the expression of miR172 in S. pimpinellifolium L3708 infected with Phytophthora infestans. Overexpression of miR172a and b in S. lycopersicum Zaofen No. 2 conferred greater resistance to P. infestans infection, as evidenced by decreased disease index, lesion sizes, and P. infestans abundance. The SOD and POD play important roles in scavenging late massive ROS in plant-pathogen interaction. Malonaldehyde (MDA) is widely recognized as an indicator of lipid peroxidation. Membrane damage in plants can be estimated by measuring leakage of electrolytes, which is evaluated by determining relative electrolyte leakage (REL). Less H₂O₂ and O₂^-, higher activities of POD and SOD, less MDA content and REL, and higher chlorophyll content and photosynthetic rate were also shown in transgenic plants after inoculation with P. infestans. Our results constitute the first step towards further investigations into the biological function and molecular mechanism of miR172-mediated silencing of AP2/ERF transcription factors in S. lycopersicum-P. infestans interaction and provide a candidate gene for breeding to enhance biotic stress-resistance in S. lycopersicum.

A Novel LncRNA, MuLnc1, Associated With Environmental Stress in Mulberry (Morus multicaulis)

Environmental stresses are major constraints that limit the leaf productivity and quality of mulberry. LncRNAs have emerged as important regulators in response to biotic and abiotic stresses in plants. However, the functions and mechanisms of most lncRNAs remain largely unknown. A novel lncRNA designated as MuLnc1 was found to be cleaved by mul-miR3954 and produce secondary siRNAs in a 21 nt phase in mulberry. It was demonstrated that one of the siRNAs produced, si161579, can silence the expression of the calmodulin-like protein gene CML27 of mulberry (MuCML27). When MuCML27 was heterologously expressed in Arabidopsis, the transgenic plants exhibited enhanced resistance to Botrytis cinerea and Pseudomonas syringae pv tomato DC3000. In addition, the transgenic MuCML27-overexpressing Arabidopsis plants are more tolerant to salt and drought stresses. Furthermore, the network of mul-miR3954-MuLnc1-siRNAs-mRNAs was modeled to elucidate the interaction between lncRNAs and sRNAs with mRNAs. All of these, taken together, suggest that MuLnc1 was associated with environmental stress in mulberry and may be considered as a potential genetic improvement target gene of mulberry. The information provided may shed light on the complicated gene expression regulatory mechanisms in mulberry stress responses.

Characterization of microRNAs of Beta macrocarpa and their responses to Beet necrotic yellow vein virus infection

Plant microRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development, defense, and symptom development. Here, 547 known miRNAs representing 129 miRNA families, and 282 potential novel miRNAs were identified in Beta macrocarpa using small RNA deep sequencing. A phylogenetic analysis was performed, and 8 Beta lineage-specific miRNAs were identified. Through a differential expression analysis, miRNAs associated with Beet necrotic yellow vein virus (BNYVV) infection were identified and confirmed using a microarray analysis and stem-loop RT-qPCR. In total, 103 known miRNAs representing 38 miRNA families, and 45 potential novel miRNAs were differentially regulated, with at least a two-fold change, in BNYVV-infected plants compared with that of the mock-inoculated control. Targets of these differentially expressed miRNAs were also predicted by degradome sequencing. These differentially expressed miRNAs were involved in hormone biosynthesis and signal transduction pathways, and enhanced axillary bud development and plant defenses. This work is the first to describe miRNAs of the plant genus Beta and may offer a reference for miRNA research in other species in the genus. It provides valuable information on the pathogenicity mechanisms of BNYVV.

Comparative Proteomic Analysis of Paulownia fortunei Response to Phytoplasma Infection with Dimethyl Sulfate Treatment

Paulownia fortunei is a widely cultivated economic forest tree species that is susceptible to infection with phytoplasma, resulting in Paulownia witches' broom (PaWB) disease. Diseased P. fortunei is characterized by stunted growth, witches' broom, shortened internodes, and etiolated and smaller leaves. To understand the molecular mechanism of its pathogenesis, we applied isobaric tags for relative and absolute quantitation (iTRAQ) and liquid chromatography coupled with tandem mass spectrometry approaches to study changes in the proteomes of healthy P. fortunei, PaWB-infected P. fortunei, and PaWB-infected P. fortunei treated with 15 mg·L⁻¹ or 75 mg·L⁻¹ dimethyl sulfate. We identified 2969 proteins and 104 and 32 differentially abundant proteins that were phytoplasma infection responsive and dimethyl sulfate responsive, respectively. Based on our analysis of the different proteomes, 27 PaWB-related proteins were identified. The protein-protein interactions of these 27 proteins were analyzed and classified into four groups (photosynthesis-related, energy-related, ribosome-related, and individual proteins). These PaWB-related proteins may help in developing a deeper understanding of how PaWB affects the morphological characteristics of P. fortunei and further establish the mechanisms involved in the response of P. fortunei to phytoplasma.

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Quantitative proteome-level analysis of paulownia witches' broom disease with methyl methane sulfonate assistance reveals diverse metabolic changes during the infection and recovery processes

Paulownia witches' broom (PaWB) disease caused by phytoplasma is a fatal disease that leads to considerable economic losses. Although there are a few reports describing studies of PaWB pathogenesis, the molecular mechanisms underlying phytoplasma pathogenicity in Paulownia trees remain uncharacterized. In this study, after building a transcriptome database containing 67,177 sequences, we used isobaric tags for relative and absolute quantification (iTRAQ) to quantify and analyze the proteome-level changes among healthy P. fortunei (PF), PaWB-infected P. fortunei (PFI), and PaWB-infected P. fortunei treated with 20 mg L^-1 or 60 mg L^-1 methyl methane sulfonate (MMS) (PFI-20 and PFI-60, respectively). A total of 2,358 proteins were identified. We investigated the proteins profiles in PF vs. PFI (infected process) and PFI-20 vs. PFI-60 (recovered process), and further found that many of the MMS-response proteins mapped to "photosynthesis" and "ribosome" pathways. Based on our comparison scheme, 36 PaWB-related proteins were revealed. Among them, 32 proteins were classified into three functional groups: (1) carbohydrate and energy metabolism, (2) protein synthesis and degradation, and (3) stress resistance. We then investigated the PaWB-related proteins involved in the infected and recovered processes, and discovered that carbohydrate and energy metabolism was inhibited, and protein synthesis and degradation decreased, as the plant responded to PaWB. Our observations may be useful for characterizing the proteome-level changes that occur at different stages of PaWB disease. The data generated in this study may serve as a valuable resource for elucidating the pathogenesis of PaWB disease during phytoplasma infection and recovery stages.

Genome-wide identification and characterization of the SPL gene family in Ziziphus jujuba

SQUAMOSA Promoter-Binding Protein-Likes (SPLs) are plant specific transcription factors playing important roles in plant growth and development. The SPL gene family has been studied in various plant species; however, there is no report about SPLs in Zizyphus jujuba. In this study, we identified 18 putative ZjSPL genes in Z. jujuba using a genome-wide analysis. Sequence features, gene structures, conserved domains and motifs were analyzed. The phylogenetic relationships of SPLs in Z. jujuba and A. thaliana were revealed. A total of 5 pairs of ZjSPLs were identified, suggesting the importance of gene duplication in SPL gene expansion in Z. jujuba. In addition, 11 of the 18 ZjSPLs, belonging to G1, G2 and G5 subgroups, were found to be targets of miR156, suggesting the conservation of miR156-mediated posttranscriptional regulation in plants. Expression analysis revealed that eight ZjSPL genes were responsive to the infection of witches'-broom phytoplasma. Our results provide a basis for the further elucidation of the biological function of ZjSPLs and their regulation in witches'-broom disease.

Genetic diversity among mulberry genotypes from seven countries

Mulberry (Morus alba) is an economically and ecologically important, widespread woody plant. It has served human beings for over hundreds of years, and it is still widely used in pharmaceuticals, food industry and farming nowadays. Using modern techniques, deeper understanding in classification and conservation resources of mulberry leads to higher-efficiency hybrids among populations. Genetic polymorphisms among 42 mulberry genotypes from seven countries over Asia and South America were detected using 17 inter simple sequence repeat (ISSR) primers. A dendrogram was constructed using the similarity matrix among genotypes and a principal component analysis (PCA) was carried out to further identify and cluster the mulberry genotypes. In the 42 genotypes, 175 distinct bands were displayed, among which 169 were polymorphic bands (96.57%). The polymorphic information content of 17 ISSR primers ranged from 0.2921 to 0.3746 with the mean of 0.3494. And Nei's index and Shanon's information index averaged 0.116 and 0.174, respectively, indicating low diversity of mulberry. For further study, cluster analysis and PCA were carried out and the results were similar. 42 genotypes were grouped, showing some hybridized combinations. Additionally, a connection between mulberry diseases and their genotypes was noted, which indicates possible application for ISSR in studying disease resistance of mulberry.

Genome-Wide Identification and Analysis of MicroRNAs Involved in Witches'-Broom Phytoplasma Response in Ziziphus jujuba

MicroRNAs (miRNAs) play an important role in responding to biotic and abiotic stresses in plants. Jujube witches'-broom a phytoplasma disease of Ziziphus jujuba is prevalent in China and is a serious problem to the industry. However, the molecular mechanism of the disease is poorly understood. In this study, genome-wide identification and analysis of microRNAs in response to witches'-broom was performed. A total of 85 conserved miRNA unique sequences belonging to 32 miRNA families and 24 novel miRNA unique sequences, including their complementary miRNA* strands were identified from small RNA libraries derived from a uninfected and witches'-broom infected Z. jujuba plant. Differentially expressed miRNAs associated with Jujube witches'-broom disease were investigated between the two libraries, and 12 up-regulated miRNAs and 10 down- regulated miRNAs identified with more than 2 fold changes. Additionally, 40 target genes of 85 conserved miRNAs and 49 target genes of 24 novel miRNAs were predicted and their putative functions assigned. Using the modified 5'-RACE method, we confirmed that SPL and MYB were cleaved by miR156 and miR159, respectively. Our results provide insight into the molecular mechanisms of witches'-broom disease in Z. jujuba.

Transcriptome, microRNA, and degradome analyses of the gene expression of Paulownia with phytoplamsa

BACKGROUND

Paulownia witches' broom (PaWB) is a fatal disease of Paulownia caused by a phytoplasma. In previous studies, we found that plants with PaWB symptoms would revert to a healthy morphology after methyl methane sulfonate (MMS) treatment. To completely understand the gene expression profiles of the Paulownia-phytoplasma interaction, three high-throughput sequencing technologies were used to investigate changes of gene expression and microRNAs (miRNAs) in healthy Paulownia tomentosa plantlets, PaWB-infected plantlets, and PaWB-infected plantlets treated with 60 mg · L(-1) MMS.

METHODS

Transcriptome, miRNAs and degradome sequencing were performed to explore the global gene expression profiles in the process of Paulownia tomentosa with phytoplasma infection.

RESULTS

A total of 98,714 all-unigenes, 62 conserved miRNAs, and 35 novel miRNAs were obtained, among which 902 differentially expressed genes (DEGs) and 24 miRNAs were found to be associated with PaWB disease. Subsequently, the target genes of these miRNAs were predicted by degradome sequencing. Interestingly, we found that 19 target genes of these differentially expressed miRNAs were among the 902 DEGs. The targets of pau-miR156g, pau-miR403, and pau-miR166c were significantly up-regulated in the P. tomentosa plantlets infected with phytoplasma. Interaction of miRNA -target genes mediated gene expression related to PaWB were identified.

CONCLUSIONS

The results elucidated the possible roles of the regulation of genes and miRNAs in the Paulownia-phytoplasma interaction, which will enrich our understanding of the mechanisms of PaWB disease in this plant.

Plant-Pathogen Interaction-Related MicroRNAs and Their Targets Provide Indicators of Phytoplasma Infection in Paulownia tomentosa × Paulownia fortunei

Paulownia witches’ broom (PaWB) caused by a phytoplasma, has caused extensive losses in the yields of paulownia timber and resulted in significant economic losses. However, the molecular mechanisms in Paulownia that underlie the phytoplasma stress are poorly characterized. In this study, we use an Illumina platform to sequence four small RNA libraries and four degradome sequencing libraries derived from healthy, PaWB-infected, and PaWB-infected 15 mg·L⁻¹ and 30 mg·L⁻¹ methyl methane sulfonate (MMS)-treated plants. In total, 125 conserved and 118 novel microRNAs (miRNAs) were identified and 33 miRNAs responsive to PaWB disease were discovered. Furthermore, 166 target genes for 18 PaWB disease-related miRNAs were obtained, and found to be involved in plant-pathogen interaction and plant hormone signal transduction metabolic pathways. Eleven miRNAs and target genes responsive to PaWB disease were examined by a quantitative real-time PCR approach. Our findings will contribute to studies on miRNAs and their targets in Paulownia, and provide new insights to further understand plant-phytoplasma interactions.

Discovery of microRNAs and transcript targets related to witches' broom disease in Paulownia fortunei by high-throughput sequencing and degradome approach

Paulownia witches' broom (PaWB) caused by the phytoplasma is a devastating disease of Paulownia trees. It has caused heavy yield losses to Paulownia production worldwide. However, knowledge of the transcriptional and post-transcriptional regulation of gene expression by microRNAs (miRNAs), especially miRNAs responsive to PaWB disease stress, is still rudimentary. In this study, to identify miRNAs and their transcript targets that are responsive to PaWB disease stress, six sequencing libraries were constructed from healthy (PF), PaWB-infected (PFI), and PaWB-infected, 20 mg L(-1) methyl methane sulfonate-treated (PFI20) P. fortunei seedlings. As a result, 95 conserved miRNAs belonging to 18 miRNA families, as well as 122 potential novel miRNAs, were identified. Most of them were found to be a response to PaWB disease-induced stress, and the expression levels of these miRNAs were validated by quantitative real-time PCR analysis. The study simultaneously identified 109 target genes from the P. fortunei for 14 conserved miRNA families and 24 novel miRNAs by degradome sequencing. Furthermore, the functions of the miRNA targets were annotated based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The results presented here provide the groundwork for further analysis of miRNAs and target genes responsive to the PaWB disease stress, and could be also useful for addressing new questions to better understand the mechanisms of plant infection by phytoplasma in the future.

RNA-Seq profile of flavescence dorée phytoplasma in grapevine

BACKGROUND

The phytoplasma-borne disease flavescence dorée is still a threat to European viticulture, despite mandatory control measures and prophylaxis against the leafhopper vector. Given the economic importance of grapevine, it is essential to find alternative strategies to contain the spread, in order to possibly reduce the current use of harmful insecticides. Further studies of the pathogen, the vector and the mechanisms of phytoplasma-host interactions could improve our understanding of the disease. In this work, RNA-Seq technology followed by three de novo assembly strategies was used to provide the first comprehensive transcriptomics landscape of flavescence dorée phytoplasma (FD) infecting field-grown Vitis vinifera leaves.

RESULTS

With an average of 8300 FD-mapped reads per library, we assembled 347 sequences, corresponding to 215 annotated genes, and identified 10 previously unannotated genes, 15 polycistronic transcripts and three genes supposedly localized in the gaps of the FD92 draft genome. Furthermore, we improved the annotation of 44 genes with the addition of 5'/3' untranslated regions. Functional classification revealed that the most expressed genes were either related to translation and protein biosynthesis or hypothetical proteins with unknown function. Some of these hypothetical proteins were predicted to be secreted, so they could be bacterial effectors with a potential role in modulating the interaction with the host plant. Interestingly, qRT-PCR validation of the RNA-Seq expression values confirmed that a group II intron represented the FD genomic region with the highest expression during grapevine infection. This mobile element may contribute to the genomic plasticity that is necessary for the phytoplasma to increase its fitness and endorse host-adaptive strategies.

CONCLUSIONS

The RNA-Seq technology was successfully applied for the first time to analyse the FD global transcriptome profile during grapevine infection. Our results provided new insights into the transcriptional organization and gene structure of FD. This may represent the starting point for the application of high-throughput sequencing technologies to study differential expression in FD and in other phytoplasmas with an unprecedented resolution.