Assaying polymorphism at DNA level for genetic diversity diagnostics of the safflower (Carthamus tinctorius L.) world germplasm resources

The Carthamus tinctorius L. genome sequence provides insights into synthesis of unsaturated fatty acids

Domesticated safflower (Carthamus tinctorius L.) is a widely cultivated edible oil crop. However, despite its economic importance, the genetic basis underlying key traits such as oil content, resistance to biotic and abiotic stresses, and flowering time remains poorly understood. Here, we present the genome assembly for C. tinctorius variety Jihong01, which was obtained by integrating Oxford Nanopore Technologies (ONT) and BGI-SEQ500 sequencing results. The assembled genome was 1,061.1 Mb, and consisted of 32,379 protein-coding genes, 97.71% of which were functionally annotated. Safflower had a recent whole genome duplication (WGD) event in evolution history and diverged from sunflower approximately 37.3 million years ago. Through comparative genomic analysis at five seed development stages, we unveiled the pivotal roles of fatty acid desaturase 2 (FAD2) and fatty acid desaturase 6 (FAD6) in linoleic acid (LA) biosynthesis. Similarly, the differential gene expression analysis further reinforced the significance of these genes in regulating LA accumulation. Moreover, our investigation of seed fatty acid composition at different seed developmental stages unveiled the crucial roles of FAD2 and FAD6 in LA biosynthesis. These findings offer important insights into enhancing breeding programs for the improvement of quality traits and provide reference resource for further research on the natural properties of safflower.

Assessing the Genetic Divergence of Onion (Allium Cepa L.) through Morpho-Physiological and Molecular Markers

Onion is one of the most widely cultivated high-value vegetables on the globe. The study of variability and genetic divergence in this crop will facilitate judicious parent selection for future breeding. The present study incorporated 20 onion genotypes to determine the genetic variability and diversity through morpho-physiological and molecular markers. The genotypes were tested for two consecutive years to determine the effects on yield-attributing and bulb-quality traits. Diversity at a molecular level with 31 simple sequence repeat (SSR) markers was determined. Adequate genetic variability existed among the genotypes, with a predominance of additive gene action regarding the inheritance of most of the characteristics. Principal component (PC) analysis detected seven PCs contributing to more than 87% of the genetic diversity in the genotypes. Total yield (kg/plot), Grade ‘C’ bulb (equatorial diameter of the bulb between 3.5 and 4.5 cm), individual bulb weight (g) and phenol content were detected as characteristics positively associated with the first two PCs. A total of 16 primers were amplified, and they produced 32 alleles in the genotypes. Genetic diversity at the morphological and molecular levels indicated moderate diversity among the tested genotypes, with the dissolution of the genotypes into six and three clusters, respectively. The present study on onion facilitates the selection of diverse parents that can generate desirable segregants in future breeding programmes.

Genetic and Genomic Diversity in a Tarwi (Lupinus mutabilis Sweet) Germplasm Collection and Adaptability to Mediterranean Climate Conditions

Lupinus mutabilis (tarwi) is a species of Andean origin with high protein and oil content and regarded as a potential crop in Europe. The success in the introduction of this crop depends in part on in depth knowledge of the intra-specific genetic variability of the collections, enabling the establishment of breeding and conservation programs. In this study, we used morphological traits, Inter-Simple Sequence Repeat markers and genome size to assess genetic and genomic diversity of 23 tarwi accessions under Mediterranean conditions. Phenotypic analyses and yield component studies point out accession LM268 as that achieving the highest seed production, producing large seeds and efficiently using primary branches as an important component of total yield, similar to the L. albus cultivars used as controls. By contrast, accession JKI-L295 presents high yield concentrated on the main stem, suggesting a semi-determinate development pattern. Genetic and genomic analyses revealed important levels of diversity, however not relatable to phenotypic diversity, reflecting the recent domestication of this crop. This is the first study of genome size diversity within L. mutabilis, revealing an average size of 2.05 pg/2C (2001 Mbp) with 9.2% variation (1897–2003 Mbp), prompting further studies for the exploitation of this diversity.

Genetic structure and diversity analysis of tall fescue populations by EST-SSR and ISSR markers

