Transcriptome Analyses Throughout Chili Pepper Fruit Development Reveal Novel Insights into the Domestication Process

Characterization and Expression Analysis of Sugar Transporters through Partial Least Square Structural Equation Model (PLS-SEM) Revealed Their Role in Pepper (Capsicum annuum L.)

Pepper (Capsicum annuum L.) is one of the most important economic crops in the world. By controlling the transport and distribution of photosynthetic products between cells and organs, sugar transporters are widely involved in growth and development, environmental adaptation, and microbial interactions. The present study was carried out at the genome-wide level to systematically characterize sugar transporters. As a result, 50 MST, 3 SUT, and 29 SWEET genes were identified and classified. The expression pattern of sugar transporters in pepper was analyzed by transcriptomic data. The expression properties of sugar transporters were further explored in pepper varieties with significant differences in weight, shape, and pungency. It was shown that the pepper sugar transporter genes had obvious spatiotemporal specific expression characteristics and exhibited variety-specific expression preferences. We focus on analyzing candidate genes that may be involved in fruit development and expansion. We further explore the response of pepper sugar transporters to adversity stress using a structural equation model. Finally, we found that the MST, SUT, and SWEET families are collectively involved in balancing pepper resistance to abiotic stress by coordinating the expression strengths of different family members. Our study may contribute to the functional study of pepper sugar transporter genes and create the prospect of utilizing sugar transporter gene resources to improve pepper variety.

Estimating Transcriptome Diversity and Specialization in Capsicum annuum L

Chili pepper fruits of the genus Capsicum represent excellent experimental models to study the growth, development, and ripening processes in a non-climacteric species at the physiological, biochemical, and molecular levels. Fruit growth, development, and ripening involve a complex, harmonious, and finely controlled regulation of gene expression. The purpose of this study was to estimate the changes in transcriptome diversity and specialization, as well as gene specificities during fruit development in this crop, and to illustrate the advantages of estimating these parameters. To achieve these aims, we programmed and made publicly available an R package. In this study, we applied these methods to a set of 179 RNA-Seq libraries from a factorial experiment that includes 12 different genotypes at various stages of fruit development. We found that the diversity of the transcriptome decreases linearly from the flower to the mature fruit, while its specialization follows a complex and non-linear behavior during this process. Additionally, by defining sets of genes with different degrees of specialization and applying Gene Ontology enrichment analysis, we identified processes, functions, and components that play a central role in particular fruit development stages. In conclusion, the estimation of diversity, specialization, and specificity summarizes the global properties of the transcriptomes, providing insights that are difficult to achieve by other means.

Dynamic shifts in primary metabolism across fruit development stages in Capsicum chinense (cv. Habanero)

The development of fleshy fruits involves changes in size and mass, followed by cell differentiation, which is associated with anatomical and histological changes. Parallel to these changes, metabolic alterations lead to the production of osmolytes and energy that modify cell turgor pressure, thereby promoting cell expansion and fruit growth. Detailed information is known about these processes in climacteric fruits (e.g. tomato); however, the regulation of metabolism and its association with anatomical changes in non-climacteric fruit development are poorly understood. In this study, we used detailed anatomical and histological analyses to define three developmental phases of chili pepper (Capsicum chinense cv. Habanero): cell division, cell expansion, and ripening. We showed that each was marked by distinct metabolic profiles, underpinning the switches in energy metabolism to support cellular processes. Interestingly, mitochondrial activity was high in the early stages of development and declined over time, with a modest increase in O2 consumption by pericarp tissues at the beginning of the ripening stage. This respiratory-like burst was associated with the degradation of starch and malate, which are the sources of energy and carbon required for other processes associated with fruit maturation.

Genetic Analyses of Seed Longevity in Capsicum annuum L. in Cold Storage Conditions

