Isolation and characterization of Ligustrum micranthum (Oleaceae) microsatellite loci using paired‐end Illumina reads

Twenty‐six microsatellite loci were developed and characterized for Ligustrum micranthum, a species endemic to the Ogasawara Islands, Japan. The genetic structure of this species must be clarified in order to restore the island's ecosystem. A total of 8511 primer pairs were designed from de novo sequencing. Of the 48 primer pairs selected, amplification and polymorphisms were tested using one population each from the Chichijima and Hahajima Islands of the Ogasawara Islands. Twenty‐six microsatellite loci were successfully amplified and the number of alleles for these loci ranged from five to 31 per locus, and the mean expected heterozygosities were 0.858 and 0.849, respectively. No significant deviation from the Hardy–Weinberg equilibrium was observed in either population, and no significant linkage disequilibrium was detected between any locus pair. The microsatellite loci reported in this study can be used in future studies to evaluate the genetic structure and mating system of L. micranthum.     

Isolation and characterization of tri‐ and tetra‐repeat microsatellite loci in the white‐spotted charr Salvelinus leucomaenis (Salmonidae)

Tri‐ and tetra‐motif repeat microsatellite marker loci were developed for the white‐spotted charr Salvelinus leucomaenis. The 454 pyrosequencing was used to discover repeat arrays, and eight microsatellite‐primer sets, available for the estimation of polymorphisms, were identified. The number of alleles in a wild population ranged from two to four and the observed and expected heterozygosities were 0·180–0·600 and 0·188–0·599, respectively.     

Development and characterization of polymorphic microsatellite markers for Neolitsea sericea using Illumina paired‐end draft sequencing data

Simple sequence repeat (SSR) markers were developed and characterized for Neolitsea sericea (Bl.) Koidz. (Lauraceae). Out of 196 designed primer pairs, a total of 144 pairs showed amplification, of which 44 had clear and stable chromatograms. Polymorphism of these 44 loci was tested using 32 individuals sampled from a single population of N. sericea. The number of alleles and the polymorphism information content varied from 3 to 12 and 0.271 to 0.853, respectively. A significant departure from the Hardy‐Weinberg equilibrium was observed in one of the 44 loci. These SSR markers are useful for population genetic studies and parentage analysis in N. sericea, which is one of the most common evergreen species in coastal Pinus thunbergii forests in central‐western Japan.     

DNA barcodes from century‐old type specimens using next‐generation sequencing

Type specimens have high scientific importance because they provide the only certain connection between the application of a Linnean name and a physical specimen. Many other individuals may have been identified as a particular species, but their linkage to the taxon concept is inferential. Because type specimens are often more than a century old and have experienced conditions unfavourable for DNA preservation, success in sequence recovery has been uncertain. This study addresses this challenge by employing next‐generation sequencing (NGS) to recover sequences for the barcode region of the cytochrome c oxidase 1 gene from small amounts of template DNA. DNA quality was first screened in more than 1800 century‐old type specimens of Lepidoptera by attempting to recover 164‐bp and 94‐bp reads via Sanger sequencing. This analysis permitted the assignment of each specimen to one of three DNA quality categories – high (164‐bp sequence), medium (94‐bp sequence) or low (no sequence). Ten specimens from each category were subsequently analysed via a PCR‐based NGS protocol requiring very little template DNA. It recovered sequence information from all specimens with average read lengths ranging from 458 bp to 610 bp for the three DNA categories. By sequencing ten specimens in each NGS run, costs were similar to Sanger analysis. Future increases in the number of specimens processed in each run promise substantial reductions in cost, making it possible to anticipate a future where barcode sequences are available from most type specimens.     

Multiplex restriction amplicon sequencing: a novel next‐generation sequencing‐based marker platform for high‐throughput genotyping

To enable rapid selection of traits in marker‐assisted breeding, markers must be technically simple, low‐cost, high‐throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3’‐ends, preceded by 6‐10 bases of specific or degenerate nucleotide sequences and then by a unique M13‐tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next‐generation sequencing‐based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.     

amplisas: a web server for multilocus genotyping using next‐generation amplicon sequencing data

