Rapid Genome-Wide Location-Specific Polymorphic SSR Marker Discovery in Black Pepper by GBS Approach

Black pepper (Piper nigrum), the “King of Spices,” is an economically important spice in India and is known for its medicinal and cultural values. SSRs, the tandem repeats of small DNA sequences, are often polymorphic in nature with diverse applications. For population structure, QTL/gene discovery, MAS, and diversity analysis, it is imperative to have their location specificity. The existing PinigSSRdb catalogs ~70K putative SSR markers but these are anonymous (unknown chromosomal location), based on 916 scaffolds rather than 26 chromosomes. Under this study, we generated ddRAD sequence data of 29 black pepper genotypes from all over India, being low-cost and most efficient technique for the identification of polymorphic markers. The major limitation of ddRAD with compromised/non-uniform coverage has been successfully overcome by taking advantage of chromosome-wise data availability. The latest black pepper genome assembly was used to extract genome-wide SSRs. A total of 276,230 genomic SSRs were mined distributed over 26 chromosomes, with relative density of 362.88 SSRs/Mb and average distance of 2.76 Kb between two SSRs. This assembly was also used to find the polymorphic SSRs in the generated GBS data of 29 black pepper genotypes utilizing rapid and cost-effective method giving 3,176 polymorphic SSRs, out of which 2015 were found to be hypervariable. The developed web-genomic resource, BlackP2MSATdb (http://webtom.cabgrid.res.in/blackp2msatdb/), is the largest and first reported web resource for genomic and polymorphic SSRs of black pepper, which is useful to develop varietal signature, coreset, physical map, QTL/gene identification, and MAS in endeavor of black pepper production.

Introduction

Black pepper, one of the oldest spices globally used, belongs to the family Piperaceae (scientific name: Piper nigrum; 2n =52). It is widely known as “King of Spices” and is economically important spice in India, hence also termed as “Black Gold.” It is valued for its characteristic pungency and flavor (which is due to the alkaloid, piperine) along with many health benefits. It is perennial woody vine that gained its world-wide economic and cultural potential due to which its production and trade widened globally. Black pepper is known to have its origin in the Western Ghats of Southern India due to the humid climate, but now, it is also grown in tropical and sub-tropical regions of countries like Ethiopia, Vietnam, Brazil, Indonesia, India, Tajikistan, Sri Lanka, China, Malaysia, and Mexico. It is cultivated on total area of 259 thousand hectares with production of 65 thousand tons in India (Spice Board of India, Ministry of Agriculture and Farmers Welfare, Govt. of India), while globally it is cultivated on 749 thousand hectares with production of 1,103 thousand tons (FAO, 2019).

Microsatellites or simple sequence repeats (SSRs) are tandem repeats of small DNA sequences, often polymorphic in nature which serve in genetic analysis, mapping, gene tagging, and marker-assisted selection of linked desired economically important traits. Several attempts have been made in previous studies in development and characterization of black pepper microsatellite markers (Raghavan et al., 2010; Joy et al., 2011; Gordo et al., 2012; Hu et al., 2015; Jagtap et al., 2016; Wu et al., 2016; Jose et al., 2017) using various data sets including RNA-seq data. The first whole-genome SSR mining in black pepper was performed by Kumari et al. (2019) and presented in the form of a database, PinigSSRdb. This very first black pepper SSR database contains ~70 thousand putative SSR markers, but all are anonymous loci (unknown chromosomal location) as they are based on 916 scaffolds rather than 26 chromosomes. For population structure and diversity analysis, it is imperative to have SSR loci with known physical location (two–three SSRs per chromosome) having polymorphism information content (PIC) value >0.6 (four−12 alleles). The number of SSR loci can be drastically minimized, if markers are selected uniformly by taking the advantage of their known physical location over chromosome. Such selected markers can be a better diversity calculus for genomic and population variability with respect to random anonymous SSR loci (You et al., 2005). These chromosomes and location-specific SSR loci can be further used to trace the inheritance of particular chromosomal regions in molecular breeding programs from the foundation genotypes (Pestsova and Röder, 2002). Hence, the role of SSR markers in black pepper is immense including the study of its diversity, development of the germplasm core set, removal of duplicates, development of varietal signature, and low-cost polymorphic markers.

The latest genome assembly of black pepper provided by Hu et al. (2019) can be used prudently to cater the above critical needs. We used this latest assembly of black pepper to extract whole-genome SSRs with improved characterization and SSR-based assessment of genetic diversity among 29 genotypes of black pepper in India utilizing a rapid and cost-effective method to extract polymorphic SSRs. This was supplemented with the development of Black Pepper Polymorphic Microsatellite Database: BlackP2MSATdb, which provides the public access to the largest and first database of genomic and polymorphic SSRs of black pepper to be utilized by scientific community in the studies involving molecular breeding along with DNA fingerprinting, varietal identification, construction of physical map, identification of QTL, map-based gene cloning, marker-assisted selection, and evolutionary studies in black pepper and its related species.

Materials and Methods

Genomic DNA Extraction and GBS Data Generation

A total of 29 black pepper genotypes (released varieties, important germplasm accessions, and wild relatives) were maintained at Indian Council of Agricultural Research-Indian Institute of Spices Research, Kozhikode, Kerala, India (11.2986° N, 75.8405° E) (Supplementary File 1). The leaf samples were collected for DNA isolation from these 29 genotypes. The ddRAD genotyping-by-sequencing (GBS) data for these 29 genotypes of black pepper were generated for this study. The leaf samples were kept in liquid nitrogen. Total genomic DNA was isolated using CTAB extraction method (Saghai-Maroof MA). Quality and quantity checks of DNA were done by Qubit Fluorometer (Thermo Fisher Scientific Inc.). For determining A260/280 ratio, 1 μl of sample was loaded in Nanodrop (Thermo Fisher Scientific Inc.). The quality of samples was checked in 0.8% agarose gel electrophoresis along with Hind III marker for the presence of intact bands. In gDNA (100 ng), barcoded adapters (~0.1 pM) were added and double-digested with EcoRI and MseI enzyme for 2 h at 37 and 65°C for 20 min. The barcoded and common adapters were ligated to the sticky ends of the digested DNA with T4 DNA Ligase enzyme at 16°C overnight. The ligated products were pooled together. The pool of ligated products was size-selected to 150–650 bp using ampure beads. The pool was PCR-amplified to generate the final library. The library pool was analyzed in Bioanalyzer 2100 (Agilent Technologies) using high-sensitivity (HS) DNA chip as per manufacturer's instructions. The library pool was sequenced on Illumina NextSeq platform. Finally, the 29 RAD-based genotyping-by-sequencing (GBS) data libraries with pair-end reads were separated from library pool and passed through quality check.

Quality Check and de novo Assembly of GBS Data Libraries

All the GBS libraries of 29 genotypes were passed through quality check using FastQC (Wingett and Andrews, 2018) keeping the default parameters, viz., Phred score ≥30 and GC distribution >40%. De novo assembly was performed for each library using SPAdes-3.13.0 (Prjibelski et al., 2020) with default options. The de novo assembly for 29 samples of generated GBS data was done to find the polymorphic SSRs using the black pepper genome as reference (Hu et al., 2019). Simply mapping was not done to avoid missing of many SSRs from genome that are not represented in RefSeq.

Genome-Wide Mining of SSR Markers and Their Primer Designing

In order to mine location-specific, genome-wide SSRs, black pepper reference genome was downloaded from http://cotton.hzau.edu.cn/EN/download.php (Hu et al., 2019). The black pepper genome has a total of 45 scaffolds, out of which 26 were considered as pseudo-chromosomes (Pn1–Pn26), while the remaining 19 scaffolds have been termed as PnU throughout the manuscript.

The SSRs were mined from the black pepper reference genome (Hu et al., 2019) and from GBS libraries of 29 genotypes using the perl script of MISA (MIcroSAtellite identification) tool (Beier et al., 2017). The criteria utilized for the identification were as follows: mononucleotide repeats motif with at least 10 repeats, dinucleotide with six, tri-, tetra-, penta- and hexa-nucleotide with five repeats (Thiel et al., 2003). Compound microsatellites were defined as those with the interval between two repeats motifs <=100 nucleotides as previous reports (Zhao et al., 2017).

For primer designing, SSRs mined from MISA along with their locations were processed using p3_in.pl perl script of MISA. Backward and forward primers from flanking sequences of each identified SSR were also extracted using Primer3 executable (Untergasser et al., 2012). The parameters used for primer designing are as follows: 18–27 bp primer length, 57–63°C melting temperature, 30–70% GC content, and 100–300 bp product size. The output obtained from Primer3 was processed using p3_out.pl perl script to extract the final set of primer sequences.

Extraction of Genomic Frequency and Distribution of SSRs

Genomic frequencies per Mb and Kb of SSRs throughout genome and on all individual chromosomes were calculated considering size of the genome and sizes of all individual chromosomes, respectively. Genomic distribution into various genomic regions as exon, intron + intergenic regions, promoter (−1 Kb−100 bp), and TTS—transcription termination sites (−100 bp−1 Kb), of all black pepper genomic SSRs was extracted using HOMER (Heinz et al., 2010) utilizing black pepper reference genome annotation file.

Polymorphic SSRs and Their Genomic Distribution

SSR mining was performed for all the 29 black pepper genotypes' libraries using MISA (Beier et al., 2017) along with primer designing using Primer3. The SSRs extracted from each genotype were compared to see the polymorphism with respect to the genome-wide SSRs extracted from the reference black pepper whole genome. The monomorphic SSRs along with polymorphic SSRs were extracted considering similar 5′ and 3′ 20 nucleotide flanking sequences along with variable number of SSR motifs using Perl scripts. Genomic distribution of polymorphic SSRs was extracted using HOMER (Heinz et al., 2010) utilizing black pepper reference genome annotation file. Hypervariable polymorphic SSR markers (those having ≥20 nucleotides) (Temnykh et al., 2001) were also extracted.

Validation of Extracted SSRs

All the wet-lab validated anonymous SSRs from published literature were collected to be validated in our identified list of SSRs using e-PCR technique. For this validation of identified genomic SSRs, wet-lab validated 116 black pepper SSR primer pairs were collected (Menezes et al., 2009; Raghavan et al., 2010; Joy et al., 2011; Jagtap et al., 2016; Wu et al., 2016; Kumari et al., 2019; Meilawati et al., 2020). A total of 3682 primers of 736 SSRs extracted from black pepper transcriptome were used (Hu et al., 2015). These primers were mapped with 200 bp up- and down-stream of each of genomic SSRs of black pepper to find exact matches, taking cue from the concept of PrimerBLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi). In addition, this validation will add the location to the anonymous SSRs in black pepper genome.

