Pepper EST database: comprehensive in silico tool for analyzing the chili pepper (Capsicum annuum) transcriptome

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

An Ontology-Driven Personalized Faceted Search for Exploring Knowledge Bases of Capsicum

Capsicum is a genus of flowering plants in the Solanaceae family in which the members are well known to have a high economic value. The Capsicum fruits, which are popularly known as peppers or chili, have been widely used by people worldwide. It serves as a spice and raw material for many products such as sauce, food coloring, and medicine. For many years, scientists have studied this plant to optimize its production. A tremendous amount of knowledge has been obtained and shared, as reflected in multiple knowledge-based systems, databases, or information systems. An approach to knowledge-sharing is through the adoption of a common ontology to eliminate knowledge understanding discrepancy. Unfortunately, most of the knowledge-sharing solutions are intended for scientists who are familiar with the subject. On the other hand, there are groups of potential users that could benefit from such systems but have minimal knowledge of the subject. For these non-expert users, finding relevant information from a less familiar knowledge base would be daunting. More than that, users have various degrees of understanding of the available content in the knowledge base. This understanding discrepancy raises a personalization problem. In this paper, we introduce a solution to overcome this challenge. First, we developed an ontology to facilitate knowledge-sharing about Capsicum to non-expert users. Second, we developed a personalized faceted search algorithm that provides multiple structured ways to explore the knowledge base. The algorithm addresses the personalization problem by identifying the degree of understanding about the subject from each user. In this way, non-expert users could explore a knowledge base of Capsicum efficiently. Our solution characterized users into four groups. As a result, our faceted search algorithm defines four types of matching mechanisms, including three ranking mechanisms as the core of our solution. In order to evaluate the proposed method, we measured the predictability degree of produced list of facets. Our findings indicated that the proposed matching mechanisms could tolerate various query types, and a high degree of predictability can be achieved by combining multiple ranking mechanisms. Furthermore, it demonstrates that our approach has a high potential contribution to biodiversity science in general, where many knowledge-based systems have been developed with limited access to users outside of the domain.

Global gene expression profiling for fruit organs and pathogen infections in the pepper, Capsicum annuum L

Hot pepper (Capsicum annuum) is one of the most consumed vegetable crops in the world and useful to human as it has many nutritional and medicinal values. Genomic resources of pepper are publically available since the pepper genomes have been completed and massive data such as transcriptomes have been deposited. Nevertheless, global transcriptome profiling is needed to identify molecular mechanisms related to agronomic traits in pepper, but limited analyses are published. Here, we report the comprehensive analysis of pepper transcriptomes during fruit ripening and pathogen infection. For the ripening, transcriptome data were obtained from placenta and pericarp at seven developmental stages. To reveal global transcriptomic landscapes during infection, leaves at six time points post-infection by one of three pathogens (Phytophthora infestans, Pepper mottle virus, and Tobacco mosaic virus P0 strain) were profiled. The massive parallel transcriptome profiling in this study will serve as a valuable resource for detection of molecular networks of fruit development and disease resistance in Capsicum annuum.

Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae)

The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.

Whole Genome Resequencing of Capsicum baccatum and Capsicum annuum to Discover Single Nucleotide Polymorphism Related to Powdery Mildew Resistance

The present study deals with genome wide identification of single-nucleotide polymorphism (SNP) markers related to powdery mildew (PM) resistance in two pepper varieties. Capsicum baccatum (PRH1- a PM resistant line) and Capsicum annuum (Saengryeg- a PM susceptible line), were resequenced to develop SNP markers. A total of 6,213,009 and 6,840,889 SNPs for PRH1 and Saengryeg respectively have been discovered. Among the SNPs, majority were classified as homozygous type SNPs, particularly in the resistant line. Moreover, the SNPs were differentially distributed among the chromosomes in both the resistant and susceptible lines. In total, 4,887,031 polymorphic SNP loci were identified between the two lines and 306,871 high-resolution melting (HRM) marker primer sets were designed. In order to understand the SNPs associated with the vital genes involved in diseases resistance and stress associated processes, chromosome-wise gene ontology analysis was performed. The results revealed the occurrence that SNPs related to diseases resistance genes were predominantly distributed in chromosome 4. In addition, 6281 SNPs associated with 46 resistance genes were identified. Among the lines, PRH1 consisted of maximum number of polymorphic SNPs related to NBS-LRR genes. The SNP markers were validated using HRM assay in 45 F4 populations and correlated with the phenotypic disease index.

