Cross-species testing and utility of microsatellite loci in Indirana frogs
© Nair et al.; licensee BioMed Central Ltd. 2012
Received: 16 May 2012
Accepted: 13 July 2012
Published: 29 July 2012
Microsatellite loci are widely used in population and conservation genetic studies of amphibians, but the availability of such markers for tropical and subtropical taxa is currently very limited. In order to develop resources for conservation genetic studies in the genus Indirana, we tested amplification success and polymorphism in 62 previously developed microsatellite loci, in eight Indirana species - including new candidate species. Developing genomic resources for this amphibian taxon is particularly important as it is endemic to the Western Ghats biodiversity hotspot, and harbours several endangered species.
The cross-species amplification success rate varied from 11.3 % to 29.0 % depending on the species, with 29 - 80 % of the amplifying loci being polymorphic. A strong negative correlation between cross-species amplification success (and polymorphism) and genetic distance separating target from source species was observed.
Our results provide additional genetic support for the existence of genetically divergent cryptic species within the genus Indirana. The tested markers should be useful for population and conservation genetic studies in this genus, and in particular, for species closely related to the source species, I. beddomii.
The fauna of the Indian Western Ghats biodiversity hotspot is well known for its diversity and high level of endemism . In particular, there are many families and genera of amphibians that are unique to the region [2, 3], with roughly 132 endemic species . However, the amphibian diversity in this region remains inadequately characterized . This is reflected by numerous taxonomic uncertainties and ambiguities , and in the fact that many species still await proper taxonomic description . Frogs belonging to the endemic genus Indirana (Ranixalidae; ) are among the poorly studied amphibian genera from the Western Ghats . At present, information regarding their interspecific and intraspecific differentiation and variability across the diverse and fragmented habitats of Western Ghats is very limited .
An earlier study based on the examination of interspecific mitochondrial and nuclear sequence variability revealed a great deal of cryptic diversity within the Indirana genus; morphologically similar species displayed a high degree (4.2–17.1 %) of genetic divergence . In order to further explore this differentiation, we have recently developed 62 polymorphic microsatellite loci for one of the Indirana species (I. beddomii; ). Here we report the results of cross-species amplification tests for these microsatellite loci in eight other species (including candidate species) within the Indirana genus. Apart from providing information about the cross-species utility of these markers, we also investigated how the amplification success and levels of polymorphism relate to the evolutionary divergence between target and source species.
The samples were collected during field surveys in the southern Western Ghats between 2008 and 2010, under licence from National Biodiversity Authority, India (licence #NBA/TECH Appl/9/85/34/08/08-09/682).
The 62 microsatellite markers to be tested were developed using samples of I. beddomii from southern Kerala (08°45’59”N, 77°06’34”E; ). However, specimens thought to be I. beddomii from sites in northern and central Kerala (Aralam 11°55’54”N, 75°50’09”E; and Periyar 09°29’27”N, 77°08’10”E) and Kudremukh in Karnataka (13°12’38”N, 75°11’19”E) have previously shown high genetic divergence (4.2-12.5 %) from I. beddomii frogs from southern Kerala, and from each other . This high genetic divergence suggested that these were cryptic species albeit morphologically similar to each other and to I. beddomii. Therefore, I. beddomii as currently recognized appears not to be a single species, but a complex of at least four distinct species . The same also seems to apply to another species, I. diplosticta. We tested for cross-species microsatellite amplification in four putative species, along with three known species in the genus (I. semipalmata I. diplosticta I. leptodactyla), and in an additional unknown species (Indirana sp., Additional file 1). Each primer pair developed for I. beddomii (the “source” species) was tested against all of the “target” species and scored for amplification success and polymorphism. Cross-species amplification tests were performed using 5–10 individuals for each of the eight species/new candidate species within the genus Indirana (Additional file 1).
DNA was extracted from toe clips using a DNeasy Blood and Tissue kit (QIAGEN). Forward primers were labelled with fluorescent dyes FAM, HEX and TET (DNA Technology A/S), and a 5'-GTTT “tail” was added to every reverse primer to facilitate accurate genotyping . PCR reactions were performed for each primer pair in a total reaction volume of 10μL consisting of 1× QIAGEN Multiplex PCR solution, 0.2-0.3 μM of primers, dH2O and 10–20 ng of template DNA. The following PCR cycling conditions were used for amplification: 95 °C for 15 min, followed by 30 cycles of 95 °C for 30 s, 55 °C for 1.5 min, 72 °C for 1 min, and then a final extension step at 60 °C for 10 min. The PCR products were diluted 1:100 and electrophoresed on a MegaBACE 1000 capillary sequencer with MegaBACE ET550-R size standard. Genotypes were scored using Fragment Profiler (ver. 1.2; GE Healthcare Life Sciences). The amplification was considered positive if one or two alleles were observed with little or no stutter in comparison to a positive control. Decreasing the PCR annealing temperature is known to increase amplification success and polymorphism detection in cross-species amplification tests . Therefore, markers that did not amplify initially were re-tested for amplification at 52 °C. Annealing temperatures lower than 52 °C were not used, as they have not proved to be useful in identifying polymorphic loci in earlier cross-species tests .
The mitochondrial 16S ribosomal RNA gene sequences (466 bp) of Indirana species were retrieved from Genbank [Genbank: JQ596642-44, JQ596648-85]. We used 3–5 sequences per species to estimate the genetic distance between species (i.e. source-target species divergence) based on Kimura two parameter model  using MEGA5 . All statistical tests were performed using the program SPSS 15 (SPSS Inc., Chicago IL, USA). The species identity of the individuals used in the study was ascertained using 16S sequence as barcodes, following methods described in Nair et al. .
Cross-species amplification success rate of the 62 microsatellite loci tested on Indirana species
No. of individuals tested
No. of amplifying loci
No. of polymorphic loci
% of amplifying loci which are polymorphic
I. beddomii * (Aralam)
I. beddomii* (Kudremukh)
I. diplosticta* (Vellarimala)
I. beddomii * (Periyar)
It is well known that the utility of the microsatellites developed for a particular species is a negative function of the genetic distance separating the target and source species [13, 17, 18], which has been shown in birds, cetaceans, frogs  and other amphibians (e.g. ). Our results from the Indirana genus conform to these general patterns, and the observed cross-species amplification success rate of 11.3-29.0 % (mean = 21.2 %) is comparable to the within-genus amplification success rate of approximately 21 % reported for Ranid frogs . The fact that we found relatively low cross-species amplification success rates of 21.0 %, 25.8 %, and 25.8 % for the three cryptic species that were previously thought to be I. beddomii, supports the idea – based on analyses of multiple nuclear and mitochondrial genes  – that the Aralam, Kudremukh and Periyar populations of what was thought to be I. beddomii are indeed distinct species .
One of the main reasons for lowered cross-species amplification success is thought to be the lack of conservation of priming sites between highly divergent species . Additionally, increased genome size (C-value) is also known to have a negative effect on cross-species amplification success [22, 23]: a decrease in the ratio of target to non-target DNA causes a reduction in amplification efficiency [24, 25]. Regardless of the proximate cause, the close correspondence between the proportion of amplifying loci and 16S divergence among the tested taxa (Figure 1a) – as well as the proportion of polymorphic loci and 16S divergence (Figure 1b & c) – suggests comparable rates of divergence of mitochondrial and nuclear genomes in these frogs. While these patterns are admittedly driven by low amplification success in the two highly divergent species, I. diplosticta and I. leptodactyla (14.7-17.5 %), the data nonetheless suggests that the cross-species utility of these markers declines with increasing evolutionary distance among taxa. For example, only 2–3 out the 62 tested loci were found to be polymorphic in the two most divergent taxa, limiting the utility of these markers in studies of divergent Indirana species. Interestingly, one of the loci (IND107) was highly polymorphic in all taxa studied, with similar numbers of alleles observed both in the source and most divergent species (Additional file 1). This locus appears to be conserved in all species within Indirana and could be associated with some functionally important gene. Conservation of microsatellite loci residing within or close to functionally important genes has been reported from other highly divergent species (e.g. [26, 27]).
Taken together, the results of this study show that the cross-species amplification success in Indirana frogs depends on the degree of evolutionary divergence between the source and target species, with the success rate and polymorphism declining rapidly with increasing divergence between the taxa. However, given the relatively high levels of microsatellite polymorphism in some of the target species, these markers may provide useful tools for future conservation genetic studies aiming to address taxonomic uncertainties, or to study genetic variability and differentiation of Indirana populations.
AN, SG and JM conceived the study. SVG, SG and KSK collected all the samples. AN performed all the experiments, analyses and had a major role in writing the manuscript together with JM. All authors contributed to and approved the final manuscript.
We thank Marika Karjalainen and Kirsi Kähkönen for assistance in laboratory and Jacquelin DeFaveri for language editing. We are grateful to Kerala and Karnataka Forest Departments for permissions for sample collections. Thanks are due to Vinodlal K. S., Vijayan P. and Biju S. Nair for help with the field work, and Amber Teacher and Takahito Shikano for discussions and helpful comments on earlier versions of this manuscript. This study was funded by Academy of Finland (grant #125092 to JM) and Department of Biotechnology, Government of India (BT/IN/Finnish/11/SG/2008 to SG).
- Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J: Biodiversity hotspots for conservation priorities. Nature. 2000, 403: 853-858. 10.1038/35002501.PubMedView ArticleGoogle Scholar
- Biju SD, Bossuyt F: New frog family from India reveals an ancient biogeographical link with the Seychelles. Nature. 2003, 425: 711-714. 10.1038/nature02019.PubMedView ArticleGoogle Scholar
- Daniels RJR: Amphibians of Peninsular India. 2005, Universities Press (India) Private Ltd, HyderabadGoogle Scholar
- Dinesh KP, Radhakrishnan C, Gururaja KV, Bhatt G: An annotated checklist of Amphibia of India with some insights into the patterns of species discoveries, distribution and endemism. Rec Zool Surv India. 2009, 302: 1-133.Google Scholar
- Inger RF: Distribution of amphibians in southern Asia and adjacent islands. Patterns of Distribution of Amphibians: A global perspective. Edited by: Duellman WE. 1999, The Johns Hopkins University Press, London, 445-482.Google Scholar
- Daniels RJR: Taxonomic uncertainties and conservation assessment of the Western Ghats. Curr Sci. 1997, 73: 169-170.Google Scholar
- Biju SD: A synopsis to the frog fauna of the Western Ghats. India. Occas Publ Indian Soc Conserv Biol. 2001, 1: 1-24.Google Scholar
- Bocxlaer IV, Roelants K, Biju SD, Nagaraju J, Bossuyt F: Late Cretaceous vicariance in Gondwanan amphibians. PLoS One. 2006, 1: e74-10.1371/journal.pone.0000074.PubMedPubMed CentralView ArticleGoogle Scholar
- Nair A, Gopalan SV, George S, Kumar KS, Teacher AGF, Merilä J: Endemic Indirana frogs of the Western Ghats biodiversity hotspot. Ann Zool Fenn, in pressGoogle Scholar
- Nair A, Gopalan SV, George S, Kumar KS, Teacher AGF, Merilä J: High cryptic diversity of endemic Indirana frogs in the Western Ghats biodiversity hotspot. Anim Conserv. doi: 10.1111/j.1469-1795.2012.00539.x.Google Scholar
- Nair A, Kumar KS, George S, Gopalan SV, Li M-H, Leder E, Merilä J: Sixty two new microsatellite markers for an endemic frog Indirana beddomii from the Western Ghats biodiversity hotspot. Conserv Genet Resour. 2011, 3: 167-171. 10.1007/s12686-010-9315-1.View ArticleGoogle Scholar
- Brownstein MJ, Carpten JD, Smith JR: Modulation of non-templated nucleotide addition by TaqDNA polymerase: primer modification that facilitate genotyping. BioTechiques. 1996, 20: 1004-1010.Google Scholar
- Primmer CR, Painter JN, Koskinen MT, Palo JU, Merilä J: Factors affecting avian cross-species microsatellite amplification. J Avian Biol. 2005, 36: 348-360. 10.1111/j.0908-8857.2005.03465.x.View ArticleGoogle Scholar
- Liu Z, Crooijmans RPMA, van der Poel JJ, Groenen MAM: Use of chicken microsatellite markers in turkey: A pessimistic view. Anim Genet. 1996, 27: 191-193.View ArticleGoogle Scholar
- Kimura MA: A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Steker G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Moore SS, Sargeant LL, King TJ, Mattick JS, Georges M, Hetzel DJS: The conservation of dinucleotide microsatellites among mammalian genomes allows the use of heterologous PCR primer pairs in closely related species. Genomics. 1991, 10: 654-660. 10.1016/0888-7543(91)90448-N.PubMedView ArticleGoogle Scholar
- Carreras-Carbonell J, Macpherson E, Pascual M: Utility of pairwise mtDNA genetic distances for predicting cross-species microsatellite amplification and polymorphism success in fishes. Conserv Genet. 2008, 9: 181-190. 10.1007/s10592-007-9322-2.View ArticleGoogle Scholar
- Hendrix R, Hauswaldt S, Veith M, Steinfartz S: Strong correlation between cross-amplification success and genetic distance across all members of “True Salamanders” (Amphibia: Salamandridae) revealed by Salamandra salamandra-specific microsatellite loci. Mol Ecol Resour. 2010, 10: 1038-1047. 10.1111/j.1755-0998.2010.02861.x.PubMedView ArticleGoogle Scholar
- Primmer CR, Merilä J: A low rate of cross-species microsatellite amplification success in ranid frogs. Conserv Genet. 2002, 3: 445-449. 10.1023/A:1020519218836.View ArticleGoogle Scholar
- Balloux F, Ecoffey E, Fumagalli L, Goudet J, Wyttenbach A, Hausser F: Microsatellite conservation, polymorphism, and GC content in shrews of the genus Sorex (Insectivora, Mammalia). Mol Biol Evol. 1998, 15: 473-475. 10.1093/oxfordjournals.molbev.a025945.PubMedView ArticleGoogle Scholar
- Garner TWJ, Schmidt BR, Hoeck P, Van Buskirk J: Di- and tetranucleotide microsatellite markers for the Alpine newt (Triturus alpestris): characterization and cross-priming in five congeners. Mol Ecol Notes. 2003, 3: 186-188. 10.1046/j.1471-8286.2003.00394.x.View ArticleGoogle Scholar
- Barbará T, Palma-Silva C, Paggi GM, Bered F, Fay MF, Lexer C: Cross-species transfer of nuclear microsatellite markers: potential and limitations. Mol Ecol. 2007, 16: 3759-3767. 10.1111/j.1365-294X.2007.03439.x.PubMedView ArticleGoogle Scholar
- Farrelly V, Rainey FA, Stackebrandt E: Effect of genome size and rrn gene copy number on PCR amplification of 16 S rRNA genes from a mixture of bacterial species. Appl Environ Microbiol. 1995, 61: 2798-2801.PubMedPubMed CentralGoogle Scholar
- Garner TWJ: Genome size and microsatellites: the effect of nuclear size on amplification potential. Genome. 2002, 45: 212-215. 10.1139/g01-113.PubMedView ArticleGoogle Scholar
- Moore SS, Hale P, Bryne K: NCAM: a polymorphic microsatellite locus conserved across eutherian mammal species. Anim Genet. 1998, 29: 33-36. 10.1046/j.1365-2052.1998.00234.x.PubMedView ArticleGoogle Scholar
- Shikano T, Ramadevi J, Shimada Y, Merilä J: Utility of sequenced genomes for microsatellite marker development in non-model organisms: a case study of functionally important genes in nine-spined sticklebacks (Pungitius pungitius). BMC Genomics. 2010, 11: 334-10.1186/1471-2164-11-334.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.