Seungill Kim

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

Hot pepper (Capsicum annuum), one of the oldest domesticated crops in the Americas, is the most widely grown spice crop in the world. We report whole-genome sequencing and assembly of the hot pepper (Mexican landrace of Capsicum annuum cv. CM334) at 186.6× coverage. We also report resequencing of two cultivated peppers and de novo sequencing of the wild species Capsicum chinense. The genome size of the hot pepper was approximately fourfold larger than that of its close relative tomato, and the genome showed an accumulation of Gypsy and Caulimoviridae family elements. Integrative genomic and transcriptomic analyses suggested that change in gene expression and neofunctionalization of capsaicin synthase have shaped capsaicinoid biosynthesis. We found differential molecular patterns of ripening regulators and ethylene synthesis in hot pepper and tomato. The reference genome will serve as a platform for improving the nutritional and medicinal values of Capsicum species.

RNA-seq pinpoints a Xanthomonas TAL-effector activated resistance gene in a large-crop genome

Transcription activator-like effector (TALE) proteins of the plant pathogenic bacterial genus Xanthomonas bind to and transcriptionally activate host susceptibility genes, promoting disease. Plant immune systems have taken advantage of this mechanism by evolving TALE binding sites upstream of resistance (R) genes. For example, the pepper Bs3 and rice Xa27 genes are hypersensitive reaction plant R genes that are transcriptionally activated by corresponding TALEs. Both R genes have a hallmark expression pattern in which their transcripts are detectable only in the presence and not the absence of the corresponding TALE. By transcriptome profiling using next-generation sequencing (RNA-seq), we tested whether we could avoid laborious positional cloning for the isolation of TALE-induced R genes. In a proof-of-principle experiment, RNA-seq was used to identify a candidate for Bs4C, an R gene from pepper that mediates recognition of the Xanthomonas TALE protein AvrBs4. We identified one major Bs4C candidate transcript by RNA-seq that was expressed exclusively in the presence of AvrBs4. Complementation studies confirmed that the candidate corresponds to the Bs4C gene and that an AvrBs4 binding site in the Bs4C promoter directs its transcriptional activation. Comparison of Bs4C with a nonfunctional allele that is unable to recognize AvrBs4 revealed a 2-bp polymorphism within the TALE binding site of the Bs4C promoter. Bs4C encodes a structurally unique R protein and Bs4C-like genes that are present in many solanaceous genomes seem to be as tightly regulated as pepper Bs4C. These findings demonstrate that TALE-specific R genes can be cloned from large-genome crops with a highly efficient RNA-seq approach.

Comparative analysis of pepper and tomato reveals euchromatin expansion of pepper genome caused by differential accumulation of Ty3/Gypsy-like elements

BackgroundAmong the Solanaceae plants, the pepper genome is three times larger than that of tomato. Although the gene repertoire and gene order of both species are well conserved, the cause of the genome-size difference is not known. To determine the causes for the expansion of pepper euchromatic regions, we compared the pepper genome to that of tomato.ResultsFor sequence-level analysis, we generated 35.6 Mb of pepper genomic sequences from euchromatin enriched 1,245 pepper BAC clones. The comparative analysis of orthologous gene-rich regions between both species revealed insertion of transposons exclusively in the pepper sequences, maintaining the gene order and content. The most common type of the transposon found was the LTR retrotransposon. Phylogenetic comparison of the LTR retrotransposons revealed that two groups of Ty3/Gypsy-like elements (Tat and Athila) were overly accumulated in the pepper genome. The FISH analysis of the pepper Tat elements showed a random distribution in heterochromatic and euchromatic regions, whereas the tomato Tat elements showed heterochromatin-preferential accumulation.ConclusionsCompared to tomato pepper euchromatin doubled its size by differential accumulation of a specific group of Ty3/Gypsy-like elements. Our results could provide an insight on the mechanism of genome evolution in the Solanaceae family.

The Hot Pepper (Capsicum annuum) MicroRNA Transcriptome Reveals Novel and Conserved Targets: A Foundation for Understanding MicroRNA Functional Roles in Hot Pepper

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nt in length which play important roles in regulating gene expression in plants. Although many miRNA studies have focused on a few model plants, miRNAs and their target genes remain largely unknown in hot pepper (Capsicum annuum), one of the most important crops cultivated worldwide. Here, we employed high-throughput sequencing technology to identify miRNAs in pepper extensively from 10 different libraries, including leaf, stem, root, flower, and six developmental stage fruits. Based on a bioinformatics pipeline, we successfully identified 29 and 35 families of conserved and novel miRNAs, respectively. Northern blot analysis was used to validate further the expression of representative miRNAs and to analyze their tissue-specific or developmental stage-specific expression patterns. Moreover, we computationally predicted miRNA targets, many of which were experimentally confirmed using 5′ rapid amplification of cDNA ends analysis. One of the validated novel targets of miR-396 was a domain rearranged methyltransferase, the major de novo methylation enzyme, involved in RNA-directed DNA methylation in plants. This work provides the first reliable draft of the pepper miRNA transcriptome. It offers an expanded picture of pepper miRNAs in relation to other plants, providing a basis for understanding the functional roles of miRNAs in pepper.

