Genome-wide transcriptional changes of ramie (Boehmeria nivea L. Gaud) in response to root-lesion nematode infection

A Meta-analysis of Interactions Between Insect Herbivores and Plant-Parasitic Nematodes

Insect herbivores and plant-parasitic nematodes are global, economically devastating pests that are present in nearly every crop and natural system worldwide. Although they may be spatially separated, they indirectly interact with each other by altering both plant chemical defense and nutrition. However, the outcome of these interactions is highly variable across different focal species. We performed a meta-analysis to determine how plant and nematode traits influence insect herbivore growth and reproduction, as well as nematode abundance and reproduction. We investigated how interactions between plant-parasitic nematodes and insect herbivores influence plant biomass, carbon, and nitrogen in the roots and shoots. We found no overall effect of nematodes on insect herbivores or insect herbivores on nematodes. However, while phloem-feeding insect reproduction was not affected by nematode feeding guild or plant family, chewing insect growth increased in the presence of cyst nematodes and decreased in the presence of gall nematodes. The effect of nematodes on chewing insect herbivore growth was also affected by the focal plant family. Nematode presence did not alter plant biomass when plants were exposed to aboveground insect herbivory, but carbon and nitrogen were higher in roots and nitrogen was higher in shoots of plants with nematodes and insects compared to plants with insects alone. Our results indicate that the mechanisms driving the outcome of aboveground-belowground interactions are still unclear, but those chewing insects may have more variable responses to nematode damage than phloem-feeders.

Comparative transcriptome analysis of R3a and Avr3a-mediated defense responses in transgenic tomato

Late blight caused by Phytophthora infestans is one of the most devastating diseases in potatoes and tomatoes. At present, several late blight resistance genes have been mapped and cloned. To better understand the transcriptome changes during the incompatible interaction process between R3a and Avr3a, in this study, after spraying DEX, the leaves of MM-R3a-Avr3a and MM-Avr3a transgenic plants at different time points were used for comparative transcriptome analysis. A total of 7,324 repeated DEGs were detected in MM-R3a-Avr3a plants at 2-h and 6-h, and 729 genes were differentially expressed at 6-h compared with 2-h. Only 1,319 repeated DEGs were found in MM-Avr3a at 2-h and 6-h, of which 330 genes have the same expression pattern. Based on GO, KEGG and WCGNA analysis of DEGs, the phenylpropanoid biosynthesis, plant-pathogen interaction, and plant hormone signal transduction pathways were significantly up-regulated. Parts of the down-regulated DEGs were enriched in carbon metabolism and the photosynthesis process. Among these DEGs, most of the transcription factors, such as WRKY, MYB, and NAC, related to disease resistance or endogenous hormones SA and ET pathways, as well as PR, CML, SGT1 gene were also significantly induced. Our results provide transcriptome-wide insights into R3a and Avr3a-mediated incompatibility interaction.

Transcriptome Analysis of High-NUE (T29) and Low-NUE (T13) Genotypes Identified Different Responsive Patterns Involved in Nitrogen Stress in Ramie (Boehmeria nivea (L.) Gaudich)

Nitrogen-use efficiency (NUE) has significant impacts on plant growth and development. NUE in plants differs substantially in physiological resilience to nitrogen stress; however, the molecular mechanisms underlying enhanced resilience of high-NUE plants to nitrogen deficiency remains unclear. We compared transcriptome-wide gene expression between high-NUE and low-NUE ramie (Boehmeria nivea (L.) Gaudich) genotypes under nitrogen (N)-deficient and normal conditions to identify the transcriptomic expression patterns that contribute to ramie resilience to nitrogen deficiency. Two ramie genotypes with contrasting NUE were used in the study, including T29 (NUE = 46.01%) and T13 (NUE = 15.81%). Our results showed that high-NUE genotypes had higher gene expression under the control condition across 94 genes, including frontloaded genes such as GDSL esterase and lipase, gibberellin, UDP-glycosyltransferase, and omega-6 fatty acid desaturase. Seventeen stress-tolerance genes showed lower expression levels and varied little in response to N-deficiency stress in high-NUE genotypes. In contrast, 170 genes were upregulated under N deficiency in high-NUE genotypes but downregulated in low-NUE genotypes compared with the controls. Furthermore, we identified the potential key genes that enable ramie to maintain physiological resilience under N-deficiency stress, and categorized these genes into three groups based on the transcriptome and their expression patterns. The transcriptomic and clustering analysis of these nitrogen-utilization-related genes could provide insight to better understand the mechanism of linking among the three gene classes that enhance resilience in high-NUE ramie genotypes.

