Genetic diversity and population structure of Perilla frutescens collected from Korea and China based on simple sequence repeats (SSRs)

Insights into structure, codon usage, repeats, and RNA editing of the complete mitochondrial genome of Perilla frutescens (Lamiaceae)

Perilla frutescens (L.) Britton, a member of the Lamiaceae family, stands out as a versatile plant highly valued for its unique aroma and medicinal properties. Additionally, P. frutescens seeds are rich in Îś-linolenic acid, holding substantial economic importance. While the nuclear and chloroplast genomes of P. frutescens have already been documented, the complete mitochondrial genome sequence remains unreported. To this end, the sequencing, annotation, and assembly of the entire Mitochondrial genome of P. frutescens were hereby conducted using a combination of Illumina and PacBio data. The assembled P. frutescens mitochondrial genome spanned 299,551 bp and exhibited a typical circular structure, involving a GC content of 45.23%. Within the genome, a total of 59 unique genes were identified, encompassing 37 protein-coding genes, 20 tRNA genes, and 2 rRNA genes. Additionally, 18 introns were observed in 8 protein-coding genes. Notably, the codons of the P. frutescens mitochondrial genome displayed a notable A/T bias. The analysis also revealed 293 dispersed repeat sequences, 77 simple sequence repeats (SSRs), and 6 tandem repeat sequences. Moreover, RNA editing sites preferentially produced leucine at amino acid editing sites. Furthermore, 70 sequence fragments (12,680 bp) having been transferred from the chloroplast to the mitochondrial genome were identified, accounting for 4.23% of the entire mitochondrial genome. Phylogenetic analysis indicated that among Lamiaceae plants, P. frutescens is most closely related to Salvia miltiorrhiza and Platostoma chinense. Meanwhile, inter-species Ka/Ks results suggested that Ka/Ks < 1 for 28 PCGs, indicating that these genes were evolving under purifying selection. Overall, this study enriches the mitochondrial genome data for P. frutescens and forges a theoretical foundation for future molecular breeding research.

Integrated analysis of miRNA, transcriptome, and degradome sequencing provides new insights into lipid metabolism in perilla seed

MicroRNAs (miRNA) are small noncoding RNAs that play a crucial as molecular regulators in lipid metabolism in various oil crops. Perilla (Perilla frutescens) is a specific oil crop known for its high alpha-linolenic acid (C18:3n3, ALA) content (>65 %) in their seed oils. In view of the regulatory mechanism of miRNAs in perilla remains unclear, we conducted miRNAs and transcriptome sequencing in two cultivars with distinct lipid compositions. A total of 525 unique miRNAs, including 142 differentially expressed miRNAs was identified in perilla seeds. The 318 miRNAs targeted 7,761 genes. Furthermore, we identified 112 regulated miRNAs and their 610 target genes involved in lipid metabolism. MiR159b and miR167a as the core nodes to regulate the expression of genes in oil biosynthesis (e.g., KAS, FATB, GPAT, FAD, DGK, LPAAT) and key regulatory TFs (e.g., MYB, ARF, DOF, SPL, NAC, TCP, and bHLH). The 1,219 miRNA-mRNA regulation modules were confirmed through degradome sequencing. Notably, pf-miR159b-MYBs and pf-miR167a-ARFs regulation modules were confirmed. They exhibited significantly different expression levels in two cultivars and believed to play important roles in oil biosynthesis in perilla seeds. This provides valuable insights into the functional analysis of miRNA-regulated lipid metabolism in perilla seeds.

Characterization of a Group of 2,3-Oxidosqualene Cyclase Genes Involved in the Biosynthesis of Diverse Triterpenoids of Perilla frutescens

Perilla frutescens (L.), a traditional edible and medicinal crop, contains diverse triterpenes with multiple pharmacological properties. However, the biosynthesis of triterpenes in perilla remains rarely revelation. In this study, nine putative 2,3-oxidosqualene cyclase (OSC) genes (PfOSC1-9) were screened from the P. frutescens genome and functionally characterized by heterologous expression. Camelliol C, a triterpenol with pharmacological effect, was first identified as abundant in perilla seeds, and the camelliol C synthase (PfOSC7) was first identified in P. frutescens utilizing a yeast system. In addition, PfOSC2, PfOSC4, and PfOSC9 were identified as cycloartenol, lupeol, and β-amyrin synthase, respectively. Molecular docking and site-directed mutagenesis revealed that changes in Leu253 of PfOSC4, Ala480 of PfOSC7, and Trp257 of PfOSC9 might lead to variations of catalytic specificity or efficiency. These results will provide key insights into the biosynthetic pathways of triterpenoids and have great significance for germplasm breeding in P. frutescens.

