Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species

A HapMap leads to a Capsicum annuum SNP infinium array: a new tool for pepper breeding

The Capsicum genus (Pepper) is a part of the Solanacae family. It has been important in many cultures worldwide for its key nutritional components and uses as spices, medicines, ornamentals and vegetables. Worldwide population growth is associated with demand for more nutritionally valuable vegetables while contending with decreasing resources and available land. These conditions require increased efficiency in pepper breeding to deal with these imminent challenges. Through resequencing of inbred lines we have completed a valuable haplotype map (HapMap) for the pepper genome based on single-nucleotide polymorphisms (SNP). The identified SNPs were annotated and classified based on their gene annotation in the pepper draft genome sequence and phenotype of the sequenced inbred lines. A selection of one marker per gene model was utilized to create the PepperSNP16K array, which simultaneously genotyped 16 405 SNPs, of which 90.7% were found to be informative. A set of 84 inbred and hybrid lines and a mapping population of 90 interspecific F2 individuals were utilized to validate the array. Diversity analysis of the inbred lines shows a distinct separation of bell versus chile/hot pepper types and separates them into five distinct germplasm groups. The interspecific population created between Tabasco (C. frutescens chile type) and P4 (C. annuum blocky type) produced a linkage map with 5546 markers separated into 1361 bins on twelve 12 linkage groups representing 1392.3 cM. This publically available genotyping platform can be used to rapidly assess a large number of markers in a reproducible high-throughput manner for pepper. As a standardized tool for genetic analyses, the PepperSNP16K can be used worldwide to share findings and analyze QTLs for important traits leading to continued improvement of pepper for consumers. Data and information on the array are available through the Solanaceae Genomics Network.

Genome-Wide Identification and Expression Analysis of WRKY Gene Family in Capsicum annuum L

The WRKY family of transcription factors is one of the most important families of plant transcriptional regulators with members regulating multiple biological processes, especially in regulating defense against biotic and abiotic stresses. However, little information is available about WRKYs in pepper (Capsicum annuum L.). The recent release of completely assembled genome sequences of pepper allowed us to perform a genome-wide investigation for pepper WRKY proteins. In the present study, a total of 71 WRKY genes were identified in the pepper genome. According to structural features of their encoded proteins, the pepper WRKY genes (CaWRKY) were classified into three main groups, with the second group further divided into five subgroups. Genome mapping analysis revealed that CaWRKY were enriched on four chromosomes, especially on chromosome 1, and 15.5% of the family members were tandemly duplicated genes. A phylogenetic tree was constructed depending on WRKY domain' sequences derived from pepper and Arabidopsis. The expression of 21 selected CaWRKY genes in response to seven different biotic and abiotic stresses (salt, heat shock, drought, Phytophtora capsici, SA, MeJA, and ABA) was evaluated by quantitative RT-PCR; Some CaWRKYs were highly expressed and up-regulated by stress treatment. Our results will provide a platform for functional identification and molecular breeding studies of WRKY genes in pepper.

Genomics-based strategies for the use of natural variation in the improvement of crop metabolism

Next-generation genomics holds great potential in the study of plant phenotypic variation. With several crop reference genomes now available, the affordable costs of de novo genome assembly or target resequencing offer the opportunity to mine the enormous amount of genetic diversity hidden in crop wild relatives. Wide introgressions from these wild ancestors species or land races represent a possible strategy to improve cultivated varieties. In this review, we discuss the mechanisms underlying metabolic diversity within plant species and the possible strategies (and barriers) to introgress novel metabolic traits into cultivated varieties. We show how deep genomic surveys uncover various types of structural variants from extended gene pools of major crops and highlight how this variation may be used for the improvement of crop metabolism.

Genome-Wide Divergence and Linkage Disequilibrium Analyses for Capsicum baccatum Revealed by Genome-Anchored Single Nucleotide Polymorphisms

Principal component analysis (PCA) with 36,621 polymorphic genome-anchored single nucleotide polymorphisms (SNPs) identified collectively for Capsicum annuum and Capsicum baccatum was used to characterize population structure and species domestication of these two important incompatible cultivated pepper species. Estimated mean nucleotide diversity (π) and Tajima's D across various chromosomes revealed biased distribution toward negative values on all chromosomes (except for chromosome 4) in cultivated C. baccatum, indicating a population bottleneck during domestication of C. baccatum. In contrast, C. annuum chromosomes showed positive π and Tajima's D on all chromosomes except chromosome 8, which may be because of domestication at multiple sites contributing to wider genetic diversity. For C. baccatum, 13,129 SNPs were available, with minor allele frequency (MAF) ≥0.05; PCA of the SNPs revealed 283 C. baccatum accessions grouped into 3 distinct clusters, for strong population structure. The fixation index (F_ST ) between domesticated C. annuum and C. baccatum was 0.78, which indicates genome-wide divergence. We conducted extensive linkage disequilibrium (LD) analysis of C. baccatum var. pendulum cultivars on all adjacent SNP pairs within a chromosome to identify regions of high and low LD interspersed with a genome-wide average LD block size of 99.1 kb. We characterized 1742 haplotypes containing 4420 SNPs (range 9-2 SNPs per haplotype). Genome-wide association study (GWAS) of peduncle length, a trait that differentiates wild and domesticated C. baccatum types, revealed 36 significantly associated genome-wide SNPs. Population structure, identity by state (IBS) and LD patterns across the genome will be of potential use for future GWAS of economically important traits in C. baccatum peppers.

New Insights on Eggplant/Tomato/Pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL

Eggplant, pepper, and tomato are the most exploited berry-producing vegetables within the Solanaceae family. Their genomes differ in size, but each has 12 chromosomes which have undergone rearrangements causing a redistribution of loci. The genome sequences of all three species are available but differ in coverage, assembly quality and percentage of anchorage. Determining their syntenic relationship and QTL orthology will contribute to exploit genomic resources and genetic data for key agronomic traits. The syntenic analysis between tomato and pepper based on the alignment of 34,727 tomato CDS to the pepper genome sequence, identified 19,734 unique hits. The resulting synteny map confirmed the 14 inversions and 10 translocations previously documented, but also highlighted 3 new translocations and 4 major new inversions. Furthermore, each of the 12 chromosomes exhibited a number of rearrangements involving small regions of 0.5-0.7 Mbp. Due to high fragmentation of the publicly available eggplant genome sequence, physical localization of most eggplant QTL was not possible, thus, we compared the organization of the eggplant genetic map with the genome sequence of both tomato and pepper. The eggplant/tomato syntenic map confirmed all the 10 translocations but only 9 of the 14 known inversions; on the other hand, a newly detected inversion was recognized while another one was not confirmed. The eggplant/pepper syntenic map confirmed 10 translocations and 8 inversions already detected and suggested a putative new translocation. In order to perform the assessment of eggplant and pepper QTL orthology, the eggplant and pepper sequence-based markers located in their respective genetic map were aligned onto the pepper genome. GBrowse in pepper was used as reference platform for QTL positioning. A set of 151 pepper QTL were located as well as 212 eggplant QTL, including 76 major QTL (PVE ≥ 10%) affecting key agronomic traits. Most were confirmed to cluster in orthologous chromosomal regions. Our results highlight that the availability of genome sequences for an increasing number of crop species and the development of "ultra-dense" physical maps provide new and key tools for detailed syntenic and orthology studies between related plant species.

Genome-Wide Characterization of the MADS-Box Gene Family in Radish (Raphanus sativus L.) and Assessment of Its Roles in Flowering and Floral Organogenesis

