Association Mapping for Important Agronomic Traits in Safflower (Carthamus tinctorius L.) Core Collection Using Microsatellite Markers

Genetic diversity, clinical uses, and phytochemical and pharmacological properties of safflower (Carthamus tinctorius L.): an important medicinal plant

Safflower (Carthamus tinctorius L.), a member of the Asteraceae family, is widely used in traditional herbal medicine. This review summarized agronomic conditions, genetic diversity, clinical application, and phytochemicals and pharmacological properties of safflower. The genetic diversity of the plant is rich. Abundant in secondary metabolites like flavonoids, phenols, alkaloids, polysaccharides, fatty acids, polyacetylene, and other bioactive components, the medicinal plant is effective for treating cardiovascular diseases, neurodegenerative diseases, and respiratory diseases. Especially, Hydroxysafflor yellow A (HYSA) has a variety of pharmacological effects. In terms of treatment and prevention of some space sickness in space travel, safflower could be a potential therapeutic agent. Further studies are still required to support the development of safflower in medicine. Our review indicates that safflower is an important medicinal plant and research prospects regarding safflower are very broad and worthy of further investigation.

Genome-wide association studies identifies genetic loci related to fatty acid and branched-chain amino acid metabolism and histone modifications under varying nitrogen treatments in safflower (Carthamus tinctorius)

Effective identification and usage of genetic variation are prerequisites for developing nutrient-efficient cultivars. A collection of 94 safflower (Carthamus tinctorius ) genotypes (G) was investigated for important morphological and photosynthetic traits at four nitrogen (N) treatments. We found significant variation for all the studied traits except chlorophyll b (chl b ) among safflower genotypes, nitrogen treatments and G×N interaction. The examined traits showed a 2.82-50.00% increase in response to N application. Biological yield (BY) reflected a significantly positive correlation with fresh shoot weight (FSW), root length (RL), fresh root weight (FRW) and number of leaves (NOL), while a significantly positive correlation was also observed among carotenoids (C), chlorophyll a (chl a ), chl b and total chlorophyll content (CT) under all treatments. Superior genotypes with respect to plant height (PH), FSW, NOL, RL, FRW and BY were clustered into Group 3, while genotypes with better mean performance regarding chl a , chl b C and CT were clustered into Group 2 as observed in principal component analysis. The identified eight best-performing genotypes could be useful to develop improved nitrogen efficient cultivars. Genome-wide association analysis resulted in 32 marker-trait associations (MTAs) under four treatments. Markers namely DArT-45481731 , DArT-17812864 , DArT-15670279 and DArT-45482737 were found consistent. Protein-protein interaction networks of loci associated with MTAs were related to fatty acid and branched-chain amino acid metabolism and histone modifications.

Agro-Morphological and Molecular Characterization Reveal Deep Insights in Promising Genetic Diversity and Marker-Trait Associations in Fagopyrum esculentum and Fagopyrum tataricum

Characterisation of genetic diversity is critical to adequately exploit the potential of germplasm collections and identify important traits for breeding programs and sustainable crop improvement. Here, we characterised the phenotypic and genetic diversity of a global collection of the two cultivated buckwheat species Fagopyrum esculentum and Fagopyrum tataricum (190 and 51 accessions, respectively) using 37 agro-morphological traits and 24 SSR markers. A wide range of variation was observed in both species for most of the traits analysed. The two species differed significantly in most traits, with traits related to seeds and flowering contributing most to differentiation. The accessions of each species were divided into three major phenoclusters with no clear geographic clustering. At the molecular level, the polymorphic SSR markers were highly informative, with an average polymorphic information content (PIC) of over 0.65 in both species. Genetic diversity, as determined by Nei's expected heterozygosity (He), was high (He = 0.77 and He = 0.66, respectively) and differed significantly between species (p = 0.03) but was homogeneously distributed between regions, confirming the lack of genetic structure as determined by clustering approaches. The weak genetic structure revealed by the phenotypic and SSR data and the low fixation indices in both species suggested frequent seed exchange and extensive cultivation and selection. In addition, 93 and 140 significant (p < 0.05) marker-trait associations (MTAs) were identified in both species using a general linear model and a mixed linear model, most of which explained >20% of the phenotypic variation in associated traits. Core collections of 23 and 13 phenotypically and genetically diverse accessions, respectively, were developed for F. esculentum and F. tataricum. Overall, the data analysed provided deep insights into the agro-morphological and genetic diversity and genetic relationships among F. esculentum and F. tataricum accessions and pointed to future directions for genomics-based breeding programs and germplasm management.

