Current trends in microsatellite genotyping

PIDS: A User-Friendly Plant DNA Fingerprint Database Management System

The high variability and somatic stability of DNA fingerprints can be used to identify individuals, which is of great value in plant breeding. DNA fingerprint databases are essential and important tools for plant molecular research because they provide powerful technical and information support for crop breeding, variety quality control, variety right protection, and molecular marker-assisted breeding. Building a DNA fingerprint database involves the production of large amounts of heterogeneous data for which storage, analysis, and retrieval are time and resource consuming. To process the large amounts of data generated by laboratories and conduct quality control, a database management system is urgently needed to track samples and analyze data. We developed the plant international DNA-fingerprinting system (PIDS) using an open source web server and free software that has automatic collection, storage, and efficient management functions based on merging and comparison algorithms to handle massive microsatellite DNA fingerprint data. PIDS also can perform genetic analyses. This system can match a corresponding capillary electrophoresis image on each primer locus as fingerprint data to upload to the server. PIDS provides free customization and extension of back-end functions to meet the requirements of different laboratories. This system can be a significant tool for plant breeders and can be applied in forensic science for human fingerprint identification, as well as in virus and microorganism research.

Genomic Encryption of Digital Data Stored in Synthetic DNA

Today, we can read human genomes and store digital data robustly in synthetic DNA. Herein, we report a strategy to intertwine these two technologies to enable the secure storage of valuable information in synthetic DNA, protected with personalized keys. We show that genetic short tandem repeats (STRs) contain sufficient entropy to generate strong encryption keys, and that only one technology, DNA sequencing, is required to simultaneously read the key and the data. Using this approach, we experimentally generated 80 bit strong keys from human DNA, and used such a key to encrypt 17 kB of digital information stored in synthetic DNA. Finally, the decrypted information was recovered perfectly from a single massively parallel sequencing run.

A 13-plex of tetra- and penta-STRs to identify zebrafish

The zebrafish species Danio rerio has become one of the major vertebrate model organisms used in biomedical research. However, there are aspects of the model that need to be improved. One of these is the ability to identify individual fish and fish lines by DNA profiling. Although many dinucleotide short tandem repeat (diSTR) markers are available for this and similar purposes, they have certain disadvantages such as an excessive polymerase slippage ("stutter") that causes difficulties in automated genotyping and cross-laboratory comparisons. Here we report on the development of a 13-plex of tetranucleotide and pentanucleotide STRs (tetraSTRs and pentaSTRs, respectively) that have low stutter. The system uses an inexpensive universal primer labelling system, which can easily be converted to a direct labeling system if desired. This 13-plex was examined in three zebrafish lines (NHGRI-1, kca33Tg, and kca66Tg, originally obtained from ZIRC). The average observed heterozygosity (Ho) and expected heterozygosity (He) in these highly inbred lines were 0.291 and 0.359, respectively, which is very similar to what has been found with diSTRs. The probability of identity (PI) for all fish tested was 2.1 × 10^-5 and the PI for siblings (PIsib) was 6.4 × 10^-3, as calculated by the Genalex package. Ninety percent of the fish tested were correctly identified with their respective strains. It is also demonstrated that this panel can be used to confirm doubled-haploid cell lines. This multiplex should find multiple uses for improving the accuracy and reproducibility of studies using the zebrafish model.

Examining multiple paternity in the raccoon dog (Nyctereutes procyonoides) in Japan using microsatellite analysis

We analyzed the genotypes of three pregnant females and their litters to investigate the phenomenon of multiple paternity in wild raccoon dogs (Nyctereutes procyonoides) using 17 microsatellite markers. If a female has mated with only one male during estrus, then the maximum number of paternal alleles will not exceed two among littermates with the same father. The results revealed two out of three litters had three or four paternal alleles at one or five microsatellite loci. Therefore, the female had mated with more than one male during estrus. To the best of our knowledge, the present study is the first to report the possibility of multiple paternity in wild raccoon dogs.

Development of genome-wide polymorphic microsatellite markers for Trichinella spiralis

Background

Trichinella nematodes are globally distributed food-borne pathogens, in which Trichinella spiralis is the most common species in China. Microsatellites are a powerful tool in population genetics and phylogeographic analysis. However, only a few microsatellite markers were reported in T. spiralis. Thus, there is a need to develop and validate genome-wide microsatellite markers for T. spiralis.

Methods

Microsatellites were selected from shotgun genomic sequences using MIcroSAtellite identification tool (MISA). The identified markers were validated in 12 isolates of T. spiralis in China.

Results

A total of 93,140 microsatellites were identified by MISA from 9267 contigs in T. spiralis genome sequences, in which 16 polymorphic loci were selected for validation by PCR with single larvae from 12 isolates of T. spiralis in China. There were 7–19 alleles per locus (average 11.25 alleles per locus). The observed heterozygosity (H_O) and expected heterozygosity (H_E) ranged from 0.325 to 0.750 and 0.737 to 0.918, respectively. The polymorphism information content (PIC) ranged from 0.719 to 0.978 (average 0.826). Among the 16 loci, markers for 10 loci could be amplified from all 12 international standard strains of Trichinella spp.

Conclusions

Sixteen highly polymorphic markers were selected and validated for T. spiralis. Primary phylogenetic analysis showed that these markers might serve as a useful tool for genetic studies of Trichinella parasites.

Evaluation of Cross-Species Transferability of SSR Markers in Foeniculum vulgare

Fennel (Foeniculum vulgare) is a species belonging to the Apiaceae family, well known for its nutritional and pharmacological properties. Despite the economic and agricultural relevance, its genomic and transcriptomic data remain poor. Microsatellites-also known as simple sequence repeats (SSRs)-are codominant markers widely used to perform cross-amplification tests starting from markers developed in related species. SSRs represent a powerful tool, especially for those species lacking genomic information. In this study, a set of primers previously designed in Daucus carota for polymorphic SSR loci was tested in commercial varieties and breeding lines of fennel in order to: (i) test their cross-genera transferability, (ii) look at their efficiency in assessing genetic diversity, and (iii) identify their usefulness for marker-assisted selection (MAS) in breeding programs. Thirty-nine SSR markers from carrot were selected and tested for their transferability score, and only 23% of them resulted suitable for fennel. The low rate of SSR transferability between the two species evidences the difficulties of the use of genomic SSR in cross-genera transferability.

