Rapid Development of Microsatellite Markers with 454 Pyrosequencing in a Vulnerable Fish, the Mottled Skate, Raja pulchra

The mottled skate, Raja pulchra, is an economically valuable fish. However, due to a severe population decline, it is listed as a vulnerable species by the International Union for Conservation of Nature. To analyze its genetic structure and diversity, microsatellite markers were developed using 454 pyrosequencing. A total of 17,033 reads containing dinucleotide microsatellite repeat units (mean, 487 base pairs) were identified from 453,549 reads. Among 32 loci containing more than nine repeat units, 20 primer sets (62%) produced strong PCR products, of which 14 were polymorphic. In an analysis of 60 individuals from two R. pulchra populations, the number of alleles per locus ranged from 1-10, and the mean allelic richness was 4.7. No linkage disequilibrium was found between any pair of loci, indicating that the markers were independent. The Hardy-Weinberg equilibrium test showed significant deviation in two of the 28 single-loci after sequential Bonferroni's correction. Using 11 primer sets, cross-species amplification was demonstrated in nine related species from four families within two classes. Among the 11 loci amplified from three other Rajidae family species; three loci were polymorphic. A monomorphic locus was amplified in all three Rajidae family species and the Dasyatidae family. Two Rajidae polymorphic loci amplified monomorphic target DNAs in four species belonging to the Carcharhiniformes class, and another was polymorphic in two Carcharhiniformes species.     

Isolation and characterization of cross-amplification microsatellite panels for species of procapra (bovidae; antilopinae)

The three Procapra species, Tibetan gazelle (P. picticaudata), Mongolian gazelle (P. gutturosa) and Przewalski's gazelle (P. przewalskii) are endemic to Asia. Several intraspecific genetic issues have been studied with species-specific microsatellite loci in these Asian gazelles. However, cross-species microsatellite panels are absent, which inhibits comparative conservation and evolutionary studies of the Procapra. In this study, we isolated 20 cross-species microsatellite loci for Procapra from both related species and the genomic library of P. przewalskii. Fifty-three samples of the three gazelles were used to characterize the markers. Allele numbers ranged from three to 20, with a mean of 7.93 per locus. Observed heterozygosity (H(O)) averaged 0.680 and expected heterozygosity (H(E)) 0.767. The mean polymorphic information content (PIC) was 0.757 for P. picticaudata, 0.803 for P. gutturosa and 0.590 for P. przewalskii. Nine loci were significantly deviated from Hardy-Weinberg (H-W) equilibrium in the three species. Significant linkage disequilibrium was detected in four pairs of loci in P. przewalskii, five pairs in P. gutturosa and 51 pairs in P. picticaudata. Considering the abundance of published loci and their high success rates of cross-amplification, testing and utilization of loci from related species is efficient for wild species of Bovidae. The cross-species microsatellite loci we developed will facilitate further interspecies genetic studies in Procapra.     

Isolation and Characterization of Polymorphic Microsatellite Loci from Metapenaeopsis barbata Using PCR-Based Isolation of Microsatellite Arrays (PIMA)

The red-spot prawn, Metapenaeopsis barbata, is a commercially important, widely distributed demersal species in the Indo-West Pacific Ocean. Overfishing has made its populations decline in the past decade. To study conservation genetics, eight polymorphic microsatellite loci were isolated. Genetic characteristics of the SSR (simple sequence repeat) fingerprints were estimated in 61 individuals from adjacent seas of Taiwan and China. The number of alleles, ranging from 2 to 4, as well as observed and expected heterozygosities in populations, ranging from 0.048 to 0.538, and 0.048 and 0.654, respectively, were detected. No deviation from Hardy-Weinberg expectations was detected at either locus. No significant linkage disequilibrium was detected in locus pairs. The polymorphic microsatellite loci will be useful for investigations of the genetic variation, population structure, and conservation genetics of this species.     

Eleven Novel Polymorphic Microsatellite Loci for Oval Squid Sepioteuthis Lessoniana (Shiro-Ika Type)