The MADS-box gene family is an important transcription factor (TF) family that is involved in various aspects of plant growth and development, especially flowering time and floral organogenesis. Although it has been reported in many plant species, the systematic identification and characterization of MADS-box TF family is still limited in radish (Raphanus sativus L.). In the present study, a comprehensive analysis of MADS-box genes was performed, and a total of 144 MADS-box family members were identified from the whole radish genome. Meanwhile, a detailed list of MADS-box genes from other 28 plant species was also investigated. Through the phylogenetic analysis between radish and Arabidopsis thaliana, all the RsMADS genes were classified into two groups including 68 type I (31 Mα, 12 Mβ and 25Mγ) and 76 type II (70 MIKC^C and 6 MIKC^∗). Among them, 41 (28.47%) RsMADS genes were located in nine linkage groups of radish from R1 to R9. Moreover, the homologous MADS-box gene pairs were identified among radish, A. thaliana, Chinese cabbage and rice. Additionally, the expression profiles of RsMADS genes were systematically investigated in different tissues and growth stages. Furthermore, quantitative real-time PCR analysis was employed to validate expression patterns of some crucial RsMADS genes. These results could provide a valuable resource to explore the potential functions of RsMADS genes in radish, and facilitate dissecting MADS-box gene-mediated molecular mechanisms underlying flowering and floral organogenesis in root vegetable crops.

Autophagy, a Conserved Mechanism for Protein Degradation, Responds to Heat, and Other Abiotic Stresses in Capsicum annuum L

Abiotic stresses negatively affect plants growth and development by inducing protein denaturation, and autophagy degrades the damaged proteins to alleviate their toxicity, however, little is known about the involvement of autophagy in pepper (Capsicum annuum L.) tolerances to abiotic stresses. In this study, we identified autophagy-related gene (ATG) members in the whole genome of pepper by HMM method and analyzed their expression profiles in response to heat and other abiotic stresses by quantitative real-time PCR. The results showed that the CaATG contained 15 core ATG members including 29 ATG proteins with their respective conserved functional domains, involving the whole process of autophagy. Under normal environmental condition, the expression of CaATG genes showed tissue- and developmental stage-specific patterns, while under abiotic stresses of salt, drought, heat, cold and carbohydrate starvation, the accumulation of autophagosome punctate increased and the expression level of CaATG genes changed with stress type-dependent pattern, which indicates the linkage of autophagy in pepper response to abiotic stresses. After treated with heat stress, both the number of up-regulated CaATG genes and the increment of autophagosome punctate were higher in pepper thermotolerant line R9 than those in thermosensitive line B6, implying an association of autophagy with heat tolerance. In addition, CaATG6 was predicted to interact with CaHSP90 family members. Our study suggests that autophagy is connected to pepper tolerances to heat and other abiotic stresses.

Genome-Wide Identification and Expression Profile of Dof Transcription Factor Gene Family in Pepper (Capsicum annuum L.)

Dof (DNA-binding One Zinc Finger) transcription factor family is unique to plants and has diverse roles associated with plant-specific phenomena, such as light, phytohormone and defense responses as well as seed development and germination. Although, genome-wide analysis of this family has been performed in many species, information regarding Dof genes in the pepper, Capsicum annuum L., is extremely limited. In this study, exhaustive searches of pepper genome revealed 33 potential CaDofs that were phylogenetically clustered into four subgroups. Twenty-nine of the 33 Dof genes could be mapped on 11 chromosomes, except for chromosome 7. The intron/exon organizations and conserved motif compositions of these genes were also analyzed. Additionally, phylogenetic analysis and classification of the Dof transcription factor family in eight plant species revealed that S. lycopersicum and C. annuum as well as O. sativa and S. bicolor Dof proteins may have evolved conservatively. Moreover, comprehensive expression analysis of CaDofs using a RNA-seq atlas and quantitative real-time polymerase chain reaction (qRT-PCR) revealed that these genes exhibit a variety of expression patterns. Most of the CaDofs were expressed in at least one of the tissues tested, whereas several genes were identified as being highly responsive to heat and salt stresses. Overall, this study describes the first genome-wide analysis of the pepper Dof family, whose genes exhibited different expression patterns in all primary fruit developmental stages and tissue types, as in response to abiotic stress. In particular, some Dof genes might be used as biomarkers for heat and salt stress. The results could expand our understanding of the roles of Dof genes in pepper.

Major Quantitative Trait Loci and Putative Candidate Genes for Powdery Mildew Resistance and Fruit-Related Traits Revealed by an Intraspecific Genetic Map for Watermelon (Citrullus lanatus var. lanatus)

An intraspecific genetic map for watermelon was constructed using an F₂ population derived from ‘Arka Manik’ × ‘TS34’ and transcript sequence variants and quantitative trait loci (QTL) for resistance to powdery mildew (PMR), seed size (SS), and fruit shape (FS) were analyzed. The map consists of 14 linkage groups (LGs) defined by 174 cleaved amplified polymorphic sequences (CAPS), 2 derived-cleaved amplified polymorphic sequence markers, 20 sequence-characterized amplified regions, and 8 expressed sequence tag-simple sequence repeat markers spanning 1,404.3 cM, with a mean marker interval of 6.9 cM and an average of 14.6 markers per LG. Genetic inheritance and QTL analyses indicated that each of the PMR, SS, and FS traits is controlled by an incompletely dominant effect of major QTLs designated as pmr2.1, ss2.1, and fsi3.1, respectively. The pmr2.1, detected on chromosome 2 (Chr02), explained 80.0% of the phenotypic variation (LOD = 30.76). This QTL was flanked by two CAPS markers, wsb2-24 (4.00 cM) and wsb2-39 (13.97 cM). The ss2.1, located close to pmr2.1 and CAPS marker wsb2-13 (1.00 cM) on Chr02, explained 92.3% of the phenotypic variation (LOD = 68.78). The fsi3.1, detected on Chr03, explained 79.7% of the phenotypic variation (LOD = 31.37) and was flanked by two CAPS, wsb3-24 (1.91 cM) and wsb3-9 (7.00 cM). Candidate gene-based CAPS markers were developed from the disease resistance and fruit shape gene homologs located on Chr.02 and Chr03 and were mapped on the intraspecific map. Colocalization of these markers with the major QTLs indicated that watermelon orthologs of a nucleotide-binding site-leucine-rich repeat class gene containing an RPW8 domain and a member of SUN containing the IQ67 domain are candidate genes for pmr2.1 and fsi3.1, respectively. The results presented herein provide useful information for marker-assisted breeding and gene cloning for PMR and fruit-related traits.

Host genetic resistance to root-knot nematodes, Meloidogyne spp., in Solanaceae: from genes to the field

Root-knot nematodes (RKNs) heavily damage most solanaceous crops worldwide. Fortunately, major resistance genes are available in a number of plant species, and their use provides a safe and economically relevant strategy for RKN control. From a structural point of view, these genes often harbour NBS-LRR motifs (i.e. a nucleotide binding site and a leucine rich repeat region near the carboxy terminus) and are organised in syntenic clusters in solanaceous genomes. Their introgression from wild to cultivated plants remains a challenge for breeders, although facilitated by marker-assisted selection. As shown with other pathosystems, the genetic background into which the resistance genes are introgressed is of prime importance to both the expression of the resistance and its durability, as exemplified by the recent discovery of quantitative trait loci conferring quantitative resistance to RKNs in pepper. The deployment of resistance genes at a large scale may result in the emergence and spread of virulent nematode populations able to overcome them, as already reported in tomato and pepper. Therefore, careful management of the resistance genes available in solanaceous crops is crucial to avoid significant reduction in the duration of RKN genetic control in the field. From that perspective, only rational management combining breeding and cultivation practices will allow the design and implementation of innovative, sustainable crop production systems that protect the resistance genes and maintain their durability.

Identification of Nicotiana benthamiana microRNAs and their targets using high throughput sequencing and degradome analysis

BACKGROUND

Nicotiana benthamiana is a widely used model plant species for research on plant-pathogen interactions as well as other areas of plant science. It can be easily transformed or agroinfiltrated, therefore it is commonly used in studies requiring protein localization, interaction, or plant-based systems for protein expression and purification. To discover and characterize the miRNAs and their cleaved target mRNAs in N. benthamiana, we sequenced small RNA transcriptomes and degradomes of two N. benthamiana accessions and validated them by Northern blots.

RESULTS

