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

Identification of Optimal Reference Genes for Normalization of qPCR Analysis during Pepper Fruit Development

Due to its high sensitivity and reproducibility, quantitative real-time PCR (qPCR) is practiced as a useful research tool for targeted gene expression analysis. For qPCR operations, the normalization with suitable reference genes (RGs) is a crucial step that eventually determines the reliability of the obtained results. Although pepper is considered an ideal model plant for the study of non-climacteric fruit development, at present no specific RG have been developed or validated for the qPCR analyses of pepper fruit. Therefore, this study aimed to identify stably expressed genes for their potential use as RGs in pepper fruit studies. Initially, a total of 35 putative RGs were selected by mining the pepper transcriptome data sets derived from the PGP (Pepper Genome Platform) and PGD (Pepper Genome Database). Their expression stabilities were further measured in a set of pepper (Capsicum annuum L. var. 007e) fruit samples, which represented four different fruit developmental stages (IM: Immature; MG: Mature green; B: Break; MR: Mature red) using the qPCR analysis. Then, based on the qPCR results, three different statistical algorithms, namely geNorm, Normfinder, and boxplot, were chosen to evaluate the expression stabilities of these putative RGs. It should be noted that nine genes were proven to be qualified as RGs during pepper fruit development, namely CaREV05 (CA00g79660); CaREV08 (CA06g02180); CaREV09 (CA06g05650); CaREV16 (Capana12g002666); CaREV21 (Capana10g001439); CaREV23 (Capana05g000680); CaREV26 (Capana01g002973); CaREV27 (Capana11g000123); CaREV31 (Capana04g002411); and CaREV33 (Capana08g001826). Further analysis based on geNorm suggested that the application of the two most stably expressed genes (CaREV05 and CaREV08) would provide optimal transcript normalization in the qPCR experiments. Therefore, a new and comprehensive strategy for the identification of optimal RGs was developed. This strategy allowed for the effective normalization of the qPCR analysis of the pepper fruit development at the whole pepper genome level. This study not only explored the optimal RGs specific for studying pepper fruit development, but also introduced a referable strategy of RG mining which could potentially be implicated in other plant species.

Characterization of BiP Genes from Pepper (Capsicum annuum L.) and the Role of CaBiP1 in Response to Endoplasmic Reticulum and Multiple Abiotic Stresses

Adverse environmental conditions have a detrimental impact on crop growth and development, and cause protein denaturation or misfolding. The binding protein (BiP) plays an important protective role by alleviating endoplasmic reticulum (ER) stress induced by misfolded proteins. In this study, we characterized three BiP genes (CaBiP1, CaBiP2, and CaBiP3) in pepper, an economically important vegetable and spice species. The role of CaBiP1 in plant tolerance to ER stress and adverse environmental conditions (including heat, salinity, osmotic and drought stress) were investigated. All the expected functional and signaling domains were detected in three BiP proteins, but the motifs and exon-intron distribution differed slightly in CaBiP3. CaBiP1 and CaBiP2 were constitutively expressed in all the tested tissues under both normal and stressed conditions, whereas CaBiP3 was mainly expressed following stress. Silencing of CaBiP1 reduced pepper tolerance to ER stress and various environment stresses, and was accompanied by increased H₂O₂ accumulation, MDA content, relative electric leakage (REL), water loss rate, and a reduction in soluble protein content and relative water content (RWC) in the leaves. Conversely, overexpression of CaBiP1 in Arabidopsis enhanced tolerance to ER stress and multiple environment stresses, as demonstrated by an increase in germination rate, root length, survival rate, RWC, the unfolded protein response (UPR) pathway, and a decrease in water loss rate. Our results suggest that CaBiP1 may contribute to plant tolerance to abiotic stresses by reducing ROS accumulation, increasing the water-retention ability, and stimulating UPR pathways and expression of stress-related genes.

Divergent evolution of multiple virus-resistance genes from a progenitor in Capsicum spp

Plants have evolved hundreds of nucleotide-binding and leucine-rich domain proteins (NLRs) as potential intracellular immune receptors, but the evolutionary mechanism leading to the ability to recognize specific pathogen effectors is elusive. Here, we cloned Pvr4 (a Potyvirus resistance gene in Capsicum annuum) and Tsw (a Tomato spotted wilt virus resistance gene in Capsicum chinense) via a genome-based approach using independent segregating populations. The genes both encode typical NLRs and are located at the same locus on pepper chromosome 10. Despite the fact that these two genes recognize completely different viral effectors, the genomic structures and coding sequences of the two genes are strikingly similar. Phylogenetic studies revealed that these two immune receptors diverged from a progenitor gene of a common ancestor. Our results suggest that sequence variations caused by gene duplication and neofunctionalization may underlie the evolution of the ability to specifically recognize different effectors. These findings thereby provide insight into the divergent evolution of plant immune receptors.

Molecular Mapping of PMR1, a Novel Locus Conferring Resistance to Powdery Mildew in Pepper (Capsicum annuum)

Powdery mildew, caused by Leveillula taurica, is a major fungal disease affecting greenhouse-grown pepper (Capsicum annuum). Powdery mildew resistance has a complex mode of inheritance. In the present study, we investigated a novel powdery mildew resistance locus, PMR1, using two mapping populations: 102 'VK515' F_2:3 families (derived from a cross between resistant parental line 'VK515R' and susceptible parental line 'VK515S') and 80 'PM Singang' F₂ plants (derived from the F₁ 'PM Singang' commercial hybrid). Genetic analysis of the F_2:3 'VK515' and F₂ 'PM Singang' populations revealed a single dominant locus for inheritance of the powdery mildew resistance trait. Genetic mapping showed that the PMR1 locus is located on syntenic regions of pepper chromosome 4 in a 4-Mb region between markers CZ2_11628 and HRM4.1.6 in 'VK515R'. Six molecular markers including one SCAR marker and five SNP markers were localized to a region 0 cM from the PMR1 locus. Two putative nucleotide-binding site leucine-rich repeat (NBS-LRR)-type disease resistance genes were identified in this PMR1 region. Genotyping-by-sequencing (GBS) and genetic mapping analysis revealed suppressed recombination in the PMR1 region, perhaps due to alien introgression. In addition, a comparison of species-specific InDel markers as well as GBS-derived SNP markers indicated that C. baccatum represents a possible source of such alien introgression of powdery mildew resistance into 'VK515R'. The molecular markers developed in this study will be especially helpful for marker-assisted selection in pepper breeding programs for powdery mildew resistance.

Genome sequence and genetic diversity of European ash trees

Ash trees (genus Fraxinus, family Oleaceae) are widespread throughout the Northern Hemisphere, but are being devastated in Europe by the fungus Hymenoscyphus fraxineus, causing ash dieback, and in North America by the herbivorous beetle Agrilus planipennis. Here we sequence the genome of a low-heterozygosity Fraxinus excelsior tree from Gloucestershire, UK, annotating 38,852 protein-coding genes of which 25% appear ash specific when compared with the genomes of ten other plant species. Analyses of paralogous genes suggest a whole-genome duplication shared with olive (Olea europaea, Oleaceae). We also re-sequence 37 F. excelsior trees from Europe, finding evidence for apparent long-term decline in effective population size. Using our reference sequence, we re-analyse association transcriptomic data, yielding improved markers for reduced susceptibility to ash dieback. Surveys of these markers in British populations suggest that reduced susceptibility to ash dieback may be more widespread in Great Britain than in Denmark. We also present evidence that susceptibility of trees to H. fraxineus is associated with their iridoid glycoside levels. This rapid, integrated, multidisciplinary research response to an emerging health threat in a non-model organism opens the way for mitigation of the epidemic.