Transcriptome analysis of near-isogenic lines for glume hairiness of wheat

Hairiness, which is a phenotypic trait common among land plants, primarily affects the stem, leaf, and floral organs. Plant hairiness is associated with complex functions. For example, glume hairiness in wheat is related to the resistance to biotic and abiotic stresses, and may also influence human health. In the present study, two pairs of near-isogenic lines (NILs) for glume hairiness, which were derived from a cross between a Tibetan semi-wild wheat accession (Triticum aestivum ssp. tibetanum Q1028) and a common wheat cultivar (T. aestivum 'Zhengmai 9023'), underwent a glume transcriptome analysis. We detected 27,935 novel genes, of which 18,027 were annotated. Additionally, 488 and 600 differentially expressed genes (DEGs) were detected in NIL1 and NIL2, respectively, with 37 DEGs detected in both NIL pairs. Moreover, 987 and 1584 single nucleotide polymorphisms (SNPs) were detected in NIL1 and NIL2, respectively, with 39 SNPs detected in both NIL pairs, of which most were located in the Hairy glume (Hg) gene region on chromosome arm 1AS. The annotation of the DEGs with gene ontology terms revealed that genes associated with hairiness in Arabidopsis and rice were similarly enriched. The possible functions of these genes related to glume hairiness were examined. The study results provide useful information for identifying candidate genes and the fine-mapping of Hg in the wheat genome.

Selection of Reference Genes for qPCR Analyses of Gene Expression in Ramie Leaves and Roots across Eleven Abiotic/Biotic Treatments

Quantitative real-time PCR (qPCR) is commonly used for deciphering gene functions. For effective qPCR analyses, suitable reference genes are needed for normalization. The objective of this study is to identify the appropriate reference gene(s) for qPCR analyses of the leaves and roots of ramie (Boehmeria nivea L.), an important natural fiber crop. To accomplish this goal, we investigated the expression patterns of eight common plant qPCR reference genes in ramie leaves and roots under five abiotic stresses, five hormonal treatments, and one biotic stress. The relative expression stabilities of the eight genes were evaluated using four common but different approaches: geNorm, NormFinder, BestKeeper, and RefFinder. Across the 11 tested conditions, ACT1 was the most stably expressed among the eight genes while GAPDH displayed the biggest variation. Overall, while variations in the suggested reference genes were found for different tissue x treatment combinations, our analyses revealed that together, genes ACT1, CYP2, and UBQ can provide robust references for gene expression studies of ramie leaves under most conditions, while genes EF-1α, TUB, and ACT1 can be used for similar studies of ramie roots. Our results should help future functional studies of the genes in ramie genome across tissues and environmental conditions.

A genomic resource derived from the integration of genome sequences, expressed transcripts and genetic markers in ramie

Background

The redundancy of genomic resources, including transcript and molecular markers, and their uncertain position in the genome have dramatically hindered the study of traits in ramie, an important natural fiber crop.

Results

We obtained a high-quality transcriptome consisting of 30,591 non-redundant transcripts using single-molecule long-read sequencing and proposed it as a universal ramie transcriptome. Additionally, 55,882 single nucleotide polymorphisms (SNPs) were identified and a high-density genetic map was developed. Based on this genetic map, 181.7 Mb ramie genome sequences were assembled into 14 chromosomes. For the convenient use of these resources, 29,286 (~ 95.7%) of the transcripts and all 55,882 SNPs, along with 1827 previously reported sequence repeat markers (SSRs), were mapped into the ramie genome, and 22,343 (~ 73.0%) transcripts, 50,154 (~ 89.7%) SNPs, and 1466 (~ 80.3%) SSRs were assigned to a specific location in the corresponding chromosome.

Conclusion

