RNA Sequencing Analyses for Deciphering Potato Molecular Responses

Investigating the genetic components of tuber bruising in a breeding population of tetraploid potatoes

BACKGROUND

Tuber bruising in tetraploid potatoes (Solanum tuberosum) is a trait of economic importance, as it affects tubers' fitness for sale. Understanding the genetic components affecting tuber bruising is a key step in developing potato lines with increased resistance to bruising. As the tetraploid setting renders genetic analyses more complex, there is still much to learn about this complex phenotype. Here, we used capture sequencing data on a panel of half-sibling populations from a breeding programme to perform a genome-wide association analysis (GWAS) for tuber bruising. In addition, we collected transcriptomic data to enrich the GWAS results. However, there is currently no satisfactory method to represent both GWAS and transcriptomics analysis results in a single visualisation and to compare them with existing knowledge about the biological system under study.

RESULTS

When investigating population structure, we found that the STRUCTURE algorithm yielded greater insights than discriminant analysis of principal components (DAPC). Importantly, we found that markers with the highest (though non-significant) association scores were consistent with previous findings on tuber bruising. In addition, new genomic regions were found to be associated with tuber bruising. The GWAS results were backed by the transcriptomics differential expression analysis. The differential expression notably highlighted for the first time the role of two genes involved in cellular strength and mechanical force sensing in tuber resistance to bruising. We proposed a new visualisation, the HIDECAN plot, to integrate the results from the genomics and transcriptomics analyses, along with previous knowledge about genomic regions and candidate genes associated with the trait.

CONCLUSION

This study offers a unique genome-wide exploration of the genetic components of tuber bruising. The role of genetic components affecting cellular strength and resistance to physical force, as well as mechanosensing mechanisms, was highlighted for the first time in the context of tuber bruising. We showcase the usefulness of genomic data from breeding programmes in identifying genomic regions whose association with the trait of interest merit further investigation. We demonstrate how confidence in these discoveries and their biological relevance can be increased by integrating results from transcriptomics analyses. The newly proposed visualisation provides a clear framework to summarise of both genomics and transcriptomics analyses, and places them in the context of previous knowledge on the trait of interest.

Transcriptional deregulation of stress-growth balance in Nicotiana benthamiana biofactories producing insect sex pheromones

Plant biofactories are a promising platform for sustainable production of high-value compounds, among which are insect sex pheromones, a green alternative to conventional insecticides in agriculture. Recently, we have constructed transgenic Nicotiana benthamiana plants ("Sexy Plants", SxP) that successfully produce a blend of moth (Lepidoptera) sex pheromone compounds (Z)-11-hexadecen-1-ol and (Z)-11-hexadecenyl acetate. However, efficient biosynthesis of sex pheromones resulted in growth and developmental penalty, diminishing the potential for commercial use of SxP in biomanufacturing. To gain insight into the underlying molecular responses, we analysed the whole-genome transcriptome and evaluated it in relation to growth and pheromone production in low- and high-producing transgenic plants of v1.0 and v1.2 SxP lines. In our study, high-producing SxPv1.2 plants accumulated the highest amounts of pheromones but still maintained better growth compared to v1.0 high producers. For an in-depth biological interpretation of the transcriptomic data, we have prepared a comprehensive functional N. benthamiana genome annotation as well as gene translations to Arabidopsis thaliana, enabling functional information transfer by using Arabidopsis knowledge networks. Differential gene expression analysis, contrasting pheromone producers to wild-type plants, revealed that while only a few genes were differentially regulated in low-producing plants, high-producing plants exhibited vast transcriptional reprogramming. They showed signs of stress-like response, manifested as downregulation of photosynthesis-related genes and significant differences in expression of hormonal signalling and secondary metabolism-related genes, the latter presumably leading to previously reported volatilome changes. Further network analyses confirmed stress-like response with activation of jasmonic acid and downregulation of gibberellic acid signalling, illuminating the possibility that the observed growth penalty was not solely a consequence of a higher metabolic burden imposed upon constitutive expression of a heterologous biosynthetic pathway, but rather the result of signalling pathway perturbation. Our work presents an example of comprehensive transcriptomic analyses of disadvantageous stress signalling in N. benthamiana biofactory that could be applied to other bioproduction systems.