Genetic diversity and population structure of Camellia huana (Theaceae), a limestone species with narrow geographic range, based on chloroplast DNA sequence and microsatellite markers

Genetic diversity evaluation of Luculia yunnanensis, a vulnerable species endemic to Yunnan, Southwestern China based on morphological traits and EST-SSR markers

Luculia yunnanensis is a vulnerable species endemic to Yunnan Province, Southwestern China, which has high ornamental value. Its wild population has not been fully protected and utilized for a long time, which is not conducive to the long-term stable development of this species. Genetic diversity assessment is the basis and prerequisite for the conservation of rare species. In this study, 21 phenotypic traits and 17 highly polymorphic EST-SSR markers were used to analyze the genetic diversity and genetic structure of 164 individuals from six L. yunnanensis populations. The coefficient of variation of 21 phenotypic traits ranged from 11.76% to 52.58% (mean=21.72%), and the coefficient of variation of 18 traits was less than 30%. The average values of Ne, I, Ho and He were 1.710, 0.619, 0.384, and 0.352, respectively. The genetic diversity of LLO (Ho = 0.476 and He = 0.426) and LCM (Ho = 0.424 and He = 0.381) populations in Lushui County was highest. The GDX populations (Ho = 0.335 and He = 0.269) isolated by Gaoligong Mountain had the lowest genetic diversity. The AMOVA results showed that 13.04% of the genetic variation was among populations and 86.96% was within populations. The average phenotypic differentiation coefficient of phenotypic traits among populations was 18.69%. The results of phenotypic and genetic variation analysis were consistent, indicating that the most of variation exists within population. Genetic structure, UPGMA clustering and PCA analysis results showed that the populations of L. yunnanensis had obvious geographical divisions, and the populations distributed in the southern region and distributed in the northern region of the Nujiang River clustered into one group respectively. Combining the results of phenotypic and molecular markers, we recommend that give priority to the protection of LLO, LCM and GDX population, in order to ensure the sustainable utilization of L. yunnanensis germplasm resources.

Comparative analysis of chloroplast genomes of Pulsatilla species reveals evolutionary and taxonomic status of newly discovered endangered species Pulsatilla saxatilis

BACKGROUND

Pulsatilla saxatilis, a new species of the genus Pulsatilla has been discovered. The morphological information of this species has been well described, but its chloroplast genome characteristics and comparison with species of the same genus remain to be reported.

RESULTS

Our results showed that the total length of chloroplast (cp.) genome of P. saxatilis is 162,659 bp, with a GC content of 37.5%. The cp. genome contains 134 genes, including 90 known protein-coding genes, 36 tRNA genes, and 8 rRNA genes. P. saxatilis demonstrated similar characteristics to other species of genus Pulsatilla. Herein, we compared cp. genomes of 10 species, including P. saxatilis, and found that the cp. genomes of the genus Pulsatilla are extremely similar, with a length of 162,322-163,851 bp. Furthermore, The SSRs of Pulsatilla ranged from 10 to 22 bp in length. Among the four structural regions of the cp. genome, most long repeats and SSRs were detected in the LSC region, followed by that in the SSC region, and least in IRA/ IRB regions. The most common types of long repeats were forward and palindromic repeats, followed by reverse repeats, and only a few complementary repeats were found in 10 cp. genomes. We also analyzed nucleotide diversity and identified ccsA_ndhD, rps16_trnK-UUU, ccsA, and rbcL, which could be used as potential molecular markers for identification of Pulsatilla species. The results of the phylogenetic tree constructed by connecting the sequences of high variation regions were consistent with those of the cp. gene phylogenetic tree, and the species more closely related to P. saxatilis was identified as the P. campanella.

CONCLUSION

It was determined that the closest species to P. saxatilis is P. campanella, which is the same as the conclusion based on pollen grain characteristics, but different from the P. chinensis determined based on morphological characteristics. By revealing information on the chloroplast characteristics, development, and evolution of the cp. genome and the potential molecular markers, this study provides effective molecular data regarding the evolution, genetic diversity, and species identification of the genus Pulsatilla.

Population Structure and Genetic Diversity of the Toona ciliata (Meliaceae) Complex Assayed with Chloroplast DNA Markers