The oval squid Sepioteuthis lessoniana is one of the most economically important squid species in Japan; however, its population structure is poorly understood due to the lack of hypervariable markers. Such information is critical for managing sustainable fisheries, as well as for ensuring the existence of wild S. lessoniana stocks. Eleven candidate microsatellite loci were isolated from a small insert genomic DNA library. Polymorphisms in these 11 loci were screened in 24 wild individuals. The number of alleles per locus was found to range from 5 to 19 alleles, and the observed heterozygosity ranged from 0.292 to 0.958. No evidence for linkage disequilibrium was detected among all the loci. The genotypic proportions conformed to Hardy-Weinberg equilibrium, except at one locus. In conclusion, these polymorphic microsatellite loci may be used to develop a genetic framework to manage S. lessoniana in the future.     

Isolation and Characterization of Sixteen Polymorphic Microsatellite Loci in the Golden Apple Snail Pomacea canaliculata

We report the characterization of 16 polymorphic microsatellite markers in the golden apple snail, Pomacea canaliculata, a pest registered in the list of "100 of the world's worst invasive alien species". The fast isolation by AFLP (Amplified Fragment Length Polymorphism) of sequences containing repeats (FIASCO) method was used to isolate microsatellite loci, and polymorphism was explored with 29 individuals collected in an invasive region from China. These primers showed a number of alleles per locus ranging from three to 13. The ranges of observed and expected heterozygosity were 0.310-0.966 and 0.523-0.898, respectively. These microsatellite markers described here will be useful for population genetic studies of P. canaliculata.     

Polymorphic Microsatellite Loci Isolated from the Squalidus argentatus Using PCR-Based Isolation of Microsatellite Arrays (PIMA)

Squalidus argentatus (Sauvage and Dabry de Thiersant 1874) is a small-sized freshwater fish which is distributed in Mainland China, Hainan Island and Taiwan. The populations of S. argentatus have dropped sharply probably due to overharvesting and water pollution recently. Eleven polymorphic microsatellite markers were developed for the cyprinid fish S. argentatus. These new markers were tested on 43 individuals collected from Yangtze River and Qiantang River. The number of alleles, observed and expected heterozygosity per locus, in two populations ranged from 3 to 14, from 0.333 to 0.954 and from 0.480 to 0.928, respectively. Only two loci are significantly deviated from Hardy-Weinberg expectations due to the heterozygote deficiency. No significant linkage disequilibrium was detected between the pairwise comparisons of these loci. These polymorphic microsatellite loci will enable us to study the genetic variation, population structure, and conservation genetics of this species in the future.     

Isolation and Characterization of Eleven Polymorphic Microsatellite Loci for the Valuable Medicinal Plant Dendrobium huoshanense and Cross-Species Amplification

Dendrobium huoshanense (Orchidaceae) is a perennial herb and a widely used medicinal plant in Traditional Chinese medicine (TCM) endemic to Huoshan County town in Anhui province in Southeast China. A microsatellite-enriched genomic DNA library of D. huoshanense was developed and screened to identify marker loci. Eleven polymorphic loci were isolated and analyzed by screening 25 individuals collected from a natural population. The number of alleles per locus ranged from 2 to 5. The observed and expected heterozygosities ranged from 0.227 to 0.818 and from 0.317 to 0.757, respectively. Two loci showed significant deviations from Hardy-Weinberg equilibrium and four of the pairwise comparisons of loci revealed linkage disequilibrium (p < 0.05). These microsatellite loci were cross-amplified for five congeneric species and seven loci can be amplified in all species. These simple sequence repeats (SSR) markers are useful in genetic studies of D. huoshanense and other related species and in conservation decision-making.     

Nuclear Microsatellite Primers for the Endangered Relict Fir,Abies pinsapo (Pinaceae) and Cross-Amplification in Related Mediterranean Species

Twelve nuclear microsatellite primers (nSSR) were developed for the endangered species Abies pinsapo Boiss. to enable the study of gene flow and genetic structure in the remaining distribution areas. Microsatellite primers were developed using next-generation sequencing (454) data from a single Abies pinsapo individual. Primers were applied to thirty individuals from the three extant localities. The number of alleles per locus ranged from one to four. Cross-amplification was tested for other Abies species from the Mediterranean Basin, and most of the loci showed higher polymorphisms in the Mediterranean species than in A. pinsapo. These microsatellite markers provide tools for conservation genetic studies in Abies pinsapo as well other Abies species from the Mediterranean Basin.     

