Genome-Wide Characterization of Simple Sequence Repeat (SSR) Loci in Chinese Jujube and Jujube SSR Primer Transferability

Comprehensive genome-wide identification and transferability of chromosome-specific highly variable microsatellite markers from citrus species

Citrus species among the most important and widely consumed fruit in the world due to Vitamin C, essential oil glands, and flavonoids. Highly variable simple sequence repeats (SSR) markers are one of the most informative and versatile molecular markers used in perennial tree genetic research. SSR survey of Citrus sinensis and Citrus maxima were identified perfect SSRs spanning nine chromosomes. Furthermore, we categorized all SSR motifs into three major classes based on their tract lengths. We designed and validated a class I SSRs in the C. sinensis and C. maxima genome through electronic polymerase chain reaction (ePCR) and found 83.89% in C. sinensis and 78.52% in C. maxima SSRs producing a single amplicon. Then, we selected extremely variable SSRs (> 40 nt) from the ePCR-verified class I SSRs and in silico validated across seven draft genomes of citrus, which provided us a subset of 84.74% in C. sinensis and 77.53% in C. maxima highly polymorphic SSRs. Out of these, 129 primers were validated on 24 citrus genotypes through wet-lab experiment. We found 127 (98.45%) polymorphic HvSSRs on 24 genotypes. The utility of the developed HvSSRs was demonstrated by analysing genetic diversity of 181 citrus genotypes using 17 HvSSRs spanning nine citrus chromosomes and were divided into 11 main groups through 17 HvSSRs. These chromosome-specific SSRs will serve as a powerful genomic tool used for future QTL mapping, molecular breeding, investigation of population genetic diversity, comparative mapping, and evolutionary studies among citrus and other relative genera/species.

Genotyping of Jujube (Ziziphus spp.) Germplasm in New Mexico and Southwestern Texas

Since the early 19th century, a substantial amount of jujube (Ziziphus spp.) germplasm has been introduced from China and Europe into the United States. However, due to a lack of passport data, cultivar mislabeling is common and the genetic background of the introduced germplasm remains unknown. In the present study, a low-density SNP array was employed to genotype 204 jujube trees sampled from multiple locations in New Mexico, Texas, Missouri, and Kentucky. Multilocus matching of SNP profiles revealed a significant rate of genetic redundancy among these jujube samples. A total of 14 synonymous groups were detected, comprising 48 accessions. Bayesian clustering analysis and neighbor-joining tree partitioned the US jujube germplasm into two major clusters. The first cluster included cultivated genotypes (Ziziphus jujuba Mill.), whereas the other major cluster comprised the wild/sour jujube (Ziziphus spinosa Hu.). The results also revealed a unique jujube population at Fabens/Tornillo, Texas, and a semi-naturalized population at Tucumcari, NM. These findings will provide valuable guidance to jujube growers and researchers on the effective utilization of jujube germplasm in the horticultural industry.

Transcriptomics and gas chromatography-mass spectrometry metabolomics reveal the mechanism of heat shock combined with 1-methylcyclopropene to regulate the cuticle wax of jujube fruit during storage

Cuticle wax is closely related to fruit quality during storage. In this study, changes in epidermal wax morphology, composition, and genes regulation induced by heat shock (HT), 1-methylcyclopropene (1-MCP) or their combination (HT + 1-MCP) were investigated in jujube fruit during cold storage. HT, 1-MCP, or HT + 1-MCP caused a smoother wax layer and fewer micro-cracks compared to the control (CK) during cold storage. It was confirmed that acids and terpenoids were the main wax components by gas chromatography-mass spectrometry. HT + 1-MCP and 1-MCP treatments could significantly increase (p < 0.05) the wax content at 45 d of cold storage. The transcriptomics results indicated that HT + 1-MCP treatment up-regulated FATB, FATB, FAB2, FAD2 and CYP716A, and maintained the wax content of jujube fruit during cold storage. These results could provide new perspective for regulating the cuticle characteristics to extend the shelf life of jujube fruit.

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.

Characterization of the Genetic Variability within Ziziphus nummularia Genotypes by Phenotypic Traits and SSR Markers with Special Reference to Geographic Distribution