Next-Generation Sequencing Approaches in Genome-Wide Discovery of Single Nucleotide Polymorphism Markers Associated with Pungency and Disease Resistance in Pepper

Pepper is an economically important horticultural plant that has been widely used for its pungency and spicy taste in worldwide cuisines. Therefore, the domestication of pepper has been carried out since antiquity. Owing to meet the growing demand for pepper with high quality, organoleptic property, nutraceutical contents, and disease tolerance, genomics assisted breeding techniques can be incorporated to develop novel pepper varieties with desired traits. The application of next-generation sequencing (NGS) approaches has reformed the plant breeding technology especially in the area of molecular marker assisted breeding. The availability of genomic information aids in the deeper understanding of several molecular mechanisms behind the vital physiological processes. In addition, the NGS methods facilitate the genome-wide discovery of DNA based markers linked to key genes involved in important biological phenomenon. Among the molecular markers, single nucleotide polymorphism (SNP) indulges various benefits in comparison with other existing DNA based markers. The present review concentrates on the impact of NGS approaches in the discovery of useful SNP markers associated with pungency and disease resistance in pepper. The information provided in the current endeavor can be utilized for the betterment of pepper breeding in future.

Cloning and expression analysis of a new anther-specific gene CaMF4 in Capsicum annuum

Our previous study on the genic male sterile-fertile line 114AB of Capsicum annuum indicated a diversity of differentially expressed cDNA fragments in fertile and sterile lines. In this study, a transcript-derived fragment (TDF), male fertile 4 (CaMF4) was chosen for further investigation to observe that this specific fragment accumulates in the flower buds of the fertile line. The full genomic DNA sequence of CaMF4 was 894 bp in length, containing two exons and one intron, and the complete coding sequence encoded a putative 11.53 kDa protein of 109 amino acids. The derived protein of CaMF4 shared similarity with the members of PGPS/D3 protein family. The expression of CaMF4 was detected in both the flower buds at stage 8 and open flowers of the male fertile line. In contrast to this observation, expression of CaMF4 was not detected in any organs of the male sterile line. Further analysis revealed that CaMF4 was expressed particularly in anthers of the fertile line. Our results suggest that CaMF4 is an anther-specific gene and might be indispensable for anther or pollen development in C. annuum.

Fine mapping of Restorer-of-fertility in pepper (Capsicum annuum L.) identified a candidate gene encoding a pentatricopeptide repeat (PPR)-containing protein

Using fine mapping techniques, the genomic region co-segregating with Restorer - of - fertility ( Rf ) in pepper was delimited to a region of 821 kb in length. A PPR gene in this region, CaPPR6 , was identified as a strong candidate for Rf based on expression pattern and characteristics of encoding sequence. Cytoplasmic-genic male sterility (CGMS) has been used for the efficient production of hybrid seeds in peppers (Capsicum annuum L.). Although the mitochondrial candidate genes that might be responsible for cytoplasmic male sterility (CMS) have been identified, the nuclear Restorer-of-fertility (Rf) gene has not been isolated. To identify the genomic region co-segregating with Rf in pepper, we performed fine mapping using an Rf-segregating population consisting of 1068 F2 individuals, based on BSA-AFLP and a comparative mapping approach. Through six cycles of chromosome walking, the co-segregating region harboring the Rf locus was delimited to be within 821 kb of sequence. Prediction of expressed genes in this region based on transcription analysis revealed four candidate genes. Among these, CaPPR6 encodes a pentatricopeptide repeat (PPR) protein with PPR motifs that are repeated 14 times. Characterization of the CaPPR6 protein sequence, based on alignment with other homologs, showed that CaPPR6 is a typical Rf-like (RFL) gene reported to have undergone diversifying selection during evolution. A marker developed from a sequence near CaPPR6 showed a higher prediction rate of the Rf phenotype than those of previously developed markers when applied to a panel of breeding lines of diverse origin. These results suggest that CaPPR6 is a strong candidate for the Rf gene in pepper.

Genome-wide analysis of Dof transcription factors reveals functional characteristics during development and response to biotic stresses in pepper