This is the first study to characterize the ramie transcriptome by long-read sequencing, and the substantial number of transcripts of significant length obtained will accelerate our understanding of ramie growth and development. This integration of genome sequences, expressed transcripts, and genetic markers will provide an extremely useful resource for genetic, molecular, and breeding studies of ramie.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5878-8) contains supplementary material, which is available to authorized users.

Comparative transcriptomics reveals the selection patterns of domesticated ramie

Although domestication has dramatically altered the phenotype, physiology, and life history of ramie (Boehmeria nivea) plants, few studies have investigated the effects of domestication on the structure and expression pattern of genes in this fiber crop. To investigate the selective pattern and genetic relationships among a cultivated variety of ramie (BNZ: B. nivea, ZZ1) and four wild species, BNT (B. nivea var. tenacissima), BNN (B. nivea var. nipononivea), BNW (B. nivea var. nivea), and BAN (B. nivea var. viridula), in the section Tilocnide, we performed an RNA sequencing analysis of these ramie species. The de novo assembly of the "all-ramie" transcriptome yielded 119,114 unigenes with an average length of 633 bp, and a total of 7,084 orthologous gene pairs were identified. The phylogenetic tree showed that the cultivar BNZ clustered with BAN in one group, BNW was closely related to BNT, and BNN formed a separate group. Introgression analysis indicated that gene flow occurred from BNZ to BNN and BAN, and between BAN and BNN. Among these orthologs, 2,425 and 269 genes underwent significant purifying and positive selection, respectively. For these positively selected genes, oxidation-reduction process (GO:0055114) and stress response pathways (GO:0006950) were enriched, indicating that modulation of the cellular redox status was important during both ramie fiber evolution and improvement. Two genes related to the suppression of flowering and one gene annotated as a flowering-promoting factor were subjected to positive selection, probably caused by human manipulation. Additionally, five genes were homologs of those involved in abiotic stress tolerance and disease resistance, with higher expression levels in the cultivar BNZ than in the wild species. Collectively, the results of this study indicated that domestication has resulted in the upregulation of many genes involved in the abiotic and biotic stress responses, fiber yield, and plant growth of ramie.

Identification of candidate effector genes of Pratylenchus penetrans

Pratylenchus penetrans is one of the most important species of root lesion nematodes (RLNs) because of its detrimental and economic impact in a wide range of crops. Similar to other plant-parasitic nematodes (PPNs), P. penetrans harbours a significant number of secreted proteins that play key roles during parasitism. Here, we combined spatially and temporally resolved next-generation sequencing datasets of P. penetrans to select a list of candidate genes aimed at the identification of a panel of effector genes for this species. We determined the spatial expression of transcripts of 22 candidate effectors within the oesophageal glands of P. penetrans by in situ hybridization. These comprised homologues of known effectors of other PPNs with diverse putative functions, as well as novel pioneer effectors specific to RLNs. It is noteworthy that five of the pioneer effectors encode extremely proline-rich proteins. We then combined in situ localization of effectors with available genomic data to identify a non-coding motif enriched in promoter regions of a subset of P. penetrans effectors, and thus a putative hallmark of spatial expression. Expression profiling analyses of a subset of candidate effectors confirmed their expression during plant infection. Our current results provide the most comprehensive panel of effectors found for RLNs. Considering the damage caused by P. penetrans, this information provides valuable data to elucidate the mode of parasitism of this nematode and offers useful suggestions regarding the potential use of P. penetrans-specific target effector genes to control this important pathogen.

QTL analysis of four main stem bark traits using a GBS-SNP-based high-density genetic map in ramie

Ramie fiber extracted from ramie stem bark (RSB) is a highly important natural fiber, and therefore, RSB is an economically important plant organ. The genetic basis of RSB traits is poorly understood. In the present study, fiber yield and three RSB traits (bark thickness, bark weight, and fiber output ratio) were subject to quantitative trait locus (QTL) analysis using an F₂ agamous line population derived from two ramie varieties (Qingdaye and Zhongzhu 1). A total of 4338 high-quality single nucleotide polymorphisms were identified using the genotyping-by-sequencing technique and were subsequently used to construct a high-density genetic map spanning 1942.9 cM. Thereafter, QTL analysis identified five, two, four, and four QTLs for bark thickness, bark weight, fiber output ratio, and fiber yield, respectively. A 5.1 cM region that corresponded to a QTL for bark thickness (qBT4a) contained 106 candidate genes, and the Zhongzhu 1 allele of one of the genes, a putative MYB gene (evm.

MODEL

scaffold7373.133_D1), included a 760-bp insertion that caused premature termination, thereby producing a protein that lacked part of the MYB domain. Because MYB transcription factors play central roles in regulating the development of secondary cellular walls and fiber biosynthesis, we propose evm.