Understanding the impacts and constraints of climate change on Ziziphus nummularia's geographical distribution is crucial for its future sustainability. In this study, we analyze information obtained from the field investigation, the distribution and response of climatic changes of Ziziphus nummularia by the use of ArcGIS analysis. The genetic diversity of 180 genotypes from three populations was studied by morphological attributes and simple sequence repeat (SSR). The results showed that the significant bioclimatic variable limiting the distribution of Z. nummularia was the mean temperature (bio 10_18.tif and bio19). Under the current climatic change, the suitable growth region of Z. nummularia is Swat (35.22° N, 72.42° E), while the future distribution would be Buner (34.39° N, 72.61° E), respectively. A total of 11 phenotypic traits were noted and had significant phenotypic variation among the traits. A total of 120 alleles were amplified. The alleles per locus ranged from 2 to 6, averaging 4.42, whereas PIC ranged from 0.33 to 0.79. Within a mean value of 0.67 per locus, expected heterozygosity was 0.57, observed heterozygosity was 0.661, and average gene diversity was 0.49. Flow estimates (6.41) indicated frequent gene flow within genotypes. The clustering, STRUCTURE, and PCoA analysis indicated Swat and Buner migration routes and evolution as well. The results indicated the prevalence of genetic variability and relationships among Z. nummularia across geographical boundaries had retained unique alleles. This may facilitate the development of agronomically desirable cultivars. However, climate change has impacted species distributions, requiring strategies to conserve genetic resources in different areas.

Genomic constitution and relationships of Ziziphus species collected from Malakand Division, Pakistan

Polyploidization has played a major role in plant evolution and can alter plant morphology, phenology, and ecology within only one or a few generations. Ziziphus species are economically as well as nutritionally important fruit-yielding trees. Identification of genotypes with unique traits or those with higher ploidy levels or a broader genetic base could lead to further improvements within the species. The current study has assessed the ploidy levels in the Ziziphus species (Ziziphus jujuba Mill. and Ziziphus nummularia (Burm. f. Wight & Arn) with phenotypic traits, flow cytometry, and chromosomal count as well as with SSRs markers. Morphological traits were inferred to be the most important drivers of trait variations among the investigated genotypes. The total sugar, total cAMPs, titratable acid, and chlorophyll (a, b, and total) were also significantly different in contrast with diploid plants, which showed that tetraploid Ziziphus had the potential to increase nutritional contents. Out of twenty (20), five (5) Z. jujuba genotypes (ZJL-9, ZJL-12, ZJL-17, ZJL-18, and ZJL-19) were found tetraploid 2n = 4x = 48, with genome size ranging from 965.9 to1238.8 Mb that was significantly higher than the tetraploid Z. jujuba Mill. variety Dongzao. Similarly, Z. nummularia ZNL-07 to ZNL-15 have found tetraploid 2n = 4x = 72 with genomic sizes ranging from 1152.2 to 1746.8 Mb respectively. Simple sequence repeats (SSRs) marker was applied to assess the genetic relationship within Ziziphus genotypes. To the best of our understanding, this is the first report on the identification of naturalized random tetraploids within the Pakistani Ziziphus species. This study provides important insights into the genomic architecture of Ziziphus species with implications for classification, conservation, and improvements of Ziziphus germplasm resources.

Genetic Diversity of Juglans mandshurica Populations in Northeast China Based on SSR Markers

Juglans mandshurica is a native tree species in Northeast China. Due to habitat destruction and human disturbance, its population size has sharply decreased. Currently, information on molecular markers of J. mandshurica is limited and cannot meet the needs of germplasm resource evaluation and molecular marker-assisted breeding of J. mandshurica. Based on transcriptomic data from three tissues (leaves, bark, and fruit pericarp), we developed expressed sequence tag-simple sequence repeats (EST-SSRs) for J. mandshurica, and 15 polymorphic EST-SSR primers were initially selected. The average number of alleles (Na), expected heterozygosity (He), and the polymorphic information content (PIC) at different loci were 18.27, 0.670, and 0.797, respectively. Population genetic diversity analysis revealed that the average Na, He, and Shannon information indices (I) for 15 J. mandshurica populations were 6.993, 0.670, and 1.455, respectively. Among them, population Hunchun exhibited the highest genetic diversity (Na = 7.933, He = 0.723, and I = 1.617), while population Heihe exhibited the lowest genetic diversity (Na = 4.200, He = 0.605, and I = 1.158). STRUCTURE analysis, neighbor-joining method cluster analysis, and principal coordinate analysis showed that the 343 individuals of J. mandshurica from 15 populations were clustered into three categories. Category 1 (green) had 147 individuals from eight populations in Qingyuan, Caohekou, Jian, Ningan, Yongji, Baishishan, Helong, and Maoershan; category 2 (blue) had 81 individuals from three populations in Hulin, Boli, and Sanchazi; and category 3 (red) had 115 individuals from four populations in Heihe, Hunchun, Fangzheng, and Liangshui. Analysis of molecular variance (AMOVA) showed that genetic variations among and within individuals accounted for 16.22% and 21.10% of the total genetic variation, respectively, indicating that genetic variations within populations were greater than genetic variations among populations. The average genetic differentiation coefficient (Fst) and gene flow (Nm) between different populations were 0.109 and 4.063, respectively, implying moderate levels of genetic differentiation and gene flow. Based on the genetic diversity characteristics of different populations, we proposed various genetic conservation strategies for J. mandshurica.

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis

BACKGROUND

Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections.

RESULTS

All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny.

CONCLUSIONS

Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

Heat shock treatment maintains the quality attributes of postharvest jujube fruits and delays their senescence process during cold storage

The effects of heat shock (HT), 1-methylcyclopropene (1-MCP), or their combination (HT + 1-MCP) on the quality of fresh jujube fruits during cold storage were studied. Among them, HT showed the best preservation effect on jujube fruits, which was more effective than others in inhibiting the increase of red index, decay incidence, and weight loss and delaying the decrease of firmness, soluble solids content (SSC), titratable acidity (TA), and ascorbic acid (AsA) content. Besides, it could delay the degradation rate of the cell wall to maintain the integrity of cell membrane, and keep the high activity of active oxygen scavenging enzymes. During cold storage, malondialdehyde (MDA) content and relative electrolyte leakage (REL) of the HT group were significantly lower than those of the control group, 1-MCP, and HT + 1-MCP group (p < .05), while superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) activities were significantly higher than those of other groups (p < .05). It was concluded that the postharvest HT treatment could effectively delay the senescence and decay of jujube fruits. PRACTICAL APPLICATIONS: Jujube fruits have high nutritional value used for food and medicine. However, they are not tolerant to storage after harvest, resulting in high economic losses. Therefore, it is of great significance to find a suitable method to maintain the quality of jujube fruits. Our results revealed the effect of HT, 1-MCP, and their combination on the quality maintenance of jujube fruits, and found that HT could effectively maintain the quality of them, which could be used as an effective method for keeping jujube fruits fresh.