The DNA-binding with one zinc finger proteins (Dofs) are a plant-specific family of transcription factors. The Dofs are involved in a variety of biological processes such as phytohormone production, seed development, and environmental adaptation. Dofs have been previously identified in several plants, but not in pepper. We identified 33 putative Dof genes in pepper (CaDofs). To gain an overview of the CaDofs, we analyzed phylogenetic relationships, protein motifs, and evolutionary history. We divided the 33 CaDofs, containing 25 motifs, into four major groups distributed on eight chromosomes. We discovered an expansion of the CaDofs dated to a recent duplication event. Segmental duplication that occurred before the speciation of the Solanaceae lineages was predominant among the CaDofs. The global gene-expression profiling of the CaDofs by RNA-seq analysis showed distinct temporal and pathogen-specific variation during development and response to biotic stresses (two TMV strains, PepMoV, and Phytophthora capsici), suggesting functional diversity among the CaDofs. These results will provide the useful clues into the responses of Dofs in biotic stresses and promote a better understanding of their multiple function in pepper and other species.

Bioinformatics resources for pollen

Bioinformatics for Pollen. Pollen plays a key role in crop production, and its development is the most delicate phase in reproduction. Different metabolic pathways are involved in pollen development, and changes in the level of some metabolites, as well as responses to stress, are correlated with the reduction in pollen viability, leading consequently to a decrease in the fruit production. However, studies on pollen may be hard because gamete development and fertilization are complex processes that occur during a short window of time. The rise of the so-called -omics sciences provided key strategies to promote molecular research in pollen tissues, starting from model organisms and moving to increasing number of species. An integrated multi-level approach based on investigations from genomics, transcriptomics, proteomics and metabolomics appears now feasible to clarify key molecular processes in pollen development and viability. To this aim, bioinformatics has a fundamental role for data production and analysis, contributing varied and ad hoc methodologies, endowed with different sensitivity and specificity, necessary for extracting added-value information from the large amount of molecular data achievable. Bioinformatics is also essential for data management, organization, distribution and integration in suitable resources. This is necessary to catch the biological features of the pollen tissues and to design effective approaches to identifying structural or functional properties, enabling the modeling of the major involved processes in normal or in stress conditions. In this review, we provide an overview of the available bioinformatics resources for pollen, ranging from raw data collections to complete databases or platforms, when available, which include data and/or results from -omics efforts on the male gametophyte. Perspectives in the fields will also be described.

Expressed Sequence Tags Analysis and Design of Simple Sequence Repeats Markers from a Full-Length cDNA Library in Perilla frutescens (L.)

Perilla frutescens is valuable as a medicinal plant as well as a natural medicine and functional food. However, comparative genomics analyses of P. frutescens are limited due to a lack of gene annotations and characterization. A full-length cDNA library from P. frutescens leaves was constructed to identify functional gene clusters and probable EST-SSR markers via analysis of 1,056 expressed sequence tags. Unigene assembly was performed using basic local alignment search tool (BLAST) homology searches and annotated Gene Ontology (GO). A total of 18 simple sequence repeats (SSRs) were designed as primer pairs. This study is the first to report comparative genomics and EST-SSR markers from P. frutescens will help gene discovery and provide an important source for functional genomics and molecular genetic research in this interesting medicinal plant.

Dynamics of the chili pepper transcriptome during fruit development

BACKGROUND

The set of all mRNA molecules present in a cell constitute the transcriptome. The transcriptome varies depending on cell type as well as in response to internal and external stimuli during development. Here we present a study of the changes that occur in the transcriptome of chili pepper fruit during development and ripening.

RESULTS

RNA-Seq was used to obtain transcriptomes of whole Serrano-type chili pepper fruits (Capsicum annuum L.; 'Tampiqueño 74') collected at 10, 20, 40 and 60 days after anthesis (DAA). 15,550,468 Illumina MiSeq reads were assembled de novo into 34,066 chili genes. We classified the expression patterns of individual genes as well as genes grouped into Biological Process ontologies and Metabolic Pathway categories using statistical criteria. For the analyses of gene groups we added the weighted expression of individual genes. This method was effective in interpreting general patterns of expression changes and increased the statistical power of the analyses. We also estimated the variation in diversity and specialization of the transcriptome during chili pepper development. Approximately 17% of genes exhibited a significant change of expression in at least one of the intervals sampled. In contrast, significant differences in approximately 63% of the Biological Processes and 80% of the Metabolic Pathways studied were detected in at least one interval. Confirming previous reports, genes related to capsaicinoid and ascorbic acid biosynthesis were significantly upregulated at 20 DAA while those related to carotenoid biosynthesis were highly expressed in the last period of fruit maturation (40-60 DAA). Our RNA-Seq data was validated by examining the expression of nine genes involved in carotenoid biosynthesis by qRT-PCR.