MODEL

scaffold7373.133_D1 as a likely candidate gene for qBT4a.

Biodegradation of ramie stalk by Flammulina velutipes: mushroom production and substrate utilization

In the textile industry, ramie stalk is byproducts with a low economic value. The potential use of this leftover as a substrate ingredient for Flammulina velutipes (F. velutipe) cultivation was evaluated. The degradation and utilization of ramie stalk by F. velutipes was evaluated through mushroom production, lignocelluloses degradation and lignocellulolytic enzymes activity. The best substrate mixture for F. velutipes cultivation comprised 50% ramie stalk, 20% cottonseed hulls, 25% wheat bran, 4% cornstarch and 2% CaCO₃. The highest biological efficiency of fruiting bodies was reached 119.7%. F. velutipes appears to degrade 12.7–32.0% lignin, 14.4–30.2% cellulose and 9.3–25.7% hemicellulose during cultivation on the different substrates. The results of enzymes activities showed that laccase and peroxidase were higher before fruiting; while cellulase and hemicellulase showed higher activities after fruiting. The biological efficiency of fruiting bodies was positively correlated with the activities of cellulase, hemicellulase and ligninolytic enzyme. The results of this study demonstrate that ramie stalk can be used as an effective supplement for increasing mushroom yield in F. velutipes.

Characterization of Lilium longiflorum cv. 'Nellie White' Infection with Root-lesion Nematode Pratylenchus penetrans by Bright-field and Transmission Electron Microscopy

Lilium longiflorum cv. Nellie White, commonly known as Easter lily, is an important floral crop with an annual wholesale value of over $26 million in the United States. The root-lesion nematode, Pratylenchus penetrans, is a major pest of lily due to the significant root damage it causes. In this study, we investigated the cytological aspects of this plant-nematode interaction using bright-field and transmission electron microscopy. We took advantage of an in vitro culture method to multiply lilies and follow the nematode infection over time. Phenotypic reactions of roots inoculated with P. penetrans were evaluated from 0 to 60 d after nematode infection. Symptom development progressed from initial randomly distributed discrete necrotic areas to advanced necrosis along entire roots of each inoculated plant. A major feature characterizing this susceptible host response to nematode infection was the formation of necrosis, browning, and tissue death involving both root epidermis and cortical cells. Degradation of consecutive cell walls resulted in loss of cell pressure, lack of cytoplasmic integrity, followed by cell death along the intracellular path of the nematode's migration. Pratylenchus penetrans was never seen in the vascular cylinder as the layer of collapsed endodermal cells presumably blocked the progression of nematodes into this area of the roots. This study presents the first detailed cytological characterization of P. penetrans infection of Easter lily plants.

Transcriptome Analysis of the Cf-12-Mediated Resistance Response to Cladosporium fulvum in Tomato

Cf-12 is an effective gene for resisting tomato leaf mold disease caused by Cladosporium fulvum (C. fulvum). Unlike many other Cf genes such as Cf-2, Cf-4, Cf-5, and Cf-9, no physiological races of C. fulvum that are virulent to Cf-12 carrying plant lines have been identified. In order to better understand the molecular mechanism of Cf-12 gene resistance response, RNA-Seq was used to analyze the transcriptome changes at three different stages of C. fulvum infection (0, 4, and 8 days post infection [dpi]). A total of 9100 differentially expressed genes (DEGs) between 4 and 0 dpi, 8643 DEGs between 8 and 0 dpi and 2547 DEGs between 8 and 4 dpi were identified. In addition, we found that 736 DEGs shared among the above three groups, suggesting the presence of a common core of DEGs in response to C. fulvum infection. These DEGs were significantly enriched in defense-signaling pathways such as the calcium dependent protein kinases pathway and the jasmonic acid signaling pathway. Additionally, we found that many transcription factor genes were among the DEGs, indicating that transcription factors play an important role in C. fulvum defense response. Our study provides new insight on the molecular mechanism of Cf resistance to C. fulvum, especially the unique features of Cf-12 in responding to C. fulvum infection.