Isolation and Characterization of Microsatellite Loci for Hibiscus aridicola (Malvaceae), an Endangered Plant Endemic to the Dry-Hot Valleys of Jinsha River in Southwest China

Hibiscus aridicola (Malvaceae) is an endangered ornamental shrub endemic to the dry-hot valleys of Jinsha River in southwest China. Only four natural populations of H. aridicola exist in the wild according to our field investigation. It can be inferred that H. aridicola is facing a very high risk of extinction in the wild and an urgent conservation strategy is required. By using a modified biotin-streptavidin capture method, a total of 40 microsatellite markers were developed and characterized in H. aridicola for the first time. Polymorphisms were evaluated in 39 individuals from four natural populations. Fifteen of the markers showed polymorphisms with two to six alleles per locus; the observed heterozygosity ranged from 0.19 to 0.72. These microsatellite loci would be useful tools for population genetics studies on H. aridicola and other con-generic species which are important to the conservation and development of endangered species.     

Isolation and Characterization of New 24 Microsatellite DNA Markers for Golden Cuttlefish (Sepia esculenta)

Twenty-four microsatellite DNA markers were isolated and characterized for golden cuttlefish (Sepia esculenta) from a (GT)(13)-enriched genomic library. Loci were tested in 48 individuals from Jiaozhou bay of China. The numbers of alleles per locus ranged from two to 25 with an average of 10.3. The observed and expected heterozygosities ranged from 0.063 to 0.896 and from 0.137 to 0.953, with averages of 0.519 and 0.633, respectively. Six loci significantly deviated from Hardy-Weinberg equilibrium after Bonferroni's correction and no significant linkage disequilibrium between loci pairs was detected. These microsatellite markers would be useful for analyzing the population genetic structure to make conservation and management decisions for S. esculenta.     

Rapid Development of Microsatellite Markers for the Endangered Fish Schizothorax biddulphi (Günther) Using Next Generation Sequencing and Cross-Species Amplification

Tarim schizothoracin (Schizothorax biddulphi) is an endemic fish species native to the Tarim River system of Xinjiang and has been classified as an extremely endangered freshwater fish species in China. Here, we used a next generation sequencing platform (ion torrent PGM™) to obtain a large number of microsatellites for S. biddulphi, for the first time. A total of 40577 contigs were assembled, which contained 1379 SSRs. In these SSRs, the number of dinucleotide repeats were the most frequent (77.08%) and AC repeats were the most frequently occurring microsatellite, followed by AG, AAT and AT. Fifty loci were randomly selected for primer development; of these, 38 loci were successfully amplified and 29 loci were polymorphic across panels of 30 individuals. The H_o ranged from 0.15 to 0.83, and H_e ranged from 0.15 to 0.85, with 3.5 alleles per locus on average. Cross-species utility indicated that 20 of these markers were successfully amplified in a related, also an endangered fish species, S. irregularis. This study suggests that PGM™ sequencing is a rapid and cost-effective tool for developing microsatellite markers for non-model species and the developed microsatellite markers in this study would be useful in Schizothorax genetic analysis.     

Development of Novel Microsatellite Markers in the Omei Treefrog (Rhacophorus omeimontis)

Eleven novel microsatellite markers were developed and characterized for the Omei treefrog (Rhacophorus omeimontis) using the fast isolation by AFLP of sequences containing repeats method. Polymorphism of each locus was tested in 24 individuals from two wild populations. The number of alleles per locus ranged from 4 to 15, the average observed and expected heterozygosity per locus ranged from 0.250 to 0.839 and from 0.562 to 0.914, respectively. Two of the 11 microsatellite loci showed significant deviations from Hardy-Weinberg equilibrium. Two locus pairs showed significant linkage disequilibrium. Neither evidence of scoring error due to stuttering nor evidence of large allele dropout was found at all of the 11 loci, but evidence of null alleles was indicated at two loci because of general excess of homozygotes for most allele size classes. These polymorphic loci will be useful markers in studying mate choice of the Omei treefrog.     