We used a comprehensive molecular approach to detect and to experimentally validate N. benthamiana miRNAs and their target mRNAs from various tissues. We identified 40 conserved miRNA families and 18 novel microRNA candidates and validated their target mRNAs with a genomic scale approach. The accumulation of thirteen novel miRNAs was confirmed by Northern blot analysis. The conserved and novel miRNA targets were found to be involved in various biological processes including transcription, RNA binding, DNA modification, signal transduction, stress response and metabolic process. Among the novel miRNA targets we found the mRNA of REPRESSOR OF SILENCING (ROS1). Regulation of ROS1 by a miRNA provides a new regulatory layer to reinforce transcriptional gene silencing by a post-transcriptional repression of ROS1 activity.

CONCLUSIONS

The identified conserved and novel miRNAs along with their target mRNAs also provides a tissue specific atlas of known and new miRNA expression and their cleaved target mRNAs of N. benthamiana. Thus this study will serve as a valuable resource to the plant research community that will be beneficial well into the future.

Evolutionarily Dynamic, but Robust, Targeting of Resistance Genes by the miR482/2118 Gene Family in the Solanaceae

Plants are exposed to pathogens around the clock. A common resistance response in plants upon pathogen detection is localized cell death. Given the irreversible nature of this response, multiple layers of negative regulation are present to prevent the untimely or misexpression of resistance genes. One layer of negative regulation is provided by a recently discovered microRNA (miRNA) gene family, miR482/2118. This family targets the transcripts of resistance genes in plants. We investigated the evolutionary history and specificity of this miRNA gene family within the Solanaceae. This plant family includes many important crop species, providing a set of well-defined resistance gene repertoires. Across 14 species from the Solanaceae, we identified eight distinct miR482/2118 gene family members. Our studies show conservation of miRNA type and number in the group of wild tomatoes and, to a lesser extent, throughout the Solanaceae. The eight orthologous miRNA gene clusters evolved under different evolutionary constraints, allowing for individual subfunctionalization of the miRNAs. Despite differences in the predicted targeting behavior of each miRNA, the miRNA-R-gene network is robust due to its high degree of interconnectivity and redundant targeting. Our data suggest that the miR482/2118 gene family acts as an evolutionary buffer for R-gene sequence diversity.

Metabolomic Characterization of Hot Pepper (Capsicum annuum "CM334") during Fruit Development

Non-targeted metabolomic analysis of hot pepper (Capsicum annuum "CM334") was performed at six development stages [16, 25, 36, 38, 43, and 48 days post-anthesis (DPA)] to analyze biochemical changes. Distinct distribution patterns were observed in the changes of metabolites, gene expressions, and antioxidant activities by early (16-25 DPA), breaker (36-38 DPA), and later (43-48 DPA) stages. In the early stages, glycosides of luteolin, apigenin, and quercetin, shikimic acid, γ-aminobutyric acid (GABA), and putrescine were highly distributed but gradually decreased over the breaker stage. At later stages, leucine, isoleucine, proline, phenylalanine, capsaicin, dihydrocapsaicin, and kaempferol glycosides were significantly increased. Pathway analysis revealed metabolite-gene interactions in the biosynthesis of amino acids, capsaicinoids, fatty acid chains, and flavonoids. The changes in antioxidant activity were highly reflective of alterations in metabolites. The present study could provide useful information about nutrient content at each stage of pepper cultivation.

Genomics and Evolution in Traditional Medicinal Plants: Road to a Healthier Life

Medicinal plants have long been utilized in traditional medicine and ethnomedicine worldwide. This review presents a glimpse of the current status of and future trends in medicinal plant genomics, evolution, and phylogeny. These dynamic fields are at the intersection of phytochemistry and plant biology and are concerned with the evolution mechanisms and systematics of medicinal plant genomes, origin and evolution of the plant genotype and metabolic phenotype, interaction between medicinal plant genomes and their environment, the correlation between genomic diversity and metabolite diversity, and so on. Use of the emerging high-end genomic technologies can be expanded from crop plants to traditional medicinal plants, in order to expedite medicinal plant breeding and transform them into living factories of medicinal compounds. The utility of molecular phylogeny and phylogenomics in predicting chemodiversity and bioprospecting is also highlighted within the context of natural-product-based drug discovery and development. Representative case studies of medicinal plant genome, phylogeny, and evolution are summarized to exemplify the expansion of knowledge pedigree and the paradigm shift to the omics-based approaches, which update our awareness about plant genome evolution and enable the molecular breeding of medicinal plants and the sustainable utilization of plant pharmaceutical resources.

High genetic abundance of Rpi-blb2/Mi-1.2/Cami gene family in Solanaceae

Background

Three NBS-LRR genes, Rpi-blb2, Mi-1.2, and Cami, constitute a very special plant resistance gene family. These genes confer resistance against 4 distantly related pathogen species in 3 different Solanaceae hosts. To characterize this noted resistance, we conducted a series of studies on this gene family.

Results

First, homologs of this gene family were identified in the pepper, tomato and potato genomes. This revealed a large variation in copy number within this gene family among species and a great divergence was found both between and within species. To gain more information pertaining to gene resistance within this family, 121 LRR regions were cloned in 16 different wild/cultivated potato accessions. Again, frequent copy number variations and a high level of divergence between homolog were observed common among accessions. The divergence within species was so high that it reaches the level of divergence between species. Also, frequent frameshift mutations and abundant gene conversion events were identified in these LRR regions.

Conclusions

Our findings suggest that this family harbors an unusually high level of genetic abundance, making it of particular interest. Together with other reported examples, our study also provides evidence that multi-resistance is a common trait in R gene families like this.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0493-z) contains supplementary material, which is available to authorized users.

Screening Genetic Resources of Capsicum Peppers in Their Primary Center of Diversity in Bolivia and Peru

For most crops, like Capsicum, their diversity remains under-researched for traits of interest for food, nutrition and other purposes. A small investment in screening this diversity for a wide range of traits is likely to reveal many traditional varieties with distinguished values. One objective of this study was to demonstrate, with Capsicum as model crop, the application of indicators of phenotypic and geographic diversity as effective criteria for selecting promising genebank accessions for multiple uses from crop centers of diversity. A second objective was to evaluate the expression of biochemical and agromorphological properties of the selected Capsicum accessions in different conditions. Four steps were involved: 1) Develop the necessary diversity by expanding genebank collections in Bolivia and Peru; 2) Establish representative subsets of ~100 accessions for biochemical screening of Capsicum fruits; 3) Select promising accessions for different uses after screening; and 4) Examine how these promising accessions express biochemical and agromorphological properties when grown in different environmental conditions. The Peruvian Capsicum collection now contains 712 accessions encompassing all five domesticated species (C. annuum, C. chinense, C. frutescens, C. baccatum, and C. pubescens). The collection in Bolivia now contains 487 accessions, representing all five domesticates plus four wild taxa (C. baccatum var. baccatum, C. caballeroi, C. cardenasii, and C. eximium). Following the biochemical screening, 44 Bolivian and 39 Peruvian accessions were selected as promising, representing wide variation in levels of antioxidant capacity, capsaicinoids, fat, flavonoids, polyphenols, quercetins, tocopherols, and color. In Peru, 23 promising accessions performed well in different environments, while each of the promising Bolivian accessions only performed well in a certain environment. Differences in Capsicum diversity and local contexts led to distinct outcomes in each country. In Peru, mild landraces with high values in health-related attributes were of interest to entrepreneurs. In Bolivia, wild Capsicum have high commercial demand.

Molecular marker development and genetic diversity exploration by RNA-seq in Platycodon grandiflorum