Functional Annotation of Genes Associated With Polymorphic SSRs

First, genes associated with polymorphic SSRs were extracted by Perl script using annotation file of black pepper reference genome. Then, functional annotation of extracted genes was performed using BLAST2GO and KEGG pathway analyses (Kanehisa, 2000).

Development of BlackP2MSATdb Database

The Black Pepper Polymorphic Microsatellite database (BlackP2MSATdb) is a three-tier architecture-based relational database having client, middle, and database tier (Figure 1). The analysis results obtained in this study are cataloged in BlackP2MSATdb. A database was prepared in MySQL database, and its web interface was prepared in PHP and HTML, which includes the following steps of data retrieval: (a) A request generates from user's system to web server, (b) a query sends to MySQL database, (c) a database response generates and sends to web interface, and (d) finally, a web server response sends to user's system. Web hosting of this database was done by Apache2 server. BlackP2MSATdb includes details of SSR markers and polymorphic SSRs of black pepper obtained from GBS data analysis of 29 black pepper genotypes. The output of SSR markers gives the information of each marker in terms of its chromosome number, marker type, marker size, start position, end position, forward primer, reverse primer, and the link to genome.

FIGURE 1

Figure 1. Schematic diagram of database development and data retrieval.

Results

Quality Check and de novo Assembly of GBS Data Libraries

The pair-end reads of all the 29 GBS libraries were subjected to the quality check. The number of paired-end reads, GC %, assembled scaffolds, and the number of extracted SSR markers for each of the GBS libraries of 29 genotypes are shown in Table 1. The average number of reads was 1,842,097, while average GC %, scaffolds, and SSR markers were 39.74, 452,843, and 3,907, respectively. Genotype 4226 was found to have the largest number of reads (3,560,292) with 6,153 SSRs, while genotype KS had the largest number of mined SSRs (9,309). Genotype P2 was found to have least number of reads (213,631) and SSRs (1,007).