Genetic Differentiation and Clonal Expansion of Chinese Botrytis cinerea Populations from Tomato and Other Crops in China

Botrytis cinerea is an important pathogen of vegetable and fruit crops but little is known about its population structure and genetics in China. We hypothesized that the geographic populations of B. cinerea in China would be genetically differentiated by host, geographic location, and/or year. In this study, we collected 393 B. cinerea isolates representing 28 populations from tomato, cherry, and nectarine from 2006 to 2014 in China. The isolates were analyzed using 14 microsatellite markers, including six new markers that provided more genotyping power than the eight previously published loci. We also investigated the B. cinerea population structure and inferred its mode of reproduction and dispersal based on genotype data. High genotypic diversity was detected in all populations, and clonal reproduction was dominant. Southern China populations harbored more genotypes than northern populations. Differentiation by host plant was evident. Between 2011 and 2012, genotypes changed only slightly among years for Liaoning populations, but they changed substantially among years for the Shanghai and Fujian populations. Clonal dispersal was detected and the farthest dispersal distance was estimated to be about 1,717 km. Two high-frequency genotypes were widely distributed in more than 10 populations and across several years. Our results provide useful, novel information for plant breeding programs and control of B. cinerea in China.

Molecular genetic diversity and differentiation of Nile tilapia (Oreochromis niloticus, L. 1758) in East African natural and stocked populations

BACKGROUND

The need for enhancing the productivity of fisheries in Africa triggered the introduction of non-native fish, causing dramatic changes to local species. In East Africa, the extensive translocation of Nile tilapia (Oreochromis niloticus) is one of the major factors in this respect. Using 40 microsatellite loci with SSR-GBS techniques, we amplified a total of 664 individuals to investigate the genetic structure of O. niloticus from East Africa in comparison to Ethiopian and Burkina Faso populations.

RESULTS

All three African regions were characterized by independent gene-pools, however, the Ethiopian population from Lake Tana was genetically more divergent (Fst = 2.1) than expected suggesting that it might be a different sub-species. In East Africa, the genetic structure was congruent with both geographical location and anthropogenic activities (Isolation By Distance for East Africa, R2 = 0.67 and Uganda, R2 = 0.24). O. niloticus from Lake Turkana (Kenya) was isolated, while in Uganda, despite populations being rather similar to each other, two main natural catchments were able to be defined. We show that these two groups contributed to the gene-pool of different non-native populations. Moreover, admixture and possible hybridization with other tilapiine species may have contributed to the genetic divergence found in some populations such as Lake Victoria. We detected other factors that might be affecting Nile tilapia genetic variation. For example, most of the populations have gone through a reduction in genetic diversity, which can be a consequence of bottleneck (G-W, < 0.5) caused by overfishing, genetic erosion due to fragmentation or founder effect resulting from stocking activities.

CONCLUSIONS

The anthropogenic activities particularly in the East African O. niloticus translocations, promoted artificial admixture among Nile Tilapia populations. Translocations may also have triggered hybridization with the native congenerics, which needs to be further studied. These events may contribute to outbreeding depression and hence compromising the sustainability of the species in the region.

Polymorphic markers for identification of parasite population in Plasmodium malariae

Background

Molecular genotyping in Plasmodium serves many aims including providing tools for studying parasite population genetics and distinguishing recrudescence from reinfection. Microsatellite typing, insertion-deletion (INDEL) and single nucleotide polymorphisms is used for genotyping, but only limited information is available for Plasmodium malariae, an important human malaria species. This study aimed to provide a set of genetic markers to facilitate the study of P. malariae population genetics.

Methods

Markers for microsatellite genotyping and pmmsp1 gene polymorphisms were developed and validated in symptomatic P. malariae field isolates from Myanmar (N = 37). Fragment analysis was used to determine allele sizes at each locus to calculate multiplicity of infections (MOI), linkage disequilibrium, heterozygosity and construct dendrograms. Nucleotide diversity (π), number of haplotypes, and genetic diversity (H_d) were assessed and a phylogenetic tree was constructed. Genome-wide microsatellite maps with annotated regions of newly identified markers were constructed.

Results

Six microsatellite markers were developed and tested in 37 P. malariae isolates which showed sufficient heterozygosity (0.530–0.922), and absence of linkage disequilibrium (I_A^S=0.03, p value > 0.05) (N = 37). In addition, a tandem repeat (VNTR)-based pmmsp1 INDEL polymorphisms marker was developed and assessed in 27 P. malariae isolates showing a nucleotide diversity of 0.0976, haplotype gene diversity of 0.698 and identified 14 unique variants. The size of VNTR consensus repeat unit adopted as allele was 27 base pairs. The markers Pm12_426 and pmmsp1 showed greatest diversity with heterozygosity scores of 0.920 and 0.835, respectively. Using six microsatellites markers, the likelihood that any two parasite strains would have the same microsatellite genotypes was 8.46 × 10⁻⁴ and was further reduced to 1.66 × 10⁻⁴ when pmmsp1 polymorphisms were included.

Conclusions

Six novel microsatellites genotyping markers and a set of pmmsp1 VNTR-based INDEL polymorphisms markers for P. malariae were developed and validated. Each marker could be independently or in combination employed to access genotyping of the parasite. The newly developed markers may serve as a useful tool for investigating parasite diversity, population genetics, molecular epidemiology and for distinguishing recrudescence from reinfection in drug efficacy studies.

Genome- and transcriptome-derived microsatellite loci in lumpfish Cyclopterus lumpus: molecular tools for aquaculture, conservation and fisheries management

