Whole-genome resequencing analysis of the medicinal plant Gardenia jasminoides.

Background: Gardenia jasminoides is a species of Chinese medicinal plant, which has high medicinal and economic value and rich genetic diversity, but the study on its genetic diversity is far not enough. Methods: In this study, one wild and one cultivated gardenia materials were resequenced using IlluminaHiSeq sequencing platform and the data were evaluated to understand the genomic characteristics of G. jasminoides. Results: After data analysis, the results showed that clean data of 11.77G, Q30 reached 90.96%. The average comparison rate between the sample and reference genome was 96.08%, the average coverage depth was 15X, and the genome coverage was 85.93%. The SNPs of FD and YP1 were identified, and 3,087,176 and 3,241,416 SNPs were developed, respectively. In addition, SNP non-synonymous mutation, InDel mutation, SV mutation and CNV mutation were also detected between the sample and the reference genome, and KEGG, GO and COG database annotations were made for genes with DNA level variation. The structural gene variation in the biosynthetic pathway of crocin and gardenia, the main medicinal substance of G. jasminoides was further explored, which provided basic data for molecular breeding and genetic diversity of G. jasminoides in the future.

Development and validation of a custom panel including 114 InDels using massively parallel sequencing for forensic application.

Insertion/deletion polymorphisms (InDels) have particular characteristics, such as a relatively low mutation rate, small amplicon size, and no stutter artifacts when genotyped via the capillary electrophoresis platform. It would be an important complementary tool for individual identification and certain kinship analyses. At present, massively parallel sequencing (MPS) has shown excellent application value in forensic studies. Therefore, in this study, we developed a custom MPS InDel panel that contains 114 InDels [77 autosomal InDels (A-InDels), 32 X-chromosomal InDels (X-InDels), and 5 Y-chromosomal InDels) based on previous studies. To assess this panel's performance, several validation experiments were performed, including sensitivity, inhibitor, degraded DNA testing, species specificity, concordance, repeatability, case-type samples, and population studies. The results showed that the lowest DNA input was 0.25 ng. All genotypes were obtained in the presence of 80 ng/µL humic acid, 2000 µmol/L calcium, 3000 µmol/L EDTA and indigo. In degraded DNA testing, 90% of loci could be detected for 16-day-old formalin-fixed hearts. In addition, this panel has good species specificity. The values of combined power of discrimination and the combined power of exclusion for 77 A-InDels were 1-3.9951 × 10-32 and 1-4.2956 × 10-7 , respectively. The combined mean exclusion chance for 32 X-InDels was 0.99999 in trios and 0.99904 in duos. The validation results indicate that this newly developed MPS multiplex system is a robust tool for forensic applications.

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica.

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

A Genomic Evaluation of Six Selected Inbred Lines of the Naturalized Plants of Milk Thistle (Silybum marianum L. Gaertn.) in Korea.

Milk thistle (Silybum marianum) belongs to the Asteraceae family and is a medicinal plant native to the Mediterranean Basin. Silymarin in achene is a widely used herbal product for chronic liver disease. There is growing interest in natural medicine using milk thistle in Korea, but the raw material completely relies on imports. Despite its economic importance, phenotypic evaluations of native resources of milk thistle in Korea have not been carried out. In addition, genomic research and molecular marker development are very limited in milk thistle. In this study, we evaluated 220 milk thistle resources consisting of 172 accessions collected from the domestic market, and 48 accessions isolated from 6 accessions distributed by the National Agrobiodiversity Center in Korea. Six plant characteristics (height, seed weight, number of flowers, seed weight per flower, spine length, and color at harvest) were measured, and six samples (M01-M06) were selected to represent the genetic diversity of the population for genomic research. To develop PCR-based and co-dominant insertion/deletion (InDel) markers, we performed genome-wide InDel detection by comparing the whole-genome resequencing data of the six selected accessions with the reference genome sequence (GCA_001541825). As a result, 177 InDel markers with high distinguishability and reproducibility were selected from the 30,845 InDel variants. Unknowingly imported alien plant resources could easily be genetically mixed, and jeopardized seed purity can cause continuous difficulties in the development of high value-added agricultural platforms utilizing natural products. The selected plant materials and 177 validated InDel markers developed via whole-genome resequencing analysis could be valuable resources for breeding, conservation, and ecological studies of natives to Korea, along with acceleration of Silybum marianum industrialization.