Platycodon grandiflorum, generally known as the bellflower or balloon flower, is the only species in the genus Platycodon of the family Campanulaceae. Platycodon plants have been traditionally used as a medicinal crop in East Asia for their antiphlogistic, antitussive, and expectorant properties. Despite these practical uses, marker-assisted selection and molecular breeding in platycodons have lagged due to the lack of genetic information on this genus. In this study, we performed RNA-seq analysis of three platycodon accessions to develop molecular markers and explore genetic diversity. First, genic simple sequence repeats (SSRs) were retrieved and compared; dinucleotide motifs were the most abundant repeats (39%-40%) followed by trinucleotide (25%-31%), tetranucleotide (1.5%-1.9%), and pentanucleotide (0.3%-1.0%) repeats. The result of in silico SSR analysis, three SSR markers were detected and showed possibility to distinguish three platycodon accessions. After several filtering procedures, 180 single nucleotide polymorphisms (SNPs) were used to design 40 cleaved amplified polymorphic sequence (CAPS) markers. Twelve of these PCR-based markers were validated as highly polymorphic and utilized to investigate genetic diversity in 21 platycodon accessions collected from various regions of South Korea. Collectively, the 12 markers yielded 35 alleles, with an average of 3 alleles per locus. Polymorphism information content (PIC) values ranged from 0.087 to 0.693, averaging 0.373 per locus. Since platycodon genetics have not been actively studied, the sequence information and the DNA markers generated from our research have the potential to contribute to further genetic improvements, genomic studies, and gene discovery in this genus.

Ultra-High Density, Transcript-Based Genetic Maps of Pepper Define Recombination in the Genome and Synteny Among Related Species

Our ability to assemble complex genomes and construct ultradense genetic maps now allows the determination of recombination rates, translocations, and the extent of genomic collinearity between populations, species, and genera. We developed two ultradense genetic linkage maps for pepper from single-position polymorphisms (SPPs) identified de novo with a 30,173 unigene pepper genotyping array. The Capsicum frutescens × C. annuum interspecific and the C. annuum intraspecific genetic maps were constructed comprising 16,167 and 3,878 unigene markers in 2108 and 783 genetic bins, respectively. Accuracies of marker groupings and orders are validated by the high degree of collinearity between the two maps. Marker density was sufficient to locate the chromosomal breakpoint resulting in the P1/P8 translocation between C. frutescens and C. annuum to a single bin. The two maps aligned to the pepper genome showed varying marker density along the chromosomes. There were extensive chromosomal regions with suppressed recombination and reduced intraspecific marker density. These regions corresponded to the pronounced nonrecombining pericentromeric regions in tomato, a related Solanaceous species. Similar to tomato, the extent of reduced recombination appears to be more pronounced in pepper than in other plant species. Alignment of maps with the tomato and potato genomes shows the presence of previously known translocations and a translocation event that was not observed in previous genetic maps of pepper.

Functional genomics of tomato: opportunities and challenges in post-genome NGS era

The Tomato Genome Sequencing Project represented a landmark venture in the history of sequencing projects where both Sanger's and next-generation sequencing (NGS) technologies were employed, and a highly accurate and one of the best assembled plant genomes along with a draft of the wild relative, Solanum pimpinellifolium, were released in 2012. However, the functional potential of the major portion of this newly generated resource is still undefined. The very first challenge before scientists working on tomato functional biology is to exploit this high-quality reference sequence for tapping of the wealth of genetic variants for improving agronomic traits in cultivated tomatoes. The sequence data generated recently by 150 Tomato Genome Consortium would further uncover the natural alleles present in different tomato genotypes. Therefore, we found it relevant to have a fresh outlook on tomato functional genomics in the context of application of NGS technologies in its post-genome sequencing phase. Herein, we provide an overview how NGS technologies vis-a-vis available reference sequence have assisted each other for their mutual improvement and how their combined use could further facilitate the development of other 'omics' tools, required to propel the Solanaceae research. Additionally, we highlight the challenges associated with the application of these cutting-edge technologies.

A call for benchmarking transposable element annotation methods

DNA derived from transposable elements (TEs) constitutes large parts of the genomes of complex eukaryotes, with major impacts not only on genomic research but also on how organisms evolve and function. Although a variety of methods and tools have been developed to detect and annotate TEs, there are as yet no standard benchmarks-that is, no standard way to measure or compare their accuracy. This lack of accuracy assessment calls into question conclusions from a wide range of research that depends explicitly or implicitly on TE annotation. In the absence of standard benchmarks, toolmakers are impeded in improving their tools, annotators cannot properly assess which tools might best suit their needs, and downstream researchers cannot judge how accuracy limitations might impact their studies. We therefore propose that the TE research community create and adopt standard TE annotation benchmarks, and we call for other researchers to join the authors in making this long-overdue effort a success.

VIGS approach reveals the modulation of anthocyanin biosynthetic genes by CaMYB in chili pepper leaves

The purple coloration of pepper leaves arises from the accumulation of anthocyanin. Three regulatory and 12 structural genes have been characterized for their involvement in the anthocyanin biosynthesis. Examination of the abundance of these genes in leaves showed that the majority of them differed between anthocyanin pigmented line Z1 and non-pigmented line A3. Silencing of the R2R3-MYB transcription factor CaMYB in pepper leaves of Z1 resulted in the loss of anthocyanin accumulation. Moreover, the expression of multiple genes was altered in the silenced leaves. The expression of MYC was significantly lower in CaMYB-silenced leaves, whereas WD40 showed the opposite pattern. Most structural genes including CHS, CHI, F3H, F3'5'H, DFR, ANS, UFGT, ANP, and GST were repressed in CaMYB-silenced foliage with the exception of PAL, C4H, and 4CL. These results indicated that MYB plays an important role in the regulation of anthocyanin biosynthetic related genes. Besides CaMYB silenced leaves rendered more sporulation of Phytophthora capsici Leonian indicating that CaMYB might be involved in the defense response to pathogens.

Molecular characterization of Pvr9 that confers a hypersensitive response to Pepper mottle virus (a potyvirus) in Nicotiana benthamiana

There are some R genes against potyviruses which were mapped in pepper. However, none of them has been characterized at the molecular level. In this study, we characterized Pvr9 which is an Rpi-blb2 ortholog from pepper and confers a hypersensitive response to Pepper mottle virus (PepMoV) in a transient expression system in Nicotiana benthamiana. This gene putatively encoded for 1298 amino acids and is located on pepper chromosome 6. PepMoV NIb was the elicitor of the Pvr9-mediated hypersensitive response. NIb from several other potyviruses also elicited the hypersensitive response. Inoculation of pepper with PepMoV resulted in a minor increase in Pvr9 transcription in the resistant cultivar CM334 and a slight down-regulation in the susceptible cultivar Floral Gem. The 5' upstream region of Pvr9 from cultivar CM334 had higher transcription activity than the region from cultivar Floral Gem. The cultivars CM334 and Floral Gem had non-functional Pvr9 homologs with loss-of-function mutations.

A Recently Evolved Alternative Splice Site in the BRANCHED1a Gene Controls Potato Plant Architecture

Amplification and diversification of transcriptional regulators that control development is a driving force of morphological evolution. A major source of protein diversity is alternative splicing, which leads to the generation of different isoforms from a single gene. The mechanisms and timing of intron evolution nonetheless remain unclear, and the functions of alternative splicing-generated protein isoforms are rarely studied. In Solanum tuberosum, the BRANCHED1a (BRC1a) gene encodes a TCP transcription factor that controls lateral shoot outgrowth. Here, we report the recent evolution in Solanum of an alternative splice site in BRC1a that leads to the generation of two BRC1a protein isoforms with distinct C-terminal regions, BRC1a(Long) and BRC1a(Short), encoded by unspliced and spliced mRNA, respectively. The BRC1a(Long) C-terminal region has a strong activation domain, whereas that of BRC1a(S) lacks an activation domain and is predicted to form an amphipathic helix, the H domain, which prevents protein nuclear targeting. BRC1a(Short) is thus mainly cytoplasmic, while BRC1a(Long) is mainly nuclear. BRC1a(Long) functions as a transcriptional activator, whereas BRC1a(Short) appears to have no transcriptional activity. Moreover, BRC1a(Short) can heterodimerize with BRC1a(Long) and act as a dominant-negative factor; it increases BRC1a(Long) concentration in cytoplasm and reduces its transcriptional activity. This alternative splicing mechanism is regulated by hormones and external stimuli that control branching. The evolution of a new alternative splicing site and a novel protein domain in Solanum BRC1a led to a multi-level mechanism of post-transcriptional and post-translational BRC1a regulation that effectively modulates its branch suppressing activity in response to environmental and endogenous cues.