The lumpfish Cyclopterus lumpus is commercially exploited in numerous areas of its range in the North Atlantic Ocean, and is important in salmonid aquaculture as a biological agent for controlling sea lice. Despite the economic importance, few genetic resources for downstream applications, such as linkage mapping, parentage analysis, marker-assisted selection (MAS), quantitative trait loci (QTL) analysis, and assessing adaptive genetic diversity are currently available for the species. Here, we identify both genome- and transcriptome-derived microsatellites loci from C. lumpus to facilitate such applications. Across 2,346 genomic contigs, we detected a total of 3,067 microsatellite loci, of which 723 were the most suitable ones for primer design. From 116,555 transcriptomic unigenes, we identified a total of 231,556 microsatellite loci, which may indicate a high coverage of the available STRs. Out of these, primer pairs could only be designed for 6,203 loci. Dinucleotide repeats accounted for 89 percent and 52 percent of the genome- and transcriptome-derived microsatellites, respectively. The genetic composition of the dominant repeat motif types showed differences from other investigated fish species. In the genome-derived microsatellites AC/GT (67.8 percent), followed by AG/CT (15.1 percent) and AT/AT (5.6 percent) were the major motifs. Transcriptome-derived microsatellites showed also most dominantly the AC/GT repeat motif (33 percent), followed by A/T (26.6 percent) and AG/CT (11 percent). Functional annotation of microsatellite-containing transcriptomic sequences showed that the majority of the expressed sequence tags encode proteins involved in cellular and metabolic processes, binding activity and catalytic reactions. Importantly, STRs linked to genes involved in immune system process, growth, locomotion and reproduction were discovered in the present study. The extensive genomic marker information reported here will facilitate molecular ecology studies, conservation initiatives and will benefit many aspects of the breeding programmes of C. lumpus.

Microsatellite loci development for three catfish species from northwestern South America

ABSTRACT The Neotropical catfish species Ageneiosus pardalis, Pimelodus grosskopfii, and Sorubim cuspicaudus are important fishery resources in Colombia that show historical declines in their capture. This study used next-generation sequencing with 454 FLX technology (Roche Applied Science) and bioinformatics analysis to develop between 18 and 24 microsatellite loci for these species. The novel microsatellite loci showed high values of polymorphic information content -PIC (A. pardalis: 0.601-0.903, P. grosskopfii: 0.748-0.946 and S. cuspicaudus: 0.383-0.876), and the average number of alleles/locus ranged from 7-15 for A. pardalis, 9-30 for P. grosskopfii and 5-14 for S. cuspicaudus. The average observed and expected heterozygosities were respectively, 0.757 ± 0.035 and 0.834 ± 0.015 for A. pardalis; 0.596 ± 0.040 and 0.881 ± 0.009 for P. grosskopfii; and 0.747 ± 0.031 and 0.757 ± 0.025 for S. cuspicaudus. For future studies, these loci can be useful to estimate the genetic diversity and population structure in these three Neotropical catfishes.

Self-(In)compatibility Systems: Target Traits for Crop-Production, Plant Breeding, and Biotechnology

Self-incompatibility (SI) mechanisms prevent self-fertilization in flowering plants based on specific discrimination between self- and non-self pollen. Since this trait promotes outcrossing and avoids inbreeding it is a widespread mechanism of controlling sexual plant reproduction. Growers and breeders have effectively exploited SI as a tool for manipulating domesticated crops for thousands of years. However, only within the past thirty years have studies begun to elucidate the underlying molecular features of SI. The specific S-determinants and some modifier factors controlling SI have been identified in the sporophytic system exhibited by Brassica species and in the two very distinct gametophytic systems present in Papaveraceae on one side and in Solanaceae, Rosaceae, and Plantaginaceae on the other. Molecular level studies have enabled SI to SC transitions (and vice versa) to be intentionally manipulated using marker assisted breeding and targeted approaches based on transgene integration, silencing, and more recently CRISPR knock-out of SI-related factors. These scientific advances have, in turn, provided a solid basis to implement new crop production and plant breeding practices. Applications of self-(in)compatibility include widely differing objectives such as crop yield and quality improvement, marker-assisted breeding through SI genotyping, and development of hybrids for overcoming intra- and interspecific reproductive barriers. Here, we review scientific progress as well as patented applications of SI, and also highlight future prospects including further elucidation of SI systems, deepening our understanding of SI-environment relationships, and new perspectives on plant self/non-self recognition.

Seventeen new microsatellites for Tamarix gallica and cross-amplification in Tamarix species

PREMISE

Microsatellite markers were developed for the western Mediterranean tree Tamarix gallica (Tamaricaceae) as part of a study of its genetic diversity and structure.

METHODS AND RESULTS

Seventeen microsatellite markers were developed for T. gallica, 14 of which were polymorphic. These microsatellites have di-, tri-, and tetranucleotide repeats with 1-13 alleles per locus and population. Levels of observed and expected heterozygosity ranged from 0.000 to 0.900 and from 0.000 to 0.863, respectively. Six microsatellites showed significant deviations from Hardy-Weinberg equilibrium in at least one population. Cross-amplification in 19 Tamarix species showed a wide transferability to other species of the genus.

CONCLUSIONS

The 14 new polymorphic microsatellite markers will be used to assess the genetic diversity and population genetic structure of T. gallica. Additionally, the successful cross-species amplification suggests their potential usefulness for investigating species delimitation and population genetics in the genus Tamarix.

The roles of vicariance and dispersal in the differentiation of two species of the Rhinella marina species complex

The high levels of Neotropical biodiversity are commonly associated with the intense Neogene-Quaternary geological events and climate dynamics. Here, we investigate the evolutionary history of two species of Neotropical closely related amphibians (R. horribilis and R. marina). We combine published data with new mitochondrial DNA sequences and multiple nuclear markers, including 12 microsatellites. The phylogenetic analyses showed support for grouping the samples in two main clades; R. horribilis (Central America and Mexico) and R. marina (South America east of the Andes). However, the phylogenetic inferences also show an evident mito-nuclear discordance. We use Approximate Bayesian Computation (ABC) to test the role of different events in the diversification between the two groups recovered. We found that both species were affected primarily by a recent Pleistocene divergence, which was similar to the divergence estimate revealed by the Isolation-with-Migration model, under persistent bidirectional gene flow through time. We provide the first evidence that R. horribilis is differentiated from the South American R. marina at the nuclear level supporting the taxonomic status of R. horribilis, which has been controversial for more than a century.

New EST-SSR Markers for Individual Genotyping of Opium Poppy Cultivars (Papaver somniferum L.)