Seed longevity is the most important trait in the genebank management system. No seed can remain infinitely viable. There are 1241 accessions of Capsicum annuum L. available at the German Federal ex situ genebank at IPK Gatersleben. C. annuum (Capsicum) is the most economically important species of the genus Capsicum. So far, there is no report that has addressed the genetic basis of seed longevity in Capsicum. Here, we convened a total of 1152 Capsicum accessions that were deposited in Gatersleben over forty years (from 1976 to 2017) and assessed their longevity by analyzing the standard germination percentage after 5-40 years of storage at -15/-18 °C. These data were used to determine the genetic causes of seed longevity, along with 23,462 single nucleotide polymorphism (SNP) markers covering all of the 12 Capsicum chromosomes. Using the association-mapping approach, we identified a total of 224 marker trait associations (MTAs) (34, 25, 31, 35, 39, 7, 21 and 32 MTAs after 5-, 10-, 15-, 20-, 25-, 30-, 35- and 40-year storage intervals) on all the Capsicum chromosomes. Several candidate genes were identified using the blast analysis of SNPs, and these candidate genes are discussed.

Gene Functional Networks from Time Expression Profiles: A Constructive Approach Demonstrated in Chili Pepper (Capsicum annuum L.)

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).

Integrated analysis of DNA methylation, transcriptome, and global metabolites in interspecific heterotic Capsicum F1 hybrid

Hybrid breeding is one of the efficacious methods of crop improvement. Here, we report our work towards understanding the molecular basis of F1 hybrid heterosis from Capsicum chinense and C. frutescens cross. Bisulfite sequencing identified a total of 70597 CG, 108797 CHG, and 38418 CHH differentially methylated regions (DMRs) across F₁ hybrid and parents, and of these, 4891 DMRs showed higher methylation in F₁ compared to the mid-parental methylation values (MPMV). Transcriptome analysis showed higher expression of 46–55% differentially expressed genes (DE-Gs) in the F₁ hybrid. The qRT-PCR analysis of 24 DE-Gs with negative promoter methylation revealed 91.66% expression similarity with the transcriptome data. A few metabolites and 65–72% enriched genes in metabolite biosynthetic pathways showed overall increased expression in the F₁ hybrid compared to parents. These findings, taken together, provided insights into the integrated role of DNA methylation, and genes and metabolites expression in the manifestation of heterosis in Capsicum.

•

Global methylation identified significantly different proportions of mCs in hybrid

•

Of common DMRs, 33.08% showed different methylation in hybrid from the mid-parental value

•

Negatively correlated DEG pDMR-genes were enriched in metabolic pathways

•

Significant higher expression of metabolites and DE-Gs were identified in the F₁ hybrid

Putative Transcription Factor Genes Associated with Regulation of Carotenoid Biosynthesis in Chili Pepper Fruits Revealed by RNA-Seq Coexpression Analysis

During the ripening process, the pericarp of chili pepper (Capsicum spp.) fruits accumulates large amounts of carotenoids. Although the carotenoid biosynthesis pathway in the Capsicum genus has been widely studied from different perspectives, the transcriptional regulation of genes encoding carotenoid biosynthetic enzymes has not been elucidated in this fruit. We analyzed RNA-Seq transcriptomic data from the fruits of 12 accessions of Capsicum annuum during the growth, development, and ripening processes using the R package named Salsa. We performed coexpression analyses between the standardized expression of genes encoding carotenoid biosynthetic enzymes (target genes (TGs)) and the genes of all expressed transcription factors (TFs). Additionally, we analyzed the promoter region of each biosynthetic gene to identify putative binding sequences for each selected TF candidate. We selected 83 TFs as putative regulators of the carotenogenic structural genes. From them, putative binding sites in the promoters of the carotenoid-biosynthesis-related structural genes were found for only 54 TFs. These results could guide the search for transcription factors involved in the regulation of the carotenogenic pathway in chili pepper fruits and might facilitate the collection of corresponding experimental evidence to corroborate their participation in the regulation of this biosynthetic pathway in Capsicum spp.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

Chile Pepper (Capsicum) Breeding and Improvement in the "Multi-Omics" Era