Toona ciliata is a deciduous or semi-deciduous tree species and belongs to the Toona genus of the Meliaceae family. Owing to low natural regeneration and over-exploitation, the species is listed as an endangered species at level II in China and its conservation has received increasing concern. Here, we sampled 447 individuals from 29 populations across the range-wide distribution of the T. ciliata complex in China and assessed their genetic variation using two chloroplast DNA markers. The results showed that the overall haplotype diversity and nucleotide diversity per site were high at h = 0.9767 and π = 0.0303 for the psbA-trnH fragment and h= 0.8999 and π = 0.0189 for the trnL-trnL fragment. Phylogenetic analysis supported the division of the natural distribution of T. ciliata complex into western and eastern regions. The genetic diversity was higher in the western region than in the eastern region, showing significant phylogeographic structure. Genetic differentiation among populations was moderate (Φst=42.87%), and the effects of isolation by distance (IBD) were significant. A neutrality test and mismatch distribution analysis indicated that the distribution of the T. ciliata complex generally did not expand, although a few local populations could likely expand after bottleneck effects. The overall results were complementary to and consolidated previous studies using mitochondrial and nuclear DNA markers. We finally discussed strategies for the genetic conservation of the T. ciliata complex.

Whole chloroplast genome-specific non-synonymous SNPs reveal the presence of substantial diversity in the pigeonpea mini-core collection

To unravel the plastid genome diversity among the cultivated groups of the pigeonpea germplasm, we characterized the SNP occurrence and distribution of 142 pigeonpea mini-core collections based on their reference-based assembly of the chloroplast genome. A total of 8921 SNPs were found, which were again filtered and finally 3871 non-synonymous SNPs were detected and used for diversity estimates. These 3871 SNPs were classified into 12 groups and were present in only 44 of the 125 genes, demonstrating the presence of a precise mechanism for maintaining the whole chloroplast genome throughout evolution. The Acetyl-CoA carboxylase D gene possesses the maximum number of SNPs (12.29%), but the Adenosine Tri-Phosphate synthatase cluster genes (atpA, atpB, atpE, atpF, atpH, and atpI) altogether bear 43.34% of the SNPs making them most diverse. Various diversity estimates, such as the number of effective alleles (1.013), Watterson's estimate (0.19), Tajima's D ( - 3.15), Shannon's information index (0.036), suggest the presence of less diversity in the cultivated gene pool of chloroplast genomes. The genetic relatedness estimates based on pairwise correlations were also in congruence with these diversity descriptors and indicate the prevalence of rare alleles in the accessions. Interestingly, no stratification was observed either through STRUCTURE, PCoA, or phylogenetic analysis, indicating the common origin of the chloroplast in all the accessions used, irrespective of their geographical distribution. Further 6194 Cleaved Amplified Polymorphic Sequences (CAPS) markers for 531 SNPs were developed and validated in a selected set of germplasm. Based on these results, we inferred that all of the cultivated gene pools of pigeonpea have a common origin for the chloroplast genome and they possess less diversity in protein-coding regions, indicating a stable and evolved plastid genome. At the same time, all diversity analysis indicates the occurrence of rare alleles, suggesting the suitability of the mini-core collection in future pigeonpea improvement programs. In addition, the development of chloroplast genome-based CAPS markers would have utility in pigeonpea breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03785-8.

Plastid phylogenomics and species discrimination in the "Chinese" clade of Roscoea (Zingiberaceae)

Roscoea is an alpine or subalpine genus from the pan-tropical family Zingiberaceae, which consists of two disjunct groups in geography, namely the "Chinese" clade and the "Himalayan" clade. Despite extensive research on the genus, Roscoea species remain poorly defined and relationships between these species are not well resolved. In this study, we used plastid genomes of nine species and one variety to resolve phylogenetic relationships within the "Chinese" clade of Roscoea and as DNA super barcodes for species discrimination. We found that Roscoea plastid genomes ranged in length from 163,063 to 163,796 bp, and encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, four rRNA genes. In addition, expansion and contraction of the IR regions showed obvious infraspecific conservatism and interspecific differentiation. Plastid phylogenomics revealed that species belonging to the "Chinese" clade of Roscoea can be divided into four distinct subclades. Furthermore, our analysis supported the independence of R. cautleoides var. pubescens, the recovery of Roscoea pubescens Z.Y. Zhu, and a close relationship between R. humeana and R. cautloides. When we used the plastid genome as a super barcode, we found that it possessed strong discriminatory power (90%) with high support values. Intergenic regions provided similar resolution, which was much better than that of protein-coding regions, hypervariable regions, and DNA universal barcodes. However, plastid genomes could not completely resolve Roscoea phylogeny or definitively discriminate species. These limitations are likely related to the complex history of Roscoea speciation, poorly defined species within the genus, and the maternal inheritance of plastid genomes.