QTL mapping and candidate gene analysis for yield and grain weight/size in Tartary buckwheat.

BACKGROUND: Grain weight/size influences not only grain yield (GY) but also nutritional and appearance quality and consumer preference in Tartary buckwheat. The identification of quantitative trait loci (QTLs)/genes for grain weight/size is an important objective of Tartary buckwheat genetic research and breeding programs. RESULTS: Herein, we mapped the QTLs for GY, 1000-grain weight (TGW), grain length (GL), grain width (GW) and grain length-width ratio (L/W) in four environments using 221 recombinant inbred lines (XJ-RILs) derived from a cross of 'Xiaomiqiao × Jinqiaomai 2'. In total, 32 QTLs, including 7 for GY, 5 for TGW, 6 for GL, 11 for GW and 3 for L/W, were detected and distributed in 24 genomic regions. Two QTL clusters, qClu-1-3 and qClu-1-5, located on chromosome Ft1, were revealed to harbour 7 stable major QTLs for GY (qGY1.2), TGW (qTGW1.2), GL (qGL1.1 and qGL1.4), GW (qGW1.7 and qGW1.10) and L/W (qL/W1.2) repeatedly detected in three and above environments. A total of 59 homologues of 27 known plant grain weight/size genes were found within the physical intervals of qClu-1-3 and qClu-1-5. Six homologues, FtBRI1, FtAGB1, FtTGW6, FtMADS1, FtMKK4 and FtANT, were identified with both non-synonymous SNP/InDel variations and significantly differential expression levels between the two parents, which may play important roles in Tatary buckwheat grain weight/size control and were chosen as core candidate genes for further investigation. CONCLUSIONS: Two stable major QTL clusters related to grain weight/size and six potential key candidate genes were identified by homology comparison, SNP/InDel variations and qRT‒qPCR analysis between the two parents. Our research provides valuable information for improving grain weight/size and yield in Tartary buckwheat breeding.

Functional inhibition of the StERF3 gene by dual targeting through CRISPR/Cas9 enhances resistance to the late blight disease in Solanum tuberosum L.

BACKGROUND: Disease-resistant cultivars are the best solution to get their maximum yield potential and avoid fungicide application. There is no doubt about the contribution, and use of R genes (resistance genes) in resistance development in plants, while S genes (susceptibility genes) also hold a strong position in pathogenesis by resistance repression, and their loss of function contributes to enhanced resistance. Hence, we attempted to knock out the function of the StERF3 gene in potatoes through CRISPR/Cas9-based genome editing and investigated the CRISPR/Cas9 approach as strategic control against late blight disease in potato plants. METHODS AND RESULTS: The StERF3 gene was edited in late blight susceptible cv. Lady Rosetta. Full allelic edited plants were identified through DnpI, and N1aIV mediated restriction digestion and then further analyzed through Indel Detection by Amplicon Analysis. Sequence analysis of targeted plants for indel identification showed full allelic editing. The detached leaf assay of full allelic edited plants demonstrated the role of the StERF3 gene in susceptibility to late blight in potatoes. In planta disease assay also showed reduced, slowed, and delayed disease progression in StERF3-loss-of-function mutants compared to wild-type (control) plants. Less fungal biomass was quantified in knockouts through Real-time qPCR that supported less susceptibility of edited plants to late blight. Besides, relatively high expression of pathogens-related genes, StPR1, and StNPR1, were also observed in StERF3-loss-of-function mutants compared to the corresponding control. CONCLUSION: The results showed the functional inhibition of StERF3 genes using the CRISPR/Cas9 approach. The functional knockouts (StERF3 gene-edited potato plants) revealed enhanced resistance against Phytophthora infestans, thereby demonstrating the best strategic control for late blight disease in potato plants.

