Characterization of the complete mitochondrial genome and a set of polymorphic microsatellite markers through next-generation sequencing for the brown brocket deer Mazama gouazoubira

The complete mitochondrial genome of the brown brocket deer Mazama gouazoubira and a set of polymorphic microsatellite markers were identified by 454-pyrosequencing. De novo genome assembly recovered 98% of the mitochondrial genome with a mean coverage of 9-fold. The mitogenome consisted of 16,356 base pairs that included 13 protein-coding genes, two ribosomal subunit genes, 22 transfer RNAs and the control region, as found in other deer. The genetic divergence between the mitogenome described here and a previously published report was ∼0.5%, with the control region and ND5 gene showing the highest intraspecific variation. Seven polymorphic loci were characterized using 15 unrelated individuals; there was moderate genetic variation across most loci (mean of 5.6 alleles/locus, mean expected heterozygosity = 0.70), with only one locus deviating significantly from Hardy-Weinberg equilibrium, probably because of null alleles. Marker independence was confirmed with tests for linkage disequilibrium. The genetic variation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the phylogeography and population genetic patterns in M. gouazoubira, particularly in the context of habitat fragmentation in South America.     

Characterisation of 12 microsatellite loci in the Vietnamese commercial clam <Emphasis Type="Italic">Lutraria rhynchaena</Emphasis> Jonas 1844 (Heterodonta: Bivalvia: Mactridae) through next-generation sequencing

The marine clam Lutraria rhynchaena is gaining popularity as an aquaculture species in Asia. Lutraria populations are present in the wild throughout Vietnam and several stocks have been established and translocated for breeding and aquaculture grow-out purposes. In this study, we demonstrate the feasibility of utilising Illumina next-generation sequencing technology to streamline the identification and genotyping of microsatellite loci from this clam species. Based on an initial partial genome scan, 48 microsatellite markers with similar melting temperatures were identified and characterised. The 12 most suitable polymorphic loci were then genotyped using 51 individuals from a population in Quang Ninh Province, North Vietnam. Genetic variation was low (mean number of alleles per locus = 2.6; mean expected heterozygosity = 0.41). Two loci showed significant deviation from Hardy–Weinberg equilibrium (HWE) and the presence of null alleles, but there was no evidence of linkage disequilibrium among loci. Three additional populations were screened (n = 7–36) to test the geographic utility of the 12 loci, which revealed 100 % successful genotyping in two populations from central Vietnam (Nha Trang). However, a second population from north Vietnam (Co To) could not be successfully genotyped and morphological evidence and mitochondrial variation suggests that this population represents a cryptic species of Lutraria. Comparisons of the Qang Ninh and Nha Trang populations, excluding the 2 loci out of HWE, revealed statistically significant allelic variation at 4 loci. We reported the first microsatellite loci set for the marine clam Lutraria rhynchaena and demonstrated its potential in differentiating clam populations. Additionally, a cryptic species population of Lutraria rhynchaena was identified during initial loci development, underscoring the overlooked diversity of marine clam species in Vietnam and the need to genetically characterise population representatives prior to microsatellite development. The rapid identification and validation of microsatellite loci using next-generation sequencing technology warrant its integration into future microsatellite loci development for key aquaculture species in Vietnam and more generally, aquaculture countries in the South East Asia region.     

Social structure and landscape genetics of the endemic New Caledonian ant <Emphasis Type="Italic">Leptomyrmex pallens</Emphasis> Emery, 1883 (Hymenoptera: Formicidae: Dolichoderinae), in the context of fire-induced rainforest fragmentation