Molecular marker development and genetic diversity exploration in Medicago polymorpha

Medicago polymorpha L. (bur clover), an invasive plant species of the genus Medicago, has been traditionally used in China as an edible vegetable crop because of its high nutritive value. However, few molecular markers for M. polymorpha have been identified. Using the recently published high-quality reference genome of M. polymorpha, we performed a specific-locus amplified fragment sequencing (SLAF-seq) analysis of 10 M. polymorpha accessions to identify molecular markers and explore genetic diversity. A total of 52,237 high-quality single nucleotide polymorphisms (SNPs) were developed. These SNPs were mostly distributed on pseudochromosome 3, least distributed on pseudochromosome 7, and relatively evenly distributed on five other pseudochromosomes of M. polymorpha. Phenotypic analysis showed that there was a great difference in phenotypic traits among different M. polymorpha accessions. Moreover, clustering all M. polymorpha accessions based on their phenotypic traits revealed three groups. Both phylogenetic analysis and principal component analysis (PCA) of all M. polymorpha accessions based on SNP markers consistently indicated that all M. polymorpha accessions could be divided into three distinct groups (I, II, and III). Subsequent genetic diversity analysis for the 10 M. polymorpha accessions validated the effectiveness of the M. polymorpha germplasm molecular markers in China. Additionally, SSR mining analysis was also performed to identify polymorphic SSR motifs, which could provide valuable candidate markers for the further breeding of M. polymorpha. Since M. polymorpha genetics have not been actively studied, the molecular markers generated from our research will be useful for further research on M. polymorpha resource utilization and marker-assisted breeding.

Occurrence of simple sequence repeats in cDNA sequences of safflower (Carthamus tinctorius) reveals the importance of SSR-containing genes for cell biology and dynamic response to environmental cues

Safflower (Carthamus tinctorius) is a diploid crop plant belonging to the family Asteraceae and is well known as one of important oilseed crops due to edible oil containing unsaturated fatty acids. In recent years it is gaining increased attention for food, pharmaceutical and industrial uses, and hence the updating its breeding methods is necessary. Genic simple sequence repeats (SSRs) in addition of being desire molecular markers, are supposed to influence gene function and the respective phenotype. This study aimed to identify SSRs in cDNA sequences and further analysis of the functional features of the SSR-containing genes to elucidate their role in biological and cellular processes. We identified 1,841 SSR regions in 1,667 cDNA sequences. Among all types of repeats, trinucleotide repeats were the most abundant (35.7%), followed by hexanucleotide (29.6%) and dinucleotide repeats (22.0%). Thirty five SSR primer pairs were validated by PCR reaction, detected a high rate of polymorphism (>57%) among safflower accessions, physically mapped on safflower genome and could clearly discriminate the cultivated accessions from wild relatives. The cDNA-derived SSR markers are suitable for evaluation of genetic diversity, linkage and association mapping studies and genome-based breeding programmes. Occurrence of SSR repeats in biologically-important classes of proteins such as kinases, transferases and transcription factors was inferred from functional analyses, which along with variability of their repeat copies, can endow the cell and whole organism the flexibility of facing with continuously changing environment, and indicate a structure-based evolution mechanism of the genome which acts as an up-to-dating tool for the cell and whole origanism, which is realized in GO terms such as involvement of most SSR-containing genes in biological, cellular and metabolic processes, especially in response to stimulus, response to stress, interaction to other organisms and defense responses.

Genomic markers analysis associated with resistance to Alternaria alternata (fr.) keissler-tomato pathotype, Solanum lycopersicum L

