De novo foliar transcriptome of Chenopodium amaranticolor and analysis of its gene expression during virus-induced hypersensitive response

Identification and verification of seed development related miRNAs in kernel almond by small RNA sequencing and qPCR

Many studies have investigated the role of miRNAs on the yield of various plants, but so far, no report is available on the identification and role of miRNAs in fruit and seed development of almonds. In this study, preliminary analysis by high-throughput sequencing of short RNAs of kernels from the crosses between almond cultivars 'Sefid' × 'Mamaee' (with small and large kernels, respectively) and 'Sefid' × 'P. orientalis' (with small kernels) showed that the expressions of several miRNAs such as Pdu-miR395a-3p, Pdu-miR8123-5p, Pdu-miR482f, Pdu-miR6285, and Pdu-miR396a were significantly different. These miRNAs targeted genes encoding different proteins such as NYFB-3, SPX1, PGSIP3 (GUX2), GH3.9, and BEN1. The result of RT-qPCR revealed that the expression of these genes showed significant differences between the crosses and developmental stages of the seeds, suggesting that these genes might be involved in controlling kernel size because the presence of these miRNAs had a negative effect on their target genes. Pollen source can influence kernel size by affecting hormonal signaling and metabolic pathways through related miRNAs, a phenomenon known as xenia.

Transcriptional Regulatory Networks Associate with Early Stages of Potato Virus X Infection of Solanum tuberosum

Potato virus X (PVX) belongs to genus Potexvirus. This study characterizes the cellular transcriptome responses to PVX infection in Russet potato at 2 and 3 days post infection (dpi). Among the 1242 differentially expressed genes (DEGs), 268 genes were upregulated, and 37 genes were downregulated at 2 dpi while 677 genes were upregulated, and 265 genes were downregulated at 3 dpi. DEGs related to signal transduction, stress response, and redox processes. Key stress related transcription factors were identified. Twenty-five pathogen resistance gene analogs linked to effector triggered immunity or pathogen-associated molecular pattern (PAMP)-triggered immunity were identified. Comparative analysis with Arabidopsis unfolded protein response (UPR) induced DEGs revealed genes associated with UPR and plasmodesmata transport that are likely needed to establish infection. In conclusion, this study provides an insight on major transcriptional regulatory networked involved in early response to PVX infection and establishment.

Response of bitter and sweet Chenopodium quinoa varieties to cucumber mosaic virus: Transcriptome and small RNASeq perspective

Saponins are secondary metabolites with antiviral properties. Low saponin (sweet) varieties of quinoa (Chenopodium quinoa) have been developed because seeds high in saponins taste bitter. The aim of this study was to elucidate the role of saponin in resistance of quinoa to Cucumber mosaic virus (CMV). Differential gene expression was studied in time-series study of CMV infection. High-throughput transcriptome sequence data were obtained from 36 samples (3 varieties × +/- CMV × 1 or 4 days after inoculation × 3 replicates). Translation, lipid, nitrogen, amino acid metabolism, and mono- and sesquiterpenoid biosynthesis genes were upregulated in CMV infections. In 'Red Head' (bitter), CMV-induced systemic symptoms were concurrent with downregulation of a key saponin biosynthesis gene, TSARL1, four days after inoculation. In local lesion responses (sweet and semi-sweet), TSARL1 levels remained up-regulated. Known microRNAs (miRNA) (81) from 11 miR families and 876 predicted novel miRNAs were identified. Differentially expressed miRNA and short interfering RNA clusters (24nt) induced by CMV infection are predicted to target genomic and intergenic regions enriched in repetitive elements. This is the first report of integrated RNASeq and sRNASeq data in quinoa-virus interactions and provides comprehensive understanding of involved genes, non-coding regions, and biological pathways in virus resistance.

Sequencing the Plastid Genome of Giant Ragweed (Ambrosia trifida, Asteraceae) From a Herbarium Specimen

