Newly developed SSR markers reveal genetic diversity and geographical clustering in spinach (Spinacia oleracea)

Genome-wide identification, characterization and primer designing of simple sequence repeats across Leguminosae family

Legumes are important clade of commercially important family Leguminosae that mainly include medicinal, flowering and edible plants. Although the genomic sequence of legumes is accessible, only the limited number of effective simple sequence repeat markers has been identified by prior research. Additional polymorphic simple sequence repeats marker discovery will aid in the genetics and breeding of legumes. In this study, 13 complete genome sequences were screened for the identification of chromosome-wise simple sequence repeats (SSRs) and 1,866,861 SSRs were identified. Based on the study, it was observed that the number of SSRs in non-coding region was more as compared to coding region and frequency of mononucleotides was highest followed by di-nucleotides while penta- and hexa-nucleotide repeats were least frequent one. The identified genome-wide SSRs and newly developed SSR markers, primers and their mapping will provide a powerful means for genetic researches across Leguminosae plants, including genetic diversity and evolutionary origin analysis, fingerprinting, QTL mapping and marker-assisted selection for breeding as well as comparative genomic analysis studies.

Genome-wide characterization and development of SSR markers for genetic diversity analysis in northwestern Himalayas Walnut (Juglans regia L.)

In the present study, we designed and validated genome-wide polymorphic SSR markers (110 SSRs) by mining the walnut genome. A total of 198,924 SSR loci were identified. Among these, successful primers were designed for 162,594 (81.73%) SSR loci. Dinucleotides were the most predominant accounting for 88.40% (175,075) of total SSRs. The SSR frequency was 377.312 SSR/Mb and it showed a decreasing trend from dinucleotide to octanucleotide motifs. We identified 20 highly polymorphic SSR markers and used them to genotype 72 walnut accessions. Over all, we obtained 118 alleles that ranged from 2 to 12 with an average value of 5.9. The higher SSR PIC values indicate their robustness in discriminating walnut genotypes. Heat map, PCA, and population structure categorized 72 walnut genotypes into 2 distinct clusters. The genetic variation within population was higher than among population as inferred by analysis of molecular variance (AMOVA). For walnut improvement, it is necessary to have a large repository of SSRs with high discriminative power. The present study reports 150,000 SSRs, which is the largest SSR repository for this important nut crop. Scientific communities may use this repository for walnut improvement such as QTL mapping, genetic studies, linkage map construction, and marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03563-6.

The Generation of Nitric Oxide from Aldehyde Dehydrogenase-2: The Role of Dietary Nitrates and Their Implication in Cardiovascular Disease Management

Reduced bioavailability of the nitric oxide (NO) signaling molecule has been associated with the onset of cardiovascular disease. One of the better-known and effective therapies for cardiovascular disorders is the use of organic nitrates, such as glyceryl trinitrate (GTN), which increases the concentration of NO. Unfortunately, chronic use of this therapy can induce a phenomenon known as "nitrate tolerance", which is defined as the loss of hemodynamic effects and a reduction in therapeutic effects. As such, a higher dosage of GTN is required in order to achieve the same vasodilatory and antiplatelet effects. Mitochondrial aldehyde dehydrogenase 2 (ALDH2) is a cardioprotective enzyme that catalyzes the bio-activation of GTN to NO. Nitrate tolerance is accompanied by an increase in oxidative stress, endothelial dysfunction, and sympathetic activation, as well as a loss of the catalytic activity of ALDH2 itself. On the basis of current knowledge, nitrate intake in the diet would guarantee a concentration of NO such as to avoid (or at least reduce) treatment with GTN and the consequent onset of nitrate tolerance in the course of cardiovascular diseases, so as not to make necessary the increase in GTN concentrations and the possible inhibition/alteration of ALDH2, which aggravates the problem of a positive feedback mechanism. Therefore, the purpose of this review is to summarize data relating to the introduction into the diet of some natural products that could assist pharmacological therapy in order to provide the NO necessary to reduce the intake of GTN and the phenomenon of nitrate tolerance and to ensure the correct catalytic activity of ALDH2.

In Silico Screening and Development of Microsatellite Markers for Genetic Analysis in Perca fluviatilis

