Development and characterization of microsatellite loci in the Centricnemus leucogrammus weevil

Centricnemus leucogrammus is a weevil characteristic of European xerothermic habitats and steppes. The species was probably more widespread during the Pleistocene glaciations, while its current distribution is limited to “warm-stage refugia.” It may be regarded as a typical representative of flightless xerothermophilous beetles. Previous studies concentrated on its genetic variation using mitochondrial genes. Here, we identified, tested and characterized 24 polymorphic microsatellite loci with the use of 454 sequencing of microsatellite enriched genomic libraries. The new set of loci will be used in studies on the population structure of this weevil and may provide valuable information for its conservation.     

Genomic Resources Notes accepted 1 October 2013 – 30 November 2013

This article documents the public availability of (i) genomic sequence data and 43 microsatellite loci for the bat species, Lasiurus borealis and Lasiurus cinereus, and (ii) complete mitochondrial and partial nuclear genomes for two jack species, Caranx ignobilis, Caranx melampygus.     

Targeted multiplex next‐generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.     

Extracting ‘legacy loci’ from an invertebrate sequence capture data set

Sequence capture studies result in rich data sets comprising hundreds to thousands of targeted genomic regions that are superseding Sanger-based data sets comprised of a few well-known loci with historical uses in phylogenetics (‘legacy loci’). However, integrating sequence capture and Sanger-based data sets is of interest as legacy loci can include different types of loci (e.g. mitochondrial and nuclear) across a potentially larger sample of species from past studies. Sequence capture data sets include nontargeted sequences, and there has been recent interest in extracting legacy loci from invertebrate data sets. Here, we use published legacy data from leaf-footed bugs (Hemiptera: Coreoidea) to recover 15 mitochondrial and seven nuclear legacy loci from off-target sequences in a sequence capture data set, explore approaches to improve legacy locus recovery, and combine these loci with sequence capture data for phylogenetic analysis. Two nuclear loci were determined to already be targeted by sequence capture baits. Most of the remaining loci were successfully recovered from off-target sequences, but this recovery varied greatly. Additionally, complementing complete mitogenomes with additional reference mitochondrial sequences from a genetic depository did not offer improvement for most of our taxa; however, supplementing these reference sequences with extracted legacy loci offered ≥6% improvement across taxa for a given mitochondrial locus (negligible improvement for nuclear loci). Phylogenetic analysis of legacy and sequence capture data produced a topology generally congruent with recent studies, but support was lower. Thus, future studies may employ the approaches used in this study to integrate legacy data with newly generated sequence capture data sets without added expenses.     

Somaclonal variation of haploid in vitro tissue culture obtained from Siberian larch (Larix sibirica Ledeb.) megagametophytes for whole genome de novo sequencing

The objective of this study was to obtain a genetically stable haploid in vitro-derived line from Siberian larch (Larix sibirica Ledeb.) using megagametophyte explants, which then could be used for different molecular genetic studies, including whole genome de novo sequencing. However, cytogenetic analysis and genotyping of 11 microsatellite loci showed high levels of genomic instability and a high frequency of mutation in the obtained megagametophyte-derived callus cultures. All cultures contained new mutations in one or more microsatellite loci.     

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.     

Characterisation of the complete mitochondrial genome and 13 microsatellite loci through next-generation sequencing for the New Caledonian spider-ant Leptomyrmex pallens

The complete mitochondrial genome and a set of polymorphic microsatellite markers were identified by 454 pyrosequencing (1/16th of a plate) for the New Caledonian rainforest spider-ant Leptomyrmex pallens. De novo genome assembly recovered the entire mitochondrial genome with mean coverage of 8.9-fold (range 1–27). The mitogenome consists of 15,591 base pairs including 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The genome arrangement is typical of insect taxa and very similar to the only other published ant mitogenome from the Solenopsis genus, with the main differences consisting of translocations and inversions of tRNAs. A total of 13 polymorphic loci were also characterized using 41 individuals from a single population in the Aoupinié region, corresponding to workers from 21 nests and 16 foraging workers. We observed moderate genetic variation across most loci (mean number of alleles per locus = 4.50; mean expected heterozygosity = 0.53) with evidence of only two loci deviating significantly from Hardy–Weinberg equilibrium due to null alleles. Marker independence was confirmed with tests for linkage disequilibrium. Most loci cross amplified for three additional Leptomyrmex species. The annotation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the colony structure, population genetic patterns, and dispersal strategy of L. pallens in the context of rainforest fragmentation in New Caledonia. Furthermore, this paper confirms a recent line of evidence that comprehensive mitochondrial data can be obtained relatively easily from small next-generation sequencing analyses. Greater synthesis of next-generation sequencing data will play a significant role in expanding the taxonomic representation of mitochondrial genome sequences.     