Heterogeneous Genetic Diversity Estimation of a Promising Domestication Medicinal Motherwort Leonurus Cardiaca Based on Chloroplast Genome Resources.

Leonurus cardiaca has a long history of use in western herbal medicine and is applied for the treatment of gynaecological conditions, anxiety, and heart diseases. Because of its botanical relationship to the primary Chinese species, L. japonicus, and extensive medical indications that go beyond the traditional indications for the Chinese species, it is a promising medicinal resource. Therefore, the features of genetic diversity and variability in the species have been prioritized. To explore these issues, we sequenced the chloroplast genomes of 22 accessions of L. cardiaca from different geographical locations worldwide using high-throughput sequencing. The results indicate that L. cardiaca has a typical quadripartite structure and range from 1,51,236 bp to 1,51,831 bp in size, forming eight haplotypes. The genomes all contain 114 distinct genes, including 80 protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes. Comparative analysis showed abundant diversity of single nucleotide polymorphisms (SNPs), indels, simple sequence repeats (SSRs) in 22 accessions. Codon usage showed highly similar results for L. cardiaca species. The phylogenetic and network analysis indicated 22 accessions forming four clades that were partly related to the geographical distribution. In summary, our study highlights the advantage of chloroplast genome with large data sets in intraspecific diversity evaluation and provides a new tool to facilitate medicinal plant conservation and domestication.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea.

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

Chloroplast genome variation and phylogenetic relationships of Atractylodes species.

BACKGROUND: Atractylodes DC is the basic original plant of the widely used herbal medicines "Baizhu" and "Cangzhu" and an endemic genus in East Asia. Species within the genus have minor morphological differences, and the universal DNA barcodes cannot clearly distinguish the systemic relationship or identify the species of the genus. In order to solve these question, we sequenced the chloroplast genomes of all species of Atractylodes using high-throughput sequencing. RESULTS: The results indicate that the chloroplast genome of Atractylodes has a typical quadripartite structure and ranges from 152,294 bp (A. carlinoides) to 153,261 bp (A. macrocephala) in size. The genome of all species contains 113 genes, including 79 protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes. Four hotspots, rpl22-rps19-rpl2, psbM-trnD, trnR-trnT(GGU), and trnT(UGU)-trnL, and a total of 42-47 simple sequence repeats (SSR) were identified as the most promising potentially variable makers for species delimitation and population genetic studies. Phylogenetic analyses of the whole chloroplast genomes indicate that Atractylodes is a clade within the tribe Cynareae; Atractylodes species form a monophyly that clearly reflects the relationship within the genus. CONCLUSIONS: Our study included investigations of the sequences and structural genomic variations, phylogenetics and mutation dynamics of Atractylodes chloroplast genomes and will facilitate future studies in population genetics, taxonomy and species identification.

Comparative Analysis of Actaea Chloroplast Genomes and Molecular Marker Development for the Identification of Authentic Cimicifugae Rhizoma.

Actaea (Ranunculaceae; syn. Cimicifuga) is a controversial and complex genus. Dried rhizomes of Actaea species are used as Korean traditional herbal medicine. Although Actaea species are valuable, given their taxonomic classification and medicinal properties, sequence information of Actaea species is limited. In this study, we determined the complete chloroplast (cp) genome sequences of three Actaea species, including A. simplex, A. dahurica, and A. biternata. The cp genomes of these species varied in length from 159,523 to 159,789 bp and contained 112 unique functional genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene order, orientation, and content were well conserved in the three cp genomes. Comparative sequence analysis revealed the presence of hotspots, including ndhC-trnV-UAC, in Actaea cp genomes. High-resolution phylogenetic relationships were established among Actaea species based on cp genome sequences. Actaea species were clustered into each Actaea section, consistent with the Angiosperm Phylogeny Group (APG) IV system of classification. We also developed a novel indel marker, based on copy number variation of tandem repeats, to facilitate the authentication of the herbal medicine Cimicifugae Rhizoma. The availability Actaea cp genomes will provide abundant information for the taxonomic and phylogenetic analyses of Actaea species, and the Actaea (ACT) indel marker will be useful for the authentication of the herbal medicine.