Next‐generation sequencing (NGS) technologies are revolutionizing the fields of biology and medicine as powerful tools for amplicon sequencing (AS). Using combinations of primers and barcodes, it is possible to sequence targeted genomic regions with deep coverage for hundreds, even thousands, of individuals in a single experiment. This is extremely valuable for the genotyping of gene families in which locus‐specific primers are often difficult to design, such as the major histocompatibility complex (MHC). The utility of AS is, however, limited by the high intrinsic sequencing error rates of NGS technologies and other sources of error such as polymerase amplification or chimera formation. Correcting these errors requires extensive bioinformatic post‐processing of NGS data. Amplicon Sequence Assignment (amplisas ) is a tool that performs analysis of AS results in a simple and efficient way, while offering customization options for advanced users. amplisas is designed as a three‐step pipeline consisting of (i) read demultiplexing, (ii) unique sequence clustering and (iii) erroneous sequence filtering. Allele sequences and frequencies are retrieved in excel spreadsheet format, making them easy to interpret. amplisas performance has been successfully benchmarked against previously published genotyped MHC data sets obtained with various NGS technologies.     

exonsampler: a computer program for genome‐wide and candidate gene exon sampling for targeted next‐generation sequencing

The computer program exonsampler automates the sampling of thousands of exon sequences from publicly available reference genome sequences and gene annotation databases. It was designed to provide exon sequences for the efficient, next‐generation gene sequencing method called exon capture. The exon sequences can be sampled by a list of gene name abbreviations (e.g. IFNG, TLR1), or by sampling exons from genes spaced evenly across chromosomes. It provides a list of genomic coordinates (a bed file), as well as a set of sequences in fasta format. User‐adjustable parameters for collecting exon sequences include a minimum and maximum acceptable exon length, maximum number of exonic base pairs (bp) to sample per gene, and maximum total bp for the entire collection. It allows for partial sampling of very large exons. It can preferentially sample upstream (5 prime) exons, downstream (3 prime) exons, both external exons, or all internal exons. It is written in the Python programming language using its free libraries. We describe the use of exonsampler to collect exon sequences from the domestic cow (Bos taurus) genome for the design of an exon‐capture microarray to sequence exons from related species, including the zebu cow and wild bison. We collected ~10% of the exome (~3 million bp), including 155 candidate genes, and ~16 000 exons evenly spaced genomewide. We prioritized the collection of 5 prime exons to facilitate discovery and genotyping of SNPs near upstream gene regulatory DNA sequences, which control gene expression and are often under natural selection.     

High‐throughput SNP discovery in the rabbit (Oryctolagus cuniculus) genome by next‐generation semiconductor‐based sequencing

The European rabbit (Oryctolagus cuniculus) is a domesticated species with one of the broadest ranges of economic and scientific applications and fields of investigation. Rabbit genome information and assembly are available (oryCun2.0), but so far few studies have investigated its variability, and massive discovery of polymorphisms has not been published yet for this species. Here, we sequenced two reduced representation libraries (RRLs) to identify single nucleotide polymorphisms (SNPs) in the rabbit genome. Genomic DNA of 10 rabbits belonging to different breeds was pooled and digested with two restriction enzymes (HaeIII and RsaI) to create two RRLs which were sequenced using the Ion Torrent Personal Genome Machine. The two RRLs produced 2 917 879 and 4 046 871 reads, for a total of 280.51 Mb (248.49 Mb with quality >20) and 417.28 Mb (360.89 Mb with quality >20) respectively of sequenced DNA. About 90% and 91% respectively of the obtained reads were mapped on the rabbit genome, covering a total of 15.82% of the oryCun2.0 genome version. The mapping and ad hoc filtering procedures allowed to reliably call 62 491 SNPs. SNPs in a few genomic regions were validated by Sanger sequencing. The Variant Effect Predictor Web tool was used to map SNPs on the current version of the rabbit genome. The obtained results will be useful for many applied and basic research programs for this species and will contribute to the development of cost‐effective solutions for high‐throughput SNP genotyping in the rabbit.     

A targeted next‐generation sequencing toolkit for exon‐based cichlid phylogenomics

