Whole transcriptome sequencing of Pseudomonas syringae pv. actinidiae-infected kiwifruit plants reveals species-specific interaction between long non-coding RNA and coding genes

An improved bacterial mRNA enrichment strategy in dual RNA sequencing to unveil the dynamics of plant-bacterial interactions

BACKGROUND

Dual RNA sequencing is a powerful tool that enables a comprehensive understanding of the molecular dynamics underlying plant-microbe interactions. RNA sequencing (RNA-seq) poses technical hurdles in the transcriptional analysis of plant-bacterial interactions, especially in bacterial transcriptomics, owing to the presence of abundant ribosomal RNA (rRNA), which potentially limits the coverage of essential transcripts. Therefore, to achieve cost-effective and comprehensive sequencing of the bacterial transcriptome, it is imperative to devise efficient methods for eliminating rRNA and enhancing the proportion of bacterial mRNA. In this study, we modified a strand-specific dual RNA-seq method with the goal of enriching the proportion of bacterial mRNA in the bacteria-infected plant samples. The enriched method involved the sequential separation of plant mRNA by poly A selection and rRNA removal for bacterial mRNA enrichment followed by strand specific RNA-seq library preparation steps. We assessed the efficiency of the enriched method in comparison to the conventional method by employing various plant-bacterial interactions, including both host and non-host resistance interactions with pathogenic bacteria, as well as an interaction with a beneficial rhizosphere associated bacteria using pepper and tomato plants respectively.

RESULTS

In all cases of plant-bacterial interactions examined, an increase in mapping efficiency was observed with the enriched method although it produced a lower read count. Especially in the compatible interaction with Xanthmonas campestris pv. Vesicatoria race 3 (Xcv3), the enriched method enhanced the mapping ratio of Xcv3-infected pepper samples to its own genome (15.09%; 1.45-fold increase) and the CDS (8.92%; 1.49-fold increase). The enriched method consistently displayed a greater number of differentially expressed genes (DEGs) than the conventional RNA-seq method at all fold change threshold levels investigated, notably during the early stages of Xcv3 infection in peppers. The Gene Ontology (GO) enrichment analysis revealed that the DEGs were predominantly enriched in proteolysis, kinase, serine type endopeptidase and heme binding activities.

CONCLUSION

The enriched method demonstrated in this study will serve as a suitable alternative to the existing RNA-seq method to enrich bacterial mRNA and provide novel insights into the intricate transcriptomic alterations within the plant-bacterial interplay.

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Harnessing the potential of non-coding RNA: An insight into its mechanism and interaction in plant biotic stress

Plants have developed diverse physical and chemical defence mechanisms to ensure their continued growth and well-being in challenging environments. Plants also have evolved intricate molecular mechanisms to regulate their responses to biotic stress. Non-coding RNA (ncRNA) plays a crucial role in this process that affects the expression or suppression of target transcripts. While there have been numerous reviews on the role of molecules in plant biotic stress, few of them specifically focus on how plant ncRNAs enhance resistance through various mechanisms against different pathogens. In this context, we explored the role of ncRNA in exhibiting responses to biotic stress endogenously as well as cross-kingdom regulation of transcript expression. Furthermore, we address the interplay between ncRNAs, which can act as suppressors, precursors, or regulators of other ncRNAs. We also delve into the regulation of ncRNAs in response to attacks from different organisms, such as bacteria, viruses, fungi, nematodes, oomycetes, and insects. Interestingly, we observed that diverse microorganisms interact with distinct ncRNAs. This intricacy leads us to conclude that each ncRNA serves a specific function in response to individual biotic stimuli. This deeper understanding of the molecular mechanisms involving ncRNAs in response to biotic stresses enhances our knowledge and provides valuable insights for future research in the field of ncRNA, ultimately leading to improvements in plant traits.

Two transcription factors, AcREM14 and AcC3H1, enhance the resistance of kiwifruit Actinidiachinensis var. chinensis to Pseudomonas syringae pv. actinidiae

Kiwifruit bacterial canker is a global disease caused by Pseudomonas syringae pv. actinidiae (Psa), which poses a major threat to kiwifruit production worldwide. Despite the economic importance of Actinidia chinensis var. chinensis, only a few resistant varieties have been identified to date. In this study, we screened 44 kiwifruit F1 hybrid lines derived from a cross between two A. chinensis var. chinensis lines and identified two offspring with distinct resistance to Psa: resistant offspring RH12 and susceptible offspring SH14. To identify traits associated with resistance, we performed a comparative transcriptomic analysis of these two lines. We identified several highly differentially expressed genes (DEGs) associated with flavonoid synthesis, pathogen interactions, and hormone signaling pathways, which play essential roles in disease resistance. Additionally, using weighted gene co-expression network analysis, we identified six core transcription factors. Moreover, qRT-PCR results demonstrated the high expression of AcC3H1 and AcREM14 in Psa-induced highly resistant hybrid lines. Ultimately, Overexpression of AcC3H1 and AcREM14 in kiwifruit enhanced disease resistance, and this was associated with upregulation of enzymatic activity and gene expression in the salicylic acid (SA) signaling pathway. Our study elucidates a molecular mechanism underlying disease resistance in kiwifruit and contributes to the advancement of research on kiwifruit breeding.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