Alternaria alternata, the causal pathogen of early blight (EB) disease, is one of the most important diseases in tomato, and other solanaceae family. We analyzed 35 tomato genotypes for quantitative/qualitative traits and biomass growth parameters, as well as the extent and structure of genetic variation associated with EB resistance. Phenotypic comparisons displayed significant differences in leaf blade width (24.95%), stem thickness (30.28%), foliage density (18.88%), and plant size (18.89%), with significant positive correlations with EB resistance (0.18-0.75). Correlation analysis showed that mature fruit size, thickness of fruit pericarp, and leaf type were significantly and negatively correlated with EB resistance (up to -0.41). The susceptible tomato seedlings represented significant reductions in biomass parameters. According to ISSR analysis, the highest resolving power (≥0.79) and heterozygosity (≥0.24) values revealed the presence of high genetic variability among the tomato genotypes. Bayesian model-based STRUCTURE analysis assembled the genotypes into 4 (best ΔK = 4) genetic groups. Combined phenotypic and molecular markers proved to be significantly useful for genetic diversity assessment associated with EB disease resistance.

Assessment of genetic diversity among 131 safflower (Carthamus tinctorius L.) accessions using peroxidase gene polymorphism (POGP) markers

BACKGROUND

Safflower (Carthamus tinctorius L.) is an old oilseed crop with a 1.4 GB genome size and its flowers are used for food coloring, dyes and pharmaceutical industries. It was domesticated from its putative wild ancestor Carthamus palestinus about forty-five hundred years ago in the fertile crescent region.The current study was aimed to determine the genetic diversity, population structure and to check the applicability of iPBS-retrotransposons markers.

METHODS AND RESULTS

Eleven POGP primers yielded 70 bands of which 61 were highly polymorphic with 87.14% polymorphism. A great level of genetic variation was examined with higher values of overall gene diversity (0.27), genetic distance (0.53), number of effective alleles (1.46), Shannon's information index (0.41) and polymorphism information contents (0.71). Analysis of molecular variance revealed high genetic variation with 79% within the population. The STRUCTURE, PCoA and Neighbor-joining analysis separated the safflower germplasm into 2 major populations and 1 un-classified population. The accessions which were from Asian countries i.e., China, Afghanistan, Turkey, Iran and Pakistan were genetically similar and clustered together in both populations A and B. The maximum genetic distance was measured 0.88 between Pakistan 26 x Pakistan 24.

CONCLUSION

Findings of this research such as maximum diversity indices, higher PIC values showed the effectiveness and utility of POGP markers for the evaluation of genetic relationships among safflower accessions. The results of this study also showed that POGP markers are less effective compared to ISSRs, iPBS-retrotransposons and DArTSeq markers. AMOVA showed high genetic variation (79%) within a population and maximum genetic distance was found between the accessions Pakistan 26- Pakistan 24 and may be suggested as candidate parents for future breeding activities of safflower. The accessions from the fertile crescent region were clustered together and proved the origin of safflower domestication. This study highlights genetic variation among safflower germplasm and could be helpfull for parental selection and planning for future breeding programs.

Identification of significant marker-trait associations for Fusarium wilt resistance in a genetically diverse core collection of safflower using AFLP and SSR markers

Safflower (Carthamus tinctorius L.), an oilseed crop, is severely affected by Fusarium oxysporum f. sp. carthami (Foc), a fungus causing Fusarium wilt (FW) resulting in up to 80% yield loss. In the present study, we used a panel of 84 diverse accessions from the composite core collection to perform association mapping for FW-resistance. Hydroponics-based screening resulted in categorization of 84 accessions as 31 immune, 19 highly resistant, 9 moderately resistant, 4 moderately susceptible, and 21 highly susceptible. Genotyping with a combination of 155 AFLP and 144 SSR markers revealed substantial genetic differentiation and structure analysis identified three main subpopulations (K = 3) with nearly 35% of admixtures in the panel. Kinship analysis at individual and population level revealed absence of or weak relatedness between the accessions. Association mapping with General Linear Model and Mixed Linear Model identified 4 marker-trait associations (MTAs) significantly linked with the FW-resistance trait. Of these, 3 robust MTAs identified in both the models exhibited phenotypic variance ranging from 4.09 to 6.45%. Locus-128 showing a low P-value and high phenotypic variance was identified as a promising marker-trait association that will facilitate marker-assisted breeding for FW-resistance in safflower.