Genome-wide analysis, expression profile of heat shock factor gene family (CaHsfs) and characterisation of CaHsfA2 in pepper (Capsicum annuum L.)

BACKGROUND

Heat shock factors (Hsfs) play crucial roles in plant developmental and defence processes. The production and quality of pepper (Capsicum annuum L.), an economically important vegetable crop, are severely reduced by adverse environmental stress conditions, such as heat, salt and osmotic stress. Although the pepper genome has been fully sequenced, the characterization of the Hsf gene family under abiotic stress conditions remains incomplete.

RESULTS

A total of 25 CaHsf members were identified in the pepper genome by bioinformatics analysis and PCR assays. They were grouped into three classes, CaHsfA, B and C, based on highly conserved Hsf domains, were distributed over 11 of 12 chromosomes, with none found on chromosome 11, and all of them, except CaHsfA5, formed a protein-protein interaction network. According to the RNA-seq data of pepper cultivar CM334, most CaHsf members were expressed in at least one tissue among root, stem, leaf, pericarp and placenta. Quantitative real-time PCR assays showed that all of the CaHsfs responded to heat stress (40 °C for 2 h), except CaHsfC1 in thermotolerant line R9 leaves, and that the expression patterns were different from those in thermosensitive line B6. Many CaHsfs were also regulated by salt and osmotic stresses, as well as exogenous Ca(2+), putrescine, abscisic acid and methyl jasmonate. Additionally, CaHsfA2 was located in the nucleus and had transcriptional activity, consistent with the typical features of Hsfs. Time-course expression profiling of CaHsfA2 in response to heat stress revealed differences in its expression level and pattern between the pepper thermosensitive line B6 and thermotolerant line R9.

CONCLUSIONS

Twenty-five Hsf genes were identified in the pepper genome and most of them responded to heat, salt, osmotic stress, and exogenous substances, which provided potential clues for further analyses of CaHsfs functions in various kinds of abiotic stresses and of corresponding signal transduction pathways in pepper.

CaAP2 transcription factor is a candidate gene for a flowering repressor and a candidate for controlling natural variation of flowering time in Capsicum annuum

The APETALA2 transcription factor homolog CaAP2 is a candidate gene for a flowering repressor in pepper, as revealed by induced-mutation phenotype, and a candidate underlying a major QTL controlling natural variation in flowering time. To decipher the genetic control of transition to flowering in pepper (Capsicum spp.) and determine the extent of gene function conservation compared to model species, we isolated and characterized several ethyl methanesulfonate (EMS)-induced mutants that vary in their flowering time compared to the wild type. In the present study, we report on the isolation of an early-flowering mutant that flowers after four leaves on the primary stem compared to nine leaves in the wild-type 'Maor'. By genetic mapping and sequencing of putative candidate genes linked to the mutant phenotype, we identified a member of the APETALA2 (AP2) transcription factor family, CaAP2, which was disrupted in the early-flowering mutant. CaAP2 is a likely ortholog of AP2 that functions as a repressor of flowering in Arabidopsis. To test whether CaAP2 has an effect on controlling natural variation in the transition to flowering in pepper, we performed QTL mapping for flowering time in a cross between early and late-flowering C. annuum accessions. We identified a major QTL in a region of chromosome 2 in which CaAP2 was the most significant marker, explaining 52 % of the phenotypic variation of the trait. Sequence comparison of the CaAP2 open reading frames in the two parents used for QTL mapping did not reveal significant variation. In contrast, significant differences in expression level of CaAP2 were detected between near-isogenic lines that differ for the flowering time QTL, supporting the putative function of CaAP2 as a major repressor of flowering in pepper.

Antagonistic Activities of Bacillus spp. Strains Isolated from Tidal Flat Sediment Towards Anthracnose Pathogens Colletotrichum acutatum and C. gloeosporioides in South Korea

Anthracnose is a fungal disease caused by Colletotrichum species that is detrimental to numerous plant species. Anthracnose control with fungicides has both human health and environmental safety implications. Despite increasing public concerns, fungicide use will continue in the absence of viable alternatives. There have been relatively less efforts to search antagonistic bacteria from mudflats harboring microbial diversity. A total of 420 bacterial strains were isolated from mudflats near the western sea of South Korea. Five bacterial strains, LB01, LB14, HM03, HM17, and LB15, were characterized as having antifungal properties in the presence of C. acutatum and C. gloeosporioides. The three Bacillus atrophaeus strains, LB14, HM03, and HM17, produced large quantities of chitinase and protease enzymes, whereas the B. amyloliquefaciens strain LB01 produced protease and cellulase enzymes. Two important antagonistic traits, siderophore production and solubilization of insoluble phosphate, were observed in the three B. atrophaeus strains. Analyses of disease suppression revealed that LB14 was most effective for suppressing the incidence of anthracnose symptoms on pepper fruits. LB14 produced antagonistic compounds and suppressed conidial germination of C. acutatum and C. gloeosporioides. The results from the present study will provide a basis for developing a reliable alternative to fungicides for anthracnose control.

Evolutionary analysis of RB/Rpi-blb1 locus in the Solanaceae family

Late blight caused by the oomycete Phytophthora infestans is one of the most severe threats to potato production worldwide. Numerous studies suggest that the most effective protective strategy against the disease would be to provide potato cultivars with durable resistance (R) genes. However, little is known about the origin and evolutional history of these durable R-genes in potato. Addressing this might foster better understanding of the dynamics of these genes in nature and provide clues for identifying potential candidate R-genes. Here, a systematic survey was executed at RB/Rpi-blb1 locus, an exclusive broad-spectrum R-gene locus in potato. As indicated by synteny analysis, RB/Rpi-blb1 homologs were identified in all tested genomes, including potato, tomato, pepper, and Nicotiana, suggesting that the RB/Rpi-blb1 locus has an ancient origin. Two evolutionary patterns, similar to those reported on RGC2 in Lactuca, and Pi2/9 in rice, were detected at this locus. Type I RB/Rpi-blb1 homologs have frequent copy number variations and sequence exchanges, obscured orthologous relationships, considerable nucleotide divergence, and high non-synonymous to synonymous substitutions (Ka/Ks) between or within species, suggesting rapid diversification and balancing selection in response to rapid changes in the oomycete pathogen genomes. These characteristics may serve as signatures for cloning of late blight resistance genes.

Characterization of the heterotrimeric G-protein family and its transmembrane regulator from capsicum (Capsicum annuum L.)

Throughout evolution, organisms have created numerous mechanisms to sense and respond to their environment. One such highly conserved mechanism involves regulation by heterotrimeric G-protein complex comprised of alpha (Gα), beta (Gβ) and gamma (Gγ) subunits. In plants, these proteins play important roles in signal transduction pathways related to growth and development including response to biotic and abiotic stresses and consequently affect yield. In this work, we have identified and characterized the complete heterotrimeric G-protein repertoire in the Capsicum annuum (Capsicum) genome which consists of one Gα, one Gβ and three Gγ genes. We have also identified one RGS gene in the Capsicum genome that acts as a regulator of the G-protein signaling. Biochemical activities of the proteins were confirmed by assessing the GTP-binding and GTPase activity of the recombinant Gα protein and its regulation by the GTPase acceleration activity of the RGS protein. Interaction between different subunits was established using yeast- and plant-based analyses. Gene and protein expression profiles of specific G-protein components revealed interesting spatial and temporal regulation patterns, especially during root development and during fruit development and maturation. This research thus details the characterization of the first heterotrimeric G-protein family from a domesticated, commercially important vegetable crop.

The WRKY transcription factor genes in eggplant (Solanum melongena L.) and Turkey Berry (Solanum torvum Sw.)