Cichlid fishes (family Cichlidae) are models for evolutionary and ecological research. Massively parallel sequencing approaches have been successfully applied to study relatively recent diversification in groups of African and Neotropical cichlids, but such technologies have yet to be used for addressing larger‐scale phylogenetic questions of cichlid evolution. Here, we describe a process for identifying putative single‐copy exons from five African cichlid genomes and sequence the targeted exons for a range of divergent (>tens of millions of years) taxa with probes designed from a single reference species (Oreochromis niloticus, Nile tilapia). Targeted sequencing of 923 exons across 10 cichlid species that represent the family's major lineages and geographic distribution resulted in a complete taxon matrix of 564 exons (649 549 bp), representing 559 genes. Maximum likelihood and Bayesian analyses in both species tree and concatenation frameworks yielded the same fully resolved and highly supported topology, which matched the expected backbone phylogeny of the major cichlid lineages. This work adds to the body of evidence that it is possible to use a relatively divergent reference genome for exon target design and successful capture across a broad phylogenetic range of species. Furthermore, our results show that the use of a third‐party laboratory coupled with accessible bioinformatics tools makes such phylogenomics projects feasible for research groups that lack direct access to genomic facilities. We expect that these resources will be used in further cichlid evolution studies and hope the protocols and identified targets will also be useful for phylogenetic studies of a wider range of organisms.     

Characterization of newly developed microsatellite loci in the stemmed yellow violet,Viola pubescens (Violaceae)

Viola pubescens is a perennial herb with a mixed mating system involving the production of both open, chasmogamous flowers and closed, cleistogamous flowers. Using a nonradioactive protocol, 12 microsatellite primers that amplified loci with perfect and imperfect repeats were developed in the species. As examined in five populations of V. pubescens using multiplex reactions, all loci exhibited polymorphism with an average of 5.4 alleles. Most primers also amplified in 12 additional Viola species and one cultivar, revealing the possibility of wide applicability across the genus. This is the first known report of microsatellite loci developed in Viola and the Violaceae.     

Efficient and sensitive identification and quantification of airborne pollen using next‐generation DNA sequencing

Pollen monitoring is an important and widely used tool in allergy research and creation of awareness in pollen‐allergic patients. Current pollen monitoring methods are microscope‐based, labour intensive and cannot identify pollen to the genus level in some relevant allergenic plant groups. Therefore, a more efficient, cost‐effective and sensitive method is needed. Here, we present a method for identification and quantification of airborne pollen using DNA sequencing. Pollen is collected from ambient air using standard techniques. DNA is extracted from the collected pollen, and a fragment of the chloroplast gene trnL is amplified using PCR. The PCR product is subsequently sequenced on a next‐generation sequencing platform (Ion Torrent). Amplicon molecules are sequenced individually, allowing identification of different sequences from a mixed sample. We show that this method provides an accurate qualitative and quantitative view of the species composition of samples of airborne pollen grains. We also show that it correctly identifies the individual grass genera present in a mixed sample of grass pollen, which cannot be achieved using microscopic pollen identification. We conclude that our method is more efficient and sensitive than current pollen monitoring techniques and therefore has the potential to increase the throughput of pollen monitoring.     

Development and characterization of 24 microsatellite markers in Primula tosaensis,an endangered primrose,using MiSeq

Primula tosaensis (Primulaceae) is an endangered primrose endemic to Japan. In this study, 24 novel microsatellite markers were developed using Illumina MiSeq sequencing to facilitate conservation of this endangered species. The genetic diversity and polymorphisms of these novel markers were measured in 32 individuals from a wild P. tosaensis population. The number of alleles and expected heterozygosities ranged from 2 to 5 (mean = 2.8) and from 0.119 to 0.724 (mean = 0.395), respectively. All loci were in Hardy–Weinberg equilibrium. The markers developed in this study will provide a powerful and practical tool for investigating the population structure and genetic diversity of P. tosaensis.     

Discovering and sequencing new plant viral genomes by next‐generation sequencing: description of a practical pipeline

Small‐scale sequencing has improved substantially in recent decades, culminating in the development of next‐generation sequencing (NGS) technologies. Modern NGS methods have helped the discovery of many new plant viruses. Nevertheless, there is still a need to establish solid assembly pipelines targeting small genomes characterised by low identities to known viral sequences. Here, we describe and discuss the fundamental steps required for discovering and sequencing new plant viral genomes by NGS. A practical pipeline and standard alternative tools used in NGS analysis are presented.     