CONCLUSIONS

In general, more profound changes in the chili fruit transcriptome were observed in the intervals between 10 to 20 and 40 to 60 DAA. The last interval, between 40 to 60 DAA, included 49% of all significant changes detected, and was characterized predominantly by a global decrease in gene expression. This period signals the end of maturation and the beginning of senescence of chili pepper fruit. The transcriptome at 60 DAA was the most specialized and least diverse of the four states sampled.

Conserved changes in the dynamics of metabolic processes during fruit development and ripening across species

Computational analyses of molecular phenotypes traditionally aim at identifying biochemical components that exhibit differential expression under various scenarios (e.g. environmental and internal perturbations) in a single species. High-throughput metabolomics technologies allow the quantification of (relative) metabolite levels across developmental stages in different tissues, organs, and species. Novel methods for analyzing the resulting multiple data tables could reveal preserved dynamics of metabolic processes across species. The problem we address in this study is 2-fold. (1) We derive a single data table, referred to as a compromise, which captures information common to the investigated set of multiple tables containing data on different fruit development and ripening stages in three climacteric (i.e. peach [Prunus persica] and two tomato [Solanum lycopersicum] cultivars, Ailsa Craig and M82) and two nonclimacteric (i.e. strawberry [Fragaria × ananassa] and pepper [Capsicum chilense]) fruits; in addition, we demonstrate the power of the method to discern similarities and differences between multiple tables by analyzing publicly available metabolomics data from three tomato ripening mutants together with two tomato cultivars. (2) We identify the conserved dynamics of metabolic processes, reflected in the data profiles of the corresponding metabolites that contribute most to the determined compromise. Our analysis is based on an extension to principal component analysis, called STATIS, in combination with pathway overenrichment analysis. Based on publicly available metabolic profiles for the investigated species, we demonstrate that STATIS can be used to identify the metabolic processes whose behavior is similarly affected during fruit development and ripening. These findings ultimately provide insights into the pathways that are essential during fruit development and ripening across species.

De novo transcriptome assembly and novel microsatellite marker information in Capsicum annuum varieties Saengryeg 211 and Saengryeg 213

BACKGROUND

Pepper, Capsicum annuum L., Solanaceae, is a major staple economically important vegetable crop worldwide. Limited functional genomics resources and whole genome association studies could be substantially improved through the application of molecular approach for the characterization of gene content and identification of molecular markers. The massive parallel pyrosequencing of two pepper varieties, the highly pungent, Saengryeg 211, and the non-pungent, Saengryeg 213, including de novo transcriptome assembly, functional annotation, and in silico discovery of potential molecular markers is described. We performed 454 GS-FLX Titanium sequencing of polyA-selected and normalized cDNA libraries generated from a single pool of transcripts obtained from mature fruits of two pepper varieties.

RESULTS

A single 454 pyrosequencing run generated 361,671 and 274,269 reads totaling 164.49 and 124.60 Mb of sequence data (average read length of 454 nucleotides), which assembled into 23,821 and 17,813 isotigs and 18,147 and 15,129 singletons for both varieties, respectively. These reads were organized into 20,352 and 15,781 'isogroups' for both varieties. Assembled sequences were functionally annotated based on homology to genes in multiple public databases and assigned with Gene Ontology (GO) terms. Sequence variants analyses identified a total of 3,766 and 2,431 potential (Simple Sequence Repeat) SSR motifs for microsatellite analysis for both varieties, where trinucleotide was the most common repeat unit (84%), followed by di (9.9%), hexa (4.1%) and pentanucleotide repeats (2.1%). GAA repeat (8.6%) was the most frequent repeat motif, followed by TGG (7.2%), TTC (6.5%), and CAG (6.2%).

CONCLUSIONS

High-throughput transcriptome assembly, annotation and large scale of SSR marker discovery has been achieved using next generation sequencing (NGS) of two pepper varieties. These valuable informations for functional genomics resource shall help to further improve the pepper breeding efforts with respect to genetic linkage maps, QTL mapping and marker-assisted trait selection.