Genome-wide association studies dissect the G × E interaction for agronomic traits in a worldwide collection of safflowers (Carthamus tinctorius L.)

Genome-wide association studies were conducted using a globally diverse safflower (Carthamus tinctorius L.) Genebank collection for grain yield (YP), days to flowering (DF), plant height (PH), 500 seed weight (SW), seed oil content (OL), and crude protein content (PR) in four environments (sites) that differed in water availability. Phenotypic variation was observed for all traits. YP exhibited low overall genetic correlations (rGoverall) across sites, while SW and OL had high rGoverall and high pairwise genetic correlations (rGij) across all pairwise sites. In total, 92 marker-trait associations (MTAs) were identified using three methods, single locus genome-wide association studies (GWAS) using a mixed linear model (MLM), the Bayesian multi-locus method (BayesR), and meta-GWAS. MTAs with large effects across all sites were detected for OL, SW, and PR, and MTAs specific for the different water stress sites were identified for all traits. Five MTAs were associated with multiple traits; 4 of 5 MTAs were variously associated with the three traits of SW, OL, and PR. This study provided insights into the phenotypic variability and genetic architecture of important safflower agronomic traits under different environments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01295-8.

Whole genome re-sequencing reveals the genetic diversity and evolutionary patterns of Eucommia ulmoides

Eucommia ulmoides (E. ulmoides) is a deciduous perennial tree belonging to the order Garryales, and is known as "living fossil" plant, along with ginkgo (Ginkgo biloba), metaspaca (Metasequoia glyptostroboides) and dove tree (Davidia involucrata Baill). However, the genetic diversity and population structure of E. ulmoides are still  ambiguous nowdays. In this study, we re-sequenced the genomes of 12 E. ulmoides accessions from different major climatic geography regions in China to elucidate the genetic diversity, population structure and evolutionary pattern. By integration of phylogenetic analysis, principal component analysis and population structure analysis based on a number of high-quality SNPs, a total of 12 E. ulmoides accessions were clustered into four different groups. This result is consistent with their geographical location except for group samples from Shanghai and Hunan province. E. ulmoides accessions from Hunan province exhibited a closer genetic relationship with E. ulmoides accessions from Shanghai in China compared with other regions, which is also supported by the result of population structure analyses. Genetic diversity analysis further revealed that E. ulmoides samples in Shanghai and Hunan province were with higher genetic diversity than those in other regions in this study. In addition, we treated the E. ulmoides materials from Shanghai and Hunan province as group A, and the other materials from other places as group B, and then analyzed the evolutionary pattern of E. ulmoides. The result showed the significant differentiation (Fst = 0.1545) between group A and group B. Some candidate highly divergent genome regions were identified in group A by selective sweep analyses, and the function analysis of candidate genes in these regions showed that biological regulation processes could be correlated with the Eu-rubber biosynthesis. Notably, nine genes were identified from selective sweep regions. They were involved in the Eu-rubber biosynthesis and expressed in rubber containing tissues. The genetic diversity research and evolution model of E. ulmoides were preliminarily explored in this study, which laid the foundation for the protection of germplasm resources and the development and utilization of multipurpose germplasm resources in the future.

Identification of a Diverse Core Set Panel of Rice From the East Coast Region of India Using SNP Markers