Chile pepper (Capsicum spp.) is a major culinary, medicinal, and economic crop in most areas of the world. For more than hundreds of years, chile peppers have "defined" the state of New Mexico, USA. The official state question, "Red or Green?" refers to the preference for either red or the green stage of chile pepper, respectively, reflects the value of these important commodities. The presence of major diseases, low yields, decreased acreages, and costs associated with manual labor limit production in all growing regions of the world. The New Mexico State University (NMSU) Chile Pepper Breeding Program continues to serve as a key player in the development of improved chile pepper varieties for growers and in discoveries that assist plant breeders worldwide. Among the traits of interest for genetic improvement include yield, disease resistance, flavor, and mechanical harvestability. While progress has been made, the use of conventional breeding approaches has yet to fully address producer and consumer demand for these traits in available cultivars. Recent developments in "multi-omics," that is, the simultaneous application of multiple omics approaches to study biological systems, have allowed the genetic dissection of important phenotypes. Given the current needs and production constraints, and the availability of multi-omics tools, it would be relevant to examine the application of these approaches in chile pepper breeding and improvement. In this review, we summarize the major developments in chile pepper breeding and present novel tools that can be implemented to facilitate genetic improvement. In the future, chile pepper improvement is anticipated to be more data and multi-omics driven as more advanced genetics, breeding, and phenotyping tools are developed.

Inheritance of gene expression throughout fruit development in chili pepper

Gene expression is the primary molecular phenotype and can be estimated in specific organs or tissues at particular times. Here we analyzed genome-wide inheritance of gene expression in fruits of chili pepper (Capsicum annuum L.) in reciprocal crosses between a domesticated and a wild accession, estimating this parameter during fruit development. We defined a general hierarchical schema to classify gene expression inheritance which can be employed for any quantitative trait. We found that inheritance of gene expression is affected by both, the time of fruit development as well as the direction of the cross, and propose that such variations could be common in many developmental processes. We conclude that classification of inheritance patterns is important to have a better understanding of the mechanisms underlying gene expression regulation, and demonstrate that sets of genes with specific inheritance pattern at particular times of fruit development are enriched in different biological processes, molecular functions and cell components. All curated data and functions for analysis and visualization are publicly available as an R package.

Transcriptomic analysis of a wild and a cultivated varieties of Capsicum annuum over fruit development and ripening

Chili pepper (Capsicum annuum) is one of the most important crops worldwide. Its fruits contain metabolites produced over the maturation process like capsaicinoids and carotenoids. This metabolic process produces internal changes in flavor, color, texture, and aroma in fruits to make them more attractive for seed dispersal organisms. The chiltepin (C. annuum L. var. glabriusculum) is a wild variety of the C. annuum L. species that is considered a source of genetic resources that could be used to improve the current chili crops. In this study, we performed a transcriptomic analysis on two fruit maturation stages: immature stage (green fruit) and mature stage (red fruit) of a wild and a cultivated pepper variety. We found 19,811 genes expressed, and 1,008 genes differentially expressed (DEGs) in at least one of the five contrast used; 730 DEGs were found only in one contrast, and most DEGs in all contrasts were downregulated. GO enrichment analysis showed that the majority of DEGs are related to stress responses. KEGG enrichment analysis detected differences in expression patterns in metabolic pathways related to phenylpropanoid biosynthesis, secondary metabolites, plant hormone signal transduction, carotenoid biosynthesis and sesquiterpenoid and triterpenoid biosynthesis. We selected 105 tomato fruit ripening-related genes, and found 53 pepper homologs differentially expressed related to shape, size, and secondary metabolite biosynthesis. According to the transcriptome analysis, the two peppers showed very similar gene expression patterns; differences in expression patterns of genes related to shape, size, ethylene and secondary metabolites biosynthesis suggest that changes produced by domestication of chilli pepper could be very specific to the expression of genes related to traits desired in commercial fruits.

A method to analyze time expression profiles demonstrated in a database of chili pepper fruit development

RNA-Seq experiments allow genome-wide estimation of relative gene expression. Estimation of gene expression at different time points generates time expression profiles of phenomena of interest, as for example fruit development. However, such profiles can be complex to analyze and interpret. We developed a methodology that transforms original RNA-Seq data from time course experiments into standardized expression profiles, which can be easily interpreted and analyzed. To exemplify this methodology we used RNA-Seq data obtained from 12 accessions of chili pepper (Capsicum annuum L.) during fruit development. All relevant data, as well as functions to perform analyses and interpretations from this experiment, were gathered into a publicly available R package: “Salsa”. Here we explain the rational of the methodology and exemplify the use of the package to obtain valuable insights into the multidimensional time expression changes that occur during chili pepper fruit development. We hope that this tool will be of interest for researchers studying fruit development in chili pepper as well as in other angiosperms.