Utilization of Novel Perilla SSR Markers to Assess the Genetic Diversity of Native Perilla Germplasm Accessions Collected from South Korea

The Perilla crop is highly regarded in South Korea, both as a health food and traditional food. However, there is still a lack of Perilla SSR primer sets (PSPSs) for studying genetic variation among accessions of cultivated and weedy types of Perilla crop (CWTPC) from South Korea. In this study, 30 PSPSs were newly developed based on transcriptome contigs in P. frutescens, and 17 of these PSPSs were used to study the genetic diversity, phylogenetic relationships and structure population among 90 accessions of the CWTPC collected from South Korea. A total of 100 alleles were detected from selected 17 PSPSs, with an average of 5.9 alleles per locus. The gene diversity (GD) ranged from 0.164 to 0.831, with an average of 0.549. The average GD values from the cultivated var. frutescens, weedy var. frutescens, and weedy var. crispa, were 0.331, 0.588, and 0.389 respectively. In addition, most variance shown by Perilla SSR markers was within a population (73%). An analysis of the population structure and phylogenetic relationships showed that the genetic relationship among accessions of the weedy var. frutescens and weedy var. crispa is closer than that for the accessions of the cultivated var. frutescens. Based on association analysis between 17 PSPSs and three seed traits in 90 Perilla accessions, we detected 11 PSPSs that together were associated with the seed size and seed hardness traits. Therefore, the newly developed PSPSs will be useful for analyzing genetic variation among accessions of the CWTPC, association mapping, and selection of important morphological traits in Perilla crop breeding programs.

Genetic variation of seed oil characteristics in native Korean germplasm of Perilla crop (Perilla frutescens L.) using SSR markers

BACKGROUND

In order to maximize the use of valuable native Perilla germplasm in South Korea, knowledge of the Perilla seed oil content and genetic variation among native Perilla germplasm resources is very important for the conservation and development of new Perilla seed oil varieties using the native Perilla germplasm accessions preserved from the Rural Development Administration Genebank (RDA-Genebank) collection from South Korea.

OBJECTIVES

In this study, we studied population structure and association mapping to identify Perilla SSR markers (PSMs) associated with the five fatty acid contents and two seed characteristics of the native Korean Perilla germplasm accessions of cultivated var. frutescens of the RDA-Genebank collected in South Korea.

METHODS

For an association mapping analysis to find PSMs associated with the five fatty acid contents and two seed characteristics of the Perilla germplasm accessions of cultivated var. frutescens, we evaluated the content of five fatty acids of 280 native Korean Perilla germplasm accessions and used 29 Perilla SSR primer sets to measure the genetic diversity and relationships, population structure, and association mapping of the native Korean Perilla germplasm accessions of the RDA-Genebank collected in South Korea.

RESULTS

Five fatty acids of 280 native Korean Perilla accessions were identified as follows: palmitic acid (PA) (5.30-8.66%), stearic acid (SA) (1.60-4.19%), oleic acid (OA) (9.60-22.5%), linoleic acid (LA) (8.38-25.4%), and linolenic acid (LNA) (52.7-76.4%). In a correlation analysis among the five fatty acids and two seed characteristics of the 280 Perilla accessions, the combinations of PA and SA (0.794**) and SA and OA (0.724**) showed a particularly high positive correlation coefficients compare to other combinations. By using an association analysis of the 29 PSMs and the five fatty acids in the 280 Perilla accessions, we found 17 PSMs (KNUPF1, KNUPF2, KNUPF4, KNUPF10, KNUPF16, KNUPF25, KNUPF26, KNUPF28, KNUPF37, KNUPF55, KNUPF62, KNUPF71, KNUPF74, KNUPF77, KNUPF85, KNUPF89, and KNUPF118) associated with the content of the five fatty acid components and two seed characteristics.

CONCLUSIONS

These PSMs are considered to be useful molecular markers related to five fatty acid components and two seed characteristics for selecting accessions from the germplasm accessions of the Perilla crop and their related weedy types through association mapping analysis and marker-assisted selection (MAS) breeding programs.

Genetic Diversity Study on Geographical Populations of the Multipurpose Species Elsholtzia stauntonii Using Transferable Microsatellite Markers