Advances in Understanding the Molecular Mechanisms of Root Lesion Nematode Host Interactions

Root lesion nematodes (RLNs) are one of the most economically important groups of plant nematodes. As migratory endoparasites, their presence in roots is less obvious than infestations of sedentary endoparasites; nevertheless, in many instances, they are the major crop pests. With increasing molecular information on nematode parasitism, available data now reflect the differences and, in particular, similarities in lifestyle between migratory and sedentary endoparasites. Far from being unsophisticated compared with sedentary endoparasites, migratory endoparasites are exquisitely suited to their parasitic lifestyle. What they lack in effectors required for induction of permanent feeding sites, they make up for with their versatile host range and their ability to move and feed from new host roots and survive adverse conditions. In this review, we summarize the current molecular data available for RLNs and highlight differences and similarities in effectors and molecular mechanisms between migratory and sedentary endoparasitic nematodes.

Transcriptome Analysis of Ramie (Boehmeria nivea L. Gaud.) in Response to Ramie Moth (Cocytodes coerulea Guenée) Infestation

The ramie moth Cocytodes coerulea Guenée (RM) is an economically important pest that seriously impairs the yield of ramie, an important natural fiber crop. The molecular mechanisms that underlie the ramie-pest interactions are unclear up to date. Therefore, a transcriptome profiling analysis would aid in understanding the ramie defense mechanisms against RM. In this study, we first constructed two cDNA libraries derived from RM-challenged (CH) and unchallenged (CK) ramie leaves. The subsequent sequencing of the CH and CK libraries yielded 40.2 and 62.8 million reads, respectively. Furthermore, de novo assembling of these reads generated 26,759 and 29,988 unigenes, respectively. An integrated assembly of data from these two libraries resulted in 46,533 unigenes, with an average length of 845 bp per unigene. Among these genes, 24,327 (52.28%) were functionally annotated by predicted protein function. A comparative analysis of the CK and CH transcriptome profiles revealed 1,980 differentially expressed genes (DEGs), of which 750 were upregulated and 1,230 were downregulated. A quantitative real-time PCR (qRT-PCR) analysis of 13 random selected genes confirmed the gene expression patterns that were determined by Illumina sequencing. Among the DEGs, the expression patterns of transcription factors, protease inhibitors, and antioxidant enzymes were studied. Overall, these results provide useful insights into the defense mechanism of ramie against RM.

Transcriptome Sequencing in Response to Salicylic Acid in Salvia miltiorrhiza

Salvia miltiorrhiza is a traditional Chinese herbal medicine, whose quality and yield are often affected by diseases and environmental stresses during its growing season. Salicylic acid (SA) plays a significant role in plants responding to biotic and abiotic stresses, but the involved regulatory factors and their signaling mechanisms are largely unknown. In order to identify the genes involved in SA signaling, the RNA sequencing (RNA-seq) strategy was employed to evaluate the transcriptional profiles in S. miltiorrhiza cell cultures. A total of 50,778 unigenes were assembled, in which 5,316 unigenes were differentially expressed among 0-, 2-, and 8-h SA induction. The up-regulated genes were mainly involved in stimulus response and multi-organism process. A core set of candidate novel genes coding SA signaling component proteins was identified. Many transcription factors (e.g., WRKY, bHLH and GRAS) and genes involved in hormone signal transduction were differentially expressed in response to SA induction. Detailed analysis revealed that genes associated with defense signaling, such as antioxidant system genes, cytochrome P450s and ATP-binding cassette transporters, were significantly overexpressed, which can be used as genetic tools to investigate disease resistance. Our transcriptome analysis will help understand SA signaling and its mechanism of defense systems in S. miltiorrhiza.

Identification of Ramie Genes in Response to Pratylenchus coffeae Infection Challenge by Digital Gene Expression Analysis