Putative WRKYs associated with regulation of fruit ripening revealed by detailed expression analysis of the WRKY gene family in pepper

WRKY transcription factors play important roles in plant development and stress responses. Here, global expression patterns of pepper CaWRKYs in various tissues as well as response to environmental stresses and plant hormones were systematically analyzed, with an emphasis on fruit ripening. The results showed that most CaWRKYs were expressed in at least two of the tissues tested. Group I, a subfamily of the entire CaWRKY gene family, had a higher expression level in vegetative tissues, whereas groups IIa and III showed relatively lower expression levels. Comparative analysis showed that the constitutively highly expressed WRKY genes were conserved in tomato and pepper, suggesting potential functional similarities. Among the identified 61 CaWRKYs, almost 60% were expressed during pepper fruit maturation, and the group I genes were in higher proportion during the ripening process, indicating an as-yet unknown function of group I in the fruit maturation process. Further analysis suggested that many CaWRKYs expressed during fruit ripening were also regulated by abiotic stresses or plant hormones, indicating that these CaWRKYs play roles in the stress-related signaling pathways during fruit ripening. This study provides new insights to the current research on CaWRKY and contributes to our knowledge about the global regulatory network in pepper fruit ripening.

Dissimilar evolutionary histories of two resistance gene families in the genus Solanum

Genomic analyses have shown that most genes in eukaryotic lineages belong to families. Gene families vary in terms of number of members, nucleotide similarity, gene integrity, expression, and function. Often, the members of gene families are arranged in clusters, which contribute to maintaining similarity among gene copies and also to generate duplicates through replication errors. Gene families offer us an opportunity to examine the forces involved in the evolution of the genomes and to study recombination events and genomic rearrangements. In this work, we focused on the evolution of two plant resistance gene families, Sw5 and Mi-1, and analyzed the completely sequenced nuclear genomes of potato and tomato. We first noticed that the potato genome carries larger resistance gene families than tomato, but all gene copies are pseudogenes. Second, phylogenetic analyses indicated that Sw5 and Mi-1 gene families had dissimilar evolutionary histories. In contrast to Sw5, Mi-1 homologues suffered repeated gene conversion events among the gene copies, particularly in the tomato genome.

Comparative Analysis of Fruit Metabolites and Pungency Candidate Genes Expression between Bhut Jolokia and Other Capsicum Species

Bhut jolokia, commonly known as Ghost chili, a native Capsicum species found in North East India was recorded as the naturally occurring hottest chili in the world by the Guinness Book of World Records in 2006. Although few studies have reported variation in pungency content of this particular species, no study till date has reported detailed expression analysis of candidate genes involved in capsaicinoids (pungency) biosynthesis pathway and other fruit metabolites. Therefore, the present study was designed to evaluate the diversity of fruit morphology, fruiting habit, capsaicinoids and other metabolite contents in 136 different genotypes mainly collected from North East India. Significant intra and inter-specific variations for fruit morphological traits, fruiting habits and 65 fruit metabolites were observed in the collected Capsicum germplasm belonging to three Capsicum species i.e., Capsicum chinense (Bhut jolokia, 63 accessions), C. frutescens (17 accessions) and C. annuum (56 accessions). The pungency level, measured in Scoville Heat Unit (SHU) and antioxidant activity measured by 2, 2-diphenyl-1-picrylhydrazyl (DPPH) free radical scavenging assay showed maximum levels in C. chinense accessions followed by C. frutescens accessions, while C. annuum accessions showed the lowest value for both the traits. The number of different fruit metabolites detected did not vary significantly among the different species but the metabolite such as benzoic acid hydroxyl esters identified in large percentage in majority of C. annuum genotypes was totally absent in the C. chinense genotypes and sparingly present in few genotypes of C. frutescens. Significant correlations were observed between fruit metabolites capsaicin, dihydrocapsaicin, hexadecanoic acid, cyclopentane, α-tocopherol and antioxidant activity. Furthermore, comparative expression analysis (through qRT-PCR) of candidate genes involved in capsaicinoid biosynthesis pathway revealed many fold higher expression of majority of the genes in C. chinense compared to C. frutescens and C. annuum suggesting that the possible reason for extremely high pungency might be due to the higher level of candidate gene(s) expression although nucleotide variation in pungency related genes may also be involved in imparting variations in level of pungency.

The historical role of species from the Solanaceae plant family in genetic research

This article evaluates the main contributions of tomato, tobacco, petunia, potato, pepper and eggplant to classical and molecular plant genetics and genomics since the beginning of the twentieth century. Species from the Solanaceae family form integral parts of human civilizations as food sources and drugs since thousands of years, and, more recently, as ornamentals. Some Solanaceous species were subjects of classical and molecular genetic research over the last 100 years. The tomato was one of the principal models in twentieth century classical genetics and a pacemaker of genome analysis in plants including molecular linkage maps, positional cloning of disease resistance genes and quantitative trait loci (QTL). Besides that, tomato is the model for the genetics of fruit development and composition. Tobacco was the major model used to establish the principals and methods of plant somatic cell genetics including in vitro propagation of cells and tissues, totipotency of somatic cells, doubled haploid production and genetic transformation. Petunia was a model for elucidating the biochemical and genetic basis of flower color and development. The cultivated potato is the economically most important Solanaceous plant and ranks third after wheat and rice as one of the world's great food crops. Potato is the model for studying the genetic basis of tuber development. Molecular genetics and genomics of potato, in particular association genetics, made valuable contributions to the genetic dissection of complex agronomic traits and the development of diagnostic markers for breeding applications. Pepper and eggplant are horticultural crops of worldwide relevance. Genetic and genomic research in pepper and eggplant mostly followed the tomato model. Comparative genome analysis of tomato, potato, pepper and eggplant contributed to the understanding of plant genome evolution.

Reverse Genetics and High Throughput Sequencing Methodologies for Plant Functional Genomics

In the post-genomic era, increasingly sophisticated genetic tools are being developed with the long-term goal of understanding how the coordinated activity of genes gives rise to a complex organism. With the advent of the next generation sequencing associated with effective computational approaches, wide variety of plant species have been fully sequenced giving a wealth of data sequence information on structure and organization of plant genomes. Since thousands of gene sequences are already known, recently developed functional genomics approaches provide powerful tools to analyze plant gene functions through various gene manipulation technologies. Integration of different omics platforms along with gene annotation and computational analysis may elucidate a complete view in a system biology level. Extensive investigations on reverse genetics methodologies were deployed for assigning biological function to a specific gene or gene product. We provide here an updated overview of these high throughout strategies highlighting recent advances in the knowledge of functional genomics in plants.

Genome-wide Diversity and Association Mapping for Capsaicinoids and Fruit Weight in Capsicum annuum L