Tall fescue is a perennial cool-season grass with economic importance especially in temperate regions of the northern hemisphere. This study was done to assess the genetic diversity and population structure of 90 tall fescue populations and cultivars using ISSR and EST-SSR markers in order to categorize valuable populations for breeding programs and to construct the core collection of tall fescue collection in Iran. The 10 EST-SSR primer pairs amplified 92 alleles. The allele numbers varied from 4 to 13 alleles per locus with an average of 9.2 alleles, of which 84 (90.6%) were polymorphic with an average of 8.4 polymorphic bands per primer. The 39 ISSR primers totally produced 387 scorable bands, of which 335 (86.6%) were polymorphic with an average of 8.6 polymorphic bands per primer. The amplified markers by ISSR primers varied from 6 to 14 markers per primer with an average of 9.92 markers per primer. The 90 tall fescue populations using both EST-SSR and ISSR data were classified into two clusters by UPGMA method that was coincide with PCA and structure analysis results. The turf-type and forage-type populations were clearly separated. Based on the results, the Iranian populations provide a valuable and novel germplasm to employ in tall fescue varietal improvement programs for both forage and turf-type applications. This progression is an important step to introduce this collection for development of a core collection of tall fescue germplasm in Iran.

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.

Genetic diversity and population structure of Haloxylon salicornicum moq. in Kuwait by ISSR markers