Development and characterization of eleven polymorphic microsatellite loci from South China field mouse (Apodemus draco)

One population distributed in Yunnan of China was regarded as a new species based on mitochondrial cytochrome b sequences. However, the usefulness of mitochondrial sequence data in determining species boundaries is not universally agreed upon, the frequency data from multiple nuclear gene loci is necessary in determining species boundaries. So, we describe in this paper the isolation and characterization of eleven microsatellite loci in the South China field mouse from genomic DNA-enriched libraries. The eleven loci were tested in 24 individuals from two populations in Southwest China. These loci were highly polymorphic with numbers of alleles per locus ranging from 9 to 24 and expected heterozygosities from 0.898 to 0.967. Eight loci followed Hardy–Weinberg expectations after Bonferroni correction for multiple comparisons. No significant linkage association was found among all these loci. The eleven polymorphic microsatellite loci will be useful in determining species boundaries of the South China field mouse.     

Isolation,characterization and cross‐species amplification of eight microsatellite DNA loci in the migratory locust (Locusta migratoria)

Eight polymorphic di‐ and trinucleotide microsatellite loci suitable for population genetic analysis were developed in Locusta migratoria from a partial phagemid genomic library enriched for microsatellite inserts. The expected heterozygosity at these loci ranges from 0.45 to 0.97, with the observed allele numbers varying between nine and 45. The overall microsatellite cloning efficiency in L. migratoria is 14%, suggesting that in migratory locusts, microsatellite sequences are abundant and should provide a valuable and easily accessible source of nuclear markers for genetic studies. These microsatellite loci were highly Locusta‐specific, with only very limited cross‐species applicability.     

Comprehensive whole genome survey analyses of male and female brown-spotted flathead fish Platycephalus sp.1

The flathead fish Platycephalus sp.1 is an ecologically and commercially important marine fish in the northwestern Pacific with notable sexual differences in growth and development. Yet the genomic data of this species is lacking. In the present study, whole genome sequencing of two individuals (one male and one female) of Platycephalus sp.1 were conducted to provide fundamental genomic information. The genome sizes were estimated to be 674.96 Mb (male) and 684.15 Mb (female) by using k-mer analyses. The heterozygosity and repeat ratios suggested possible male heterogamety of this species. The draft genome sequences were initially assembled and genome-wide microsatellite motifs were identified. Besides, the complete mitochondrial genome sequences were assembled and the phylogenetic analyses genetically supported the validation of Platycephalus sp.1. The reported genomic data and genetic markers in this study could be useful in future comparative genomics and evolutionary biology studies.     

Comprehensive Draft Genome Analyses of Three Rockfishes (Scorpaeniformes,Sebastiscus) via Genome Survey Sequencing

Sebastiscus species, marine rockfishes, are of essential economic value. However, the genomic data of this genus is lacking and incomplete. Here, whole genome sequencing of all species of Sebastiscus was conducted to provide fundamental genomic information. The genome sizes were estimated to be 802.49 Mb (S. albofasciatus), 786.79 Mb (S. tertius), and 776.00 Mb (S. marmoratus) by using k-mer analyses. The draft genome sequences were initially assembled, and genome-wide microsatellite motifs were identified. The heterozygosity, repeat ratios, and numbers of microsatellite motifs all suggested possibly that S. tertius is more closely related to S. albofasciatus than S. marmoratus at the genetic level. Moreover, the complete mitochondrial genome sequences were assembled from the whole genome data and the phylogenetic analyses genetically supported the validation of Sebastiscus species. This study provides an important genome resource for further studies of Sebastiscus species.     

Genetic diversity of the Pacific walrus (<Emphasis Type="Italic">Odobenus rosmarus divergens</Emphasis>) in the western part of the Chukchi Sea