WRKY transcription factors, which play critical roles in stress responses, have not been characterized in eggplant or its wild relative, turkey berry. The recent availability of RNA-sequencing data provides the opportunity to examine WRKY genes from a global perspective. We identified 50 and 62 WRKY genes in eggplant (SmelWRKYs) and turkey berry (StorWRKYs), respectively, all of which could be classified into three groups (I–III) based on the WRKY protein structure. The SmelWRKYs and StorWRKYs contain ~76% and ~95% of the number of WRKYs found in other sequenced asterid species, respectively. Positive selection analysis revealed that different selection constraints could have affected the evolution of these groups. Positively-selected sites were found in Groups IIc and III. Branch-specific selection pressure analysis indicated that most WRKY domains from SmelWRKYs and StorWRKYs are conserved and have evolved at low rates since their divergence. Comparison to homologous WRKY genes in Arabidopsis revealed several potential pathogen resistance-related SmelWRKYs and StorWRKYs, providing possible candidate genetic resources for improving stress tolerance in eggplant and probably other Solanaceae plants. To our knowledge, this is the first report of a genome-wide analyses of the SmelWRKYs and StorWRKYs.

Functional studies of transcription factors involved in plant defenses in the genomics era

Plant transcription factors (TFs) play roles in diverse biological processes including defense responses to pathogens. Here, we provide an overview of recent studies of plant TFs with regard to defense responses. TFs play roles in plant innate immunity by regulating genes related to pathogen-associated molecular pattern-triggered immunity, effector-triggered immunity, hormone signaling pathways and phytoalexin synthesis. Currently, genome-wide phylogenetic and transcriptomic analyses are as important as functional analyses in the study of plant TFs. The integration of genomics information with the knowledge obtained from functional studies provides new insights into the regulation of plant defense mechanisms as well as engineering crops with improved resistance to invading pathogens.

Advances in the understanding of cuticular waxes in Arabidopsis thaliana and crop species

The aerial parts of plants are covered with a cuticle, a hydrophobic layer consisting of cutin polyester and cuticular waxes that protects them from various environmental stresses. Cuticular waxes mainly comprise very long chain fatty acids and their derivatives such as aldehydes, alkanes, secondary alcohols, ketones, primary alcohols, and wax esters that are also important raw materials for the production of lubricants, adhesives, cosmetics, and biofuels. The major function of cuticular waxes is to control non-stomatal water loss and gas exchange. In recent years, the in planta roles of many genes involved in cuticular wax biosynthesis have been characterized not only from model organisms like Arabidopsis thaliana and saltwater cress (Eutrema salsugineum), but also crop plants including maize, rice, wheat, tomato, petunia, Medicago sativa, Medicago truncatula, rapeseed, and Camelina sativa through genetic, biochemical, molecular, genomic, and cell biological approaches. In this review, we discuss recent advances in the understanding of the biological functions of genes involved in cuticular wax biosynthesis, transport, and regulation of wax deposition from Arabidopsis and crop species, provide information on cuticular wax amounts and composition in various organs of nine representative plant species, and suggest the important issues that need to be investigated in this field of study.

Development of molecular markers tightly linked to Pvr4 gene in pepper using next-generation sequencing

It is imperative to identify highly polymorphic and tightly linked markers of a known trait for molecular marker-assisted selection. Potyvirus resistance 4 (Pvr4) locus in pepper confers resistance to three pathotypes of potato virus Y and to pepper mottle virus. We describe the use of next-generation sequencing technology to generate molecular markers tightly linked to Pvr4. Initially, comparative genomics was carried out, and a syntenic region of tomato on chromosome ten was used to generate PCR-based markers and map Pvr4. Subsequently, the genomic sequence of pepper was used, and more than 5000 single-nucleotide variants (SNVs) were identified within the interval. In addition, we identified nucleotide binding site-leucine-rich repeat-type disease resistance genes within the interval. Several of these SNVs were converted to molecular markers desirable for large-scale molecular breeding programmes.

RNA-dependent RNA polymerase (NIb) of the potyviruses is an avirulence factor for the broad-spectrum resistance gene Pvr4 in Capsicum annuum cv. CM334

Potyviruses are one of the most destructive viral pathogens of Solanaceae plants. In Capsicum annuum landrace CM334, a broad-spectrum gene, Pvr4 is known to be involved in resistance against multiple potyviruses, including Pepper mottle virus (PepMoV), Pepper severe mosaic virus (PepSMV), and Potato virus Y (PVY). However, a potyvirus avirulence factor against Pvr4 has not been identified. To identify the avirulence factor corresponding to Pvr4 in potyviruses, we performed Agrobacterium-mediated transient expressions of potyvirus protein coding regions in potyvirus-resistant (Pvr4) and -susceptible (pvr4) pepper plants. Hypersensitive response (HR) was observed only when a RNA-dependent RNA polymerase (NIb) of PepMoV, PepSMV, or PVY was expressed in Pvr4-bearing pepper leaves in a genotype-specific manner. In contrast, HR was not observed when the NIb of Tobacco etch virus (TEV), a virulent potyvirus, was expressed in Pvr4-bearing pepper leaves. Our results clearly demonstrate that NIbs of PepMoV, PepSMV, and PVY serve as avirulence factors for Pvr4 in pepper plants.

NexGenEx-Tom: a gene expression platform to investigate the functionalities of the tomato genome

BACKGROUND

Next Generation Sequencing technologies (NGS) unexpectedly pushed forward the capability of solving genome organization and of widely depicting gene expression. However, although the flourishing of tools to process the NGS data, versatile and user-friendly computational environments for integrative and comparative analyses of the results from the increasing amount of collections are still required.

DESCRIPTION

Here we present the architecture and the facilities of NexGenEx-, a web based platform that offers processed NGS transcriptome collections and enables immediate analyses of the results. The platform allows gene expression investigations, profiling and comparisons, and exploits different resources.

CONCLUSION

In the current version, NexGenEx-Tom includes processed and normalized NGS expression data from three collections covering several tissue/stages from different genotypes. Beyond providing a user-friendly interface, the platform was designed with the aim to easily be expanded to include other NGS based transcriptome collections. It can also integrate different genome releases, possibly from different cultivars or genotypes, but even from different species. The platform is proposed as an example effort in tomato, and is described as a profitable approach for the exploitation of these challenging and precious datasets.

Current status and prospects for the study of Nicotiana genomics, genetics, and nicotine biosynthesis genes

Nicotiana, a member of the Solanaceae family, is one of the most important research model plants, and of high agricultural and economic value worldwide. To better understand the substantial and rapid research progress with Nicotiana in recent years, its genomics, genetics, and nicotine gene studies are summarized, with useful web links. Several important genetic maps, including a high-density map of N. tabacum consisting of ~2,000 markers published in 2012, provide tools for genetics research. Four whole genome sequences are from allotetraploid species, including N. benthamiana in 2012, and three N. tabacum cultivars (TN90, K326, and BX) in 2014. Three whole genome sequences are from diploids, including progenitors N. sylvestris and N. tomentosiformis in 2013 and N. otophora in 2014. These and additional studies provide numerous insights into genome evolution after polyploidization, including changes in gene composition and transcriptome expression in N. tabacum. The major genes involved in the nicotine biosynthetic pathway have been identified and the genetic basis of the differences in nicotine levels among Nicotiana species has been revealed. In addition, other progress on chloroplast, mitochondrial, and NCBI-registered projects on Nicotiana are discussed. The challenges and prospects for genomic, genetic and application research are addressed. Hence, this review provides important resources and guidance for current and future research and application in Nicotiana.

An InDel-based linkage map of hot pepper (Capsicum annuum)