Cross-Species,Amplifiable Microsatellite Markers for Neoverrucid Barnacles from Deep-Sea Hydrothermal Vents Developed Using Next-Generation Sequencing

Barnacles of the genus Neoverruca are abundant near deep-sea hydrothermal vents of the northwestern Pacific Ocean, and are useful for understanding processes of population formation and maintenance of deep-sea vent faunas. Using next-generation sequencing, we isolated 12 polymorphic microsatellite loci from Neoverruca sp., collected in the Okinawa Trough. These microsatellite loci revealed 2–19 alleles per locus. The expected and observed heterozygosities ranged from 0.286 to 1.000 and 0.349 to 0.935, respectively. Cross-species amplification showed that 9 of the 12 loci were successfully amplified for Neoverruca brachylepadoformis in the Mariana Trough. A pairwise F_ST value calculated using nine loci showed significant genetic differentiation between the two species. Consequently, the microsatellite markers we developed will be useful for further population genetic studies to elucidate genetic diversity, differentiation, classification, and evolutionary processes in the genus Neoverruca.     

Development of 30 Novel Polymorphic Expressed Sequence Tags (EST)-Derived Microsatellite Markers for the Miiuy Croaker, Miichthys miiuy

Expressed sequence tags (ESTs) can be used to identify microsatellite markers. We developed 30 polymorphic microsatellite markers from 5053 ESTs of the Miichthys miiuy. Out of 123 EST derived microsatellites for which PCR primers were designed, 30 loci were polymorphic in 30 individuals from a single natural population with 2-13 alleles per locus. The observed and expected heterozygosities were from 0.1024 to 0.7917 and from 0.2732 to 0.8845, respectively. Nine loci deviated from the Hardy-Weinberg equilibrium, and linkage disequilibrium was significant between 22 pairs of loci. These polymorphic microsatellite loci will be useful for genetic diversity analysis and molecule-assisted breeding for M. miiuy.     

Isolation and Characterization of Microsatellite Loci in the Chinese Cobra Naja atra (Elapidae)

We characterize thirteen polymorphic microsatellite loci isolated from Naja atra genomic libraries, which were enriched for AC-motif microsatellites. The thirteen loci were screened on a group of 48 individuals from two populations, one in Yong'an and the other in Ganzhou. These markers revealed a relatively high degree of genetic diversity (4-12 alleles per locus) and heterozygosity (Ho ranged from 0.213-0.854 and He ranged from 0.301-0.838). Tests for departure from Hardy-Weinberg equilibrium and for linkage disequilibrium were conducted for each of the two populations separately. After sequential Bonferroni correction, none of the 13 loci showed significant departures from Hardy-Weinberg equilibrium. Hierarchical analysis of molecular variance indicated that a small but significant (P < 0.001) proportion (16.0%) of the total variation in the microsatellite DNA data were attributable to differences among populations, indicating geographical structuring and restricted gene flow. It could be attributable to the Wuyi mountains in the area having a sufficiently isolating effect to significantly reduce gene flow. Our microsatellite data also showed a low N(m) (1.31) value in the two populations from mainland China. Thus, the Yong'an and Ganzhou populations could be treated as distinct evolutionarily significant units (ESUs). The high level of polymorphism revealed by these microsatellite markers will be useful for the study of gene flow, population structure and evolutionary history of N. atra.     

Development of Microsatellite Markers for the Korean Mussel, Mytilus coruscus (Mytilidae) Using Next-Generation Sequencing