The transcriptomic response to heat stress of a jujube (Ziziphus jujuba Mill.) cultivar is featured with changed expression of long noncoding RNAs

Long non-coding RNA (lncRNA) of plant species undergoes dynamic regulation and acts in developmental and stress regulation. Presently, there is little information regarding the identification of lncRNAs in jujube (Ziziphus jujuba Mill.), and it is uncertain whether the lncRNAs could respond to heat stress (HS) or not. In our previous study, a cultivar (Hqing1-HR) of Z. jujuba were treated by HS (45°C) for 0, 1, 3, 5 and 7 days, and it was found that HS globally changed the gene expression by RNA sequencing (RNA-seq) experiments and informatics analyses. In the current study, 8260 lncRNAs were identified successfully from the previous RNA-seq data, and it indicated that lncRNAs expression was also altered globally, suggesting that the lncRNAs might play vital roles in response to HS. Furthermore, bioinformatics analyses of potential target mRNAs of lncRNAs with cis-acting mechanism were performed, and it showed that multiple differentially expressed (DE) mRNAs co-located with DElncRNAs were highly enriched in pathways associated with response to stress and regulation of metabolic process. Taken together, these findings not only provide a comprehensive identification of lncRNAs but also useful clues for molecular mechanism response to HS in jujube.

Genome survey of Zanthoxylum bungeanum and development of genomic-SSR markers in congeneric species

Zanthoxylum bungeanum, a spice and medicinal plant, is cultivated in many parts of China and some countries in Southeast Asia; however, data on its genome are lacking. In the present study, we performed a whole-genome survey and developed novel genomic-SSR markers of Z. bungeanum. Clean data (∼197.16 Gb) were obtained and assembled into 11185221 scaffolds with an N50 of 183 bp. K-mer analysis revealed that Z. bungeanum has an estimated genome size of 3971.92 Mb, and the GC content, heterozygous rate, and repeat sequence rate are 37.21%, 1.73%, and 86.04%, respectively. These results indicate that the genome of Z. bungeanum is complex. Furthermore, 27153 simple sequence repeat (SSR) loci were identified from 57288 scaffolds with a minimum length > 1 kb. Mononucleotide repeats (19706) were the most abundant type, followed by dinucleotide repeats (5154). The most common motifs were A/T, followed by AT/AT; these SSRs accounted for 71.42% and 11.84% of all repeats, respectively. A total of 21243 non-repeating primer pairs were designed, and 100 were randomly selected and validated by PCR analysis using DNA from 10 Z. bungeanum individuals and 5 Zanthoxylum armatum individuals. Finally, 36 polymorphic SSR markers were developed with polymorphism information content (PIC) values ranging from 0.16 to 0.75. Cluster analysis revealed that Z. bungeanum and Z. armatum could be divided into two major clusters, suggesting that these newly developed SSR markers are useful for genetic diversity and germplasm resource identification in Z. bungeanum and Z. armatum.

Genomic survey sequencing, development and characterization of single- and multi-locus genomic SSR markers of Elymus sibiricus L

BACKGROUND

Siberian wildrye (Elymus sibiricus L.) attracts considerable interest for grassland establishment and pasture recovery in the Qinghai-Tibet Plateau (QTP) due to its excellence in strong stress tolerance, high nutritional value and ease to cultivate. However, the lack of genomic information of E. sibiricus hampers its genetics study and breeding process.

