- Research note
- Open Access
Characterization of 29 polymorphic microsatellite markers developed by genomic screening of Sumatran rhinoceros (Dicerorhinus sumatrensis)
BMC Research Notes volume 14, Article number: 119 (2021)
The Sumatran rhinoceros is critically endangered, with fewer than 100 individuals surviving across its current range. Accurate census estimates of the remaining populations are essential for development and implementation of conservation plans. In order to enable molecular censusing, we here develop microsatellite markers with amplicon sizes of short length, appropriate for non-invasive fecal sampling.
Due to limited sample quantity and potential lack of genome-wide diversity, Illumina sequence reads were generated from two Sumatran rhinoceros samples. Genomic screening identified reads with short tandem repeats and loci that were polymorphic within the dataset. Twenty-nine novel polymorphic microsatellite markers were characterized (A = 2.4; HO = 0.30). These were sufficient to distinguish among individuals (PID < 0.0001), and to distinguish among siblings (PID(sib) < 0.0001). Among rhinos in Indonesia, almost all markers were established as polymorphic and effective for genotyping DNA from fecal samples. Notably, the markers amplified and displayed microsatellite polymorphisms using DNA extracted from 11 fecal samples collected non-invasively from wild Sumatran rhinoceros. These microsatellite markers provide an important resource for a census and genetic studies of wild Sumatran rhinos.
In the past two decades the population of Sumatran rhinoceros (Dicerorhinus sumatrensis) has declined by more than 50% ; with less than one percent of its former geographic range occupied and < 100 individuals of this critically endangered species surviving in isolated populations in Sumatra and Borneo. Despite the need for an accurate census and determination of the remaining genetic diversity for the conservation and management of surviving populations, an accurate census of the remaining populations has been elusive.
Molecular tools to amplify DNA from dung would be of great utility for a critical census of Sumatran rhinoceros. With rigorous fecal collection protocols followed, DNA from dung samples has been used to estimate population sizes in wildlife [2, 3]. Microsatellite markers are ideal for genetic profiling of populations due to high rates of evolution leading to intra-species polymorphisms, well-understood mutation dynamics, and proven utility in wildlife management and conservation of non-model taxa [4, 5]. Dung sample collection is relatively straightforward, may allow for a large proportion of the population to be sampled, and involves no stress from handling or direct observation of individuals . Further, dung sampling allows for real time monitoring of changes in genetic diversity and population dynamics .
The evaluation of microsatellite DNA profiles from dung samples poses two challenges. First, dung samples contain a low quantity and quality of host DNA, as well as inhibitors and contaminants [8,9,10]. Both of these lead to reduced microsatellite amplification success with DNA from dung when compared to DNA from blood or tissue samples. Second, there are a limited number of Sumatran rhinoceros microsatellite markers, and these were not optimized for use with dung DNA. More specifically, the amplicon length had not been minimized .
To generate more precise and accurate census estimates from in situ Sumatran rhinoceros populations, here we report optimized new microsatellite markers that target short regions of Sumatran rhinoceros DNA. We provide evidence of the amplification reliability and the presence of polymorphisms in Sumatran rhinoceros genotypes from dung. Specifically, we report the (a) identification of 29 novel short amplicon polymorphic Sumatran rhinoceros microsatellite markers, (b) characterization of variability of these markers using high quality DNA, (c) optimization of amplification success and verification of polymorphisms of these markers using DNA from dung, and (d) first field results from use of these markers on 11 dung samples from free-ranging Sumatran rhinoceros in Indonesia. Our results suggest that these new tools will be of value in the conservation and management of extant Sumatran rhinoceros.
Methods and materials
Sample collection and DNA extraction
We used several types of Sumatran rhino samples (Table 1): previously extracted high quality DNA (n = 6) from the San Diego Zoo Institute for Conservation Research (ICR) and the Royal Ontario Museum (ROM); blood samples from rhinos kept at the Cincinnati Zoo (n = 2) or at the Sumatran Rhino Sanctuary (SRS) within Way Kambas National Park (WK) in Sumatra (n = 3); fecal samples from captive rhinoceros at the Cincinnati Zoo (n = 2) and SRS (n = 3); and fecal samples (n = 11) from an unknown number of wild rhino individuals at the Burkit Barisan Selatan National Park (BBS) in Sumatra. Across these samples, there were paired blood-fecal samples (each pair from the same individual) for three of the rhinos (two at the Cincinnati Zoo and one at SRS).
At the University of Illinois at Urbana-Champaign (UIUC), high quality DNA was extracted from the Cincinnati Zoo blood samples (n = 2) using the Qiagen DNeasy Blood and Tissue Kit. At the Eijkman Institute for Molecular Biology (EIMB), DNA was extracted from the SRS rhino blood samples (n = 3) using a salting out procedure . At both UIUC and EIMB, fecal DNA was extracted from Cincinnati Zoo (n = 2), SRS (n = 3) and BBS (n = 11) using the QIAmp DNA Stool Kit (Qiagen) with a modified protocol .
Bioinformatic identification of polymorphic microsatellite loci
Whole genome sequences (Illumina MiSeq v3) were generated from two high quality DNA samples from Sumatran rhinoceros Dsu-33 and Dsu-35, both originally from Sumatra (Table 1). Reads from both rhinos were combined, and reads with short tandem repeat (STR) sequences that corresponded to the same locus were identified. Only those combined reads that exhibited polymorphism at an STR locus (across the four chromosomes sequenced) were used to design primers, using MsatCommander v1.0.8 . To increase PCR success with the degraded Sumatran rhino DNA available from dung, primers were designed to amplify a short target product with a maximum length of 150 bp (inclusive of the primer lengths which total 36 to 44 bp). After running in silico PCR with the IPCRESS program , loci showing repetitive elements, monomorphism, a very broad allele size range, or sequences closely matching human were all excluded from consideration.
PCR amplification and locus polymorphism
Candidate primer pairs were evaluated at UIUC (Additional file 1) for amplification reliability and accuracy using available high quality (non-fecal) DNA from six Sumatran rhinos that had not been used in marker development, from the ICR and ROM (Table 1). A standard PCR mix and amplification protocol (Additional file 1) were used to test each of the candidate primer pairs. An agarose gel was used to verify a single amplicon of the correct size.
Estimates of microsatellite variability
Variability at the loci was characterized using six high quality DNA extracts (rhinos Dsu-28, 29, 44, 63, 64 and 66; Table 1). FSTAT, v220.127.116.11  and GenAlEx, v6.1 [17, 18] were used to calculate the number of alleles per locus (A), expected heterozygosity (HE), and observed heterozygosity (HO). CERVUS v3.0.7  was used to calculate the probabilities of individual identity, i.e., PID and PID(sib) values (for unrelated individuals and full siblings, respectively) for each marker, and cumulative PID and PID(sib) across multiple markers, using the equations of Waits et al. 2001 . Cumulative PID and PID(sib) were calculated by multiplying PID and PID(sib) across loci. Using GENEPOP (https://genepop.curtin.edu.au/), significant linkage disequilibrium was not detected among loci. Genome scaffolds were identified for loci (Additional file 1: Table S1).
Amplification and genotyping of DNA from fecal samples of rhinos in Indonesia
To ensure that reliable and accurate dung genotyping of wild Sumatran rhinoceros can be conducted in the range country of Sumatran rhinos, we completed similar evaluations of DNA extraction protocols and PCR optimization steps at UIUC and EIMB. At EIMB, the Qiagen Multiplex PCR kit proved useful for amplifying low-quality DNA from dung samples collected at SRS and BBS (Additional file 1: Figure S1), with the PCR results evaluated for each marker for: shape of peaks, secondary peaks, DNA slippage bands, intensity in relative fluorescent units, and ease of allele identification. These observations were considered in compiling a set of markers that is more highly recommended for use, listed first in Table 2.
Our bioinformatics pipeline saved time and resources, because the candidate microsatellite loci chosen for further evaluation were polymorphic in two sequenced Sumatran rhinos. Of a total 30,556,224 Illumina sequencing reads, 176,357 reads (2.4%) from rhino Dsu-33 and 167,849 reads (2.8%) from Dsu-35 contained microsatellite motifs. Our bioinformatics method identified 861 polymorphic di-, tri-, and tetra-nucleotide repeat loci within the two sequenced Sumatran rhinos of which 229 had suitable priming regions. Loci were then excluded if they displayed < 6 tandem repeats, had too broad a size range, showed evidence of being in repetitive elements, or closely matched sequences of human DNA (a potential contaminating factor).
The above resulted in 55 polymorphic candidate loci for UIUC laboratory PCR optimization using six test samples (Dsu-28, 29, 44, 63, 64 and 66; Table 1). From the optimization efforts, 53/55 primer pairs produced amplicons within the expected size range when tested using at least two samples; while 24/53 were excluded for having secondary bands outside of the expected size range, for genotypes that were difficult to score, or for being monomorphic in the test samples. The monomorphic loci in the test samples suggest that some polymorphisms detected in the two Illumina-sequenced rhinos (different individuals from the six rhinos used to characterize variability) were due to low-frequency alleles. The polymorphism of the 29 remaining markers (Table 2; Additional file 1: Table S2) was characterized across the test rhinos: average allelic richness (A) = 2.4; number of alleles ranged from two to four; average expected heterozygosity (HE) = 0.45; and average observed heterozygosity (HO) = 0.30. Overall, FIS = 0.44, likely due to population structure among geographically separated populations, or to inbreeding. Using the most informative markers, i.e., those with the lowest PID values, as few as seven loci could cumulatively distinguish individual identity (cumulative PID < 0.0001 ) (Table 2). Using a more conservative standard, PID(sib), as few as 19 optimized loci could confidently distinguish individual identity among siblings (cumulative PID(sib) < 0.0001 ; Table 2).
DNA from the dung of the two rhinos at the Cincinnati Zoo was used at UIUC to optimize PCR mix and cycling conditions (Additional file 1). At EIMB additional optimization was completed for the 29 markers (Additional file 1: Figure S1), with polymorphisms examined across three blood and three fecal samples from five SRS rhinos (for one SRS rhino, both types of sample were obtained). To conserve DNA resources, not all samples were used to test each marker. Yet 28/29 markers successfully amplified, and 27 of these were polymorphic in initial testing. For 13 markers randomly chosen among the 27, all 6 SRS samples were used to test marker quality and polymorphism. For some of these 13 markers, genotypes could be established for all of the SRS rhino samples (Additional file 1: Table S3). For 4/6 SRS samples, all 13 of the markers generated genotypes; for the other two (both fecal) samples, eight and twelve markers were successful (Additional file 1: Table S3). To examine the accuracy of genotypes from fecal samples, three paired blood-fecal samples were available (Table 1). While the number is low, comparing genotypes based on the two types of samples (Additional file 1) indicated that fecal samples yield accurate genotypes.
To further test their utility for a census using dung collected from wild Sumatran rhinos, the 13 randomly chosen new markers were tested on 11 fecal samples collected from an unknown number of wild rhinoceros in BBS (Table 1). Their identity as rhino samples was established by mtDNA sequencing because wild rhino and Asian tapir (Tapirus indicus) fecal samples are sometimes difficult to distinguish in the field. The specificity of primer pairs was tested using Asian tapir DNA and human DNA (a conceivable contaminant). One or both of these species amplified for 12/29 markers, although their amplicon sizes were almost always substantially different from those of rhino (Additional file 1: Table S4).
After species identity was established for 11 fecal samples from wild rhinos, the 13 randomly chosen markers were tested. For nine of the 13 markers, DNA from least nine of the 11 fecal samples amplified successfully and could be scored (Additional file 1: Table S3). For DNA from nine of the 11 fecal samples from wild Sumatran rhinos, at least 9 of the markers amplified successfully. This survey thus established proof of principle for the utility of the markers for fecal censusing of wild Sumatran rhinoceros.
Recent estimates suggest that fewer than 100 Sumatran rhinoceros individuals persist across Sumatra, with few in Borneo, although current census estimates have a large degree of uncertainty . Since an accurate census could guide management decisions, it is crucial to have reliable methods for non-invasively estimating population size. Overall, we generated a panel of polymorphic microsatellite markers useful for genetically distinguishing among individuals. They are appropriate for fecal DNA given their specificity for rhinos and short amplicon lengths. As almost all markers successfully amplified DNA from fecal samples, they will be useful for informing conservation managers about the population size and genetic characteristics of wild Sumatran rhinoceros.
As additional fecal samples are collected from wild rhinoceros, PCR conditions and marker choice may be further optimized for some primer pairs (Table 2). Although initial results suggest that the genotyping error rate is low, further testing of these markers using a larger number of paired blood and fecal samples is required to better estimate genotyping error (i.e., allelic dropout and false alleles) .
Availability of data and materials
Raw sequencing reads were deposited in the Sequence Read Archive at the NCBI: accession number SRR12548982 for Dsu-33 (BioSample accessions: SAMN15944577); and SRR12548981 for Dsu-35 (BioSample accessions: SAMN15944578). The datasets used and/or analyzed by the current study are available from the corresponding author on reasonable request. Sequences for each locus are shown in Additional file 1: Figure S2.
Burkit Barisan Selatan National Park
Eijkman Institute for Molecular Biology
Institute for Conservation Research (San Diego Zoo)
Royal Ontario Museum
Sumatran Rhino Sanctuary
University of Illinois at Urbana-Champaign
Way Kambas National Park
Pusparini W, Sievert PR, Fuller TK, Randhir TO, Andayani N. Rhinos in the parks: an island-wide survey of the last wild population of the Sumatran rhinoceros. PLoS ONE. 2015;10:9.
Guschanski K, Vigilant L, McNeilage A, Gray M, Kagoda E, Robbins MM. Counting elusive animals: comparing field and genetic census of the entire mountain gorilla population of Bwindi Impenetrable National Park Uganda. Biol Conserv. 2009;142(2):290–300.
Goode MJ, Beaver JT, Muller LI, Clark JD, van Manen FT, Harper CA, Basinger PS. Capture-recapture of white-tailed deer using DNA from fecal pellet groups. Wildlife Biol. 2014;20(5):270–8.
Goldstein DB, Schlötterer C. Microsatellites: Evolution and Applications. Oxford: Oxford University Press; 1999.
Ishida Y, Oleksyk TK, Georgiadis NJ, David VA, Zhao K, Stephens RM, Kolokotronis SO, Roca AL. Reconciling apparent conflicts between mitochondrial and nuclear phylogenies in African elephants. PLoS ONE. 2011;6(6):e20642.
Zhan X, Li M, Zhang Z, Goossens B, Chen Y, Wang H, Bruford MW, Wei F. Molecular censusing doubles giant panda population estimate in a key nature reserve. Curr Biol. 2006;16(12):R451-452.
Arandjelovic M, Vigilant L. Non-invasive genetic censusing and monitoring of primate populations. Am J Primatol. 2018;80(3):e22743.
Ishida Y, Demeke Y. Van Coeverden de Groot PJ, Georgiadis NJ, Leggett KEA, Fox VE, Roca AL: Distinguishing forest and savanna African elephants using short nuclear DNA sequences. J Hered. 2011;102(5):610–6.
Ishida Y, Demeke Y. Short amplicon microsatellite markers for low quality elephant DNA. Conserv Genet Resour. 2012;4(2):491–4.
Taberlet P, Waits LP, Luikart G. Noninvasive genetic sampling: look before you leap. Trends Ecol Evol. 1999;14(8):323–7.
Scott C, Foose T, Morales JC, Fernando P, Melnick DJ, Boag PT, Davila JA. Optimization of novel polymorphic microsatellites in the endangered Sumatran rhinoceros (Dicerorhinus sumatrensis). Mol Ecol Notes. 2004;4:194–6.
Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215.
Archie EA, Moss CJ, Alberts SC. Characterization of tetranucleotide microsatellite loci in the African savannah elephant (Loxodonta africana africana). Mol Ecol Notes. 2003;3(2):244–6.
Faircloth BC. MSATCOMMANDER: detection of microsatellite repeat arrays and automated, locus-specific primer design. Mol Ecol Resour. 2008;8(1):92–4.
Slater GS, Birney E. Automated generation of heuristics for biological sequence comparison. Bmc Bioinform. 2005;6:1.
Goudet J. FSTAT (Version 1.2): A computer program to calculate F-statistics. J Heredity. 1995;86(6):485–6.
Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel Population genetic software for teaching and research-an update. Bioinformatics. 2012;28(19):2537–9.
Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16(5):1099–106.
Waits LP, Luikart G, Taberlet P. Estimating the probability of identity among genotypes in natural populations: cautions and guidelines. Mol Ecol. 2001;10(1):249–56.
Havmoller RG, Payne J, Ramono W, Ellis S, Yoganand K, Long B, Dinerstein E, Williams AC, Putra RH, Gawi J, et al. Will current conservation responses save the Critically Endangered Sumatran rhinoceros Dicerorhinus sumatrensis? Oryx. 2016;50(2):355–9.
McKelvey KS, Schwartz MK. Providing reliable and accurate genetic capture-mark-recapture estimates in a cost-effective way. J Wildl Manag. 2004;68(3):453–6.
We thank the governments of Indonesia and Malaysia for enabling this research. We thank the Association and Zoos & Aquariums Rhino Advisory Group and the Rhino Research Council for enabling the collection of samples, and the USFWS inspectors at the Memphis Wildlife Inspection Office for guidance during the import process from Canada to the USA. For help in obtaining samples, we thank (in addition to co-authors) the Directorate General, Conservation of Natural Resources and Environment (KSDAE), Ministry of Environment and Forestry of the Republic of Indonesia; Indonesia Institute of Science; Burkit Barisan Selatan National Park Management Authority; Dr. Terri Roth (Cincinnati Zoo); Dr. Cynthia Steiner, Marlys Houck, Suellen Charter, and Julie Fronczek (San Diego Zoo ICR); Dr. Burton K. Lim and Jacqueline Miller (ROM). We thank the UIUC High-Throughput Sequencing and Genotyping Unit.
For funding we thank the US Fish and Wildlife Service Rhinoceros and Tiger Conservation Fund (F15AP00758 and F18AP00853); the International Rhino Foundation; the Ministry of Research and Technology/National Research and Innovation Agency of the Republic of Indonesia; the World Wildlife Fund; the National Science and Engineering Research Council (Canada); and the UIUC ACES Office of International Program (Seed Grant and Graduate Student International Research Grant). JRB was supported by a UIUC Graduate College Dissertation Completion Fellowship.
Ethics approval and consent to participate
All research was conducted with UIUC IACUC approval (protocol # 15053). Samples were imported under appropriate permits required by CITES and national governments; CITES/ESA Permit 14US84465A/9, CITES COSE Permit 12US757718/9. Samples from Sumatra were transferred within Indonesia under appropriate permits required by Directorate General, Conservation of Natural Resources, and Environment of Lampung Province; 103/SATS-DN/BKSDA.L-09/XII/2012.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Brandt, J.R., Saidah, S.H., Zhao, K. et al. Characterization of 29 polymorphic microsatellite markers developed by genomic screening of Sumatran rhinoceros (Dicerorhinus sumatrensis). BMC Res Notes 14, 119 (2021). https://doi.org/10.1186/s13104-021-05522-x
- Dicerorhinus sumatrensis
- Short tandem repeats
- Non-invasive sampling