Integrated Analysis of lncRNA–mRNA Regulatory Networks Related to Lipid Metabolism in High-Oleic-Acid Rapeseed

A high oleic acid content is considered an essential characteristic in the breeding of high-quality rapeseed in China. Long-chain non-coding RNA (lncRNA) molecules play an important role in the plant’s growth and its response to stress. To better understand the role of lncRNAs in regulating plant reproductive development, we analyzed whole-transcriptome and physiological data to characterize the dynamic changes in lncRNA expression during the four representative times of seed development of high- and low-oleic-acid rapeseed in three regions. We identified 21 and 14 lncRNA and mRNA modules, respectively. These modules were divided into three types related to region, development stages, and material. Next, we analyzed the key modules related to the oil content and the oleic acid, linoleic acid, and linolenic acid contents with physiological data and constructed the key functional network analysis on this basis. Genes related to lipid metabolism, such as 3-ketoacyl-CoA synthase 16 (KCS16) and acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1), were present in the co-expression network, suggesting that the effect of these genes on lipid metabolism might be embodied by the expression of these lncRNAs. Our results provide a fresh insight into region-, development-stage-, and material-biased changes in lncRNA expression in the seeds of Brassica napus. Some of these lncRNAs may participate in the regulatory network of lipid accumulation and metabolism, together with regulated genes. These results may help elucidate the regulatory system of lncRNAs in the lipid metabolism of high-oleic-acid rapeseed seeds.

Survival of Pseudomonas syringae pv. actinidiae in detached kiwifruit leaves at different environmental conditions

Pseudomonas syringae pv. actinidiae (Psa) is the causal agent of kiwifruit canker, a serious threat to commercial kiwifruit production worldwide. Studies of the movement path and the survival time of Psa in the host are crucial for integrated management programs. Hence, we used Psa with GFPuv gene (Psa-GFPuv) strain to investigate the movement path of Psa in leaves and branches, and the survival time of Psa in leaves under different environmental conditions. We found that the pathogen Psa spread longitudinally in the branches and leaves rather than transverse path. Additionally, the survival time of bacteria in fallen leaves under different environmental conditions were simulated by the way of Psa infecting the detached kiwifruit leaves. Psa survives the longest, up to 43 days in detached kiwifruit leaves with high humidity (above 80%) at 5 °C, and up to 32 days with low humidity (20%). At 15 °C, the Psa can survive in detached kiwifruit leaves for 20-30 days with increasing humidity. At 25 °C, it can only survive for 3 days with low humidity (20%) and 15 days with high humidity (above 80%). Furthermore, the population growth experiments showed that bacterial growth of Psa was more favorable in detached kiwifruit leaves with above 80% humidity at 5 °C. These results suggest that the survival condition of Psa in detached kiwifruit leaves is significantly affected by environmental conditions, and provide the basis for the control timing and technology of kiwifruit canker.

From molecular basics to agronomic benefits: Insights into noncoding RNA-mediated gene regulation in plants

The development of plants is largely dependent on their growth environment. To better adapt to a particular habitat, plants have evolved various subtle regulatory mechanisms for altering gene expression. Non coding RNAs (ncRNAs) constitute a major portion of the transcriptomes of eukaryotes. Various ncRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycles of plants. In this review, we summarize the current understanding of the biogenesis and contributions of small nucle olar RNA (snoRNA)- and regulatory long non coding RNA (lncRNA)-mediated gene regulation in plant development and environmental responses. Many regulatory ncRNAs appear to be associated with increased yield, quality and disease resistance of various species and cultivars. These ncRNAs may potentially be used as genetic resources for improving agronomic traits and for molecular breeding. The challenges in understanding plant ncRNA biology and the possibilities to make better use of these valuable gene resources in the future are discussed in this review.

Genomic Variation and Host Interaction among Pseudomonas syringae pv. actinidiae Strains in Actinidia chinensis ‘Hongyang’

Kiwifruit bacterial canker is a recent epidemic disease caused by Pseudomonas syringae pv. actinidiae (Psa), which has undergone worldwide expansion in a short time and resulted in significant economic losses. ‘Hongyang’ (Actinidia chinensis), a widely grown cultivar because of its health-beneficial nutrients and appreciated red-centered inner pericarp, is highly sensitive to Psa. In this work, ten Psa strains were isolated from ‘Hongyang’ and sequenced for genome analysis. The results indicated divergences in pathogenicity and pathogenic-related genes among the Psa strains. Significantly, the interruption at the 596 bp of HrpR in two low-pathogenicity strains reemphasized this gene, expressing a transcriptional regulator for the effector secretion system, as an important pathogenicity-associated locus of Psa. The transcriptome analysis of ‘Hongyang’ infected with different Psa strains was performed by RNA-seq of stem tissues locally (at the inoculation site) and systemically. Psa infection re-programmed the host genes expression, and the susceptibility to Psa might be attributed to the down-regulation of several genes involved in plant-pathogen interactions, especially calcium signaling transduction, as well as fatty acid elongation. This suppression was found in both low- and high-pathogenicity Psa inoculated tissues, but the effect was stronger with more virulent strains. Taken together, the divergences of P. syringae pv. actinidiae in pathogenicity, genome, and resulting transcriptomic response of A. chinensis provide insights into unraveling the molecular mechanism of Psa-kiwifruit interactions and resistance improvement in the kiwifruit crop.