Habitat fragmentation is a major threat to biodiversity, as it can alter ecological processes at various spatial and trophic scales. At the species level, fragmentation leading to the isolation of populations can trigger reductions in genetic diversity, potentially having detrimental effects on population fitness, adaptability and ultimately population persistence. Leptomyrmex pallens is a widespread rainforest ant endemic to New Caledonia but now confined to habitat patches that have been fragmented by anthropogenic fire regimes over the last 200 years. We investigated the social structure of L. pallens in the Aoupinié region (c.a. 4900 ha), and assessed the impacts of habitat fragmentation on its population genetic structure. Allele frequencies at 13 polymorphic microsatellite loci were compared among 411 worker ants from 21 nests distributed across the region. High within-nest relatedness (r = 0.70 ± 0.02), and a single queen found in 38 % of the nests by pedigree analysis indicate that the species is monogynous to weakly polygynous. Estimates of gene flow and genetic structure across the region were subsequently determined using a combined dataset of single workers per nest and of unrelated foraging workers. These estimates coupled with a comprehensive landscape genetic analysis revealed no evidence of significant population structure or habitat effects, suggesting that the Aoupinié region harbours a single panmictic population. In contrast, analyses of mitochondrial DNA sequence data revealed a high degree of genetic structuring, indicating limited maternal gene flow and suggesting that gene flow among nests is driven primarily by winged males. Overall these findings suggest that fire-induced habitat fragmentation has had little impact on the population dynamics of L. pallens. Additional studies of less mobile species should therefore be conducted to gain further insights into fire related disturbances on the unique biodiversity and function of New Caledonian ecosystems.     

Isolation and characterization of microsatellite markers in a highland fish, <Emphasis Type="BoldItalic">Pareuchiloglanis sinensis</Emphasis> (Siluriformes: Sisoridae) by next-generation sequencing

Pareuchiloglanis sinensis (Siluriformes: Sisoridae) is an endemic and highland fish species which occurs only in some rivers of south-west China. In this study, the isolation and characterization of polymorphic microsatellite loci of this fish species by next-generation sequencing is described. A total of 9471 simple-sequence repeats (SSRs) were observed from RNA-seq data. One hundred and twenty primer pairs were chosen randomly and validated across 48 P. sinensis individuals collected from the Dadu river (a branch of the Yangtze river) of which 28 polymorphic microsatellite loci were detected. The number of alleles ranged from 2 to 14, with an average of seven alleles per locus. Twenty loci exhibited high polymorphism with the polymorphism information content (PIC) larger than 0.5. The mean observed and expected heterozygosity varied from 0.104 to 0.958 and 0.157 to 0.844, with an average of 0.583 and 0.613, respectively. The microsatellite markers characterized in the current study serve as a useful tool for the conversation genetic studies and population evaluation of P. sinensis.     

Development of 13 microsatellite markers for <Emphasis Type="Italic">Castanopsis tribuloides</Emphasis> (Fagaceae) using next-generation sequencing

Catanopsis tribuloides is a climax tree species commonly distributed in evergreen forests and has been used to restore degraded areas in northern Thailand. To aid in study of genetic diversity of the species, microsatellite markers, which are specific to C. tribuloides, were developed using whole genome sequencing by next-generation sequencing technology. The primers for microsatellite were developed and screened for optimal annealing temperature by PCR assay. The loci primers specific with C. tribuloides, 13 polymorphic microsatellite primers were successfully developed. The results from genetic information analyzing showed the number of alleles presented were between 2 and 24. Accordingly, the expected and observed heterozygosity obtained were between 0.298 and 0.920 and 0.364 to 1.000, respectively. Null allele frequency was presented 0.000–0.199. Genetic information was generated 10 loci primers significantly deviated from Hardy–Weinberg Equilibrium. All 13 primer pairs of loci were not significant with linkage disequilibrium. A set of microsatellite markers in this study could be applied to gene flow, genetic structure and population genetic studies in the future.     

Isolation and characterization of 20 polynucleotide microsatellite markers in a vulnerable Korean snail, <Emphasis Type="Italic">Ellobium chinense,</Emphasis> using 454 pyrosequencing