Genetic diversity, population genetic structure and gene flow in the rare and endangered wild plant Cypripedium macranthos revealed by genotyping-by-sequencing

BACKGROUND

Genetic diversity, genetic structure, and gene flow in plant populations and their influencing factors are important in conservation biology. Cypripedium macranthos is one of the few wild orchids with high ornamental value in northern China. However, over the past decade, excessive collection, trading, tourism development, habitat fragmentation, deceptive pollination, and seed germination difficulties have all caused a sharp decline in the number of C. macranthos individuals and its population. In order to propose a scientific and effective conservation strategy, the genetic diversity, genetic structure and gene flow of the current CM population are urgent scientific issues to be clarified.

RESULTS

Here, 99 individuals of C. macranthos from north and northeast China were analyzed to evaluate the genetic diversity, gene flow among populations, and genetic structure by genotyping-by-sequencing. More than 68.44 Gb high-quality clean reads and 41,154 SNPs were obtained. Our data based on bioinformatics methods revealed that C. macranthos has lower genetic diversity, high levels of historical gene flow, and moderate-to-high genetic differentiation between populations. The gene migration model revealed that the direction of gene flow was mainly from northeast populations to north populations in China. The results of genetic structure analysis showed that 11 C. macranthos populations can be considered as two groups, and further divided into four subgroups. Moreover, the Mantel test detected no significant "Isolation by Distance" between populations.

CONCLUSIONS

Our study demonstrates that the present genetic diversity and genetic structure of C. macranthos populations were mainly caused by biological characteristics, human interference, habitat fragmentation, and restricted gene flow. Finally, constructive measures, which can provide a basis for the proposal of conservation strategies, have been suggested.

Genome-wide association study of leaf-related traits in tea plant in Guizhou based on genotyping-by-sequencing

BACKGROUND

Studying the genetic characteristics of tea plant (Camellia spp.) leaf traits is essential for improving yield and quality through breeding and selection. Guizhou Plateau, an important part of the original center of tea plants, has rich genetic resources. However, few studies have explored the associations between tea plant leaf traits and single nucleotide polymorphism (SNP) markers in Guizhou.

RESULTS

In this study, we used the genotyping-by-sequencing (GBS) method to identify 100,829 SNP markers from 338 accessions of tea germplasm in Guizhou Plateau, a region with rich genetic resources. We assessed population structure based on high-quality SNPs, constructed phylogenetic relationships, and performed genome-wide association studies (GWASs). Four inferred pure groups (G-I, G-II, G-III, and G-IV) and one inferred admixture group (G-V), were identified by a population structure analysis, and verified by principal component analyses and phylogenetic analyses. Through GWAS, we identified six candidate genes associated with four leaf traits, including mature leaf size, texture, color and shape. Specifically, two candidate genes, located on chromosomes 1 and 9, were significantly associated with mature leaf size, while two genes, located on chromosomes 8 and 11, were significantly associated with mature leaf texture. Additionally, two candidate genes, located on chromosomes 1 and 2 were identified as being associated with mature leaf color and mature leaf shape, respectively. We verified the expression level of two candidate genes was verified using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and designed a derived cleaved amplified polymorphism (dCAPS) marker that co-segregated with mature leaf size, which could be used for marker-assisted selection (MAS) breeding in Camellia sinensis.

CONCLUSIONS

In the present study, by using GWAS approaches with the 338 tea accessions population in Guizhou, we revealed a list of SNPs markers and candidate genes that were significantly associated with four leaf traits. This work provides theoretical and practical basis for the genetic breeding of related traits in tea plant leaves.