TABLE 1

Table 1. Reads, quality, scaffolds, and extracted SSRs in each of GBS libraries of 29 black pepper genotypes.

Genome-Wide Mining of SSR Markers and Their Primer Designing

The SSRs (mono- to hexa-nucleotide and compound SSRs) were mined from the black pepper reference genome (genome size ~761 Mb having 45 scaffolds) and GBS libraries of 29 genotypes using MISA tool. A total of 276,230 SSRs were identified, out of which 275,966 SSRs were found in 26 pseudo-chromosomes (termed as Pn1–26) and only 264 SSRs were found in the remaining 19 scaffolds. The genomic distribution of SSRs showed the abundance of SSRs in intergenic + intronic region (175,641), followed by promoter (70,396), TTS (28,861), and exonic regions (1,331). The average distance between two SSRs in black pepper genome was 2.76 Kb, and the frequency of SSRs per Mb was 362.88. A total of 457 types of mono- to hexa-nucleotide SSR motifs were among all detected SSRs. Out of these SSR motifs, 4, 12, 60, 136, 104, and 141 types belong to mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide repeats, respectively (Table 2). Among 276,230 SSRs, the mononucleotide repeats were the most abundant in proportion and frequency per Mb (137,675, 49.84%) followed by the proportion and frequency of dinucleotide (62,083, 22.48%), compound (38,132, 13.80%), tri-nucleotide (29,560, 10.7%), tetra-nucleotide (7,297, 02.64%), penta-nucleotide (562, 0.20%), and hexa-nucleotide (920, 0.33%). Frequencies of a particular motif type (di-, tri-, tetra-, and penta-nucleotide) are represented in Supplementary Table S1. This table shows the frequency and pattern of di- (including all motifs), tri- (including motifs with ≥141 frequency), tetra- (including motifs with ≥10 frequency), penta- (including motifs with ≥3 frequency), and hexa-nucleotide (including motifs with ≥2 frequency) SSR motifs. A total of 12 and 60 SSR motifs extracted were di- and tri-types, respectively (Figures 2A,B), while 36 out of 136 tetra-nucleotide motifs were having ≥3 frequency (Figure 2C). A total of 52 out of 104 penta-nucleotide motifs had the frequency ≥ 5 (Figure 2D). Abundance of TA motif type (27,311) was observed, followed by AT (24,802), AAT (4,219), and TTA (4,132).

TABLE 2

Table 2. Summary of SSRs in black pepper genome.

FIGURE 2

Figure 2. Graphical representation of (A) di- (included all motifs), (B) tri- (included all motifs), (C) tetra- (included motifs with ≥3 frequency), and (D) penta-nucleotide (included motifs with ≥5 frequency) SSR motifs with their frequencies and patterns.

Extraction of Genomic Frequency and Distribution of SSRs

A comprehensive overview of the mined SSRs was studied. Based on the overall size of black pepper genome, the overall frequency of occurrence of SSRs per Mb within the genome was found to be 362.88 SSRs. The average distance between two SSRs was 2.76 Kb.

The stacked bar graph shows the proportions of SSRs for each chromosome with SSR motifs stacked as mono-, di-, tri-, tetra-, penta-, hexa-, and compound types (Figure 3), while Figure 4 delineates the chromosome-wise numbers of mono- to hexa-nucleotides SSRs in three dimensions. For each of the pseudo-chromosomes of black pepper, the frequency distribution of SSRs per Mb and the average distance between two SSRs in Kb were found (Table 3). It was seen that the pseudo-chromosome Pn1 had the largest number of SSRs (17501), followed by Pn2, Pn5, Pn3, Pn4, Pn7, Pn8, Pn9, Pn13, Pn10, Pn15, Pn12, Pn11, Pn14, and Pn6. The minimum number of SSRs was detected in PnU (280). The highest frequency of SSR markers was observed on Pn7 (398.4), while the least was on Pn6. The average distance between two SSRs is the maximum (3.73) in Pn6 and the least in Pn7.

FIGURE 3

Figure 3. Chromosome-wise distribution of mono- to hexa-nucleotide SSRs and compound SSRs.

FIGURE 4

Figure 4. Distribution of SSR motif repeat number from mono- to hexa-nucleotide along with compound SSRs.

TABLE 3

Table 3. Distribution of SSRs along with frequency per Mb and distance between two SSRs in each of pseudo-chromosomes of black pepper.

Polymorphic Markers, Their Genomic Distribution, and Hypervariable Polymorphic Markers

For the extraction of polymorphic SSRs, genomic SSRs extracted from black pepper reference genome were compared with the SSRs extracted from GBS libraries of 29 genotypes. A total of 5,949 genomic SSRs were found common in both black pepper reference genome and 29 genotypes collectively, out of which 2,773 genomic SSRs were monomorphic and 3,176 genomic SSRs were polymorphic. The number of monomorphic and polymorphic SSRs in each of the GBS libraries of 29 black pepper genotypes along with their genomic distributions is shown in Table 4. Figure 5 shows chromosome-wise distribution of polymorphic SSRs which obviates the abundance of polymorphic SSRs on chromosome Pn1, followed by Pn2, Pn4, and Pn3. Figure 6 shows the chromosome-wise distribution of mono-, di-, tri-, tetra-, penta-, hexa-nucleotide, and compound motifs of polymorphic SSRs found in each of the 29 genotypes. Out of all polymorphic SSRs, 2015 were hypervariable in reference genome and/or in the particular genotype where it is showing polymorphism. Table 5 shows mono-, di-, tri-, tetra-, penta-, hexa-nucleotide, and compound motif-wise distribution of hypervariable polymorphic SSRs. Chromosome-wise distribution of hypervariable polymorphic SSRs along with frequency per Mb and distance in Kb between two hypervariable SSRs is represented in Table 6.

TABLE 4

Table 4. Monomorphic and polymorphic SSRs in each of the GBS libraries of 29 black pepper genotypes and their genomic distributions.

FIGURE 5

Figure 5. Distribution of polymorphic SSRs on different chromosomes of black pepper genome.

FIGURE 6