Elsholtzia stauntonii Benth. (Lamiaceae) is an economically important ornamental, medicinal and aromatic plant species. To meet the increasing market demand for E. stauntonii, it is necessary to assess genetic diversity within the species to accelerate the process of genetic improvement. Analysis of the transferability of simple sequence repeat (SSR) markers from related species or genera is a fast and economical method to evaluate diversity, and can ensure the availability of molecular markers in crops with limited genomic resources. In this study, the cross-genera transferability of 497 SSR markers selected from other members of the Lamiaceae (Salvia L., Perilla L., Mentha L., Hyptis Jacq., Leonurus L., Pogostemon Desf., Rosmarinus L., and Scutella L.) to E. stauntonii was 9.05% (45 primers). Among the 45 transferable markers, 10 markers revealed relatively high polymorphism in E. stauntonii. The genetic variation among 825 individuals from 18 natural populations of E. stauntonii in Hebei Province of China was analyzed using the 10 polymorphic SSR markers. On the basis of the SSR data, the average number of alleles (N _A), expected heterozygosity (H _E), and Shannon's information index (I) of the 10 primers pairs were 7.000, 0.478, and 0.688, respectively. Lower gene flow (N _m = 1.252) and high genetic differentiation (F _st = 0.181) were detected in the populations. Analysis of molecular variance (AMOVA) revealed that most of the variation (81.47%) was within the populations. Integrating the results of STRUCTURE, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) clustering, and principal coordinate analysis, the 825 samples were grouped into two clusters associated with geographical provenance (southwestern and northeastern regions), which was consistent with the results of a Mantel test (r = 0.56, p < 0.001). Overall, SSR markers developed in related genera were effective to study the genetic structure and genetic diversity in geographical populations of E. stauntonii. The results provide a theoretical basis for conservation of genetic resources, genetic improvement, and construction of a core collection for E. stauntonii.

Construction of a core collection of native Perilla germplasm collected from South Korea based on SSR markers and morphological characteristics

The leaves and seed oil of Perilla crop (Perilla frutescens L.) have attracted interest as health foods in East Asia. This crop has been traditionally cultivated and used for a long time as a folk plant, especially in Korea. In our study, the 22 SSR markers and eight morphological traits were used to assess the genetic diversity and population structure, to select a core collection of 400 Perilla accessions conserved in the RDA-Genebank of South Korea. A total of 173 alleles were detected and the number of alleles per locus ranged from 4 to 15 (average = 7.9). Gene diversity and polymorphic information content ranged from 0.138 to 0.868 (average = 0.567) and 0.134 to 0.853 (average = 0.522), respectively. The 400 accessions were not clearly distinguished geographically by STRUCTURE and UPGMA analyses. A core collection (44 accessions) was selected from the entire collection by using PowerCore. The core collection accounted for 11.0% of the entire Perilla collection, including 100% of the number of alleles maintained in the whole collection and with similar or greater Shannon-Weaver and Nei diversity indices than the whole collection. The core collection selected by SSR markers was evenly distributed in three clusters on a scatter plot by eight morphological traits. The first core collection of Perilla accessions was constructed, and it maintained allelic richness. Further modification of the core collection is expected with the continuous addition of new accessions of the two cultivated types of Perilla crop and their weedy types.

Identifying SSR Markers Related to Seed Fatty Acid Content in Perilla Crop (Perilla frutescens L.)

Perilla seed oil has been attracting attention in South Korea as a health food. Five fatty acids of 100 Perilla accessions were identified as follows: palmitic acid (PA) (5.10-9.13%), stearic acid (SA) (1.70-3.99%), oleic acid (OA) (11.1-21.9%), linoleic acid (LA) (10.2-23.4%), and linolenic acid (LNA) (54.3-75.4%). Additionally, the 100 Perilla accessions were divided into two groups (high or low) based on the total fatty acid content (TFAC). By using an association analysis of 40 simple sequence repeat (SSR) markers and the six Perilla seed oil traits in the 100 Perilla accessions, we detected four SSR markers associated with TFAC, five SSR markers associated with LNA, one SSR marker associated with LA, two SSR markers each associated with OA and PA, and four SSR markers associated with SA. Among these SSR markers, four SSR markers (KNUPF14, KNUPF62, KNUPF72, KNUPF85) were all associated with TFAC and LNA. Moreover, two SSR markers (KNUPF62, KNUPF85) were both associated with TFAC, LNA, and OA. Therefore, these SSR markers are considered to be useful molecular markers for selecting useful accessions related to fatty acid contents in Perilla germplasm and for improving the seed oil quality of Perilla crop through marker-assisted selection (MAS) breeding programs.

Identifying SSR markers associated with seed characteristics in Perilla (Perilla frutescens L.)