Root lesion disease, caused by Pratylenchus coffeae, seriously impairs the growth and yield of ramie, an important natural fiber crop. The ramie defense mechanism against P. coffeae infection is poorly understood, which hinders efforts to improve resistance via breeding programs. In this study, the transcriptome of the resistant ramie cultivar Qingdaye was characterized using Illumina sequence technology. About 46.3 million clean pair end (PE) reads were generated and assembled into 40,826 unigenes with a mean length of 830 bp. Digital gene expression (DGE) analysis was performed on both the control roots (CK) and P. coffeae-challenged roots (CH), and the differentially expressed genes (DEGs) were identified. Approximately 10.16 and 8.07 million cDNA reads in the CK and CH cDNA libraries were sequenced, respectively. A total of 137 genes exhibited different transcript abundances between the two libraries. Among them, the expressions of 117 and 20 DEGs were up- and down-regulated in P. coffeae-challenged ramie, respectively. The expression patterns of 15 candidate genes determined by qRT-PCR confirmed the results of DGE analysis. Time-course expression profiles of eight defense-related genes in susceptible and resistant ramie cultivars were different after P. coffeae inoculation. The differential expression of protease inhibitors, pathogenesis-related proteins (PRs), and transcription factors in resistant and susceptible ramie during P. coffeae infection indicated that cystatin likely plays an important role in nematode resistance.

Transcriptome comparison reveals the patterns of selection in domesticated and wild ramie (Boehmeria nivea L. Gaud)

Ramie is an old fiber crop, cultivated for thousands of years in China. The cultivar ramie evolved from the wild species Qingyezhuma (QYZM, Boehmeria nivea var. tenacissima). However, the mechanism of domestication of this old fiber crop is poorly understood. In order to characterize the selective pattern in ramie domestication, orthologous genes between the transcriptomes of domesticated ramie variety Zhongzhu 1 (ZZ1) and wild QYZM were assessed using bidirectional best-hit method and ratio of non-synonymous (Ka) to synonymous (Ks) nucleotide substitutions was estimated. Sequence comparison of 56,932 and 59,246 unigenes from the wild QYZM and domesticated ZZ1, respectively, helped identify 10,745 orthologous unigene pairs with a total orthologous length of 10.18 Mb. Among these unigenes, 85 and 13 genes were found to undergo significant purifying and positive selection, respectively. Most of the selected genes were homologs of those involved in abiotic stress tolerance or disease resistance in other plants, suggesting that abiotic and biotic stresses were important selective pressures in ramie domestication. Two genes probably related to the fiber yield of ramie were subjected to positive selection, which may be caused by human manipulation. Thus, our results show the pervasive effects of artificial and natural selections on the accelerated domestication of ramie from its wild relative.

Identification of 32 full-length NAC transcription factors in ramie (Boehmeria nivea L. Gaud) and characterization of the expression pattern of these genes

NAM, ATAF, and CUC (NAC) genes are plant-specific transcription factors (TFs) that play key roles in plant growth, development, and stress tolerance. To date, none of the ramie NAC (BnNAC) genes had been identified, even though ramie is one of the most important natural fiber crops. In order to mine the BnNAC TFs and identify their potential function, the search for BnNAC genes against two pools of unigenes de novo assembled from the RNA-seq in our two previous studies was performed, and a total of 32 full-length BnNAC genes were identified in this study. Forty-seven function-known NAC proteins published in other species, in concert with these 32 BnNAC proteins were subjected to phylogenetic analysis, and the result showed that all the 79 NAC proteins can be divided into eight groups (NAC-I-VIII). Among the 32 BnNAC genes, 24, 2, and 1 gene showed higher expression in stem xylem, leaf, and flower, respectively. Furthermore, the expression of 14, 11 and 4 BnNAC genes was regulated by drought, cadmium stress, and infection by root lesion nematode, respectively. Interestingly, there were five BnNAC TFs which showed high homology with the NAC TFs of other species involved in regulating the secondary wall synthesis, and their expressions were not regulated by drought and cadmium stress. These results suggested that the BnNAC family might have a functional diversity. The identification of these 32 full-length BnNAC genes and the characterization of their expression pattern provide a basis for future clarification of their functions in ramie growth and development.