Evolution of the large genome in Capsicum annuum occurred through accumulation of single‐type long terminal repeat retrotransposons and their derivatives

Although plant genome sizes are extremely diverse, the mechanism underlying the expansion of huge genomes that did not experience whole-genome duplication has not been elucidated. The pepper, Capsicum annuum, is an excellent model for studies of genome expansion due to its large genome size (2700 Mb) and the absence of whole genome duplication. As most of the pepper genome structure has been identified as constitutive heterochromatin, we investigated the evolution of this region in detail. Our findings show that the constitutive heterochromatin in pepper was actively expanded 20.0-7.5 million years ago through a massive accumulation of single-type Ty3/Gypsy-like elements that belong to the Del subgroup. Interestingly, derivatives of the Del elements, such as non-autonomous long terminal repeat retrotransposons and long-unit tandem repeats, played important roles in the expansion of constitutive heterochromatic regions. This expansion occurred not only in the existing heterochromatic regions but also into the euchromatic regions. Furthermore, our results revealed a repeat of unit length 18-24 kb. This repeat was found not only in the pepper genome but also in the other solanaceous species, such as potato and tomato. These results represent a characteristic mechanism for large genome evolution in plants.

Integrative structural annotation of de novo RNA-Seq provides an accurate reference gene set of the enormous genome of the onion (Allium cepa L.)

The onion (Allium cepa L.) is one of the most widely cultivated and consumed vegetable crops in the world. Although a considerable amount of onion transcriptome data has been deposited into public databases, the sequences of the protein-coding genes are not accurate enough to be used, owing to non-coding sequences intermixed with the coding sequences. We generated a high-quality, annotated onion transcriptome from de novo sequence assembly and intensive structural annotation using the integrated structural gene annotation pipeline (ISGAP), which identified 54,165 protein-coding genes among 165,179 assembled transcripts totalling 203.0 Mb by eliminating the intron sequences. ISGAP performed reliable annotation, recognizing accurate gene structures based on reference proteins, and ab initio gene models of the assembled transcripts. Integrative functional annotation and gene-based SNP analysis revealed a whole biological repertoire of genes and transcriptomic variation in the onion. The method developed in this study provides a powerful tool for the construction of reference gene sets for organisms based solely on de novo transcriptome data. Furthermore, the reference genes and their variation described here for the onion represent essential tools for molecular breeding and gene cloning in Allium spp.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

BackgroundTransposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.ResultsWe report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.ConclusionsOur study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

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.

De Novo Transcriptome Analysis to Identify Anthocyanin Biosynthesis Genes Responsible for Tissue-Specific Pigmentation in Zoysiagrass (Zoysia japonica Steud.).

Zoysiagrass (Zoysia japonica Steud.) is commonly found in temperate climate regions and widely used for lawns, in part, owing to its uniform green color. However, some zoysiagrass cultivars accumulate red to purple pigments in their spike and stolon tissues, thereby decreasing the aesthetic value. Here we analyzed the anthocyanin contents of two zoysiagrass cultivars ‘Anyang-jungji’ (AJ) and ‘Greenzoa’ (GZ) that produce spikes and stolons with purple and green colors, respectively, and revealed that cyanidin and petunidin were primarily accumulated in the pigmented tissues. In parallel, we performed a de novo transcriptome assembly and identified differentially expressed genes between the two cultivars. We found that two anthocyanin biosynthesis genes encoding anthocyanidin synthase (ANS) and dihydroflavonol 4-reductase (DFR) were preferentially upregulated in the purple AJ spike upon pigmentation. Both ANS and DFR genes were also highly expressed in other zoysiagrass cultivars with purple spikes and stolons, but their expression levels were significantly low in the cultivars with green tissues. We observed that recombinant ZjDFR1 and ZjANS1 proteins successfully catalyze the conversions of dihydroflavonols into leucoanthocyanidins and leucoanthocyanidins into anthocyanidins, respectively. These findings strongly suggest that upregulation of ANS and DFR is responsible for tissue-specific anthocyanin biosynthesis and differential pigmentation in zoysiagrass. The present study also demonstrates the feasibility of a de novo transcriptome analysis to identify the key genes associated with specific traits, even in the absence of reference genome information.

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.