Development and Application of InDel Markers for Authentication of the Korean Herbs Zanthoxylum schinifolium and Zanthoxylum piperitum.

Zanthoxylum schinifolium and Zanthoxylum piperitum are the sources of the well-known traditional Korean herbal medicines "sancho" (prickly ash) and "chopi" (Korean pepper), respectively. Sancho and chopi are often indiscriminately mixed due to the similar appearance of the herbal materials when used as spices and herbal medicines. Moreover, commercial sancho and chopi products often contain adulterants, which is insufficient to ensure food efficacy and safety. In this study, we developed hypervariable insertion/deletion (InDel) markers to distinguish between sancho and chopi products by comparing the complete chloroplast genome sequences of four Zanthoxylum species deposited in the National Center for Biotechnology Information (NCBI) GenBank. Comparative analyses of the nucleotide diversity (Pi) of these Zanthoxylum genomes revealed four hypervariable divergent sites (trnH-psbA, psbZ-trnG, trnfM-rps14, and trnF-ndhK) with Pi > 0.025 among 520 windows. Of these four regions, including two genic and two intergenic regions, only psbZ-trnG yielded accurate PCR amplification results between commercial sancho and chopi products from the Korean herbal medicine market. We therefore selected psbZ-trnG, an InDel-variable locus with high discriminatory powers, as a candidate DNA barcode locus. This InDel marker could be used as a valuable, simple, and efficient tool for identifying these medicinal herbs, thereby increasing the safety of these spices and herbal materials in the food market.

Characterization of genome-wide genetic variations between two varieties of tea plant (Camellia sinensis) and development of InDel markers for genetic research.

BACKGROUND: Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) are the major genetic variations and are distributed extensively across the whole plant genome. However, few studies of these variations have been conducted in the long-lived perennial tea plant. RESULTS: In this study, we investigated the genome-wide genetic variations between Camellia sinensis var. sinensis 'Shuchazao' and Camellia sinensis var. assamica 'Yunkang 10', identified 7,511,731 SNPs and 255,218 InDels based on their whole genome sequences, and we subsequently analyzed their distinct types and distribution patterns. A total of 48 InDel markers that yielded polymorphic and unambiguous fragments were developed when screening six tea cultivars. These markers were further deployed on 46 tea cultivars for transferability and genetic diversity analysis, exhibiting information with an average 4.02 of the number of alleles (Na) and 0.457 of polymorphism information content (PIC). The dendrogram showed that the phylogenetic relationships among these tea cultivars are highly consistent with their genetic backgrounds or original places. Interestingly, we observed that the catechin/caffeine contents between 'Shuchazao' and 'Yunkang 10' were significantly different, and a large number of SNPs/InDels were identified within catechin/caffeine biosynthesis-related genes. CONCLUSION: The identified genome-wide genetic variations and newly-developed InDel markers will provide a valuable resource for tea plant genetic and genomic studies, especially the SNPs/InDels within catechin/caffeine biosynthesis-related genes, which may serve as pivotal candidates for elucidating the molecular mechanism governing catechin/caffeine biosynthesis.

The Complete Chloroplast Genomes of Six Ipomoea Species and Indel Marker Development for the Discrimination of Authentic Pharbitidis Semen (Seeds of I. nil or I. purpurea).