The genome sequencing and comparative analysis of a wild kiwifruit Actinidia eriantha

The current kiwifruit industry is mainly based on the cultivars derived from the species Actinidia chinensis (Ac) which may bring risks such as canker disease. Introgression of desired traits from wild relatives is an important method for improving kiwifruit cultivars. Actinidia eriantha (Ae) is a particularly important taxon used for hybridization or introgressive breeding of new kiwifruit cultivars because of its valued species-specific traits. Here, we assembled a chromosome-scale high-quality genome of a Ae sample which was directly collected from its wild populations. Our analysis revealed that 41.3% of the genome consists of repetitive elements, comparable to the percentage in Ac and Ae cultivar "White" genomes. The genomic structural variation, including the presence/absence-variation (PAV) of genes, is distinct between Ae and Ac, despite both sharing the same two kiwifruit-specific whole genome duplication (WGD) events. This suggests that a post-WGD divergence mechanism occurred during their evolution. We further investigated genes involved in ascorbic acid biosynthesis and disease-resistance of Ae, and we found introgressive genome could contribute to the complex relationship between Ae and other representative kiwifruit taxa. Collectively, the Ae genome offers valuable genetic resource to accelerate kiwifruit breeding applications.

Genome-wide identification and characterization of Fusarium circinatum-responsive lncRNAs in Pinus radiata

Background

One of the most promising strategies of Pine Pitch Canker (PPC) management is the use of reproductive plant material resistant to the disease. Understanding the complexity of plant transcriptome that underlies the defence to the causal agent Fusarium circinatum, would greatly facilitate the development of an accurate breeding program. Long non-coding RNAs (lncRNAs) are emerging as important transcriptional regulators under biotic stresses in plants. However, to date, characterization of lncRNAs in conifer trees has not been reported. In this study, transcriptomic identification of lncRNAs was carried out using strand-specific paired-end RNA sequencing, from Pinus radiata samples inoculated with F. circinatum at an early stage of infection.

Results

Overall, 13,312 lncRNAs were predicted through a bioinformatics approach, including long intergenic non-coding RNAs (92.3%), antisense lncRNAs (3.3%) and intronic lncRNAs (2.9%). Compared with protein-coding RNAs, pine lncRNAs are shorter, have lower expression, lower GC content and harbour fewer and shorter exons. A total of 164 differentially expressed (DE) lncRNAs were identified in response to F. circinatum infection in the inoculated versus mock-inoculated P. radiata seedlings. The predicted cis-regulated target genes of these pathogen-responsive lncRNAs were related to defence mechanisms such as kinase activity, phytohormone regulation, and cell wall reinforcement. Co-expression network analysis of DE lncRNAs, DE protein-coding RNAs and lncRNA target genes also indicated a potential network regulating pectinesterase activity and cell wall remodelling.

Conclusions

This study presents the first comprehensive genome-wide analysis of P. radiata lncRNAs and provides the basis for future functional characterizations of lncRNAs in relation to pine defence responses against F. circinatum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08408-9.

Genome-wide identification and characterization of the TIFY gene family in kiwifruit

Background

The TIFY gene family is a group of plant-specific transcription factors involved in regulation of plant growth and development and a variety of stress responses. However, the TIFY family has not yet been well characterized in kiwifruit, a popular fruit with important nutritional and economic value.

Results

A total of 27 and 21 TIFY genes were identified in the genomes of Actinidia eriantha and A. chinensis, respectively. Phylogenetic analyses showed that kiwifruit TIFY genes could be classified into four major groups, JAZ, ZML, TIFY and PPD, and the JAZ group could be further clustered into six subgroups (JAZ I to JAZ VI). Members within the same group or subgroup have similar exon-intron structures and conserved motif compositions. The kiwifruit TIFY genes are unevenly distributed on the chromosomes, and the segmental duplication events played a vital role in the expansion of the TIFY genes in kiwifruit. Syntenic analyses of TIFY genes between kiwifruit and other five plant species (including Arabidopsis thaliana, Camellia sinensis, Oryza sativa, Solanum lycopersicum and Vitis vinifera) and between the two kiwifruit species provided valuable clues for understanding the potential evolution of the kiwifruit TIFY family. Molecular evolutionary analysis showed that the evolution of kiwifruit TIFY genes was primarily constrained by intense purifying selection. Promoter cis-element analysis showed that most kiwifruit TIFY genes possess multiple cis-elements related to stress-response, phytohormone signal transduction and plant growth and development. The expression pattern analyses indicated that TIFY genes might play a role in different kiwifruit tissues, including fruit at specific development stages. In addition, several TIFY genes with high expression levels during Psa (Pseudomonas syringae pv. actinidiae) infection were identified, suggesting a role in the process of Pas infection.