Accumulated capsaicinoid content and increased fruit size are traits resulting from Capsicum annuum domestication. In this study, we used a diverse collection of C. annuum to generate 66,960 SNPs using genotyping by sequencing. The study identified 1189 haplotypes containing 3413 SNPs. Length of individual linkage disequilibrium (LD) blocks varied along chromosomes, with regions of high and low LD interspersed with an average LD of 139 kb. Principal component analysis (PCA), Bayesian model based population structure analysis and an Euclidean tree built based on identity by state (IBS) indices revealed that the clustering pattern of diverse accessions are in agreement with capsaicin content (CA) and fruit weight (FW) classifications indicating the importance of these traits in shaping modern pepper genome. PCA and IBS were used in a mixed linear model of capsaicin and dihydrocapsaicin content and fruit weight to reduce spurious associations because of confounding effects of subpopulations in genome-wide association study (GWAS). Our GWAS results showed SNPs in Ankyrin-like protein, IKI3 family protein, ABC transporter G family and pentatricopeptide repeat protein are the major markers for capsaicinoids and of 16 SNPs strongly associated with FW in both years of the study, 7 are located in known fruit weight controlling genes.

Genome-wide SNP discovery and population structure analysis in pepper (Capsicum annuum) using genotyping by sequencing

Background

Knowledge on population structure and genetic diversity in vegetable crops is essential for association mapping studies and genomic selection. Genotyping by sequencing (GBS) represents an innovative method for large scale SNP detection and genotyping of genetic resources. Herein we used the GBS approach for the genome-wide identification of SNPs in a collection of Capsicum spp. accessions and for the assessment of the level of genetic diversity in a subset of 222 cultivated pepper (Capsicum annum) genotypes.

Results

GBS analysis generated a total of 7,568,894 master tags, of which 43.4% uniquely aligned to the reference genome CM334. A total of 108,591 SNP markers were identified, of which 105,184 were in C. annuum accessions. In order to explore the genetic diversity of C. annuum and to select a minimal core set representing most of the total genetic variation with minimum redundancy, a subset of 222 C. annuum accessions were analysed using 32,950 high quality SNPs. Based on Bayesian and Hierarchical clustering it was possible to divide the collection into three clusters. Cluster I had the majority of varieties and landraces mainly from Southern and Northern Italy, and from Eastern Europe, whereas clusters II and III comprised accessions of different geographical origins. Considering the genome-wide genetic variation among the accessions included in cluster I, a second round of Bayesian (K = 3) and Hierarchical (K = 2) clustering was performed. These analysis showed that genotypes were grouped not only based on geographical origin, but also on fruit-related features.

Conclusions

GBS data has proven useful to assess the genetic diversity in a collection of C. annuum accessions. The high number of SNP markers, uniformly distributed on the 12 chromosomes, allowed the accessions to be distinguished according to geographical origin and fruit-related features. SNP markers and information on population structure developed in this study will undoubtedly support genome-wide association mapping studies and marker-assisted selection programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3297-7) contains supplementary material, which is available to authorized users.

Genome sequence and analysis of the Japanese morning glory Ipomoea nil

Ipomoea is the largest genus in the family Convolvulaceae. Ipomoea nil (Japanese morning glory) has been utilized as a model plant to study the genetic basis of floricultural traits, with over 1,500 mutant lines. In the present study, we have utilized second- and third-generation-sequencing platforms, and have reported a draft genome of I. nil with a scaffold N50 of 2.88 Mb (contig N50 of 1.87 Mb), covering 98% of the 750 Mb genome. Scaffolds covering 91.42% of the assembly are anchored to 15 pseudo-chromosomes. The draft genome has enabled the identification and cataloguing of the Tpn1 family transposons, known as the major mutagen of I. nil, and analysing the dwarf gene, CONTRACTED, located on the genetic map published in 1956. Comparative genomics has suggested that a whole genome duplication in Convolvulaceae, distinct from the recent Solanaceae event, has occurred after the divergence of the two sister families.

The evolution of inflorescence diversity in the nightshades and heterochrony during meristem maturation

One of the most remarkable manifestations of plant evolution is the diversity for floral branching systems. These “inflorescences” arise from stem cell populations in shoot meristems that mature gradually to reproductive states in response to environmental and endogenous signals. The morphology of the shoot meristem maturation process is conserved across distantly related plants, raising the question of how diverse inflorescence architectures arise from seemingly common maturation programs. In tomato and related nightshades (Solanaceae), inflorescences range from solitary flowers to highly branched structures bearing hundreds of flowers. Since reproductive barriers between even closely related Solanaceae have precluded a genetic dissection, we captured and compared meristem maturation transcriptomes from five domesticated and wild species reflecting the evolutionary continuum of inflorescence complexity. We find these divergent species share hundreds of dynamically expressed genes, enriched for transcription factors. Meristem stages are defined by distinct molecular states and point to modified maturation schedules underlying architectural variation. These modified schedules are marked by a peak of transcriptome expression divergence during the reproductive transition, driven by heterochronic shifts of dynamic genes, including transcriptional regulators with known roles in flowering. Thus, evolutionary diversity in Solanaceae inflorescence complexity is determined by subtle modifications of transcriptional programs during a critical transitional window of meristem maturation, which we propose underlies similar cases of plant architectural variation. More broadly, our findings parallel the recently described transcriptome “inverse hourglass” model for animal embryogenesis, suggesting both plant and animal morphological variation is guided by a mid-development period of transcriptome divergence.

Comparative analyses of ubiquitin-like ATG8 and cysteine protease ATG4 autophagy genes in the plant lineage and cross-kingdom processing of ATG8 by ATG4

Autophagy is important for degradation and recycling of intracellular components. In a diversity of genera and species, orthologs and paralogs of the yeast Atg4 and Atg8 proteins are crucial in the biogenesis of double-membrane autophagosomes that carry the cellular cargoes to vacuoles and lysosomes. Although many plant genome sequences are available, the ATG4 and ATG8 sequence analysis is limited to some model plants. We identified 28 ATG4 and 116 ATG8 genes from the available 18 different plant genome sequences. Gene structures and protein domain sequences of ATG4 and ATG8 are conserved in plant lineages. Phylogenetic analyses classified ATG8s into 3 subgroups suggesting divergence from the common ancestor. The ATG8 expansion in plants might be attributed to whole genome duplication, segmental and dispersed duplication, and purifying selection. Our results revealed that the yeast Atg4 processes Arabidopsis ATG8 but not human LC3A (HsLC3A). In contrast, HsATG4B can process yeast and plant ATG8s in vitro but yeast and plant ATG4s cannot process HsLC3A. Interestingly, in Nicotiana benthamiana plants the yeast Atg8 is processed compared to HsLC3A. However, HsLC3A is processed when coexpressed with HsATG4B in plants. Molecular modeling indicates that lack of processing of HsLC3A by plant and yeast ATG4 is not due to lack of interaction with HsLC3A. Our in-depth analyses of ATG4 and ATG8 in the plant lineage combined with results of cross-kingdom ATG8 processing by ATG4 further support the evolutionarily conserved maturation of ATG8. Broad ATG8 processing by HsATG4B and lack of processing of HsLC3A by yeast and plant ATG4s suggest that the cross-kingdom ATG8 processing is determined by ATG8 sequence rather than ATG4.

Genome-wide analysis of the Hsp70 family genes in pepper (Capsicum annuum L.) and functional identification of CaHsp70-2 involvement in heat stress