Ipomoea L. is the largest genus within the Convolvulaceae and contains 600-700 species. Ipomoea species (morning glories) are economically valuable as horticultural species and scientifically valuable as ecological model plants to investigate mating systems, molecular evolution, and both plant-herbivore and plant-parasite interactions. Furthermore, the dried seeds of I. nil or I. purpurea are used in Korean traditional herbal medicines. In this study, chloroplast (cp) genomes were sequenced from six Ipomoea species, namely, I. nil and I. purpurea and, for the first time, I. triloba, I. lacunosa, I. hederacea, and I. hederacea var. integriuscula. The cp genomes were 161,354-161,750 bp in length and exhibited conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 transfer RNA genes, and 4 ribosomal RNA genes. The gene order, content, and orientation of the six Ipomoea cp genomes were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the six Ipomoea cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (petN-psbM, trnI-CAU-ycf2, ndhH-ndhF, psbC-trnS, and ccsA-ndhD). In addition, the protein-coding genes accD, cemA, and ycf2 exhibited high sequence variability and were under positive selection (Ka/Ks > 1), indicating adaptive evolution to the environment within the Ipomoea genus. Phylogenetic analysis of the six Ipomoea species revealed that these species clustered according to the APG IV system. In particular, I. nil and I. hederacea had monophyletic positions, with I. purpurea as a sister. I. triloba and I. lacunosa in the section Batatas and I. hederacea and I. hederacea var. integriuscula in the section Quamoclit were supported in this study with strong bootstrap values and posterior probabilities. We uncovered high-resolution phylogenetic relationships between Ipomoeeae. Finally, indel markers (IPOTY and IPOYCF) were developed for the discrimination of the important herbal medicine species I. nil and I. purpurea. The cp genomes and analyses in this study provide useful information for taxonomic, phylogenetic, and evolutionary analysis of the Ipomoea genome, and the indel markers will be useful for authentication of herbal medicines.

Molecular Discrimination of Cynanchum wilfordii and Cynanchum auriculatum by InDel Markers of Chloroplast DNA.

The tuber of Cynanchum wilfordii (Baekshuoh Radix in Korean) is an important medicinal herb in Korea and China; however, it is difficult to differentiate C. wilfordii from a related medicinal herb, C. auriculatum (Baishouwu Radix in Chinese). We sought to develop a molecular method that could be used to distinguish between the tubers of C. wilfordii and C. auriculatum. We aligned the chloroplast genome sequences (available in the NCBI database) of the two species and identified three species-specific insertion and deletion (InDel) sites in the trnQ-psbK, rps2-rpoC2, and psaJ-rpl33 intergenic spacer (IGS) regions. To confirm the presence of these three InDels and validate their use as markers, we designed three primer pairs to amplify the trnQ-psbK, rps2-rpoC2, and psaJ-rpl33 IGS regions. Polymerase chain reaction (PCR) amplification of the trnQ-psbK IGS region yielded a 249 bp fragment for C. wilfordii, and 419 bp fragment for C. auriculatum, whereas the rps2-rpoC2 IGS primers produced a 629 bp fragment from C. wilfordii and a 282 bp fragment from C. auriculatum. In the psaJ-rpl33 IGS region, allele fragments of 342 and 360 bp in length were amplified from C. wilfordii, whereas 249 and 250 bp fragment were amplified from C. auriculatum. We propose these three InDel markers as a valuable, simple, and efficient tool for identifying these medicinal herbs and will thus reduce adulteration of these herbal materials in commercial markets.

[Development of indel markers for molecular authentication of Panax ginseng and P. quinquefolius].

Panax ginseng and P. quinquefolius are two kinds of important medicinal herbs. They are morphologically similar but have different pharmacological effects. Therefore, botanical origin authentication of these two ginsengs is of great importance for ensuring pharmaceutical efficacy and food safety. Based on the fact that intron position in orthologous genes is highly conserved across plant species, intron length polymorphisms were exploited from unigenes of ginseng. Specific primers were respectively designed for these two species based on their insertion/deletion sequences of cytochrome P450 and glyceraldehyde 3-phosphate dehydrogenase, and multiplex PCR was conducted for molecular authentication of P.ginseng and P. quinquefolius. The results showed that the developed multiplex PCR assay was effective for molecular authentication of P.ginseng and P. quinquefolius without strict PCR condition and the optimization of reaction system.This study provides a preferred ideal marker system for molecular authentication of ginseng,and the presented method can be employed in origin authentication of other herbal preparations.