Two independent pepper (Capsicum annuum) genomes were published recently, opening a new era of molecular genetics research on pepper. However, pepper molecular marker technologies are still mainly focusing on the simple sequence repeats derived from public database or genomic library. The development and application of the third generation marker system such as single nucleotide polymorphisms, structure variations as well as insertion/deletion polymorphisms (InDels) is still in its infancy. In the present study, we developed InDel markers for pepper genetic mapping with the convenience of two whole-genome re-sequenced inbred lines BA3 (C. annuum) and B702 (C. annuum). A total of 154,519 and 149,755 InDel (1-5 bp) sites were identified for BA3 and B702, respectively, by the alignment of re-sequencing reads to Zunla-1 reference genome. Then, 14,498 InDel sites (only 4 and 5 bp) that are different between BA3 and B702 were predicted. Finally, within a random set of 1,000 primer pairs, 251 InDel markers were validated and mapped onto a linkage map using F₂ population derived from the intraspecific cross BA3 × B702. The first InDel-based map, named as BB-InDel map, consisted of 12 linkage groups, covered a genetic distance of 1,178.01 cM and the average distance between bin markers was 5.01 cM. Compared to the Zunla-1 reference physical map, high consistency was observed on all 12 chromosomes, and the total length of scaffold anchored and physical distance covered by this map was 299.66 and 2,558.68 Mb, respectively, which accounted for 8.95 and 76.38 % of the Zunla-1 reference genome (3.35 Gb), respectively. Furthermore, 37 scaffolds (total length of 36.21 Mb) from the pseudo-chromosome (P0) of the current genome assembly were newly assigned to the corresponding chromosomes by 40 InDel markers. Thus, this map provided good genome coverage and would be useful for basic and applied research in pepper.

Genome-wide analysis of the CaHsp20 gene family in pepper: comprehensive sequence and expression profile analysis under heat stress

The Hsp20 genes are present in all plant species and play important roles in alleviating heat stress and enhancing plant thermotolerance by preventing the irreversible aggregation of denaturing proteins. However, very little is known about the CaHsp20 gene family in pepper (Capsicum annuum L.), an important vegetable crop with character of temperate but thermosensitive. In this study, a total of 35 putative pepper Hsp20 genes (CaHsp20s) were identified and renamed on the basis of their molecular weight, and then their gene structure, genome location, gene duplication, phylogenetic relationship, and interaction network were also analyzed. The expression patterns of CaHsp20 genes in four different tissues (root, stem, leaf, and flower) from the thermotolerant line R9 under heat stress condition were measured using semi-quantitative RT-PCR. The transcripts of most CaHsp20 genes maintained a low level in all of the four tissues under normal temperature condition, but were highly induced by heat stress, while the expression of CaHsp16.6b, 16.7, and 23.8 were only detected in specific tissues and were not so sensitive to heat stress like other CaHsp20 genes. In addition, compared to those in thermotolerant line R9, the expression peak of most CaHsp20 genes in thermosensitive line B6 under heat stress was hysteretic, and several CaHsp20 genes (CaHsp16.4, 18.2a, 18.7, 21.2, 22.0, 25.8, and 25.9) showed higher expression levels in both line B6 and R9. These data suggest that the CaHsp20 genes may be involved in heat stress and defense responses in pepper, which provides the basis for further functional analyses of CaHsp20s in the formation of pepper acquired thermotoleance.

Next Generation Sequencing Technologies: The Doorway to the Unexplored Genomics of Non-Model Plants

Non-model plants i.e., the species which have one or all of the characters such as long life cycle, difficulty to grow in the laboratory or poor fecundity, have been schemed out of sequencing projects earlier, due to high running cost of Sanger sequencing. Consequently, the information about their genomics and key biological processes are inadequate. However, the advent of fast and cost effective next generation sequencing (NGS) platforms in the recent past has enabled the unearthing of certain characteristic gene structures unique to these species. It has also aided in gaining insight about mechanisms underlying processes of gene expression and secondary metabolism as well as facilitated development of genomic resources for diversity characterization, evolutionary analysis and marker assisted breeding even without prior availability of genomic sequence information. In this review we explore how different Next Gen Sequencing platforms, as well as recent advances in NGS based high throughput genotyping technologies are rewarding efforts on de-novo whole genome/transcriptome sequencing, development of genome wide sequence based markers resources for improvement of non-model crops that are less costly than phenotyping.

Molecular and cellular control of cell death and defense signaling in pepper

Pepper (Capsicum annuum L.) provides a good experimental system for studying the molecular and functional genomics underlying the ability of plants to defend themselves against microbial pathogens. Cell death is a genetically programmed response that requires specific host cellular factors. Hypersensitive response (HR) is defined as rapid cell death in response to a pathogen attack. Pepper plants respond to pathogen attacks by activating genetically controlled HR- or disease-associated cell death. HR cell death, specifically in incompatible interactions between pepper and Xanthomonas campestris pv. vesicatoria, is mediated by the molecular genetics and biochemical machinery that underlie pathogen-induced cell death in plants. Gene expression profiles during the HR-like cell death response, virus-induced gene silencing and transient and transgenic overexpression approaches are used to isolate and identify HR- or disease-associated cell death genes in pepper plants. Reactive oxygen species, nitric oxide, cytosolic calcium ion and defense-related hormones such as salicylic acid, jasmonic acid, ethylene and abscisic acid are involved in the execution of pathogen-induced cell death in plants. In this review, we summarize recent molecular and cellular studies of the pepper cell death-mediated defense response, highlighting the signaling events of cell death in disease-resistant pepper plants. Comprehensive knowledge and understanding of the cellular functions of pepper cell death response genes will aid the development of novel practical approaches to enhance disease resistance in pepper, thereby helping to secure the future supply of safe and nutritious pepper plants worldwide.

Genome-wide analysis of the CaHsp20 gene family in pepper: comprehensive sequence and expression profile analysis under heat stress

The Hsp20 genes are present in all plant species and play important roles in alleviating heat stress and enhancing plant thermotolerance by preventing the irreversible aggregation of denaturing proteins. However, very little is known about the CaHsp20 gene family in pepper (Capsicum annuum L.), an important vegetable crop with character of temperate but thermosensitive. In this study, a total of 35 putative pepper Hsp20 genes (CaHsp20s) were identified and renamed on the basis of their molecular weight, and then their gene structure, genome location, gene duplication, phylogenetic relationship, and interaction network were also analyzed. The expression patterns of CaHsp20 genes in four different tissues (root, stem, leaf, and flower) from the thermotolerant line R9 under heat stress condition were measured using semi-quantitative RT-PCR. The transcripts of most CaHsp20 genes maintained a low level in all of the four tissues under normal temperature condition, but were highly induced by heat stress, while the expression of CaHsp16.6b, 16.7, and 23.8 were only detected in specific tissues and were not so sensitive to heat stress like other CaHsp20 genes. In addition, compared to those in thermotolerant line R9, the expression peak of most CaHsp20 genes in thermosensitive line B6 under heat stress was hysteretic, and several CaHsp20 genes (CaHsp16.4, 18.2a, 18.7, 21.2, 22.0, 25.8, and 25.9) showed higher expression levels in both line B6 and R9. These data suggest that the CaHsp20 genes may be involved in heat stress and defense responses in pepper, which provides the basis for further functional analyses of CaHsp20s in the formation of pepper acquired thermotoleance.

Genome-wide identification and transcriptional expression analysis of mitogen-activated protein kinase and mitogen-activated protein kinase kinase genes in Capsicum annuum

The tripartite mitogen-activated protein kinase (MAPK) signaling cascades have been implicated in plant growth, development, and environment adaptation, but a comprehensive understanding of MAPK signaling at genome-wide level is limited in Capsicum annuum. Herein, genome-wide identification and transcriptional expression analysis of MAPK and MAPK kinase (MAPKK) were performed in pepper. A total of 19 pepper MAPK (CaMAPKs) genes and five MAPKK (CaMAPKKs) genes were identified. Phylogenetic analysis indicated that CaMAPKs and CaMAPKKs could be classified into four groups and each group contains similar exon-intron structures. However, significant divergences were also found. Notably, five members of the pepper MAPKK family were much less conserved than those found in Arabidopsis, and 9 Arabidopsis MAPKs did not have orthologs in pepper. Additionally, 7 MAPKs in Arabidopsis had either two or three orthologs in the pepper genome, and six pepper MAPKs and one MAPKK differing in sequence were found in three pepper varieties. Quantitative real-time RT-PCR analysis showed that the majority of MAPK and MAPKK genes were ubiquitously expressed and transcriptionally modified in pepper leaves after treatments with heat, salt, and Ralstonia solanacearum inoculation as well as exogenously applied salicylic acid, methyl jasmonate, ethephon, and abscisic acid. The MAPKK-MAPK interactome was tested by yeast two-hybrid assay, the results showed that one MAPKK might interact with multiple MAPKs, one MAPK might also interact with more than one MAPKKs, constituting MAPK signaling networks which may collaborate in transmitting upstream signals into appropriate downstream cellular responses and processes. These results will facilitate future functional characterization of MAPK cascades in pepper.