Hsp70s function as molecular chaperones and are encoded by a multi-gene family whose members play a crucial role in plant response to stress conditions, and in plant growth and development. Pepper (Capsicum annuum L.) is an important vegetable crop whose genome has been sequenced. Nonetheless, no overall analysis of the Hsp70 gene family is reported in this crop plant to date. To assess the functionality of Capsicum annuum Hsp70 (CaHsp70) genes, pepper genome database was analyzed in this research. A total of 21 CaHsp70 genes were identified and their characteristics were also described. The promoter and transcript expression analysis revealed that CaHsp70s were involved in pepper growth and development, and heat stress response. Ectopic expression of a cytosolic gene, CaHsp70-2, regulated expression of stress-related genes and conferred increased thermotolerance in transgenic Arabidopsis. Taken together, our results provide the basis for further studied to dissect CaHsp70s' function in response to heat stress as well as other environmental stresses.

Resequencing of Capsicum annuum parental lines (YCM334 and Taean) for the genetic analysis of bacterial wilt resistance

BACKGROUND

Bacterial wilt (BW) is a widespread plant disease that affects a broad range of dicot and monocot hosts and is particularly harmful for solanaceous plants, such as pepper, tomato, and eggplant. The pathogen responsible for BW is the soil-borne bacterium, Ralstonia solanacearum, which can adapt to diverse temperature conditions and is found in climates ranging from tropical to temperate. Resistance to BW has been detected in some pepper plant lines; however, the genomic loci and alleles that mediate this are poorly studied in this species.

RESULTS

We resequenced the pepper cultivars YCM344 and Taean, which are parental recombinant inbred lines (RIL) that display differential resistance phenotypes against BW, with YCM344 being highly resistant to infection with this pathogen. We identified novel single nucleotide polymorphisms (SNPs) and insertions/deletions (Indels) that are only present in both parental lines, as compared to the reference genome and further determined variations that distinguish these two cultivars from one another. We then identified potentially informative SNPs that were found in genes related to those that have been previously associated with disease resistance, such as the R genes and stress response genes. Moreover, via comparative analysis, we identified SNPs located in genomic regions that have homology to known resistance genes in the tomato genomes.

CONCLUSIONS

From our SNP profiling in both parental lines, we could identify SNPs that are potentially responsible for BW resistance, and practically, these may be used as markers for assisted breeding schemes using these populations. We predict that our analyses will be valuable for both better understanding the YCM334/Taean-derived populations, as well as for enhancing our knowledge of critical SNPs present in the pepper genome.

Discovery of putative capsaicin biosynthetic genes by RNA-Seq and digital gene expression analysis of pepper

The Indian pepper ‘Guijiangwang’ (Capsicum frutescens L.), one of the world’s hottest chili peppers, is rich in capsaicinoids. The accumulation of the alkaloid capsaicin and its analogs in the epidermal cells of the placenta contribute to the pungency of Capsicum fruits. To identify putative genes involved in capsaicin biosynthesis, RNA-Seq was used to analyze the pepper’s expression profiles over five developmental stages. Five cDNA libraries were constructed from the total RNA of placental tissue and sequenced using an Illumina HiSeq 2000. More than 19 million clean reads were obtained from each library, and greater than 50% of the reads were assignable to reference genes. Digital gene expression (DGE) profile analysis using Solexa sequencing was performed at five fruit developmental stages and resulted in the identification of 135 genes of known function; their expression patterns were compared to the capsaicin accumulation pattern. Ten genes of known function were identified as most likely to be involved in regulating capsaicin synthesis. Additionally, 20 new candidate genes were identified related to capsaicin synthesis. We use a combination of RNA-Seq and DGE analyses to contribute to the understanding of the biosynthetic regulatory mechanism(s) of secondary metabolites in a nonmodel plant and to identify candidate enzyme-encoding genes.

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 Identification of BAHD Acyltransferases and In vivo Characterization of HQT-like Enzymes Involved in Caffeoylquinic Acid Synthesis in Globe Artichoke

Globe artichoke (Cynara cardunculus L. var. scolymus) is a rich source of compounds promoting human health (phytonutrients), among them caffeoylquinic acids (CQAs), mainly represented by chlorogenic acid (CGA), and dicaffeoylquinic acids (diCQAs). The enzymes involved in their biosynthesis belong to the large family of BAHD acyltransferases. Following a survey of the globe artichoke genome, we identified 69 BAHD proteins carrying the catalytic site (HXXXD). Their phylogenetic analysis together with another 43 proteins, from 21 species, representative of the BAHD family, highlighted their grouping in seven major clades. Nine globe artichoke acyltransferases clustered in a sub-group of Clade V, with 3 belonging to hydroxycinnamoyl-CoA:quinate hydroxycinnamoyl transferase (HQT) and 2 to hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase (HCT) like proteins. We focused our attention on the former, HQT1, HQT2, and HQT3, as they are known to play a key role in CGA biosynthesis. The expression of genes coding for the three HQTs and correlation of expression with the CQA content is reported for different globe artichoke tissues. For the first time in the globe artichoke, we developed and applied the virus-induced gene silencing approach with the goal of assessing in vivo the effect of HQT1 silencing, which resulted in a marked reduction of both CGA and diCQAs. On the other hand, when the role of the three HQTs was assessed in leaves of Nicotiana benthamiana through their transient overexpression, significant increases in mono- and diCQAs content were observed. Using transient GFP fusion proteins expressed in N. benthamiana leaves we also established the sub-cellular localization of these three enzymes.

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.

Development of a SNP array and its application to genetic mapping and diversity assessment in pepper (Capsicum spp.)

The development and application of single nucleotide polymorphisms (SNPs) is in its infancy for pepper. Here, a set of 15,000 SNPs were chosen from the resequencing data to develop an array for pepper with 12,720 loci being ultimately synthesized. Of these, 8,199 (~64.46%) SNPs were found to be scorable and covered ~81.18% of the whole genome. With this array, a high-density interspecific genetic map with 5,569 SNPs was constructed using 297 F2 individuals, and genetic diversity of a panel of 399 pepper elite/landrace lines was successfully characterized. Based on the genetic map, one major QTL, named Up12.1, was detected for the fruit orientation trait. A total of 65 protein-coding genes were predicted within this QTL region based on the current annotation of the Zunla-1 genome. In summary, the thousands of well-validated SNP markers, high-density genetic map and genetic diversity information will be useful for molecular genetics and innovative breeding in pepper. Furthermore, the mapping results lay foundation for isolating the genes underlying variation in fruit orientation of Capsicum.

Genome-Wide Identification, Expression Diversication of Dehydrin Gene Family and Characterization of CaDHN3 in Pepper (Capsicum annuum L.)

Dehydrins (DHNs) play a crucial role in enhancing abiotic stress tolerance in plants. Although DHNs have been identified and characterized in many plants, there is little known about Capsicum annuum L., one of the economically important vegetable crops. In this study, seven CaDHNs in the pepper genome were identified, which could be divided into two classes: YnSKn- and SKn-type, based on their highly conserved domains. Quantitative real-time PCR (qRT-PCR) results showed that the seven DHN genes were expressed in all tissues and might be involved in the growth and development of pepper. The gene expression profiles analysis suggested that most of the CaDHN genes were induced by various stresses (low temperature, salt and mannitol) and signaling molecules (ABA, SA and MeJA). Furthermore, the CaDHN3 (YSK2)-silenced pepper plants showed obvious lower resistance to abiotic stresses (cold, salt and mannitol) than the control plants (TRV2:00). So the CaDHN3 might act as a positive role in resisting abiotic stresses. This study lays the foundation for further studies into the regulation of their expression under various conditions.

