PepperHub, an Informatics Hub for the Chili Pepper Research Community
Feng Liu,Huiyang Yu,Yingtian Deng,Jingyuan Zheng,Minglei Liu,Lijun Ou,Bozhi Yang,Xiongze Dai,Yanqing Ma,Shengyu Feng,Shuang He,Xuefeng Li,Zhuqing Zhang,Wenchao Chen,Shudong Zhou,Rong Chen,Minmin Liu,Sha Yang,Ruimin Wei,Huadong Li
DOI: https://doi.org/10.1016/j.molp.2017.03.005
IF: 27.5
2017-01-01
Molecular Plant
Abstract:Pepper belongs to the Solanaceae family, which includes many important vegetable crops such as tomato, potato, and eggplant. Not only widely used as vegetables and spicy ingredients, pepper also has diverse applications in pharmaceutics, natural coloring agents, cosmetics, defense repellents, and as ornamental plants (Kim et al., 2014Kim S. Park M. Yom S.I. Kim Y.M. Lee J.M. Lee H.A. Seo E. Choi J. Cheong K. Kim K.T. et al.Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species.Nat. Genet. 2014; 46: 270-278Crossref PubMed Scopus (646) Google Scholar, Qin et al., 2014Qin C. Yu C. Shen Y. Fang X. Chen L. Min J. Cheng J. Zhao S. Xu M. Luo Y. et al.Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization.Proc. Natl. Acad. Sci. USA. 2014; 111: 5135-5140Crossref PubMed Scopus (506) Google Scholar). Pepper is among the most widely cultivated and consumed vegetables in the world, with annual production reaching to 38 million tons in 2011 (www.fao.org). Pepper fruits have significant diversity in morphology and color, and they provide good models for fruit developmental biology (Paran and van der Knaap, 2007Paran I. van der Knaap E. Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper.J. Exp. Bot. 2007; 58: 3841-3852Crossref PubMed Scopus (248) Google Scholar, Rivera et al., 2016Rivera A. Monteagudo A.B. Igartua E. Taboada A. Garcia-Ulloa A. Pomar F. Riveiro-Leira M. Silvar C. Assessing genetic and phenotypic diversity in pepper (Capsicum annuum L.) landraces from North-West Spain.Sci. Hortic. (Amsterdam). 2016; 203: 1-11Crossref Scopus (32) Google Scholar). Like all other crops, pepper plants are often confronted with different pathogens and pests (Pernezny et al., 2003Pernezny K. Roberts P.D. Murphy J.F. Goldberg N.P. Compendium of Pepper Diseases. APS Press, St. Paul, MN2003Google Scholar), and diverse abiotic stress conditions, which necessitate basic studies on the mechanisms of pepper plants responding to various stimuli to facilitate breeding efforts for tolerant cultivars. Centralized and specialized informatic webservers, such as TAIR (www.arabidopsis.org) (Lamesch et al., 2012Lamesch P. Berardini T.Z. Li D. Swarbreck D. Wilks C. Sasidharan R. Muller R. Dreher K. Alexander D.L. Garcia-Hernandez M. et al.The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools.Nucleic Acids Res. 2012; 40: D1202-D1210Crossref PubMed Scopus (1366) Google Scholar) and Solgenomics (www.solgenomics.net) (Fernandez-Pozo et al., 2015Fernandez-Pozo N. Menda N. Edwards J.D. Saha S. Tecle I.Y. Strickler S.R. Bombarely A. Fisher-York T. Pujar A. Forester H. et al.The Sol Genomics Network (SGN)–from genotype to phenotype to breeding.Nucleic Acids Res. 2015; 43: D1036-D1041Crossref PubMed Scopus (343) Google Scholar), have integrated the genomic, transcriptomic, and proteomic data for their targeted plant species, provided convenient public platforms for the research communities, and greatly promoted related research. Despite its significant agricultural importance and interesting fruit biology, development of genetic resources and informatics platforms for pepper research has lagged far behind that of other crop plants such as rice and tomato. The release of three draft pepper genomes constituted a milestone in the development of genetic resources for the pepper research community (Kim et al., 2014Kim S. Park M. Yom S.I. Kim Y.M. Lee J.M. Lee H.A. Seo E. Choi J. Cheong K. Kim K.T. et al.Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species.Nat. Genet. 2014; 46: 270-278Crossref PubMed Scopus (646) Google Scholar, Qin et al., 2014Qin C. Yu C. Shen Y. Fang X. Chen L. Min J. Cheng J. Zhao S. Xu M. Luo Y. et al.Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization.Proc. Natl. Acad. Sci. USA. 2014; 111: 5135-5140Crossref PubMed Scopus (506) Google Scholar). Although some of the data are accessible from public databases, a comprehensive informatics platform remains to be developed to cover different aspects of omics data, and to be friendly to biologists without bioinformatics expertise. Here we present the Pepper Informatics Hub (PepperHub), which consists of five main modules, including Genome (pepper genome database), Transcriptome (pepper transcriptome database), sRNome (pepper small RNA database), Variome (pepper genetic variation), and Proteome (pepper proteome database) (Figure 1A). The Genome module hosts the reference genomes of zunla and CM334 and provides Gbrowse and BLAST functions. Gbrowse allows users to browse the gene structure and sequence information and BLAST allows users to query and retrieve the genomic DNA, cDNA, and protein sequences of both Zunla and CM334 genomes (Kim et al., 2014Kim S. Park M. Yom S.I. Kim Y.M. Lee J.M. Lee H.A. Seo E. Choi J. Cheong K. Kim K.T. et al.Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species.Nat. Genet. 2014; 46: 270-278Crossref PubMed Scopus (646) Google Scholar, Qin et al., 2014Qin C. Yu C. Shen Y. Fang X. Chen L. Min J. Cheng J. Zhao S. Xu M. Luo Y. et al.Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization.Proc. Natl. Acad. Sci. USA. 2014; 111: 5135-5140Crossref PubMed Scopus (506) Google Scholar). The Variome module allows users to browse or search for genetic variation data among the pepper re-sequencing population (Figure 1B). The SNP and INDEL search functions allow users to retrieve SNP or INDEL in a table format for a given chromosome region or a region encompassing a given gene; the Accession section allows users to search the allele information of SNP/INDEL in each re-sequenced accession; the Gbrowse function allows users to browse the SNP and INDEL in the genome browser. The Proteome module is designed to access pepper protein-protein interaction data (see Supplemental Materials and Methods) and allows users to search and visualize interacting partners for a given gene product with the option of showing a one- or two-layer interaction network (Figure 1C). The sRNome module hosts the annotated pepper miRNA and published small RNA databases (Kim et al., 2014Kim S. Park M. Yom S.I. Kim Y.M. Lee J.M. Lee H.A. Seo E. Choi J. Cheong K. Kim K.T. et al.Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species.Nat. Genet. 2014; 46: 270-278Crossref PubMed Scopus (646) Google Scholar, Qin et al., 2014Qin C. Yu C. Shen Y. Fang X. Chen L. Min J. Cheng J. Zhao S. Xu M. Luo Y. et al.Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization.Proc. Natl. Acad. Sci. USA. 2014; 111: 5135-5140Crossref PubMed Scopus (506) Google Scholar, Liu et al., 2017Liu Z. Zhang Y. Ou L. Kang L. Liu Y. Lv J. Wei G. Yang B. Yang S. Chen W. et al.Identification and characterization of novel microRNAs for fruit development and quality in hot pepper (Capsicum annuum L.).Gene. 2017; 608: 66-72Crossref PubMed Scopus (34) Google Scholar) and allows users to analyze small RNA accumulation from miRNA and other sRNA producing loci. A target search function was added to the sRNome module for users to search targets of annotated miRNAs (Liu et al., 2017Liu Z. Zhang Y. Ou L. Kang L. Liu Y. Lv J. Wei G. Yang B. Yang S. Chen W. et al.Identification and characterization of novel microRNAs for fruit development and quality in hot pepper (Capsicum annuum L.).Gene. 2017; 608: 66-72Crossref PubMed Scopus (34) Google Scholar). The Transcriptome module includes a large volume of new transcriptome data resulting from high-throughput mRNA sequencing using an elite pepper breeding line 6421 (see Supplemental Materials and Methods and Supplemental Figure 1). Messenger RNAs from 188 samples of different organs/tissues during successive developmental stages or different stressed plants at consecutive time points were sequenced in triplicate and over 3 Tb of transcriptome data were produced (Supplemental Data 1; Supplemental Tables 1 and 2). The Transcriptome web module was constructed for users to retrieve gene expression data, and to visualize the co-expression network for all pepper genes (Figure 1D; Supplemental Figure 2 and Supplemental Data 1). A default user name “123” and password “123” are used for security purposes. Here, we use pepper MYB transcription factors’ gene expression as a case study to demonstrate how to obtain useful information using the Transcriptome module of PepperHub. To this end, a comprehensive list of MYB transcription factor genes containing the characteristic DNA-binding domain was identified (see the Supplemental Materials and Methods and Supplemental Table 3). First, a heatmap was generated for the expression levels of all pepper MYB genes in various experiments (Supplemental Figure 3) with the ProfileHeatmap function on the webserver (Supplemental Figure 2C). From these results, a few MYB genes were identified, which showed interesting patterns during pericarp development, ABA treatment, and other biological processes (Figure 1D and Supplemental Figure 4). To identify potential MYB transcription factors involved in the regulation of capsanthin biosynthesis, co-expressed genes were retrieved for MYB transcription factors expressed in pericarp (Figure 1D) with the CoExpNetwork function on the webserver (Supplemental Figure 2E and Supplemental Table 4). A comprehensive list of pepper genes involved in catalysis of capsanthin was identified by querying the annotated pepper genes with known enzymes in the pathway using BLAST and KEGG (Goldstein and Brown, 1990Goldstein J.L. Brown M.S. Regulation of the mevalonate pathway.Nature. 1990; 343: 425-430Crossref PubMed Scopus (4544) Google Scholar, Guzman et al., 2010Guzman I. Hamby S. Romero J. Bosland P.W. O'Connell M.A. Variability of carotenoid biosynthesis in orange colored Capsicum spp.Plant Sci. 2010; 179: 49-59Crossref PubMed Scopus (108) Google Scholar) (Supplemental Table 5 and Supplemental Figure 5A). By comparing the MYB co-expressed gene list and the capsanthin biosynthesis gene list, several MYB transcription factors were identified. Five of them co-expressed with genes encoding enzymes catalyzing mevalonate synthesis, while the sixth one co-expressed with several genes encoding enzymes catalyzing downstream reactions from geranylgeranyl diphosphate to lutein, capsanthin, and capsorubin (Supplemental Figure 5A). Further analyses of the expression pattern of these MYB genes using the ProfileCartoon and ProfileChart functions (Supplemental Figure 2A and 2D) revealed that two of them, Capana12g002172 and Capana03g000766, showed preferential expression in pericarp (Figure 1E and Supplemental Figure 5B) and good co-expression pattern with capsanthin biosynthesis genes (Figure 1F and 1G). These results indicated that pepper MYB transcription factors Capana12g002172 and Capana03g000766 may play a potential role in regulating capsanthin biosynthesis during pepper fruit development through activation of different structural genes in the capsanthin biosynthesis pathway. Toward building a central public data platform for the pepper research community, we constructed PepperHub, which not only provides access to published pepper genome and small RNA data but also offers a large volume of newly generated high-quality transcriptome datasets (Supplemental Table 1) with in-depth coverage of the annotated pepper genes (Supplemental Table 2). Besides providing access to public data for pepper research, PepperHub also aims to be user friendly, particularly to molecular biologists without bioinformatics expertise. Various web-based tools are provided for users to retrieve and visualize different types of data. For vegetable plants, transcriptome datasets can be found and downloaded from public databases such as TED (Fei et al., 2006Fei Z. Tang X. Alba R. Giovannoni J. Tomato Expression Database (TED): a suite of data presentation and analysis tools.Nucleic Acids Res. 2006; 34: D766-D770Crossref PubMed Scopus (64) Google Scholar), and the eFP Browser provides analysis tools for one or two genes (Winter et al., 2007Winter D. Vinegar B. Nahal H. Ammar R. Wilson G.V. Provart N.J. An “Electronic Fluorescent Pictograph” browser for exploring and analyzing large-scale biological data sets.PLoS One. 2007; 2: e718Crossref PubMed Scopus (1856) Google Scholar). However, there is no data server providing tools for high-throughput analyses of vegetable transcriptome data. The PepperHub transcriptome module has integrated these functionalities and provides more option for users to visualize expression of median to large sets of genes using line charts and heatmaps, which allow users to screen for genes of special interest from a long list of genes. This function was demonstrated by our case study of the gene expression of the pepper MYB family. Our analyses identified two pericarp-specific MYB transcription factors, Capana12g002172 and Capana03g000766, for which expression increased right at the onset (G5 in Figure 1F and 1G) of the fruit color break stage (Supplemental Figure 1A and 2B), suggesting that they may play a potential role in the regulation of fruit color break. In summary, PepperHub provides an integrative public platform for sharing and analyzing research data. Undoubtedly, PepperHub will accelerate research in pepper functional genomics and would serve as a valuable resource for studying fruit developmental biology and stress responses in general. PepperHub is available at http://www.hnivr.org/pepperhub. This work was supported by the National Key Research and Development Program of China (2016YFD0101704), National Science Foundation of China (31470105), and Huazhong Agricultural University startup fund (2013RC001).