Figure 6. Distribution of polymorphic SSR motif repeat number from mono- to hexa-nucleotide along with compound SSRs (termed as Mix in figure). The y-axis shows the abundance of polymorphic SSRs that have different motif repeat numbers on x-axis in different genotypes (29 genotypes of black pepper) on z-axis which are discriminated by legends of different colors.

TABLE 5

Table 5. Frequency of hypervariable SSRs with di-, tri-, tetra-, penta-, and hexa-nucleotide motifs along with frequency per Mb and distance in Kb between two hypervariable SSRs.

TABLE 6

Table 6. Chromosome-wise frequency of hypervariable SSRs along with frequency per Mb and distance in Kb between two hypervariable SSRs.

Validation of Extracted SSRs

For the validation of genomic-wide SSRs mined from black pepper reference genome, wet-lab validated and previously reported primers from eight studies were utilized as mentioned under methodology section. A total of 3,798 primer sequences (116 wet-lab validated SSRs from literature with one primer pair plus 736 SSRs with 5 primer pairs (forward and reverse) for most of SSRs from Hu et al., 2015) of 852 validated SSR repeats from various published literature were used for SSR validation. The validated microsatellites (both in vivo and in silico) used for in silico re-validation were taken from eight different studies represented by countries like China, Brazil, and India. Results obtained are shown in Table 7 along with SSR motifs and references of re-validated primers of 360 SSR repeats (primer IDs provided in Supplementary Table S2), out of which 330 SSRs were from Hu et al. (2015).

TABLE 7

Table 7. Validated black pepper genome-wide SSR markers and their method of SSR discovery (The bold ones are the polymorphic SSRs, which could not be determined in case of SSRs listed in Kumari et al., 2019).

All the anonymous in vitro /wet-lab validated SSRs from published literature were validated in silico in our identified list of SSRs using e-PCR technique which is a kind of re-validation of wet-lab validated markers. In addition, this validation will add the location to those anonymous SSRs in black pepper genome.

Functional Annotation of Genes Associated With Polymorphic SSRs

A total of 835 genes were found associated with polymorphic SSR markers extracted from annotation file of black pepper reference genome. The functional annotation of these gene was performed using Blast2GO and KEGG pathway analyses. A total of 30, 36, and 55 GO terms were found in cellular component, molecular function, and biological process classes, respectively. Figure 7 shows GO terms which were associated with ≥10 genes in all three classes. A total of 67 KEGG pathways were found associated with polymorphic SSRs, out of which pathways associated with ≥3 genes are shown in Figure 8.

FIGURE 7

Figure 7. GO terms of genes associated with polymorphic SSRs in black pepper (each term with ≥10 such genes).

FIGURE 8

Figure 8. KEGG pathways associated with polymorphic SSRs in black pepper (each pathway with ≥3 such genes).

Development of BlackP2MSATdb Database

The Black Pepper Polymorphic Microsatellite Database, BlackP2MSATdb genomic resource, is successfully developed using three-tier architecture containing 276,320 SSR loci which are mined from whole-genome sequence data black pepper (Hu et al., 2015) and the GBS data of the 29 black pepper genotypes. It also catalogs 3,176 polymorphic markers which were found polymorphic in 29 genotypes. It has four tabs, viz., Home, Statistics, Data, and Team. The “Home” page has the brief introduction about database. The “Statistics” page included two graphical representation of the proportion of all included data as a sum of all SSRs along with all 3,176 polymorphic SSRs in 29 different black pepper genotypes. The “Data” page navigates the users to get the chromosome-wise SSRs and the motif types. It also provides the polymorphic markers of the 29 black pepper genotypes with respect to the black pepper genome. The “Team” page included the team member name and link to the profile page of each member. The various interfaces of BlackP2MSATdb genomic resource are shown in Figure 9. BlackP2MSATdb is freely available at http://webtom.cabgrid.res.in/blackp2msatdb/.

FIGURE 9

Figure 9. Screenshot of interfaces of BlackP2MSATdb genomic resource.

Discussion

This is the first comprehensive study for genome-wide microsatellite mining from black pepper latest reference genome recently sequenced (Hu et al., 2019). ddRAD data of 29 black pepper genotypes including released varieties, important germplasm accessions, and wild relatives have been successfully used for discovery of polymorphic markers by computational approach. Similar approach of efficient SSR mining has been reported in other plants also (Xiao et al., 2016). This study is the first in this species, utilizing genomic SSRs for rapid extraction of polymorphic SSRs in other 29 black pepper genotypes to study genetic diversity. Black pepper has immense economic importance world-wide, yet no study till date provides physical localization of black pepper SSRs on 26 chromosomes. Present finding in the form of genomic resources can be used in rapid genetic dissection of complex traits including identification of QTLs for various desired traits, resistance to biotic and abiotic stress, genetic enhancement, and varietal development (Guzmán et al., 2020; Nalbandyan et al., 2020; Younis et al., 2020).

The previous study involving whole-genome SSR mining gave 69126 putative SSR markers utilizing the draft genome comprising of 916 scaffolds and the genome coverage of 80X (Kumari et al., 2019). The total mined SSR markers in the present study are much larger (~4-fold) in number as compared to the previously reported black pepper SSR markers by Kumari et al. (2019), which is due to more than 2-fold depth of the genome assembly being utilized. Around 99.87% of identified SSR markers were localized on 26 chromosomes. Such chromosome and location-specific SSR loci can be further used to trace the inheritance of particular chromosomal regions in molecular breeding programs from the foundation genotypes (Saghai-Maroof et al., 1984; Ganie et al., 2016). SSR markers selected uniformly by taking the advantage of their known physical location over chromosome can be a better diversity calculus for genomic and population variability (You et al., 2005) with respect to random anonymous SSR loci as provided in previous studies (Kumari et al., 2019).

Genomic frequency of identified SSR per Mb and distance between two SSRs suggests that the SSRs are distributed ubiquitously over the entire chromosomes (~362 SSRs/Mb), which can prove to be a better representative for the analysis of genome variability as suggested by Guo et al. (2009). It has the advantage being specific to the location over the targeted chromosome and hence can be used as the marker of any specific genes. Similar density of 341–671 SSRs/Mb has been reported in case of sugar beet (Iquebal et al., 2015).