Genome-Wide Comparative Analyses Reveal the Dynamic Evolution of Nucleotide-Binding Leucine-Rich Repeat Gene Family among Solanaceae Plants

Plants have evolved an elaborate innate immune system against invading pathogens. Within this system, intracellular nucleotide-binding leucine-rich repeat (NLR) immune receptors are known play critical roles in effector-triggered immunity (ETI) plant defense. We performed genome-wide identification and classification of NLR-coding sequences from the genomes of pepper, tomato, and potato using fixed criteria. We then compared genomic duplication and evolution features. We identified intact 267, 443, and 755 NLR-encoding genes in tomato, potato, and pepper genomes, respectively. Phylogenetic analysis and classification of Solanaceae NLRs revealed that the majority of NLR super family members fell into 14 subgroups, including a TIR-NLR (TNL) subgroup and 13 non-TNL subgroups. Specific subgroups have expanded in each genome, with the expansion in pepper showing subgroup-specific physical clusters. Comparative analysis of duplications showed distinct duplication patterns within pepper and among Solanaceae plants suggesting subgroup- or species-specific gene duplication events after speciation, resulting in divergent evolution. Taken together, genome-wide analysis of NLR family members provide insights into their evolutionary history in Solanaceae. These findings also provide important foundational knowledge for understanding NLR evolution and will empower broader characterization of disease resistance genes to be used for crop breeding.

Resistance to Tospoviruses in Vegetable Crops: Epidemiological and Molecular Aspects

During the past three decades, the economic impact of tospoviruses has increased, causing high yield losses in a variety of crops and ornamentals. Owing to the difficulty in combating thrips vectors with insecticides, the best way to limit/prevent tospovirus-induced diseases involves a management strategy that includes virus resistance. This review briefly presents current tospovirus taxonomy, diversity, molecular biology, and cytopathology as an introduction to a more extensive description of the two main resistance genes employed against tospoviruses: the Sw5 gene in tomato and the Tsw in pepper. Natural and experimental resistance-breaking (RB) isolates allowed the identification of the viral avirulence protein triggering each of the two resistance gene products; epidemiology of RB isolates is discussed to reinforce the need for allelic variants and the need to search for new/alternative resistance genes. Ongoing efforts for alternative resistance strategies are described not only for Tomato spotted wilt virus (TSWV) in pepper and tomato but also for other vegetable crops heavily impacted by tospoviruses.

Capsaicin: From Plants to a Cancer-Suppressing Agent

Capsaicinoids are plant secondary metabolites, capsaicin being the principal responsible for the pungency of chili peppers. It is biosynthesized through two pathways involved in phenylpropanoid and fatty acid metabolism. Plant capsaicin concentration is mainly affected by genetic, environmental and crop management factors. However, its synthesis can be enhanced by the use of elicitors. Capsaicin is employed as food additive and in pharmaceutical applications. Additionally, it has been found that capsaicin can act as a cancer preventive agent and shows wide applications against various types of cancer. This review is an approach in contextualizing the use of controlled stress on the plant to increase the content of capsaicin, highlighting its synthesis and its potential use as anticancer agent.

Transcriptome Analysis of Capsicum Chlorosis Virus-Induced Hypersensitive Resistance Response in Bell Capsicum

BACKGROUND

Capsicum chlorosis virus (CaCV) is an emerging pathogen of capsicum, tomato and peanut crops in Australia and South-East Asia. Commercial capsicum cultivars with CaCV resistance are not yet available, but CaCV resistance identified in Capsicum chinense is being introgressed into commercial Bell capsicum. However, our knowledge of the molecular mechanisms leading to the resistance response to CaCV infection is limited. Therefore, transcriptome and expression profiling data provide an important resource to better understand CaCV resistance mechanisms.

METHODOLOGY/PRINCIPAL FINDINGS

We assembled capsicum transcriptomes and analysed gene expression using Illumina HiSeq platform combined with a tag-based digital gene expression system. Total RNA extracted from CaCV/mock inoculated CaCV resistant (R) and susceptible (S) capsicum at the time point when R line showed a strong hypersensitive response to CaCV infection was used in transcriptome assembly. Gene expression profiles of R and S capsicum in CaCV- and buffer-inoculated conditions were compared. None of the genes were differentially expressed (DE) between R and S cultivars when mock-inoculated, while 2484 genes were DE when inoculated with CaCV. Functional classification revealed that the most highly up-regulated DE genes in R capsicum included pathogenesis-related genes, cell death-associated genes, genes associated with hormone-mediated signalling pathways and genes encoding enzymes involved in synthesis of defense-related secondary metabolites. We selected 15 genes to confirm DE expression levels by real-time quantitative PCR.

CONCLUSION/SIGNIFICANCE

DE transcript profiling data provided comprehensive gene expression information to gain an understanding of the underlying CaCV resistance mechanisms. Further, we identified candidate CaCV resistance genes in the CaCV-resistant C. annuum x C. chinense breeding line. This knowledge will be useful in future for fine mapping of the CaCV resistance locus and potential genetic engineering of resistance into CaCV-susceptible crops.

Phylogenetic relationships, diversification and expansion of chili peppers (Capsicum, Solanaceae)

BACKGROUND AND AIMS

Capsicum (Solanaceae), native to the tropical and temperate Americas, comprises the well-known sweet and hot chili peppers and several wild species. So far, only partial taxonomic and phylogenetic analyses have been done for the genus. Here, the phylogenetic relationships between nearly all taxa of Capsicum were explored to test the monophyly of the genus and to obtain a better knowledge of species relationships, diversification and expansion.

METHODS

Thirty-four of approximately 35 Capsicum species were sampled. Maximum parsimony and Bayesian inference analyses were performed using two plastid markers (matK and psbA-trnH) and one single-copy nuclear gene (waxy). The evolutionary changes of nine key features were reconstructed following the parsimony ancestral states method. Ancestral areas were reconstructed through a Bayesian Markov chain Monte Carlo analysis.

KEY RESULTS

Capsicum forms a monophyletic clade, with Lycianthes as a sister group, following both phylogenetic approaches. Eleven well-supported clades (four of them monotypic) can be recognized within Capsicum, although some interspecific relationships need further analysis. A few features are useful to characterize different clades (e.g. fruit anatomy, chromosome base number), whereas some others are highly homoplastic (e.g. seed colour). The origin of Capsicum is postulated in an area along the Andes of western to north-western South America. The expansion of the genus has followed a clockwise direction around the Amazon basin, towards central and south-eastern Brazil, then back to western South America, and finally northwards to Central America.

CONCLUSIONS

New insights are provided regarding interspecific relationships, character evolution, and geographical origin and expansion of Capsicum A clearly distinct early-diverging clade can be distinguished, centred in western-north-western South America. Subsequent rapid speciation has led to the origin of the remaining clades. The diversification of Capsicum has culminated in the origin of the main cultivated species in several regions of South to Central America.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

The rubber tree genome shows expansion of gene family associated with rubber biosynthesis