Haloxylon salicornicum moq. Bunge ex Boiss (Rimth) is one of the native plants of Kuwait, extensively depleting through the anthropogenic activities. It is important to conserve Haloxylon community in Kuwait as it can tolerate extreme adverse conditions of drought and salinity to be potentially used in the desert and urban revegetation and greenery national programs. Therefore, a set of 16 inter simple sequence repeat (ISSR) markers were used to assess genetic diversity and population structure of 108 genotypes from six locations in Kuwait. The ISSR primers produced 195 unambiguous and reproducible bands out of which 167 bands were polymorphic (86.5%) with a mean PIC value of 0.31. The overall average values of Nei’s gene diversity (h') and Shannon’s diversity index (I) were 0.254 and 0.375, respectively. Results of AMOVA revealed high genetic variations within populations (77.8%) and low among populations (22%). The values of Fixation index (F_ST = 0.22; P = 0.0), Genetic differentiation (G_{ST =} 0.262; G’_ST = 0.327; D = 0.335 and Gene flow (N_M = 0.880) were indicative of heterozygous populations. The results of STRUCTURE and split decomposition analysis suggested that the Rimth accessions of Kuwait can be grouped into five and six subpopulations, respectively. Principal coordinate analysis (PCoA) grouped them into three clusters. The pairwise Nei’s genetic distances (D_S) among populations demonstrated a narrow range from 0.047 to 0.187 (Scale-0.0 to 1.0). The Mantel’s test revealed a weak correlation (r²- 0.188; P-0.013) between the genetic distance and geographic distances. Our results suggest that the narrowly distributed Haloxylon community in Kuwait demonstrated a high genetic diversity within the populations however the overall population structure was weak.

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

Carthamus tinctorius L. (safflower) is an important oilseed crop producing seed oil rich in unsaturated fatty acids. Scarcity of identified marker-trait associations is a major limitation toward development of successful marker-assisted breeding programs in safflower. In the present study, a safflower panel (CartAP) comprising 124 accessions derived from two core collections was assayed for its suitability for association mapping. Genotyping of CartAP using microsatellite markers revealed significant genetic diversity indicated by Shannon information index (H = 0.7537) and Nei's expected heterozygosity (I = 0.4432). In Principal Coordinate Analysis, the CartAP accessions were distributed homogeneously in all quadrants indicating their diverse nature. Distance-based Neighbor Joining analysis did not delineate the CartAP accessions in consonance with their geographical origin. Bayesian analysis of population structure of CartAP demonstrated the unstructured nature of the association panel. Kinship analysis at population (G_ij ) and individual level (F_ij ) revealed absence of or weak relatedness between the CartAP accessions. The above parameters established the suitability of CartAP for association mapping. We performed association mapping using phenotypic data for eight traits of agronomic value (viz., seed oil content, oleic acid, linoleic acid, plant height, number of primary branches, number of capitula per plant, 100-seed weight and days to 50% flowering) available for two growing seasons (2011-2012 and 2012-2013) through General Linear Model and Mixed Linear Model. Our study identified ninety-six significant marker-trait associations (MTAs; P < 0.05) of which, several MTAs with correlation coefficient (R²) > 10% were consistently represented in both models and in both seasons for traits viz., oil content, oleic acid content, linoleic acid content and number of primary branches. Several MTAs with high R²-values were detected either in a majority or in some environments (models and/or seasons). Many MTAs were also common between traits (viz., oleic/linoleic acid content; plant height/days to 50% flowering; number of primary branches/number of capitula per plant) that showed positive or negative correlation in their phenotypic values. The marker-trait associations identified in this study will facilitate marker-assisted breeding and identification of genetic determinants of trait variability.

Exploration of Genetic and Genomic Resources for Abiotic and Biotic Stress Tolerance in Pearl Millet

Pearl millet is one of the most important small-grained C4 Panicoid crops with a large genome size (∼2352 Mb), short life cycle and outbreeding nature. It is highly resilient to areas with scanty rain and high temperature. Pearl millet is a nutritionally superior staple crop for people inhabiting hot, drought-prone arid and semi-arid regions of South Asia and Africa where it is widely grown and used for food, hay, silage, bird feed, building material, and fuel. Having excellent nutrient composition and exceptional buffering capacity against variable climatic conditions and pathogen attack makes pearl millet a wonderful model crop for stress tolerance studies. Pearl millet germplasm show a large range of genotypic and phenotypic variations including tolerance to abiotic and biotic stresses. Conventional breeding for enhancing abiotic and biotic stress resistance in pearl millet have met with considerable success, however, in last few years various novel approaches including functional genomics and molecular breeding have been attempted in this crop for augmenting yield under adverse environmental conditions, and there is still a lot of scope for further improvement using genomic tools. Discovery and use of various DNA-based markers such as EST-SSRs, DArT, CISP, and SSCP-SNP in pearl millet not only help in determining population structure and genetic diversity but also prove to be important for developing strategies for crop improvement at a faster rate and greater precision. Molecular marker-based genetic linkage maps and identification of genomic regions determining yield under abiotic stresses particularly terminal drought have paved way for marker-assisted selection and breeding of pearl millet cultivars. Reference collections and marker-assisted backcrossing have also been used to improve biotic stress resistance in pearl millet specifically to downy mildew. Whole genome sequencing of pearl millet genome will give new insights for processing of functional genes and assist in crop improvement programs through molecular breeding approaches. This review thus summarizes the exploration of pearl millet genetic and genomic resources for improving abiotic and biotic stress resistance and development of cultivars superior in stress tolerance.

Development of Genomic Microsatellite Markers in Carthamus tinctorius L. (Safflower) Using Next Generation Sequencing and Assessment of Their Cross-Species Transferability and Utility for Diversity Analysis

BACKGROUND

Safflower (Carthamus tinctorius L.), an Asteraceae member, yields high quality edible oil rich in unsaturated fatty acids and is resilient to dry conditions. The crop holds tremendous potential for improvement through concerted molecular breeding programs due to the availability of significant genetic and phenotypic diversity. Genomic resources that could facilitate such breeding programs remain largely underdeveloped in the crop. The present study was initiated to develop a large set of novel microsatellite markers for safflower using next generation sequencing.

PRINCIPAL FINDINGS

Low throughput genome sequencing of safflower was performed using Illumina paired end technology providing ~3.5X coverage of the genome. Analysis of sequencing data allowed identification of 23,067 regions harboring perfect microsatellite loci. The safflower genome was found to be rich in dinucleotide repeats followed by tri-, tetra-, penta- and hexa-nucleotides. Primer pairs were designed for 5,716 novel microsatellite sequences with repeat length ≥ 20 bases and optimal flanking regions. A subset of 325 microsatellite loci was tested for amplification, of which 294 loci produced robust amplification. The validated primers were used for assessment of 23 safflower accessions belonging to diverse agro-climatic zones of the world leading to identification of 93 polymorphic primers (31.6%). The numbers of observed alleles at each locus ranged from two to four and mean polymorphism information content was found to be 0.3075. The polymorphic primers were tested for cross-species transferability on nine wild relatives of cultivated safflower. All primers except one showed amplification in at least two wild species while 25 primers amplified across all the nine species. The UPGMA dendrogram clustered C. tinctorius accessions and wild species separately into two major groups. The proposed progenitor species of safflower, C. oxyacantha and C. palaestinus were genetically closer to cultivated safflower and formed a distinct cluster. The cluster analysis also distinguished diploid and tetraploid wild species of safflower.

CONCLUSION

Next generation sequencing of safflower genome generated a large set of microsatellite markers. The novel markers developed in this study will add to the existing repertoire of markers and can be used for diversity analysis, synteny studies, construction of linkage maps and marker-assisted selection.

Efficacy of RAPD, ISSR and DAMD markers in assessment of genetic variability and population structure of wild Musa acuminata colla

North east India is considered as one of the major biodiversity hotspots worldwide and centre of origin of several plant species including Musa. Musa acuminata Colla is known to be one of the wild progenitors of cultivated bananas and plantains. Three single primer based DNA marker techniques viz., random amplified polymorphic DNA (RAPD), inter simple sequence repeat (ISSR) and directed amplification of minisatellites DNA (DAMD) were used for diversity diagnostics among 25 genotypes of wild M. acuminata collected from Meghalaya province of north east India. A total of 58 primers (26-RAPD, 21-ISSR, and11-DAMD) yielded 451 DNA fragments, of which 395 (87.58 %) were found to be polymorphic in nature. The polymorphic information content (PIC) values were almost identical for each marker system. The resolving power of the marker system was found to be highest in RAPD (3.96) whereas ISSR resolved highest marker index (16.39) in the study. Selected amplicon data obtained through single primer amplification reactions were utilized for determination of diversity within and among the populations of M. acuminata. Nei's genetic differentiation (Gst) value (0.451) indicated higher proportion of the genetic variation within the populations which is supported by the AMOVA analysis (88 %). The study provides insight into the efficacy of RAPD, ISSR and DAMD to analyse the genetic variation existing in the wild Musa germplasm, which can further be exploited for quality trait improvement and domestication of such important horticultural crops. The genetic diversity based population structure may shed light on the genetic basis of speciation and evolution of various species within the genus Musa.

Genetic assessment of safflower (Carthamus tinctorius L.) collection with microsatellite markers acquired via pyrosequencing method

A genetic evaluation of safflower germplasm collections derived from different geographical regions and countries will provide useful information for sustainable conservation and the utilization of genetic diversity. However, the molecular marker information is limited for evaluation of genetic diversity of safflower germplasm. In this study, we acquired 509 putative genomic SSR markers for sufficient genome coverage using next-generation sequencing methods and characterized thirty polymorphic SSRs in safflower collection composed of 100 diverse accessions. The average allele number and expected heterozygosity were 2.8 and 0.386, respectively. Analysis of population structure and phylogeny based on thirty SSR profiles revealed genetic admixture between geographical regions contrary to genetic clustering. However, the accessions from Korea were genetically conserved in distinctive groups in contrast to other safflower gene pool. In conclusion, these new genomic SSRs will facilitate valuable studies to clarify genetic relationships as well as conduct population structure analyses, genetic map construction and association analysis for safflower.

Sequence analysis of KpnI repeat sequences to revisit the phylogeny of the Genus Carthamus L

Repetitive DNA sequences constitute a significant proportion of eukaryotic genomes. Knowledge about the distribution of repetitive DNA sequences is necessary in order to gain insights into the organization, evolution and behavior of eukaryotic genomes. Therefore, we used two repetitive DNA sequences pCtKpnI-I and pCtKpnI-II, earlier reported in Carthamus tinctorius L. to study the phylogeny and to revise the taxonomic status of the taxa belonging to the genus. The study unraveled two major lines within the genus Carthamus; one line included all the diploid taxa (2n = 24) and the other line comprised the taxa with 2n = 20 and the polyploid taxa (2n = 44 and 64). The results of the present study will prove useful in molecular breeding for improving some targeted agronomic traits in genus Carthamus.

Integration of gene-based markers in a pearl millet genetic map for identification of candidate genes underlying drought tolerance quantitative trait loci

BACKGROUND

Identification of genes underlying drought tolerance (DT) quantitative trait loci (QTLs) will facilitate understanding of molecular mechanisms of drought tolerance, and also will accelerate genetic improvement of pearl millet through marker-assisted selection. We report a map based on genes with assigned functional roles in plant adaptation to drought and other abiotic stresses and demonstrate its use in identifying candidate genes underlying a major DT-QTL.

RESULTS

Seventy five single nucleotide polymorphism (SNP) and conserved intron spanning primer (CISP) markers were developed from available expressed sequence tags (ESTs) using four genotypes, H 77/833-2, PRLT 2/89-33, ICMR 01029 and ICMR 01004, representing parents of two mapping populations. A total of 228 SNPs were obtained from 30.5 kb sequenced region resulting in a SNP frequency of 1/134 bp. The positions of major pearl millet linkage group (LG) 2 DT-QTLs (reported from crosses H 77/833-2 × PRLT 2/89-33 and 841B × 863B) were added to the present consensus function map which identified 18 genes, coding for PSI reaction center subunit III, PHYC, actin, alanine glyoxylate aminotransferase, uridylate kinase, acyl-CoA oxidase, dipeptidyl peptidase IV, MADS-box, serine/threonine protein kinase, ubiquitin conjugating enzyme, zinc finger C- × 8-C × 5-C × 3-H type, Hd3, acetyl CoA carboxylase, chlorophyll a/b binding protein, photolyase, protein phosphatase1 regulatory subunit SDS22 and two hypothetical proteins, co-mapping in this DT-QTL interval. Many of these candidate genes were found to have significant association with QTLs of grain yield, flowering time and leaf rolling under drought stress conditions.

CONCLUSIONS

We have exploited available pearl millet EST sequences to generate a mapped resource of seventy five new gene-based markers for pearl millet and demonstrated its use in identifying candidate genes underlying a major DT-QTL in this species. The reported gene-based markers represent an important resource for identification of candidate genes for other mapped abiotic stress QTLs in pearl millet. They also provide a resource for initiating association studies using candidate genes and also for comparing the structure and function of distantly related plant genomes such as other Poaceae members.

Single primer amplification reaction (SPAR) reveals inter- and intra-specific natural genetic variation in five species of Cymbidium (Orchidaceae)

A total of 53 primers belonging to three SPAR methods, viz. RAPD, ISSR and DAMD, collectively produced 456 polymorphic amplicons with 96.6% polymorphism at inter-specific level in five species of Cymbidium, viz. C. aloifolium, C. mastersii, C. elegans, C. eburneum and C. tigrinum, whereas at intra-specific level, the observed polymorphism ranged from 51.2% to 77.1% among them. Three SPARs collectively revealed 25 unique species-specific amplicons; most of them were amplified with RAPD and DAMD primers besides few bands which were either missed (absent) or lost (heterozygosity). UPGMA clustering evidently distinguished the representatives of C. aloifolium and C. tigrinum, with distinct genetic distance, which may be due to their entirely different habitats as well as discrete morphological characteristics. Upon analysis of the data generated, all the three SPAR methods, either independently and/or in combination, revealed wide range of genetic variation between and within five species of Cymbidium. Comparison of matrix of individual SPAR method revealed that analysis of natural genetic variation using combination of SPAR methods, rather than an isolated approach, is highly effective. The critical analyses of the amplicon data are indicative of DAMD as the most powerful SPAR method by showing highest resolving power (Rp) followed by ISSR and RAPD. Alternatively, the total polymorphic information content was highest in case of RAPD followed by other two SPAR methods. Thus, the present investigation for the first time provides a valuable baseline data for genetic variation at inter- and intra-specific levels in horticultural Cymbidiums and also addresses conservation concerns.

Development of molecular markers and linkage maps for the Carthamus species C. tinctorius and C. oxyacanthus

A set of SSR and RFLP markers for safflower ( Carthamus tinctorius ) and jeweled distaff thistle ( C. oxyacanthus ) was generated from cDNA and genomic libraries and by mining public and proprietary sequence databases. In total, 1412 PCR-based markers and 75 RFLP markers were screened and polymorphic loci were mapped in an intraspecific F2 population of C. tinctorius and an interspecific BC1 population of C. tinctorius × C. oxyacanthus. The two populations shared one common parent and the resulting linkage maps could be compared for synteny. The level of polymorphism was low in both populations and only 8.2% and 13.7% of the analyzed markers could be mapped in the intraspecific and interspecific maps, respectively. The two maps showed significant colinearity of markers in several regions and an apparent translocation or inversion event on one linkage group. Noticeable segregation distortion was found on one linkage group of the C. tinctorius map and dense clustering of loci occurred on several linkage groups of the C. oxyacanthus map. The two maps represent the first major linkage analysis of Carthamus species. The molecular tools will be useful for a variety of genetic and genomic applications in safflower and its related species and have been used in our laboratory to map a flower color gene in C. tinctorius.