Identified genomic SSRs with mononucleotide repeats were the most abundant in proportion, followed by dinucleotide and compound repeats. The abundance of mono-repeat was observed in all the genotypes which might be due to the inherent limitation of the next-generation sequencing (NGS) chemistry used for the data generation as reported by Haseneyer et al. (2011). The next abundant dinucleotide motif type has also been widely reported in other crops like wheat, watermelon, P. dactylifera, E. oleifera, and E. guineensis (Kariin and Burge, 1995; Shioiri and Takahata, 2001; Kumpatla and Mukhopadhyay, 2005; Deng et al., 2016; Zhu et al., 2016).

In the next part of this study, genetic diversity assessment based on black pepper genomic SSRs was performed in 29 black pepper genotypes by rapid extraction of polymorphic SSRs. A total of 3160 genomic SSRs were found to be polymorphic in these 29 genotypes, most of which were intergenic along with promoter and TTS. It is in concordance with the fact that intergenic SSRs have higher utility in DNA fingerprinting and varietal identification in comparison with genic SSRs (Kalia et al., 2011). Similarly, EST-SSRs extracted from expression data were useful for genetic analysis and were found to be relatively low polymorphic and concentrated in genic regions reported in other studies of black pepper (Gordo et al., 2012; Hu et al., 2015). Thus, the present study further strengthens the fact of abundant occurrence of polymorphic SSRs in intergenic region. Functional annotation of genes associated with polymorphic SSRs suggests that the most of these genes were involved in molecular function of various metabolic pathways, as reported in Curcuma alismatifolia (Taheri et al., 2019).

The provided method of rapid extraction of polymorphic SSRs could be helpful to extract putative polymorphic SSRs in large group of genotypes. In order to overcome the challenge of traditional way of differentiating varieties with identical morphological characteristics based on phenotypic observation, the revolutionary tool of SSRs has proved to be very effective. They play an important role in identification and development of variety, improvement of trait of interest, QTL mapping, and other diversity analysis. They also aid in the product traceability, the evolutionary, and genetic studies including taxonomic comparison (Zietkiewicz et al., 1994; Iquebal et al., 2013; Yu et al., 2017). The phylogeny of nine black peppers varieties is reported to be determined using the morphological and molecular markers which can be important for breeding scheme (Meilawati et al., 2020).

The catalog of putative SSR markers in BlackP2MSATdb can be a very useful web resource for the varietal differentiation since it also houses the ready-to-use primers required for the rapid genotyping. Literature mentions such approach for varietal differentiation as in case of variety of crops like sugarcane, sesame, capsicum, tea, barley, and eggplant (Karakousis et al., 2003; Manigbas and Villegas, 2004; Stàgel et al., 2008; Sahu et al., 2012; Shirasawa et al., 2013; Bhattacharyya et al., 2014). The web-genomic resource, BlackP2MSATdb, of these extracted genomic and polymorphic SSR markers was prepared to provide access to scientific community to utilize them for further studies involving molecular breeding along with construction of physical map, map-based gene cloning, marker-assisted selection, and evolutionary studies in black pepper and its related species. Also, the developed BlackP2MSATdb resource can be useful regarding the phenomena like “microsatellite mediated enhancement” (MME) of gene expression (Kumar and Bhatia, 2016).

Such chromosome and location-specific SSR loci can be further used to trace the inheritance of particular chromosomal regions in molecular breeding programs from the foundation genotypes (Ganie et al., 2016). The ddRAD method is advantageous in the sense that it is very low-cost and most efficient technique to identify polymorphic SSR markers (Razak et al., 2021). The major limitation of ddRAD, that is, non-uniform coverage of the entire genome, has been successfully overcome in the present work by mapping the SSRs on latest available reference sequence of black pepper (Aballay et al., 2021). The present genomic resource, BlackP2MSATdb having chromosome-wise, location-specific SSR markers, can be of immense use in making of coreset and mapping of genes/QTLs for germplasm management and improvement. The obtained >3 thousand polymorphic SSR markers in this study are very valuable as they are distributed over all the chromosomes. Such markers are imperatively needed for efficient and accurate diversity study and development of varietal signature (Korir et al., 2013). From this work, user can now select equal number of loci from each chromosome for diversity analysis.

We found high degree of polymorphism in all the 29 genotypes with respect to the number of SSRs discovered in each genotype. The obtained >3K polymorphic SSR markers in this study are very valuable as they are distributed over all the chromosomes. Such markers are imperatively needed for efficient and accurate diversity study and development of varietal signature (Bohra et al., 2017). This work clearly demonstrates that ddRAD data are still advantageous in terms of not only rapidity and computational ease but also per unit cost of polymorphic SSR marker discovery which is economical than the whole-genome sequencing. Present finding in the form of genomic resources can be used in rapid genetic dissection of complex traits including identification of QTLs for various desired traits, resistance to biotic and abiotic stress, genetic enhancement, and varietal development (Guzmán et al., 2020; Nalbandyan et al., 2020; Younis et al., 2020).

It was also found in this study that a good number of polymorphic SSRs were hypervariable in reference genome and/or in black pepper genotypes, which is in corroboration of the fact that the hypervariable microsatellites provide a general source of polymorphic DNA markers (Powell et al., 1995). These hypervariable polymorphic SSRs could be utilized in various studies involving diversity analysis, hybrid purity testing, and trait mapping (Marker trait association-MTA) and have been utilized in various previous studies (Singh et al., 2012; Varshney et al., 2012; Dutta et al., 2013; Bohra et al., 2017) in case of pigeon pea.

Conclusion

The ddRAD-based GBS approach was used for the rapid discovery of genome-wide location-specific polymorphic SSR markers in black pepper. These polymorphic markers were obtained from 29 diverse black pepper genotypes covering the released varieties, important germplasm accessions, and wild relatives. The first comprehensive web-genomic resource of black pepper, BlackP2MSATdb, was developed that catalogs genome-wide deep mining of SSR markers from the genome assembly of black pepper with 26 chromosomes at a single platform. These markers can be searched chromosome-wise along with their physical location. The validation of these markers was done by e-PCR technique by considering all the anonymous wet-lab validated SSRs from published literature. BlackP2MSATdb catalogs a total of 276,230 putative SSR markers identified from the whole genome of black pepper, with relative density 362.88 SSRs per Mb and average distance between two SSRs being 2.76 Kb. A total of 3,176 polymorphic markers were common to black pepper reference genome and 29 genotypes collectively, out of which 2015 were hypervariable. This web-genomic resource can be of immense use in molecular breeding, construction of physical map, identification of QTL, map-based gene cloning, marker-assisted selection, and evolutionary studies in black pepper and its related species.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found at: National Center for Biotechnology Information (NCBI) BioProject database under accession number PRJNA489154.