RESULTS

In this study, we performed a genome survey and developed a set of SSR markers for E. sibiricus based on Next-generation sequencing (NGS). We generated 469.17 Gb clean sequence which is 58.64× of the 6.86 Gb estimated genome size. We assembled a draft genome of 4.34 Gb which has 73.23% repetitive elements, a heterozygosity ratio of 0.01% and GC content of 45.68%. Based on the gnomic sequences we identified 67,833 SSR loci and from which four hundred were randomly selected to develop markers. Finally, 30 markers exhibited polymorphism between accessions and ten were identified as single-locus SSR. These newly developed markers along with previously reported 30 ones were applied to analyze genetic polymorphism among 27 wild E. sibiricus accessions. We found that single-locus SSRs are superior to multi-loci SSRs in effectiveness.

CONCLUSIONS

This study provided insights into further whole genome sequencing of E. sibiricus in strategy selection. The novel developed SSR markers will facilitate genetics study and breeding for Elymus species.

Genome survey sequencing and genetic diversity of cultivated Akebia trifoliata assessed via phenotypes and SSR markers

Akebia trifoliata (Lardizabalaceae) is an important medicinal plant with multiple pharmacological effects. However, the lack of genomic information had limited the further excavation and utilization of this plant. An initial survey of the genome A. trifoliata was performed by next-generation sequencing, and then the genome size was inferred by flow cytometry. The whole genome survey of A. trifoliata generated 61.90 Gb of sequence data with approximately 95.51 × coverage. The genome size, heterozygosity and GC content obtained by k-mer analysis were almost 648.07 Mb, 0.72% and 36.11%, respectively. The genome size calculated by flow cytometry was 685.77 Mb, which was consistent with the results of genome survey. A total of 851,957 simple sequence repeats (SSR) were identified in the A. trifoliata genome. Twenty-eight phenotypic traits and thirty pairs of SSR primers were selected for the analysis of the genetic diversity of 43 accessions of cultivated A. trifoliata. The results showed that 216 bands were generated by 30 pairs of SSR primers, of which 189 (87.5%) were polymorphic. In addition, the phenotypes and SSR markers were used for cluster analysis of 43 cultivated accessions. The results of the two clustering methods were partially consistent. The genome survey of A. trifoliata demonstrated that the genome size of this plant was about 648.07 Mb. In the present study, the size and characteristics of the genome of A. trifoliata were reported for the first time, which greatly enriched the genomic resources of A. trifoliata for the further research and utilization.

Genome Survey Sequencing of In Vivo Mother Plant and In Vitro Plantlets of Mikania cordata

Mikania cordata, the only native congener of the invasive weed Mikania micrantha in China, is an ideal species for comparative study to reveal the invasion mechanism. However, its genome resources are lagging far behind its congener, which limits the comparative genomic analysis. Our goal is to characterize the genome of M. cordata by next-generation sequencing and propose a scheme for long-read genome sequencing. Previous studies have shown that the genomic resources of the host plant would be affected by the endophytic microbial DNA. An aseptic sample of M. cordata will ensure the proper genome in downstream analysis. Because endophytes are ubiquitous in the greenhouse-grown M. cordata, the in vitro culture with cefotaxime or timentin treatment was undertaken to obtain the aseptic plantlets. The in vivo mother plant and in vitro plantlets were used to survey the genome. The microbial contamination in M. cordata was recognized by blast search and eliminated from the raw reads. The decontaminated sequencing reads were used to predict the genome size, heterozygosity, and repetitive rate. The in vivo plant was so contaminated that microbes occupied substantial sequencing resources and misled the scaffold assembly. Compared with cefotaxime, treatment with timentin performed better in cultivating robust in vitro plantlets. The survey result from the in vitro plantlets was more accurate due to low levels of contamination. The genome size was estimated to be 1.80 Gb with 0.50% heterozygosity and 78.35% repetitive rate. Additionally, 289,831 SSRs were identified in the genome. The genome is heavily contaminated and repetitive; therefore, the in vitro culture technique and long-read sequencing technology are recommended to generate a high-quality and highly contiguous genome.

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

Simple Summary

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

Abstract

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

The historical and current research progress on jujube-a superfruit for the future

Jujube (Ziziphus jujuba Mill.), or Chinese date, is the most important species of Rhamnaceae, a large cosmopolitan family, and is one of the oldest cultivated fruit trees in the world. It originates from the middle and lower reaches of the Yellow River, the 'mother river' of the Chinese people. It is distributed in at least 48 countries on all continents except Antarctica and is becoming increasingly important, especially in arid and semiarid marginal lands. Based on a systematic analysis of the unique characteristics of jujube, we suggest that it deserves to be recognized as a superfruit. We summarized historical research achievements from the past 3000 years and reviewed recent research advances since 1949 in seven fields, including genome sequencing and application, germplasm resources and systematic taxonomy, breeding and genetics, cultivation theory and techniques, pest control, postharvest physiology and techniques, and nutrition and processing. Based on the challenges facing the jujube industry, we discuss eight research aspects to be focused on in the future.