Conclusions

In this study, the kiwifruit TIFY genes were identified from two assembled kiwifruit genomes. In addition, their basic physiochemical properties, chromosomal localization, phylogeny, gene structures and conserved motifs, synteny analyses, promoter cis-elements and expression patters were systematically examined. The results laid a foundation for further understanding the function of TIFY genes in kiwifruit, and provided a new potential approach for the prevention and treatment of Psa infection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08398-8.

Defence-related pathways, phytohormones and primary metabolism are key players in kiwifruit plant tolerance to Pseudomonas syringae pv. actinidiae

The reasons underlying the differential tolerance of Actinidia spp. to the pandemic pathogen Pseudomonas syringae pv. actinidiae (Psa) have not yet been elucidated. We hypothesized that differential plant-defence strategies linked to transcriptome regulation, phytohormones and primary metabolism might be key and that Actinidia chinensis susceptibility results from an inefficient activation of defensive mechanisms and metabolic impairments shortly following infection. Here, 48 h postinoculation bacterial density was 10-fold higher in A. chinensis var. deliciosa than in Actinidia arguta, accompanied by significant increases in glutamine, ornithine, jasmonic acid (JA) and salicylic acid (SA) (up to 3.2-fold). Actinidia arguta showed decreased abscisic acid (ABA) (0.7-fold), no changes in primary metabolites, and 20 defence-related genes that were only differentially expressed in this species. These include GLOX1, FOX1, SN2 and RBOHA, which may contribute to its higher tolerance. Results suggest that A. chinensis' higher susceptibility to Psa is due to an inefficient activation of plant defences, with the involvement of ABA, JA and SA, leading to impairments in primary metabolism, particularly the ammonia assimilation cycle. A schematic overview on the interaction between Psa and genotypes with distinct tolerance is provided, highlighting the key transcriptomic and metabolomic aspects contributing to the different plant phenotypes after infection.

Long Non-Coding RNAs as Emerging Regulators of Pathogen Response in Plants

Long non-coding RNAs (lncRNAs) are transcripts without protein-coding potential that contain more than 200 nucleotides that play important roles in plant survival in response to different stresses. They interact with molecules such as DNA, RNA, and protein, and play roles in the regulation of chromatin remodeling, RNA metabolism, and protein modification activities. These lncRNAs regulate the expression of their downstream targets through epigenetic changes, at the level of transcription and post-transcription. Emerging information from computational biology and functional characterization of some of them has revealed their diverse mechanisms of action and possible roles in biological processes such as flowering time, reproductive organ development, as well as biotic and abiotic stress responses. In this review, we have mainly focused on the role of lncRNAs in biotic stress response due to the limited availability of knowledge in this domain. We have discussed the available molecular mechanisms of certain known lncRNAs against specific pathogens. Further, considering that fungal, viral, and bacterial diseases are major factors in the global food crisis, we have highlighted the importance of lncRNAs against pathogen responses and the progress in plant research to develop a better understanding of their functions and molecular mechanisms.

Analysis of long non-coding RNAs and mRNAs in harvested kiwifruit in response to the yeast antagonist, Wickerhamomyces anomalus

•

W. anomalus exhibits good biocontrol activity against blue and gray mold on kiwifruit.

•

LncRNAs in kiwifruit may be involved in activating plant hormone signal transduction pathways in response to the biocontrol yeast.

•

LncRNAs in kiwifruit may modulate the production of related TFs and secondary metabolites.

•

The expression of downstream defense-related genes in kiwifruit increases in response to the application of the biocontrol yeast.

Biological control utilizing antagonistic yeasts is an effective method for controlling postharvest diseases. Long non-coding RNAs (lncRNAs) have been found to be involved in a variety of plant growth and development processes, including those associated with plant disease resistance. In the present study, the yeast antagonist, Wickerhamomyces anomalus, was found to strongly inhibit postharvest blue mold (Penicillium expansum) and gray mold (Botrytis cinerea) decay of kiwifruit. Additionally, lncRNA high-throughput sequencing and bioinformatic analysis was used to identify lncRNAs in W. anomalus-treated wounds in kiwifruit and predict their function based on putative target genes. Our results indicate that lncRNAs may be involved in increasing ethylene (ET), jasmonic acid (JA), abscisic acid (ABA), and auxin (IAA) levels, as well as activating signal transduction pathways that regulate the expression of several transcription factors (WRKY72, WRKY53, JUB1AP2). These transcription factors (TFs) then mediate the expression of downstream, defense-related genes (ZAR1, PAD4, CCR4, NPR4) and the synthesis of secondary metabolites, thus, potentially enhancing disease resistance. Notably, by stimulating the accumulation of antifungal compounds, such as phenols and lignin, disease resistance in kiwifruit was enhanced. Our study provides new information on the mechanism underlying the induction of disease resistance in kiwifruit by W. anomalus, as well as a new disease resistance strategy that can be used to enhance the defense response of fruit to pathogenic fungi.