High-quality simple sequence repeat (SSR) markers are invaluable tools for revealing genetic variability which could be utilized for many purposes, such as breeding new varieties or the identifying current ones, among other applications. Based on the analysis of 3.7 million EST sequences and 15 genomic sequences from bacterial artificial chromosome (BAC) libraries, 200 trinucleotide genic (EST)-SSR and three genomic (gSSR) markers were tested, where 17 of them fulfilled all criteria for quality markers. Moreover, the reproducibility of these new markers was verified by two genetics laboratories, with a mean error rate per allele and per locus equal to 0.17%. These markers were tested on 38 accessions of Papaver somniferum and nine accessions of another five species of the Papaver and Argemone genera. In total, 118 alleles were detected for all accessions (median = 7; three to ten alleles per locus) and 88 alleles (median = 5; three to nine alleles per locus) within P. somniferum alone. Multivariate methods and identity analysis revealed high resolution capabilities of the new markers, where all but three pair accessions (41 out of 47) had a unique profile and opium poppy was distinguished from other species.

Genetic diversity and population genetic structure of six dromedary camel (camelus dromedarius) populations in Saudi Arabia

Camels are an integral and essential component of the Saudi Arabian heritage. The genetic diversity and population genetic structure of dromedary camels are poorly documented in Saudi Arabia so this study was carried out to investigate the genetic diversity of both local and exotic camel breeds. The genetic diversity was evaluated within and among camel populations using 21 microsatellite loci. Hair and blood samples were collected from 296 unrelated animals representing 4 different local breeds, namely Majaheem (MG), Maghateer (MJ), Sofr (SO), and Shaul (SH), and two exotic breeds namely Sawahli (SL) and Somali (SU). Nineteen out of 21 microsatellite loci generated multi-locus fingerprints for the studied camel individuals, with an average of 13.3 alleles per locus. Based on the genetic analyses, the camels were divided into two groups: one contained the Saudi indigenous populations (MG, MJ, SH and SO) and the other contained the non-Saudi ones (SU and SL). There was very little gene flow occurring between the two groups. The African origin of SU and SL breeds may explain their close genetic relationship. It is anticipated that the genetic diversity assessment is important to preserve local camel genetic resources and develop future breeding programs to improve camel productivity.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves, Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early-breeding house sparrows and hybrids and rural late-breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z-chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land-use changes in a mosaic landscape.

Complete chloroplast genome sequence and phylogenetic analysis of wasabi (Eutrema japonicum) and its relatives

In Japan, two Eutrema species, wasabi (Eutrema japonicum, the important traditional Japanese condiment) and yuriwasabi (E. tenue), have been recognized as endemic species. We sequenced complete chloroplast (cp) genomes of seven wasabi and yuriwasabi accessions from Japan to study their phylogeny and evolution, using molecular dating of species divergence. Phylogenetic analyses of the complete cp DNA of these two Japanese species and five other Eurasian Eutrema species revealed that wasabi and yuriwasabi did not form a monophyletic group. One yuriwasabi accession (Gifu) formed a clade with E. yunnanense from China, indicating that this accession should be considered as a different species from the other yuriwasabi accessions. We reveal that Japanese Eutrema species diverged from the ‘E. yunnanense–yuriwasabi (Gifu)’ clade approximately 1.3 million years ago (Mya), suggesting that the connection between Japan and the Eurasian continent has existed more recently than the Quaternary period. The abundance of cp sequence data in this study also allowed the detection of genetic differentiation among wasabi cultivars. The two polymorphic sites detected between ‘Fujidaruma’ and ‘Shimane No.3’ were used to develop genotyping markers. The cp genome information provided here will thus inform the evolutionary histories of Japanese Eutrema species and help in genotyping wasabi cultivars.

Genetically Independent Tetranucleotide to Hexanucleotide Core Motif SSR Markers for Identifying Lentinula edodes Cultivars

For the purpose of protecting the rights of Lentinula edodes breeders, we developed a new simple sequence repeat (SSR) marker set consisting only of genetically independent tetranucleotide or longer core motifs. Using available genome sequences for five L. edodes strains, we designed primers for 13 SSR markers that amplified polymorphic sequences in 20 L. edodes cultivars. We evaluated the independence of every possible marker pair based on genotype data. Consequently, eight genetically independent markers were selected. The polymorphic information content values of the markers ranged from 0.269 to 0.764, with an average of 0.409. The markers could distinguish among 20 L. edodes cultivars and produced highly repeatable and reproducible results. The markers developed in this study will enable the precise identification of L. edodes cultivars, and may be useful for protecting breeders' rights.

A high-density exome capture genotype-by-sequencing panel for forestry breeding in Pinus radiata

Development of genome-wide resources for application in genomic selection or genome-wide association studies, in the absence of full reference genomes, present a challenge to the forestry industry, where longer breeding cycles could benefit from the accelerated selection possible through marker-based breeding value predictions. In particular, large conifer megagenomes require a strategy to reduce complexity, whilst ensuring genome-wide coverage is achieved. Using a transcriptome-based reference template, we have successfully developed a high density exome capture genotype-by-sequencing panel for radiata pine (Pinus radiata D.Don), capable of capturing in excess of 80,000 single nucleotide polymorphism (SNP) markers with a minor allele frequency above 0.03 in the population tested. This represents approximately 29,000 gene models from a core set of 48,914 probes. A set of 704 SNP markers capable of pedigree reconstruction and differentiating individual genotypes were tested within two full-sib mapping populations. While as few as 70 markers could reconstruct parentage in almost all cases, the impact of missing genotypes was noticeable in several offspring. Therefore, 60 sets of 110 randomly selected SNP markers were compared for both parentage reconstruction and clone differentiation. The performance in parentage reconstruction showed little variation over 60 iterations. However, there was notable variation in discriminatory power between closely related individuals, indicating a higher density SNP marker panel may be required to elucidate hidden relationships in complex pedigrees.

Isolation and characterisation of 14 novel microsatellite markers through Next Generation Sequencing for the commercial Atlantic seabob shrimp Xiphopenaeus kroyeri

Assessing population genetic structure is a crucial step to support fisheries and conservation management. DNA microsatellite molecular markers are a widely used tool in population genotyping. In the present study, we characterised and developed 14 novel polymorphic microsatellite markers for a decapod crustacean, the Atlantic seabob shrimp Xiphopenaeus kroyeri (Heller, 1862), through rapid and cost-effective Illumina shotgun sequencing and a Galaxy-based bioinformatic pipeline. We genotyped 60 individuals from 2 populations with the newly developed microsatellites, resulting in the detection of 3 to 29 alleles per locus. Four loci deviated from Hardy-Weinberg equilibrium. Cross-amplification in a cryptic congeneric species was successful for eight loci (57%). The microsatellite loci developed in this study will be highly relevant for genetic and evolutionary studies of X. kroyeri, and for the stock management of this commercially exploited species.