Perca fluviatilis is an economically important species of freshwater fish. To understand the genetic structure of P. fluviatilis in China, 268 samples were collected from Wulungu Lake (WL), Jili Lake (JL), the Wulungu River (WR), and the Kalaeerqisi River (KR). These samples were then analyzed using microsatellite markers. A total of 98,425 microsatellite markers were developed based on the genomic data, and 29 polymorphic microsatellite markers were selected to analyze genetic diversity in this study. The number of alleles (Na) and observed heterozygosity (Ho) per population ranged from 4.621 (KR) to 11.172 (WL) and from 0.510 (KR) to 0.716 (JL), respectively. The results of the polymorphic information content (PIC) showed that the WL, JL, and WR populations were highly polymorphic (PIC≥ 0.5) and that the KR population was moderately polymorphic (0.25 ≤ PIC < 0.5). The genetic differentiation coefficient (Fst) among the four P. fluviatilis populations was 0.074, indicating moderate genetic differentiation among the populations in Xinjiang. The reason for the significant difference between the rivers and lakes could be the presence of a dam blocking the flow of P. fluviatilis. The development of microsatellite markers provides support for population genetics in the future. The evaluation of the genetic structure of P. fluviatilis in Xinjiang provides a reference for the reproduction and conservation of P. fluviatilis.

Genetic Diversity, Structure, and Selective Sweeps in Spinacia turkestanica Associated With the Domestication of Cultivated Spinach

Genotype-by-sequencing (GBS) was used to explore the genetic diversity and structure of Spinacia turkestanica, and the selective sweeps involved in domestication of cultivated spinach, S. oleracea, from S. turkestanica. A total 7,065 single nucleotide polymorphisms (SNPs) generated for 16 Spinacia oleracea and 76 S. turkestanica accessions placed the S. oleracea accessions in one group, Q1, and the 76 S. turkestanica accessions, which originated from Central Asia, in two distinct groups, Q2 and Q3. The Q2 group shared greater genetic identity with the S. oleracea accessions, Q1, than the Q3 S. turkestanica group. Likewise, the S. oleracea Q1 group had a smaller Fst (0.008) with the Q2 group than with the Q3 group (Fst = 0.012), and a greater gene flow (Nm = 30.13) with the Q2 group than with the Q3 group (Nm = 21.83). The Q2 accessions originated primarily from Uzbekistan while the Q3 accessions originated mostly from Tajikistan. The Zarafshan Mountain Range appears to have served as a physical barrier that largely separated members of the Q2 and Q3 groups of S. turkestanica. Accessions with admixtures of Q2 and Q3 were collected primarily from lower elevations at the southern end of the Zarafshan Mountain Range in Uzbekistan. Selective sweep regions identified at 32, 49, and 52 Mb on chromosomes 1, 2, and 3, respectively, appear to have played a vital role in the domestication of S. oleracea as they are correlated with important domestication traits, including day length sensitivity for bolting (flowering). High XP-CLR scores at the 52 Mb genomic region of chromosome three suggest that a selective sweep at this region was responsible for early differentiation of S. turkestanica into two groups in Central Asia.

Quantitative trait loci (QTL) analysis of leaf related traits in spinach (Spinacia oleracea L.)

BACKGROUND

Spinach (Spinacia oleracea L.) is an important leafy vegetable crop, and leaf-related traits including leaf length, leaf width, and petiole length, are important commercial traits. However, the underlying genes remain unclear. The objective of the study was to conduct QTL mapping of leaf-related traits in spinach.

RESULTS

A BC₁ population was used to construct the linkage map and for QTL mapping of leaf length, leaf width, petiole length, and the ratio of leaf length to width in 2015 and 2019. Two genetic linkage maps were constructed by specific locus amplified fragment sequencing (SLAF-seq), and kompetitive allele specific PCR (KASP) technology, respectively using BC₁ population in 2015. Based on the results of 2015, the specific linkage groups (LG) detected QTLs were generated using BC₁ population in 2019. A total of 13 QTLs were detected for leaf-related traits, only five QTLs being repeatedly detected in multiple years or linkage maps. Interestingly, the major QTLs of leaf length, petiole length, and the ratio of leaf length to width were highly associated with the same SNP markers (KM3102838, KM1360385 and KM2191098). A major QTL of leaf width was mapped on chromosome 1 from 41.470-42.045 Mb. And 44 genes were identified within the region. Based on the GO analysis, these genes were significantly enriched on ribonuclease, lyase activity, phosphodiester bond hydrolysis process, and cell wall component, thus it might change cell size to determine leaves shape.

CONCLUSIONS

Five QTLs for leaf-related traits were repeatedly detected at least two years or linkage maps. The major QTLs of leaf length, petiole length, and the ratio of leaf length to width were mapped on the same loci. And three genes (Spo10792, Spo21018, and Spo21019) were identified as important candidate genes for leaf width.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Mining the red deer genome (CerEla1.0) to develop X-and Y-chromosome-linked STR markers

