Chili pepper (Capsicum spp.) is one of the most important horticultural crops worldwide, and its unique organoleptic properties and health benefits have been established for centuries. However, there is little knowledge about how metabolites are distributed throughout fruit parts. This work focuses on the use of liquid chromatography coupled with high resolution mass spectrometry (UHPLC-ESI-HRMS) to estimate the global metabolite profiles of the pericarp, placenta, and seeds of Tabasco pepper fruits (Capsicum frutescens L.) at the red mature stage of ripening. Our main results putatively identified 60 differential compounds between these tissues and seeds. Firstly, we found that pericarp has a higher content of glycosides, showing on average a fold change of 5 and a fold change of 14 for terpenoids when compared with other parts of the fruit. While placenta was the richest tissue in capsaicinoid-related compounds, alkaloids, and tocopherols, with a 35, 3, and 7 fold change, respectively. However, the seeds were richer in fatty acids and saponins with fold changes of 86 and 224, respectively. Therefore, our study demonstrates that a non-targeted metabolomic approach may help to improve our understanding of unexplored areas of plant metabolism and also may be the starting point for a detailed analysis in complex plant parts, such as fruits.

Chili pepper (Capsicum spp.) is an important crop, as well as a model for fruit development studies and domestication. Here, we performed a time-course experiment to estimate standardized gene expression profiles with respect to fruit development for six domesticated and four wild chili pepper ancestors. We sampled the transcriptomes every 10 days from flowering to fruit maturity, and found that the mean standardized expression profiles for domesticated and wild accessions significantly differed. The mean standardized expression was higher and peaked earlier for domesticated vs. wild genotypes, particularly for genes involved in the cell cycle that ultimately control fruit size. We postulate that these gene expression changes are driven by selection pressures during domestication and show a robust network of cell cycle genes with a time shift in expression, which explains some of the differences between domesticated and wild phenotypes.

Gene co-expression networks are powerful tools to understand functional interactions between genes. However, large co-expression networks are difficult to interpret and do not guarantee that the relations found will be true for different genotypes. Statistically verified time expression profiles give information about significant changes in expressions through time, and genes with highly correlated time expression profiles, which are annotated in the same biological process, are likely to be functionally connected. A method to obtain robust networks of functionally related genes will be useful to understand the complexity of the transcriptome, leading to biologically relevant insights. We present an algorithm to construct gene functional networks for genes annotated in a given biological process or other aspects of interest. We assume that there are genome-wide time expression profiles for a set of representative genotypes of the species of interest. The method is based on the correlation of time expression profiles, bound by a set of thresholds that assure both, a given false discovery rate, and the discard of correlation outliers. The novelty of the method consists in that a gene expression relation must be repeatedly found in a given set of independent genotypes to be considered valid. This automatically discards relations particular to specific genotypes, assuring a network robustness, which can be set a priori. Additionally, we present an algorithm to find transcription factors candidates for regulating hub genes within a network. The algorithms are demonstrated with data from a large experiment studying gene expression during the development of the fruit in a diverse set of chili pepper genotypes. The algorithm is implemented and demonstrated in a new version of the publicly available R package “Salsa” (version 1.0).

Chili pepper ( Capsicum spp.) is both an important crop and a model for domestication studies. Here we performed a time course experiment to estimate standardized gene expression profiles across fruit development for six domesticated and four wild chili pepper ancestors. We sampled the transcriptome every 10 days, from flower to fruit at 60 Days After Anthesis (DAA), and found that the mean standardized expression profile for domesticated and wild accessions significantly differed. The mean standardized expression was higher and peaked earlier for domesticated vs. wild genotypes, particularly for genes involved in the cell cycle that ultimately control fruit size. We postulate that these gene expression changes are driven by selection pressures during domestication and show a robust network of cell cycle genes with a time-shift in expression which explains some of the differences between domesticated and wild phenotypes.