We study 117 Pacific walrus samples from three rookeries within the western part of Chukchi Sea (Cape Vankarem, Cape Serdtse-Kamen, and Kolyuchin Island). We analyze the variability of nuclear (20 microsatellite loci) and mitochondrial DNA (three fragments). Two microsatellite loci which are described as microsatellites for the first time are used in this study: repeated sequences within introns of Coro1c and Plod2 genes. A high degree of genetic diversity is demonstrated for both nuclear and mitochondrial markers compared to Atlantic walrus. A high degree of genetic diversity is preserved within populations of Pacific walrus, despite a strong decline in the recent past. We discover the absence of significant differentiation for microsatellite loci and the presence of weak differentiation for mtDNA (mainly for a D-loop fragment). Walrus specimens that use the rookeries of the western part of Chukchi Sea are thought to belong to a single reproductive group.     

Automated screening and primer design of fish microsatellite DNA loci on pyrosequencing data

In recent years, massive sequencing approaches have allowed us to determine genomic structures of various organisms rapidly, raising novel applicability of the high-throughput sequence data obtained to various fields of biological studies. We present here a pipeline to search for microsatellite DNA and design PCR primers encompassing the microsatellites on genomic sequence data produced by 454 pyrosequencing. We tested this pipeline, called ‘Auto-primer’, on several fish genomic sequences and obtained many and various candidates for microsatellite DNA loci useful for detecting intraspecies genetic variability. This in silico search for microsatellite DNA is superior to conventional cloning methods, since any sequence patterns of repeat unit can be screened.     

Contrasting levels of genetic differentiation among populations of wolverines (<Emphasis Type="Italic">Gulo gulo</Emphasis>) from northern Canada revealed by nuclear and mitochondrial loci

Habitat loss, fragmentation, overharvest, and other anthropogenic factors have resulted in population and distribution declines for North American wolverines (Gulo gulo). Currently, wolverines east of the Hudson Bay are endangered and possibly extinct, whereas the status of wolverines throughout the remaining Holarctic is vulnerable. Three previous studies using nuclear loci have detected little to no significant structuring among wolverines sampled across northern Canada. Based on these results it has been suggested that wolverines in northern Canada represent a single, panmictic population. However, as has been shown in numerous studies, in cases of female site fidelity, it is possible to have demographically autonomous populations even with male-biased gene flow. To better assess the genetic structure of wolverines in northern Canada, we examined nine microsatellite loci and DNA sequence variation from a 200 bp fragment of the mitochondrial (mtDNA) control region for 270 wolverines from nine collecting areas representing three regions of northern Canada. In agreement with previous studies, microsatellite analyses revealed a lack of significant population substructure (F _ST=0.0004). However, analysis of molecular variance, comparisons of pairwise F _ST values and nested-clade analysis of the mtDNA data revealed considerable genetic structuring among samples of wolverines from these three regions of northern Canada. These mitochondrial data provide evidence that wolverines in Canada are genetically structured due to female philopatry. The contrasting patterns of genetic differentiation based on nuclear and mitochondrial data highlight the importance of examining both nuclear and mitochondrial loci when attempting to elucidate patterns of genetic structure.     

Hybridisation,paternal leakage and mitochondrial DNA linearization in three anomalous fish (Scombridae)

Using mitochondrial DNA for species identification and population studies assumes that the genome is maternally inherited, circular, located in the cytoplasm and lacks recombination. This study explores the mitochondrial genomes of three anomalous mackerel. Complete mitochondrial genome sequencing plus nuclear microsatellite genotyping of these fish identified them as Scomberomorus munroi (spotted mackerel). Unlike normal S. munroi, these three fish also contained different linear, mitochondrial genomes of Scomberomorus semifasciatus (grey mackerel). The results are best explained by hybridisation, paternal leakage and mitochondrial DNA linearization. This unusual observation may provide an explanation for mtDNA outliers in animal population studies.     

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.     

Rapid identification of thousands of copperhead snake (Agkistrodon contortrix) microsatellite loci from modest amounts of 454 shotgun genome sequence