Microsatellites are widely applied in population and forensic genetics, wildlife studies and parentage testing in animal breeding, among others, and recently, high-throughput sequencing technologies have greatly facilitated the identification of microsatellite markers. In this study the genomic data of Cervus elaphus (CerEla1.0) was exploited, in order to identify microsatellite loci along the red deer genome and for designing the cognate primers. The bioinformatics pipeline identified 982,433 microsatellite motifs genome-wide, assorted along the chromosomes, from which 45,711 loci mapped to the X- and 1096 to the Y-chromosome. Primers were successfully designed for 170,873 loci, and validated with an independently developed autosomal tetranucleotide STR set. Ten X- and five Y-chromosome-linked microsatellites were selected and tested by two multiplex PCR setups on genomic DNA samples of 123 red deer stags. The average number of alleles per locus was 3.3, and the average gene diversity value of the markers was 0.270. The overall observed and expected heterozygosities were 0.755 and 0.832, respectively. Polymorphic Information Content (PIC) ranged between 0.469 and 0.909 per locus with a mean value of 0.813. Using the X- and Y-chromosome linked markers 19 different Y-chromosome and 72 X-chromosome lines were identified. Both the X- and the Y-haplotypes split to two distinct clades each. The Y-chromosome clades correlated strongly with the geographic origin of the haplotypes of the samples. Segregation and admixture of subpopulations were demonstrated by the use of the combination of nine autosomal and 16 sex chromosomal STRs concerning southwestern and northeastern Hungary. In conclusion, the approach demonstrated here is a very efficient method for developing microsatellite markers for species with available genomic sequence data, as well as for their use in individual identifications and in population genetics studies.

Development of a Large Gene-Associated SSR Marker Set and in-Depth Genetic Characterization in Scarlet Sage

Salvia splendens, scarlet or tropical sage, is a tender perennial herbaceous flowering plant popularly grown in public and private gardens all over the world. In this study, we developed a set of simple sequence repeats (SSRs) from genome-wide sequences to assess the genetic diversity and population structure among 112 cultivars. We obtained 364,379 SSRs by mining scarlet sage’s recently published whole genome sequence; 14,545 gene-associated SSR loci were identified in 2 kb gene flanking regions. Among the 768 gene-associated SSR primer sets we screened, 576 loci successfully amplified in DNA pools of 3–4 different cultivars, of which 271 remained polymorphic when tested across eight individual plants. We searched for the related gene functions attributable to these gene-associated SSRs using diverse databases, resulting in 259 Non-redundant matching sequences, 205 individual Gene Ontology (GO) terms, 236 assigned to eukaryotic orthologous groups, and 67 KEGG-annotated (Kyoto Encyclopedia of Genes and Genomes) sequences. We finally selected 41 polymorphic SSR loci to infer genetic diversity and population structure among 112 S. splendens accessions. Based on the developed gene-associated SSRs, clustering analyses consistently revealed two distinct genetic groups within the core collection of S. splendens cultivars. This work developed and characterized an exhaustive set of genome-wide gene-associated SSR markers for scarlet sage. These SSRs can provide species identification, genetic diversity and population structure information for S. splendens, and will therefore be important tools for the management and protection of S. splendens germplasm.

Genetic diversity of Ziziphus mauritiana germplasm based on SSR markers and ploidy level estimation

A set of reliable SSR markers were developed for Ziziphus mauritiana. The genetic relationship of Z. mauritiana germplasms was generally consistent with their geographical origin, and low diversity in the maternal lineage was revealed. Ziziphus mauritiana, known as Indian jujube, is an important fruit crop that is native to southern Asia and eastern Africa. There is a variety of germplasm resources, and particularly many new cultivars were selected and introduced into wide tropical regions in recent years. However, there are few practical molecular markers for cultivar authentication and genetic analysis. In this study, we developed 55 polymorphic nuclear SSR markers based on restriction-site associated DNA sequences and transcriptome sequencing. We selected 14 robust nSSR markers for further analysis of 117 Z. mauritiana accessions from four countries (45 from China, 39 from Vietnam, 25 from Pakistan and 8 from Myanmar). In total, 137 alleles were detected and DNA fingerprints for each accession were constructed. Cluster analysis based on the unweighted pair group method with arithmetic mean displayed that most accessions clustered consistently with their geographic origin. In addition, there was common and high degree polyploidization based on nSSR and flow cytometry analyses. Only two of the 50 SSR loci in noncoding regions from the chloroplast genome had polymorphisms, and 5 haplotypes in total were identified among the 117 accessions. Haplotype C with 89 accessions was the most dominant haplotype and presented in four countries. This indicates low diversity in the maternal lineage of tested Z. mauritiana germplasm. Our research provides reliable marker resources for cultivar authentication and new insights into the genetic diversity, polyploidization and domestication of Z. mauritiana.

De Novo Transcriptome Analysis of Dalbergia odorifera and Transferability of SSR Markers Developed from the Transcriptome