Genome-wide identification and expression analysis of calcium-dependent protein kinase and its closely related kinase genes in Capsicum annuum

As Ca2+ sensors and effectors, calcium-dependent protein kinases (CDPKs) play important roles in plant growth, development, and response to environmental cues. However, no CDPKs have been characterized in Capsicum annuum thus far. Herein, a genome wide comprehensive analysis of genes encoding CDPKs and CDPK-related protein kinases (CRKs) was performed in pepper, a total of 31 CDPK genes and five closely related kinase genes were identified, which were phylogenetically divided into four distinct subfamilies and unevenly distributed across nine chromosomes. Conserved sequence and exon-intron structures were found to be shared by pepper CDPKs within the same subfamily, and the expansion of the CDPK family in pepper was found to be due to segmental duplication events. Five CDPKs in the C. annuum variety CM334 were found to be mutated in the Chiltepin variety, and one CDPK present in CM334 was lost in Chiltepin. The majority of CDPK and CRK genes were expressed in different pepper tissues and developmental stages, and 10, 12, and 8 CDPK genes were transcriptionally modified by salt, heat, and Ralstonia solanacearum stresses, respectively. Furthermore, these genes were found to respond specifically to one stress as well as respond synergistically to two stresses or three stresses, suggesting that these CDPK genes might be involved in the specific or synergistic response of pepper to salt, heat, and R. solanacearum. Our results lay the foundation for future functional characterization of pepper CDPK and its closely related gene families.

New insights into Capsicum spp relatedness and the diversification process of Capsicum annuum in Spain

The successful exploitation of germplasm banks, harbouring plant genetic resources indispensable for plant breeding, will depend on our ability to characterize their genetic diversity. The Vegetable Germplasm Bank of Zaragoza (BGHZ) (Spain) holds an important Capsicum annuum collection, where most of the Spanish pepper variability is represented, as well as several accessions of other domesticated and non-domesticated Capsicum spp from all over the five continents. In the present work, a total of 51 C. annuum landraces (mainly from Spain) and 51 accessions from nine Capsicum species maintained at the BGHZ were evaluated using 39 microsatellite (SSR) markers spanning the whole genome. The 39 polymorphic markers allowed the detection of 381 alleles, with an average of 9.8 alleles per locus. A sizeable proportion of alleles (41.2%) were recorded as specific alleles and the majority of these were present at very low frequencies (rare alleles). Multivariate and model-based analyses partitioned the collection in seven clusters comprising the ten different Capsicum spp analysed: C. annuum, C. chinense, C. frutescens, C. pubescens, C. bacatum, C. chacoense and C. eximium. The data clearly showed the close relationships between C. chinense and C. frutescens. C. cardenasii and C. eximium were indistinguishable as a single, morphologically variable species. Moreover, C. chacoense was placed between C. baccatum and C. pubescens complexes. The C. annuum group was structured into three main clusters, mostly according to the pepper fruit shape, size and potential pungency. Results suggest that the diversification of C. annuum in Spain may occur from a rather limited gene pool, still represented by few landraces with ancestral traits. This ancient population would suffer from local selection at the distinct geographical regions of Spain, giving way to pungent and elongated fruited peppers in the South and Center, while sweet blocky and triangular types in Northern Spain.

Endogenous pararetrovirus sequences associated with 24 nt small RNAs at the centromeres of Fritillaria imperialis L. (Liliaceae), a species with a giant genome

Endogenous pararetroviral sequences are the most commonly found virus sequences integrated into angiosperm genomes. We describe an endogenous pararetrovirus (EPRV) repeat in Fritillaria imperialis, a species that is under study as a result of its exceptionally large genome (1C = 42 096 Mbp, approximately 240 times bigger than Arabidopsis thaliana). The repeat (FriEPRV) was identified from Illumina reads using the RepeatExplorer pipeline, and exists in a complex genomic organization at the centromere of most, or all, chromosomes. The repeat was reconstructed into three consensus sequences that formed three interconnected loops, one of which carries sequence motifs expected of an EPRV (including the gag and pol domains). FriEPRV shows sequence similarity to members of the Caulimoviridae pararetrovirus family, with phylogenetic analysis indicating a close relationship to Petuvirus. It is possible that no complete EPRV sequence exists, although our data suggest an abundance that exceeds the genome size of Arabidopsis. Analysis of single nucleotide polymorphisms revealed elevated levels of C→T and G→A transitions, consistent with deamination of methylated cytosine. Bisulphite sequencing revealed high levels of methylation at CG and CHG motifs (up to 100%), and 15-20% methylation, on average, at CHH motifs. FriEPRV's centromeric location may suggest targeted insertion, perhaps associated with meiotic drive. We observed an abundance of 24 nt small RNAs that specifically target FriEPRV, potentially providing a signature of RNA-dependent DNA methylation. Such signatures of epigenetic regulation suggest that the huge genome of F. imperialis has not arisen as a consequence of a catastrophic breakdown in the regulation of repeat amplification.

The Sol Genomics Network (SGN)--from genotype to phenotype to breeding

The Sol Genomics Network (SGN, http://solgenomics.net) is a web portal with genomic and phenotypic data, and analysis tools for the Solanaceae family and close relatives. SGN hosts whole genome data for an increasing number of Solanaceae family members including tomato, potato, pepper, eggplant, tobacco and Nicotiana benthamiana. The database also stores loci and phenotype data, which researchers can upload and edit with user-friendly web interfaces. Tools such as BLAST, GBrowse and JBrowse for browsing genomes, expression and map data viewers, a locus community annotation system and a QTL analysis tools are available. A new tool was recently implemented to improve Virus-Induced Gene Silencing (VIGS) constructs called the SGN VIGS tool. With the growing genomic and phenotypic data in the database, SGN is now advancing to develop new web-based breeding tools and implement the code and database structure for other species or clade-specific databases.

Insights into transcriptional regulation of β-D-N-acetylhexosaminidase, an N-glycan-processing enzyme involved in ripening-associated fruit softening

Tomato (Solanum lycopersicum) fruit ripening-specific N-glycan processing enzyme, β-D-N-acetylhexosaminidase (β-Hex), plays an important role in the ripening-associated fruit-softening process. However, the regulation of fruit ripening-specific expression of β-Hex is not well understood. We have identified and functionally characterized the fruit ripening-specific promoter of β-Hex and provided insights into its transcriptional regulation during fruit ripening. Our results demonstrate that RIPENING INHIBITOR (RIN), a global fruit ripening regulator, and ABSCISIC ACID STRESS RIPENING 1 (SlASR1), a poorly characterized ripening-related protein, are the transcriptional regulators of β-Hex. Both RIN and SlASR1 directly bound to the β-Hex promoter fragments containing CArG and C₂₋₃(C/G)A cis-acting elements, the binding sites for RIN and SlASR1, respectively. Moreover, β-Hex expression/promoter activity in tomato fruits was downregulated once expression of either RIN or SlASR1 was suppressed; indicating that RIN and SlASR1 positively regulate the transcription of β-Hex during fruit ripening. Interestingly, RIN could also bind to the SlASR1 promoter, which contains several CArG cis-acting elements, and SlASR1 expression was suppressed in rin mutant fruits, indicating that RIN also acts as a positive regulator of SlASR1 expression during fruit ripening. Taken together, these results suggest that RIN, both directly and indirectly, through SlASR1, regulates the transcription of β-Hex during fruit ripening. The fruit ripening-specific promoter of β-Hex could be a useful tool in regulating gene expression during fruit ripening.