Transcriptome analysis of Capsicum annuum varieties Mandarin and Blackcluster: assembly, annotation and molecular marker discovery

Next generation sequencing technologies have proven to be a rapid and cost-effective means to assemble and characterize gene content and identify molecular markers in various organisms. Pepper (Capsicum annuum L., Solanaceae) is a major staple vegetable crop, which is economically important and has worldwide distribution. High-throughput transcriptome profiling of two pepper cultivars, Mandarin and Blackcluster, using 454 GS-FLX pyrosequencing yielded 279,221 and 316,357 sequenced reads with a total 120.44 and 142.54Mb of sequence data (average read length of 431 and 450 nucleotides). These reads resulted from 17,525 and 16,341 'isogroups' and were assembled into 19,388 and 18,057 isotigs, and 22,217 and 13,153 singletons for both the cultivars, respectively. Assembled sequences were annotated functionally based on homology to genes in multiple public databases. Detailed sequence variant analysis identified a total of 9701 and 12,741 potential SNPs which eventually resulted in 1025 and 1059 genotype specific SNPs, for both the varieties, respectively, after examining SNP frequency distribution for each mapped unigenes. These markers for pepper will be highly valuable for marker-assisted breeding and other genetic studies.

A pepper (Capsicum annuum L.) metacaspase 9 (Camc9) plays a role in pathogen-induced cell death in plants

Metacaspases, which belong to the cysteine-type C14 protease family, are most structurally similar to mammalian caspases than any other caspase-like protease in plants. Atmc9 (Arabidopsis thaliana metacaspase 9) has a unique domain structure, and distinct biochemical characteristics, such as Ca²⁺ binding, pH, redox status, S-nitrosylation and specific protease inhibitors. However, the biological roles of Atmc9 in plant-pathogen interactions remain largely unknown. In this study, a metacaspase gene present as a single copy in the pepper genome, and sharing 54% amino acid sequence identity with Atmc9, was isolated and named Capsicum annuum metacaspase 9 (Camc9). Camc9 encodes a 318-amino-acid polypeptide with an estimated molecular weight of 34.6 kDa, and shares approximately 40% amino acid sequence identity with known type II metacaspases in plants. Quantitative reverse transcription-polymerase chain reaction analyses revealed that the expression of Camc9 was induced by infections of Xanthomonas campestris pv. vesicatoria race 1 and race 3 and treatment with methyl jasmonate. Suppression of Camc9 expression using virus-induced gene silencing enhanced disease resistance and suppressed cell death symptom development following infection with virulent bacterial pathogens. By contrast, overexpression of Camc9 by transient or stable transformation enhanced disease susceptibility and pathogen-induced cell death by regulation of reactive oxygen species production and defence-related gene expression. These results suggest that Camc9 is a possible member of the metacaspase gene family and plays a role as a positive regulator of pathogen-induced cell death in the plant kingdom.

Multiple classes of immune-related proteases associated with the cell death response in pepper plants

Proteases regulate a large number of biological processes in plants, such as metabolism, physiology, growth, and defense. In this study, we carried out virus-induced gene silencing assays with pepper cDNA clones to elucidate the biological roles of protease superfamilies. A total of 153 representative protease genes from pepper cDNA were selected and cloned into a Tobacco rattle virus-ligation independent cloning vector in a loss-of-function study. Silencing of 61 proteases resulted in altered phenotypes, such as the inhibition of shoot growth, abnormal leaf shape, leaf color change, and lethality. Furthermore, the silencing experiments revealed that multiple proteases play a role in cell death and immune response against avirulent and virulent pathogens. Among these 153 proteases, 34 modulated the hypersensitive cell death response caused by infection with an avirulent pathogen, and 16 proteases affected disease symptom development caused by a virulent pathogen. Specifically, we provide experimental evidence for the roles of multiple protease genes in plant development and immune defense following pathogen infection. With these results, we created a broad sketch of each protease function. This information will provide basic information for further understanding the roles of the protease superfamily in plant growth, development, and defense.

Characterization of Capsicum annuum genetic diversity and population structure based on parallel polymorphism discovery with a 30K unigene Pepper GeneChip