A small and air-breathing snail, Ellobium chinense (Ellobiidae), is a vulnerable species by International Union for Conservation of Nature (IUCN). To protect and manage habitat and population of E. chinense, microsatellite markers were developed using 454 pyrosequencing and 20 polymorphic microsatellite markers were identified. A total of 146,704 sequences containing a minimum of four repeat motifs (mean, 631 base pairs) were identified from 499,505 reads. Among 80 loci containing more than nine repeat units, 34 primer sets (42.5 %) produced strong PCR products, of which 20 were polymorphic among 48 samples of E. chinense. All loci exhibited high genetic variability, with an average of 18.9 alleles per locus, and the mean observed and expected heterozygosities were 0.65 and 0.90, respectively. In addition, cross-amplification was tested for all 20 loci in the same family species, Melampus sincaporensis. 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 20 newly developed microsatellites will be useful in future conservation genetic studies of this species.     

Characterization of 11 novel polymorphic microsatellite loci in the threatened Korean loach,Iksookimia koreensis,isolated using a next-generation sequencing method

To investigate the genetic diversity and population structure of Iksookimia koreensis, we characterized 11 novel polymorphic microsatellite markers developed using next-generation sequencing. The number of alleles per locus ranged from 4 to 10 (average = 6.26). Observed and expected heterozygosities ranged from 0.333 to 0.866 and from 0.375 to 0.866, respectively. No loci showed significant deviation from Hardy–Weinberg equilibrium (HWE). These loci were also used successfully to study the genetic diversity of a closely related species, Iksookimia longicorpa. Four of the 11 loci amplified in the two species showed different allele frequency and distribution, indicating deep genetic divergence between I. koreensis and I. longicorpa. The newly developed microsatellite markers reported here will provide a useful tool for examining gene flow, population genetic structure, and genetic diversity of these species.     

Isolation and characterization of 19 novel microsatellite loci in Rhododendron aureum and Rhododendron brachycarpum (Ericaceae)

Nineteen polymorphic microsatellite loci are described for Rhododendron aureum Georgi (Ericaceae), an endangered species in Korea, and the closely related Rhododendron brachycarpum. The sequences containing repeat motifs were identified using low-depth next generation sequencing and 148 microsatellite loci were examined for amplification success and the detection of polymorphisms. All 19 loci were polymorphic and the number of alleles for these markers varied from 4 to 25, with an average of 15 alleles per locus. Three R. aureum loci showed departure from Hardy–Weinberg equilibrium (HWE) and one R. brachycarpum locus deviated significantly from HWE. These microsatellite loci will be useful for investigating the genetic diversity and population structure of both species.     

Complete mitochondrial genome of the darkling beetle Gonocephalum outreyi (Coleoptera: Tenebrionidae) with phylogenetic implications

The complete mitochondrial genome (mitogenome) of Gonocephalum outreyi was determined by using next-generation sequencing approach. The full length of this mitogenome is 15,836?bp, which consists of 37 typical metazoan mitochondrial genes with an identical genome organization to ancestral insects. The majority of the protein-coding genes begin with the codon ATN, except for cox1 and cox2 with AAT and AAA, respectively. To elucidate the phylogenetic position of G. outreyi, we used various sequence coding schemes for protein-coding genes and the combined nucleotide sequences of all mitochondrial genes for tree building under the Bayesian and Maximum Likelihood inferences. The phylogenetic results consistently supported G. outreyi as a member of the family Tenebrionidae. The monophyly of both Tenebrionoidea and Tenebrionidae were strongly supported. The Scraptiidae and Melandryidae were recovered to be non-monophyletic in regards to the Osphya. Within Tenebrionidae, the subfamilies Diaperinae and Tenebrioninae were found to be non-monophyletic.     

A set of microsatellite markers for population genetics of leopard cat (Prionailurus bengalensis) and cross-species amplification in other felids