Transcriptomic and Proteomic Profiling Reveal the Key Role of AcMYB16 in the Response of Pseudomonas syringae pv. actinidiae in Kiwifruit

Kiwifruit bacterial canker caused by Pseudomonas syringae pv. actinidiae (Psa), is an important disease of kiwifruit (Actinidia Lind.). Plant hormones may induce various secondary metabolites to resist pathogens via modulation of hormone-responsive transcription factors (TFs), as reported in past studies. In this study, we showed that JA accumulated in the susceptible cultivar Actinidia chinensis 'Hongyang' but decreased in the resistant cultivar of A. chinensis var. deliciosa 'Jinkui' in response to Psa. Integrated transcriptomic and proteomic analyses were carried out using the resistant cultivar 'Jinkui'. A total of 5,045 differentially expressed genes (DEGs) and 1,681 differentially expressed proteins (DEPs) were identified after Psa infection. Two pathways, 'plant hormone signal transduction' and 'phenylpropanoid biosynthesis,' were activated at the protein and transcript levels. In addition, a total of 27 R2R3-MYB transcription factors (TFs) were involved in the response to Psa of 'Jinkui,' including the R2R3-MYB TF subgroup 4 gene AcMYB16, which was downregulated in 'Jinkui' but upregulated in 'Hongyang.' The promoter region of AcMYB16 has a MeJA responsiveness cis-acting regulatory element (CRE). Transient expression of the AcMYB16 gene in the leaves of 'Jinkui' induced Psa infection. Together, these data suggest that AcMYB16 acts as a repressor to regulate the response of kiwifruit to Psa infection. Our work will help to unravel the processes of kiwifruit resistance to pathogens and will facilitate the development of varieties with resistance against bacterial pathogens.

Expanding the application of a UV-visible reporter for transient gene expression and stable transformation in plants

Green fluorescent protein (GFP) has been widely used for monitoring gene expression and protein localization in diverse organisms. However, highly sensitive imaging equipment, like fluorescence microscope, is usually required for the visualization of GFP, limitings its application to fixed locations in samples. A reporter that can be visualized in real-time regardless the shape, size and location of the target samples will increase the flexibility and efficiency of research work. Here, we report the application of a GFP-like protein, called eYGFPuv, in both transient expression and stable transformation, in two herbaceous plant species (Arabidopsis and tobacco) and two woody plant species (poplar and citrus). We observed bright fluorescence under UV light in all of the four plant species without any effects on plant growth or development. eYGFPuv was shown to be effective for imaging transient expression in leaf and root tissues. With a focus on in vitro transformation, we demonstrated that the transgenic events expressing 1x eYGFPuv could be easily identified visually during the callus stage and the shoot stage, enabling early and efficient selection of transformants. Furthermore, whole-plant level visualization of eYGFPuv revealed its ubiquitous stability in transgenic plants. In addition, our transformation experiments showed that eYGFPuv can also be used to select transgenic plants without antibiotics. This work demonstrates the feasibility of utilizing 1x eYGFPuv in studies of gene expression and plant transformation in diverse plants.

Genome-wide analysis of long noncoding RNAs affecting floral bud dormancy in pears in response to cold stress

The versatile role of long noncoding RNAs (lncRNAs) in plant growth and development has been established, but a systematic identification and analysis of lncRNAs in the pear has not been reported. Bud dormancy is a crucial and complicated protective mechanism for plants in winter. The roles of lncRNAs in the dormancy process remain largely unclear. In this study, we induced pear floral buds to enter into different dormant statuses by simulating four different chilling accumulation conditions. Then, a time series of RNA-seq analysis was performed and we identified 7594 lncRNAs in Pyrus pyrifolia (Burm. F.) Nakai that have not been identified. The sequence and expression of the lncRNAs were confirmed by PCR analysis. In total, 6253 lncRNAs were predicted to target protein-coding genes including 692 cis-regulated pairs (596 lncRNAs) and 13,158 trans-regulated pairs (6181 lncRNAs). Gene Ontology analysis revealed that most of lncRNAs' target genes were involved in catalytic activity, metabolic processes and cellular processes. In the trend analysis, 124 long-term cold response lncRNAs and 80 short-term cold response lncRNAs were predicted. Regarding the lncRNA-miRNA regulatory networks, 59 lncRNAs were identified as potential precursors for miRNA members of 20 families, 586 lncRNAs were targets of 261 pear miRNAs and 53 lncRNAs were endogenous target mimics for 26 miRNAs. In addition, three cold response lncRNAs, two miRNAs and their target genes were selected for expression confirmed. The trend of their expression was consistent with the predicted relationships among them and suggested possible roles of lncRNAs in ABA metabolic pathway. Our findings not only suggest the potential roles of lncRNAs in regulating the dormancy of pear floral buds but also provide new insights into the lncRNA-miRNA-mRNA regulatory network in plants.