The Tsushima leopard cat exhibits extremely low genetic diversity compared with the Korean Amur leopard cat: Implications for conservation

We examined genetic diversity of the wild Tsushima leopard cat-a regional population of the Amur leopard cat-using microsatellite markers. In addition, we compared genetic diversity of the Tsushima leopard cat with that of the Korean population of Amur leopard cat. Although bias should be considered when applying cross-species amplification, the Tsushima leopard cat showed a lower index of molecular genetic diversity than did the Korean population. These results were consistent with those obtained using other genetic markers, such as mitochondrial DNA and Y chromosome sequences. This low genetic diversity of the wild Tsushima leopard cat may be derived from the founding population. Furthermore, our results suggest that the captive populations held in Japanese zoos may show extremely low genetic diversity, leading to difficulties in genetic management of the Tsushima leopard cat. Moreover, the two regional populations were clearly separated using these marker sets. In the present study, we demonstrated that the genetic diversity of the Tsushima leopard cat is extremely low compared with that of the continental regional population. Importantly, the Japanese captive population for ex situ conservation was derived from a founding population with extremely low genetic diversity; hence, we assume that both the captive and wild populations showed extremely low genetic diversities. Our findings emphasize the need to develop carefully considered management strategies for genetic conservation.

Human mediated translocation of Pacific paper mulberry [Broussonetia papyrifera (L.) L'Hér. ex Vent. (Moraceae)]: Genetic evidence of dispersal routes in Remote Oceania

Paper mulberry, Broussonetia papyrifera (L.) L’Hér. ex Vent. (Moraceae), a dioecious species, was transported by humans from Taiwan to the islands of Remote Oceania. Its introduction and cultivation in Remote Oceania was intentional due to its cultural importance as a fiber source for barkcloth textiles. The aim of this study was to explore the genetic diversity and structure of paper mulberry populations within Remote Oceania in order to infer dispersal patterns that may reflect past human interaction among island groups. We present the integrated analysis of 380 samples (313 contemporary and 67 herbarium specimens) collected in Near and Remote Oceania. Genetic characterization was based on a set of ten microsatellites developed for B. papyrifera and complemented with the analysis of the ribosomal internal transcribed spacer ITS-1 sequence, a sex marker and the chloroplast ndhF–rpl32 intergenic spacer. Microsatellite data identify a total of 64 genotypes, despite this being a clonally propagated crop, and show three major dispersal hubs within Remote Oceania, centered on the islands of Fiji, Tonga, and Pitcairn. Of 64 genotypes identified, 55 correspond to genotypes associated to female-sexed plants that probably descend from plants introduced by the prehistoric Austronesian-speaking voyagers. The ratio of accessions to genotypes between herbarium and contemporary samples, suggests recent loss of genetic diversity. In addition to the chloroplast haplotypes described previously, we detected two new haplotypes within Remote Oceania both originating in Taiwan. This is the first study of a commensal species to show genetic structuring within Remote Oceania. In spite of the genetic bottleneck, the presence of only one sex, a timespan of less than 5000 years, and asexual propagation of this crop in Remote Oceania, we detect genetic diversity and regional structuring. These observations suggest specific migration routes between island groups within Remote Oceania.

Identification and evaluation of a core microsatellite panel for use in white-tailed deer (Odocoileus virginianus)

Background

Microsatellite loci have been used extensively over the past two decades to study the genetic characteristics of non-model species. The ease of microsatellite development and ability to adapt markers from related species has led to the proliferation of available markers for many commonly studied species. Because it is often infeasible to genotype individuals across all available loci, researchers generally rely on subsets of markers. Marker choice can bias inferences made using disparate suites of loci. This has been a primary motivation for efforts to identify uniform marker panels. Here, we use the geographic distribution of previous studies to identify microsatellite loci for white-tailed deer (Odocoileus virginianus) with the potential for widespread use, and we evaluate the effectiveness of this panel in a portion of the range where few previous studies have been conducted. The purpose was to consolidate the numerous genetic resources for this species into a manageable panel and to provide a uniform methodology that improves comparisons between past and future studies.

Results

We reviewed microsatellite panels from 58 previous or ongoing projects and identified 106 candidate loci. We developed a multiplex protocol and evaluated the efficacy of 17 of the most commonly used loci using 720 DNA samples collected from the Mid-Atlantic region of the United States of America. Amplification errors were detected in six of these loci. The 11 remaining loci were highly polymorphic, exhibited low frequencies of null alleles, and were easy to interpret with the aid of allele binning software.

Conclusions

The development of broadly-applicable, core microsatellite panels has the potential to improve repeatability and comparative ability for commonly studied species. The properties of the consolidated 11 microsatellite panel suggest that they are applicable for many common research objectives for white-tailed deer. The geographic distribution of previous studies using these markers provides a greater degree of confidence regarding the robustness to common sources of error related to amplification anomalies, such as null alleles, relative to loci with more limited use. While this does not replace further evaluation of genotyping errors, it does provide a common platform that benefits future research studies.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0750-z) contains supplementary material, which is available to authorized users.

Markers in time and space: A review of the last decade of plant phylogeographic approaches

Plant studies comprise a relatively small proportion of the phylogeographic literature, likely as a consequence of the fundamental challenges posed by the complex genomic structures and life history strategies of these organisms. Comparative plastomics (i.e., comparisons of mutation rates within and among regions of the chloroplast genome) across plant lineages has led to an increased understanding of which markers are likely to provide the most information at low taxonomic levels. However, the extent to which the results of such work have influenced the literature has not been fully assessed, nor has the extent to which plant phylogeographers explicitly analyse markers in time and space, both of which are integral components of the field. Here, we reviewed more than 400 publications from the last decade of plant phylogeography to specifically address the following questions: (i) What is the phylogenetic breadth of studies to date? (ii) What molecular markers have been used, and why were they chosen? (iii) What kinds of markers are most frequently used and in what combinations? (iv) How frequently are divergence time estimation and ecological niche modelling used in plant phylogeography? Our results indicate that chloroplast DNA sequence data remain the primary tool of choice, followed distantly by nuclear DNA sequences and microsatellites. Less than half (42%) of all studies use divergence time estimation, while even fewer use ecological niche modelling (14%). We discuss the implications of our findings, as well as the need for community standards on data reporting.