We describe the isolation and characterization of novel microsatellite loci from the leopard cat, Prionailurus bengalensis Kerr, 1792 (Family Felidae). Using Illumina HiSeq2500 sequencing technology, we sequenced the leopard cat genome and identified 1.5 million loci of simple sequence repeats with di- to deca-nucleotide motifs. We developed twelve polymorphic markers with tetra-nucleotide motif types after screening 35 loci for amplification and polymorphism. The observed and expected heterozygosities of the markers were 0.438 and 0.423, respectively. The number of alleles per locus ranged from 2 to 7, with a mean polymorphism information content of 0.383. Eleven loci were at Hardy-Weinberg equilibrium and no linkage disequilibrium was detected among any pairs of loci. We tested cross-species amplification of these markers across five other felids (Panthera tigris, P. pardus, P. onca, Acinonyx jubatus, and Felis catus). All loci were transferable to at least one other feline species and four amplified all five species. The microsatellite markers developed in this study will be valuable for estimating ecological parameters of populations and to establish conservation and management strategies for feline species.     

Development of 21 polymorphic microsatellite markers for the black-banded sea krait, <Emphasis Type="Italic">Laticauda semifasciata</Emphasis> (Elapidae: Laticaudinae), and cross-species amplification for two other congeneric species

The genus Laticauda (Reptilia: Elapidae), commonly known as sea kraits, is venomous marine amphibious snakes distributed throughout the south and southeast Asian islands and mostly found in coastal waters. To facilitate genetic studies, we have developed microsatellite loci for L. semifasciata using the 454 GS-FLX pyrosequencing technique. A total of 65,680 sequences containing a minimum of five repeat motifs were identified from 451,659 reads. Among 80 loci containing more than nine repeat units, 34 primer sets (42.5%) produced strong PCR products, of which 21 were polymorphic among 36 samples of L. semifasciata. All loci exhibited high genetic variability, with an average of 7.38 alleles per locus, and the mean observed and expected heterozygosities were 0.73 and 0.76, respectively. The cross-species amplification of these loci in two laticaudine species, L. colubrina and L. laticaudata, revealed a high transferability (78.6%) and polymorphism (59.5%) of the loci. Our work demonstrated the utility of next-generation 454 sequencing as the rapid and cost-effective method for development of microsatellite markers. The high level of polymorphism in these microsatellite loci will be useful for the detection of population subdivision and the study of migration, gene flow, relatedness and philopatry of L. semifasciata and other laticaudine species.     

Complete Taiwanese Macaque (Macaca cyclopis) Mitochondrial Genome: Reference-Assisted de novo Assembly with Multiple k-mer Strategy

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.     

Development of microsatellite markers for <Emphasis Type="Italic">Myracrodruon urundeuva</Emphasis> (F.F. &amp; M.F. Allemão), a highly endangered species from tropical forest based on next-generation sequencing

Myracrodruon urundeuva is a tree species of high economic importance due the strength and durability of its wood. Threatened of extinction in Brazil, it is present only in a few forest remnants, mostly in conservation units. Currently, there is little information on the genetic diversity of natural populations in Brazil and even less information about the genome of this species. Here, new species-specific microsatellite loci were developed based on next-generation sequencing (Illumina). More than 100,000 loci were identified in the run, with di- to hexanucleotides motifs. Of these, 20 loci were selected for validation in 30 individuals, with 15 successfully polymorphic loci detected. The number of alleles ranged among loci from 3 to 16, with an average of 7.73, expected (H _e ) and observed (H _o ) heterozygosity ranged from 0.246 to 0.902 and from 0.103 to 0.867, respectively. These results point out that these new set of markers has a great potential for use in population genetic studies for genetic conservation of the species.     

The development of 10 novel polymorphic microsatellite markers through next generation sequencing and a preliminary population genetic analysis for the endangered Glenelg spiny crayfish, Euastacus bispinosus