The widely cultivated pepper, Capsicum spp., important as a vegetable and spice crop world-wide, is one of the most diverse crops. To enhance breeding programs, a detailed characterization of Capsicum diversity including morphological, geographical and molecular data is required. Currently, molecular data characterizing Capsicum genetic diversity is limited. The development and application of high-throughput genome-wide markers in Capsicum will facilitate more detailed molecular characterization of germplasm collections, genetic relationships, and the generation of ultra-high density maps. We have developed the Pepper GeneChip® array from Affymetrix for polymorphism detection and expression analysis in Capsicum. Probes on the array were designed from 30,815 unigenes assembled from expressed sequence tags (ESTs). Our array design provides a maximum redundancy of 13 probes per base pair position allowing integration of multiple hybridization values per position to detect single position polymorphism (SPP). Hybridization of genomic DNA from 40 diverse C. annuum lines, used in breeding and research programs, and a representative from three additional cultivated species (C. frutescens, C. chinense and C. pubescens) detected 33,401 SPP markers within 13,323 unigenes. Among the C. annuum lines, 6,426 SPPs covering 3,818 unigenes were identified. An estimated three-fold reduction in diversity was detected in non-pungent compared with pungent lines, however, we were able to detect 251 highly informative markers across these C. annuum lines. In addition, an 8.7 cM region without polymorphism was detected around Pun1 in non-pungent C. annuum. An analysis of genetic relatedness and diversity using the software Structure revealed clustering of the germplasm which was confirmed with statistical support by principle components analysis (PCA) and phylogenetic analysis. This research demonstrates the effectiveness of parallel high-throughput discovery and application of genome-wide transcript-based markers to assess genetic and genomic features among Capsicum annuum.

Isolation and characterization of CaMF3, an anther-specific gene in Capsicum annuum L

Previous work on gene expression analysis based on RNA sequencing identified a variety of differentially expressed cDNA fragments in the genic male sterile-fertile line 114AB of Capsicum annuum L. In this work, we examined the accumulation of one of the transcript-derived fragments (TDFs), CaMF3 (male fertile 3), in the flower buds of a fertile line. The full genomic DNA sequence of CaMF3 was 1,951 bp long and contained 6 exons and 5 introns, with the complete sequence encoding a putative 25.89 kDa protein of 234 amino acids. The predicted protein of CaMF3 shared sequence similarity with members of the isoamyl acetate-hydrolyzing esterase (IAH1) protein family. CaMF3 expression was detected only in flower buds at stages 7 and 8 and in open flowers of a male fertile line; no expression was observed in any organs of a male sterile line. Fine expression analysis revealed that CaMF3 was expressed specifically in anthers of the fertile line. These results suggest that CaMF3 is an anther-specific gene that may be essential for anther or pollen development in C. annuum.

A common plant cell-wall protein HyPRP1 has dual roles as a positive regulator of cell death and a negative regulator of basal defense against pathogens

Although hybrid proline-rich proteins (HyPRPs) are ubiquitous in plants, little is known about their roles other than as cell-wall structural proteins. We identified the gene HyPRP1 in Capsicum annuum and Nicotiana benthamiana, which encodes a protein containing proline-rich domain and eight-cysteine motif (8CM) that is constitutively expressed in various organs, mostly in the root, but is down-regulated upon inoculation with either incompatible or compatible pathogens. Ectopic expression of HyPRP1 in plants accelerated cell death, showing developmental abnormality with down-regulation of ROS-scavenging genes, and enhanced pathogen susceptibility suppressing expression of defense-related genes. Conversely, silencing of HyPRP1 suppressed pathogen-induced cell death, but enhanced disease resistance, with up-regulation of defense-related genes and inhibition of in planta growth of bacterial pathogens independently of signal molecule-mediated pathways. Furthermore, the secreted 8CM was sufficient for these HyPRP1 functions. Together, our results suggest that a common plant cell-wall structural protein, HyPRP1, performs distinct dual roles in positive regulation of cell death and negative regulation of basal defense against pathogen.

The capsicum transcriptome DB: a "hot" tool for genomic research