We report the first plastome sequence of giant ragweed (Ambrosia trifida); with this new genome information, we assessed the phylogeny of Asteraceae and the transcriptional profiling against glyphosate resistance in giant ragweed. Assembly and genic features show a normal angiosperm quadripartite plastome structure with no signatures of deviation in gene directionality. Comparative analysis revealed large inversions across the plastome of giant ragweed and the previously sequenced members of the plant family. Asteraceae plastid genomes contain two inversions of 22.8 and 3.3 kb; the former is located between trnS-GCU and trnG-UCC genes, and the latter between trnE-UUC and trnT-GGU genes. The plastid genome sequences of A. trifida and the related species, Ambrosia artemisiifolia, are identical in gene content and arrangement, but they differ in length. The phylogeny is well-resolved and congruent with previous hypotheses about the phylogenetic relationship of Asteraceae. Transcriptomic analysis revealed divergence in the relative expressions at the exonic and intronic levels, providing hints toward the ecological adaptation of the genus. Giant ragweed shows various levels of glyphosate resistance, with introns displaying higher expression patterns at resistant time points after the assumed herbicide treatment.

Characterization of the genome of a phylogenetically distinct tospovirus and its interactions with the local lesion-induced host Chenopodium quinoa by whole-transcriptome analyses

Chenopodium quinoa is a natural local lesion host of numerous plant viruses, including tospoviruses (family Bunyaviridae). Groundnut chlorotic fan-spot tospovirus (GCFSV) has been shown to consistently induce local lesions on the leaves of C. quinoa 4 days post-inoculation (dpi). To reveal the whole genome of GCFSV and its interactions with C. quinoa, RNA-seq was performed to determine the transcriptome profiles of C. quinoa leaves. The high-throughput reads from infected C. quinoa leaves were used to identify the whole genome sequence of GCFSV and its single nucleotide polymorphisms. Our results indicated that GCFSV is a phylogenetically distinct tospovirus. Moreover, 27,170 coding and 29,563 non-coding sequences of C. quinoa were identified through de novo assembly, mixing reads from mock and infected samples. Several key genes involved in the modulation of hypersensitive response (HR) were identified. The expression levels of 4,893 deduced complete genes annotated using the Arabidopsis genome indicated that several HR-related orthologues of pathogenesis-related proteins, transcription factors, mitogen-activated protein kinases, and defense proteins were significantly expressed in leaves that formed local lesions. Here, we also provide new insights into the replication progression of a tospovirus and the molecular regulation of the C. quinoa response to virus infection.

Making Plants Break a Sweat: the Structure, Function, and Evolution of Plant Salt Glands

Salt stress is a complex trait that poses a grand challenge in developing new crops better adapted to saline environments. Some plants, called recretohalophytes, that have naturally evolved to secrete excess salts through salt glands, offer an underexplored genetic resource for examining how plant development, anatomy, and physiology integrate to prevent excess salt from building up to toxic levels in plant tissue. In this review we examine the structure and evolution of salt glands, salt gland-specific gene expression, and the possibility that all salt glands have originated via evolutionary modifications of trichomes. Salt secretion via salt glands is found in more than 50 species in 14 angiosperm families distributed in caryophyllales, asterids, rosids, and grasses. The salt glands of these distantly related clades can be grouped into four structural classes. Although salt glands appear to have originated independently at least 12 times, they share convergently evolved features that facilitate salt compartmentalization and excretion. We review the structural diversity and evolution of salt glands, major transporters and proteins associated with salt transport and secretion in halophytes, salt gland relevant gene expression regulation, and the prospect for using new genomic and transcriptomic tools in combination with information from model organisms to better understand how salt glands contribute to salt tolerance. Finally, we consider the prospects for using this knowledge to engineer salt glands to increase salt tolerance in model species, and ultimately in crops.