In India, rice (Oryza sativa L.) is cultivated under a variety of climatic conditions. Due to the fragility of the coastal ecosystem, rice farming in these areas has lagged behind. Salinity coupled with floods has added to this trend. Hence, to prevent genetic erosion, conserving and characterizing the coastal rice, is the need of the hour. This work accessed the genetic variation and population structure among 2,242 rice accessions originating from India’s east coast comprising Andhra Pradesh, Orissa, and Tamil Nadu, using 36 SNP markers, and have generated a core set (247 accessions) as well as a mini-core set (30 accessions) of rice germplasm. All the 36 SNP loci were biallelic and 72 alleles found with average two alleles per locus. The genetic relatedness of the total collection was inferred using the un-rooted neighbor-joining tree, which grouped all the genotypes (2,242) into three major clusters. Two groups were obtained with a core set and three groups obtained with a mini core set. The mean PIC value of total collection was 0.24, and those of the core collection and mini core collection were 0.27 and 0.32, respectively. The mean heterozygosity and gene diversity of the overall collection were 0.07 and 0.29, respectively, and the core set and mini core set revealed 0.12 and 0.34, 0.20 and 0.40 values, respectively, representing 99% of distinctiveness in the core and mini core sets. Population structure analysis showed maximum population at K = 4 for total collection and core collection. Accessions were distributed according to their population structure confirmed by PCoA and AMOVA analysis. The identified small and diverse core set panel will be useful in allele mining for biotic and abiotic traits and managing the genetic diversity of the coastal rice collection. Validation of the 36-plex SNP assay was done by comparing the genetic diversity parameters across two different rice core collections, i.e., east coast and northeast rice collection. The same set of SNP markers was found very effective in deciphering diversity at different genetic parameters in both the collections; hence, these marker sets can be utilized for core development and diversity analysis studies.

Genetic diversity and evolutionary patterns of Taraxacum kok-saghyz Rodin

Taraxacum kok-saghyz Rodin (TKS) is an important potential alternative source of natural inulin and rubber production, which has great significance for the production of industrial products. In this study, we sequenced 58 wild TKS individuals collected from four different geography regions worldwide to elucidate the population structure, genetic diversity, and the patterns of evolution. Also, the first flowering time, crown diameter, morphological characteristics of leaf, and scape of all TKS individuals were measured and evaluated statistically. Phylogenetic analysis based on SNPs and cluster analysis based on agronomic traits showed that all 58 TKS individuals could be roughly divided into three distinct groups: (a) Zhaosu County in Xinjiang (population AB, including a few individuals from population C and D); (b) Tekes County in Xinjiang (population C); and (c) Tuzkol lake in Kazakhstan (population D). Population D exhibited a closer genetic relationship with population C compared with population AB. Genetic diversity analysis further revealed that population expansion from C and D to AB occurred, as well as gene flow between them. Additionally, some natural selection regions were identified in AB population. Function annotation of candidate genes identified in these regions revealed that they mainly participated in biological regulation processes, such as transporter activity, structural molecule activity, and molecular function regulator. We speculated that the genes identified in selective sweep regions may contribute to TKS adaptation to the Yili River Valley of Xinjiang. In general, this study provides new insights in clarifying population structure and genetic diversity analysis of TKS using SNP molecular markers and agronomic traits.

Genetic dissection of flowering time in flax (Linum usitatissimum L.) through single- and multi-locus genome-wide association studies

In a rapidly changing climate, flowering time (FL) adaptation is important to maximize seed yield in flax (Linum usitatissimum L.). However, our understanding of the genetic mechanism underlying FL in this multipurpose crop remains limited. With the aim of dissecting the genetic architecture of FL in flax, a genome-wide association study (GWAS) was performed on 200 accessions of the flax core collection evaluated in four environments. Two single-locus and six multi-locus models were applied using 70,935 curated single nucleotide polymorphism (SNP) markers. A total of 40 quantitative trait nucleotides (QTNs) associated with 27 quantitative trait loci (QTL) were identified in at least two environments. The number of QTL with positive-effect alleles in accessions was significantly correlated with FL (r = 0.77 to 0.82), indicating principally additive gene actions. Nine QTL were significant in at least three of the four environments accounting for 3.06-14.71% of FL variation. These stable QTL spanned regions that harbored 27 Arabidopsis thaliana and Oryza sativa FL-related orthologous genes including FLOWERING LOCUS T (Lus10013532), FLOWERING LOCUS D (Lus10028817), transcriptional regulator SUPERMAN (Lus10021215), and gibberellin 2-beta-dioxygenase 2 (Lus10037816). In silico gene expression analysis of the 27 FL candidate gene orthologous suggested that they might play roles in the transition from vegetative to reproductive phase, flower development and fertilization. Our results provide new insights into the QTL architecture of flowering time in flax, identify potential candidate genes for further studies, and demonstrate the effectiveness of combining different GWAS models for the genetic dissection of complex traits.