Development and Application of SSR Markers Related to Genes Involved in Leaf Adaxial-Abaxial Polarity Establishment in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

In Chinese cabbage (Brassica rapa L. ssp. pekinensis), leaf adaxial-abaxial (ad-ab) polarity is tightly related to leaf incurvature, an essential factor for the formation of leafy heads. Therefore, identification of the genes responsible for leaf ad-ab polarity and studying their genetic variation may clarify the mechanism of leafy head formation. By comparing the sequences of the genes regulating leaf ad-ab polarity development in Arabidopsis thaliana (A. thaliana), 41 candidate genes distributed on 10 chromosomes were found to be responsible for the establishment of ad-ab polarity in Chinese cabbage. Orthologous genes, including 10 single copies, 14 double copies, and one triple copies, were detected in the Chinese cabbage. The gene structure and conserved domain analyses showed that the number of exons of the 41 candidate genes range from one to 25, and that most genes share the conserved motifs 1, 6, and 10. Based on the 41 candidate genes, 341 simple sequence repeats (SSRs) were detected, including five replicated types: single, double, triple, quintuple, and sextuple nucleotide replications. Among these sequence repeat (SSR) loci, 323 loci were used to design 969 specific primers, and 362 primer pairs were selected randomly and evaluated using 12 Chinese cabbage accessions with different heading types. 23 primer pairs resulting with clear, polymorphic bands, combined with other 127 markers, was used to construct a linkage map by using an F₂ population containing 214 lines derived from the hybrid of the overlapping heading Chinese cabbage “14Q-141” and the outward curling heading Chinese cabbage “14Q-279.” The result showed that the sequences of markers in the genetic linkage map and the physical map was consistent in general. Our study could help to accelerate the breeding process of leafy head quality in Chinese cabbage.

Characterization and Development of Genomic SSRs in Pecan (Carya illinoinensis)

Research Highlights: The distribution of simple sequence repeat (SSR) motifs in two draft genomes of pecan was evaluated. Sixty-six SSR loci were validated by PCR amplification in pecan. Twenty-two new development markers can be used for genetic study in genus Carya. Background and Objectives: Pecan has good nutritional and health benefits and is an important crop worldwide. However, the genetic research in this species is insufficient. One of the main reasons for this is the lack of enough accurate, convenient, and economical molecular markers. Among different marker types, SSR loci are enormously useful in genetic studies. However, the number of SSRs in C. illinoinensis (Wangenh.) K. Koch is limited. Materials and Methods: The distribution of SSR motifs in the pecan genome was analyzed. Then, the primers for each SSR were designed. To evaluate their availability, 74 SSR loci were randomly selected and amplified in pecan. Finally, 22 new SSRs and eight former ones were picked to evaluate the genetic diversity in 60 pecan genotypes and to determine their transferability in other Carya species. Results: 145,714 and 143,041 SSR motifs were obtained from two draft genomes of ‘87MX3-2’ and ‘Pawnee’, respectively. In total, 9145 candidate primers were obtained. Sixty-six (89.19%) primers amplified the target products. Among the 30 SSRs, 29 loci showed polymorphism in 60 pecan genotypes. The polymorphic information content (PIC) values ranged from 0.012 to 0.906. In total, 26, 25, and 22 SSRs can be used in C. cathayensis Sarg., C. dabieshanensis W. C. Cheng & R. H. Chang, and C. hunanensis W.C. Liu, respectively. Finally, the dendrogram of all individuals was constructed. The results agree with the geographic origin of the four species and the pedigree relationships between different pecan cultivars. Conclusions: The characterization of SSRs in the pecan genome and the new SSRs will promote the progress of genetic study and breeding in pecan, as well as other species of genus Carya.

Detection and application of genome-wide variations in peach for association and genetic relationship analysis

BACKGROUND

Peach (Prunus persica L.) is a diploid species and model plant of the Rosaceae family. In the past decade, significant progress has been made in peach genetic research via DNA markers, but the number of these markers remains limited.

RESULTS

In this study, we performed a genome-wide DNA markers detection based on sequencing data of six distantly related peach accessions. A total of 650,693~1,053,547 single nucleotide polymorphisms (SNPs), 114,227~178,968 small insertion/deletions (InDels), 8386~12,298 structure variants (SVs), 2111~2581 copy number variants (CNVs) and 229,357~346,940 simple sequence repeats (SSRs) were detected and annotated. To demonstrate the application of DNA markers, 944 SNPs were filtered for association study of fruit ripening time and 15 highly polymorphic SSRs were selected to analyze the genetic relationship among 221 accessions.