Author Contributions

SJ, MI, JK, and DK conceived theme of the study and designed experiment. JK and PU collected the samples. KS, AN, MI, SJ, and UA were involved in computational analysis and development of web resources. AR collected resources. SJ, MI, DK, and AR supervised the study. AN, KS, SJ, and MI were involved in writing the original draft. DK and AR reviewed and edited the manuscript. All authors read and approved the final manuscript and contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors are thankful to CABin grant, Indian Council of Agricultural Research, Ministry of Agriculture and Farmers' Welfare, Govt. of India (F. no. Agril. Edn. 4–1/2013-A&P) for providing financial and infrastructural support to carry out this research and for creation of Advanced Super Computing Hub for Omics Knowledge in Agriculture (ASHOKA) facility where the work was carried out. The grant of IARI Merit scholarship to AN is duly acknowledged.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.846937/full#supplementary-material

Supplementary File 1. Black pepper variety details.

Supplementary Table S1. SSR motifs and their frequency for di- (including all motifs), tri- (including motifs with frequency ≥141), tetra- (including motifs with frequency ≥10), penta- (including motifs with frequency ≥3), and hexa-nucleotide (including motifs with frequency ≥2) type.

Supplementary Table S2. Set of validated microsatellites used for in silico re-validation.

Abbreviations

SSR, simple sequence repeats; QTL, quantitative trait loci; DNA, deoxyribonucleic acid; GBS, genotyping-by-sequencing; CTAB, cetyltrimethylammonium bromide; HS, high sensitivity; MISA, MIcroSAtellite identification; GC, guanine–cytosine; TTS, transcription termination sites; HOMER, Hypergeometric Optimization of Motif EnRichment; BLAST, Basic Local Alignment Search Tool; KEGG, Kyoto Encyclopedia of Genes and Genomes; NGS, next-generation sequencing.

References

Aballay, M. M., Aguirre, N. C., Filippi, C. V., Valentini, G. H., and Sánchez, G. (2021). Fine-tuning the performance of ddRAD-seq in the peach genome. Sci. Rep. 11, 1–13. doi: 10.1038/s41598-021-85815-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Beier, S., Thiel, T., Münch, T., Scholz, U., and Mascher, M. (2017). MISA-web: a web server for microsatellite prediction. Bioinformatics 33, 2583–2585. doi: 10.1093/bioinformatics/btx198

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhattacharyya, U., Pandey, S. K., and Dasgupta, T. (2014). Identification of EST-SSRs and FDM in sesame (Sesamum indicum L.) through data mining. Sch. J. Agric. Sci. 4, 60–69. Available online at: http://www.scholarly-journals.com/sjas/archive/2014/February/pdf/Bhattacharyya%20et%20al.pdf

Google Scholar

Bohra, A., Jha, R., Pandey, G., et al. (2017). New hypervariable SSR markers for diversity analysis, hybrid purity testing and trait mapping in pigeonpea [Cajanus cajan (L.) Millspaugh]. Front. Plant Sci. 8:377. doi: 10.3389/fpls.2017.00377

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, P., Wang, M., Feng, K., Cui, L., Tong, W., Song, W., et al. (2016). Genome-wide characterization of microsatellites in Triticeae species: abundance, distribution and evolution. Sci. Rep. 6, 1–13. doi: 10.1038/srep32224

PubMed Abstract | CrossRef Full Text | Google Scholar

Dutta, S., Mahato, A. K., Sharma, P., Ranjeet, S. R., Sharma, T. R., and Singh, N. K. (2013). Highly variable Arhar simple sequence repeat markers for molecular diversity and phylogenetic studies in pigeonpea (Cajanus cajan (L.) Millisp.). Plant Breed. 132, 1439–0523. doi: 10.1111/pbr.12014

CrossRef Full Text | Google Scholar

FAO (2019). The State of Food and Agriculture 2019. Moving Forward on Food Loss and Waste Reduction. Rome.

Ganie, S. A., Borgohain, M. J., Kritika, K., Talukdar, A., Pani, D. R., and Mondal, T. K. (2016). Assessment of genetic diversity of Saltol QTL among the rice (Oryza sativa L.) genotypes. Physiol. Mol. Biol. Plants 22, 107–114. doi: 10.1007/s12298-016-0342-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Gordo, S., Pinheiro, D. G., Moreira, E. C., Rodrigues, S. M., Poltronieri, M. C., de Lemos, O. F., et al. (2012). High-throughput sequencing of black pepper root transcriptome. BMC Plant Biol. 12, 1–9. doi: 10.1186/1471-2229-12-168

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, W. J., Ling, J., and Li, P. (2009). Consensus features of microsatellite distribution: microsatellite contents are universally correlated with recombination rates and are preferentially depressed by centromeres in multicellular eukaryotic genomes. Genomics 93, 323–331. doi: 10.1016/j.ygeno.2008.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Guzmán, F. A., Moore, S., de Vicente, M. C., and Jahn, M. M. (2020). Microsatellites to enhance characterization, conservation and breeding value of Capsicum germplasm. Genet. Resour. Crop Evol. 67, 569–585. doi: 10.1007/s10722-019-00801-w

CrossRef Full Text | Google Scholar

Haseneyer, G., Schmutzer, T., Seidel, M., Zhou, R., Mascher, M., Schön, C. C., et al. (2011). From RNA-seq to large-scale genotyping-genomics resources for rye (Secale cereale L.). BMC Plant Biol. 11, 1–13. doi: 10.1186/1471-2229-11-131

PubMed Abstract | CrossRef Full Text | Google Scholar

