Microsatellite based analysis of genetic diversity of popular black pepper genotypes in South India

Genetic diversity and population structure of Piper nigrum (black pepper) accessions based on next-generation SNP markers

Despite the economic importance of Piper nigrum (black pepper), a highly valued crop worldwide, development and utilization of genomic resources have remained limited, with diversity assessments often relying on only a few samples or DNA markers. Here we employed restriction-site associated DNA sequencing to analyze 175 P. nigrum accessions from eight main black pepper growing regions in Sri Lanka. The sequencing effort resulted in 1,976 million raw reads, averaging 11.3 million reads per accession, revealing 150,356 high-quality single nucleotide polymorphisms (SNPs) distributed across 26 chromosomes. Population structure analysis revealed two subpopulations (K = 2): a dominant group consisting of 152 accessions sourced from both home gardens and large-scale cultivations, and a smaller group comprising 23 accessions exclusively from native collections in home gardens. This clustering was further supported by principal component analysis, with the first two principal components explaining 35.2 and 12.1% of the total variation. Genetic diversity analysis indicated substantial gene flow (Nm = 342.21) and a low fixation index (FST = 0.00073) between the two subpopulations, with no clear genetic differentiation among accessions from different agro-climatic regions. These findings demonstrate that most current black pepper genotypes grown in Sri Lanka share a common genetic background, emphasizing the necessity to broaden the genetic base to enhance resilience to biotic and abiotic stresses. This study represents the first attempt at analyzing black pepper genetic diversity using high-resolution SNP markers, laying the foundation for future genome-wide association studies for SNP-based gene discovery and breeding.

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.

Chemical analysis and classification of black pepper (Piper nigrum L.) based on their country of origin using mass spectrometric methods and chemometrics

The current study applied gas chromatography-mass spectrometry (GC-MS), liquid chromatography-mass spectrometry (LC-MS), and thermal desorption direct analysis in real-time mass spectrometry (TD-DART-MS) methods to the analysis of black pepper (Piper nigrum L.) samples from different countries. The black pepper powder samples were analyzed directly by TD-DART-MS without any extraction, but for GC-MS and LC-MS methods, a methanol extraction procedure was employed before the analysis. Various compounds, such as piperamides and terpenes, were detected. Partial least squares-discriminant analysis (PLS-DA) was used to classify black pepper samples based on their origins. Total ion mass spectrum (TMS) data profiles from GC-MS, LC-MS, and TD-DART-MS methods were constructed and evaluated for the performance of classification. A cubic-root data transformation was tested in the data preprocessing and found to be effective for improving the classification rates. The average classification rates of PLS-DA models with GC-MS-cubic-root-TMS, LC-MS-cubic-root-TMS, and DART-MS-cubic-root-TMS data representations were 94.1 ± 0.6%, 87.7 ± 0.6%, and 97.0 ± 0.3% respectively, for 100-time bootstrapped-Latin-partition cross-validation. This study presents for the first time the analysis of plant-based food materials by using TD-DART-MS, and it has been demonstrated as a simple and high-throughput method for classification studies.

Evaluation of DNA markers for molecular identification of three Piper species from Brazilian Atlantic Rainforest

Piper is one of two large genera in the Piperaceae, and with ca. 2600 species, is one of the largest plant genera in the world. Species delimitation and evaluation of genetic diversity among populations are important requisites for conservation and adequate exploitation of economically important species. DNA barcoding has been used as a powerful tool and a practical method for species characterization and delimitation. The present work aims to evaluate molecular markers for barcoding three Piper species native to Brazil: P. gaudichaudianum (“jaborandi” or “pariparoba”), P. malacophyllum (“pariparoba-murta”) and P. regnellii (“caapeba” or “pariparoba”). A reference DNA barcode library was developed using sequences of three candidate regions: ITS2, trnH-psbA and rbcL. Transferability of the microsatellite (SSR) primers Psol 3, Psol 6 and Psol 10, designed originally for Piper solmsianum, to the three Piper species was also evaluated. The discriminatory power of the markers was based on the determination of inter- and intraspecific distances, phylogenetic reconstruction, and clustering analysis, as well as BLASTn comparison. Sequences of ITS2 enabled efficient species identification by means of the BLASTn procedure. Based on these sequences, intraspecific divergence was lower than interspecific variation. Maximum Parsimony analyses based on ITS2 sequences provided three resolved clades, each corresponding to one of the three analysed species. Sequences of trnH-psbA and rbcL had lower discriminatory value. Analyses combining sequences of these regions were less effective toward the attainment of resolved and strongly supported clades of all species. In summary, robustly supported clades of P. regnellii were obtained in most of the analyses, based either on isolated or combined sequences. The SSRs primers Psol 3, Psol 6 and Psol 10 were shown to be transferable to P. gaudichaudianum and P. regnellii, but not to P. malacophyllum. Preliminary cluster analyses based on the polymorphism of the amplified products suggested that Psol 3 has lower potential than Psol 6 and Psol 10 for discrimination of Piper species.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

Sensitive detection of Piper nigrum L. adulterants by a novel screening approach based on qPCR