CONCLUSIONS

The results showed that the use of high-throughput sequencing to develop DNA markers is fast and effective. Comprehensive identification of DNA markers, including SVs and SSRs, would be of benefit to genetic diversity evaluation, genetic mapping, and molecular breeding of peach.

Genome-wide characterization and analysis of microsatellite sequences in camelid species

Microsatellites or simple sequence repeats (SSRs) are among the genetic markers most widely utilized in research. This includes applications in numerous fields such as genetic conservation, paternity testing, and molecular breeding. Though ordered draft genome assemblies of camels have been announced, including for the Arabian camel, systemic analysis of camel SSRs is still limited. The identification and development of informative and robust molecular SSR markers are essential for marker assisted breeding programs and paternity testing. Here we searched and compared perfect SSRs with 1–6 bp nucleotide motifs to characterize microsatellites for draft genome sequences of the Camelidae. We analyzed and compared the occurrence, relative abundance, relative density, and guanine-cytosine (GC) content in four taxonomically different camelid species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. A total of 546762, 544494, 547974, and 437815 SSRs were mined, respectively. Mononucleotide SSRs were the most frequent in the four genomes, followed in descending order by di-, tetra-, tri-, penta-, and hexanucleotide SSRs. GC content was highest in dinucleotide SSRs and lowest in mononucleotide SSRs. Our results provide further evidence that SSRs are more abundant in noncoding regions than in coding regions. Similar distributions of microsatellites were found in all four species, which indicates that the pattern of microsatellites is conserved in family Camelidae.

Comparative analyses of simple sequence repeats (SSRs) in 23 mosquito species genomes: Identification, characterization and distribution (Diptera: Culicidae)

Simple sequence repeats (SSRs) exist in both eukaryotic and prokaryotic genomes and are the most popular genetic markers, but the SSRs of mosquito genomes are still not well understood. In this study, we identified and analyzed the SSRs in 23 mosquito species using Drosophila melanogaster as reference at the whole-genome level. The results show that SSR numbers (33 076-560 175/genome) and genome sizes (574.57-1342.21 Mb) are significantly positively correlated (R² = 0.8992, P < 0.01), but the correlation in individual species varies in these mosquito species. In six types of SSR, mono- to trinucleotide SSRs are dominant with cumulative percentages of 95.14%-99.00% and densities of 195.65/Mb-787.51/Mb, whereas tetra- to hexanucleotide SSRs are rare with 1.12%-4.22% and 3.76/Mb-40.23/Mb. The (A/T)n, (AC/GT)n and (AGC/GCT)n are the most frequent motifs in mononucleotide, dinucleotide and trinucleotide SSRs, respectively, and the motif frequencies of tetra- to hexanucleotide SSRs appear to be species-specific. The 10-20 bp length of SSRs are dominant with the number of 110 561 ± 93 482 and the frequency of 87.25% ± 5.73% on average, and the number and frequency decline with the increase of length. Most SSRs (83.34% ± 7.72%) are located in intergenic regions, followed by intron regions (11.59% ± 5.59%), exon regions (3.74% ± 1.95%), and untranslated regions (1.32% ± 1.39%). The mono-, di- and trinucleotide SSRs are the main SSRs in both gene regions (98.55% ± 0.85%) and exon regions (99.27% ± 0.52%). An average of 42.52% of total genes contains SSRs, and the preference for SSR occurrence in different gene subcategories are species-specific. The study provides useful insights into the SSR diversity, characteristics and distribution in 23 mosquito species of genomes.

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.

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.

Variability and population genetic structure in Achyrocline flaccida (Weinm.) DC., a species with high value in folk medicine in South America

Better knowledge of medicinal plant species and their conservation is an urgent need worldwide. Decision making for conservation strategies can be based on the knowledge of the variability and population genetic structure of the species and on the events that may influence these genetic parameters. Achyrocline flaccida (Weinm.) DC. is a native plant from the grassy fields of South America with high value in folk medicine. In spite of its importance, no genetic and conservation studies are available for the species. In this work, microsatellite and ISSR (inter-simple sequence repeat) markers were used to estimate the genetic variability and structure of seven populations of A. flaccida from southern Brazil. The microsatellite markers were inefficient in A. flaccida owing to a high number of null alleles. After the evaluation of 42 ISSR primers on one population, 10 were selected for further analysis of seven A. flaccida populations. The results of ISSR showed that the high number of exclusive absence of loci might contribute to the inter-population differentiation. Genetic variability of the species was high (Nei’s diversity of 0.23 and Shannon diversity of 0.37). AMOVA indicated higher genetic variability within (64.7%) than among (33.96%) populations, and the variability was unevenly distributed (F_ST 0.33). Gene flow among populations ranged from 1.68 to 5.2 migrants per generation, with an average of 1.39. The results of PCoA and Bayesian analyses corroborated and indicated that the populations are structured. The observed genetic variability and population structure of A. flaccida are discussed in the context of the vegetation formation history in southern Brazil, as well as the possible anthropogenic effects. Additionally, we discuss the implications of the results in the conservation of the species.