Integrated analysis of lncRNA and mRNA transcriptomes reveals the potential regulatory role of lncRNA in kiwifruit ripening and softening

Kiwifruit has gained increasing attention worldwide for its unique flavor and high nutritional value. Rapid softening after harvest greatly shortens its shelf-life and reduces the commercial value. Therefore, it is imperative and urgent to identify and clarify its softening mechanism. This study aimed to analyze and compare the long noncoding RNA (lncRNA) and mRNA expression patterns in ABA-treated (ABA) and room temperature (RT)-stored fruits with those in freshly harvested fruits (CK) as control. A total of 697 differentially expressed genes (DEGs) and 81 differentially expressed lncRNAs (DELs) were identified while comparing ABA with CK, and 458 DEGs and 143 DELs were detected while comparing RT with CK. The Kyoto Encyclopedia of Genes and Genomes analysis of the identified DEGs and the target genes of DELs revealed that genes involved in starch and sucrose metabolism, brassinosteroid biosynthesis, plant hormone signal transduction, and flavonoid biosynthesis accounted for a large part. The co-localization networks, including 38 DEGs and 31 DELs in ABA vs. CK, and 25 DEGs and 25 DELs in RT vs. CK, were also performed. Genes related to fruit ripening, such as genes encoding β-galactosidase, mannan endo-1,4-β-mannosidase, pectinesterase/pectinesterase inhibitor, and NAC transcription factor, were present in the co-localization network, suggesting that lncRNAs were involved in regulating kiwifruit ripening. Notably, several ethylene biosynthesis- and signaling-related genes, including one 1-aminocyclopropane-1-carboxylic acid oxidase gene and three ethylene response factor genes, were found in the co-localization network of ABA vs. CK, suggesting that the promoting effect of ABA on ethylene biosynthesis and fruit softening might be embodied by increasing the expression of these lncRNAs. These results may help understand the regulatory mechanism of lncRNAs in ripening and ABA-induced fruit softening of kiwifruit.

Kiwifruit Genome Database (KGD): a comprehensive resource for kiwifruit genomics