Transcriptional Responses of Chilean Quinoa (Chenopodium quinoa Willd.) Under Water Deficit Conditions Uncovers ABA-Independent Expression Patterns

HIGHLIGHTS R49 genotype displayed best performance on selected physiological parameters and highest tolerance to drought.R49 drought over-represented transcripts has exhibited 19% of genes (306 contigs) that presented no homology to published databases.Expression pattern for canonical responses to drought such as ABA biosynthesis and other genes induced in response to drought were assessed by qPCR. Global freshwater shortage is one of the biggest challenges of our time, often associated to misuse, increased consumption demands and the effects of climate change, paralleled with the desertification of vast areas. Chenopodium quinoa (Willd.) represents a very promising species, due to both nutritional content and cultivation under water constraint. We characterized drought tolerance of three Chilean genotypes and selected Genotype R49 (Salares ecotype) based upon Relative Water Content (RWC), Electrolyte Leakage (EL) and maximum efficiency of photosystem II (Fv/Fm) after drought treatment, when compared to another two genotypes. Exploratory RNA-Seq of R49 was generated by Illumina paired-ends method comparing drought and control irrigation conditions. We obtained 104.8 million reads, with 54 million reads for control condition and 51 million reads for drought condition. Reads were assembled in 150,952 contigs, were 31,523 contigs have a reading frame of at least 300 nucleotides (100 aminoacids). BLAST2GO annotation showed a 15% of genes without homology to NCBI proteins, but increased to 19% (306 contigs) when focused into drought-induced genes. Expression pattern for canonical drought responses such as ABA biosynthesis and other genes induced were assessed by qPCR, suggesting novelty of R49 drought responses.

Transcriptome analysis of woodland strawberry (Fragaria vesca) response to the infection by Strawberry vein banding virus (SVBV)

BACKGROUND

Woodland strawberry (Fragaria vesca) infected with Strawberry vein banding virus (SVBV) exhibits chlorotic symptoms along the leaf veins. However, little is known about the molecular mechanism of strawberry disease caused by SVBV.

METHODS

We performed the next-generation sequencing (RNA-Seq) study to identify gene expression changes induced by SVBV in woodland strawberry using mock-inoculated plants as a control.

RESULTS

Using RNA-Seq, we have identified 36,850 unigenes, of which 517 were differentially expressed in the virus-infected plants (DEGs). The unigenes were annotated and classified with Gene Ontology (GO), Clusters of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The KEGG pathway analysis of these genes suggested that strawberry disease caused by SVBV may affect multiple processes including pigment metabolism, photosynthesis and plant-pathogen interactions.

CONCLUSIONS

Our research provides comprehensive transcriptome information regarding SVBV infection in strawberry.

Differentially expressed genes of Chenopodium amaranticolor in response to cymbidium mosaic virus infection

Cymbidium mosaic virus (CymMV)-induced expressed sequence tag (EST) clones from Chenopodium amaranticolor were identified. CymMV was mechanically inoculated onto C. amaranticolor, and local lesion symptoms were observed. Inoculated leaves were collected on serial days post inoculation (dpi) to identify activated or suppressed genes. mRNA isolation and suppression subtractive hybridization (SSH) were then performed to identify differentially expressed genes related to the local lesion response. Fifty-three ESTs, including genes related to defense and stress responses (e.g., lipoxygenase, jasmonate-induced protein, and heat shock protein), were generated. In addition, a large proportion of the ESTs were found to be involved in photosynthesis, as determined by their functional categories. Expression levels of several EST genes were observed using quantitative real-time reverse transcription-polymerase chain reaction, and the evaluated genes showed varying levels of expression during the experimental period. In this study, differentially expressed sequences via SSH were identified from CymMV-infected C. amaranticolor, and profiling and annotation were carried out to determine the expression pattern of CymMV and its interaction with C. amaranticolor.

Gene Expression Changes during the Gummosis Development of Peach Shoots in Response to Lasiodiplodia theobromae Infection Using RNA-Seq