The Glenelg spiny crayfish, Euastacus bispinosus, is an iconic freshwater invertebrate of south eastern Australia and listed as ‘endangered’ under the Environment Protection and Biodiversity Conservation Act 1999, and ‘vulnerable’ under the International Union for Conservation of Nature’s Red List. The species has suffered major population declines as a result of over-fishing, low environmental flows, the introduction of invasive fish species and habitat degradation. In order to develop an effective conservation strategy, patterns of gene flow, genetic structure and genetic diversity across the species distribution need to be clearly understood. In this study we develop a suite of polymorphic microsatellite markers by next generation sequencing. A total of 15 polymorphic loci were identified and 10 characterized using 22 individuals from the lower Glenelg River. We observed low to moderate genetic variation across most loci (mean number of alleles per locus = 2.80; mean expected heterozygosity = 0.36) with no evidence of individual loci deviating significantly from Hardy–Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. Individuals from two additional sites (Crawford River, Victoria; Ewens Ponds Conservation Park, South Australia) were genotyped at all 10 loci and a preliminary investigation of genetic diversity and population structure was undertaken. Analyses indicate high levels of genetic differentiation among sample locations (F _ST  = 0.49), while the Ewens Ponds population is genetically homogeneous, indicating a likely small founder group and ongoing inbreeding. Management actions will be needed to restore genetic diversity in this and possibly other at risk populations. These markers will provide a valuable resource for future population genetic assessments so that an effective framework can be developed for implementing conservation strategies for E. bispinosus.     

Rapid microsatellite marker development in the endangered neotropical freshwater turtle Podocnemis lewyana (Testudines: Podocnemididae) using 454 sequencing

Using high throughput sequencing we obtained a large number of microsatellites from Podocnemis lewyana, an endemic turtle from northwestern South America. We used 454 Genome Sequence FLX platform of sheared genomic DNA from randomly sampling approximately 17% of the haploid genome. We identified 86,501 reads (8.1% of all reads) that contained our definition of microsatellite loci. AC and TC were the most abundant motifs in the P. lewyana genome. TGC and AAAC were most abundant tri and tetra-nucleotide motifs respectively. 72.7% of microsatellite reads had flanking sequence regions suitable for primer design and PCR amplification. We validated the identified potentially amplifiable loci (PAL) and tested for polymorphism by selecting 15 loci corresponding to tetranucleotides. Twelve loci showed polymorphism in eight individuals. These findings demonstrates that microsatellite detection using next-generation sequencing is an efficient way of getting a lot of loci for listed taxa and in turn will have a large impact on future genetic studies aiming to understand and implement conservation plans for this highly threatened freshwater turtle.     

Isolation and characterization of microsatellite marker loci in the Wagner’s mustached bat <Emphasis Type="BoldItalic">Pteronotus psilotis</Emphasis> (Chiroptera: Mormoopidae) and cross-amplification in other related species

Pteronotus psilotis, a mormoopid bat, is an insectivorous, gregarious and strict cave-dwelling species that is found areas between the sea level and an elevation of about 1000 masl. This species is present in diverse habitats ranging from rain forest to dry deciduous forest. Nine microsatellite loci were developed for Wagner’s mustached bat, Pteronotus psilotis using the next-generation sequencing approach, and their utility for population genetics studies was assessed. All loci were polymorphic (7–15 alleles) and characterized in 30 individuals from three P. psilotis populations, with the levels of observed and expected heterozygosity ranging from 0.280 to 0.867 and 0.584 to 0.842, respectively. One locus showed significant departures from Hardy–Weinberg expectations after Bonferroni correction. Cross-amplification in 11 other bat species was tested, for which eight microsatellites were successfully amplified, and of these seven were polymorphic. The development of these new microsatellite loci will contribute to investigations of genetic population structure, genetic diversity and gene flow in P. psilotis populations, as well as in other closely related bat species.     

Development of Genomic Microsatellite Markers in Carthamus tinctorius L. (Safflower) Using Next Generation Sequencing and Assessment of Their Cross-Species Transferability and Utility for Diversity Analysis

Background