Genetic Diversity and Association Analysis among Germplasms of Diospyros kaki in Zhejiang Province Based on SSR Markers

In subtropical to temperate regions, persimmon (Diospyros kaki Thunb.) is an economically important fruit crop cultivated for its edible fruits. Persimmons are distributed abundantly and widely in Zhejiang Province, representing a valuable resource for the breeding of new cultivars and studying the origin and evolution of persimmon. In this study, we elucidated the genetic structures and diversity patterns of 179 persimmon germplasms from 16 different ecologic populations in Zhejiang Province based on the analysis of 17 SSR markers. The results show that there was a medium degree of genetic diversity for persimmon found in Zhejiang Province. With the exception of the Tiantai Mountain and Xin’an River populations, we found extensive gene exchange had occurred among the other populations. The 179 D. kaki germplasms from the 16 populations could be separated into three distinct clusters (I, II, and III) with a higher mean pairwise genetic differentiation index (FST) (0.2714). Nearly all samples of Cluster-I were distributed inland. Cluster-II and Cluster-III contained samples that were widely distributed throughout Zhejiang Province including all samples from the coastal populations and the Northeast Plain populations. In addition, we performed association mapping with nine traits (fruit crude fiber content, fruit calcium content, fruit water content, fruit longitudinal diameter, fruit aspect ratio, seed width, seed length, leaf aspect ratio, and number of lateral veins) using these markers. This led to the identification of 13 significant marker–trait associations (MTAs; p < 0.00044, 0.1/228) using a general linear model, of which, six MTAs with a correlation coefficient (R2) >10% were consistently represented in the general linear model with p < 0.00044 in the two models. The genetic structures and diversity patterns of the persimmon germplasms revealed in this study will provide a reference for the efficient conservation and further utilization of persimmon germplasms. The MTAs identified in this study will be useful for future marker-assisted breeding of persimmon.

Designing a Mini-Core Collection Effectively Representing 3004 Diverse Rice Accessions

Genetic diversity provides the foundation for plant breeding and genetic research. Over 3000 rice genomes were recently sequenced as part of the 3K Rice Genome (3KRG) Project. We added four additional Indian rice accessions to create a panel of 3004 accessions. However, such a large collection of germplasm is difficult to preserve and evaluate. The construction of core and mini-core collections is an efficient method for the management of genetic resources. In this study, we developed a mini-core comprising 520 accessions that captured most of the SNPs and represented all of the phenotypes and geographic regions from the original panel. The mini-core was validated using different statistical analyses and contained representatives from all major rice groups, including japonica, indica, aus/boro, and aromatic/basmati. Genome-wide association analyses of the mini-core panel efficiently reproduced the marker-trait associations identified in the original panel. Haplotype analysis validated the utility of the mini-core panel. In the current era with many ongoing large-scale sequencing projects, such a strategy for mini-core design should be useful in many crops. The rice mini-core collection developed in this study would be valuable for agronomic trait evaluation and useful for rice improvement via marker-assisted molecular breeding.

Genetic Diversity, Population Structure and Marker-Trait Association for 100-Seed Weight in International Safflower Panel Using SilicoDArT Marker Information