Lasiodiplodia theobromae is a causal agent of peach (Prunus persica L.) tree gummosis, a serious disease affecting peach cultivation and production. However, the molecular mechanism underlying the pathogenesis remains unclear. RNA-Seq was performed to investigate gene expression in peach shoots inoculated or mock-inoculated with L. theobromae. A total of 20772 genes were detected in eight samples; 4231, 3750, 3453, and 3612 differentially expressed genes were identified at 12, 24, 48, and 60 h after inoculation, respectively. Furthermore, 920 differentially co-expressed genes (515 upregulated and 405 downregulated) were found, respectively. Gene ontology annotation revealed that phenylpropanoid biosynthesis and metabolism, uridine diphosphate-glucosyltransferase activity, and photosynthesis were the most differentially regulated processes during gummosis development. Significant differences were also found in the expression of genes involved in glycometabolism and in ethylene and jasmonic acid biosynthesis and signaling. These data illustrate the dynamic changes in gene expression in the inoculated peach shoots at the transcriptome level. Overall, gene expression in defense response and glycometabolism might result in the gummosis of peach trees induced by L. theobromae.

Digital Gene Expression Analysis to Screen Disease Resistance-Relevant Genes from Leaves of Herbaceous Peony (Paeonia lactiflora Pall.) Infected by Botrytis cinerea

Herbaceous peony (Paeonia lactiflora Pall.) is a well-known traditional flower in China and is widely used for landscaping and garden greening due to its high ornamental value. However, disease spots usually appear after the flowering of the plant and may result in the withering of the plant in severe cases. This study examined the disease incidence in an herbaceous peony field in the Yangzhou region, Jiangsu Province. Based on morphological characteristics and molecular data, the disease in this area was identified as a gray mold caused by Botrytis cinerea. Based on previously obtained transcriptome data, eight libraries generated from two herbaceous peony cultivars ‘Zifengyu’ and ‘Dafugui’ with different susceptibilities to the disease were then analyzed using digital gene expression profiling (DGE). Thousands of differentially expressed genes (DEGs) were screened by comparing the eight samples, and these genes were annotated using the Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) database. The pathways related to plant-pathogen interaction, secondary metabolism synthesis and antioxidant system were concentrated, and 51, 76, and 13 disease resistance-relevant candidate genes were identified, respectively. The expression patterns of these candidate genes differed between the two cultivars: their expression of the disease-resistant cultivar ‘Zifengyu’ sharply increased during the early stages of infection, while it was relatively subdued in the disease-sensitive cultivar ‘Dafugui’. A selection of ten candidate genes was evaluated by quantitative real-time PCR (qRT-PCR) to validate the DGE data. These results revealed the transcriptional changes that took place during the interaction of herbaceous peony with B. cinerea, providing insight into the molecular mechanisms of host resistance to gray mold.

Transcriptome Analysis of Beta macrocarpa and Identification of Differentially Expressed Transcripts in Response to Beet Necrotic Yellow Vein Virus Infection

BACKGROUND

Rhizomania is one of the most devastating diseases of sugar beet. It is caused by Beet necrotic yellow vein virus (BNYVV) transmitted by the obligate root-infecting parasite Polymyxa betae. Beta macrocarpa, a wild beet species widely used as a systemic host in the laboratory, can be rub-inoculated with BNYVV to avoid variation associated with the presence of the vector P. betae. To better understand disease and resistance between beets and BNYVV, we characterized the transcriptome of B. macrocarpa and analyzed global gene expression of B. macrocarpa in response to BNYVV infection using the Illumina sequencing platform.

RESULTS