The spice made from the fruits of Piper nigrum L. (Piperaceae) has high economic value since the beginnings of international trade. Because its price has been increasing, adulterations with papaya seeds, cayenne pepper and maize flour were reported. These have been screened by methodologies dedicated to the detection of single adulterants lacking sensitivity and specificity. Herein we propose a specific, highly-sensitive, high-throughput and affordable qPCR-based methodology for the detection of P. nigrum contaminants (Carica papaya, Zea mays and Capsicum annuum) using plant DNA barcodes trnL and psbA-trnH. The method enables the specific detection of contaminants in a short time with low limits of detection (LOD₆ values of 1, 2 and 10 Haploid Genome Equivalents). A market survey (29 samples) revealed 41% of samples contaminated, though about ¾ at very low levels indicating accidental contamination. The proposed tool will contribute to the improvement of quality of this much traded spice.

Isolation and characterization of genomic microsatellite markers for small cardamom (Elettaria cardamomum Maton) for utility in genetic diversity analysis

Microsatellite markers in small cardamom (Elettaria cardamomum Maton) were developed using the selective hybridization enrichment method. A total of 140 microsatellite repeats were identified from 270 clones. Primers were designed for 58 microsatellites and 44 primer pairs amplified products of expected size in cardamom. These markers were used for studying the diversity of 20 important small cardamom genotypes, and six markers were found to be polymorphic. The number of alleles ranged from 2 to 7 with an average of 3.6 per locus. Polymorphic information content values ranged from 0.14 to 0.38 based on dominant scoring. The two markers ECM 47a and ECMG 28 generated specific banding patterns for the genotypes MCC7 (Pink tiller) and APG434 (MA18) respectively. Dendrogram illustrated the genetic similarity between different genotypes of Kerala and Karnataka regions. It differentiated the closely related genotypes and released varieties into separate groups. Principal coordinate analysis revealed PV1 and ICRI 1 as the most divergent genotypes. The study demonstrated that these markers are informative and can be further utilized for generating reliable molecular data for assisting the crop improvement of small cardamom. Cross generic transferability (71.4 %) of the developed primers proved that they are useful for phylogenetic studies in the family Zingiberaceae. This is the first report of de novo isolation, characterisation and utilization of microsatellite markers for the genetic diversity analysis of small cardamom.

De novo transcriptome sequencing reveals a considerable bias in the incidence of simple sequence repeats towards the downstream of 'Pre-miRNAs' of black pepper

Next generation sequencing has an advantageon transformational development of species with limited available sequence data as it helps to decode the genome and transcriptome. We carried out the de novo sequencing using illuminaHiSeq™ 2000 to generate the first leaf transcriptome of black pepper (Piper nigrum L.), an important spice variety native to South India and also grown in other tropical regions. Despite the economic and biochemical importance of pepper, a scientifically rigorous study at the molecular level is far from complete due to lack of sufficient sequence information and cytological complexity of its genome. The 55 million raw reads obtained, when assembled using Trinity program generated 2,23,386 contigs and 1,28,157 unigenes. Reports suggest that the repeat-rich genomic regions give rise to small non-coding functional RNAs. MicroRNAs (miRNAs) are the most abundant type of non-coding regulatory RNAs. In spite of the widespread research on miRNAs, little is known about the hair-pin precursors of miRNAs bearing Simple Sequence Repeats (SSRs). We used the array of transcripts generated, for the in silico prediction and detection of ‘43 pre-miRNA candidates bearing different types of SSR motifs’. The analysis identified 3913 different types of SSR motifs with an average of one SSR per 3.04 MB of thetranscriptome. About 0.033% of the transcriptome constituted ‘pre-miRNA candidates bearing SSRs’. The abundance, type and distribution of SSR motifs studied across the hair-pin miRNA precursors, showed a significant bias in the position of SSRs towards the downstream of predicted ‘pre-miRNA candidates’. The catalogue of transcripts identified, together with the demonstration of reliable existence of SSRs in the miRNA precursors, permits future opportunities for understanding the genetic mechanism of black pepper and likely functions of ‘tandem repeats’ in miRNAs.

Identification of an miRNA candidate reflects the possible significance of transcribed microsatellites in the hairpin precursors of black pepper

Plant miRNAs (18-24nt) are generated by the RNase III-type Dicer endonuclease from the endogenous hairpin precursors ('pre-miRNAs') with significant regulatory functions. The transcribed regions display a higher frequency of microsatellites, when compared to other regions of the genomic DNA. Simple sequence repeats (SSRs) resulting from replication slippage occurring in transcripts affect the expression of genes. The available experimental evidence for the incidence of SSRs in the miRNA precursors is limited. Considering the potential significance of SSRs in the miRNA genes, we carried out a preliminary analysis to verify the presence of SSRs in the pri-miRNAs of black pepper (Piper nigrum L.). We isolated a (CT) dinucleotide SSR bearing transcript using SMART strategy. The transcript was predicted to be a 'pri-miRNA candidate' with Dicer sites based on miRNA prediction tools and MFOLD structural predictions. The presence of this 'miRNA candidate' was confirmed by real-time TaqMan assays. The upstream sequence of the 'miRNA candidate' by genome walking when subjected to PlantCARE showed the presence of certain promoter elements, and the deduced amino acid showed significant similarity with NAP1 gene, which affects the transcription of many genes. Moreover the hairpin-like precursor overlapped the neighbouring NAP1 gene. In silico analysis revealed distinct putative functions for the 'miRNA candidate', of which majority were related to growth. Hence, we assume that this 'miRNA candidate' may get activated during transcription of NAP gene, thereby regulating the expression of many genes involved in developmental processes.