Substantial differences exist in seed dormancy between cultivated crops and their wild progenitors. The purpose of this study was to identify simple sequence repeat (SSR) markers associated with seed characteristics in cultivated and weedy types of Perilla crop. By using an association analysis of 29 SSR markers and three seed traits in 38 Perilla accessions, we detected six SSR markers associated with the seed germination rate (SGR), eight SSR markers associated with seed hardness (SH), and seven SSR markers associated with seed size (SS). Among these SSR markers, three (KNUPF3, KNUPF25, KNUPF60) were associated with the SGR, SH, and SS traits. Correlation analysis among the three seed traits of the 38 Perilla accessions showed a positive correlation coefficient for the combination of SGR and SS (0.811**) and a negative correlation coefficient for the combinations of SGR and SH (- 0.706**), and SS and SH (- 0.899**). A phylogenetic tree constructed using the unweighted pair group method with arithmetic mean (UPGMA) revealed that accessions of cultivated P. frutescens var. frutescens could be distinguished from weedy accessions of P. frutescens var. frutescens and P. frutescens var. crispa using the 29 SSR markers. Selected SSR markers related to the three seed traits distinguished accessions of cultivated and weedy types. Therefore, these results are very important for understanding the seed characteristics of cultivated and weedy types of Perilla crop. It will further help for improving the seed quality of Perilla crop through marker-assisted selection (MAS) breeding programs.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s12298-021-00933-3) contains supplementary material, which is available to authorized users.

Assessment of genetic diversity and population structure among a collection of Korean Perilla germplasms based on SSR markers

BACKGROUND

Information on the genetic variation of genetic resource collections is very important for both the conservation and utilization of crop germplasms in genebanks. Var. frutescens of Perilla crop is extensively cultivated in South Korea as both an oil crop and a vegetable crop.

OBJECTIVES

We used SSR markers to evaluate the genetic diversity, genetic relationships, and population structure of 155 accessions of var. frutescens that have been selected as genetic resources for the development of leaf vegetable cultivars and preserved in the RDA-Genebank collection from South Korea.

METHODS

A total of 155 accessions of var. frutescens of Perilla crop collected in South Korea were obtained from the RDA-Genebank of the Republic of Korea. We selected 20 SSR markers representing the polymorphism of and adequately amplifying all the Perilla accessions.

RESULTS

The average GD and PIC values were 0.642 and 0.592, respectively, with ranges of 0.244-0.935 and 0.232- 0.931. The genetic variability in the southern region of South Korea was higher than that in the central region. The clustering patterns were not clearly distinguished between the accessions of var. frutescens from the central and southern regions of South Korea.

CONCLUSION

These results regarding the genetic diversity and population structure of the 155 accessions of var. frutescens of South Korea provide useful information for understanding the genetic variability of this crop and selecting and managing core germplasm sets in the RDA-Genebank of the Republic of Korea.

Genome-Wide Characterization and Comparative Analyses of Simple Sequence Repeats among Four Miniature Pig Breeds

Simple Summary

Simple sequence repeats (SSRs) are present at high densities in regulatory elements, suggesting that they may affect gene function and phenotypic traits. Therefore, SSRs can be exploited in marker-assisted selection. In addition, they can be widely used as molecular markers to study genetic diversity, population structure, and evolution. While SSRs have been widely studied in many mammalian species, very little research has focused on genome-wide SSRs of miniature pigs, a small but special group of pigs that express the dwarf phenotype. Based on the SSR-enriched library building and sequencing, about 30,000 novel polymorphic SSRs for four miniature pig breeds were mapped to the Duroc pig reference genome. The four miniature pig breeds had different numbers and types of SSRs and distributions of repeat units. There were 2518 polymorphic SSRs in the intron or exon regions that were common to all four breeds and functional analyses revealed 17 genes that were associated with body size and other genes that were associated with growth and development. In conclusion, the SSRs detected in the miniature pigs in this study may provide useful genetic markers for the selection of farm animals and the polymorphic SSRs provide valuable insights into the determination of mature body size, as well as the immunity, growth and development of animals.

Abstract

Simple sequence repeats (SSRs) are commonly used as molecular markers in research on genetic diversity and discrimination among taxa or breeds because polymorphisms in these regions contribute to gene function and phenotypically important traits. In this study, we investigated genome-wide characteristics, repeat units, and polymorphisms of SSRs using sequencing data from SSR-enriched libraries created from Wuzhishan (WZS), Bama (BM), inbred Luchuan (LC) and Zangxiang (ZX) miniature pig breeds. The numbers and types of SSRs, distributions of repeat units and polymorphic SSRs varied among the four breeds. Compared to the Duroc pig reference genome, 2518 polymorphic SSRs were unique and common to all four breeds and functional annotation revealed that they may affect the coding and regulatory regions of genes. Several examples, such as FGF23, MYF6, IGF1R, and LEPROT, are associated with growth and development in pigs. Three of the polymorphic SSRs were selected to confirm the polymorphism and the corresponding alleles through fluorescence polymerase chain reaction (PCR) and capillary electrophoresis. Together, this study provides useful insights into the discovery, characteristics and distribution of SSRs in four pig breeds. The polymorphic SSRs, especially those common and unique to all four pig breeds, might affect associated genes and play important roles in growth and development.