Hevea brasiliensis Muell. Arg, a member of the family Euphorbiaceae, is the sole natural resource exploited for commercial production of high-quality natural rubber. The properties of natural rubber latex are almost irreplaceable by synthetic counterparts for many industrial applications. A paucity of knowledge on the molecular mechanisms of rubber biosynthesis in high yield traits still persists. Here we report the comprehensive genome-wide analysis of the widely planted H. brasiliensis clone, RRIM 600. The genome was assembled based on ~155-fold combined coverage with Illumina and PacBio sequence data and has a total length of 1.55 Gb with 72.5% comprising repetitive DNA sequences. A total of 84,440 high-confidence protein-coding genes were predicted. Comparative genomic analysis revealed strong synteny between H. brasiliensis and other Euphorbiaceae genomes. Our data suggest that H. brasiliensis's capacity to produce high levels of latex can be attributed to the expansion of rubber biosynthesis-related genes in its genome and the high expression of these genes in latex. Using cap analysis gene expression data, we illustrate the tissue-specific transcription profiles of rubber biosynthesis-related genes, revealing alternative means of transcriptional regulation. Our study adds to the understanding of H. brasiliensis biology and provides valuable genomic resources for future agronomic-related improvement of the rubber tree.

Plant Genetic Background Increasing the Efficiency and Durability of Major Resistance Genes to Root-knot Nematodes Can Be Resolved into a Few Resistance QTLs

With the banning of most chemical nematicides, the control of root-knot nematodes (RKNs) in vegetable crops is now based essentially on the deployment of single, major resistance genes (R-genes). However, these genes are rare and their efficacy is threatened by the capacity of RKNs to adapt. In pepper, several dominant R-genes are effective against RKNs, and their efficacy and durability have been shown to be greater in a partially resistant genetic background. However, the genetic determinants of this partial resistance were unknown. Here, a quantitative trait loci (QTL) analysis was performed on the F2:3 population from the cross between Yolo Wonder, an accession considered partially resistant or resistant, depending on the RKN species, and Doux Long des Landes, a susceptible cultivar. A genetic linkage map was constructed from 130 F2 individuals, and the 130 F3 families were tested for resistance to the three main RKN species, Meloidogyne incognita, M. arenaria, and M. javanica. For the first time in the pepper-RKN pathosystem, four major QTLs were identified and mapped to two clusters. The cluster on chromosome P1 includes three tightly linked QTLs with specific effects against individual RKN species. The fourth QTL, providing specific resistance to M. javanica, mapped to pepper chromosome P9, which is known to carry multiple NBS-LRR repeats, together with major R-genes for resistance to nematodes and other pathogens. The newly discovered cluster on chromosome P1 has a broad spectrum of action with major additive effects on resistance. These data highlight the role of host QTLs involved in plant-RKN interactions and provide innovative potential for the breeding of new pepper cultivars or rootstocks combining quantitative resistance and major R-genes, to increase both the efficacy and durability of RKN control by resistance genes.

A high-quality carrot genome assembly provides new insights into carotenoid accumulation and asterid genome evolution

We report a high-quality chromosome-scale assembly and analysis of the carrot (Daucus carota) genome, the first sequenced genome to include a comparative evolutionary analysis among members of the euasterid II clade. We characterized two new polyploidization events, both occurring after the divergence of carrot from members of the Asterales order, clarifying the evolutionary scenario before and after radiation of the two main asterid clades. Large- and small-scale lineage-specific duplications have contributed to the expansion of gene families, including those with roles in flowering time, defense response, flavor, and pigment accumulation. We identified a candidate gene, DCAR_032551, that conditions carotenoid accumulation (Y) in carrot taproot and is coexpressed with several isoprenoid biosynthetic genes. The primary mechanism regulating carotenoid accumulation in carrot taproot is not at the biosynthetic level. We hypothesize that DCAR_032551 regulates upstream photosystem development and functional processes, including photomorphogenesis and root de-etiolation.

Short interspersed nuclear elements (SINEs) are abundant in Solanaceae and have a family-specific impact on gene structure and genome organization

Short interspersed nuclear elements (SINEs) are highly abundant non-autonomous retrotransposons that are widespread in plants. They are short in size, non-coding, show high sequence diversity, and are therefore mostly not or not correctly annotated in plant genome sequences. Hence, comparative studies on genomic SINE populations are rare. To explore the structural organization and impact of SINEs, we comparatively investigated the genome sequences of the Solanaceae species potato (Solanum tuberosum), tomato (Solanum lycopersicum), wild tomato (Solanum pennellii), and two pepper cultivars (Capsicum annuum). Based on 8.5 Gbp sequence data, we annotated 82 983 SINE copies belonging to 10 families and subfamilies on a base pair level. Solanaceae SINEs are dispersed over all chromosomes with enrichments in distal regions. Depending on the genome assemblies and gene predictions, 30% of all SINE copies are associated with genes, particularly frequent in introns and untranslated regions (UTRs). The close association with genes is family specific. More than 10% of all genes annotated in the Solanaceae species investigated contain at least one SINE insertion, and we found genes harbouring up to 16 SINE copies. We demonstrate the involvement of SINEs in gene and genome evolution including the donation of splice sites, start and stop codons and exons to genes, enlargement of introns and UTRs, generation of tandem-like duplications and transduction of adjacent sequence regions.

Structurally distinct Arabidopsis thaliana NLR immune receptors recognize tandem WY domains of an oomycete effector

Nucleotide-binding leucine-rich repeat (NB-LRR, or NLR) receptors mediate pathogen recognition. The Arabidopsis thaliana NLR RPP1 recognizes the tandem WY-domain effector ATR1 from the oomycete Hyaloperonospora arabidopsidis through direct association with C-terminal LRRs. We isolated and characterized homologous NLR genes RPP1-EstA and RPP1-ZdrA from two Arabidopsis ecotypes, Estland (Est-1) and Zdarec (Zdr-1), responsible for recognizing a novel spectrum of ATR1 alleles. RPP1-EstA and -ZdrA encode nearly identical NLRs that are phylogenetically distinct from known immunity-activating RPP1 homologs and possess greatly expanded LRR domains. Site-directed mutagenesis and truncation analysis of ATR1 suggests that these homologs recognize a novel surface of the 2(nd) WY domain of ATR1, partially specified by a C-terminal region of the LRR domain. Synteny comparison with RPP1 loci involved in hybrid incompatibility suggests that these functions evolved independently. Closely related RPP1 homologs have diversified their recognition spectra through LRR expansion and sequence variation, allowing them to detect multiple surfaces of the same pathogen effector. Engineering NLR receptor specificity may require a similar combination of repeat expansion and tailored amino acid variation.

Genome-Wide Identification and Analysis of the SBP-Box Family Genes under Phytophthora capsici Stress in Pepper (Capsicum annuum L.)