Kiwifruit (Actinidia spp.) plants produce economically important fruits containing abundant, balanced phytonutrients with extraordinarily high vitamin C contents. Since the release of the first kiwifruit reference genome sequence in 2013, large volumes of genome and transcriptome data have been rapidly accumulated for a handful of kiwifruit species. To efficiently store, analyze, integrate, and disseminate these large-scale datasets to the research community, we constructed the Kiwifruit Genome Database (KGD; http://kiwifruitgenome.org/). The database currently contains all publicly available genome and gene sequences, gene annotations, biochemical pathways, transcriptome profiles derived from public RNA-Seq datasets, and comparative genomic analysis results such as syntenic blocks and homologous gene pairs between different kiwifruit genome assemblies. A set of user-friendly query interfaces, analysis tools and visualization modules have been implemented in KGD to facilitate translational and applied research in kiwifruit, which include JBrowse, a popular genome browser, and the NCBI BLAST sequence search tool. Other notable tools developed within KGD include a genome synteny viewer and tools for differential gene expression analysis as well as gene ontology (GO) term and pathway enrichment analysis.

Early Pathogen Recognition and Antioxidant System Activation Contributes to Actinidia arguta Tolerance Against Pseudomonas syringae Pathovars actinidiae and actinidifoliorum

Actinidia chinensis and A. arguta have distinct tolerances to Pseudomonas syringae pv. actinidiae (Psa), but the reasons underlying the inter-specific variation remain unclear. This study aimed to integrate the metabolic and molecular responses of these two kiwifruit species against the highly pathogenic Psa and the less pathogenic P. syringae pv. actinidifoliorum (Pfm) bacterial strains. Disease development was monitored weekly till 21 days post inoculation (dpi), analysing a broad number and variety of parameters including: colony forming units (CFU), foliar symptoms, total chlorophylls, lipid peroxidation, soluble polyphenols, lignin and defense-related gene expression. At the end of the experimental period A. chinensis inoculated with Psa presented the highest endophytic bacterial population, whereas A. arguta inoculated with Pfm showed the lowest values, also resulting in a lower extent of leaf symptoms. Metabolic responses to infection were also more pronounced in A. chinensis with decreased total chlorophylls (up to 55%) and increased lipid peroxidation (up to 53%), compared with non-inoculated plants. Moreover, at 14 dpi soluble polyphenols and lignin concentrations were significantly higher (112 and 26%, respectively) in Psa-inoculated plants than in controls, while in A. arguta no significant changes were observed in those metabolic responses, except for lignin concentration which was, in general, significantly higher in Psa-inoculated plants (by at least 22%), comparing with control and Pfm-inoculated plants. Genes encoding antioxidant enzymes (SOD, APX and CAT) were upregulated at an earlier stage in Psa-inoculated A. arguta than in A. chinensis. In contrast, genes related with phenylpropanoids (LOX1) and ethylene (SAM) pathways were downregulated in A. arguta, but upregulated in A. chinensis in the later phases of infection. Expression of Pto3, responsible for pathogen recognition, occurred 2 dpi in A. arguta, but only 14 dpi in A. chinensis. In conclusion, we found that A. arguta is more tolerant to Psa and Pfm infection than A. chinensis and its primary and secondary metabolism is less impacted. A. arguta higher tolerance seems to be related with early pathogen recognition, the activation of plant antioxidant system, and to the suppression of ET and JA pathways from an earlier moment after infection.

Genome-wide identification and integrated analysis of lncRNAs in rice backcross introgression lines (BC2F12)

BACKGROUND

Distant hybridization is an important way to create interspecific genetic variation and breed new varieties in rice. A lot of backcross introgression lines (BILs) had been constructed for the scientific issues in rice. However, studies on the critical regulatory factor lncRNA in cultivated rice, wild rice and their BIL progenies were poorly reported.

RESULTS

Here, high-throughput RNA sequencing technology was used to explore the functional characteristics and differences of lncRNAs in O. sativa, O. longistaminata and their three BC₂F₁₂ progenies. A total of 1254 lncRNAs were screened out, and the number of differentially expressed lncRNAs between progenies and O. sativa were significantly less than that between progenies and O. longistaminata. Some lncRNAs regulated more than one mRNA, and 89.5% of lncRNAs regulated the expression of target genes through cis-acting. A total of 78 lncRNAs and 271 mRNAs were targeted by 280 miRNAs, and 22 lncRNAs were predicted to be the precursor of 20 microRNAs. Some miRNAs were found to target their own potential precursor lncRNAs. Over 50% of lncRNAs showed parental expression level dominance (ELD) in all three progenies, and most lncRNAs showed ELD-O. sativa rather than ELD-O. longistaminata. Further analysis showed that lncRNAs might regulate the expression of plant hormone-related genes and the adaptability of O. sativa, O. longistaminata and their progenies.

CONCLUSIONS

Taken together, the above results provided valuable clues for elucidating the functional features and expression differences of lncRNAs between O. sativa, O. longistaminata and their BIL progenies, and expanded our understanding about the biological functions of lncRNAs in rice.

Screening and characterization of long noncoding RNAs involved in the albinism of Ananas comosus var. bracteatus leaves

Long noncoding RNAs (lncRNAs) have been reported to play key regulatory roles in plant growth, development, and biotic and abiotic stress physiology. Revealing the mechanism of lncRNA regulation in the albino portions of leaves is important for understanding the development of chimeric leaves in Ananas comosus var. bracteatus. In this study, a total of 3,543 candidate lncRNAs were identified, among which 1,451 were differentially expressed between completely green (CGr) and completely white (CWh) leaves. LncRNAs tend to have shorter transcripts, lower expression levels, and greater expression specificity than protein-coding genes. Predicted lncRNA targets were functionally annotated by the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A lncRNA-mRNA interaction network was constructed, and 36 target mRNAs related to chlorophyll metabolism were predicted to interact with 86 lncRNAs. Among these, 25 significantly differentially expressed lncRNAs putatively interacted with 16 target mRNAs. Based on an expression pattern analysis of the lncRNAs and their target mRNAs, the lncRNAs targeting magnesium chelatase subunit H (ChlH), protochlorophyllide oxidoreductase (POR), and heme o synthase (COX10) were suggested as key regulators of chlorophyll metabolism. This study provides the first lncRNA database for A. comosus var. bracteatus and contributes greatly to understanding the mechanism of epigenetic regulation of leaf albinism.

A Systematic View Exploring the Role of Chloroplasts in Plant Abiotic Stress Responses

Chloroplasts are intracellular semiautonomous organelles central to photosynthesis and are essential for plant growth and yield. The significance of the function of chloroplast-related genes in response to climate change has not been well studied in crops. In the present study, the initial focus was on genes that were predicted to be located in the chloroplast genome in rice, a model crop plant, with genes either preferentially expressed in the leaf or ubiquitously expressed in all organs. The characteristics were analyzed by Gene Ontology (GO) enrichment and MapMan functional classification tools. It was then identified that 110 GO terms (45 for leaf expression and 65 for ubiquitous expression) and 1,695 genes mapped to MapMan overviews were strongly associated with chloroplasts. In particular, the MapMan cellular response overview revealed a close association between heat stress response and chloroplast-related genes in rice. Moreover, features of these genes in response to abiotic stress were analyzed using a large-scale publicly available transcript dataset. Consequently, the expression of 215 genes was found to be upregulated in response to high temperature stress. Conversely, genes that responded to other stresses were extremely limited. In other words, chloroplast-related genes were found to affect abiotic stress response mainly through high temperature response, with little effect on response to drought and salinity stress. These results suggest that genes involved in diurnal rhythm in the leaves participate in the reaction to recognize temperature changes in the environment. Furthermore, the predicted protein–protein interaction network analysis associated with high temperature stress is expected to provide a very important basis for the study of molecular mechanisms by which chloroplasts will respond to future climate changes.

Identification and functional prediction of cold-related long non-coding RNA (lncRNA) in grapevine

Plant long non-coding RNA (lncRNA) undergoes dynamic regulation and acts in developmental and stress regulation. In this study, we surveyed the expression dynamics of lncRNAs in grapevine (Vitis vinifera L.) under cold stress using high-throughput sequencing. Two-hundred and three known lncRNAs were significantly up-regulated and 144 known lncRNAs were significantly down-regulated in cold-treated grapevine. In addition, 2 088 novel lncRNA transcripts were identified in this study, with 284 novel lncRNAs significantly up-regulated and 182 novel lncRNAs significantly down-regulated in cold-treated grapevine. Two-hundred and forty-two differentially expressed grapevine lncRNAs were predicted to target 326 protein-coding genes in a cis-regulatory relationship. Many differentially expressed grapevine lncRNAs targeted stress response-related genes, such as CBF4 transcription factor genes, late embryogenesis abundant protein genes, peroxisome biogenesis protein genes, and WRKY transcription factor genes. Sixty-two differentially expressed grapevine lncRNAs were predicted to target 100 protein-coding genes in a trans-regulatory relationship. The expression of overall target genes in both cis and trans-regulatory relationships were positively related to the expression of lncRNAs in grapevines under cold stress. We identified 31 known lncRNAs as 34 grapevine micro RNA (miRNA) precursors and some miRNAs may be derived from multiple lncRNAs. We found 212 lncRNAs acting as targets of miRNAs in grapevines, involving 150 miRNAs; additionally, 120 grapevine genes were predicted as targets of grapevine miRNAs and lncRNAs. We found one gene cluster that was up-regulated and showed the same expression trend. In this cluster, many genes may be involved in abiotic stress response such as WRKY, Hsf, and NAC transcription factor genes.

Multiple quantitative trait loci contribute to resistance to bacterial canker incited by Pseudomonas syringae pv. actinidiae in kiwifruit (Actinidia chinensis)

Pseudomonas syringae pv. actinidiae (Psa) biovar 3, a virulent, canker-inducing pathogen is an economic threat to the kiwifruit (Actinidia spp.) industry worldwide. The commercially grown diploid (2×) A. chinensis var. chinensis is more susceptible to Psa than tetraploid and hexaploid kiwifruit. However information on the genetic loci modulating Psa resistance in kiwifruit is not available. Here we report mapping of quantitative trait loci (QTLs) regulating resistance to Psa in a diploid kiwifruit population, derived from a cross between an elite Psa-susceptible 'Hort16A' and a resistant male breeding parent P1. Using high-density genetic maps and intensive phenotyping, we identified a single QTL for Psa resistance on Linkage Group (LG) 27 of 'Hort16A' revealing 16-19% phenotypic variance and candidate alleles for susceptibility and resistance at this loci. In addition, six minor QTLs were identified in P1 on distinct LGs, exerting 4-9% variance. Resistance in the F1 population is improved by additive effects from 'Hort16A' and P1 QTLs providing evidence that divergent genetic pathways interact to combat the virulent Psa strain. Two different bioassays further identified new QTLs for tissue-specific responses to Psa. The genetic marker at LG27 QTL was further verified for association with Psa resistance in diploid Actinidia chinensis populations. Transcriptome analysis of Psa-resistant and susceptible genotypes in field revealed hallmarks of basal defense and provided candidate RNA-biomarkers for screening for Psa resistance in greenhouse conditions.

Genome-wide identification and characterization of mRNAs and lncRNAs involved in cold stress in the wild banana (Musa itinerans)

Cold stress seriously affects banana growth, yield and fruit quality. Long noncoding RNAs (lncRNAs) have been demonstrated as key regulators of biotic and abiotic stress in plants, but the identification and prediction of cold responsive mRNAs and lncRNAs in wild banana remains unexplored. In present study, a cold resistant wild banana line from China was used to profile the cold-responsive mRNAs and lncRNAs by RNA-seq under cold stress conditions, i.e. 13°C (critical growth temperature), 4°C (chilling temperature), 0°C (freezing temperature) and normal growing condition, i.e. 28°C (control group). A total of 12,462 lncRNAs were identified in cold-stressed wild banana. In mRNA, much more alternative splicing events occurred in wild banana under the cold stress conditions compared with that in the normal growing condition. The GO analysis of differential expression genes (DEGs) showed the biochemical processes and membrane related genes responded positively to the cold stress. The KEGG pathway enrichment analysis of the DEGs showed that the pathways of photosynthesis, photosynthesis–antenna proteins, circadian rhythm–plant, glutathione metabolism, starch and sucrose metabolism, cutin/suberine/biosynthesis were altered or affected by the cold stress conditions. Our analyses of the generated transcriptome and lncRNAs provide new insights into regulating expression of genes and lncRNAs that respond to cold stress in the wild banana.