Application of high-throughput amplicon sequencing-based SSR genotyping in genetic background screening

Background

Host genetic backgrounds affect gene functions. The genetic backgrounds of genetically engineered organisms must be identified to confirm their genetic backgrounds identity with those of recipients. Marker-assisted backcrossing (MAB), transgenesis and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (CRISPR/Cas9) editing are three commonly used genetic engineering techniques. However, methods for genetic background screening between genetically engineered organisms and corresponding recipients suffer from low efficiency, low accuracy or high cost.

Results

Here, we improved our previously reported AmpSeq-SSR method, an amplicon sequencing-based simple sequence repeat (SSR) genotyping method, by selecting SSR loci with high polymorphism among varieties. Ultimately, a set of 396 SSRs was generated and applied to evaluate the genetic backgrounds identity between rice lines developed through MAB, transgenesis, and CRISPR/Cas9 editing and the respective recipient rice. We discovered that the percentage of different SSRs between the MAB-developed rice line and its recipient was as high as 23.5%. In contrast, only 0.8% of SSRs were different between the CRISPR/Cas9-system-mediated rice line and its recipient, while no SSRs showed different genotypes between the transgenic rice line and its recipient. Furthermore, most differential SSRs induced by MAB technology were located in non-coding regions (62.9%), followed by untranslated regions (21.0%) and coding regions (16.1%). Trinucleotide repeats were the most prevalent type of altered SSR. Most importantly, all altered SSRs located in coding regions were trinucleotide repeats.

Conclusions

This method is not only useful for the background evaluation of genetic resources but also expands our understanding of the unintended effects of different genetic engineering techniques. While the work we present focused on rice, this method can be readily extended to other organisms.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5800-4) contains supplementary material, which is available to authorized users.

Population Structure of Pythium ultimum from Greenhouse Floral Crops in Michigan

Pythium ultimum causes seedling damping-off and root and crown rot in greenhouse ornamental plants. To understand the population dynamics and assess population structure of P. ultimum in Michigan floriculture crops, simple sequence repeats (SSRs) were developed using the previously published P. ultimum predicted transcriptome. A total of 166 isolates sampled from 2011 to 2013 from five, one, and three greenhouses in Kalamazoo, Kent, and Wayne Counties, respectively, were analyzed using six polymorphic and fluorescently labeled SSR markers. The average unbiased Simpson's index (λ_u, 0.95), evenness (E₅, 0.56), and recovery of 12 major clones out of the 65 multilocus genotypes obtained, suggests that P. ultimum is not a recent introduction into Michigan greenhouses. Analyses revealed a clonal population, with limited differentiation among seasons, hosts, and counties sampled. Results also indicated the presence of common genotypes among years, suggesting that sanitation measures should be enhanced to eradicate resident P. ultimum populations. Finally, the presence of common genotypes among counties suggests that there is an exchange of infected plant material among greenhouse facilities, or that there is a common source of inoculum coming to the region. Continued monitoring of pathogen populations will enhance our understanding of population dynamics of P. ultimum in Michigan and facilitate improvement of control strategies.

A trait-based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life-history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species-level life-history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G'_ST , and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G'_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.

Quantitative trait loci mapping for the shear force value in breast muscle of F2 chickens

The shear force value is one of the major traits that determine meat quality. In the present study, we performed QTL analysis for chicken breast muscle shear force value at 7 wk of age using 545 single nucleotide polymorphism (SNP) markers developed via restriction-site associated DNA sequencing (RAD-seq). An F2 resource family was generated by mating Oh-Shamo, a native Japanese chicken breed, and the White Plymouth Rock chicken breed. A total of 215 F2 birds were produced. Simple interval mapping revealed one significant main-effect QTL between 6.28 and 8.10 Mb SNPs on the chromosome Z with a logarithm of odds score of 5.53 at the genome-wide 5% level. At this QTL, the confidence interval, phenotypic variance explained, and additive effect were 26 cM, 12.24%, and -0.31 in males and -0.34 in females, respectively. No QTL with epistatic interaction effects were detected. To our knowledge, this is the first report on a QTL affecting the shear force value in the chicken breast muscle, using SNP markers derived from RAD-seq.

No evidence for an intragenomic arms race under paternal genome elimination in Planococcus mealybugs

Genomic conflicts arising during reproduction might play an important role in shaping the striking diversity of reproductive strategies across life. Among these is paternal genome elimination (PGE), a form of haplodiploidy which has independently evolved several times in arthropods. PGE males are diploid but transmit maternally inherited chromosomes only, whereas paternal homologues are excluded from sperm. Mothers thereby effectively monopolize the parentage of sons, at the cost of the father's reproductive success. This creates striking conflict between the sexes that could result in a co-evolutionary arms race between paternal and maternal genomes over gene transmission, yet empirical evidence that such an arms race indeed takes place under PGE is scarce. This study addresses this by testing whether PGE is complete when paternal genotypes are exposed to divergent maternal backgrounds in intraspecific and hybrid crosses of the citrus mealybug, Planococcus citri, and the closely related Planococcus ficus. We determined whether males can transmit genetic information through their sons by tracking inheritance of two traits in a three-generation pedigree: microsatellite markers and sex-specific pheromone preferences. Our results suggest leakages of single paternal chromosomes through males occurring at a low frequency, but we find no evidence for transmission of paternal pheromone preferences from fathers to sons. The absence of differences between hybrid and intraspecific crosses in leakage rate of paternal alleles suggests that a co-evolutionary arms race cannot be demonstrated on this evolutionary timescale, but we conclude that there is scope for intragenomic conflict between parental genomes in mealybugs. Finally, we discuss how these paternal escapes can occur and what these findings may reveal about the evolutionary dynamics of this bizarre genetic system.

Polyploidy in Fruit Tree Crops of the Genus Annona (Annonaceae)