Development of SSR molecular markers and genetic diversity analysis of Clematis acerifolia from Taihang Mountains

Investigating the genetic diversity and population structure is important in conserving narrowly distributed plants. In this study, 90 Clematis acerifolia (C. acerifolia) plants belonging to nine populations were collected from the Taihang Mountains in Beijing, Hebei, and Henan. Twenty-nine simple sequence repeats (SSR) markers developed based on RAD-seq data were used to analyze the genetic diversity and population structure of C. acerifolia. The mean PIC value for all markers was 0.2910, indicating all SSR markers showed a moderate degree of polymorphism. The expected heterozygosity of the whole populations was 0.3483, indicating the genetic diversity of both C. acerifolia var. elobata and C. acerifolia were low. The expected heterozygosity of C. acerifolia var. elobata (He = 0.2800) was higher than that of C. acerifolia (He = 0.2614). Genetic structure analysis and principal coordinate analysis demonstrated that C. acerifolia and C. acerifolia var. elobata showed great genetic differences. Molecular variance analysis (AMOVA) demonstrated that within-population genetic variation (68.31%) was the main contributor to the variation of the C. acerifolia populations. Conclusively, C. acerifolia var. elobata had higher genetic diversity than C. acerifolia, and there are significant genetic differences between C. acerifolia and C. acerifolia var. elobata, and small genetic variations within the C. acerifolia populations. Our results provide a scientific and rational basis for the conservation of C. acerifolia and provide a reference for the conservation of other cliff plants.

Complete Plastome of Physalis angulata var. villosa, Gene Organization, Comparative Genomics and Phylogenetic Relationships among Solanaceae

Physalis angulata var. villosa, rich in withanolides, has been used as a traditional Chinese medicine for many years. To date, few extensive molecular studies of this plant have been conducted. In the present study, the plastome of P. angulata var. villosa was sequenced, characterized and compared with that of other Physalis species, and a phylogenetic analysis was conducted in the family Solanaceae. The plastome of P. angulata var. villosa was 156,898 bp in length with a GC content of 37.52%, and exhibited a quadripartite structure typical of land plants, consisting of a large single-copy (LSC, 87,108 bp) region, a small single-copy (SSC, 18,462 bp) region and a pair of inverted repeats (IR: IRA and IRB, 25,664 bp each). The plastome contained 131 genes, of which 114 were unique and 17 were duplicated in IR regions. The genome consisted of 85 protein-coding genes, eight rRNA genes and 38 tRNA genes. A total of 38 long, repeat sequences of three types were identified in the plastome, of which forward repeats had the highest frequency. Simple sequence repeats (SSRs) analysis revealed a total of 57 SSRs, of which the T mononucleotide constituted the majority, with most of SSRs being located in the intergenic spacer regions. Comparative genomic analysis among nine Physalis species revealed that the single-copy regions were less conserved than the pair of inverted repeats, with most of the variation being found in the intergenic spacer regions rather than in the coding regions. Phylogenetic analysis indicated a close relationship between Physalis and Withania. In addition, Iochroma, Dunalia, Saracha and Eriolarynx were paraphyletic, and clustered together in the phylogenetic tree. Our study published the first sequence and assembly of the plastome of P. angulata var. villosa, reported its basic resources for evolutionary studies and provided an important tool for evaluating the phylogenetic relationship within the family Solanaceae.

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

Genetic diversity and phylogeography of the endemic species Chimonobambusa utilis growing in southwest China: Chloroplast DNA sequence and microsatellite marker analyses