SQUAMOSA promoter binding protein (SBP)-box genes encode plant-specific transcription factors that are extensively involved in many physiological and biochemical processes, including growth, development, and signal transduction. However, pepper (Capsicum annuum L.) SBP-box family genes have not been well characterized. We investigated SBP-box family genes in the pepper genome and characterized these genes across both compatible and incompatible strain of Phytophthora capsici, and also under different hormone treatments. The results indicated that total 15 members were identified and distributed on seven chromosomes of pepper. Phylogenetic analysis showed that SBP-box genes of pepper can be classified into six groups. In addition, duplication analysis within pepper genome, as well as between pepper and Arabidopsis genomes demonstrated that there are four pairs of homology of SBP-box genes in the pepper genome and 10 pairs between pepper and Arabidopsis genomes. Tissue-specific expression analysis of the CaSBP genes demonstrated their diverse spatiotemporal expression patterns. The expression profiles were similarly analyzed following exposure to P. capsici inoculation and hormone treatments. It was shown that nine of the CaSBP genes (CaSBP01, 02, 03, 04, 05, 06, 11, 12, and 13) exhibited a dramatic up-regulation after compatible HX-9 strain (P. capsici) inoculation, while CaSBP09 and CaSBP15 were down-regulated. In case of PC strain (P. capsici) infection six of the CaSBP genes (CaSBP02, 05, 06, 11, 12, and 13) were arose while CaSBP14 was down regulated. Furthermore, Salicylic acid, Methyl jasmonate and their biosynthesis inhibitors treatment indicated that some of the CaSBP genes are potentially involved in these hormone regulation pathways. This genome-wide identification, as well as characterization of evolutionary relationships and expression profiles of the pepper CaSBP genes, will help to improve pepper stress tolerance in the future.

Neofunctionalization of Chromoplast Specific Lycopene Beta Cyclase Gene (CYC-B) in Tomato Clade

The ancestor of tomato underwent whole genome triplication ca. 71 Myr ago followed by widespread gene loss. However, few of the triplicated genes are retained in modern day tomato including lycopene beta cyclase that mediates conversion of lycopene to β-carotene. The fruit specific β-carotene formation is mediated by a chromoplast-specific paralog of lycopene beta cyclase (CYC-B) gene. Presently limited information is available about how the variations in CYC-B gene contributed to its neofunctionalization. CYC-B gene in tomato clade contained several SNPs and In-Dels in the coding sequence (33 haplotypes) and promoter region (44 haplotypes). The CYC-B gene coding sequence in tomato appeared to undergo purifying selection. The transit peptide sequence of CYC-B protein was predicted to have a stronger plastid targeting signal than its chloroplast specific paralog indicating a possible neofunctionalization. In promoter of two Bog (Beta old gold) mutants, a NUPT (nuclear plastid) DNA fragment of 256 bp, likely derived from a S. chilense accession, was present. In transient expression assay, this promoter was more efficient than the "Beta type" promoter. CARGATCONSENSUS box sequences are required for the binding of the MADS-box regulatory protein RIPENING INHIBITOR (RIN). The loss of CARGATCONSENSUS box sequence from CYC-B promoter in tomato may be related to attenuation of its efficiency to promote higher accumulation of β-carotene than lycopene during fruit ripening.

Assembly of the draft genome of buckwheat and its applications in identifying agronomically useful genes

Buckwheat (Fagopyrum esculentum Moench; 2n = 2x = 16) is a nutritionally dense annual crop widely grown in temperate zones. To accelerate molecular breeding programmes of this important crop, we generated a draft assembly of the buckwheat genome using short reads obtained by next-generation sequencing (NGS), and constructed the Buckwheat Genome DataBase. After assembling short reads, we determined 387,594 scaffolds as the draft genome sequence (FES_r1.0). The total length of FES_r1.0 was 1,177,687,305 bp, and the N50 of the scaffolds was 25,109 bp. Gene prediction analysis revealed 286,768 coding sequences (CDSs; FES_r1.0_cds) including those related to transposable elements. The total length of FES_r1.0_cds was 212,917,911 bp, and the N50 was 1,101 bp. Of these, the functions of 35,816 CDSs excluding those for transposable elements were annotated by BLAST analysis. To demonstrate the utility of the database, we conducted several test analyses using BLAST and keyword searches. Furthermore, we used the draft genome as a reference sequence for NGS-based markers, and successfully identified novel candidate genes controlling heteromorphic self-incompatibility of buckwheat. The database and draft genome sequence provide a valuable resource that can be used in efforts to develop buckwheat cultivars with superior agronomic traits.

Large-Scale Analyses of Angiosperm Nucleotide-Binding Site-Leucine-Rich Repeat Genes Reveal Three Anciently Diverged Classes with Distinct Evolutionary Patterns

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.

Mapping of a Novel Race Specific Resistance Gene to Phytophthora Root Rot of Pepper (Capsicum annuum) Using Bulked Segregant Analysis Combined with Specific Length Amplified Fragment Sequencing Strategy

Phytophthora root rot caused by Phytophthora capsici (P. capsici) is a serious limitation to pepper production in Southern China, with high temperature and humidity. Mapping PRR resistance genes can provide linked DNA markers for breeding PRR resistant varieties by molecular marker-assisted selection (MAS). Two BC1 populations and an F2 population derived from a cross between P. capsici-resistant accession, Criollo de Morelos 334 (CM334) and P. capsici-susceptible accession, New Mexico Capsicum Accession 10399 (NMCA10399) were used to investigate the genetic characteristics of PRR resistance. PRR resistance to isolate Byl4 (race 3) was controlled by a single dominant gene, PhR10, that was mapped to an interval of 16.39Mb at the end of the long arm of chromosome 10. Integration of bulked segregant analysis (BSA) and Specific Length Amplified Fragment sequencing (SLAF-seq) provided an efficient genetic mapping strategy. Ten polymorphic Simple Sequence Repeat (SSR) markers were found within this region and used to screen the genotypes of 636 BC1 plants, delimiting PhR10 to a 2.57 Mb interval between markers P52-11-21 (1.5 cM away) and P52-11-41 (1.1 cM). A total of 163 genes were annotated within this region and 31 were predicted to be associated with disease resistance. PhR10 is a novel race specific gene for PRR, and this paper describes linked SSR markers suitable for marker-assisted selection of PRR resistant varieties, also laying a foundation for cloning the resistance gene.

FISH-mapping of the 5S rDNA locus in chili peppers (Capsicum-Solanaceae)

We present here the physical mapping of the 5S rDNA locus in six wild and five cultivated taxa of Capsicum by means of a genus-specific FISH probe. In all taxa, a single 5S locus per haploid genome that persistently mapped onto the short arm of a unique metacentric chromosome pair at intercalar position, was found. 5S FISH signals of almost the same size and brightness intensity were observed in all the analyzed taxa. This is the first cytological characterization of the 5S in wild taxa of Capsicum by using a genus-derived probe, and the most exhaustive and comprehensive in the chili peppers up to now. The information provided here will aid the cytomolecular characterization of pepper germplasm to evaluate variability and can be instrumental to integrate physical, genetic and genomic maps already generated in the genus.

Evolutionary Profiling of Group II Pyridoxal-Phosphate-Dependent Decarboxylases Suggests Expansion and Functional Diversification of Histidine Decarboxylases in Tomato

Pyridoxal phosphate (PLP)-dependent enzymes are one of the most important enzymes involved in plant N metabolism. Here, we explored the evolution of group II PLP-dependent decarboxylases (PLP_deC), including aromatic L-amino acid decarboxylase, glutamate decarboxylase, and histidine decarboxylase in the plant lineage. Gene identification analysis revealed a higher number of genes encoding PLP_deC in higher plants than in lower plants. Expression profiling of PLP_deC orthologs and syntelogs in (L.) Heynh., pepper ( L.), and tomato ( L.) pointed toward conserved as well as distinct roles in developmental processes such as fruit maturation and ripening and abiotic stress responses. We further characterized a putative promoter of tomato ripening-associated gene () operating in a complex regulatory circuit. Our analysis provides a firm basis for further in-depth exploration of the PLP_deC gene family, particularly in the economically important Solanaceae family.