Genome-Wide Organization and Expression Profiling of the SBP-Box Gene Family in Chinese Jujube (Ziziphus jujuba Mill.)

Transcription factors play vital roles in the developmental processes of plants. The SQUAMOSA promoter binding protein (SBP) genes encode a family of plant-specific transcription factors and plays many crucial roles in plant development. In this study, 16 SBP-box gene family members were identified in Ziziphus jujuba Mill. Dongzao (Dongzao), which were distributed over 8 chromosomes. They were classified into seven groups according to their phylogenetic relationships with other SBP-box gene families. Within each group, genes shared similar exon-intron structures and motif locations. The number of exons varied among the groups. We identified 12 homologous gene pairs between Dongzao and Arabidopsis. Expression profiling revealed that ZjSBP02 and ZjSBP14 expressed highly in mature fruits, ZjSBP01 expressed higher in mature leaves than other tissues and the expression level of ZjSBP12 was much higher in the flowers. The transcriptome analysis indicated that ZjSBPs had different expression patterns in various tissues. This study represents the first systematic analysis of the SBP-box gene family in Z. jujuba. The data presented here provides a foundation for understanding the crucial roles of ZjSBP genes in plant growth and development.

Genetic Diversity, Population Structure, and Linkage Disequilibrium of a Core Collection of Ziziphus jujuba Assessed with Genome-wide SNPs Developed by Genotyping-by-sequencing and SSR Markers

Chinese jujube (Ziziphus jujuba Mill) is an economically important fruit species native to China with high nutritious and medicinal value. Genotyping-by-sequencing was used to detect and genotype single nucleotide polymorphisms (SNPs) in a core collection of 150 Chinese jujube accessions and further to characterize their genetic diversity, population structure, and linkage disequilibrium (LD). A total of 4,680 high-quality SNPs were identified, of which 38 sets of tri-allelic SNPs were detected. The average polymorphism information content (PIC) values based on bi-allelic SNPs and tri-allelic SNPs were 0.27 and 0.38, respectively. STRUCTURE and principal coordinate analyses based on SNPs revealed that the 150 accessions could be clustered into two groups. However, neighbor-joining trees indicated the accessions should be grouped into three major clusters. Our data confirm that the resolving power for genetic diversity was similar for the SSRs and SNPs. In contrast, regarding population structure, the resolving power was higher for SSRs than for SNPs. The LD pattern in Chinese jujube was investigated for the first time. We observed a relatively rapid LD decay with a short range (∼10 kb) for all pseudo-chromosomes and for individual pseudo-chromosomes. Our findings provide important information for future genome-wide association analyses and marker-assisted selective breeding of Chinese jujube.

Genome survey of pistachio (Pistacia vera L.) by next generation sequencing: Development of novel SSR markers and genetic diversity in Pistacia species

BACKGROUND

Pistachio (Pistacia vera L.) is one of the most important nut crops in the world. There are about 11 wild species in the genus Pistacia, and they have importance as rootstock seed sources for cultivated P. vera and forest trees. Published information on the pistachio genome is limited. Therefore, a genome survey is necessary to obtain knowledge on the genome structure of pistachio by next generation sequencing. Simple sequence repeat (SSR) markers are useful tools for germplasm characterization, genetic diversity analysis, and genetic linkage mapping, and may help to elucidate genetic relationships among pistachio cultivars and species.

RESULTS

To explore the genome structure of pistachio, a genome survey was performed using the Illumina platform at approximately 40× coverage depth in the P. vera cv. Siirt. The K-mer analysis indicated that pistachio has a genome that is about 600 Mb in size and is highly heterozygous. The assembly of 26.77 Gb Illumina data produced 27,069 scaffolds at N50 = 3.4 kb with a total of 513.5 Mb. A total of 59,280 SSR motifs were detected with a frequency of 8.67 kb. A total of 206 SSRs were used to characterize 24 P. vera cultivars and 20 wild Pistacia genotypes (four genotypes from each five wild Pistacia species) belonging to P. atlantica, P. integerrima, P. chinenesis, P. terebinthus, and P. lentiscus genotypes. Overall 135 SSR loci amplified in all 44 cultivars and genotypes, 41 were polymorphic in six Pistacia species. The novel SSR loci developed from cultivated pistachio were highly transferable to wild Pistacia species.

CONCLUSIONS

The results from a genome survey of pistachio suggest that the genome size of pistachio is about 600 Mb with a high heterozygosity rate. This information will help to design whole genome sequencing strategies for pistachio. The newly developed novel polymorphic SSRs in this study may help germplasm characterization, genetic diversity, and genetic linkage mapping studies in the genus Pistacia.