Heinz, S., Benner, C., Spann, N., Bertolino, E., Lin, Y. C., Laslo, P., et al. (2010). Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576–589. doi: 10.1016/j.molcel.2010.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, L., Hao, C., Fan, R., Wu, B., Tan, L., and Wu, H. (2015). De novo assembly and characterization of fruit transcriptome in black pepper (Piper nigrum). PLoS ONE. 10, e0129822. doi: 10.1371/journal.pone.0129822

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, L., Xu, Z., Wang, M., Fan, R., Yuan, D., Wu, B., et al. (2019). The chromosome-scale reference genome of black pepper provides insight into piperine biosynthesis. Nat. Commun. 10, 1–11. doi: 10.1038/s41467-019-12607-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Iquebal, M. A., Arora, V., Verma, N., Rai, A., and Kumar, D. (2013). First whole genome based microsatellite DNA marker database of tomato for mapping and variety identification. BMC Plant Biol. 13, 1–7. doi: 10.1186/1471-2229-13-197

CrossRef Full Text | Google Scholar

Iquebal, M. A., Jaiswal, S., Angadi, U. B., Sablok, G., Arora, V., Kumar, S., et al. (2015). SBMDb: first whole genome putative microsatellite DNA marker database of sugarbeet for bioenergy and industrial applications. Database 15, 1–10. doi: 10.1093/database/bav111

PubMed Abstract | CrossRef Full Text | Google Scholar

Jagtap, A. B., Sujatha, R., Nazeem, P. A., Meena, O. P., and Pathania, S. (2016). Morpho-molecular characterization of putative interspecific crosses in black pepper ('Piper nigrum' L. and 'Piper colubrinum'). Plant Omics. 9, 73–80. doi: 10.3316/informit.887902476345816

PubMed Abstract | CrossRef Full Text | Google Scholar

Jose, S., Sujatha, R., and Deeshma, K. P. (2017). Novel EST-SSR marker development and validation in black pepper cultivars and varieties. J. Trop. Agric. 55, 175–179.

Google Scholar

Joy, N., Prasanth, V. P., and Soniya, E. V. (2011). Microsatellite based analysis of genetic diversity of popular black pepper genotypes in South India. Genetica 139, 1033–1043. doi: 10.1007/s10709-011-9605-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalia, R. K., Rai, M. K., Kalia, S., Singh, R., and Dhawan, A. K. (2011). Microsatellite markers: an overview of the recent progress in plants. Euphytica 177, 309–334. doi: 10.1007/s10681-010-0286-9

CrossRef Full Text | Google Scholar

Kanehisa, M. (2000). Post-genome Informatics. Oxford: Oxford University Press.

Google Scholar

Karakousis, A., Barr, A. R., Chalmers, K. J., Ablett, G. A., Holton, T. A., Henry, R. J., et al. (2003). Potential of SSR markers for plant breeding and variety identification in Australian barley germplasm. Aust. J. Agric. Res. 54, 1197–1210. doi: 10.1071/AR02178

CrossRef Full Text | Google Scholar

Kariin, S., and Burge, C. (1995). Dinucleotide relative abundance extremes: a genomic signature. Trends Genet. 11, 283–290. doi: 10.1016/S0168-9525(00)89076-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Korir, N. K., Han, J., Shangguan, L., Wang, C., Kayesh, E., Zhang, Y., et al. (2013). Plant variety and cultivar identification: advances and prospects. Crit. Rev. Biotechnol. 33, 111–125. doi: 10.3109/07388551.2012.675314

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., and Bhatia, S. (2016). A polymorphic (GA/CT) n-SSR influences promoter activity of Tryptophan decarboxylase gene in Catharanthus roseus L. Don. Sci. Rep. 6, 1–12. doi: 10.1038/srep33280

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, R., Wankhede, D. P., Bajpai, A., Maurya, A., Prasad, K., Gautam, D., et al. (2019). Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): a valuable resource for boosting genomics applications. PLoS ONE 14, e0226002. doi: 10.1371/journal.pone.0226002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumpatla, S. P., and Mukhopadhyay, S. (2005). Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species. Genome 48, 985–998. doi: 10.1139/g05-060

PubMed Abstract | CrossRef Full Text | Google Scholar

Manigbas, N. L., and Villegas, L. C. (2004). Microsatellite markers in hybridity tests to identify true hybrids of sugarcane. Philipp J. Crop Sci. 29, 23–32. Available online at: https://www.cabdirect.org/cabdirect/abstract/20093019321

Google Scholar

Meilawati, N. L. W., Susilowati, M., and Bermawie, N. (2020). “Phyllogenetic of nine superior black pepper (Piper nigrum L.) varieties based on morphological and molecular markers,” in IOP Conference Series: Earth and Environmental Science, Vol. 418 (IOP Publishing), 012056. Available online at: https://iopscience.iop.org/article/10.1088/1755-1315/418/1/012056/pdf

Google Scholar

Menezes, I. C., Cidade, F. W., Souza, A. P., and Sampaio, I. C. (2009). Isolation and characterization of microsatellite loci in the black pepper, Piper nigrum L. (piperaceae). Conserv. Genet. Resourc. 1, 209–212. doi: 10.1007/s12686-009-9051-6

CrossRef Full Text | Google Scholar

Nalbandyan, A. A., Hussein, A. S., Fedulova, T. P., Cherepukhina, I. V., Kryukova, T. I., Rudenko, T. S., et al. (2020). Differentiation of sugar beet varieties using SSR markers: a tool to create promising hybrids. Russ. Agric. Sci. 46, 442–446. doi: 10.3103/S1068367420050146

CrossRef Full Text | Google Scholar

Pestsova, E., and Röder, M. (2002). Microsatellite analysis of wheat chromosome 2D allows the reconstruction of chromosomal inheritance in pedigrees of breeding programmes. Theor. Appl. Genet. 106, 84–91. doi: 10.1007/s00122-002-0998-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Powell, W., Morgante, M., Andre, C., McNicol, J. W., Machray, G. C., Doyle, J. J., et al. (1995). Hypervariable microsatellites provide a general source of polymorphic DNA markers for the chloroplast genome. Curr. Biol. 5, 1023–1029. doi: 10.1016/S0960-9822(95)00206-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Prjibelski, A., Antipov, D., Meleshko, D., Lapidus, A., and Korobeynikov, A. (2020). Using SPAdes de novo assembler. Curr. Protoc. Bioinformatics 70, e102. doi: 10.1002/cpbi.102

PubMed Abstract | CrossRef Full Text | Google Scholar

Raghavan, R., Elumalai, S., Babu, K. N., and Hittalmani, S. (2010). Molecular characterization of black pepper (Piper nigrum) using RAPD and SSR markers. Biosci. Biotechnol. Res. Asia. 7, 1011–1015. Available online at: https://pdfs.semanticscholar.org/1559/7311fd0492d75e54ff0b5fced542c3b27dbd.pdf