Evolution of plant genome architecture

We have witnessed an explosion in our understanding of the evolution and structure of plant genomes in recent years. Here, we highlight three important emergent realizations: (1) that the evolutionary history of all plant genomes contains multiple, cyclical episodes of whole-genome doubling that were followed by myriad fractionation processes; (2) that the vast majority of the variation in genome size reflects the dynamics of proliferation and loss of lineage-specific transposable elements; and (3) that various classes of small RNAs help shape genomic architecture and function. We illustrate ways in which understanding these organism-level and molecular genetic processes can be used for crop plant improvement.

CaCDPK15 positively regulates pepper responses to Ralstonia solanacearum inoculation and forms a positive-feedback loop with CaWRKY40 to amplify defense signaling

CaWRKY40 is a positive regulator of pepper (Capsicum annum) response to Ralstonia solanacearum inoculation (RSI), but the underlying mechanism remains largely unknown. Here, we functionally characterize CaCDPK15 in the defense signaling mediated by CaWRKY40. Pathogen-responsive TGA, W, and ERE boxes were identified in the CaCDPK15 promoter (pCaCDPK15), and pCaCDPK15-driven GUS expression was significantly enhanced in response to RSI and exogenously applied salicylic acid, methyl jasmonate, abscisic acid, and ethephon. Virus-induced gene silencing (VIGS) of CaCDPK15 significantly increased the susceptibility of pepper to RSI and downregulated the immunity-associated markers CaNPR1, CaPR1, and CaDEF1. By contrast, transient CaCDPK15 overexpression significantly activated hypersensitive response associated cell death, upregulated the immunity-associated marker genes, upregulated CaWRKY40 expression, and enriched CaWRKY40 at the promoters of its targets genes. Although CaCDPK15 failed to interact with CaWRKY40, the direct binding of CaWRKY40 to pCaCDPK15 was detected by chromatin immunoprecipitation, which was significantly potentiated by RSI in pepper plants. These combined results suggest that RSI in pepper induces CaCDPK15 and indirectly activates downstream CaWRKY40, which in turn potentiates CaCDPK15 expression. This positive-feedback loop would amplify defense signaling against RSI and efficiently activate strong plant immunity.

Dramatic Number Variation of R Genes in Solanaceae Species Accounted for by a Few R Gene Subfamilies

Most disease resistance genes encode nucleotide-binding-site (NBS) and leucine-rich-repeat (LRR) domains, and the NBS-LRR encoding genes are often referred to as R genes. Using newly developed approach, 478, 485, 1,194, 1,665, 2,042 and 374 R genes were identified from the genomes of tomato Heinz1706, wild tomato LA716, potato DM1-3, pepper Zunla-1 and wild pepper Chiltepin and tobacco TN90, respectively. The majority of R genes from Solanaceae were grouped into 87 subfamilies, including 16 TIR-NBS-LRR (TNL) and 71 non-TNL subfamilies. Each subfamily was annotated manually, including identification of intron/exon structure and intron phase. Interestingly, TNL subfamilies have similar intron phase patterns, while the non-TNL subfamilies have diverse intron phase due to frequent gain of introns. Prevalent presence/absence polymorphic R gene loci were found among Solanaceae species, and an integrated map with 427 R loci was constructed. The pepper genome (2,042 in Chiltepin) has at least four times of R genes as in tomato (478 in Heinz1706). The high number of R genes in pepper genome is due to the amplification of R genes in a few subfamilies, such as the Rpi-blb2 and BS2 subfamilies. The mechanism underlying the variation of R gene number among different plant genomes is discussed.

The genome sequence of the outbreeding globe artichoke constructed de novo incorporating a phase-aware low-pass sequencing strategy of F1 progeny

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F₁ individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.

A New Ethylene-Responsive Factor CaPTI1 Gene of Pepper (Capsicum annuum L.) Involved in the Regulation of Defense Response to Phytophthora capsici

Ethylene-responsive factors (ERF) are usually considered to play diverse roles in plant response to biotic and abiotic stresses. In this study, an ERF gene CaPTI1 was isolated from pepper transcriptome database. CaPTI1 contains an open reading frame (ORF) of 543 bp, which encodes a putative polypeptide of 180 amino acids with a theoretical molecular weight of 20.30 kDa. Results of expression profile showed that CaPTI1 had a highest expression level in roots and this gene could not only response to the infection of Phytophthora capsici and the stresses of cold and drought, but also be induced by the signaling molecule (salicylic acid, Methyl Jasmonate, Ethephon, and hydogen peroxide). Furthermore, virus-induce gene silencing (VIGS) of CaPTI1 in pepper weakened the defense response significantly by reducing the expression of defense related genes CaPR1, CaDEF1 and CaSAR82 and also the root activity. These results suggested that CaPTI1 is involved in the regulation of defense response to P. capsici in pepper.

A comprehensive characterization of simple sequence repeats in pepper genomes provides valuable resources for marker development in Capsicum

The sequences of the full set of pepper genomes including nuclear, mitochondrial and chloroplast are now available for use. However, the overall of simple sequence repeats (SSR) distribution in these genomes and their practical implications for molecular marker development in Capsicum have not yet been described. Here, an average of 868,047.50, 45.50 and 30.00 SSR loci were identified in the nuclear, mitochondrial and chloroplast genomes of pepper, respectively. Subsequently, systematic comparisons of various species, genome types, motif lengths, repeat numbers and classified types were executed and discussed. In addition, a local database composed of 113,500 in silico unique SSR primer pairs was built using a homemade bioinformatics workflow. As a pilot study, 65 polymorphic markers were validated among a wide collection of 21 Capsicum genotypes with allele number and polymorphic information content value per marker raging from 2 to 6 and 0.05 to 0.64, respectively. Finally, a comparison of the clustering results with those of a previous study indicated the usability of the newly developed SSR markers. In summary, this first report on the comprehensive characterization of SSR motifs in pepper genomes and the very large set of SSR primer pairs will benefit various genetic studies in Capsicum.

An ultra-high-density bin map facilitates high-throughput QTL mapping of horticultural traits in pepper (Capsicum annuum)

Most agricultural traits are controlled by quantitative trait loci (QTLs); however, there are few studies on QTL mapping of horticultural traits in pepper (Capsicum spp.) due to the lack of high-density molecular maps and the sequence information. In this study, an ultra-high-density map and 120 recombinant inbred lines (RILs) derived from a cross between C. annuum ‘Perennial’ and C. annuum ‘Dempsey’ were used for QTL mapping of horticultural traits. Parental lines and RILs were resequenced at 18× and 1× coverage, respectively. Using a sliding window approach, an ultra-high-density bin map containing 2,578 bins was constructed. The total map length of the map was 1,372 cM, and the average interval between bins was 0.53 cM. A total of 86 significant QTLs controlling 17 horticultural traits were detected. Among these, 32 QTLs controlling 13 traits were major QTLs. Our research shows that the construction of bin maps using low-coverage sequence is a powerful method for QTL mapping, and that the short intervals between bins are helpful for fine-mapping of QTLs. Furthermore, bin maps can be used to improve the quality of reference genomes by elucidating the genetic order of unordered regions and anchoring unassigned scaffolds to linkage groups.