Genome duplication or polyploidy is one of the main factors of speciation in plants. It is especially frequent in hybrids and very valuable in many crops. The genus Annona belongs to the Annonaceae, a family that includes several fruit tree crops, such as cherimoya (Annona cherimola), sugar apple (Annona squamosa), their hybrid atemoya (A. cherimola × A. squamosa) or pawpaw (Asimina triloba). In this work, genome content was evaluated in several Annona species, A. triloba and atemoya. Surprisingly, while the hybrid atemoya has been reported as diploid, flow cytometry analysis of a progeny obtained from an interspecific cross between A. cherimola and A. squamosa showed an unusual ploidy variability that was also confirmed karyotype analysis. While the progeny from intraspecific crosses of A. cherimola showed polyploid genotypes that ranged from 2.5 to 33%, the hybrid atemoyas from the interspecific cross showed 35% of triploids from a total of 186 genotypes analyzed. With the aim of understanding the possible implications of the production of non-reduced gametes, pollen performance, pollen size and frequency distribution of pollen grains was quantified in the progeny of this cross and the parents. A large polymorphism in pollen grain size was found within the interspecific progeny with higher production of unreduced pollen in triploids (38%) than in diploids (29%). Moreover, using PCR amplification of selected microsatellite loci, while 13.7% of the pollen grains from the diploids showed two alleles, 41.28% of the grains from the triploids amplified two alleles and 5.63% showed up to three alleles. This suggests that the larger pollen grains could correspond to diploid and, in a lower frequency, to triploid pollen. Pollen performance was also affected with lower pollen germination in the hybrid triploids than in both diploid parents. The results confirm a higher percentage of polyploids in the interspecific cross, affecting pollen grain size and pollen performance. The occurrence of unreduced gametes in A. cherimola, A. squamosa and their interspecific progeny that may result in abnormalities of ploidy such as the triploids and tetraploids observed in this study, opens an interesting opportunity to study polyploidy in Annonaceae.

The future of parentage analysis: From microsatellites to SNPs and beyond

Parentage analysis is a cornerstone of molecular ecology that has delivered fundamental insights into behaviour, ecology and evolution. Microsatellite markers have long been the king of parentage, their hypervariable nature conferring sufficient power to correctly assign offspring to parents. However, microsatellite markers have seen a sharp decline in use with the rise of next-generation sequencing technologies, especially in the study of population genetics and local adaptation. The time is ripe to review the current state of parentage analysis and see how it stands to be affected by the emergence of next-generation sequencing approaches. We find that single nucleotide polymorphisms (SNPs), the typical next-generation sequencing marker, remain underutilized in parentage analysis but are gaining momentum, with 58 SNP-based parentage analyses published thus far. Many of these papers, particularly the earlier ones, compare the power of SNPs and microsatellites in a parentage context. In virtually every case, SNPs are at least as powerful as microsatellite markers. As few as 100-500 SNPs are sufficient to resolve parentage completely in most situations. We also provide an overview of the analytical programs that are commonly used and compatible with SNP data. As the next-generation parentage enterprise grows, a reliance on likelihood and Bayesian approaches, as opposed to strict exclusion, will become increasingly important. We discuss some of the caveats surrounding the use of next-generation sequencing data for parentage analysis and conclude that the future is bright for this important realm of molecular ecology.

Comparing RADseq and microsatellites for estimating genetic diversity and relatedness - Implications for brown trout conservation

The conservation and management of endangered species requires information on their genetic diversity, relatedness and population structure. The main genetic markers applied for these questions are microsatellites and single nucleotide polymorphisms (SNPs), the latter of which remain the more resource demanding approach in most cases. Here, we compare the performance of two approaches, SNPs obtained by restriction-site-associated DNA sequencing (RADseq) and 16 DNA microsatellite loci, for estimating genetic diversity, relatedness and genetic differentiation of three, small, geographically close wild brown trout (Salmo trutta) populations and a regionally used hatchery strain. The genetic differentiation, quantified as F ST, was similar when measured using 16 microsatellites and 4,876 SNPs. Based on both marker types, each brown trout population represented a distinct gene pool with a low level of interbreeding. Analysis of SNPs identified half- and full-siblings with a higher probability than the analysis based on microsatellites, and SNPs outperformed microsatellites in estimating individual-level multilocus heterozygosity. Overall, the results indicated that moderately polymorphic microsatellites and SNPs from RADseq agreed on estimates of population genetic structure in moderately diverged, small populations, but RADseq outperformed microsatellites for applications that required individual-level genotype information, such as quantifying relatedness and individual-level heterozygosity. The results can be applied to other small populations with low or moderate levels of genetic diversity.

An RNA Sequencing Transcriptome Analysis and Development of EST-SSR Markers in Chinese Hawthorn through Illumina Sequencing

Chinese hawthorn (Crataegus pinnatifida) is an important ornamental and economic horticultural plant. However, the lack of molecular markers has limited the development and utilization of hawthorn germplasm resources. Simple sequence repeats (SSRs) derived from expressed sequence tags (ESTs) allow precise and effective cultivar characterization and are routinely used for genetic diversity analysis. Thus, we first reported the development of polymorphic EST-SSR markers in C. pinnatifida with perfect repeats using Illumina RNA-Seq technique. In total, we investigated 14,364 unigenes, from which 5091 EST-SSR loci were mined. Di-nucleotides (2012, 39.52%) were the most abundant SSRs, followed by mono- (1989, 39.07%), and tri-nucleotides (1024, 20.11%). On the basis of these EST-SSRs, a total of 300 primer pairs were designed and used for polymorphism analysis in 70 accessions collected from different geographical regions of China. Of 239 (79.67%) pairs of primer-generated amplification products, 163 (54.33%) pairs of primers showed polymorphism. Finally, 33 primers with high polymorphism were selected for genetic diversity analysis and tested on 70 individuals with low-cost fluorescence-labeled M13 primers using capillary electrophoresis genotyping platform. A total of 108 alleles were amplified by 33 SSR markers, with the number of alleles (Na) ranging from 2 to 14 per locus (mean: 4.939), and the effective number of alleles (Ne) ranging from 1.258 to 3.214 (mean: 2.221). The mean values of gene diversity (He), observed heterozygosity (Ho), and polymorphism information content (PIC) were 0.524 (range 0.205–0.689), 0.709 (range 0.132–1.000), and 0.450 (range 0.184–0.642), respectively. Furthermore, the dendrogram constructed based on the EST-SSR separated the cultivars into two main clusters. In sum, our study was the first comprehensive study on the development and analysis of a large set of SSR markers in hawthorn. The results suggested that the use of NGS techniques for SSR development represented a powerful tool for genetic studies. Additionally, fluorescence-labeled M13 markers proved to be a valuable method for genotyping. All of these EST-SSR markers have agronomic potential and constitute a scientific basis for future studies on the identification, classification, and innovation of hawthorn germplasms.