Google Scholar

Razak, S. A., Ghazalli, M. N., Azman, N. H. E. N., Majid, A. M. A., and Ismail, S.N. (2021). RAD sequencing for the development of microsatellite markers for identification of Malaysian taro cultivars. Biotechnol. Biotechnol. Equip. 35, 1284–1290. doi: 10.1080/13102818.2021.1969278

CrossRef Full Text | Google Scholar

Saghai-Maroof, M. A., Soliman, K. M., Jorgensen, R. A., and Allard, R. (1984). Ribosomal DNA spacer-length polymorphisms in barley: Mendelian inheritance, chromosomal location, and population dynamics. Proc. Nat. Acad. Sci. U.S.A. 81, 8014–8018. doi: 10.1073/pnas.81.24.8014

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahu, J., Sarmah, R., Dehury, B., Sarma, K., Sahoo, S., Sahu, M., et al. (2012). Mining for SSRs and FDMs from expressed sequence tags of Camellia sinensis. Bioinformation 8, 260. doi: 10.6026/97320630008260

PubMed Abstract | CrossRef Full Text | Google Scholar

Shioiri, C., and Takahata, N. (2001). Skew of mononucleotide frequencies, relative abundance of dinucleotides, and DNA strand asymmetry. J. Mol. Evol. 53, 364–376. doi: 10.1007/s002390010226

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirasawa, K., Ishii, K., Kim, C., Ban, T., Suzuki, M., Ito, T., et al. (2013). Development of Capsicum EST-SSR markers for species identification and in silico mapping onto the tomato genome sequence. Mol. Breed. 31, 101–110. doi: 10.1007/s11032-012-9774-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, N. K., Gupta, D. K., Jayaswal, P. K., Mahato, A. K., Dutta, S., Singh, S., et al. (2012). The first draft of the pigeonpea genome sequence. J. Plant Biochem. Biotechnol. 21, 98–112. doi: 10.1007/s13562-011-0088-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Stàgel, A., Portis, E., Toppino, L., Rotino, G. L., and Lanteri, S. (2008). Gene-based microsatellite development for mapping and phylogeny studies in eggplant. BMC Genom. 9, 1–14. doi: 10.1186/1471-2164-9-357

PubMed Abstract | CrossRef Full Text | Google Scholar

Taheri, S., Abdullah, T. L., Rafii, M. Y., Harikrishna, J. A., Werbrouck, S. P., Teo, C. H., et al. (2019). De novo assembly of transcriptomes, mining, and development of novel EST-SSR markers in Curcuma alismatifolia (Zingiberaceae family) through illumina sequencing. Sci. Rep. 9, 1–14. doi: 10.1038/s41598-019-39944-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Temnykh, S., DeClerck, G., Lukashova, A., Lipovich, L., Cartinhour, S., and McCouch, S. (2001). Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genet. Res. 11, 1441–1452. doi: 10.1101/gr.184001

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiel, T., Michalek, W., Varshney, R., and Graner, A. (2003). Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor. Appl. Genet. 106, 411–422. doi: 10.1007/s00122-002-1031-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M., et al. (2012). Primer3-new capabilities and interfaces. Nucleic Acids Res. 40, e115. doi: 10.1093/nar/gks596

PubMed Abstract | CrossRef Full Text | Google Scholar

Varshney, R. K., Chen, W., Li, Y., Bharti, A. K., Saxena, R. K., Schlueter, J. A., et al. (2012). Draft genome sequence of pigeon pea (Cajanus cajan), an orphan legume crop of resource-poor farmers. Nat. Biotechnol. 30, 83–89. doi: 10.1038/nbt.2022

PubMed Abstract | CrossRef Full Text | Google Scholar

Wingett, S. W., and Andrews, S. (2018). FastQ Screen: a tool for multi-genome mapping and quality control. F1000Research 7:1338. doi: 10.12688/f1000research.15931.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B. D., Fan, R., Hu, L. S., Wu, H. S., and Hao, C. Y. (2016). Genetic diversity in the germplasm of black pepper determined by EST-SSR markers. Genet. Mol. Res. 15, 1–9. doi: 10.4238/gmr.15018099

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, Y., Xia, W., Ma, J., Mason, A. S., Fan, H., Shi, P., et al. (2016). Genome-wide identification and transferability of microsatellite markers between Palmae species. Front. Plant Sci. 7:1578. doi: 10.3389/fpls.2016.01578

PubMed Abstract | CrossRef Full Text | Google Scholar

You, G. X., Zhang, X. Y., and Wang, L. F. (2005). An estimation of the minimum number of SSR loci needed to reveal genetic relationships in wheat varieties: information from 96 random accessions with maximized genetic diversity. Mol. Breed. 14, 397–406. doi: 10.1007/s11032-005-0285-z

CrossRef Full Text | Google Scholar

Younis, A., Ramzan, F., Ramzan, Y., Zulfiqar, F., Ahsan, M., and K. B. (2020). Molecular markers improve abiotic stress tolerance in crops: a review. Plants 9, 1374. doi: 10.3390/plants9101374

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J., Dossa, K., Wang, L., Zhang, Y., Wei, X., Liao, B., et al. (2017). PMDBase: a database for studying microsatellite DNA and marker development in plants. Nucleic Acids Res. 45, D1046–D1053. doi: 10.1093/nar/gkw906

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Qiu, J., Agarwal, G., Wang, J., Ren, X., Xia, H., et al. (2017). Genome-wide discovery of microsatellite markers from diploid progenitor species, Arachis duranensis and A. ipaensis, and their application in cultivated peanut (A. hypogaea). Front. Plant Sci. 8:1209. doi: 10.3389/fpls.2017.01209

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H., Song, P., Koo, D. H., Guo, L., Li, Y., Sun, S., et al. (2016). Genome wide characterization of simple sequence repeats in watermelon genome and their application in comparative mapping and genetic diversity analysis. BMC Genom. 17, 1–17. doi: 10.1186/s12864-016-2870-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zietkiewicz, E., Rafalski, A., and Labuda, D. (1994). Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification. Genomics 20, 176–183. doi: 10.1006/geno.1994.1151

PubMed Abstract | CrossRef Full Text | Google Scholar