Chili pepper (Capsicum annuum) is an economically important crop with no available public genome sequence. We describe a genomic resource to facilitate Capsicum annuum research. A collection of Expressed Sequence Tags (ESTs) derived from five C. annuum organs (root, stem, leaf, flower and fruit) were sequenced using the Sanger method and multiple leaf transcriptomes were deeply sampled using with GS-pyrosequencing. A hybrid assembly of 1,324,516 raw reads yielded 32,314 high quality contigs as validated by coverage and identity analysis with existing pepper sequences. Overall, 75.5% of the contigs had significant sequence similarity to entries in nucleic acid and protein databases; 23% of the sequences have not been previously reported for C. annuum and expand sequence resources for this species. A MySQL database and a user-friendly Web interface were constructed with search-tools that permit queries of the ESTs including sequence, functional annotation, Gene Ontology classification, metabolic pathways, and assembly information. The Capsicum Transcriptome DB is free available from http://www.bioingenios.ira.cinvestav.mx:81/Joomla/

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

Identification of reference genes for reverse transcription quantitative real-time PCR normalization in pepper (Capsicum annuum L.)

Reverse transcription quantitative real-time polymerase chain reaction (qRT-PCR) is a commonly used technology for gene expression and transcriptome analysis. Normalization is a process that is necessary to accurately analyze qRT-PCR data. Stability of reference gene expression is required for this process. Due to the large variation in expression levels of reference genes obtained from different experimental conditions, gene expression stabilities must be evaluated and identified in all experimental systems. In the present paper, the stability of the expression levels of seven potential reference genes in pepper are assessed using qRT-PCR analysis to determine optimal reference genes. These reference genes are evaluated in different pepper tissues, abiotic stress, and hormonal treatment samples. Three common statistical algorithms, geNorm, NormFinder, and BestKeeper, are used to identify expression stability and provide an accurate selection of reference genes. Two reference genes, beta tubulin and ubiquitin-conjugating protein (UBI-3), showed high stability in sample pools with abiotic stress and hormonal treatments. Among the sample pools tested, UBI-3 and glyceraldehyde-3-phosphate dehydrogenase expression levels were the most stable in different tissues. Therefore, these reference genes are selected for qRT-PCR analysis under the experimental conditions tested in pepper. In contrast, ubiquitin-conjugating enzyme and actin genes are identified as the least stable reference genes in all the groups tested, confirming that they are not suitable for normalization. Validation of these candidate genes could provide useful guidelines for reference gene selection in qRT-PCR studies in pepper.

Transcriptome profiling and molecular marker discovery in red pepper, Capsicum annuum L. TF68

Transcriptome from high throughput sequencing-by-synthesis is a good resource of molecular markers. In this study, we present utility of massively parallel sequencing by synthesis for profiling the transcriptome of red pepper (Capsicum annuum L. TF68) using 454 GS-FLX pyrosequencing. Through the generation of approximately 30.63 megabases (Mb) of expressed sequence tag (EST) data with the average length of 375 base pairs (bp), 9,818 contigs and 23,712 singletons were obtained by raw reads assembly. Using BLAST alignment against NCBI non-redundant and a UniProt protein database, 30% of the tentative consensus sequences were assigned to specific function annotation, while 24% returned alignments of unknown function, leaving up to 46% with no alignment. Functional classification using FunCat revealed that sequences with putative known function were distributed cross 18 categories. All unigenes have an approximately equal distribution on chromosomes by aligning with tomato (Solanum lycopersicum) pseudomolecules. Furthermore, 1,536 high quality single nucleotide discrepancies were discovered using the Bukang mature fruit cDNA collection (dbEST ID: 23667) as a reference. Moreover, 758 simple sequence repeat (SSR) motif loci were mined from 614 contigs, from which 572 primer sets were designed. The SSR motifs corresponded to di- and tri- nucleotide motifs (27.03 and 61.92%, respectively). These molecular markers may be of great value for application in linkage mapping and association mapping research.

In silico identification and characterization of microRNAs and their putative target genes in Solanaceae plants

MicroRNAs (miRNAs) are a class of small, single-stranded, noncoding RNAs ranging from 19 to 25 nucleotides. The miRNA control various cellular functions by negatively regulating gene expression at the post-transcriptional level. The miRNA regulation over their target genes has a central role in regulating plant growth and development; however, only a few reports have been published on the function of miRNAs in the family Solanaceae. We identified Solanaceae miRNAs and their target genes by analyzing expressed sequence tag (EST) data from five different Solanaceae species. A comprehensive bioinformatic analysis of EST data of Solanaceae species revealed the presence of at least 11 miRNAs and 54 target genes in pepper (Capsicum annuum L.), 22 miRNAs and 221 target genes in potato (Solanum tuberosum L.), 12 miRNAs and 417 target genes in tomato (Solanum lycopersicum L.), 46 miRNAs and 60 target genes in tobacco (Nicotiana tabacum L.), and 7 miRNAs and 28 target genes in Nicotiana benthamiana. The identified Solanaceae miRNAs and their target genes were deposited in the SolmiRNA database, which is freely available for academic research only at http://genepool.kribb.re.kr/SolmiRNA. Our data indicate that the Solanaceae family has both conserved and specific miRNAs and that their target genes may play important roles in growth and development of Solanaceae plants.