Chimonobambusa utilis (Keng) Keng F is an endemic species distributed only in the Daluoshan Mountains, southwest China. Ch. utilis is popular due to its unique flavor and deliciousness and plays an important role in the industrial revolution in many counties in China. A total of 20 natural populations were sampled from the entire distribution range of Ch. utilis. In the present study, we used five EST-SSR molecular markers, three chloroplast DNA (trnH-psbA, atpF-atpH, and psbK-psbI), and one ITS molecular marker to elucidate the genetic diversity and phylogeography analyses of these Ch. utilis populations. The results exhibited that Ch. utilis populations showed lower genetic diversity than other angiosperms (H_T = 0.752, H_S = 0.364, and F_ST = 0.05021 for EST-SSR; H_T = 0.956, Hs = 0.507, and F_ST = 0.70121 for cpDNA; H_T = 0.868, Hs = 0.495, and F_ST = 0.70121 for nrDNA). A total of 40 alleles were detected for five polymorphic loci. We detected 20 polymorphic sites and 11 haplotypes within 1,398 bp of cpDNA and 59 polymorphic sites and 32 haplotypes within the 589 bp of the ITS sequence. Based on the haplotype distribution, we infer that there were at least two glacial refuges of Ch. utilis populations during the Quaternary Ice Age. The genetic and geographic distance were correlated (p < 0.05), indicating that narrow distribution might be the primary cause of the low genetic differentiation of Ch. utilis populations. Based on the genetic diversity of Ch. utilis populations, we recommend implementing effective genetic resource management and sustainable utilization.

Genome Survey and SSR Analysis of Camellia nitidissima Chi (Theaceae)

Camellia nitidissima Chi (CNC), a species of golden Camellia, is well known as "the queen of camellias." It is an ornamental, medicinal, and edible plant grown in China. In this study, we conducted a genome survey sequencing analysis and simple sequence repeat (SSR) identification of CNC using the Illumina sequencing platform. The 21-mer analysis predicted its genome size to be 2,778.82 Mb, with heterozygosity and repetition rates of 1.42% and 65.27%, respectively. The CNC genome sequences were assembled into 9,399,197 scaffolds, covering ∼2,910 Mb and an N50 of 869 base pair. Its genomic characteristics were found to be similar to those of Camellia oleifera. In addition, 1,940,616 SSRs were identified from the genome data, including mono-(61.85%), di-(28.71%), tri-(6.51%), tetra-(1.85%), penta-(0.57%), and hexanucleotide motifs (0.51%). We believe these data will provide a useful foundation for the development of novel molecular markers for CNC as well as for further whole-genome sequencing of CNC.

Population Genetic Structure Analysis Reveals Significant Genetic Differentiation of the Endemic Species Camellia chekiangoleosa Hu. with a Narrow Geographic Range

In order to protect and utilize the germplasm resource better, it is highly necessary to carry out a study on the genetic diversity of Camellia chekiangoleosa Hu. However, systematic research on population genetics analysis of the species is comparatively rare. Herein, 16 highly variable simple sequence repeat (SSR) markers were used for genetic structure assessment in 12 natural C. chekiangoleosa populations. The genetic diversity of C. chekiangoleosa was low (h = 0.596), within which, central populations (such as Damaoshan (DMS), Sanqingshan (SQS), and Gutianshan (GTS)) at the junction of four main mountain ranges presented high diversity and represented the center of the C. chekiangoleosa diversity distribution; the Hengshan (HS) population in the west showed the lowest diversity, and the diversity of the eastern and coastal populations was intermediate. C. chekiangoleosa exhibited a high level of genetic differentiation, and the variation among populations accounted for approximately 24% of the total variation. The major reasons for this situation are the small population scale and bottleneck effects in some populations (HS and Lingshan (LS)), coupled with inbreeding within the population and low gene flow among populations (Nm = 0.796). To scientifically protect the genetic diversity of C. chekiangoleosa, in situ conservation measures should be implemented for high-diversity populations, while low-diversity populations should be restored by reintroduction.

Chloroplast and Nuclear Genetic Diversity Explain the Limited Distribution of Endangered and Endemic Thuja sutchuenensis in China

Narrow-ranged species face challenges from natural disasters and human activities, and to address why species distributes only in a limited region is of great significance. Here we investigated the genetic diversity, gene flow, and genetic differentiation in six wild and three cultivated populations of Thuja sutchuenensis, a species that survive only in the Daba mountain chain, using chloroplast simple sequence repeats (cpSSR) and nuclear restriction site-associated DNA sequencing (nRAD-seq). Wild T. sutchuenensis populations were from a common ancestral population at 203 ka, indicating they reached the Daba mountain chain before the start of population contraction at the Last Interglacial (LIG, ∼120-140 ka). T. sutchuenensis populations showed relatively high chloroplast but low nuclear genetic diversity. The genetic differentiation of nRAD-seq in any pairwise comparisons were low, while the cpSSR genetic differentiation values varied with pairwise comparisons of populations. High gene flow and low genetic differentiation resulted in a weak isolation-by-distance effect. The genetic diversity and differentiation of T. sutchuenensis explained its survival in the Daba mountain chain, while its narrow ecological niche from the relatively isolated and unique environment in the "refugia" limited its distribution.