Development of indel markers for molecular authentication of Panax ginseng and P. quinquefolius / 中国中药杂志

Panax ginseng and P. quinquefolius are two kinds of important medicinal herbs. They are morphologically similar but have different pharmacological effects. Therefore, botanical origin authentication of these two ginsengs is of great importance for ensuring pharmaceutical efficacy and food safety. Based on the fact that intron position in orthologous genes is highly conserved across plant species, intron length polymorphisms were exploited from unigenes of ginseng. Specific primers were respectively designed for these two species based on their insertion/deletion sequences of cytochrome P450 and glyceraldehyde 3-phosphate dehydrogenase, and multiplex PCR was conducted for molecular authentication of P.ginseng and P. quinquefolius. The results showed that the developed multiplex PCR assay was effective for molecular authentication of P.ginseng and P. quinquefolius without strict PCR condition and the optimization of reaction system.This study provides a preferred ideal marker system for molecular authentication of ginseng,and the presented method can be employed in origin authentication of other herbal preparations.

Viability and genetic stability of potato spindle tuber viroid mutants with indels in specific loops of the rod-like secondary structure.

Maintenance of the rod-like structure of potato spindle tuber viroid (PSTVd), which contains over 20 loops and bulges between double-stranded helices, is important for viroid biology. To study tolerance to modifications of the stem-loop structures and PSTVd capacity for mutation repair, we have created 6 mutants carrying 3-4 nucleotides deletions or insertions at three unique restriction sites, EagI, StyI and AvaII. Differences in the infectivity of these in vitro generated PSTVd mutants can result from where the mutations map, as well as from the extent to which the secondary structure of the molecule is affected. Deletion or insertion of 4 nucleotides at the EagI and StyI sites led to loss of infectivity. However, mutants with deletion (PSTVd-Ava-del) or insertion (PSTVd-Ava-in) of 3 nucleotides (221GAC223), at the AvaII site (loop 20) were viable but not genetically stable. In all analyzed plants, reversion to the wild type PSTVd-S23 sequence was observed for the PSTVd-Ava-in mutant a few weeks after agroinfiltration. Analysis of PSTVd-Ava-del progeny allowed the identification of 10 new sequence variants carrying various modifications, some of them having retained the original three nucleotide deletion at the AvaII site. Interestingly, other variants gained three nucleotides in the deletion site but did not revert to the original wild type sequence. The genetic stability of the progeny PSTVd-Ava-del sequence variants was evaluated in tomato leaves (early infection) and in both leaves and roots (late infection), respectively.

Development of User-Friendly Method to Distinguish Subspecies of the Korean Medicinal Herb Perilla frutescens Using Multiplex-PCR.

Perilla (Perilla frutescens) is an economically and culturally important plant in East Asia. Plant breeding between cultivars has enhanced the genetic diversity of perilla overall, but means that functionally diverse subspecies are more difficult to identify and distinguish. In this study, we developed gene-based DNA markers to distinguish between the Korean herbal medicinal perilla varieties. We identified informative simple sequence repeat (SSR) regions on the promoter regions of the Myb-P1 and dihydroflavonol 4-reductase (DFR) genes, as well as a large insertion-deletion (indel) region in the limonene synthase (LS) gene, and developed markers to characterize the distinct subspecies differences (PfMyb-P1pro, PfDFRpro, and PfLS, respectively). Using the PfLS primers, a 430-bp region could be amplified from P. frutescens var. acuta, crispa, and f. viridis (known as Jasoyeop, Jureum-soyeop, and Chungsoyeop, respectively), but not from P. frutescens var. japonica (Dlggae). The PfMybpro primers resulted in PCR products of 314 or 316, 330, 322, and 315 bp from Dlggae, Jasoyeop, Jureum-soyeop, and Chungsoyeop, respectively, and the PfDFRpro primers resulted in products of 189 or 202, 187 or 189, 185 or 189, and 193bp, respectively, for the four perilla subspecies. Combining these three reactions into a single multiplex PCR approach resulted in subspecies-specific PCR band patterns for six common types of commercial perilla, distinguishing between three varieties of Dlggae (Cham-Dlggae, Ip-Dlggae, and Bora-Dlggae), as well as identifying Jasoyeop, Jureum-soyeop, and Chungsoyeop. These user-friendly markers will be valuable as a simple and efficient method for identifying the Korean medicinal herb Jasoyeop, as well as distinguishing between other functionally distinct subspecies, which may have broad applications in the Korean herbal industry.