Safflower is an important oilseed crop mainly grown in the arid and semi-arid regions of the world. The aim of this study was to explore phenotypic and genetic diversity, population structure, and marker-trait association for 100-seed weight in 94 safflower accessions originating from 26 countries using silicoDArT markers. Analysis of variance revealed statistically significant genotypic effects (p < 0.01), while Turkey samples resulted in higher 100-seed weight compared to Pakistan samples. A Constellation plot divided the studied germplasm into two populations on the basis of their 100-seed weight. Various mean genetic diversity parameters including observed number of alleles (1.99), effective number of alleles (1.54), Shannon's information index (0.48), expected heterozygosity (0.32), and unbiased expected heterozygosity (0.32) for the entire population exhibited sufficient genetic diversity using 12232 silicoDArT markers. Analysis of molecular variance (AMOVA) revealed that most of the variations (91%) in world safflower panel are due to differences within country groups. A model-based structure grouped the 94 safflower accessions into populations A, B, C and an admixture population upon membership coefficient. Neighbor joining analysis grouped the safflower accessions into two populations (A and B). Principal coordinate analysis (PCoA) also clustered the safflower accessions on the basis of geographical origin. Three accessions; Egypt-5, Egypt-2, and India-2 revealed the highest genetic distance and hence might be recommended as candidate parental lines for safflower breeding programs. The mixed linear model i.e., the Q + K model, demonstrated that two DArTseq markers (DArT-45483051 and DArT-15672391) had significant association (p < 0.01) for 100-seed weight. We envisage that identified DArTseq markers associated with 100-seed weight will be helpful to develop high-yielding cultivars of safflower through marker-assisted breeding in the near future.

Development of core-collections for Guizhou tea genetic resources and GWAS of leaf size using SNP developed by genotyping-by-sequencing

An accurate depiction of the genetic relationship, the development of core collection, and genome-wide association analysis (GWAS) are key for the effective exploitation and utilization of genetic resources. Here, genotyping-by-sequencing (GBS) was used to characterize 415 tea accessions mostly collected from the Guizhou region in China. A total of 30,282 high-quality SNPs was used to estimate the genetic relationships, develop core collections, and perform GWAS. We suggest 198 and 148 accessions to represent the core set and mini-core set, which consist of 47% and 37% of the whole collection, respectively, and contain 93–95% of the total SNPs. Furthermore, the frequencies of all alleles and genotypes in the whole set were very well retained in the core set and mini-core set. The 415 accessions were clustered into 14 groups and the core and the mini-core collections contain accessions from each group, species, cultivation status and growth habit. By analyzing the significant SNP markers associated with multiple traits, nine SNPs were found to be significantly associated with four leaf size traits, namely MLL, MLW, MLA and MLSI (P < 1.655E−06). This study characterized the genetic distance and relationship of tea collections, suggested the core collections, and established an efficient GWAS analysis of GBS result.

In-depth genome diversity, population structure and linkage disequilibrium analysis of worldwide diverse safflower (Carthamus tinctorius L.) accessions using NGS data generated by DArTseq technology

Safflower (Carthamus tinctorius L.) is one of the most important oilseed crops for its seed oil rich in unsaturated fatty acids. Precise utilization of diverse genetic resources is fundamental in breeding programs to improve high yield genotypes with desirable traits. In this study, for the first time we report successful application of DArTseq technology; an efficient genotyping-by-sequencing (NGS); to analysis genetic diversity and population structure of 89 safflower accessions from worldwide origins. Totally, 19,639 DArTseq markers (10,130 SilicoDArTs and 9509 SNPs) generated through DArTseq genotyping. After filtering the data, 3431 polymorphic DArTseq markers (1136 SilicoDArTs and 2295 SNPs) used for genetic diversity, population structure and linkage disequilibrium analysis in safflower genotypes. All the SilicoDArT and SNP markers showed high reproducibility and call rate. Polymorphism information content (PIC) values ranged from 0.1 to 0.5, while ≥ 0.50% of SilicoDArTs and ≥ 0.64% SNPs showed PIC values more than median. Genotypes grouping using DArTseq markers resulted in three distinct clusters. Results showed weak correlation between safflower diversity pattern and origins. Analysis of molecular variance revealed that the majority of genetic variation was attributed to the differences among varieties within cluster populations and there was no significant molecular variance between origins. However, safflower of accessions belonged to Iran, Turkey, Pakistan and India indeed appear to be genetically similar and grouped close in referred cluster, while the accessions from Near East (Afghanistan, China) being distinct. Our results were in agreement with hypothesis that safflower domesticated in somewhere west of Fertile Crescent and then expanded through Africa and Europe. Present study using a panel of globally diverse safflower accessions and large number of DArTseq markers set the stage for future analysis of safflower domestication using large germplasm from proposed domestication centers. Also, studied germplasm in this study can be used as a valuable source for future genomic studies in safflower for mapping desirable traits through genome-wide association mapping studies.