The overall de novo assembly of cDNA sequence data generated 75,917 unigenes, with an average length of 1054 bp. Based on a BLASTX search (E-value ≤ 10-5) against the non-redundant (NR, NCBI) protein, Swiss-Prot, the Gene Ontology (GO), Clusters of Orthologous Groups of proteins (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, there were 39,372 unigenes annotated. In addition, 4,834 simple sequence repeats (SSRs) were also predicted, which could serve as a foundation for various applications in beet breeding. Furthermore, comparative analysis of the two transcriptomes revealed that 261 genes were differentially expressed in infected compared to control plants, including 128 up- and 133 down-regulated genes. GO analysis showed that the changes in the differently expressed genes were mainly enrichment in response to biotic stimulus and primary metabolic process.

CONCLUSION

Our results not only provide a rich genomic resource for beets, but also benefit research into the molecular mechanisms of beet- BNYV Vinteraction.

Analysis of gene expression changes in peach leaves in response to Plum pox virus infection using RNA-Seq

Differences in gene expression were studied after Plum pox virus (PPV, sharka disease) infection in peach GF305 leaves with and without sharka symptoms using RNA-Seq. For each sample, more than 80% of 100-nucleotide paired-end (PE) Illumina reads were aligned on the peach reference genome. In the symptomatic sample, a significant proportion of reads were mapped to PPV reference genomes (1.04% compared with 0.00002% in non-symptomatic leaves), allowing for the ultra-deep assembly of the complete genome of the PPV isolate used (9775 nucleotides, missing only 11 nucleotides at the 5' genome end). In addition, significant alternative splicing events were detected in 359 genes and 12 990 single nucleotide polymorphisms (SNPs) were identified, 425 of which could be annotated. Gene ontology annotation revealed that the high-ranking mRNA target genes associated with the expression of sharka symptoms are mainly related to the response to biotic stimuli, to lipid and carbohydrate metabolism and to the negative regulation of catalytic activity. A greater number of differentially expressed genes were observed in the early asymptomatic phase of PPV infection in comparison with the symptomatic phase. These early infection events were associated with the induction of genes related to pathogen resistance, such as jasmonic acid, chitinases, cytokinin glucosyl transferases and Lys-M proteins. Once the virus had accumulated, the overexpression of Dicer protein 2a genes suggested a gene silencing plant response that was suppressed by the virus HCPro and P1 proteins. These results illustrate the dynamic nature of the peach-PPV interaction at the transcriptome level and confirm that sharka symptom expression is a complex process that can be understood on the basis of changes in plant gene expression.

Resistance to Plum pox virus strain C in Arabidopsis thaliana and Chenopodium foetidum involves genome-linked viral protein and other viral determinants and might depend on compatibility with host translation initiation factors

Research performed on model herbaceous hosts has been useful to unravel the molecular mechanisms that control viral infections. The most common Plum pox virus (PPV) strains are able to infect Nicotiana species as well as Chenopodium and Arabidopsis species. However, isolates belonging to strain C (PPV-C) that have been adapted to Nicotiana spp. are not infectious either in Chenopodium foetidum or in Arabidopsis thaliana. In order to determine the mechanism underlying this interesting host-specific behavior, we have constructed chimerical clones derived from Nicotiana-adapted PPV isolates from the D and C strains, which differ in their capacity to infect A. thaliana and C. foetidum. With this approach, we have identified the nuclear inclusion a protein (VPg+Pro) as the major pathogenicity determinant that conditions resistance in the presence of additional secondary determinants, different for each host. Genome-linked viral protein (VPg) mutations similar to those involved in the breakdown of eIF4E-mediated resistance to other potyviruses allow some PPV chimeras to infect A. thaliana. These results point to defective interactions between a translation initiation factor and the viral VPg as the most probable cause of host-specific incompatibility, in which other viral factors also participate, and suggest that complex interactions between multiple viral proteins and translation initiation factors not only define resistance to potyviruses in particular varieties of susceptible hosts but also contribute to establish nonhost resistance.

Tobacco rattle virus-based virus-induced gene silencing in Nicotiana benthamiana

Tobacco rattle virus (TRV)-based virus-induced gene silencing (VIGS) is widely used in various plant species to downregulate the expression of a target plant gene. TRV is a bipartite, positive-strand RNA virus with the TRV1 and TRV2 genomes. To induce post-transcriptional gene silencing (PTGS), the TRV2 genome is genetically modified to carry a fragment of the target gene and delivered into the plant (along with the TRV1 genome) by agroinoculation. TRV1- and TRV2-carrying Agrobacterium strains are then co-inoculated into 3-week-old plant leaves by one of three methods: a needleless syringe, the agrodrench method or by pricking with a toothpick. Target gene silencing occurs in the newly developed noninoculated leaves within 2-3 weeks of TRV inoculation. The TRV-VIGS protocol described here takes only 4 weeks to implement, and it is faster and easier to perform than other gene silencing techniques that are currently available. Although we use Nicotiana benthamiana as an example, the protocol is adaptable to other plant species.

Presence of poly(A) and poly(A)-rich tails in a positive-strand RNA virus known to lack 3׳ poly(A) tails

Here we show that Tobacco mosaic virus (TMV), a positive-strand RNA virus known to end with 3׳ tRNA-like structures, does possess a small fraction of gRNA bearing polyadenylate tails. Particularly, many tails are at sites corresponding to the 3׳ end of near full length gRNA, and are composed of poly(A)-rich sequences containing the other nucleotides in addition to adenosine, resembling the degradation-stimulating poly(A) tails observed in all biological kingdoms. Further investigations demonstrate that these polyadenylated RNA species are not enriched in chloroplasts. Silencing of cpPNPase, a chloroplast-localized polynucleotide polymerase known to not only polymerize the poly(A)-rich tails but act as a 3׳ to 5׳ exoribonuclease, does not change the profile of polyadenylate tails associated with TMV RNA. Nevertheless, because similar tails were also detected in other phylogenetically distinct positive-strand RNA viruses lacking poly(A) tails, such kind of polyadenylation may reflect a common but as-yet-unknown interface between hosts and viruses.

Comparative transcriptome analysis of two rice varieties in response to rice stripe virus and small brown planthoppers during early interaction

Rice stripe, a virus disease, transmitted by a small brown planthopper (SBPH), has greatly reduced production of japonica rice in East Asia, especially in China. Although we have made great progress in mapping resistance genes, little is known about the mechanism of resistance.

By de novo transcriptome assembling, we gained sufficient transcript data to analyze changes in gene expression of early interaction in response to SBPH and RSV infection in rice. Respectively 648 and 937 DEGs were detected from the disease-resistant (Liaonong 979) and the susceptible (Fengjin) varieties, most of which were up-regulated. We found 37 genes related to insect resistance, which mainly included genes for jasmonate-induced protein, TIFY protein, lipoxygenase, as well as trypsin inhibitor genes and transcription factor genes. In the interaction process between RSV and rice, 87 genes were thought to be related to RSV resistance; these primarily included 12 peroxidase biosynthesis genes, 12 LRR receptor-like protein kinase genes, 6 genes coding pathogenesis-related proteins, 4 glycine-rich cell wall structural protein genes, 2 xyloglucan hydrolase genes and a cellulose synthase. The results indicate that the rice-pathogen interaction happened both in disease-resistant and susceptible varieties, and some genes related to JA biosynthesis played key roles in the interaction between SBPHs and rice. When rice was infected by RSV a hypersensitive reaction (HR) in the disease-resistant variety was suppressed, which resulted from an increase in peroxidase expression and down-regulation of LRR receptor-like protein kinase and pathogenesis-related proteins, while, the changes of peroxidase biosynthesis, glycine-rich cell wall structural protein, cellulose synthase and xyloglucan endotransglucosylase/hydrolase could lead to the strengthening of physical barriers of rice, which may be an important resistance mechanism to RSV in rice.