Unlocking the secondary gene‐pool of barley with next‐generation sequencing

Crop wild relatives (CWR) provide an important source of allelic diversity for any given crop plant species for counteracting the erosion of genetic diversity caused by domestication and elite breeding bottlenecks. Hordeum bulbosum L. is representing the secondary gene pool of the genus Hordeum. It has been used as a source of genetic introgressions for improving elite barley germplasm (Hordeum vulgare L.). However, genetic introgressions from H. bulbosum have yet not been broadly applied, due to a lack of suitable molecular tools for locating, characterizing, and decreasing by recombination and marker‐assisted backcrossing the size of introgressed segments. We applied next‐generation sequencing (NGS) based strategies for unlocking genetic diversity of three diploid introgression lines of cultivated barley containing chromosomal segments of its close relative H. bulbosum. Firstly, exome capture‐based (re)‐sequencing revealed large numbers of single nucleotide polymorphisms (SNPs) enabling the precise allocation of H. bulbosum introgressions. This SNP resource was further exploited by designing a custom multiplex SNP genotyping assay. Secondly, two‐enzyme‐based genotyping‐by‐sequencing (GBS) was employed to allocate the introgressed H. bulbosum segments and to genotype a mapping population. Both methods provided fast and reliable detection and mapping of the introgressed segments and enabled the identification of recombinant plants. Thus, the utilization of H. bulbosum as a resource of natural genetic diversity in barley crop improvement will be greatly facilitated by these tools in the future.     

Isolation and characterization of nine microsatellite markers from Coffea arabica L., showing wide cross‐species amplifications

Genetic improvement of coffee (Coffea arabica L.) is constrained by low genetic diversity and lack of genetic markers, suitable screening tools, information on the genetic make‐up of available gene pool and long generation time. In this context, use of DNA markers such as microsatellites that provide high genetic‐resolution becomes highly desirable. Here, we report the development of nine new microsatellite markers from partial genomic library of an elite variety of Coffea arabica. The developed microsatellites revealed robust cross‐species amplifications in 17 related species of coffee, and their Polymorphic Information Content varied from 0 to 0.6, 0 to 0.78 and 0.67 to 0.90 for the arabica, robusta genotypes and species representatives, respectively. The data thus suggest their potential use as genetic markers for assessment of germplasm diversity and linkage analysis of coffee.     

Microsatellite markers developed using a next‐generation sequencing technique for Neotrogla spp. (Psocodea: Prionoglarididae), cave dwelling insects with sex‐reversed genitalia

The genus Neotrogla (Psocodea: Prinoglarididae) comprises four named species from Brazil. Females of this cave‐dwelling insect are characterized by a conspicuous penis‐like intromittent organ, termed a gynosome, which is inserted into the vagina‐like male genitalia during copulation. Another evolutionarily novel structure, the spermathecal plate, enables a female to simultaneously store two freshly deposited spermatophores (consisting of sperm and possibly nutritious substances) in her sperm storage organ (spermatheca). It is unknown whether the two spermatophores are derived from two different males. To investigate the mating ecology and population genetic structures of these insects with sex‐reversed genitalia, 16 novel highly polymorphic microsatellite loci were isolated and characterized based on ~2,275 Mbp genomic sequences from an undescribed Neotrogla species. Our first screening detected 99,888 candidate loci. Similar to other hemipteroid insects studied thus far, AAT motif microsatellites were conspicuously dominant. We further screened 99 sequences, for which 50 pairs of polymerase chain reaction primers were successfully designed. Sixteen of these primers successfully amplified products of the expected size in the 11 Neotrogla sp. individuals collected from two caves. The number of alleles per loci varied from two to nine, with no significant deviation from the Hardy–Weinberg equilibrium in either population. Although the caves sampled were only approximately 1 km apart, significant genetic differentiation was detected between the two populations. In total, 13, 12, 13 and 11 loci were cross‐amplified in N. aurora, N. brasiliensis, N. curvata and N. truncata, respectively, indicating the applicability of these microsatellite loci for metapopulation genetic studies in multiple Neotrogla species.     

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.