Genetic Diversity and Population Structure Analysis of Dalbergia Odorifera Germplasm and Development of a Core Collection Using Microsatellite Markers

Dalbergia odorifera T. Chen (Fabaceae) is a woody tree species indigenous to Hainan Island in China. Due to its high medicinal and commercial value, this tree species has been planted over 3500 ha² in southern China. There is an urgent need for improvement of the D. odorifera germplasm, however, limited information on germplasm collection, conservation, and assessment of genetic resources is available. Therefore, we have built a database of 251 individuals collected across the whole of southern China, which included 42 wild trees and 210 cultivated trees, with the following objectives. (1) Evaluate genetic diversity and population structure of the database using 19 microsatellite markers and (2) develop a core collection for improvement and breeding programs. Totally, the 19 microsatellite markers harbored 77 alleles across the database with the polymorphic information content (PIC) ranging from 0.03 to 0.66. Medium genetic diversity level was inferred by Nei’s gene diversity (0.38), Shannon’s information index (0.65), and observed (0.33) and expected heterozygosity (0.38). Structure analysis showed that four was the optimum cluster size using the model-based Bayesian procedure, and the 251 D. odorifera individuals were grouped into five populations including four pure ones (RP1-4) and one mixed one (MIX) based on their maximum membership coefficients. Among these populations, the expected heterozygosity varied from 0.30 (RP3) to 0.38 (RP4). Analysis of molecular variance (AMOVA) showed 11% genetic variation existed among populations, and moderate population differentiation was inferred by the matrix of pairwise Fst (genetic differentiation among populations), which was in the range of 0.031 to 0.095. Moreover, a core collection of 31 D. odorifera individuals including six wild and 25 cultivated trees was developed, which was only 12.4% of the database but conserved the whole genetic diversity. The results of this study provided additional insight into the genetic structure of the large D. odorifera germplasm, and the core collection will be useful for the efficient and sustainable utilization of genetic resources, as well as efficient improvement in breeding programs.

Mobile genomic element diversity in world collection of safflower (Carthamus tinctorius L.) panel using iPBS-retrotransposon markers

Safflower (Carthamus tinctorius L.) is a multipurpose crop of dry land yielding very high quality of edible oil. Present study was aimed to investigate the genetic diversity and population structure of 131 safflower accessions originating from 28 different countries using 13 iPBS-retrotransposon markers. A total of 295 iPBS bands were observed among which 275 (93.22%) were found polymorphic. Mean Polymorphism information content (0.48) and diversity parameters including mean effective number of alleles (1.33), mean Shannon’s information index (0.33), overall gene diversity (0.19), Fstatistic (0.21), and inbreeding coefficient (1.00) reflected the presence of sufficient amount of genetic diversity in the studied plant materials. Analysis of molecular variance (AMOVA) showed that more than 40% of genetic variation was derived from populations. Model-based structure, principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA) algorithms clustered the 131 safflower accessions into four main populations A, B, C, D and an unclassified population, with no meaningful geographical origin. Most diverse accessions originated from Asian countries including Afghanistan, Pakistan, China, Turkey, and India. Four accessions, Turkey3, Afghanistan4, Afghanistan2, and Pakistan24 were found most genetically distant and might be recommended as a candidate parents for breeding purposes. The findings of this study are most probably supported by the seven similarity centers hypothesis of safflower. This is a first study to explore the genetic diversity and population structure in safflower accessions using the iPBS-retrotransposon markers. The information provided in this work will therefore be helpful for scientists interested in safflower breeding.