Identifying the genetic diversity, genetic structure and a core collection of Ziziphus jujuba Mill. var. jujuba accessions using microsatellite markers

Ziziphus is a genus of spiny shrubs and small trees in the Rhamnaceae family. This group has a controversial taxonomy, with more than 200 species described, including Chinese jujube (Ziziphus jujuba Mill. var. jujuba) and Indian jujube (Z. mauritiana), as well as several other important cultivated fruit crops. Using 24 SSR markers distributed across the Chinese jujube genome, 962 jujube accessions from the two largest germplasm repositories were genotyped with the aim of analyzing the genetic diversity and structure and constructing a core collection that retain high genetic diversity. A molecular profile comparison revealed 622 unique genotypes, among which 123 genotypes were genetically identical to at least one other accessions. STRUCTURE analysis and multivariate analyses (Cluster and PCoA) roughly divided the accessions into three major groups, with some admixture among groups. A simulated annealing algorithm and a heuristic algorithm were chosen to construct the core collection. A final core of 150 accessions was selected, comprising 15.6% of the analyzed accessions and retaining more than 99.5% of the total alleles detected. We found no significant differences in allele frequency distributions or in genetic diversity parameters between the chosen core accessions and the 622 genetically unique accessions. This work contributes to the understanding of Chinese jujube diversification and the protection of important germplasm resources.

Identification of SNP and SSR Markers in Finger Millet Using Next Generation Sequencing Technologies

Finger millet is an important cereal crop in eastern Africa and southern India with excellent grain storage quality and unique ability to thrive in extreme environmental conditions. Since negligible attention has been paid to improving this crop to date, the current study used Next Generation Sequencing (NGS) technologies to develop both Simple Sequence Repeat (SSR) and Single Nucleotide Polymorphism (SNP) markers. Genomic DNA from cultivated finger millet genotypes KNE755 and KNE796 was sequenced using both Roche 454 and Illumina technologies. Non-organelle sequencing reads were assembled into 207 Mbp representing approximately 13% of the finger millet genome. We identified 10,327 SSRs and 23,285 non-homeologous SNPs and tested 101 of each for polymorphism across a diverse set of wild and cultivated finger millet germplasm. For the 49 polymorphic SSRs, the mean polymorphism information content (PIC) was 0.42, ranging from 0.16 to 0.77. We also validated 92 SNP markers, 80 of which were polymorphic with a mean PIC of 0.29 across 30 wild and 59 cultivated accessions. Seventy-six of the 80 SNPs were polymorphic across 30 wild germplasm with a mean PIC of 0.30 while only 22 of the SNP markers showed polymorphism among the 59 cultivated accessions with an average PIC value of 0.15. Genetic diversity analysis using the polymorphic SNP markers revealed two major clusters; one of wild and another of cultivated accessions. Detailed STRUCTURE analysis confirmed this grouping pattern and further revealed 2 sub-populations within wild E. coracana subsp. africana. Both STRUCTURE and genetic diversity analysis assisted with the correct identification of the new germplasm collections. These polymorphic SSR and SNP markers are a significant addition to the existing 82 published SSRs, especially with regard to the previously reported low polymorphism levels in finger millet. Our results also reveal an unexploited finger millet genetic resource that can be included in the regional breeding programs in order to efficiently optimize productivity.

Characterization of SSR genomic abundance and identification of SSR markers for population genetics in Chinese jujube (Ziziphus jujuba Mill.)

Chinese jujube (Ziziphus jujuba Mill. [Rhamnaceae]), native to China, is a major dried fruit crop in Asia. Although many simple sequence repeat (SSR) markers are available for phylogenetic analysis of jujube cultivars, few of these are validated on the level of jujube populations. In this study, we first examined the abundance of jujube SSRs with repeated unit lengths of 1–6 base pairs, and compared their distribution with those in Arabidopsis thaliana. We identified 280,596 SSRs in the assembled genome of jujube. The density of SSRs in jujube was 872.60 loci/Mb, which was much higher than in A. thaliana (221.78 loci/Mb). (A+ T)-rich repeats were dominant in the jujube genome. We then randomly selected 100 SSRs in the jujube genome with long repeats and used them to successfully design 70 primer pairs. After screening using a series of criteria, a set of 20 fluorescently labeled primer pairs was further selected and screened for polymorphisms among three jujube populations. The average number of alleles per locus was 12.8. Among the three populations, mean observed and expected heterozygosities ranged from 0.858 to 0.967 and 0.578 to 0.844, respectively. After testing in three populations, all SSRs loci were in Hardy-Weinberg equilibrium (HWE) in at least one population. Finally, removing high null allele frequency loci and linked loci, a set of 17 unlinked loci was in HWE. These markers will facilitate the study of jujube genetic structure and help elucidate the evolutionary history of this important fruit crop.