Mytilus coruscus (family Mytilidae) is one of the most important marine shellfish species in Korea. During the past few decades, this species has become endangered due to the loss of habitats and overfishing. Despite this species' importance, information on its genetic background is scarce. In this study, we developed microsatellite markers for M. coruscus using next-generation sequencing. A total of 263,900 raw reads were obtained from a quarter-plate run on the 454 GS-FLX titanium platform, and 176,327 unique sequences were generated with an average length of 381 bp; 2569 (1.45%) sequences contained a minimum of five di- to tetra-nucleotide repeat motifs. Of the 51 loci screened, 46 were amplified successfully, and 22 were polymorphic among 30 individuals, with seven of trinucleotide repeats and three of tetranucleotide repeats. All loci exhibited high genetic variability, with an average of 17.32 alleles per locus, and the mean observed and expected heterozygosities were 0.67 and 0.90, respectively. In addition, cross-amplification was tested for all 22 loci in another congener species, M. galloprovincialis. None of the primer pairs resulted in effective amplification, which might be due to their high mutation rates. Our work demonstrated the utility of next-generation 454 sequencing as a method for the rapid and cost-effective identification of microsatellites. The high degree of polymorphism exhibited by the 22 newly developed microsatellites will be useful in future conservation genetic studies of this species.     

Isolation and Characterization of Microsatellite Loci in Pistacia weinmannifolia (Anacardiaceae)

Fourteen polymorphic microsatellite loci were isolated from the genomic DNA of Pistacia weinmannifolia, using the Fast Isolation by AFLP of Sequences Containing repeats (FIASCO) method, and screened on 12 individuals from each of two wild populations. The 14 polymorphic loci had an average of 4.1 alleles per locus varying from 1 to 9. The observed (H(o)) and expected (H(e)) heterozygosities across the two populations ranged from 0.000 to 0.933 and from 0.000 to 0.906, respectively. Tests for departure from Hardy-Weinberg equilibrium (HWE) and genotypic linkage disequilibrium (LD) were conducted for each of the two populations separately. It was found that no locus significantly deviated from HWE proportions and no significant LD was detected between loci (p < 0.001). In the test of cross-species utility, we successfully amplified nine (64.2%) of 14 loci in P. chinensis and four (28.6%) in P. mexicana. The relatively high level of polymorphism for these markers will facilitate further studies of gene flow, population structure and evolutionary history of P. weinmannifolia and its congeners.     

A Set of Novel Microsatellite Markers Developed for a Distylous Species Luculia gratissima (Rubiaceae)

Luculia gratissima (Wall.) Sweet (Rubiaceae) is a perennial shrub distributed in the southeast margin of the Tibetan Plateau in southwest China and adjacent region of Nepal and Myanmar. The plant is a distylous species with reciprocally placed stigmas and anthers in each floral morph. By using the Fast Isolation by Amplified Fragment Length Polymorphism (AFLP) of Sequences Containing (FIASCO) repeats protocol, 19 primer sets were identified in two wild populations. Of these primers, 10 displayed polymorphisms and nine were monomorphic. The number of alleles per locus ranged from two to five, values for observed and expected heterozygosities ranged from 0.000 to 1.000 and from 0.289 to 0.760, with averages of 0.303 and 0.555, respectively. These microsatellite loci will facilitate further studies on breeding system, gene flow patterns, and population structure of L. gratissima and its allied species.     

Isolation and Characterisation of 11 Polymorphic Microsatellite Markers in Papaver rhoeas L. (Corn Poppy), a Major Annual Plant Species from Cultivated Areas

Papaver rhoeas, an annual plant species in the Papaveraceae family, is part of the biodiversity of agricultural ecosystems and also a noxious agronomic weed. We developed microsatellite markers to study the genetic diversity of P. rhoeas, using an enriched microsatellite library coupled with 454 next-generation sequencing. A total of 13,825 sequences were obtained that yielded 1795 microsatellite loci. After discarding loci with less than six repeats of the microsatellite motif, automated primer design was successful for 598 loci. We tested 74 of these loci for amplification with a total of 97 primer pairs. Thirty loci passed our tests and were subsequently tested for polymorphism using 384 P. rhoeas plants originating from 12 populations from France. Of the 30 loci, 11 showed reliable polymorphism not affected by the presence of null alleles. The number of alleles and the expected heterozygosity ranged from 3 to 7.4 and from 0.27 to 0.73, respectively. A low but significant genetic differentiation among populations was observed (F_ST = 0.04; p < 0.001). The 11 validated polymorphic microsatellite markers developed in this work will be useful in studies of genetic diversity and population structure of P. rhoeas, assisting in designing management strategies for the control or the conservation of this species.