Optimal integration of next-generation sequencing into mainstream research requires re-evaluation of how problems can be reasonably overcome and what questions can be asked. One potential application is the rapid acquisition of genomic information to identify microsatellite loci for evolutionary, population genetic and chromosome linkage mapping research on non-model and not previously sequenced organisms. Here, we report on results using high-throughput sequencing to obtain a large number of microsatellite loci from the venomous snake Agkistrodon contortrix, the copperhead. We used the 454 Genome Sequencer FLX next-generation sequencing platform to sample randomly ∼27 Mbp (128 773 reads) of the copperhead genome, thus sampling about 2% of the genome of this species. We identified microsatellite loci in 11.3% of all reads obtained, with 14 612 microsatellite loci identified in total, 4564 of which had flanking sequences suitable for polymerase chain reaction primer design. The random sequencing-based approach to identify microsatellites was rapid, cost-effective and identified thousands of useful microsatellite loci in a previously unstudied species.     

Microsatellite Analysis of Genetic Variation and Population Subdivision for the Black Muntjac, Muntiacus crinifrons

The black muntjac (Muntiacus crinifrons) is a rare deer found only in a restricted region in east China. Recent studies of mitochondrial DNA diversity have shown a markedly low level of nucleotide diversity for the species, and the Suichang population was genetically differentiated from the two other populations, in Huangshan and Tianmushan mountains. In this study, we extended the analysis of genetic diversity and population subdivision for the black muntjac using data from 11 highly polymorphic nuclear DNA microsatellite loci. Contrary to the results based on mtDNA data, the microsatellite loci revealed that the black muntjac retained a rather high nuclear genetic diversity (overall average H (E) = 0.78). Nevertheless, both types of markers supported the idea that the extant black muntjac population is genetically disrupted (overall phi (ST) = 0.16 for mtDNA and overall F (ST) = 0.053 for microsatellite, both P < 0.001). The correlation between genetic differentiation and geographic distance was not significant (Mantel test; P > 0.05), implying that the patterns of genetic differentiation observed in this study might result from recent habitat fragmentation or loss. Based on the results from the mtDNA and nuclear DNA data sets, two management units were defined for the species, Huangshan/Tianmushan and Suichang. We also recommend that a new captive population be established with individuals from the Suichang region as a founder source.     

High-throughput identification of informative nuclear Loci for shallow-scale phylogenetics and phylogeography

One of the major challenges for researchers studying phylogeography and shallow-scale phylogenetics is the identification of highly variable and informative nuclear loci for the question of interest. Previous approaches to locus identification have generally required extensive testing of anonymous nuclear loci developed from genomic libraries of the target taxon, testing of loci of unknown utility from other systems, or identification of loci from the nearest model organism with genomic resources. Here, we present a fast and economical approach to generating thousands of variable, single-copy nuclear loci for any system using next-generation sequencing. We performed Illumina paired-end sequencing of three reduced-representation libraries (RRLs) in chorus frogs (Pseudacris) to identify orthologous, single-copy loci across libraries and to estimate sequence divergence at multiple taxonomic levels. We also conducted PCR testing of these loci across the genus Pseudacris and outgroups to determine whether loci developed for phylogeography can be extended to deeper phylogenetic levels. Prior to sequencing, we conducted in silico digestion of the most closely related reference genome (Xenopus tropicalis) to generate expectations for the number of loci and degree of coverage for a particular experimental design. Using the RRL approach, we: (i) identified more than 100,000 single-copy nuclear loci, 6339 of which were obtained for divergent conspecifics and 904 of which were obtained for heterospecifics; (ii) estimated average nuclear sequence divergence at 0.1% between alleles within an individual, 1.1% between conspecific individuals that represent two different clades, and 1.8% between species; and (iii) determined from PCR testing that 53% of the loci successfully amplify within-species and also many amplify to the genus-level and deeper in the phylogeny (16%). Our study effectively identified nuclear loci present in the genome that have levels of sequence divergence on par with mitochondrial loci commonly used in phylogeography. Specifically, we estimated that ～7% of loci in the chorus frog genome are >3% divergent within species; this translates to a prediction of approximately 50,000 single-copy loci in the genome with >3% divergence. Moreover, successful amplification of many loci at deeper phylogenetic levels indicates that the RRL approach represents an efficient method for rapid identification of informative loci for both phylogenetics and phylogeography. We conclude by making recommendations for minimizing the cost and maximizing the efficiency of locus identification for future studies in this field.