Safflower (Carthamus tinctorius L.), an Asteraceae member, yields high quality edible oil rich in unsaturated fatty acids and is resilient to dry conditions. The crop holds tremendous potential for improvement through concerted molecular breeding programs due to the availability of significant genetic and phenotypic diversity. Genomic resources that could facilitate such breeding programs remain largely underdeveloped in the crop. The present study was initiated to develop a large set of novel microsatellite markers for safflower using next generation sequencing.

Principal Findings

Low throughput genome sequencing of safflower was performed using Illumina paired end technology providing ~3.5X coverage of the genome. Analysis of sequencing data allowed identification of 23,067 regions harboring perfect microsatellite loci. The safflower genome was found to be rich in dinucleotide repeats followed by tri-, tetra-, penta- and hexa-nucleotides. Primer pairs were designed for 5,716 novel microsatellite sequences with repeat length ≥ 20 bases and optimal flanking regions. A subset of 325 microsatellite loci was tested for amplification, of which 294 loci produced robust amplification. The validated primers were used for assessment of 23 safflower accessions belonging to diverse agro-climatic zones of the world leading to identification of 93 polymorphic primers (31.6%). The numbers of observed alleles at each locus ranged from two to four and mean polymorphism information content was found to be 0.3075. The polymorphic primers were tested for cross-species transferability on nine wild relatives of cultivated safflower. All primers except one showed amplification in at least two wild species while 25 primers amplified across all the nine species. The UPGMA dendrogram clustered C. tinctorius accessions and wild species separately into two major groups. The proposed progenitor species of safflower, C. oxyacantha and C. palaestinus were genetically closer to cultivated safflower and formed a distinct cluster. The cluster analysis also distinguished diploid and tetraploid wild species of safflower.

Conclusion

Next generation sequencing of safflower genome generated a large set of microsatellite markers. The novel markers developed in this study will add to the existing repertoire of markers and can be used for diversity analysis, synteny studies, construction of linkage maps and marker-assisted selection.     

The mitogenome of Onchocerca volvulus from the Brazilian Amazonia focus

We report here the first complete mitochondria genome of Onchocerca volvulus from a focus outside of Africa. An O. volvulus mitogenome from the Brazilian Amazonia focus was obtained using a combination of high-throughput and Sanger sequencing technologies. Comparisons made between this mitochondrial genome and publicly available mitochondrial sequences identified 46 variant nucleotide positions and suggested that our Brazilian mitogenome is more closely related to Cameroon-origin mitochondria than West African-origin mitochondria. As well as providing insights into the origins of Latin American onchocerciasis, the Brazilian Amazonia focus mitogenome may also have value as an epidemiological resource.     

Characterization of 13 polymorphic microsatellite loci in the Japanese land leech

We developed 13 polymorphic microsatellite loci of the Japanese land leech (Haemadipsa japonica; Haemadipsidea) using an Illumina MiSeq sequencing approach. A total of 42,064 nuclear DNA contigs were filtered for microsatellite motifs, among which 30,873 simple sequence repeat loci were identified. From these sequences, we selected 30 primer sets, and 13 of these loci were successfully amplified. Polymorphism of the 13 loci was tested using 16 individuals sampled from sixteen populations across Japan. The number of alleles and polymorphism information content varied from 5 to 17 and 0.335 to 0.883, respectively, and observed and expected heterozygosity values ranged from 0.143 to 0.875 and 0.349 to 0.893, respectively, indicating that these loci are polymorphic. Furthermore, we established useful multiplex PCR using these loci. The 13 microsatellite loci described in this paper are the first nuclear microsatellite markers for a land leech species.     

Isolation and characterization of microsatellite loci from Macadamia

Thirty‐three microsatellite loci were isolated for the Australian rainforest tree Macadamia integrifolia. Genotyping across a test panel of 43 commercial cultivars generated an average polymorphic information content of 0.480. Five loci showed no polymorphism across cultivars. Significant linkage disequilibrium was detected in 10 pairwise comparisons, including two pairs of loci identified from the same clone sequence. The 33 microsatellite loci represent a significant tool for genome mapping and population genetic studies.