Dalbergia odorifera T. Chen (Fabaceae), indigenous to Hainan Island, is a precious rosewood (Hainan hualimu) in China. However, only limited genomic information is available which has resulted in a lack of molecular markers, limiting the development and utilization of the germplasm resources. In this study, we aim to enrich genomic information of D. odorifera, and develop a series of transferable simple sequence repeat (SSR) markers for Dalbergia species. Therefore, we performed transcriptome sequencing for D. odorifera by pooling leaf tissues from three trees. A dataset of 138,516,418 reads was identified and assembled into 115,292 unigenes. Moreover, 35,774 simple sequence repeats (SSRs) were identified as potential SSR markers. A set of 19 SSR markers was successfully transferred across species of Dalbergia odorifera T. Chen, Dalbergia tonkinensis Prain, and Dalbergia cochinchinensis Pierre ex Laness. In total, 112 alleles (3–13 alleles/locus) were presented among 60 Dalbergia trees, and polymorphic information content ranged from 0.38 to 0.75. The mean observed and mean expected heterozygosity was 0.34 and 0.40 in D. odorifera, 0.27 and 0.32 in D. tonkinensis, and 0.29 and 0.33 in D. cochinchinensis, respectively. The cluster analysis classified these 60 trees into three major groups according to the three Dalbergia species based on the genetic similarity coefficients, indicating these newly developed transferable markers can be used to explore the relationships among Dalbergia species and assist genetic research. All these unigenes and SSR markers will be useful for breeding programs in the future.

The landscape of transposable elements and satellite DNAs in the genome of a dioecious plant spinach (Spinacia oleracea L.)

Background

Repetitive sequences, including transposable elements (TEs) and satellite DNAs, occupy a considerable portion of plant genomes. Analysis of the repeat fraction benefits the understanding of genome structure and evolution. Spinach (Spinacia oleracea L.), an important vegetable crop, is also a model dioecious plant species for studying sex determination and sex chromosome evolution. However, the repetitive sequences of the spinach genome have not been fully investigated.

Results

We extensively analyzed the repetitive components of draft spinach genome, especially TEs and satellites, by different strategies. A total of 16,002 full-length TEs were identified. Among the most abundant long terminal repeat (LTR) retrotransposons (REs), Copia elements were overrepresented compared with Gypsy ones. Angela was the most dominating Copia lineage; Ogre/Tat was the most abundant Gypsy lineage. The mean insertion age of LTR-REs was 1.42 million years; approximately 83.7% of these elements were retrotransposed during the last two million years. RepeatMasker totally masked about 64.05% of the spinach genome, with LTR-REs, non-LTR-REs, and DNA transposons occupying 49.2, 2.4, and 5.6%, respectively. Fluorescence in situ hybridization (FISH) analysis showed that most LTR-REs dispersed all over the chromosomes, by contrast, elements of CRM lineage were distributed at the centromeric region of all chromosomes. In addition, Ogre/Tat lineage mainly accumulated on sex chromosomes, and satellites Spsat2 and Spsat3 were exclusively located at the telomeric region of the short arm of sex chromosomes.

Conclusions

We reliably annotated the TE fraction of the draft genome of spinach. FISH analysis indicates that Ogre/Tat lineage and the sex chromosome-specific satellites DNAs might participate in sex chromosome formation and evolution. Based on FISH signals of microsatellites, together with 45S rDNA, a fine karyotype of spinach was established. This study improves our knowledge of repetitive sequence organization in spinach genome and aids in accurate spinach karyotype construction.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0147-6) contains supplementary material, which is available to authorized users.

Genome-wide characterization of microsatellites and genetic diversity assessment of spinach in the Chinese germplasm collection

Spinach is a nutritional leafy green vegetable, and it also serves as a model species for studying sex chromosome evolution. Genetic marker development and genome structure analysis are important in breeding practice and theoretical evolution studies of spinach. In this study, the frequency and distribution of different microsatellites in the recently released draft spinach genome were characterized. A total of 261,002 perfect microsatellites were identified (estimated frequency: ~262.1 loci/Mbp). The most abundant microsatellites were tetranucleotide and trinucleotide, accounting for 33.2% and 27.7% of the total number of microsatellites, respectively. A total of 105 primer pairs were designed and screened, and 34 were polymorphic among the detected spinach cultivars. Combined with seven primer sets developed previously, 41 primer pairs were used to investigate genetic diversity among 43 spinach cultivars in China. The average polymorphism information content value of the 41 markers was 0.43, representing an intermediate level. The spinach cultivars had a low genetic diversity, and no detectable common factors were shared by each group in the UPGMA dendrogram. This study's findings facilitate further investigations on the organization of the microsatellites in spinach genome and provide clues for future breeding applications of spinach in China.