Complete sequencing and comparative analyses of the pepper (Capsicum annuum L.) plastome revealed high frequency of tandem repeats and large insertion/deletions on pepper plastome

Plants in the family Solanaceae are used as model systems in comparative and evolutionary genomics. The complete chloroplast genomes of seven solanaceous species have been sequenced, including tobacco, potato and tomato, but not peppers. We analyzed the complete chloroplast genome sequence of the hot pepper, Capsicum annuum. The pepper chloroplast genome was 156,781 bp in length, including a pair of inverted repeats (IR) of 25,783 bp. The content and the order of 133 genes in the pepper chloroplast genome were identical to those of other solanaceous plastomes. To characterize pepper plastome sequence, we performed comparative analysis using complete plastome sequences of pepper and seven solanaceous plastomes. Frequency and contents of large indels and tandem repeat sequences and distribution pattern of genome-wide sequence variations were investigated. In addition, a phylogenetic analysis using concatenated alignments of coding sequences was performed to determine evolutionary position of pepper in Solanaceae. Our results revealed two distinct features of pepper plastome compared to other solanaceous plastomes. Firstly, large indels, including insertions on accD and rpl20 gene sequences, were predominantly detected in the pepper plastome compared to other solanaceous plastomes. Secondly, tandem repeat sequences were particularly frequent in the pepper plastome. Taken together, our study represents unique features of evolution of pepper plastome among solanaceous plastomes.

SolEST database: a "one-stop shop" approach to the study of Solanaceae transcriptomes

BACKGROUND

Since no genome sequences of solanaceous plants have yet been completed, expressed sequence tag (EST) collections represent a reliable tool for broad sampling of Solanaceae transcriptomes, an attractive route for understanding Solanaceae genome functionality and a powerful reference for the structural annotation of emerging Solanaceae genome sequences.

DESCRIPTION

We describe the SolEST database http://biosrv.cab.unina.it/solestdb which integrates different EST datasets from both cultivated and wild Solanaceae species and from two species of the genus Coffea. Background as well as processed data contained in the database, extensively linked to external related resources, represent an invaluable source of information for these plant families. Two novel features differentiate SolEST from other resources: i) the option of accessing and then visualizing Solanaceae EST/TC alignments along the emerging tomato and potato genome sequences; ii) the opportunity to compare different Solanaceae assemblies generated by diverse research groups in the attempt to address a common complaint in the SOL community.

CONCLUSION

Different databases have been established worldwide for collecting Solanaceae ESTs and are related in concept, content and utility to the one presented herein. However, the SolEST database has several distinguishing features that make it appealing for the research community and facilitates a "one-stop shop" for the study of Solanaceae transcriptomes.

Basic properties and information theory of Audic-Claverie statistic for analyzing cDNA arrays

BACKGROUND

The Audic-Claverie method 1 has been and still continues to be a popular approach for detection of differentially expressed genes in the SAGE framework. The method is based on the assumption that under the null hypothesis tag counts of the same gene in two libraries come from the same but unknown Poisson distribution. The problem is that each SAGE library represents only a single measurement. We ask: Given that the tag count samples from SAGE libraries are extremely limited, how useful actually is the Audic-Claverie methodology? We rigorously analyze the A-C statistic that forms a backbone of the methodology and represents our knowledge of the underlying tag generating process based on one observation.

RESULTS

We show that the A-C statistic and the underlying Poisson distribution of the tag counts share the same mode structure. Moreover, the K-L divergence from the true unknown Poisson distribution to the A-C statistic is minimized when the A-C statistic is conditioned on the mode of the Poisson distribution. Most importantly, the expectation of this K-L divergence never exceeds 1/2 bit.

CONCLUSION

A rigorous underpinning of the Audic-Claverie methodology has been missing. Our results constitute a rigorous argument supporting the use of Audic-Claverie method even though the SAGE libraries represent very sparse samples.