Complete chloroplast genome sequences of Solanum commersonii and its application to chloroplast genotype in somatic hybrids with Solanum tuberosum.

KEY MESSAGE: Chloroplast genome of Solanum commersonii and S olanum tuberosum were completely sequenced, and Indel markers were successfully applied to distinguish chlorotypes demonstrating the chloroplast genome was randomly distributed during protoplast fusion. Somatic hybridization has been widely employed for the introgression of resistance to several diseases from wild Solanum species to overcome sexual barriers in potato breeding. Solanum commersonii is a major resource used as a parent line in somatic hybridization to improve bacterial wilt resistance in interspecies transfer to cultivated potato (S. tuberosum). Here, we sequenced the complete chloroplast genomes of Lz3.2 (S. commersonii) and S. tuberosum (PT56), which were used to develop fusion products, then compared them with those of five members of the Solanaceae family, S. tuberosum, Capsicum annum, S. lycopersicum, S. bulbocastanum and S. nigrum and Coffea arabica as an out-group. We then developed Indel markers for application in chloroplast genotyping. The complete chloroplast genome of Lz3.2 is composed of 155,525 bp, which is larger than the PT56 genome with 155,296 bp. Gene content, order and orientation of the S. commersonii chloroplast genome were highly conserved with those of other Solanaceae species, and the phylogenetic tree revealed that S. commersonii is located within the same node of S. tuberosum. However, sequence alignment revealed nine Indels between S. commersonii and S. tuberosum in their chloroplast genomes, allowing two Indel markers to be developed. The markers could distinguish the two species and were successfully applied to chloroplast genotyping (chlorotype) in somatic hybrids and their progenies. The results obtained in this study confirmed the random distribution of the chloroplast genome during protoplast fusion and its maternal inheritance and can be applied to select proper plastid genotypes in potato breeding program.

Resequencing at ≥40-Fold Depth of the Parental Genomes of a Solanum lycopersicum × S. pimpinellifolium Recombinant Inbred Line Population and Characterization of Frame-Shift InDels That Are Highly Likely to Perturb Protein Function.

A recombinant in-bred line population derived from a cross between Solanum lycopersicum var. cerasiforme (E9) and S. pimpinellifolium (L5) has been used extensively to discover quantitative trait loci (QTL), including those that act via rootstock genotype, however, high-resolution single-nucleotide polymorphism genotyping data for this population are not yet publically available. Next-generation resequencing of parental lines allows the vast majority of polymorphisms to be characterized and used to progress from QTL to causative gene. We sequenced E9 and L5 genomes to 40- and 44-fold depth, respectively, and reads were mapped to the reference Heinz 1706 genome. In L5 there were three clear regions on chromosome 1, chromosome 4, and chromosome 8 with increased rates of polymorphism. Two other regions were highly polymorphic when we compared Heinz 1706 with both E9 and L5 on chromosome 1 and chromosome 10, suggesting that the reference sequence contains a divergent introgression in these locations. We also identified a region on chromosome 4 consistent with an introgression from S. pimpinellifolium into Heinz 1706. A large dataset of polymorphisms for the use in fine-mapping QTL in a specific tomato recombinant in-bred line population was created, including a high density of InDels validated as simple size-based polymerase chain reaction markers. By careful filtering and interpreting the SnpEff prediction tool, we have created a list of genes that are predicted to have highly perturbed protein functions in the E9 and L5 parental lines.