A Comprehensive Guide to Potato Transcriptome Assembly.

We have witnessed a rapid advancement in high-throughput genome sequencing and the maturation of long-read technologies. However, an accurate assembly of polyploid potato genomes still remains challenging. Sequencing the double-monoploid genome of Solanum tuberosum Group Phureja (Xu et al., Nature 475:189-195, 2011) has enabled functional studies of polyploid potato cultivars using RNA sequencing (RNA-Seq) technologies, although with the limitation of not covering cultivar-specific gene expression. The accumulated RNA-Seq datasets from these cultivars can be leveraged to assemble tetraploid potato transcriptomes that enable the analysis of genes that are not limited to reference genome annotations. To increase transcriptomes' quality, short-read assemblies are nowadays complemented with full-length transcriptome sequencing using Pacific Biosciences or Oxford Nanopore platforms. In this chapter we give a detailed guide on a pipeline for de novo transcriptome assembly of polyploid potato genotypes and their integration into a pan-transcriptome.

Methodologies for Discovery and Quantitative Profiling of sRNAs in Potato.

Small RNAs (sRNAs) are short noncoding RNAs involved in the regulation of a wide range of biological processes in plants. Advances in high-throughput sequencing and development of new computational tools had facilitated the discovery of different classes of sRNAs, their quantification, and elucidation of their functional role in gene expression regulation by target transcript predictions. The workflow presented here allows identification of different sRNA species: known and novel potato miRNAs, and their sequence variants (isomiRs), as well as identification of phased small interfering RNAs (phasiRNAs). Moreover, it includes steps for differential expression analysis to search for regulated sRNAs across different tested biological conditions. In addition, it describes two different methods for predicting sRNA targets, in silico prediction, and degradome sequencing data analysis. All steps of the workflow are written in a clear and user-friendly way; thus they can be followed also by the users with minimal bioinformatics knowledge. We also included several in-house scripts together with valuable notes to facilitate data (pre)processing steps and to reduce the analysis time.

Precision transcriptomics of viral foci reveals the spatial regulation of immune-signaling genes and identifies RBOHD as an important player in the incompatible interaction between potato virus Y and potato.

Whereas the activation of resistance (R) proteins has been intensively studied, the downstream signaling mechanisms leading to the restriction of the pathogen remain mostly unknown. We studied the immunity network response conditioned by the potato Ny-1 gene against potato virus Y. We analyzed the processes in the cell death zone and surrounding tissue on the biochemical and gene expression levels in order to reveal the spatiotemporal regulation of the immune response. We show that the transcriptional response in the cell death zone and surrounding tissue is dependent on salicylic acid (SA). For some genes the spatiotemporal regulation is completely lost in the SA-deficient line, whereas other genes show a different response, indicating multiple connections between hormonal signaling modules. The induction of NADPH oxidase RBOHD expression occurs specifically on the lesion border during the resistance response. In plants with silenced RBOHD, the functionality of the resistance response is perturbed and the spread of the virus is not arrested at the site of infection. RBOHD is required for the spatial accumulation of SA, and conversely RBOHD is under the transcriptional regulation of SA. Using spatially resolved RNA-seq, we also identified spatial regulation of an UDP-glucosyltransferase, another component in feedback activation of SA biosynthesis, thus deciphering a novel aspect of resistance signaling.

Cultivar-specific transcriptome and pan-transcriptome reconstruction of tetraploid potato.

Although the reference genome of Solanum tuberosum Group Phureja double-monoploid (DM) clone is available, knowledge on the genetic diversity of the highly heterozygous tetraploid Group Tuberosum, representing most cultivated varieties, remains largely unexplored. This lack of knowledge hinders further progress in potato research. In conducted investigation, we first merged and manually curated the two existing partially-overlapping DM genome-based gene models, creating a union of genes in Phureja scaffold. Next, we compiled available and newly generated RNA-Seq datasets (cca. 1.5 billion reads) for three tetraploid potato genotypes (cultivar Désirée, cultivar Rywal, and breeding clone PW363) with diverse breeding pedigrees. Short-read transcriptomes were assembled using several de novo assemblers under different settings to test for optimal outcome. For cultivar Rywal, PacBio Iso-Seq full-length transcriptome sequencing was also performed. EvidentialGene redundancy-reducing pipeline complemented with in-house developed scripts was employed to produce accurate and complete cultivar-specific transcriptomes, as well as to attain the pan-transcriptome. The generated transcriptomes and pan-transcriptome represent a valuable resource for potato gene variability exploration, high-throughput omics analyses, and breeding programmes.