Insights into evolving global populations of Phytophthora infestans via new complementary mtDNA haplotype markers and nuclear SSRs

In many parts of the world the damaging potato late blight pathogen, Phytophthora infestans, is spread as a succession of clonal lineages. The discrimination of genetic diversity within such evolving populations provides insights into the processes generating novel lineages and the pathways and drivers of pathogen evolution and dissemination at local and global scales. This knowledge, in turn, helps optimise management practices. Here we combine two key methods for dissecting mitochondrial and nuclear diversity and resolve intra and inter-lineage diversity of over 100 P. infestans isolates representative of key clonal lineages found globally. A novel set of PCR primers that amplify five target regions are provided for mitochondrial DNA sequence analysis. These five loci increased the number of mtDNA haplotypes resolved from four with the PCR RFLP method to 37 (17, 6, 8 and 4 for Ia, Ib, IIa, and IIb haplotypes, respectively, plus 2 Herb-1 haplotypes). As with the PCR RFLP method, two main lineages, I and II were defined. Group I contained 25 mtDNA haplotypes that grouped broadly according to the Ia and Ib types and resolved several sub-clades amongst the global sample. Group II comprised two distinct clusters with four haplotypes corresponding to the RFLP type IIb and eight haplotypes resolved within type IIa. The 12-plex SSR assay revealed 90 multilocus genotypes providing accurate discrimination of dominant clonal lineages and other genetically diverse isolates. Some association of genetic diversity and geographic region of contemporary isolates was observed; US and Mexican isolates formed a loose grouping, distinct from isolates from Europe, South America and other regions. Diversity within clonal lineages was observed that varied according to the age of the clone. In combination, these fine-scale nuclear and maternally inherited mitochondrial markers enabled a greater level of discrimination among isolates than previously available and provided complementary perspectives on evolutionary questions relating to the diversity, phylogeography and the origins and spread of clonal lineages of P. infestans.

Microsatellite development from genome skimming and transcriptome sequencing: comparison of strategies and lessons from frog species

Background

Even though microsatellite loci frequently have been isolated using recently developed next-generation sequencing (NGS) techniques, this task is still difficult because of the subsequent polymorphism screening requires a substantial amount of time. Selecting appropriate polymorphic microsatellites is a critical issue for ecological and evolutionary studies. However, the extent to which assembly strategy, read length, sequencing depth, and library layout produce a measurable effect on microsatellite marker development remains unclear. Here, we use six frog species for genome skimming and two frog species for transcriptome sequencing to develop microsatellite markers, and investigate the effect of different isolation strategies on the yield of microsatellites.

Results

The results revealed that the number of isolated microsatellites increases with increased data quantity and read length. Assembly strategy could influence the yield and the polymorphism of microsatellite development. Larger k-mer sizes produced fewer total number of microsatellite loci, but these loci had a longer repeat length, suggesting greater polymorphism. However, the proportion of each type of nucleotide repeats was not affected; dinucleotide repeats were always the dominant type. Finally, the transcriptomic microsatellites displayed lower levels of polymorphisms and were less abundant than genomic microsatellites, but more likely to be functionally linked loci.

Conclusions

These observations provide deep insight into the evolution and distribution of microsatellites and how different isolation strategies affect microsatellite development using NGS.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5329-y) contains supplementary material, which is available to authorized users.

Advances in genotyping microsatellite markers through sequencing and consequences of scoring methods for Ceratonia siliqua (Leguminosae)

PREMISE OF THE STUDY

Simple sequence repeat (SSR) or microsatellite markers have been used in a broad range of studies mostly scoring alleles on the basis of amplicon size as a proxy for the number of repeat units of an SSR motif. However, additional sources of variation within the SSR or in the flanking regions have largely remained undetected.

METHODS

In this study, we implemented a next-generation sequencing-based genotyping approach in a newly characterized set of 18 nuclear SSR markers for the carob tree, Ceratonia siliqua. Our aim was to evaluate the effect of three different methods of scoring molecular variation present within microsatellite markers on the genetic diversity and structure results.

RESULTS

The analysis of the sequences of 77 multilocus genotypes from four populations revealed SSR variation and additional sources of polymorphism in 87% of the loci analyzed (42 single-nucleotide polymorphisms and five insertion/deletion polymorphisms), as well as divergent paralog copies in two loci. Ignoring sequence variation under standard amplicon size genotyping resulted in incorrect identification of 69% of the alleles, with important effects on the genetic diversity and structure estimates.

DISCUSSION

Next-generation sequencing allows the detection and scoring of SSRs, single-nucleotide polymorphisms, and insertion/deletion polymorphisms to increase the resolution of population genetic studies.

Isolation and characterization of 21 polymorphic microsatellite loci for the rockpool shrimp Palaemon elegans using Illumina MiSeq sequencing

The rockpool shrimp Palaemon elegans is considered an important crustacean species within the European coastline fauna. This species is experiencing an ongoing geographical expansion beyond its native distribution range due to unintentional human introductions. A better knowledge of the genetic diversity, geographic structure and connectivity of its populations is necessary. In the present study, microsatellite loci were isolated using the Illumina MiSeq platform. The microsatellite-enriched library sequencing produced 3.9 million raw reads. Reads were processed and primer pairs were designed for microsatellite sequences amplification. Ninety-six microsatellite loci were preliminary screened in individuals from Atlantic and Mediterranean localities. From them, 21 loci exhibited reliable polymorphism and were thoroughly characterized in 30 individuals from a Cantabrian locality (Spain). No linkage disequilibrium between pairs of loci was detected. Number of alleles per locus ranged from 2 to 12. Observed and expected heterozygosities ranged from 0.033 to 0.833 and from 0.033 to 0.869 respectively. No significant departure from the Hardy-Weinberg equilibrium was detected in most of loci. This is the first time that microsatellite markers have been developed for P. elegans. This characterized microsatellite suite provides new suitable tools for further analyses, facilitating the understanding of population genetics both in natural and introduced populations.