Genetic relationships and low diversity among the tea-oil Camellia species in Sect. Oleifera, a bulk woody oil crop in China

Tea-oil Camellia is one of the four woody oil crops in the world and has high ecological, economic and medicinal values. However, there are great differences in the classification and merging of tea-oil Camellia Sect. Oleifera species, which brings difficulties to the innovative utilization and production of tea-oil Camellia resources. Here, ISSR, SRAP and chloroplast sequence markers were analyzed in 18 populations of tea-oil Camellia Sect. Oleifera species to explore their phylogenetic relationships and genetic diversity. The results showed that their genetic diversity were low, with mean H and π values of 0.16 and 0.00140, respectively. There was high among-population genetic differentiation, with ISSR and SRAP markers showing an Fst of 0.38 and a high Nm of 1.77 and cpDNA markers showing an Fst of 0.65 and a low Nm of 0.27. The C. gauchowensis, C. vietnamensis and Hainan Island populations formed a single group, showing the closest relationships, and supported being the same species for them with the unifying name C. drupifera and classifying the resources on Hainan Island as C. drupifera. The tea-oil Camellia resources of Hainan Island should be classified as a special ecological type or variety of C. drupifera. However, cpDNA marker-based STRUCTURE analysis showed that the genetic components of the C. osmantha population formed an independent, homozygous cluster; hence, C. osmantha should be a new species in Sect. Oleifera. The C. oleifera var. monosperma and C. oleifera populations clustered into two distinct clades, and the C. oleifera var. monosperma populations formed an independent cluster, accounting for more than 99.00% of its genetic composition; however, the C. oleifera populations contained multiple different cluster components, indicating that C. oleifera var. monosperma significantly differs from C. oleifera and should be considered the independent species C. meiocarpa. Haplotype analysis revealed no rapid expansion in the tested populations, and the haplotypes of C. oleifera, C. meiocarpa and C. osmantha evolved from those of C. drupifera. Our results support the phylogenetic classification of Camellia subgenera, which is highly significant for breeding and production in tea-oil Camellia.

Revealing Distinctions in Genetic Diversity and Adaptive Evolution Between Two Varieties of Camellia sinensis by Whole-Genome Resequencing

Camellia sinensis var. sinensis (CSS) and C. sinensis var. assamica (CSA) are the two most economically important tea varieties. They have different characteristics and geographical distribution. Their genetic diversity and differentiation are unclear. Here, we identified 18,903,625 single nucleotide polymorphisms (SNPs) and 7,314,133 insertion-deletion mutations (indels) by whole-genome resequencing of 30 cultivated and three wild related species. Population structure and phylogenetic tree analyses divided the cultivated accessions into CSS and CSA containing 6,440,419 and 6,176,510 unique variations, respectively. The CSS subgroup possessed higher genetic diversity and was enriched for rare alleles. The CSA subgroup had more non-synonymous mutations and might have experienced a greater degree of balancing selection. The evolution rate (dN/dS) and KEGG enrichment indicated that genes involved in the synthesis and metabolism of flavor substances were positively selected in both CSS and CSA subpopulations. However, there are extensive genome differentiation regions (2959 bins and approximately 148 M in size) between the two subgroups. Compared with CSA (141 selected regions containing 124 genes), the CSS subgroup (830 selected regions containing 687 genes) displayed more selection regions potentially related to environmental adaptability. Fifty-three pairs of polymorphic indel markers were developed. Some markers were located in hormone-related genes with distinct alleles in the two cultivated subgroups. These identified variations and selected regions provide clues for the differentiation and adaptive evolution of tea varieties. The newly developed indel markers will be valuable in further genetic research on tea plants.