Abstract
Habanero peppers constantly face biotic and abiotic stresses such as pathogen/pest infections, extreme temperature, drought and UV radiation. In addition, the fruit cutin lipid composition plays an important role in post-harvest water loss rates, which in turn causes shriveling and reduced fruit quality and storage. In this study, we integrated metabolome and transcriptome profiling pertaining to cutin in two habanero genotypes: PI 224448 and PI 257145. The fruits were selected by the waxy or glossy phenotype on their surfaces. Metabolomics analysis showed a significant variation in cutin composition, with about 6-fold higher cutin in PI 257145 than PI 224448. It also revealed that 10,16-dihydroxy hexadecanoic acid is the most abundant monomer in PI 257145. Transcriptomic analysis of high-cutin PI 257145 and low-cutin PI 224448 resulted in the identification of 2703 statistically significant differentially expressed genes, including 1693 genes upregulated and 1010 downregulated in high-cutin PI 257145. Genes and transcription factors such as GDSL lipase, glycerol-3 phosphate acyltransferase 6, long-chain acyltransferase 2, cytochrome P450 86A/77A, SHN1, ANL2 and HDG1 highly contributed to the high cutin content in PI 257145. We predicted a putative cutin biosynthetic pathway for habanero peppers based on deep transcriptome analysis. This is the first study of the transcriptome and metabolome pertaining to cutin in habanero peppers. These analyses improve our knowledge of the molecular mechanisms regulating the accumulation of cutin in habanero pepper fruits. These resources can be built on for developing cultivars with high cutin content that show resistance to biotic and abiotic stresses with superior postharvest appearance.